Focused Exome Sequencing Gives a High Diagnostic Yield in the Indian Subcontinent.

The genetically isolated yet heterogeneous and highly consanguineous Indian population has shown a higher prevalence of rare genetic disorders. However, there is a significant socioeconomic burden for genetic testing to be accessible to the general population. In the current study, we analyzed next-generation sequencing data generated through focused exome sequencing from individuals with different phenotypic manifestations referred for genetic testing to achieve a molecular diagnosis. Pathogenic or likely pathogenic variants are reported in 280 of 833 cases with a diagnostic yield of 33.6%. Homozygous sequence and copy number variants were found as positive diagnostic findings in 131 cases (15.7%) because of the high consanguinity in the Indian population. No relevant findings related to reported phenotype were identified in 6.2% of the cases. Patients referred for testing due to metabolic disorder and neuromuscular disorder had higher diagnostic yields. Carrier testing of asymptomatic individuals with a family history of the disease, through focused exome sequencing, achieved positive diagnosis in 54 of 118 cases tested. Copy number variants were also found in trans with single-nucleotide variants and mitochondrial variants in a few of the cases. The diagnostic yield and the findings from this study signify that a focused exome test is a good lower-cost alternative for whole-exome and whole-genome sequencing and as a first-tier approach to genetic testing.

DrugMechDB: A Curated Database of Drug Mechanisms.

Computational drug repositioning methods have emerged as an attractive and effective solution to find new candidates for existing therapies, reducing the time and cost of drug development. Repositioning methods based on biomedical knowledge graphs typically offer useful supporting biological evidence. This evidence is based on reasoning chains or subgraphs that connect a drug to a disease prediction. However, there are no databases of drug mechanisms that can be used to train and evaluate such methods. Here, we introduce the Drug Mechanism Database (DrugMechDB), a manually curated database that describes drug mechanisms as paths through a knowledge graph. DrugMechDB integrates a diverse range of authoritative free-text resources to describe 4,583 drug indications with 32,249 relationships, representing 14 major biological scales. DrugMechDB can be employed as a benchmark dataset for assessing computational drug repositioning models or as a valuable resource for training such models.

Combined sequence and copy number analysis improves diagnosis of limb girdle and other myopathies.

OBJECTIVE: Clinical and genetic heterogeneities make diagnosis of limb-girdle muscular dystrophy (LGMD) and other overlapping disorders of muscle weakness complicated and expensive. We aimed to develop a comprehensive next generation sequence-based multi-gene panel ("The Lantern Focused Neuromuscular Panel") to detect both sequence variants and copy number variants in one assay. METHODS: Patients with clinical diagnosis of LGMD or other overlapping muscular dystrophies in the United States were tested by PerkinElmer Genomics in 2018-2021 via "The Lantern Project," a sponsored diagnostic testing program. Sixty-six genes related to LGMD subtypes- and other myopathies were investigated. Main outcomes were diagnostic yield, gene-variant spectrum, and LGMD subtypes' prevalence. RESULTS: Molecular diagnosis was established in 19.6% (1266) of 6473 cases. Major genes contributing to LGMD were identified including CAPN3 (5.4%, 68), DYSF (4.0%, 51), GAA (3.7%, 47), ANO5 (3.6%, 45), and FKRP (2.7%, 34). Genes of other overlapping MD subtypes identified included PABPN1 (10.5%, 133), VCP (2.2%, 28), MYOT (1.2% 15), LDB3 (1.0%, 13), COL6A1 (1.5%, 19), FLNC (1.1%, 14), and DNAJB6 (0.8%, 10). Different sizes of copy number variants including single exon, multi-exon, and whole genes were identified in 7.5% (95) cases in genes including DMD, EMD, CAPN3, ANO5, SGCG, COL6A2, DOK7, and LAMA2. INTERPRETATION: "The Lantern Focused Neuromuscular Panel" enables identification of LGMD subtypes and other myopathies with overlapping clinical features. Prevalence of some MD subtypes was higher than previously reported. Widespread deployment of this comprehensive NGS panel has the potential to ensure early, accurate diagnosis as well as re-define MD epidemiology.

