Determinants of Power in Gene-Based Burden Testing for Monogenic Disorders

Digenic Inheritance in Rare Disorders and Mitochondrial Disease-Crossing the Frontier to a More Comprehensive Understanding of Etiology

Our understanding of rare disease genetics has been shaped by a monogenic disease model. While the traditional monogenic disease model has been successful in identifying numerous disease-associated genes and significantly enlarged our knowledge in the field of human genetics, it has limitations in explaining phenomena like phenotypic variability and reduced penetrance. Widening the perspective beyond Mendelian inheritance has the potential to enable a better understanding of disease complexity in rare disorders. Digenic inheritance is the simplest instance of a non-Mendelian disorder, characterized by the functional interplay of variants in two disease-contributing genes. Known digenic disease causes show a range of pathomechanisms underlying digenic interplay, including direct and indirect gene product interactions as well as epigenetic modifications. This review aims to systematically explore the background of digenic inheritance in rare disorders, the approaches and challenges when investigating digenic inheritance, and the current evidence for digenic inheritance in mitochondrial disorders.

PCDHA9 as a candidate gene for amyotrophic lateral sclerosis

Amyotrophic lateral sclerosis (ALS) is a devastating neurodegenerative disease. To identify additional genetic factors, we analyzed exome sequences in a large cohort of Chinese ALS patients and found a homozygous variant (p.L700P) in PCDHA9 in three unrelated patients. We generated Pcdhα9 mutant mice harboring either orthologous point mutation or deletion mutation. These mice develop progressive spinal motor loss, muscle atrophy, and structural/functional abnormalities of the neuromuscular junction, leading to paralysis and early lethality. TDP-43 pathology is detected in the spinal motor neurons of aged mutant mice. Mechanistically, we demonstrate that Pcdha9 mutation causes aberrant activation of FAK and PYK2 in aging spinal cord, and dramatically reduced NKA-α1 expression in motor neurons. Our single nucleus multi-omics analysis reveals disturbed signaling involved in cell adhesion, ion transport, synapse organization, and neuronal survival in aged mutant mice. Together, our results present PCDHA9 as a potential ALS gene and provide insights into its pathogenesis.

Whole exome sequencing reveals genetic landscape associated with left ventricular outflow tract obstruction in Chinese Han population

Left ventricular outflow tract obstruction (LVOTO), a major form of outflow tract malformation, accounts for a substantial portion of congenital heart defects (CHDs). Unlike its prevalence, the genetic architecture of LVOTO remains largely unknown. To unveil the genetic mutations and risk genes potentially associated with LVOTO, we enrolled a cohort of 106 LVOTO patients and 100 healthy controls and performed a whole-exome sequencing (WES). 71,430 rare deleterious mutations were found in LVOTO patients. By using gene-based burden testing, we further found 32 candidate genes enriched in LVOTO patient including known pathological genes such as GATA5 and GATA6. Most variants of 32 risk genes occur simultaneously rather exclusively suggesting polygenic inherence of LVOTO and 14 genes out of 32 risk genes interact with previously discovered CHD genes. Single cell RNA-seq further revealed dynamic expressions of GATA5, GATA6, FOXD3 and MYO6 in endocardium and neural crest lineage indicating the mutations of these genes lead to LVOTO possibly through different lineages. These findings uncover the genetic architecture of LVOTO which advances the current understanding of LVOTO genetics.

Learning the kernel for rare variant genetic association test

Introduction: Compared to Genome-Wide Association Studies (GWAS) for common variants, single-marker association analysis for rare variants is underpowered. Set-based association analyses for rare variants are powerful tools that capture some of the missing heritability in trait association studies. Methods: We extend the convex-optimized SKAT (cSKAT) test set procedure which learns from data the optimal convex combination of kernels, to the full Generalised Linear Model (GLM) setting with arbitrary non-genetic covariates. We call this extended cSKAT (ecSKAT) and show that the resulting optimization problem is a quadratic programming problem that can be solved with no additional cost compared to cSKAT. Results: We show that a modified objective is related to an upper bound for the p-value through a decreasing exponential term in the objective function, indicating that optimizing this objective function is a principled way of learning the combination of kernels. We evaluate the performance of the proposed method on continuous and binary traits using simulation studies and illustrate its application using UK Biobank Whole Exome Sequencing data on hand grip strength and systemic lupus erythematosus rare variant association analysis. Discussion: Our proposed ecSKAT method enables correcting for important confounders in association studies such as age, sex or population structure for both quantitative and binary traits. Simulation studies showed that ecSKAT can recover sensible weights and achieve higher power across different sample sizes and misspecification settings. Compared to the burden test and SKAT method, ecSKAT gives a lower p-value for the genes tested in both quantitative and binary traits in the UKBiobank cohort.

Rare variant association on unrelated individuals in case-control studies using aggregation tests: existing methods and current limitations

