Search

Wild castor grows in the high-altitude tropical desert of the African Plateau, a region known for high ultraviolet radiation, strong light, and extremely dry condition. To investigate the potential genetic basis of adaptation to both highland and tropical deserts, we generated a chromosome-level genome sequence assembly of the wild castor accession WT05, with a genome size of 316 Mb, a scaffold N50 of 31.93 Mb, and a contig N50 of 8.96 Mb, respectively. Compared with cultivated castor and other Euphorbiaceae species, the wild castor exhibits positive selection and gene family expansion for genes involved in DNA repair, photosynthesis, and abiotic stress responses. Genetic variations associated with positive selection were identified in several key genes, such as LIG1, DDB2, and RECG1, involved in nucleotide excision repair. Moreover, a study of genomic diversity among wild and cultivated accessions revealed genomic regions containing selection signatures associated with the adaptation to extreme environments. The identification of the genes and alleles with selection signatures provides insights into the genetic mechanisms underlying the adaptation of wild castor to the high-altitude tropical desert and would facilitate direct improvement of modern castor varieties.

Oil produced by castor (Ricinus communis) has broad industrial applications. However, knowledge on the genetic diversity, especially genetic alterations that occurred during domestication and subsequent traits selection, of this oil crop is limited. Here, our population genomics analyses show that the Chinese castors have developed a geographic pattern, classified into the southern-, the middle-, and the northern-China groups. We detect a number of candidate genomic loci that are associated with the selection signals during the geographical differentiation and domestication. Using genome-wide association analysis, we identify candidate genes associated with nine agronomically important traits. One of the candidate genes encoding a glycosyltransferase related to cellulose and lignin biosynthesis is associated with both capsule dehiscence and endocarp thickness. We hypothesize that the abundance of cellulose or lignin in endocarp is an important factor for capsule dehiscence. Our results provide foundation for castor breeding and genetic study., Castor is an important industrial oil crop, but knowledge on its genetic diversity is limited. Here, Fan et al. show geographic pattern of Chinese castors that have developed during domestication by population genetic analyses, and reveal candidate genes associated with agronomically important traits.

Wild castor grows in the high-altitude tropical desert of the African Plateau, a region known for high ultraviolet radiations, strong lights and extremely dry conditions. To investigate the potential genetic basis of adaptation to both highland and tropical deserts, we generated a chromosome-level genome sequence of the wild castor accession WT05, with genome size of 316 Mb and scaffold and contig N50 sizes of 31.93 Mb and 8.96 Mb, respectively. Compared with cultivated castor and other Euphorbiaceae species, the wild castor exhibits positive selection and gene family expansion for genes involved in DNA repair, photosynthesis and abiotic stress responses. Genetic variations associated with positive selection were identified in several key genes, such as LIG1, DDB2, and RECG1, involved in nucleotide excision repair. Moreover, a study of genomic diversity among wild and cultivated accessions revealed genomic regions containing selection signatures associated with the adaptation to extreme environments. The identification of the genes and alleles with selection signatures provides insights into the genetic mechanisms underlying the adaptation of wild castor to the high-altitude tropical desert and would facilitate direct improvement of modern castor varieties.

RNA editing is a post-transcriptional or cotranscriptional process that changes the sequence of the precursor transcript by substitutions, insertions, or deletions. Almost all of the land plants undergo RNA editing in organelles (plastids and mitochondria). Although several software tools have been developed to identify RNA editing events, there has been a great challenge to distinguish true RNA editing events from genome variation, sequencing errors, and other factors. Here we introduce REDO, a comprehensive application tool for identifying RNA editing events in plant organelles based on variant call format files from RNA-sequencing data. REDO is a suite of Perl scripts that illustrate a bunch of attributes of RNA editing events in figures and tables. REDO can also detect RNA editing events in multiple samples simultaneously and identify the significant differential proportion of RNA editing loci. Comparing with similar tools, such as REDItools, REDO runs faster with higher accuracy, and more specificity at the cost of slightly lower sensitivity. Moreover, REDO annotates each RNA editing site in RNAs, whereas REDItools reports only possible RNA editing sites in genome, which need additional steps to obtain RNA editing profiles for RNAs. Overall, REDO can identify potential RNA editing sites easily and provide several functions such as detailed annotations, statistics, figures, and significantly differential proportion of RNA editing sites among different samples.

