Homozygous and hemizygous CNV detection from exome sequencing data in a Mendelian disease cohort

Genomic Balancing Act: deciphering DNA rearrangements in the complex chromosomal aberration involving 5p15.2, 2q31.1, and 18q21.32

Despite extensive research into the genetic underpinnings of neurodevelopmental disorders (NDD), many clinical cases remain unresolved. We studied a female proband with a NDD, mildly dysmorphic facial features, and brain stem hypoplasia on neuroimaging. Comprehensive genomic analyses revealed a terminal 5p loss and a terminal 18q gain in the proband while a diploid copy number for chromosomes 5 and 18 in both parents. Genomic investigations in the proband identified an unbalanced translocation t(5;18) with additional genetic material from chromosome 2 (2q31.3) inserted at the breakpoint, pointing to a complex chromosomal rearrangement (CCR) involving 5p15.2, 2q31.3, and 18q21.32. Breakpoint junction analyses enabled by long-read genome sequencing unveiled the presence of four distinct junctions in the father, who is a carrier of a balanced CCR. The proband inherited from the father both the abnormal chromosome 5 resulting in segmental aneusomies of chr5 (loss) and chr18 (gain) and a der(2) homologue. Evidences suggest a chromoplexy mechanism for this CCR derivation, involving double-strand breaks (DSBs) repaired by non-homologous end joining (NHEJ) or alternative end joining (alt-EJ). The complexity of the CCR and the segregation of homologues elucidate the genetic model for this family. This study demonstrates the importance of combining multiple genomic technologies to uncover genetic causes of complex neurodevelopmental syndromes and to better understand genetic disease mechanisms.

Optical genome mapping unveils hidden structural variants in neurodevelopmental disorders

While short-read sequencing currently dominates genetic research and diagnostics, it frequently falls short of capturing certain structural variants (SVs), which are often implicated in the etiology of neurodevelopmental disorders (NDDs). Optical genome mapping (OGM) is an innovative technique capable of capturing SVs that are undetectable or challenging-to-detect via short-read methods. This study aimed to investigate NDDs using OGM, specifically focusing on cases that remained unsolved after standard exome sequencing. OGM was performed in 47 families using ultra-high molecular weight DNA. Single-molecule maps were assembled de novo, followed by SV and copy number variant calling. We identified 7 variants of interest, of which 5 (10.6%) were classified as likely pathogenic or pathogenic, located in BCL11A, OPHN1, PHF8, SON, and NFIA. We also identified an inversion disrupting NAALADL2, a gene which previously was found to harbor complex rearrangements in two NDD cases. Variants in known NDD genes or candidate variants of interest missed by exome sequencing mainly consisted of larger insertions (> 1kbp), inversions, and deletions/duplications of a low number of exons (1-4 exons). In conclusion, in addition to improving molecular diagnosis in NDDs, this technique may also reveal novel NDD genes which may harbor complex SVs often missed by standard sequencing techniques.

A Novel Subtype of Acquired Generalized Lipodystrophy Associated With Subcutaneous Panniculitis-Like T-cell Lymphoma

Acquired generalized lipodystrophy (AGL) is an extremely rare disease that is characterized by loss of body fat affecting nearly all parts of the body. It is often associated with autoimmune diseases or panniculitis, whereas in other patients the underlying etiology is unclear. We report a 52-year-old male individual who was diagnosed with subcutaneous panniculitis-like T-cell lymphoma (SPTCL) that spontaneously went into remission. Years later he developed new subcutaneous nodules most concerning for relapse SPTCL or lupus panniculitis, followed by onset of hemophagocytic lymphohistiocytosis (HLH) that was treated with allogeneic stem cell transplantation. Notably, around the same time, he also developed generalized subcutaneous fat loss of both upper and lower extremities, chest, abdomen, and face that persisted after treatment of the HLH. Whole exome sequencing was performed to search for pathogenic variants that are associated with SPTCL, including those in hepatitis A virus cellular receptor 2 (HAVCR2), but did not detect any potential disease-causing variant. Our report brings to the attention a novel subtype of panniculitis-variety of AGL. Whether generalized loss of subcutaneous fat in this patient is due to lymphoma-associated panniculitis or due to development of adipose tissue-directed autoantibodies as a paraneoplastic "autoimmune" manifestation of SPTCL remains unclear.

NODAL variants are associated with a continuum of laterality defects from simple D-transposition of the great arteries to heterotaxy

BACKGROUND

NODAL signaling plays a critical role in embryonic patterning and heart development in vertebrates. Genetic variants resulting in perturbations of the TGF-β/NODAL signaling pathway have reproducibly been shown to cause laterality defects in humans. To further explore this association and improve genetic diagnosis, the study aims to identify and characterize a broader range of NODAL variants in a large number of individuals with laterality defects.

METHODS

We re-analyzed a cohort of 321 proband-only exomes of individuals with clinically diagnosed laterality congenital heart disease (CHD) using family-based, rare variant genomic analyses. To this cohort we added 12 affected subjects with known NODAL variants and CHD from institutional research and clinical cohorts to investigate an allelic series. For those with candidate contributory variants, variant allele confirmation and segregation analysis were studied by Sanger sequencing in available family members. Array comparative genomic hybridization and droplet digital PCR were utilized for copy number variants (CNV) validation and characterization. We performed Human Phenotype Ontology (HPO)-based quantitative phenotypic analyses to dissect allele-specific phenotypic differences.

RESULTS

Missense, nonsense, splice site, indels, and/or structural variants of NODAL were identified as potential causes of heterotaxy and other laterality defects in 33 CHD cases. We describe a recurrent complex indel variant for which the nucleic acid secondary structure predictions implicate secondary structure mutagenesis as a possible mechanism for formation. We identified two CNV deletion alleles spanning NODAL in two unrelated CHD cases. Furthermore, 17 CHD individuals were found (16/17 with known Hispanic ancestry) to have the c.778G > A:p.G260R NODAL missense variant which we propose reclassification from variant of uncertain significance (VUS) to likely pathogenic. Quantitative HPO-based analyses of the observed clinical phenotype for all cases with p.G260R variation, including heterozygous, homozygous, and compound heterozygous cases, reveal clustering of individuals with biallelic variation. This finding provides evidence for a genotypic-phenotypic correlation and an allele-specific gene dosage model.

CONCLUSION

Our data further support a role for rare deleterious variants in NODAL as a cause for sporadic human laterality defects, expand the repertoire of observed anatomical complexity of potential cardiovascular anomalies, and implicate an allele specific gene dosage model.

Incorporating CNV analysis improves the yield of exome sequencing for rare monogenic disorders-an important consideration for resource-constrained settings

Exome sequencing (ES) is a recommended first-tier diagnostic test for many rare monogenic diseases. It allows for the detection of both single-nucleotide variants (SNVs) and copy number variants (CNVs) in coding exonic regions of the genome in a single test, and this dual analysis is a valuable approach, especially in limited resource settings. Single-nucleotide variants are well studied; however, the incorporation of copy number variant analysis tools into variant calling pipelines has not been implemented yet as a routine diagnostic test, and chromosomal microarray is still more widely used to detect copy number variants. Research shows that combined single and copy number variant analysis can lead to a diagnostic yield of up to 58%, increasing the yield with as much as 18% from the single-nucleotide variant only pipeline. Importantly, this is achieved with the consideration of computational costs only, without incurring any additional sequencing costs. This mini review provides an overview of copy number variant analysis from exome data and what the current recommendations are for this type of analysis. We also present an overview on rare monogenic disease research standard practices in resource-limited settings. We present evidence that integrating copy number variant detection tools into a standard exome sequencing analysis pipeline improves diagnostic yield and should be considered a significantly beneficial addition, with relatively low-cost implications. Routine implementation in underrepresented populations and limited resource settings will promote generation and sharing of CNV datasets and provide momentum to build core centers for this niche within genomic medicine.

Expanding the allelic spectrum of ELOVL4-related autosomal recessive neuro-ichthyosis

BACKGROUND

Very long-chain fatty acids (VLCFAs) composed of more than 20 carbon atoms are essential in the biosynthesis of cell membranes in the brain, skin, and retina. VLCFAs are elongated beyond 28 carbon atoms by ELOVL4 enzyme. Variants in ELOVL4 are associated with three Mendelian disorders: autosomal dominant (AD) Stargardt-like macular dystrophy type 3, AD spinocerebellar ataxia, and autosomal recessive disorder congenital ichthyosis, spastic quadriplegia and impaired intellectual development (ISQMR). Only seven subjects from five unrelated families with ISQMR have been described, all of which have biallelic single-nucleotide variants.

