Genome Research current issue

Nascent transcript analysis of glucocorticoid crosstalk with TNF defines primary and cooperative inflammatory repression [RESEARCH]

2019-11-01T08:06:42-07:00

The glucocorticoid receptor (NR3C1, also known as GR) binds to specific DNA sequences and directly induces transcription of anti-inflammatory genes that contribute to cytokine repression, frequently in cooperation with NF-kB. Whether inflammatory repression also occurs through local interactions between GR and inflammatory gene regulatory elements has been controversial. Here, using global run-on sequencing (GRO-seq) in human airway epithelial cells, we show that glucocorticoid signaling represses transcription within 10 min. Many repressed regulatory regions reside within "hyper-ChIPable" genomic regions that are subject to dynamic, yet nonspecific, interactions with some antibodies. When this artifact was accounted for, we determined that transcriptional repression does not require local GR occupancy. Instead, widespread transcriptional induction through canonical GR binding sites is associated with reciprocal repression of distal TNF-regulated enhancers through a chromatin-dependent process, as evidenced by chromatin accessibility and motif displacement analysis. Simultaneously, transcriptional induction of key anti-inflammatory effectors is decoupled from primary repression through cooperation between GR and NF-kB at a subset of regulatory regions. Thus, glucocorticoids exert bimodal restraints on inflammation characterized by rapid primary transcriptional repression without local GR occupancy and secondary anti-inflammatory effects resulting from transcriptional cooperation between GR and NF-kB.

A-to-I RNA editing contributes to the persistence of predicted damaging mutations in populations [RESEARCH]

Mai, T.-L., Chuang, T.-J. — 2019-11-01T08:06:42-07:00

Adenosine-to-inosine (A-to-I) RNA editing is a very common co-/posttranscriptional modification that can lead to A-to-G changes at the RNA level and compensate for G-to-A genomic changes to a certain extent. It has been shown that each healthy individual can carry dozens of missense variants predicted to be severely deleterious. Why strongly detrimental variants are preserved in a population and not eliminated by negative natural selection remains mostly unclear. Here, we ask if RNA editing correlates with the burden of deleterious A/G polymorphisms in a population. Integrating genome and transcriptome sequencing data from 447 human lymphoblastoid cell lines, we show that nonsynonymous editing activities (prevalence/level) are negatively correlated with the deleteriousness of A-to-G genomic changes and positively correlated with that of G-to-A genomic changes within the population. We find a significantly negative correlation between nonsynonymous editing activities and allele frequency of A within the population. This negative editing-allele frequency correlation is particularly strong when editing sites are located in highly important genes/loci. Examinations of deleterious missense variants from the 1000 Genomes Project further show a significantly higher proportion of rare missense mutations for G-to-A changes than for other types of changes. The proportion for G-to-A changes increases with increasing deleterious effects of the changes. Moreover, the deleteriousness of G-to-A changes is significantly positively correlated with the percentage of editing enzyme binding motifs at the variants. Overall, we show that nonsynonymous editing is associated with the increased burden of G-to-A missense mutations in healthy individuals, expanding RNA editing in pathogenomics studies.

Global analyses of the dynamics of mammalian microRNA metabolism [RESEARCH]

Kingston, E. R., Bartel, D. P. — 2019-11-01T08:06:42-07:00

