Transcriptome-directed analysis for Mendelian disease diagnosis overcomes limitations of conventional genomic testing

Rapid genome sequencing for pediatrics

The advancements made in next-generation sequencing (NGS) technology over the past two decades have transformed our understanding of genetic variation in humans and had a profound impact on our ability to diagnose patients with rare genetic diseases. In this review, we discuss the recently developed application of rapid NGS techniques, used to diagnose pediatric patients with suspected rare diseases who are critically ill. We highlight the challenges associated with performing such clinical diagnostics tests in terms of the laboratory infrastructure, bioinformatic analysis pipelines, and the ethical considerations that need to be addressed. We end by looking at what future developments in this field may look like and how they can be used to augment the genetic data to further improve the diagnostic rates for these high-priority patients.

Transcriptome analysis provides critical answers to the "variants of uncertain significance" conundrum

While whole-genome and exome sequencing have transformed our collective understanding of genetics' role in disease pathogenesis, there are certain conditions and populations for whom DNA-level data fails to identify the underlying genetic etiology. Specifically, patients of non-White race and non-European ancestry are disproportionately affected by "variants of unknown/uncertain significance" (VUS), limiting the scope of precision medicine for minority patients and perpetuating health disparities. VUS often include deep intronic and splicing variants which are difficult to interpret from DNA data alone. RNA analysis can illuminate the consequences of VUS, thereby allowing for their reclassification as pathogenic versus benign. Here we review the critical role transcriptome analysis plays in clarifying VUS in both neoplastic and non-neoplastic diseases.

OMIXCARE: OMICS technologies solved about 33% of the patients with heterogeneous rare neuro-developmental disorders and negative exome sequencing results and identified 13% additional candidate variants

Purpose: Patients with rare or ultra-rare genetic diseases, which affect 350 million people worldwide, may experience a diagnostic odyssey. High-throughput sequencing leads to an etiological diagnosis in up to 50% of individuals with heterogeneous neurodevelopmental or malformation disorders. There is a growing interest in additional omics technologies in translational research settings to examine the remaining unsolved cases. Methods: We gathered 30 individuals with malformation syndromes and/or severe neurodevelopmental disorders with negative trio exome sequencing and array comparative genomic hybridization results through a multicenter project. We applied short-read genome sequencing, total RNA sequencing, and DNA methylation analysis, in that order, as complementary translational research tools for a molecular diagnosis. Results: The cohort was mainly composed of pediatric individuals with a median age of 13.7 years (4 years and 6 months to 35 years and 1 month). Genome sequencing alone identified at least one variant with a high level of evidence of pathogenicity in 8/30 individuals (26.7%) and at least a candidate disease-causing variant in 7/30 other individuals (23.3%). RNA-seq data in 23 individuals allowed two additional individuals (8.7%) to be diagnosed, confirming the implication of two pathogenic variants (8.7%), and excluding one candidate variant (4.3%). Finally, DNA methylation analysis confirmed one diagnosis identified by genome sequencing (Kabuki syndrome) and identified an episignature compatible with a BAFopathy in a patient with a clinical diagnosis of Coffin-Siris with negative genome and RNA-seq results in blood. Conclusion: Overall, our integrated genome, transcriptome, and DNA methylation analysis solved 10/30 (33.3%) cases and identified a strong candidate gene in 4/30 (13.3%) of the patients with rare neurodevelopmental disorders and negative exome sequencing results.

Can tandem alternative splicing and evasion of premature termination codon surveillance contribute to attenuated Peutz-Jeghers syndrome?

Alternative use of short distance tandem sites such as NAGN_n AG are a common mechanism of alternative splicing; however, single nucleotide variants are rarely reported as likely to generate or to disrupt tandem splice sites. We identify a pathogenic intron 5 STK11 variant (NM_000455.4:c.[735-6A>G];[=]) segregating with the mucocutaneous features but not the hamartomatous polyps of Peutz-Jeghers syndrome in two individuals. By RNAseq analysis of peripheral blood mRNA, this variant was shown to generate a novel and preferentially used tandem proximal splice acceptor (AAGTGAAG). The variant transcript (NM_000455.4:c.734_734 + 1insTGAAG), which encodes a frameshift (p.[Tyr246Glufs*43]) constituted 36%-43% of STK11 transcripts suggesting partial escape from nonsense mediated mRNA decay and translation of a truncated protein. A review of the ClinVar database identified other similar variants. We suggest that nucleotide changes creating or disrupting tandem alternative splice sites are a pertinent disease mechanism and require contextualization for clinical reporting. Additionally, we hypothesize that some pathogenic STK11 variants cause an attenuated phenotype.

