Genetic diversity fuels gene discovery for tobacco and alcohol use.

Tobacco and alcohol use are heritable behaviours associated with 15% and 5.3% of worldwide deaths, respectively, due largely to broad increased risk for disease and injury1-4. These substances are used across the globe, yet genome-wide association studies have focused largely on individuals of European ancestries5. Here we leveraged global genetic diversity across 3.4 million individuals from four major clines of global ancestry (approximately 21% non-European) to power the discovery and fine-mapping of genomic loci associated with tobacco and alcohol use, to inform function of these loci via ancestry-aware transcriptome-wide association studies, and to evaluate the genetic architecture and predictive power of polygenic risk within and across populations. We found that increases in sample size and genetic diversity improved locus identification and fine-mapping resolution, and that a large majority of the 3,823 associated variants (from 2,143 loci) showed consistent effect sizes across ancestry dimensions. However, polygenic risk scores developed in one ancestry performed poorly in others, highlighting the continued need to increase sample sizes of diverse ancestries to realize any potential benefit of polygenic prediction.

Sequencing of 53,831 diverse genomes from the NHLBI TOPMed Program.

The Trans-Omics for Precision Medicine (TOPMed) programme seeks to elucidate the genetic architecture and biology of heart, lung, blood and sleep disorders, with the ultimate goal of improving diagnosis, treatment and prevention of these diseases. The initial phases of the programme focused on whole-genome sequencing of individuals with rich phenotypic data and diverse backgrounds. Here we describe the TOPMed goals and design as well as the available resources and early insights obtained from the sequence data. The resources include a variant browser, a genotype imputation server, and genomic and phenotypic data that are available through dbGaP (Database of Genotypes and Phenotypes)1. In the first 53,831 TOPMed samples, we detected more than 400 million single-nucleotide and insertion or deletion variants after alignment with the reference genome. Additional previously undescribed variants were detected through assembly of unmapped reads and customized analysis in highly variable loci. Among the more than 400 million detected variants, 97% have frequencies of less than 1% and 46% are singletons that are present in only one individual (53% among unrelated individuals). These rare variants provide insights into mutational processes and recent human evolutionary history. The extensive catalogue of genetic variation in TOPMed studies provides unique opportunities for exploring the contributions of rare and noncoding sequence variants to phenotypic variation. Furthermore, combining TOPMed haplotypes with modern imputation methods improves the power and reach of genome-wide association studies to include variants down to a frequency of approximately 0.01%.

Rare variants in long non-coding RNAs are associated with blood lipid levels in the TOPMed whole-genome sequencing study.

Long non-coding RNAs (lncRNAs) are known to perform important regulatory functions in lipid metabolism. Large-scale whole-genome sequencing (WGS) studies and new statistical methods for variant set tests now provide an opportunity to assess more associations between rare variants in lncRNA genes and complex traits across the genome. In this study, we used high-coverage WGS from 66,329 participants of diverse ancestries with measurement of blood lipids and lipoproteins (LDL-C, HDL-C, TC, and TG) in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) program to investigate the role of lncRNAs in lipid variability. We aggregated rare variants for 165,375 lncRNA genes based on their genomic locations and conducted rare-variant aggregate association tests using the STAAR (variant-set test for association using annotation information) framework. We performed STAAR conditional analysis adjusting for common variants in known lipid GWAS loci and rare-coding variants in nearby protein-coding genes. Our analyses revealed 83 rare lncRNA variant sets significantly associated with blood lipid levels, all of which were located in known lipid GWAS loci (in a ±500-kb window of a Global Lipids Genetics Consortium index variant). Notably, 61 out of 83 signals (73%) were conditionally independent of common regulatory variation and rare protein-coding variation at the same loci. We replicated 34 out of 61 (56%) conditionally independent associations using the independent UK Biobank WGS data. Our results expand the genetic architecture of blood lipids to rare variants in lncRNAs.

Rare coding variants in 35 genes associate with circulating lipid levels-A multi-ancestry analysis of 170,000 exomes.

