Analysis Of Whole Exome Sequencing Data Research Articles

Copy Number Analysis of Whole Exome Sequencing Data

Abstract Background Genetic alterations implicated in constitutional disorders range in size from single nucleotide substitutions to losses and gains of millions of bases. Whole exome sequencing (WES) captures many clinically relevant variants, but an analytic pipeline optimized for detection of small substitutions, insertions, and deletions may miss pathogenic copy number variants (CNVs). Historically, our laboratory performed separate cytogenetic tests to evaluate for larger genomic imbalances. Addition of copy number analysis to our WES pipeline could improve the assay’s diagnostic yield. This study describes part of our strategy to validate copy number analysis of our constitutional whole exome sequencing data. Methods We identified 18 CNV specimens that underwent both chromosomal microarray analysis (CMA) and WES between 2017 and 2023. Variants included 13 deletions (size range 14Kb-2Mb) and 5 duplications (size range 325Kb-10Mb). The standard method of constitutional copy number analysis in our clinical laboratory is CMA using the Affymetrix Cytoscan HD Array and analysis software. CNVs identified on CMA were reviewed with reference to the Database of Genomic Variants and multiple academic hospital databases. WES was performed on an Illumina sequencing system following capture-based library enrichment with the Integrated DNA Technologies xGen Exome Research, CNV Backbone, Human mtDNA Research, and Human ID Research Panel probes. CNVs identified on WES were called using the Broad Institute’s Genomic Analysis Toolkit (GATK) with a quality score cutoff of 120. Whole exome sequencing data was visualized on the Integrative Genomics Viewer. Results 14 out of 18 CNVs reported on CMA were also detected on WES with a quality score at or above 120 (78%). This includes 10 out of 13 deletions (77%) and 3 out of 5 duplications (60%). For both deletions and duplications, the size ranges of detected variants at or above the cutoff score and undetected variants overlapped. One duplication reported on CMA was detected on WES with a quality score below 120. Visualization of whole exome sequencing data showed variation in probe density that may have affected CNV calling. Conclusions The majority of CNVs reported after CMA (14/18, 78%) were also detected on WES data using a GATK quality score cutoff of 120. The data does not suggest an association between CNV size and detection by WES, although interpretation is limited by the small number of cases. Data visualization may be required to identify some CNVs. Future steps to characterize the limitations of copy number analysis of WES data include assessments of variance in probe density and read depth across the exome. Copy number analysis of WES data from cases with normal CMA results will also be performed for assay validation.

Genetic link between KIF1A mutations and amyotrophic lateral sclerosis: evidence from whole-exome sequencing.

Genetics have been shown to have a substantial impact on amyotrophic lateral sclerosis (ALS). The ALS process involves defects in axonal transport and cytoskeletal dynamics. It has been identified that KIF1A, responsible for encoding a kinesin-3 motor protein that carries synaptic vesicles, is considered a genetic predisposing factor for ALS. The analysis of whole-exome sequencing data from 1,068 patients was conducted to examine the genetic link between ALS and KIF1A. For patients with KIF1A gene mutations and a family history, we extended the analysis to their families and reanalyzed them using Sanger sequencing for cosegregation analysis. In our cohort, the KIF1A mutation frequency was 1.31% (14/1,068). Thirteen nonsynonymous variants were detected in 14 ALS patients. Consistent with the connection between KIF1A and ALS, the missense mutation p.A1083T (c.3247G>A) was shown to cosegregate with disease. The mutations related to ALS in our study were primarily located in the cargo-binding region at the C-terminal, as opposed to the mutations of motor domain at the N-terminal of KIF1A which were linked to hereditary peripheral neuropathy and spastic paraplegia. We observed high clinical heterogeneity in ALS patients with missense mutations in the KIF1A gene. KIF5A is a more frequent determinant of ALS in the European population, while KIF1A accounts for a similar proportion of ALS in both the European and Chinese populations. Our investigation revealed that mutations in the C-terminus of KIF1A could increase the risk of ALS, support the pathogenic role of KIF1A in ALS and expand the phenotypic and genetic spectrum of KIF1A-related ALS.

FBLN2 is associated with Goldenhar syndrome and is essential for cranial neural crest cell development.