DrugMechDB: A Curated Database of Drug Mechanisms.

Computational drug repositioning methods have emerged as an attractive and effective solution to find new candidates for existing therapies, reducing the time and cost of drug development. Repositioning methods based on biomedical knowledge graphs typically offer useful supporting biological evidence. This evidence is based on reasoning chains or subgraphs that connect a drug to disease predictions. However, there are no databases of drug mechanisms that can be used to train and evaluate such methods. Here, we introduce the Drug Mechanism Database (DrugMechDB), a manually curated database that describes drug mechanisms as paths through a knowledge graph. DrugMechDB integrates a diverse range of authoritative free-text resources to describe 4,583 drug indications with 32,249 relationships, representing 14 major biological scales. DrugMechDB can be employed as a benchmark dataset for assessing computational drug repurposing models or as a valuable resource for training such models.

Pompe disease ascertained through The Lantern Project, 2018-2021: Next-generation sequencing and enzymatic testing to overcome obstacles to diagnosis.

The Lantern Project is an ongoing complimentary diagnostic program for patients in the United States sponsored by Sanofi and implemented by PerkinElmer Genomics. It combines specific enzymatic, biomarker, and genetic testing to facilitate rapid, accurate laboratory diagnosis of Pompe disease and several other lysosomal storage diseases, and a multigene next-generation sequencing panel including Pompe disease, LGMD, and other neuromuscular disorders. This article reports data for Pompe disease collected from October 2018 through December 2021, including acid α-glucosidase (GAA) enzyme assay and GAA sequencing (standard or expedited for positive newborn screening [NBS] to rule out infantile-onset Pompe disease [IOPD]) and the Focused Neuromuscular Panel, which includes GAA. One hundred forty patients (12 received only GAA enzyme testing, 128 had GAA sequencing alone or in addition to enzyme assay) have been confirmed with Pompe disease in this project. Eight of the 140 had a variant of unknown significance, but GAA activity ≤2.10 µmol/L/h, thus were confirmed with Pompe disease. Three diagnosed patients 0-2 years old had cross-reactive immunologic material (CRIM)-negative GAA variants and thus IOPD. One additional infant with presumptive IOPD had a homozygous frameshift c.1846del, likely CRIM-negative; symptoms were not provided. Among the 128 patients with molecular results, the c.-32-13T>G splice variant was homozygous in 11, compound-heterozygous in 98, and absent in 19. Proximal muscle weakness (58 patients) was the most common sign reported at testing; elevated creatine kinase (29 patients) was the most common laboratory result. The most common symptom categories were muscular (73 patients), musculoskeletal (13 patients), and respiratory (23 patients). Clinical information was not available for 42 samples, and 17 infants had only "abnormal NBS" or "low GAA" reported. Cardiac symptoms in 7 included potentially age-related conditions in five c.-32-13T>G-compound-heterozygous adults (myocardial infarction, heart murmur/palpitations, congestive heart failure: 1 each; 2 with atrial fibrillation) and hypertrophic cardiomyopathy in 2 children (1 and 2 years old) with presumptive IOPD. One novel GAA variant was observed in a patient with enzyme activity 0.31 µmol/L/h: c.1853_1854ins49, a frameshift pathogenic variant. The Lantern Project demonstrates the combinatorial utility of enzyme assay, targeted single-gene testing, and a focused neuromuscular next-generation sequencing panel in diagnosing Pompe disease.

Beyond the single average tumor: Understanding IO combinations using a clinical QSP model that incorporates heterogeneity in patient response.