Rates of production and degradation together specify microRNA (miRNA) abundance and dynamics. Here, we used approach-to-steady-state metabolic labeling to assess these rates for 176 miRNAs in contact-inhibited mouse embryonic fibroblasts (MEFs), 182 miRNAs in dividing MEFs, and 127 miRNAs in mouse embryonic stem cells (mESCs). MicroRNA duplexes, each comprising a mature miRNA and its passenger strand, are produced at rates as fast as 110 ± 50 copies/cell/min, which exceeds rates reported for any mRNAs. These duplexes are rapidly loaded into Argonaute, with <30 min typically required for duplex loading and silencing-complex maturation. Within Argonaute, guide strands have stabilities that vary by 100-fold. Half-lives also vary globally between cell lines, with median values ranging from 11 to 34 h in mESCs and contact-inhibited MEFs, respectively. Moreover, relative half-lives for individual miRNAs vary between cell types, implying the influence of cell-specific factors in dictating turnover rate. The apparent influence of miRNA regions most important for targeting, together with the effect of one target on miR-7 accumulation, suggest that targets fulfill this role. Analysis of the tailing and trimming of miRNA 3' termini showed that the flux was typically greatest through the isoform tailed with a single uridine, although changes in this flux did not correspond to changes in stability, which suggested that the processes of tailing and trimming might be independent from that of decay. Together, these results establish a framework for describing the dynamics and regulation of miRNAs throughout their life cycle.

Cotargeting among microRNAs in the brain [RESEARCH]

Cherone, J. M., Jorgji, V., Burge, C. B. — 2019-11-01T08:06:42-07:00

MicroRNAs (miRNAs) play roles in diverse developmental and disease processes. Distinct miRNAs have hundreds to thousands of conserved mRNA binding sites but typically direct only modest repression via single sites. Cotargeting of individual mRNAs by different miRNAs could potentially achieve stronger and more complex patterns of repression. By comparing target sets of different miRNAs, we identified hundreds of pairs of miRNAs that share more mRNA targets than expected (often by twofold or more) relative to stringent controls. Genetic perturbations revealed a functional overlap in neuronal differentiation for the cotargeting pair miR-138/miR-137. Clustering of all cotargeting pairs revealed a group of nine predominantly brain-enriched miRNAs that share many targets. In reporter assays, subsets of these miRNAs together repressed gene expression by five- to 10-fold, often showing cooperative repression. Together, our results uncover an unexpected pattern in which combinations of miRNAs collaborate to robustly repress cotargets, and suggest important developmental roles for cotargeting.

The subgenomes show asymmetric expression of alleles in hybrid lineages of Megalobrama amblycephala x Culter alburnus [RESEARCH]

2019-11-01T08:06:42-07:00

Hybridization drives rapid speciation by shaping novel genotypic and phenotypic profiles. Genomic incompatibility and transcriptome shock have been observed in hybrids, although this is rarer in animals than in plants. Using the newly sequenced genomes of the blunt snout bream (Megalobrama amblycephala [BSB]) and the topmouth culter (Culter alburnus [TC]), we focused on the sequence variation and gene expression changes in the reciprocal intergeneric hybrid lineages (F₁–F₃) of BSB x TC. A genome-wide transcriptional analysis identified 145–974 expressed recombinant genes in the successive generations of hybrid fish, suggesting the rapid emergence of allelic variation following hybridization. Some gradual changes of gene expression with additive and dominance effects and various cis and trans regulations were observed from F₁ to F₃ in the two hybrid lineages. These asymmetric patterns of gene expression represent the alternative strategies for counteracting deleterious effects of the subgenomes and improving adaptability of novel hybrids. Furthermore, we identified positive selection and additive expression patterns in transforming growth factor, beta 1b (tgfb1b), which may account for the morphological variations of the pharyngeal jaw in the two hybrid lineages. Our current findings provide insights into the evolution of vertebrate genomes immediately following hybridization.

Gene expression profiling of single cells from archival tissue with laser-capture microdissection and Smart-3SEQ [METHOD]

2019-11-01T08:06:42-07:00

RNA sequencing (RNA-seq) is a sensitive and accurate method for quantifying gene expression. Small samples or those whose RNA is degraded, such as formalin-fixed paraffin-embedded (FFPE) tissue, remain challenging to study with nonspecialized RNA-seq protocols. Here, we present a new method, Smart-3SEQ, that accurately quantifies transcript abundance even with small amounts of total RNA and effectively characterizes small samples extracted by laser-capture microdissection (LCM) from FFPE tissue. We also obtain distinct biological profiles from FFPE single cells, which have been impossible to study with previous RNA-seq protocols, and we use these data to identify possible new macrophage phenotypes associated with the tumor microenvironment. We propose Smart-3SEQ as a highly cost-effective method to enable large gene expression profiling experiments unconstrained by sample size and tissue availability. In particular, Smart-3SEQ's compatibility with FFPE tissue unlocks an enormous number of archived clinical samples; combined with LCM it allows unprecedented studies of small cell populations and single cells isolated by their in situ context.

