Whole exome sequencing analysis identifies genes for alcohol consumption.

Alcohol consumption is a heritable behavior seriously endangers human health. However, genetic studies on alcohol consumption primarily focuses on common variants, while insights from rare coding variants are lacking. Here we leverage whole exome sequencing data across 304,119 white British individuals from UK Biobank to identify protein-coding variants associated with alcohol consumption. Twenty-five variants are associated with alcohol consumption through single variant analysis and thirteen genes through gene-based analysis, ten of which have not been reported previously. Notably, the two unreported alcohol consumption-related genes GIGYF1 and ANKRD12 show enrichment in brain function-related pathways including glial cell differentiation and are strongly expressed in the cerebellum. Phenome-wide association analyses reveal that alcohol consumption-related genes are associated with brain white matter integrity and risk of digestive and neuropsychiatric diseases. In summary, this study enhances the comprehension of the genetic architecture of alcohol consumption and implies biological mechanisms underlying alcohol-related adverse outcomes.

Exome sequencing identifies novel genetic variants associated with varicose veins.

BACKGROUND: Varicose veins (VV) are one of the common human diseases, but the role of genetics in its development is not fully understood. METHODS: We conducted an exome-wide association study of VV using whole-exome sequencing data from the UK Biobank, and focused on common and rare variants using single-variant association analysis and gene-level collapsing analysis. FINDINGS: A total of 13,823,269 autosomal genetic variants were obtained after quality control. We identified 36 VV-related independent common variants mapping to 34 genes by single-variant analysis and three rare variant genes (PIEZO1, ECE1, FBLN7) by collapsing analysis, and most associations between genes and VV were replicated in FinnGen. PIEZO1 was the closest gene associated with VV (P = 5.05 × 10-31), and it was found to reach exome-wide significance in both single-variant and collapsing analyses. Two novel rare variant genes (ECE1 and METTL21A) associated with VV were identified, of which METTL21A was associated only with females. The pleiotropic effects of VV-related genes suggested that body size, inflammation, and pulmonary function are strongly associated with the development of VV. CONCLUSIONS: Our findings highlight the importance of causal genes for VV and provide new directions for treatment.

Exome sequencing identifies novel genes underlying primary congenital glaucoma in the National Birth Defects Prevention Study.

BACKGROUND: Primary congenital glaucoma (PCG) affects approximately 1 in 10,000 live born infants in the United States (U.S.). PCG has a autosomal recessive inheritance pattern, and variable expressivity and reduced penetrance have been reported. Likely causal variants in the most commonly mutated gene, CYP1B1, are less prevalent in the U.S., suggesting that alternative genes may contribute to the condition. This study utilized exome sequencing to investigate the genetic architecture of PCG in the U.S. and to identify novel genes and variants. METHODS: We studied 37 family trios where infants had PCG and were part of the National Birth Defects Prevention Study (births 1997-2011), a U.S. multicenter study of birth defects. Samples underwent exome sequencing and sequence reads were aligned to the human reference sample (NCBI build 37/hg19). Variant filtration was conducted under de novo and Mendelian inheritance models using GEMINI. RESULTS: Among candidate variants, CYP1B1 was most represented (five trios, 13.5%). Twelve probands (32%) had potentially pathogenic variants in other genes not previously linked to PCG but important in eye development and/or to underlie Mendelian conditions with potential phenotypic overlap (e.g., CRYBB2, RXRA, GLI2). CONCLUSION: Variation in the genes identified in this population-based study may help to further explain the genetics of PCG.

Exome sequencing in undiagnosed congenital myopathy reveals new genes and refines genes-phenotypes correlations.