METHODS

We performed clinical exome sequencing on probands from four unrelated families with neuro-ichthyosis.

RESULTS

We identified three novel homozygous ELOVL4 variants. Two of the families originated from the same Saudi tribe and had the exact homozygous exonic deletion in ELOVL4, while the third and fourth probands had two different novel homozygous missense variants. Seven out of the eight affected subjects had profound developmental delay, epilepsy, axial hypotonia, peripheral hypertonia, and ichthyosis. Delayed myelination and corpus callosum hypoplasia were seen in two of five subjects with brain magnetic rosonance imaging and cerebral atrophy in three.

CONCLUSION

Our study expands the allelic spectrum of ELOVL4-related ISQMR. The detection of the same exonic deletion in two unrelated Saudi family from same tribe suggests a tribal founder mutation.

Homozygous 22q11.2 distal type II microdeletion is associated with syndromic neurodevelopmental delay

Genomic disorders result from heterozygous copy number variants (CNVs). Homozygous deletions spanning numerous genes are rare, despite the potential contribution of consanguinity to such instances. CNVs in the 22q11.2 region are mediated by nonallelic homologous recombination between pairs of low copy repeats (LCRs), from amongst eight LCRs designated A-H. Heterozygous distal type II deletions (LCR-E to LCR-F) have incomplete penetrance and variable expressivity, and can lead to neurodevelopmental issues, minor craniofacial anomalies, and congenital abnormalities. We report siblings with global developmental delay, hypotonia, minor craniofacial anomalies, ocular abnormalities, and minor skeletal issues, in whom chromosomal microarray identified a homozygous distal type II deletion. The deletion was brought to homozygosity as a result of a consanguineous marriage between two heterozygous carriers of the deletion. The phenotype of the children was strikingly more severe and complex than that of the parents. This report suggests that the distal type II deletion harbors a dosage-sensitive gene or regulatory element, which leads to a more severe phenotype when deleted on both chromosomes.

Rare variant enrichment analysis supports GREB1L as a contributory driver gene in the etiology of Mayer-Rokitansky-Küster-Hauser syndrome

Mayer-Rokitansky-Küster-Hauser (MRKH) syndrome is characterized by aplasia of the female reproductive tract; the syndrome can include renal anomalies, absence or dysgenesis, and skeletal anomalies. While functional models have elucidated several candidate genes, only WNT4 (MIM: 603490) variants have been definitively associated with a subtype of MRKH with hyperandrogenism (MIM: 158330). DNA from 148 clinically diagnosed MRKH probands across 144 unrelated families and available family members from North America, Europe, and South America were exome sequenced (ES) and by family-based genomics analyzed for rare likely deleterious variants. A replication cohort consisting of 442 Han Chinese individuals with MRKH was used to further reproduce GREB1L findings in diverse genetic backgrounds. Proband and OMIM phenotypes annotated using the Human Phenotype Ontology were analyzed to quantitatively delineate the phenotypic spectrum associated with GREB1L variant alleles found in our MRKH cohort and those previously published. This study reports 18 novel GREB1L variant alleles, 16 within a multiethnic MRKH cohort and two within a congenital scoliosis cohort. Cohort-wide analyses for a burden of rare variants within a single gene identified likely damaging variants in GREB1L (MIM: 617782), a known disease gene for renal hypoplasia and uterine abnormalities (MIM: 617805), in 16 of 590 MRKH probands. GREB1L variant alleles, including a CNV null allele, were found in 8 MRKH type 1 probands and 8 MRKH type II probands. This study used quantitative phenotypic analyses in a worldwide multiethnic cohort to identify and strengthen the association of GREB1L to isolated uterine agenesis (MRKH type I) and syndromic MRKH type II.

Study of pathogenic genes in a pedigree with familial dilated cardiomyopathy

BACKGROUND

Dilated cardiomyopathy (DCM) is a genetically heterogeneous cardiac disorder characterized by left ventricular dilation and contractile dysfunction. The substantial genetic heterogeneity evident in patients with DCM contributes to variable disease severity and complicates overall prognosis, which can be very poor.

AIM

To identify pathogenic genes in DCM through pedigree analysis.

METHODS

Our research team identified a patient with DCM in the clinic. Through investigation, we found that the family of this patient has a typical DCM pedigree. High-throughput sequencing technology, next-generation sequencing, was used to sequence the whole exomes of seven samples in the pedigree.

RESULTS

A novel and potentially pathogenic gene mutation-ANK2p.F3067L-was discovered. The mutation was completely consistent with the clinical information for this DCM pedigree. Sanger sequencing was used to further verify the locus of the mutation in pedigree samples. These results were consistent with those of high-throughput sequencing.

CONCLUSIONS

ANK2p.F3067L is considered a novel and potentially pathogenic gene mutation in DCM.

TRAPPC9-related neurodevelopmental disorder: Report of a homozygous deletion in TRAPPC9 due to paternal uniparental isodisomy

TRAPPC9 loss-of-function biallelic variants are associated with an autosomal recessive intellectual disability syndrome (Online Mendelian Inheritance of Man no. 613192), also characterized by microcephaly, hypertelorism, obesity, growth delay, and behavioral differences. Here, we describe an 8-year-old Hispanic female with neurodevelopmental disorder, partial epilepsy, microcephaly, bilateral cleft lip and alveolus, growth delay, and dysmorphic features. She had abnormal myelination, mega cisterna magna, and colpocephaly on brain magnetic resonance imaging (MRI). Microarray showed a single ~146 Mb region of homozygosity (ROH) encompassing all of Chromosome 8, consistent with uniparental isodisomy (UPD). Exome sequencing performed in-house did not identify single nucleotide variants to explain her phenotype. Algorithms developed in-house and further evaluation of BAM files revealed a homozygous deletion overlapping Exon 2 in TRAPPC9 within the ROH. Subsequent del/dup analyses with exon-level oligo array confirmed a likely pathogenic deletion in TRAPPC9 (NM_031466.5): arr[GRCh37] 8q24.3(141460661_141461780)x0. Our case highlights the implications of downstream analyses from UPD/ROH given the increased risk for AR conditions, the strengths of combining orthologous molecular methods to establish a diagnosis and further delineates the TRAPPC9-related phenotype in an individual of Hispanic ancestry.

A biallelic frameshift indel in PPP1R35 as a cause of primary microcephaly

Protein phosphatase 1 regulatory subunit 35 (PPP1R35) encodes a centrosomal protein required for recruiting microtubule-binding elongation machinery. Several proteins in this centriole biogenesis pathway correspond to established primary microcephaly (MCPH) genes, and multiple model organism studies hypothesize PPP1R35 as a candidate MCPH gene. Here, using exome sequencing (ES) and family-based rare variant analyses, we report a homozygous, frameshifting indel deleting the canonical stop codon in the last exon of PPP1R35 [Chr7: c.753_*3delGGAAGCGTAGACCinsCG (p.Trp251Cysfs*22)]; the variant allele maps in a 3.7 Mb block of absence of heterozygosity (AOH) in a proband with severe MCPH (-4.3 SD at birth, -6.1 SD by 42 months), pachygyria, and global developmental delay from a consanguineous Turkish kindred. Droplet digital PCR (ddPCR) confirmed mutant mRNA expression in fibroblasts. In silico prediction of the translation of mutant PPP1R35 is expected to be elongated by 18 amino acids before encountering a downstream stop codon. This complex indel allele is absent in public databases (ClinVar, gnomAD, ARIC, 1000 genomes) and our in-house database of 14,000+ exomes including 1800+ Turkish exomes supporting predicted pathogenicity. Comprehensive literature searches for PPP1R35 variants yielded two probands affected with severe microcephaly (-15 SD and -12 SD) with the same homozygous indel from a single, consanguineous, Iranian family from a cohort of 404 predominantly Iranian families. The lack of heterozygous cases in two large cohorts representative of the genetic background of these two families decreased our suspicion of a founder allele and supports the contention of a recurrent mutation. We propose two potential secondary structure mutagenesis models for the origin of this variant allele mediated by hairpin formation between complementary GC rich segments flanking the stop codon via secondary structure mutagenesis.