FFPEcap-seq: a method for sequencing capped RNAs in formalin-fixed paraffin-embedded samples [METHOD]

2019-11-01T08:06:42-07:00

The majority of clinical cancer specimens are preserved as formalin-fixed paraffin-embedded (FFPE) samples. For clinical molecular tests to have wide-reaching impact, they must be applicable to FFPE material. Accurate quantitative measurements of RNA derived from FFPE specimens is challenging because of low yields and high amounts of degradation. Here, we present FFPEcap-seq, a method specifically designed for sequencing capped 5' ends of RNA derived from FFPE samples. FFPEcap-seq combines enzymatic enrichment of 5' capped RNAs with template switching to create sequencing libraries. We find that FFPEcap-seq can faithfully capture mRNA expression levels in FFPE specimens while also detecting enhancer RNAs that arise from distal regulatory regions. FFPEcap-seq is a fast and straightforward method for making high-quality 5' end RNA-seq libraries from FFPE-derived RNA.

Identification and dynamic quantification of regulatory elements using total RNA [METHOD]

Duttke, S. H., Chang, M. W., Heinz, S., Benner, C. — 2019-11-01T08:06:42-07:00

The spatial and temporal regulation of transcription initiation is pivotal for controlling gene expression. Here, we introduce capped-small RNA-seq (csRNA-seq), which uses total RNA as starting material to detect transcription start sites (TSSs) of both stable and unstable RNAs at single-nucleotide resolution. csRNA-seq is highly sensitive to acute changes in transcription and identifies an order of magnitude more regulated transcripts than does RNA-seq. Interrogating tissues from species across the eukaryotic kingdoms identified unstable transcripts resembling enhancer RNAs, pri-miRNAs, antisense transcripts, and promoter upstream transcripts in multicellular animals, plants, and fungi spanning 1.6 billion years of evolution. Integration of epigenomic data from these organisms revealed that histone H3 trimethylation (H3K4me3) was largely confined to TSSs of stable transcripts, whereas H3K27ac marked nucleosomes downstream from all active TSSs, suggesting an ancient role for posttranslational histone modifications in transcription. Our findings show that total RNA is sufficient to identify transcribed regulatory elements and capture the dynamics of initiated stable and unstable transcripts at single-nucleotide resolution in eukaryotes.

SiCloneFit: Bayesian inference of population structure, genotype, and phylogeny of tumor clones from single-cell genome sequencing data [METHOD]

Zafar, H., Navin, N., Chen, K., Nakhleh, L. — 2019-11-01T08:06:42-07:00

Accumulation and selection of somatic mutations in a Darwinian framework result in intra-tumor heterogeneity (ITH) that poses significant challenges to the diagnosis and clinical therapy of cancer. Identification of the tumor cell populations (clones) and reconstruction of their evolutionary relationship can elucidate this heterogeneity. Recently developed single-cell DNA sequencing (SCS) technologies promise to resolve ITH to a single-cell level. However, technical errors in SCS data sets, including false-positives (FP) and false-negatives (FN) due to allelic dropout, and cell doublets, significantly complicate these tasks. Here, we propose a nonparametric Bayesian method that reconstructs the clonal populations as clusters of single cells, genotypes of each clone, and the evolutionary relationship between the clones. It employs a tree-structured Chinese restaurant process as the prior on the number and composition of clonal populations. The evolution of the clonal populations is modeled by a clonal phylogeny and a finite-site model of evolution to account for potential mutation recurrence and losses. We probabilistically account for FP and FN errors, and cell doublets are modeled by employing a Beta-binomial distribution. We develop a Gibbs sampling algorithm comprising partial reversible-jump and partial Metropolis-Hastings updates to explore the joint posterior space of all parameters. The performance of our method on synthetic and experimental data sets suggests that joint reconstruction of tumor clones and clonal phylogeny under a finite-site model of evolution leads to more accurate inferences. Our method is the first to enable this joint reconstruction in a fully Bayesian framework, thus providing measures of support of the inferences it makes.