BACKGROUND: Congenital myopathies are severe genetic diseases with a strong impact on patient autonomy and often on survival. A large number of patients do not have a genetic diagnosis, precluding genetic counseling and appropriate clinical management. Our objective was to find novel pathogenic variants and genes associated with congenital myopathies and to decrease diagnostic odysseys and dead-end. METHODS: To identify pathogenic variants and genes implicated in congenital myopathies, we established and conducted the MYOCAPTURE project from 2009 to 2018 to perform exome sequencing in a large cohort of 310 families partially excluded for the main known genes. RESULTS: Pathogenic variants were identified in 156 families (50%), among which 123 families (40%) had a conclusive diagnosis. Only 44 (36%) of the resolved cases were linked to a known myopathy gene with the corresponding phenotype, while 55 (44%) were linked to pathogenic variants in a known myopathy gene with atypical signs, highlighting that most genetic diagnosis could not be anticipated based on clinical-histological assessments in this cohort. An important phenotypic and genetic heterogeneity was observed for the different genes and for the different congenital myopathy subtypes, respectively. In addition, we identified 14 new myopathy genes not previously associated with muscle diseases (20% of all diagnosed cases) that we previously reported in the literature, revealing novel pathomechanisms and potential therapeutic targets. CONCLUSIONS: Overall, this approach illustrates the importance of massive parallel gene sequencing as a comprehensive tool for establishing a molecular diagnosis for families with congenital myopathies. It also emphasizes the contribution of clinical data, histological findings on muscle biopsies, and the availability of DNA samples from additional family members to the diagnostic success rate. This study facilitated and accelerated the genetic diagnosis of congenital myopathies, improved health care for several patients, and opened novel perspectives for either repurposing of existing molecules or the development of novel treatments.

Clinical exome analysis and targeted gene repair of the c.1354dupT variant in iPSC lines from patients with PROM1-related retinopathies exhibiting diverse phenotypes.

BACKGROUND: Inherited retinal dystrophies (IRD) are one of the main causes of incurable blindness worldwide. IRD are caused by mutations in genes that encode essential proteins for the retina, leading to photoreceptor degeneration and loss of visual function. IRD generates an enormous global financial burden due to the lack of understanding of a significant part of its pathophysiology, molecular diagnosis, and the near absence of non-palliative treatment options. Patient-derived induced pluripotent stem cells (iPSC) for IRD seem to be an excellent option for addressing these questions, serving as exceptional tools for in-depth studies of IRD pathophysiology and testing new therapeutic approaches. METHODS: From a cohort of 8 patients with PROM1-related IRD, we identified 3 patients carrying the same variant (c.1354dupT) but expressing three different IRD phenotypes: Cone and rod dystrophy (CORD), Retinitis pigmentosa (RP), and Stargardt disease type 4 (STGD4). These three target patients, along with one healthy relative from each, underwent comprehensive ophthalmic examinations and their genetic panel study was expanded through clinical exome sequencing (CES). Subsequently, non-integrative patient-derived iPSC were generated and fully characterized. Correction of the c.1354dupT mutation was performed using CRISPR/Cas9, and the genetic restoration of the PROM1 gene was confirmed through flow cytometry and western blotting in the patient-derived iPSC lines. RESULTS: CES revealed that 2 target patients with the c.1354dupT mutation presented monoallelic variants in genes associated with the complement system or photoreceptor differentiation and peroxisome biogenesis disorders, respectively. The pluripotency and functionality of the patient-derived iPSC lines were confirmed, and the correction of the target mutation fully restored the capability of encoding Prominin-1 (CD133) in the genetically repaired patient-derived iPSC lines. CONCLUSIONS: The c.1354dupT mutation in the PROM1 gene is associated to three distinct AR phenotypes of IRD. This pleotropic effect might be related to the influence of monoallelic variants in other genes associated with retinal dystrophies. However, further evidence needs to be provided. Future experiments should include gene-edited patient-derived iPSC due to its potential as disease modelling tools to elucidate this matter in question.

Genetic Evaluation of the Patients with Clinically Diagnosed Inborn Errors of Immunity by Whole Exome Sequencing: Results from a Specialized Research Center for Immunodeficiency in Türkiye.

