We describe the application of a microarray platform, which combines information from exon body and splice-junction probes, to perform a quantitative analysis of tissue-specific alternative splicing (AS) for thousands of exons in mammalian cells. Through this system, we have analyzed global features of AS in major mouse tissues. The results provide numerous inferences for the functions of tissue-specific AS, insights into how the evolutionary history of exons can impact on their inclusion levels, and also information on how global regulatory properties of AS define tissue type. Like global transcription profiles, global AS profiles reflect tissue identity. Interestingly, we find that transcription and AS act independently on different sets of genes in order to define tissue-specific expression profiles. These results demonstrate the utility of our quantitative microarray platform and data for revealing important global regulatory features of AS.

Although subcellular mRNA trafficking has been demonstrated as a mechanism to control protein distribution, it is generally believed that most protein localization occurs subsequent to translation. To address this point, we developed and employed a high-resolution fluorescent in situ hybridization procedure to comprehensively evaluate mRNA localization dynamics during early Drosophila embryogenesis. Surprisingly, of the 3370 genes analyzed, 71% of those expressed encode subcellularly localized mRNAs. Dozens of new and striking localization patterns were observed, implying an equivalent variety of localization mechanisms. Tight correlations between mRNA distribution and subsequent protein localization and function, indicate major roles for mRNA localization in nucleating localized cellular machineries. A searchable web resource documenting mRNA expression and localization dynamics has been established and will serve as an invaluable tool for dissecting localization mechanisms and for predicting gene functions and interactions.

In animals, egg activation triggers a cascade of posttranscriptional events that act on maternally synthesized RNAs. We show that, in Drosophila, the PAN GU (PNG) kinase sits near the top of this cascade, triggering translation of SMAUG (SMG), a multifunctional posttranscriptional regulator conserved from yeast to humans. Although PNG is required for cytoplasmic polyadenylation of smg mRNA, it regulates translation via mechanisms that are independent of its effects on the poly(A) tail. Analyses of mutants suggest that PNG relieves translational repression by PUMILIO (PUM) and one or more additional factors, which act in parallel through the smg mRNA's 3' untranslated region (UTR). Microarray-based gene expression profiling shows that SMG is a major regulator of maternal transcript destabilization. SMG-dependent mRNAs are enriched for gene ontology annotations for function in the cell cycle, suggesting a possible causal relationship between failure to eliminate these transcripts and the cell cycle defects in smg mutants.

Identification of non-coding mutations driving tumorigenesis requires alternative approaches to coding mutations. Enriched associations between mutated regulatory elements and altered cis-regulation in tumors are a promising approach to stratify candidate non-coding driver mutations. Here we provide a bioinformatics pipeline to mine data from the Cancer Genomic Commons (GDC) for such associations. The pipeline integrates RNA and whole-genome sequencing with genotyping data to reveal putative non-coding driver mutations by cancer type. For complete information on the generation and use of this protocol, please refer to Cheng et al. (2021).

Identifying associations between genotypes and gene expression levels using microarrays has enabled systematic interrogation of regulatory variation underlying complex phenotypes. This approach has vast potential for functional characterization of disease states, but its prohibitive cost, given hundreds to thousands of individual samples from populations have to be genotyped and expression profiled, has limited its widespread application.

Oct4 is a widely recognized pluripotency factor as it maintains Embryonic Stem (ES) cells in a pluripotent state, and, in vivo, prevents the inner cell mass (ICM) in murine embryos from differentiating into trophectoderm. However, its function in somatic tissue after this developmental stage is not well characterized. Using a tamoxifen-inducible Cre recombinase and floxed alleles of Oct4, we investigated the effect of depleting Oct4 in mouse embryos between the pre-streak and headfold stages, ~E6.0-E8.0, when Oct4 is found in dynamic patterns throughout the embryonic compartment of the mouse egg cylinder. We found that depletion of Oct4 ~E7.5 resulted in a severe phenotype, comprised of craniorachischisis, random heart tube orientation, failed turning, defective somitogenesis and posterior truncation. Unlike in ES cells, depletion of the pluripotency factors Sox2 and Oct4 after E7.0 does not phenocopy, suggesting that ~E7.5 Oct4 is required within a network that is altered relative to the pluripotency network. Oct4 is not required in extraembryonic tissue for these processes, but is required to maintain cell viability in the embryo and normal proliferation within the primitive streak. Impaired expansion of the primitive streak occurs coincident with Oct4 depletion ∼E7.5 and precedes deficient convergent extension which contributes to several aspects of the phenotype.

During early embryonic development, one of the two X chromosomes in mammalian female cells is inactivated to compensate for a potential imbalance in transcript levels with male cells, which contain a single X chromosome. Here, we use mouse female embryonic stem cells (ESCs) with non-random X chromosome inactivation (XCI) and polymorphic X chromosomes to study the dynamics of gene silencing over the inactive X chromosome by high-resolution allele-specific RNA-seq.

