Small RNA profiling of Xenopus embryos reveals novel miRNAs and a new class of small RNAs derived from intronic transposable elements

Revised Annotation and Characterization of Novel Aedes albopictus miRNAs and Their Potential Functions in Dengue Virus Infection

Simple Summary

Aedes albopictus (Ae. albopictus) is an important vector of the dengue virus. Genetics and molecular studies of virus infection in mosquito vectors are important to uncover the basic biology of the virus. It has been reported that miRNAs are important and possess functional roles in virus infection in Ae. albopictus. Here, we report a comprehensive catalog of miRNAs using the latest genome version of Ae. albopictus. We discovered a total of 72 novel mature miRNAs, 44 of which were differentially expressed in C6/36 cells infected with the dengue virus. Target prediction analysis revealed that the differentially expressed miRNAs were involved in lipid metabolism and protein processing in the endoplasmic reticulum. Results from this study provide a valuable resource for researchers to study miRNAs in this mosquito vector, especially in host–virus interactions.

Abstract

The Asian tiger mosquito, Ae. albopictus, is a highly invasive species that transmits several arboviruses including dengue (DENV), Zika (ZIKV), and chikungunya (CHIKV). Although several studies have identified microRNAs (miRNAs) in Ae. albopictus, it is crucial to extend and improve current annotations with both the newly improved genome assembly and the increased number of small RNA-sequencing data. We combined our high-depth sequence data and 26 public datasets to re-annotate Ae. albopictus miRNAs and found a total of 72 novel mature miRNAs. We discovered that the expression of novel miRNAs was lower than known miRNAs. Furthermore, compared to known miRNAs, novel miRNAs are prone to expression in a stage-specific manner. Upon DENV infection, a total of 44 novel miRNAs were differentially expressed, and target prediction analysis revealed that miRNA-target genes were involved in lipid metabolism and protein processing in endoplasmic reticulum. Taken together, the miRNA annotation profile provided here is the most comprehensive to date. We believed that this would facilitate future research in understanding virus–host interactions, particularly in the role of miRNAs.

Computational meta-analysis of ribosomal RNA fragments: potential targets and interaction mechanisms

The most abundant cellular RNA species, ribosomal RNA (rRNA), appears to be a source of massive amounts of non-randomly generated fragments. We found rRNA fragments (rRFs) in immunoprecipitated Argonaute (Ago-IP) complexes in human and mouse cells and in small RNA sequencing datasets. In human Ago1-IP, guanine-rich rRFs were preferentially cut in single-stranded regions of mature rRNAs between pyrimidines and adenosine, and non-randomly paired with cellular transcripts in crosslinked chimeras. Numerous identical rRFs were found in the cytoplasm and nucleus in mouse Ago2-IP. We report specific interaction motifs enriched in rRF-target pairs. Locations of such motifs on rRFs were compatible with the Ago structural features and patterns of the Ago-RNA crosslinking in both species. Strikingly, many of these motifs may bind to double-stranded regions on target RNAs, suggesting a potential pathway for regulating translation by unwinding mRNAs. Occurring on either end of rRFs and matching intronic, untranslated or coding regions in targets, such interaction sites extend the concept of microRNA seed regions. Targeting both borders of certain short introns, rRFs may be involved in their biogenesis or function, facilitated by Ago. Frequently dismissed as noise, rRFs are poised to greatly enrich the known functional spectrum of small RNA regulation.

Transposable element-derived sequences in vertebrate development

Transposable elements (TEs) are major components of all vertebrate genomes that can cause deleterious insertions and genomic instability. However, depending on the specific genomic context of their insertion site, TE sequences can sometimes get positively selected, leading to what are called "exaptation" events. TE sequence exaptation constitutes an important source of novelties for gene, genome and organism evolution, giving rise to new regulatory sequences, protein-coding exons/genes and non-coding RNAs, which can play various roles beneficial to the host. In this review, we focus on the development of vertebrates, which present many derived traits such as bones, adaptive immunity and a complex brain. We illustrate how TE-derived sequences have given rise to developmental innovations in vertebrates and how they thereby contributed to the evolutionary success of this lineage.

