Search

Background: The aim of this study was to evaluate the diagnostic yield of routine exome sequencing (ES) in fetuses with ultrasound anomalies., Methods: We performed a retrospective analysis of the ES results of 629 fetuses with isolated or multiple anomalies referred in 2019-2022. Variants in a gene panel consisting of approximately 3400 genes associated with multiple congenital anomalies and/or intellectual disability were analyzed. We used trio analysis and filtering for de novo variants, compound heterozygous variants, homozygous variants, X-linked variants, variants in imprinted genes, and known pathogenic variants., Results: Pathogenic and likely pathogenic variants (class five and four, respectively) were identified in 14.0% (88/629, 95% CI 11.5%-16.9%) of cases. In the current cohort, the probability of detecting a monogenetic disorder was ∼1:7 (88/629, 95% CI 1:8.7-1:5.9), ranging from 1:9 (49/424) in cases with one major anomaly to 1:5 (32/147) in cases with multiple system anomalies., Conclusions: Our results indicate that a notable number of fetuses (1:7) with ultrasound anomalies and a normal chromosomal microarray have a (likely) pathogenic variant that can be detected through prenatal ES. These results warrant implementation of exome sequencing in selected cases, including those with an isolated anomaly on prenatal ultrasound., (© 2024 The Author(s). Prenatal Diagnosis published by John Wiley & Sons Ltd.)

Objective: To evaluate the diagnostic yield of exome sequencing (ES) in fetuses and neonates with prenatally detected congenital diaphragmatic hernia (CDH) and normal copy number variant (CNV) analysis., Methods: We conducted a retrospective cohort study of prenatally diagnosed CDH cases seen between 2019 and 2022. All cases who underwent prenatal or postnatal genetic testing were reviewed. The results from the ES analysis that identified pathogenic or likely pathogenic single nucleotide variants are described., Results: In total, 133 fetuses with CDH were seen, of whom 98 (74%) had an isolated CDH and 35 (26%) had a complex CDH (associated structural anomalies) on prenatal examination. ES was performed in 68 cases, and eight pathogenic or likely pathogenic variants were found, accounting for a 12% diagnostic yield (10% [5/50] in isolated cases and 17% [3/18] in complex CDH)., Conclusions: In 12% of fetuses and neonates with CDH and normal CNV analysis results, pathogenic or likely pathogenic variants were identified with ES. These data indicate that there is a substantial diagnostic yield when offering ES in prenatally detected CDH, both in complex and isolated cases., (© 2024 The Author(s). Prenatal Diagnosis published by John Wiley & Sons Ltd.)

Objectives: Determine the incremental diagnostic yield of prenatal exome sequencing (pES) over chromosome microarray (CMA) or G-banding karyotype in fetuses with central nervous system (CNS) abnormalities., Methods: Data were collected via electronic searches from January 2010 to April 2022 in MEDLINE, Cochrane, Web of Science and EMBASE. The NHS England prenatal exome cohort was also included. Incremental yield was calculated as a pooled value using a random-effects model., Results: Thirty studies were included (n = 1583 cases). The incremental yield with pES for any CNS anomaly was 32% [95%CI 27%-36%; I 2  = 72%]. Subgroup analysis revealed apparent incremental yields in; (a) isolated CNS anomalies; 27% [95%CI 19%-34%; I 2  = 74%]; (b) single CNS anomaly; 16% [95% CI 10%-23%; I 2  = 41%]; (c) more than one CNS anomaly; 31% [95% Cl 21%-40%; I 2  = 56%]; and (d) the anatomical subtype with the most optimal yield was Type 1 malformation of cortical development, related to abnormal cell proliferation or apoptosis, incorporating microcephalies, megalencephalies and dysplasia; 40% (22%-57%; I 2  = 68%). The commonest syndromes in isolated cases were Lissencephaly 3 and X-linked hydrocephalus., Conclusions: Prenatal exome sequencing provides a high incremental diagnostic yield in fetuses with CNS abnormalities with optimal yields in cases with multiple CNS anomalies, particularly those affecting the midline, posterior fossa and cortex., (© 2023 The Authors. Prenatal Diagnosis published by John Wiley & Sons Ltd.)

Purpose: This study presents a case of familial transmission of thyroxine-binding globulin (TBG) deficiency. The SERPINA7-gene which codes for TBG is located on the X-chromosome (Xq21-22). More than 45 mutations have been reported to cause TBG- deficiency from various countries, but none from India so far. Genetic analysis of SERPINA7 gene was carried out to determine the cause of low TBG levels in one family., Methods: DNA samples of the propositus and the family members were subjected to Polymerase Chain Reaction (PCR) followed by direct sequencing. Allele-specific PCR and Next-gen sequencing (NGS) were employed to confirm the site of the mutation. Thyroid function tests were estimated by Radioimmunoassay (RIA) and Immunoradiometric assay (IRMA) kits. X-chromosomal inactivation status was analyzed in the female members harboring the mutation., Results: A mutational screening in this family revealed a novel frame-shift mutation S353Q, 354fs3X in the exon 4 of the SERPINA7 gene which will be referred to as TBG-complete deficiency-India (TBG-CD-Ind). One out of four female family members harboring the mutation showed selective X-chromosomal inactivation. The affected family members were clinically euthyroid initially, showed changes in the thyroid function when tested after a long time span. However, the changes in the thyroid function in the affected family members had an autoimmune etiology., Conclusion: This study presents the first report of TBG-CD from India wherein a novel frameshift mutation referred to as TBG-CD-Ind (S353Q, 354fs3X) in the SERPINA7 gene was detected. No apparent association was identified between thyroid function and the TBG-mutation in the affected subjects. A detailed biochemical and genomic testing to determine the exact cause of discordant TFT in the patients would certainly aid in the unequivocal diagnosis of the thyroid function and for the precise individualized treatment., (© 2021. Italian Society of Endocrinology (SIE).)

Background: Bone morphogenetic protein 15 (BMP15) is expressed in oocytes and plays a crucial role in the reproduction of mono-ovulating species. In humans, BMP15 gene mutations lead to imperfect protein function and premature ovarian insufficiency. Here we investigated the BMP15 gene variants in a population of Iranian women with premature ovarian insufficiency. We conducted predictive bioinformatics analysis to further study the outcomes of BMP15 gene alterations., Methods: Twenty-four well-diagnosed premature ovarian insufficiency cases with normal karyotype participated in this study. The entire coding sequence and exon-intron junctions of the BMP15 gene were analyzed by direct sequencing. In-silico analysis was applied using various pipelines integrated into the Ensembl Variant Effect Predictor online tool. The clinical interpretation was performed based on the approved guidelines., Results: By gene screening of BMP15, we discovered p.N103K, p.A180T, and p.M184T heterozygous variants in 3 unrelated patients. The p.N103K and p.M184T were not annotated on gnomAD, 1000 Genome and/or dbSNP. These mutations were not identified in 800 Iranians whole-exome sequencing that is recorded on Iranom database. We identified the p.N103K variant in a patient with secondary amenorrhea at the age of 17, elevated FSH and atrophic ovaries. The p.M184T was detected in a sporadic case with atrophic ovaries and very high FSH who developed secondary amenorrhea at the age of 31., Conclusions: Here we newly identified p.N103K and p.M184T mutation in the BMP15 gene associated with idiopathic premature ovarian insufficiency. Both mutations have occurred in the prodomain region of protein. Despite prodomain cleavage through dimerization, it is actively involved in the mature protein function. Further studies elucidating the roles of prodomain would lead to a better understanding of the disease pathogenesis., (© 2022. The Author(s).)

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., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2022 The Authors. Published by Elsevier Inc. All rights reserved.)

