Coffea arabica, an allotetraploid hybrid of Coffea eugenioides and Coffea canephora, is the source of approximately 60% of coffee products worldwide, and its cultivated accessions have undergone several population bottlenecks. We present chromosome-level assemblies of a di-haploid C. arabica accession and modern representatives of its diploid progenitors, C. eugenioides and C. canephora. The three species exhibit largely conserved genome structures between diploid parents and descendant subgenomes, with no obvious global subgenome dominance. We find evidence for a founding polyploidy event 350,000–610,000 years ago, followed by several pre-domestication bottlenecks, resulting in narrow genetic variation. A split between wild accessions and cultivar progenitors occurred ~30.5 thousand years ago, followed by a period of migration between the two populations. Analysis of modern varieties, including lines historically introgressed with C. canephora, highlights their breeding histories and loci that may contribute to pathogen resistance, laying the groundwork for future genomics-based breeding of C. arabica.
Circulating nucleic acids are found in free form in body fluids and may serve as minimally invasive tools for cancer diagnosis and prognosis. Only a few studies have investigated the potential application of circulating mRNAs and microRNAs (miRNAs) in prostate cancer (PCa). The Cancer Genome Atlas (TCGA) database was used for an in silico analysis to identify circulating mRNA and miRNA as potential markers of PCa. A total of 2,267 genes and 49 miRNAs were differentially expressed between normal and tumor samples. The prediction analyses of target genes and integrative analysis of mRNA and miRNA expression revealed eleven genes and eight miRNAs which were validated by RT-qPCR in plasma samples from 102 untreated PCa patients and 50 cancer-free individuals. Two genes, OR51E2 and SIM2, and two miRNAs, miR-200c and miR-200b, showed significant association with PCa. Expression levels of these transcripts distinguished PCa patients from controls (67% sensitivity and 75% specificity). PCa patients and controls with prostate-specific antigen (PSA) ≤ 4.0 ng/mL were discriminated based on OR51E2 and SIM2 expression levels. The miR-200c expression showed association with Gleason score and miR-200b, with bone metastasis, bilateral tumor, and PSA > 10.0 ng/mL. The combination of circulating mRNA and miRNA was useful for the diagnosis and prognosis of PCa.
Long noncoding RNAs (lncRNAs) correspond to a eukaryotic noncoding RNA class that gained great attention in the past years as a higher layer of regulation for gene expression in cells. There is, however, a lack of specific computational approaches to reliably predict lncRNA in plants, which contrast the variety of prediction tools available for mammalian lncRNAs. This distinction is not that obvious, given that biological features and mechanisms generating lncRNAs in the cell are likely different between animals and plants. Considering this, we present a machine learning analysis and a classifier approach called RNAplonc (https://github.com/TatianneNegri/RNAplonc/) to identify lncRNAs in plants. Our feature selection analysis considered 5468 features, and it used only 16 features to robustly identify lncRNA with the REPTree algorithm. That was the base to create the model and train it with lncRNA and mRNA data from five plant species (thale cress, cucumber, soybean, poplar and Asian rice). After an extensive comparison with other tools largely used in plants (CPC, CPC2, CPAT and PLncPRO), we found that RNAplonc produced more reliable lncRNA predictions from plant transcripts with 87.5% of the best result in eight tests in eight species from the GreeNC database and four independent studies in monocotyledonous (Brachypodium) and eudicotyledonous (Populus and Gossypium) species.
Transposable elements (TEs) are the most represented sequences occurring in eukaryotic genomes. Few methods provide the classification of these sequences into deeper levels, such as superfamily level, which could provide useful and detailed information about these sequences. Most methods that classify TE sequences use handcrafted features such as k-mers and homology-based search, which could be inefficient for classifying non-homologous sequences. Here we propose an approach, called transposable elements pepresentation learner (TERL), that preprocesses and transforms one-dimensional sequences into two-dimensional space data (i.e., image-like data of the sequences) and apply it to deep convolutional neural networks. This classification method tries to learn the best representation of the input data to classify it correctly. We have conducted six experiments to test the performance of TERL against other methods. Our approach obtained macro mean accuracies and F1-score of 96.4% and 85.8% for superfamilies and 95.7% and 91.5% for the order sequences from RepBase, respectively. We have also obtained macro mean accuracies and F1-score of 95.0% and 70.6% for sequences from seven databases into superfamily level and 89.3% and 73.9% for the order level, respectively. We surpassed accuracy, recall and specificity obtained by other methods on the experiment with the classification of order level sequences from seven databases and surpassed by far the time elapsed of any other method for all experiments. Therefore, TERL can learn how to predict any hierarchical level of the TEs classification system and is about 20 times and three orders of magnitude faster than TEclass and PASTEC, respectively https://github.com/muriloHoracio/TERL. Contact:murilocruz@alunos.utfpr.edu.br
Large-scale transcriptome projects have shown that the number of RNA transcripts not coding for proteins (non-coding RNAs) is much larger than previously recognized. High-throughput technologies, coupled with bioinformatics approaches, have produced increasing amounts of data, highlighting the role of non-coding RNAs (ncRNAs) in biological processes. Data generated by these studies include diverse non-coding RNA classes from organisms of different kingdoms, which were obtained using different experimental and computational assays. This has led to a rapid increase of specialized RNA databases. The fast growth in the number of available databases makes integration of stored information a difficult task. We present here NRDR, a Non-coding RNA Databases Resource for information retrieval on ncRNA databases (www.ncrnadatabases.org). We performed a survey of 102 public databases on ncRNAs and we have introduced four categorizations to classify these databases and to help researchers quickly search and find the information they need: RNA family, information source, information content and available search mechanisms. NRDR is a useful databases searching tool that will facilitate research on ncRNAs.
8 September 2018 in Functional & Integrative Genomics
Coffea arabica L. is an important agricultural commodity, accounting for 60% of traded coffee worldwide. Nitrogen (N) is a macronutrient that is usually limiting to plant yield; however, molecular mechanisms of plant acclimation to N limitation remain largely unknown in tropical woody crops. In this study, we investigated the transcriptome of coffee roots under N starvation, analyzing poly-A+ libraries and small RNAs. We also evaluated the concentration of selected amino acids and N-source preferences in roots. Ammonium was preferentially taken up over nitrate, and asparagine and glutamate were the most abundant amino acids observed in coffee roots. We obtained 34,654 assembled contigs by mRNA sequencing, and validated the transcriptional profile of 12 genes by RT-qPCR. Illumina small RNA sequencing yielded 8,524,332 non-redundant reads, resulting in the identification of 86 microRNA families targeting 253 genes. The transcriptional pattern of eight miRNA families was also validated. To our knowledge, this is the first catalog of differentially regulated amino acids, N sources, mRNAs, and sRNAs in Arabica coffee roots.