An integrative framework for clinical diagnosis and knowledge discovery from exome sequencing data

Non-silent single nucleotide genetic variants, like nonsense changes and insertion-deletion variants, that affect protein function and length substantially are prevalent and are frequently misclassified. The low sensitivity and specificity of existing variant effect predictors for nonsense and indel variations restrict their use in clinical applications. We propose the Pathogenic Mutation Prediction (PMPred) method to predict the pathogenicity of single nucleotide variations, which impair protein function by prematurely terminating a protein's elongation during its synthesis. The prediction starts by monitoring functional effects (Gene Ontology annotation changes) of the change in sequence, using an existing ensemble machine learning model (UniGOPred). This, in turn, reveals the mutations that significantly deviate functionally from the wild-type sequence. We have identified novel harmful mutations in patient data and present them as motivating case studies. We also show that our method has increased sensitivity and specificity compared to state-of-the-art, especially in single nucleotide variations that produce large functional changes in the final protein. As further validation, we have done a comparative docking study on such a variation that is misclassified by existing methods and, using the altered binding affinities, show how PMPred can correctly predict the pathogenicity when other tools miss it. PMPred is freely accessible as a web service at https://pmpred.kansil.org/, and the related code is available at https://github.com/kansil/PMPred.

DNA methylome, R-loop and clinical exome profiling of patients with sporadic amyotrophic lateral sclerosis

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disorder characterized by the death of motor neurons, the aetiology of which is essentially unknown. Here, we present an integrative epigenomic study in blood samples from seven clinically characterised sporadic ALS patients to elucidate molecular factors associated with the disease. We used clinical exome sequencing (CES) to study DNA variants, DNA-RNA hybrid immunoprecipitation sequencing (DRIP-seq) to assess R-loop distribution, and reduced representation bisulfite sequencing (RRBS) to examine DNA methylation changes. The above datasets were combined to create a comprehensive repository of genetic and epigenetic changes associated with the ALS cases studied. This repository is well-suited to unveil new correlations within individual patients and across the entire patient cohort. The molecular attributes described here are expected to guide further mechanistic studies on ALS, shedding light on the underlying genetic causes and facilitating the development of new epigenetic therapies to combat this life-threatening disease.

Human whole-exome genotype data for Alzheimer's disease

The heterogeneity of the whole-exome sequencing (WES) data generation methods present a challenge to a joint analysis. Here we present a bioinformatics strategy for joint-calling 20,504 WES samples collected across nine studies and sequenced using ten capture kits in fourteen sequencing centers in the Alzheimer's Disease Sequencing Project. The joint-genotype called variant-called format (VCF) file contains only positions within the union of capture kits. The VCF was then processed specifically to account for the batch effects arising from the use of different capture kits from different studies. We identified 8.2 million autosomal variants. 96.82% of the variants are high-quality, and are located in 28,579 Ensembl transcripts. 41% of the variants are intronic and 1.8% of the variants are with CADD > 30, indicating they are of high predicted pathogenicity. Here we show our new strategy can generate high-quality data from processing these diversely generated WES samples. The improved ability to combine data sequenced in different batches benefits the whole genomics research community.

TG468: a text graph convolutional network for predicting clinical response to immune checkpoint inhibitor therapy

Enhancing cancer treatment efficacy remains a significant challenge in human health. Immunotherapy has witnessed considerable success in recent years as a treatment for tumors. However, due to the heterogeneity of diseases, only a fraction of patients exhibit a positive response to immune checkpoint inhibitor (ICI) therapy. Various single-gene-based biomarkers and tumor mutational burden (TMB) have been proposed for predicting clinical responses to ICI; however, their predictive ability is limited. We propose the utilization of the Text Graph Convolutional Network (GCN) method to comprehensively assess the impact of multiple genes, aiming to improve the predictive capability for ICI response. We developed TG468, a Text GCN model framing drug response prediction as a text classification task. By combining natural language processing (NLP) and graph neural network techniques, TG468 effectively handles sparse and high-dimensional exome sequencing data. As a result, TG468 can distinguish survival time for patients who received ICI therapy and outperforms single gene biomarkers, TMB and some classical machine learning models. Additionally, TG468's prediction results facilitate the identification of immune status differences among specific patient types in the Cancer Genome Atlas dataset, providing a rationale for the model's predictions. Our approach represents a pioneering use of a GCN model to analyze exome data in patients undergoing ICI therapy and offers inspiration for future research using NLP technology to analyze exome sequencing data.