Whole exome sequencing (WES) is a powerful tool for the cataloguing of population-specific genetic diseases. Within this proof-of-concept study we evaluated whether analysis of a small number of individual exomes is capable of identifying recurrent pathogenic alleles. We considered 106 exomes of subjects of Russian origin and revealed 13 genetic variants, which occurred more than twice and fulfilled the criteria for pathogenicity. All these alleles turned out to be indeed recurrent, as revealed by the analysis of 1045 healthy Russian donors. Eight of these variants (NAGA c.973G>A, ACADM c.985A>C, MPO c.2031-2A>C, SLC3A1 c.1400T>C, LRP2 c.6160G>A, BCHE c.293A>G, MPO c.752T>C, FCN3 c.349delC) are non-Russian-specific, as their high prevalence was previously demonstrated in other European populations. The remaining five disease-associated alleles appear to be characteristic for subjects of Russian origin and include CLCN1 c.2680C>T (myotonia congenita), DHCR7 c.453G>A (Smith-Lemli-Opitz syndrome), NUP93 c.1162C>T (steroid-resistant nephrotic syndrome, type 12), SLC26A2 c.1957T>A (multiple epiphyseal dysplasia) and EIF3F c.694T>G (mental retardation). These recessive disease conditions may be of particular relevance for the Russian Federation and other countries with a significant Slavic population., (Copyright © 2022 Elsevier Masson SAS. All rights reserved.)

Background: Gallbladder carcinoma (GBC) is a relatively rare but highly aggressive cancer with late clinical detection and a poor prognosis. However, the lack of models with features consistent with human gallbladder tumours has hindered progress in pathogenic mechanisms and therapies., Methods: We established organoid lines derived from human GBC as well as normal gallbladder and benign gallbladder adenoma (GBA) tissues. The histopathology signatures of organoid cultures were identified by H&E staining, immunohistochemistry and immunofluorescence. The genetic and transcriptional features of organoids were analysed by whole-exome sequencing and RNA sequencing. A set of compounds targeting the most active signalling pathways in GBCs were screened for their ability to suppress GBC organoids. The antitumour effects of candidate compounds, CUDC-101 and CUDC-907, were evaluated in vitro and in vivo., Results: The established organoids were cultured stably for more than 6 months and closely recapitulated the histopathology, genetic and transcriptional features, and intratumour heterogeneity of the primary tissues at the single-cell level. Notably, expression profiling analysis of the organoids revealed a set of genes that varied across the three subtypes and thus may participate in the malignant progression of gallbladder diseases. More importantly, we found that the dual PI3K/HDAC inhibitor CUDC-907 significantly restrained the growth of various GBC organoids with minimal toxicity to normal gallbladder organoids., Conclusions: Patient-derived organoids are potentially a useful platform to explore molecular pathogenesis of gallbladder tumours and discover personalized drugs., (© 2022 The Authors. Clinical and Translational Medicine published by John Wiley & Sons Australia, Ltd on behalf of Shanghai Institute of Clinical Bioinformatics.)

Purpose: Previous studies have shown an approximately two-fold elevation in the relative risk of urinary bladder cancer (UBC) among people with a family history that could not be entirely explained by shared environmental exposures, thus suggesting a genetic component in its predisposition. Multiple genome-wide association studies and recent gene panel sequencing studies identified several genetic loci that are associated with UBC risk; however, the list of UBC-associated variants and genes is incomplete., Materials and Methods: We exome sequenced eight patients from three multiplex UBC pedigrees and a group of 77 unrelated familial UBC cases matched to 241 cancer-free controls. In addition, we examined pathogenic germline variation in 444 candidate genes in 392 The Cancer Genome Atlas UBC cases., Results: In the pedigrees, segregating variants were family-specific although the identified genes clustered in common pathways, most notably DNA repair ( MLH1 and MSH2 ) and cellular metabolism ( IDH1 and ME1 ). In the familial UBC group, the proportion of pathogenic and likely pathogenic variants was significantly higher in cases compared with controls ( P = .003). Pathogenic and likely pathogenic variant load was also significantly increased in genes involved in cilia biogenesis ( P = .001). In addition, a pathogenic variant in CHEK2 (NM_007194.4:c.1100del; p.T367Mfs*15) was over-represented in cases (variant frequency = 2.6%; 95% CI, 0.71 to 6.52) compared with controls (variant frequency = 0.21%; 95% CI, 0.01 to 1.15), but was not statistically significant., Conclusion: These results point to a complex polygenic predisposition to UBC. Despite heterogeneity, the genes cluster in several biologically relevant pathways and processes, for example, DNA repair, cilia biogenesis, and cellular metabolism. Larger studies are required to determine the importance of CHEK2 in UBC etiology., Competing Interests: Talia Wegman-Ostrosky Travel, Accommodations, Expenses: Pfizer Bin Zhu Employment: GeneDx BioReference Douglas R. Stewart Employment: Genome Medical No other potential conflicts of interest were reported. Talia Wegman-Ostrosky Travel, Accommodations, Expenses: Pfizer Bin Zhu Employment: GeneDx BioReference Douglas R. Stewart Employment: Genome Medical No other potential conflicts of interest were reported., (© 2021 by American Society of Clinical Oncology.)

Background & Aims: Our previous genomic whole-exome sequencing (WES) data identified the key ErbB pathway mutations that play an essential role in regulating the malignancy of gallbladder cancer (GBC). Herein, we tested the hypothesis that individual cellular components of the tumor microenvironment (TME) in GBC function differentially to participate in ErbB pathway mutation-dependent tumor progression., Methods: We engaged single-cell RNA-sequencing to reveal transcriptomic heterogeneity and intercellular crosstalk from 13 human GBCs and adjacent normal tissues. In addition, we performed WES analysis to reveal the genomic variations related to tumor malignancy. A variety of bulk RNA-sequencing, immunohistochemical staining, immunofluorescence staining and functional experiments were employed to study the difference between tissues with or without ErbB pathway mutations., Results: We identified 16 cell types from a total of 114,927 cells, in which epithelial cells, M2 macrophages, and regulatory T cells were predominant in tumors with ErbB pathway mutations. Furthermore, epithelial cell subtype 1, 2 and 3 were mainly found in adenocarcinoma and subtype 4 was present in adenosquamous carcinoma. The tumors with ErbB pathway mutations harbored larger populations of epithelial cell subtype 1 and 2, and expressed higher levels of secreted midkine (MDK) than tumors without ErbB pathway mutations. Increased MDK resulted in an interaction with its receptor LRP1, which is expressed by tumor-infiltrating macrophages, and promoted immunosuppressive macrophage differentiation. Moreover, the crosstalk between macrophage-secreted CXCL10 and its receptor CXCR3 on regulatory T cells was induced in GBC with ErbB pathway mutations. Elevated MDK was correlated with poor overall survival in patients with GBC., Conclusions: This study has provided valuable insights into transcriptomic heterogeneity and the global cellular network in the TME, which coordinately functions to promote the progression of GBC with ErbB pathway mutations; thus, unveiling novel cellular and molecular targets for cancer therapy., Lay Summary: We employed single-cell RNA-sequencing and functional assays to uncover the transcriptomic heterogeneity and intercellular crosstalk present in gallbladder cancer. We found that ErbB pathway mutations reduced anti-cancer immunity and led to cancer development. ErbB pathway mutations resulted in immunosuppressive macrophage differentiation and regulatory T cell activation, explaining the reduced anti-cancer immunity and worse overall survival observed in patients with these mutations., Competing Interests: Conflict of interest The authors declare no conflicts of interest that pertain to this work. Please refer to the accompanying ICMJE disclosure forms for further details., (Copyright © 2021 The Author(s). Published by Elsevier B.V. All rights reserved.)

Background: Uniparental disomy (UPD) is the inheritance of two homologous chromosomes from the same parent. UPD may result in clinical phenotypes when occurring on chromosomes with specific imprinting pattern, when leading to homozygosity of a deleterious recessive allele inherited from one carrier parent, or when associated with a mosaic aneuploidy. Due to the importance of UPD in genetic disease etiology, UPD analysis has started to be implemented in the context of exome sequencing (ES) or genome sequencing., Methods: We developed an in-house algorithm TRIPS (Trio Parentage/UPD Studies) to identify UPD events in trio ES cases. This method identifies regions with uniparental inheritance by utilizing the trio genotyping data obtained from the concurrent SNP array to delineate the parental origin of the SNPs in the proband., Results: We identified 16 UPD events from 2675 ES trios. Among those, four events led to imprinting disorders, seven unmasked a pathogenic/likely pathogenic variant in a recessive disease gene, and two were consistent with a mosaic genome wide paternal UPD pattern. Twelve of these UPD events directly contributed to the molecular diagnosis of the patients., Conclusion: Our study demonstrated the contribution of UPD to the molecular diagnosis in one clinical ES cohort, thus UPD analysis should be incorporated into routine clinical ES interpretation., (© 2021 The Authors. Molecular Genetics & Genomic Medicine published by Wiley Periodicals LLC.)

