Where genotype is not predictive of phenotype: towards an understanding of the molecular basis of reduced penetrance in human inherited disease

Novel RB1 germline mutation in a healthy man

BACKGROUND

Retinoblastoma (Rb) most frequently presents as a unilateral sporadic disease up to 40% of cases, however, arise from a monoallelic germline pathogenic variant. Only 10% of the germline mutations are inherited, and high penetrance is seen in up to 90% of these cases. As an effort to optimize counseling and screening, mutations are classified according to inheritance patterns. However, RB1 spectrum is highly heterogeneous, and information for unaffected carriers remains scarce.

MATERIALS AND METHODS

The Mexican family of a 5-month-old patient diagnosed with Rb was studied. The family consisted of five individuals (father, mother, and three siblings). Genetic testing using a next-generation sequencing assay targeting RB1 with oligonucleotide baits designed to capture its exons and 20 bases flanking intronic sequences was performed in every family member. Clinical history and a complete ophthalmological examination (best-corrected visual acuity, slit-lamp biomicroscopy, macular optical coherence tomography, fundus autofluorescence, optical coherence tomography angiography, and electrophysiological studies) were performed in members testing positive to RB1 mutation.

RESULTS

The father and her five-month-old daughter tested positive for a non-synonymous RB1 mutation c.459del (p.Lys154Serfs*21). The girl presented with bilateral retinoblastoma, successfully treated with cryotherapy and intravenous chemotherapy. The father had no relevant findings on imaging studies or ophthalmologic evaluation.

CONCLUSIONS

This report describes a rare case of a novel low-penetrance RB1 germline mutation. Long-term follow-up of the father will include periodic evaluation of the eyes and orbits, and surveillance for systemic sarcoma and secondary malignancies. Implications for unaffected individuals need to be further studied.

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.

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.

Loss of POGZ alters neural differentiation of human embryonic stem cells

POGZ is a pogo transposable element derived protein with multiple zinc finger domains. Many de novo loss-of-function (LoF) variants of the POGZ gene are associated with autism and other neurodevelopmental disorders. However, the role of POGZ in human cortical development remains poorly understood. Here we generated multiple POGZ LoF lines in H9 human embryonic stem cells (hESCs) using CRISPR/CAS9 genome editing. These lines were then differentiated into neural structures, similar to those found in early to mid-fetal human brain, a critical developmental stage for studying disease mechanisms of neurodevelopmental disorders. We found that the loss of POGZ reduced neural stem cell proliferation in excitatory cortex-patterned neural rosettes, structures analogous to the cortical ventricular zone in human fetal brain. As a result, fewer intermediate progenitor cells and early born neurons were generated. In addition, neuronal migration from the apical center to the basal surface of neural rosettes was perturbed due to the loss of POGZ. Furthermore, cortical-like excitatory neurons derived from multiple POGZ homozygous knockout lines exhibited a more simplified dendritic architecture compared to wild type lines. Our findings demonstrate how POGZ regulates early neurodevelopment in the context of human cells, and provide further understanding of the cellular pathogenesis of neurodevelopmental disorders associated with POGZ variants.

Lifetime risk of autosomal recessive neurodegeneration with brain iron accumulation (NBIA) disorders calculated from genetic databases

Background

Neurodegeneration with brain iron accumulation (NBIA) are a group of clinically and genetically heterogeneous diseases characterized by iron overload in basal ganglia and progressive neurodegeneration. Little is known about the epidemiology of NBIA disorders. In the absence of large-scale population-based studies, obtaining reliable epidemiological data requires innovative approaches.

Methods

All pathogenic variants were collected from the 13 genes associated with autosomal recessive NBIA (PLA2G6, PANK2, COASY, ATP13A2, CP, AP4M1, FA2H, CRAT, SCP2, C19orf12, DCAF17, GTPBP2, REPS1). The allele frequencies of these disease-causing variants were assessed in exome/genome collections: the Genome Aggregation Database (gnomAD) and our in-house database. Lifetime risks were calculated from the sum of allele frequencies in the respective genes under assumption of Hardy-Weinberg equilibrium.

Findings

The combined estimated lifetime risk of all 13 investigated NBIA disorders is 0.88 (95% confidence interval 0.70–1.10) per 100,000 based on the global gnomAD dataset (n = 282,912 alleles), 0.92 (0.65–1.29) per 100,000 in the European gnomAD dataset (n = 129,206), and 0.90 (0.48–1.62) per 100,000 in our in-house database (n = 44,324). Individually, the highest lifetime risks (>0.15 per 100,000) are found for disorders caused by variants in PLA2G6, PANK2 and COASY.

Interpretation

This population-genetic estimation on lifetime risks of recessive NBIA disorders reveals frequencies far exceeding previous population-based numbers. Importantly, our approach represents lifetime risks from conception, thus including prenatal deaths. Understanding the true lifetime risk of NBIA disorders is important in estimating disease burden, allocating resources and targeting specific interventions.

Clinical exome sequencing - mistakes and caveats

Massive parallel sequencing technology has become the predominant technique for genetic diagnostics and research. Many genetic laboratories have wrestled with the challenges of setting up genetic testing workflows based on a completely new technology. The learning curve we went through as a laboratory was accompanied by growing pains while we gained new knowledge and expertise. Here we discuss some important mistakes that have been made in our laboratory through 10 years of clinical exome sequencing but that have given us important new insights on how to adapt our working methods. We provide these examples and the lessons that we learned to help other laboratories avoid to make the same mistakes.

Rare disorders have many faces: in silico characterization of rare disorder spectrum

Background

The diagnostic journey for many rare disease patients remains challenging despite use of latest genetic technological advancements. We hypothesize that some patients remain undiagnosed due to more complex diagnostic scenarios that are currently not considered in genome analysis pipelines. To better understand this, we characterized the rare disorder (RD) spectrum using various bioinformatics resources (e.g., Orphanet/Orphadata, Human Phenotype Ontology, Reactome pathways) combined with custom-made R scripts.

Results

Our in silico characterization led to identification of 145 borderline-common, 412 rare and 2967 ultra-rare disorders. Based on these findings and point prevalence, we would expect that approximately 6.53%, 0.34%, and 0.30% of individuals in a randomly selected population have a borderline-common, rare, and ultra-rare disorder, respectively (equaling to 1 RD patient in 14 people). Importantly, our analyses revealed that (1) a higher proportion of borderline-common disorders were caused by multiple gene defects and/or other factors compared with the rare and ultra-rare disorders, (2) the phenotypic expressivity was more variable for the borderline-common disorders than for the rarer disorders, and (3) unique clinical characteristics were observed across the disorder categories forming the spectrum.

Conclusions

Recognizing that RD patients who remain unsolved even after genome sequencing might belong to the more common end of the RD spectrum support the usage of computational pipelines that account for more complex genetic and phenotypic scenarios.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13023-022-02217-9.

Implementation of Early Next-Generation Sequencing for Inborn Errors of Immunity: A Prospective Observational Cohort Study of Diagnostic Yield and Clinical Implications in Dutch Genome Diagnostic Centers

Objective

Inborn errors of immunity (IEI) are a heterogeneous group of disorders, affecting different components of the immune system. Over 450 IEI related genes have been identified, with new genes continually being recognized. This makes the early application of next-generation sequencing (NGS) as a diagnostic method in the evaluation of IEI a promising development. We aimed to provide an overview of the diagnostic yield and time to diagnosis in a cohort of patients suspected of IEI and evaluated by an NGS based IEI panel early in the diagnostic trajectory in a multicenter setting in the Netherlands.

Study Design

We performed a prospective observational cohort study. We collected data of 165 patients with a clinical suspicion of IEI without prior NGS based panel evaluation that were referred for early NGS using a uniform IEI gene panel. The diagnostic yield was assessed in terms of definitive genetic diagnoses, inconclusive diagnoses and patients without abnormalities in the IEI gene panel. We also assessed time to diagnosis and clinical implications.

Results

For children, the median time from first consultation to diagnosis was 119 days versus 124 days for adult patients (U=2323; p=0.644). The median turn-around time (TAT) of genetic testing was 56 days in pediatric patients and 60 days in adult patients (U=1892; p=0.191). A definitive molecular diagnosis was made in 25/65 (24.6%) of pediatric patients and 9/100 (9%) of adults. Most diagnosed disorders were identified in the categories of immune dysregulation (n=10/25; 40%), antibody deficiencies (n=5/25; 20%), and phagocyte diseases (n=5/25; 20%). Inconclusive outcomes were found in 76/165 (46.1%) patients. Within the patient group with a genetic diagnosis, a change in disease management occurred in 76% of patients.