Techniques for Profiling the Cellular Immune Response and Their Implications for Interventional Oncology

In recent years there has been increased interest in using the immune contexture of the primary tumors to predict the patient's prognosis. The tumor microenvironment of patients with cancers consists of different types of lymphocytes, tumor-infiltrating leukocytes, dendritic cells, and others. Different technologies can be used for the evaluation of the tumor microenvironment, all of which require a tissue or cell sample. Image-guided tissue sampling is a cornerstone in the diagnosis, stratification, and longitudinal evaluation of therapeutic efficacy for cancer patients receiving immunotherapies. Therefore, interventional radiologists (IRs) play an essential role in the evaluation of patients treated with systemically administered immunotherapies. This review provides a detailed description of different technologies used for immune assessment and analysis of the data collected from the use of these technologies. The detailed approach provided herein is intended to provide the reader with the knowledge necessary to not only interpret studies containing such data but also design and apply these tools for clinical practice and future research studies.

Clinical Study of 30 Novel KCNQ2 Variants/Deletions in KCNQ2-Related Disorders

Background

KCNQ2-related disorder is typically characterized as neonatal onset seizure and epileptic encephalopathy. The relationship between its phenotype and genotype is still elusive. This study aims to provide clinical features, management, and prognosis of patients with novel candidate variants of the KCNQ2 gene.

Methods

We enrolled patients with novel variants in the KCNQ2 gene from the China Neonatal Genomes Project between January 2018 and January 2021. All patients underwent next-generation sequencing tests and genetic data were analyzed by an in-house pipeline. The pathogenicity of variants was classified according to the guideline of the American College of Medical Genetics. Each case was evaluated by two geneticists back to back. Patients' information was acquired from clinical records.

Results

A total of 30 unrelated patients with novel variants in the KCNQ2 gene were identified, including 19 patients with single-nucleotide variants (SNVs) and 11 patients with copy number variants (CNVs). For the 19 SNVs, 12 missense variants and 7 truncating variants were identified. Of them, 36.8% (7/19) of the KCNQ2 variants were located in C-terminal regions, 15.7% (3/19) in segment S2, and 15.7% (3/19) in segment S4. Among them, 18 of 19 patients experienced seizures in the early neonatal period. However, one patient presented neurodevelopmental delay (NDD) as initial phenotype when he was 2 months old, and he had severe NDD when he was 3 years old. This patient did not present seizure but had abnormal electrographic background activity and brain imaging. Moreover, for the 11 patients with CNVs, 20q13.3 deletions involving EEF1A2, KCNQ2, and CHRNA4 genes were detected. All of them presented neonatal-onset seizures, responded to antiepileptic drugs, and had normal neurological development.

Conclusion

In this study, patients with novel KCNQ2 variants have variable phenotypes, whereas patients with 20q13.3 deletion involving EEF1A2, KCNQ2, and CHRNA4 genes tend to have normal neurological development.

Use of Next-Generation Sequencing for Identifying Mitochondrial Disorders

Mitochondria are major contributors to ATP synthesis, generating more than 90% of the total cellular energy production through oxidative phosphorylation (OXPHOS): metabolite oxidation, such as the β-oxidation of fatty acids, and the Krebs's cycle. OXPHOS inadequacy due to large genetic lesions in mitochondrial as well as nuclear genes and homo- or heteroplasmic point mutations in mitochondrially encoded genes is a characteristic of heterogeneous, maternally inherited genetic disorders known as mitochondrial disorders that affect multisystemic tissues and organs with high energy requirements, resulting in various signs and symptoms. Several traditional diagnostic approaches, including magnetic resonance imaging of the brain, cardiac testing, biochemical screening, variable heteroplasmy genetic testing, identifying clinical features, and skeletal muscle biopsies, are associated with increased risks, high costs, a high degree of false-positive or false-negative results, or a lack of precision, which limits their diagnostic abilities for mitochondrial disorders. Variable heteroplasmy levels, mtDNA depletion, and the identification of pathogenic variants can be detected through genetic sequencing, including the gold standard Sanger sequencing. However, sequencing can be time consuming, and Sanger sequencing can result in the missed recognition of larger structural variations such as CNVs or copy-number variations. Although each sequencing method has its own limitations, genetic sequencing can be an alternative to traditional diagnostic methods. The ever-growing roster of possible mutations has led to the development of next-generation sequencing (NGS). The enhancement of NGS methods can offer a precise diagnosis of the mitochondrial disorder within a short period at a reasonable expense for both research and clinical applications.

Novel pathogenic variants and quantitative phenotypic analyses of Robinow syndrome: WNT signaling perturbation and phenotypic variability

Robinow syndrome (RS) is a genetically heterogeneous disorder with six genes that converge on the WNT/planar cell polarity (PCP) signaling pathway implicated (DVL1, DVL3, FZD2, NXN, ROR2, and WNT5A). RS is characterized by skeletal dysplasia and distinctive facial and physical characteristics. To further explore the genetic heterogeneity, paralog contribution, and phenotypic variability of RS, we investigated a cohort of 22 individuals clinically diagnosed with RS from 18 unrelated families. Pathogenic or likely pathogenic variants in genes associated with RS or RS phenocopies were identified in all 22 individuals, including the first variant to be reported in DVL2. We retrospectively collected medical records of 16 individuals from this cohort and extracted clinical descriptions from 52 previously published cases. We performed Human Phenotype Ontology (HPO) based quantitative phenotypic analyses to dissect allele-specific phenotypic differences. Individuals with FZD2 variants clustered into two groups with demonstrable phenotypic differences between those with missense and truncating alleles. Probands with biallelic NXN variants clustered together with the majority of probands carrying DVL1, DVL2, and DVL3 variants, demonstrating no phenotypic distinction between the NXN-autosomal recessive and dominant forms of RS. While phenotypically similar diseases on the RS differential matched through HPO analysis, clustering using phenotype similarity score placed RS-associated phenotypes in a unique cluster containing WNT5A, FZD2, and ROR2 apart from non-RS-associated paralogs. Through human phenotype analyses of this RS cohort and OMIM clinical synopses of Mendelian disease, this study begins to tease apart specific biologic roles for non-canonical WNT-pathway proteins.

VarGenius-HZD Allows Accurate Detection of Rare Homozygous or Hemizygous Deletions in Targeted Sequencing Leveraging Breadth of Coverage

Homozygous deletions (HDs) may be the cause of rare diseases and cancer, and their discovery in targeted sequencing is a challenging task. Different tools have been developed to disentangle HD discovery but a sensitive caller is still lacking. We present VarGenius-HZD, a sensitive and scalable algorithm that leverages breadth-of-coverage for the detection of rare homozygous and hemizygous single-exon deletions (HDs). To assess its effectiveness, we detected both real and synthetic rare HDs in fifty exomes from the 1000 Genomes Project obtaining higher sensitivity in comparison with state-of-the-art algorithms that each missed at least one event. We then applied our tool on targeted sequencing data from patients with Inherited Retinal Dystrophies and solved five cases that still lacked a genetic diagnosis. We provide VarGenius-HZD either stand-alone or integrated within our recently developed software, enabling the automated selection of samples using the internal database. Hence, it could be extremely useful for both diagnostic and research purposes.

Targeting the Non-Coding Genome for the Diagnosis of Disorders of Sex Development

Disorders of sex development (DSD) are a complex group of conditions with highly variable clinical phenotypes, most often caused by failure of gonadal development. DSD are estimated to occur in around 1.7% of all live births. Whilst the understanding of genes involved in gonad development has increased exponentially, approximately 50% of patients with a DSD remain without a genetic diagnosis, possibly implicating non-coding genomic regions instead. Here, we review how variants in the non-coding genome of DSD patients can be identified using techniques such as array comparative genomic hybridization (CGH) to detect copy number variants (CNVs), and more recently, whole genome sequencing (WGS). Once a CNV in a patient's non-coding genome is identified, putative regulatory elements such as enhancers need to be determined within these vast genomic regions. We will review the available online tools and databases that can be used to refine regions with potential enhancer activity based on chromosomal accessibility, histone modifications, transcription factor binding site analysis, chromatin conformation, and disease association. We will also review the current in vitro and in vivo techniques available to demonstrate the functionality of the identified enhancers. The review concludes with a clinical update on the enhancers linked to DSD.

High prevalence of multilocus pathogenic variation in neurodevelopmental disorders in the Turkish population

