Enhancing clinical genomic accuracy with panelGC: a novel metric and tool for quantifying and monitoring GC biases in hybridization capture panel sequencing.

Accurate assessment of fragment abundance within a genome is crucial in clinical genomics applications such as the analysis of copy number variation (CNV). However, this task is often hindered by biased coverage in regions with varying guanine-cytosine (GC) content. These biases are particularly exacerbated in hybridization capture sequencing due to GC effects on probe hybridization and polymerase chain reaction (PCR) amplification efficiency. Such GC content-associated variations can exert a negative impact on the fidelity of CNV calling within hybridization capture panels. In this report, we present panelGC, a novel metric, to quantify and monitor GC biases in hybridization capture sequencing data. We establish the efficacy of panelGC, demonstrating its proficiency in identifying and flagging potential procedural anomalies, even in situations where instrument and experimental monitoring data may not be readily accessible. Validation using real-world datasets demonstrates that panelGC enhances the quality control and reliability of hybridization capture panel sequencing.

cSurvival: a web resource for biomarker interactions in cancer outcomes and in cell lines.

Survival analysis is a technique for identifying prognostic biomarkers and genetic vulnerabilities in cancer studies. Large-scale consortium-based projects have profiled >11 000 adult and >4000 pediatric tumor cases with clinical outcomes and multiomics approaches. This provides a resource for investigating molecular-level cancer etiologies using clinical correlations. Although cancers often arise from multiple genetic vulnerabilities and have deregulated gene sets (GSs), existing survival analysis protocols can report only on individual genes. Additionally, there is no systematic method to connect clinical outcomes with experimental (cell line) data. To address these gaps, we developed cSurvival (https://tau.cmmt.ubc.ca/cSurvival). cSurvival provides a user-adjustable analytical pipeline with a curated, integrated database and offers three main advances: (i) joint analysis with two genomic predictors to identify interacting biomarkers, including new algorithms to identify optimal cutoffs for two continuous predictors; (ii) survival analysis not only at the gene, but also the GS level; and (iii) integration of clinical and experimental cell line studies to generate synergistic biological insights. To demonstrate these advances, we report three case studies. We confirmed findings of autophagy-dependent survival in colorectal cancers and of synergistic negative effects between high expression of SLC7A11 and SLC2A1 on outcomes in several cancers. We further used cSurvival to identify high expression of the Nrf2-antioxidant response element pathway as a main indicator for lung cancer prognosis and for cellular resistance to oxidative stress-inducing drugs. Altogether, these analyses demonstrate cSurvival's ability to support biomarker prognosis and interaction analysis via gene- and GS-level approaches and to integrate clinical and experimental biomedical studies.

eVITTA: a web-based visualization and inference toolbox for transcriptome analysis.

Transcriptome profiling is essential for gene regulation studies in development and disease. Current web-based tools enable functional characterization of transcriptome data, but most are restricted to applying gene-list-based methods to single datasets, inefficient in leveraging up-to-date and species-specific information, and limited in their visualization options. Additionally, there is no systematic way to explore data stored in the largest transcriptome repository, NCBI GEO. To fill these gaps, we have developed eVITTA (easy Visualization and Inference Toolbox for Transcriptome Analysis; https://tau.cmmt.ubc.ca/eVITTA/). eVITTA provides modules for analysis and exploration of studies published in NCBI GEO (easyGEO), detailed molecular- and systems-level functional profiling (easyGSEA), and customizable comparisons among experimental groups (easyVizR). We tested eVITTA on transcriptomes of SARS-CoV-2 infected human nasopharyngeal swab samples, and identified a downregulation of olfactory signal transducers, in line with the clinical presentation of anosmia in COVID-19 patients. We also analyzed transcriptomes of Caenorhabditis elegans worms with disrupted S-adenosylmethionine metabolism, confirming activation of innate immune responses and feedback induction of one-carbon cycle genes. Collectively, eVITTA streamlines complex computational workflows into an accessible interface, thus filling the gap of an end-to-end platform capable of capturing both broad and granular changes in human and model organism transcriptomes.

Mediator subunit MDT-15/MED15 and Nuclear Receptor HIZR-1/HNF4 cooperate to regulate toxic metal stress responses in Caenorhabditis elegans.

Zinc is essential for cellular functions as it is a catalytic and structural component of many proteins. In contrast, cadmium is not required in biological systems and is toxic. Zinc and cadmium levels are closely monitored and regulated as their excess causes cell stress. To maintain homeostasis, organisms induce metal detoxification gene programs through stress responsive transcriptional regulatory complexes. In Caenorhabditis elegans, the MDT-15 subunit of the evolutionarily conserved Mediator transcriptional coregulator is required to induce genes upon exposure to excess zinc and cadmium. However, the regulatory partners of MDT-15 in this response, its role in cellular and physiological stress adaptation, and the putative role for mammalian MED15 in the metal stress responses remain unknown. Here, we show that MDT-15 interacts physically and functionally with the Nuclear Hormone Receptor HIZR-1 to promote molecular, cellular, and organismal adaptation to cadmium and excess zinc. Using gain- and loss-of-function mutants and qRT-PCR and reporter analysis, we find that mdt-15 and hizr-1 cooperate to induce zinc and cadmium responsive genes. Moreover, the two proteins interact physically in yeast-two-hybrid assays and this interaction is enhanced by the addition of zinc or cadmium, the former a known ligand of HIZR-1. Functionally, mdt-15 and hizr-1 mutants show defective storage of excess zinc in the gut and are hypersensitive to zinc-induced reductions in egg-laying. Furthermore, mdt-15 but not hizr-1 mutants are hypersensitive to cadmium-induced reductions in egg-laying, suggesting potential divergence of regulatory pathways. Lastly, mammalian MDT-15 orthologs bind genomic regulatory regions of metallothionein and zinc transporter genes in a cadmium and zinc-stimulated fashion, and human MED15 is required to induce a metallothionein gene in lung adenocarcinoma cells exposed to cadmium. Collectively, our data show that mdt-15 and hizr-1 cooperate to regulate cadmium detoxification and zinc storage and that this mechanism is at least partially conserved in mammals.

A "developmental hourglass" in fungi.

The "developmental hourglass" concept suggests that intermediate developmental stages are most resistant to evolutionary changes and that differences between species arise through divergence later in development. This high conservation during middevelopment is illustrated by the "waist" of the hourglass and it represents a low probability of evolutionary change. Earlier molecular surveys both on animals and on plants have shown that the genes expressed at the waist stage are more ancient and more conserved in their expression. The existence of such a developmental hourglass has not been explored in fungi, another eukaryotic kingdom. In this study, we generated a series of transcriptomic data covering the entire lifecycle of a model mushroom-forming fungus, Coprinopsis cinerea, and we observed a molecular hourglass over its development. The "young fruiting body" is the stage that expresses the evolutionarily oldest (lowest transcriptome age index) transcriptome and gives the strongest signal of purifying selection (lowest transcriptome divergence index). We also demonstrated that all three kingdoms-animals, plants, and fungi-display high expression levels of genes in "information storage and processing" at the waist stages, whereas the genes in "metabolism" become more highly expressed later. Besides, the three kingdoms all show underrepresented "signal transduction mechanisms" at the waist stages. The synchronic existence of a molecular "hourglass" across the three kingdoms reveals a mutual strategy for eukaryotes to incorporate evolutionary innovations.

Transcriptional coactivator MED15 is required for beta cell maturation.

Mediator, a co-regulator complex required for RNA Polymerase II activity, interacts with tissue-specific transcription factors to regulate development and maintain homeostasis. We observe reduced Mediator subunit MED15 expression in endocrine hormone-producing pancreatic islets isolated from people living with type 2 diabetes and sought to understand how MED15 and Mediator control gene expression programs important for the function of insulin-producing ß-cells. Here we show that Med15 is expressed during mouse ß-cell development and maturation. Knockout of Med15 in mouse ß-cells causes defects in ß-cell maturation without affecting ß-cell mass or insulin expression. ChIP-seq and co-immunoprecipitation analyses found that Med15 binds ß-cell transcription factors Nkx6-1 and NeuroD1 to regulate key ß-cell maturation genes. In support of a conserved role during human development, human embryonic stem cell-derived ß-like cells, genetically engineered to express high levels of MED15, express increased levels of maturation markers. We provide evidence of a conserved role for Mediator in ß-cell maturation and demonstrate an additional layer of control that tunes ß-cell transcription factor function.

Nuclear hormone receptor NHR-49 acts in parallel with HIF-1 to promote hypoxia adaptation in Caenorhabditis elegans.

The response to insufficient oxygen (hypoxia) is orchestrated by the conserved hypoxia-inducible factor (HIF). However, HIF-independent hypoxia response pathways exist that act in parallel with HIF to mediate the physiological hypoxia response. Here, we describe a hypoxia response pathway controlled by Caenorhabditis elegans nuclear hormone receptor NHR-49, an orthologue of mammalian peroxisome proliferator-activated receptor alpha (PPARα). We show that nhr-49 is required for animal survival in hypoxia and is synthetic lethal with hif-1 in this context, demonstrating that these factors act in parallel. RNA-seq analysis shows that in hypoxia nhr-49 regulates a set of genes that are hif-1-independent, including autophagy genes that promote hypoxia survival. We further show that nuclear hormone receptor nhr-67 is a negative regulator and homeodomain-interacting protein kinase hpk-1 is a positive regulator of the NHR-49 pathway. Together, our experiments define a new, essential hypoxia response pathway that acts in parallel with the well-known HIF-mediated hypoxia response.

Discovery of microRNA-like RNAs during early fruiting body development in the model mushroom Coprinopsis cinerea.

Coprinopsis cinerea is a model mushroom particularly suited for the study of fungal fruiting body development and the evolution of multicellularity in fungi. While microRNAs (miRNAs) have been extensively studied in animals and plants for their essential roles in post-transcriptional regulation of gene expression, miRNAs in fungi are less well characterized and their potential roles in controlling mushroom development remain unknown. To identify miRNA-like RNAs (milRNAs) in C. cinerea and explore their expression patterns during the early developmental transition of mushroom development, small RNA libraries of vegetative mycelium and primordium were generated and putative milRNA candidates were identified following the standards of miRNA prediction in animals and plants. Two out of 22 novel predicted milRNAs, cci-milR-12c and cci-milR-13e-5p, were validated by northern blot and stem-loop reverse transcription real-time PCR. Cci-milR-12c was differentially expressed whereas the expression levels of cci-milR-13e-5p were similar in the two developmental stages. Target prediction of the validated milRNAs resulted in genes associated with fruiting body development, including pheromone, hydrophobin, cytochrome P450, and protein kinase. Essential genes for miRNA biogenesis, including three coding for Dicer-like (DCL), one for Argonaute (AGO), one for AGO-like and one for quelling deficient-2 (QDE-2) proteins, were also identified in the C. cinerea genome. Phylogenetic analysis showed that the DCL and AGO proteins of C. cinerea were more closely related to those in other basidiomycetes and ascomycetes than to those in animals and plants. Taken together, our findings provided the first evidence for milRNAs in the model mushroom and their potential roles in regulating fruiting body development. New information on the evolutionary relationship of milRNA biogenesis proteins across kingdoms has also provided new insights for guiding further functional and evolutionary studies of miRNAs.

Population genomic analysis uncovers environmental stress-driven selection and adaptation of Lentinula edodes population in China.

The elucidation of genome-wide variations could help reveal aspects of divergence, domestication, and adaptation of edible mushrooms. Here, we resequenced the whole genomes of 39 wild and 21 cultivated strains of Chinese Lentinula edodes, the shiitake mushroom. We identified three distinct genetic groups in the Chinese L. edodes population with robust differentiation. Results of phylogenetic and population structure analyses suggest that the cultivated strains and most of the wild trains of L. edodes in China possess different gene pools and two outlier strains show signatures of hybridization between groups. Eighty-four candidate genes contributing to population divergence were detected in outlier analysis, 18 of which are involved in response to environmental stresses. Gene enrichment analysis of group-specific single nucleotide polymorphisms showed that the cultivated strains were genetically diversified in biological processes related to stress response. As the formation of fruiting bodies is a stress-response process, we postulate that environment factors, such as temperature, drove the population divergence of L. edodes in China by natural or artificial selection. We also found phenotypic variations between groups and identified some wild strains that have potential to diversify the genetic pool for improving agricultural traits of L. edodes cultivars in China.

Loss of ARID1A expression leads to sensitivity to ROS-inducing agent elesclomol in gynecologic cancer cells.

Inactivating mutations in ARID1A are found in a broad spectrum of cancer types, with the highest frequency in gynecologic cancers. However, therapeutic strategies targeting ARID1A-mutant cancer cells remain limited. In this study, we aimed to identify drugs sensitivities in ARID1A-mutant cancer cell lines. By analyzing the Genomics of Drug Sensitivity in Cancer database, we found that ARID1A-mutant cancer cell lines were more sensitive to treatment with the reactive oxygen species (ROS)-inducing agent elesclomol. In a panel of 14 gynecologic cancer cell lines, treatment with elesclomol inhibited growth and induced apoptosis more potently in ARID1A-mutant cells. Knockdown of ARID1A in RMG1 and OVCA432 ovarian cancer cells resulted in increased sensitivity to elesclomol, whereas restoration of ARID1A expression in TOV21G ovarian cancer cells resulted in increased resistance to elesclomol. Furthermore, we found that knockdown of ARID1A expression resulted in increased intracellular ROS levels. In ovarian clear cell carcinoma patient samples, low expression of ARID1A correlated with high expression of 8-hydroxyguanosine, a marker for oxidative stress. In summary, we demonstrate for the first time that loss of ARID1A leads to accumulation of ROS and suggest that elesclomol may be used to target ARID1A-mutant gynecologic cancer cells.

High-density transcriptional initiation signals underline genomic islands in bacteria.

Genomic islands (GIs), frequently associated with the pathogenicity of bacteria and having a substantial influence on bacterial evolution, are groups of "alien" elements which probably undergo special temporal-spatial regulation in the host genome. Are there particular hallmark transcriptional signals for these "exotic" regions? We here explore the potential transcriptional signals that underline the GIs beyond the conventional views on basic sequence composition, such as codon usage and GC property bias. It showed that there is a significant enrichment of the transcription start positions (TSPs) in the GI regions compared to the whole genome of Salmonella enterica and Escherichia coli. There was up to a four-fold increase for the 70% GIs, implying high-density TSPs profile can potentially differentiate the GI regions. Based on this feature, we developed a new sliding window method GIST, Genomic-island Identification by Signals of Transcription, to identify these regions. Subsequently, we compared the known GI-associated features of the GIs detected by GIST and by the existing method Islandviewer to those of the whole genome. Our method demonstrates high sensitivity in detecting GIs harboring genes with biased GI-like function, preferred subcellular localization, skewed GC property, shorter gene length and biased "non-optimal" codon usage. The special transcriptional signals discovered here may contribute to the coordinate expression regulation of foreign genes. Finally, by using GIST, we detected many interesting GIs in the 2011 German E. coli O104:H4 outbreak strain TY-2482, including the microcin H47 system and gene cluster ycgXEFZ-ymgABC that activates the production of biofilm matrix. The aforesaid findings highlight the power of GIST to predict GIs with distinct intrinsic features to the genome. The heterogeneity of cumulative TSPs profiles may not only be a better identity for "alien" regions, but also provide hints to the special evolutionary course and transcriptional regulation of GI regions.