Objectives: Pancreatic cancer (PC) is one of the most lethal cancers. Currently, the most commonly used PC cell lines were all derived from White individuals, and few models can be used to study in Asian populations. In China, the incidence of PC is increasing every year, highlighting the need to develop new Chinese PC cell lines for cancer research., Methods: A total of 203 patients diagnosed with PC, enrolled from January 2012 to December 2017 in our hospital (Department of Hepatobiliary and Pancreatic Surgery), had received surgery. Primary PC patient tumor samples were harvested sterilely from surgery, minced, and digested in collagenase. Cells were cultured in plates precoated with l-polylysine., Results: Two primary PC cell lines, Si-Liang 187 and Si-Liang 188, were established, both of which grew as adherent monolayers and demonstrated distinguished phenotypes. Si-Liang 187 shows the typical mesenchymal pattern, whereas Si-Liang 188 shows epithelial phenotype. Moreover, they carry different genetic backgrounds according to whole exome sequencing., Conclusions: Two new Chinese PC cell lines were established, both of which were characterized and confirmed with high tumorigenicity. They may serve as useful tools for pathogenesis research when evaluating new treatment strategies in Asian patients., Competing Interests: The authors declare no conflict of interest., (Copyright © 2022 Wolters Kluwer Health, Inc. All rights reserved.)

Objective: To investigate the genetic etiology of skeletal dysplasia in highly selected fetuses during the first and second trimesters using deep phenotyping and exome sequencing (ES)., Method: Fetuses with short femurs were identified using the established prenatal diagnostic approach. A multidisciplinary team reviewed fetal phenotypic information (prenatal ultrasound findings, fetal postmortem, and radiographs) in a cohort of highly selected fetuses with skeletal dysplasia during the first and second trimesters. The affected families underwent multiplatform genetic tests., Results: Of the 27 affected fetuses, 21 (77.8%) had pathogenic or potential pathogenic variations in the following genes: COL1A1, FGFR3, COL2A1, COL1A2, FLNB, DYNC2LI1, and TRIP11. Two fetuses had compound heterozygous mutations in DYNC2LI1 and TRIP11, respectively, and the other 19 carried de novo autosomal dominant variants. Novel variants were identified in COL1A1, COL2A1, COL1A2, DYNC2LI1, and TRIP11 in 11 fetuses. We also included the first description of the phenotype of odontochondrodysplasia in a prenatal setting., Conclusions: ES or panel sequencing offers a high diagnostic yield for fetal skeletal dysplasia during the first and second trimesters. Comprehensive and complete phenotypic information is indispensable for genetic analysis and the expansion of genotype-phenotype correlations in fetal skeletal abnormalities., (© 2021 John Wiley & Sons Ltd.)

Introduction: Nonimmune hydrops fetalis (NIHF) has varied etiology. We assessed the etiological spectrum and evaluated the utility of fetal whole exome sequencing (fWES) for the diagnosis of NIHF., Methods: In this prospective cohort study, we evaluated antenatally diagnosed fetuses with NIHF between July 2018 and December 2019 according to the routine diagnostic algorithm. Fetuses that remained undiagnosed after routine NIHF workup were subjected to fetal chromosomal microarray and/or WES. Pregnancies were followed up for clinical outcomes., Results: Of the 45 fetuses, consanguinity and recurrent hydrops fetalis were observed in 13.3% (6/45) and 28.8% (13/45), respectively. Overall, an etiological diagnosis was possible in 75.5% (34/45) of fetuses, while the cause remained unknown in 24.4% (11/45). A genetic etiology was identified in 46.6% (21/45): aneuploidy and monogenic disorders in 28.8% (13/45) and 17.8% (8/45), respectively. fWES on 19 fetuses detected disease-causing variants in 42.1% (8/19). Nine novel variants were detected in RAPSN, ASCC1, NEB, PKD1L1, GUSB, and PIEZO1. Only 8.8% (4/45) of the cohort survived without morbidity., Conclusions: This study describes the etiological spectrum and the disease-causing variants in an Indian cohort of hydropic fetuses., (© 2021 John Wiley & Sons Ltd.)

Objective: To determine the incremental yield of exome sequencing (ES) over chromosomal microarray analysis (CMA) or karyotyping in prenatally diagnosed non-immune hydrops fetalis (NIHF)., Methods: A prospective cohort study (comprising an extended group of the Prenatal Assessment of Genomes and Exomes (PAGE) study) was performed which included 28 cases of prenatally diagnosed NIHF undergoing trio ES following negative CMA or karyotyping. These cases were combined with data from a systematic review of the literature. MEDLINE, EMBASE, CINAHL and ClinicalTrials.gov databases were searched electronically (January 2000 to October 2020) for studies reporting on the incremental yield of ES over CMA or karyotyping in fetuses with prenatally detected NIHF. Inclusion criteria for the systematic review were: (i) at least two cases of NIHF undergoing sequencing; (ii) testing initiated based on prenatal ultrasound-based phenotype; and (iii) negative CMA or karyotyping result. The incremental diagnostic yield of ES was assessed in: (i) all cases of NIHF; (ii) isolated NIHF; (iii) NIHF associated with an additional fetal structural anomaly; and (iv) NIHF according to severity (i.e. two vs three or more cavities affected)., Results: In the extended PAGE study cohort, the additional diagnostic yield of ES over CMA or karyotyping was 25.0% (7/28) in all NIHF cases, 21.4% (3/14) in those with isolated NIHF and 28.6% (4/14) in those with non-isolated NIHF. In the meta-analysis, the pooled incremental yield based on 21 studies (306 cases) was 29% (95% CI, 24-34%; P < 0.00001; I 2  = 0%) in all NIHF, 21% (95% CI, 13-30%; P < 0.00001; I 2  = 0%) in isolated NIHF and 39% (95% CI, 30-49%; P < 0.00001; I 2  = 1%) in NIHF associated with an additional fetal structural anomaly. In the latter group, congenital limb contractures were the most prevalent additional structural anomaly associated with a causative pathogenic variant, occurring in 17.3% (19/110) of cases. The incremental yield did not differ significantly according to hydrops severity. The most common genetic disorders identified were RASopathies, occurring in 30.3% (27/89) of cases with a causative pathogenic variant, most frequently due to a PTPN11 variant (44.4%; 12/27). The predominant inheritance pattern in causative pathogenic variants was autosomal dominant in monoallelic disease genes (57.3%; 51/89), with most being de novo (86.3%; 44/51)., Conclusions: Use of prenatal next-generation sequencing in both isolated and non-isolated NIHF should be considered in the development of clinical pathways. Given the wide range of potential syndromic diagnoses and heterogeneity in the prenatal phenotype of NIHF, exome or whole-genome sequencing may prove to be a more appropriate testing approach than a targeted gene panel testing strategy. © 2021 The Authors. Ultrasound in Obstetrics & Gynecology published by John Wiley & Sons Ltd on behalf of International Society of Ultrasound in Obstetrics and Gynecology., (© 2021 The Authors. Ultrasound in Obstetrics & Gynecology published by John Wiley & Sons Ltd on behalf of International Society of Ultrasound in Obstetrics and Gynecology.)