Neurodevelopmental disorders (NDDs) are clinically and genetically heterogenous; many such disorders are secondary to perturbation in brain development and/or function. The prevalence of NDDs is > 3%, resulting in significant sociocultural and economic challenges to society. With recent advances in family-based genomics, rare-variant analyses, and further exploration of the Clan Genomics hypothesis, there has been a logarithmic explosion in neurogenetic "disease-associated genes" molecular etiology and biology of NDDs; however, the majority of NDDs remain molecularly undiagnosed. We applied genome-wide screening technologies, including exome sequencing (ES) and whole-genome sequencing (WGS), to identify the molecular etiology of 234 newly enrolled subjects and 20 previously unsolved Turkish NDD families. In 176 of the 234 studied families (75.2%), a plausible and genetically parsimonious molecular etiology was identified. Out of 176 solved families, deleterious variants were identified in 218 distinct genes, further documenting the enormous genetic heterogeneity and diverse perturbations in human biology underlying NDDs. We propose 86 candidate disease-trait-associated genes for an NDD phenotype. Importantly, on the basis of objective and internally established variant prioritization criteria, we identified 51 families (51/176 = 28.9%) with multilocus pathogenic variation (MPV), mostly driven by runs of homozygosity (ROHs) - reflecting genomic segments/haplotypes that are identical-by-descent. Furthermore, with the use of additional bioinformatic tools and expansion of ES to additional family members, we established a molecular diagnosis in 5 out of 20 families (25%) who remained undiagnosed in our previously studied NDD cohort emanating from Turkey.

Next-Generation Sequencing in the Field of Primary Immunodeficiencies: Current Yield, Challenges, and Future Perspectives

Primary immunodeficiencies comprise a group of inborn errors of immunity that display significant clinical and genetic heterogeneity. Next-generation sequencing techniques and predominantly whole exome sequencing have revolutionized the understanding of the genetic and molecular basis of genetic diseases, thereby also leading to a sharp increase in the discovery of new genes associated with primary immunodeficiencies. In this review, we discuss the current diagnostic yield of this generic diagnostic approach by evaluating the studies that have employed next-generation sequencing techniques in cohorts of patients with primary immunodeficiencies. The average diagnostic yield for primary immunodeficiencies is determined to be 29% (range 10-79%) and 38% specifically for whole-exome sequencing (range 15-70%). The significant variation between studies is mainly the result of differences in clinical characteristics of the studied cohorts but is also influenced by varying sequencing approaches and (in silico) gene panel selection. We further discuss other factors contributing to the relatively low yield, including the inherent limitations of whole-exome sequencing, challenges in the interpretation of novel candidate genetic variants, and promises of exploring the non-coding part of the genome. We propose strategies to improve the diagnostic yield leading the way towards expanded personalized treatment in PIDs.

Deep clinicopathological phenotyping identifies a previously unrecognized pathogenic EMD splice variant

Exome sequencing (ES) has revolutionized rare disease management, yet only ~25%-30% of patients receive a molecular diagnosis. A limiting factor is the quality of available phenotypic data. Here, we describe how deep clinicopathological phenotyping yielded a molecular diagnosis for a 19-year-old proband with muscular dystrophy and negative clinical ES. Deep phenotypic analysis identified two critical data points: (1) the absence of emerin protein in muscle biopsy and (2) clinical features consistent with Emery-Dreifuss muscular dystrophy. Sequencing data analysis uncovered an ultra-rare, intronic variant in EMD, the gene encoding emerin. The variant, NM_000117.3: c.188-6A > G, is predicted to impact splicing by in silico tools. This case thus illustrates how better integration of clinicopathologic data into ES analysis can enhance diagnostic yield with implications for clinical practice.

Herpes simplex encephalitis in a patient with a distinctive form of inherited IFNAR1 deficiency

Inborn errors of TLR3-dependent IFN-α/β- and IFN-λ-mediated immunity in the CNS can underlie herpes simplex virus 1 (HSV-1) encephalitis (HSE). The respective contributions of IFN-α/β and IFN-λ are unknown. We report a child homozygous for a genomic deletion of the entire coding sequence and part of the 3'-UTR of the last exon of IFNAR1, who died of HSE at the age of 2 years. An older cousin died following vaccination against measles, mumps, and rubella at 12 months of age, and another 17-year-old cousin homozygous for the same variant has had other, less severe, viral illnesses. The encoded IFNAR1 protein is expressed on the cell surface but is truncated and cannot interact with the tyrosine kinase TYK2. The patient's fibroblasts and EBV-B cells did not respond to IFN-α2b or IFN-β, in terms of STAT1, STAT2, and STAT3 phosphorylation or the genome-wide induction of IFN-stimulated genes. The patient's fibroblasts were susceptible to viruses, including HSV-1, even in the presence of exogenous IFN-α2b or IFN-β. HSE is therefore a consequence of inherited complete IFNAR1 deficiency. This viral disease occurred in natural conditions, unlike those previously reported in other patients with IFNAR1 or IFNAR2 deficiency. This experiment of nature indicates that IFN-α/β are essential for anti-HSV-1 immunity in the CNS.

Benchmarking germline CNV calling tools from exome sequencing data

Whole-exome sequencing is an attractive alternative to microarray analysis because of the low cost and potential ability to detect copy number variations (CNV) of various sizes (from 1-2 exons to several Mb). Previous comparison of the most popular CNV calling tools showed a high portion of false-positive calls. Moreover, due to a lack of a gold standard CNV set, the results are limited and incomparable. Here, we aimed to perform a comprehensive analysis of tools capable of germline CNV calling available at the moment using a single CNV standard and reference sample set. Compiling variants from previous studies with Bayesian estimation approach, we constructed an internal standard for NA12878 sample (pilot National Institute of Standards and Technology Reference Material) including 110,050 CNV or non-CNV exons. The standard was used to evaluate the performance of 16 germline CNV calling tools on the NA12878 sample and 10 correlated exomes as a reference set with respect to length distribution, concordance, and efficiency. Each algorithm had a certain range of detected lengths and showed low concordance with other tools. Most tools are focused on detection of a limited number of CNVs one to seven exons long with a false-positive rate below 50%. EXCAVATOR2, exomeCopy, and FishingCNV focused on detection of a wide range of variations but showed low precision. Upon unified comparison, the tools were not equivalent. The analysis performed allows choosing algorithms or ensembles of algorithms most suitable for a specific goal, e.g. population studies or medical genetics.

Detection of homozygous and hemizygous complete or partial exon deletions by whole-exome sequencing

The detection of copy number variations (CNVs) in whole-exome sequencing (WES) data is important, as CNVs may underlie a number of human genetic disorders. The recently developed HMZDelFinder algorithm can detect rare homozygous and hemizygous (HMZ) deletions in WES data more effectively than other widely used tools. Here, we present HMZDelFinder_opt, an approach that outperforms HMZDelFinder for the detection of HMZ deletions, including partial exon deletions in particular, in WES data from laboratory patient collections that were generated over time in different experimental conditions. We show that using an optimized reference control set of WES data, based on a PCA-derived Euclidean distance for coverage, strongly improves the detection of HMZ complete exon deletions both in real patients carrying validated disease-causing deletions and in simulated data. Furthermore, we develop a sliding window approach enabling HMZDelFinder_opt to identify HMZ partial deletions of exons that are undiscovered by HMZDelFinder. HMZDelFinder_opt is a timely and powerful approach for detecting HMZ deletions, particularly partial exon deletions, in WES data from inherently heterogeneous laboratory patient collections.

Genetic identification of pathogenic variations of the DMD gene: a retrospective study from 10,481 neonatal patients based on next-generation sequencing data

Background

An elevated level of creatine kinase (CK) is usually the primary screening marker for Duchenne muscular dystrophy (DMD)/Becker muscular dystrophy (BMD). This study investigated the clinical application of next-generation sequencing (NGS) in newborns with a possible diagnosis of DMD/BMD in the neonatal intensive care unit (NICU).

Methods

NGS data from the NICU between June 1, 2016, and June 30, 2020, were reanalyzed by an in-house pipeline. Other methods confirmed the genetic findings, and clinical follow-up was performed until August 1, 2020.

Results