Over the past years, progress made in next-generation sequencing technologies and bioinformatics have sparked a surge in association studies. Especially, genome-wide association studies (GWASs) have demonstrated their effectiveness in identifying disease associations with common genetic variants. Yet, rare variants can contribute to additional disease risk or trait heterogeneity. Because GWASs are underpowered for detecting association with such variants, numerous statistical methods have been recently proposed. Aggregation tests collapse multiple rare variants within a genetic region (e.g. gene, gene set, genomic loci) to test for association. An increasing number of studies using such methods successfully identified trait-associated rare variants and led to a better understanding of the underlying disease mechanism. In this review, we compare existing aggregation tests, their statistical features and scope of application, splitting them into the five classical classes: burden, adaptive burden, variance-component, omnibus and other. Finally, we describe some limitations of current aggregation tests, highlighting potential direction for further investigations.

Excalibur: A new ensemble method based on an optimal combination of aggregation tests for rare-variant association testing for sequencing data

The development of high-throughput next-generation sequencing technologies and large-scale genetic association studies produced numerous advances in the biostatistics field. Various aggregation tests, i.e. statistical methods that analyze associations of a trait with multiple markers within a genomic region, have produced a variety of novel discoveries. Notwithstanding their usefulness, there is no single test that fits all needs, each suffering from specific drawbacks. Selecting the right aggregation test, while considering an unknown underlying genetic model of the disease, remains an important challenge. Here we propose a new ensemble method, called Excalibur, based on an optimal combination of 36 aggregation tests created after an in-depth study of the limitations of each test and their impact on the quality of result. Our findings demonstrate the ability of our method to control type I error and illustrate that it offers the best average power across all scenarios. The proposed method allows for novel advances in Whole Exome/Genome sequencing association studies, able to handle a wide range of association models, providing researchers with an optimal aggregation analysis for the genetic regions of interest.

An exome-wide study of renal operational tolerance

Background

Renal operational tolerance is a rare and beneficial state of prolonged renal allograft function in the absence of immunosuppression. The underlying mechanisms are unknown. We hypothesized that tolerance might be driven by inherited protein coding genetic variants with large effect, at least in some patients.

Methods

We set up a European survey of over 218,000 renal transplant recipients and collected DNAs from 40 transplant recipients who maintained good allograft function without immunosuppression for at least 1 year. We performed an exome-wide association study comparing the distribution of moderate to high impact variants in 36 tolerant patients, selected for genetic homogeneity using principal component analysis, and 192 controls, using an optimal sequence-kernel association test adjusted for small samples.

Results

We identified rare variants of HOMER2 (3/36, FDR 0.0387), IQCH (5/36, FDR 0.0362), and LCN2 (3/36, FDR 0.102) in 10 tolerant patients vs. 0 controls. One patient carried a variant in both HOMER2 and LCN2. Furthermore, the three genes showed an identical variant in two patients each. The three genes are expressed at the primary cilium, a key structure in immune responses.

Conclusion

Rare protein coding variants are associated with operational tolerance in a sizable portion of patients. Our findings have important implications for a better understanding of immune tolerance in transplantation and other fields of medicine.ClinicalTrials.gov, identifier: NCT05124444.

Genetic background of idiopathic neurodevelopmental delay patients with significant brain deviation volume

BACKGROUND

Significant brain volume deviation is an essential phenotype in children with neurodevelopmental delay (NDD), but its genetic basis has not been fully characterized. This study attempted to analyze the genetic factors associated with significant whole-brain deviation volume (WBDV).

METHODS

We established a reference curve based on 4222 subjects ranging in age from the first postnatal day to 18 years. We recruited only NDD patients without acquired etiologies or positive genetic results. Cranial magnetic resonance imaging (MRI) and clinical exome sequencing (2742 genes) data were acquired. A genetic burden test was performed, and the results were compared between patients with and without significant WBDV. Literature review analyses and BrainSpan analysis based on the human brain developmental transcriptome were performed to detect the potential role of genetic risk factors in human brain development.

RESULTS

We recruited a total of 253 NDD patients. Among them, 26 had significantly decreased WBDV (<-2 standard deviations [SDs]), and 14 had significantly increased WBDV (>+2 SDs). NDD patients with significant WBDV had higher rates of motor development delay (49.8% [106/213] vs . 75.0% [30/40], P  = 0.003) than patients without significant WBDV. Genetic burden analyses found 30 genes with an increased allele frequency of rare variants in patients with significant WBDV. Analyses of the literature further demonstrated that these genes were not randomly identified: burden genes were more related to the brain development than background genes ( P  = 1.656e -9 ). In seven human brain regions related to motor development, we observed burden genes had higher expression before 37-week gestational age than postnatal stages. Functional analyses found that burden genes were enriched in embryonic brain development, with positive regulation of synaptic growth at the neuromuscular junction, positive regulation of deoxyribonucleic acid templated transcription, and response to hormone, and these genes were shown to be expressed in neural progenitors. Based on single cell sequencing analyses, we found TUBB2B gene had elevated expression levels in neural progenitor cells, interneuron, and excitatory neuron and SOX15 had high expression in interneuron and excitatory neuron.