Objective: To evaluate the utility of clinical exome sequencing (ES)-based carrier screening in Chinese consanguineous couples., Methods: Consanguineous couples were screened for autosomal recessive (AR) disorders using the clinical ES of 5000 genes associated with human diseases., Results: We recruited 14 couples who elected to have sequencing. One couple was related as first cousins and 13 as second cousins. Both partners carrying the same pathogenic variant were detected in four couples. One couple was found in which one partner carried a splice variant, and the other had a missence variant of the same gene. These five couples were identified as being at risk of having a child affected by an AR disorder., Conclusion: Our study demonstrates that ES-based preconception screening yields a clinical value for Chinese consanguineous couples. It enables to detect at-risk couples for rare AR diseases., (© 2021 John Wiley & Sons Ltd.)

Exome sequencing (ES) is an important diagnostic tool for individuals with neurodevelopmental disorders (NDD) and/or multiple congenital anomalies (MCA). However, the cost of ES limits the test's accessibility for many patients. We evaluated the yield of publicly funded clinical ES, performed at a tertiary center in Israel, over a 3-year period (2018-2020). Probands presented with (1) moderate-to-profound global developmental delay (GDD)/intellectual disability (ID); or (2) mild GDD/ID with epilepsy or congenital anomaly; and/or (3) MCA. Subjects with normal chromosomal microarray analysis who met inclusion criteria were included, totaling 280 consecutive cases. Trio ES (proband and parents) was the default option. In 252 cases (90.0%), indication of NDD was noted. Most probands were males (62.9%), and their mean age at ES submission was 9.3 years (range 1 month to 51 years). Molecular diagnosis was reached in 109 probands (38.9%), mainly due to de novo variants (91/109, 83.5%). Disease-causing variants were identified in 92 genes, 15 of which were implicated in more than a single case. Male sex, families with multiple-affected members and premature birth were significantly associated with lower ES yield (p < 0.05). Other factors, including MCA and coexistence of epilepsy, autism spectrum disorder, microcephaly or abnormal brain magnetic resonance imaging findings, were not associated with the yield. To conclude, our findings support the utility of clinical ES in a real-world setting, as part of a publicly funded genetic workup for individuals with GDD/ID and/or MCA., (© 2021. The Author(s).)

Background: Prenatal whole-exome sequencing (WES) is becoming increasingly used when karyotype and microarray tests are not diagnostic of fetal malformations. Although the value of WES clearly emerges in terms of higher diagnostic rates, the limitations of prenatal phenotyping together with the counseling challenges for variants of uncertain significance and incidental results suggest that the routine application of prenatal WES is not yet easy., Methods: Structurally abnormal fetuses with a mean gestational age of 24 weeks (range 13-38 weeks) were recruited from the Chong Qing Health Center for Women and Children. We performed a retrospective WES investigation in 85 fetuses, using DNA from amniotic fluid (66 samples, 77.6%), umbilical cord blood (10 samples, 11.8%), and fetal tissues (9 samples, 10.6%). Parental DNA was extracted from peripheral blood., Results: Molecular diagnosis was obtained in 16 of the 85 fetuses (18.8%). According to the variant segregation mode and family history, 7 fetuses (43.75%) were affected by an autosomal dominant condition (6 variants were de novo and 1 variant was inherited from an unknowingly affected father), 7 fetuses (43.75%) had an autosomal recessive syndrome always associated with compound heterozygosity, and 2 fetuses (12.5%) had an X-linked condition (one mother was a carrier). In addition, the highest diagnostic rate was observed in fetuses with multisystem abnormalities (38.9%, 7/18). A variant of uncertain significance was detected in 16 samples (18.8%, 16/85)., Conclusion: Our study confirms that prenatal WES is an efficient tool for studying fetal abnormalities, although further improvements are needed to establish stronger fetal genotype-phenotype correlations., (Copyright © 2021. Published by Elsevier Masson SAS.)

Single-nucleotide polymorphisms (SNPs) are the most common genetic variations for various complex human diseases, including cancers. Genome-wide association studies (GWAS) have identified numerous SNPs that increase cancer risks, such as breast cancer, colorectal cancer, and leukemia. These SNPs were cataloged for scientific use. However, GWAS are often conducted on certain populations in which the Orang Asli and Malays were not included. Therefore, we have developed a bioinformatic pipeline to mine the whole-genome sequence databases of the Orang Asli and Malays to determine the presence of pathogenic SNPs that might increase the risks of cancers among them. Five different in silico tools, SIFT, PROVEAN, Poly-Phen-2, Condel, and PANTHER, were used to predict and assess the functional impacts of the SNPs. Out of the 80 cancer-related nsSNPs from the GWAS dataset, 52 nsSNPs were found among the Orang Asli and Malays. They were further analyzed using the bioinformatic pipeline to identify the pathogenic variants. Three nsSNPs; rs1126809 (TYR), rs10936600 (LRRC34), and rs757978 (FARP2), were found as the most damaging cancer pathogenic variants. These mutations alter the protein interface and change the allosteric sites of the respective proteins. As TYR, LRRC34, and FARP2 genes play important roles in numerous cellular processes such as cell proliferation, differentiation, growth, and cell survival; therefore, any impairment on the protein function could be involved in the development of cancer. rs1126809, rs10936600, and rs757978 are the important pathogenic variants that increase the risks of cancers among the Orang Asli and Malays. The roles and impacts of these variants in cancers will require further investigations using in vitro cancer models., (© 2021. The Author(s).)

In childhood cancer, the frequency of cancer-associated germline variants and their inheritance patterns are not thoroughly investigated. Moreover, the identification of children carrying a genetic predisposition by clinical means remains challenging. In this single-center study, we performed trio whole-exome sequencing and comprehensive clinical evaluation of a prospectively enrolled cohort of 160 children with cancer and their parents. We identified in 11/160 patients a pathogenic germline variant predisposing to cancer and a further eleven patients carried a prioritized VUS with a strong association to the cancerogenesis of the patient. Through clinical screening, 51 patients (31.3%) were identified as suspicious for an underlying cancer predisposition syndrome (CPS), but only in ten of those patients a pathogenic variant could be identified. In contrast, one patient with a classical CPS and ten patients with prioritized VUS were classified as unremarkable in the clinical work-up. Taken together, a monogenetic causative variant was detected in 13.8% of our patients using WES. Nevertheless, the still unclarified clinical suspicious cases emphasize the need to consider other genetic mechanisms including new target genes, structural variants, or polygenic interactions not previously associated with cancer predisposition., (© 2021. The Author(s).)

Exome association studies to date have generally been underpowered to systematically evaluate the phenotypic impact of very rare coding variants. We leveraged extensive haplotype sharing between 49,960 exome-sequenced UK Biobank participants and the remainder of the cohort (total n ≈ 500,000) to impute exome-wide variants with accuracy R 2  > 0.5 down to minor allele frequency (MAF) ~0.00005. Association and fine-mapping analyses of 54 quantitative traits identified 1,189 significant associations (P < 5 × 10 -8 ) involving 675 distinct rare protein-altering variants (MAF < 0.01) that passed stringent filters for likely causality. Across all traits, 49% of associations (578/1,189) occurred in genes with two or more hits; follow-up analyses of these genes identified allelic series containing up to 45 distinct 'likely-causal' variants. Our results demonstrate the utility of within-cohort imputation in population-scale genome-wide association studies, provide a catalog of likely-causal, large-effect coding variant associations and foreshadow the insights that will be revealed as genetic biobank studies continue to grow., (© 2021. The Author(s), under exclusive licence to Springer Nature America, Inc.)

Background: In patients with advanced melanoma (MM), genomic profiling may guide treatment decisions in the frontline setting and beyond as specific tumor mutations can be treated with targeted therapy (TT). The range of panel sizes used to identify targetable mutations (TM) can range from a few dozen to whole exome sequencing (WES)., Aim: We investigated the impact of panel size and mutation status on first-line treatment selection and outcomes in MM., Methods and Results: We analyzed data for 1109 MM patients from three cohorts: 169 patients at NYULH and profiled with the 50 gene Ion Torrent panel (IT), 195 patients at MSKCC, profiled with the 400-gene MSK-IMPACT panel (MSK-I) and 745 patients at seven different sites profiled with WES. Data for cohorts 2 and 3 were extrapolated from the publicly available cBioPortal. Treatment information was available for 100%, 25%, and 0% of patients in cohort 1, 2, and 3, respectively. BRAF and NRAS were among the top five most commonly mutated genes in the IT and MSK-I, whereas for WES only BRAF was a top five mutation. There was no significant difference in OS for BRAF MUT patients treated with immune checkpoint inhibitors (ICI) vs TT in cohort 1 (P = .19), nor for BRAF MUT patients from cohort 1 treated with ICI vs those from cohort 2 treated with TT (P = .762)., Conclusion: Public datasets provide population-level data; however, the heterogeneity of reported clinical information limits their value and calls for data standardization. Without evidence of clear clinical benefit of a larger panel size, there is a rationale for adopting smaller, more cost effective panels in MM., (© 2021 The Authors. Cancer Reports published by Wiley Periodicals LLC.)