Of the 10,481 newborns, 19 (0.18%, 19/10,481) cases with pathogenic variations of the DMD gene were identified, including 13 (68.4%, 13/19) deletions, 4 (21.1%, 4/19) duplications, and 2 (10.5%, 2/19) nonsense mutations. Eight of the cases were diagnosed with DMD. Therapeutic strategies were modified for these patients. Six cases were diagnosed with BMD. Five patients except for 1 deceased patient were further followed-up, and clinical management was adjusted based on the clinical symptoms. The remaining 5 cases were indeterminate for DMD and BMD. Genetic counseling and further follow-up were performed or suggested.

Conclusions

Our study showed that DMD/BMD could be diagnosed earlier in the neonatal stage before the typical clinical symptoms appear. Early diagnosis may provide an opportunity for guiding the care and treatment of patients. However, ethical issues need to be kept in mind in the process of genetic counseling.

Phenotype and genotype spectra of a Chinese cohort with nephronophthisis-related ciliopathy

BACKGROUND

Nephronophthisis-related ciliopathies (NPHP-RC) account for the majority of cases of monogenetically caused end-stage renal disease (ESRD) in children. Exploring the correlation between the phenotype and genotype of NPHP-RC is helpful for early diagnosis and management. We investigated the phenotype and genotype spectra of NPHP-RC in a Chinese multicentre cohort.

METHODS

Crosss-ectional and longitudinal data of 60 patients from 57 families with pathogenic NPHP-RC gene mutations distributed in 22 regions of China were collected into a unified, anonymous database. The mean observation time of this cohort was 3.5±3.1 years.

RESULTS

Mutations in NPHP1 and NPHP3 were the most common genetic defects. Overall, 45% of patients presented with isolated nephronophthisis (NPH), and 55% exhibited the extrarenal phenotype, which frequently involved the liver (41.7%, n=25), central nervous system (26.7%, n=16), eyes (26.7%, n=16) and skeletal system (11.7%, n=7). Accidental detection of elevated serum creatinine and non-specific symptoms caused by chronic kidney disease occurred in 65% of patients. Patients carrying NPHP1 mutations mainly presented with isolated NPH (90%, 18/20) and progressed to ESRD at a mean age of 12.9±0.5 years. The mean age of ESRD onset in the non-NPHP1 group was lower than that in the NPHP1 group (6.2±1.4 years, p<0.001), especially for patients carrying NPHP3 mutations (3.1±1.2 years), showing a heterogeneous phenotype characterised by Bardet-Biedl syndrome (12.5%, n=5), Joubert syndrome (7.5%, n=3), COACH syndrome (2.5%, n=1), Mainzer-Saldino syndrome (2.5%, n=1), short-rib thoracic dysplasia (2.5%, n=1) and unclassified symptoms (32.5%, n=13).

CONCLUSIONS

The Chinese Children Genetic Kidney Disease Database registry characterised the spectrum of the phenotype and genotype of NPHP-RC in the Chinese population. NPHP1 and NPHP3 were the most common pathogenic genes. Rapid progression to ESRD and liver involvement were noted in patients with NPHP3 mutations.

A flexible computational pipeline for research analyses of unsolved clinical exome cases

Exome sequencing has enabled molecular diagnoses for rare disease patients but often with initial diagnostic rates of ~25-30%. Here we develop a robust computational pipeline to rank variants for reassessment of unsolved rare disease patients. A comprehensive web-based patient report is generated in which all deleterious variants can be filtered by gene, variant characteristics, OMIM disease and Phenolyzer scores, and all are annotated with an ACMG classification and links to ClinVar. The pipeline ranked 21/34 previously diagnosed variants as top, with 26 in total ranked ≤7th, 3 ranked ≥13th; 5 failed the pipeline filters. Pathogenic/likely pathogenic variants by ACMG criteria were identified for 22/145 unsolved cases, and a previously undefined candidate disease variant for 27/145. This open access pipeline supports the partnership between clinical and research laboratories to improve the diagnosis of unsolved exomes. It provides a flexible framework for iterative developments to further improve diagnosis.

Next Generation Sequencing and Bioinformatics Analysis of Family Genetic Inheritance

Mendelian and complex genetic trait diseases continue to burden and affect society both socially and economically. The lack of effective tests has hampered diagnosis thus, the affected lack proper prognosis. Mendelian diseases are caused by genetic mutations in a singular gene while complex trait diseases are caused by the accumulation of mutations in either linked or unlinked genomic regions. Significant advances have been made in identifying novel diseases associated mutations especially with the introduction of next generation and third generation sequencing. Regardless, some diseases are still without diagnosis as most tests rely on SNP genotyping panels developed from population based genetic analyses. Analysis of family genetic inheritance using whole genomes, whole exomes or a panel of genes has been shown to be effective in identifying disease-causing mutations. In this review, we discuss next generation and third generation sequencing platforms, bioinformatic tools and genetic resources commonly used to analyze family based genomic data with a focus on identifying inherited or novel disease-causing mutations. Additionally, we also highlight the analytical, ethical and regulatory challenges associated with analyzing personal genomes which constitute the data used for family genetic inheritance.

Inborn errors of type I IFN immunity in patients with life-threatening COVID-19

Clinical outcome upon infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) ranges from silent infection to lethal coronavirus disease 2019 (COVID-19). We have found an enrichment in rare variants predicted to be loss-of-function (LOF) at the 13 human loci known to govern Toll-like receptor 3 (TLR3)- and interferon regulatory factor 7 (IRF7)-dependent type I interferon (IFN) immunity to influenza virus in 659 patients with life-threatening COVID-19 pneumonia relative to 534 subjects with asymptomatic or benign infection. By testing these and other rare variants at these 13 loci, we experimentally defined LOF variants underlying autosomal-recessive or autosomal-dominant deficiencies in 23 patients (3.5%) 17 to 77 years of age. We show that human fibroblasts with mutations affecting this circuit are vulnerable to SARS-CoV-2. Inborn errors of TLR3- and IRF7-dependent type I IFN immunity can underlie life-threatening COVID-19 pneumonia in patients with no prior severe infection.

Clinical genomics and contextualizing genome variation in the diagnostic laboratory

INTRODUCTION

The human genome contains the instructions for the development and biological homeostasis of the human organism and the genetic transmission of traits. Genome variation in human populations is the basis of evolution; individual or personal genomes vary tremendously, making each of us truly unique.

AREAS COVERED

Assaying this individual variation using genomic technologies has many applications in clinical medicine, from elucidating the biology of disease to designing strategies to ameliorate perturbations from homeostasis. Detecting pathogenic rare variation in a genome may provide a molecular diagnosis that can be informative for patient management and family healthcare.

EXPERT OPINION

Despite the increasing clinical use of unbiased genomic testing, including chromosome microarray analysis (CMA) with array comparative genomic hybridization (aCGH) or SNP arrays, clinical exome sequencing (cES), and whole-genome sequencing (WGS), to survey genome-wide for molecular aberrations, clinical acumen paired with an understanding of the limitations of each testing type will be needed to achieve molecular diagnoses. Potential opportunities for improving case solved rates, functionally annotating the majority of genes in the human genome, and further understanding genetic contributions to disease will empower clinical genomics and the precision medicine initiative.

The Impact of Next-Generation Sequencing on the Diagnosis, Treatment, and Prevention of Hereditary Neuromuscular Disorders

The impact of high-throughput sequencing in genetic neuromuscular disorders cannot be overstated. The ability to rapidly and affordably sequence multiple genes simultaneously has enabled a second golden age of Mendelian disease gene discovery, with flow-on impacts for rapid genetic diagnosis, evidence-based treatment, tailored therapy development, carrier-screening, and prevention of disease recurrence in families. However, there are likely many more neuromuscular disease genes and mechanisms to be discovered. Many patients and families remain without a molecular diagnosis following targeted panel sequencing, clinical exome sequencing, or even genome sequencing. Here we review how massively parallel, or next-generation, sequencing has changed the field of genetic neuromuscular disorders, and anticipate future benefits of recent technological innovations such as RNA-seq implementation and detection of tandem repeat expansions from short-read sequencing.

Clinical Utility of Next-Generation Sequencing for Developmental Disorders in the Rehabilitation Department: Experiences from a Single Chinese Center