CONCLUSION

Idiopathic NDD patients with significant brain volume changes detected by MRI had an increased prevalence of motor development delay, which could be explained by the genetic differences characterized herein.

Whole exome sequencing study identifies candidate loss of function variants and locus heterogeneity in familial cholesteatoma

Cholesteatoma is a rare progressive disease of the middle ear. Most cases are sporadic, but some patients report a positive family history. Identifying functionally important gene variants associated with this disease has the potential to uncover the molecular basis of cholesteatoma pathology with implications for disease prevention, surveillance, or management. We performed an observational WES study of 21 individuals treated for cholesteatoma who were recruited from ten multiply affected families. These family studies were complemented with gene-level mutational burden analysis. We also applied functional enrichment analyses to identify shared properties and pathways for candidate genes and their products. Filtered data collected from pairs and trios of participants within the ten families revealed 398 rare, loss of function (LOF) variants co-segregating with cholesteatoma in 389 genes. We identified six genes DENND2C, DNAH7, NBEAL1, NEB, PRRC2C, and SHC2, for which we found LOF variants in two or more families. The parallel gene-level analysis of mutation burden identified a significant mutation burden for the genes in the DNAH gene family, which encode products involved in ciliary structure. Functional enrichment analyses identified common pathways for the candidate genes which included GTPase regulator activity, calcium ion binding, and degradation of the extracellular matrix. The number of candidate genes identified and the locus heterogeneity that we describe within and between multiply affected families suggest that the genetic architecture for familial cholesteatoma is complex.

Human pluripotent stem cells for the modelling of retinal pigment epithelium homeostasis and disease: A review

Human pluripotent stem cells (hPSCs), which include induced pluripotent stem cells and embryonic stem cells, are powerful tools for studying human development, physiology and disease, including those affecting the retina. Cells from selected individuals, or specific genetic backgrounds, can be differentiated into distinct cell types allowing the modelling of diseases in a dish for therapeutic development. hPSC‐derived retinal cultures have already been used to successfully model retinal pigment epithelium (RPE) degeneration for various retinal diseases including monogenic conditions and complex disease such as age‐related macular degeneration. Here, we will review the current knowledge gained in understanding the molecular events involved in retinal disease using hPSC‐derived retinal models, in particular RPE models. We will provide examples of various conditions to illustrate the scope of applications associated with the use of hPSC‐derived RPE models.

Whole-Exome Sequencing in Congenital Hypothyroidism Due to Thyroid Dysgenesis

Context: Congenital hypothyroidism due to thyroid dysgenesis (CHTD) is a predominantly sporadic and nonsyndromic (NS) condition of unknown etiology. NS-CHTD shows a 40-fold increase in relative risk among first-degree relatives (1 in 100 compared with a birth prevalence of 1 in 4000 in the general population), but a discordance rate between monozygotic (MZ) twins of 92%. This suggests a two-hit mechanism, combining a genetic predisposition (incomplete penetrance of inherited variants) with postzygotic events (accounting for MZ twin discordance). Objective: To evaluate whether whole-exome sequencing (WES) allows to identify new predisposing genes in NS-CHTD. Methods: We performed a case-control study by comparing the whole exome of 36 nonconsanguineous cases of NS-CHTD (33 with lingual thyroid ectopy and 3 with athyreosis, based on technetium pertechnetate scintigraphy at diagnosis) with that of 301 unaffected controls to assess for enrichment in rare protein-altering variants. We performed an unbiased approach using a gene-based burden with a false discovery rate correction. Moreover, we identified all rare pathogenic and likely pathogenic variants, based on in silico prediction tools, in 27 genes previously associated with congenital hypothyroidism (CH) (thyroid dysgenesis [TD] and dyshormonogenesis). Results: After correction for multiple testing, no enrichment in rare protein-altering variants was observed in NS-CHTD. Pathogenic or likely pathogenic variants (21 variants in 12 CH genes) were identified in 42% of cases. Eight percent of cases had variants in more than one gene (oligogenic group); these were not more severely affected than monogenic cases. Moreover, cases with protein-altering variants in dyshormonogenesis-related genes were not more severely affected than those without. Conclusions: No new predisposing genes were identified following an unbiased analysis of WES data in a well-characterized NS-CHTD cohort. Nonetheless, the discovery rate of rare pathogenic or likely pathogenic variants was 42%. Eight percent of the cases harbored multiple variants in genes associated with TD or dyshormonogenesis, but these variants did not explain the variability of hypothyroidism observed in dysgenesis. WES did not identify a genetic cause in NS-CHTD cases, confirming the complex etiology of this disease. Additional studies in larger cohorts and/or novel discovery approaches are required.

Genetic Determinants of Sudden Unexpected Death in Pediatrics

PURPOSE

This study aimed to evaluate genetic contributions to sudden unexpected death in pediatrics (SUDP).