Goldenhar syndrome, a rare craniofacial malformation, is characterized by developmental anomalies in the first and second pharyngeal arches. Its etiology is considered to be heterogenous, including both genetic and environmental factors that remain largely unknown. To further elucidate the genetic cause in a five-generation Goldenhar syndrome pedigree and exploit the whole-exome sequencing (WES) data of this pedigree, we generated collapsed haplotype pattern markers based on WES and employed rare variant nonparametric linkage analysis. FBLN2 was identified as a candidate gene via analysis of WES data across the significant linkage region. A fbln2 knockout zebrafish line was established by CRISPR/Cas9 to examine the gene's role in craniofacial cartilage development. fbln2 was expressed specifically in the mandible during the zebrafish early development, while fbln2 knockout zebrafish exhibited craniofacial malformations with abnormal chondrocyte morphologies. Functional studies revealed that fbln2 knockout caused abnormal chondrogenic differentiation, apoptosis, and proliferation of cranial neural crest cells (CNCCs), and downregulated the bone morphogenic protein (BMP) signaling pathway in the zebrafish model. This study demonstrates the role of FBLN2 in CNCC development and BMP pathway regulation, and highlights FBLN2 as a candidate gene for Goldenhar syndrome, which may have implications for the selection of potential screening targets and the development of treatments for conditions like microtia-atresia.

Discovery of Pathogenic Variants Associated with Idiopathic Recurrent Pregnancy Loss Using Whole-Exome Sequencing.

Idiopathic recurrent pregnancy loss (RPL) is defined as at least two pregnancy losses before 20 weeks of gestation. Approximately 5% of pregnant couples experience idiopathic RPL, which is a heterogeneous disease with various causes including hormonal, chromosomal, and intrauterine abnormalities. Although how pregnancy loss occurs is still unknown, numerous biological factors are associated with the incidence of pregnancy loss, including genetic variants. Whole-exome sequencing (WES) was conducted on blood samples from 56 Korean patients with RPL and 40 healthy controls. The WES data were aligned by means of bioinformatic analysis, and the detected variants were annotated using machine learning tools to predict the pathogenicity of protein alterations. Each indicated variant was confirmed using Sanger sequencing. A replication study was also conducted in 112 patients and 114 controls. The Variant Effect Scoring Tool, Combined Annotation Dependent Depletion tool, Sorting Intolerant from Tolerant annotation tool, and various databases detected 10 potential variants previously associated with spontaneous abortion genes in patients by means of a bioinformatic analysis of WES data. Several variants were detected in more than one patient. Interestingly, several of the detected genes were functionally clustered, including some with a secretory function (mucin 4; MUC4; rs200737893 G>A and hyaluronan-binding protein 2; HABP2; rs542838125 G>T), in which growth arrest-specific 2 Like 2 (GAS2L2; rs140842796 C>T) and dynamin 2 (DNM2; rs763894364 G>A) are functionally associated with cell protrusion and the cytoskeleton. ATP Binding Cassette Subfamily C Member 6 (ABCC6) was the only gene with two variants. HABP2 (rs542838125 G>T), MUC4 (rs200737893 G>A), and GAS2L2 (rs140842796 C>T) were detected in only the patient group in the replication study. The combination of WES and machine learning tools is a useful method to detect potential variants associated with RPL. Using bioinformatic tools, we found 10 potential variants in 9 genes. WES data from patients are needed to better understand the causes of RPL.

Loss of GPR75 protects against non-alcoholic fatty liver disease and body fat accumulation

Approximately 1 in 4 people worldwide have non-alcoholic fatty liver disease (NAFLD); however, there are currently no medications to treat this condition. This study investigated the role of adiposity-associated orphan G protein-coupled receptor 75 (GPR75) in liver lipid accumulation. We profiled Gpr75 expression and report that it is most abundant in the brain. Next, we generated the first single-cell-level analysis of Gpr75 and identified a subpopulation co-expressed with key appetite-regulating hypothalamic neurons. CRISPR-Cas9-deleted Gpr75 mice fed a palatable western diet high in fat adjusted caloric intake to remain in energy balance, thereby preventing NAFLD. Consistent with mouse results, analysis of whole-exome sequencing data from 428,719 individuals (UK Biobank) revealed that variants in GPR75 are associated with a reduced likelihood of hepatic steatosis. Here, we provide a significant advance in understanding of the expression and function of GPR75, demonstrating that it is a promising pharmaceutical target for NAFLD treatment.