Molecular diagnosis of inborn errors of immunity (IEI) plays a critical role in determining patients' long-term prognosis, treatment options, and genetic counseling. Over the past decade, the broader utilization of next-generation sequencing (NGS) techniques in both research and clinical settings has facilitated the evaluation of a significant proportion of patients for gene variants associated with IEI. In addition to its role in diagnosing known gene defects, the application of high-throughput techniques such as targeted, exome, and genome sequencing has led to the identification of novel disease-causing genes. However, the results obtained from these different methods can vary depending on disease phenotypes or patient characteristics. In this study, we conducted whole-exome sequencing (WES) in a sizable cohort of IEI patients, consisting of 303 individuals from 21 different clinical immunology centers in Türkiye. Our analysis resulted in likely genetic diagnoses for 41.1% of the patients (122 out of 297), revealing 52 novel variants and uncovering potential new IEI genes in six patients. The significance of understanding outcomes across various IEI cohorts cannot be overstated, and we believe that our findings will make a valuable contribution to the existing literature and foster collaborative research between clinicians and basic science researchers.

The role of distinct APOBEC/ADAR mRNA levels in mutational signatures linked to aging and ultraviolet radiation.

The APOBEC/AID family is known for its mutator activity, and recent evidence also supports the potential impact of ADARs. Furthermore, the mutator impacts of APOBEC/ADAR mutations have not yet been investigated. Assessment of pancancer TCGA exomes identified enriched somatic variants among exomes with nonsynonymous APOBEC1, APOBEC3B, APOBEC3C, ADAR, and ADARB1 mutations, compared to exomes with synonymous ones. Principal component (PC) analysis reduced the number of potential players to eight in cancer exomes/genomes, and to five in cancer types. Multivariate regression analysis was used to assess the impact of the PCs on each COSMIC mutational signature among pancancer exomes/genomes and particular cancers, identifying several novel links, including SBS17b, SBS18, and ID7 mainly determined by APOBEC1 mRNA levels; SBS40, ID1, and ID2 by age; SBS3 and SBS16 by APOBEC3A/APOBEC3B mRNA levels; ID5 and DBS9 by DNA repair/replication (DRR) defects; and SBS7a-d, SBS38, ID4, ID8, ID13, and DBS1 by ultraviolet (UV) radiation/ADARB1 mRNA levels. APOBEC/ADAR mutations appeared to potentiate the impact of DRR defects on several mutational signatures, and some factors seemed to inversely affect certain signatures. These findings potentially implicate certain APOBEC/ADAR mutations/mRNA levels in distinct mutational signatures, particularly APOBEC1 mRNA levels in aging-related signatures and ADARB1 mRNA levels in UV radiation-related signatures.

Systematic rare variant analyses identify RAB32 as a susceptibility gene for familial Parkinson's disease.

Despite substantial progress, causal variants are identified only for a minority of familial Parkinson's disease (PD) cases, leaving high-risk pathogenic variants unidentified1,2. To identify such variants, we uniformly processed exome sequencing data of 2,184 index familial PD cases and 69,775 controls. Exome-wide analyses converged on RAB32 as a novel PD gene identifying c.213C > G/p.S71R as a high-risk variant presenting in ~0.7% of familial PD cases while observed in only 0.004% of controls (odds ratio of 65.5). This variant was confirmed in all cases via Sanger sequencing and segregated with PD in three families. RAB32 encodes a small GTPase known to interact with LRRK2 (refs. 3,4). Functional analyses showed that RAB32 S71R increases LRRK2 kinase activity, as indicated by increased autophosphorylation of LRRK2 S1292. Here our results implicate mutant RAB32 in a key pathological mechanism in PD-LRRK2 kinase activity5-7-and thus provide novel insights into the mechanistic connections between RAB family biology, LRRK2 and PD risk.

Exome sequence analysis identifies rare coding variants associated with a machine learning-based marker for coronary artery disease.

Coronary artery disease (CAD) exists on a spectrum of disease represented by a combination of risk factors and pathogenic processes. An in silico score for CAD built using machine learning and clinical data in electronic health records captures disease progression, severity and underdiagnosis on this spectrum and could enhance genetic discovery efforts for CAD. Here we tested associations of rare and ultrarare coding variants with the in silico score for CAD in the UK Biobank, All of Us Research Program and BioMe Biobank. We identified associations in 17 genes; of these, 14 show at least moderate levels of prior genetic, biological and/or clinical support for CAD. We also observed an excess of ultrarare coding variants in 321 aggregated CAD genes, suggesting more ultrarare variant associations await discovery. These results expand our understanding of the genetic etiology of CAD and illustrate how digital markers can enhance genetic association investigations for complex diseases.