Discovery of MicroRNAs from Batrachuperus yenyuanensis Using Deep Sequencing and Prediction of Their Targets

MicroRNAs (miRNAs), a family of ∼22-nucleotide non-coding single-stranded RNA molecules, are considered as key post-transcriptional regulators of gene expression that regulate various biological processes in living organism. Many miRNAs have been identified in animals; however, few have been reported in Hynobiidae species. The present study is aimed to identify a full repertoire of miRNAs in Batrachuperus yenyuanensis (Yenyuan stream salamander), which would significantly increase our knowledge of miRNAs in amphibians. A small RNA library was constructed from B. yenyuanensis and sequenced using deep sequencing. As a result, 1,717,751 clean reads were obtained, representing 356 known and 80 novel miRNAs. Additionally, expression levels of eight randomly selected miRNAs in B. yenyuanensis were confirmed using the stem-loop quantitative real-time reverse transcription PCR. In addition, 13,972 targets were predicted for these identified miRNAs, although the physiological functions of many of these targets remain unknown. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis suggested that the predicted targets are involved in a variety of physiological regulatory functions in B. yenyuanensis. These results provide useful information for further research on the miRNAs involved in the growth and development of B. yenyuanensis, as well as adaptation of this species to its high-altitude habitats.

Identification and differential regulation of microRNAs during thyroid hormone-dependent metamorphosis in Microhyla fissipes

BACKGROUND

Anuran metamorphosis, which is obligatorily initiated and sustained by thyroid hormone (TH), is a dramatic example of extensive morphological, biochemical and cellular changes occurring during post-embryonic development. Thus, it provides an ideal model to understand the actions of the hormone and molecular mechanisms underlying these developmental and apoptotic processes. In addition to transcriptional factors, microRNAs (miRNAs) play key roles in diverse biological processes via post-transcriptional repression of mRNAs. However, the possible role of miRNAs in anuran metamorphosis is not well understood. Screening and identification of TH-responding miRNAs are required to reveal the integrated regulatory mechanisms of TH during metamorphosis. Given the specific role of TRs during M. fissipes metamorphosis and the characteristics of M. fissipes as an ideal model, Illumina sequencing technology was employed to get a full scope of miRNA in M. fissipes metamorphosis treated by T3.

RESULTS

Morphological and histological analysis revealed that 24 h T3 treatment M. fissipes tadpoles resembled that at the climax of natural metamorphosis. Thus, small RNA libraries were constructed from control and 24 h T3 treatment groups. A total of 164 conserved miRNAs and 36 predicted novel miRNAs were characterized. Furthermore, 5' first and ninth nucleotides of miRNAs were significantly enriched in U in our study. In all, 21 miRNAs were differentially expressed between the T3 and control groups (p < 0.01). A total of 10,206 unigenes were identified as target genes of these differentially expressed miRNAs. KEGG pathway analysis indicated that the most overrepresented miRNA target genes were enriched in the "PI3k-Akt signaling pathway". In addition, a network associated with the TH signaling pathway provides an opportunity to further understand the complex biological processes that occur in metamorphosis.

CONCLUSIONS

We identified a large number of miRNAs during M. fissipes metamorphosis, and 21 of them were differentially expressed in the two groups that represented two different metamorphic stages. These miRNAs may play important roles during metamorphosis. The study gives us clues for further studies of the mechanisms of anuran metamorphosis and provides a model to study the mechanism of TH-affected biological processes in humans.

microRNAs associated with early neural crest development in Xenopus laevis

BACKGROUND

The neural crest (NC) is a class of transitory stem cell-like cells unique to vertebrate embryos. NC cells arise within the dorsal neural tube where they undergo an epithelial to mesenchymal transition in order to migrate and differentiate throughout the developing embryo. The derivative cell types give rise to multiple tissues, including the craniofacial skeleton, peripheral nervous system and skin pigment cells. Several well-studied gene regulatory networks underpin NC development, which when disrupted can lead to various neurocristopathies such as craniofrontonasal dysplasia, DiGeorge syndrome and some forms of cancer. Small RNAs, such as microRNAs (miRNAs) are non-coding RNA molecules important in post-transcriptional gene silencing and critical for cellular regulation of gene expression.

RESULTS