Diagnostic genome sequencing improves diagnostic yield: a prospective single-centre study in 1000 patients with inherited eye diseases

PURPOSE

Genome sequencing (GS) is expected to reduce the diagnostic gap in rare disease genetics. We aimed to evaluate a scalable framework for genome-based analyses 'beyond the exome' in regular care of patients with inherited retinal degeneration (IRD) or inherited optic neuropathy (ION).

METHODS

PCR-free short-read GS was performed on 1000 consecutive probands with IRD/ION in routine diagnostics. Complementary whole-blood RNA-sequencing (RNA-seq) was done in a subset of 74 patients. An open-source bioinformatics analysis pipeline was optimised for structural variant (SV) calling and combined RNA/DNA variation interpretation.

RESULTS

A definite genetic diagnosis was established in 57.4% of cases. For another 16.7%, variants of uncertain significance were identified in known IRD/ION genes, while the underlying genetic cause remained unresolved in 25.9%. SVs or alterations in non-coding genomic regions made up for 12.7% of the observed variants. The RNA-seq studies supported the classification of two unclear variants.

CONCLUSION

GS is feasible in clinical practice and reliably identifies causal variants in a substantial proportion of individuals. GS extends the diagnostic yield to rare non-coding variants and enables precise determination of SVs. The added diagnostic value of RNA-seq is limited by low expression levels of the major IRD disease genes in blood.

Genomics and tumor microenvironment of breast mucoepidermoid carcinoma based on whole-exome and RNA sequencing

Mammary mucoepidermoid carcinoma (MEC) is a rare entity. The molecular characteristics of breast MEC have not been fully investigated due to its rarity. We performed a retrospective study among 1000 patients with breast carcinomas and identified four cases of breast MEC. Clinical and demographic data were collected. Immunohistochemistry panels which were used to diagnose salivary gland MEC and breast carcinomas were also performed. MAML2 rearrangements were detected by FISH and fusion partners were identified by RNA sequencing. Whole-exome sequencing (WES) was used to reveal the genomes of these four breast MEC. Then, the biological functions and features of breast MEC were further compared with those of invasive breast carcinomas and salivary gland MEC.According to Ellis and Auclair's methods, these four breast MEC could be classified as low-grade breast MEC. All the patients were alive, and disease-free survival (PFS) ranged from 20 months to 67 months. Among these four breast MEC, two cases were triple-negative, and the other two cases were found to be ER positive, with one also showing HER2 equivocal by immunohistochemical staining, but no amplification in FISH. FISH analysis confirmed the presence of the MAML2 translocation in three of four tumors, and CRTC1-MAML2 fusion was confirmed in two of them by RNA-sequencing. The average coverage size of WES for the tumor mutation burden estimation was 32 Mb. MUC4, RP1L1 and QRICH2 mutations were identified in at least three tumors, and these mutation also existed in breast invasive carcinoma databases (TCGA, Cell 2015; TCGA, Nature 2012). The results showed that there were many genes in breast MEC overlapping with the breast invasive carcinoma databases mentioned above, range from 5 to 63 genes (median:21 genes). Next, we assessed immune cell infiltration levels in these tumors. In all these tumors, M2 macrophages and plasma cell were in the high infiltration group. Our breast MEC showed different results from the salivary gland MEC, whose plasma cells were in the low infiltration group. Overall, we first analyzed the genomics and tumor microenvironment of breast mucoepidermoid carcinoma and proposed our hypothesis that although MECs arising in the breast resemble their salivary gland counterparts phenotypically, our findings indicate that breast MECs probably resemble invasive breast carcinomas at the genetic level and immune cell infiltration levels. More cases and in deep research need to be done to further understand this rare carcinoma.

Whole-Exome Screening and Analysis of Signaling Pathways in Multiple Endocrine Neoplasia Type 1 Patients with Different Outcomes: Insights into Cellular Mechanisms and Possible Functional Implications