Genomic imprinting is an epigenetic process that restricts gene expression to either the maternally or paternally inherited allele. Many theories have been proposed to explain its evolutionary origin, but understanding has been limited by a paucity of data mapping the breadth and dynamics of imprinting within any organism. We generated an atlas of imprinting spanning 33 mouse and 45 human developmental stages and tissues. Nearly all imprinted genes were imprinted in early development and either retained their parent-of-origin expression in adults or lost it completely. Consistent with an evolutionary signature of parental conflict, imprinted genes were enriched for coexpressed pairs of maternally and paternally expressed genes, showed accelerated expression divergence between human and mouse, and were more highly expressed than their non-imprinted orthologs in other species. Our approach demonstrates a general framework for the discovery of imprinting in any species and sheds light on the causes and consequences of genomic imprinting in mammals.

We developed PolyA-seq, a strand-specific and quantitative method for high-throughput sequencing of 3' ends of polyadenylated transcripts, and used it to globally map polyadenylation (polyA) sites in 24 matched tissues in human, rhesus, dog, mouse, and rat. We show that PolyA-seq is as accurate as existing RNA sequencing (RNA-seq) approaches for digital gene expression (DGE), enabling simultaneous mapping of polyA sites and quantitative measurement of their usage. In human, we confirmed 158,533 known sites and discovered 280,857 novel sites (FDR < 2.5%). On average 10% of novel human sites were also detected in matched tissues in other species. Most novel sites represent uncharacterized alternative polyA events and extensions of known transcripts in human and mouse, but primarily delineate novel transcripts in the other three species. A total of 69.1% of known human genes that we detected have multiple polyA sites in their 3'UTRs, with 49.3% having three or more. We also detected polyadenylation of noncoding and antisense transcripts, including constitutive and tissue-specific primary microRNAs. The canonical polyA signal was strongly enriched and positionally conserved in all species. In general, usage of polyA sites is more similar within the same tissues across different species than within a species. These quantitative maps of polyA usage in evolutionarily and functionally related samples constitute a resource for understanding the regulatory mechanisms underlying alternative polyadenylation.

Accurate identification of novel, functional noncoding (nc) RNA features in genome sequence has proven more difficult than for exons. Current algorithms identify and score potential RNA secondary structures on the basis of thermodynamic stability, conservation, and/or covariance in sequence alignments. Neither the algorithms nor the information gained from the individual inputs have been independently assessed. Furthermore, due to issues in modelling background signal, it has been difficult to gauge the precision of these algorithms on a genomic scale, in which even a seemingly small false-positive rate can result in a vast excess of false discoveries.

Candida albicans is the most common cause of systemic fungal infections in humans and is considerably more virulent than its closest known relative, Candida dubliniensis. To investigate this difference, we constructed interspecies hybrids and quantified mRNA levels produced from each genome in the hybrid. This approach systematically identified expression differences in orthologous genes arising from cis-regulatory sequence changes that accumulated since the two species last shared a common ancestor, some 10 million y ago. We documented many orthologous gene-expression differences between the two species, and we pursued one striking observation: All 15 genes coding for the enzymes of glycolysis showed higher expression from the C. albicans genome than the C. dubliniensis genome in the interspecies hybrid. This pattern requires evolutionary changes to have occurred at each gene; the fact that they all act in the same direction strongly indicates lineage-specific natural selection as the underlying cause. To test whether these expression differences contribute to virulence, we created a C. dubliniensis strain in which all 15 glycolysis genes were produced at modestly elevated levels and found that this strain had significantly increased virulence in the standard mouse model of systemic infection. These results indicate that small expression differences across a deeply conserved set of metabolism enzymes can play a significant role in the evolution of virulence in fungal pathogens.

The idea that most morphological adaptations can be attributed to changes in the cis-regulation of gene expression levels has been gaining increasing acceptance, despite the fact that only a handful of such cases have so far been demonstrated. Moreover, because each of these cases involves only one gene, we lack any understanding of how natural selection may act on cis-regulation across entire pathways or networks. Here we apply a genome-wide test for selection on cis-regulation to two subspecies of the mouse Mus musculus. We find evidence for lineage-specific selection at over 100 genes involved in diverse processes such as growth, locomotion, and memory. These gene sets implicate candidate genes that are supported by both quantitative trait loci and a validated causality-testing framework, and they predict a number of phenotypic differences, which we confirm in all four cases tested. Our results suggest that gene expression adaptation is widespread and that these adaptations can be highly polygenic, involving cis-regulatory changes at numerous functionally related genes. These coordinated adaptations may contribute to divergence in a wide range of morphological, physiological, and behavioral phenotypes.