To uncover novel small RNAs in NC development we used high definition adapters and next generation sequencing of libraries derived from ectodermal explants of Xenopus laevis embryos induced to form neural and NC tissue. Ectodermal and blastula animal pole (blastula) stage tissues were also sequenced. We show that miR-427 is highly abundant in all four tissue types though in an isoform specific manner and we define a set of 11 miRNAs that are enriched in the NC. In addition, we show miR-301a and miR-338 are highly expressed in both the NC and blastula suggesting a role for these miRNAs in maintaining the stem cell-like phenotype of NC cells.

CONCLUSION

We have characterised the miRNAs expressed in Xenopus embryonic explants treated to form ectoderm, neural or NC tissue. This has identified novel tissue specific miRNAs and highlighted differential expression of miR-427 isoforms.

Transcriptomic analysis of maternally provisioned cues for phenotypic plasticity in the annual killifish, Austrofundulus limnaeus

Background

Genotype and environment can interact during development to produce novel adaptive traits that support life in extreme conditions. The development of the annual killifish Austrofundulus limnaeus is unique among vertebrates because the embryos have distinct cell movements that separate epiboly from axis formation during early development, can enter into a state of metabolic dormancy known as diapause and can survive extreme environmental conditions. The ability to enter into diapause can be maternally programmed, with young females producing embryos that do not enter into diapause. Alternately, embryos can be programmed to “escape” from diapause and develop directly by both maternal factors and embryonic incubation conditions. Thus, maternally packaged gene products are hypothesized to regulate developmental trajectory and perhaps the other unique developmental characters in this species.

Results

Using high-throughput RNA sequencing, we generated transcriptomic profiles of mRNAs, long non-coding RNAs and small non-coding RNAs (sncRNAs) in 1–2 cell stage embryos of A. limnaeus. Transcriptomic analyses suggest maternal programming of embryos through alternatively spliced mRNAs and antisense sncRNAs. Comparison of these results to those of comparable studies on zebrafish and other fishes reveals a surprisingly high abundance of transcripts involved in the cellular response to stress and a relatively lower expression of genes required for rapid transition through the cell cycle.

Conclusions

Maternal programming of developmental trajectory is unlikely accomplished by differential expression of diapause-specific genes. Rather, evidence suggests a role for trajectory-specific splice variants of genes expressed in both phenotypes. In addition, based on comparative studies with zebrafish, the A. limnaeus 1–2 cell stage transcriptome is unique in ways that are consistent with their unique life history. These results not only impact our understanding of the genetic mechanisms that regulate entrance into diapause, but also provide insight into the epigenetic regulation of gene expression during development.

Electronic supplementary material

The online version of this article (doi:10.1186/s13227-017-0069-7) contains supplementary material, which is available to authorized users.

A gene regulatory program controlling early Xenopus mesendoderm formation: Network conservation and motifs

Germ layer formation is among the earliest differentiation events in metazoan embryos. In triploblasts, three germ layers are formed, among which the endoderm gives rise to the epithelial lining of the gut tube and associated organs including the liver, pancreas and lungs. In frogs (Xenopus), where early germ layer formation has been studied extensively, the process of endoderm specification involves the interplay of dozens of transcription factors. Here, we review the interactions between these factors, summarized in a transcriptional gene regulatory network (GRN). We highlight regulatory connections conserved between frog, fish, mouse, and human endodermal lineages. Especially prominent is the conserved role and regulatory targets of the Nodal signaling pathway and the T-box transcription factors, Vegt and Eomes. Additionally, we highlight network topologies and motifs, and speculate on their possible roles in development.

Mechanisms of Vertebrate Germ Cell Determination

Two unique characteristics of the germ line are the ability to persist from generation to generation and to retain full developmental potential while differentiating into gametes. How the germ line is specified that allows it to retain these characteristics within the context of a developing embryo remains unknown and is one focus of current research. Germ cell specification proceeds through one of two basic mechanisms: cell autonomous or inductive. Here, we discuss how germ plasm driven germ cell specification (cell autonomous) occurs in both zebrafish and the frog Xenopus. We describe the segregation of germ cells during embryonic development of solitary and colonial ascidians to provide an evolutionary context to both mechanisms. We conclude with a discussion of the inductive mechanism as exemplified by both the mouse and axolotl model systems. Regardless of mechanism, several general themes can be recognized including the essential role of repression and posttranscriptional regulation of gene expression.