Multiple endocrine neoplasia type 1 (MEN1) is a syndrome characterized by tumors in multiple organs. Although being a dominantly inherited monogenic disease, disease phenotypes are unpredictable and differ even among members of the same family. There is growing evidence for the role of modifier genes in the alteration of the course of this disease. However, genome-wide screening data are still lacking. In our study, we addressed the different outcomes of the disease, focusing on pituitary and adrenocortical tumors. By means of exome sequencing we identified the affected signaling pathways that segregated with those symptoms. Most significantly, we identified damaging alterations in numerous structural genes responsible for cell adhesion and migration. Additionally, in the case of pituitary tumors, genes related to neuronal function, survival, and morphogenesis were repeatedly identified, while in patients with adrenocortical tumors, TLR10, which is involved in the regulation of the innate immunity, was commonly modified. Our data show that using exome screening, it is possible to find signatures which correlate with the given clinical MEN1 outcomes, providing evidence that studies addressing modifier effects in MEN1 are reasonable.

Comprehensive review and expanding the genetic landscape of Cornelia-de-Lange spectrum: insights from novel mutations and skin biopsy in exome-negative cases

BACKGROUND

Cornelia de Lange Syndrome (CdLS) is a rare genetic disorder characterized by a range of physical, cognitive, and behavioral abnormalities. This study aimed to perform a comprehensive review of the literature on CdLS and investigate two cases of CdLS with distinct phenotypes that underwent WES to aid in their diagnosis.

METHODS

We conducted a comprehensive review of the literature on CdLS along with performing whole-exome sequencing on two CdLS patients with distinct phenotypes, followed by Sanger sequencing validation and in-silico analysis.

RESULTS

The first case exhibited a classic CdLS phenotype, but the initial WES analysis of blood-derived DNA failed to identify any mutations in CdLS-related genes. However, a subsequent WES analysis of skin-derived DNA revealed a novel heterozygous mutation in the NIPBL gene (NM_133433.4:c.6534_6535del, p.Met2178Ilefs*8). The second case was presented with a non-classic CdLS phenotype, and WES analysis of blood-derived DNA identified a heterozygous missense variant in the SMC1A gene (NM_006306.4:c.2320G>A, p.Asp774Asn).

CONCLUSIONS

The study shows the importance of considering mosaicism in classic CdLS cases and the value of WES for identifying genetic defects. These findings contribute to our understanding of CdLS genetics and underscore the need for comprehensive genetic testing to enhance the diagnosis and management of CdLS patients.

Using multi-scale genomics to associate poorly annotated genes with rare diseases

BACKGROUND

Next-generation sequencing (NGS) has significantly transformed the landscape of identifying disease-causing genes associated with genetic disorders. However, a substantial portion of sequenced patients remains undiagnosed. This may be attributed not only to the challenges posed by harder-to-detect variants, such as non-coding and structural variations but also to the existence of variants in genes not previously associated with the patient's clinical phenotype. This study introduces EvORanker, an algorithm that integrates unbiased data from 1,028 eukaryotic genomes to link mutated genes to clinical phenotypes.

METHODS

EvORanker utilizes clinical data, multi-scale phylogenetic profiling, and other omics data to prioritize disease-associated genes. It was evaluated on solved exomes and simulated genomes, compared with existing methods, and applied to 6260 knockout genes with mouse phenotypes lacking human associations. Additionally, EvORanker was made accessible as a user-friendly web tool.

RESULTS

In the analyzed exomic cohort, EvORanker accurately identified the "true" disease gene as the top candidate in 69% of cases and within the top 5 candidates in 95% of cases, consistent with results from the simulated dataset. Notably, EvORanker outperformed existing methods, particularly for poorly annotated genes. In the case of the 6260 knockout genes with mouse phenotypes, EvORanker linked 41% of these genes to observed human disease phenotypes. Furthermore, in two unsolved cases, EvORanker successfully identified DLGAP2 and LPCAT3 as disease candidates for previously uncharacterized genetic syndromes.

CONCLUSIONS

We highlight clade-based phylogenetic profiling as a powerful systematic approach for prioritizing potential disease genes. Our study showcases the efficacy of EvORanker in associating poorly annotated genes to disease phenotypes observed in patients. The EvORanker server is freely available at https://ccanavati.shinyapps.io/EvORanker/ .

HAT: de novo variant calling for highly accurate short-read and long-read sequencing data

MOTIVATION

de novo variants (DNVs) are variants that are present in offspring but not in their parents. DNVs are both important for examining mutation rates as well as in the identification of disease-related variation. While efforts have been made to call DNVs, calling of DNVs is still challenging from parent-child sequenced trio data. We developed Hare And Tortoise (HAT) as an automated DNV detection workflow for highly accurate short-read and long-read sequencing data. Reliable detection of DNVs is important for human genomics and HAT addresses this need.