Conclusion

In this cohort, the highest yields of NGS based evaluation for IEI early in the diagnostic trajectory were found in pediatric patients, and in the disease categories immune dysregulation and phagocyte diseases. In cases where a definitive diagnosis was made, this led to important disease management implications in a large majority of patients. More research is needed to establish a uniform diagnostic pathway for cases with inconclusive diagnoses, including variants of unknown significance.

Autism genes converge on asynchronous development of shared neuron classes

Genetic risk for autism spectrum disorder (ASD) is associated with hundreds of genes spanning a wide range of biological functions1-6. The alterations in the human brain resulting from mutations in these genes remain unclear. Furthermore, their phenotypic manifestation varies across individuals7,8. Here we used organoid models of the human cerebral cortex to identify cell-type-specific developmental abnormalities that result from haploinsufficiency in three ASD risk genes-SUV420H1 (also known as KMT5B), ARID1B and CHD8-in multiple cell lines from different donors, using single-cell RNA-sequencing (scRNA-seq) analysis of more than 745,000 cells and proteomic analysis of individual organoids, to identify phenotypic convergence. Each of the three mutations confers asynchronous development of two main cortical neuronal lineages-γ-aminobutyric-acid-releasing (GABAergic) neurons and deep-layer excitatory projection neurons-but acts through largely distinct molecular pathways. Although these phenotypes are consistent across cell lines, their expressivity is influenced by the individual genomic context, in a manner that is dependent on both the risk gene and the developmental defect. Calcium imaging in intact organoids shows that these early-stage developmental changes are followed by abnormal circuit activity. This research uncovers cell-type-specific neurodevelopmental abnormalities that are shared across ASD risk genes and are finely modulated by human genomic context, finding convergence in the neurobiological basis of how different risk genes contribute to ASD pathology.

Retrospective analysis of a clinical exome sequencing cohort reveals the mutational spectrum and identifies candidate disease-associated loci for BAFopathies

PURPOSE

BRG1/BRM-associated factor (BAF) complex is a chromatin remodeling complex that plays a critical role in gene regulation. Defects in the genes encoding BAF subunits lead to BAFopathies, a group of neurodevelopmental disorders with extensive locus and phenotypic heterogeneity.

METHODS

We retrospectively analyzed data from 16,243 patients referred for clinical exome sequencing (ES) with a focus on the BAF complex. We applied a genotype-first approach, combining predicted genic constraints to propose candidate BAFopathy genes.

RESULTS

We identified 127 patients carrying pathogenic variants, likely pathogenic variants, or de novo variants of unknown clinical significance in 11 known BAFopathy genes. Those include 34 patients molecularly diagnosed using ES reanalysis with new gene-disease evidence (n = 21) or variant reclassifications in known BAFopathy genes (n = 13). We also identified de novo or predicted loss-of-function variants in 4 candidate BAFopathy genes, including ACTL6A, BICRA (implicated in Coffin-Siris syndrome during this study), PBRM1, and SMARCC1.

CONCLUSION

We report the mutational spectrum of BAFopathies in an ES cohort. A genotype-driven and pathway-based reanalysis of ES data identified new evidence for candidate genes involved in BAFopathies. Further mechanistic and phenotypic characterization of additional patients are warranted to confirm their roles in human disease and to delineate their associated phenotypic spectrums.

Trypsinogen (PRSS1 and PRSS2) gene dosage correlates with pancreatitis risk across genetic and transgenic studies: a systematic review and re-analysis

Trypsinogen (PRSS1, PRSS2) copy number gains and regulatory variants have both been proposed to elevate pancreatitis risk through a gene dosage effect (i.e., by increasing the expression of wild-type protein). However, to date, their impact on pancreatitis risk has not been thoroughly evaluated whilst the underlying pathogenic mechanisms remain to be explicitly investigated in mouse models. Genetic studies of the rare trypsinogen duplication and triplication copy number variants (CNVs), and the common rs10273639C variant, were collated from PubMed and/or ClinVar. Mouse studies that analyzed the influence of a transgenically expressed wild-type human PRSS1 or PRSS2 gene on the development of pancreatitis were identified from PubMed. The genetic effects of the different risk genotypes, in terms of odds ratios, were calculated wherever appropriate. The genetic effects of the rare trypsinogen duplication and triplication CNVs were also evaluated by reference to their associated disease subtypes. We demonstrate a positive correlation between increased trypsinogen gene dosage and pancreatitis risk in the context of the rare duplication and triplication CNVs, and between the level of trypsinogen expression and disease risk in the context of the heterozygous and homozygous rs10273639C-tagged genotypes. We retrospectively identify three mouse transgenic studies that are informative in relation to the pathogenic mechanism underlying the trypsinogen gene dosage effect in pancreatitis. Trypsinogen gene dosage correlates with pancreatitis risk across genetic and transgenic studies, highlighting the fundamental role of dysregulated expression of wild-type trypsinogen in the etiology of pancreatitis. Specifically downregulating trypsinogen expression in the pancreas may serve as a potential therapeutic and/or prevention strategy for pancreatitis.

Computational Genomics in the Era of Precision Medicine: Applications to Variant Analysis and Gene Therapy

Rapid methodological advances in statistical and computational genomics have enabled researchers to better identify and interpret both rare and common variants responsible for complex human diseases. As we continue to see an expansion of these advances in the field, it is now imperative for researchers to understand the resources and methodologies available for various data types and study designs. In this review, we provide an overview of recent methods for identifying rare and common variants and understanding their roles in disease etiology. Additionally, we discuss the strategy, challenge, and promise of gene therapy. As computational and statistical approaches continue to improve, we will have an opportunity to translate human genetic findings into personalized health care.

Population-Based Penetrance of Deleterious Clinical Variants

IMPORTANCE

Population-based assessment of disease risk associated with gene variants informs clinical decisions and risk stratification approaches.

OBJECTIVE

To evaluate the population-based disease risk of clinical variants in known disease predisposition genes.

DESIGN, SETTING, AND PARTICIPANTS

This cohort study included 72 434 individuals with 37 780 clinical variants who were enrolled in the BioMe Biobank from 2007 onwards with follow-up until December 2020 and the UK Biobank from 2006 to 2010 with follow-up until June 2020. Participants had linked exome and electronic health record data, were older than 20 years, and were of diverse ancestral backgrounds.

EXPOSURES

Variants previously reported as pathogenic or predicted to cause a loss of protein function by bioinformatic algorithms (pathogenic/loss-of-function variants).

MAIN OUTCOMES AND MEASURES

The primary outcome was the disease risk associated with clinical variants. The risk difference (RD) between the prevalence of disease in individuals with a variant allele (penetrance) vs in individuals with a normal allele was measured.

RESULTS

Among 72 434 study participants, 43 395 were from the UK Biobank (mean [SD] age, 57 [8.0] years; 24 065 [55%] women; 2948 [7%] non-European) and 29 039 were from the BioMe Biobank (mean [SD] age, 56 [16] years; 17 355 [60%] women; 19 663 [68%] non-European). Of 5360 pathogenic/loss-of-function variants, 4795 (89%) were associated with an RD less than or equal to 0.05. Mean penetrance was 6.9% (95% CI, 6.0%-7.8%) for pathogenic variants and 0.85% (95% CI, 0.76%-0.95%) for benign variants reported in ClinVar (difference, 6.0 [95% CI, 5.6-6.4] percentage points), with a median of 0% for both groups due to large numbers of nonpenetrant variants. Penetrance of pathogenic/loss-of-function variants for late-onset diseases was modified by age: mean penetrance was 10.3% (95% CI, 9.0%-11.6%) in individuals 70 years or older and 8.5% (95% CI, 7.9%-9.1%) in individuals 20 years or older (difference, 1.8 [95% CI, 0.40-3.3] percentage points). Penetrance of pathogenic/loss-of-function variants was heterogeneous even in known disease predisposition genes, including BRCA1 (mean [range], 38% [0%-100%]), BRCA2 (mean [range], 38% [0%-100%]), and PALB2 (mean [range], 26% [0%-100%]).

CONCLUSIONS AND RELEVANCE