piRNAs and their diverse roles: a transposable element-driven tactic for gene regulation?

P-element-induced wimpy testis (PIWI)-interacting RNAs (piRNAs) are small, noncoding RNAs known for silencing transposable elements (TEs) in the germline of animals. Most genomes host TEs, which are notorious for mobilizing themselves and endangering survival of the host if not controlled. By silencing TEs in the germline, piRNAs prevent harmful mutations from being passed on to the next generation. How piRNAs are generated and how they silence TEs were the focus of researchers ever since their discovery. Now a spate of recent papers are beginning to tell us that piRNAs can play roles beyond TE silencing and are involved in diverse cellular processes from mRNA regulation to development or genome rearrangement. In this review, we discuss some of these recently reported roles. Data on these new roles are often rudimentary, and the involvement of piRNAs in these processes is yet to be definitely established. What is interesting is that the reports are on animals widely separated on the phylogenetic tree of life and that piRNAs were also found outside the gonadal tissues. Some of these piRNAs map to TE sequences, prompting us to hypothesize that genomes may have co-opted the TE-derived piRNA system for their own regulation.-Sarkar, A., Volff, J.-N., Vaury, C. piRNAs and their diverse roles: a transposable element-driven tactic for gene regulation?

Heterochromatic histone modifications at transposons in Xenopus tropicalis embryos

Transposable elements are parasitic genomic elements that can be deleterious for host gene function and genome integrity. Heterochromatic histone modifications are involved in the repression of transposons. However, it remains unknown how these histone modifications mark different types of transposons during embryonic development. Here we document the variety of heterochromatic epigenetic signatures at parasitic elements during development in Xenopus tropicalis, using genome-wide ChIP-sequencing data and ChIP-qPCR analysis. We show that specific subsets of transposons in various families and subfamilies are marked by different combinations of the heterochromatic histone modifications H4K20me3, H3K9me2/3 and H3K27me3. Many DNA transposons are marked at the blastula stage already, whereas at retrotransposons the histone modifications generally accumulate at the gastrula stage or later. Furthermore, transposons marked by H3K9me3 and H4K20me3 are more prominent in gene deserts. Using intra-subfamily divergence as a proxy for age, we show that relatively young DNA transposons are preferentially marked by early embryonic H4K20me3 and H3K27me3. In contrast, relatively young retrotransposons are marked by increasing H3K9me3 and H4K20me3 during development, and are also linked to piRNA-sized small non-coding RNAs. Our results implicate distinct repression mechanisms that operate in a transposon-selective and developmental stage-specific fashion.

The Influence of Genotype and Environment on Small RNA Profiles in Grapevine Berry

Understanding the molecular mechanisms involved in the interaction between the genetic composition and the environment is crucial for modern viticulture. We approached this issue by focusing on the small RNA transcriptome in grapevine berries of the two varieties Cabernet Sauvignon and Sangiovese, growing in adjacent vineyards in three different environments. Four different developmental stages were studied and a total of 48 libraries of small RNAs were produced and sequenced. Using a proximity-based pipeline, we determined the general landscape of small RNAs accumulation in grapevine berries. We also investigated the presence of known and novel miRNAs and analyzed their accumulation profile. The results showed that the distribution of small RNA-producing loci is variable between the two cultivars, and that the level of variation depends on the vineyard. Differently, the profile of miRNA accumulation mainly depends on the developmental stage. The vineyard in Riccione maximizes the differences between the varieties, promoting the production of more than 1000 specific small RNA loci and modulating their expression depending on the cultivar and the maturation stage. In total, 89 known vvi-miRNAs and 33 novel vvi-miRNA candidates were identified in our samples, many of them showing the accumulation profile modulated by at least one of the factors studied. The in silico prediction of miRNA targets suggests their involvement in berry development and in secondary metabolites accumulation such as anthocyanins and polyphenols.

Identification of microRNAs and microRNA targets in Xenopus gastrulae: The role of miR-26 in the regulation of Smad1

