A harmonized public resource of deeply sequenced diverse human genomes

Genetic variant classification by predicted protein structure: A case study on IRF6

Next-generation genome sequencing has revolutionized genetic testing, identifying numerous rare disease-associated gene variants. However, to impute pathogenicity, computational approaches remain inadequate and functional testing of gene variant is required to provide the highest level of evidence. The emergence of AlphaFold2 has transformed the field of protein structure determination, and here we outline a strategy that leverages predicted protein structure to enhance genetic variant classification. We used the gene IRF6 as a case study due to its clinical relevance, its critical role in cleft lip/palate malformation, and the availability of experimental data on the pathogenicity of IRF6 gene variants through phenotype rescue experiments in irf6-/- zebrafish. We compared results from over 30 pathogenicity prediction tools on 37 IRF6 missense variants. IRF6 lacks an experimentally derived structure, so we used predicted structures to explore associations between mutational clustering and pathogenicity. We found that among these variants, 19 of 37 were unanimously predicted as deleterious by computational tools. Comparing in silico predictions with experimental findings, 12 variants predicted as pathogenic were experimentally determined as benign. Even with the recently published AlphaMissense model, 15/18 (83%) of the predicted pathogenic variants were experimentally determined as benign. In comparison, mapping variants to the protein revealed deleterious mutation clusters around the protein binding domain, whereas N-terminal variants tend to be benign, suggesting the importance of structural information in determining pathogenicity of mutations in this gene. In conclusion, incorporating gene-specific structural features of known pathogenic/benign mutations may provide meaningful insights into pathogenicity predictions in a gene-specific manner and facilitate the interpretation of variant pathogenicity.

Evaluation of genotype imputation using Glimpse tools on low coverage ancient DNA

Ancient DNA provides a unique frame for directly studying human population genetics in time and space. Still, since most of the ancient genomic data is low coverage, analysis is confronted with a low number of SNPs, genotype uncertainties, and reference-bias. Here, we for the first time benchmark the two distinct versions of Glimpse tools on 120 ancient human genomes from Eurasia including those largely from previously under-evaluated regions and compare the performance of genotype imputation with de facto analysis approaches for low coverage genomic data analysis. We further investigate the impact of two distinct reference panels on imputation accuracy for low coverage genomic data. We compute accuracy statistics and perform PCA and f4-statistics to explore the behaviour of genotype imputation on low coverage data regarding (i)two versions of Glimpse, (ii)two reference panels, (iii)four post-imputation filters and coverages, as well as (iv)data type and geographical origin of the samples on the analyses. Our results reveal that even for 0.1X coverage ancient human genomes, genotype imputation using Glimpse-v2 is suitable. Additionally, using the 1000 Genomes merged with Human Genome Diversity Panel improves the accuracy of imputation for the rare variants with low MAF, which might be important not only for ancient genomics but also for modern human genomic studies based on low coverage data and for haplotype-based analysis. Most importantly, we reveal that genotype imputation of low coverage ancient human genomes reduces the genetic affinity of the samples towards human reference genome. Through solving one of the most challenging biases in data analysis, so-called reference bias, genotype imputation using Glimpse v2 is promising for low coverage ancient human genomic data analysis and for rare-variant-based and haplotype-based analysis.

Genetic admixture predictors of fetal alcohol spectrum disorders (FASD) in a South African population

Ancestrally admixed populations are underrepresented in genetic studies of complex diseases, which are still dominated by European-descent populations. This is relevant not only from a representation standpoint but also because of admixed populations' unique features, including being enriched for rare variants, for which effect sizes are disproportionately larger than common polymorphisms. Furthermore, results from these populations may be generalizable to other populations. The South African Cape Coloured (SACC) population is genetically admixed, with one of the highest prevalences of fetal alcohol spectrum disorders (FASD) worldwide. We profiled its admixture and examined associations between ancestry profiles and FASD outcomes using two longitudinal birth cohorts (N=308 mothers, 280 children) designed to examine effects of prenatal alcohol exposure on development. Participants were genotyped via MEGA-ex array to capture common and rare variants. Rare variants were overrepresented in our SACC cohorts, with numerous polymorphisms being monomorphic in other reference populations (e.g., ∼30,000 and ∼ 221,000 variants in gnomAD European and Asian populations, respectively). The cohorts showed global African (51 %; Bantu and San); European (26 %; Northern/Western); South Asian (18 %); and East Asian (5 %; largely Southern regions) ancestries. The cohorts exhibited high rates of homozygosity (6 %), with regions of homozygosity harboring more deleterious variants when lying within African local-ancestry genomic segments. Both maternal and child ancestry profiles were associated with higher FASD risk, and maternal and child ancestry-by-prenatal alcohol exposure interaction effects were seen on child cognition. Our findings indicate that the SACC population may be a valuable asset to identify novel disease-associated genetic loci for FASD and other diseases.