Epigenomic Approaches for the Diagnosis of Rare Diseases

Rare diseases affect more than 300 million people worldwide. Diagnosing rare diseases is a major challenge as they have different causes and etiologies. Careful assessment of clinical symptoms often leads to the testing of the most common genetic alterations that could explain the disease. Patients with negative results for these tests frequently undergo whole exome or genome sequencing, leading to the identification of the molecular cause of the disease in 50% of patients at best. Therefore, a significant proportion of patients remain undiagnosed after sequencing their genome. Recently, approaches based on functional aspects of the genome, including transcriptomics and epigenomics, are beginning to emerge. Here, we will review these approaches, including studies that have successfully provided diagnoses for complex undiagnosed cases.

The Analysis of the Human Megakaryocyte and Platelet Coding Transcriptome in Healthy and Diseased Subjects

Platelets are generated and released into the bloodstream from their precursor cells, megakaryocytes that reside in the bone marrow. Though platelets have no nucleus or DNA, they contain a full transcriptome that, during platelet formation, is transported from the megakaryocyte to the platelet. It has been described that transcripts in platelets can be translated into proteins that influence platelet response. The platelet transcriptome is highly dynamic and has been extensively studied using microarrays and, more recently, RNA sequencing (RNA-seq) in relation to diverse conditions (inflammation, obesity, cancer, pathogens and others). In this review, we focus on bulk and single-cell RNA-seq studies that have aimed to characterize the coding transcriptome of healthy megakaryocytes and platelets in humans. It has been noted that bulk RNA-seq has limitations when studying in vitro-generated megakaryocyte cultures that are highly heterogeneous, while single-cell RNA-seq has not yet been applied to platelets due to their very limited RNA content. Next, we illustrate how these methods can be applied in the field of inherited platelet disorders for gene discovery and for unraveling novel disease mechanisms using RNA from platelets and megakaryocytes and rare disease bioinformatics. Next, future perspectives are discussed on how this field of coding transcriptomics can be integrated with other next-generation technologies to decipher unexplained inherited platelet disorders in a multiomics approach.

Advances in studying human gametogenesis and embryonic development in China

Reproductive medicine in China has developed rapidly since 1988 due to the support from the government and scientific exploration. However, the success rate of assisted reproduction technology (ART) is around 30-40% and many unknown "black boxes" in gametogenesis and embryo development are still present. With the development of single-cell and low-input sequencing technologies, the network of transcriptome and epigenetic regulation (DNA methylation, chromatin accessibility, and histone modifications) during the development of human primordial germ cells (PGCs), gametes and embryos has been investigated in depth. Furthermore, pre-implantation genetic testing (PGT) has also rapidly developed. In this review, we summarize and analyze China's outstanding progress in these fields.

Short amplicon reverse transcription-polymerase chain reaction detects aberrant splicing in genes with low expression in blood missed by ribonucleic acid sequencing analysis for clinical diagnosis

Use of blood RNA sequencing (RNA-seq) as a splicing analysis tool for clinical interpretation of variants of uncertain significance (VUSs) found via whole-genome and exome sequencing can be difficult for genes that have low expression in the blood due to insufficient read count coverage aligned to specific genes of interest. Here, we present a short amplicon reverse transcription-polymerase chain reaction(RT-PCR) for the detection of genes with low blood expression. Short amplicon RT-PCR, is designed to span three exons where an exon harboring a variant is flanked by one upstream and one downstream exon. We tested short amplicon RT-PCRs for genes that have median transcripts per million (TPM) values less than one according to the genotype-tissue expression database. Median TPM values of genes analyzed in this study are SYN1 = 0.8549, COL1A1 = 0.6275, TCF4 = 0.4009, DSP = .2894, TTN = 0.2851, COL5A2 = 0.1036, TERT = 0.04452, NTRK2 = 0.0344, ABCA4 = 0.00744, PRPH = 0, and WT1 = 0. All these genes show insufficient exon-spanning read coverage in our RNA-seq data to allow splicing analysis. We successfully detected all genes tested except PRPH and WT1. Aberrant splicing was detected in SYN1, TCF4, NTRK2, TTN, and TERT VUSs. Therefore, our results show short amplicon RT-PCR is a useful alternative for the analysis of splicing events in genes with low TPM in blood RNA for clinical diagnostics.