A quantitative systems pharmacology model for metastatic melanoma was developed for immuno-oncology with the goal of predicting efficacy of combination checkpoint therapy with pembrolizumab and ipilimumab. This literature-based model is developed at multiple scales: (i) tumor and immune cell interactions at a lesion level; (ii) multiple heterogeneous target lesions, nontarget lesion growth, and appearance of new metastatic lesion at a patient level; and (iii) interpatient differences at a population level. The model was calibrated to pembrolizumab and ipilimumab monotherapy in patients with melanoma from Robert et al., specifically, waterfall plot showing target lesion response and overall response rate (Response Evaluation Criteria in Solid Tumors [RECIST] version 1.1), which additionally considers nontarget lesion growth and appearance of new metastatic lesions. We then used the model to predict waterfall and RECIST version 1.1 for combination treatment reported in Long et al. A key insight from this work was that nontarget lesions growth and appearance of new metastatic lesion contributed significantly to disease progression, despite reduction in target lesions. Further, the lesion level simulations of combination therapy show substantial efficacy in warm lesions (intermediary immunogenicity) but limited advantage of combination in both cold and hot lesions (low and high immunogenicity). Because many patients with metastatic disease are expected to have a mixture of these lesions, disease progression in such patients may be driven by a subset of cold lesions that are unresponsive to checkpoint inhibitors. These patients may benefit more from the combinations which include therapies to target cold lesions than double checkpoint inhibitors.

A single NGS-based assay covering the entire genomic sequence of the DMD gene facilitates diagnostic and newborn screening confirmatory testing.

Molecular diagnosis for Duchenne and Becker muscular dystrophies (DMD/BMD) involves a two-tiered approach for detection of deletions/duplications using MLPA or array CGH, followed by sequencing of coding and flanking intronic regions to detect sequence variants, which is time-consuming and expensive. We have developed a comprehensive next-generation sequencing (NGS)-based single-step assay to sequence the entire 2.2 Mb of the DMD gene to detect all copy number and sequence variants in both index males and carrier females. Assay validation was 100% concordant with other methodologies. A total of 772 samples have been tested, of which 62% (N = 480) were index cases with a clinical suspicion of DMD. Carrier testing females account for 38% (N = 292). Molecular diagnosis was confirmed in 86% (N = 413) of the index cases. Intragenic deletions and duplications (single-exon or multi-exon) were detected in 60% (N = 247) and 14% (N = 58) of the index cases, respectively. Full-sequence analysis of the entire gene allows for detection of deep intronic pathogenic variants and accurate breakpoint detection of CNVs involving similar exons, which could have an impact on the outcome of clinical trials. This comprehensive assay is highly sensitive for diagnostic testing for DMD and is also suitable for confirmatory testing for newborn screening for DMD.

Cadmium exposure upregulates SNAIL through miR-30 repression in human lung epithelial cells.

Cadmium (Cd) is a known human lung carcinogen. In addition, Cd exposure is associated with several lung diseases including emphysema, chronic obstructive pulmonary disease (COPD), asthma and fibrosis. Although earlier studies have identified several processes dysregulated by Cd exposure, the underlying mechanisms remain unclear. Here, we examined the transcriptome of lung epithelial cells exposed to Cd to understand the molecular basis of Cd-induced diseases. Computational analysis of the transcriptome predicted a significant number of Cd-upregulated genes to be targets of miR-30 family miRNAs. Experimental validation showed downregulation of all the miR-30 family members in Cd exposed cells. We found SNAIL, an EMT master regulator, to be the most upregulated among the miR-30 targets. Furthermore, we found decrease in the levels of epithelial marker E- cadherin (CDH1) and increase in the levels of mesenchymal markers, ZEB1 and vimentin. This suggested induction of EMT in Cd exposed cells. Luciferase reporter assays showed that miR-30 repressed SNAIL by directly targeting its 3' UTR. Over expression of miR-30e and transfection of miR-30e mimics reduced Cd-induced SNAIL upregulation. Our results suggest that miR-30 negatively regulates SNAIL in lung epithelial cells and that Cd-induced downregulation of miR-30 relieves this repression, resulting in SNAIL upregulation and EMT induction. EMT plays a major role in many diseases associated with Cd exposure including fibrosis, COPD, and cancer and metastasis. Therefore, our identification of miR-30 downregulation in Cd exposed cells and the consequent activation of SNAIL provides important mechanistic insights into lung diseases associated with Cd exposure.