In 2 large biobank cohorts, the estimated penetrance of pathogenic/loss-of-function variants was variable but generally low. Further research of population-based penetrance is needed to refine variant interpretation and clinical evaluation of individuals with these variant alleles.

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.

Analysis of the Frequency of 10 Polymorphic Markers of CDKN2A and RB1 Genes in Russian Populations

The study of population frequencies of rare clinically significant alleles is a prerequisite of the development of personalized medicine. We performed genotyping of 1785 DNA samples from representatives of Russian populations according to 10 benign polymorphic markers of two genes involved in oncogenesis: 3 variants of the CDKN2A gene (rs3731249, rs116150891, and rs6413464) and 7 markers of the RB1 gene (rs149800437, rs147754935, rs183898408, rs146897002, rs4151539, rs187912365, and rs144668210). Genotyping was performed using the Illumina biochip test system. The sample covered 28 populations of the Russian Federation and neighboring countries, which were later combined into 3 groups (Asian, European, and Caucasian). The information from the ALFA (NCBI) project was used as reference for the frequencies of these polymorphisms in the Asian and European populations. It was shown that rare alleles in 8 of 10 studied polymorphic markers are presented in Russian populations of European and Caucasian origin with frequencies that are tens and hundreds of times higher than the available data for Western European populations, and in Russian Asian populations, alternative alleles of 5 markers absent in the Asian population of the ALFA project were found. In the subpopulation of Astrakhan Tatars, exceptionally high frequencies of rare alleles were identified; this requires further study.

Unraveling the genotype by environment interaction in a thermosensitive fish with a polygenic sex determination system

In most animals, sex determination occurs at conception, when sex chromosomes are segregated following Mendelian laws. However, in multiple reptiles and fishes, this genetic sex can be overridden by external factors after fertilization or birth. In some species, the genetic sex may also be governed by multiple genes, further limiting our understanding of sex determination in such species. We used the European sea bass (Dicentrarchus labrax) as a model and combined genomic (using a single nucleotide polymorphism chip) and transcriptomic (RNA-Sequencing) approaches to thoroughly depict this polygenic sex determination system and its interaction with temperature. We estimated genetic sex tendency (eGST), defined as the estimated genetic liability to become a given sex under a liability threshold model for sex determination, which accurately predicts the future phenotypic sex. We found evidence that energetic pathways, concerning the regulation of lipids and glucose, are involved in sex determination and could explain why females tend to exhibit higher energy levels and improved growth compared to males. Besides, early exposure to high-temperature up-regulated sox3, followed by sox9a in individuals with intermediate eGST, but not in individuals showing highly female-biased eGST, providing the most parsimonious explanation for temperature-induced masculinization. This gonadal state was maintained likely by DNA methylation and the up-regulation of several genes involved in histone modifications, including jmjd1c Overall, we describe a sex determination system resulting from continuous genetic and environmental influences in an animal. Our results provide significant progress in our understanding of the mechanisms underlying temperature-induced masculinization in fish.

Gastric cancer genetic predisposition and clinical presentations: Established heritable causes and potential candidate genes

Tumour risk syndromes (TRS) are characterized by an increased risk of early-onset cancers in a familial context. High cancer risk is mostly driven by loss-of-function variants in a single cancer-associated gene. Presently, predisposition to diffuse gastric cancer (DGC) is explained by CDH1 and CTNNA1 pathogenic and likely pathogenic variants (P/LP), causing Hereditary Diffuse Gastric Cancer (HDGC); while APC promoter 1B single nucleotide variants predispose to Gastric Adenocarcinoma and Proximal Polyposis of the Stomach (GAPPS). Familial Intestinal Gastric Cancer (FIGC), recognized as a GC-predisposing disease, remains understudied and genetically unsolved. GC can also occur in the spectrum of other TRS. Identification of heritable causes allows defining diagnostic testing criteria, helps to clinically classify GC families into the appropriate TRS, and allows performing pre-symptomatic testing identifying at-risk individuals for downstream surveillance, risk reduction and/or treatment. However, most of HDGC, some GAPPS, and most FIGC patients/families remain unsolved, expecting a heritable factor to be discovered. The missing heritability in GC-associated tumour risk syndromes (GC-TRS) is likely explained not by a single major gene, but by a diversity of genes, some, predisposing to other TRS. This would gain support if GC-enriched small families or apparently isolated early-onset GC cases were hiding a family history compatible with another TRS. Herein, we revisited current knowledge on GC-TRS, and searched in the literature for individuals/families bearing P/LP variants predisposing for other TRS, but whose probands display a clinical presentation and/or family history also fitting GC-TRS criteria. We found 27 families with family history compatible with HDGC or FIGC, harbouring 28 P/LP variants in 16 TRS-associated genes, mainly associated with DNA repair. PALB2 or BRCA2 were the most frequently mutated candidate genes in individuals with family history compatible with HDGC and FIGC, respectively. Consolidation of PALB2 and BRCA2 as HDGC- or FIGC-associated genes, respectively, holds promise and worth additional research. This analysis further highlighted the influence, that proband's choice and small or unreported family history have, for a correct TRS diagnosis, genetic screening, and disease management. In this review, we provide a rational for identification of particularly relevant candidate genes in GC-TRS.

Bilirubin: The yellow hormone?

Bilirubin is a tetrapyrrolic compound originating from heme catabolism. Although originally considered only a potentially dangerous waste product, it has become increasingly evident that this molecule represents an important modulator of various biological functions in the human body. Bilirubin appears to have versatile functions, from cell signaling (behaving almost like a "real" hormonal substance), modulation of metabolism, to immune regulation, affecting biological activities with apparent clinical and even therapeutic consequences. These activities may be the reason for the lower incidence of diseases of civilisation (cardiovascular diseases, arterial hypertension, diabetes, obesity, metabolic syndrome, certain cancers, autoimmune, and neurodegenerative diseases) observed in individuals with a chronic mild unconjugated hyperbilirubinemia, a typical sign of Gilbert's syndrome. While higher serum concentrations of unconjugated bilirubin may serve as an important protective factor against these diseases, low levels of bilirubin are associated with the opposite effect.

Penetrance estimation of PRRT2 variants in paroxysmal kinesigenic dyskinesia and infantile convulsions

Proline-rich transmembrane protein 2 (PRRT2) is the leading cause of paroxysmal kinesigenic dyskinesia (PKD), benign familial infantile epilepsy (BFIE), and infantile convulsions with choreoathetosis (ICCA). Reduced penetrance of PRRT2 has been observed in previous studies, whereas the exact penetrance has not been evaluated well. The objective of this study was to estimate the penetrance of PRRT2 and determine its influencing factors. We screened 222 PKD index patients and their available relatives, identified 39 families with pathogenic or likely pathogenic (P/LP) PRRT2 variants via Sanger sequencing, and obtained 184 PKD/BFIE/ICCA families with P/LP PRRT2 variants from the literature. Penetrance was estimated as the proportion of affected variant carriers. PRRT2 penetrance estimate was 77.6% (95% confidence interval (CI) 74.5%-80.7%) in relatives and 74.5% (95% CI 70.2%-78.8%) in obligate carriers. In addition, we first observed that penetrance was higher in truncated than in non-truncated variants (75.8% versus 50.0%, P = 0.01), higher in Asian than in Caucasian carriers (81.5% versus 68.5%, P = 0.004), and exhibited no difference in gender or parental transmission. Our results are meaningful for genetic counseling, implying that approximately three-quarters of PRRT2 variant carriers will develop PRRT2-related disorders, with patients from Asia or carrying truncated variants at a higher risk.

The Genetics of Inherited Cholestatic Disorders in Neonates and Infants: Evolving Challenges

Many inherited conditions cause cholestasis in the neonate or infant. Next-generation sequencing methods can facilitate a prompt diagnosis in some of these cases; application of these methods in patients with liver diseases of unknown cause has also uncovered novel gene-disease associations and improved our understanding of physiological bile secretion and flow. By helping to define the molecular basis of certain cholestatic disorders, these methods have also identified new targets for therapy as well patient subgroups more likely to benefit from specific therapies. At the same time, sequencing methods have presented new diagnostic challenges, such as the interpretation of single heterozygous genetic variants. This article discusses those challenges in the context of neonatal and infantile cholestasis, focusing on difficulties in predicting variant pathogenicity, the possibility of other causal variants not identified by the genetic screen used, and phenotypic variability among patients with variants in the same genes. A prospective, observational study performed between 2010-2013, which sequenced six important genes (ATP8B1, ABCB11, ABCB4, NPC1, NPC2 and SLC25A13) in an international cohort of 222 patients with infantile liver disease, is given as an example of potential benefits and challenges that clinicians could face having received a complex genetic result. Further studies including large cohorts of patients with paediatric liver disease are needed to clarify the spectrum of phenotypes associated with, as well as appropriate clinical response to, single heterozygous variants in cholestasis-associated genes.