PhISCS: a combinatorial approach for subperfect tumor phylogeny reconstruction via integrative use of single-cell and bulk sequencing data [METHOD]

2019-11-01T08:06:42-07:00

Available computational methods for tumor phylogeny inference via single-cell sequencing (SCS) data typically aim to identify the most likely perfect phylogeny tree satisfying the infinite sites assumption (ISA). However, the limitations of SCS technologies including frequent allele dropout and variable sequence coverage may prohibit a perfect phylogeny. In addition, ISA violations are commonly observed in tumor phylogenies due to the loss of heterozygosity, deletions, and convergent evolution. In order to address such limitations, we introduce the optimal subperfect phylogeny problem which asks to integrate SCS data with matching bulk sequencing data by minimizing a linear combination of potential false negatives (due to allele dropout or variance in sequence coverage), false positives (due to read errors) among mutation calls, and the number of mutations that violate ISA (real or because of incorrect copy number estimation). We then describe a combinatorial formulation to solve this problem which ensures that several lineage constraints imposed by the use of variant allele frequencies (VAFs, derived from bulk sequence data) are satisfied. We express our formulation both in the form of an integer linear program (ILP) and—as a first in tumor phylogeny reconstruction—a Boolean constraint satisfaction problem (CSP) and solve them by leveraging state-of-the-art ILP/CSP solvers. The resulting method, which we name PhISCS, is the first to integrate SCS and bulk sequencing data while accounting for ISA violating mutations. In contrast to the alternative methods, typically based on probabilistic approaches, PhISCS provides a guarantee of optimality in reported solutions. Using simulated and real data sets, we demonstrate that PhISCS is more general and accurate than all available approaches.

Quantitative mitochondrial DNA copy number determination using droplet digital PCR with single-cell resolution [METHOD]

O'Hara, R., Tedone, E., Ludlow, A., Huang, E., Arosio, B., Mari, D., Shay, J. W. — 2019-11-01T08:06:42-07:00

Mitochondria are involved in a number of diverse cellular functions, including energy production, metabolic regulation, apoptosis, calcium homeostasis, cell proliferation, and motility, as well as free radical generation. Mitochondrial DNA (mtDNA) is present at hundreds to thousands of copies per cell in a tissue-specific manner. mtDNA copy number also varies during aging and disease progression and therefore might be considered as a biomarker that mirrors alterations within the human body. Here, we present a new quantitative, highly sensitive droplet digital PCR (ddPCR) method, droplet digital mitochondrial DNA measurement (ddMDM), to measure mtDNA copy number not only from cell populations but also from single cells. Our developed assay can generate data in as little as 3 h, is optimized for 96-well plates, and also allows the direct use of cell lysates without the need for DNA purification or nuclear reference genes. We show that ddMDM is able to detect differences between samples whose mtDNA copy number was close enough as to be indistinguishable by other commonly used mtDNA quantitation methods. By utilizing ddMDM, we show quantitative changes in mtDNA content per cell across a wide variety of physiological contexts including cancer progression, cell cycle progression, human T cell activation, and human aging.

Single-pollen-cell sequencing for gamete-based phased diploid genome assembly in plants [METHOD]

2019-11-01T08:06:42-07:00

