Search

The National Genomics Data Center (NGDC), part of the China National Center for Bioinformation (CNCB), provides a family of database resources to support global academic and industrial communities. With the explosive accumulation of multi-omics data generated at an unprecedented rate, CNCB-NGDC constantly expands and updates core database resources by big data archive, integrative analysis and value-added curation. In the past year, efforts have been devoted to integrating multiple omics data, synthesizing the growing knowledge, developing new resources and upgrading a set of major resources. Particularly, several database resources are newly developed for infectious diseases and microbiology (MPoxVR, KGCoV, ProPan), cancer-trait association (ASCancer Atlas, TWAS Atlas, Brain Catalog, CCAS) as well as tropical plants (TCOD). Importantly, given the global health threat caused by monkeypox virus and SARS-CoV-2, CNCB-NGDC has newly constructed the monkeypox virus resource, along with frequent updates of SARS-CoV-2 genome sequences, variants as well as haplotypes. All the resources and services are publicly accessible at https://ngdc.cncb.ac.cn.

Piwi-interacting RNAs are a type of small noncoding RNA that have various functions. piRBase is a manually curated resource focused on assisting piRNA functional analysis. piRBase release v3.0 is committed to providing more comprehensive piRNA related information. The latest release covers >181 million unique piRNA sequences, including 440 datasets from 44 species. More disease-related piRNAs and piRNA targets have been collected and displayed. The regulatory relationships between piRNAs and targets have been visualized. In addition to the reuse and expansion of the content in the previous version, the latest version has additional new content, including gold standard piRNA sets, piRNA clusters, piRNA variants, splicing-junction piRNAs, and piRNA expression data. In addition, the entire web interface has been redesigned to provide a better experience for users. piRBase release v3.0 is free to access, browse, search, and download at http://bigdata.ibp.ac.cn/piRBase.

Figs and fig pollinators are one of the few classic textbook examples of obligate pollination mutualism. The specific dependence of fig pollinators on the relatively safe living environment with sufficient food sources in the enclosed fig syconia implies that they are vulnerable to habitat changes. However, there is still no extensive genomic evidence to reveal the evolutionary footprint of this long-term mutually beneficial symbiosis in fig pollinators. In fig syconia, there are also non-pollinator species. The non-pollinator species differ in their evolutionary and life histories from pollinators. We conducted comparative analyses on 11 newly sequenced fig wasp genomes and one previously published genome. The pollinators colonized the figs approximately 66.9 million years ago, consistent with the origin of host figs. Compared with non-pollinators, many more genes in pollinators were subject to relaxed selection. Seven genes were absent in pollinators in response to environmental stress and immune activation. Pollinators had more streamlined gene repertoires in the innate immune system, chemosensory toolbox, and detoxification system. Our results provide genomic evidence for the differentiation between pollinators and nonpollinators. The data suggest that owing to the long-term adaptation to the fig, some genes related to functions no longer required are absent in pollinators.

The National Genomics Data Center (NGDC), part of the China National Center for Bioinformation (CNCB), provides a suite of database resources to support worldwide research activities in both academia and industry. With the explosive growth of multi-omics data, CNCB-NGDC is continually expanding, updating and enriching its core database resources through big data deposition, integration and translation. In the past year, considerable efforts have been devoted to 2019nCoVR, a newly established resource providing a global landscape of SARS-CoV-2 genomic sequences, variants, and haplotypes, as well as Aging Atlas, BrainBase, GTDB (Glycosyltransferases Database), LncExpDB, and TransCirc (Translation potential for circular RNAs). Meanwhile, a series of resources have been updated and improved, including BioProject, BioSample, GWH (Genome Warehouse), GVM (Genome Variation Map), GEN (Gene Expression Nebulas) as well as several biodiversity and plant resources. Particularly, BIG Search, a scalable, one-stop, cross-database search engine, has been significantly updated by providing easy access to a large number of internal and external biological resources from CNCB-NGDC, our partners, EBI and NCBI. All of these resources along with their services are publicly accessible at https://bigd.big.ac.cn.