Genetic and Bioinformatic Strategies to Improve Diagnosis in Three Inherited Bleeding Disorders in Bogotá, Colombia

Inherited bleeding disorders (IBDs) are the most frequent congenital diseases in the Colombian population; three of them are hemophilia A (HA), hemophilia B (HB), and von Willebrand Disease (VWD). Currently, diagnosis relies on multiple clinical laboratory assays to assign a phenotype. Due to the lack of accessibility to these tests, patients can receive an incomplete diagnosis. In these cases, genetic studies reinforce the clinical diagnosis. The present study characterized the molecular genetic basis of 11 HA, three HB, and five VWD patients by sequencing the F8, F9, or the VWF gene. Twelve variations were found in HA patients, four in HB patients, and 19 in WVD patients. From these variations a total of 25 novel variations were found. Disease-causing variations were used as positive controls for validation of the high-resolution melting (HRM) variant-scanning technique. This approach is a low-cost genetic diagnostic method proposed to be incorporated in developing countries. For the data analysis, we developed an accessible open-source code in Python that improves HRM data analysis with better sensitivity of 95% and without bias when using different HRM equipment and software. Analysis of amplicons with a length greater than 300 bp can be performed by implementing an analysis by denaturation domains.

A Protocol for Preconceptional Screening of Consanguineous Couples Using Whole Exome Sequencing

Genetic studies performed in consanguineous couples suggest that the reproductive risk that distinguish them from other couples in the general population is related to autosomal recessive (AR) diseases. This risk is scattered among the thousands of known and potential AR diseases. Thus, for effective preconceptional screening of consanguineous couples it is necessary a test that encompasses the largest number of genes possible. For that reason, we decided to create a protocol based on whole exome sequencing (WES). We sequenced completely the exomes of 39 consanguineous couples at high coverage (∼100×). Applying bioinformatics filters, we could detect genetic variants that were simultaneously present in both members of the couple in all genes listed in the Clinical Genomics Database as causally related to AR diseases. Shared variants were then assessed for pathogenicity. For non-truncating variants (missense and in-frame indels) we considered as pathogenic or likely pathogenic only the variants included as such in the ClinVar database. Shared truncating variants (frameshift, non-sense, and canonical splice variants) were considered likely pathogenic when loss-of-function was a known mechanism of disease. The 39 consanguineous cases included two couples with a coefficient of genetic relationship (CGR) of 0.25, 26 couples with a CGR of 0.125, three couples with a CGR of 0.0625 and eight couples with a CGR of 0.03125. In 21 of the 39 couples (53.8%) we ascertained sharing of heterozygosity for at least one variant considered pathogenic or likely pathogenic for an AR disease. In eight couples we found sharing of heterozygosity for at least two pathogenic variants. Once the specific pathogenic variant was identified, it became possible for the couple to undergo prenatal diagnosis or, if desired, preimplantation genetic diagnosis (PGD) involving in vitro fertilization and embryo screening. In conclusion, our results demonstrate that preconceptional screening by WES is a useful new procedure that should be incorporated in the genetic counseling of all consanguineous couples.

The Italian National Registry for FSHD: an enhanced data integration and an analytics framework towards Smart Health Care and Precision Medicine for a rare disease

Background

The Italian Clinical network for FSHD (ICNF) has established the Italian National Registry for FSHD (INRF), collecting data from patients affected by Facioscapulohumeral dystrophy (FSHD) and their relatives. The INRF has gathered data from molecular analysis, clinical evaluation, anamnestic information, and family history from more than 3500 participants.

Methods

A data management framework, called Mediator Environment for Multiple Information Sources (MOMIS) FSHD Web Platform, has been developed to provide charts, maps and search tools customized for specific needs. Patients’ samples and their clinical information derives from the Italian Clinical network for FSHD (ICNF), a consortium consisting of fourteen neuromuscular clinics distributed across Italy. The tools used to collect, integrate, and visualize clinical, molecular and natural history information about patients affected by FSHD and their relatives are described.

Results

The INRF collected the molecular data regarding FSHD diagnosis conducted on 7197 subjects and identified 3362 individuals carrying a D4Z4 Reduced Allele (DRA): 1634 were unrelated index cases. In 1032 cases the molecular testing has been extended to 3747 relatives, 1728 carrying a DRA. Since 2009 molecular analysis has been accompanied by clinical evaluation based standardized evaluation protocols. In the period 2009–2020, 3577 clinical forms have been collected, 2059 follow the Comprehensive Clinical Evaluation form (CCEF). The integration of standardized clinical information and molecular data has made possible to demonstrate the wide phenotypic variability of FSHD. The MOMIS (Mediator Environment for Multiple Information Sources) data integration framework allowed performing genotype–phenotype correlation studies, and generated information of medical importance either for clinical practice or genetic counseling.

Conclusion

The platform implemented for the FSHD Registry data collection based on OpenClinica meets the requirement to integrate patient/disease information, as well as the need to adapt dynamically to security and privacy concerns. Our results indicate that the quality of data collection in a multi-integrated approach is fundamental for clinical and epidemiological research in a rare disease and may have great value in allowing us to redefine diagnostic criteria and disease markers for FSHD. By extending the use of the MOMIS data integration framework to other countries and the longitudinal systematic collection of standardized clinical data will facilitate the understanding of disease natural history and offer valuable inputs towards trial readiness. This approach is of high significance to FSHD medical community and also to rare disease research in general.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13023-021-02100-z.

Targeted Sequencing of 242 Clinically Important Genes in the Russian Population From the Ivanovo Region

We performed a targeted sequencing of 242 clinically important genes mostly associated with cardiovascular diseases in a representative population sample of 1,658 individuals from the Ivanovo region northeast of Moscow. Approximately 11% of 11,876 detected variants were not found in the Single Nucleotide Polymorphism Database (dbSNP) or reported earlier in the Russian population. Most novel variants were singletons and doubletons in our sample, and virtually no novel alleles presumably specific for the Russian population were able to reach the frequencies above 0.1-0.2%. The overwhelming majority (99.3%) of variants detected in this study in three or more copies were shared with other populations. We found two dominant and seven recessive known pathogenic variants with allele frequencies significantly increased compared to those in the gnomAD non-Finnish Europeans. Of the 242 targeted genes, 28 were in the list of 59 genes for which the American College of Medical Genetics and Genomics (ACMG) recommended the reporting of incidental findings. Based on the number of variants detected in the sequenced subset of ACMG59 genes, we approximated the prevalence of known pathogenic and novel or rare protein-truncating variants in the complete set of ACMG59 genes in the Ivanovo population at 1.4 and 2.8%, respectively. We analyzed the available clinical data and observed the incomplete penetrance of known pathogenic variants in the 28 ACMG59 genes: only 1 individual out of 12 with such variants had the phenotype most likely related to the variant. When known pathogenic and novel or rare protein-truncating variants were considered together, the overall rate of confirmed phenotypes was about 19%, with maximum in the subset of novel protein-truncating variants. We report three novel protein truncating variants in APOB and one in MYH7 observed in individuals with hypobetalipoproteinemia and hypertrophic cardiomyopathy, respectively. Our results provide a valuable reference for the clinical interpretation of gene sequencing in Russian and other populations.

AML with germline DDX41 variants is a clinicopathologically distinct entity with an indolent clinical course and favorable outcome

Germline DDX41 variants in myeloid neoplasms (MNs) are not uncommon, and we explored the prevalence and characterized the clinical and pathologic features in a cohort of 3132 unrelated adult MN patients. By targeted next-generation sequencing, we identified 28 patients (20 men and 8 women) with pathogenic germline DDX41 variants who developed acute myeloid leukemia (AML), in which only 3 (11%) had a family history (FH) of MNs. A subacute clinical course of cytopenia (mean duration of 11.2 months, range 0-72 months) prior to the initial AML diagnosis was accompanied by a low blast count (median at 30%, range 20-70%) in hypocellular marrows (93% of all patients), in vast contrast to the typical proliferative subtypes of AML in the elderly. Most patients had a normal karyotype (75%) and acquired a second DDX41 variant (69%). A favorable overall survival (OS) was observed in comparison to that of common subtypes of AML with wild-type DDX41 in age-matched patients. Our study demonstrated that the frequent germline pathogenic DDX41 variants characterized a clinically distinct AML entity. Features characteristic of DDX41-mutated AML include male predominance, often lack of FH, indolent course, low proliferative potential, frequent somatic DDX41 variants, and a favorable OS.