Genome assemblies from diploid organisms create mosaic sequences alternating between parental alleles, which can create erroneous gene models and other problems. In animals, a popular strategy to generate haploid genome-resolved assemblies has been the sampling of (haploid) gametes, and the advent of single-cell sequencing has further advanced such methods. However, several challenges for the isolation and amplification of DNA from plant gametes have limited such approaches in plants. Here, we combined a new approach for pollen protoplast isolation with a single-cell DNA amplification technique and then used a "barcoding" bioinformatics strategy to incorporate haploid-specific sequence data from 12 pollen cells, ultimately enabling the efficient and accurate phasing of the pear genome into its A and B haploid genomes. Beyond revealing that 8.12% of the genes in the pear reference genome feature mosaic assemblies and enabling a previously impossible analysis of allelic affects in pear gene expression, our new haploid genome assemblies provide high-resolution information about recombination during meiosis in pollen. Considering that outcrossing pear is an angiosperm species featuring very high heterozygosity, our method for rapidly phasing genome assemblies is potentially applicable to several yet-unsequenced outcrossing angiosperm species in nature.

Dynamics of microRNA expression during mouse prenatal development [RESOURCES]

2019-11-01T08:06:42-07:00

MicroRNAs (miRNAs) play a critical role as posttranscriptional regulators of gene expression. The ENCODE Project profiled the expression of miRNAs in an extensive set of organs during a time-course of mouse embryonic development and captured the expression dynamics of 785 miRNAs. We found distinct organ-specific and developmental stage–specific miRNA expression clusters, with an overall pattern of increasing organ-specific expression as embryonic development proceeds. Comparative analysis of conserved miRNAs in mouse and human revealed stronger clustering of expression patterns by organ type rather than by species. An analysis of messenger RNA expression clusters compared with miRNA expression clusters identifies the potential role of specific miRNA expression clusters in suppressing the expression of mRNAs specific to other developmental programs in the organ in which these miRNAs are expressed during embryonic development. Our results provide the most comprehensive time-course of miRNA expression as part of an integrated ENCODE reference data set for mouse embryonic development.

Genes essential for embryonic stem cells are associated with neurodevelopmental disorders [RESOURCES]

Shohat, S., Shifman, S. — 2019-11-01T08:06:42-07:00

Mouse embryonic stem cells (mESCs) are key components in generating mouse models for human diseases and performing basic research on pluripotency, yet the number of genes essential for mESCs is still unknown. We performed a genome-wide screen for essential genes in mESCs and compared it to screens in human cells. We found that essential genes are enriched for basic cellular functions, are highly expressed in mESCs, and tend to lack paralog genes. We discovered that genes that are essential specifically in mESCs play a role in pathways associated with their pluripotent state. We show that 29.5% of human genes intolerant to loss-of-function mutations are essential in mouse or human ESCs, and that the human phenotypes most significantly associated with genes essential for ESCs are neurodevelopmental. Our results provide insights into essential genes in the mouse, the pathways which govern pluripotency, and suggest that many genes associated with neurodevelopmental disorders are essential at very early embryonic stages.

A chromosome-level assembly of the Atlantic herring genome--detection of a supergene and other signals of selection [RESOURCES]

2019-11-01T08:06:42-07:00

The Atlantic herring is a model species for exploring the genetic basis for ecological adaptation, due to its huge population size and extremely low genetic differentiation at selectively neutral loci. However, such studies have so far been hampered because of a highly fragmented genome assembly. Here, we deliver a chromosome-level genome assembly based on a hybrid approach combining a de novo Pacific Biosciences (PacBio) assembly with Hi-C-supported scaffolding. The assembly comprises 26 autosomes with sizes ranging from 12.4 to 33.1 Mb and a total size, in chromosomes, of 726 Mb, which has been corroborated by a high-resolution linkage map. A comparison between the herring genome assembly with other high-quality assemblies from bony fishes revealed few inter-chromosomal but frequent intra-chromosomal rearrangements. The improved assembly facilitates analysis of previously intractable large-scale structural variation, allowing, for example, the detection of a 7.8-Mb inversion on Chromosome 12 underlying ecological adaptation. This supergene shows strong genetic differentiation between populations. The chromosome-based assembly also markedly improves the interpretation of previously detected signals of selection, allowing us to reveal hundreds of independent loci associated with ecological adaptation.