Claudia Kutter

Five-Vertebrate ChIP-seq Reveals the Evolutionary Dynamics of Transcription Factor Binding

Subtle Variation Despite vast phenotypic differences, vertebrates have many readily recognizable specific cell types, like liver hepatocytes. The gene expression that defines specific cells depends on evolutionarily conserved orthologous transcription factors. Schmidt et al. (p. 1036, published online 8 April) studied the conservation and divergence in the genome-wide binding of two such transcription factors, CEBPA and HNF4A, in livers from human, dog, mouse, short-tailed opossum, and chicken. Although the sequence bound by orthologous transcription factors was similar, the vast majority of binding events were unique to each species. Binding of two liver-specific transcription factors in several vertebrate species reveals complex regulatory evolution. Transcription factors (TFs) direct gene expression by binding to DNA regulatory regions. To explore the evolution of gene regulation, we used chromatin immunoprecipitation with high-throughput sequencing (ChIP-seq) to determine experimentally the genome-wide occupancy of two TFs, CCAAT/enhancer-binding protein alpha and hepatocyte nuclear factor 4 alpha, in the livers of five vertebrates. Although each TF displays highly conserved DNA binding preferences, most binding is species-specific, and aligned binding events present in all five species are rare. Regions near genes with expression levels that are dependent on a TF are often bound by the TF in multiple species yet show no enhanced DNA sequence constraint. Binding divergence between species can be largely explained by sequence changes to the bound motifs. Among the binding events lost in one lineage, only half are recovered by another binding event within 10 kilobases. Our results reveal large interspecies differences in transcriptional regulation and provide insight into regulatory evolution.

The Evolutionary Landscape of Alternative Splicing in Vertebrate Species

Whence Species Variation? Vertebrates have widely varying phenotypes that are at odds with their much more limited proteincoding genotypes and conserved messenger RNA expression patterns. Genes with multiple exons and introns can undergo alternative splicing, potentially resulting in multiple protein isoforms (see the Perspective by Papasaikas and Valcárcel). Barbosa-Morais et al. (p. 1587) and Merkin et al. (p. 1593) analyzed alternative splicing across the genomes of a variety of vertebrates, including human, primates, rodents, opossum, platypus, chicken, lizard, and frog. The findings suggest that the evolution of alternative splicing has for the most part been very rapid and that alternative splicing patterns of most organs more strongly reflect the identity of the species rather than the organ type. Species-classifying alternative splicing can affect key regulators, often in disordered regions of proteins that may influence protein-protein interactions, or in regions involved in protein phosphorylation. The patterns and complexity of messenger RNA splicing across vertebrates cluster by species rather than by organ. How species with similar repertoires of protein-coding genes differ so markedly at the phenotypic level is poorly understood. By comparing organ transcriptomes from vertebrate species spanning ~350 million years of evolution, we observed significant differences in alternative splicing complexity between vertebrate lineages, with the highest complexity in primates. Within 6 million years, the splicing profiles of physiologically equivalent organs diverged such that they are more strongly related to the identity of a species than they are to organ type. Most vertebrate species-specific splicing patterns are cis-directed. However, a subset of pronounced splicing changes are predicted to remodel protein interactions involving trans-acting regulators. These events likely further contributed to the diversification of splicing and other transcriptomic changes that underlie phenotypic differences among vertebrate species.

Waves of Retrotransposon Expansion Remodel Genome Organization and CTCF Binding in Multiple Mammalian Lineages

Summary CTCF-binding locations represent regulatory sequences that are highly constrained over the course of evolution. To gain insight into how these DNA elements are conserved and spread through the genome, we defined the full spectrum of CTCF-binding sites, including a 33/34-mer motif, and identified over five thousand highly conserved, robust, and tissue-independent CTCF-binding locations by comparing ChIP-seq data from six mammals. Our data indicate that activation of retroelements has produced species-specific expansions of CTCF binding in rodents, dogs, and opossum, which often functionally serve as chromatin and transcriptional insulators. We discovered fossilized repeat elements flanking deeply conserved CTCF-binding regions, indicating that similar retrotransposon expansions occurred hundreds of millions of years ago. Repeat-driven dispersal of CTCF binding is a fundamental, ancient, and still highly active mechanism of genome evolution in mammalian lineages. PaperClip

Molecular characterization of geminivirus-derived small RNAs in different plant species