Recently, RNA-seq has become widely used technology for transcriptome profiling due to its single-base accuracy and high-throughput speciality. In this study, we applied a computational approach on an integrated RNA-seq dataset across 15 normal mouse tissues, and consequently assigned 8408 house-keeping (HK) genes and 2581 tissue-specific (TS) genes among UCSC RefGene annotation. Apart from some basic genomic features, we also performed expression, function and pathway analysis with clustering, DAVID and Ingenuity Pathway Analysis, indicating the physiological connections (tissues) and diverse biological roles of HK genes (fundamental processes) and TS genes (tissue-corresponding processes). Moreover, we used RT-PCR method to test 18 candidate HK genes and finally identified a novel list of highly stable internal control genes: Ywhae, Ddb 1, Eif4h, etc. In summary, this study provides a new HK gene and TS gene resource for further genetic and evolution research and helps us better understand morphogenesis and biological diversity in mouse.

Background: Date palm (Phoenix dactylifera L.) is a cultivated woody plant with agricultural and economic importance in many countries around the world. With the advantages of next generation sequencing technologies, genome sequences for many date palm cultivars have been released recently. Short sequence repeat (SSR) and single nucleotide polymorphism (SNP) can be identified from these genomic data, and have been proven to be very useful biomarkers in plant genome analysis and breeding. Results: Here, we first improved the date palm genome assembly using 130X of HiSeq data generated in our lab. Then 246,445 SSRs (214,901 SSRs and 31,544 compound SSRs) were annotated in this genome assembly; among the SSRs, mononucleotide SSRs (58.92%) were the most abundant, followed by di- (29.92%), tri- (8.14%), tetra- (2.47%), penta- (0.36%), and hexa-nucleotide SSRs (0.19%). The high-quality PCR primer pairs were designed for most (174,497; 70.81% out of total) SSRs. We also annotated 6,375,806 SNPs with raw read depth≥3 in 90% cultivars. To further reduce false positive SNPs, we only kept 5,572,650 (87.40% out of total) SNPs with at least 20% cultivars support for downstream analyses. The high-quality PCR primer pairs were also obtained for 4,177,778 (65.53%) SNPs. We reconstructed the phylogenetic relationships among the 62 cultivars using these variants and found that they can be divided into three clusters, namely North Africa, Egypt – Sudan, and Middle East – South Asian, with Egypt – Sudan being the admixture of North Africa and Middle East – South Asian cultivars; we further confirmed these clusters using principal component analysis. Moreover, 34,346 SSRs and 4,177,778 SNPs with PCR primers were assigned to shared cultivars for cultivar classification and diversity analysis. All these SSRs, SNPs and their classification are available in our database, and can be used for cultivar identification, comparison, and molecular breeding. Conclusion: DRDB is a comprehensive genomic resource database of date palm. It can serve as a bioinformatics platform for date palm genomics, genetics, and molecular breeding. DRDB is freely available at http://drdb.big.ac.cn/home.

Background Exon recognition and splicing precisely and efficiently by spliceosome is the key to generate mature mRNAs. About one third or a half of disease-related mutations affect RNA splicing. Software PVAAS has been developed to identify variants associated with aberrant splicing by directly using RNA-seq data. However, it bases on the assumption that annotated splicing site is normal splicing, which is not true in fact. Results We develop the ISVASE, a tool for specifically identifying sequence variants associated with splicing events (SVASE) by using RNA-seq data. Comparing with PVAAS, our tool has several advantages, such as multi-pass stringent rule-dependent filters and statistical filters, only using split-reads, independent sequence variant identification in each part of splicing (junction), sequence variant detection for both of known and novel splicing event, additional exon-exon junction shift event detection if known splicing events provided, splicing signal evaluation, known DNA mutation and/or RNA editing data supported, higher precision and consistency, and short running time. Using a realistic RNA-seq dataset, we performed a case study to illustrate the functionality and effectiveness of our method. Moreover, the output of SVASEs can be used for downstream analysis such as splicing regulatory element study and sequence variant functional analysis. Conclusions ISVASE is useful for researchers interested in sequence variants (DNA mutation and/or RNA editing) associated with splicing events. The package is freely available at https://sourceforge.net/projects/isvase/. Electronic supplementary material The online version of this article (doi:10.1186/s12859-017-1732-7) contains supplementary material, which is available to authorized users.