Exome-wide evidence of compound heterozygous effects across common phenotypes in the UK Biobank.

The phenotypic impact of compound heterozygous (CH) variation has not been investigated at the population scale. We phased rare variants (MAF ∼0.001%) in the UK Biobank (UKBB) exome-sequencing data to characterize recessive effects in 175,587 individuals across 311 common diseases. A total of 6.5% of individuals carry putatively damaging CH variants, 90% of which are only identifiable upon phasing rare variants (MAF < 0.38%). We identify six recessive gene-trait associations (p < 1.68 × 10-7) after accounting for relatedness, polygenicity, nearby common variants, and rare variant burden. Of these, just one is discovered when considering homozygosity alone. Using longitudinal health records, we additionally identify and replicate a novel association between bi-allelic variation in ATP2C2 and an earlier age at onset of chronic obstructive pulmonary disease (COPD) (p < 3.58 × 10-8). Genetic phase contributes to disease risk for gene-trait pairs: ATP2C2-COPD (p = 0.000238), FLG-asthma (p = 0.00205), and USH2A-visual impairment (p = 0.0084). We demonstrate the power of phasing large-scale genetic cohorts to discover phenome-wide consequences of compound heterozygosity.

Clinical exome sequencing unravels the diverse spectrum of genetic heterogeneity and genotype-phenotype correlations in hypertrophic cardiomyopathy.

BACKGROUND: Catalogues of pathogenic genetic mutations in hypertrophic cardiomyopathy (HCM) are disproportionately small when compared to that of the size of the population with South Asian ancestry and their collective increased risk of heart disease. METHODS: We conducted clinical exome sequencing of 200 HCM patients to identified cardiomyopathy-associated genetic mutations. The clinical and echocardiographic characteristics of genotype-positive and genotype-negative patients were compared, and the likelihood of detecting a positive genetic test result was evaluated. Allelic burden analysis was done to compare the minor allele frequencies (MAF) of the pathogenic or likely pathogenic (P/LP) variants and variants of uncertain significance (VUSs) identified in the cohort against various population genomics databases. RESULTS: The genetic yield was 40% for P/LP variants, with MYBPC3 and MYH7 as the predominant sarcomere genes. Younger age-at-diagnosis, family history of HCM, asymmetric hypertrophic (ASH) pattern, the ratio of the interventricular septum to posterior wall thickness (IVS/PW ratio), left atrial (LA) dimensions, severe mitral regurgitation grade (MR grade), late gadolinium enhancement (LGE) detected fibrosis and absence of hypertension were associated with an increased likelihood of HCM-associated variants. Patients who experienced ventricular tachycardia and premature cardiovascular death were significantly likely to carry MYBPC3 or loss-of-function variants. LA and interventricular septal (IVS) dimensions were associated with MYH7 variants. The rare variant burden for P/LP variants and VUSs was significantly enriched in HCM cases compared to population controls. CONCLUSION: Our study provides a comprehensive evaluation of HCM-associated genetic mutations from an Indian population. The identified genotype-phenotype associations could improve the yield of targeted genetic testing in HCM.

Genome Sequencing for Diagnosing Rare Diseases.