METHODS

We phenotyped and performed exome sequencing for 352 SUDP cases. We analyzed variants in 294 "SUDP genes" with mechanisms plausibly related to sudden death. In a subset of 73 cases with parental data (trios), we performed exome-wide analyses and conducted cohort-wide burden analyses.

RESULTS

In total, we identified likely contributory variants in 37 of 352 probands (11%). Analysis of SUDP genes identified pathogenic/likely pathogenic variants in 12 of 352 cases (SCN1A, DEPDC5 [2], GABRG2, SCN5A [2], TTN [2], MYBPC3, PLN, TNNI3, and PDHA1) and variants of unknown significance-favor-pathogenic in 17 of 352 cases. Exome-wide analyses of the 73 cases with family data additionally identified 4 de novo pathogenic/likely pathogenic variants (SCN1A [2], ANKRD1, and BRPF1) and 4 de novo variants of unknown significance-favor-pathogenic. Comparing cases with controls, we demonstrated an excess burden of rare damaging SUDP gene variants (odds ratio, 2.94; 95% confidence interval, 2.37-4.21) and of exome-wide de novo variants in the subset of 73 with trio data (odds ratio, 3.13; 95% confidence interval, 1.91-5.16).

CONCLUSION

We provide strong evidence for a role of genetic factors in SUDP, involving both candidate genes and novel genes for SUDP and expanding phenotypes of disease genes not previously associated with sudden death.

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.

Constitutively active SARM1 variants that induce neuropathy are enriched in ALS patients

Background

In response to injury, neurons activate a program of organized axon self-destruction initiated by the NAD⁺ hydrolase, SARM1. In healthy neurons SARM1 is autoinhibited, but single amino acid changes can abolish autoinhibition leading to constitutively active SARM1 enzymes that promote degeneration when expressed in cultured neurons.

Methods

To investigate whether naturally occurring human variants might disrupt SARM1 autoinhibition and potentially contribute to risk for neurodegenerative disease, we assayed the enzymatic activity of all 42 rare SARM1 alleles identified among 8507 amyotrophic lateral sclerosis (ALS) patients and 9671 controls. We then intrathecally injected mice with virus expressing SARM1 constructs to test the capacity of an ALS-associated constitutively active SARM1 variant to promote neurodegeneration in vivo.

Results

Twelve out of 42 SARM1 missense variants or small in-frame deletions assayed exhibit constitutive NADase activity, including more than half of those that are unique to the ALS patients or that occur in multiple patients. There is a > 5-fold enrichment of constitutively active variants among patients compared to controls. Expression of constitutively active ALS-associated SARM1 alleles in cultured dorsal root ganglion (DRG) neurons is pro-degenerative and cytotoxic. Intrathecal injection of an AAV expressing the common SARM1 reference allele is innocuous to mice, but a construct harboring SARM1^V184G, the constitutively active variant found most frequently among the ALS patients, causes axon loss, motor dysfunction, and sustained neuroinflammation.

Conclusions

These results implicate rare hypermorphic SARM1 alleles as candidate genetic risk factors for ALS and other neurodegenerative conditions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13024-021-00511-x.

Association of apolipoprotein E ε4 allele and amyotrophic lateral sclerosis in Chinese population

Background: Amyotrophic lateral sclerosis (ALS) has a complex genetic origin, and how immune dysregulation may contribute to ALS etiology remain unclear. Given the roles played by apolipoprotein E (APOE) signaling in neuroinflammation and neurodegeneration, an improved knowledge of the association between APOE genotypes and ALS risk in Chinese population may help to understand the underlying etiology of the disease. Methods: A retrospective case-control study with participants of Chinese ancestry was conducted, with a total of 683 ALS patients and 369 healthy controls analyzed for APOE genotypes using Sanger sequencing. In addition, 282 of these patients were further analyzed for known ALS risk variants and rare deleterious variants related to immune disorders via whole exome sequencing. Results: Among the 683 ALS patients analyzed (346 males, 337 females; mean age at onset [SD]: 51.9 [10.9]), 145 patients (21.1%) carried ε4, the proportion of which was significantly higher than 16.0% in controls (59/369; OR, 1.42; 95%CI, 1.02-1.98; p = 0.02). There is no evidence supporting the association between APOE genotypes and disease phenotypes. We also didn't find any enrichment of currently known ALS risk variants or variants in genes related to immune abnormality in specific APOE genotypes. Conclusion: Our study highlighted the importance of trans-ethnic studies in identifying genetic risk factors, and the relevance of APOE in ALS etiopathogenesis in Chinese population.

VINYL: Variant prIoritizatioN by survivaL analysis

MOTIVATION

Clinical applications of genome re-sequencing technologies typically generate large amounts of data that need to be carefully annotated and interpreted to identify genetic variants potentially associated with pathological conditions. In this context, accurate and reproducible methods for the functional annotation and prioritization of genetic variants are of fundamental importance.

RESULTS