The red palm weevil (RPW), Rhynchophorus ferrugineus (Coleoptera: Curculionidae), is an invasive, concealed and destructive tissue borer, and it becomes a lethal pest of the palm family of plants and has been reported to attack 20 palm species around the globe. Here we report a systematic transcriptomic study on embryogenesis of RPW, where we analyze the transcriptomes across five developmental stages of RPW embryogenesis, involving four embryonic stages (E1, E2, E3 and E4) and one larval stage (L1). Using the RNA-seq and next-generation platforms, we generated 80 to 91 million reads for each library and assemble 22 532 genes that are expressed at different embryonic stages. Among the total transcripts from the five embryonic development stages, we found that 30.45 % are differentially expressed, 10.10 % show stage-specificity and even a larger fraction, 62.88 %, exhibit constitutive expression in all the stages. We also analyzes the expression dynamics of several conserved signaling pathways (such as Hedgehog, JAK-STAT, Notch, TGF-β, Ras/MAPK and Wnt), as well as key developmental genes, including those related to apoptosis, axis formation, Hox complex, neurogenesis and segmentation. The datasets provide an essential resource for gene annotation and RPW functional genomics, including studies by using tools and concepts from multiple disciplines, such as development, physiology, biochemistry, molecular biology and genetics.

MicroRNAs (miRNAs) are small endogenous noncoding RNAs. Plant miRNAs are known to play important regulatory roles in homeostasis, stress response, and diverse developmental processes. Here, we describe the identification of conserved miRNAs in date palm (Phoenix dactylifera L.) based on transcriptomic data acquired across multistage fruit development and genome sequences, which include 238 plant conserved miRNAs and 276 novel P. dactylifera-specific miRNAs.

The red palm weevil (RPW; Rhynchophorus ferrugineus) is a devastating pest of palms, prevalent in the Middle East as well as many other regions of the world. Here, we report a large-scale de novo complementary DNA (cDNA) sequencing effort that acquired ∼5 million reads and assembled them into 26 765 contigs from 12 libraries made from samples of different RPW developmental stages based on the Roche/454 GS FLX platform. We annotated these contigs based on the publically available known insect genes and the Tribolium castaneum genome assembly. We find that over 80% of coding sequences (CDS) from the RPW contigs have high-identity homologs to known proteins with complete CDS. Gene expression analysis shows that the pupa and larval stages have the highest and lowest expression levels, respectively. In addition, we also identified more than 60 000 single nucleotide polymorphisms and 1 200 simple sequence repeat markers. This study provides the first large-scale cDNA dataset for RPW, a much-needed resource for future molecular studies.

Coconut (Cocos nucifera L.), a member of the palm family (Arecaceae), is one of the most economically important crops in tropics, serving as an important source of food, drink, fuel, medicine, and construction material. Here we report an assembly of the coconut (C. nucifera, Oman local Tall cultivar) mitochondrial (mt) genome based on next-generation sequencing data. This genome, 678,653bp in length and 45.5% in GC content, encodes 72 proteins, 9 pseudogenes, 23 tRNAs, and 3 ribosomal RNAs. Within the assembly, we find that the chloroplast (cp) derived regions account for 5.07% of the total assembly length, including 13 proteins, 2 pseudogenes, and 11 tRNAs. The mt genome has a relatively large fraction of repeat content (17.26%), including both forward (tandem) and inverted (palindromic) repeats. Sequence variation analysis shows that the Ti/Tv ratio of the mt genome is lower as compared to that of the nuclear genome and neutral expectation. By combining public RNA-Seq data for coconut, we identify 734 RNA editing sites supported by at least two datasets. In summary, our data provides the second complete mt genome sequence in the family Arecaceae, essential for further investigations on mitochondrial biology of seed plants.