BACKGROUND: Genetic variants that cause rare disorders may remain elusive even after expansive testing, such as exome sequencing. The diagnostic yield of genome sequencing, particularly after a negative evaluation, remains poorly defined. METHODS: We sequenced and analyzed the genomes of families with diverse phenotypes who were suspected to have a rare monogenic disease and for whom genetic testing had not revealed a diagnosis, as well as the genomes of a replication cohort at an independent clinical center. RESULTS: We sequenced the genomes of 822 families (744 in the initial cohort and 78 in the replication cohort) and made a molecular diagnosis in 218 of 744 families (29.3%). Of the 218 families, 61 (28.0%) - 8.2% of families in the initial cohort - had variants that required genome sequencing for identification, including coding variants, intronic variants, small structural variants, copy-neutral inversions, complex rearrangements, and tandem repeat expansions. Most families in which a molecular diagnosis was made after previous nondiagnostic exome sequencing (63.5%) had variants that could be detected by reanalysis of the exome-sequence data (53.4%) or by additional analytic methods, such as copy-number variant calling, to exome-sequence data (10.8%). We obtained similar results in the replication cohort: in 33% of the families in which a molecular diagnosis was made, or 8% of the cohort, genome sequencing was required, which showed the applicability of these findings to both research and clinical environments. CONCLUSIONS: The diagnostic yield of genome sequencing in a large, diverse research cohort and in a small clinical cohort of persons who had previously undergone genetic testing was approximately 8% and included several types of pathogenic variation that had not previously been detected by means of exome sequencing or other techniques. (Funded by the National Human Genome Research Institute and others.).

Genetic drivers and cellular selection of female mosaic X chromosome loss.

Mosaic loss of the X chromosome (mLOX) is the most common clonal somatic alteration in leukocytes of female individuals1,2, but little is known about its genetic determinants or phenotypic consequences. Here, to address this, we used data from 883,574 female participants across 8 biobanks; 12% of participants exhibited detectable mLOX in approximately 2% of leukocytes. Female participants with mLOX had an increased risk of myeloid and lymphoid leukaemias. Genetic analyses identified 56 common variants associated with mLOX, implicating genes with roles in chromosomal missegregation, cancer predisposition and autoimmune diseases. Exome-sequence analyses identified rare missense variants in FBXO10 that confer a twofold increased risk of mLOX. Only a small fraction of associations was shared with mosaic Y chromosome loss, suggesting that distinct biological processes drive formation and clonal expansion of sex chromosome missegregation. Allelic shift analyses identified X chromosome alleles that are preferentially retained in mLOX, demonstrating variation at many loci under cellular selection. A polygenic score including 44 allelic shift loci correctly inferred the retained X chromosomes in 80.7% of mLOX cases in the top decile. Our results support a model in which germline variants predispose female individuals to acquiring mLOX, with the allelic content of the X chromosome possibly shaping the magnitude of clonal expansion.

A deep catalogue of protein-coding variation in 983,578 individuals.

Rare coding variants that substantially affect function provide insights into the biology of a gene1-3. However, ascertaining the frequency of such variants requires large sample sizes4-8. Here we present a catalogue of human protein-coding variation, derived from exome sequencing of 983,578 individuals across diverse populations. In total, 23% of the Regeneron Genetics Center Million Exome (RGC-ME) data come from individuals of African, East Asian, Indigenous American, Middle Eastern and South Asian ancestry. The catalogue includes more than 10.4 million missense and 1.1 million predicted loss-of-function (pLOF) variants. We identify individuals with rare biallelic pLOF variants in 4,848 genes, 1,751 of which have not been previously reported. From precise quantitative estimates of selection against heterozygous loss of function (LOF), we identify 3,988 LOF-intolerant genes, including 86 that were previously assessed as tolerant and 1,153 that lack established disease annotation. We also define regions of missense depletion at high resolution. Notably, 1,482 genes have regions that are depleted of missense variants despite being tolerant of pLOF variants. Finally, we estimate that 3% of individuals have a clinically actionable genetic variant, and that 11,773 variants reported in ClinVar with unknown significance are likely to be deleterious cryptic splice sites. To facilitate variant interpretation and genetics-informed precision medicine, we make this resource of coding variation from the RGC-ME dataset publicly accessible through a variant allele frequency browser.

Expedited Exome Reanalysis Following Deep Phenotyping and Muscle Biopsy in Suspected Mitochondrial Disorder.