Objectives: To study the utility of clinical exome sequencing (CES) using next generation sequencing (NGS) in evaluating neonates with suspected genetic conditions., Methods: This is an observational study conducted in a tertiary care neonatal unit. We included neonates with suspected genetic conditions, for whom CES were done either by direct sampling or from stored DNA. Data was collected from the Sri Ramachandra centre of excellence in perinatal health (SCOPE) case records of 2016-2019. Yield of CES, percentage of pathogenic, non-pathogenic and variant of uncertain significance (VUS) and associated disorders were studied., Results: CES was done in 36 neonates. Variants were detected in 78% (28/36). However, significant variants with clinical correlation were present in 20 (56%) babies. Test was carried out from the stored sample in 10 (28%) babies. Mean turn-around time was 39 ± 7 days. Specialist was involved in 1 and treatment changes were done in 5 neonates. Five out of 8 VUS were clinically correlating. Inborn errors of metabolism were the commonest (60%). Two VUS were ascertained as likely pathogenic after parental segregation analysis., Conclusion: CES has a definite role in evaluation of suspected genetic conditions for diagnosis and prognostication. It also helps scientific society to build in additional evidence so that the "VUS" could be asserted as "likely pathogenic" . Our experience reiterates the importance of storing and archiving DNA of the affected child., (Copyright © 2021 Elsevier Masson SAS. All rights reserved.)

Objectives/hypothesis: To assess the relative cost-effectiveness of exome sequencing for isolated congenital deafness compared with standard care., Study Design: Incremental cost-effectiveness and cost-benefit analyses were undertaken from the perspective of the Australian healthcare system using an 18-year time horizon., Methods: A decision tree was used to model the costs and outcomes associated with exome sequencing and standard care for infants presenting with isolated congenital deafness., Results: Exome sequencing resulted in an incremental cost of AU$1,000 per child and an additional 30 diagnoses per 100 children tested. The incremental cost-effectiveness ratio was AU$3,333 per additional diagnosis. The mean societal willingness to pay for exome sequencing was estimated at AU$4,600 per child tested relative to standard care, resulting in a positive net benefit of AU$3,600. Deterministic and probabilistic sensitivity analyses confirmed the cost-effectiveness of exome sequencing., Conclusions: Our findings demonstrate the cost-effectiveness of exome sequencing in congenital hearing loss, through increased diagnostic rate and consequent improved process of care by reducing or ceasing diagnostic investigation or facilitating targeted further investigation. We recommend equitable funding for exome sequencing in infants presenting with isolated congenital hearing loss., Level of Evidence: N/A. Laryngoscope, 131:E2371-E2377, 2021., (© 2020 The American Laryngological, Rhinological and Otological Society, Inc.)

Background: There is limited real-world information on use of tumor mutational burden (TMB) testing and characteristics of patients receiving it. Materials & methods: Patients ≥18 years old and diagnosed with advanced solid tumors between 1 January 2015 and 31 January 2019 with TMB testing (TMB cohort) and without it (non-TMB) were included in this retrospective, observational study. Results: The TMB cohort (n = 202) was younger than non-TMB (n = 212) (mean age: 62.1 vs 65.6 at diagnosis; p = 0.005). There were more Black patients in the TMB cohort (21.3 vs 11.8% in non-TMB; p = 0.004). Clinical characteristics were comparable between the two cohorts; however, systemic anticancer treatment was higher among TMB cohort (91.6 vs 77.8% in non-TMB). Conclusion: Notable differences were observed between patients receiving TMB test and those not receiving it.

Objective: Laboratories performing prenatal exome sequencing (ES) frequently limit analysis to predetermined gene lists. We used a diagnostic postnatal ES cohort to assess how many of the genes diagnosed are not included in a number of select fixed lists used for prenatal diagnosis., Methods: Of 601 postnatal ES tests, pathogenic variants related to neurodevelopmental disorders were detected in 138 probands. We evaluated if causative genes were present in the following: (1) Developmental Disorders Genotype-Phenotype database list, (2) a commercial laboratory list for prenatal ES, (3) the PanelApp fetal anomalies panel, and (4) a published list used for prenatal diagnosis by ES (Prenatal Assessment of Genomes and Exomes study)., Results: The percentages of cases where the diagnosed gene was not included in the selected four lists were; 11.6%, 17.24%, 23.2%, and 10.9%, respectively. In 13/138 (9.4%) cases, the causative gene was not included in any of the lists; in 4/13 (∼30%) cases noninclusion was explained by a relatively recent discovery of gene-phenotype association., Conclusions: A significant number of genes related to neurocognitive phenotypes are not included in some of the lists used for prenatal ES data interpretation. These are not only genes related to recently discovered disorders, but also genes with well-established gene-phenotype., (© 2021 John Wiley & Sons Ltd.)

Dozens of genes are associated with idiopathic hypogonadotropic hypogonadism (IHH) and an oligogenic etiology has been suggested. However, the associated genes may account for only approximately 50% cases. In addition, a genomic systematic pedigree analysis is still lacking. Here, we conducted whole exome sequencing (WES) on 18 unrelated men affected by IHH and their corresponding parents. Notably, one reported and 10 novel variants in eight known IHH causative genes (AXL, CCDC141, CHD7, DMXL2, FGFR1, PNPLA6, POLR3A, and PROKR2), nine variants in nine recently reported candidate genes (DCAF17, DCC, EGF, IGSF10, NOTCH1, PDE3A, RELN, SLIT2, and TRAPPC9), and four variants in four novel candidate genes for IHH (CCDC88C, CDON, GADL1, and SPRED3) were identified in 77.8% (14/18) of IHH cases. Among them, eight (8/18, 44.4%) cases carried more than one variant in IHH-related genes, supporting the oligogenic model. Interestingly, we found that those variants tended to be maternally inherited (maternal with n = 17 vs paternal with n = 7; P = 0.028). Our further retrospective investigation of published reports replicated the maternal bias (maternal with n = 46 vs paternal with n = 28; P = 0.024). Our study extended a variant spectrum for IHH and provided the first evidence that women are probably more tolerant to variants of IHH-related genes than men., Competing Interests: None

Objectives: To conduct qualitative interviews with healthcare providers working in different countries to understand their experiences of dealing with uncertain results from prenatal chromosome microarray analysis (CMA) and exome sequencing (ES)., Methods: Semi-structured interviews with 31 healthcare providers who report or return prenatal CMA and/or ES results (clinicians, genetic counsellors and clinical scientists) in six countries with differing healthcare systems; Australia (4), Denmark (5), Netherlands (6), Singapore (4), Sweden (6) and United Kingdom (6). The topic guide explored the main sources of uncertainty and their management., Results: There was variation in reporting practices both between and across countries for variants of uncertain significance, however, there was broad agreement on reporting practices for incidental findings. There was also variation in who decides what results are reported (clinical scientists or clinicians). Technical limitations and lack of knowledge (to classify variants and of prenatal phenotypes) were significant challenges, as were turnaround times and lack of guidelines., Conclusion: Health professionals around the globe are dealing with similar sources of uncertainty, but managing them in different ways, Continued dialogue with international colleagues on ways of managing uncertain results is important to compare and contrast the benefits and limitations of the different approaches., (© 2021 The Authors. Prenatal Diagnosis published by John Wiley & Sons Ltd.)