Comorbidity of bathing suit ichthyosis and limb-girdle muscular dystrophy type 2 A in a Tunisian patient revealed by Whole Exome Sequencing

Elevated rates of consanguinity and inbreeding are responsible for the high prevalence of recessively inherited diseases among inbred populations including Tunisia. In addition, the co-occurrence of two of these conditions, called also comorbidity, within the same individual or in members of the same family are often described in Tunisia which is challenging for diagnosis. The high throughput sequencing has improved the diagnosis of inherited diseases. We report here on a 32-year-old woman born to consanguineous parents. She presented with congenital ichthyosis and muscular dystrophy. She was primarily suspected as suffering from Chanarin-Dorfman syndrome (CDS) with unusual form. Screening of founder mutations allowed only the elucidation of the molecular etiology of Ichthyosis. As the result was inconclusive, Whole Exome Sequencing (WES) was conducted. WES data analysis led to the identification of a mutation in the CAPN3 gene underlying limb-girdle muscular dystrophy type 2A (LGMD2A). Sanger sequencing confirmed the familial segregation of mutations. This work presents the first case worldwide of individual comorbidity of bathing suit ichthyosis and LGMD2A. The co-occurrence of two diseases should be systematically considered when establishing a diagnosis especially in consanguineous populations. WES is a powerful tool for molecular diagnosis in particular for revealing comorbidities and rectifying the diagnosis.

Discovery of a Novel Shared Variant Among RTEL1 Gene and RTEL1-TNFRSF6B lncRNA at Chromosome 20q13.33 in Familial Progressive Myoclonus Epilepsy.

Background: Progressive myoclonus epilepsy (PME) is a neurodegenerative disorder marked by recurrent seizures and progressive myoclonus. To date, based on the phenotypes and causal genes, more than 40 subtypes of PMEs have been identified, and more remain to be characterized. Our study is aimed at identifying the aberrant gene(s) possibly associated with PMEs in two siblings born to asymptomatic parents, in the absence of known genetic mutations. Methods: Clinical assessments and molecular analyses, such as the repeat expansion test for CSTB; SCA1, 2, 3, 6, and 7; whole exome sequencing (WES); and mitochondrial genome sequencing coupled with computational analysis, were performed. Results: A family-based segregation analysis of WES data was performed to identify novel genes associated with PMEs. The potassium channel, KCNH8 [c.298T>C; (p.Tyr100His)], a DNA repair gene, regulator of telomere elongation helicase 1 (RTEL1) [c.691G>T; (p.Asp231Tyr)] and long noncoding RNA, RTEL1-TNFRSF6B [chr20:62298898_G>T; NR_037882.1, hg19] were among the candidate genes that were found to be associated with PMEs. These homozygous variations in siblings belong to genes with a loss-of-function intolerant (pLI) score of ≤ 0.86, expected to be detrimental by multiple computational analyses, and were heterozygous in parents. Additionally, computational analysis and the expression of RTEL1 and RTEL1-TNFRSF6B revealed that RTEL1-TNFRSF6B may modulate RTEL1 via hsa-miR-3529-3p. In the patient with the severe phenotype, a further deleterious mutation in SLC22A17 was identified. No de novo variants specific to these probands were identified in the mitochondrial genome. Conclusions: Our study is the first to report variants in KCNH8, RTEL1, and RTEL1-TNFRSF6B among PME cases. These genes when characterized fully may shed light on pathogenicity and have the potential to be used in the diagnosis of PME.

Identification and bioinformatics analysis of a novel mutation in PLA2G6 gene in a patient with neurodegenerative disorder

