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

Impact of genome build on RNA-seq interpretation and diagnostics

Transcriptomics is a powerful tool for unraveling the molecular effects of genetic variants and disease diagnosis. Prior studies have demonstrated that choice of genome build impacts variant interpretation and diagnostic yield for genomic analyses. To identify the extent genome build also impacts transcriptomics analyses, we studied the effect of the hg19, hg38, and CHM13 genome builds on expression quantification and outlier detection in 386 rare disease and familial control samples from both the Undiagnosed Diseases Network and Genomics Research to Elucidate the Genetics of Rare Disease Consortium. Across six routinely collected biospecimens, 61% of quantified genes were not influenced by genome build. However, we identified 1,492 genes with build-dependent quantification, 3,377 genes with build-exclusive expression, and 9,077 genes with annotation-specific expression across six routinely collected biospecimens, including 566 clinically relevant and 512 known OMIM genes. Further, we demonstrate that between builds for a given gene, a larger difference in quantification is well correlated with a larger change in expression outlier calling. Combined, we provide a database of genes impacted by build choice and recommend that transcriptomics-guided analyses and diagnoses are cross referenced with these data for robustness.

Impact of genome build on RNA-seq interpretation and diagnostics

Transcriptomics is a powerful tool for unraveling the molecular effects of genetic variants and disease diagnosis. Prior studies have demonstrated that choice of genome build impacts variant interpretation and diagnostic yield for genomic analyses. To identify the extent genome build also impacts transcriptomics analyses, we studied the effect of the hg19, hg38, and CHM13 genome builds on expression quantification and outlier detection in 386 rare disease and familial control samples from both the Undiagnosed Diseases Network (UDN) and Genomics Research to Elucidate the Genetics of Rare Disease (GREGoR) Consortium. We identified 2,800 genes with build-dependent quantification across six routinely-collected biospecimens, including 1,391 protein-coding genes and 341 known rare disease genes. We further observed multiple genes that only have detectable expression in a subset of genome builds. Finally, we characterized how genome build impacts the detection of outlier transcriptomic events. Combined, we provide a database of genes impacted by build choice, and recommend that transcriptomics-guided analyses and diagnoses are cross-referenced with these data for robustness.

RNA Sequencing of Urine-Derived Cells for the Characterization and Diagnosis of Osteogenesis Imperfecta

DNA sequencing is a reliable tool for identifying genetic variants in osteogenesis imperfecta (OI) but cannot always establish pathogenicity, particularly in variants altering splicing. RNA sequencing can provide functional evidence of the effect of a variant on the transcript but requires cells expressing the relevant genes. Here, we used urine-derived cells (UDC) to characterize genetic variants in patients with suspected or confirmed OI and provide evidence on the pathogenicity of variants of uncertain significance (VUS). Urine samples were obtained from 45 children and adolescents; UDC culture was successful in 40 of these participants (age range 4-20 years, 21 females), including 18 participants with OI or suspected OI who had a candidate variant or VUS on DNA sequencing. RNA was extracted from UDC and sequenced on an Illumina NextSeq550 device. Principal component analysis showed that the gene expression profiles of UDC and fibroblasts (based on Genotype Tissue Expression [GTEx] Consortium data) clustered close together and had less variability than those of whole blood cells. Transcript abundance was sufficient for analysis by RNA sequencing (defined as a median gene expression level of ≥10 transcripts per million) for 25 of the 32 bone fragility genes (78%) that were included in our diagnostic DNA sequencing panel. These results were similar to GTEx data for fibroblasts. Abnormal splicing was identified in 7 of the 8 participants with pathogenic or likely pathogenic variants in the splice region or deeper within the intron. Abnormal splicing was also observed in 2 VUS (COL1A1 c.2829+5G>A and COL1A2 c.693+6T>G), but no splice abnormality was observed in 3 other VUS. Abnormal deletions and duplications could also be observed in UDC transcripts. In conclusion, UDC are suitable for RNA transcript analysis in patients with suspected OI and can provide functional evidence for pathogenicity, in particular of variants affecting splicing. © 2023 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Building the foundation for a community-generated national research blueprint for inherited bleeding disorders: research priorities for ultra-rare inherited bleeding disorders

BACKGROUND

Ultra-rare inherited bleeding disorders (BDs) present important challenges for generating a strong evidence foundation for optimal diagnosis and management. Without disorder-appropriate treatment, affected individuals potentially face life-threatening bleeding, delayed diagnosis, suboptimal management of invasive procedures, psychosocial distress, pain, and decreased quality-of-life.

RESEARCH DESIGN AND METHODS

The National Hemophilia Foundation (NHF) and the American Thrombosis and Hemostasis Network identified the priorities of people with inherited BDs and their caregivers, through extensive inclusive community consultations, to inform a blueprint for future decades of research. Multidisciplinary expert Working Group (WG) 3 distilled highly feasible transformative ultra-rare inherited BD research opportunities from the community-identified priorities.

RESULTS