A novel, de novo intronic variant in POGZ causes White-Sutton syndrome

White-Sutton syndrome (WHSUS), which is caused by heterozygous pathogenic variants in POGZ, is characterized by a spectrum of intellectual disabilities and global developmental delay with or without features of autism spectrum disorder. Additional features may include hypotonia, behavioral abnormalities, ophthalmic abnormalities, hearing loss, sleep apnea, microcephaly, dysmorphic facial features, and rarely, congenital diaphragmatic hernia (CDH). We present a 6-year-old female with features of WHSUS, including CDH, but with nondiagnostic clinical trio exome sequencing. Exome sequencing reanalysis revealed a heterozygous, de novo, intronic variant in POGZ (NM_015100.3:c.2546-20T>A). RNA sequencing revealed that this intronic variant leads to skipping of exon 18. This exon skipping event results in a frameshift with a predicted premature stop codon in the last exon and escape from nonsense-mediated mRNA decay (NMD). To our knowledge, this case is the first case of WHSUS caused by a de novo, intronic variant that is not near a canonical splice site within POGZ. These findings emphasize the limitations of standard clinical exome filtering algorithms and the importance of research reanalysis of exome data together with RNA sequencing to confirm a suspected diagnosis of WHSUS. As the sixth reported case of CDH with heterozygous pathogenic variants in POGZ and features consistent with WHSUS, this report supports the conclusion that WHSUS should be considered in the differential diagnosis for patients with syndromic CDH.

PRUNE1 c.933G>A synonymous variant induces exon 7 skipping, disrupts the DHHA2 domain, and leads to an atypical NMIHBA syndrome presentation: Case report and review of the literature