Nickel exposure induces persistent mesenchymal phenotype in human lung epithelial cells through epigenetic activation of ZEB1.

Nickel (Ni) is an environmental and occupational carcinogen, and exposure to Ni is associated with lung and nasal cancers in humans. Furthermore, Ni exposure is implicated in several lung diseases including chronic inflammatory airway diseases, asthma, and fibrosis. However, the mutagenic potential of Ni is low and does not correlate with its potent toxicity and carcinogenicity. Therefore, mechanisms underlying Ni exposure-associated diseases remain poorly understood. Since the health risks of environmental exposures often continue post exposure, understanding the exposure effects that persist after the termination of exposure could provide mechanistic insights into diseases. By examining the persistent effects of Ni exposure, we report that Ni induces epithelial-mesenchymal transition (EMT) and that the mesenchymal phenotype remains irreversible even after the termination of exposure. Ni-induced EMT was dependent on the irreversible upregulation of ZEB1, an EMT master regulator, via resolution of its promoter bivalency. ZEB1, upon activation, downregulated its repressors as well as the cell-cell adhesion molecule, E-cadherin, resulting in the cells undergoing EMT and switching to persistent mesenchymal status. ZEB1 depletion in cells exposed to Ni attenuated Ni-induced EMT. Moreover, Ni exposure did not induce EMT in ZEB1-depleted cells. Activation of EMT, during which the epithelial cells lose cell-cell adhesion and become migratory and invasive, plays a major role in asthma, fibrosis, and cancer and metastasis, lung diseases associated with Ni exposure. Therefore, our finding of irreversible epigenetic activation of ZEB1 by Ni exposure and the acquisition of persistent mesenchymal phenotype would have important implications in understanding Ni-induced diseases.

Identification of a unique gene expression signature in mercury and 2,3,7,8-tetrachlorodibenzo-p-dioxin co-exposed cells.

Mercury (Hg) and 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) are major environmental contaminants that commonly co-occur in the environment. Both Hg and TCDD are associated with a number of human diseases including cancers. While the individual toxicological effects of Hg and TCDD have been extensively investigated, studies on co-exposure are limited to a few genes and pathways. Therefore, a significant knowledge gap exists in the understanding of the deleterious effects of co-exposure to Hg and TCDD. Due to the prevalence of Hg and TCDD co-contamination in the environment and the major human health hazards they pose, it is important to obtain a fuller understanding of genome-wide effects of Hg and TCDD co-exposure. In this study, by performing a comprehensive transcriptomic analysis of human bronchial epithelial cells (BEAS-2B) exposed to Hg and TCDD individually and in combination, we have uncovered a subset of genes with altered expression only in the co-exposed cells. We also identified the additive as well as antagonistic effects of Hg and TCDD on gene expression. Moreover, we found that co-exposure impacted several biological and disease processes not affected by Hg or TCDD individually. Our studies show that the consequences of Hg and TCDD co-exposure on the transcriptional program and biological processes could be substantially different from single exposures, thus providing new insights into the co-exposure-specific pathogenic processes.

Systems Biology Genetic Approach Identifies Serotonin Pathway as a Possible Target for Obstructive Sleep Apnea: Results from a Literature Search Review.