WG3 identified three focus areas with the potential to advance the needs of all people with ultra-rare inherited BDs and scored the feasibility, impact, and risk of priority initiatives, including 13 in systems biology and mechanistic science; 2 in clinical research, data collection, and research infrastructure; and 5 in the regulatory process for novel therapeutics and required data collection.

CONCLUSIONS

Centralization and expansion of expertise and resources, flexible innovative research and regulatory approaches, and inclusion of all people with ultra-rare inherited BDs and their health care professionals will be essential to capitalize on the opportunities outlined herein.

Novel compound heterozygous variants in EMC1 associated with global developmental delay: a lesson from a non-silent synonymous exonic mutation

Background

The endoplasmic reticulum-membrane protein complex (EMC) as a molecular chaperone is required for the proper synthesis, folding and traffic of several transmembrane proteins. Variants in the subunit 1 of EMC (EMC1) have been implicated in neurodevelopmental disorders.

Methods

Whole exome sequencing (WES) with Sanger sequencing validation was performed for a Chinese family, including the proband (a 4-year-old girl who displayed global developmental delay, severe hypotonia and visual impairment), her affected younger sister and her non-consanguineous parents. RT-PCR assay and Sanger sequencing were used to detect abnormal RNA splicing.

Results

Novel compound heterozygous variants in EMC1, including the maternally inherited chr1: 19566812_1956800delinsATTCTACTT[hg19];NM_015047.3:c.765_777delins ATTCTACTT;p.(Leu256fsTer10) and the paternally inherited chr1:19549890G> A[hg19];NM_015047.3:c.2376G>A;p.(Val792=) are identified in the proband and her affected sister. RT-PCR assay followed by Sanger sequencing reveals that the c.2376G>A variant leads to aberrant splicing, with retention of intron 19 (561bp) in the mature mRNA, which is presumed to introduce a premature translational termination codon (p.(Val792fsTer31)).

Conclusion

Novel compound heterozygous variants in EMC1 have been identified in individuals with global developmental delay. Non-silent synonymous mutations should be kept in mind in genetic analysis.

De Novo Somatic Mosaicism of CYBB Caused by Intronic LINE-1 Element Insertion Resulting in Chronic Granulomatous Disease

Chronic granulomatosis disease (CGD) is a rare inborn error of immunity, characterized by phagocytic respiratory outbreak dysfunction. Mutations causing CGD occur in CYBB on the X chromosome and in the autosomal genes CYBA, NCF1, NCF2, NCF4, RAC2, and CYBC1. Nevertheless, some patients are clinically diagnosed with CGD, due to abnormal respiratory outbursts, while the pathogenic gene mutation is unidentified. Here, we report a patient with CGD who first presented with Bacillus Calmette-Guérin disease and had recurrent pneumonia. He was diagnosed with CGD by nitro blue tetrazolium and respiratory burst tests. Detailed assessment of neutrophil activity revealed that patient neutrophils were almost entirely nonfunctional. Sanger sequencing detected a 6-kb insertion of a LINE-1 transposable element in the third intron of CYBB, leading to abnormal splicing and pseudoexon insertion, as well as introduction of a premature termination codon, resulting in predicted protein truncation. Clonal analysis demonstrated that the patient had somatic mosaicism, and the phagocytes were almost all variant CYBB, while the mosaicism rate of PBMC was about 65%. Finally, deep RNA sequencing and gp91^phox expression analysis confirmed the pathogenicity of the mutation. In conclusion, we demonstrate that insertion of a LINE-1 transposon in a CYBB intron was responsible for CGD in our patient. Intron LINE-1 transposon element insertion should be examined in CGD patients without any known disease-causing gene mutation, in addition to identification of new genes.

The Utility of Small Fishes for the Genetic Study of Human Age-Related Disorders

Animal models have been used to model human diseases, and among them, small fishes have been highlighted for their usefulness in various ways, such as the low cost of maintenance, ease of genetic modification, small size for easy handling, and strength in imaging studies due to their relative transparency. Recently, the use of turquoise killifish, Nothobranchius furzeri, which is known to exhibit various aging phenotypes in a short period, has attracted attention in research on aging and age-related diseases. However, when using animal models, it is important to keep their genetic background and interspecies differences in mind for translating them into human diseases. In this article, we obtained the gene symbols of protein-coding genes of turquoise killifish, medaka, zebrafish, and humans from NCBI datasets and extracted common shared genes among four species to explore the potential of interspecies translational research and to apply small fish models for human age-related disorders. Common shared protein-coding genes were analyzed with the Reactome Pathway Database to determine the coverage of these genes in each pathway in humans. We applied common shared genes to the Orphanet database to establish a list of human diseases that contain common shared genes among the four species. As examples, the senescence-related pathways and some pathways of human age-related diseases, such as Alzheimer’s disease, Parkinson’s disease, frontotemporal dementia, nonalcoholic fatty liver disease, progeria, hepatocellular carcinoma, and renal cell carcinoma, were extracted from the curated pathway and disease list to discuss the further utility of fish models for human age-related disorders.