RESULTS

HAT is a computational workflow that begins with aligned read data (i.e. CRAM or BAM) from a parent-child sequenced trio and outputs DNVs. HAT detects high-quality DNVs from Illumina short-read whole-exome sequencing, Illumina short-read whole-genome sequencing, and highly accurate PacBio HiFi long-read whole-genome sequencing data. The quality of these DNVs is high based on a series of quality metrics including number of DNVs per individual, percent of DNVs at CpG sites, and percent of DNVs phased to the paternal chromosome of origin.

AVAILABILITY AND IMPLEMENTATION

https://github.com/TNTurnerLab/HAT.

Carrier frequency estimation of pathogenic variants of autosomal recessive and X-linked recessive mendelian disorders using exome sequencing data in 1,642 Thais

BACKGROUND

People with autosomal recessive disorders often were born without awareness of the carrier status of their parents. The American College of Medical Genetics and Genomics (ACMG) recommends screening 113 genes known to cause autosomal recessive and X-linked conditions in couples seeking to learn about their risk of having children with these disorders to have an appropriate reproductive plan.

METHODS

We analyzed the exome sequencing data of 1,642 unrelated Thai individuals to identify the pathogenic variant (PV) frequencies in genes recommended by ACMG.

RESULTS

In the 113 ACMG-recommended genes, 165 PV and likely PVs in 60 genes of 559 exomes (34%, 559/1642) were identified. The carrier rate was increased to 39% when glucose-6-phosphate dehydrogenase (G6PD) was added. The carrier rate was still as high as 14.7% when thalassemia and hemoglobinopathies were excluded. In addition to thalassemia, hemoglobinopathies, and G6PD deficiency, carrier frequencies of > 1% were found for Gaucher disease, primary hyperoxaluria, Pendred syndrome, and Wilson disease. Nearly 2% of the couples were at risk of having offsprings with the tested autosomal recessive conditions.

CONCLUSIONS

Based on the study samples, the expanded carrier screening, which specifically targeted common autosomal recessive conditions in Thai individuals, will benefit clinical outcomes, regarding preconception/prenatal genetic carrier screening.

Phenotype-driven reanalysis reveals five novel pathogenic variants in 40 exome-negative families with Charcot-Marie-Tooth Disease

BACKGROUND

To identify genetic causes in 40 whole exome sequencing (WES)-negative Charcot-Marie-Tooth (CMT) families and provide a summary of the clinical and genetic features of the diagnosed patients.

METHODS

The clinical information and sequencing data of 40 WES-negative families out of 131 CMT families were collected, and phenotype-driven reanalysis was conducted using the Exomiser software.

RESULTS

The molecular diagnosis was regained in 4 families, increasing the overall diagnosis rate by 3.0%. One family with adolescent-onset pure CMT1 was diagnosed [POLR3B: c.2810G>A (p.R937Q)] due to the novel genotype-phenotype association. One infantile-onset, severe CMT1 family with deep sensory disturbance was diagnosed by screening the BAM file and harbored c.1174C>T (p.R392*) and 875_927delinsCTGCCCACTCTGCCCACTCTGCCCACTCTG (p.V292Afs53) of PRX. Two families were diagnosed due to characteristic phenotypes, including an infantile-onset ICMT family with renal dysfunction harboring c.213_233delinsGAGGAGCA (p.S72Rfs34) of INF2 and an adolescent-onset CMT2 family with optic atrophy harboring c.560C>T (p.P187L) and c.616A>G (p.K206E) of SLC25A46. According to the American College of Medical Genetics and Genomics (ACMG) guidelines, the variants of POLR3B and SLC25A46 were classified as likely pathogenic, and the variants of INF2 and PRX were pathogenic. All these variants were first reported worldwide except for p.R392* of PRX.

CONCLUSIONS

We identified five novel pathogenic variants in POLR3B, PRX, INF2, and SLC25A46, which broaden their phenotypic and genotypic spectrums. Regular phenotype-driven reanalysis is a powerful strategy for increasing the diagnostic yield of WES-negative CMT patients, and long-term follow-up and screening BAM files for contiguous deletion and missense variants are both essential for reanalysis.

Whole exome germline sequencing in early-onset prostate cancer patients: Genomic findings and clinical outcomes