RATIONALE: Overall validity of existing genetic biomarkers in the diagnosis of obstructive sleep apnea (OSA) remains unclear. The objective of this systematic genetic study is to identify "novel" biomarkers for OSA using systems biology approach. METHODS: Candidate genes for OSA were extracted from PubMed, MEDLINE, and Embase search engines and DisGeNET database. The gene ontology (GO) analyses and candidate genes prioritization were performed using Enrichr tool. Genes pertaining to the top 10 pathways were extracted and used for Ingenuity Pathway Analysis. RESULTS: In total, we have identified 153 genes. The top 10 pathways associated with OSA include (i) serotonin receptor interaction, (ii) pathways in cancer, (iii) AGE-RAGE signaling in diabetes, (iv) infectious diseases, (v) serotonergic synapse, (vi) inflammatory bowel disease, (vii) HIF-1 signaling pathway, (viii) PI3-AKT signaling pathway, (ix) regulation lipolysis in adipocytes, and (x) rheumatoid arthritis. After removing the overlapping genes, we have identified 23 candidate genes, out of which >30% of the genes were related to the genes involved in the serotonin pathway. Among these 4 serotonin receptors SLC6A4, HTR2C, HTR2A, and HTR1B were strongly associated with OSA. CONCLUSIONS: This preliminary report identifies several potential candidate genes associated with OSA and also describes the possible regulatory mechanisms.

Oxidative stress under ambient and physiological oxygen tension in tissue culture.

Oxygen (O2) levels range from 2-9% in vivo. However, cell culture experiments are performed at atmospheric O2 levels (21%). Oxidative stress due to generation of reactive oxygen species (ROS) in cells cultured at higher than physiological levels is implicated in multitude of deleterious effects including DNA damage, genomic instability and senescence. In addition, oxidative stress activates redox sensitive transcription factors related to inflammatory signaling and apoptotic signaling. Furthermore, several chromatin-modifying enzymes are affected by ROS, potentially impacting epigenetic regulation of gene expression. While primary cells are cultured at lower O2 levels due to their inability to grow at higher O2, the immortalized cells, which display no such apparent growth difficulties, are typically cultured at 21% O2. This review will provide an overview of issues associated with increased oxygen levels in in vitro cell culture and point out the benefits of using lower levels of oxygen tension even for immortalized cells.

Nuclear Factor κB1/RelA Mediates Inflammation in Human Lung Epithelial Cells at Atmospheric Oxygen Levels.

Oxygen levels range from 2% to 9% in vivo. Atmospheric O2 levels (21%) are known to induce cell proliferation defects and cellular senescence in primary cell cultures. However, the mechanistic basis of the deleterious effects of higher O2 levels is not fully understood. On the other hand, immortalized cells including cancer cell lines, which evade cellular senescence are normally cultured at 21% O2 and the effects of higher O2 on these cells are understudied. Here, we addressed this problem by culturing immortalized human bronchial epithelial (BEAS-2B) cells at ambient atmospheric, 21% O2 and lower, 10% O2. Our results show increased inflammatory response at 21% O2 but not at 10% O2. We found higher RelA binding at the NF-κB1/RelA target gene promoters as well as upregulation of several pro-inflammatory cytokines in cells cultured at 21% O2. RelA knockdown prevented the upregulation of the pro-inflammatory cytokines at 21% O2, suggesting NF-κB1/RelA as a major mediator of inflammatory response in cells cultured at 21% O2. Interestingly, unlike the 21% O2 cultured cells, exposure of 10% O2 cultured cells to H2O2 did not elicit inflammatory response, suggesting increased ability to tolerate oxidative stress in cells cultured at lower O2 levels.

A mycobacterial phosphoribosyltransferase promotes bacillary survival by inhibiting oxidative stress and autophagy pathways in macrophages and zebrafish.