Large-scale gene sequencing studies for complex traits have the potential to identify causal genes with therapeutic implications. We performed gene-based association testing of blood lipid levels with rare (minor allele frequency < 1%) predicted damaging coding variation by using sequence data from >170,000 individuals from multiple ancestries: 97,493 European, 30,025 South Asian, 16,507 African, 16,440 Hispanic/Latino, 10,420 East Asian, and 1,182 Samoan. We identified 35 genes associated with circulating lipid levels; some of these genes have not been previously associated with lipid levels when using rare coding variation from population-based samples. We prioritize 32 genes in array-based genome-wide association study (GWAS) loci based on aggregations of rare coding variants; three (EVI5, SH2B3, and PLIN1) had no prior association of rare coding variants with lipid levels. Most of our associated genes showed evidence of association among multiple ancestries. Finally, we observed an enrichment of gene-based associations for low-density lipoprotein cholesterol drug target genes and for genes closest to GWAS index single-nucleotide polymorphisms (SNPs). Our results demonstrate that gene-based associations can be beneficial for drug target development and provide evidence that the gene closest to the array-based GWAS index SNP is often the functional gene for blood lipid levels.

A framework for detecting noncoding rare-variant associations of large-scale whole-genome sequencing studies.

Large-scale whole-genome sequencing studies have enabled analysis of noncoding rare-variant (RV) associations with complex human diseases and traits. Variant-set analysis is a powerful approach to study RV association. However, existing methods have limited ability in analyzing the noncoding genome. We propose a computationally efficient and robust noncoding RV association detection framework, STAARpipeline, to automatically annotate a whole-genome sequencing study and perform flexible noncoding RV association analysis, including gene-centric analysis and fixed window-based and dynamic window-based non-gene-centric analysis by incorporating variant functional annotations. In gene-centric analysis, STAARpipeline uses STAAR to group noncoding variants based on functional categories of genes and incorporate multiple functional annotations. In non-gene-centric analysis, STAARpipeline uses SCANG-STAAR to incorporate dynamic window sizes and multiple functional annotations. We apply STAARpipeline to identify noncoding RV sets associated with four lipid traits in 21,015 discovery samples from the Trans-Omics for Precision Medicine (TOPMed) program and replicate several of them in an additional 9,123 TOPMed samples. We also analyze five non-lipid TOPMed traits.

Whole genome sequence analysis of platelet traits in the NHLBI Trans-Omics for Precision Medicine (TOPMed) initiative.

Platelets play a key role in thrombosis and hemostasis. Platelet count (PLT) and mean platelet volume (MPV) are highly heritable quantitative traits, with hundreds of genetic signals previously identified, mostly in European ancestry populations. We here utilize whole genome sequencing (WGS) from NHLBI's Trans-Omics for Precision Medicine initiative (TOPMed) in a large multi-ethnic sample to further explore common and rare variation contributing to PLT (n = 61 200) and MPV (n = 23 485). We identified and replicated secondary signals at MPL (rs532784633) and PECAM1 (rs73345162), both more common in African ancestry populations. We also observed rare variation in Mendelian platelet-related disorder genes influencing variation in platelet traits in TOPMed cohorts (not enriched for blood disorders). For example, association of GP9 with lower PLT and higher MPV was partly driven by a pathogenic Bernard-Soulier syndrome variant (rs5030764, p.Asn61Ser), and the signals at TUBB1 and CD36 were partly driven by loss of function variants not annotated as pathogenic in ClinVar (rs199948010 and rs571975065). However, residual signal remained for these gene-based signals after adjusting for lead variants, suggesting that additional variants in Mendelian genes with impacts in general population cohorts remain to be identified. Gene-based signals were also identified at several genome-wide association study identified loci for genes not annotated for Mendelian platelet disorders (PTPRH, TET2, CHEK2), with somatic variation driving the result at TET2. These results highlight the value of WGS in populations of diverse genetic ancestry to identify novel regulatory and coding signals, even for well-studied traits like platelet traits.

Whole-genome sequencing in diverse subjects identifies genetic correlates of leukocyte traits: The NHLBI TOPMed program.

Many common and rare variants associated with hematologic traits have been discovered through imputation on large-scale reference panels. However, the majority of genome-wide association studies (GWASs) have been conducted in Europeans, and determining causal variants has proved challenging. We performed a GWAS of total leukocyte, neutrophil, lymphocyte, monocyte, eosinophil, and basophil counts generated from 109,563,748 variants in the autosomes and the X chromosome in the Trans-Omics for Precision Medicine (TOPMed) program, which included data from 61,802 individuals of diverse ancestry. We discovered and replicated 7 leukocyte trait associations, including (1) the association between a chromosome X, pseudo-autosomal region (PAR), noncoding variant located between cytokine receptor genes (CSF2RA and CLRF2) and lower eosinophil count; and (2) associations between single variants found predominantly among African Americans at the S1PR3 (9q22.1) and HBB (11p15.4) loci and monocyte and lymphocyte counts, respectively. We further provide evidence indicating that the newly discovered eosinophil-lowering chromosome X PAR variant might be associated with reduced susceptibility to common allergic diseases such as atopic dermatitis and asthma. Additionally, we found a burden of very rare FLT3 (13q12.2) variants associated with monocyte counts. Together, these results emphasize the utility of whole-genome sequencing in diverse samples in identifying associations missed by European-ancestry-driven GWASs.