Parental Origin of the RB1 Gene Mutations in Families with Low Penetrance Hereditary Retinoblastoma

Simple Summary

Some families with hereditary retinoblastoma exhibit mild phenotype with low penetrance and variable expressivity, including complete absence of clinical signs of the disease in some carriers of the germline RB1 mutation. The identification of low-penetrance mutations in the RB1 gene and the study of their inheritance in pedigrees is contributing to understanding the mechanisms underlying the development of retinoblastoma with low penetrance. It is important both for further expansion of knowledge in the field of molecular genetics of retinoblastoma, and for competent genetic counseling and subsequent clinical management of families with this form of the disease. Our results support an assumption that parental origin of an RB1 mutation influences the likelihood of developing retinoblastoma. We also revealed a relatively high frequency of asymptomatic carriage of the RB1 mutations among the parents of retinoblastoma patients, highlighting the utmost necessity for molecular analysis among the probands’ relatives irrespective of their clinical status and family history of retinoblastoma.

Abstract

Our aim was to identify RB1 alterations causing hereditary low penetrance retinoblastoma and to evaluate how the parental origin of an RB1 mutation affects its phenotypic expression. By NGS and MLPA, RB1 mutations were found in 191 from 332 unrelated retinoblastoma patients. Among patients with identified RB1 mutations but without clinical family history of retinoblastoma, 7% (12/175) were found to have hereditary disease with one of the parents being an asymptomatic carrier of an RB1 mutation. Additionally, in two families with retinoblastoma history, mutations were inherited by probands from unaffected parents. Overall, nine probands inherited RB1 mutations from clinically unaffected fathers and five, from mothers. Yet, we gained explanations of maternal “unaffectedness” in most cases, either as somatic mosaicism or as clinical presentation of retinomas in involution, rendering the proportion of paternal to maternal truly asymptomatic mutation carriers as 9:1 (p = 0.005). This observation supports an assumption that parental origin of an RB1 mutation influences the likelihood of developing retinoblastoma. Additionally, our study revealed a relatively high frequency of asymptomatic carriage of the RB1 mutations among the parents of retinoblastoma patients, highlighting the utmost necessity of molecular analysis among the probands’ relatives irrespective of their clinical status and family history of retinoblastoma.

Heritability of Locomotor Stereotypies in Chilean Horses

The Chilean horse is a breed of closed registry. Stall-walking and weaving are locomotor stereotypies that affect this breed, and genetic predisposition has been suggested for both conditions. The objective of the present study was to estimate heritability of stall-walking and weaving in Chilean horses. Owners of 2,098 horses registered in the Chilean horse Stud Book, which were or had been stabled for at least 1 year, were asked to provide for identification data of the animal and presence or absence of stall-walking and/or weaving. The Chilean Horse Stud Book was accessed online, to collect information on name and registration number of the sire and dam of each horse. The prevalence of stall-walking and weaving was calculated from the sample (n = 2,098). The database to estimate the heritability (h²) included all the sample horses (n = 2,098) and their sires and dams. Additionally, all the ancestors available in the Chilean Horse Stud Book were included in the database for 297 sampled horses including those with locomotor stereotypies. The genealogical database consisted of 7,187 individuals. The prevalence of stall-walking was 2.05% and prevalence of weaving was 1.43%, being more frequent in males (P < .05) for stall-walking. Heritability of stall-walking was low (h² = 0.213 ± 0.08) and moderate for weaving (h² = 0.435 ± 0.06). Heritability and prevalence of locomotor stereotypies found in the Chilean horse suggest that they are hereditary disorders, highly influenced by environmental factors.

Integration of multiple imaging platforms to uncover cardiovascular defects in adult zebrafish

Aims

Mammalian models have been instrumental in investigating adult heart function and human disease. However, electrophysiological differences with human hearts and high costs motivate the need for non-mammalian models. The zebrafish is a well-established genetic model to study cardiovascular development and function; however, analysis of cardiovascular phenotypes in adult specimens is particularly challenging as they are opaque.

Methods and results

Here, we optimized and combined multiple imaging techniques including echocardiography, magnetic resonance imaging, and micro-computed tomography to identify and analyse cardiovascular phenotypes in adult zebrafish. Using alk5a/tgfbr1a mutants as a case study, we observed morphological and functional cardiovascular defects that were undetected with conventional approaches. Correlation analysis of multiple parameters revealed an association between haemodynamic defects and structural alterations of the heart, as observed clinically.

Conclusion

We report a new, comprehensive, and sensitive platform to identify otherwise indiscernible cardiovascular phenotypes in adult zebrafish.

Decoding 'Unnecessary Complexity': A Law of Complexity and a Concept of Hidden Variation Behind "Missing Heritability" in Precision Medicine

The high hopes for the Human Genome Project and personalized medicine were not met because the relationship between genotypes and phenotypes turned out to be more complex than expected. In a previous study we laid the foundation of a theory of complexity and showed that because of the blind nature of evolution, and molecular and historical contingency, cells have accumulated unnecessary complexity, complexity beyond what is necessary and sufficient to describe an organism. Here we provide empirical evidence and show that unnecessary complexity has become integrated into the genome in the form of redundancy and is relevant to molecular evolution of phenotypic complexity. Unnecessary complexity creates uncertainty between molecular and phenotypic complexity, such that phenotypic complexity (CP) is higher than molecular complexity (CM), which is higher than DNA complexity (CD). The qualitative inequality in complexity is based on the following hierarchy: CP > CM > CD. This law-like relationship holds true for all complex traits, including complex diseases. We present a hypothesis of two types of variation, namely open and closed (hidden) systems, show that hidden variation provides a hitherto undiscovered "third source" of phenotypic variation, beside genotype and environment, and argue that "missing heritability" for some complex diseases is likely to be a case of "diluted heritability". There is a need for radically new ways of thinking about the principles of genotype-phenotype relationship. Understanding how cells use hidden, pathway variation to respond to stress can shed light on why two individuals who share the same risk factors may not develop the same disease, or how cancer cells escape death.

Dynapaenia and sarcopaenia in chronic haemodialysis patients: do muscle weakness and atrophy similarly influence poor outcome?

BACKGROUND

Sarcopaenia, defined as a decline in both muscle mass and function, has been recognized as a major determinant of poor outcome in haemodialysis (HD) patients. It is generally assumed that sarcopaenia is driven by muscle atrophy related to protein-energy wasting. However, dynapaenia, defined as weakness without atrophy, has been characterized by a different disease phenotype from sarcopaenia. The aim of this study was to compare the characteristics and prognosis of sarcopaenic and dynapaenic patients among a prospective cohort of chronic HD (CHD) patients.

METHODS

Two hundred and thirty-two CHD patients were enrolled from January to July 2016 and then followed prospectively until December 2018. At inclusion, weakness and atrophy were, respectively, evaluated by maximal voluntary force (MVF) and creatinine index (CI). Sarcopaenia was defined as the association of weakness and atrophy (MVF and CI below the median) while dynapaenia was defined as weakness not related to atrophy (MVF below the median, and CI above the median).

RESULTS

From a total of 187 prevalent CHD patients [65% of men, age 65.3 (49.7-82.0) years], 44 died during the follow-up period of 23.7 (12.4-34.9) months. Sarcopaenia and dynapaenia were observed in 33.7 and 16% of the patients, respectively. Compared with patients with sarcopaenia, patients with dynapaenia were younger and with a lower Charlson score. In contrast, mortality rate was similar in both groups (38 and 27%, respectively). After adjustment for age, sex, lean tissue index, serum albumin, high-sensitivity C-reactive protein (hs-CRP), haemoglobin (Hb), normalized protein catabolic rate (nPCR), dialysis vintage and Charlson score, only patients with dynapaenia were at increased risk of death [hazard ratio (HR) = 2.99, confidence interval 1.18-7.61; P = 0.02].

CONCLUSIONS