Molecular and clinical characterization of a founder mutation causing G6PC3 deficiency

G6PC3 deficiency is a monogenic immunometabolic disorder that causes syndromic congenital neutropenia. Patients display heterogeneous extra-hematological manifestations, contributing to delayed diagnosis. Here, we investigated the origin and functional consequence of the G6PC3 c.210delC variant found in patients of Mexican origin. Based on the shared haplotypes amongst carriers of the c.210delC mutation, we estimated that this variant originated from a founder effect in a common ancestor. Furthermore, by ancestry analysis, we concluded that it originated in the indigenous Mexican population. At the protein level, we showed that this frameshift mutation leads to an aberrant protein expression in overexpression and patient-derived cells. G6PC3 pathology is driven by the intracellular accumulation of the metabolite 1,5-anhydroglucitol-6-phosphate (1,5-AG6P) that inhibits glycolysis. We characterized how the variant c.210delC impacts glycolysis by performing extracellular flux assays on patient-derived cells. When treated with 1,5-anhydroglucitol (1,5-AG), the precursor to 1,5-AG6P, patient-derived cells exhibited markedly reduced engagement of glycolysis. Finally, we compared the clinical presentation of patients with the mutation c.210delC and all other G6PC3 deficient patients reported in the literature to date, and we found that c.210delC carriers display all prominent clinical features observed in prior G6PC3 deficient patients. In conclusion, G6PC3 c.210delC is a loss-of-function mutation that arose from a founder effect in the indigenous Mexican population. These findings may facilitate the diagnosis of additional patients in this geographical area. Moreover, the in vitro 1,5-AG-dependent functional assay used in our study could be employed to assess the pathogenicity of additional G6PC3 variants.

Biostatistical Aspects of Whole Genome Sequencing Studies: Preprocessing and Quality Control

Rapid advances in high-throughput DNA sequencing technologies have enabled large-scale whole genome sequencing (WGS) studies. Before performing association analysis between phenotypes and genotypes, preprocessing and quality control (QC) of the raw sequence data need to be performed. Because many biostatisticians have not been working with WGS data so far, we first sketch Illumina's short-read sequencing technology. Second, we explain the general preprocessing pipeline for WGS studies. Third, we provide an overview of important QC metrics, which are applied to WGS data: on the raw data, after mapping and alignment, after variant calling, and after multisample variant calling. Fourth, we illustrate the QC with the data from the GENEtic SequencIng Study Hamburg-Davos (GENESIS-HD), a study involving more than 9000 human whole genomes. All samples were sequenced on an Illumina NovaSeq 6000 with an average coverage of 35× using a PCR-free protocol. For QC, one genome in a bottle (GIAB) trio was sequenced in four replicates, and one GIAB sample was successfully sequenced 70 times in different runs. Fifth, we provide empirical data on the compression of raw data using the DRAGEN original read archive (ORA). The most important quality metrics in the application were genetic similarity, sample cross-contamination, deviations from the expected Het/Hom ratio, relatedness, and coverage. The compression ratio of the raw files using DRAGEN ORA was 5.6:1, and compression time was linear by genome coverage. In summary, the preprocessing, joint calling, and QC of large WGS studies are feasible within a reasonable time, and efficient QC procedures are readily available.

Germline variation contributes to false negatives in CRISPR-based experiments with varying burden across ancestries