Reproducibility and Repeatability of Intravoxel Incoherent Motion MRI Acquisition Methods in Liver.

BACKGROUND: Intravoxel incoherent motion (IVIM) diffusion weighted MRI (DWI) has potential for evaluating hepatic fibrosis but image acquisition technique influence on diffusion parameter estimation bears investigation. PURPOSE: To minimize variability and maximize repeatably in abdominal DWI in terms of IVIM parameter estimates. STUDY TYPE: Prospective test-retest and image quality comparison. SUBJECTS: Healthy volunteers (3F/7M, 29.9 ± 12.9 years) and Family Study subjects (18F/12M, 51.7 ± 16.7 years), without and with liver steatosis. FIELD STRENGTH/SEQUENCE: Abdominal single-shot echo-planar imaging (EPI) and simultaneous multi-slice (SMS) DWI sequences with respiratory triggering (RT), breath-holding (BH), and navigator echo (NE) at 3 Tesla. ASSESSMENT: SMS-BH, EPI-NE, and SMS-RT data from twice-scanned healthy volunteers were analyzed using 6 × b-values (0-800 s⋅mm-2 ) and lower (LO) and higher (HI) b-value ranges. Family Study subjects were scanned using SMS and standard EPI sequences. The biexponential IVIM model was used to estimate fast-diffusion coefficient (Df ), fraction of fast diffusion (f), and slow-diffusion coefficient (Ds ). Scan time, estimated signal-to-noise ratio (eSNR), eSNR per acquisition, and distortion ratio were compared. STATISTICAL TESTS: Coefficients of variation (CoV) and Bland Altman analyses were performed for test-retest repeatability. Interclass correlation coefficient (ICC) assessed interobserver agreement with P < 0.05 deemed significant. RESULTS: Within-subject CoVs among volunteers (N = 10) for f and Ds were lowest in EPI-NE-LO (11.6%) and SMS-RT-HI (11.1%). Inter-observer ICCs for f and Ds were highest for EPI-NE-LO (0.63) and SMS-RT-LO (0.76). Df could not be estimated for most subjects. Estimated eSNR (EPI = 21.9, SMS = 4.7) and eSNR time (EPI = 6.7, SMS = 16.6) were greater for SMS, with less distortion in the liver region (DR-PE: EPI = 23.6, SMS = 13.1). DATA CONCLUSION: Simultaneous multislice acquisitions had significantly less variability and higher ICCs of Ds , higher eSNR, less distortion, and reduced scan time compared to EPI. EVIDENCE LEVEL: 2 TECHNICAL EFFICACY: Stage 1.

Genome-wide significant risk loci for mood disorders in the Old Order Amish founder population.

Genome-wide association studies (GWAS) of mood disorders in large case-control cohorts have identified numerous risk loci, yet pathophysiological mechanisms remain elusive, primarily due to the very small effects of common variants. We sought to discover risk variants with larger effects by conducting a genome-wide association study of mood disorders in a founder population, the Old Order Amish (OOA, n = 1,672). Our analysis revealed four genome-wide significant risk loci, all of which were associated with >2-fold relative risk. Quantitative behavioral and neurocognitive assessments (n = 314) revealed effects of risk variants on sub-clinical depressive symptoms and information processing speed. Network analysis suggested that OOA-specific risk loci harbor novel risk-associated genes that interact with known neuropsychiatry-associated genes via gene interaction networks. Annotation of the variants at these risk loci revealed population-enriched, non-synonymous variants in two genes encoding neurodevelopmental transcription factors, CUX1 and CNOT1. Our findings provide insight into the genetic architecture of mood disorders and a substrate for mechanistic and clinical studies.

Genetic correlations between liver fat content, metabolic health, and adiposity distribution in the Fels Longitudinal Study.