BACKGROUND

Whole exome sequencing (WES) furthered our understanding of various tumors. We assessed the occurrence of germline likely pathogenic/pathogenic (LP/P) variants, disease features, and clinical outcomes in early-onset prostate cancer.

METHODS

This retrospective study (N = 134) included consecutive prostate cancer patients who donated blood samples for research purposes to the Kahn-Sagol-Maccabi biobank. Patients diagnosed at >65 years were excluded. Clinical characteristics were extracted from the medical records. Germline WES was performed with analysis reporting on oncogenetic, two immunogenic, and a secondary minimum list panels (121, 468, 76, and 59 genes, respectively).

RESULTS

Median age at diagnosis was 61 (range 46-65) years; 131 (98%) were diagnosed with local disease. The median follow-up time from diagnosis was 14 (range <1-25) years. Of the patients with local disease, 32 (24%) and 10 (8%) had biochemical and distant recurrences, respectively. Twenty-five patients (19%) had ≥1 additional cancer (excluding non-melanoma skin cancer), most frequently bladder (6), colorectal (5), and lymphoma (5). Seven (5%) deaths were reported, with only one related to prostate cancer. LP/P variants were identified in 8 patients (6%), all in genes from the oncogenetic panel: ATM, BRCA1 (in two patients), BRCA2 (in two patients), HOXB13, MUTYH, and MYH7. Of these eight patients, with a median follow-up of 7 years (range <1-15), two (25%) had biochemical recurrences, one had (12.5%) distant recurrence, and no deaths were reported.

CONCLUSIONS

In this cohort of 134 early-onset prostate cancer patients, we identified germline LP/P variants in an oncogenetic panel in 6% of participants, with no unique clinical outcome.

Intestinal-Type Adenocarcinoma in Head and Neck: Dissecting Oncogenic Gene Alterations Through Whole Transcriptome and Exome Analysis

Adenocarcinomas of the nasal/paranasal sinuses are uncommon, but intestinal-type adenocarcinomas (ITACs) are important. Due to the rarity of these tumors, their molecular profile is not well known. To further investigate the molecular profile and find potential oncogenic drivers, we compared the whole transcriptome and exome of ITACs at different anatomic locations in the head and neck. Twenty-one head and neck adenocarcinomas were used in this study, divided into 10 sinonasal adenocarcinomas (SNT) and 11 extrasinonasal (T) head and neck adenocarcinomas according to anatomic location and histology. Tumor samples along with normal mucosa were microdissected from formalin-fixed, paraffin-embedded samples, and RNA and DNA were subjected to whole-transcriptome and -exome shotgun sequencing. Analysis of ITACs at sinonasal locations showed 410 subtype-specific differentially expressed (DE) genes and noncoding transcripts compared with the group of other anatomic locations, with 2909 subtype-specific DE genes. The groups shared 872 genes, with 17 highly different or opposing DE genes. Whole-exome mutation analysis revealed the gene MLL3 (KMT2C) to be exhibiting the most frequent loss-of-function mutations in all adenocarcinomas investigated. The results suggest that the head and neck ITACs investigated were mainly caused by loss-of-function mutations in MLL3 that disabled chromatin methylation and remodeling of all MLL3-targeted enhancers in the tumors. This changed the activity of multiple genes/gene clusters, supporting oncogenicity mostly via pathways of signaling, dedifferentiation, proliferation, migration, and immune and inflammatory deregulation, indicating a truly epigenetic event as the root cause for the heterogenous diversity of these enteric types of cancer. The data of this study form the basis for understanding cell fate determination and cellular homeostasis in the normal respiratory mucosa at different anatomic sites and show the contribution of different mucosal components to the etiology/molecular pathology of ITAC.

Whole exome sequencing reveals novel candidate variants for endometriosis utilizing multiple affected members in a single family

BACKGROUND

Endometriosis is an estrogen-dependent, chronic inflammatory disease that affects 10% of women during the reproductive ages. Despite the estimated 50% heritability for the condition, only 26% was associated with common genetic variants. Thus, necessity of identifying rare variants for the missing heritability is implicated in the literature. Therefore, our study aimed to identify novel rare genetic variants involved in the pathogenesis of endometriosis utilizing a family of multiple affected members.

METHODS