DNA geminiviruses are thought to be targets of RNA silencing. Here, we characterize small interfering (si) RNAs—the hallmarks of silencing—associated with Cabbage leaf curl begomovirus in Arabidopsis and African cassava mosaic begomovirus in Nicotiana benthamiana and cassava. We detected 21, 22 and 24 nt siRNAs of both polarities, derived from both the coding and the intergenic regions of these geminiviruses. Genetic evidence showed that all the 24 nt and a substantial fraction of the 22 nt viral siRNAs are generated by the dicer-like proteins DCL3 and DCL2, respectively. The viral siRNAs were 5′ end phosphorylated, as shown by phosphatase treatments, and methylated at the 3′-nucleotide, as shown by HEN1 miRNA methylase-dependent resistance to β-elimination. Similar modifications were found in all types of endogenous and transgene-derived siRNAs tested, but not in a major fraction of siRNAs from a cytoplasmic RNA tobamovirus. We conclude that several distinct silencing pathways are involved in DNA virus-plant interactions.

Rapid turnover of long noncoding RNAs and the evolution of gene expression.

A large proportion of functional sequence within mammalian genomes falls outside protein-coding exons and can be transcribed into long RNAs. However, the roles in mammalian biology of long noncoding RNA (lncRNA) are not well understood. Few lncRNAs have experimentally determined roles, with some of these being lineage-specific. Determining the extent by which transcription of lncRNA loci is retained or lost across multiple evolutionary lineages is essential if we are to understand their contribution to mammalian biology and to lineage-specific traits. Here, we experimentally investigated the conservation of lncRNA expression among closely related rodent species, allowing the evolution of DNA sequence to be uncoupled from evolution of transcript expression. We generated total RNA (RNAseq) and H3K4me3-bound (ChIPseq) DNA data, and combined both to construct catalogues of transcripts expressed in the adult liver of Mus musculus domesticus (C57BL/6J), Mus musculus castaneus, and Rattus norvegicus. We estimated the rate of transcriptional turnover of lncRNAs and investigated the effects of their lineage-specific birth or death. LncRNA transcription showed considerably greater gain and loss during rodent evolution, compared with protein-coding genes. Nucleotide substitution rates were found to mirror the in vivo transcriptional conservation of intergenic lncRNAs between rodents: only the sequences of noncoding loci with conserved transcription were constrained. Finally, we found that lineage-specific intergenic lncRNAs appear to be associated with modestly elevated expression of genomically neighbouring protein-coding genes. Our findings show that nearly half of intergenic lncRNA loci have been gained or lost since the last common ancestor of mouse and rat, and they predict that such rapid transcriptional turnover contributes to the evolution of tissue- and lineage-specific gene expression.

MicroRNA-mediated regulation of stomatal development in Arabidopsis.

The proper number and distribution of stomata are essential for the efficient exchange of gases between the atmosphere and the aerial parts of plants. We show that the density and development of stomatal complexes on the epidermis of Arabidopsis thaliana leaves depend, in part, on the microRNA-mediated regulation of Agamous-like16 (AGL16), which is a member of the MADS box protein family. AGL16 mRNA is targeted for sequence-specific degradation by miR824, a recently evolved microRNA conserved in the Brassicaceae and encoded at a single genetic locus. Primary stomatal complexes can give rise to higher-order complexes derived from satellite meristemoids. Expression of a miR824-resistant AGL16 mRNA, but not the wild-type AGL16 mRNA, in transgenic plants increased the incidence of stomata in higher-order complexes. By contrast, reduced expression of AGL16 mRNA in the agl16-1 deficiency mutant and in transgenic lines overexpressing miR824 decreased the incidence of stomata in higher-order complexes. These findings and the nonoverlapping patterns of AGL16 mRNA and miR824 localization led us to propose that the miR824/AGL16 pathway functions in the satellite meristemoid lineage of stomatal development.

Aberrant methylation of tRNAs links cellular stress to neuro-developmental disorders

Mutations in the cytosine‐5 RNA methyltransferase NSun2 cause microcephaly and other neurological abnormalities in mice and human. How post‐transcriptional methylation contributes to the human disease is currently unknown. By comparing gene expression data with global cytosine‐5 RNA methylomes in patient fibroblasts and NSun2‐deficient mice, we find that loss of cytosine‐5 RNA methylation increases the angiogenin‐mediated endonucleolytic cleavage of transfer RNAs (tRNA) leading to an accumulation of 5′ tRNA‐derived small RNA fragments. Accumulation of 5′ tRNA fragments in the absence of NSun2 reduces protein translation rates and activates stress pathways leading to reduced cell size and increased apoptosis of cortical, hippocampal and striatal neurons. Mechanistically, we demonstrate that angiogenin binds with higher affinity to tRNAs lacking site‐specific NSun2‐mediated methylation and that the presence of 5′ tRNA fragments is sufficient and required to trigger cellular stress responses. Furthermore, the enhanced sensitivity of NSun2‐deficient brains to oxidative stress can be rescued through inhibition of angiogenin during embryogenesis. In conclusion, failure in NSun2‐mediated tRNA methylation contributes to human diseases via stress‐induced RNA cleavage.