Mycobacterium tuberculosis employs various strategies to modulate host immune responses to facilitate its persistence in macrophages. The M. tuberculosis cell wall contains numerous glycoproteins with unknown roles in pathogenesis. Here, by using Concanavalin A and LC-MS analysis, we identified a novel mannosylated glycoprotein phosphoribosyltransferase, encoded by Rv3242c from M. tuberculosis cell walls. Homology modeling, bioinformatic analyses, and an assay of phosphoribosyltransferase activity in Mycobacterium smegmatis expressing recombinant Rv3242c (MsmRv3242c) confirmed the mass spectrometry data. Using Mycobacterium marinum-zebrafish and the surrogate MsmRv3242c infection models, we proved that phosphoribosyltransferase is involved in mycobacterial virulence. Histological and infection assays showed that the M. marinum mimG mutant, an Rv3242c orthologue in a pathogenic M. marinum strain, was strongly attenuated in adult zebrafish and also survived less in macrophages. In contrast, infection with wild type and the complemented ΔmimG:Rv3242c M. marinum strains showed prominent pathological features, such as severe emaciation, skin lesions, hemorrhaging, and more zebrafish death. Similarly, recombinant MsmRv3242c bacteria showed increased invasion in non-phagocytic epithelial cells and longer intracellular survival in macrophages as compared with wild type and vector control M. smegmatis strains. Further mechanistic studies revealed that the Rv3242c- and mimG-mediated enhancement of intramacrophagic survival was due to inhibition of autophagy, reactive oxygen species, and reduced activities of superoxide dismutase and catalase enzymes. Infection with MsmRv3242c also activated the MAPK pathway, NF-κB, and inflammatory cytokines. In summary, we show that a novel mycobacterial mannosylated phosphoribosyltransferase acts as a virulence and immunomodulatory factor, suggesting that it may constitute a novel target for antimycobacterial drugs.

Epigenetic dysregulation by nickel through repressive chromatin domain disruption.

Investigations into the genomic landscape of histone modifications in heterochromatic regions have revealed histone H3 lysine 9 dimethylation (H3K9me2) to be important for differentiation and maintaining cell identity. H3K9me2 is associated with gene silencing and is organized into large repressive domains that exist in close proximity to active genes, indicating the importance of maintenance of proper domain structure. Here we show that nickel, a nonmutagenic environmental carcinogen, disrupted H3K9me2 domains, resulting in the spreading of H3K9me2 into active regions, which was associated with gene silencing. We found weak CCCTC-binding factor (CTCF)-binding sites and reduced CTCF binding at the Ni-disrupted H3K9me2 domain boundaries, suggesting a loss of CTCF-mediated insulation function as a potential reason for domain disruption and spreading. We furthermore show that euchromatin islands, local regions of active chromatin within large H3K9me2 domains, can protect genes from H3K9me2-spreading-associated gene silencing. These results have major implications in understanding H3K9me2 dynamics and the consequences of chromatin domain disruption during pathogenesis.

Role of glycans and glycoproteins in disease development by Mycobacterium tuberculosis.

Glycoproteins play a critical role in host-pathogen interactions, antigenicity, and virulence determination, and are therefore, considered as potential drug targets. The cell wall of Mycobacterium tuberculosis (Mtb), the causative agent of tuberculosis (TB), dominantly contains sugars and lipids. Despite the efforts taken by the World Health Organization to reduce the incidence rate, the prevalence of TB is increasing in certain regions. This is mainly attributed to the emergence of multidrug-resistant bacteria. Factors that contribute to Mtb virulence and antigenicity remain elusive. However, several studies have shown that sugars and lipids are mainly responsible for Mtb pathogenesis and resistance to numerous drugs. This review gives insight into the role of glycoproteins in mycobacterium pathogenesis, disease development, and its implications in drug development.

Bio-informatics based analysis of genes implicated in alcohol mediated liver injury.

Alcohol induced liver injury has been studied extensively. Using literature search and bioinformatics tools, the present study characterizes the genes involved in alcohol induced liver injury. The cellular and metabolic processes in which genes involved in alcohol induced liver injury are implicated are also discussed. The genes related to alcohol induced liver injury are also involved in affecting certain molecular functions and metabolism of drugs, besides being associated with diseases. In conclusion, the changes in regulation of genes implicated in alcohol induced liver injury apart from causing alcohol mediated hepatic dysfunction may affect other vital processes in the body.