This study investigated the clinical and genetic characteristics of developmental disorders (DDs) in children attending a rehabilitation department. A total of 94 children with suspected rare and undiagnosed DDs were included in this study. All patients were subjected to next-generation sequencing by means of proband single whole-exome sequencing (Pro-WES) or trio whole-exome sequencing (Trio-WES). To investigate the copy number variations (CNVs), 63 patients were subjected to the trio strategy, and 17 cases were subjected to the proband single strategy. The patients developed early and suffered from severe symptoms. WES reached a high diagnostic rate (48.7%, 46/94), and de novo (48.3%, 28/58) was the main pathogenic form. Most identified single-nucleotide variations (SNVs)/small insertions and deletions (indels) were found only in one patient. The number of uncertain significant locus in the patients taking Trio-WES was significantly lower than that in patients taking Pro-WES (2.1% vs 2.8%). Compared with hereditary mutations passed from parents, pathogenicity was more obvious in de novo mutations. The diagnostic rate of WES accompanied by CNVseq (57.5%, 46/80) was significantly higher (p = 0.016) than WES alone. Next-generation sequencing exhibited a satisfactory diagnostic rate for DDs patients in the rehabilitation department. Compared with the proband-only model, the family trio strategy should be employed more frequently because it can reduce the number of uncertain significant sites and help to identify de novo pathogenic mutations.

Integrated sequencing and array comparative genomic hybridization in familial Parkinson disease

Objective

To determine how single nucleotide variants (SNVs) and copy number variants (CNVs) contribute to molecular diagnosis in familial Parkinson disease (PD), we integrated exome sequencing (ES) and genome-wide array-based comparative genomic hybridization (aCGH) and further probed CNV structure to reveal mutational mechanisms.

Methods

We performed ES on 110 subjects with PD and a positive family history; 99 subjects were also evaluated using genome-wide aCGH. We interrogated ES and aCGH data for pathogenic SNVs and CNVs at Mendelian PD gene loci. We confirmed SNVs via Sanger sequencing and further characterized CNVs with custom-designed high-density aCGH, droplet digital PCR, and breakpoint sequencing.

Results

Using ES, we discovered individuals with known pathogenic SNVs in GBA (p.Glu365Lys, p.Thr408Met, p.Asn409Ser, and p.Leu483Pro) and LRRK2 (p.Arg1441Gly and p.Gly2019Ser). Two subjects were each double heterozygotes for variants in GBA and LRRK2. Based on aCGH, we additionally discovered cases with an SNCA duplication and heterozygous intragenic GBA deletion. Five additional subjects harbored both SNVs (p.Asn52Metfs*29, p.Thr240Met, p.Pro437Leu, and p.Trp453*) and likely disrupting CNVs at the PRKN locus, consistent with compound heterozygosity. In nearly all cases, breakpoint sequencing revealed microhomology, a mutational signature consistent with CNV formation due to DNA replication errors.

Conclusions

Integrated ES and aCGH yielded a genetic diagnosis in 19.3% of our familial PD cohort. Our analyses highlight potential mechanisms for SNCA and PRKN CNV formation, uncover multilocus pathogenic variation, and identify novel SNVs and CNVs for further investigation as potential PD risk alleles.

CNVs cause autosomal recessive genetic diseases with or without involvement of SNV/indels

PURPOSE

Improved resolution of molecular diagnostic technologies enabled detection of smaller sized exonic level copy-number variants (CNVs). The contribution of CNVs to autosomal recessive (AR) conditions may be better recognized using a large clinical cohort.

METHODS

We retrospectively investigated the CNVs' contribution to AR conditions in cases subjected to chromosomal microarray analysis (CMA, N = ~70,000) and/or clinical exome sequencing (ES, N = ~12,000) at Baylor Genetics; most had pediatric onset neurodevelopmental disorders.

RESULTS

CNVs contributed to biallelic variations in 87 cases, including 81 singletons and three affected sibling pairs. Seventy cases had CNVs affecting both alleles, and 17 had a CNV and a single-nucleotide variant (SNV)/indel in trans. In total, 94.3% of AR-CNVs affected one gene; among these 41.4% were single-exon and 35.0% were multiexon partial-gene events. Sixty-nine percent of homozygous AR-CNVs were embedded in homozygous genomic intervals. Five cases had large deletions unmasking an SNV/indel on the intact allele for a recessive condition, resulting in multiple molecular diagnoses.

CONCLUSIONS

AR-CNVs are often smaller in size, transmitted through generations, and underrecognized due to limitations in clinical CNV detection methods. Our findings from a large clinical cohort emphasized integrated CNV and SNV/indel analyses for precise clinical and molecular diagnosis especially in the context of genomic disorders.

Clinical exome sequencing as the first-tier test for diagnosing developmental disorders covering both CNV and SNV: a Chinese cohort

Background

Developmental disorders (DDs) are early onset disorders affecting 5%–10% of children worldwide. Chromosomal microarray analysis detecting CNVs is currently recommended as the first-tier test for DD diagnosis. However, this analysis omits a high percentage of disease-causing single nucleotide variations (SNVs) that warrant further sequencing. Currently, next-generation sequencing can be used in clinical scenarios detecting CNVs, and the use of exome sequencing in the DD cohort ahead of the microarray test has not been evaluated.

Methods

Clinical exome sequencing (CES) was performed on 1090 unrelated Chinese DD patients who were classified into five phenotype subgroups. CNVs and SNVs were both detected and analysed based on sequencing data.

Results

An overall diagnostic rate of 41.38% was achieved with the combinational analysis of CNV and SNV. Over 12.02% of patients were diagnosed based on CNV, which was comparable with the published CMA diagnostic rate, while 0.74% were traditionally elusive cases who had dual diagnosis or apparently homozygous mutations that were clarified. The diagnostic rates among subgroups ranged from 21.82% to 50.32%. The top three recurrent cytobands with diagnostic CNVs were 15q11.2-q13.1, 22q11.21 and 7q11.23. The top three genes with diagnostic SNVs were: MECP2, SCN1A and SCN2A. Both the diagnostic rate and spectrums of CNVs and SNVs showed differences among the phenotype subgroups.

Conclusion

With a higher diagnostic rate, more comprehensive observation of variations and lower cost compared with conventional strategies, simultaneous analysis of CNVs and SNVs based on CES showed potential as a new first-tier choice to diagnose DD.

Role of carotenoids and retinoids during heart development

The nutritional requirements of the developing embryo are complex. In the case of dietary vitamin A (retinol, retinyl esters and provitamin A carotenoids), maternal derived nutrients serve as precursors to signaling molecules such as retinoic acid, which is required for embryonic patterning and organogenesis. Despite variations in the composition and levels of maternal vitamin A, embryonic tissues need to generate a precise amount of retinoic acid to avoid congenital malformations. Here, we summarize recent findings regarding the role and metabolism of vitamin A during heart development and we survey the association of genes known to affect retinoid metabolism or signaling with various inherited disorders. A better understanding of the roles of vitamin A in the heart and of the factors that affect retinoid metabolism and signaling can help design strategies to meet nutritional needs and to prevent birth defects and disorders associated with altered retinoid metabolism. This article is part of a Special Issue entitled Carotenoids recent advances in cell and molecular biology edited by Johannes von Lintig and Loredana Quadro.

Sialidosis type II: Expansion of phenotypic spectrum and identification of a common mutation in seven patients

Sialidosis, an autosomal recessive disorder, is characterized by progressive lysosomal storage of sialylated glycopeptides and oligosaccharides. It occurs as a result of biallelic mutations in the NEU1 gene. Sialidosis is traditionally classified as a milder, late-onset type I and a severe early-onset type II disease. The presence of a cherry-red spot is a well-established ophthalmological clue to the disorder. We present a clinical-radiological report of seven unrelated patients with molecularly confirmed sialidosis type II. To the best of our knowledge, This is the largest reported series of patients with Sialidosis type II. A novel, previously unreported ophthalmic phenotype of bulls-eye maculopathy, is described. All seven phenotypically heterogeneous patients had the same pathogenic variant (c.679G > A; p.Gly227Arg) at a homozygous level in the NEU1 gene. We propose that this is a common mutation in north Indians for this rare disorder. We also observed an overlap of symptoms and a continuum of phenotypes in type I and II Sialidosis.

Deficiencies in vesicular transport mediated by TRAPPC4 are associated with severe syndromic intellectual disability