Spodoptera litura is a threatening agricultural insect in tropical and subtropical areas and accounts for tremendous annual crop losses. As seen in virtually all insect species, serine proteases (SPs) are crucial to S. litura. The expression pattern of SPs from the midgut of S. litura was studied through expressed sequence tags (ESTs) analysis. One of SP (SlSP1) was chosen for detailed study, because the expression of the gene was midgut and larvae specific. SlSP1 was conducted as a model of its evolution, structure, and potential binding activity with corresponding substrates. SlSP1 is composed of 255 amino acids including a signal peptide at N-terminal followed by a putative activation peptide and the mature protein along with five putative phosphorylation sites, three disulphide bridges, and two N-glycosylation positions. At least nine conserved motifs were obtained in multiple sequence alignments. Some conserved residues, such as the catalytic triad His84, Asp127, and Ser229 as well as six cysteines at position 66, 82, 194, 211, 223, and 247, were examined. After homology modeling and molecular dynamics simulation, the resultant three-dimensional (3D) structure of SlSP1 was docked with the substrates 2PTC-Arg and 2PTC-Lys, respectively. Molecular Mechanic/PoissonBoltzmann surface area analysis was applied to anticipate optimal binding mode and crucial active sites of this enzyme. The residues Trp28, Gly187, Aso188, Arg249, Ile250, Lys246, and Lys278 are crucial for the substrate binding and molecule process. This information can be used in logical design of SPs inhibitors. New inhibitors may be a basis for development of a new pest control technology.

The date palm (Phoenix dactylifera L.), famed for its sugar-rich fruits (dates) and cultivated by humans since 4,000 B.C., is an economically important crop in the Middle East, Northern Africa, and increasingly other places where climates are suitable. Despite a long history of human cultivation, the understanding of P. dactylifera genetics and molecular biology are rather limited, hindered by lack of basic data in high quality from genomics and transcriptomics. Here we report a large-scale effort in generating gene models (assembled expressed sequence tags or ESTs and mapped to a genome assembly) for P. dactylifera, using the long-read pyrosequencing platform (Roche/454 GS FLX Titanium) in high coverage. We built fourteen cDNA libraries from different P. dactylifera tissues (cultivar Khalas) and acquired 15,778,993 raw sequencing reads—about one million sequencing reads per library—and the pooled sequences were assembled into 67,651 non-redundant contigs and 301,978 singletons. We annotated 52,725 contigs based on the plant databases and 45 contigs based on functional domains referencing to the Pfam database. From the annotated contigs, we assigned GO (Gene Ontology) terms to 36,086 contigs and KEGG pathways to 7,032 contigs. Our comparative analysis showed that 70.6 % (47,930), 69.4 % (47,089), 68.4 % (46,441), and 69.3 % (47,048) of the P. dactylifera gene models are shared with rice, sorghum, Arabidopsis, and grapevine, respectively. We also assigned our gene models into house-keeping and tissue-specific genes based on their tissue specificity. Electronic supplementary material The online version of this article (doi:10.1007/s11103-012-9924-z) contains supplementary material, which is available to authorized users.