The success of genome-wide association studies (GWAS) in identifying common, low-penetrance variant-cancer associations for the past decade is undisputed. However, discovering additional high-penetrance cancer mutations in unknown cancer predisposing genes requires detection of variant-cancer association of ultra-rare coding variants. Consequently, large-scale next-generation sequence data with associated phenotype information are needed. Here, we used genotype data on 166,281 Icelanders, of which, 49,708 were whole-genome sequenced and 408,595 individuals from the UK Biobank, of which, 41,147 were whole-exome sequenced, to test for association between loss-of-function burden in autosomal genes and basal cell carcinoma (BCC), the most common cancer in Caucasians. A total of 25,205 BCC cases and 683,058 controls were tested. Rare germline loss-of-function variants in PTPN14 conferred substantial risks of BCC (OR, 8.0; P = 1.9 × 10 -12 ), with a quarter of carriers getting BCC before age 70 and over half in their lifetime. Furthermore, common variants at the PTPN14 locus were associated with BCC, suggesting PTPN14 as a new, high-impact BCC predisposition gene. A follow-up investigation of 24 cancers and three benign tumor types showed that PTPN14 loss-of-function variants are associated with high risk of cervical cancer (OR, 12.7, P = 1.6 × 10 -4 ) and low age at diagnosis. Our findings, using power-increasing methods with high-quality rare variant genotypes, highlight future prospects for new discoveries on carcinogenesis. SIGNIFICANCE: This study identifies the tumor-suppressor gene PTPN14 as a high-impact BCC predisposition gene and indicates that inactivation of PTPN14 by germline sequence variants may also lead to increased risk of cervical cancer., (©2021 American Association for Cancer Research.)

Consanguinity, commonplace in many regions around the globe, is associated with an increased risk of autosomal recessive (AR) genetic disorders. Consequently, consanguineous couples undergoing preimplantation genetic diagnosis (PGD) for one Mendelian disorder may be at increased risk for a child with a second, unrelated AR genetic disorder. We examined the yield of exome analysis for carrier screening of additional AR disorders, beyond the primary diagnosis, amongst consanguineous vs. non-consanguineous populations. Parental samples from trio exomes of 102 consanguineous families and 105 non-consanguineous controls were evaluated for shared carrier status, after disregarding the primary molecular diagnosis. Results were sub-classified according to disease severity. Secondary shared carrier status for pathogenic and likely pathogenic variants leading to AR disorders of moderate to profound severity was identified in 10/102 (9.8%) consanguineous couples, as compared to 1/105 (0.95%) non-consanguineous couples (χ 2  = 8.0565, p value < 0.005). Higher inbreeding coefficient values, calculated from individual exomes, correlated with secondary shared carrier status for diseases of moderate to profound severity (r = 0.17, p value < 0.0125). Our results indicate that consanguineous couples undergoing PGD are at increased risk for a second genetic disease of moderate to profound severity. This study represents an underestimate of the rate of secondary shared carrier status due to inability to detect deep intronic variants, no assessment of copy number variants, and false negative results stemming from stringent variant interpretation. False positive results may result from inaccuracies in public databases. Additional studies in consanguineous populations will determine whether exome-based carrier screening should be recommended to all couples undergoing PGD.

Rare variants are thought to play an important role in the etiology of complex diseases and may explain a significant fraction of the missing heritability in genetic disease studies. Next-generation sequencing facilitates the association of rare variants in coding or regulatory regions with complex diseases in large cohorts at genome-wide scale. However, rare variant association studies (RVAS) still lack power when cohorts are small to medium-sized and if genetic variation explains a small fraction of phenotypic variance. Here we present a novel Bayesian rare variant Association Test using Integrated Nested Laplace Approximation (BATI). Unlike existing RVAS tests, BATI allows integration of individual or variant-specific features as covariates, while efficiently performing inference based on full model estimation. We demonstrate that BATI outperforms established RVAS methods on realistic, semi-synthetic whole-exome sequencing cohorts, especially when using meaningful biological context, such as functional annotation. We show that BATI achieves power above 70% in scenarios in which competing tests fail to identify risk genes, e.g. when risk variants in sum explain less than 0.5% of phenotypic variance. We have integrated BATI, together with five existing RVAS tests in the 'Rare Variant Genome Wide Association Study' (rvGWAS) framework for data analyzed by whole-exome or whole genome sequencing. rvGWAS supports rare variant association for genes or any other biological unit such as promoters, while allowing the analysis of essential functionalities like quality control or filtering. Applying rvGWAS to a Chronic Lymphocytic Leukemia study we identified eight candidate predisposition genes, including EHMT2 and COPS7A., Competing Interests: The authors have declared that no competing interests exist.

Despite clear technical superiority of genome sequencing (GS) over other diagnostic methods such as exome sequencing (ES), few studies are available regarding the advantages of its clinical application. We analyzed 1007 consecutive index cases for whom GS was performed in a diagnostic setting over a 2-year period. We reported pathogenic and likely pathogenic (P/LP) variants that explain the patients' phenotype in 212 of the 1007 cases (21.1%). In 245 additional cases (24.3%), a variant of unknown significance (VUS) related to the phenotype was reported. We especially investigated patients which had had ES with no genetic diagnosis (n = 358). For this group, GS diagnostic yield was 14.5% (52 patients with P/LP out of 358). GS should be especially indicated for ES-negative cases since up to 29.6% of them  could benefit from GS testing (14.5% with P/LP, n = 52 and 15.1% with VUS, n = 54). Genetic diagnoses in most of the ES-negative/GS-positive cases were determined by technical superiority of GS, i.e., access to noncoding regions and more uniform coverage. Importantly, we reported 79 noncoding variants, of which, 41 variants were classified as P/LP. Interpretation of noncoding variants remains challenging, and in many cases, complementary methods based on direct enzyme assessment, biomarker testing and RNA analysis are needed for variant classification and diagnosis. We present the largest cohort of patients with GS performed in a clinical setting to date. The results of this study should direct the decision for GS as standard second-line, or even first-line stand-alone test.

The main purpose of pathway or gene set analysis methods is to provide mechanistic insight into the large amount of data produced in high-throughput studies. These tools were developed for gene expression analyses, but they have been rapidly adopted by other high-throughput techniques, becoming one of the foremost tools of omics research.Currently, according to different biological questions and data, we can choose among a vast plethora of methods and databases. Here we use two published examples of RNAseq datasets to approach multiple analyses of gene sets, networks and pathways using freely available and frequently updated software. Finally, we conclude this chapter by presenting a survival pathway analysis of a multiomics dataset. During this overview of different methods, we focus on visualization, which is a fundamental but challenging step in this computational field.

A complete RNA-Seq analysis involves the use of several different tools, with substantial software and computational requirements. The Galaxy platform simplifies the execution of such bioinformatics analyses by embedding the needed tools in its web interface, while also providing reproducibility. Here, we describe how to perform a reference-based RNA-Seq analysis using Galaxy, from data upload to visualization and functional enrichment analysis of differentially expressed genes.

Genetic research often utilizes or generates information that is potentially sensitive to individuals, families, or communities. For these reasons, genetic research may warrant additional scrutiny from investigators and governmental regulators, compared to other types of biomedical research. The informed consent process should address the range of social and psychological issues that may arise in genetic research. This article addresses a number of these issues, including recruitment of participants, disclosure of results, psychological impact of results, insurance and employment discrimination, community engagement, consent for tissue banking, and intellectual property issues. Points of consideration are offered to assist in the development of protocols and consent processes in light of contemporary debates on a number of these issues. © 2020 Wiley Periodicals LLC., (© 2020 Wiley Periodicals LLC.)

Small non-coding RNAs (ncRNAs) are short non-coding sequences involved in gene regulation in many biological processes and diseases. The lack of a complete comprehension of their biological functionality, especially in a genome-wide scenario, has demanded new computational approaches to annotate their roles. It is widely known that secondary structure is determinant to know RNA function and machine learning based approaches have been successfully proven to predict RNA function from secondary structure information. Here we show that RNA function can be predicted with good accuracy from a lightweight representation of sequence information without the necessity of computing secondary structure features which is computationally expensive. This finding appears to go against the dogma of secondary structure being a key determinant of function in RNA. Compared to recent secondary structure based methods, the proposed solution is more robust to sequence boundary noise and reduces drastically the computational cost allowing for large data volume annotations. Scripts and datasets to reproduce the results of experiments proposed in this study are available at: https://github.com/bioinformatics-sannio/ncrna-deep., Competing Interests: The authors have declared that no competing interests exist.

