Analysis of 589,306 genomes identifies individuals resilient to severe Mendelian childhood diseases

A PLURIPOTENT STEM CELL PLATFORM FOR IN VITRO SYSTEMS GENETICS STUDIES OF MOUSE DEVELOPMENT

The directed differentiation of pluripotent stem cells (PSCs) from panels of genetically diverse individuals is emerging as a powerful experimental system for characterizing the impact of natural genetic variation on developing cell types and tissues. Here, we establish new PSC lines and experimental approaches for modeling embryonic development in a genetically diverse, outbred mouse stock (Diversity Outbred mice). We show that a range of inbred and outbred PSC lines can be stably maintained in the primed pluripotent state (epiblast stem cells -- EpiSCs) and establish the contribution of genetic variation to phenotypic differences in gene regulation and directed differentiation. Using pooled in vitro fertilization, we generate and characterize a genetic reference panel of Diversity Outbred PSCs (n = 230). Finally, we demonstrate the feasibility of pooled culture of Diversity Outbred EpiSCs as "cell villages", which can facilitate the differentiation of large numbers of EpiSC lines for forward genetic screens. These data can complement and inform similar efforts within the stem cell biology and human genetics communities to model the impact of natural genetic variation on phenotypic variation and disease-risk.

Matrix metalloproteinase-9 deficiency confers resilience in fibrodysplasia ossificans progressiva in a man and mice

Single case studies of extraordinary disease resilience may provide therapeutic insight into conditions for which no definitive treatments exist. An otherwise healthy 35-year-old man (patient-R) with the canonical pathogenic ACVR1R206H variant and the classic congenital great toe malformation of fibrodysplasia ossificans progressiva (FOP) had extreme paucity of post-natal heterotopic ossification (HO) and nearly normal mobility. We hypothesized that patient-R lacked a sufficient post-natal inflammatory trigger for HO. A plasma biomarker survey revealed a reduction in total matrix metalloproteinase-9 (MMP-9) compared to healthy controls and individuals with quiescent FOP. Whole exome sequencing identified compound heterozygous variants in MMP-9 (c.59C > T, p.A20V and c.493G > A, p.D165N). Structural analysis of the D165N variant predicted both decreased MMP-9 secretion and activity that were confirmed by enzyme-linked immunosorbent assay and gelatin zymography. Further, human proinflammatory M1-like macrophages expressing either MMP-9 variant produced significantly less Activin A, an obligate ligand for HO in FOP, compared to wildtype controls. Importantly, MMP-9 inhibition by genetic, biologic, or pharmacologic means in multiple FOP mouse models abrogated trauma-induced HO, sequestered Activin A in the extracellular matrix (ECM), and induced regeneration of injured skeletal muscle. Our data suggest that MMP-9 is a druggable node linking inflammation to HO, orchestrates an existential role in the pathogenesis of FOP, and illustrates that a single patient's clinical phenotype can reveal critical molecular mechanisms of disease that unveil novel treatment strategies.

Sterol O-Acyltransferase 1 (SOAT1): A Genetic Modifier of Niemann-Pick Disease, Type C1

Niemann-Pick disease type C1 (NPC1) is a lysosomal disorder due to impaired intracellular cholesterol transport out of the endolysosomal compartment.. Marked heterogeneity has been observed in individuals with the same NPC1 genotype, thus suggesting a significant effect of modifier genes. Prior work demonstrated that decreased SOAT1 activity decreased disease severity in an NPC1 mouse model. Thus, we hypothesized that a polymorphism associated with decreased SOAT1 expression might influence the NPC1 phenotype. Phenotyping and genomic sequencing of 117 individuals with NPC1 was performed as part of a Natural History trial. Phenotyping included determination of disease severity and disease burden. Significant clinical heterogeneity is present in individuals homozygous for the NPC1I1061T variant and in siblings. Analysis of the SOAT1 polymorphism, rs1044925 (A>C), showed a significant association of the C-allele with earlier age of neurological onset. The C-allele may be associated with a higher Annualized Severity Index Score as well as increased frequency of liver disease and seizures. A polymorphism associated with decreased expression of SOAT1 appears to be a genetic modifier of the NPC1 phenotype. This finding is consistent with prior data showing decreased phenotypic severity in Npc1-/-:Soat1-/- mice and supports efforts to investigate the potential of SOAT1 inhibitors as a potential therapy for NPC1.

Pleiotropy, epistasis and the genetic architecture of quantitative traits

Pleiotropy (whereby one genetic polymorphism affects multiple traits) and epistasis (whereby non-linear interactions between genetic polymorphisms affect the same trait) are fundamental aspects of the genetic architecture of quantitative traits. Recent advances in the ability to characterize the effects of polymorphic variants on molecular and organismal phenotypes in human and model organism populations have revealed the prevalence of pleiotropy and unexpected shared molecular genetic bases among quantitative traits, including diseases. By contrast, epistasis is common between polymorphic loci associated with quantitative traits in model organisms, such that alleles at one locus have different effects in different genetic backgrounds, but is rarely observed for human quantitative traits and common diseases. Here, we review the concepts and recent inferences about pleiotropy and epistasis, and discuss factors that contribute to similarities and differences between the genetic architecture of quantitative traits in model organisms and humans.

PICOTEES: a privacy-preserving online service of phenotype exploration for genetic-diagnostic variants from Chinese children cohorts