In this paper, we present VINYL, a flexible and fully automated system for the functional annotation and prioritization of genetic variants. Extensive analyses of both real and simulated datasets suggest that VINYL can identify clinically relevant genetic variants in a more accurate manner compared to equivalent state of the art methods, allowing a more rapid and effective prioritization of genetic variants in different experimental settings. As such we believe that VINYL can establish itself as a valuable tool to assist healthcare operators and researchers in clinical genomics investigations.

AVAILABILITY

VINYL is available at http://beaconlab.it/VINYL and https://github.com/matteo14c/VINYL.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Novel Mutations and Genes That Impact on Growth in Short Stature of Undefined Aetiology: The EPIGROW Study

BACKGROUND

Children with short stature of undefined aetiology (SS-UA) may have undiagnosed genetic conditions.

PURPOSE

To identify mutations causing short stature (SS) and genes related to SS, using candidate gene sequence data from the European EPIGROW study.

METHODS

First, we selected exonic single nucleotide polymorphisms (SNPs), in cases and not controls, with minor allele frequency (MAF) < 2%, whose carriage fitted the mode of inheritance. Known mutations were identified using Ensembl and gene-specific databases. Variants were classified as pathogenic, likely pathogenic, or variant of uncertain significance using criteria from the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. If predicted by ≥ 5/10 algorithms (eg, Polyphen2) to be deleterious, this was considered supporting evidence of pathogenicity. Second, gene-based burden testing determined the difference in SNP frequencies between cases and controls across all and then rare SNPs. For genotype/phenotype relationships, we used PLINK, based on haplotype, MAF > 2%, genotype present in > 75%, and Hardy Weinberg equilibrium P > 10-4.

RESULTS

First, a diagnostic yield of 10% (27/263) was generated by 2 pathogenic (nonsense in ACAN) and a further 25 likely pathogenic mutations, including previously known missense mutations in FANCB, IGFIR, MMP13, NPR2, OBSL1, and PTPN11. Second, genes related to SS: all methods identified PEX2. Another 7 genes (BUB1B, FANCM, CUL7, FANCA, PTCH1, TEAD3, BCAS3) were identified by both gene-based approaches and 6 (A2M, EFEMP1, PRKCH, SOS2, RNF135, ZBTB38) were identified by gene-based testing for all SNPs and PLINK.

CONCLUSIONS

Such panels improve diagnosis in SS-UA, extending known disease phenotypes. Fourteen genes related to SS included some known to cause growth disorders as well as novel targets.

Germline Genetic Risk Variants for Progressive Multifocal Leukoencephalopathy

Progressive multifocal leukoencephalopathy (PML) is a rare demyelinating disorder of the brain caused by reactivation of the JC virus (JCV), a polyomavirus that infects at least 60% of the population but is asymptomatic or results in benign symptoms in most people. PML occurs as a secondary disease in a variety of disorders or as a serious adverse event from immunosuppressant agents, but is mainly found in three groups: HIV-infected patients, patients with hematological malignancies, or multiple sclerosis (MS) patients on the immunosuppressant therapy natalizumab. It is severely debilitating and is deadly in ~50% HIV cases, ~90% of hematological malignancy cases, and ~24% of MS-natalizumab cases. A PML risk prediction test would have clinical utility in all at risk patient groups but would be particularly beneficial in patients considering therapy with immunosuppressant agents known to cause PML, such as natalizumab, rituximab, and others. While a JC antibody test is currently used in the clinical decision process for natalizumab, it is suboptimal because of its low specificity and requirement to periodically retest patients for seroconversion or to assess if a patient's JCV index has increased. Whereas a high specificity genetic risk prediction test comprising host genetic risk variants (i.e., germline variants occurring at higher frequency in PML patients compared to the general population) could be administered one time to provide clinicians with additional risk prediction information that is independent of JCV serostatus. Prior PML case reports support the hypothesis that PML risk is greater in patients with a genetically caused immunodeficiency disorder. To identify germline PML risk variants, we performed exome sequencing on 185 PML cases (70 in a discovery cohort and 115 in a replication cohort) and used the gnomAD variant database for interpretation. Our study yielded 19 rare variants (maximum allele frequency of 0.02 in gnomAD ethnically matched populations) that impact 17 immune function genes (10 are known to cause inborn errors of immunity). Modeling of these variants in a PML genetic risk test for MS patients considering natalizumab treatment indicates that at least a quarter of PML cases may be preventable.

Clinical genetic testing in endocrinology: Current concepts and contemporary challenges