Background: Congenital muscular dystrophy type 1A (MDC1A), also termed merosin-deficient congenital muscular dystrophy (CMD), is a severe form of CMD caused by mutations in the laminin α2 gene (LAMA2). Of the more than 300 likely pathogenic variants found in the Leiden Open Variant Database, the majority are truncating mutations leading to complete LAMA2 loss of function, but multiple copy number variants (CNVs) have also been reported with variable frequency., Methods: We collected a cohort of individuals diagnosed with likely MDC1A and sought to identify both single nucleotide variants and small and larger CNVs via exome sequencing by extending the analysis of sequencing data to detect splicing changes and CNVs., Results: Standard exome analysis identified multiple novel LAMA2 variants in our cohort, but only four cases carried biallelic variants. Since likely truncating LAMA2 variants are often found in heterozygosity without a second allele, we performed additional splicing and CNV analysis on exome data and identified one splice change outside of the canonical sequences and three CNVs, in the remaining four cases., Conclusions: Our findings support the expectation that a portion of MDC1A cases may be caused by at least one CNV allele and show how these changes can be effectively identified by additional analysis of existing exome data., (© 2020 The Authors. Molecular Genetics & Genomic Medicine published by Wiley Periodicals LLC.)

Background: Exome sequencing has recently become more readily available, and more information about incidental findings has been disclosed. However, data from East Asia are scarce. We studied the application of exome sequencing to the identification of pathogenic/likely pathogenic variants in the ACMG 59 gene list and the frequency of these variants in the Taiwanese population., Methods: This study screened 161 Taiwanese exomes for variants from the ACMG 59 gene list. The identified variants were reviewed based on information from different databases and the available literature and classified according to the ACMG standard guidelines., Results: We identified seven pathogenic/likely pathogenic variants in eight individuals, with five participants with autosomal recessive variants in one allele and three participants with autosomal dominant variants. Approximately 1.86% (3/161) of the Taiwanese individuals had a reportable pathogenic/likely pathogenic variant as determined by whole-exome sequencing (WES), which was comparable to the proportions published previously in other countries. We further investigated the high carrier rate of rare variants in the ATP7B gene, which might indicate a founder effect in our population., Conclusion: This study was the first to provide Taiwanese population data of incidental findings and emphasized a high carrier rate of candidate pathogenic/likely pathogenic variants in the ATP7B gene., (© 2020 The Authors. Molecular Genetics & Genomic Medicine published by Wiley Periodicals LLC.)

Objective: We aimed to investigate the value of whole-exome sequencing (WES) in fetuses with congenital anomalies of the kidney and urinary tract (CAKUT) with or without other structural anomalies but with normal findings upon karyotyping and chromosome microarray analysis (CMA)., Methods: Cases with CAKUT with or without other structural anomalies were screened for eligibility. Fetuses with abnormal karyotyping or CMA results were excluded. We performed WES on DNA samples from eligible fetus-parental trios and identified diagnostic genetic variants based on ultrasonographic features., Results: A total of 163 eligible fetus-parental trios were successfully analyzed by WES. We found 26 likely pathogenic or pathogenic variants in 18 genes from 20 fetuses, with a total proportion of diagnostic genetic variants of 12.3% (20/163). Genetic variants were significantly more frequently detected in fetuses with multisystem anomalies (27.0%, 10/37), enlarged kidney/echogenic kidney (20%, 4/20), and multicystic dysplastic kidney (11.1%, 4/36). Pregnancy outcome data showed that 88 (94.6%, 88/93) of the surviving cases with negative WES results had a good prognosis in early childhood., Conclusions: Our study is the largest to use WES prenatally for CAKUT and shows that WES can be used diagnostically to define the molecular defects that underlie unexplained CAKUT., (© 2020 John Wiley & Sons, Ltd.)

Background: Mutations in the exonuclease domain of POLE, a DNA polymerase associated with DNA replication and repair, lead to cancers with ultra-high mutation rates. Most studies focus on intestinal and uterine cancers with POLE mutations. These cancers exhibit a significant immune cell infiltrate and favorable prognosis. We questioned whether loss of function of other DNA polymerases can cooperate to POLE to generate the ultramutator phenotype., Methods: We used cases and data from 15 cancer types in The Cancer Genome Atlas to investigate mutation frequencies of 14 different DNA polymerases. We tested whether tumor mutation burden, patient outcome (disease-free survival) and immune cell infiltration measured by ESTIMATE can be attributed to mutations in POLQ and POLZ/REV3L., Results: Thirty six percent of colorectal, stomach and endometrial cancers with POLE mutations carried additional mutations in POLQ (E/Q), POLZ/REV3L (E/Z) or both DNA polymerases (E/Z/Q). The mutation burden in these tumors was significantly greater compared to POLE-only (E) mutant tumors (p < 0.001). In addition, E/Q, E/Z, and E/Q/Z mutant tumors possessed an increased frequency of mutations in the POLE exonuclease domain (p = 0.013). Colorectal, stomach and endometrial E/Q, E/Z, and E/Q/Z mutant tumors within TCGA demonstrated 100% disease-free survival, even if the POLE mutations occurred outside the exonuclease domain (p = 0.003). However, immune scores in these tumors were related to microsatellite instability (MSI) and not POLE mutation status. This suggests that the host immune response may not be the sole mechanism for prolonged disease-free survival of ultramutated tumors in this cohort., Conclusion: Results in this study demonstrate that mutations in POLQ and REV3L in POLE mutant tumors should undergo further investigation to determine whether POLQ and REV3L mutations contribute to the ultramutator phenotype and favorable outcome of patients with POLE mutant tumors.

A high rate of consanguinity leads to a high prevalence of autosomal recessive disorders in inbred populations. One example of inbred populations is the Arab communities in Israel and the Palestinian Authority. In the Palestinian Authority in particular, due to limited access to specialized medical care, most patients do not receive a genetic diagnosis and can therefore neither receive genetic counseling nor possibly specific treatment. We used whole-exome sequencing as a first-line diagnostic tool in 83 Palestinian and Israeli Arab families with suspected neurogenetic disorders and were able to establish a probable genetic diagnosis in 51% of the families (42 families). Pathogenic, likely pathogenic or highly suggestive candidate variants were found in the following genes extending and refining the mutational and phenotypic spectrum of these rare disorders: ACO2, ADAT3, ALS2, AMPD2, APTX, B4GALNT1, CAPN1, CLCN1, CNTNAP1, DNAJC6, GAMT, GPT2, KCNQ2, KIF11, LCA5, MCOLN1, MECP2, MFN2, MTMR2, NT5C2, NTRK1, PEX1, POLR3A, PRICKLE1, PRKN, PRX, SCAPER, SEPSECS, SGCG, SLC25A15, SPG11, SYNJ1, TMCO1, and TSEN54. Further, this cohort has proven to be ideal for prioritization of new disease genes. Two separately published candidate genes (WWOX and PAX7) were identified in this study. Analyzing the runs of homozygosity (ROHs) derived from the Exome sequencing data as a marker for the rate of inbreeding, revealed significantly longer ROHs in the included families compared with a German control cohort. The total length of ROHs correlated with the detection rate of recessive disease-causing variants. Identification of the disease-causing gene led to new therapeutic options in four families.