Screening for muscle functionality is highly warranted to identify patients with muscle functional impairment without muscle atrophy. In contrast to sarcopaenia, dynapaenia should appear as a phenotype induced by uraemic milieu, characterized by young patients with low Charlson score and poor prognosis outcome independently of serum albumin, hs-CRP, Hb, nPCR and dialysis vintage.

Protocadherin 19 Clustering Epilepsy and Neurosteroids: Opportunities for Intervention

Steroids yield great influence on neurological development through nuclear hormone receptor (NHR)-mediated gene regulation. We recently reported that cell adhesion molecule protocadherin 19 (encoded by the PCDH19 gene) is involved in the coregulation of steroid receptor activity on gene expression. PCDH19 variants cause early-onset developmental epileptic encephalopathy clustering epilepsy (CE), with altered steroidogenesis and NHR-related gene expression being identified in these individuals. The implication of hormonal pathways in CE pathogenesis has led to the investigation of various steroid-based antiepileptic drugs in the treatment of this disorder, with mixed results so far. Therefore, there are many unmet challenges in assessing the antiseizure targets and efficiency of steroid-based therapeutics for CE. We review and assess the evidence for and against the implication of neurosteroids in the pathogenesis of CE and in view of their possible clinical benefit.

The genetic structure of the Turkish population reveals high levels of variation and admixture

The construction of population-based variomes has contributed substantially to our understanding of the genetic basis of human inherited disease. Here, we investigated the genetic structure of Turkey from 3,362 unrelated subjects whose whole exomes (n = 2,589) or whole genomes (n = 773) were sequenced to generate a Turkish (TR) Variome that should serve to facilitate disease gene discovery in Turkey. Consistent with the history of present-day Turkey as a crossroads between Europe and Asia, we found extensive admixture between Balkan, Caucasus, Middle Eastern, and European populations with a closer genetic relationship of the TR population to Europeans than hitherto appreciated. We determined that 50% of TR individuals had high inbreeding coefficients (≥0.0156) with runs of homozygosity longer than 4 Mb being found exclusively in the TR population when compared to 1000 Genomes Project populations. We also found that 28% of exome and 49% of genome variants in the very rare range (allele frequency < 0.005) are unique to the modern TR population. We annotated these variants based on their functional consequences to establish a TR Variome containing alleles of potential medical relevance, a repository of homozygous loss-of-function variants and a TR reference panel for genotype imputation using high-quality haplotypes, to facilitate genome-wide association studies. In addition to providing information on the genetic structure of the modern TR population, these data provide an invaluable resource for future studies to identify variants that are associated with specific phenotypes as well as establishing the phenotypic consequences of mutations in specific genes.

Augmenting and directing long-range CRISPR-mediated activation in human cells

Epigenetic editing is an emerging technology that uses artificial transcription factors (aTFs) to regulate expression of a target gene. Although human genes can be robustly upregulated by targeting aTFs to promoters, the activation induced by directing aTFs to distal transcriptional enhancers is substantially less robust and consistent. Here we show that long-range activation using CRISPR-based aTFs in human cells can be made more efficient and reliable by concurrently targeting an aTF to the target gene promoter. We used this strategy to direct target gene choice for enhancers capable of regulating more than one promoter and to achieve allele-selective activation of human genes by targeting aTFs to single-nucleotide polymorphisms embedded in distally located sequences. Our results broaden the potential applications of the epigenetic editing toolbox for research and therapeutics.

Homologous Recombination Deficiency: Cancer Predispositions and Treatment Implications

Homologous recombination (HR) is a highly accurate DNA repair mechanism. Several HR genes are established cancer susceptibility genes with clinically actionable pathogenic variants (PVs). Classically, BRCA1 and BRCA2 germline PVs are associated with significant breast and ovarian cancer risks. Patients with BRCA1 or BRCA2 PVs display worse clinical outcomes but respond better to platinum-based chemotherapies and poly-ADP ribose polymerase inhibitors, a trait termed "BRCAness." With the advent of whole-exome sequencing and multigene panels, PVs in other HR genes are increasingly identified among familial cancers. As such, several genes such as PALB2 are reclassified as cancer predisposition genes. But evidence for cancer risks remains unclear for many others. In this review, we will discuss cancer predispositions and treatment implications beyond BRCA1 and BRCA2, with a focus on 24 HR genes: 53BP1, ATM, ATR, ATRIP, BARD1, BLM, BRIP1, DMC1, MRE11A, NBN, PALB2, RAD50, RAD51, RAD51B, RAD51C, RAD51D, RIF1, RMI1, RMI2, RPA1, TOP3A, TOPBP1, XRCC2, and XRCC3. IMPLICATIONS FOR PRACTICE: This review provides a comprehensive reference for readers to quickly identify potential cancer predisposing homologous recombination (HR) genes, and to generate research questions for genes with inconclusive evidence. This review also evaluates the "BRCAness" of each HR member. Clinicians can refer to these discussions to identify potential candidates for future clinical trials.

Variants in the ethylmalonyl-CoA decarboxylase (ECHDC1) gene: a novel player in ethylmalonic aciduria?

Ethylmalonic acid (EMA) is a major and potentially cytotoxic metabolite associated with short-chain acyl-CoA dehydrogenase (SCAD) deficiency, a condition whose status as a disease is uncertain. Unexplained high EMA is observed in some individuals with complex neurological symptoms, who carry the SCAD gene (ACADS) variants, c.625G>A and c.511C>T. The variants have a high allele frequency in the general population, but are significantly overrepresented in individuals with elevated EMA. This has led to the idea that these variants need to be associated with variants in other genes to cause hyperexcretion of ethylmalonic acid and possibly a diseased state. Ethylmalonyl-CoA decarboxylase (ECHDC1) has been described and characterized as an EMA metabolite repair enzyme, however, its clinical relevance has never been investigated. In this study, we sequenced the ECHDC1 gene (ECHDC1) in 82 individuals, who were reported with unexplained high EMA levels due to the presence of the common ACADS variants only. Three individuals with ACADS c.625G>A variants were found to be heterozygous for ECHDC1 loss-of-function variants. Knockdown experiments of ECHDC1, in healthy human cells with different ACADS c.625G>A genotypes, showed that ECHDC1 haploinsufficiency and homozygosity for the ACADS c.625G>A variant had a synergistic effect on cellular EMA excretion. This study reports the first cases of ECHDC1 gene defects in humans and suggests that ECHDC1 may be involved in elevated EMA excretion in only a small group of individuals with the common ACADS variants. However, a direct link between ECHDC1/ACADS deficiency, EMA and disease could not be proven.

Single-cell multimodal analysis in a case with reduced penetrance of Progranulin-Frontotemporal Dementia

We identified an autosomal dominant progranulin mutation carrier without symptoms of dementia in her lifetime (Reduced Penetrance Mutation Carrier, RedPenMC). This resistance to develop expected pathology presents a unique opportunity to interrogate neurodegenerative mechanisms. We performed multimodal single-nuclei analyses of post-mortem frontal cortex from RedPenMC, including transcriptomics and global levels of chromatin marks. RedPenMC had an increased ratio of GRN-expressing microglia, higher levels of activating histone mark H3k4me3 in microglia and lower levels of the repressive chromatin marks H3k9me1 and H3k9me3 in the frontal cortex than her affected mutation carrier son and evidence of higher protein levels of progranulin in both plasma and brain homogenates. Although the study is limited to one case, the results support that restoring brain progranulin levels may be sufficient to escape neurodegeneration and FTD. In addition to previously identified modifier genes, it is possible that epigenetic marks may contribute to the increased progranulin expression in cases of reduced penetrance. These findings may stimulate similar follow-up studies and new therapeutic approaches.

Estimating the Posttest Probability of Long QT Syndrome Diagnosis for Rare KCNH2 Variants

BACKGROUND

The proliferation of genetic profiling has revealed many associations between genetic variations and disease. However, large-scale phenotyping efforts in largely healthy populations, coupled with DNA sequencing, suggest variants currently annotated as pathogenic are more common in healthy populations than previously thought. In addition, novel and rare variants are frequently observed in genes associated with disease both in healthy individuals and those under suspicion of disease. This raises the question of whether these variants can be useful predictors of disease. To answer this question, we assessed the degree to which the presence of a variant in the cardiac potassium channel gene KCNH2 was diagnostically predictive for the autosomal dominant long QT syndrome.

METHODS