Recent advances in DNA sequencing technology have led to an unprecedented period of disease-gene discovery offering many new opportunities for genetic testing in the clinical setting. Endocrinology has seen a rapid expansion in the taxonomy of monogenic disorders, which can be detected by an expanding portfolio of genetic tests in both diagnostic and predictive settings. Successful testing relies on many factors including the ability to identify those at increased risk of genetic disease in the busy clinic as well as a working knowledge of the various testing platforms and their limitations. The clinical utility of a given test is dependent upon many factors, which include the reliability of the genetic testing platform, the accuracy of the test result interpretation and knowledge of disease penetrance and expression. The increasing adoption of "high-content" genetic testing based on next-generation sequencing (NGS) to diagnose hereditary endocrine disorders brings a number of challenges including the potential for uncertain test results and/or genetic findings unrelated to the indication for testing. Therefore, it is increasingly important that the clinician is aware of the current evolution in genetic testing, and understands the different settings in which it may be employed. This review provides an overview of the genetic testing workflow, focusing on each of the major components required for successful testing in adult and paediatric endocrine settings. In addition, the challenges of variant interpretation are highlighted, as are issues related to informed consent, prenatal diagnosis and predictive testing. Finally, the future directions of genetic testing relevant to endocrinology are discussed.

Rare-variant collapsing analyses for complex traits: guidelines and applications

The first phase of genome-wide association studies (GWAS) assessed the role of common variation in human disease. Advances optimizing and economizing high-throughput sequencing have enabled a second phase of association studies that assess the contribution of rare variation to complex disease in all protein-coding genes. Unlike the early microarray-based studies, sequencing-based studies catalogue the full range of genetic variation, including the evolutionarily youngest forms. Although the experience with common variants helped establish relevant standards for genome-wide studies, the analysis of rare variation introduces several challenges that require novel analysis approaches.

Revisiting Non-BRCA1/2 Familial Whole Exome Sequencing Datasets Implicates NCK1 as a Cancer Gene

Through linkage and candidate gene screening, many breast cancer (BC) predisposition genes have been identified in the past 20 years. However, the majority of genetic risks that contribute to familial BC remains undetermined. In this study, we revisited whole exome sequencing datasets from non-BRCA1/2 familial BC patients, to search for novel BC predisposition genes. Based on the infinite mutation model, we supposed that rare non-silent variants that cooccurred between familial and TCGA-germline datasets, might play a predisposition contributing role. In our analysis, we not only identified novel potential pathogenic variants from known cancer predisposition genes, such as MRE11, CTR9 but also identified novel candidate predisposition genes, such as NCK1. According to the TCGA mRNA expression dataset of BC, NCK1 was significantly upregulated in basal-like subtypes and downregulated in luminal subtypes. In vitro, NCK1 mutants (D73H and R42Q) transfected MCF7 cell lines, which attributed to the luminal subtype, were much more viable and invasive than the wild type. On the other side, our results also showed that overall survival and disease-free survival of patients with NCK1 variations might be dependent on the genomic context. In conclusion, genetic heterogeneity exists among non-BRCA1/2 BC pedigrees and NCK1 could be a novel BC predisposition gene.

Exome-Based Rare-Variant Analyses in CKD

BACKGROUND

Studies have identified many common genetic associations that influence renal function and all-cause CKD, but these explain only a small fraction of variance in these traits. The contribution of rare variants has not been systematically examined.

METHODS

We performed exome sequencing of 3150 individuals, who collectively encompassed diverse CKD subtypes, and 9563 controls. To detect causal genes and evaluate the contribution of rare variants we used collapsing analysis, in which we compared the proportion of cases and controls carrying rare variants per gene.

RESULTS

The analyses captured five established monogenic causes of CKD: variants in PKD1, PKD2, and COL4A5 achieved study-wide significance, and we observed suggestive case enrichment for COL4A4 and COL4A3. Beyond known disease-associated genes, collapsing analyses incorporating regional variant intolerance identified suggestive dominant signals in CPT2 and several other candidate genes. Biallelic mutations in CPT2 cause carnitine palmitoyltransferase II deficiency, sometimes associated with rhabdomyolysis and acute renal injury. Genetic modifier analysis among cases with APOL1 risk genotypes identified a suggestive signal in AHDC1, implicated in Xia-Gibbs syndrome, which involves intellectual disability and other features. On the basis of the observed distribution of rare variants, we estimate that a two- to three-fold larger cohort would provide 80% power to implicate new genes for all-cause CKD.

CONCLUSIONS

This study demonstrates that rare-variant collapsing analyses can validate known genes and identify candidate genes and modifiers for kidney disease. In so doing, these findings provide a motivation for larger-scale investigation of rare-variant risk contributions across major clinical CKD categories.

The Genetic Landscape of Diamond-Blackfan Anemia