The conserved transport protein particle (TRAPP) complexes regulate key trafficking events and are required for autophagy. TRAPPC4, like its yeast Trs23 orthologue, is a core component of the TRAPP complexes and one of the essential subunits for guanine nucleotide exchange factor activity for Rab1 GTPase. Pathogenic variants in specific TRAPP subunits are associated with neurological disorders. We undertook exome sequencing in three unrelated families of Caucasian, Turkish and French-Canadian ethnicities with seven affected children that showed features of early-onset seizures, developmental delay, microcephaly, sensorineural deafness, spastic quadriparesis and progressive cortical and cerebellar atrophy in an effort to determine the genetic aetiology underlying neurodevelopmental disorders. All seven affected subjects shared the same identical rare, homozygous, potentially pathogenic variant in a non-canonical, well-conserved splice site within TRAPPC4 (hg19:chr11:g.118890966A>G; TRAPPC4: NM_016146.5; c.454+3A>G). Single nucleotide polymorphism array analysis revealed there was no haplotype shared between the tested Turkish and Caucasian families suggestive of a variant hotspot region rather than a founder effect. In silico analysis predicted the variant to cause aberrant splicing. Consistent with this, experimental evidence showed both a reduction in full-length transcript levels and an increase in levels of a shorter transcript missing exon 3, suggestive of an incompletely penetrant splice defect. TRAPPC4 protein levels were significantly reduced whilst levels of other TRAPP complex subunits remained unaffected. Native polyacrylamide gel electrophoresis and size exclusion chromatography demonstrated a defect in TRAPP complex assembly and/or stability. Intracellular trafficking through the Golgi using the marker protein VSVG-GFP-ts045 demonstrated significantly delayed entry into and exit from the Golgi in fibroblasts derived from one of the affected subjects. Lentiviral expression of wild-type TRAPPC4 in these fibroblasts restored trafficking, suggesting that the trafficking defect was due to reduced TRAPPC4 levels. Consistent with the recent association of the TRAPP complex with autophagy, we found that the fibroblasts had a basal autophagy defect and a delay in autophagic flux, possibly due to unsealed autophagosomes. These results were validated using a yeast trs23 temperature sensitive variant that exhibits constitutive and stress-induced autophagic defects at permissive temperature and a secretory defect at restrictive temperature. In summary we provide strong evidence for pathogenicity of this variant in a member of the core TRAPP subunit, TRAPPC4 that associates with vesicular trafficking and autophagy defects. This is the first report of a TRAPPC4 variant, and our findings add to the growing number of TRAPP-associated neurological disorders.

A diagnostic ceiling for exome sequencing in cerebellar ataxia and related neurological disorders

Genetic ataxias are associated with mutations in hundreds of genes with high phenotypic overlap complicating the clinical diagnosis. Whole-exome sequencing (WES) has increased the overall diagnostic rate considerably. However, the upper limit of this method remains ill-defined, hindering efforts to address the remaining diagnostic gap. To further assess the role of rare coding variation in ataxic disorders, we reanalyzed our previously published exome cohort of 76 predominantly adult and sporadic-onset patients, expanded the total number of cases to 260, and introduced analyses for copy number variation and repeat expansion in a representative subset. For new cases (n = 184), our resulting clinically relevant detection rate remained stable at 47% with 24% classified as pathogenic. Reanalysis of the previously sequenced 76 patients modestly improved the pathogenic rate by 7%. For the combined cohort (n = 260), the total observed clinical detection rate was 52% with 25% classified as pathogenic. Published studies of similar neurological phenotypes report comparable rates. This consistency across multiple cohorts suggests that, despite continued technical and analytical advancements, an approximately 50% diagnostic rate marks a relative ceiling for current WES-based methods and a more comprehensive genome-wide assessment is needed to identify the missing causative genetic etiologies for cerebellar ataxia and related neurodegenerative diseases.

Comprehensive Outline of Whole Exome Sequencing Data Analysis Tools Available in Clinical Oncology

Whole exome sequencing (WES) enables the analysis of all protein coding sequences in the human genome. This technology enables the investigation of cancer-related genetic aberrations that are predominantly located in the exonic regions. WES delivers high-throughput results at a reasonable price. Here, we review analysis tools enabling utilization of WES data in clinical and research settings. Technically, WES initially allows the detection of single nucleotide variants (SNVs) and copy number variations (CNVs), and data obtained through these methods can be combined and further utilized. Variant calling algorithms for SNVs range from standalone tools to machine learning-based combined pipelines. Tools for CNV detection compare the number of reads aligned to a dedicated segment. Both SNVs and CNVs help to identify mutations resulting in pharmacologically druggable alterations. The identification of homologous recombination deficiency enables the use of PARP inhibitors. Determining microsatellite instability and tumor mutation burden helps to select patients eligible for immunotherapy. To pave the way for clinical applications, we have to recognize some limitations of WES, including its restricted ability to detect CNVs, low coverage compared to targeted sequencing, and the missing consensus regarding references and minimal application requirements. Recently, Galaxy became the leading platform in non-command line-based WES data processing. The maturation of next-generation sequencing is reinforced by Food and Drug Administration (FDA)-approved methods for cancer screening, detection, and follow-up. WES is on the verge of becoming an affordable and sufficiently evolved technology for everyday clinical use.

A Patient with Berardinelli-Seip Syndrome, Novel AGPAT2 Splicesite Mutation and Concomitant Development of Non-diabetic Polyneuropathy

Primary polyneuropathy in the context of Seip-Berardinelli type 1 seipinopathy, or congenital generalized lipodystrophy type 1 (CGL1) has not been previously reported. We report the case history of a 27 year old female CGL1 patient presenting with an unusual additional development of non-diabetic peripheral neuropathy and learning disabilities in early adolescence. Whole exome sequencing (WES) of the patient genome identified a novel variant, homozygous for a 52 bp intronic deletion in the AGPAT2 locus, coding for 1-acylglycerol-3-phosphate O-acyltransferase 2, which is uniquely associated with CGL1 seipinopathies, with no molecular evidence for dual diagnosis. Functional studies using RNA isolated from patient peripheral blood leucocytes showed abnormal RNA splicing resulting in the loss of 25 amino acids from the patient AGPAT2 protein coding sequence. Stability and transcription levels for the misspliced AGPAT2 mRNA in our patient nonetheless remained normal. Any AGPAT2 protein produced in our patient is therefore likely to be dysfunctional. However, formal linkage of this deletion to the neuropathy observed remains to be shown. The classical clinical presentation of a patient with AGPAT2-associated lipodystrophy shows normal cognition and no development of polyneuropathy. Cognitive disabilities and polyneuropathy are features associated exclusively with clinical CGL type 2 arising from seipin (BSCL2) gene mutations. This case study suggests that in some genetic contexts, AGPAT2 mutations can also produce phenotypes with primary polyneuropathy.

SeQuiLa-cov: A fast and scalable library for depth of coverage calculations

BACKGROUND

Depth of coverage calculation is an important and computationally intensive preprocessing step in a variety of next-generation sequencing pipelines, including the analysis of RNA-sequencing data, detection of copy number variants, or quality control procedures.

RESULTS

Building upon big data technologies, we have developed SeQuiLa-cov, an extension to the recently released SeQuiLa platform, which provides efficient depth of coverage calculations, reaching >100× speedup over the state-of-the-art tools. The performance and scalability of our solution allow for exome and genome-wide calculations running locally or on a cluster while hiding the complexity of the distributed computing with Structured Query Language Application Programming Interface.

CONCLUSIONS

SeQuiLa-cov provides significant performance gain in depth of coverage calculations streamlining the widely used bioinformatic processing pipelines.

Exome Sequencing of a Primary Ovarian Insufficiency Cohort Reveals Common Molecular Etiologies for a Spectrum of Disease

CONTEXT

Primary ovarian insufficiency (POI) encompasses a spectrum of premature menopause, including both primary and secondary amenorrhea. For 75% to 90% of individuals with hypergonadotropic hypogonadism presenting as POI, the molecular etiology is unknown. Common etiologies include chromosomal abnormalities, environmental factors, and congenital disorders affecting ovarian development and function, as well as syndromic and nonsyndromic single gene disorders suggesting POI represents a complex trait.

OBJECTIVE

To characterize the contribution of known disease genes to POI and identify molecular etiologies and biological underpinnings of POI.

DESIGN, SETTING, AND PARTICIPANTS

We applied exome sequencing (ES) and family-based genomics to 42 affected female individuals from 36 unrelated Turkish families, including 31 with reported parental consanguinity.

RESULTS

