Wanfei Liu

Genome sequence of the date palm Phoenix dactylifera L

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.4u2009Mb in size and covers >90% of the genome (~671u2009Mb) 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.

wapRNA: a web-based application for the processing of RNA sequences.

SUMMARYnmRNA/miRNA-seq technology is becoming the leading technology to globally profile gene expression and elucidate the transcriptional regulation mechanisms in living cells. Although there are many tools available for analyzing RNA-seq data, few of them are available as easy accessible online web tools for processing both mRNA and miRNA data for the RNA-seq based user community. As such, we have developed a comprehensive web application tool for processing mRNA-seq and miRNA-seq data. Our web tool wapRNA includes four different modules: mRNA-seq and miRNA-seq sequenced from SOLiD or Solexa platform and all the modules were tested on previously published experimental data. We accept raw sequence data with an optional reads filter, followed by mapping and gene annotation or miRNA prediction. wapRNA also integrates downstream functional analyses such as Gene Ontology, KEGG pathway, miRNA targets prediction and comparison of genes or miRNAs different expression in different samples. Moreover, we provide the executable packages for installation on users local server.nnnAVAILABILITYnwapRNA is freely available for use at http://waprna.big.ac.cn. The executable packages and the instruction for installation can be downloaded from our web site.nnnCONTACTnhusn@big.ac.cn; songshh@big.ac.cn.nnnSUPPLEMENTARY INFORMATIONnSupplementary data are available at Bioinformatics online.

Complete genome sequence of the pathogenic bacterium Riemerella anatipestifer strain RA-GD.

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.

Comparative analyses of H3K4 and H3K27 trimethylations between the mouse cerebrum and testis.

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.

Large-scale collection and annotation of gene models for date palm (Phoenix dactylifera, L.)

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.6xa0% (47,930), 69.4xa0% (47,089), 68.4xa0% (46,441), and 69.3xa0% (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.

Thousands of Novel Transcripts Identified in Mouse Cerebrum, Testis, and ES Cells Based on ribo-minus RNA Sequencing

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.

Discovery, Identification and Comparative Analysis of Non-Specific Lipid Transfer Protein (nsLtp) Family in Solanaceae

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.

Profiling microRNA expression during multi-staged date palm (Phoenix dactylifera L.) fruit development.

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.

The association between H3K4me3 and antisense transcription.

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.

A large-scale gene discovery for the red palm weevil Rhynchophorus ferrugineus （Coleoptera： Curculionidae）

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 26u2009765 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 60u2009000 single nucleotide polymorphisms and 1u2009200 simple sequence repeat markers. This study provides the first large‐scale cDNA dataset for RPW, a much‐needed resource for future molecular studies.