Diamond-Blackfan anemia (DBA) is a rare bone marrow failure disorder that affects 7 out of 1,000,000 live births and has been associated with mutations in components of the ribosome. In order to characterize the genetic landscape of this heterogeneous disorder, we recruited a cohort of 472 individuals with a clinical diagnosis of DBA and performed whole-exome sequencing (WES). We identified relevant rare and predicted damaging mutations for 78% of individuals. The majority of mutations were singletons, absent from population databases, predicted to cause loss of function, and located in 1 of 19 previously reported ribosomal protein (RP)-encoding genes. Using exon coverage estimates, we identified and validated 31 deletions in RP genes. We also observed an enrichment for extended splice site mutations and validated their diverse effects using RNA sequencing in cell lines obtained from individuals with DBA. Leveraging the size of our cohort, we observed robust genotype-phenotype associations with congenital abnormalities and treatment outcomes. We further identified rare mutations in seven previously unreported RP genes that may cause DBA, as well as several distinct disorders that appear to phenocopy DBA, including nine individuals with biallelic CECR1 mutations that result in deficiency of ADA2. However, no new genes were identified at exome-wide significance, suggesting that there are no unidentified genes containing mutations readily identified by WES that explain >5% of DBA-affected case subjects. Overall, this report should inform not only clinical practice for DBA-affected individuals, but also the design and analysis of rare variant studies for heterogeneous Mendelian disorders.

Burden Testing of Rare Variants Identified through Exome Sequencing via Publicly Available Control Data

The genetic causes of many Mendelian disorders remain undefined. Factors such as lack of large multiplex families, locus heterogeneity, and incomplete penetrance hamper these efforts for many disorders. Previous work suggests that gene-based burden testing-where the aggregate burden of rare, protein-altering variants in each gene is compared between case and control subjects-might overcome some of these limitations. The increasing availability of large-scale public sequencing databases such as Genome Aggregation Database (gnomAD) can enable burden testing using these databases as controls, obviating the need for additional control sequencing for each study. However, there exist various challenges with using public databases as controls, including lack of individual-level data, differences in ancestry, and differences in sequencing platforms and data processing. To illustrate the approach of using public data as controls, we analyzed whole-exome sequencing data from 393 individuals with idiopathic hypogonadotropic hypogonadism (IHH), a rare disorder with significant locus heterogeneity and incomplete penetrance against control subjects from gnomAD (n = 123,136). We leveraged presumably benign synonymous variants to calibrate our approach. Through iterative analyses, we systematically addressed and overcame various sources of artifact that can arise when using public control data. In particular, we introduce an approach for highly adaptable variant quality filtering that leads to well-calibrated results. Our approach "re-discovered" genes previously implicated in IHH (FGFR1, TACR3, GNRHR). Furthermore, we identified a significant burden in TYRO3, a gene implicated in hypogonadotropic hypogonadism in mice. Finally, we developed a user-friendly software package TRAPD (Test Rare vAriants with Public Data) for performing gene-based burden testing against public databases.

Cancer genetics, precision prevention and a call to action

More than 15 years have passed since the identification, through linkage, of 'first-wave' susceptibility genes for common cancers (BRCA1, BRCA2, MLH1 and MSH2). These genes have strong frequency-penetrance profiles, such that the associated clinical utility probably remains relevant regardless of the context of ascertainment. 'Second-wave' genes, not tractable by linkage, were subsequently identified by mutation screening of candidate genes (PALB2, ATM, CHEK2, BRIP1, RAD51C and RAD51D). Their innately weaker frequency-penetrance profiles have rendered delineation of cancer associations, risks and variant pathogenicity challenging, thereby compromising their clinical application. Early germline exome-sequencing endeavors for common cancers did not yield the long-anticipated slew of 'next-wave' genes but instead implied a highly polygenic genomic architecture requiring much larger experiments to make any substantive inroads into gene discovery. As such, the 'genetic economics' of frequency penetrance clearly indicates that focused identification of carriers of first-wave-gene mutations is most impactful for cancer control. With screening, prevention and early detection at the forefront of the cancer management agenda, we propose that the time is nigh for the initiation of national population-testing programs to identify carriers of first-wave gene mutation carriers. To fully deliver a precision prevention program, long-term, large-scale mutation studies that capture longitudinal clinical data and serial biosamples are required.

Molecular characterization of exonic rearrangements and frame shifts in the dystrophin gene in Duchenne muscular dystrophy patients in a Saudi community

BACKGROUND

In individuals with Duchenne muscular dystrophy (DMD), exon skipping treatment to restore a wild-type phenotype or correct the frame shift of the mRNA transcript of the dystrophin (DMD) gene are mutation-specific. To explore the molecular characterization of DMD rearrangements and predict the reading frame, we simultaneously screened all 79 DMD gene exons of 45 unrelated male DMD patients using a multiplex ligation-dependent probe amplification (MLPA) assay for deletion/duplication patterns. Multiplex PCR was used to confirm single deletions detected by the MLPA.

RESULTS

There was an obvious diagnostic delay, with an extremely statistically significant difference between the age at initial symptoms and the age of clinical evaluation of DMD cases (t value, 10.3; 95% confidence interval 5.95-8.80, P < 0.0001); the mean difference between the two groups was 7.4 years. Overall, we identified 147 intragenic rearrangements: 46.3% deletions and 53.7% duplications. Most of the deletions (92.5%) were between exons 44 and 56, with exon 50 being the most frequently involved (19.1%). Eight new rearrangements, including a mixed deletion/duplication and double duplications, were linked to seven cases with DMD. Of all the cases, 17.8% had duplications with no hot spots. In addition, confirmation of the reading frame hypothesis helped account for new DMD rearrangements in this study. We found that 81% of our Saudi patients would potentially benefit from exon skipping, of which 42.9% had a mutation amenable to skipping of exon 51.