BACKGROUND AND AIMS: Hepatic steatosis is known to be heritable, but its genetic basis is mostly uncharacterized. Steatosis is associated with metabolic and adiposity features; recent studies hypothesize that shared genetic effects between these traits could account for some of the unexplained heritability. This study aimed to quantify these genetic associations in a family-based sample of non-Hispanic white adults. METHODS AND RESULTS: 704 participants (18-95 years, 55.8% female) from the Fels Longitudinal Study with an MRI assessment of liver fat were included. Quantitative genetic analyses estimated the age- and sex-adjusted heritability of individual traits and the genetic correlations within trait pairs. Mean liver fat was 5.95% (SE = 0.23) and steatosis (liver fat >5.56%) was present in 29.8% of participants. Heritability (h2± SE) of steatosis was 0.72 ± 0.17 (p = 6.80e-6). All other traits including liver enzymes, fasting glucose, HOMA-IR, visceral and subcutaneous adipose tissue (VAT, SAT), body mass index, body fat percent, waist circumference, lipids and blood pressure were also heritable. Significant genetic correlations were found between liver fat and all traits except aspartate aminotransferase (AST), and among most trait pairs. Highest genetic correlations were between liver fat and HOMA-IR (0.85 ± 0.08, p = 1.73e-8), fasting glucose and ALT (0.89 ± 0.26, p = 6.68e-5), and HOMA-IR with: waist circumference (0.81 ± 0.12, p = 3.76e-6), body fat percent (0.78 ± 0.12 p = 2.42e-5) and VAT (0.73 ± 0.07, p = 6.37e-8). CONCLUSIONS: Common genes may exist between liver fat accumulation, metabolic features and adiposity phenotypes.

Admixture mapping of cognitive function in diverse Hispanic and Latino adults: Results from the Hispanic Community Health Study/Study of Latinos.

INTRODUCTION: We conducted admixture mapping and fine-mapping analyses to identify ancestry-of-origin loci influencing cognitive abilities. METHODS: We estimated the association of local ancestry intervals across the genome with five neurocognitive measures in 7140 diverse Hispanic and Latino adults (mean age 55 years). We prioritized genetic variants in associated loci and tested them for replication in four independent cohorts. RESULTS: We identified nine local ancestry-associated regions for the five neurocognitive measures. There was strong biological support for the observed associations to cognitive function at all loci and there was statistical evidence of independent replication at 4q12, 9p22.1, and 13q12.13. DISCUSSION: Our study identified multiple novel loci harboring genes implicated in cognitive functioning and dementia, and uncovered ancestry-relevant genetic variants. It adds to our understanding of the genetic architecture of cognitive function in Hispanic and Latino adults and demonstrates the power of admixture mapping to discover unique haplotypes influencing cognitive function, complementing genome-wide association studies. HIGHLIGHTS: We identified nine ancestry-of-origin chromosomal regions associated with five neurocognitive traits. In each associated region, we identified single nucleotide polymorphisms (SNPs) that explained, at least in part, the admixture signal and were tested for replication in independent samples of Black, non-Hispanic White, and Hispanic/Latino adults with the same or similar neurocognitive tests. Statistical evidence of independent replication of the prioritized SNPs was observed for three of the nine associations, at chr4q12, chr9p22.1, and chr13q12.13. At all loci, there was strong biological support for the observed associations to cognitive function and dementia, prioritizing genes such as KIT, implicated in autophagic clearance of neurotoxic proteins and on mast cell and microglial-mediated inflammation; SLC24A2, implicated in synaptic plasticity associated with learning and memory; and MTMR6, implicated in phosphoinositide lipids metabolism.

Rapid, Phase-free Detection of Long Identity-by-Descent Segments Enables Effective Relationship Classification.