miR393 and Secondary siRNAs Regulate Expression of the TIR1/AFB2 Auxin Receptor Clade and Auxin-Related Development of Arabidopsis Leaves

The phytohormone auxin is a key regulator of plant growth and development that exerts its functions through F-box receptors. Arabidopsis (Arabidopsis thaliana) has four partially redundant of these receptors that comprise the TRANSPORT INHIBITOR RESPONSE1/AUXIN SIGNALING F-BOX1 auxin receptor (TAAR) clade. Recent studies have shown that the microRNA miR393 regulates the expression of different sets of TAAR genes following pathogen infection or nitrate treatment. Here we report that miR393 helps regulate auxin-related development of leaves. We found that AtMIR393B is the predominant source for miR393 in all aerial organs and that miR393 down-regulates all four TAAR genes by guiding the cleavage of their mRNAs. A mutant unable to produce miR393 shows developmental abnormalities of leaves and cotyledons reminiscent of enhanced auxin perception by TAARs. Interestingly, miR393 initiates the biogenesis of secondary siRNAs from the transcripts of at least two of the four TAAR genes. Our results indicate that these siRNAs, which we call siTAARs, help regulate the expression of TAAR genes as well as several unrelated genes by guiding the cleavage of their mRNAs. Thus, miR393 and possibly siTAARs regulate auxin perception and certain auxin-related aspects of leaf development.

Epigenetic conservation at gene regulatory elements revealed by non-methylated DNA profiling in seven vertebrates

Two-thirds of gene promoters in mammals are associated with regions of non-methylated DNA, called CpG islands (CGIs), which counteract the repressive effects of DNA methylation on chromatin. In cold-blooded vertebrates, computational CGI predictions often reside away from gene promoters, suggesting a major divergence in gene promoter architecture across vertebrates. By experimentally identifying non-methylated DNA in the genomes of seven diverse vertebrates, we instead reveal that non-methylated islands (NMIs) of DNA are a central feature of vertebrate gene promoters. Furthermore, NMIs are present at orthologous genes across vast evolutionary distances, revealing a surprising level of conservation in this epigenetic feature. By profiling NMIs in different tissues and developmental stages we uncover a unifying set of features that are central to the function of NMIs in vertebrates. Together these findings demonstrate an ancient logic for NMI usage at gene promoters and reveal an unprecedented level of epigenetic conservation across vertebrate evolution. DOI: http://dx.doi.org/10.7554/eLife.00348.001

Role of cis-12-Oxo-Phytodienoic Acid in Tomato Embryo Development

Oxylipins including jasmonates are signaling compounds in plant growth, development, and responses to biotic and abiotic stresses. In Arabidopsis (Arabidopsis thaliana) most mutants affected in jasmonic acid (JA) biosynthesis and signaling are male sterile, whereas the JA-insensitive tomato (Solanum lycopersicum) mutant jai1 is female sterile. The diminished seed formation in jai1 together with the ovule-specific accumulation of the JA biosynthesis enzyme allene oxide cyclase (AOC), which correlates with elevated levels of JAs, suggest a role of oxylipins in tomato flower/seed development. Here, we show that 35S::SlAOC-RNAi lines with strongly reduced AOC in ovules exhibited reduced seed set similarly to the jai1 plants. Investigation of embryo development of wild-type tomato plants showed preferential occurrence of AOC promoter activity and AOC protein accumulation in the developing seed coat and the embryo, whereas 12-oxo-phytodienoic acid (OPDA) was the dominant oxylipin occurring nearly exclusively in the seed coat tissues. The OPDA- and JA-deficient mutant spr2 was delayed in embryo development and showed an increased programmed cell death in the developing seed coat and endosperm. In contrast, the mutant acx1a, which accumulates preferentially OPDA and residual amount of JA, developed embryos similar to the wild type, suggesting a role of OPDA in embryo development. Activity of the residual amount of JA in the acx1a mutant is highly improbable since the known reproductive phenotype of the JA-insensitive mutant jai1 could be rescued by wound-induced formation of OPDA. These data suggest a role of OPDA or an OPDA-related compound for proper embryo development possibly by regulating carbohydrate supply and detoxification.