Search

Abstract Machine learning can be used to define subtypes of psychiatric conditions based on shared biological foundations of mental disorders. Here we analyzed cross-sectional brain images from 4,222 individuals with schizophrenia and 7038 healthy subjects pooled across 41 international cohorts from the ENIGMA, non-ENIGMA cohorts and public datasets. Using the Subtype and Stage Inference (SuStaIn) algorithm, we identify two distinct neurostructural subgroups by mapping the spatial and temporal ‘trajectory’ of gray matter change in schizophrenia. Subgroup 1 was characterized by an early cortical-predominant loss with enlarged striatum, whereas subgroup 2 displayed an early subcortical-predominant loss in the hippocampus, striatum and other subcortical regions. We confirmed the reproducibility of the two neurostructural subtypes across various sample sites, including Europe, North America and East Asia. This imaging-based taxonomy holds the potential to identify individuals with shared neurobiological attributes, thereby suggesting the viability of redefining existing disorder constructs based on biological factors.

Abstract Genome-wide association studies (GWAS) have become well-powered to detect loci associated with telomere length. However, no prior work has validated genes nominated by GWAS to examine their role in telomere length regulation. We conducted a multi-ancestry meta-analysis of 211,369 individuals and identified five novel association signals. Enrichment analyses of chromatin state and cell-type heritability suggested that blood/immune cells are the most relevant cell type to examine telomere length association signals. We validated specific GWAS associations by overexpressing KBTBD6 or POP5 and demonstrated that both lengthened telomeres. CRISPR/Cas9 deletion of the predicted causal regions in K562 blood cells reduced expression of these genes, demonstrating that these loci are related to transcriptional regulation of KBTBD6 and POP5. Our results demonstrate the utility of telomere length GWAS in the identification of telomere length regulation mechanisms and validate KBTBD6 and POP5 as genes affecting telomere length regulation.

Summary: Background: Metabolic ageing biomarkers may capture the age-related shifts in metabolism, offering a precise representation of an individual’s overall metabolic health. Methods: Utilising comprehensive lipidomic datasets from two large independent population cohorts in Australia (n = 14,833, including 6630 males, 8203 females), we employed different machine learning models, to predict age, and calculated metabolic age scores (mAge). Furthermore, we defined the difference between mAge and age, termed mAgeΔ, which allow us to identify individuals sharing similar age but differing in their metabolic health status. Findings: Upon stratification of the population into quintiles by mAgeΔ, we observed that participants in the top quintile group (Q5) were more likely to have cardiovascular disease (OR = 2.13, 95% CI = 1.62–2.83), had a 2.01-fold increased risk of 12-year incident cardiovascular events (HR = 2.01, 95% CI = 1.45–2.57), and a 1.56-fold increased risk of 17-year all-cause mortality (HR = 1.56, 95% CI = 1.34–1.79), relative to the individuals in the bottom quintile group (Q1). Survival analysis further revealed that men in the Q5 group faced the challenge of reaching a median survival rate due to cardiovascular events more than six years earlier and reaching a median survival rate due to all-cause mortality more than four years earlier than men in the Q1 group. Interpretation: Our findings demonstrate that the mAge score captures age-related metabolic changes, predicts health outcomes, and has the potential to identify individuals at increased risk of metabolic diseases. Funding: The specific funding of this article is provided in the acknowledgements section.