Identity-by-descent (IBD) segments are a useful tool for applications ranging from demographic inference to relationship classification, but most detection methods rely on phasing information and therefore require substantial computation time. As genetic datasets grow, methods for inferring IBD segments that scale well will be critical. We developed IBIS, an IBD detector that locates long regions of allele sharing between unphased individuals, and benchmarked it with Refined IBD, GERMLINE, and TRUFFLE on 3,000 simulated individuals. Phasing these with Beagle 5 takes 4.3 CPU days, followed by either Refined IBD or GERMLINE segment detection in 2.9 or 1.1 h, respectively. By comparison, IBIS finishes in 6.8 min or 7.8 min with IBD2 functionality enabled: speedups of 805-946× including phasing time. TRUFFLE takes 2.6 h, corresponding to IBIS speedups of 20.2-23.3×. IBIS is also accurate, inferring ≥7 cM IBD segments at quality comparable to Refined IBD and GERMLINE. With these segments, IBIS classifies first through third degree relatives in real Mexican American samples at rates meeting or exceeding other methods tested and identifies fourth through sixth degree pairs at rates within 0.0%-2.0% of the top method. While allele frequency-based approaches that do not detect segments can infer relationship degrees faster than IBIS, the fastest are biased in admixed samples, with KING inferring 30.8% fewer fifth degree Mexican American relatives correctly compared with IBIS. Finally, we ran IBIS on chromosome 2 of the UK Biobank dataset and estimate its runtime on the autosomes to be 3.3 days parallelized across 128 cores.

Imaging local genetic influences on cortical folding.

Recent progress in deciphering mechanisms of human brain cortical folding leave unexplained whether spatially patterned genetic influences contribute to this folding. High-resolution in vivo brain MRI can be used to estimate genetic correlations (covariability due to shared genetic factors) in interregional cortical thickness, and biomechanical studies predict an influence of cortical thickness on folding patterns. However, progress has been hampered because shared genetic influences related to folding patterns likely operate at a scale that is much more local (<1 cm) than that addressed in prior imaging studies. Here, we develop methodological approaches to examine local genetic influences on cortical thickness and apply these methods to two large, independent samples. We find that such influences are markedly heterogeneous in strength, and in some cortical areas are notably stronger in specific orientations relative to gyri or sulci. The overall, phenotypic local correlation has a significant basis in shared genetic factors and is highly symmetric between left and right cortical hemispheres. Furthermore, the degree of local cortical folding relates systematically with the strength of local correlations, which tends to be higher in gyral crests and lower in sulcal fundi. The relationship between folding and local correlations is stronger in primary sensorimotor areas and weaker in association areas such as prefrontal cortex, consistent with reduced genetic constraints on the structural topology of association cortex. Collectively, our results suggest that patterned genetic influences on cortical thickness, measurable at the scale of in vivo MRI, may be a causal factor in the development of cortical folding.

Rare coding variants in RCN3 are associated with blood pressure.

BACKGROUND: While large genome-wide association studies have identified nearly one thousand loci associated with variation in blood pressure, rare variant identification is still a challenge. In family-based cohorts, genome-wide linkage scans have been successful in identifying rare genetic variants for blood pressure. This study aims to identify low frequency and rare genetic variants within previously reported linkage regions on chromosomes 1 and 19 in African American families from the Trans-Omics for Precision Medicine (TOPMed) program. Genetic association analyses weighted by linkage evidence were completed with whole genome sequencing data within and across TOPMed ancestral groups consisting of 60,388 individuals of European, African, East Asian, Hispanic, and Samoan ancestries. RESULTS: Associations of low frequency and rare variants in RCN3 and multiple other genes were observed for blood pressure traits in TOPMed samples. The association of low frequency and rare coding variants in RCN3 was further replicated in UK Biobank samples (N = 403,522), and reached genome-wide significance for diastolic blood pressure (p = 2.01 × 10- 7). CONCLUSIONS: Low frequency and rare variants in RCN3 contributes blood pressure variation. This study demonstrates that focusing association analyses in linkage regions greatly reduces multiple-testing burden and improves power to identify novel rare variants associated with blood pressure traits.

Lipidomic profiling in the Strong Heart Study identified American Indians at risk of chronic kidney disease.

Dyslipidemia associates with and usually precedes the onset of chronic kidney disease (CKD), but a comprehensive assessment of molecular lipid species associated with risk of CKD is lacking. Here, we sought to identify fasting plasma lipids associated with risk of CKD among American Indians in the Strong Heart Family Study, a large-scale community-dwelling of individuals, followed by replication in Mexican Americans from the San Antonio Family Heart Study and Caucasians from the Australian Diabetes, Obesity and Lifestyle Study. We also performed repeated measurement analysis to examine the temporal relationship between the change in the lipidome and change in kidney function between baseline and follow-up of about five years apart. Network analysis was conducted to identify differential lipid classes associated with risk of CKD. In the discovery cohort, we found that higher baseline level of multiple lipid species, including glycerophospholipids, glycerolipids and sphingolipids, was significantly associated with increased risk of CKD, independent of age, sex, body mass index, diabetes and hypertension. Many lipid species were replicated in at least one external cohort at the individual lipid species and/or the class level. Longitudinal change in the plasma lipidome was significantly associated with change in the estimated glomerular filtration rate after adjusting for covariates, baseline lipids and the baseline rate. Network analysis identified distinct lipidomic signatures differentiating high from low-risk groups. Thus, our results demonstrated that disturbed lipid metabolism precedes the onset of CKD. These findings shed light on the mechanisms linking dyslipidemia to CKD and provide potential novel biomarkers for identifying individuals with early impaired kidney function at preclinical stages.