Prune exopolyphosphatase-1 (PRUNE1) encodes a member of the aspartic acid-histidine-histidine (DHH) phosphodiesterase superfamily that regulates cell migration and proliferation during brain development. In 2015, biallelic PRUNE1 loss-of-function variants were identified to cause the neurodevelopmental disorder with microcephaly, hypotonia, and variable brain abnormalities (NMIHBA, OMIM#617481). NMIHBA is characterized by the namesake features and structural brain anomalies including thinning of the corpus callosum, cerebral and cerebellar atrophy, and delayed myelination. To date, 47 individuals have been reported in the literature, but the phenotypic spectrum of PRUNE1-related disorders and their causative variants remains to be characterized fully. Here, we report a novel homozygous PRUNE1 NM_021222.2:c.933G>A synonymous variant identified in a 6-year-old boy with intellectual and developmental disabilities, hypotonia, and spastic diplegia, but with the absence of microcephaly, brain anomalies, or seizures. Fibroblast RNA sequencing revealed that the PRUNE1 NM_021222.1:c.933G>A variant resulted in an in-frame skipping of the penultimate exon 7, removing 53 amino acids from an important protein domain. This case represents the first synonymous variant and the third pathogenic variant known to date affecting the DHH-associated domain (DHHA2 domain). These findings extend the genotypic and phenotypic spectrums in PRUNE1-related disorders and highlight the importance of considering synonymous splice site variants in atypical presentations.

Integrating Genetic Structural Variations and Whole-Genome Sequencing Into Clinical Neurology

Advances in genome sequencing technologies have unlocked new possibilities in identifying disease-associated and causative genetic markers, which may in turn enhance disease diagnosis and improve prognostication and management strategies. With the capability of examining genetic variations ranging from single-nucleotide mutations to large structural variants, whole-genome sequencing (WGS) is an increasingly adopted approach to dissect the complex genetic architecture of neurologic diseases. There is emerging evidence for different structural variants and their roles in major neurologic and neurodevelopmental diseases. This review first describes different structural variants and their implicated roles in major neurologic and neurodevelopmental diseases, and then discusses the clinical relevance of WGS applications in neurology. Notably, WGS-based detection of structural variants has shown promising potential in enhancing diagnostic power of genetic tests in clinical settings. Ongoing WGS-based research in structural variations and quantifying mutational constraints can also yield clinical benefits by improving variant interpretation and disease diagnosis, while supporting biomarker discovery and therapeutic development. As a result, wider integration of WGS technologies into health care will likely increase diagnostic yields in difficult-to-diagnose conditions and define potential therapeutic targets or intervention points for genome-editing strategies.

Case Report: Exome and RNA Sequencing Identify a Novel de novo Missense Variant in HNRNPK in a Chinese Patient With Au-Kline Syndrome

Au-Kline syndrome is a severe multisystemic syndrome characterized by several congenital defects, including intellectual disability. Loss-of-function and missense variants in the HNRNPK gene are associated with a range of dysmorphic features. This report describes an eleven-year-old Chinese boy with intellectual disability and developmental delays. Family-based whole-exome and Sanger sequencing identified a de novo missense variant in HNRNPK (NM_002140.3: c.143T > A, p. Leu48Val). In silico analysis predicted that this variant would be damaged in a highly conserved residue in the K homology 1 (KH1) domain. Bioinformatic analysis showed that the affinity change (ΔΔG) caused by this variant was -0.033 kcal/mol, indicating that it would have reduced affinity for RNA binding. Transcript analysis of the peripheral blood from this case found 42 aberrantly expressed and 86 aberrantly spliced genes (p-value <0.01). Functional enrichment analysis confirmed that the biological functions of these genes, including protein binding and transcriptional regulation, are associated with HNRNPK. In summary, this study identifies the first Chinese patient with a novel de novo heterozygous HNRNPK gene variant that contributes to Au-Kline syndrome and expands current knowledge of the clinical spectrum of HNRNPK variants.

Clinical implementation of RNA sequencing for Mendelian disease diagnostics

BACKGROUND

Lack of functional evidence hampers variant interpretation, leaving a large proportion of individuals with a suspected Mendelian disorder without genetic diagnosis after whole genome or whole exome sequencing (WES). Research studies advocate to further sequence transcriptomes to directly and systematically probe gene expression defects. However, collection of additional biopsies and establishment of lab workflows, analytical pipelines, and defined concepts in clinical interpretation of aberrant gene expression are still needed for adopting RNA sequencing (RNA-seq) in routine diagnostics.

METHODS

We implemented an automated RNA-seq protocol and a computational workflow with which we analyzed skin fibroblasts of 303 individuals with a suspected mitochondrial disease that previously underwent WES. We also assessed through simulations how aberrant expression and mono-allelic expression tests depend on RNA-seq coverage.

RESULTS

We detected on average 12,500 genes per sample including around 60% of all disease genes-a coverage substantially higher than with whole blood, supporting the use of skin biopsies. We prioritized genes demonstrating aberrant expression, aberrant splicing, or mono-allelic expression. The pipeline required less than 1 week from sample preparation to result reporting and provided a median of eight disease-associated genes per patient for inspection. A genetic diagnosis was established for 16% of the 205 WES-inconclusive cases. Detection of aberrant expression was a major contributor to diagnosis including instances of 50% reduction, which, together with mono-allelic expression, allowed for the diagnosis of dominant disorders caused by haploinsufficiency. Moreover, calling aberrant splicing and variants from RNA-seq data enabled detecting and validating splice-disrupting variants, of which the majority fell outside WES-covered regions.

CONCLUSION

Together, these results show that streamlined experimental and computational processes can accelerate the implementation of RNA-seq in routine diagnostics.

Genome sequencing reveals novel noncoding variants in PLA2G6 and LMNB1 causing progressive neurologic disease

Neurodegenerative disorders and leukodystrophies are progressive neurologic conditions that can occur following the disruption of intricately coordinated patterns of gene expression. Exome sequencing has been adopted as an effective diagnostic tool for determining the underlying genetic etiology of Mendelian neurologic disorders, however genome sequencing offer advantages in its ability to identify and characterize copy number, structural, and sequence variants in noncoding regions. Genome sequencing from peripheral leukocytes was performed on two patients with progressive neurologic disease of unknown etiology following negative genetic investigations including exome sequencing. RNA sequencing from peripheral blood was performed to determine gene expression patterns in one of the patients. Potential causative variants were matched to the patients' clinical presentation. The first proband was found to be heterozygous for a likely pathogenic missense variant in PLA2G6 (c.386T>C; p.Leu129Pro) and have an additional deep intronic variant in PLA2G6 (c.2035-926G>A). RNA sequencing indicated this latter variant created a splice acceptor site leading to the incorporation of a pseudo-exon introducing a premature termination codon. The second proband was heterozygous for a 261 kb deletion upstream of LMNB1 that included an enhancer region. Previous reports of copy number variants spanning this region of cis-acting regulatory elements corroborated its pathogenicity. When combined with clinical presentations, these findings led to a definitive diagnosis of autosomal recessive infantile neuroaxonal dystrophy and autosomal dominant adult-onset demyelinating leukodystrophy, respectively. In patients with progressive neurologic disease of unknown etiology, genome sequencing with the addition of RNA analysis where appropriate should be considered for the identification of causative noncoding pathogenic variants.

RNA sequencing role and application in clinical diagnostic

Although whole-exome sequencing and whole-genome sequencing has tremendously improved our understanding of the genetic etiology of human disorders, about half of the patients still do not receive a molecular diagnosis. The high fraction of variants with uncertain significance and the challenges of interpretation of noncoding variants have urged scientists to implement RNA sequencing (RNA-seq) in the diagnostic approach as a high throughput assay to complement genomic data with functional evidence. RNA-seq data can be used to identify aberrantly spliced genes, detect allele-specific expression, and identify gene expression outliers. Amongst eight studies utilizing RNA-seq, a mean diagnostic uplift of 15% has been reported. Here, we provide an overview of how RNA-seq has been implemented to aid in identifying the causal variants of Mendelian disorders.

MRSD: A quantitative approach for assessing suitability of RNA-seq in the investigation of mis-splicing in Mendelian disease

Variable levels of gene expression between tissues complicates the use of RNA sequencing of patient biosamples to delineate the impact of genomic variants. Here, we describe a gene- and tissue-specific metric to inform the feasibility of RNA sequencing. This overcomes limitations of using expression values alone as a metric to predict RNA-sequencing utility. We have derived a metric, minimum required sequencing depth (MRSD), that estimates the depth of sequencing required from RNA sequencing to achieve user-specified sequencing coverage of a gene, transcript, or group of genes. We applied MRSD across four human biosamples: whole blood, lymphoblastoid cell lines (LCLs), skeletal muscle, and cultured fibroblasts. MRSD has high precision (90.1%-98.2%) and overcomes transcript region-specific sequencing biases. Applying MRSD scoring to established disease gene panels shows that fibroblasts, of these four biosamples, are the optimum source of RNA for 63.1% of gene panels. Using this approach, up to 67.8% of the variants of uncertain significance in ClinVar that are predicted to impact splicing could be assayed by RNA sequencing in at least one of the biosamples. We demonstrate the utility and benefits of MRSD as a metric to inform functional assessment of splicing aberrations, in particular in the context of Mendelian genetic disorders to improve diagnostic yield.

Toward transcriptomics as a primary tool for rare disease investigation

In the past 5 years transcriptome or RNA-sequencing (RNA-seq) has steadily emerged as a complementary assay for rare disease diagnosis and discovery. In this perspective, we summarize several recent developments and challenges in the use of RNA-seq for rare disease investigation. Using an accessible patient sample, such as blood, skin, or muscle, RNA-seq enables the assay of expressed RNA transcripts. Analysis of RNA-seq allows the identification of aberrant or outlier gene expression and alternative splicing as functional evidence to support rare disease study and diagnosis. Further, many types of variant effects can be profiled beyond coding variants, as the consequences of noncoding variants that impact gene expression and splicing can be directly observed. This is particularly apparent for structural variants that disproportionately underlie outlier gene expression and for splicing variants in which RNA-seq can both measure aberrant canonical splicing and detect deep intronic effects. However, a major potential limitation of RNA-seq in rare disease investigation is the developmental and cell type specificity of gene expression as a pathogenic variant's effect may be limited to a specific spatiotemporal context and access to a patient's tissue sample from the relevant tissue and timing of disease expression may not be possible. We speculate that as advances in computational methods and emerging experimental techniques overcome both developmental and cell type specificity, there will be broadening use of RNA sequencing and multiomics in rare disease diagnosis and delivery of precision health.

What Has the Undiagnosed Diseases Network Taught Us About the Clinical Applications of Genomic Testing?

Genetic testing has undergone a revolution in the last decade, particularly with the advent of next-generation sequencing and its associated reductions in costs and increases in efficiencies. The Undiagnosed Diseases Network (UDN) has been a leader in the application of such genomic testing for rare disease diagnosis. This review discusses the current state of genomic testing performed within the UDN, with a focus on the strengths and limitations of whole-exome and whole-genome sequencing in clinical diagnostics and the importance of ongoing data reanalysis. The role of emerging technologies such as RNA and long-read sequencing to further improve diagnostic rates in the UDN is also described. This review concludes with a discussion of the challenges faced in insurance coverage of comprehensive genomic testing as well as the opportunities for a larger role of testing in clinical medicine.

Analysis of Pathogenic Pseudoexons Reveals Novel Mechanisms Driving Cryptic Splicing

Understanding pre-mRNA splicing is crucial to accurately diagnosing and treating genetic diseases. However, mutations that alter splicing can exert highly diverse effects. Of all the known types of splicing mutations, perhaps the rarest and most difficult to predict are those that activate pseudoexons, sometimes also called cryptic exons. Unlike other splicing mutations that either destroy or redirect existing splice events, pseudoexon mutations appear to create entirely new exons within introns. Since exon definition in vertebrates requires coordinated arrangements of numerous RNA motifs, one might expect that pseudoexons would only arise when rearrangements of intronic DNA create novel exons by chance. Surprisingly, although such mutations do occur, a far more common cause of pseudoexons is deep-intronic single nucleotide variants, raising the question of why these latent exon-like tracts near the mutation sites have not already been purged from the genome by the evolutionary advantage of more efficient splicing. Possible answers may lie in deep intronic splicing processes such as recursive splicing or poison exon splicing. Because these processes utilize intronic motifs that benignly engage with the spliceosome, the regions involved may be more susceptible to exonization than other intronic regions would be. We speculated that a comprehensive study of reported pseudoexons might detect alignments with known deep intronic splice sites and could also permit the characterisation of novel pseudoexon categories. In this report, we present and analyse a catalogue of over 400 published pseudoexon splice events. In addition to confirming prior observations of the most common pseudoexon mutation types, the size of this catalogue also enabled us to suggest new categories for some of the rarer types of pseudoexon mutation. By comparing our catalogue against published datasets of non-canonical splice events, we also found that 15.7% of pseudoexons exhibit some splicing activity at one or both of their splice sites in non-mutant cells. Importantly, this included seven examples of experimentally confirmed recursive splice sites, confirming for the first time a long-suspected link between these two splicing phenomena. These findings have the potential to improve the fidelity of genetic diagnostics and reveal new targets for splice-modulating therapies.

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.

RNA sequencing combining with whole exome sequencing reveals a compound heterozygous variant in ATM in a girl with atypical ataxia-telangiectasia

Ataxia-telangiectasia (A-T) is an infrequent autosomal recessive multisystem disorder characterized by progressive cerebellar ataxia, oculo-cutaneous telangiectasia, a tendency to malignancies and variable immunodeficiency. Here we described a 5-year-old girl with atypical A-T symptoms. And 2 different ATM variants c.5939_5948del in exon 40 and c.2639-384A > G in intron 17 were detected by whole exome sequencing (WES) combined with RNA sequencing (RNA-seq). The variant spectrum of ATM was expanded. RNA-seq makes up for deficiencies of WES. We proposed a new approach, a dual-omics that combines RNA-seq with WES, for the diagnosis of genetic diseases. Moreover, our study discussed the phenotypic heterogeneity of A-T among family members as well as individuals. For children with recurrent infections and immunodeficiency, we suggested focusing on A-T after the exclusion of other potential diseases.

Targeted long-read sequencing identifies missing disease-causing variation

Despite widespread clinical genetic testing, many individuals with suspected genetic conditions lack a precise diagnosis, limiting their opportunity to take advantage of state-of-the-art treatments. In some cases, testing reveals difficult-to-evaluate structural differences, candidate variants that do not fully explain the phenotype, single pathogenic variants in recessive disorders, or no variants in genes of interest. Thus, there is a need for better tools to identify a precise genetic diagnosis in individuals when conventional testing approaches have been exhausted. We performed targeted long-read sequencing (T-LRS) using adaptive sampling on the Oxford Nanopore platform on 40 individuals, 10 of whom lacked a complete molecular diagnosis. We computationally targeted up to 151 Mbp of sequence per individual and searched for pathogenic substitutions, structural variants, and methylation differences using a single data source. We detected all genomic aberrations-including single-nucleotide variants, copy number changes, repeat expansions, and methylation differences-identified by prior clinical testing. In 8/8 individuals with complex structural rearrangements, T-LRS enabled more precise resolution of the mutation, leading to changes in clinical management in one case. In ten individuals with suspected Mendelian conditions lacking a precise genetic diagnosis, T-LRS identified pathogenic or likely pathogenic variants in six and variants of uncertain significance in two others. T-LRS accurately identifies pathogenic structural variants, resolves complex rearrangements, and identifies Mendelian variants not detected by other technologies. T-LRS represents an efficient and cost-effective strategy to evaluate high-priority genes and regions or complex clinical testing results.

Genetic Testing and Counseling and Child Neurology

Advances in genetic technology have decreased the cost and increased the accessibility of genetic testing, and introduced new therapeutic options for many genetic conditions. With new treatments available for previously untreatable neurogenetic conditions, identifying a genetic diagnosis has become of great importance. This article provides a review of basic genetic concepts, ethical and counseling considerations with genetic testing, and genetic testing strategies, and highlights a series of clinical care pearls.

Demystifying emerging bulk RNA-Seq applications: the application and utility of bioinformatic methodology

Significant innovations in next-generation sequencing techniques and bioinformatics tools have impacted our appreciation and understanding of RNA. Practical RNA sequencing (RNA-Seq) applications have evolved in conjunction with sequence technology and bioinformatic tools advances. In most projects, bulk RNA-Seq data is used to measure gene expression patterns, isoform expression, alternative splicing and single-nucleotide polymorphisms. However, RNA-Seq holds far more hidden biological information including details of copy number alteration, microbial contamination, transposable elements, cell type (deconvolution) and the presence of neoantigens. Recent novel and advanced bioinformatic algorithms developed the capacity to retrieve this information from bulk RNA-Seq data, thus broadening its scope. The focus of this review is to comprehend the emerging bulk RNA-Seq-based analyses, emphasizing less familiar and underused applications. In doing so, we highlight the power of bulk RNA-Seq in providing biological insights.

Splicing in the Diagnosis of Rare Disease: Advances and Challenges

Mutations which affect splicing are significant contributors to rare disease, but are frequently overlooked by diagnostic sequencing pipelines. Greater ascertainment of pathogenic splicing variants will increase diagnostic yields, ending the diagnostic odyssey for patients and families affected by rare disorders, and improving treatment and care strategies. Advances in sequencing technologies, predictive modeling, and understanding of the mechanisms of splicing in recent years pave the way for improved detection and interpretation of splice affecting variants, yet several limitations still prohibit their routine ascertainment in diagnostic testing. This review explores some of these advances in the context of clinical application and discusses challenges to be overcome before these variants are comprehensively and routinely recognized in diagnostics.

BRCA1 and BRCA2 whole cDNA analysis in unsolved hereditary breast/ovarian cancer patients

Germline pathogenic variants in BRCA1 and BRCA2 genes (BRCA1/2) explain an important fraction of hereditary breast/ovarian cancer (HBOC) cases. Genetic testing generally involves examining coding regions and exon/intron boundaries, thus the frequency of deleterious variants in non-coding regions is unknown. Here we analysed BRCA1/2 whole cDNA in a large cohort of 320 unsolved high-risk HBOC cases in order to identify potential splicing alterations explained by variants in BRCA1/2 deep intronic regions. Whole RNA splicing profiles were analysed by RT-PCR using Sanger sequencing or high-resolution electrophoresis in a QIAxcel instrument. Known predominant BRCA1/2 alternative splicing events were detected, together with two novel events BRCA1 ▼21 and BRCA2 Δ18q_27p. BRCA2 exon 3 skipping was detected in one patient (male) affected with breast cancer, caused by a known Portuguese founder mutation (c.156_157insAluYa5). An altered BRCA2 splicing pattern was detected in three patients, consisting in the up-regulation of ▼20A, Δ22 and ▼20A+Δ22 transcripts. In silico analysis and semi-quantitative data identified the polymorphism BRCA2 c.8755-66T>C as a potential modifier of Δ22 levels. Our findings suggest that mRNA alterations in BRCA1/2 caused by deep intronic variants are rare in Spanish population. However, RNA analysis complements DNA-based strategies allowing the identification of alterations that could go undetected by conventional testing.

Biallelic Pathogenic Variants in TNNT3 Associated With Congenital Myopathy

OBJECTIVE

Pathogenic variants in TNNT3, the gene encoding fast skeletal muscle troponin T, were first described in autosomal dominant distal arthrogryposis type 2B2. Recently, a homozygous splice site variant, c.681+1G>A, was identified in a patient with nemaline myopathy and distal arthrogryposis. Here, we describe the second individual with congenital myopathy associated with biallelic TNNT3 variants.

METHODS

Clinical exome sequencing data from a patient with molecularly undiagnosed congenital myopathy underwent research reanalysis. Clinical and histopathologic data were collected and compared with the single reported patient with TNNT3-related congenital myopathy.

RESULTS

A homozygous TNNT3 variant, c.481-1G>A, was identified. This variant alters a consensus splice acceptor and is predicted to affect splicing by multiple in silico prediction tools. Both the patient reported here and the previously published patient exhibited limb, bulbar, and respiratory muscle weakness from birth, which improved over time. Other shared features include history of polyhydramnios, hypotonia, scoliosis, and high-arched palate. Distal arthrogryposis and nemaline rods, findings reported in the first patient with TNNT3-related congenital myopathy, were not observed in the patient reported here.

CONCLUSIONS

This report provides further evidence for the association of biallelic TNNT3 variants with severe recessive congenital myopathy with or without nemaline rods and distal arthrogryposis. TNNT3 sequencing and copy number analysis should be incorporated into the workup of congenital myopathies.

How Machine Learning and Statistical Models Advance Molecular Diagnostics of Rare Disorders Via Analysis of RNA Sequencing Data

Rare diseases, although individually rare, collectively affect approximately 350 million people worldwide. Currently, nearly 6,000 distinct rare disorders with a known molecular basis have been described, yet establishing a specific diagnosis based on the clinical phenotype is challenging. Increasing integration of whole exome sequencing into routine diagnostics of rare diseases is improving diagnostic rates. Nevertheless, about half of the patients do not receive a genetic diagnosis due to the challenges of variant detection and interpretation. During the last years, RNA sequencing is increasingly used as a complementary diagnostic tool providing functional data. Initially, arbitrary thresholds have been applied to call aberrant expression, aberrant splicing, and mono-allelic expression. With the application of RNA sequencing to search for the molecular diagnosis, the implementation of robust statistical models on normalized read counts allowed for the detection of significant outliers corrected for multiple testing. More recently, machine learning methods have been developed to improve the normalization of RNA sequencing read count data by taking confounders into account. Together the methods have increased the power and sensitivity of detection and interpretation of pathogenic variants, leading to diagnostic rates of 10-35% in rare diseases. In this review, we provide an overview of the methods used for RNA sequencing and illustrate how these can improve the diagnostic yield of rare diseases.

Genetic and phenotypic analysis of 101 patients with developmental delay or intellectual disability using whole-exome sequencing

Whole-exome sequencing (WES) enables identification of pathogenic variants, including copy number variants (CNVs). In this study, we performed WES in 101 Japanese patients with unexplained developmental delay (DD) or intellectual disability (ID) (63 males and 38 females), 98 of them with trio-WES. Pathogenic variants were identified in 54 cases (53.5%), including four cases with pathogenic CNVs. In one case, a pathogenic variant was identified by reanalysis of exome data; and in two cases, two molecular diagnoses were identified. Among 58 pathogenic variants, 49 variants occurred de novo in 48 patients, including two somatic variants. The accompanying autism spectrum disorder and external ear anomalies were associated with detection of pathogenic variants with odds ratios of 11.88 (95% confidence interval [CI] 2.52-56.00) and 3.46 (95% CI 1.23-9.73), respectively. These findings revealed the importance of reanalysis of WES data and detection of CNVs and somatic variants in increasing the diagnostic yield for unexplained DD/ID. In addition, genetic testing is recommended when patients suffer from the autism spectrum disorder or external ear anomalies, which potentially suggests the involvement of genetic factors associated with gene expression regulation.

Detection of aberrant gene expression events in RNA sequencing data

RNA sequencing (RNA-seq) has emerged as a powerful approach to discover disease-causing gene regulatory defects in individuals affected by genetically undiagnosed rare disorders. Pioneering studies have shown that RNA-seq could increase the diagnosis rates over DNA sequencing alone by 8-36%, depending on the disease entity and tissue probed. To accelerate adoption of RNA-seq by human genetics centers, detailed analysis protocols are now needed. We present a step-by-step protocol that details how to robustly detect aberrant expression levels, aberrant splicing and mono-allelic expression in RNA-seq data using dedicated statistical methods. We describe how to generate and assess quality control plots and interpret the analysis results. The protocol is based on the detection of RNA outliers pipeline (DROP), a modular computational workflow that integrates all the analysis steps, can leverage parallel computing infrastructures and generates browsable web page reports.