Mitochondrial disorders (MDs) are genetic ailments affecting all age groups. Epidemiological data and frequencies of gene mutations in pediatric patients in China are scarce. This retrospective study assessed 101 patients with suspected MDs treated at the Neurology Department of Children's Hospital, Fudan University, in 2011-2017. Mitochondrial (mtDNA) and nuclear (nDNA) samples were assessed by long-range polymerase chain reaction (PCR)-based whole mtDNA sequencing and whole exome sequencing (WES) for identifying pathogenic mutations. Muscle samples underwent various staining protocols and immunofluorescence for detecting selected proteins. Seventeen mutations in the MT-TL1, MT-COX2, MT-ND4, MT, tRNA TRNE, MT-TN, MT-TK, MT-ATP6, MT-ND6, MT-ND3 and MT-CO3 genes were identified in 39 patients, of which m.3243A > G, m.3303C > T, m.8993T > C/G, m.9176T > C, and m.10191T > C were most common. Mitochondrial myopathy and MELAS were most common for m.3243A > G mutation. Four novel mutations were detected, including m.9478insT, m.5666T > C, m.8265T > C, and m.8380-13600 deletion mutations related to Leigh syndrome, mitochondrial myopathy and KSS, respectively. Thirty-three mutations in the TK2, POLG, IBA57, HADHB, FBXL4, ALDH5A1, FOXRED1, TPK1, NDUFAF5, NDUFAF7, NDUFV1, CARS2, PDHA1, and HIBCH genes were identified in 19 patients, including 23 currently unknown. Higher rates of TK2, POLG, IBA57, and HADHB mutations were found in nDNA-mutated MD compared with the remaining individuals. Besides, IBA57 c.286T > C (p.Y96H), TK2 c.497A > T (p.D166V) founder mutations critically contributed to MDs. Comprehensive genomic analysis plays a critical role in pediatric MD diagnosis. These data summarize the relative frequencies of different gene mutations in a large Chinese population, and identified 23 novel MD-associated nDNA and 4 novel mtDNA mutations., Competing Interests: Declaration of competing interest The authors declare that they have no competing interests., (Copyright © 2020. Published by Elsevier Masson SAS.)

As single-cell RNA sequencing technologies mature, massive gene expression profiles can be obtained. Consequently, cell clustering and annotation become two crucial and fundamental procedures affecting other specific downstream analyses. Most existing single-cell RNA-seq (scRNA-seq) data clustering algorithms do not take into account the available cell annotation results on the same tissues or organisms from other laboratories. Nonetheless, such data could assist and guide the clustering process on the target dataset. Identifying marker genes through differential expression analysis to manually annotate large amounts of cells also costs labor and resources. Therefore, in this paper, we propose a novel end-to-end cell supervised clustering and annotation framework called scAnCluster, which fully utilizes the cell type labels available from reference data to facilitate the cell clustering and annotation on the unlabeled target data. Our algorithm integrates deep supervised learning, self-supervised learning and unsupervised learning techniques together, and it outperforms other customized scRNA-seq supervised clustering methods in both simulation and real data. It is particularly worth noting that our method performs well on the challenging task of discovering novel cell types that are absent in the reference data.

Objective: To investigate the genetic causes and clinical outcomes of nonimmune hydrops fetalis (NIHF)., Methods: Cohort of cases of NIHF between July 2013 and December 2018. Initial genetic testing included quantitative fluorescence polymerase chain reaction for aneuploidies, karyotyping and chromosomal microarray analysis (CMA). In negative results, whole exome sequencing (WES) of the fetuses and parents was performed. Clinical post-natal follow-up assessments were conducted., Results: One hundred and nine patients fulfilled the study inclusion criteria and were sequentially genetically assessed by karyotype, CMA and WES. Among them, 24.8% (27/109) had a clinically significant genetic abnormality: 21 (19%) had abnormal karyotypes; 3/72 had pathogenic/likely pathogenic copy number variants (additional yield = 4.2%); and 3 had single gene disorders. The pregnancy termination and live birth rates of the cases with positive genetic testing results were significantly different from those with negative results (92.6% vs 53.7% and 3.7% vs 31.7%, respectively, P < .05 for both). During clinical follow-up of the survivors, 3/23 (13.0%) children developed an additional phenotype., Conclusion: This study improves our understanding of the diagnostic yield of CMA and WES for NIHF. A genetic diagnosis of NIHF can help determine the fetal prognosis and recurrence risk and influence pregnancy decision-making., (© 2020 John Wiley & Sons, Ltd.)

Heterogeneity is a fundamental feature of complex phenotypes. So far, genomic screenings have profiled thousands of samples providing insights into the transcriptome of the cell. However, disentangling the heterogeneity of these transcriptomic Big Data to identify defective biological processes remains challenging. Here we present GSECA, a method exploiting the bimodal behavior of RNA-sequencing gene expression profiles to identify altered gene sets in heterogeneous patient cohorts. Using simulated and experimental RNA-sequencing data sets, we show that GSECA provides higher performances than other available algorithms in detecting truly altered biological processes in large cohorts. Applied to 5941 samples from 14 different cancer types, GSECA correctly identified the alteration of the PI3K/AKT signaling pathway driven by the somatic loss of PTEN and verified the emerging role of PTEN in modulating immune-related processes. In particular, we showed that, in prostate cancer, PTEN loss appears to establish an immunosuppressive tumor microenvironment through the activation of STAT3, and low PTEN expression levels have a detrimental impact on patient disease-free survival. GSECA is available at https://github.com/matteocereda/GSECA., (© The Author(s) 2019. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Advantages and diagnostic effectiveness of the two most widely used resequencing approaches, whole exome (WES) and whole genome (WGS) sequencing, are often debated. WES dominated large-scale resequencing projects because of lower cost and easier data storage and processing. Rapid development of 3 rd generation sequencing methods and novel exome sequencing kits predicate the need for a robust statistical framework allowing informative and easy performance comparison of the emerging methods. In our study we developed a set of statistical tools to systematically assess coverage of coding regions provided by several modern WES platforms, as well as PCR-free WGS. We identified a substantial problem in most previously published comparisons which did not account for mappability limitations of short reads. Using regression analysis and simple machine learning, as well as several novel metrics of coverage evenness, we analyzed the contribution from the major determinants of CDS coverage. Contrary to a common view, most of the observed bias in modern WES stems from mappability limitations of short reads and exome probe design rather than sequence composition. We also identified the ~ 500 kb region of human exome that could not be effectively characterized using short read technology and should receive special attention during variant analysis. Using our novel metrics of sequencing coverage, we identified main determinants of WES and WGS performance. Overall, our study points out avenues for improvement of enrichment-based methods and development of novel approaches that would maximize variant discovery at optimal cost.

The MLL3 gene has been shown to be recurrently mutated in many malignancies including in families with acute myeloid leukemia. We demonstrate that many MLL3 variant calls made by exome sequencing are false positives due to misalignment to homologous regions, including a region on chr21, and can only be validated by long-range PCR. Numerous other recurrently mutated genes reported in COSMIC and TCGA databases have pseudogenes and cannot also be validated by conventional short read-based sequencing approaches. Genome-wide identification of pseudogene regions demonstrates that frequency of these homologous regions is increased with sequencing read lengths below 200 bps. To enable identification of poor quality sequencing variants in prospective studies, we generated novel genome-wide maps of regions with poor mappability that can be used in variant calling algorithms. Taken together, our findings reveal that pseudogene regions are a source of false-positive mutations in cancers.

Complex diseases such as inflammatory bowel disease (IBD), which consists of ulcerative colitis and Crohn's disease, are a significant medical burden-70 000 new cases of IBD are diagnosed in the United States annually. In this review, we examine the history of genetic variant discovery in complex disease with a focus on IBD. We cover methods that have been applied to microsatellite, common variant, targeted resequencing and whole-exome and -genome data, specifically focusing on the progression of technologies towards rare-variant discovery. The inception of these methods combined with better availability of population level variation data has led to rapid discovery of IBD-causative and/or -associated variants at over 200 loci; over time, these methods have grown exponentially in both power and ascertainment to detect rare variation. We highlight rare-variant discoveries critical to the elucidation of the pathogenesis of IBD, including those in NOD2, IL23R, CARD9, RNF186 and ADCY7. We additionally identify the major areas of rare-variant discovery that will evolve in the coming years. A better understanding of the genetic basis of IBD and other complex diseases will lead to improved diagnosis, prognosis, treatment and surveillance., (© The Author(s) 2019. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)