The growth in biomedical data resources has raised potential privacy concerns and risks of genetic information leakage. For instance, exome sequencing aids clinical decisions by comparing data through web services, but it requires significant trust between users and providers. To alleviate privacy concerns, the most commonly used strategy is to anonymize sensitive data. Unfortunately, studies have shown that anonymization is insufficient to protect against reidentification attacks. Recently, privacy-preserving technologies have been applied to preserve application utility while protecting the privacy of biomedical data. We present the PICOTEES framework, a privacy-preserving online service of phenotype exploration for genetic-diagnostic variants (https://birthdefectlab.cn:3000/). PICOTEES enables privacy-preserving queries of the phenotype spectrum for a single variant by utilizing trusted execution environment technology, which can protect the privacy of the user's query information, backend models, and data, as well as the final results. We demonstrate the utility and performance of PICOTEES by exploring a bioinformatics dataset. The dataset is from a cohort containing 20,909 genetic testing patients with 3,152,508 variants from the Children's Hospital of Fudan University in China, dominated by the Chinese Han population (>99.9%). Our query results yield a large number of unreported diagnostic variants and previously reported pathogenicity.

Rare copy-number variants as modulators of common disease susceptibility

BACKGROUND

Copy-number variations (CNVs) have been associated with rare and debilitating genomic disorders (GDs) but their impact on health later in life in the general population remains poorly described.

METHODS

Assessing four modes of CNV action, we performed genome-wide association scans (GWASs) between the copy-number of CNV-proxy probes and 60 curated ICD-10 based clinical diagnoses in 331,522 unrelated white British UK Biobank (UKBB) participants with replication in the Estonian Biobank.

RESULTS

We identified 73 signals involving 40 diseases, all of which indicating that CNVs increased disease risk and caused earlier onset. We estimated that 16% of these associations are indirect, acting by increasing body mass index (BMI). Signals mapped to 45 unique, non-overlapping regions, nine of which being linked to known GDs. Number and identity of genes affected by CNVs modulated their pathogenicity, with many associations being supported by colocalization with both common and rare single-nucleotide variant association signals. Dissection of association signals provided insights into the epidemiology of known gene-disease pairs (e.g., deletions in BRCA1 and LDLR increased risk for ovarian cancer and ischemic heart disease, respectively), clarified dosage mechanisms of action (e.g., both increased and decreased dosage of 17q12 impacted renal health), and identified putative causal genes (e.g., ABCC6 for kidney stones). Characterization of the pleiotropic pathological consequences of recurrent CNVs at 15q13, 16p13.11, 16p12.2, and 22q11.2 in adulthood indicated variable expressivity of these regions and the involvement of multiple genes. Finally, we show that while the total burden of rare CNVs-and especially deletions-strongly associated with disease risk, it only accounted for ~ 0.02% of the UKBB disease burden. These associations are mainly driven by CNVs at known GD CNV regions, whose pleiotropic effect on common diseases was broader than anticipated by our CNV-GWAS.

CONCLUSIONS

Our results shed light on the prominent role of rare CNVs in determining common disease susceptibility within the general population and provide actionable insights for anticipating later-onset comorbidities in carriers of recurrent CNVs.

Functional analysis of germline VANGL2 variants using rescue assays of vangl2 knockout zebrafish

Developmental studies have shown that the evolutionarily conserved Wnt Planar Cell Polarity (PCP) pathway is essential for the development of a diverse range of tissues and organs including the brain, spinal cord, heart and sensory organs, as well as establishment of the left-right body axis. Germline mutations in the highly conserved PCP gene VANGL2 in humans have only been associated with central nervous system malformations, and functional testing to understand variant impact has not been performed. Here we report three new families with missense variants in VANGL2 associated with heterotaxy and congenital heart disease p.(Arg169His), non-syndromic hearing loss p.(Glu465Ala) and congenital heart disease with brain defects p.(Arg135Trp). To test the in vivo impact of these and previously described variants, we have established clinically-relevant assays using mRNA rescue of the vangl2 mutant zebrafish. We show that all variants disrupt Vangl2 function, although to different extents and depending on the developmental process. We also begin to identify that different VANGL2 missense variants may be haploinsufficient and discuss evidence in support of pathogenicity. Together, this study demonstrates that zebrafish present a suitable pipeline to investigate variants of unknown significance and suggests new avenues for investigation of the different developmental contexts of VANGL2 function that are clinically meaningful.

Meta-analysis of dispensable essential genes and their interactions with bypass suppressors

Genes have been historically classified as essential or non-essential based on their requirement for viability. However, genomic mutations can sometimes bypass the requirement for an essential gene, challenging the binary classification of gene essentiality. Such dispensable essential genes represent a valuable model for understanding the incomplete penetrance of loss-of-function mutations often observed in natural populations. Here, we compiled data from multiple studies on essential gene dispensability in Saccharomyces cerevisiae to comprehensively characterize these genes. In analyses spanning different evolutionary timescales, dispensable essential genes exhibited distinct phylogenetic properties compared with other essential and non-essential genes. Integration of interactions with suppressor genes that can bypass the gene essentiality revealed the high functional modularity of the bypass suppression network. Furthermore, dispensable essential and bypass suppressor gene pairs reflected simultaneous changes in the mutational landscape of S. cerevisiae strains. Importantly, species in which dispensable essential genes were non-essential tended to carry bypass suppressor mutations in their genomes. Overall, our study offers a comprehensive view of dispensable essential genes and illustrates how their interactions with bypass suppressors reflect evolutionary outcomes.

Chargaff's second parity rule lies at the origin of additive genetic interactions in quantitative traits to make omnigenic selection possible

Background

Francis Crick's central dogma provides a residue-by-residue mechanistic explanation of the flow of genetic information in living systems. However, this principle may not be sufficient for explaining how random mutations cause continuous variation of quantitative highly polygenic complex traits. Chargaff's second parity rule (CSPR), also referred to as intrastrand DNA symmetry, defined as near-exact equalities G ≈ C and A ≈ T within a single DNA strand, is a statistical property of cellular genomes. The phenomenon of intrastrand DNA symmetry was discovered more than 50 years ago; at present, it remains unclear what its biological role is, what the mechanisms are that force cellular genomes to comply strictly with CSPR, and why genomes of certain noncellular organisms have broken intrastrand DNA symmetry. The present work is aimed at studying a possible link between intrastrand DNA symmetry and the origin of genetic interactions in quantitative traits.

Methods

Computational analysis of single-nucleotide polymorphisms in human and mouse populations and of nucleotide composition biases at different codon positions in bacterial and human proteomes.

Results

The analysis of mutation spectra inferred from single-nucleotide polymorphisms observed in murine and human populations revealed near-exact equalities of numbers of reverse complementary mutations, indicating that random genetic variations obey CSPR. Furthermore, nucleotide compositions of coding sequences proved to be statistically interwoven via CSPR because pyrimidine bias at the 3rd codon position compensates purine bias at the 1st and 2nd positions.

Conclusions

According to Fisher's infinitesimal model, we propose that accumulation of reverse complementary mutations results in a continuous phenotypic variation due to small additive effects of statistically interwoven genetic variations. Therefore, additive genetic interactions can be inferred as a statistical entanglement of nucleotide compositions of separate genetic loci. CSPR challenges the neutral theory of molecular evolution-because all random mutations participate in variation of a trait-and provides an alternative solution to Haldane's dilemma by making a gene function diffuse. We propose that CSPR is symmetry of Fisher's infinitesimal model and that genetic information can be transferred in an implicit contactless manner.

Model Organism Modifier (MOM): a user-friendly Galaxy workflow to detect modifiers from genome sequencing data using Caenorhabditis elegans

Genetic modifiers are variants modulating phenotypic outcomes of a primary detrimental variant. They contribute to rare diseases phenotypic variability, but their identification is challenging. Genetic screening with model organisms is a widely used method for demystifying genetic modifiers. Forward genetics screening followed by whole genome sequencing allows the detection of variants throughout the genome but typically produces thousands of candidate variants making the interpretation and prioritization process very time-consuming and tedious. Despite whole genome sequencing is more time and cost-efficient, usage of computational pipelines specific to modifier identification remains a challenge for biological-experiment-focused laboratories doing research with model organisms. To facilitate a broader implementation of whole genome sequencing in genetic screens, we have developed Model Organism Modifier or MOM, a pipeline as a user-friendly Galaxy workflow. Model Organism Modifier analyses raw short-read whole genome sequencing data and implements tailored filtering to provide a Candidate Variant List short enough to be further manually curated. We provide a detailed tutorial to run the Galaxy workflow Model Organism Modifier and guidelines to manually curate the Candidate Variant Lists. We have tested Model Organism Modifier on published and validated Caenorhabditis elegans modifiers screening datasets. As whole genome sequencing facilitates high-throughput identification of genetic modifiers in model organisms, Model Organism Modifier provides a user-friendly solution to implement the bioinformatics analysis of the short-read datasets in laboratories without expertise or support in Bioinformatics.

Dystonin modifiers of junctional epidermolysis bullosa and models of epidermolysis bullosa simplex without dystonia musculorum

The Lamc2jeb junctional epidermolysis bullosa (EB) mouse model has been used to demonstrate that significant genetic modification of EB symptoms is possible, identifying as modifiers Col17a1 and six other quantitative trait loci, several with strong candidate genes including dystonin (Dst/Bpag1). Here, CRISPR/Cas9 was used to alter exon 23 in mouse skin specific isoform Dst-e (Ensembl GRCm38 transcript name Dst-213, transcript ID ENSMUST00000183302.5, protein size 2639AA) and validate a proposed arginine/glutamine difference at amino acid p1226 in B6 versus 129 mice as a modifier of EB. Frame shift deletions (FSD) in mouse Dst-e exon 23 (Dst-eFSD/FSD) were also identified that cause mice carrying wild-type Lamc2 to develop a phenotype similar to human EB simplex without dystonia musculorum. When combined, Dst-eFSD/FSD modifies Lamc2jeb/jeb (FSD+jeb) induced disease in unexpected ways implicating an altered balance between DST-e (BPAG1e) and a rarely reported rodless DST-eS (BPAG1eS) in epithelium as a possible mechanism. Further, FSD+jeb mice with pinnae removed are found to provide a test bed for studying internal epithelium EB disease and treatment without severe skin disease as a limiting factor while also revealing and accelerating significant nasopharynx symptoms present but not previously noted in Lamc2jeb/jeb mice.

Stochastic gene expression and environmental stressors trigger variable somite segmentation phenotypes

Mutations of several genes cause incomplete penetrance and variable expressivity of phenotypes, which are usually attributed to modifier genes or gene-environment interactions. Here, we show stochastic gene expression underlies the variability of somite segmentation defects in embryos mutant for segmentation clock genes her1 or her7. Phenotypic strength is further augmented by low temperature and hypoxia. By performing live imaging of the segmentation clock reporters, we further show that groups of cells with higher oscillation amplitudes successfully form somites while those with lower amplitudes fail to do so. In unfavorable environments, the number of cycles with high amplitude oscillations and the number of successful segmentations proportionally decrease. These results suggest that individual oscillation cycles stochastically fail to pass a threshold amplitude, resulting in segmentation defects in mutants. Our quantitative methodology is adaptable to investigate variable phenotypes of mutant genes in different tissues.

Systematic evaluation of genome sequencing for the diagnostic assessment of autism spectrum disorder and fetal structural anomalies

Short-read genome sequencing (GS) holds the promise of becoming the primary diagnostic approach for the assessment of autism spectrum disorder (ASD) and fetal structural anomalies (FSAs). However, few studies have comprehensively evaluated its performance against current standard-of-care diagnostic tests: karyotype, chromosomal microarray (CMA), and exome sequencing (ES). To assess the clinical utility of GS, we compared its diagnostic yield against these three tests in 1,612 quartet families including an individual with ASD and in 295 prenatal families. Our GS analytic framework identified a diagnostic variant in 7.8% of ASD probands, almost 2-fold more than CMA (4.3%) and 3-fold more than ES (2.7%). However, when we systematically captured copy-number variants (CNVs) from the exome data, the diagnostic yield of ES (7.4%) was brought much closer to, but did not surpass, GS. Similarly, we estimated that GS could achieve an overall diagnostic yield of 46.1% in unselected FSAs, representing a 17.2% increased yield over karyotype, 14.1% over CMA, and 4.1% over ES with CNV calling or 36.1% increase without CNV discovery. Overall, GS provided an added diagnostic yield of 0.4% and 0.8% beyond the combination of all three standard-of-care tests in ASD and FSAs, respectively. This corresponded to nine GS unique diagnostic variants, including sequence variants in exons not captured by ES, structural variants (SVs) inaccessible to existing standard-of-care tests, and SVs where the resolution of GS changed variant classification. Overall, this large-scale evaluation demonstrated that GS significantly outperforms each individual standard-of-care test while also outperforming the combination of all three tests, thus warranting consideration as the first-tier diagnostic approach for the assessment of ASD and FSAs.

Diagnostic implications of pitfalls in causal variant identification based on 4577 molecularly characterized families

Despite large sequencing and data sharing efforts, previously characterized pathogenic variants only account for a fraction of Mendelian disease patients, which highlights the need for accurate identification and interpretation of novel variants. In a large Mendelian cohort of 4577 molecularly characterized families, numerous scenarios in which variant identification and interpretation can be challenging are encountered. We describe categories of challenges that cover the phenotype (e.g. novel allelic disorders), pedigree structure (e.g. imprinting disorders masquerading as autosomal recessive phenotypes), positional mapping (e.g. double recombination events abrogating candidate autozygous intervals), gene (e.g. novel gene-disease assertion) and variant (e.g. complex compound inheritance). Overall, we estimate a probability of 34.3% for encountering at least one of these challenges. Importantly, our data show that by only addressing non-sequencing-based challenges, around 71% increase in the diagnostic yield can be expected. Indeed, by applying these lessons to a cohort of 314 cases with negative clinical exome or genome reports, we could identify the likely causal variant in 54.5%. Our work highlights the need to have a thorough approach to undiagnosed diseases by considering a wide range of challenges rather than a narrow focus on sequencing technologies. It is hoped that by sharing this experience, the yield of undiagnosed disease programs globally can be improved.

Genetic control of mRNA splicing as a potential mechanism for incomplete penetrance of rare coding variants

Exonic variants present some of the strongest links between genotype and phenotype. However, these variants can have significant inter-individual pathogenicity differences, known as variable penetrance. In this study, we propose a model where genetically controlled mRNA splicing modulates the pathogenicity of exonic variants. By first cataloging exonic inclusion from RNA-sequencing data in GTEx V8, we find that pathogenic alleles are depleted on highly included exons. Using a large-scale phased whole genome sequencing data from the TOPMed consortium, we observe that this effect may be driven by common splice-regulatory genetic variants, and that natural selection acts on haplotype configurations that reduce the transcript inclusion of putatively pathogenic variants, especially when limiting to haploinsufficient genes. Finally, we test if this effect may be relevant for autism risk using families from the Simons Simplex Collection, but find that splicing of pathogenic alleles has a penetrance reducing effect here as well. Overall, our results indicate that common splice-regulatory variants may play a role in reducing the damaging effects of rare exonic variants.

Incomplete penetrance and variable expressivity in monogenic diabetes; a challenge but also an opportunity

Monogenic Forms of Diabetes (MFD) account for about 3% of all diabetes, and their accurate diagnosis often results in life-changing therapeutic reassignment for the patients. Like other Mendelian diseases, reduced penetrance and variable expressivity are often seen in several different types of MFD, where symptoms develop only in a portion of the persons who carry the pathogenic variant or vary widely in symptom severity and age of onset. This complicates diagnosis and disease management in MFD. In addition to its clinical importance, knowledge of genetic modifiers that confer penetrance and expressivity variability opens possibilities to identify protective genetic variants which may help probe the mechanisms of more common forms of diabetes and shed light in new therapeutic strategies. In this review, we will mainly address penetrance and expressivity variation in different types of MFD, factors that confer such variations and opportunities that come with such knowledge. Related literature was searched in PubMed, Medline and Embase. Papers with publication year from 1974 to 2023 are included. Data are either sourced from literatures or from OMIM, Clinvar and 1000 genome browser.

Seven naturally variant loci serve as genetic modifiers of Lamc2jeb induced non-Herlitz junctional Epidermolysis Bullosa in mice

Epidermolysis Bullosa (EB) is a group of rare genetic disorders that compromise the structural integrity of the skin such that blisters and subsequent erosions occur after minor trauma. While primary genetic risk of all subforms of EB adhere to Mendelian patterns of inheritance, their clinical presentations and severities can vary greatly, implying genetic modifiers. The Lamc2jeb mouse model of non-Herlitz junctional EB (JEB-nH) demonstrated that genetic modifiers can contribute substantially to the phenotypic variability of JEB and likely other forms of EB. The innocuous changes in an 'EB related gene', Col17a1, have shown it to be a dominant modifier of Lamc2jeb. This work identifies six additional Quantitative Trait Loci (QTL) that modify disease in Lamc2jeb/jeb mice. Three QTL include other known 'EB related genes', with the strongest modifier effect mapping to a region including the epidermal hemi-desmosomal structural gene dystonin (Dst-e/Bpag1-e). Three other QTL map to intervals devoid of known EB-associated genes. Of these, one contains the nuclear receptor coactivator Ppargc1a as its primary candidate and the others contain related genes Pparg and Igf1, suggesting modifier pathways. These results, demonstrating the potent disease modifying effects of normally innocuous genetic variants, greatly expand the landscape of genetic modifiers of EB and therapeutic approaches that may be applied.

Phenotypic effects of genetic variants associated with autism

While over 100 genes have been associated with autism, little is known about the prevalence of variants affecting them in individuals without a diagnosis of autism. Nor do we fully appreciate the phenotypic diversity beyond the formal autism diagnosis. Based on data from more than 13,000 individuals with autism and 210,000 undiagnosed individuals, we estimated the odds ratios for autism associated to rare loss-of-function (LoF) variants in 185 genes associated with autism, alongside 2,492 genes displaying intolerance to LoF variants. In contrast to autism-centric approaches, we investigated the correlates of these variants in individuals without a diagnosis of autism. We show that these variants are associated with a small but significant decrease in fluid intelligence, qualification level and income and an increase in metrics related to material deprivation. These effects were larger for autism-associated genes than in other LoF-intolerant genes. Using brain imaging data from 21,040 individuals from the UK Biobank, we could not detect significant differences in the overall brain anatomy between LoF carriers and non-carriers. Our results highlight the importance of studying the effect of the genetic variants beyond categorical diagnosis and the need for more research to understand the association between these variants and sociodemographic factors, to best support individuals carrying these variants.

Insights into the genetic architecture of congenital heart disease from animal modeling

Congenital heart disease (CHD) is observed in up to 1% of live births and is one of the leading causes of mortality from birth defects. While hundreds of genes have been implicated in the genetic etiology of CHD, their role in CHD pathogenesis is still poorly understood. This is largely a reflection of the sporadic nature of CHD, as well as its variable expressivity and incomplete penetrance. We reviewed the monogenic causes and evidence for oligogenic etiology of CHD, as well as the role of de novo mutations, common variants, and genetic modifiers. For further mechanistic insight, we leveraged single-cell data across species to investigate the cellular expression characteristics of genes implicated in CHD in developing human and mouse embryonic hearts. Understanding the genetic etiology of CHD may enable the application of precision medicine and prenatal diagnosis, thereby facilitating early intervention to improve outcomes for patients with CHD.

Widespread epistasis among beneficial genetic variants revealed by high-throughput genome editing

The phenotypic effect of any genetic variant can be altered by variation at other genomic loci. Known as epistasis, these genetic interactions shape the genotype-phenotype map of every species, yet their origins remain poorly understood. To investigate this, we employed high-throughput genome editing to measure the fitness effects of 1,826 naturally polymorphic variants in four strains of Saccharomyces cerevisiae. About 31% of variants affect fitness, of which 24% have strain-specific fitness effects indicative of epistasis. We found that beneficial variants are more likely to exhibit genetic interactions and that these interactions can be mediated by specific traits such as flocculation ability. This work suggests that adaptive evolution will often involve trade-offs where a variant is only beneficial in some genetic backgrounds, potentially explaining why many beneficial variants remain polymorphic. In sum, we provide a framework to understand the factors influencing epistasis with single-nucleotide resolution, revealing widespread epistasis among beneficial variants.

Clinical Features of Genetic Resilience in Chronic Obstructive Pulmonary Disease

Introduction

In the personalized risk quantification of chronic obstructive pulmonary disease (COPD), genome-wide association studies and polygenic risk scores (PRS) complement traditional risk factors, such as age and cigarette smoking. However, despite being at considerable levels of risk, some individuals do not develop COPD. Research on COPD resilience remains largely unexplored.

Methods

We applied the previously published COPD PRS to whole genome sequencing data from non-Hispanic white and African American individuals in the COPDGene study. We defined genetic resilience as individuals unaffected by COPD with a polygenic risk score above the 90 ^th percentile. We defined risk-matched case individuals as those with COPD (i.e., FEV ₁ /FVC < 0.70) and a PRS above the 90 ^th percentile. We defined low risk individuals without COPD (i.e., FEV ₁ /FVC > 0.70) as a polygenic risk score below the 10 ^th percentile. We compared genetically resilient individuals to risk-matched individuals with COPD and low risk individuals by demographics, lung function, respiratory symptoms, co-morbidities, and chest CT scan measurements. We also performed survival analyses, differential expression analysis, and matching for sensitivity analyses.

Results

We identified 211 resilient individuals without COPD, 605 genetic risk-matched individuals with COPD, and 527 low-risk individuals without COPD. Resilient individuals had higher FEV ₁ % predicted and lower percent emphysema. In contrast, resilient individuals had higher airway wall thickness compared to low-risk unaffected individuals. While there was no difference in survival between low-risk and resilient individuals, resilient individuals had higher survival compared to risk matched cases. We also identified two genes that were differentially expressed between low-risk unaffected individuals and resilient individuals.

Conclusion

Genetically resilient individuals had a reduced burden of COPD disease-related measures compared to risk-matched cases but had subtly increased measures compared to low-risk unaffected individuals. Further genetic studies will be needed to illuminate the underlying pathobiology of our observations.

Noncancer-related Secondary Findings in a Cohort of 231 Children With Cancer and Their Parents

Application of next-generation sequencing may lead to the detection of secondary findings (SF) not related to the initially analyzed disease but to other severe medically actionable diseases. However, the analysis of SFs is not yet routinely performed. We mined whole-exome sequencing data of 231 pediatric cancer patients and their parents who had been treated in our center for the presence of SFs. By this approach, we identified in 6 children (2.6%) pathogenic germline variants in 5 of the noncancer-related genes on the American College of Medical Genetics and Genomics (ACMG) SF v3.0 list, of which the majority were related to cardiovascular diseases ( RYR2 , MYBPC3 , KCNQ1 ). Interestingly, only the patient harboring the KCNQ1 variant showed at the time point of the analysis signs of the related Long QT syndrome. Moreover, we report 3 variants of unknown significance which, although not classified as pathogenic, have been reported in the literature to occur in individuals with the respective disease. While the frequency of patients with SFs is low, the impact of such findings on the patients' life is enormous, with regard to the potential prevention of life-threatening diseases. Hence, we are convinced that such actionable SF should be routinely analyzed.

Causal relationship between gut microbiota and cancers: a two-sample Mendelian randomisation study

BACKGROUND

Evidence from observational studies and clinical trials suggests that the gut microbiota is associated with cancer. However, the causal association between gut microbiota and cancer remains to be determined.

METHODS

We first identified two sets of gut microbiota based on phylum, class, order, family, and genus level information, and cancer data were obtained from the IEU Open GWAS project. We then performed two-sample Mendelian randomisation (MR) to determine whether the gut microbiota is causally associated with eight cancer types. Furthermore, we performed a bi-directional MR analysis to examine the direction of the causal relations.

RESULTS

We identified 11 causal relationships between genetic liability in the gut microbiome and cancer, including those involving the genus Bifidobacterium. We found 17 strong associations between genetic liability in the gut microbiome and cancer. Moreover, we found 24 associations between genetic liability in the gut microbiome and cancer using multiple datasets.

CONCLUSIONS

Our MR analysis revealed that the gut microbiota was causally associated with cancers and may be useful in providing new insights for further mechanistic and clinical studies of microbiota-mediated cancer.

Incidental molecular diagnoses and heterozygous risk alleles in a carrier screening cohort

PURPOSE

Expanded pan-ethnic carrier screening is an effective tool for the management of reproductive risk. However, growth in the number of conditions screened, in combination with increasingly more comprehensive test methodologies, can lead to the detection of genetic findings that may affect the health of the tested individual. The objective of this study was to investigate the frequency of pathogenic genotypes in a presumed healthy carrier screening cohort to facilitate broader discussions regarding disclosure of genetic information from carrier screening.

METHODS

A retrospective analysis of 73,755 targeted carrier screens was performed to identify individuals with pathogenic genotypes and heterozygous risk alleles.

RESULTS

In this study, we identified 79 individuals (0.11%) with pathogenic genotypes associated with moderate to profound autosomal recessive or X-linked conditions. In addition, 10 cases had chromosome X dosage abnormalities suggestive of a sex chromosome abnormality. Heterozygote risk alleles represented the majority of ancillary findings in this cohort, including 280 female carriers of FMR1 premutation alleles, 15 heterozygous females with pathogenic DMD variants, and 174 heterozygotes with pathogenic variants in genes that may confer increased risk for somatic malignancies in the heterozygous state.

CONCLUSION

These data suggest that nearly 1% of individuals undergoing carrier screening will have a finding that may require clinical evaluation or surveillance.

Genetic control of mRNA splicing as a potential mechanism for incomplete penetrance of rare coding variants

Exonic variants present some of the strongest links between genotype and phenotype. However, these variants can have significant inter-individual pathogenicity differences, known as variable penetrance. In this study, we propose a model where genetically controlled mRNA splicing modulates the pathogenicity of exonic variants. By first cataloging exonic inclusion from RNA-seq data in GTEx v8, we find that pathogenic alleles are depleted on highly included exons. Using a large-scale phased WGS data from the TOPMed consortium, we observe that this effect may be driven by common splice-regulatory genetic variants, and that natural selection acts on haplotype configurations that reduce the transcript inclusion of putatively pathogenic variants, especially when limiting to haploinsufficient genes. Finally, we test if this effect may be relevant for autism risk using families from the Simons Simplex Collection, but find that splicing of pathogenic alleles has a penetrance reducing effect here as well. Overall, our results indicate that common splice-regulatory variants may play a role in reducing the damaging effects of rare exonic variants.

MED13 and glycolysis are conserved modifiers of α-synuclein-associated neurodegeneration

α-Synuclein (α-syn) is important in synucleinopathies such as Parkinson's disease (PD). While genome-wide association studies (GWASs) of synucleinopathies have identified many risk loci, the underlying genes have not been shown for most loci. Using Drosophila, we screened 3,471 mutant chromosomes for genetic modifiers of α-synuclein and identified 12 genes. Eleven modifiers have human orthologs associated with diseases, including MED13 and CDC27, which lie within PD GWAS loci. Drosophila Skd/Med13 and glycolytic enzymes are co-upregulated by α-syn-associated neurodegeneration. While elevated α-syn compromises mitochondrial function, co-expressing skd/Med13 RNAi and α-syn synergistically increase the ratio of oxidized-to-reduced glutathione. The resulting neurodegeneration can be suppressed by overexpressing a glycolytic enzyme or treatment with deferoxamine, suggesting that compensatory glycolysis is neuroprotective. In addition, the functional relationship between α-synuclein, MED13, and glycolytic enzymes is conserved between flies and mice. We propose that hypoxia-inducible factor and MED13 are part of a druggable pathway for PD.

Whole exome sequencing in dense families suggests genetic pleiotropy amongst Mendelian and complex neuropsychiatric syndromes

Whole Exome Sequencing (WES) studies provide important insights into the genetic architecture of serious mental illness (SMI). Genes that are central to the shared biology of SMIs may be identified by WES in families with multiple affected individuals with diverse SMI (F-SMI). We performed WES in 220 individuals from 75 F-SMI families and 60 unrelated controls. Within pedigree prioritization employed criteria of rarity, functional consequence, and sharing by ≥ 3 affected members. Across the sample, gene and gene-set-wide case-control association analysis was performed with Sequence Kernel Association Test (SKAT). In 14/16 families with ≥ 3 sequenced affected individuals, we identified a total of 78 rare predicted deleterious variants in 78 unique genes shared by ≥ 3 members with SMI. Twenty (25%) genes were implicated in monogenic CNS syndromes in OMIM (OMIM-CNS), a fraction that is a significant overrepresentation (Fisher's Exact test OR = 2.47, p = 0.001). In gene-set SKAT, statistically significant association was noted for OMIM-CNS gene-set (SKAT-p = 0.005) but not the synaptic gene-set (SKAT-p = 0.17). In this WES study in F-SMI, we identify private, rare, protein altering variants in genes previously implicated in Mendelian neuropsychiatric syndromes; suggesting pleiotropic influences in neurodevelopment between complex and Mendelian syndromes.

Scientific opportunities in resilience research for cardiovascular health and wellness. Report from a National Heart, Lung, and Blood Institute workshop

Exposure of biological systems to acute or chronic insults triggers a host of molecular and physiological responses to either tolerate, adapt, or fully restore homeostasis; these responses constitute the hallmarks of resilience. Given the many facets, dimensions, and discipline-specific focus, gaining a shared understanding of "resilience" has been identified as a priority for supporting advances in cardiovascular health. This report is based on the working definition: "Resilience is the ability of living systems to successfully maintain or return to homeostasis in response to physical, molecular, individual, social, societal, or environmental stressors or challenges," developed after considering many factors contributing to cardiovascular resilience through deliberations of multidisciplinary experts convened by the National Heart, Lung, and Blood Institute during a workshop entitled: "Enhancing Resilience for Cardiovascular Health and Wellness." Some of the main emerging themes that support the possibility of enhancing resilience for cardiovascular health include optimal energy management and substrate diversity, a robust immune system that safeguards tissue homeostasis, and social and community support. The report also highlights existing research challenges, along with immediate and long-term opportunities for resilience research. Certain immediate opportunities identified are based on leveraging existing high-dimensional data from longitudinal clinical studies to identify vascular resilience measures, create a 'resilience index,' and adopt a life-course approach. Long-term opportunities include developing quantitative cell/organ/system/community models to identify resilience factors and mechanisms at these various levels, designing experimental and clinical interventions that specifically assess resilience, adopting global sharing of resilience-related data, and cross-domain training of next-generation researchers in this field.

Developmental disruption and restoration of brain synaptome architecture in the murine Pax6 neurodevelopmental disease model

Neurodevelopmental disorders of genetic origin delay the acquisition of normal abilities and cause disabling phenotypes. Nevertheless, spontaneous attenuation and even complete amelioration of symptoms in early childhood and adolescence can occur in many disorders, suggesting that brain circuits possess an intrinsic capacity to overcome the deficits arising from some germline mutations. We examined the molecular composition of almost a trillion excitatory synapses on a brain-wide scale between birth and adulthood in mice carrying a mutation in the homeobox transcription factor Pax6, a neurodevelopmental disorder model. Pax6 haploinsufficiency had no impact on total synapse number at any age. By contrast, the molecular composition of excitatory synapses, the postnatal expansion of synapse diversity and the acquisition of normal synaptome architecture were delayed in all brain regions, interfering with networks and electrophysiological simulations of cognitive functions. Specific excitatory synapse types and subtypes were affected in two key developmental age-windows. These phenotypes were reversed within 2-3 weeks of onset, restoring synapse diversity and synaptome architecture to the normal developmental trajectory. Synapse subtypes with rapid protein turnover mediated the synaptome remodeling. This brain-wide capacity for remodeling of synapse molecular composition to recover and maintain the developmental trajectory of synaptome architecture may help confer resilience to neurodevelopmental genetic disorders.

Clinical genomics and precision medicine

Precision Medicine emerges from the genomic paradigm of health and disease. For precise molecular diagnoses of genetic diseases, we must analyze the Whole Exome (WES) or the Whole Genome (WGS). By not needing exon capture, WGS is more powerful to detect single nucleotide variants and copy number variants. In healthy individuals, we can observe monogenic highly penetrant variants, which may be causally responsible for diseases, and also susceptibility variants, associated with common polygenic diseases. But there is the major problem of penetrance. Thus, there is the question of whether it is worthwhile to perform WGS in all healthy individuals as a step towards Precision Medicine. The genetic architecture of disease is consistent with the fact that they are all polygenic. Moreover, ancestry adds another layer of complexity. We are now capable of obtaining Polygenic Risk Scores for all complex diseases using only data from new generation sequencing. Yet, review of available evidence does not at present favor the idea that WGS analyses are sufficiently developed to allow reliable predictions of the risk components for monogenic and polygenic hereditary diseases in healthy individuals. Probably, it is still better for WGS to remain reserved for the diagnosis of pathogenic variants of Mendelian diseases.

Loss-of-function mutation survey revealed that genes with background-dependent fitness are rare and functionally related in yeast

In different individuals, the same mutation can lead to different phenotypes due to genetic background effects. This is commonly observed in various systems, including many human diseases. While isolated examples of such background effects have been observed, a systematic view across a large number of individuals is still lacking. Here, we surveyed genetic background effects associated with gene loss-of-function mutations across a population of natural isolates of the yeast Saccharomyces cerevisiae. We found that ∼15% of genes can display a background-dependent fitness change. Genes related to mitochondrial functions are significantly overrepresented, and showed reversed patterns of fitness gain or loss with genes involved in transcription and chromatin remodeling as well as in nuclear–cytoplasmic transport, suggesting a potential functional rewiring.

In natural populations, the same mutation can lead to different phenotypic outcomes due to the genetic variation that exists among individuals. Such genetic background effects are commonly observed, including in the context of many human diseases. However, systematic characterization of these effects at the species level is still lacking to date. Here, we sought to comprehensively survey background-dependent traits associated with gene loss-of-function (LoF) mutations in 39 natural isolates of Saccharomyces cerevisiae using a transposon saturation strategy. By analyzing the modeled fitness variability of a total of 4,469 genes, we found that 15% of them, when impacted by a LoF mutation, exhibited a significant gain- or loss-of-fitness phenotype in certain natural isolates compared with the reference strain S288C. Out of these 632 genes with predicted background-dependent fitness effects, around 2/3 impact multiple backgrounds with a gradient of predicted fitness change while 1/3 are specific to a single genetic background. Genes related to mitochondrial function are significantly overrepresented in the set of genes showing a continuous variation and display a potential functional rewiring with other genes involved in transcription and chromatin remodeling as well as in nuclear–cytoplasmic transport. Such rewiring effects are likely modulated by both the genetic background and the environment. While background-specific cases are rare and span diverse cellular processes, they can be functionally related at the individual level. All genes with background-dependent fitness effects tend to have an intermediate connectivity in the global genetic interaction network and have shown relaxed selection pressure at the population level, highlighting their potential evolutionary characteristics.

T. A genome sequencing system for universal newborn screening, diagnosis, and precision medicine for severe genetic diseases

Newborn screening (NBS) dramatically improves outcomes in severe childhood disorders by treatment before symptom onset. In many genetic diseases, however, outcomes remain poor because NBS has lagged behind drug development. Rapid whole-genome sequencing (rWGS) is attractive for comprehensive NBS because it concomitantly examines almost all genetic diseases and is gaining acceptance for genetic disease diagnosis in ill newborns. We describe prototypic methods for scalable, parentally consented, feedback-informed NBS and diagnosis of genetic diseases by rWGS and virtual, acute management guidance (NBS-rWGS). Using established criteria and the Delphi method, we reviewed 457 genetic diseases for NBS-rWGS, retaining 388 (85%) with effective treatments. Simulated NBS-rWGS in 454,707 UK Biobank subjects with 29,865 pathogenic or likely pathogenic variants associated with 388 disorders had a true negative rate (specificity) of 99.7% following root cause analysis. In 2,208 critically ill children with suspected genetic disorders and 2,168 of their parents, simulated NBS-rWGS for 388 disorders identified 104 (87%) of 119 diagnoses previously made by rWGS and 15 findings not previously reported (NBS-rWGS negative predictive value 99.6%, true positive rate [sensitivity] 88.8%). Retrospective NBS-rWGS diagnosed 15 children with disorders that had been undetected by conventional NBS. In 43 of the 104 children, had NBS-rWGS-based interventions been started on day of life 5, the Delphi consensus was that symptoms could have been avoided completely in seven critically ill children, mostly in 21, and partially in 13. We invite groups worldwide to refine these NBS-rWGS conditions and join us to prospectively examine clinical utility and cost effectiveness.

Polygenic Resilience Modulates the Penetrance of Parkinson Disease Genetic Risk Factors

OBJECTIVE

The aim of the current study is to understand why some individuals avoid developing Parkinson disease (PD) despite being at relatively high genetic risk, using the largest datasets of individual-level genetic data available.

METHODS

We calculated polygenic risk score to identify controls and matched PD cases with the highest burden of genetic risk for PD in the discovery cohort (International Parkinson's Disease Genomics Consortium, 7,204 PD cases and 9,412 controls) and validation cohorts (Comprehensive Unbiased Risk Factor Assessment for Genetics and Environment in Parkinson's Disease, 8,968 cases and 7,598 controls; UK Biobank, 2,639 PD cases and 14,301 controls; Accelerating Medicines Partnership-Parkinson's Disease Initiative, 2,248 cases and 2,817 controls). A genome-wide association study meta-analysis was performed on these individuals to understand genetic variation associated with resistance to disease. We further constructed a polygenic resilience score, and performed multimarker analysis of genomic annotation (MAGMA) gene-based analyses and functional enrichment analyses.

RESULTS

A higher polygenic resilience score was associated with a lower risk for PD (β = -0.054, standard error [SE] = 0.022, p = 0.013). Although no single locus reached genome-wide significance, MAGMA gene-based analyses nominated TBCA as a putative gene. Furthermore, we estimated the narrow-sense heritability associated with resilience to PD (h2  = 0.081, SE = 0.035, p = 0.0003). Subsequent functional enrichment analysis highlighted histone methylation as a potential pathway harboring resilience alleles that could mitigate the effects of PD risk loci.

INTERPRETATION

The present study represents a novel and comprehensive assessment of heritable genetic variation contributing to PD resistance. We show that a genetic resilience score can modify the penetrance of PD genetic risk factors and therefore protect individuals carrying a high-risk genetic burden from developing PD. ANN NEUROL 2022;92:270-278.

Common genetic variation associated with Mendelian disease severity revealed through cryptic phenotype analysis

Clinical heterogeneity is common in Mendelian disease, but small sample sizes make it difficult to identify specific contributing factors. However, if a disease represents the severely affected extreme of a spectrum of phenotypic variation, then modifier effects may be apparent within a larger subset of the population. Analyses that take advantage of this full spectrum could have substantially increased power. To test this, we developed cryptic phenotype analysis, a model-based approach that infers quantitative traits that capture disease-related phenotypic variability using qualitative symptom data. By applying this approach to 50 Mendelian diseases in two cohorts, we identify traits that reliably quantify disease severity. We then conduct genome-wide association analyses for five of the inferred cryptic phenotypes, uncovering common variation that is predictive of Mendelian disease-related diagnoses and outcomes. Overall, this study highlights the utility of computationally-derived phenotypes and biobank-scale cohorts for investigating the complex genetic architecture of Mendelian diseases.

Inherited variants in CHD3 show variable expressivity in Snijders Blok-Campeau syndrome

PURPOSE

Common diagnostic next-generation sequencing strategies are not optimized to identify inherited variants in genes associated with dominant neurodevelopmental disorders as causal when the transmitting parent is clinically unaffected, leaving a significant number of cases with neurodevelopmental disorders undiagnosed.

METHODS

We characterized 21 families with inherited heterozygous missense or protein-truncating variants in CHD3, a gene in which de novo variants cause Snijders Blok-Campeau syndrome.

RESULTS

Computational facial and Human Phenotype Ontology-based comparisons showed that the phenotype of probands with inherited CHD3 variants overlaps with the phenotype previously associated with de novo CHD3 variants, whereas heterozygote parents are mildly or not affected, suggesting variable expressivity. In addition, similarly reduced expression levels of CHD3 protein in cells of an affected proband and of healthy family members with a CHD3 protein-truncating variant suggested that compensation of expression from the wild-type allele is unlikely to be an underlying mechanism. Notably, most inherited CHD3 variants were maternally transmitted.

CONCLUSION

Our results point to a significant role of inherited variation in Snijders Blok-Campeau syndrome, a finding that is critical for correct variant interpretation and genetic counseling and warrants further investigation toward understanding the broader contributions of such variation to the landscape of human disease.

Using Stem Cell Models to Explore the Genetics Underlying Psychiatric Disorders: Linking Risk Variants, Genes, and Biology in Brain Disease

There is an urgent and unmet need to advance our ability to translate genetic studies of psychiatric disorders into clinically actionable information, which could transform diagnostics and even one day lead to novel (and potentially presymptomatic) therapeutic interventions. Today, although there are hundreds of significant loci associated with psychiatric disorders, resolving the target gene(s) and pathway(s) impacted by each is a major challenge. Integrating human induced pluripotent stem cell-based approaches with CRISPR-mediated genomic engineering strategies makes it possible to study the impact of patient-specific variants within the cell types of the brain. As the scale and scope of functional genomic studies expands, so does our ability to resolve the complex interplay of the many risk variants linked to psychiatric disorders. In this review, the author discusses some of the technological advances that make it possible to ask exciting questions that are fundamental to our understanding of psychiatric disorders. How do distinct risk variants converge and interact with each other (and the environment) across the diverse cell types that comprise the human brain? Can clinical trajectories and/or therapeutic response be predicted from genetic profiles? Just as critically, by spreading the message that genetic risk for psychiatric disorders is biological and fundamentally no different than for other human conditions, we can dispel the stigma associated with mental illness.

The individual and global impact of copy-number variants on complex human traits

The impact of copy-number variations (CNVs) on complex human traits remains understudied. We called CNVs in 331,522 UK Biobank participants and performed genome-wide association studies (GWASs) between the copy number of CNV-proxy probes and 57 continuous traits, revealing 131 signals spanning 47 phenotypes. Our analysis recapitulated well-known associations (e.g., 1q21 and height), revealed the pleiotropy of recurrent CNVs (e.g., 26 and 16 traits for 16p11.2-BP4-BP5 and 22q11.21, respectively), and suggested gene functionalities (e.g., MARF1 in female reproduction). Forty-eight CNV signals (38%) overlapped with single-nucleotide polymorphism (SNP)-GWASs signals for the same trait. For instance, deletion of PDZK1, which encodes a urate transporter scaffold protein, decreased serum urate levels, while deletion of RHD, which encodes the Rhesus blood group D antigen, associated with hematological traits. Other signals overlapped Mendelian disorder regions, suggesting variable expressivity and broad impact of these loci, as illustrated by signals mapping to Rotor syndrome (SLCO1B1/3), renal cysts and diabetes syndrome (HNF1B), or Charcot-Marie-Tooth (PMP22) loci. Total CNV burden negatively impacted 35 traits, leading to increased adiposity, liver/kidney damage, and decreased intelligence and physical capacity. Thirty traits remained burden associated after correcting for CNV-GWAS signals, pointing to a polygenic CNV architecture. The burden negatively correlated with socio-economic indicators, parental lifespan, and age (survivorship proxy), suggesting a contribution to decreased longevity. Together, our results showcase how studying CNVs can expand biological insights, emphasizing the critical role of this mutational class in shaping human traits and arguing in favor of a continuum between Mendelian and complex diseases.

Complex Changes in the Efficiency of the Expression of Many Genes in Monogenic Diseases, Mucopolysaccharidoses, May Arise from Significant Disturbances in the Levels of Factors Involved in the Gene Expression Regulation Processes

Monogenic diseases are primarily caused by mutations in a single gene; thus, they are commonly recognized as genetic disorders with the simplest mechanisms. However, recent studies have indicated that the molecular mechanisms of monogenic diseases can be unexpectedly complicated, and their understanding requires complex studies at the molecular level. Previously, we have demonstrated that in mucopolysaccharidoses (MPS), a group of monogenic lysosomal storage diseases, several hundreds of genes reveal significant changes in the expression of various genes. Although the secondary effects of the primary biochemical defect and the inefficient degradation of glycosaminoglycans (GAGs) might be considered, the scale of the changes in the expression of a large fraction of genes cannot be explained by a block in one biochemical pathway. Here, we demonstrate that in cellular models of 11 types of MPS, the expression of genes coding for proteins involved in the regulation of the expression of many other genes at various stages (such as signal transduction, transcription, splicing, RNA degradation, translation, and others) is significantly disturbed relative to the control cells. This conclusion was based on transcriptomic studies, supported by biochemical analyses of levels of selected proteins encoded by genes revealing an especially high level of dysregulation in MPS (EXOSC9, SRSF10, RPL23, and NOTCH3 proteins were investigated). Interestingly, the reduction in GAGs levels, through the inhibition of their synthesis normalized the amounts of EXOSC9, RPL23, and NOTCH3 in some (but not all) MPS types, while the levels of SRSF10 could not be corrected in this way. These results indicate that different mechanisms are involved in the dysregulation of the expression of various genes in MPS, pointing to a potential explanation for the inability of some therapies (such as enzyme replacement therapy or substrate reduction therapy) to fully correct the physiology of MPS patients. We suggest that the disturbed expression of some genes, which appears as secondary or tertiary effects of GAG storage, might not be reversible, even after a reduction in the amounts of the storage material.

Whole-genome sequencing of 1,171 elderly admixed individuals from São Paulo, Brazil

As whole-genome sequencing (WGS) becomes the gold standard tool for studying population genomics and medical applications, data on diverse non-European and admixed individuals are still scarce. Here, we present a high-coverage WGS dataset of 1,171 highly admixed elderly Brazilians from a census-based cohort, providing over 76 million variants, of which ~2 million are absent from large public databases. WGS enables identification of ~2,000 previously undescribed mobile element insertions without previous description, nearly 5 Mb of genomic segments absent from the human genome reference, and over 140 alleles from HLA genes absent from public resources. We reclassify and curate pathogenicity assertions for nearly four hundred variants in genes associated with dominantly-inherited Mendelian disorders and calculate the incidence for selected recessive disorders, demonstrating the clinical usefulness of the present study. Finally, we observe that whole-genome and HLA imputation could be significantly improved compared to available datasets since rare variation represents the largest proportion of input from WGS. These results demonstrate that even smaller sample sizes of underrepresented populations bring relevant data for genomic studies, especially when exploring analyses allowed only by WGS.

Whole genome sequencing (WGS) data on non-European and admixed individuals remains scarce. Here, the authors analyse WGS data from 1,171 admixed elderly Brazilians from a census cohort, characterising population-specific genetic variation and exploring the clinical utility of this expanded dataset.

Genetic resiliency associated with dominant lethal TPM1 mutation causing atrial septal defect with high heritability

Analysis of large-scale human genomic data has yielded unexplained mutations known to cause severe disease in healthy individuals. Here, we report the unexpected recovery of a rare dominant lethal mutation in TPM1, a sarcomeric actin-binding protein, in eight individuals with large atrial septal defect (ASD) in a five-generation pedigree. Mice with Tpm1 mutation exhibit early embryonic lethality with disrupted myofibril assembly and no heartbeat. However, patient-induced pluripotent-stem-cell-derived cardiomyocytes show normal beating with mild myofilament defect, indicating disease suppression. A variant in TLN2, another myofilament actin-binding protein, is identified as a candidate suppressor. Mouse CRISPR knock-in (KI) of both the TLN2 and TPM1 variants rescues heart beating, with near-term fetuses exhibiting large ASD. Thus, the role of TPM1 in ASD pathogenesis unfolds with suppression of its embryonic lethality by protective TLN2 variant. These findings provide evidence that genetic resiliency can arise with genetic suppression of a deleterious mutation.

Genetic basis of a spontaneous mutation’s expressivity

Genetic background often influences the phenotypic consequences of mutations, resulting in variable expressivity. How standing genetic variants collectively cause this phenomenon is not fully understood. Here, we comprehensively identify loci in a budding yeast cross that impact the growth of individuals carrying a spontaneous missense mutation in the nuclear-encoded mitochondrial ribosomal gene MRP20. Initial results suggested that a single large effect locus influences the mutation’s expressivity, with one allele causing inviability in mutants. However, further experiments revealed this simplicity was an illusion. In fact, many additional loci shape the mutation’s expressivity, collectively leading to a wide spectrum of mutational responses. These results exemplify how complex combinations of alleles can produce a diversity of qualitative and quantitative responses to the same mutation.

Developmental disorders caused by haploinsufficiency of transcriptional regulators: a perspective based on cell fate determination

Many human birth defects and neurodevelopmental disorders are caused by loss-of-function mutations in a single copy of transcription factor (TF) and chromatin regulator genes. Although this dosage sensitivity has long been known, how and why haploinsufficiency (HI) of transcriptional regulators leads to developmental disorders (DDs) is unclear. Here I propose the hypothesis that such DDs result from defects in cell fate determination that are based on disrupted bistability in the underlying gene regulatory network (GRN). Bistability, a crucial systems biology concept to model binary choices such as cell fate decisions, requires both positive feedback and ultrasensitivity, the latter often achieved through TF cooperativity. The hypothesis explains why dosage sensitivity of transcriptional regulators is an inherent property of fate decisions, and why disruption of either positive feedback or cooperativity in the underlying GRN is sufficient to cause disease. I present empirical and theoretical evidence in support of this hypothesis and discuss several issues for which it increases our understanding of disease, such as incomplete penetrance. The proposed framework provides a mechanistic, systems-level explanation of HI of transcriptional regulators, thus unifying existing theories, and offers new insights into outstanding issues of human disease. This article has an associated Future Leader to Watch interview with the author of the paper.

Low frequency of treatable pediatric disease alleles in gnomAD: An opportunity for future genomic screening of newborns

Hematopoietic stem cell transplant (HSCT) can prevent progression of several genetic disorders. Although a subset of these disorders are identified on newborn screening panels, others are not identified until irreversible symptoms develop. Genetic testing is an efficient methodology to ascertain pre-symptomatic children, but the penetrance of risk-associated variants in the general population is not well understood. We developed a list of 127 genes associated with disorders treatable with HSCT. We identified likely pathogenic or pathogenic (LP/P) and loss-of-function (LoF) variants in these genes in the Genome Aggregation Database (gnomAD), a dataset containing exome and genome sequencing data from 141,456 healthy adults. Within gnomAD, we identified 59 individuals with a LP/P or LoF variant in 15 genes. Genes were associated with bone marrow failure syndromes, bleeding disorders, primary immunodeficiencies, osteopetrosis, metabolic disorders, and epidermolysis bullosa. In conclusion, few ostensibly healthy adults had genotypes associated with pediatric disorders treatable with HSCTs. Given that most of these disorders do not have biomarkers that could be cheaply and universally assessed on a standard newborn screen, our data suggest that genetic testing may be a complementary approach to traditional newborn screening methodology that has the potential to improve mortality and is not expected to lead to a high burden of false-positive results.

Incomplete Penetrance and Variable Expressivity: From Clinical Studies to Population Cohorts

The same genetic variant found in different individuals can cause a range of diverse phenotypes, from no discernible clinical phenotype to severe disease, even among related individuals. Such variants can be said to display incomplete penetrance, a binary phenomenon where the genotype either causes the expected clinical phenotype or it does not, or they can be said to display variable expressivity, in which the same genotype can cause a wide range of clinical symptoms across a spectrum. Both incomplete penetrance and variable expressivity are thought to be caused by a range of factors, including common variants, variants in regulatory regions, epigenetics, environmental factors, and lifestyle. Many thousands of genetic variants have been identified as the cause of monogenic disorders, mostly determined through small clinical studies, and thus, the penetrance and expressivity of these variants may be overestimated when compared to their effect on the general population. With the wealth of population cohort data currently available, the penetrance and expressivity of such genetic variants can be investigated across a much wider contingent, potentially helping to reclassify variants that were previously thought to be completely penetrant. Research into the penetrance and expressivity of such genetic variants is important for clinical classification, both for determining causative mechanisms of disease in the affected population and for providing accurate risk information through genetic counseling. A genotype-based definition of the causes of rare diseases incorporating information from population cohorts and clinical studies is critical for our understanding of incomplete penetrance and variable expressivity. This review examines our current knowledge of the penetrance and expressivity of genetic variants in rare disease and across populations, as well as looking into the potential causes of the variation seen, including genetic modifiers, mosaicism, and polygenic factors, among others. We also considered the challenges that come with investigating penetrance and expressivity.

Driving mosaicism: somatic variants in reference population databases and effect on variant interpretation in rare genetic disease

BACKGROUND

Genetic variation databases provide invaluable information on the presence and frequency of genetic variants in the 'untargeted' human population, aggregated with the primary goal to facilitate the interpretation of clinically important variants. The presence of somatic variants in such databases can affect variant assessment in undiagnosed rare disease (RD) patients. Previously, the impact of somatic mosaicism was only considered in relation to two Mendelian disease-associated genes. Here, we expand the analyses to identify additional mosaicism-prone genes in blood-derived reference population databases.

RESULTS

To identify additional mosaicism-prone genes relevant to RDs, we focused on known/previously established ClinVar pathogenic and likely pathogenic single-nucleotide variants, residing in genes associated with early onset, severe autosomal dominant diseases. We asked whether any of these variants are present in a higher-than-expected frequency in the reference population databases and whether there is evidence of somatic origin (i.e., allelic imbalance) rather than germline heterozygosity (~ half of the reads supporting alternative allele). The mosaicism-prone genes identified were further categorized according to the processes they are involved in. Beyond the previously reported ASXL1 and DNMT3A, we identified 7 additional autosomal dominant RD-associated genes with known pathogenic single-nucleotide variants present in the reference population databases and good evidence of allelic imbalance: BRAF, CBL, FGFR3, IDH2, KRAS, PTPN11 and SETBP1. From this group of 9 genes, the majority (n = 7) was important for hematopoiesis. In addition, 4 of these genes were involved in cell proliferation. Further assessment of the known 156 hematopoietic genes led to identification of 48 genes (21 not yet associated with RDs) with at least some evidence of mosaicism detectable in reference population databases.

CONCLUSIONS

These results stress the importance of considering genes involved in hematopoiesis and cell proliferation when interpreting the presence and frequency of genetic variants in blood-derived reference population databases, both public and private. This is especially important when considering new variants of uncertain significance in known hematopoietic/cell proliferation RD genes and future novel gene-disease associations involving this class of genes.

The good, the bad, and the ugly: Evolutionary and pathological aspects of gene dosage alterations

Diploid organisms contain a maternal and a paternal genome complement that is thought to provide robustness and allow developmental progression despite genetic perturbations that occur in heterozygosity. However, changes affecting gene dosage from the chromosome down to the individual gene level possess a significant pathological potential and can lead to developmental disorders (DDs). This indicates that expression from a balanced gene complement is highly relevant for proper cellular and organismal function in eukaryotes. Paradoxically, gene and whole chromosome duplications are a principal driver of evolution, while heteromorphic sex chromosomes (XY and ZW) are naturally occurring aneuploidies important for sex determination. Here, we provide an overview of the biology of gene dosage at the crossroads between evolutionary benefit and pathogenicity during disease. We describe the buffering mechanisms and cellular responses to alterations, which could provide a common ground for the understanding of DDs caused by copy number alterations.

Non-genetic and genetic rewiring underlie adaptation to hypomorphic alleles of an essential gene

Adaptive evolution to cellular stress is a process implicated in a wide range of biological and clinical phenomena. Two major routes of adaptation have been identified: non-genetic changes, which allow expression of different phenotypes in novel environments, and genetic variation achieved by selection of fitter phenotypes. While these processes are broadly accepted, their temporal and epistatic features in the context of cellular evolution and emerging drug resistance are contentious. In this manuscript, we generated hypomorphic alleles of the essential nuclear pore complex (NPC) gene NUP58. By dissecting early and long-term mechanisms of adaptation in independent clones, we observed that early physiological adaptation correlated with transcriptome rewiring and upregulation of genes known to interact with the NPC; long-term adaptation and fitness recovery instead occurred via focal amplification of NUP58 and restoration of mutant protein expression. These data support the concept that early phenotypic plasticity allows later acquisition of genetic adaptations to a specific impairment. We propose this approach as a genetic model to mimic targeted drug therapy in human cells and to dissect mechanisms of adaptation.