PLA2G6-associated neurodegeneration disorders are a group of conditions that result from mutations in the PLA2G6 gene. These disorders can cause a range of symptoms affecting the brain and nervous system, such as movement problems, cognitive decline, and vision loss. In this case, Whole Exome Sequencing (WES) was conducted on a 9-year-old boy born to consanguineous parents, who have an increased risk of passing on recessive genetic disorders. Although the boy did not exhibit any clinical symptoms until he was seven months old, his condition deteriorated severely by the age of two, rendering him unable to sit or move around, and he did not develop the ability to speak. DNAs were extracted from the blood samples of the proband (the 9-year-old boy) and his parents. The proband's DNA was sequenced using WES to identify potential disease-causing mutations. To confirm the WES findings, PCR-Sanger sequencing analysis was performed on the DNAs obtained from the proband (the 9-year-old boy) and his parents. The effect of the mutation on the protein structure was studied using the I-TASSER web-tool, which uses computational methods to predict protein structure and function based on amino acid sequences.The analysis of WES data revealed a homozygous dinucleotide TG insertion at position c.1362–1363 in exon 10 of the PLA2G6 gene (NC_000022.11(NM_003560.4):c.1362_1363insTG). This insertion introduced a premature stop codon at position 455 (p.Met455Ter) in the protein sequence. The validation of the WES results using PCR-Sanger sequencing confirmed the presence of the mutation in the proband and revealed that both parents were carriers of the mutation. The premature stop codon, located 65 nucleotides before the exon-junction complex, is likely to be detected by the nonsense-mediated mRNA decay (NMD) pathway, leading to the degradation of the mRNA. The use of the I-TASSER web-tool further supported the findings by demonstrating significant structural and functional changes in the mutant protein. Observing a very severe phenotype in our case supports the idea that the mutation was highly deleterious.

Whole exome sequencing highlights rare variants in CTCF, DNMT1, DNMT3A, EZH2 and SUV39H1 as associated with FSHD.

Introduction: Despite the progress made in the study of Facioscapulohumeral Dystrophy (FSHD), the wide heterogeneity of disease complicates its diagnosis and the genotype-phenotype correlation among patients and within families. In this context, the present work employed Whole Exome Sequencing (WES) to investigate known and unknown genetic contributors that may be involved in FSHD and may represent potential disease modifiers, even in presence of a D4Z4 Reduced Allele (DRA). Methods: A cohort of 126 patients with clinical signs of FSHD were included in the study, which were characterized by D4Z4 sizing, methylation analysis and WES. Specific protocols were employed for D4Z4 sizing and methylation analysis, whereas the Illumina® Next-Seq 550 system was utilized for WES. The study included both patients with a DRA compatible with FSHD diagnosis and patients with longer D4Z4 alleles. In case of patients harboring relevant variants from WES, the molecular analysis was extended to the family members. Results: The WES data analysis highlighted 20 relevant variants, among which 14 were located in known genetic modifiers (SMCHD1, DNMT3B and LRIF1) and 6 in candidate genes (CTCF, DNMT1, DNMT3A, EZH2 and SUV39H1). Most of them were found together with a permissive short (4-7RU) or borderline/long DRA (8-20RU), supporting the possibility that different genes can contribute to disease heterogeneity in presence of a FSHD permissive background. The segregation and methylation analysis among family members, together with clinical findings, provided a more comprehensive picture of patients. Discussion: Our results support FSHD pathomechanism being complex with a multigenic contribution by several known (SMCHD1, DNMT3B, LRIF1) and possibly other candidate genes (CTCF, DNMT1, DNMT3A, EZH2, SUV39H1) to disease penetrance and expressivity. Our results further emphasize the importance of extending the analysis of molecular findings within the proband's family, with the purpose of providing a broader framework for understanding single cases and allowing finer genotype-phenotype correlations in FSHD-affected families.

A novel loss of function mutation in the PHD domain of the RAG2 gene, affecting zinc-binding affinity.

Severe combined immune deficiencies (SCIDs) are genetically heterogeneous disorders that lead to the absence or malfunction of adaptive immune cells, including T- and B-cells. Pathogenic variants in the RAG2 gene are associated with this disease. A couple with consanguineous marriage from the Iranian-Azeri-Turkish ethnic group was referred to the genetic lab. Two children of this family died due to SCID disease with symptoms of skin granulomas, lack of developed T- and B-cells, and intact NK cells. To infer their genotypes, DNA samples obtained from the parents were subjected to whole-exome sequencing (WES). WES data analysis revealed that both parents were carriers of a pathogenic variant, NC_000011.10 (NM_000536.4):c.1268G > C, in the RAG2 gene. This variant was absent in our cohort of 400 healthy individuals from the same ethnic group. To gain insight into the consequence of the variant on the protein function, further analysis was performed by applying bioinformatics tools. This study revealed that the replacement of cysteine with serine at the zinc-binding domain diminished the domain's affinity to zinc ion, resulting in the loss of the mutant protein's ability to bind to the recombination signal sequence (RSS). The formation of the RAG2-RSS complex is vital for T- and B-cell development. The identification of a novel pathogenic variant, c.1268G > C, revealed that this variant in the zinc-binding domain diminished the affinity of the zinc ion to the mutant protein and consequently led to the loss of its ability to bind to the RSS.