A family composed of four affected women along with their two unaffected mothers were recruited at a single gynecology and infertility clinic specialized in endometriosis. All patients presented with endometriomas, which was visualized by transvaginal ultrasonography. Two affected individuals had received laparoscopic endometrioma excision and therefore were diagnosed with recurrent disease. One mother had a history of endometrial serous adenocarcinoma (ESC) for which she underwent hysterectomy with bilateral oophorectomy. Three endometriosis cases were whole exome sequenced on Illumina NextSeq 550 platform with an average of 90% coverage. Candidate genes were confirmed by Sanger sequencing and followed-up with family segregation.

RESULTS

Novel rare variants were identified in TNFRSF1B (NM_001066.3: c.1072G>A, p.(Ala358Thr)) and GEN1 (NM_001130009.3: c.1574C>T, p.(Ser525Leu)) as possible genetic causes of endometriosis. A third novel rare variant was identified in CRABP1 (NM_004378.3:c.54G>C, p.(Glu18Asp)) only on the mother with ESC history and her daughters.

CONCLUSION

Novel candidate genetic variants that might contribute to endometriosis were suggested that need replication through independent cohorts or validation by functional studies. The family has also received genetic counseling and that the affected daughters are on clinical follow-up, accordingly.

Exome-informed formulations of food proteins enhance body growth and feed conversion efficiency in ad libitum-fed mice

Meeting requirements for dietary proteins, especially of essential amino acids (EAAs), is critical for the life-long health of living organisms. However, defining EAA targets for preparing biologically-matched nutrition that satisfies metabolic requirements for protein remains challenging. Previous research has shown the advantages of 'exome matching' in representing the specific requirement of dietary AAs, where the target dietary AA profile was derived from in silico translation of the genome of an organism, specifically responsible for protein expression (the 'exome'). However, past studies have assessed these effects in only one sex, for few parameters (body mass and composition), and have used purified diets in which protein is supplied as a mixture of individual AAs. Here, for the first time, we utilise a computational method to guide the formulation of custom protein blends and test if exome matching can be achieved at the intact protein level, through blending standard protein ingredients, ultimately leading to optimal growth, longevity and reproductive function. Mice were provided ad libitum (ad lib) access to one of the four iso-energetic protein-limited diets, two matched and two mis-matched to the mouse exome target, and fed at a fixed protein energy level of 6.2%. During or following 13-weeks of feeding, the food intake, body growth, composition and reproductive functions were measured. Compared to the two mis-matched diets, male and female animals on the exome-matched diet with protein digestibility correction applied, exhibited significantly improved growth rates and final body mass. The feed conversion efficiency in the same diet was also increased by 62% and 40% over the worst diets for males and females, respectively. Male, not female, exhibited higher accretion of lean body mass with the matched, digestibility-corrected diet. All reproductive function measures in both sexes were comparable among diets, with the exception of testicular daily sperm production in males, which was higher in the two matched diets versus the mis-matched diets. The results collectively demonstrate the pronounced advantages of exome-matching in supporting body growth and improving feed conversion efficiency in both sexes. However, the potential impact of this approach in enhancing fertility needs further investigation.

A comprehensive assessment of exome capture methods for RNA sequencing of formalin-fixed and paraffin-embedded samples

RNA-Seq analysis of Formalin-Fixed and Paraffin-Embedded (FFPE) samples has emerged as a highly effective approach and is increasingly being used in clinical research and drug development. However, the processing and storage of FFPE samples are known to cause extensive degradation of RNAs, which limits the discovery of gene expression or gene fusion-based biomarkers using RNA sequencing, particularly methods reliant on Poly(A) enrichment. Recently, researchers have developed an exome targeted RNA-Seq methodology that utilizes biotinylated oligonucleotide probes to enrich RNA transcripts of interest, which could overcome these limitations. Nevertheless, the standardization of this experimental framework, including probe designs, sample multiplexing, sequencing read length, and bioinformatic pipelines, remains an essential requirement. In this study, we conducted a comprehensive comparison of three main commercially available exome capture kits and evaluated key experimental parameters, to provide the overview of the advantages and limitations associated with the selection of library preparation protocols and sequencing platforms. The results provide valuable insights into the best practices for obtaining high-quality data from FFPE samples.

Beyond the exome: utility of long-read whole genome sequencing in exome-negative autosomal recessive diseases

BACKGROUND