Histone H3 lysine 4 trimethylation (H3K4me3) is well known to occur in the promoter region of genes for transcription activation. However, when investigating the H3K4me3 profiles in the mouse cerebrum and testis, we discovered that H3K4me3 also has a significant enrichment at the 3′ end of actively transcribed (sense) genes, named as 3′-H3K4me3. 3′-H3K4me3 is associated with ∼15% of protein-coding genes in both tissues. In addition, we examined the transcriptional initiation signals including RNA polymerase II (RNAPII) binding sites and 5′-CAGE-tag that marks transcriptional start sites. Interestingly, we found that 3′-H3K4me3 is associated with the initiation of antisense transcription. Furthermore, 3′-H3K4me3 modification levels correlate positively with the antisense expression levels of the associated sense genes, implying that 3′-H3K4me3 is involved in the activation of antisense transcription. Taken together, our findings suggest that H3K4me3 may be involved in the regulation of antisense transcription that initiates from the 3′ end of sense genes. In addition, a positive correlation was also observed between the expression of antisense and the associated sense genes with 3′-H3K4me3 modification. More importantly, we observed the 3′-H3K4me3 enrichment among genes in human, fruitfly and Arabidopsis, and found that the sequences of 3′-H3K4me3-marked regions are highly conserved and essentially indistinguishable from known promoters in vertebrate. Therefore, we speculate that these 3′-H3K4me3-marked regions may serve as potential promoters for antisense transcription and 3′-H3K4me3 appear to be a universal epigenetic feature in eukaryotes. Our results provide a novel insight into the epigenetic roles of H3K4me3 and the regulatory mechanism of antisense transcription.

Plant non-specific lipid transfer proteins (nsLtps) have been reported to be involved in plant defense activity against bacterial and fungal pathogens. In this study, we identified 135 (122 putative and 13 previously identified) Solanaceae nsLtps, which are clustered into 8 different groups. By comparing with Boutrot’s nsLtp classification, we classified these eight groups into five types (I, II, IV, IX and X). We compared Solanaceae nsLtps with Arabidopsis and Gramineae nsLtps and found that (1) Types I, II and IV are shared by Solanaceae, Gramineae and Arabidopsis; (2) Types III, V, VI and VIII are shared by Gramineae and Arabidopsis but not detected in Solanaceae so far; (3) Type VII is only found in Gramineae whereas type IX is present only in Arabidopsis and Solanaceae; (4) Type X is a new type that accounts for 52.59% Solanaceae nsLtps in our data, and has not been reported in any other plant so far. We further built and compared the three-dimensional structures of the eight groups, and found that the major functional diversification within the nsLtp family could be predated to the monocot/dicot divergence, and many gene duplications and sequence variations had happened in the nsLtp family after the monocot/dicot divergence, especially in Solanaceae.

Transcripts are expressed spatially and temporally and they are very complicated, precise and specific; however, most studies are focused on protein-coding related genes. Recently, massively parallel cDNA sequencing (RNA-seq) has emerged to be a new and promising tool for transcriptome research, and numbers of non-coding RNAs, especially lincRNAs, have been widely identified and well characterized as important regulators of diverse biological processes. In this study, we used ultra-deep RNA-seq data from 15 mouse tissues to study the diversity and dynamic of non-coding RNAs in mouse. Using our own criteria, we identified totally 16,249 non-coding genes (21,569 non-coding RNAs) in mouse. We annotated these non-coding RNAs by diverse properties and found non-coding RNAs are generally shorter, have fewer exons, express in lower level and are more strikingly tissue-specific compared with protein-coding genes. Moreover, these non-coding RNAs show significant enrichment with transcriptional initiation and elongation signals including histone modifications (H3K4me3, H3K27me3 and H3K36me3), RNAPII binding sites and CAGE tags. The gene set enrichment analysis (GSEA) result revealed several sets of lincRNAs associated with diverse biological processes such as immune effector process, muscle development and sexual reproduction. Taken together, this study provides a more comprehensive annotation of mouse non-coding RNAs and gives an opportunity for future functional and evolutionary study of mouse non-coding RNAs.