MicroRNAs (miRNAs) are known to play diverse roles in the regulation of vertebrate development. To investigate miRNA-target mRNA relationships in embryonic development, we have carried out small-RNA sequencing to identify miRNAs expressed in the early gastrula of Xenopus laevis. We identify a total of 180 miRNAs, and we have identified the locations of the miRNA precursor sequences in the X. laevis genome. Of these miRNAs, 141 represent miRs previously identified in Xenopus tropicalis. Alignment to human miRNAs led to the identification of 39 miRNAs that have not previously been described for Xenopus. We have also used a biochemical approach to isolate mRNAs that are associated with the RNA-Induced Silencing Complex (RISC) in early gastrulae and thus candidate targets of miRNA-dependent regulation. Interrogation of this RISC-associated mRNA pool by RT-PCR indicates that a number of genes essential for early patterning and specification may be under regulation by miRNAs. Smad1 transcripts are associated with the RISC; target prediction algorithms identify a single miRNA-binding site for miR-26, which is common to the 3'UTRs of Smad1a and Smad1b. Disruption of the interaction between miR-26 and the Smad1 3'UTR via a Target Protector Morpholino Oligonucleotide (TPMO) leads to a 2-fold increase in Smad1 protein accumulation, moderate increases in the expression of BMP4/Smad1 target genes, and a reduction in organizer gene expression, as well as a partially ventralized phenotype in approximately 25% of embryos. Overexpression of miR-26 resulted in moderately decreased expression of Smad1-dependent genes and an expansion of the region expressing the Organizer gene not1. Our findings indicate that interactions between miR-26 and the Smad1 3'UTR modulate Smad1 function in the establishment of axial patterning; they also establish a foundation for the functional analysis of miRNAs and their regulatory interactions during gastrulation.

Detecting endogenous retrovirus-driven tissue-specific gene transcription

Transposable elements (TEs) comprise approximately half of the human genome, and several independent lines of investigation have demonstrated their role in rewiring gene expression during development, evolution, and oncogenesis. The identification of their regulatory effects has largely been idiosyncratic, by linking activity with isolated genes. Their distribution throughout the genome raises critical questions—do these elements contribute to broad tissue- and lineage-specific regulation? If so, in what manner, as enhancers, promoters, RNAs? Here, we devise a novel approach to systematically dissect the genome-wide consequences of TE insertion on gene expression, and test the hypothesis that classes of endogenous retrovirus long terminal repeats (LTRs) exert tissue-specific regulation of adjacent genes. Using correlation of expression patterns across 18 tissue types, we reveal the tissue-specific uncoupling of gene expression due to 62 different LTR classes. These patterns are specific to the retroviral insertion, as the same genes in species without the LTRs do not exhibit the same effect. Although the LTRs can be transcribed themselves, the most highly transcribed TEs do not have the largest effects on adjacent regulation of coding genes, suggesting they function predominantly as enhancers. Moreover, the tissue-specific patterns of gene expression that are detected by our method arise from a limited number of genes, rather than as a general consequence of LTR integration. These findings identify basic principles of co-opting LTRs for genome evolution, and support the utility of our method for the analysis of TE, or other specific gene sets, in relation to the rest of the genome.

Adenylation of maternally inherited microRNAs by Wispy

Early development depends heavily on accurate control of maternally inherited mRNAs, and yet it remains unknown how maternal microRNAs are regulated during maternal-to-zygotic transition (MZT). We here find that maternal microRNAs are highly adenylated at their 3' ends in mature oocytes and early embryos. Maternal microRNA adenylation is widely conserved in fly, sea urchin, and mouse. We identify Wispy, a noncanonical poly(A) polymerase, as the enzyme responsible for microRNA adenylation in flies. Knockout of wispy abrogates adenylation and results in microRNA accumulation in eggs, whereas overexpression of Wispy increases adenylation and reduces microRNA levels in S2 cells. Wispy interacts with Ago1 through protein-protein interaction, which may allow the effective and selective adenylation of microRNAs. Thus, adenylation may contribute to the clearance of maternally deposited microRNAs during MZT. Our work provides mechanistic insights into the regulation of maternal microRNAs and illustrates the importance of RNA tailing in development.