Long-read whole genome sequencing (lrWGS) has the potential to address the technical limitations of exome sequencing in ways not possible by short-read WGS. However, its utility in autosomal recessive Mendelian diseases is largely unknown.

METHODS

In a cohort of 34 families in which the suspected autosomal recessive diseases remained undiagnosed by exome sequencing, lrWGS was performed on the Pacific Bioscience Sequel IIe platform.

RESULTS

Likely causal variants were identified in 13 (38%) of the cohort. These include (1) a homozygous splicing SV in TYMS as a novel candidate gene for lethal neonatal lactic acidosis, (2) a homozygous non-coding SV that we propose impacts STK25 expression and causes a novel neurodevelopmental disorder, (3) a compound heterozygous SV in RP1L1 with complex inheritance pattern in a family with inherited retinal disease, (4) homozygous deep intronic variants in LEMD2 and SNAP91 as novel candidate genes for neurodevelopmental disorders in two families, and (5) a promoter SNV in SLC4A4 causing non-syndromic band keratopathy. Surprisingly, we also encountered causal variants that could have been identified by short-read exome sequencing in 7 families. The latter highlight scenarios that are especially challenging at the interpretation level.

CONCLUSIONS

Our data highlight the continued need to address the interpretation challenges in parallel with efforts to improve the sequencing technology itself. We propose a path forward for the implementation of lrWGS sequencing in the setting of autosomal recessive diseases in a way that maximizes its utility.

A novel 268 kb deletion combined with a splicing variant in IL7R causes of severe combined immunodeficiency in a Chinese family: a case report

BACKGROUND

Severe combined immunodeficiency (SCID) is a group of fatal primary immunodeficiencies characterized by the severe impairment of T-cell differentiation. IL7R deficiency is a rare form of SCID that usually presents in the first months of life with severe and opportunistic infections, failure to thrive, and a high risk of mortality unless treated. Although recent improvements in early diagnosis have been achieved through newborn screening, few IL7R-related SCID patients had been reported in the Chinese population.

CASE PRESENTATION

Here, we retrospectively analyzed a case of SCID in a 5-month-old girl with symptoms, including severe T-cell depletion, recurrent fever, oral ulcers, pneumonia, hepatosplenomegaly, bone marrow hemophagocytosis, and bacterial and viral infections. Whole-exome sequencing (WES), quantitative PCR (qPCR), and chromosome microarray analysis (CMA) were performed to identify the patient's genetic etiology. We identified a 268 kb deletion and a splicing variant, c.221 + 1G > A, in the proband. These two variants of IL7R were inherited from the father and mother.

CONCLUSIONS

To our knowledge, this is the first report of whole IL7R gene deletion in combination with a pathogenic splicing variant in a patient with SCID. This deletion also expands the pathogenic variation spectrum of SCID caused by IL7R. The incorporation of exome-based copy number variant analysis makes WES a powerful molecular diagnostic technique for the clinical diagnosis of pediatric patients.

Benchmarking and improving the performance of variant-calling pipelines with RecallME

MOTIVATION

The steady increment of Whole Genome/Exome sequencing and the development of novel Next Generation Sequencing-based gene panels requires continuous testing and validation of variant calling (VC) pipelines and the detection of sequencing-related issues to be maintained up-to-date and feasible for the clinical settings. State of the art tools are reliable when used to compute standard performance metrics. However, the need for an automated software to discriminate between bioinformatic and sequencing issues and to optimize VC parameters remains unmet.

RESULTS

The aim of the current work is to present RecallME, a bioinformatic suite that tracks down difficult-to-detect variants as insertions and deletions in highly repetitive regions, thus providing the maximum reachable recall for both single nucleotide variants and small insertion and deletions and to precisely guide the user in the pipeline optimization process.

AVAILABILITY AND IMPLEMENTATION

Source code is freely available under MIT license at https://github.com/mazzalab-ieo/recallme. RecallME web application is available at https://translational-oncology-lab.shinyapps.io/recallme/. To use RecallME, users must obtain a license for ANNOVAR by themselves.

Exome analysis focusing on epilepsy-related genes in children and adults with sudden unexplained death

PURPOSE

