Model-free Estimation of Recent Genetic Relatedness

Individual differences in internalizing symptoms in late childhood: A variance decomposition into cortical thickness, genetic and environmental differences

The brain undergoes extensive development during late childhood and early adolescence. Cortical thinning is a prominent feature of this development, and some researchers have suggested that differences in cortical thickness may be related to internalizing symptoms, which typically increase during the same period. However, research has yielded inconclusive results. We utilized a new method that estimates the combined effect of individual differences in vertex-wise cortical thickness on internalizing symptoms. This approach allows for many small effects to be distributed across the cortex and avoids the necessity of correcting for multiple tests. Using a sample of 8763 children aged 8.9 to 11.1 from the ABCD study, we decomposed the total variation in caregiver-reported internalizing symptoms into differences in cortical thickness, additive genetics, and shared family environmental factors and unique environmental factors. Our results indicated that individual differences in cortical thickness accounted for less than 0.5% of the variation in internalizing symptoms. In contrast, the analysis revealed a substantial effect of additive genetics and family environmental factors on the different components of internalizing symptoms, ranging from 06% to 48% and from 0% to 34%, respectively. Overall, while this study found a minimal association between cortical thickness and internalizing symptoms, additive genetics, and familial environmental factors appear to be of importance for describing differences in internalizing symptoms in late childhood. RESEARCH HIGHLIGHTS: We utilized a new method for modelling the total contribution of vertex-wise individual differences in cortical thickness to internalizing symptoms in late childhood. The total contribution of individual differences in cortical thickness accounted for <0.5% of the variance in internalizing symptoms. Additive genetics and shared family environmental variation accounted for 17% and 34% of the variance in internalizing symptoms, respectively. Our results suggest that cortical thickness is not an important indicator for internalizing symptoms in childhood, whereas genetic and environmental differences have a substantial impact.

Whole-genome sequencing study in Koreans identifies novel loci for Alzheimer's disease

INTRODUCTION

The genetic basis of Alzheimer's disease (AD) in Koreans is poorly understood.

METHODS

We performed an AD genome-wide association study using whole-genome sequence data from 3540 Koreans (1583 AD cases, 1957 controls) and single-nucleotide polymorphism array data from 2978 Japanese (1336 AD cases, 1642 controls). Significant findings were evaluated by pathway enrichment and differential gene expression analysis in brain tissue from controls and AD cases with and without dementia prior to death.

RESULTS

We identified genome-wide significant associations with APOE in the total sample and ROCK2 (rs76484417, p = 2.71×10-8) among APOE ε4 non-carriers. A study-wide significant association was found with aggregated rare variants in MICALL1 (MICAL like 1) (p = 9.04×10-7). Several novel AD-associated genes, including ROCK2 and MICALL1, were differentially expressed in AD cases compared to controls (p < 3.33×10-3). ROCK2 was also differentially expressed between AD cases with and without dementia (p = 1.34×10-4).

DISCUSSION

Our results provide insight into genetic mechanisms leading to AD and cognitive resilience in East Asians.

HIGHLIGHTS

Novel genome-wide significant associations for AD identified with ROCK2 and MICALL1. ROCK2 and MICALL1 are differentially expressed between AD cases and controls in the brain. This is the largest whole-genome-sequence study of AD in an East Asian population.

Secure discovery of genetic relatives across large-scale and distributed genomic data sets

Finding relatives within a study cohort is a necessary step in many genomic studies. However, when the cohort is distributed across multiple entities subject to data-sharing restrictions, performing this step often becomes infeasible. Developing a privacy-preserving solution for this task is challenging owing to the burden of estimating kinship between all the pairs of individuals across data sets. We introduce SF-Relate, a practical and secure federated algorithm for identifying genetic relatives across data silos. SF-Relate vastly reduces the number of individual pairs to compare while maintaining accurate detection through a novel locality-sensitive hashing (LSH) approach. We assign individuals who are likely to be related together into buckets and then test relationships only between individuals in matching buckets across parties. To this end, we construct an effective hash function that captures identity-by-descent (IBD) segments in genetic sequences, which, along with a new bucketing strategy, enable accurate and practical private relative detection. To guarantee privacy, we introduce an efficient algorithm based on multiparty homomorphic encryption (MHE) to allow data holders to cooperatively compute the relatedness coefficients between individuals and to further classify their degrees of relatedness, all without sharing any private data. We demonstrate the accuracy and practical runtimes of SF-Relate on the UK Biobank and All of Us data sets. On a data set of 200,000 individuals split between two parties, SF-Relate detects 97% of third-degree or closer relatives within 15 h of runtime. Our work enables secure identification of relatives across large-scale genomic data sets.

Rare variant contribution to the heritability of coronary artery disease

Whole genome sequences (WGS) enable discovery of rare variants which may contribute to missing heritability of coronary artery disease (CAD). To measure their contribution, we apply the GREML-LDMS-I approach to WGS of 4949 cases and 17,494 controls of European ancestry from the NHLBI TOPMed program. We estimate CAD heritability at 34.3% assuming a prevalence of 8.2%. Ultra-rare (minor allele frequency ≤ 0.1%) variants with low linkage disequilibrium (LD) score contribute ~50% of the heritability. We also investigate CAD heritability enrichment using a diverse set of functional annotations: i) constraint; ii) predicted protein-altering impact; iii) cis-regulatory elements from a cell-specific chromatin atlas of the human coronary; and iv) annotation principal components representing a wide range of functional processes. We observe marked enrichment of CAD heritability for most functional annotations. These results reveal the predominant role of ultra-rare variants in low LD on the heritability of CAD. Moreover, they highlight several functional processes including cell type-specific regulatory mechanisms as key drivers of CAD genetic risk.

Genetic analysis of cognitive preservation in the midwestern Amish reveals a novel locus on chromosome 2

INTRODUCTION

Alzheimer's disease (AD) remains a debilitating condition with limited treatments and additional therapeutic targets needed. Identifying AD protective genetic loci may identify new targets and accelerate identification of therapeutic treatments. We examined a founder population to identify loci associated with cognitive preservation into advanced age.

METHODS

Genome-wide association and linkage analyses were performed on 946 examined and sampled Amish individuals, aged 76-95, who were either cognitively unimpaired (CU) or impaired (CI).

RESULTS

A total of 12 single nucleotide polymorphisms (SNPs) demonstrated suggestive association (P ≤ 5 × 10-4) with cognitive preservation. Genetic linkage analyses identified > 100 significant (logarithm of the odds [LOD] ≥ 3.3) SNPs, some which overlapped with the association results. Only one locus on chromosome 2 retained significance across multiple analyses.

DISCUSSION

A novel significant result for cognitive preservation on chromosome 2 includes the genes LRRTM4 and CTNNA2. Additionally, the lead SNP, rs1402906, impacts the POU3F2 transcription factor binding affinity, which regulates LRRTM4 and CTNNA2.

HIGHLIGHTS

GWAS and linkage identified over 100 loci associated with cognitive preservation. One locus on Chromosome 2 retained significance over multiple analyses. Predicted TFBSs near rs1402906 regulate genes associated with neurocognition.

MethylGenotyper: Accurate Estimation of SNP Genotypes and Genetic Relatedness from DNA Methylation Data