Impact of Rare and Common Genetic Variants on Diabetes Diagnosis by Hemoglobin A1c in Multi-Ancestry Cohorts: The Trans-Omics for Precision Medicine Program.

Hemoglobin A1c (HbA1c) is widely used to diagnose diabetes and assess glycemic control in individuals with diabetes. However, nonglycemic determinants, including genetic variation, may influence how accurately HbA1c reflects underlying glycemia. Analyzing the NHLBI Trans-Omics for Precision Medicine (TOPMed) sequence data in 10,338 individuals from five studies and four ancestries (6,158 Europeans, 3,123 African-Americans, 650 Hispanics, and 407 East Asians), we confirmed five regions associated with HbA1c (GCK in Europeans and African-Americans, HK1 in Europeans and Hispanics, FN3K and/or FN3KRP in Europeans, and G6PD in African-Americans and Hispanics) and we identified an African-ancestry-specific low-frequency variant (rs1039215 in HBG2 and HBE1, minor allele frequency (MAF) = 0.03). The most associated G6PD variant (rs1050828-T, p.Val98Met, MAF = 12% in African-Americans, MAF = 2% in Hispanics) lowered HbA1c (-0.88% in hemizygous males, -0.34% in heterozygous females) and explained 23% of HbA1c variance in African-Americans and 4% in Hispanics. Additionally, we identified a rare distinct G6PD coding variant (rs76723693, p.Leu353Pro, MAF = 0.5%; -0.98% in hemizygous males, -0.46% in heterozygous females) and detected significant association with HbA1c when aggregating rare missense variants in G6PD. We observed similar magnitude and direction of effects for rs1039215 (HBG2) and rs76723693 (G6PD) in the two largest TOPMed African American cohorts, and we replicated the rs76723693 association in the UK Biobank African-ancestry participants. These variants in G6PD and HBG2 were monomorphic in the European and Asian samples. African or Hispanic ancestry individuals carrying G6PD variants may be underdiagnosed for diabetes when screened with HbA1c. Thus, assessment of these variants should be considered for incorporation into precision medicine approaches for diabetes diagnosis.

Efficient Variant Set Mixed Model Association Tests for Continuous and Binary Traits in Large-Scale Whole-Genome Sequencing Studies.

With advances in whole-genome sequencing (WGS) technology, more advanced statistical methods for testing genetic association with rare variants are being developed. Methods in which variants are grouped for analysis are also known as variant-set, gene-based, and aggregate unit tests. The burden test and sequence kernel association test (SKAT) are two widely used variant-set tests, which were originally developed for samples of unrelated individuals and later have been extended to family data with known pedigree structures. However, computationally efficient and powerful variant-set tests are needed to make analyses tractable in large-scale WGS studies with complex study samples. In this paper, we propose the variant-set mixed model association tests (SMMAT) for continuous and binary traits using the generalized linear mixed model framework. These tests can be applied to large-scale WGS studies involving samples with population structure and relatedness, such as in the National Heart, Lung, and Blood Institute's Trans-Omics for Precision Medicine (TOPMed) program. SMMATs share the same null model for different variant sets, and a virtue of this null model, which includes covariates only, is that it needs to be fit only once for all tests in each genome-wide analysis. Simulation studies show that all the proposed SMMATs correctly control type I error rates for both continuous and binary traits in the presence of population structure and relatedness. We also illustrate our tests in a real data example of analysis of plasma fibrinogen levels in the TOPMed program (n = 23,763), using the Analysis Commons, a cloud-based computing platform.

Genetic influence on cognitive development between childhood and adulthood.