This analysis identified likely damaging, potentially contributing variants and molecular diagnoses in 16 families (44%), including 11 families with likely damaging variants in known genes and five families with predicted deleterious variants in disease genes (IGSF10, MND1, MRPS22, and SOHLH1) not previously associated with POI. Of the 16 families, 2 (13%) had evidence for potentially pathogenic variants at more than one locus. Absence of heterozygosity consistent with identity-by-descent mediated recessive disease burden contributes to molecular diagnosis in 15 of 16 (94%) families. GeneMatcher allowed identification of additional families from diverse genetic backgrounds.

CONCLUSIONS

ES analysis of a POI cohort further characterized locus heterogeneity, reaffirmed the association of genes integral to meiotic recombination, demonstrated the likely contribution of genes involved in hypothalamic development, and documented multilocus pathogenic variation suggesting the potential for oligogenic inheritance contributing to the development of POI.

Paralog Studies Augment Gene Discovery: DDX and DHX Genes

Members of a paralogous gene family in which variation in one gene is known to cause disease are eight times more likely to also be associated with human disease. Recent studies have elucidated DHX30 and DDX3X as genes for which pathogenic variant alleles are involved in neurodevelopmental disorders. We hypothesized that variants in paralogous genes encoding members of the DExD/H-box RNA helicase superfamily might also underlie developmental delay and/or intellectual disability (DD and/or ID) disease phenotypes. Here we describe 15 unrelated individuals who have DD and/or ID, central nervous system (CNS) dysfunction, vertebral anomalies, and dysmorphic features and were found to have probably damaging variants in DExD/H-box RNA helicase genes. In addition, these individuals exhibit a variety of other tissue and organ system involvement including ocular, outer ear, hearing, cardiac, and kidney tissues. Five individuals with homozygous (one), compound-heterozygous (two), or de novo (two) missense variants in DHX37 were identified by exome sequencing. We identified ten total individuals with missense variants in three other DDX/DHX paralogs: DHX16 (four individuals), DDX54 (three individuals), and DHX34 (three individuals). Most identified variants are rare, predicted to be damaging, and occur at conserved amino acid residues. Taken together, these 15 individuals implicate the DExD/H-box helicases in both dominantly and recessively inherited neurodevelopmental phenotypes and highlight the potential for more than one disease mechanism underlying these disorders.

The Genomics of Arthrogryposis, a Complex Trait: Candidate Genes and Further Evidence for Oligogenic Inheritance

Arthrogryposis is a clinical finding that is present either as a feature of a neuromuscular condition or as part of a systemic disease in over 400 Mendelian conditions. The underlying molecular etiology remains largely unknown because of genetic and phenotypic heterogeneity. We applied exome sequencing (ES) in a cohort of 89 families with the clinical sign of arthrogryposis. Additional molecular techniques including array comparative genomic hybridization (aCGH) and Droplet Digital PCR (ddPCR) were performed on individuals who were found to have pathogenic copy number variants (CNVs) and mosaicism, respectively. A molecular diagnosis was established in 65.2% (58/89) of families. Eleven out of 58 families (19.0%) showed evidence for potential involvement of pathogenic variation at more than one locus, probably driven by absence of heterozygosity (AOH) burden due to identity-by-descent (IBD). RYR3, MYOM2, ERGIC1, SPTBN4, and ABCA7 represent genes, identified in two or more families, for which mutations are probably causative for arthrogryposis. We also provide evidence for the involvement of CNVs in the etiology of arthrogryposis and for the idea that both mono-allelic and bi-allelic variants in the same gene cause either similar or distinct syndromes. We were able to identify the molecular etiology in nine out of 20 families who underwent reanalysis. In summary, our data from family-based ES further delineate the molecular etiology of arthrogryposis, yielded several candidate disease-associated genes, and provide evidence for mutational burden in a biological pathway or network. Our study also highlights the importance of reanalysis of individuals with unsolved diagnoses in conjunction with sequencing extended family members.

Genome sequencing and implications for rare disorders

The practice of genomic medicine stands to revolutionize our approach to medical care, and to realize this goal will require discovery of the relationship between rare variation at each of the ~ 20,000 protein-coding genes and their consequent impact on individual health and expression of Mendelian disease. The step-wise evolution of broad-based, genome-wide cytogenetic and molecular genomic testing approaches (karyotyping, chromosomal microarray [CMA], exome sequencing [ES]) has driven much of the rare disease discovery to this point, with genome sequencing representing the newest member of this team. Each step has brought increased sensitivity to interrogate individual genomic variation in an unbiased method that does not require clinical prediction of the locus or loci involved. Notably, each step has also brought unique limitations in variant detection, for example, the low sensitivity of ES for detection of triploidy, and of CMA for detection of copy neutral structural variants. The utility of genome sequencing (GS) as a clinical molecular diagnostic test, and the increased sensitivity afforded by addition of long-read sequencing or other -omics technologies such as RNAseq or metabolomics, are not yet fully explored, though recent work supports improved sensitivity of variant detection, at least in a subset of cases. The utility of GS will also rely upon further elucidation of the complexities of genetic and allelic heterogeneity, multilocus rare variation, and the impact of rare and common variation at a locus, as well as advances in functional annotation of identified variants. Much discovery remains to be done before the potential utility of GS is fully appreciated.

Exome sequencing in routine diagnostics: a generic test for 254 patients with primary immunodeficiencies

BACKGROUND

Diagnosis of primary immunodeficiencies (PIDs) is complex and cumbersome yet important for the clinical management of the disease. Exome sequencing may provide a genetic diagnosis in a significant number of patients in a single genetic test.

METHODS

In May 2013, we implemented exome sequencing in routine diagnostics for patients suffering from PIDs. This study reports the clinical utility and diagnostic yield for a heterogeneous group of 254 consecutively referred PID patients from 249 families. For the majority of patients, the clinical diagnosis was based on clinical criteria including rare and/or unusual severe bacterial, viral, or fungal infections, sometimes accompanied by autoimmune manifestations. Functional immune defects were interpreted in the context of aberrant immune cell populations, aberrant antibody levels, or combinations of these factors.

RESULTS

For 62 patients (24%), exome sequencing identified pathogenic variants in well-established PID genes. An exome-wide analysis diagnosed 10 additional patients (4%), providing diagnoses for 72 patients (28%) from 68 families altogether. The genetic diagnosis directly indicated novel treatment options for 25 patients that received a diagnosis (34%).

CONCLUSION

Exome sequencing as a first-tier test for PIDs granted a diagnosis for 28% of patients. Importantly, molecularly defined diagnoses indicated altered therapeutic options in 34% of cases. In addition, exome sequencing harbors advantages over gene panels as a truly generic test for all genetic diseases, including in silico extension of existing gene lists and re-analysis of existing data.

Diagnostic Yields of Trio-WES Accompanied by CNVseq for Rare Neurodevelopmental Disorders

Objective

This study is to investigate the diagnostic yield of the combination of trio whole exome sequencing (Trio-WES) and copy number variation sequencing (CNVseq) for rare neurodevelopmental disorders (NDDs).

Methods

Clinical data from consecutive pediatric patients who were diagnosed with rare NDDs that were suspected to be monogenic disorders, who were admitted to our hospital from April 2017 to March 2019, and who underwent next generation sequencing (NGS) were extracted from the medical records. Patients for whom Trio-WES and CNVseq data were available were enrolled in this study. Sanger sequencing was applied for the validation of the variants identified by Trio-WES. Sequence alignment and structural modeling were conducted for analyzing the possibility of the variants in the onset of the NDDs.

Results

In total, 54 patients were enrolled in this study, with the median age of 15 (8–26) months. A total of 242 phenotypic abnormalities belonging to 20 different systems were identified in the cohort. Twenty-four patients were diagnosed by Trio-WES, eight patients were diagnosed by CNVseq, and one case was identified by both WES and CNVseq. Compared with Trio-WES, the diagnosis rate of Trio-WES accompanied by CNVseq was significantly higher (P = 0.016). Trio-WES identified 36 variants in 26 different genes, among which 27 variants were novel. CNVseq detected four duplications and eight deletions, ranging from 310 kb to 23.27 Mb. Our case examples demonstrated the high heterogeneity of NDDs and showed the challenges of rare NDDs for physicians.

Conclusion

The significantly higher diagnosis rate of Trio-WES accompanied by CNVseq makes this strategy a potential alternative to the most widely used approaches for pediatric children with rare and undiagnosed NDDs.