CONCLUSIONS

Our study could generate considerable data on mutational rearrangements that may promote future experimental therapies in Saudi Arabia.

Molecular characterization of exonic rearrangements and frame shifts in the dystrophin gene in Duchenne muscular dystrophy patients in a Saudi community

BACKGROUND

In individuals with Duchenne muscular dystrophy (DMD), exon skipping treatment to restore a wild-type phenotype or correct the frame shift of the mRNA transcript of the dystrophin (DMD) gene are mutation-specific. To explore the molecular characterization of DMD rearrangements and predict the reading frame, we simultaneously screened all 79 DMD gene exons of 45 unrelated male DMD patients using a multiplex ligation-dependent probe amplification (MLPA) assay for deletion/duplication patterns. Multiplex PCR was used to confirm single deletions detected by the MLPA.

RESULTS

There was an obvious diagnostic delay, with an extremely statistically significant difference between the age at initial symptoms and the age of clinical evaluation of DMD cases (t value, 10.3; 95% confidence interval 5.95-8.80, P < 0.0001); the mean difference between the two groups was 7.4 years. Overall, we identified 147 intragenic rearrangements: 46.3% deletions and 53.7% duplications. Most of the deletions (92.5%) were between exons 44 and 56, with exon 50 being the most frequently involved (19.1%). Eight new rearrangements, including a mixed deletion/duplication and double duplications, were linked to seven cases with DMD. Of all the cases, 17.8% had duplications with no hot spots. In addition, confirmation of the reading frame hypothesis helped account for new DMD rearrangements in this study. We found that 81% of our Saudi patients would potentially benefit from exon skipping, of which 42.9% had a mutation amenable to skipping of exon 51.

CONCLUSIONS

Our study could generate considerable data on mutational rearrangements that may promote future experimental therapies in Saudi Arabia.

Whole-exome sequencing and gene-based rare variant association tests suggest that PLA2G4E might be a risk gene for panic disorder

Panic disorder (PD) is characterized by recurrent and unexpected panic attacks, subsequent anticipatory anxiety, and phobic avoidance. Recent epidemiological and genetic studies have revealed that genetic factors contribute to the pathogenesis of PD. We performed whole-exome sequencing on one Japanese family, including multiple patients with panic disorder, which identified seven rare protein-altering variants. We then screened these genes in a Japanese PD case-control group (384 sporadic PD patients and 571 controls), resulting in the detection of three novel single nucleotide variants as potential candidates for PD (chr15: 42631993, T>C in GANC; chr15: 42342861, G>T in PLA2G4E; chr20: 3641457, G>C in GFRA4). Statistical analyses of these three genes showed that PLA2G4E yielded the lowest p value in gene-based rare variant association tests by Efficient and Parallelizable Association Container Toolbox algorithms; however, the p value did not reach the significance threshold in the Japanese. Likewise, in a German case-control study (96 sporadic PD patients and 96 controls), PLA2G4E showed the lowest p value but again did not reach the significance threshold. In conclusion, we failed to find any significant variants or genes responsible for the development of PD. Nonetheless, our results still leave open the possibility that rare protein-altering variants in PLA2G4E contribute to the risk of PD, considering the function of this gene.

Exome-wide Association Study Identifies GREB1L Mutations in Congenital Kidney Malformations

Renal agenesis and hypodysplasia (RHD) are major causes of pediatric chronic kidney disease and are highly genetically heterogeneous. We conducted whole-exome sequencing in 202 case subjects with RHD and identified diagnostic mutations in genes known to be associated with RHD in 7/202 case subjects. In an additional affected individual with RHD and a congenital heart defect, we found a homozygous loss-of-function (LOF) variant in SLIT3, recapitulating phenotypes reported with Slit3 inactivation in the mouse. To identify genes associated with RHD, we performed an exome-wide association study with 195 unresolved case subjects and 6,905 control subjects. The top signal resided in GREB1L, a gene implicated previously in Hoxb1 and Shha signaling in zebrafish. The significance of the association, which was p = 2.0 × 10^-5 for novel LOF, increased to p = 4.1 × 10^-6 for LOF and deleterious missense variants combined, and augmented further after accounting for segregation and de novo inheritance of rare variants (joint p = 2.3 × 10^-7). Finally, CRISPR/Cas9 disruption or knockdown of greb1l in zebrafish caused specific pronephric defects, which were rescued by wild-type human GREB1L mRNA, but not mRNA containing alleles identified in case subjects. Together, our study provides insight into the genetic landscape of kidney malformations in humans, presents multiple candidates, and identifies SLIT3 and GREB1L as genes implicated in the pathogenesis of RHD.