MicroRNAs (miRNAs) play crucial roles in multiple stages of plant development and regulate gene expression at posttranscriptional and translational levels. In this study, we first identified 238 conserved miRNAs in date palm (Phoenix dactylifera) based on a high-quality genome assembly and defined 78 fruit-development-associated (FDA) miRNAs, whose expression profiles are variable at different fruit development stages. Using experimental data, we subsequently detected 276 novel P. dactylifera-specific FDA miRNAs and predicted their targets. We also revealed that FDA miRNAs function mainly in regulating genes involved in starch/sucrose metabolisms and other carbon metabolic pathways; among them, 221 FDA miRNAs exhibit negative correlation with their corresponding targets, which suggests their direct regulatory roles on mRNA targets. Our data define a comprehensive set of conserved and novel FDA miRNAs along with their expression profiles, which provide a basis for further experimentation in assigning discrete functions of these miRNAs in P. dactylifera fruit development.

Date palm (Phoenix dactylifera L.) is a cultivated woody plant species with agricultural and economic importance. Here we report a genome assembly for an elite variety (Khalas), which is 605.4 Mb in size and covers >90% of the genome (~671 Mb) and >96% of its genes (~41,660 genes). Genomic sequence analysis demonstrates that P. dactylifera experienced a clear genome-wide duplication after either ancient whole genome duplications or massive segmental duplications. Genetic diversity analysis indicates that its stress resistance and sugar metabolism-related genes tend to be enriched in the chromosomal regions where the density of single-nucleotide polymorphisms is relatively low. Using transcriptomic data, we also illustrate the date palm’s unique sugar metabolism that underlies fruit development and ripening. Our large-scale genomic and transcriptomic data pave the way for further genomic studies not only on P. dactylifera but also other Arecaceae plants., The date palm is one of the most economically important plants of the palm family. Here, the authors present a high-quality genome assembly of the date palm Phoenix dactylifera, and reveal insights into the unique sugar metabolism underlying fruit ripening.

The high-throughput next-generation sequencing technologies provide an excellent opportunity for the detection of less-abundance transcripts that may not be identifiable by previously available techniques. Here, we report a discovery of thousands of novel transcripts (mostly non-coding RNAs) that are expressed in mouse cerebrum, testis, and embryonic stem (ES) cells, through an in-depth analysis of rmRNA-seq data. These transcripts show significant associations with transcriptional start and elongation signals. At the upstream of these transcripts we observed significant enrichment of histone marks (histone H3 lysine 4 trimethylation, H3K4me3), RNAPII binding sites, and cap analysis of gene expression tags that mark transcriptional start sites. Along the length of these transcripts, we also observed enrichment of histone H3 lysine 36 trimethylation (H3K36me3). Moreover, these transcripts show strong purifying selection in their genomic loci, exonic sequences, and promoter regions, implying functional constraints on the evolution of these transcripts. These results define a collection of novel transcripts in the mouse genome and indicate their potential functions in the mouse tissues and cells.

Riemerella anatipestifer is a well-described pathogen of waterfowl and other avian species which can cause a great loss to the poultry industry. Here we obtained the complete genome sequence of R. anatipestifer strain RA-GD, which was isolated from an infected duck in Guangzhou, China, and was cultivated in our laboratory.

The global features of H3K4 and H3K27 trimethylations (H3K4me3 and H3K27me3) have been well studied in recent years, but most of these studies were performed in mammalian cell lines. In this work, we generated the genome-wide maps of H3K4me3 and H3K27me3 of mouse cerebrum and testis using ChIP-seq and their high-coverage transcriptomes using ribominus RNA-seq with SOLiD technology. We examined the global patterns of H3K4me3 and H3K27me3 in both tissues and found that modifications are closely-associated with tissue-specific expression, function and development. Moreover, we revealed that H3K4me3 and H3K27me3 rarely occur in silent genes, which contradicts the findings in previous studies. Finally, we observed that bivalent domains, with both H3K4me3 and H3K27me3, existed ubiquitously in both tissues and demonstrated an invariable preference for the regulation of developmentally-related genes. However, the bivalent domains tend towards a “winner-takes-all” approach to regulate the expression of associated genes. We also verified the above results in mouse ES cells. As expected, the results in ES cells are consistent with those in cerebrum and testis. In conclusion, we present two very important findings. One is that H3K4me3 and H3K27me3 rarely occur in silent genes. The other is that bivalent domains may adopt a “winner-takes-all” principle to regulate gene expression.