Summary: Background: Decreased levels of circulating ethanolamine plasmalogens [PE(P)], and a concurrent increase in phosphatidylethanolamine (PE) are consistently reported in various cardiometabolic conditions. Here we devised, a plasmalogen score (Pls Score) that mirrors a metabolic signal that encompasses the levels of PE(P) and PE and captures the natural variation in circulating plasmalogens and perturbations in their metabolism associated with disease, diet, and lifestyle. Methods: We utilised, plasma lipidomes from the Australian Obesity, Diabetes and Lifestyle study (AusDiab; n = 10,339, 55% women) a nationwide cohort, to devise the Pls Score and validated this in the Busselton Health Study (BHS; n = 4,492, 56% women, serum lipidome) and in a placebo-controlled crossover trial involving Shark Liver Oil (SLO) supplementation (n = 10, 100% men). We examined the association of the Pls Score with cardiometabolic risk factors, type 2 diabetes mellitus (T2DM), cardiovascular disease and all-cause mortality (over 17 years). Findings: In a model, adjusted for age, sex and BMI, individuals in the top quintile of the Pls Score (Q5) relative to Q1 had an OR of 0.31 (95% CI 0.21–0.43), 0.39 (95% CI 0.25–0.61) and 0.42 (95% CI 0.30–0.57) for prevalent T2DM, incident T2DM and prevalent cardiovascular disease respectively, and a 34% lower mortality risk (HR = 0.66; 95% CI 0.56–0.78). Significant associations between diet and lifestyle habits and Pls Score exist and these were validated through dietary supplementation of SLO that resulted in a marked change in the Pls Score. Interpretation: The Pls Score as a measure that captures the natural variation in circulating plasmalogens, was not only inversely related to cardiometabolic risk and all-cause mortality but also associate with diet and lifestyle. Our results support the potential utility of the Pls Score as a biomarker for metabolic health and its responsiveness to dietary interventions. Further research is warranted to explore the underlying mechanisms and optimise the practical implementation of the Pls Score in clinical and population settings. Funding: National Health and Medical Research Council (NHMRC grant 233200), National Health and Medical Research Council of Australia (Project grant APP1101320), Health Promotion Foundation of Western Australia, and National Health and Medical Research Council of Australia Senior Research Fellowship (#1042095).

Abstract Recent advancements in plasma lipidomic profiling methodology have significantly increased specificity and accuracy of lipid measurements. This evolution, driven by improved chromatographic and mass spectrometric resolution of newer platforms, has made it challenging to align datasets created at different times, or on different platforms. Here we present a framework for harmonising such plasma lipidomic datasets with different levels of granularity in their lipid measurements. Our method utilises elastic-net prediction models, constructed from high-resolution lipidomics reference datasets, to predict unmeasured lipid species in lower-resolution studies. The approach involves (1) constructing composite lipid measures in the reference dataset that map to less resolved lipids in the target dataset, (2) addressing discrepancies between aligned lipid species, (3) generating prediction models, (4) assessing their transferability into the targe dataset, and (5) evaluating their prediction accuracy. To demonstrate our approach, we used the AusDiab population-based cohort (747 lipid species) as the reference to impute unmeasured lipid species into the LIPID study (342 lipid species). Furthermore, we compared measured and imputed lipids in terms of parameter estimation and predictive performance, and validated imputations in an independent study. Our method for harmonising plasma lipidomic datasets will facilitate model validation and data integration efforts.

Background: Socioeconomic Status (SES) is a potent environmental determinant of health. To our knowledge, no assessment of genotype-environment interaction has been conducted to consider the joint effects of socioeconomic status and genetics on risk for metabolic disease. We analyzed data from the Mexican American Family Studies (MAFS) to evaluate the hypothesis that genotype-by-environment interaction (GxE) is an essential determinant of variation in risk factors for metabolic syndrome (MS).Methods: We employed a maximum likelihood estimation of the decomposition of variance components to detect GxE interaction. After excluding individuals with diabetes and individuals on medication for diabetes, hypertension, or dyslipidemia, we analyzed 12 MS risk factors: fasting glucose (FG), fasting insulin (FI), 2-h glucose (2G), 2-h insulin (2I), body mass index (BMI), waist circumference (WC), leptin (LP), high-density lipoprotein-cholesterol (HDL-C), triglycerides (TG), total serum cholesterol (TSC), systolic blood pressure (SBP), and diastolic blood pressure (DBP). Our SES variable used a combined score of Duncan’s socioeconomic index and education years. Heterogeneity in the additive genetic variance across the SES continuum and a departure from unity in the genetic correlation coefficient were taken as evidence of GxE interaction. Hypothesis tests were conducted using standard likelihood ratio tests.Results: We found evidence of GxE for fasting glucose, 2-h glucose, 2-h insulin, BMI, and triglycerides. The genetic effects underlying the insulin/glucose metabolism component of MS are upregulated at the lower end of the SES spectrum. We also determined that the household variance for systolic blood pressure decreased with increasing SES.Conclusion: These results show a significant change in the GxE interaction underlying the major components of MS in response to changes in socioeconomic status. Further mRNA sequencing studies will identify genes and canonical gene pathways to support our molecular-level hypotheses.

Abstract Obesity is a risk factor for type 2 diabetes and cardiovascular disease. However, a substantial proportion of patients with these conditions have a seemingly normal body mass index (BMI). Conversely, not all obese individuals present with metabolic disorders giving rise to the concept of “metabolically healthy obese”. We use lipidomic-based models for BMI to calculate a metabolic BMI score (mBMI) as a measure of metabolic dysregulation associated with obesity. Using the difference between mBMI and BMI (mBMIΔ), we identify individuals with a similar BMI but differing in their metabolic health and disease risk profiles. Exercise and diet associate with mBMIΔ suggesting the ability to modify mBMI with lifestyle intervention. Our findings show that, the mBMI score captures information on metabolic dysregulation that is independent of the measured BMI and so provides an opportunity to assess metabolic health to identify “at risk” individuals for targeted intervention and monitoring.

OBJECTIVES/GOALS: Metabolic dysfunction-associated steatotic liver disease (MASLD) is a major public health concern due to its increasing prevalence and association with type 2 diabetes mellitus. Non-invasive magnetic resonance-based biomarkers can aid in the monitoring of disease progression and identification of patients at risk for chronic liver disease. METHODS/STUDY POPULATION: Over 600 subjects will be recruited from the San Antonio Mexican American Family Study and from a second study, which consists of (i) T2DM patients diagnosed with either MASLD or metabolic dysfunction-associated steatohepatitis (MASH) or (ii) metabolically healthy controls. Hydrogen-1 MRS and diffusion-weighted MRI (DW-MRI) will be used to measure liver fat fraction and liver stiffness biomarkers, respectively. Several potential biomarkers of liver stiffness will be evaluated in vivo using the intravoxel incoherent motion (IVIM) model for DW-MRI. To further improve the diagnostic accuracy of patients with liver fibrosis, we will integrate MRI/MRS data with relevant clinical indicators of hepatic metabolism. Results will be compared to biopsy samples to evaluate the model’s diagnostic accuracy. RESULTS/ANTICIPATED RESULTS: Based on preliminary data, we predict that IVIM will be able to accurately diagnose hepatic fibrosis in patients with MASLD, allowing it to be implemented in clinics with high-field MRI units easily. Previous studies have shown correlations between IVIM estimates and fibrosis stages, however, none included additional clinical indicators of liver disease in their models. We have already found significant differences in metabolic measurements such as fasting plasma glucose and HbA1c levels. Additionally, the use of machine learning in developing these models has shown improvements in the ability to extract features from the data. The aim is to achieve high accuracy and robustness in the staging of liver fibrosis. DISCUSSION/SIGNIFICANCE: Over 100 million people in the US are affected by MASLD. Without treatment, it progresses from hepatic steatosis to MASH, fibrosis (liver stiffening), and ultimately to hepatic cirrhosis and hepatocellular carcinoma (HCC). Continued research efforts and clinical implementation of MRI and MRS are vital in combating the growing burden of MASLD.

Background: Socioeconomic status (SES) is a potent environmental determinant of health. To our knowledge, no assessment of genotype-environment interaction has been conducted to consider the joint effects of socioeconomic status and genetics on risk for cardiovascular disease (CVD). We analyzed Mexican American Family Studies (MAFS) data to evaluate the hypothesis that genotype-by-environment interaction (GxE) is an important determinant of variation in CVD risk factors.Methods: We employed a linear mixed model to investigate GxE in Mexican American extended families. We studied two proxies for CVD [Pooled Cohort Equation Risk Scores/Framingham Risk Scores (FRS/PCRS) and carotid artery intima-media thickness (CA-IMT)] in relation to socioeconomic status as determined by Duncan’s Socioeconomic Index (SEI), years of education, and household income.Results: We calculated heritability for FRS/PCRS and carotid artery intima-media thickness. There was evidence of GxE due to additive genetic variance heterogeneity and genetic correlation for FRS, PCRS, and CA-IMT measures for education (environment) but not for household income or SEI.Conclusion: The genetic effects underlying CVD are dynamically modulated at the lower end of the SES spectrum. There is a significant change in the genetic architecture underlying the major components of CVD in response to changes in education.

Summary: Background: Identifying novel epigenetic signatures associated with serum immunoglobulin E (IgE) may improve our understanding of molecular mechanisms underlying asthma and IgE-mediated diseases. Methods: We performed an epigenome-wide association study using whole blood from Framingham Heart Study (FHS; n = 3,471, 46% females) participants and validated results using the Childhood Asthma Management Program (CAMP; n = 674, 39% females) and the Genetic Epidemiology of Asthma in Costa Rica Study (CRA; n = 787, 41% females). Using the closest gene to each IgE-associated CpG, we highlighted biologically plausible pathways underlying IgE regulation and analyzed the transcription patterns linked to IgE-associated CpGs (expression quantitative trait methylation loci; eQTMs). Using prior UK Biobank summary data from genome-wide association studies of asthma and allergy, we performed Mendelian randomization (MR) for causal inference testing using the IgE-associated CpGs from FHS with methylation quantitative trait loci (mQTLs) as instrumental variables. Findings: We identified 490 statistically significant differentially methylated CpGs associated with IgE in FHS, of which 193 (39.3%) replicated in CAMP and CRA (FDR < 0.05). Gene ontology analysis revealed enrichment in pathways related to transcription factor binding, asthma, and other immunological processes. eQTM analysis identified 124 cis-eQTMs for 106 expressed genes (FDR < 0.05). MR in combination with drug-target analysis revealed CTSB and USP20 as putatively causal regulators of IgE levels (Bonferroni adjusted P

Most genetic association studies have been done on single nucleotide polymorphisms and small indels, while other types of variants have been less studied. Here, the authors use whole genome sequencing in a diverse population to identify and provide experimental evidence for associations between structural variants and blood-cell traits.

Although the common genetic variants contributing to blood lipid levels have been studied, the contribution of rare variants is less understood. Here, the authors perform a rare coding and noncoding variant association study of blood lipid levels using whole genome sequencing data.

A large portion of the heterogeneity in coronavirus disease 2019 (COVID-19) susceptibility and severity of illness (SOI) remains poorly understood. Recent evidence suggests that SARS-CoV-2 infection-associated damage to alveolar epithelial type 2 cells (AT2s) in the distal lung may directly contribute to disease severity and poor prognosis in COVID-19 patients. Our in vitro modeling of SARS-CoV-2 infection in induced pluripotent stem cell (iPSC)-derived AT2s from 10 different individuals showed interindividual variability in infection susceptibility and the postinfection cellular viral load. To understand the underlying mechanism of the AT2′s capacity to regulate SARS-CoV-2 infection and cellular viral load, a genome-wide differential gene expression analysis between the mock and SARS-CoV-2 infection-challenged AT2s was performed. The 1393 genes, which were significantly (one-way ANOVA FDR-corrected p ≤ 0.05; FC abs ≥ 2.0) differentially expressed (DE), suggest significant upregulation of viral infection-related cellular innate immune response pathways (p-value ≤ 0.05; activation z-score ≥ 3.5), and significant downregulation of the cholesterol- and xenobiotic-related metabolic pathways (p-value ≤ 0.05; activation z-score ≤ −3.5). Whilst the effect of post-SARS-CoV-2 infection response on the infection susceptibility and postinfection viral load in AT2s is not clear, interestingly, pre-infection (mock-challenged) expression of 238 DE genes showed a high correlation with the postinfection SARS-CoV-2 viral load (FDR-corrected p-value ≤ 0.05 and r2-absolute ≥ 0.57). The 85 genes whose expression was negatively correlated with the viral load showed significant enrichment in viral recognition and cytokine-mediated innate immune GO biological processes (p-value range: 4.65 × 10−10 to 2.24 × 10−6). The 153 genes whose expression was positively correlated with the viral load showed significant enrichment in cholesterol homeostasis, extracellular matrix, and MAPK/ERK pathway-related GO biological processes (p-value range: 5.06 × 10−5 to 6.53 × 10−4). Overall, our results strongly suggest that AT2s’ pre-infection innate immunity and metabolic state affect their susceptibility to SARS-CoV-2 infection and viral load.

This study of 23,000 non-diabetic individuals highlights loci associated with fasting glucose and fasting insulin in diverse cohorts with whole genome sequence data.

Dysregulation of lipid metabolism is associated with coronary artery disease (CAD). Here, the authors perform GWAS of the serum lipidome to identify variants associated with lipid species that are putatively in the mechanistic pathway to CAD.

IntroductionThe cocktail-party problem refers to the difficulty listeners face when trying to attend to relevant sounds that are mixed with irrelevant ones. Previous studies have shown that solving these problems relies on perceptual as well as cognitive processes. Previously, we showed that speech-reception thresholds (SRTs) on a cocktail-party listening task were influenced by genetic factors. Here, we estimated the degree to which these genetic factors overlapped with those influencing cognitive abilities.MethodsWe measured SRTs and hearing thresholds (HTs) in 493 listeners, who ranged in age from 18 to 91 years old. The same individuals completed a cognitive test battery comprising 18 measures of various cognitive domains. Individuals belonged to large extended pedigrees, which allowed us to use variance component models to estimate the narrow-sense heritability of each trait, followed by phenotypic and genetic correlations between pairs of traits.ResultsAll traits were heritable. The phenotypic and genetic correlations between SRTs and HTs were modest, and only the phenotypic correlation was significant. By contrast, all genetic SRT–cognition correlations were strong and significantly different from 0. For some of these genetic correlations, the hypothesis of complete pleiotropy could not be rejected.DiscussionOverall, the results suggest that there was substantial genetic overlap between SRTs and a wide range of cognitive abilities, including abilities without a major auditory or verbal component. The findings highlight the important, yet sometimes overlooked, contribution of higher-order processes to solving the cocktail-party problem, raising an important caveat for future studies aiming to identify specific genetic factors that influence cocktail-party listening.

Abstract 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.

Abstract Background Hypertriglyceridemia has emerged as a critical coronary artery disease (CAD) risk factor. Rare loss-of-function (LoF) variants in apolipoprotein C-III have been reported to reduce triglycerides (TG) and are cardioprotective in American Indians and Europeans. However, there is a lack of data in other Europeans and non-Europeans. Also, whether genetically increased plasma TG due to ApoC-III is causally associated with increased CAD risk is still unclear and inconsistent. The objectives of this study were to verify the cardioprotective role of earlier reported six LoF variants of APOC3 in South Asians and other multi-ethnic cohorts and to evaluate the causal association of TG raising common variants for increasing CAD risk. Methods We performed gene-centric and Mendelian randomization analyses and evaluated the role of genetic variation encompassing APOC3 for affecting circulating TG and the risk for developing CAD. Results One rare LoF variant (rs138326449) with a 37% reduction in TG was associated with lowered risk for CAD in Europeans (p = 0.007), but we could not confirm this association in Asian Indians (p = 0.641). Our data could not validate the cardioprotective role of other five LoF variants analysed. A common variant rs5128 in the APOC3 was strongly associated with elevated TG levels showing a p-value 2.8 × 10− 424. Measures of plasma ApoC-III in a small subset of Sikhs revealed a 37% increase in ApoC-III concentrations among homozygous mutant carriers than the wild-type carriers of rs5128. A genetically instrumented per 1SD increment of plasma TG level of 15 mg/dL would cause a mild increase (3%) in the risk for CAD (p = 0.042). Conclusions Our results highlight the challenges of inclusion of rare variant information in clinical risk assessment and the generalizability of implementation of ApoC-III inhibition for treating atherosclerotic disease. More studies would be needed to confirm whether genetically raised TG and ApoC-III concentrations would increase CAD risk.

Abstract Chronic kidney disease (CKD) is a persistent impairment of kidney function. Genome-wide association studies (GWAS) have revealed multiple genetic loci associated with CKD susceptibility but the complete genetic basis is not yet clear. Since CKD shares risk factors with cardiovascular diseases and diabetes, there may be pleiotropic loci at play but may go undetected when using single phenotype GWAS. Here, we used multi-phenotype GWAS in the Norfolk Island isolate (n = 380) to identify new loci associated with CKD. We performed a principal components analysis on different combinations of 29 quantitative traits to extract principal components (PCs) representative of multiple correlated phenotypes. GWAS of a PC derived from glomerular filtration rate, serum creatinine, and serum urea identified a suggestive peak (pmin = 1.67 × 10–7) that mapped to KCNIP4. Inclusion of other secondary CKD measurements with these three kidney function traits identified the KCNIP4 locus with GWAS significance (pmin = 1.59 × 10–9). Finally, we identified a group of two SNPs with increased minor allele frequencies as potential functional variants. With the use of genetic isolate and the PCA-based multi-phenotype GWAS approach, we have revealed a potential pleotropic effect locus for CKD. Further studies are required to assess functional relevance of this locus.

IntroductionThis project aimed to investigate the association between biometric components of metabolic syndrome (MetS) with gray matter volume (GMV) obtained with magnetic resonance imaging (MRI) from a large cohort of community-based adults (n = 776) subdivided by age and sex and employing brain regions of interest defined previously as the “Neural Signature of MetS” (NS-MetS).MethodsLipid profiles, biometrics, and regional brain GMV were obtained from the Genetics of Brain Structure (GOBS) image archive. Participants underwent T1-weighted MR imaging. MetS components (waist circumference, fasting plasma glucose, triglycerides, HDL cholesterol, and blood pressure) were defined using the National Cholesterol Education Program Adult Treatment Panel III. Subjects were grouped by age: early adult (18–25 years), young adult (26–45 years), and middle-aged adult (46–65 years). Linear regression modeling was used to investigate associations between MetS components and GMV in five brain regions comprising the NS-MetS: cerebellum, brainstem, orbitofrontal cortex, right insular/limbic cluster and caudate.ResultsIn both men and women of each age group, waist circumference was the single component most strongly correlated with decreased GMV across all NS-MetS regions. The brain region most strongly correlated to all MetS components was the posterior cerebellum.ConclusionThe posterior cerebellum emerged as the region most significantly associated with MetS individual components, as the only region to show decreased GMV in young adults, and the region with the greatest variance between men and women. We propose that future studies investigating neurological effects of MetS and its comorbidities—namely diabetes and obesity—should consider the NS-MetS and the differential effects of age and sex.

Summary: Communicating in everyday situations requires solving the cocktail-party problem, or segregating the acoustic mixture into its constituent sounds and attending to those of most interest. Humans show dramatic variation in this ability, leading some to experience real-world problems irrespective of whether they meet criteria for clinical hearing loss. Here, we estimated the genetic contribution to cocktail-party listening by measuring speech-reception thresholds (SRTs) in 425 people from large families and ranging in age from 18 to 91 years. Roughly half the variance of SRTs was explained by genes (h2 = 0.567). The genetic correlation between SRTs and hearing thresholds (HTs) was medium (ρG = 0.392), suggesting that the genetic factors influencing cocktail-party listening were partially distinct from those influencing sound sensitivity. Aging and socioeconomic status also strongly influenced SRTs. These findings may represent a first step toward identifying genes for “hidden hearing loss,” or hearing problems in people with normal HTs.

Penetrance of variants in monogenic disease and clinical utility of common polygenic variation has not been well explored on a large-scale. Here, the authors use exome sequencing data from 77,184 individuals to generate penetrance estimates and assess the utility of polygenic variation in risk prediction of monogenic variants.

This study examines the impact of G × E interaction effects on non-alcoholic fatty liver disease (NAFLD) among Mexican Americans in the Rio Grande Valley (RGV) of South Texas. We examined potential G × E interaction using variance components models and likelihood-based statistical inference in the phenotypic expression of NAFLD, including hepatic steatosis and hepatic fibrosis (identified using vibration controlled transient elastography and controlled attenuation parameter measured by the FibroScan Device). We screened for depression using the Beck Depression Inventory-II (BDI-II). We identified significant G × E interactions for hepatic fibrosis × BDI-II. These findings provide evidence that genetic factors interact with depression to influence the expression of hepatic fibrosis.

IntroductionCardiovascular disease (CVD) is the leading cause of mortality worldwide and is the leading cause of death in the US. Lipid dysregulation is a well-known precursor to metabolic diseases, including CVD. There is a growing body of literature that suggests MRI-derived epicardial fat volume, or epicardial adipose tissue (EAT) volume, is linked to the development of coronary artery disease. Interestingly, epicardial fat is also actively involved in lipid and energy homeostasis, with epicardial adipose tissue having a greater capacity for release and uptake of free fatty acids. However, there is a scarcity of knowledge on the influence of plasma lipids on EAT volume.AimThe focus of this study is on the identification of novel lipidomic species associated with CMRI-derived measures of epicardial fat in Mexican American individuals.MethodsWe performed lipidomic profiling on 200 Mexican American individuals. High-throughput mass spectrometry enabled rapid capture of precise lipidomic profiles, providing measures of 799 unique species from circulating plasma samples. Because of our extended pedigree design, we utilized a standard quantitative genetic linear mixed model analysis to determine whether lipids were correlated with EAT by formally testing for association between each lipid species and the CMRI epicardial fat phenotype.ResultsAfter correction for multiple testing using the FDR approach, we identified 135 lipid species showing significant association with epicardial fat. Of those, 131 lipid species were positively correlated with EAT, where increased circulating lipid levels were correlated with increased epicardial fat. Interestingly, the top 10 lipid species associated with an increased epicardial fat volume were from the deoxyceramide (Cer(m)) and triacylglycerol (TG) families. Deoxyceramides are atypical and neurotoxic sphingolipids. Triacylglycerols are an abundant lipid class and comprise the bulk of storage fat in tissues. Pathologically elevated TG and Cer(m) levels are related to CVD risk and, in our study, to EAT volume.ConclusionOur results indicate that specific lipid abnormalities such as enriched saturated triacylglycerols and the presence of toxic ceramides Cer(m) in plasma of our individuals could precede CVD with increased EAT volume.

The influence of X chromosome genetic variation on blood lipids and coronary heart disease (CHD) is not well understood. Here, the authors analyse X chromosome sequencing data across 65,322 multi-ancestry individuals, identifying associations of the Xq23 locus with lipid changes and reduced risk of CHD and diabetes mellitus.

Abstract Low-frequency 1q21.1 distal deletion and duplication copy number variant (CNV) carriers are predisposed to multiple neurodevelopmental disorders, including schizophrenia, autism and intellectual disability. Human carriers display a high prevalence of micro- and macrocephaly in deletion and duplication carriers, respectively. The underlying brain structural diversity remains largely unknown. We systematically called CNVs in 38 cohorts from the large-scale ENIGMA-CNV collaboration and the UK Biobank and identified 28 1q21.1 distal deletion and 22 duplication carriers and 37,088 non-carriers (48% male) derived from 15 distinct magnetic resonance imaging scanner sites. With standardized methods, we compared subcortical and cortical brain measures (all) and cognitive performance (UK Biobank only) between carrier groups also testing for mediation of brain structure on cognition. We identified positive dosage effects of copy number on intracranial volume (ICV) and total cortical surface area, with the largest effects in frontal and cingulate cortices, and negative dosage effects on caudate and hippocampal volumes. The carriers displayed distinct cognitive deficit profiles in cognitive tasks from the UK Biobank with intermediate decreases in duplication carriers and somewhat larger in deletion carriers—the latter potentially mediated by ICV or cortical surface area. These results shed light on pathobiological mechanisms of neurodevelopmental disorders, by demonstrating gene dose effect on specific brain structures and effect on cognitive function.

Abstract High concentrations of carotenoids are protective against cardiometabolic risk traits (CMTs) in adults and children. We recently showed in non-diabetic Mexican American (MA) children that serum α-carotene and β-carotene are inversely correlated with obesity measures and triglycerides and positively with HDL cholesterol and that they were under strong genetic influences. Additionally, we previously described a Pediatric Metabolic Index (PMI) that helps in the identification of children who are at risk for cardiometabolic diseases. Here, we quantified serum lycopene and β-cryptoxanthin concentrations in approximately 580 children from MA families using an ultraperformance liquid chromatography-photodiode array and determined their heritabilities and correlations with CMTs. Using response surface methodology (RSM), we determined two-way interactions of carotenoids and PMI on Matsuda insulin sensitivity index (ISI). The concentrations of lycopene and β-cryptoxanthin were highly heritable [h2 = 0.98, P = 7 × 10–18 and h2 = 0.58, P = 1 × 10–7]. We found significant (P ≤ 0.05) negative phenotypic correlations between β-cryptoxanthin and five CMTs: body mass index (− 0.22), waist circumference (− 0.25), triglycerides (− 0.18), fat mass (− 0.23), fasting glucose (− 0.09), and positive correlations with HDL cholesterol (0.29). In contrast, lycopene only showed a significant negative correlation with fasting glucose (− 0.08) and a positive correlation with HDL cholesterol (0.18). Importantly, we found that common genetic influences significantly contributed to the observed phenotypic correlations. RSM showed that increased serum concentrations of α- and β-carotenoids rather than that of β-cryptoxanthin or lycopene had maximal effects on ISI. In summary, our findings suggest that the serum carotenoids are under strong additive genetic influences and may have differential effects on susceptibility to CMTs in children.

Abstract Insulin is an essential hormone that regulates glucose homeostasis and metabolism. Insulin resistance (IR) arises when tissues fail to respond to insulin, and it leads to serious health problems including Type 2 Diabetes (T2D). Obesity is a major contributor to the development of IR and T2D. We previously showed that gene expression of alcohol dehydrogenase 1B (ADH1B) was inversely correlated with obesity and IR in subcutaneous adipose tissue of Mexican Americans. In the current study, a meta-analysis of the relationship between ADH1B expression and BMI in Mexican Americans, African Americans, Europeans, and Pima Indians verified that BMI was increased with decreased ADH1B expression. Using established human subcutaneous pre-adipocyte cell lines derived from lean (BMI

Drugs targeting cardiovascular disease (CVD) can have negative consequences for liver function. Here, the authors combine genome wide analyses on 69,479 individuals to identify loss-of-function variants with beneficial effects on CVD-related traits without negative impacts on liver function.

Polygenic risk scores (PRSs) aggregate the effects of genetic variants across the genome and are used to predict risk of complex diseases, such as obesity. Current PRSs only include common variants (minor allele frequency (MAF) ≥1%), whereas the contribution of rare variants in PRSs to predict disease remains unknown. Here, we examine whether augmenting the standard common variant PRS (PRScommon) with a rare variant PRS (PRSrare) improves prediction of obesity. We used genome-wide genotyped and imputed data on 451,145 European-ancestry participants of the UK Biobank, as well as whole exome sequencing (WES) data on 184,385 participants. We performed single variant analyses (for both common and rare variants) and gene-based analyses (for rare variants) for association with BMI (kg/m2), obesity (BMI ≥ 30 kg/m2), and extreme obesity (BMI ≥ 40 kg/m2). We built PRSscommon and PRSsrare using a range of methods (Clumping+Thresholding [C+T], PRS-CS, lassosum, gene-burden test). We selected the best-performing PRSs and assessed their performance in 36,757 European-ancestry unrelated participants with whole genome sequencing (WGS) data from the Trans-Omics for Precision Medicine (TOPMed) program. The best-performing PRScommon explained 10.1% of variation in BMI, and 18.3% and 22.5% of the susceptibility to obesity and extreme obesity, respectively, whereas the best-performing PRSrare explained 1.49%, and 2.97% and 3.68%, respectively. The PRSrare was associated with an increased risk of obesity and extreme obesity (ORobesity = 1.37 per SDPRS, Pobesity = 1.7x10-85; ORextremeobesity = 1.55 per SDPRS, Pextremeobesity = 3.8x10-40), which was attenuated, after adjusting for PRScommon (ORobesity = 1.08 per SDPRS, Pobesity = 9.8x10-6; ORextremeobesity= 1.09 per SDPRS, Pextremeobesity = 0.02). When PRSrare and PRScommon are combined, the increase in explained variance attributed to PRSrare was small (incremental Nagelkerke R2 = 0.24% for obesity and 0.51% for extreme obesity). Consistently, combining PRSrare to PRScommon provided little improvement to the prediction of obesity (PRSrare AUC = 0.591; PRScommon AUC = 0.708; PRScombined AUC = 0.710). In summary, while rare variants show convincing association with BMI, obesity and extreme obesity, the PRSrare provides limited improvement over PRScommon in the prediction of obesity risk, based on these large populations.

To assess the burden of type 2 diabetes (T2D) and its genetic profile in endogamous populations of India given the paucity of data, we aimed to determine the prevalence of T2D and estimate its heritability using family-based cohorts from three distinct Endogamous Ethnic Groups (EEGs) representing Northern (Rajasthan [Agarwals: AG]) and Southern (Tamil Nadu [Chettiars: CH] and Andhra Pradesh [Reddys: RE]) states of India. For comparison, family-based data collected previously from another North Indian Punjabi Sikh (SI) EEG was used. In addition, we examined various T2D-related cardiometabolic traits and determined their heritabilities. These studies were conducted as part of the Indian Diabetes Genetic Studies in collaboration with US (INDIGENIUS) Consortium. The pedigree, demographic, phenotypic, covariate data and samples were collected from the CH, AG, and RE EEGs. The status of T2D was defined by ADA guidelines (fasting glucose ≥ 126 mg/dl or HbA1c ≥ 6.5% and/or use of diabetes medication/history). The prevalence of T2D in CH (N = 517, families = 21, mean age = 47y, mean BMI = 27), AG (N = 530, Families = 25, mean age = 43y, mean BMI = 27), and RE (N = 500, Families = 22, mean age = 46y, mean BMI = 27) was found to be 33%, 37%, and 36%, respectively, Also, the study participants from these EEGs were found to be at increased cardiometabolic risk (e.g., obesity and prediabetes). Similar characteristics for the SI EEG (N = 1,260, Families = 324, Age = 51y, BMI = 27, T2D = 75%) were obtained previously. We used the variance components approach to carry out genetic analyses after adjusting for covariate effects. The heritability (h2) estimates of T2D in the CH, RE, SI, and AG were found to be 30%, 46%, 54%, and 82% respectively, and statistically significant (P ≤ 0.05). Other T2D related traits (e.g., BMI, lipids, blood pressure) in AG, CH, and RE EEGs exhibited strong additive genetic influences (h2 range: 17% [triglycerides/AG and hs-CRP/RE] - 86% [glucose/non-T2D/AG]). Our findings highlight the high burden of T2D in Indian EEGs with significant and differential additive genetic influences on T2D and related traits.

Hispanic/Latinos have been underrepresented in genome-wide association studies (GWAS) for anthropometric traits despite their notable anthropometric variability, ancestry proportions, and high burden of growth stunting and overweight/obesity. To address this knowledge gap, we analyzed densely imputed genetic data in a sample of Hispanic/Latino adults to identify and fine-map genetic variants associated with body mass index (BMI), height, and BMI-adjusted waist-to-hip ratio (WHRadjBMI). We conducted a GWAS of 18 studies/consortia as part of the Hispanic/Latino Anthropometry (HISLA) Consortium (stage 1, n = 59,771) and generalized our findings in 9 additional studies (stage 2, n = 10,538). We conducted a trans-ancestral GWAS with summary statistics from HISLA stage 1 and existing consortia of European and African ancestries. In our HISLA stage 1 + 2 analyses, we discovered one BMI locus, as well as two BMI signals and another height signal each within established anthropometric loci. In our trans-ancestral meta-analysis, we discovered three BMI loci, one height locus, and one WHRadjBMI locus. We also identified 3 secondary signals for BMI, 28 for height, and 2 for WHRadjBMI in established loci. We show that 336 known BMI, 1,177 known height, and 143 known WHRadjBMI (combined) SNPs demonstrated suggestive transferability (nominal significance and effect estimate directional consistency) in Hispanic/Latino adults. Of these, 36 BMI, 124 height, and 11 WHRadjBMI SNPs were significant after trait-specific Bonferroni correction. Trans-ancestral meta-analysis of the three ancestries showed a small-to-moderate impact of uncorrected population stratification on the resulting effect size estimates. Our findings demonstrate that future studies may also benefit from leveraging diverse ancestries and differences in linkage disequilibrium patterns to discover novel loci and additional signals with less residual population stratification.

CVD is the leading cause of death worldwide, and genetic investigations into the human lipidome may provide insight into CVD risk. The aim of this study was to estimate the heritability of circulating lipid species and their genetic correlation with CVD traits. Targeted lipidomic profiling was performed on 4,492 participants from the Busselton Family Heart Study to quantify the major fatty acids of 596 lipid species from 33 classes. We estimated narrow-sense heritabilities of lipid species/classes and their genetic correlations with eight CVD traits: BMI, HDL-C, LDL-C, triglycerides, total cholesterol, waist-hip ratio, systolic blood pressure, and diastolic blood pressure. We report heritabilities and genetic correlations of new lipid species/subclasses, including acylcarnitine (AC), ubiquinone, sulfatide, and oxidized cholesteryl esters. Over 99% of lipid species were significantly heritable (h2: 0.06–0.50) and all lipid classes were significantly heritable (h2: 0.14–0.50). The monohexosylceramide and AC classes had the highest median heritabilities (h2 = 0.43). The largest genetic correlation was between clinical triglycerides and total diacylglycerol (rg = 0.88). We observed novel positive genetic correlations between clinical triglycerides and phosphatidylglycerol species (rg: 0.64–0.82), and HDL-C and alkenylphosphatidylcholine species (rg: 0.45–0.74). Overall, 51% of the 4,768 lipid species-CVD trait genetic correlations were statistically significant after correction for multiple comparisons. This is the largest lipidomic study to address the heritability of lipids and their genetic correlation with CVD traits. Future work includes identifying putative causal genetic variants for lipid species and CVD using genome-wide SNP and whole-genome sequencing data.

Abstract Introduction Coronavirus disease 2019 (COVID‐19) has caused >3.5 million deaths worldwide and affected >160 million people. At least twice as many have been infected but remained asymptomatic or minimally symptomatic. COVID‐19 includes central nervous system manifestations mediated by inflammation and cerebrovascular, anoxic, and/or viral neurotoxicity mechanisms. More than one third of patients with COVID‐19 develop neurologic problems during the acute phase of the illness, including loss of sense of smell or taste, seizures, and stroke. Damage or functional changes to the brain may result in chronic sequelae. The risk of incident cognitive and neuropsychiatric complications appears independent from the severity of the original pulmonary illness. It behooves the scientific and medical community to attempt to understand the molecular and/or systemic factors linking COVID‐19 to neurologic illness, both short and long term. Methods This article describes what is known so far in terms of links among COVID‐19, the brain, neurological symptoms, and Alzheimer's disease (AD) and related dementias. We focus on risk factors and possible molecular, inflammatory, and viral mechanisms underlying neurological injury. We also provide a comprehensive description of the Alzheimer's Association Consortium on Chronic Neuropsychiatric Sequelae of SARS‐CoV‐2 infection (CNS SC2) harmonized methodology to address these questions using a worldwide network of researchers and institutions. Results Successful harmonization of designs and methods was achieved through a consensus process initially fragmented by specific interest groups (epidemiology, clinical assessments, cognitive evaluation, biomarkers, and neuroimaging). Conclusions from subcommittees were presented to the whole group and discussed extensively. Presently data collection is ongoing at 19 sites in 12 countries representing Asia, Africa, the Americas, and Europe. Discussion The Alzheimer's Association Global Consortium harmonized methodology is proposed as a model to study long‐term neurocognitive sequelae of SARS‐CoV‐2 infection. Key Points The following review describes what is known so far in terms of molecular and epidemiological links among COVID‐19, the brain, neurological symptoms, and AD and related dementias (ADRD) The primary objective of this large‐scale collaboration is to clarify the pathogenesis of ADRD and to advance our understanding of the impact of a neurotropic virus on the long‐term risk of cognitive decline and other CNS sequelae. No available evidence supports the notion that cognitive impairment after SARS‐CoV‐2 infection is a form of dementia (ADRD or otherwise). The longitudinal methodologies espoused by the consortium are intended to provide data to answer this question as clearly as possible controlling for possible confounders. Our specific hypothesis is that SARS‐CoV‐2 triggers ADRD‐like pathology following the extended olfactory cortical network (EOCN) in older individuals with specific genetic susceptibility. The proposed harmonization strategies and flexible study designs offer the possibility to include large samples of under‐represented racial and ethnic groups, creating a rich set of harmonized cohorts for future studies of the pathophysiology, determinants, long‐term consequences, and trends in cognitive aging, ADRD, and vascular disease. We provide a framework for current and future studies to be carried out within the Consortium. and offers a “green paper” to the research community with a very broad, global base of support, on tools suitable for low‐ and middle‐income countries aimed to compare and combine future longitudinal data on the topic. The Consortium proposes a combination of design and statistical methods as a means of approaching causal inference of the COVID‐19 neuropsychiatric sequelae. We expect that deep phenotyping of neuropsychiatric sequelae may provide a series of candidate syndromes with phenomenological and biological characterization that can be further explored. By generating high‐quality harmonized data across sites we aim to capture both descriptive and, where possible, causal associations.

Summary: Genetic studies on telomere length are important for understanding age-related diseases. Prior GWASs for leukocyte TL have been limited to European and Asian populations. Here, we report the first sequencing-based association study for TL across ancestrally diverse individuals (European, African, Asian, and Hispanic/Latino) from the NHLBI Trans-Omics for Precision Medicine (TOPMed) program. We used whole-genome sequencing (WGS) of whole blood for variant genotype calling and the bioinformatic estimation of telomere length in n = 109,122 individuals. We identified 59 sentinel variants (p < 5 × 10−9) in 36 loci associated with telomere length, including 20 newly associated loci (13 were replicated in external datasets). There was little evidence of effect size heterogeneity across populations. Fine-mapping at OBFC1 indicated that the independent signals colocalized with cell-type-specific eQTLs for OBFC1 (STN1). Using a multi-variant gene-based approach, we identified two genes newly implicated in telomere length, DCLRE1B (SNM1B) and PARN. In PheWAS, we demonstrated that our TL polygenic trait scores (PTSs) were associated with an increased risk of cancer-related phenotypes.

Imaging genetics analyses use neuroimaging traits as intermediate phenotypes to infer the degree of genetic contribution to brain structure and function in health and/or illness. Coefficients of relatedness (CR) summarize the degree of genetic similarity among subjects and are used to estimate the heritability – the proportion of phenotypic variance explained by genetic factors. The CR can be inferred directly from genome-wide genotype data to explain the degree of shared variation in common genetic polymorphisms (SNP-heritability) among related or unrelated subjects. We developed a central processing and graphics processing unit (CPU and GPU) accelerated Fast and Powerful Heritability Inference (FPHI) approach that linearizes likelihood calculations to overcome the ∼N2–3 computational effort dependency on sample size of classical likelihood approaches. We calculated for 60 regional and 1.3 × 105 voxel-wise traits in N = 1,206 twin and sibling participants from the Human Connectome Project (HCP) (550 M/656 F, age = 28.8 ± 3.7 years) and N = 37,432 (17,531 M/19,901 F; age = 63.7 ± 7.5 years) participants from the UK Biobank (UKBB). The FPHI estimates were in excellent agreement with heritability values calculated using Genome-wide Complex Trait Analysis software (r = 0.96 and 0.98 in HCP and UKBB sample) while significantly reducing computational (102–4 times). The regional and voxel-wise traits heritability estimates for the HCP and UKBB were likewise in excellent agreement (r = 0.63–0.76, p < 10−10). In summary, the hardware-accelerated FPHI made it practical to calculate heritability values for voxel-wise neuroimaging traits, even in very large samples such as the UKBB. The patterns of additive genetic variance in neuroimaging traits measured in a large sample of related and unrelated individuals showed excellent agreement regardless of the estimation method. The code and instruction to execute these analyses are available at www.solar-eclipse-genetics.org.

The Rio Grande Valley (RGV) in South Texas has one of the highest prevalence of obesity and type 2 diabetes (T2D) in the United States (US). We report for the first time the T2D prevalence in persons with HIV (PWH) in the RGV and the interrelationship between T2D, cardiometabolic risk factors, HIV-related indices, and antiretroviral therapies (ART). The PWH in this study received medical care at Valley AIDS Council (VAC) clinic sites located in Harlingen and McAllen, Texas. Henceforth, this cohort will be referred to as Valley AIDS Council Cohort (VACC). Cross-sectional analyses were conducted using retrospective data obtained from 1,827 registries. It included demographic and anthropometric variables, cardiometabolic traits, and HIV-related virological and immunological indices. For descriptive statistics, we used mean values of the quantitative variables from unbalanced visits across 20 months. Robust regression methods were used to determine the associations. For comparisons, we used cardiometabolic trait data obtained from HIV-uninfected San Antonio Mexican American Family Studies (SAMAFS; N = 2,498), and the Mexican American population in the National Health and Nutrition Examination Survey (HHANES; N = 5,989). The prevalence of T2D in VACC was 51% compared to 27% in SAMAFS and 19% in HHANES, respectively. The PWH with T2D in VACC were younger (4.7 years) and had lower BMI (BMI 2.43 units less) when compared to SAMAFS individuals. In contrast, VACC individuals had increased blood pressure and dyslipidemia. The increased T2D prevalence in VACC was independent of BMI. Within the VACC, ART was associated with viral load and CD4+ T cell counts but not with metabolic dysfunction. Notably, we found that individuals with any INSTI combination had higher T2D risk: OR 2.08 (95%CI 1.67, 2.6; p < 0.001). In summary, our results suggest that VACC individuals may develop T2D at younger ages independent of obesity. The high burden of T2D in these individuals necessitates rigorously designed longitudinal studies to draw potential causal inferences and develop better treatment regimens.

The de novo ceramide synthesis pathway is essential to human biology and health, but genetic influences remain unexplored. The core function of this pathway is the generation of biologically active ceramide from its precursor, dihydroceramide. Dihydroceramides have diverse, often protective, biological roles; conversely, increased ceramide levels are biomarkers of complex disease. To explore the genetics of the ceramide synthesis pathway, we searched for deleterious nonsynonymous variants in the genomes of 1,020 Mexican Americans from extended pedigrees. We identified a Hispanic ancestry-specific rare functional variant, L175Q, in delta 4-desaturase, sphingolipid 1 (DEGS1), a key enzyme in the pathway that converts dihydroceramide to ceramide. This amino acid change was significantly associated with large increases in plasma dihydroceramides. Indexes of DEGS1 enzymatic activity were dramatically reduced in heterozygotes. CRISPR/Cas9 genome editing of HepG2 cells confirmed that the L175Q variant results in a partial loss of function for the DEGS1 enzyme. Understanding the biological role of DEGS1 variants, such as L175Q, in ceramide synthesis may improve the understanding of metabolic-related disorders and spur ongoing research of drug targets along this pathway.

Sea turtle fibropapillomatosis (FP) is a tumor promoting disease that is one of several threats globally to endangered sea turtle populations. The prevalence of FP is highest in green sea turtle (Chelonia mydas) populations, and historically has shown considerable temporal growth. FP tumors can significantly affect the ability of turtles to forage for food and avoid predation and can grow to debilitating sizes. In the current study, based in South Texas, we have applied transcriptome sequencing to FP tumors and healthy control tissue to study the gene expression profiles of FP. By identifying differentially expressed turtle genes in FP, and matching these genes to their closest human ortholog we draw on the wealth of human based knowledge, specifically human cancer, to identify new insights into the biology of sea turtle FP. We show that several genes aberrantly expressed in FP tumors have known tumor promoting biology in humans, including CTHRC1 and NLRC5, and provide support that disruption of the Wnt signaling pathway is a feature of FP. Further, we profiled the expression of current targets of immune checkpoint inhibitors from human oncology in FP tumors and identified potential candidates for future studies.

AbstractBackgroundQuestions remain regarding whether genetic influences on early life psychopathology overlap with cognition and show developmental variation.MethodsUsing data from 9,421 individuals aged 8–21 from the Philadelphia Neurodevelopmental Cohort, factors of psychopathology were generated using a bifactor model of item-level data from a psychiatric interview. Five orthogonal factors were generated: anxious-misery (mood and anxiety), externalizing (attention deficit hyperactivity and conduct disorder), fear (phobias), psychosis-spectrum, and a general factor. Genetic analyses were conducted on a subsample of 4,662 individuals of European American ancestry. A genetic relatedness matrix was used to estimate heritability of these factors, and genetic correlations with executive function, episodic memory, complex reasoning, social cognition, motor speed, and general cognitive ability. Gene × Age analyses determined whether genetic influences on these factors show developmental variation.ResultsExternalizing was heritable (h2 = 0.46, p = 1 × 10−6), but not anxious-misery (h2 = 0.09, p = 0.183), fear (h2 = 0.04, p = 0.337), psychosis-spectrum (h2 = 0.00, p = 0.494), or general psychopathology (h2 = 0.21, p = 0.040). Externalizing showed genetic overlap with face memory (ρg = −0.412, p = 0.004), verbal reasoning (ρg = −0.485, p = 0.001), spatial reasoning (ρg = −0.426, p = 0.010), motor speed (ρg = 0.659, p = 1x10−4), verbal knowledge (ρg = −0.314, p = 0.002), and general cognitive ability (g)(ρg = −0.394, p = 0.002). Gene × Age analyses revealed decreasing genetic variance (γg = −0.146, p = 0.004) and increasing environmental variance (γe = 0.059, p = 0.009) on externalizing.ConclusionsCognitive impairment may be a useful endophenotype of externalizing psychopathology and, therefore, help elucidate its pathophysiological underpinnings. Decreasing genetic variance suggests that gene discovery efforts may be more fruitful in children than adolescents or young adults.

AimAcanthosis nigricans (AN) is a strong correlate of obesity and is considered a marker of insulin resistance (IR). AN is associated with various other cardiometabolic risk factors (CMRFs). However, the direct causal relationship of IR with AN in obesity has been debated. Therefore, we aimed to examine the complex causal relationships among the troika of AN, obesity, and IR in Mexican Americans (MAs).MethodsWe used data from 670 non-diabetic MA children, aged 6-17 years (49% girls). AN (prevalence 33%) severity scores (range 0-5) were used as a quasi-quantitative trait (AN-q) for analysis. We used the program SOLAR for determining phenotypic, genetic, and environmental correlations between AN-q and CMRFs (e.g., BMI, HOMA-IR, lipids, blood pressure, hs-C-reactive protein (CRP), and Harvard physical fitness score (PFS)). The genetic and environmental correlations were subsequently used in mediation analysis (AMOS program). Model comparisons were made using goodness-of-fit indexes.ResultsHeritability of AN-q was 0.75 (pConclusionOur study suggests that obesity explains the association of IR with AN, but no causal relationship between IR and AN in Mexican American children.

Abstract High-throughput platforms allow the characterization of thousands of previously known methylation sites. These platforms have great potential for investigating the epigenetic effects that are partially responsible for gene expression control. Methylation sites provide a bridge for the investigation of real-time environmental contributions on genomic events by the alteration of methylation status of those sites. Using the data provided by GAW20’s organization committee, we calculated the heritability estimates of each cytosine-phosphate-guanine (CpG) island before and after the use of fenofibrate, a lipid-control drug. Surprisingly, we detected substantially high heritability estimates before drug usage. This somewhat unexpected high sample correlation was corrected by the use of principal components and the distributions of heritability estimates before and after fenofibrate treatment, which made the distributions comparable. The methylation sites located near a gene were collected and a genetic relationship matrix estimated to represent the overall correlation between samples. We implemented a random-effect association test to screen genes whose methylation patterns partially explain the observable high-density lipoprotein (HDL) heritability. Our leading association was observed for the TMEM52 gene that encodes a transmembrane protein, and is largely expressed in the liver, had not been previously associated with HDL until this manuscript. Using a variance component decomposition framework with the linear mixed model allows the integration of data from different sources, such as methylation, gene expression, metabolomics, and proteomics. The decomposition of the genetic variance component decomposition provides a flexible analytical approach for the challenges of this new omics era.

Abstract The heritability of a phenotype is an estimation of the percent of variance in that phenotype that is attributable to additive genetic factors. Heritability is optimally estimated in family-based sample populations. Traditionally, this involves use of a pedigree-based kinship coefficient generated from the collected genealogical relationships between family members. An alternative, when dense genotype data are available, is to directly measure the empirical kinship between samples. This study compares the use of pedigree and empirical kinships in the GAW20 data set. Two phenotypes were assessed: triglyceride levels and high-density lipoprotein cholesterol (HDL-C) levels pre- and postintervention with the cholesterol-reducing drug fenofibrate. Using SOLAR (Sequential Oligogenic Linkage Analysis Routines), pedigree-based kinships and empirically calculated kinships (using IBDLD and LDAK) were used to calculate phenotype heritability. In addition, a genome-wide association study was conducted using each kinship model for each phenotype to identify genetic variants significantly associated with phenotypic variation. The variant rs247617 was significantly associated with HDL-C levels both pre- and post-fenofibrate intervention. Overall, the phenotype heritabilities calculated using pedigree based kinships or either of the empirical kinships generated using IBDLD or LDAK were comparable. Phenotype heritabilities estimated from empirical kinships generated using IBDLD were closest to the pedigree-based estimations. Given that there was not an appreciable amount of unknown relatedness between the pedigrees in this data set, a large increase in heritability in using empirical kinship was not expected, and our calculations support this. Importantly, these results demonstrate that when sufficient genotypic data are available, empirical kinship estimation is a practical alternative to using pedigree-based kinships.

Abstract Genome-wide association studies have helped us identify a wealth of genetic variants associated with complex human phenotypes. Because most variants explain a small portion of the total phenotypic variation, however, marker-based studies remain limited in their ability to predict such phenotypes. Here, we show how modern statistical genetic techniques borrowed from animal breeding can be employed to increase the accuracy of genomic prediction of complex phenotypes and the power of genetic mapping studies. Specifically, using the triglyceride data of the GAW20 data set, we apply genomic-best linear unbiased prediction (G-BLUP) methods to obtain empirical genetic values (EGVs) for each triglyceride phenotype and each individual. We then study 2 different factors that influence the prediction accuracy of G-BLUP for the analysis of human data: (a) the choice of kinship matrix, and (b) the overall level of relatedness. The resulting genetic values represent the total genetic component for the phenotype of interest and can be used to represent a trait without its environmental component. Finally, using empirical data, we demonstrate how this method can be used to increase the power of genetic mapping studies. In sum, our results show that dense genome-wide data can be used in a wider scope than previously anticipated.

Abstract We conducted a genome-wide linkage scan to detect loci that influence the levels of fasting triglycerides in plasma. Fasting triglyceride levels were available at 4 time points (visits), 2 pre- and 2 post-fenofibrate intervention. Multipoint identity-by-descent (MIBD) matrices were derived from genotypes using IBDLD. Variance-component linkage analyses were then conducted using SOLAR (Sequential Oligogenic Linkage Analysis Routines). We found evidence of linkage (logarithm of odds [LOD] ≥3) at 5 chromosomal regions with triglyceride levels in plasma. The highest LOD scores were observed for linkage to the estimated genetic value (additive genetic component) of the log-normalized triglyceride levels in plasma. Our results suggest that a chromosome 10 locus at 37 cM (LODpre = 3.01, LODpost = 3.72) influences fasting triglyceride levels in plasma regardless of the fenofibrate intervention, and that loci in chromosomes 1 at 170 cM and 4 at 24 cM ceases to affect the triglyceride levels when fenofibrate is present, while the regions in chromosomes 6 at 136 to 162 cM and 11 at 39 to 40 cM appear to influence triglyceride levels in response to fenofibrate.

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.

Most genome-wide association and fine-mapping studies to date have been conducted in individuals of European descent, and genetic studies of populations of Hispanic/Latino and African ancestry are limited. In addition, these populations have more complex linkage disequilibrium structure. In order to better define the genetic architecture of these understudied populations, we leveraged >100,000 phased sequences available from deep-coverage whole genome sequencing through the multi-ethnic NHLBI Trans-Omics for Precision Medicine (TOPMed) program to impute genotypes into admixed African and Hispanic/Latino samples with genome-wide genotyping array data. We demonstrated that using TOPMed sequencing data as the imputation reference panel improves genotype imputation quality in these populations, which subsequently enhanced gene-mapping power for complex traits. For rare variants with minor allele frequency (MAF) < 0.5%, we observed a 2.3- to 6.1-fold increase in the number of well-imputed variants, with 11-34% improvement in average imputation quality, compared to the state-of-the-art 1000 Genomes Project Phase 3 and Haplotype Reference Consortium reference panels. Impressively, even for extremely rare variants with minor allele count 86%. Subsequent association analyses of TOPMed reference panel-imputed genotype data with hematological traits (hemoglobin (HGB), hematocrit (HCT), and white blood cell count (WBC)) in ~21,600 African-ancestry and ~21,700 Hispanic/Latino individuals identified associations with two rare variants in the HBB gene (rs33930165 with higher WBC [p = 8.8x10-15] in African populations, rs11549407 with lower HGB [p = 1.5x10-12] and HCT [p = 8.8x10-10] in Hispanics/Latinos). By comparison, neither variant would have been genome-wide significant if either 1000 Genomes Project Phase 3 or Haplotype Reference Consortium reference panels had been used for imputation. Our findings highlight the utility of the TOPMed imputation reference panel for identification of novel rare variant associations not previously detected in similarly sized genome-wide studies of under-represented African and Hispanic/Latino populations.

Genome-wide association studies (GWAS) of neuroimaging data pose a significant computational burden because of the need to correct for multiple testing in both the genetic and the imaging data. Here, Ganjgahi et al. develop WLS-REML which significantly reduces computation running times in brain imaging GWAS.