Successful cognitive development between childhood and adulthood has important consequences for future mental and physical wellbeing, as well as occupational and financial success. Therefore, delineating the genetic influences underlying changes in cognitive abilities during this developmental period will provide important insights into the biological mechanisms that govern both typical and atypical maturation. Using data from the Philadelphia Neurodevelopmental Cohort (PNC), a large population-based sample of individuals aged 8 to 21 years old (n = 6634), we used an empirical relatedness matrix to establish the heritability of general and specific cognitive functions and determine if genetic factors influence cognitive maturation (i.e., Gene × Age interactions) between childhood and early adulthood. We found that neurocognitive measures across childhood and early adulthood were significantly heritable. Moreover, genetic variance on general cognitive ability, or g, increased significantly between childhood and early adulthood. Finally, we did not find evidence for decay in genetic correlation on neurocognition throughout childhood and adulthood, suggesting that the same genetic factors underlie cognition at different ages throughout this developmental period. Establishing significant Gene × Age interactions in neurocognitive functions across childhood and early adulthood is a necessary first step in identifying genes that influence cognitive development, rather than genes that influence cognition per se. Moreover, since aberrant cognitive development confers risk for several psychiatric disorders, further examination of these Gene × Age interactions may provide important insights into their etiology.

Genetic determinants of metabolic biomarkers and their associations with cardiometabolic traits in Hispanic/Latino adolescents.

BACKGROUND: Metabolic regulation plays a significant role in energy homeostasis, and adolescence is a crucial life stage for the development of cardiometabolic disease (CMD). This study aims to investigate the genetic determinants of metabolic biomarkers-adiponectin, leptin, ghrelin, and orexin-and their associations with CMD risk factors. METHODS: We characterized the genetic determinants of the biomarkers among Hispanic/Latino adolescents of the Santiago Longitudinal Study (SLS) and identified the cumulative effects of genetic variants on adiponectin and leptin using biomarker polygenic risk scores (PRS). We further investigated the direct and indirect effect of the biomarker PRS on downstream body fat percent (BF%) and glycemic traits using structural equation modeling. RESULTS: We identified putatively novel genetic variants associated with the metabolic biomarkers. A substantial amount of biomarker variance was explained by SLS-specific PRS, and the prediction was improved by including the putatively novel loci. Fasting blood insulin and insulin resistance were associated with PRS for adiponectin, leptin, and ghrelin, and BF% was associated with PRS for adiponectin and leptin. We found evidence of substantial mediation of these associations by the biomarker levels. CONCLUSIONS: The genetic underpinnings of metabolic biomarkers can affect the early development of CMD, partly mediated by the biomarkers. IMPACT: This study characterized the genetic underpinnings of four metabolic hormones and investigated their potential influence on adiposity and insulin biology among Hispanic/Latino adolescents. Fasting blood insulin and insulin resistance were associated with polygenic risk score (PRS) for adiponectin, leptin, and ghrelin, with evidence of some degree of mediation by the biomarker levels. Body fat percent (BF%) was also associated with PRS for adiponectin and leptin. This provides important insight on biological mechanisms underlying early metabolic dysfunction and reveals candidates for prevention efforts. Our findings also highlight the importance of ancestrally diverse populations to facilitate valid studies of the genetic architecture of metabolic biomarker levels.

Crossover interference and sex-specific genetic maps shape identical by descent sharing in close relatives.

Simulations of close relatives and identical by descent (IBD) segments are common in genetic studies, yet most past efforts have utilized sex averaged genetic maps and ignored crossover interference, thus omitting features known to affect the breakpoints of IBD segments. We developed Ped-sim, a method for simulating relatives that can utilize either sex-specific or sex averaged genetic maps and also either a model of crossover interference or the traditional Poisson model for inter-crossover distances. To characterize the impact of previously ignored mechanisms, we simulated data for all four combinations of these factors. We found that modeling crossover interference decreases the standard deviation of pairwise IBD proportions by 10.4% on average in full siblings through second cousins. By contrast, sex-specific maps increase this standard deviation by 4.2% on average, and also impact the number of segments relatives share. Most notably, using sex-specific maps, the number of segments half-siblings share is bimodal; and when combined with interference modeling, the probability that sixth cousins have non-zero IBD sharing ranges from 9.0 to 13.1%, depending on the sexes of the individuals through which they are related. We present new analytical results for the distributions of IBD segments under these models and show they match results from simulations. Finally, we compared IBD sharing rates between simulated and real relatives and find that the combination of sex-specific maps and interference modeling most accurately captures IBD rates in real data. Ped-sim is open source and available from https://github.com/williamslab/ped-sim.