Epigenome-wide association studies (EWAS) are susceptible to widespread confounding caused by population structure and genetic relatedness. Nevertheless, kinship estimation is challenging in EWAS without genotyping data. Here, we proposed MethylGenotyper, a method that for the first time enables accurate genotyping at thousands of single nucleotide polymorphisms (SNPs) directly from commercial DNA methylation microarrays. We modeled the intensities of methylation probes near SNPs with a mixture of three beta distributions corresponding to different genotypes and estimated parameters with an expectation-maximization algorithm. We conducted extensive simulations to demonstrate the performance of the method. When applying MethylGenotyper to the Infinium EPIC array data of 4662 Chinese samples, we obtained genotypes at 4319 SNPs with a concordance rate of 98.26%, enabling the identification of 255 pairs of close relatedness. Furthermore, we showed that MethylGenotyper allows for the estimation of both population structure and cryptic relatedness among 702 Australians of diverse ancestry. We also implemented MethylGenotyper in a publicly available R package (https://github.com/Yi-Jiang/MethylGenotyper) to facilitate future large-scale EWAS.

Obstructive sleep apnea mediates genetic risk of Diabetes Mellitus: The Hispanic Community Health Study/Study of Latinos

Objective

We sought to evaluate whether obstructive sleep apnea (OSA), and other sleep disorders, increase genetic risk of developing diabetes mellitus (DM).

Research Design and Methods

Using GWAS summary statistics from the DIAGRAM consortium and Million Veteran Program, we developed multi-ancestry Type 2 Diabetes (T2D) polygenic risk scores (T2D-PRSs) useful in admixed Hispanic/Latino individuals. We estimated the association of the T2D-PRS with cross-sectional and incident DM in the Hispanic Community Health Study/Study of Latinos (HCHS/SOL). We conducted a mediation analysis with T2D-PRSs as an exposure, incident DM as an outcome, and OSA as a mediator. Additionally, we performed Mendelian randomization (MR) analysis to assess the causal relationship between T2D and OSA.

Results

Of 12,342 HCHS/SOL participants, at baseline, 48.4% were normoglycemic, 36.6% were hyperglycemic, and 15% had diabetes, and 50.9% identified as female. Mean age was 41.5, and mean BMI was 29.4. T2D-PRSs was strongly associated with baseline DM and with incident DM. At baseline, a 1 SD increase in the primary T2D-PRS had DM adjusted odds ratio (OR) = 2.67, 95% CI [2.40; 2.97] and a higher incident DM rate (incident rate ratio (IRR) = 2.02, 95% CI [1.75; 2.33]). In a stratified analysis based on OSA severity categories the associations were stronger in individuals with mild OSA compared to those with moderate to severe OSA. Mediation analysis suggested that OSA mediates the T2D-PRS association with DM. In two-sample MR analysis, T2D-PRS had a causal effect on OSA, OR = 1.03, 95% CI [1.01; 1.05], and OSA had a causal effect on T2D, with OR = 2.34, 95% CI [1.59; 3.44].

Conclusions

OSA likely mediates genetic effects on T2D.

Analysis of C-reactive protein omics-measures associates methylation risk score with sleep health and related health outcomes

Introduction

DNA methylation (DNAm) predictors of high sensitivity C-reactive protein (CRP) offer a stable and accurate means of assessing chronic inflammation, bypassing the CRP protein fluctuations secondary to acute illness. Poor sleep health is associated with elevated inflammation (including elevated blood CRP levels) which may explain associations of sleep insufficiency with metabolic, cardiovascular and neurological diseases. Our study aims to characterize the relationships among sleep health phenotypes and CRP markers -blood, genetic, and epigenetic indicators-within the Hispanic Community Health Study/Study of Latinos (HCHS/SOL).

Methods

In HCHS/SOL, methylation risk scores (MRS)-CRP and polygenetic risk score (PRS)-CRP were constructed separately as weighted sums of methylation beta values or allele counts, respectively, for each individual. Sleep health phenotypes were measured using self-reported questionnaires and objective measurements. Survey-weighted linear regression established the association between the multiple sleep phenotypes (obstructive sleep apnea (OSA), sleep duration, insomnia and excessive sleepiness symptom), cognitive assessments, diabetes and hypertension with CRP markers while adjusting for age, sex, BMI, study center, and the first five principal components of genetic ancestry in HCHS/SOL.

Results

We included 2221 HCHS/SOL participants (age range 37-76 yrs, 65.7% female) in the analysis. Both the MRS-CRP (95% confidence interval (CI): 0.32-0.42, p = 3.3 × 10-38) and the PRS-CRP (95% CI: 0.15-0.25, p = 1 × 10-14) were associated with blood CRP level. Moreover, MRS-CRP was associated with sleep health phenotypes (OSA, long sleep duration) and related conditions (diabetes and hypertension), while PRS-CRP markers were not associated with these traits. Circulating CRP level was associated with sleep duration and diabetes. Associations between OSA traits and metabolic comorbidities weakened after adjusting for MRS-CRP, most strongly for diabetes, and least for hypertension.

Conclusions

MRS-CRP is a promising estimate for systemic and chronic inflammation as reflected by circulating CRP levels, which either mediates or serves as a common cause of the association between sleep phenotypes and related comorbidities, especially in the presence of diabetes.

Characterizing feral swine movement across the contiguous United States using neural networks and genetic data

Globalization has led to the frequent movement of species out of their native habitat. Some of these species become highly invasive and capable of profoundly altering invaded ecosystems. Feral swine (Sus scrofa × domesticus) are recognized as being among the most destructive invasive species, with populations established on all continents except Antarctica. Within the United States (US), feral swine are responsible for extensive crop damage, the destruction of native ecosystems, and the spread of disease. Purposeful human-mediated movement of feral swine has contributed to their rapid range expansion over the past 30 years. Patterns of deliberate introduction of feral swine have not been well described as populations may be established or augmented through small, undocumented releases. By leveraging an extensive genomic database of 18,789 samples genotyped at 35,141 single nucleotide polymorphisms (SNPs), we used deep neural networks to identify translocated feral swine across the contiguous US. We classified 20% (3364/16,774) of sampled animals as having been translocated and described general patterns of translocation using measures of centrality in a network analysis. These findings unveil extensive movement of feral swine well beyond their dispersal capabilities, including individuals with predicted origins >1000 km away from their sampling locations. Our study provides insight into the patterns of human-mediated movement of feral swine across the US and from Canada to the northern areas of the US. Further, our study validates the use of neural networks for studying the spread of invasive species.

Torque Teno Sus Virus 1: A Potential Surrogate Pathogen to Study Pig-Transmitted Transboundary Animal Diseases

Understanding the epidemiology and transmission dynamics of transboundary animal diseases (TADs) among wild pigs (Sus scrofa) will aid in preventing the introduction or containment of TADs among wild populations. Given the challenges associated with studying TADs in free-ranging populations, a surrogate pathogen system may predict how pathogens may circulate and be maintained within wild free-ranging swine populations, how they may spill over into domestic populations, and how management actions may impact transmission. We assessed the suitability of Torque teno sus virus 1 (TTSuV1) to serve as a surrogate pathogen for molecular epidemiological studies in wild pigs by investigating the prevalence, persistence, correlation with host health status and genetic variability at two study areas: Archbold's Buck Island Ranch in Florida and Savannah River Site in South Carolina. We then conducted a molecular epidemiological case study within Archbold's Buck Island Ranch site to determine how analysis of this pathogen could inform transmission dynamics of a directly transmitted virus. Prevalence was high in both study areas (40%, n = 190), and phylogenetic analyses revealed high levels of genetic variability within and between study areas. Our case study showed that pairwise host relatedness and geographic distance were highly correlated to pairwise viral genetic similarity. Molecular epidemiological analyses revealed a distinct pattern of direct transmission from pig to pig occurring within and between family groups. Our results suggest that TTSuV1 is highly suitable for molecular epidemiological analyses and will be useful for future studies of transmission dynamics in wild free-ranging pigs.

It's all relative: A multi-generational study using ForenSeq™ Kintelligence

The successful application of Forensic Investigative Genetic Genealogy (FIGG) to the identification of unidentified human remains and perpetrators of serious crime has led to a growing interest in its use internationally, including Australia. Routinely, FIGG has relied on the generation of high-density single nucleotide polymorphism (SNP) profiles from forensic samples using whole genome array (WGA) (∼650,000 or more SNPs) or whole genome sequencing (WGS) (millions of SNPs) for DNA segment-based comparisons in commercially available genealogy databases. To date, this approach has required DNA of a quality and quantity that is often not compatible with forensic samples. Furthermore, it requires the management of large data sets that include SNPs of medical relevance. The ForenSeq™ Kintelligence kit, comprising of 10,230 SNPs including 9867 for kinship association, was designed to overcome these challenges using a targeted amplicon sequencing-based method developed for low DNA inputs, inhibited and/or degraded forensic samples. To assess the ability of the ForenSeq™ Kintelligence workflow to correctly predict biological relationships, a comparative study comprising of 12 individuals from a family (with varying degrees of relatedness from 1st to 6th degree relatives) was undertaken using ForenSeq™ Kintelligence and a WGA approach using the Illumina Global Screening Array-24 version 3.0 Beadchip. All expected 1st, 2nd, 3rd, 4th and 5th degree relationships were correctly predicted using ForenSeq™ Kintelligence, while the expected 6th degree relationships were not detected. Given the (often) limited availability of forensic samples, findings from this study will assist Australian Law enforcement and other agencies considering the use of FIGG, to determine if the ForenSeq™ Kintelligence is suitable for existing workflows and casework sample types considered for FIGG.

Fetal Gene Regulatory Gene Deletions are Associated with Poor Cognition and Altered Cortical Morphology in Schizophrenia and Community-Based Samples

Schizophrenia spectrum disorders (SSDs) are characterized by substantial clinical and genetic heterogeneity. Multiple recurrent copy number variants (CNVs) increase risk for SSDs; however, how known risk CNVs and broader genome-wide CNVs influence clinical variability is unclear. The current study examined associations between borderline intellectual functioning or childhood-onset psychosis, known risk CNVs, and burden of deletions affecting genes in 18 previously validated neurodevelopmental gene-sets in 618 SSD individuals. CNV associations were assessed for replication in 235 SSD relatives and 583 controls, and 9,930 youth from the Adolescent Brain Cognitive Development (ABCD) Study. Known SSD- and neurodevelopmental disorder (NDD)-risk CNVs were associated with borderline intellectual functioning in SSD cases (odds ratios (OR) = 7.09 and 4.57, respectively); NDD-risk deletions were nominally associated with childhood-onset psychosis (OR = 4.34). Furthermore, deletion of genes involved in regulating gene expression during fetal brain development was associated with borderline intellectual functioning across SSD cases and non-cases (OR = 2.58), with partial replication in the ABCD cohort. Exploratory analyses of cortical morphology showed associations between fetal gene regulatory gene deletions and altered gray matter volume and cortical thickness across cohorts. Results highlight contributions of known risk CNVs to phenotypic variability in SSD and the utility of a neurodevelopmental framework for identifying mechanisms that influence phenotypic variability in SSDs, as well as the broader population, with implications for personalized medicine approaches to care.

Rare Sequence Variation Underlying Suspected Familial Cerebral Small-Vessel Disease

BACKGROUND

Cerebral small-vessel disease (cSVD) is the leading monogenic cause of stroke. Despite genetic screening in routine diagnosis, many cases remain without a known causative variant. Using a cohort with suspected familial cSVD and whole-genome sequencing, we screened for variants in genes associated with monogenic cSVD and searched for novel variants associated with the disease.

METHODS AND RESULTS

Rare variants were identified in whole-genome sequencing data from the NBR (National Institute for Health Research BioResource Rare Disease) study. Pathogenic variants in known monogenic cSVD genes were identified. Gene-based burden tests and family analysis were performed to identify novel variants associated with familial cSVD. A total of 257 suspected cSVD cases (mean ± SD age, 56.2 ± 16.1 years), and 13 086 controls with other nonstroke diseases (5874 [44.9%] men) were studied. A total of 8.9% of the cases carried a variant in known cSVD genes. Excluding these known causes, 23.6% of unrelated subjects with cSVD carried predicted deleterious variants in the Genomics England gene panel, but no association was found with cSVD in burden tests. We identified potential associations with cSVD in noncoding genes, including RP4-568F9.3 (adjusted P = 7.1 × 10-25), RP3-466I7.1 (adjusted P = 8.9 × 10-16), and ZNF209P (adjusted P = 1.0 × 10-15), and matrisomal genes (adjusted P = 5.1 × 10-6), including FAM20C, INHA, LAMC1, and VWA5B2.

CONCLUSIONS

Predicted deleterious variants in known cSVD genes were present in 23.6% of unrelated cases with cSVD, but none of the genes were associated with the disease. Rare variants in noncoding and matrisomal genes could potentially contribute to cSVD development. These genes could play a role in tissue development and brain endothelial cell function. However, further studies are needed to confirm their pathophysiological roles.

Whole genome sequencing based analysis of inflammation biomarkers in the Trans-Omics for Precision Medicine (TOPMed) consortium

Inflammation biomarkers can provide valuable insight into the role of inflammatory processes in many diseases and conditions. Sequencing based analyses of such biomarkers can also serve as an exemplar of the genetic architecture of quantitative traits. To evaluate the biological insight, which can be provided by a multi-ancestry, whole-genome based association study, we performed a comprehensive analysis of 21 inflammation biomarkers from up to 38 465 individuals with whole-genome sequencing from the Trans-Omics for Precision Medicine (TOPMed) program (with varying sample size by trait, where the minimum sample size was n = 737 for MMP-1). We identified 22 distinct single-variant associations across 6 traits-E-selectin, intercellular adhesion molecule 1, interleukin-6, lipoprotein-associated phospholipase A2 activity and mass, and P-selectin-that remained significant after conditioning on previously identified associations for these inflammatory biomarkers. We further expanded upon known biomarker associations by pairing the single-variant analysis with a rare variant set-based analysis that further identified 19 significant rare variant set-based associations with 5 traits. These signals were distinct from both significant single variant association signals within TOPMed and genetic signals observed in prior studies, demonstrating the complementary value of performing both single and rare variant analyses when analyzing quantitative traits. We also confirm several previously reported signals from semi-quantitative proteomics platforms. Many of these signals demonstrate the extensive allelic heterogeneity and ancestry-differentiated variant-trait associations common for inflammation biomarkers, a characteristic we hypothesize will be increasingly observed with well-powered, large-scale analyses of complex traits.

Cross-platform proteomics signatures of extreme old age

In previous work, we used a SomaLogic platform targeting approximately 5000 proteins to generate a serum protein signature of centenarians that we validated in independent studies that used the same technology. We set here to validate and possibly expand the results by profiling the serum proteome of a subset of individuals included in the original study using liquid chromatography tandem mass spectrometry (LC-MS/MS). Following pre-processing, the LC-MS/MS data provided quantification of 398 proteins, with only 266 proteins shared by both platforms. At 1% FDR statistical significance threshold, the analysis of LC-MS/MS data detected 44 proteins associated with extreme old age, including 23 of the original analysis. To identify proteins for which associations between expression and extreme-old age were conserved across platforms, we performed inter-study conservation testing of the 266 proteins quantified by both platforms using a method that accounts for the correlation between the results. From these tests, a total of 80 proteins reached 5% FDR statistical significance, and 26 of these proteins had concordant pattern of gene expression in whole blood generated in an independent set. This signature of 80 proteins points to blood coagulation, IGF signaling, extracellular matrix (ECM) organization, and complement cascade as important pathways whose protein level changes provide evidence for age-related adjustments that distinguish centenarians from younger individuals. The comparison with blood transcriptomics also highlights a possible role for neutrophil degranulation in aging.

Obesity risk in young adults from the Jerusalem Perinatal Study (JPS): the contribution of polygenic risk and early life exposure

BACKGROUND/OBJECTIVES

The effects of early life exposures on offspring life-course health are well established. This study assessed whether adding early socio-demographic and perinatal variables to a model based on polygenic risk score (PRS) improves prediction of obesity risk.

METHODS

We used the Jerusalem Perinatal study (JPS) with data at birth and body mass index (BMI) and waist circumference (WC) measured at age 32. The PRS was constructed using over 2.1M common SNPs identified in genome-wide association study (GWAS) for BMI. Linear and logistic models were applied in a stepwise approach. We first examined the associations between genetic variables and obesity-related phenotypes (e.g., BMI and WC). Secondly, socio-demographic variables were added and finally perinatal exposures, such as maternal pre-pregnancy BMI (mppBMI) and gestational weight gain (GWG) were added to the model. Improvement in prediction of each step was assessed using measures of model discrimination (area under the curve, AUC), net reclassification improvement (NRI) and integrated discrimination improvement (IDI).

RESULTS

One standard deviation (SD) change in PRS was associated with a significant increase in BMI (β = 1.40) and WC (β = 2.45). These associations were slightly attenuated (13.7-14.2%) with the addition of early life exposures to the model. Also, higher mppBMI was associated with increased offspring BMI (β = 0.39) and WC (β = 0.79) (p < 0.001). For obesity (BMI ≥ 30) prediction, the addition of early socio-demographic and perinatal exposures to the PRS model significantly increased AUC from 0.69 to 0.73. At an obesity risk threshold of 15%, the addition of early socio-demographic and perinatal exposures to the PRS model provided a significant improvement in reclassification of obesity (NRI, 0.147; 95% CI 0.068-0.225).

CONCLUSIONS

Inclusion of early life exposures, such as mppBMI and maternal smoking, to a model based on PRS improves obesity risk prediction in an Israeli population-sample.

A harmonized public resource of deeply sequenced diverse human genomes

Underrepresented populations are often excluded from genomic studies owing in part to a lack of resources supporting their analyses. The 1000 Genomes Project (1kGP) and Human Genome Diversity Project (HGDP), which have recently been sequenced to high coverage, are valuable genomic resources because of the global diversity they capture and their open data sharing policies. Here, we harmonized a high-quality set of 4094 whole genomes from 80 populations in the HGDP and 1kGP with data from the Genome Aggregation Database (gnomAD) and identified over 153 million high-quality SNVs, indels, and SVs. We performed a detailed ancestry analysis of this cohort, characterizing population structure and patterns of admixture across populations, analyzing site frequency spectra, and measuring variant counts at global and subcontinental levels. We also show substantial added value from this data set compared with the prior versions of the component resources, typically combined via liftOver and variant intersection; for example, we catalog millions of new genetic variants, mostly rare, compared with previous releases. In addition to unrestricted individual-level public release, we provide detailed tutorials for conducting many of the most common quality-control steps and analyses with these data in a scalable cloud-computing environment and publicly release this new phased joint callset for use as a haplotype resource in phasing and imputation pipelines. This jointly called reference panel will serve as a key resource to support research of diverse ancestry populations.

yQTL Pipeline: A structured computational workflow for large scale quantitative trait loci discovery and downstream visualization

Quantitative trait loci (QTL) denote regions of DNA whose variation is associated with variations in quantitative traits. QTL discovery is a powerful approach to understand how changes in molecular and clinical phenotypes may be related to DNA sequence changes. However, QTL discovery analysis encompasses multiple analytical steps and the processing of multiple input files, which can be laborious, error prone, and hard to reproduce if performed manually. To facilitate and automate large-scale QTL analysis, we developed the yQTL Pipeline, where the 'y' indicates the dependent quantitative variable being modeled. Prior to the association test, the pipeline supports the calculation or the direct input of pre-defined genome-wide principal components and genetic relationship matrix when applicable. User-specified covariates can also be provided. Depending on whether familial relatedness exists among the subjects, genome-wide association tests will be performed using either a linear mixed-effect model or a linear model. The options to run an ANOVA model or testing the interaction with a covariate are also available. Using the workflow management tool Nextflow, the pipeline parallelizes the analysis steps to optimize run-time and ensure results reproducibility. In addition, a user-friendly R Shiny App is developed to facilitate result visualization. It can generate Manhattan and Miami plots of phenotype traits, genotype-phenotype boxplots, and trait-QTL connection networks. We applied the yQTL Pipeline to analyze metabolomics profiles of blood serum from the New England Centenarians Study (NECS) participants. A total of 9.1M SNPs and 1,052 metabolites across 194 participants were analyzed. Using a p-value cutoff 5e-8, we found 14,983 mQTLs associated with 312 metabolites. The built-in parallelization of our pipeline reduced the run time from ~90 min to ~26 min. Visualization using the R Shiny App revealed multiple mQTLs shared across multiple metabolites. The yQTL Pipeline is available with documentation on GitHub at https://github.com/montilab/yQTLpipeline.

Partitioning variance in cortical morphometry into genetic, environmental, and subject-specific components

The relative contributions of genetic variation and experience in shaping the morphology of the adolescent brain are not fully understood. Using longitudinal data from 11,665 subjects in the ABCD Study, we fit vertex-wise variance components including family effects, genetic effects, and subject-level effects using a computationally efficient framework. Variance in cortical thickness and surface area is largely attributable to genetic influence, whereas sulcal depth is primarily explained by subject-level effects. Our results identify areas with heterogeneous distributions of heritability estimates that have not been seen in previous work using data from cortical regions. We discuss the biological importance of subject-specific variance and its implications for environmental influences on cortical development and maturation.

Data Harmonization Guidelines to Combine Multi-platform Genomic Data from Admixed Populations and Boost Power in Genome-Wide Association Studies

Data harmonization involves combining data from multiple independent sources and processing the data to produce one uniform dataset. Merging separate genotypes or whole-genome sequencing datasets has been proposed as a strategy to increase the statistical power of association tests by increasing the effective sample size. However, data harmonization is not a widely adopted strategy due to the difficulties with merging data (including confounding produced by batch effects and population stratification). Detailed data harmonization protocols are scarce and are often conflicting. Moreover, data harmonization protocols that accommodate samples of admixed ancestry are practically non-existent. Existing data harmonization procedures must be modified to ensure the heterogeneous ancestry of admixed individuals is incorporated into additional downstream analyses without confounding results. Here, we propose a set of guidelines for merging multi-platform genetic data from admixed samples that can be adopted by any investigator with elementary bioinformatics experience. We have applied these guidelines to aggregate 1544 tuberculosis (TB) case-control samples from six separate in-house datasets and conducted a genome-wide association study (GWAS) of TB susceptibility. The GWAS performed on the merged dataset had improved power over analyzing the datasets individually and produced summary statistics free from bias introduced by batch effects and population stratification. © 2024 Wiley Periodicals LLC. Basic Protocol 1: Processing separate datasets comprising array genotype data Alternate Protocol 1: Processing separate datasets comprising array genotype and whole-genome sequencing data Alternate Protocol 2: Performing imputation using a local reference panel Basic Protocol 2: Merging separate datasets Basic Protocol 3: Ancestry inference using ADMIXTURE and RFMix Basic Protocol 4: Batch effect correction using pseudo-case-control comparisons.

Genetic Associations with Placental Proteins in Maternal Serum Identify Biomarkers for Hypertension in Pregnancy

Background

Preeclampsia is a complex syndrome that accounts for considerable maternal and perinatal morbidity and mortality. Despite its prevalence, no effective disease-modifying therapies are available. Maternal serum placenta-derived proteins have been in longstanding use as markers of risk for aneuploidy and placental dysfunction, but whether they have a causal contribution to preeclampsia is unknown.

Objective

We aimed to investigate the genetic regulation of serum placental proteins in early pregnancy and their potential causal links with preeclampsia and gestational hypertension.

Study design

This study used a nested case-control design with nulliparous women enrolled in the nuMoM2b study from eight clinical sites across the United States between 2010 and 2013. The first- and second-trimester serum samples were collected, and nine proteins were measured, including vascular endothelial growth factor (VEGF), placental growth factor, endoglin, soluble fms-like tyrosine kinase-1 (sFlt-1), a disintegrin and metalloproteinase domain-containing protein 12 (ADAM-12), pregnancy-associated plasma protein A, free beta-human chorionic gonadotropin, inhibin A, and alpha-fetoprotein. This study used genome-wide association studies to discern genetic influences on these protein levels, treating proteins as outcomes. Furthermore, Mendelian randomization was used to evaluate the causal effects of these proteins on preeclampsia and gestational hypertension, and their further causal relationship with long-term hypertension, treating proteins as exposures.

Results

A total of 2,352 participants were analyzed. We discovered significant associations between the pregnancy zone protein locus and concentrations of ADAM-12 (rs6487735, P= 3.03×10 -22 ), as well as between the vascular endothelial growth factor A locus and concentrations of both VEGF (rs6921438, P= 7.94×10 -30 ) and sFlt-1 (rs4349809, P= 2.89×10 -12 ). Our Mendelian randomization analyses suggested a potential causal association between first-trimester ADAM-12 levels and gestational hypertension (odds ratio=0.78, P= 8.6×10 -4 ). We also found evidence for a potential causal effect of preeclampsia (odds ratio=1.75, P =8.3×10 -3 ) and gestational hypertension (odds ratio=1.84, P =4.7×10 -3 ) during the index pregnancy on the onset of hypertension 2-7 years later. The additional mediation analysis indicated that the impact of ADAM-12 on postpartum hypertension could be explained in part by its indirect effect through gestational hypertension (mediated effect=-0.15, P= 0.03).

Conclusions

Our study discovered significant genetic associations with placental proteins ADAM-12, VEGF, and sFlt-1, offering insights into their regulation during pregnancy. Mendelian randomization analyses demonstrated evidence of potential causal relationships between the serum levels of placental proteins, particularly ADAM-12, and gestational hypertension, potentially informing future prevention and treatment investigations.

Machine learning models for predicting blood pressure phenotypes by combining multiple polygenic risk scores

We construct non-linear machine learning (ML) prediction models for systolic and diastolic blood pressure (SBP, DBP) using demographic and clinical variables and polygenic risk scores (PRSs). We developed a two-model ensemble, consisting of a baseline model, where prediction is based on demographic and clinical variables only, and a genetic model, where we also include PRSs. We evaluate the use of a linear versus a non-linear model at both the baseline and the genetic model levels and assess the improvement in performance when incorporating multiple PRSs. We report the ensemble model's performance as percentage variance explained (PVE) on a held-out test dataset. A non-linear baseline model improved the PVEs from 28.1 to 30.1% (SBP) and 14.3% to 17.4% (DBP) compared with a linear baseline model. Including seven PRSs in the genetic model computed based on the largest available GWAS of SBP/DBP improved the genetic model PVE from 4.8 to 5.1% (SBP) and 4.7 to 5% (DBP) compared to using a single PRS. Adding additional 14 PRSs computed based on two independent GWASs further increased the genetic model PVE to 6.3% (SBP) and 5.7% (DBP). PVE differed across self-reported race/ethnicity groups, with primarily all non-White groups benefitting from the inclusion of additional PRSs. In summary, non-linear ML models improves BP prediction in models incorporating diverse populations.

Polygenic architecture of brain structure and function, behaviors, and psychopathologies in children

The human brain undergoes structural and functional changes during childhood, a critical period in cognitive and behavioral development. Understanding the genetic architecture of the brain development in children can offer valuable insights into the development of the brain, cognition, and behaviors. Here, we integrated brain imaging-genetic-phenotype data from over 8,600 preadolescent children of diverse ethnic backgrounds using multivariate statistical techniques. We found a low-to-moderate level of SNP-based heritability in most IDPs, which is lower compared to the adult brain. Using sparse generalized canonical correlation analysis (SGCCA), we identified several covariation patterns among genome-wide polygenic scores (GPSs) of 29 traits, 7 different modalities of brain imaging-derived phenotypes (IDPs), and 266 cognitive and psychological phenotype data. In structural MRI, significant positive associations were observed between total grey matter volume, left ventral diencephalon volume, surface area of right accumbens and the GPSs of cognition-related traits. Conversely, negative associations were found with the GPSs of ADHD, depression and neuroticism. Additionally, we identified a significant positive association between educational attainment GPS and regional brain activation during the N-back task. The BMI GPS showed a positive association with fractional anisotropy (FA) of connectivity between the cerebellum cortex and amygdala in diffusion MRI, while the GPSs for educational attainment and cannabis use were negatively associated with the same IDPs. Our GPS-based prediction models revealed substantial genetic contributions to cognitive variability, while the genetic basis for many mental and behavioral phenotypes remained elusive. This study delivers a comprehensive map of the relationships between genetic profiles, neuroanatomical diversity, and the spectrum of cognitive and behavioral traits in preadolescence.

The impact of assortative mating, participation bias and socioeconomic status on the polygenic risk of behavioural and psychiatric traits

To investigate assortative mating (AM), participation bias and socioeconomic status (SES) with respect to the genetics of behavioural and psychiatric traits, we estimated AM signatures using gametic phase disequilibrium and within-spouses and within-siblings polygenic risk score correlation analyses, also performing a SES conditional analysis. The cross-method meta-analysis identified AM genetic signatures for multiple alcohol-related phenotypes, bipolar disorder, major depressive disorder, schizophrenia and Tourette syndrome. Here, after SES conditioning, we observed changes in the AM genetic signatures for maximum habitual alcohol intake, frequency of drinking alcohol and Tourette syndrome. We also observed significant gametic phase disequilibrium differences between UK Biobank mental health questionnaire responders versus non-responders for major depressive disorder and alcohol use disorder. These results highlight the impact of AM, participation bias and SES on the polygenic risk of behavioural and psychiatric traits, particularly in alcohol-related traits.

Cross-platform proteomics signatures of extreme old age

In previous work we used a Somalogic platform targeting approximately 5000 proteins to generate a serum protein signature of centenarians that we validated in independent studies that used the same technology. We set here to validate and possibly expand the results by profiling the serum proteome of a subset of individuals included in the original study using liquid chromatography tandem mass spectrometry (LC-MS/MS). Following pre-processing, the LC-MS/MS data provided quantification of 398 proteins, with only 266 proteins shared by both platforms. At 1% FDR statistical significance threshold, the analysis of LC-MS/MS data detected 44 proteins associated with extreme old age, including 23 of the original analysis. To identify proteins for which associations between expression and extreme-old age were conserved across platforms, we performed inter-study conservation testing of the 266 proteins quantified by both platforms using a method that accounts for the correlation between the results. From these tests, a total of 80 proteins reached 5% FDR statistical significance, and 26 of these proteins had concordant pattern of gene expression in whole blood. This signature of 80 proteins points to blood coagulation, IGF signaling, extracellular matrix (ECM) organization, and complement cascade as important pathways whose protein level changes provide evidence for age-related adjustments that distinguish centenarians from younger individuals.

Adjusting for principal components can induce spurious associations in genome-wide association studies in admixed populations

Principal component analysis (PCA) is widely used to control for population structure in genome-wide association studies (GWAS). Top principal components (PCs) typically reflect population structure, but challenges arise in deciding how many PCs are needed and ensuring that PCs do not capture other artifacts such as regions with atypical linkage disequilibrium (LD). In response to the latter, many groups suggest performing LD pruning or excluding known high LD regions prior to PCA. However, these suggestions are not universally implemented and the implications for GWAS are not fully understood, especially in the context of admixed populations. In this paper, we investigate the impact of pre-processing and the number of PCs included in GWAS models in African American samples from the Women's Women's Health Initiative SNP Health Association Resource and two Trans-Omics for Precision Medicine Whole Genome Sequencing Project contributing studies (Jackson Heart Study and Genetic Epidemiology of Chronic Obstructive Pulmonary Disease Study). In all three samples, we find the first PC is highly correlated with genome-wide ancestry whereas later PCs often capture local genomic features. The pattern of which, and how many, genetic variants are highly correlated with individual PCs differs from what has been observed in prior studies focused on European populations and leads to distinct downstream consequences: adjusting for such PCs yields biased effect size estimates and elevated rates of spurious associations due to the phenomenon of collider bias. Excluding high LD regions identified in previous studies does not resolve these issues. LD pruning proves more effective, but the optimal choice of thresholds varies across datasets. Altogether, our work highlights unique issues that arise when using PCA to control for ancestral heterogeneity in admixed populations and demonstrates the importance of careful pre-processing and diagnostics to ensure that PCs capturing multiple local genomic features are not included in GWAS models.

Integrative Analysis of Germline Rare Variants in Clear and Non-clear Cell Renal Cell Carcinoma

Background and objective

Previous germline studies on renal cell carcinoma (RCC) have usually pooled clear and non-clear cell RCCs and have not adequately accounted for population stratification, which might have led to an inaccurate estimation of genetic risk. Here, we aim to analyze the major germline drivers of RCC risk and clinically relevant but underexplored germline variant types.

Methods

We first characterized germline pathogenic variants (PVs), cryptic splice variants, and copy number variants (CNVs) in 1436 unselected RCC patients. To evaluate the enrichment of PVs in RCC, we conducted a case-control study of 1356 RCC patients ancestry matched with 16 512 cancer-free controls using approaches accounting for population stratification and histological subtypes, followed by characterization of secondary somatic events.

Key findings and limitations

Clear cell RCC patients (n = 976) exhibited a significant burden of PVs in VHL compared with controls (odds ratio [OR]: 39.1, p = 4.95e-05). Non-clear cell RCC patients (n = 380) carried enrichment of PVs in FH (OR: 77.9, p = 1.55e-08) and MET (OR: 1.98e11, p = 2.07e-05). In a CHEK2-focused analysis with European participants, clear cell RCC (n = 906) harbored nominal enrichment of low-penetrance CHEK2 variants-p.Ile157Thr (OR: 1.84, p = 0.049) and p.Ser428Phe (OR: 5.20, p = 0.045), while non-clear cell RCC (n = 295) exhibited nominal enrichment of CHEK2 loss of function PVs (OR: 3.51, p = 0.033). Patients with germline PVs in FH, MET, and VHL exhibited significantly earlier age of cancer onset than patients without germline PVs (mean: 46.0 vs 60.2 yr, p < 0.0001), and more than half had secondary somatic events affecting the same gene (n = 10/15, 66.7%). Conversely, CHEK2 PV carriers exhibited a similar age of onset to patients without germline PVs (mean: 60.1 vs 60.2 yr, p = 0.99), and only 30.4% carried somatic events in CHEK2 (n = 7/23). Finally, pathogenic germline cryptic splice variants were identified in SDHA and TSC1, and pathogenic germline CNVs were found in 18 patients, including CNVs in FH, SDHA, and VHL.

Conclusions and clinical implications

This analysis supports the existing link between several RCC risk genes and RCC risk manifesting in earlier age of onset. It calls for caution when assessing the role of CHEK2 due to the burden of founder variants with varying population frequency. It also broadens the definition of the RCC germline landscape of pathogenicity to incorporate previously understudied types of germline variants.

Patient summary

In this study, we carefully compared the frequency of rare inherited mutations with a focus on patients' genetic ancestry. We discovered that subtle variations in genetic background may confound a case-control analysis, especially in evaluating the cancer risk associated with specific genes, such as CHEK2. We also identified previously less explored forms of rare inherited mutations, which could potentially increase the risk of kidney cancer.

Gene-environment pathways to cognitive intelligence and psychotic-like experiences in children

In children, psychotic-like experiences (PLEs) are related to risk of psychosis, schizophrenia, and other mental disorders. Maladaptive cognitive functioning, influenced by genetic and environmental factors, is hypothesized to mediate the relationship between these factors and childhood PLEs. Using large-scale longitudinal data, we tested the relationships of genetic and environmental factors (such as familial and neighborhood environment) with cognitive intelligence and their relationships with current and future PLEs in children. We leveraged large-scale multimodal data of 6,602 children from the Adolescent Brain and Cognitive Development Study. Linear mixed model and a novel structural equation modeling (SEM) method that allows estimation of both components and factors were used to estimate the joint effects of cognitive phenotypes polygenic scores (PGSs), familial and neighborhood socioeconomic status (SES), and supportive environment on NIH Toolbox cognitive intelligence and PLEs. We adjusted for ethnicity (genetically defined), schizophrenia PGS, and additionally unobserved confounders (using computational confound modeling). Our findings indicate that lower cognitive intelligence and higher PLEs are significantly associated with lower PGSs for cognitive phenotypes, lower familial SES, lower neighborhood SES, and less supportive environments. Specifically, cognitive intelligence mediates the effects of these factors on PLEs, with supportive parenting and positive school environments showing the strongest impact on reducing PLEs. This study underscores the influence of genetic and environmental factors on PLEs through their effects on cognitive intelligence. Our findings have policy implications in that improving school and family environments and promoting local economic development may enhance cognitive and mental health in children.

A harmonized public resource of deeply sequenced diverse human genomes

Underrepresented populations are often excluded from genomic studies due in part to a lack of resources supporting their analyses. The 1000 Genomes Project (1kGP) and Human Genome Diversity Project (HGDP), which have recently been sequenced to high coverage, are valuable genomic resources because of the global diversity they capture and their open data sharing policies. Here, we harmonized a high quality set of 4,094 whole genomes from HGDP and 1kGP with data from the Genome Aggregation Database (gnomAD) and identified over 153 million high-quality SNVs, indels, and SVs. We performed a detailed ancestry analysis of this cohort, characterizing population structure and patterns of admixture across populations, analyzing site frequency spectra, and measuring variant counts at global and subcontinental levels. We also demonstrate substantial added value from this dataset compared to the prior versions of the component resources, typically combined via liftover and variant intersection; for example, we catalog millions of new genetic variants, mostly rare, compared to previous releases. In addition to unrestricted individual-level public release, we provide detailed tutorials for conducting many of the most common quality control steps and analyses with these data in a scalable cloud-computing environment and publicly release this new phased joint callset for use as a haplotype resource in phasing and imputation pipelines. This jointly called reference panel will serve as a key resource to support research of diverse ancestry populations.

How Can Genomics Help or Hinder Wildlife Conservation?

Genomic data are becoming increasingly affordable and easy to collect, and new tools for their analysis are appearing rapidly. Conservation biologists are interested in using this information to assist in management and planning but are typically limited financially and by the lack of genomic resources available for non-model taxa. It is therefore important to be aware of the pitfalls as well as the benefits of applying genomic approaches. Here, we highlight recent methods aimed at standardizing population assessments of genetic variation, inbreeding, and forms of genetic load and methods that help identify past and ongoing patterns of genetic interchange between populations, including those subjected to recent disturbance. We emphasize challenges in applying some of these methods and the need for adequate bioinformatic support. We also consider the promises and challenges of applying genomic approaches to understand adaptive changes in natural populations to predict their future adaptive capacity.

Intracranial Volume Is Driven by Both Genetics and Early Life Exposures: The SOL-INCA-MRI Study

Intracranial volume (ICV) reflects maximal brain development and is associated with later-life cognitive abilities. We quantified ICV among first- and second-generation Hispanic and Latino adults from the Study of Latinos-Investigation of Cognitive Aging - MRI (SOL-INCA-MRI), estimated ICV heritability, and tested its associations with previously reported genetic variants, both individually and as a genetic risk score (GRS). We also estimated the association of ICV with early life environmental measures: nativity or age of immigration and parental education. The estimated heritability of ICV was 19% (95% CI, 0.1%-56%) in n=1781 unrelated SOL-INCA-MRI individuals. Four of 10 tested genetic variants were associated with ICV and an increase of 1 SD of the ICV-GRS was associated with an increase of 10.37 cm3 in the ICV (95% CI, 5.29-15.45). Compared to being born in the continental United States, immigrating to the United States at age 11 years or older was associated with 24 cm3 smaller ICV (95% CI, -39.97 to -8.06). Compared to both parents having less than high-school education, at least 1 parent completing high-school education was associated with 15.4 cm3 greater ICV (95% CI, 4.46-26.39). These data confirm the importance of early life health on brain development.

The genetic legacy of the first successful reintroduction of a mammal to Britain: Founder events and attempted genetic rescue in Scotland's beaver population

Conservation translocations often inherently involve a risk of genetic diversity loss, and thus loss of adaptive potential, but this risk is rarely quantified or monitored through time. The reintroduction of beavers to Scotland, via the Scottish Beaver Trial in Knapdale, is an example of a translocation that took place in the absence of genetic data for the founder individuals and resulted in a small and suspected to be genetically depauperate population. In this study we use a high-density SNP panel to assess the genetic impact of that initial translocation and the effect of subsequent reinforcement translocations using animals from a different genetic source to the original founders. We demonstrate that the initial translocation did, indeed, lead to low genetic diversity (H o = 0.052) and high mean kinship (KING-robust = 0.159) in the Knapdale population compared to other beaver populations. We also show that the reinforcement translocations have succeeded in increasing genetic diversity (H o = 0.196) and reducing kinship (KING robust = 0.028) in Knapdale. As yet, there is no evidence of admixture between the two genetic lineages that are now present in Knapdale and such admixture is necessary to realise the full genetic benefits of the reinforcement and for genetic reinforcement and then rescue to occur; future genetic monitoring will be required to assess whether this has happened. We note that, should admixture occur, the Knapdale population will harbour combinations of genetic diversity not currently seen elsewhere in Eurasian beavers, posing important considerations for the future management of this population. We consider our results in the wider context of beaver conservation throughout Scotland and the rest of Britain, and advocate for more proactive genetic sampling of all founders to allow the full integration of genetic data into translocation planning in general.

Genetic and brain similarity independently predict childhood anthropometrics and neighborhood socioeconomic conditions

Linking the developing brain with individual differences in clinical and demographic traits is challenging due to the substantial interindividual heterogeneity of brain anatomy and organization. Here we employ an integrative approach that parses individual differences in both cortical thickness and common genetic variants, and assess their effects on a wide set of childhood traits. The approach uses a linear mixed model framework to obtain the unique effects of each type of similarity, as well as their covariance. We employ this approach in a sample of 7760 unrelated children in the ABCD cohort baseline sample (mean age 9.9, 46.8% female). In general, associations between cortical thickness similarity and traits were limited to anthropometrics such as height, weight, and birth weight, as well as a marker of neighborhood socioeconomic conditions. Common genetic variants explained significant proportions of variance across nearly all included outcomes, although estimates were somewhat lower than previous reports. No significant covariance of the effects of genetic and cortical thickness similarity was found. The present findings highlight the connection between anthropometrics as well as neighborhood socioeconomic conditions and the developing brain, which appear to be independent from individual differences in common genetic variants in this population-based sample.

yQTL Pipeline: a structured computational workflow for large scale quantitative trait loci discovery and downstream visualization

Quantitative trait loci (QTL) denote regions of DNA whose variation is associated with variations in quantitative traits. QTL discovery is a powerful approach to understand how changes in molecular and clinical phenotypes may be related to DNA sequence changes. However, QTL discovery analysis encompasses multiple analytical steps and the processing of multiple input files, which can be laborious, error prone, and hard to reproduce if performed manually. In order to facilitate and automate large-scale QTL analysis, we developed the yQTL Pipeline, where the 'y' indicates the dependent quantitative variable being modeled. Prior to genome-wide association test, the pipeline supports the calculation or the direct input of pre-defined genome-wide principal components and genetic relationship matrix when applicable. User-specified covariates can also be provided. Depending on whether familial relatedness exists among the subjects, genome-wide association tests will be performed using either a linear mixed-effect model or a linear model. Using the workflow management tool Nextflow, the pipeline parallelizes the analysis steps to optimize run-time and ensure results reproducibility. In addition, a user-friendly R Shiny App is developed to facilitate result visualization. Upon uploading the result file, it can generate Manhattan plots of user-selected phenotype traits and trait-QTL connection networks based on user-specified p-value thresholds. We applied the yQTL Pipeline to analyze metabolomics profiles of blood serum from the New England Centenarians Study (NECS) participants. A total of 9.1M SNPs and 1,052 metabolites across 194 participants were analyzed. Using a p-value cutoff 5e-8, we found 14,983 mQTLs cumulatively associated with 312 metabolites. The built-in parallelization of our pipeline reduced the run time from ~90 min to ~26 min. Visualization using the R Shiny App revealed multiple mQTLs shared across multiple metabolites. The yQTL Pipeline is available with documentation on GitHub at https://github.com/montilab/yQTL-Pipeline.

Early life stress modulates the genetic influence on brain structure and cognitive function in children

The enduring influence of early life stress (ELS) on brain and cognitive development has been widely acknowledged, yet the precise mechanisms underlying this association remain elusive. We hypothesize that ELS might disrupt the genome-wide influence on brain morphology and connectivity development, consequently exerting a detrimental impact on children's cognitive ability. We analyzed the multimodal data of DNA genotypes, brain imaging (structural and diffusion MRI), and neurocognitive battery (NIH Toolbox) of 4276 children (ages 9-10 years, European ancestry) from the Adolescent Brain Cognitive Development (ABCD) study. The genome-wide influence on cognitive function was estimated using the polygenic score (GPS). By using brain morphometry and tractography, we identified the brain correlates of the cognition GPSs. Statistical analyses revealed relationships for the gene-brain-cognition pathway. The brain structural variance significantly mediated the genetic influence on cognition (indirect effect = 0.016, PFDR < 0.001). Of note, this gene-brain relationship was significantly modulated by abuse, resulting in diminished cognitive capacity (Index of Moderated Mediation = -0.007; 95 % CI = -0.012 ∼ -0.002). Our results support a novel gene-brain-cognition model likely elucidating the long-lasting negative impact of ELS on children's cognitive development.

Building a vertically integrated genomic learning health system: The biobank at the Colorado Center for Personalized Medicine

Precision medicine initiatives across the globe have led to a revolution of repositories linking large-scale genomic data with electronic health records, enabling genomic analyses across the entire phenome. Many of these initiatives focus solely on research insights, leading to limited direct benefit to patients. We describe the biobank at the Colorado Center for Personalized Medicine (CCPM Biobank) that was jointly developed by the University of Colorado Anschutz Medical Campus and UCHealth to serve as a unique, dual-purpose research and clinical resource accelerating personalized medicine. This living resource currently has more than 200,000 participants with ongoing recruitment. We highlight the clinical, laboratory, regulatory, and HIPAA-compliant informatics infrastructure along with our stakeholder engagement, consent, recontact, and participant engagement strategies. We characterize aspects of genetic and geographic diversity unique to the Rocky Mountain region, the primary catchment area for CCPM Biobank participants. We leverage linked health and demographic information of the CCPM Biobank participant population to demonstrate the utility of the CCPM Biobank to replicate complex trait associations in the first 33,674 genotyped individuals across multiple disease domains. Finally, we describe our current efforts toward return of clinical genetic test results, including high-impact pathogenic variants and pharmacogenetic information, and our broader goals as the CCPM Biobank continues to grow. Bringing clinical and research interests together fosters unique clinical and translational questions that can be addressed from the large EHR-linked CCPM Biobank resource within a HIPAA- and CLIA-certified environment.

Genetic ancestry and diagnostic yield of exome sequencing in a diverse population

It has been suggested that diagnostic yield (DY) from Exome Sequencing (ES) may be lower among patients with non-European ancestries than those with European ancestry. We examined the association of DY with estimated continental/subcontinental genetic ancestry in a racially/ethnically diverse pediatric and prenatal clinical cohort. Cases (N = 845) with suspected genetic disorders underwent ES for diagnosis. Continental/subcontinental genetic ancestry proportions were estimated from the ES data. We compared the distribution of genetic ancestries in positive, negative, and inconclusive cases by Kolmogorov-Smirnov tests and linear associations of ancestry with DY by Cochran-Armitage trend tests. We observed no reduction in overall DY associated with any genetic ancestry (African, Native American, East Asian, European, Middle Eastern, South Asian). However, we observed a relative increase in proportion of autosomal recessive homozygous inheritance versus other inheritance patterns associated with Middle Eastern and South Asian ancestry, due to consanguinity. In this empirical study of ES for undiagnosed pediatric and prenatal genetic conditions, genetic ancestry was not associated with the likelihood of a positive diagnosis, supporting the equitable use of ES in diagnosis of previously undiagnosed but potentially Mendelian disorders across all ancestral populations.

Genetic analysis of cognitive preservation in the midwestern Amish reveals a novel locus on chromosome 2

INTRODUCTION

Alzheimer disease (AD) remains a debilitating condition with limited treatments and additional therapeutic targets needed. Identifying AD protective genetic loci may identify new targets and accelerate identification of therapeutic treatments. We examined a founder population to identify loci associated with cognitive preservation into advanced age.

METHODS

Genome-wide association and linkage analyses were performed on 946 examined and sampled Amish individuals, aged 76-95, who were either cognitively unimpaired (CU) or impaired (CI).

RESULTS

12 SNPs demonstrated suggestive association (P≤5×10-4) with cognitive preservation. Genetic linkage analyses identified >100 significant (LOD≥3.3) SNPs, some which overlapped with the association results. Only one locus on chromosome 2 retained significance across multiple analyses.

DISCUSSION

A novel significant result for cognitive preservation on chromosome 2 includes the genes LRRTM4 and CTNNA2. Additionally, the lead SNP, rs1402906, impacts the POU3F2 transcription factor binding affinity, which regulates LRRTM4 and CTNNA2.

Machine learning models for blood pressure phenotypes combining multiple polygenic risk scores

We construct non-linear machine learning (ML) prediction models for systolic and diastolic blood pressure (SBP, DBP) using demographic and clinical variables and polygenic risk scores (PRSs). We developed a two-model ensemble, consisting of a baseline model, where prediction is based on demographic and clinical variables only, and a genetic model, where we also include PRSs. We evaluate the use of a linear versus a non-linear model at both the baseline and the genetic model levels and assess the improvement in performance when incorporating multiple PRSs. We report the ensemble model's performance as percentage variance explained (PVE) on a held-out test dataset. A non-linear baseline model improved the PVEs from 28.1% to 30.1% (SBP) and 14.3% to 17.4% (DBP) compared with a linear baseline model. Including seven PRSs in the genetic model computed based on the largest available GWAS of SBP/DBP improved the genetic model PVE from 4.8% to 5.1% (SBP) and 4.7% to 5% (DBP) compared to using a single PRS. Adding additional 14 PRSs computed based on two independent GWASs further increased the genetic model PVE to 6.3% (SBP) and 5.7% (DBP). PVE differed across self-reported race/ethnicity groups, with primarily all non-White groups benefitting from the inclusion of additional PRSs.

Whole-genome sequencing uncovers two loci for coronary artery calcification and identifies ARSE as a regulator of vascular calcification

Coronary artery calcification (CAC) is a measure of atherosclerosis and a well-established predictor of coronary artery disease (CAD) events. Here we describe a genome-wide association study (GWAS) of CAC in 22,400 participants from multiple ancestral groups. We confirmed associations with four known loci and identified two additional loci associated with CAC (ARSE and MMP16), with evidence of significant associations in replication analyses for both novel loci. Functional assays of ARSE and MMP16 in human vascular smooth muscle cells (VSMCs) demonstrate that ARSE is a promoter of VSMC calcification and VSMC phenotype switching from a contractile to a calcifying or osteogenic phenotype. Furthermore, we show that the association of variants near ARSE with reduced CAC is likely explained by reduced ARSE expression with the G allele of enhancer variant rs5982944. Our study highlights ARSE as an important contributor to atherosclerotic vascular calcification, and a potential drug target for vascular calcific disease.

Psychosis Endophenotypes: A Gene-Set-Specific Polygenic Risk Score Analysis

BACKGROUND AND HYPOTHESIS

Endophenotypes can help to bridge the gap between psychosis and its genetic predispositions, but their underlying mechanisms remain largely unknown. This study aims to identify biological mechanisms that are relevant to the endophenotypes for psychosis, by partitioning polygenic risk scores into specific gene sets and testing their associations with endophenotypes.

STUDY DESIGN

We computed polygenic risk scores for schizophrenia and bipolar disorder restricted to brain-related gene sets retrieved from public databases and previous publications. Three hundred and seventy-eight gene-set-specific polygenic risk scores were generated for 4506 participants. Seven endophenotypes were also measured in the sample. Linear mixed-effects models were fitted to test associations between each endophenotype and each gene-set-specific polygenic risk score.

STUDY RESULTS

After correction for multiple testing, we found that a reduced P300 amplitude was associated with a higher schizophrenia polygenic risk score of the forebrain regionalization gene set (mean difference per SD increase in the polygenic risk score: -1.15 µV; 95% CI: -1.70 to -0.59 µV; P = 6 × 10-5). The schizophrenia polygenic risk score of forebrain regionalization also explained more variance of the P300 amplitude (R2 = 0.032) than other polygenic risk scores, including the genome-wide polygenic risk scores.

CONCLUSIONS

Our finding on reduced P300 amplitudes suggests that certain genetic variants alter early brain development thereby increasing schizophrenia risk years later. Gene-set-specific polygenic risk scores are a useful tool to elucidate biological mechanisms of psychosis and endophenotypes, offering leads for experimental validation in cellular and animal models.

Review: Computational analysis of human skeletal remains in ancient DNA and forensic genetics

Degraded DNA is used to answer questions in the fields of ancient DNA (aDNA) and forensic genetics. While aDNA studies typically center around human evolution and past history, and forensic genetics is often more concerned with identifying a specific individual, scientists in both fields face similar challenges. The overlap in source material has prompted periodic discussions and studies on the advantages of collaboration between fields toward mutually beneficial methodological advancements. However, most have been centered around wet laboratory methods (sampling, DNA extraction, library preparation, etc.). In this review, we focus on the computational side of the analytical workflow. We discuss limitations and considerations to consider when working with degraded DNA. We hope this review provides a framework to researchers new to computational workflows for how to think about analyzing highly degraded DNA and prompts an increase of collaboration between the forensic genetics and aDNA fields.

A statistical framework for powerful multi-trait rare variant analysis in large-scale whole-genome sequencing studies

Large-scale whole-genome sequencing (WGS) studies have improved our understanding of the contributions of coding and noncoding rare variants to complex human traits. Leveraging association effect sizes across multiple traits in WGS rare variant association analysis can improve statistical power over single-trait analysis, and also detect pleiotropic genes and regions. Existing multi-trait methods have limited ability to perform rare variant analysis of large-scale WGS data. We propose MultiSTAAR, a statistical framework and computationally-scalable analytical pipeline for functionally-informed multi-trait rare variant analysis in large-scale WGS studies. MultiSTAAR accounts for relatedness, population structure and correlation among phenotypes by jointly analyzing multiple traits, and further empowers rare variant association analysis by incorporating multiple functional annotations. We applied MultiSTAAR to jointly analyze three lipid traits (low-density lipoprotein cholesterol, high-density lipoprotein cholesterol and triglycerides) in 61,861 multi-ethnic samples from the Trans-Omics for Precision Medicine (TOPMed) Program. We discovered new associations with lipid traits missed by single-trait analysis, including rare variants within an enhancer of NIPSNAP3A and an intergenic region on chromosome 1.

Mosaic chromosomal alterations in blood across ancestries using whole-genome sequencing

Megabase-scale mosaic chromosomal alterations (mCAs) in blood are prognostic markers for a host of human diseases. Here, to gain a better understanding of mCA rates in genetically diverse populations, we analyzed whole-genome sequencing data from 67,390 individuals from the National Heart, Lung, and Blood Institute Trans-Omics for Precision Medicine program. We observed higher sensitivity with whole-genome sequencing data, compared with array-based data, in uncovering mCAs at low mutant cell fractions and found that individuals of European ancestry have the highest rates of autosomal mCAs and the lowest rates of chromosome X mCAs, compared with individuals of African or Hispanic ancestry. Although further studies in diverse populations will be needed to replicate our findings, we report three loci associated with loss of chromosome X, associations between autosomal mCAs and rare variants in DCPS, ADM17, PPP1R16B and TET2 and ancestry-specific variants in ATM and MPL with mCAs in cis.

BRASS: Permutation methods for binary traits in genetic association studies with structured samples

In genetic association analysis of complex traits, permutation testing can be a valuable tool for assessing significance when the distribution of the test statistic is unknown or not well-approximated. This commonly arises, e.g, in tests of gene-set, pathway or genome-wide significance, or when the statistic is formed by machine learning or data adaptive methods. Existing applications include eQTL mapping, association testing with rare variants, inclusion of admixed individuals in genetic association analysis, and epistasis detection among many others. For genetic association testing in samples with population structure and/or relatedness, use of naive permutation can lead to inflated type 1 error. To address this in quantitative traits, the MVNpermute method was developed. However, for association mapping of a binary trait, the relationship between the mean and variance makes both naive permutation and the MVNpermute method invalid. We propose BRASS, a permutation method for binary traits, for use in association mapping in structured samples. In addition to modeling structure in the sample, BRASS allows for covariates, ascertainment and simultaneous testing of multiple markers, and it accommodates a wide range of test statistics. In simulation studies, we compare BRASS to other permutation and resampling-based methods in a range of scenarios that include population structure, familial relatedness, ascertainment and phenotype model misspecification. In these settings, we demonstrate the superior control of type 1 error by BRASS compared to the other 6 methods considered. We apply BRASS to assess genome-wide significance for association analyses in domestic dog for elbow dysplasia (ED) and idiopathic epilepsy (IE). For both traits we detect previously identified associations, and in addition, for ED, we detect significant association with a SNP on chromosome 35 that was not detected by previous analyses, demonstrating the potential of the method.

Decoding Genetics, Ancestry, and Geospatial Context for Precision Health

Mass General Brigham, an integrated healthcare system based in the Greater Boston area of Massachusetts, annually serves 1.5 million patients. We established the Mass General Brigham Biobank (MGBB), encompassing 142,238 participants, to unravel the intricate relationships among genomic profiles, environmental context, and disease manifestations within clinical practice. In this study, we highlight the impact of ancestral diversity in the MGBB by employing population genetics, geospatial assessment, and association analyses of rare and common genetic variants. The population structures captured by the genetics mirror the sequential immigration to the Greater Boston area throughout American history, highlighting communities tied to shared genetic and environmental factors. Our investigation underscores the potency of unbiased, large-scale analyses in a healthcare-affiliated biobank, elucidating the dynamic interplay across genetics, immigration, structural geospatial factors, and health outcomes in one of the earliest American sites of European colonization.

Rare variants in CAPN2 increase risk for isolated hypoplastic left heart syndrome

Hypoplastic left heart syndrome (HLHS) is a severe congenital heart defect (CHD) characterized by hypoplasia of the left ventricle and aorta along with stenosis or atresia of the aortic and mitral valves. HLHS represents only ∼4%-8% of all CHDs but accounts for ∼25% of deaths. HLHS is an isolated defect (i.e., iHLHS) in 70% of families, the vast majority of which are simplex. Despite intense investigation, the genetic basis of iHLHS remains largely unknown. We performed exome sequencing on 331 families with iHLHS aggregated from four independent cohorts. A Mendelian-model-based analysis demonstrated that iHLHS was not due to single, large-effect alleles in genes previously reported to underlie iHLHS or CHD in >90% of families in this cohort. Gene-based association testing identified increased risk for iHLHS associated with variation in CAPN2 (p = 1.8 × 10-5), encoding a protein involved in functional adhesion. Functional validation studies in a vertebrate animal model (Xenopus laevis) confirmed CAPN2 is essential for cardiac ventricle morphogenesis and that in vivo loss of calpain function causes hypoplastic ventricle phenotypes and suggest that human CAPN2707C>T and CAPN21112C>T variants, each found in multiple individuals with iHLHS, are hypomorphic alleles. Collectively, our findings show that iHLHS is typically not a Mendelian condition, demonstrate that CAPN2 variants increase risk of iHLHS, and identify a novel pathway involved in HLHS pathogenesis.

Polygenic prediction across populations is influenced by ancestry, genetic architecture, and methodology

Polygenic risk scores (PRSs) developed from multi-ancestry genome-wide association studies (GWASs), PRSmulti, hold promise for improving PRS accuracy and generalizability across populations. To establish best practices for leveraging the increasing diversity of genomic studies, we investigated how various factors affect the performance of PRSmulti compared with PRSs constructed from single-ancestry GWASs (PRSsingle). Through extensive simulations and empirical analyses, we showed that PRSmulti overall outperformed PRSsingle in understudied populations, except when the understudied population represented a small proportion of the multi-ancestry GWAS. Furthermore, integrating PRSs based on local ancestry-informed GWASs and large-scale, European-based PRSs improved predictive performance in understudied African populations, especially for less polygenic traits with large-effect ancestry-enriched variants. Our work highlights the importance of diversifying genomic studies to achieve equitable PRS performance across ancestral populations and provides guidance for developing PRSs from multiple studies.

Evaluation of population structure inferred by principal component analysis or the admixture model

Principal component analysis (PCA) is commonly used in genetics to infer and visualize population structure and admixture between populations. PCA is often interpreted in a way similar to inferred admixture proportions, where it is assumed that individuals belong to one of several possible populations or are admixed between these populations. We propose a new method to assess the statistical fit of PCA (interpreted as a model spanned by the top principal components) and to show that violations of the PCA assumptions affect the fit. Our method uses the chosen top principal components to predict the genotypes. By assessing the covariance (and the correlation) of the residuals (the differences between observed and predicted genotypes), we are able to detect violation of the model assumptions. Based on simulations and genome-wide human data, we show that our assessment of fit can be used to guide the interpretation of the data and to pinpoint individuals that are not well represented by the chosen principal components. Our method works equally on other similar models, such as the admixture model, where the mean of the data is represented by linear matrix decomposition.