We estimated the probability of a long QT diagnosis given the presence of each KCNH2 variant using Bayesian methods that incorporated variant features such as changes in variant function, protein structure, and in silico predictions. We call this estimate the posttest probability of disease. Our method was applied to over 4000 individuals heterozygous for 871 missense or in-frame insertion/deletion variants in KCNH2 and validated against a separate international cohort of 933 individuals heterozygous for 266 missense or in-frame insertion/deletion variants.

RESULTS

Our method was well-calibrated for the observed fraction of heterozygotes diagnosed with long QT syndrome. Heuristically, we found that the innate diagnostic information one learns about a variant from 3-dimensional variant location, in vitro functional data, and in silico predictors is equivalent to the diagnostic information one learns about that same variant by clinically phenotyping 10 heterozygotes. Most importantly, these data can be obtained in the absence of any clinical observations.

CONCLUSIONS

We show how variant-specific features can inform a prior probability of disease for rare variants even in the absence of clinically phenotyped heterozygotes.

Burden of Rare Variants in the OTOG Gene in Familial Meniere's Disease

OBJECTIVES

Meniere's disease (MD) is a rare inner ear disorder characterized by sensorineural hearing loss, episodic vertigo, and tinnitus. Familial MD has been reported in 6 to 9% of sporadic cases, and few genes including FAM136A, DTNA, PRKCB, SEMA3D, and DPT have been involved in single families, suggesting genetic heterogeneity. In this study, the authors recruited 46 families with MD to search for relevant candidate genes for hearing loss in familial MD.

DESIGN

Exome sequencing data from MD patients were analyzed to search for rare variants in hearing loss genes in a case-control study. A total of 109 patients with MD (73 familial cases and 36 early-onset sporadic patients) diagnosed according to the diagnostic criteria defined by the Barany Society were recruited in 11 hospitals. The allelic frequencies of rare variants in hearing loss genes were calculated in individuals with familial MD. A single rare variant analysis and a gene burden analysis (GBA) were conducted in the dataset selecting 1 patient from each family. Allelic frequencies from European and Spanish reference datasets were used as controls.

RESULTS

A total of 5136 single-nucleotide variants in hearing loss genes were considered for single rare variant analysis in familial MD cases, but only 1 heterozygous likely pathogenic variant in the OTOG gene (rs552304627) was found in 2 unrelated families. The gene burden analysis found an enrichment of rare missense variants in the OTOG gene in familial MD. So, 15 of 46 families (33%) showed at least 1 rare missense variant in the OTOG gene, suggesting a key role in familial MD.

CONCLUSIONS

The authors found an enrichment of multiplex rare missense variants in the OTOG gene in familial MD. This finding supports OTOG as a relevant gene in familial MD and set the groundwork for genetic testing in MD.

High shear stress enhances endothelial permeability in the presence of the risk haplotype at 9p21.3

Single nucleotide polymorphisms (SNPs) are exceedingly common in non-coding loci, and while they are significantly associated with a myriad of diseases, their specific impact on cellular dysfunction remains unclear. Here, we show that when exposed to external stressors, the presence of risk SNPs in the 9p21.3 coronary artery disease (CAD) risk locus increases endothelial monolayer and microvessel dysfunction. Endothelial cells (ECs) derived from induced pluripotent stem cells of patients carrying the risk haplotype (R/R WT) differentiated similarly to their non-risk and isogenic knockout (R/R KO) counterparts. Monolayers exhibited greater permeability and reactive oxygen species signaling when the risk haplotype was present. Addition of the inflammatory cytokine TNFα further enhanced EC monolayer permeability but independent of risk haplotype; TNFα also did not substantially alter haplotype transcriptomes. Conversely, when wall shear stress was applied to ECs in a microfluidic vessel, R/R WT vessels were more permeable at lower shear stresses than R/R KO vessels. Transcriptomes of sheared cells clustered more by risk haplotype than by patient or clone, resulting in significant differential regulation of EC adhesion and extracellular matrix genes vs static conditions. A subset of previously identified CAD risk genes invert expression patterns in the presence of high shear concomitant with altered cell adhesion genes, vessel permeability, and endothelial erosion in the presence of the risk haplotype, suggesting that shear stress could be a regulator of non-coding loci with a key impact on CAD.

An overview of germline variations in genes of primary immunodeficiences through integrative analysis of ClinVar, HGMD® and dbSNP databases

Primary immunodeficiencies (PID) are a diverse group of genetic disorders caused by inadequate development and function of immune system. Identifying genetic etiology is important for genetic counselling and treatment decisions. Clinical relevance of genetic variants is a complex problem depending on gene-specific and variant specific genotype-phenotype interactions. To address this challenge, we aimed to characterize the pathogenic landscape of PID genes by combining the analysis of germline variations reported in ClinVar and HGMD^® and identification of damaging variations available in dbSNP. We generated a joint ClinVar/HGMD database, which included 111,940 variants, among them 32,452 were classified as pathogenic/likely pathogenic. From a total of 5,415,794 bi- or multiallelic variants in PID genes recorded in dbSNP, we retrieved 38,291 high impact (HI) biallelic variants with presumably disruptive impact in the protein, of them 25,500 variants were not present in ClinVar/HGMD. Using a functional prediction algorithm, we additionally identified 28,507 deleterious and 56,016 neutral missense variants among dbSNP variants and created a collection of damaging and neutral variations in PID genes, not currently present in ClinVar/HGMD, with their allele frequencies and mappings to protein domains. The distribution of pathogenic variants from ClinVar/HGMD, HI variants and deleterious missense variants from dbSNP was analyzed in the context of hereditary pattern and gene specific metrics, such as pLI and haploinsufficiency. Our report summarized data on complex gene-specific variability in PID genes and might be useful for the identification of the most promising variants and gene regions for further study.

LONGITUDINAL STUDY OF RPE65-ASSOCIATED INHERITED RETINAL DEGENERATIONS

PURPOSE

To study the disease course of RPE65-associated inherited retinal degenerations (IRDs) as a function of the genotype, define a critical age for blindness, and identify potential modifiers.

METHODS

Forty-five patients with IRD from 33 families with biallelic RPE65 mutations, 28 stemming from a genetic isolate. We collected retrospective data from medical charts. Coexisting variants in 108 IRD-associated genes were identified with Molecular Inversion Probe analysis.

RESULTS

Most patients were diagnosed within the first years of life. Daytime visual function ranged from near-normal to blindness in the first four decades and met WHO criteria for blindness for visual acuity and visual field in the fifth decade. p.(Thr368His) was the most common variant (54%). Intrafamilial variability and interfamilial variability in disease severity and progression were observed. Molecular Inversion Probe analysis confirmed all RPE65 variants and identified one additional variant in LRAT and one in EYS in two separate patients.

CONCLUSION

All patients with RPE65-associated IRDs developed symptoms within the first year of life. Visual function in childhood and adolescence varied but deteriorated inevitably toward blindness after age 40. In this study, genotype was not predictive of clinical course. The variance in severity of disease could not be explained by double hits in other IRD genes.

Integrated Genome and Transcriptome Sequencing to Solve a Neuromuscular Puzzle: Miyoshi Muscular Dystrophy and Early Onset Primary Dystonia in Siblings of the Same Family

BACKGROUND

Neuromuscular disorders (NMD), many of which are hereditary, affect muscular function. Due to advances in high-throughput sequencing technologies, the diagnosis of hereditary NMDs has dramatically improved in recent years.

METHODS AND RESULTS

In this study, we report an family with two siblings exhibiting two different NMD, Miyoshi muscular dystrophy (MMD) and early onset primary dystonia (EOPD). Whole exome sequencing (WES) identified a novel monoallelic frameshift deletion mutation (dysferlin: c.4404delC/p.I1469Sfs∗17) in the Dysferlin gene in the index patient who suffered from MMD. This deletion was inherited from his unaffected father and was carried by his younger sister with EOPD. However, immunostaining staining revealed an absence of dysferlin expression in the proband's muscle tissue and thus suggested the presence of the second underlying mutant allele in dysferlin. Using integrated RNA sequencing (RNA-seq) and whole genome sequencing (WGS) of muscle tissue, a novel deep intronic mutation in dysferlin (dysferlin: c.5341-415A > G) was discovered in the index patient. This mutation caused aberrant mRNA splicing and inclusion of an additional pseudoexon (PE) which we termed PE48.1. This PE was inherited from his unaffected mother. PE48.1 inclusion altered the Dysferlin sequence, causing premature termination of translation.