BACKGROUND: Exome sequencing (ES) is a useful tool in diagnosing suspected mitochondrial disease but can miss pathogenic variants for several reasons. Additional testing, such as muscle biopsy or biochemical testing, can be helpful in exome-negative cases. METHODS: We report a patient who presented with repeated episodes of lactic acidosis and failure to thrive. RESULTS: ES and mitochondrial sequencing were initially negative but clinical suspicion for mitochondrial disease remained high. After muscle biopsy showed evidence of mitochondrial dysfunction, the ES was reanalyzed and revealed novel variants in AARS2. CONCLUSION: This case demonstrates the importance of muscle biopsy and biochemical testing in evaluating patients with a high suspicion of mitochondrial disease, even in the genomics era. Closed-loop communication between molecular genetics laboratories and clinical geneticists is an important step to help establish diagnosis in unsolved cases.

Analysis of somatic mutations in whole blood from 200,618 individuals identifies pervasive positive selection and novel drivers of clonal hematopoiesis.

Human aging is marked by the emergence of a tapestry of clonal expansions in dividing tissues, particularly evident in blood as clonal hematopoiesis (CH). CH, linked to cancer risk and aging-related phenotypes, often stems from somatic mutations in a set of established genes. However, the majority of clones lack known drivers. Here we infer gene-level positive selection in whole blood exomes from 200,618 individuals in UK Biobank. We identify 17 additional genes, ZBTB33, ZNF318, ZNF234, SPRED2, SH2B3, SRCAP, SIK3, SRSF1, CHEK2, CCDC115, CCL22, BAX, YLPM1, MYD88, MTA2, MAGEC3 and IGLL5, under positive selection at a population level, and validate this selection pattern in 10,837 whole genomes from single-cell-derived hematopoietic colonies. Clones with mutations in these genes grow in frequency and size with age, comparable to classical CH drivers. They correlate with heightened risk of infection, death and hematological malignancy, highlighting the significance of these additional genes in the aging process.

Tumor Antigens beyond the Human Exome.

With the advent of immunotherapeutics, a new era in the combat against cancer has begun. Particularly promising are neo-epitope-targeted therapies as the expression of neo-antigens is tumor-specific. In turn, this allows the selective targeting and killing of cancer cells whilst healthy cells remain largely unaffected. So far, many advances have been made in the development of treatment options which are tailored to the individual neo-epitope repertoire. The next big step is the achievement of efficacious "off-the-shelf" immunotherapies. For this, shared neo-epitopes propose an optimal target. Given the tremendous potential, a thorough understanding of the underlying mechanisms which lead to the formation of neo-antigens is of fundamental importance. Here, we review the various processes which result in the formation of neo-epitopes. Broadly, the origin of neo-epitopes can be categorized into three groups: canonical, noncanonical, and viral neo-epitopes. For the canonical neo-antigens that arise in direct consequence of somatic mutations, we summarize past and recent findings. Beyond that, our main focus is put on the discussion of noncanonical and viral neo-epitopes as we believe that targeting those provides an encouraging perspective to shape the future of cancer immunotherapeutics.

Characterization of the Common Genetic Variation in the Spanish Population of Navarre.

Large-scale genomic studies have significantly increased our knowledge of genetic variability across populations. Regional genetic profiling is essential for distinguishing common benign variants from disease-causing ones. To this end, we conducted a comprehensive characterization of exonic variants in the population of Navarre (Spain), utilizing whole genome sequencing data from 358 unrelated individuals of Spanish origin. Our analysis revealed 61,410 biallelic single nucleotide variants (SNV) within the Navarrese cohort, with 35% classified as common (MAF > 1%). By comparing allele frequency data from 1000 Genome Project (excluding the Iberian cohort of Spain, IBS), Genome Aggregation Database, and a Spanish cohort (including IBS individuals and data from Medical Genome Project), we identified 1069 SNVs common in Navarre but rare (MAF ≤ 1%) in all other populations. We further corroborated this observation with a second regional cohort of 239 unrelated exomes, which confirmed 676 of the 1069 SNVs as common in Navarre. In conclusion, this study highlights the importance of population-specific characterization of genetic variation to improve allele frequency filtering in sequencing data analysis to identify disease-causing variants.