Reducing disparities is vital for equitable access to precision treatments in cancer. Socioenvironmental factors are a major driver of disparities, but differences in genetic variation likely also contribute. The impact of genetic ancestry on prioritization of cancer targets in drug discovery pipelines has not been systematically explored due to the absence of pre-clinical data at the appropriate scale. Here, we analyze data from 611 genome-scale CRISPR/Cas9 viability experiments in human cell line models to identify ancestry-associated genetic dependencies essential for cell survival. Surprisingly, we find that most putative associations between ancestry and dependency arise from artifacts related to germline variants. Our analysis suggests that for 1.2-2.5% of guides, germline variants in sgRNA targeting sequences reduce cutting by the CRISPR/Cas9 nuclease, disproportionately affecting cell models derived from individuals of recent African descent. We propose three approaches to mitigate this experimental bias, enabling the scientific community to address these disparities.

Individuals with JAK1 variants are affected by syndromic features encompassing autoimmunity, atopy, colitis, and dermatitis

Inborn errors of immunity lead to autoimmunity, inflammation, allergy, infection, and/or malignancy. Disease-causing JAK1 gain-of-function (GoF) mutations are considered exceedingly rare and have been identified in only four families. Here, we use forward and reverse genetics to identify 59 individuals harboring one of four heterozygous JAK1 variants. In vitro and ex vivo analysis of these variants revealed hyperactive baseline and cytokine-induced STAT phosphorylation and interferon-stimulated gene (ISG) levels compared with wild-type JAK1. A systematic review of electronic health records from the BioME Biobank revealed increased likelihood of clinical presentation with autoimmunity, atopy, colitis, and/or dermatitis in JAK1 variant-positive individuals. Finally, treatment of one affected patient with severe atopic dermatitis using the JAK1/JAK2-selective inhibitor, baricitinib, resulted in clinically significant improvement. These findings suggest that individually rare JAK1 GoF variants may underlie an emerging syndrome with more common presentations of autoimmune and inflammatory disease (JAACD syndrome). More broadly, individuals who present with such conditions may benefit from genetic testing for the presence of JAK1 GoF variants.

The Polygenic Score Catalog: new functionality and tools to enable FAIR research