CONCLUSION

Using integrated genome and transcriptome sequencing, we discovered hereditary MMD and EOPD affecting two siblings of same family. Our results added further weight to the combined use of RNA-seq and WGS as an important method for detection of deep intronic gene mutations, and suggest that integrated sequencing assays are an effective strategy for the diagnosis of hereditary NMDs.

Shifting landscapes of human MTHFR missense-variant effects

Most rare clinical missense variants cannot currently be classified as pathogenic or benign. Deficiency in human 5,10-methylenetetrahydrofolate reductase (MTHFR), the most common inherited disorder of folate metabolism, is caused primarily by rare missense variants. Further complicating variant interpretation, variant impacts often depend on environment. An important example of this phenomenon is the MTHFR variant p.Ala222Val (c.665C>T), which is carried by half of all humans and has a phenotypic impact that depends on dietary folate. Here we describe the results of 98,336 variant functional-impact assays, covering nearly all possible MTHFR amino acid substitutions in four folinate environments, each in the presence and absence of p.Ala222Val. The resulting atlas of MTHFR variant effects reveals many complex dependencies on both folinate and p.Ala222Val. MTHFR atlas scores can distinguish pathogenic from benign variants and, among individuals with severe MTHFR deficiency, correlate with age of disease onset. Providing a powerful tool for understanding structure-function relationships, the atlas suggests a role for a disordered loop in retaining cofactor at the active site and identifies variants that enable escape of inhibition by S-adenosylmethionine. Thus, a model based on eight MTHFR variant effect maps illustrates how shifting landscapes of environment- and genetic-background-dependent missense variation can inform our clinical, structural, and functional understanding of MTHFR deficiency.

Novel missense mutation of SASH1 in a Chinese family with dyschromatosis universalis hereditaria

BACKGROUND

Dyschromatosis universalis hereditaria (DUH) is a pigmentary dermatosis characterized by generalized mottled macules with hypopigmention and hyperpigmention. ABCB6 and SASH1 are recently reported pathogenic genes related to DUH, and the aim of this study was to identify the causative mutations in a Chinese family with DUH.

METHODS

Sanger sequencing was performed to investigate the clinical manifestation and molecular genetic basis of these familial cases of DUH, bioinformatics tools and multiple sequence alignment were used to analyse the pathogenicity of mutations.

RESULTS

A novel missense mutation, c.1529G>A, in the SASH1 gene was identified, and this mutation was not found in the National Center for Biotechnology Information Database of Short Genetic Variation, Online Mendelian Inheritance in Man, ClinVar, or 1000 Genomes Project databases. All in silico predictors suggested that the observed substitution mutation was deleterious. Furthermore, multiple sequence alignment of SASH1 revealed that the p.S510N mutation was highly conserved during evolution. In addition, we reviewed the previously reported DUH-related gene mutations in SASH1 and ABCB6.

CONCLUSION

Although the affected family members had identical mutations, differences in the clinical manifestations of these family members were observed, which reveals the complexity of the phenotype-influencing factors in DUH. Our findings reveal the mutation responsible for DUH in this family and broaden the mutational spectrum of the SASH1 gene.

Challenges in the diagnosis and discovery of rare genetic disorders using contemporary sequencing technologies

Next generation sequencing (NGS) has revolutionised rare disease diagnostics. Concomitant with advancing technologies has been a rise in the number of new gene disorders discovered and diagnoses made for patients and their families. However, despite the trend towards whole exome and whole genome sequencing, diagnostic rates remain suboptimal. On average, only ~30% of patients receive a molecular diagnosis. National sequencing projects launched in the last 5 years are integrating clinical diagnostic testing with research avenues to widen the spectrum of known genetic disorders. Consequently, efforts to diagnose genetic disorders in a clinical setting are now often shared with efforts to prioritise candidate variants for the detection of new disease genes. Herein we discuss some of the biggest obstacles precluding molecular diagnosis and discovery of new gene disorders. We consider bioinformatic and analytical challenges faced when interpreting next generation sequencing data and showcase some of the newest tools available to mitigate these issues. We consider how incomplete penetrance, non-coding variation and structural variants are likely to impact diagnostic rates, and we further discuss methods for uplifting novel gene discovery by adopting a gene-to-patient-based approach.

Unraveling the genetic complexities of combined retinal dystrophy and hearing impairment

Usher syndrome, the most prevalent cause of combined hereditary vision and hearing impairment, is clinically and genetically heterogeneous. Moreover, several conditions with phenotypes overlapping Usher syndrome have been described. This makes the molecular diagnosis of hereditary deaf–blindness challenging. Here, we performed exome sequencing and analysis on 7 Mexican and 52 Iranian probands with combined retinal degeneration and hearing impairment (without intellectual disability). Clinical assessment involved ophthalmological examination and hearing loss questionnaire. Usher syndrome, most frequently due to biallelic variants in MYO7A (USH1B in 16 probands), USH2A (17 probands), and ADGRV1 (USH2C in 7 probands), was diagnosed in 44 of 59 (75%) unrelated probands. Almost half of the identified variants were novel. Nine of 59 (15%) probands displayed other genetic entities with dual sensory impairment, including Alström syndrome (3 patients), cone-rod dystrophy and hearing loss 1 (2 probands), and Heimler syndrome (1 patient). Unexpected findings included one proband each with Scheie syndrome, coenzyme Q10 deficiency, and pseudoxanthoma elasticum. In four probands, including three Usher cases, dual sensory impairment was either modified/aggravated or caused by variants in distinct genes associated with retinal degeneration and/or hearing loss. The overall diagnostic yield of whole exome analysis in our deaf–blind cohort was 92%. Two (3%) probands were partially solved and only 3 (5%) remained without any molecular diagnosis. In many cases, the molecular diagnosis is important to guide genetic counseling, to support prognostic outcomes and decisions with currently available and evolving treatment modalities.

Homozygous WNT9B variants in two families with bilateral renal agenesis/hypoplasia/dysplasia

WNT9B plays a key role in the development of the mammalian urogenital system. It is essential for the induction of mesonephric and metanephric tubules, the regulation of renal tubule morphogenesis, and the regulation of renal progenitor cell expansion and differentiation. To our knowledge, WNT9B has not been associated with renal defects in humans; however, WNT9B^-/- mice have renal agenesis/hypoplasia and reproductive tract abnormalities. We report four individuals from two unrelated consanguineous families with bilateral renal agenesis/hypoplasia/dysplasia and homozygous variants in WNT9B. The proband from Family 1 has bilateral renal cystic dysplasia and chronic kidney disease. He has two deceased siblings who presented with bilateral renal hypoplasia/agenesis. The three affected family members were homozygous for a missense variant in WNT9B (NM_003396.2: c.949G>A/p.(Gly317Arg)). The proband from Family 2 has renal hypoplasia/dysplasia, chronic kidney disease, and is homozygous for a nonsense variant in WNT9B (NM_003396.2: c.11dupC/p.(Pro5Alafs*52)). Two of her siblings died in the neonatal period, one confirmed to be in the context of oligohydramnios. The proband's unaffected brother is also homozygous for the nonsense variant in WNT9B, suggesting nonpenetrance. We propose a novel association of WNT9B and renal anomalies in humans. Further study is needed to delineate the contribution of WNT9B to genitourinary anomalies in humans.

Determinants of penetrance and variable expressivity in monogenic metabolic conditions across 77,184 exomes

Hundreds of thousands of genetic variants have been reported to cause severe monogenic diseases, but the probability that a variant carrier develops the disease (termed penetrance) is unknown for virtually all of them. Additionally, the clinical utility of common polygenetic variation remains uncertain. Using exome sequencing from 77,184 adult individuals (38,618 multi-ancestral individuals from a type 2 diabetes case-control study and 38,566 participants from the UK Biobank, for whom genotype array data were also available), we apply clinical standard-of-care gene variant curation for eight monogenic metabolic conditions. Rare variants causing monogenic diabetes and dyslipidemias display effect sizes significantly larger than the top 1% of the corresponding polygenic scores. Nevertheless, penetrance estimates for monogenic variant carriers average 60% or lower for most conditions. We assess epidemiologic and genetic factors contributing to risk prediction in monogenic variant carriers, demonstrating that inclusion of polygenic variation significantly improves biomarker estimation for two monogenic dyslipidemias.