NONCODE (http://www.noncode.org/) is a comprehensive database of collection and annotation of noncoding RNAs, especially long non-coding RNAs (lncRNAs) in animals. NONCODEV6 is dedicated to providing the full scope of lncRNAs across plants and animals. The number of lncRNAs in NONCODEV6 has increased from 548 640 to 644 510 since the last update in 2017. The number of human lncRNAs has increased from 172 216 to 173 112. The number of mouse lncRNAs increased from 131 697 to 131 974. The number of plant lncRNAs is 94 697. The relationship between lncRNAs in human and cancer were updated with transcriptome sequencing profiles. Three important new features were also introduced in NONCODEV6: (i) updated human lncRNA-disease relationships, especially cancer; (ii) lncRNA annotations with tissue expression profiles and predicted function in five common plants; iii) lncRNAs conservation annotation at transcript level for 23 plant species. NONCODEV6 is accessible through http://www.noncode.org/.

Retrotransposons are populated in vertebrate genomes, and when active, are thought to cause genome instability with potential benefit to genome evolution. Retrotransposon-derived RNAs are also known to give rise to small endo-siRNAs to help maintain heterochromatin at their sites of transcription; however, as not all heterochromatic regions are equally active in transcription, it remains unclear how heterochromatin is maintained across the genome. Here, we address these problems by defining the origins of repeat-derived RNAs and their specific chromatin locations in Drosophila S2 cells. We demonstrate that repeat RNAs are predominantly derived from active gypsy elements and processed by Dcr-2 into small RNAs to help maintain pericentromeric heterochromatin. We also show in cultured S2 cells that synthetic repeat-derived endo-siRNA mimics are sufficient to rescue Dcr-2-deficiency-induced defects in heterochromatin formation in interphase and chromosome segregation during mitosis, demonstrating that active retrotransposons are required for stable genetic inheritance.

Highly structured RNA molecules usually interact with each other, and associate with various RNA-binding proteins, to regulate critical biological processes. However, RNA structures and interactions in intact cells remain largely unknown. Here, by coupling proximity ligation mediated by RNA-binding proteins with deep sequencing, we report an RNA in situ conformation sequencing (RIC-seq) technology for the global profiling of intra- and intermolecular RNA–RNA interactions. This technique not only recapitulates known RNA secondary structures and tertiary interactions, but also facilitates the generation of three-dimensional (3D) interaction maps of RNA in human cells. Using these maps, we identify noncoding RNA targets globally, and discern RNA topological domains and trans-interacting hubs. We reveal that the functional connectivity of enhancers and promoters can be assigned using their pairwise-interacting RNAs. Furthermore, we show that CCAT1-5L—a super-enhancer hub RNA—interacts with the RNA-binding protein hnRNPK, as well as RNA derived from the MYC promoter and enhancer, to boost MYC transcription by modulating chromatin looping. Our study demonstrates the power and applicability of RIC-seq in discovering the 3D structures, interactions and regulatory roles of RNA. RNA in situ conformation sequencing (RIC-seq) enables the generation of three-dimensional interaction maps of RNA in cells, which sheds light on the interactions and regulatory functions of RNA.

Mobile element insertions (MEIs) are a major class of structural variants (SVs) and have been linked to many human genetic disorders, including hemophilia, neurofibromatosis, and various cancers. However, human MEI resources from large-scale genome sequencing are still lacking compared to those for SNPs and SVs. Here, we report a comprehensive map of 36 699 non-reference MEIs constructed from 5675 genomes, comprising 2998 Chinese samples (∼26.2×, NyuWa) and 2677 samples from the 1000 Genomes Project (∼7.4×, 1KGP). We discovered that LINE-1 insertions were highly enriched in centromere regions, implying the role of chromosome context in retroelement insertion. After functional annotation, we estimated that MEIs are responsible for about 9.3% of all protein-truncating events per genome. Finally, we built a companion database named HMEID for public use. This resource represents the latest and largest genomewide study on MEIs and will have broad utility for exploration of human MEI findings.

Chemotherapeutic agents, such as 5-fluorouracil (5-FU) and oxaliplatin (Oxi), can not only kill the cancer cell but also influence the proliferation of gut microbiota; however, the interaction between these drugs and gut microbiota remains poorly understood. In this study, we developed a powerful framework for taxonomy composition and genomic variation analysis to investigate the mutagenesis effect and proliferation influence of chemotherapeutic agents, such as 5-FU and Oxi, on gut microbiota and the interaction between these drugs and gut microbiota during chemotherapy. Using the gut microbiome data, we detected 1.45 million variations among the chemotherapy groups and found the drugs significantly affected mutation signatures of gut microbiota. Oxi notably increased transversion rate, whereas 5-FU reduced the rate. Traits related to cell division and nutrient mobilization showed evidence of strong selection pressure from chemotherapeutic agents. In addition, drug-associated bacteriome shift patterns and functional alterations were found: the metabolism changes in the 5-FU group implied that gut microbiota could provide additional nicotinamide adenine dinucleotide (NAD+) to inhibit cancer cell autophagy; in the Oxi group, the ribosome and lysine biosynthesis genes were obviously enriched. Our study provides a blueprint for characterizing the role of microbes and drug–microbe interaction in the gut microbiota response to chemotherapy.

The National Genomics Data Center (NGDC), part of the China National Center for Bioinformation (CNCB), provides a family of database resources to support global research in both academia and industry. With the explosively accumulated multi-omics data at ever-faster rates, CNCB-NGDC is constantly scaling up and updating its core database resources through big data archive, curation, integration and analysis. In the past year, efforts have been made to synthesize the growing data and knowledge, particularly in single-cell omics and precision medicine research, and a series of resources have been newly developed, updated and enhanced. Moreover, CNCB-NGDC has continued to daily update SARS-CoV-2 genome sequences, variants, haplotypes and literature. Particularly, OpenLB, an open library of bioscience, has been established by providing easy and open access to a substantial number of abstract texts from PubMed, bioRxiv and medRxiv. In addition, Database Commons is significantly updated by cataloguing a full list of global databases, and BLAST tools are newly deployed to provide online sequence search services. All these resources along with their services are publicly accessible at https://ngdc.cncb.ac.cn.

Small proteins specifically refer to proteins consisting of less than 100 amino acids translated from small open reading frames (sORFs), which were usually missed in previous genome annotation. The significance of small proteins has been revealed in current years, along with the discovery of their diverse functions. However, systematic annotation of small proteins is still insufficient. SmProt was specially developed to provide valuable information on small proteins for scientific community. Here we present the update of SmProt, which emphasizes reliability of translated sORFs, genetic variants in translated sORFs, disease-specific sORF translation events or sequences, and remarkably increased data volume. More components such as non-ATG translation initiation, function, and new sources are also included. SmProt incorporated 638,958 unique small proteins curated from 3,165,229 primary records, which were computationally predicted from 419 ribosome profiling (Ribo-seq) datasets or collected from literature and other sources from 370 cell lines or tissues in 8 species (Homo sapiens, Mus musculus, Rattus norvegicus, Drosophila melanogaster, Danio rerio, Saccharomyces cerevisiae, Caenorhabditis elegans, and Escherichia coli). In addition, small protein families identified from human microbiomes were also collected. All datasets in SmProt are free to access, and available for browse, search, and bulk downloads at http://bigdata.ibp.ac.cn/SmProt/ .

RNA-binding proteins (RBPs) have essential functions during germline and early embryo development. However, current methods are unable to identify the in vivo targets of a RBP in these low-abundance cells. Here, by coupling RBP-mediated reverse transcription termination with linear amplification of complementary DNA ends and sequencing, we present the LACE-seq method for identifying RBP-regulated RNA networks at or near the single-oocyte level. We determined the binding sites and regulatory mechanisms for several RBPs, including Argonaute 2 (Ago2), Mili, Ddx4 and Ptbp1, in mature mouse oocytes. Unexpectedly, transcriptomics and proteomics analysis of Ago2-/- oocytes revealed that Ago2 interacts with endogenous small interfering RNAs (endo-siRNAs) to repress mRNA translation globally. Furthermore, the Ago2 and endo-siRNA complexes fine-tune the transcriptome by slicing long terminal repeat retrotransposon-derived chimeric transcripts. The precise mapping of RBP-binding sites in low-input cells opens the door to studying the roles of RBPs in embryonic development and reproductive diseases.

PIWI-interacting RNAs are a class of small RNAs that is most abundantly expressed in animal germline. Substantial research is going on to reveal the functions of piRNAs in the epigenetic and post-transcriptional regulation of transposons and genes. To collect and annotate these data, we developed piRBase, a database assisting piRNA functional study. Since its launch in 2014, piRBase has integrated 264 data sets from 21 organisms, and the number of collected piRNAs has reached 173 million. The latest piRBase release (v2.0, 2018) was more focused on the comprehensive annotation of piRNA sequences, as well as the increasing number of piRNAs. In addition, piRBase release v2.0 also contained the potential information of piRNA targets and disease related piRNA. All datasets in piRBase is free to access, and available for browse, search and bulk downloads at http://www.regulatoryrna.org/database/piRNA/.

Summary The lack of haplotype reference panels and whole-genome sequencing resources specific to the Chinese population has greatly hindered genetic studies in the world’s largest population. Here, we present the NyuWa genome resource, based on deep (26.2×) sequencing of 2,999 Chinese individuals, and construct a NyuWa reference panel of 5,804 haplotypes and 19.3 million variants, which is a high-quality publicly available Chinese population-specific reference panel with thousands of samples. Compared with other panels, the NyuWa reference panel reduces the Han Chinese imputation error rate by a margin ranging from 30% to 51%. Population structure and imputation simulation tests support the applicability of one integrated reference panel for northern and southern Chinese. In addition, a total of 22,504 loss-of-function variants in coding and noncoding genes are identified, including 11,493 novel variants. These results highlight the value of the NyuWa genome resource in facilitating genetic research in Chinese and Asian populations.

The National Genomics Data Center (NGDC) provides a suite of database resources to support worldwide research activities in both academia and industry. With the rapid advancements in higher-throughput and lower-cost sequencing technologies and accordingly the huge volume of multi-omics data generated at exponential scales and rates, NGDC is continually expanding, updating and enriching its core database resources through big data integration and value-added curation. In the past year, efforts for update have been mainly devoted to BioProject, BioSample, GSA, GWH, GVM, NONCODE, LncBook, EWAS Atlas and IC4R. Newly released resources include three human genome databases (PGG.SNV, PGG.Han and CGVD), eLMSG, EWAS Data Hub, GWAS Atlas, iSheep and PADS Arsenal. In addition, four web services, namely, eGPS Cloud, BIG Search, BIG Submission and BIG SSO, have been significantly improved and enhanced. All of these resources along with their services are publicly accessible at https://bigd.big.ac.cn.

Noncoding RNAs (ncRNAs) play crucial regulatory roles in a variety of biological circuits. To document regulatory interactions between ncRNAs and biomolecules, we previously created the NPInter database (http://bigdata.ibp.ac.cn/npinter). Since the last version of NPInter was issued, a rapidly growing number of studies have reported novel interactions and accumulated numerous high-throughput interactome data. We have therefore updated NPInter to its fourth edition in which are integrated 600 000 new experimentally identified ncRNA interactions. ncRNA–DNA interactions derived from ChIRP-seq data and circular RNA interactions have been included in the database. Additionally, disease associations were annotated to the interacting molecules. The database website has also been redesigned with a more user-friendly interface and several additional functional modules. Overall, NPInter v4.0 now provides more comprehensive data and services for researchers working on ncRNAs and their interactions with other biomolecules.

A long-standing question in the field of embryogenesis is how the zygotic genome is precisely activated by maternal factors, allowing normal early embryonic development. We have previously shown that N6-methyladenine (6mA) DNA modification is highly dynamic in early Drosophila embryos and forms an epigenetic mark. However, little is known about how 6mA-formed epigenetic information is decoded. Here we report that the Fox-family protein Jumu binds 6mA-marked DNA and acts as a maternal factor to regulate the maternal-to-zygotic transition. We find that zelda encoding the pioneer factor Zelda is marked by 6mA. Our genetic assays suggest that Jumu controls the proper zygotic genome activation (ZGA) in early embryos, at least in part, by regulating zelda expression. Thus, our findings not only support that the 6mA-formed epigenetic marks can be read by specific transcription factors, but also uncover a mechanism by which the Jumu regulates ZGA partially through Zelda in early embryos., N6-methyladenine (6mA) DNA modification is a dynamic epigenetic mark in Drosophila embryos, but how 6mA is decoded is unclear. Here, the authors show that the protein Jumu binds 6mA-marked DNA to regulate the maternal-to-zygotic transition, partly through regulation the expression of the 6mA marked pioneer factor zelda.

Horizontal gene transfer (HGT) is one of the major mechanisms contributing to microbial genome diversification. A number of computational methods for finding horizontally transferred genes have been proposed in the past decades; however none of them has provided a reliable detector yet. In existing parametric approaches, only one single compositional property can participate in the detection process, or the results obtained through each single property are just simply combined. It's known that different properties may mean different information, so the single property can't sufficiently contain the information encoded by gene sequences. In addition, the class imbalance problem in the datasets, which also results in great errors for the gene detection, hasn't been considered by the published methods. Here we developed an effective classifier system (Hgtident) that used support vector machine (SVM) by combining unusual properties effectively for HGT detection.Our approach Hgtident includes the introduction of more representative datasets, optimization of SVM model, feature selection, handling of imbalance problem in the datasets and extensive performance evaluation via systematic cross-validation methods. Through feature selection, we found that JS-DN and JS-CB have higher discriminating power for HGT detection, while GC1-GC3 and k-mer (k = 1, 2, …, 7) make the least contribution. Extensive experiments indicated the new classifier could reduce Mean error dramatically, and also improve Recall by a certain level. For the testing genomes, compared with the existing popular multiple-threshold approach, on average, our Recall and Mean error was respectively improved by 2.81% and reduced by 26.32%, which means that numerous false positives were identified correctly.Hgtident introduced here is an effective approach for better detecting HGT. Combining multiple features of HGT is also essential for a wider range of HGT events detection.

Mutations in several general pre-mRNA splicing factors have been linked to myelodysplastic syndromes (MDSs) and solid tumors. These mutations have generally been assumed to cause disease by the resultant splicing defects, but different mutations appear to induce distinct splicing defects, raising the possibility that an alternative common mechanism is involved. Here we report a chain of events triggered by multiple splicing factor mutations, especially high-risk alleles in SRSF2 and U2AF1, including elevated R-loops, replication stress, and activation of the ataxia telangiectasia and Rad3-related protein (ATR)-Chk1 pathway. We further demonstrate that enhanced R-loops, opposite to the expectation from gained RNA binding with mutant SRSF2, result from impaired transcription pause release because the mutant protein loses its ability to extract the RNA polymerase II (Pol II) C-terminal domain (CTD) kinase—the positive transcription elongation factor complex (P-TEFb)—from the 7SK complex. Enhanced R-loops are linked to compromised proliferation of bone-marrow-derived blood progenitors, which can be partially rescued by RNase H overexpression, suggesting a direct contribution of augmented R-loops to the MDS phenotype.

RNAcentral is a comprehensive database of non-coding RNA (ncRNA) sequences, collating information on ncRNA sequences of all types from a broad range of organisms. We have recently added a new genome mapping pipeline that identifies genomic locations for ncRNA sequences in 296 species. We have also added several new types of functional annotations, such as tRNA secondary structures, Gene Ontology annotations, and miRNA-target interactions. A new quality control mechanism based on Rfam family assignments identifies potential contamination, incomplete sequences, and more. The RNAcentral database has become a vital component of many workflows in the RNA community, serving as both the primary source of sequence data for academic and commercial groups, as well as a source of stable accessions for the annotation of genomic and functional features. These examples are facilitated by an improved RNAcentral web interface, which features an updated genome browser, a new sequence feature viewer, and improved text search functionality. RNAcentral is freely available at https://rnacentral.org.

Whole genome sequencing technology has facilitated the discovery of a large number of somatic mutations in enhancers (SMEs), whereas the utility of SMEs in tumorigenesis has not been fully explored. Here we present Ennet, a method to comprehensively investigate SMEs enriched networks (SME-networks) in cancer by integrating SMEs, enhancer-gene interactions and gene-gene interactions. Using Ennet, we performed a pan-cancer analysis in 2004 samples from 8 cancer types and found many well-known cancer drivers were involved in the SME-networks, includingESR1,SMAD3,MYC,EGFR,BCL2andPAX5. Meanwhile, Ennet also identified many new networks with less characterization but have potentially important roles in cancer, including a large SME-network in medulloblastoma (MB), which contains genes enriched in the glutamate receptor and neural development pathways. Interestingly, SME-networks are specific across cancer types, and the vast majority of the genes identified by Ennet have few mutations in gene bodies. Collectively, our work suggests that using enhancer-only somatic mutations can be an effective way to discover potential cancer-driving networks. Ennet provides a new perspective to explore new mechanisms for tumor progression from SMEs.

Biological processes, especially developmental processes, are often dynamic. Previous BodyMap projects for human and mouse have provided researchers with portals to tissue-specific gene expression, but these efforts have not included dynamic gene expression patterns. Over the past few years, substantial progress in our understanding of the molecular mechanisms of protein-coding and long noncoding RNA (lncRNA) genes in development processes has been achieved through numerous time series RNA sequencing (RNA-seq) studies. However, none of the existing databases focuses on these time series data, thus rendering the exploration of dynamic gene expression patterns inconvenient. Here, we present Dynamic BodyMap (Dynamic-BM), a database for temporal gene expression profiles, obtained from 2203 time series of RNA-seq samples, covering >25 tissues from five species. Dynamic-BM has a user-friendly Web interface designed for browsing and searching the dynamic expression pattern of genes from different sources. It is an open resource for efficient data exploration, providing dynamic expression profiles of both protein-coding genes and lncRNAs to facilitate the generation of new hypotheses in developmental biology research. Additionally, Dynamic-BM includes a literature-based knowledgebase for lncRNAs associated with tissue development and a list of manually selected lncRNA candidates that may be involved in tissue development. Dynamic-BM is available at http://bioinfo.ibp.ac.cn/Dynamic-BM.

Eukaryotic horizontal gene transfer (HGT) events are increasingly being discovered yet few reports have summarized multiple occurrences in a wide range of species. We systematically investigated HGT events in the order Lepidoptera by employing a series of filters. Bombyx mori, Danaus plexippus and Heliconius melpomene had 13, 12 and 12 HGTs, respectively, from bacteria and fungi. These HGTs contributed a total of 64 predicted genes: 22 to B. mori, 22 to D. plexippus and 20 to H. melpomene. Several new genes were generated by post-transfer duplications. Post-transfer duplication of a suite of functional HGTs has rarely been reported in higher organisms. The distributional patterns of paralogues for certain genes differed in the three species, indicating potential independent duplication or loss events. All of these HGTs had homologues expressed in some other lepidopterans, indicating ancient transfer events. Most HGTs were involved in the metabolism of sugar and amino acids. These HGTs appeared to have experienced amelioration, purifying selection and accelerated evolution to adapt to the background genome of the recipient. The discovery of ancient, massive HGTs and duplications in lepidopterans and their adaptive evolution provides further insights into the evolutionary significance of the events from donors to multicellular host recipients.

Motivation Long non-coding RNAs (lncRNAs) are essential in many molecular pathways, and are frequently associated with disease but the mechanisms of most lncRNAs have not yet been characterized. Genetic variations, including single nucleotide polymorphisms (SNPs) and structural variations, are widely distributed in the genome, including lncRNA gene regions. As the number of studies on lncRNAs grows rapidly, it is necessary to evaluate the effects of genetic variations on lncRNAs. Results Here, we present LncVar, a database of genetic variation associated with long non-coding genes in six species. We collected lncRNAs from the NONCODE database, and evaluated their conservation. We systematically integrated transcription factor binding sites and m6A modification sites of lncRNAs and provided comprehensive effects of SNPs on transcription and modification of lncRNAs. We collected putatively translated open reading frames (ORFs) in lncRNAs, and identified both synonymous and non-synonymous SNPs in ORFs. We also collected expression quantitative trait loci of lncRNAs from the literature. Furthermore, we identified lncRNAs in CNV regions as prognostic biomarker candidates of cancers and predicted lncRNA gene fusion events from RNA-seq data from cell lines. The LncVar database can be used as a resource to evaluate the effects of the variations on the biological function of lncRNAs. Availability and Implementation LncVar is available at http://bioinfo.ibp.ac.cn/LncVar. Supplementary information Supplementary materials are available at Bioinformatics online.

Invertebrates can acquire functional genes via horizontal gene transfer (HGT) from bacteria but fishes are not known to do so. We provide the first reliable evidence of one HGT event from marine bacteria to fishes. The HGT appears to have occurred after emergence of the teleosts. The transferred gene is expressed and regulated developmentally. Its successful integration and expression may change the genetic and metabolic repertoire of fishes. In addition, this gene contains conserved domains and similar tertiary structures in fishes and their putative donor bacteria. Thus, it may function similarly in both groups. Evolutionary analyses indicate that it evolved under purifying selection, further indicating its conserved function. We document the first likely case of HGT of functional gene from prokaryote to fishes. This discovery certifies that HGT can influence vertebrate evolution.

microRNAs (miRNAs) represent a class of noncoding RNA species, believed to be regulating gene expression by binding to complementary sites in the 3'UTRs of target mRNAs. They play important regulatory roles in various metabolic pathways in most eukaryotes. The recent discovery of virus encoded miRNAs suggests that viruses may be using them to regulate host and viral gene expression. Another class of closely related small interfering RNAs (siRNAs) also has been found within the HIV-1 genome and shown to be exerting a limited impact on virus reproduction. Additionally, an additional type of viral noncoding RNAs named small noncoding RNAs (sncRNAs) ranging from a few tens to a few hundred nucleotides in length, has also been identified. sncRNAs have a wide phylogenesis and high levels of expression, suggesting they may play an important roles in different species. Here we discuss the genomic organization, expression, conservation as well as potential function of virally encoded miRNA, siRNA, and sncRNAs.

Summary: We here present BioCircos.js, an interactive and lightweight JavaScript library especially for biological data interactive visualization. BioCircos.js facilitates the development of web-based applications for circular visualization of various biological data, such as genomic features, genetic variations, gene expression and biomolecular interactions. Availability and implementation: BioCircos.js and its manual are freely available online at http://bioinfo.ibp.ac.cn/biocircos/. Contact: rschen@ibp.ac.cn Supplementary information: Supplementary data are available at Bioinformatics online.

Motivation: PIWI-interacting RNAs (piRNAs) are a class of small non-coding RNAs that are highly abundant in the germline. One important role of piRNAs is to defend genome integrity by guiding PIWI proteins to silence transposable elements (TEs), which have a high potential to cause deleterious effects on their host. The mechanism of piRNA-mediated post-transcriptional silencing was also observed to affect mRNAs, suggesting that piRNAs might play a broad role in gene expression regulation. However, there has been no systematic report with regard to how many protein-coding genes might be targeted and regulated by piRNAs. Results: We trained a support vector machine classifier based on a combination of Miwi CLIP-Seq-derived features and position-derived features to predict the potential targets of piRNAs on mRNAs in the mouse. Reanalysis of a published microarray dataset suggested that the expression level of the 2587 protein-coding genes predicted as piRNA targets showed significant upregulation as a whole after abolishing the slicer activity of Miwi, supporting the conclusion that they are subject to piRNA-mediated regulation. Availability and implementation: A web version of the method called pirnaPre as well as our results for browse is available at http://www.regulatoryrna.org/software/piRNA/piRNA_target_mRNA/index.php. Contact: crs@sun5.ibp.ac.cn or heshunmin@gmail.com Supplementary information: Supplementary data are available at Bioinformatics online.

Fig wasps exhibit extreme intraspecific morphological divergence in the wings, compound eyes, antennae, body color and size. Corresponding to this, behaviors and lifestyles between two sexes are also different: females can emerge from fig and fly to other fig tree to oviposit and pollinate, while males live inside fig for all their lifetime. Genetic regulation may drive these extreme intraspecific morphological and behavioral divergence. Transcription factors (TFs) involved in morphological development and physiological activity may exhibit sex-specific expressions. Herein, we detect 865 TFs by using genomic and transcriptomic data of the fig wasp Ceratosolen solmsi. Analyses of transcriptomic data indicated that up-regulated TFs in females show significant enrichment in development of the wing, eye and antenna in all stages, from larva to adult. Meanwhile, TFs related to the development of a variety of organs display sex-specific patterns of expression in the adults and these may contribute significantly to their sexual dimorphism. In addition, up-regulated TFs in adult males exhibit enrichment in genitalia development and circadian rhythm, which correspond with mating and protandry. This finding is consistent with their sex-specific behaviors. In conclusion, our results strongly indicate that TFs play important roles in the sexual dimorphism of fig wasps.

The piRNA machinery is known for its role in mediating epigenetic silencing of transposons. Recent studies suggest that this function also involves piRNA-guided cleavage of transposon-derived transcripts. As many piRNAs also appear to have the capacity to target diverse mRNAs, this raises the intriguing possibility that piRNAs may act extensively as siRNAs to degrade specific mRNAs. To directly test this hypothesis, we compared mouse PIWI (MIWI)-associated piRNAs with experimentally identified cleaved mRNA fragments from mouse testes, and observed cleavage sites that predominantly occur at position 10 from the 5' end of putative targeting piRNAs. We also noted strong biases for U and A residues at nucleotide positions 1 and 10, respectively, in both piRNAs and mRNA fragments, features that resemble the pattern of piRNA amplification by the 'ping-pong' cycle. Through mapping of MIWI-RNA interactions by CLIP-seq and gene expression profiling, we found that many potential piRNA-targeted mRNAs directly interact with MIWI and show elevated expression levels in the testes of Miwi catalytic mutant mice. Reporter-based assays further revealed the importance of base pairing between piRNAs and mRNA targets and the requirement for both the slicer activity and piRNA-loading ability of MIWI in piRNA-mediated target repression. Importantly, we demonstrated that proper turnover of certain key piRNA targets is essential for sperm formation. Together, these findings reveal the siRNA-like function of the piRNA machinery in mouse testes and its central requirement for male germ cell development and maturation.

Small proteins is the general term for proteins with length shorter than 100 amino acids. Identification and functional studies of small proteins have advanced rapidly in recent years, and several studies have shown that small proteins play important roles in diverse functions including development, muscle contraction and DNA repair. Identification and characterization of previously unrecognized small proteins may contribute in important ways to cell biology and human health. Current databases are generally somewhat deficient in that they have either not collected small proteins systematically, or contain only predictions of small proteins in a limited number of tissues and species. Here, we present a specifically designed web-accessible database, small proteins database (SmProt, http://bioinfo.ibp.ac.cn/SmProt), which is a database documenting small proteins. The current release of SmProt incorporates 255 010 small proteins computationally or experimentally identified in 291 cell lines/tissues derived from eight popular species. The database provides a variety of data including basic information (sequence, location, gene name, organism, etc.) as well as specific information (experiment, function, disease type, etc.). To facilitate data extraction, SmProt supports multiple search options, including species, genome location, gene name and their aliases, cell lines/tissues, ORF type, gene type, PubMed ID and SmProt ID. SmProt also incorporates a service for the BLAST alignment search and provides a local UCSC Genome Browser. Additionally, SmProt defines a high-confidence set of small proteins and predicts the functions of the small proteins.

piRNAs are a class of small RNAs that is most abundantly expressed in the animal germ line. Presently, substantial research is going on to reveal the functions of piRNAs in the epigenetic and post-transcriptional regulation of transposons and genes. A piRNA database for collection, annotation and structuring of these data will be a valuable contribution to the field, and we have therefore developed the piRBase platform which integrates various piRNA-related high-throughput data. piRBase has the largest collection of piRNAs among existing databases, and contains at present 77 million piRNA sequences from nine organisms. Repeat-derived and gene-derived piRNAs, which possibly participate in the regulation of the corresponding elements, have been given particular attention. Furthermore, epigenetic data and reported piRNA targets were also collected. To our knowledge, this is the first piRNA database that systematically integrates epigenetic and post-transcriptional regulation data to support piRNA functional analysis. We believe that piRBase will contribute to a better understanding of the piRNA functions. Database URL: http://www.regulatoryrna.org/database/piRNA/

Background In metazoans, Piwi-related Argonaute proteins play important roles in maintaining germline integrity and fertility and have been linked to a class of germline-enriched small RNAs termed piRNAs. Caenorhabditis elegans encodes two Piwi family proteins called PRG-1 and PRG-2, and PRG-1 interacts with the C. elegans piRNAs (21U-RNAs). Previous studies found that mutation of prg-1 causes a marked reduction in the expression of 21U-RNAs, temperature-sensitive defects in fertility and other phenotypic defects. Results In this study, we wanted to systematically demonstrate the function of PRG-1 in the regulation of small RNAs and their targets. By analyzing small RNAs and mRNAs with and without a mutation in prg-1 during C. elegans development, we demonstrated that (1) mutation of prg-1 leads to a decrease in the expression of 21U-RNAs, and causes 35 ~ 40% of miRNAs to be down-regulated; (2) in C. elegans, approximately 3% (6% in L4) of protein-coding genes are differentially expressed after mutating prg-1, and 60 ~ 70% of these substantially altered protein-coding genes are up-regulated; (3) the target genes of the down-regulated miRNAs and the candidate target genes of the down-regulated 21U-RNAs are enriched in the up-regulated protein-coding genes; and (4) PRG-1 regulates protein-coding genes by down-regulating small RNAs (miRNAs and 21U-RNAs) that target genes that participate in the development of C. elegans. Conclusions In prg-1-mutated C. elegans, the expression of miRNAs and 21U-RNAs was reduced, and the protein-coding targets, which were associated with the development of C. elegans, were up-regulated. This may be the mechanism underlying PRG-1 function. Electronic supplementary material The online version of this article (doi:10.1186/1471-2164-15-321) contains supplementary material, which is available to authorized users.

The NONCODE database is an integrated knowledge database designed for the analysis of non-coding RNAs (ncRNAs). Since NONCODE was first released 3 years ago, the number of known ncRNAs has grown rapidly, and there is growing recognition that ncRNAs play important regulatory roles in most organisms. In the updated version of NONCODE (NONCODE v2.0), the number of collected ncRNAs has reached 206 226, including a wide range of microRNAs, Piwi-interacting RNAs and mRNA-like ncRNAs. The improvements brought to the database include not only new and updated ncRNA data sets, but also an incorporation of BLAST alignment search service and access through our custom UCSC Genome Browser. NONCODE can be found under http://www.noncode.org or http://noncode.bioinfo.org.cn.

Background Fig pollinating wasps form obligate symbioses with their fig hosts. This mutualism arose approximately 75 million years ago. Unlike many other intimate symbioses, which involve vertical transmission of symbionts to host offspring, female fig wasps fly great distances to transfer horizontally between hosts. In contrast, male wasps are wingless and cannot disperse. Symbionts that keep intimate contact with their hosts often show genome reduction, but it is not clear if the wide dispersal of female fig wasps will counteract this general tendency. We sequenced the genome of the fig wasp Ceratosolen solmsi to address this question. Results The genome size of the fig wasp C. solmsi is typical of insects, but has undergone dramatic reductions of gene families involved in environmental sensing and detoxification. The streamlined chemosensory ability reflects the overwhelming importance of females finding trees of their only host species, Ficus hispida, during their fleeting adult lives. Despite long-distance dispersal, little need exists for detoxification or environmental protection because fig wasps spend nearly all of their lives inside a largely benign host. Analyses of transcriptomes in females and males at four key life stages reveal that the extreme anatomical sexual dimorphism of fig wasps may result from a strong bias in sex-differential gene expression. Conclusions Our comparison of the C. solmsi genome with other insects provides new insights into the evolution of obligate mutualism. The draft genome of the fig wasp, and transcriptomic comparisons between both sexes at four different life stages, provide insights into the molecular basis for the extreme anatomical sexual dimorphism of this species.