The evolution of drug resistance has a profound impact on human health. Candida glabrata is a leading human fungal pathogen that can rapidly evolve resistance to echinocandins, which target cell wall biosynthesis and are front-line therapeutics for Candida infections. Here, we provide the first global analysis of mutations accompanying the evolution of fungal drug resistance in a human host utilizing a series of C. glabrata isolates that evolved echinocandin resistance in a patient treated with the echinocandin caspofungin for recurring bloodstream candidemia. Whole genome sequencing identified a mutation in the drug target, FKS2, accompanying a major resistance increase, and 8 additional non-synonymous mutations. The FKS2-T1987C mutation was sufficient for echinocandin resistance, and associated with a fitness cost that was mitigated with further evolution, observed in vitro and in a murine model of systemic candidemia. A CDC6-A511G(K171E) mutation acquired before FKS2-T1987C(S663P), conferred a small resistance increase. Elevated dosage of CDC55, which acquired a C463T(P155S) mutation after FKS2-T1987C(S663P), ameliorated fitness. To discover strategies to abrogate echinocandin resistance, we focused on the molecular chaperone Hsp90 and downstream effector calcineurin. Genetic or pharmacological compromise of Hsp90 or calcineurin function reduced basal tolerance and resistance. Hsp90 and calcineurin were required for caspofungin-dependent FKS2 induction, providing a mechanism governing echinocandin resistance. A mitochondrial respiration-defective petite mutant in the series revealed that the petite phenotype does not confer echinocandin resistance, but renders strains refractory to synergy between echinocandins and Hsp90 or calcineurin inhibitors. The kidneys of mice infected with the petite mutant were sterile, while those infected with the HSP90-repressible strain had reduced fungal burden. We provide the first global view of mutations accompanying the evolution of fungal drug resistance in a human host, implicate the premier compensatory mutation mitigating the cost of echinocandin resistance, and suggest a new mechanism of echinocandin resistance with broad therapeutic potential.

Systematic identification and functional characterization of novel types of noncoding (nc)RNA in genomes is more difficult than it is for protein coding mRNAs, since ncRNAs typically do not possess sequence features such as splicing or translation signals, or long open reading frames. Recent "tiling" microarray studies have reported that a surprisingly larger proportion of mammalian genomes is transcribed than was previously anticipated. However, these non-genic transcripts often appear to be low in abundance, and their functional significance is not known.

We used a microarray containing probes that tile all known yeast noncoding RNAs (ncRNAs) to investigate RNA biogenesis on a global scale. The microarray verified a general loss of Box C/D snoRNAs in the TetO7-BCD1 mutant, which had previously been shown for only a handful of snoRNAs. We also monitored the accumulation of improperly processed flank sequences of pre-RNAs in strains depleted for known RNA nucleases, including RNase III, Dbr1p, Xrn1p, Rat1p and components of the exosome and RNase P complexes. Among the hundreds of aberrant RNA processing events detected, two novel substrates of Rnt1p (the RUF1 and RUF3 snoRNAs) were identified. We also identified a relationship between tRNA 5' end processing and tRNA splicing, processes that were previously thought to be independent. This analysis demonstrates the applicability of microarray technology to the study of global analysis of ncRNA synthesis and provides an extensive directory of processing events mediated by yeast ncRNA processing enzymes.

Despite the recent availability of complete genome sequences of tumors from thousands of patients, isolating disease-causing (driver) non-coding mutations from the plethora of somatic variants remains challenging, and only a handful of validated examples exist. By integrating whole-genome sequencing, genetic data, and allele-specific gene expression from TCGA, we identified 320 somatic non-coding mutations that affect gene expression in cis (FDR<0.25). These mutations cluster into 47 cis-regulatory elements that modulate expression of their subject genes through diverse molecular mechanisms. We further show that these mutations have hallmark features of non-coding drivers; namely, that they preferentially disrupt transcription factor binding motifs, are associated with a selective advantage, increased oncogene expression and decreased tumor suppressor expression.

Mapping genetic interactions in mammalian cells is limited due to technical obstacles. Here we describe a method called TCGI (tRNA-CRISPR for genetic interactions) to generate a high-efficient, barcode-free and scalable pairwise CRISPR libraries in mammalian cells for identifying genetic interactions. We have generated a genome- wide library to identify genes genetically interacting with TAZ in cell viability regulation. Validation of candidate synergistic genes reveals the screening accuracy of 85% and TAZ-MCL1 is characterized as combinational drug targets for non-small cell lung cancer treatments. TCGI has dramatically improved the current methods for mapping genetic interactions and screening drug targets for combinational therapies.

Maternal diabetes and high-fat feeding during pregnancy have been linked to later life outcomes in offspring. To investigate the effects of both maternal and paternal hyperglycemia on offspring phenotypes, we utilized an autosomal dominant mouse model of diabetes (hypoinsulinemic hyperglycemia in Akita mice). We determined metabolic and skeletal phenotypes in wildtype offspring of Akita mothers and fathers.

Vertebrates share the same general body plan and organs, possess related sets of genes, and rely on similar physiological mechanisms, yet show great diversity in morphology, habitat and behavior. Alteration of gene regulation is thought to be a major mechanism in phenotypic variation and evolution, but relatively little is known about the broad patterns of conservation in gene expression in non-mammalian vertebrates.