Genetic studies in sudden infant death syndrome (SIDS) and sudden unexplained death (SUD) cohorts have indicated that cardiovascular diseases might have contributed to sudden unexpected death in 20-35 % of autopsy-negative cases. Sudden unexpected death can also occur in people with epilepsy, termed as sudden unexpected death in epilepsy (SUDEP). The pathophysiological mechanisms of SUDEP are not well understood, but are likely multifactorial, including seizure-induced hypoventilation and arrhythmias as well as genetic risk factors. The sudden death of some of the SIDS/SUD victims might also be explained by genetic epilepsy, therefore this study aimed to expand the post-mortem genetic analysis of SIDS/SUD cases to epilepsy-related genes.

METHODS

Existing whole-exome sequencing data from our 155 SIDS and 45 SUD cases were analyzed, with a focus on 365 epilepsy-related genes. Nine of the SUD victims had a known medical history of epilepsy, seizures or other underlying neurological conditions and were therefore classified as SUDEP cases.

RESULTS

In our SIDS and SUD cohorts, we found epilepsy-related pathogenic/likely pathogenic variants in the genes OPA1, RAI1, SCN3A, SCN5A and TSC2.

CONCLUSION

Post-mortem analysis of epilepsy-related genes identified potentially disease-causing variants that might have contributed to the sudden death events in our SIDS/SUD cases. However, the interpretation of identified variants remains challenging and often changes over time as more data is gathered. Overall, this study contributes insight in potentially pathophysiological epilepsy-related mechanisms in SIDS, SUD and SUDEP victims and underlines the importance of sensible counselling on the risk and preventive measures in genetic epilepsy.

A Novel Truncating Mutation in PAX1 Gene Causes Otofaciocervical Syndrome Without Immunodeficiency

Otofaciocervical syndrome (OTFCS) is a rare genetic disorder of both autosomal recessive and autosomal dominant patterns of inheritance. It is caused by biallelic or monoallelic mutations in PAX1 or EYA1 genes, respectively. Here, we report an OTFCS2 female patient of 1st consanguineous healthy parents. She manifested facial dysmorphism, hearing loss, intellectual disability (ID), and delayed language development (DLD) as the main clinical phenotype. The novel homozygous variant c.1212dup (p.Gly405Argfs*51) in the PAX1 gene was identified by whole exome sequencing (WES), and family segregation confirmed the heterozygous status of the mutation in the parents using the Sanger sequencing. The study recorded a novel PAX1 variant representing the sixth report of OTFCS2 worldwide and the first Egyptian study expanding the geographic area where the disorder was confined.

Informing a value care model: lessons from an integrated adult neurogenomics clinic

BACKGROUND

Advances in genomics provide improved opportunities for diagnosis of complex neurogenetic disorders, yet the optimal approach to translate these benefits to the outpatient clinic is unclear.

AIMS

We retrospectively reviewed referral indications and outcomes of an integrated multidisciplinary team (MDT) clinic pathway for adults with suspected neurogenetic disorders. The associated cost implications were estimated.

METHODS

Consecutive patients who attended the neurogenomics clinic from January 2017 to April 2020 were included. The clinic comprised neurologists, clinical geneticists and genetic counsellors, who assessed each patient concurrently.

RESULTS

Ninety-nine new patients were referred spanning 45 different clinical diagnoses. Following MDT clinical assessment, 23% (23/99) of referral diagnoses were revised prior to molecular testing. Eighty-one patients (82%) underwent genetic testing, including 43 exome-based panels, 15 whole-genome sequencing, 14 single gene tests, 27 repeat-primed polymerase chain reaction testing and two chromosomal microarrays. Overall, 33/99 patients (33%) received a diagnosis, either a molecular diagnosis (n = 24, of which 22 were diagnostic and two were predictive) or a clinical diagnosis (n = 9). Of the clinical diagnosis cohort, five patients received a diagnosis without molecular testing and four patients whose negative testing (one diagnostic and three predictive) allowed exclusion of genetic differentials and, hence, confirmation of clinical diagnoses. The diagnostic rate following MDT and diagnostic testing was 30% (28/94), excluding the five predictive testing cases. MDT assessment aligned with eventual molecular diagnoses in 96% of cases. The estimated average costs were AU$1386 per patient undergoing MDT assessment and AU$4159 per diagnosis achieved.

CONCLUSIONS

We present an integrated multidisciplinary neurogenomics clinic pathway providing a diagnostic yield of 33% (30% excluding predictive testing cases), with costing implications. The relatively high diagnostic yield may be attributed to multidisciplinary input integrating accurate phenotyping of complex disorders and interpretation of genomic findings.