Polygenic scores (PGS) have transformed human genetic research and have multiple potential clinical applications, including risk stratification for disease prevention and prediction of treatment response. Here, we present a series of recent enhancements to the PGS Catalog (www.PGSCatalog.org), the largest findable, accessible, interoperable, and reusable (FAIR) repository of PGS. These include expansions in data content and ancestral diversity as well as the addition of new features. We further present the PGS Catalog Calculator (pgsc_calc, https://github.com/PGScatalog/pgsc_calc), an open-source, scalable and portable pipeline to reproducibly calculate PGS that securely democratizes equitable PGS applications by implementing genetic ancestry estimation and score normalization using reference data. With the PGS Catalog & calculator users can now quantify an individual's genetic predisposition for hundreds of common diseases and clinically relevant traits. Taken together, these updates and tools facilitate the next generation of PGS, thus lowering barriers to the clinical studies necessary to identify where PGS may be integrated into clinical practice.

Molecular and clinical characterization of a founder mutation causing G6PC3 deficiency

Background

G6PC3 deficiency is a rare genetic disorder that causes syndromic congenital neutropenia. It is driven by the intracellular accumulation of a metabolite named 1,5-anhydroglucitol-6-phosphate (1,5-AG6P) that inhibits glycolysis. Patients display heterogeneous extra-hematological manifestations, contributing to delayed diagnosis.

Objective

The G6PC3 c.210delC variant has been identified in patients of Mexican origin. We set out to study the origin and functional consequence of this mutation. Furthermore, we sought to characterize the clinical phenotypes caused by it.

Methods

Using whole-genome sequencing data, we conducted haplotype analysis to estimate the age of this allele and traced its ancestral origin. We examined how this mutation affected G6PC3 protein expression and performed extracellular flux assays on patient-derived cells to characterize how this mutation impacts glycolysis. Finally, we compared the clinical presentations of patients with the c.210delC mutation relative to other G6PC3 deficient patients published to date.

Results

Based on the length of haplotypes shared amongst ten carriers of the G6PC3 c.210delC mutation, we estimated that this variant originated in a common ancestor of indigenous American origin. The mutation causes a frameshift that introduces a premature stop codon, leading to a complete loss of G6PC3 protein expression. When treated with 1,5-anhydroglucitol (1,5-AG), the precursor to 1,5-AG6P, patient-derived cells exhibited markedly reduced engagement of glycolysis. Clinically, c.210delC carriers display all the clinical features of syndromic severe congenital neutropenia type 4 observed in prior reports of G6PC3 deficiency.

Conclusion

The G6PC3 c.210delC is a loss-of-function mutation that arose from a founder effect in the indigenous Mexican population. These findings may facilitate the diagnosis of additional patients in this geographical area. Moreover, the in vitro 1,5-AG-dependent functional assay used in our study could be employed to assess the pathogenicity of additional G6PC3 variants.

A multi-ancestry GWAS of Fuchs corneal dystrophy highlights the contributions of laminins, collagen, and endothelial cell regulation

Fuchs endothelial corneal dystrophy (FECD) is a leading indication for corneal transplantation, but its molecular etiology remains poorly understood. We performed genome-wide association studies (GWAS) of FECD in the Million Veteran Program followed by multi-ancestry meta-analysis with the previous largest FECD GWAS, for a total of 3970 cases and 333,794 controls. We confirm the previous four loci, and identify eight novel loci: SSBP3, THSD7A, LAMB1, PIDD1, RORA, HS3ST3B1, LAMA5, and COL18A1. We further confirm the TCF4 locus in GWAS for admixed African and Hispanic/Latino ancestries and show an enrichment of European-ancestry haplotypes at TCF4 in FECD cases. Among the novel associations are low frequency missense variants in laminin genes LAMA5 and LAMB1 which, together with previously reported LAMC1, form laminin-511 (LM511). AlphaFold 2 protein modeling, validated through homology, suggests that mutations at LAMA5 and LAMB1 may destabilize LM511 by altering inter-domain interactions or extracellular matrix binding. Finally, phenome-wide association scans and colocalization analyses suggest that the TCF4 CTG18.1 trinucleotide repeat expansion leads to dysregulation of ion transport in the corneal endothelium and has pleiotropic effects on renal function.

Genetic admixture predictors of fetal alcohol spectrum disorders (FASD) in the South African Cape Coloured population

Ancestrally admixed populations are underrepresented in genetic studies of complex diseases, which are still dominated by European-descent populations. This is relevant not only from a representation standpoint but also because of admixed populations' unique features, including being enriched for rare variants, for which effect sizes are disproportionately larger than common polymorphisms. Furthermore, results from these populations may be generalizable to other populations. The South African Cape Coloured (SACC) population is genetically admixed, with one of the highest prevalences of fetal alcohol spectrum disorders (FASD) worldwide. We profiled its admixture and examined associations between ancestry profiles and FASD outcomes using two longitudinal birth cohorts ( N =308 mothers, 280 children) designed to examine effects of prenatal alcohol exposure on development. Participants were genotyped via MEGA-ex array to capture common and rare variants. Rare variants were overrepresented in our SACC cohorts, with numerous polymorphisms being monomorphic in other reference populations (e.g., ∼30,000 and ∼221,000 variants in gnomAD European and Asian populations, respectively). The cohorts showed global African (51%; Bantu and San); European (26%; Northern/Western); South Asian (18%); and East Asian (5%; largely Southern regions) ancestries. The cohorts exhibited high rates of homozygosity (6%), with regions of homozygosity harboring more deleterious variants when lying within African local-ancestry genomic segments. Both maternal and child ancestry profiles were associated with FASD risk and altered severity of prenatal alcohol exposure-related cognitive deficits in the child. Our findings indicate that the SACC population may be a valuable asset to identify novel disease-associated genetic loci for FASD and other diseases.

Evaluating the Efficacy of Type 2 Diabetes Polygenic Risk Scores in an Independent European Population

Numerous type 2 diabetes (T2D) polygenic risk scores (PGSs) have been developed to predict individuals' predisposition to the disease. An independent assessment and verification of the best-performing PGS are warranted to allow for a rapid application of developed models. To date, only 3% of T2D PGSs have been evaluated. In this study, we assessed all (n = 102) presently published T2D PGSs in an independent cohort of 3718 individuals, which has not been included in the construction or fine-tuning of any T2D PGS so far. We further chose the best-performing PGS, assessed its performance across major population principal component analysis (PCA) clusters, and compared it with newly developed population-specific T2D PGS. Our findings revealed that 88% of the published PGSs were significantly associated with T2D; however, their performance was lower than what had been previously reported. We found a positive association of PGS improvement over the years (p-value = 8.01 × 10-4 with PGS002771 currently showing the best discriminatory power (area under the receiver operating characteristic (AUROC) = 0.669) and PGS003443 exhibiting the strongest association PGS003443 (odds ratio (OR) = 1.899). Further investigation revealed no difference in PGS performance across major population PCA clusters and when compared with newly developed population-specific PGS. Our findings revealed a positive trend in T2D PGS performance, consistently identifying high-T2D-risk individuals in an independent European population.

The Scalable Variant Call Representation: Enabling Genetic Analysis Beyond One Million Genomes

The Variant Call Format (VCF) is widely used in genome sequencing but scales poorly. For instance, we estimate a 150,000 genome VCF would occupy 900 TiB, making it both costly and complicated to produce and analyze. The issue stems from VCF's requirement to densely represent both reference-genotypes and allele-indexed arrays. These requirements lead to unnecessary data duplication and, ultimately, very large files. To address these challenges, we introduce the Scalable Variant Call Representation (SVCR). This representation reduces file sizes by ensuring they scale linearly with samples. SVCR achieves this by adopting reference blocks from the Genomic Variant Call Format (GVCF) and employing local allele indices. SVCR is also lossless and mergeable, allowing for N+1 and N+K incremental joint-calling. We present two implementations of SVCR: SVCR-VCF, which encodes SVCR in VCF format, and VDS, which uses Hail's native format. Our experiments confirm the linear scalability of SVCR-VCF and VDS, in contrast to the super-linear growth seen with standard VCF files. We also discuss the VDS Combiner, a scalable, open-source tool for producing a VDS from GVCFs and unique features of VDS which enable rapid data analysis. SVCR, and VDS in particular, ensure the scientific community can generate, analyze, and disseminate genetics datasets with millions of samples.

A genomic mutational constraint map using variation in 76,156 human genomes

The depletion of disruptive variation caused by purifying natural selection (constraint) has been widely used to investigate protein-coding genes underlying human disorders1-4, but attempts to assess constraint for non-protein-coding regions have proved more difficult. Here we aggregate, process and release a dataset of 76,156 human genomes from the Genome Aggregation Database (gnomAD)-the largest public open-access human genome allele frequency reference dataset-and use it to build a genomic constraint map for the whole genome (genomic non-coding constraint of haploinsufficient variation (Gnocchi)). We present a refined mutational model that incorporates local sequence context and regional genomic features to detect depletions of variation. As expected, the average constraint for protein-coding sequences is stronger than that for non-coding regions. Within the non-coding genome, constrained regions are enriched for known regulatory elements and variants that are implicated in complex human diseases and traits, facilitating the triangulation of biological annotation, disease association and natural selection to non-coding DNA analysis. More constrained regulatory elements tend to regulate more constrained protein-coding genes, which in turn suggests that non-coding constraint can aid the identification of constrained genes that are as yet unrecognized by current gene constraint metrics. We demonstrate that this genome-wide constraint map improves the identification and interpretation of functional human genetic variation.

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.

CHARR efficiently estimates contamination from DNA sequencing data

DNA sample contamination is a major issue in clinical and research applications of whole genome and exome sequencing. Even modest levels of contamination can substantially affect the overall quality of variant calls and lead to widespread genotyping errors. Currently, popular tools for estimating the contamination level use short-read data (BAM/CRAM files), which are expensive to store and manipulate and often not retained or shared widely. We propose a new metric to estimate DNA sample contamination from variant-level whole genome and exome sequence data, CHARR, Contamination from Homozygous Alternate Reference Reads, which leverages the infiltration of reference reads within homozygous alternate variant calls. CHARR uses a small proportion of variant-level genotype information and thus can be computed from single-sample gVCFs or callsets in VCF or BCF formats, as well as efficiently stored variant calls in Hail VDS format. Our results demonstrate that CHARR accurately recapitulates results from existing tools with substantially reduced costs, improving the accuracy and efficiency of downstream analyses of ultra-large whole genome and exome sequencing datasets.