Var∣Decrypt: a novel and user-friendly tool to explore and prioritize variants in whole-exome sequencing data

BackgroundHigh-throughput sequencing (HTS) offers unprecedented opportunities for the discovery of causative gene variants in multiple human disorders including cancers, and has revolutionized clinical diagnostics. However, despite more than a decade of use of HTS-based assays, extracting relevant functional information from whole-exome sequencing (WES) data remains challenging, especially for non-specialists lacking in-depth bioinformatic skills.ResultsTo address this limitation, we developed Var∣Decrypt, a web-based tool designed to greatly facilitate WES data browsing and analysis. Var∣Decrypt offers a wide range of gene and variant filtering possibilities, clustering and enrichment tools, providing an efficient way to derive patient-specific functional information and to prioritize gene variants for functional analyses. We applied Var∣Decrypt on WES datasets of 10 acute erythroid leukemia patients, a rare and aggressive form of leukemia, and recovered known disease oncogenes in addition to novel putative drivers. We additionally validated the performance of Var∣Decrypt using an independent dataset of ~ 90 multiple myeloma WES, recapitulating the identified deregulated genes and pathways, showing the general applicability and versatility of Var∣Decrypt for WES analysis.ConclusionDespite years of use of WES in human health for diagnosis and discovery of disease drivers, WES data analysis still remains a complex task requiring advanced bioinformatic skills. In that context, there is a need for user-friendly all-in-one dedicated tools for data analysis, to allow biologists and clinicians to extract relevant biological information from patient datasets. Here, we provide Var∣Decrypt (trial version accessible here: https://vardecrypt.com/app/vardecrypt), a simple and intuitive Rshiny application created to fill this gap. Source code and detailed user tutorial are available at https://gitlab.com/mohammadsalma/vardecrypt.

Universal cutoff for tumor mutational burden in predicting the efficacy of anti-PD-(L)1 therapy for advanced cancers.

Background: The US Food and Drug Administration (FDA)'s tumor-agnostic approval of pembrolizumab in high tumor mutational burden (TMB-high, i.e., TMB≥10 mut/Mb) cases, based on the data from KEYNOTE-158, has raised considerable concerns among the immuno-oncology community. This study aims to statistically infer the optimal universal cutoff in defining TMB-high that is predictive of the efficacy of anti-PD-(L) 1 therapy in advanced solid tumors. Methods: We integrated MSK-IMPACT TMB data from a public cohort and the objective response rate (ORR) for anti-PD-(L) 1 monotherapy across diverse cancer types in published trials. The optimal TMB cutoff was determined by varying the universal cutoff to define TMB-high across cancer types and examining the cancer-level correlation between objective response rate and the proportion of TMB-high cases. The utility of this cutoff in predicting overall survival (OS) benefits from anti-PD-(L) 1 therapy was then evaluated in a validation cohort of advanced cancers with coupled MSK-IMPACT TMB and OS data. In silico analysis of whole-exome sequencing data from The Cancer Genome Atlas was further employed to assess the generalizability of the identified cutoff among panels comprising several hundred genes. Results: The cancer type-level analysis identified 10 mut/Mb as the optimal cutoff for MSK-IMPACT in defining TMB-high, with the corresponding TMB-high (TMB≥10 mut/Mb) percentage strongly correlated with ORR for PD-(L) 1 blockade across cancer types [correlation coefficient, 0.72 (95% CI, 0.45-0.88)]. This cutoff was also the optimum in defining TMB-high (via MSK-IMPACT) when predicting OS benefits from anti-PD-(L) 1 therapy in the validation cohort. In this cohort, TMB≥10 mut/Mb was associated with significantly improved OS (hazard ratio, 0.58 [95% CI, 0.48-0.71]; p < 0.001). Moreover, in silico analyses revealed excellent agreement of TMB≥10 mut/Mb cases between MSK-IMPACT and the FDA-approved panels and between MSK-IMPACT and various randomly sampled panels. Conclusion: Our study demonstrates that 10 mut/Mb is the optimal, universal cutoff for TMB-high that guides the clinical application of anti-PD-(L) 1 therapy for advanced solid tumors. It also provides rigorous evidence beyond KEYNOTE-158 for the utility of TMB≥10 mut/Mb in predicting the efficacy of PD-(L) 1 blockade in broader settings, which could help to mitigate the challenges in embracing the tumor-agnostic approval of pembrolizumab in TMB-high cases.

Whole exome sequencing analysis and prenatal diagnosis in children with neurodevelopmental disorders

To explore the application value of whole exome sequencing (WES) in the diagnosis of prenatal and postnatal neurodevelopmental disorders (NDDs). A total of 70 patients diagnosed with NDDs who underwent WES at the Medical Genetics Center of the Maternal and Child Health Hospital of Hubei Province between June 2020 and July 2021 were retrospectively analyzed. Genomic DNA was extracted from peripheral blood samples and amniotic fluid. WES-based copy number variant (CNV) analysis was integrated into the routine WES data analysis pipeline. The results showed that a molecular diagnosis rate could be made in 21/70 (30%) cases. Of 21 positive cases, 14 (23%) cases were detected by single-nucleotide variant/small insertion/deletion (SNV/Indel) analysis, of which 12 variants were novel, 6 (9.8%) cases were detected by WES-based CNV analysis, and 1 (1.6%) case was detected by a combination of both. The diagnostic yield of WES combined with CNV analysis was higher than that of SNV/Indel analysis alone (30%, 21/70 vs. 20%, 14/70). Of the 28 prenatally diagnosed cases, 6 cases were found to have inherited parental variation for NDDs, 10 cases were found not to have the same pathogenic variation as the proband, and the remaining 12 cases were found to have no pathogenic or likely pathogenic variation that could explain the NDDs phenotype. Clinical follow-up showed that 5 families opted for abortion and the remaining had no current abnormalities. In conclusion, WES may be an effective method to clarify the genetic etiology and prenatal diagnosis of NDDs, which is helpful in assessing the prognosis to aid clinical management and reproductive guidance.

WWC2 expression in the testis: Implications for spermatogenesis and male fertility.

The family of WWC proteins is known to regulate cell proliferation and organ growth control via the Hippo signaling pathway. As WWC proteins share a similar domain structure and a common set of interacting proteins, they are supposed to fulfill compensatory functions in cells and tissues. While all three WWC family members WWC1, WWC2, and WWC3 are found co-expressed in most human organs including lung, brain, kidney, and liver, in the testis only WWC2 displays a relatively high expression. In this study, we investigated the testicular WWC2 expression in spermatogenesis and male fertility. We show that the Wwc2 mRNA expression level in mouse testes is increased during development in parallel with germ cell proliferation and differentiation. The cellular expression of each individual WWC family member was evaluated in published single-cell mRNA datasets of murine and human testes demonstrating a high WWC2 expression predominantly in early spermatocytes. In line with this, immunohistochemistry revealed cytosolic WWC2 protein expression in primary spermatocytes from human testes displaying full spermatogenesis. In accordance with these findings, markedly lower WWC2 expression levels were detected in testicular tissues from mice and men lacking germ cells. Finally, analysis of whole-exome sequencing data of male patients affected by infertility and unexplained severe spermatogenic failure revealed several heterozygous, rare WWC2 gene variants with a proposed damaging function and putative impact on WWC2 protein structure. Taken together, our findings provide novel insights into the testicular expression of WWC2 and show its cell-specific expression in spermatocytes. As rare WWC2 variants were identified in the background of disturbed spermatogenesis, WWC2 may be a novel candidate gene for male infertility.

Developmental and epileptic encephalopathy 89: A novel bi-allelic variant, molecular dynamics simulation, and a comprehensive clinical and molecular profile.

Gamma-aminobutyric acid (GABA), the major inhibitory neurotransmitter in the adult central nervous system, plays an important role during embryonic neural network formation. GAD67 is the rate-limiting enzyme in GABA synthesis, and its deficiency leads to developmental and epileptic encephalopathy 89 (DEE 89). Patients who suffered from this syndrome generally manifested severe to profound neurodevelopmental delay, seizures, and often congenital anomalies such as the cleft palate or/and omphalocele. Up to now, only three papers on this syndrome have been published, and our knowledge about the disease's clinical course and pathophysiology is in its infancy. We used whole-exome sequencing (WES) and multiple in-silico tools to detect a potential causative variant in a patient with severe neurodevelopmental delay and refractory epilepsy. Moreover, by molecular docking and molecular dynamics simulation, we investigate the effect of the candidate variant on the GAD67 function and structure. WES data analysis revealed a novel deleterious variant (NM_000817.3: c.850C>T; p.Leu284Phe) in the GAD1 gene, which encodes the GAD67 enzyme. Molecular docking and molecular dynamics simulation showed that this variant has deleterious effects on the structure and function of the GAD67. This study's patient, in addition to typical symptoms of the DEE89, showed microcephaly and clonus in the toe, which were novel clinical findings. Our findings expand the mutational and clinical spectrum of DEE 89. Also, by gathering clinical symptoms and genetic findings of previously reported cases, moreover providing a comprehensive clinical picture of the disease, we found that there was no common drug therapy among patients whose epilepsy was controlled. Furthermore, the comparison of clinical symptoms between patients with missense and truncating mutations did not show any significant clinical difference, except that patients with missense mutations did not show cleft palates or omphaloceles.

Determination of the Cancer Genome Atlas (TCGA) Endometrial Cancer Molecular Subtypes Using the Variant Interpretation and Clinical Decision Support Software MH Guide

The Cancer Genome Atlas (TCGA) network (United States National Cancer Institute) identified four molecular endometrial cancer (EC) subtypes using an extensive multi-method approach. The aim of this study was to determine the four TCGA EC molecular subtypes using a single-method whole-exome sequencing (WES)-based approach provided by MH Guide (Molecular Health, Heidelberg, Germany). WES and clinical data of n = 232 EC patients were obtained from TCGA. The four TCGA EC molecular subtypes designated as (i) Mutated Polymerase ε (POLE), (ii) Microsatellite Instability (MSI), (iii) Copy Number (CN) low and, (iv) CN-high were determined using the MH Guide software. The prognostic value of the subtypes determined by MH Guide were compared with the TCGA classification. Analysis of WES data using the MH Guide software led to the precise identification of the four EC molecular subtypes analogous to the TCGA classification. Both approaches displayed high concordance in terms of prognostic significance. The multi-method-based TCGA EC molecular subtypes can reliably be reproduced by the single-method-based MH Guide approach. The easy-to-implement single-method MH Guide approach represents a promising diagnostic tool.

Impact of Genomic and Clinical Factors on Outcome of Children ≥18 Months of Age with Stage 3 Neuroblastoma with Unfavorable Histology and without MYCN Amplification: A Children's Oncology Group (COG) Report.

Patients ≥18 months of age with International Neuroblastoma Staging System (INSS) stage 3 unfavorable histology (UH), MYCN-nonamplified (MYCN-NA) tumors have favorable survival rates compared with other high-risk neuroblastoma populations. The impact of select clinical and biological factors on overall survival (OS) and event-free survival (EFS) were evaluated. Patients enrolled on Children's Oncology Group (COG) A3973 (n = 34), ANBL0532 (n = 27), and/or biology protocol ANBL00B1 (n = 72) were analyzed. Tumors with available DNA (n = 65) and RNA (n = 42) were subjected to whole-exome sequencing (WES) and RNA sequencing. WES analyses and gene expression profiling were evaluated for their impact on survival. Multivariate analyses of EFS/OS using significant factors from univariate analyses were performed. 5-year EFS/OS for patients treated with high-risk therapy on A3973 and ANBL0532 were 73.0% ± 8.1%/87.9% ± 5.9% and 61.4% ± 10.2%/73.0% ± 9.2%, respectively (P = 0.1286 and P = 0.2180). In the A3973/ANBL0532 cohort, patients with less than partial response (PR; n = 5) at end-induction had poor outcomes (5-year EFS/OS: 0%/20.0% ± 17.9%. Univariate analyses of WES data revealed that subjects whose tumors had chromosome 1p or 11q loss/LOH and chromosome 5 or 9 segmental chromosomal aberrations had inferior EFS compared with those with tumors without these aberrations. Multivariate analysis revealed that 11q loss/LOH was an independent predictor of inferior OS [HR, 3.116 (95% confidence interval, 1.034-9.389), P = 0.0435]. Patients ≥18 months of age at diagnosis who had tumors with UH and MYCN-NA INSS stage 3 neuroblastoma assigned to high-risk therapy had an 81.6% ± 5.3% 5-year OS. Less than PR to induction therapy and chromosome 11q loss/LOH are independent predictors of inferior outcome and identify patients who should be eligible for future high-risk clinical trials.

Genetic Variant Overlap Analysis Identifies Established and Putative Genes Involved in Pulmonary Fibrosis.

In only around 40% of families with pulmonary fibrosis (PF) a suspected genetic cause can be found. Genetic overlap analysis of Whole Exome Sequencing (WES) data may be a powerful tool to discover new shared variants in novel genes for PF. As a proof of principle, we first selected unrelated PF patients for whom a genetic variant was detected (n = 125) in established PF genes and searched for overlapping variants. Second, we performed WES (n = 149) and identified novel potentially deleterious variants shared by at least two unrelated PF patients. These variants were genotyped in validation cohorts (n = 2748). In 125 unrelated patients, a potentially deleterious variant was detected in known PF genes of which 15 variants in six genes overlapped, involving 51 patients. Overlap analysis of WES data identified two novel variants of interest: TOM1L2 c.421T > C p.(Y141H) and TDP1c.1373dupG p.(S459fs*5), neither gene had been related to pulmonary fibrosis before. Both proteins were present in the alveolar epithelium. No apparent characteristics of telomere disease were observed. This study underlines the potential of searching for overlapping rare potentially deleterious variants to identify disease-associated variants and genes. A previously unreported variant was found in two putative new PF genes, but further research is needed to determine causality.

Two Distinct Deleterious Causative Variants in a Family with Multiple Cancer-Affected Patients.

Only 5 to 10% of cancers are hereditary, but they are particularly important since they can be passed down from generation to generation, and family members are at elevated risk. Although screening methods are one of the essential strategies for dealing with hereditary cancers, they do not have high specificity and sensitivity. The emergence of whole-exome sequencing (WES) causes a significant increase in the diagnostic rate of cancer-causing variants in at-risk families. We performed WES on the proband's DNA sample from an Iranian family with multiple cancer-affected members to identify potential causative variants. Multiple in silico tools were used to evaluate the candidate variants' pathogenicity and their effects on the protein's structure, function, and stability. Moreover, the candidate variants were co-segregated in the family with Sanger sequencing. The WES data analysis identified two pathogenic variants (CHEK2: NM_007194.4: c.538C>T, p.Arg180Cys and MLH1: NM_000249.4, c.844G>A, p.Ala282Thr). Sanger sequencing data showed each of the variants was incompletely segregated with phenotype, but both of them explained the patient's phenotype together. Also, the structural analysis demonstrated that due to the variant (c.538C>T), a salt bridge between arginine 180 and glutamic acid 149 was lost. Indeed, several protein stability tools described both variants as destabilizing. Herein, we interestingly identify two distinct deleterious causative variants (CHEK2: NM_007194.4: c.538C>T, p.Arg180Cys and MLH1: NM_000249.4, c.844G>A, p.Ala282Thr) in a family with several cancer-affected members. Furthermore, this study's findings established the utility of WES in the genetic diagnostics of cancer.

Whole-Exome Sequencing of Pakistani Consanguineous Families Identified Pathogenic Variants in Genes of Intellectual Disability.

Intellectual disability (ID) is a condition of significant limitation of cognitive functioning and adaptive behavior, with 50% of etiology attributed to genetic predisposition. We recruited two consanguineous Pakistani families manifesting severe ID and developmental delay. The probands were subjected to whole exome sequencing (WES) and variants were further prioritized based on population frequency, predicted pathogenicity and functional relevance. The WES data analysis identified homozygous pathogenic variants in genes MBOAT7 and TRAPPC9. The pathogenicity of the variants was supported by co-segregation analysis and in silico tool. The findings of this study expand mutation spectrum and provide additional evidence to the role of MBOAT7 and TRAPPC9 in causation of ID.