Crop genome sequencing: lessons and rationales

DNase I hypersensitive sites (DHSs) are hallmarks of chromatin zones containing transcriptional regulatory elements, making them critical in understanding the regulatory mechanisms of gene expression. Although large amounts of DHSs in the plant genome have been identified by high-throughput techniques, current DHSs obtained from experimental methods cover only a fraction of plant species and cell processes. Furthermore, these experimental methods are both time-consuming and expensive. Hence, it is urgent to develop automated computational means to efficiently and accurately predict DHSs in the plant genome. Recently, several methods have been proposed to predict the DHSs. However, all these methods took a lot of time to build the model, making them inappropriate for data with massive volume. In the present work, a new ensemble extreme learning machine (ELM)-based model called pDHS-ELM was proposed to predict the DHSs in the plant genome by fusing two different modes of pseudo-nucleotide composition. Here, two kinds of features including reverse complement kmer and pseudo-nucleotide composition were used to represent the DHSs. The ELM model was used to build the base classifiers. Then, an ensemble framework was employed to combine the outputs of these base classifiers. When applied to DHSs in Arabidopsis thaliana and rice (Oryza sativa) genome, the proposed method could obtain accuracies up to 88.48 and 87.58%, respectively. Compared with the state-of-the-art techniques, pDHS-ELM achieved higher sensitivity, specificity, and Matthew's correlation coefficient with much less training and test time. By employing pDHS-ELM, we identified 42,370 and 103,979 DHSs in A. thaliana and rice genome, respectively. The predicted DHSs were depleted of bulk nucleosomes and were tightly associated with transcription factors. Approximately 90% of the predicted DHSs were overlapped with transcription factors. Meanwhile, we demonstrated that the predicted DHSs were also associated with DNA methylation, nucleosome positioning/occupancy, and histone modification. This result suggests that pDHS-ELM can be considered as a new promising and powerful tool for transcriptional regulatory elements analysis. Our pDHS-ELM tool is available from the following website https://github.com/shanxinzhang/pDHS-ELM/ .

A "gene-island" sequencing strategy has been developed that expedites the targeted acquisition of orthologous gene sequences from related species for comparative genome analysis. A 152-kb bacterial artificial chromosome (BAC) clone from sorghum (Sorghum bicolor) encoding phytochrome A (PHYA) was fully sequenced, revealing 16 open reading frames with a gene density similar to many regions of the rice (Oryza sativa) genome. The sequences of genes in the orthologous region of the maize (Zea mays) and rice genomes were obtained using the gene-island sequencing method. BAC clones containing the orthologous maize and rice PHYA genes were identified, sheared, subcloned, and probed with the sorghum PHYA-containing BAC DNA. Sequence analysis revealed that approximately 75% of the cross-hybridizing subclones contained sequences orthologous to those within the sorghum PHYA BAC and less than 25% contained repetitive and/or BAC vector DNA sequences. The complete sequence of four genes, including up to 1 kb of their promoter regions, was identified in the maize PHYA BAC. Nine orthologous gene sequences were identified in the rice PHYA BAC. Sequence comparison of the orthologous sorghum and maize genes aided in the identification of exons and conserved regulatory sequences flanking each open reading frame. Within genomic regions where micro-colinearity of genes is absolutely conserved, gene-island sequencing is a particularly useful tool for comparative analysis of genomes between related species.

There is controversy as to whether gene expression is silenced in the functional centromere. The complete genomic sequences of the centromeric regions in higher eukaryotes have not been fully elucidated, because the presence of highly repetitive sequences complicates many aspects of genomic sequencing. We performed resequencing, assembly, and sequence finishing of two P1-derived artificial chromosome clones in the centromeric region of rice (Oryza sativa L.) chromosome 5 (Cen5). The pericentromeric region, where meiotic recombination is silenced, is located at the center of chromosome 5 and is 2.14 Mb long; a total of six restriction-fragment-length polymorphism markers (R448, C1388, S20487S, E3103S, C53260S, and R2059) genetically mapped at 54.6 cM were located in this region. In the pericentromeric region, 28 genes were annotated on the short arm and 45 genes on the long arm. To quantify all transcripts in this region, we performed massive parallel sequencing of mRNA. Transcriptional density (total length of transcribed region/length of the genomic region) and expression level (number of uniquely mapped reads/length of transcribed region) were calculated on the basis of the mapped reads on the rice genome. Transcriptional density and expression level were significantly lower in Cen5 than in the average of the other chromosomal regions. Moreover, transcriptional density in Cen5 was significantly lower on the short arm than on the long arm; the distribution of transcriptional density was asymmetric. The genomic sequence of Cen5 has been integrated into the most updated reference rice genome sequence constructed by the International Rice Genome Sequencing Project.

BackgroundSucrose nonfermenting-1 (SNF1)-related protein kinases (SnRKs) play important roles in regulating metabolism and stress responses in plants, providing a conduit for crosstalk between metabolic and stress signalling, in some cases involving the stress hormone, abscisic acid (ABA). The burgeoning and divergence of the plant gene family has led to the evolution of three subfamilies, SnRK1, SnRK2 and SnRK3, of which SnRK2 and SnRK3 are unique to plants. Therefore, the study of SnRKs in crops may lead to the development of strategies for breeding crop varieties that are more resilient under stress conditions. In the present study, we describe the SnRK gene family of barley (Hordeum vulgare), the widespread cultivation of which can be attributed to its good adaptation to different environments.ResultsThe barley HvSnRK gene family was elucidated in its entirety from publicly-available genome data and found to comprise 50 genes. Phylogenetic analyses assigned six of the genes to the HvSnRK1 subfamily, 10 to HvSnRK2 and 34 to HvSnRK3. The search was validated by applying it to Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa) genome data, identifying 50 SnRK genes in rice (four OsSnRK1, 11 OsSnRK2 and 35 OsSnRK3) and 39 in Arabidopsis (three AtSnRK1, 10 AtSnRK2 and 26 AtSnRK3). Specific motifs were identified in the encoded barley proteins, and multiple putative regulatory elements were found in the gene promoters, with light-regulated elements (LRE), ABA response elements (ABRE) and methyl jasmonate response elements (MeJa) the most common. RNA-seq analysis showed that many of the HvSnRK genes responded to ABA, some positively, some negatively and some with complex time-dependent responses.ConclusionsThe barley HvSnRK gene family is large, comprising 50 members, subdivided into HvSnRK1 (6 members), HvSnRK2 (10 members) and HvSnRK3 (34 members), showing differential positive and negative responses to ABA.

Plant chitinases, a class of glycosyl hydrolases, participate in various aspects of normal plant growth and development, including cell wall metabolism and disease resistance. The rice (Oryza sativa) genome encodes 37 putative chitinases and chitinase-like proteins. However, none of them has been characterized at the genetic level. In this study, we report the isolation of a brittle culm mutant, bc15, and the map-based cloning of the BC15/OsCTL1 (for chitinase-like1) gene affected in the mutant. The gene encodes the rice chitinase-like protein BC15/OsCTL1. Mutation of BC15/OsCTL1 causes reduced cellulose content and mechanical strength without obvious alterations in plant growth. Bioinformatic analyses indicated that BC15/OsCTL1 is a class II chitinase-like protein that is devoid of both an amino-terminal cysteine-rich domain and the chitinase activity motif H-E-T-T but possesses an amino-terminal transmembrane domain. Biochemical assays demonstrated that BC15/OsCTL1 is a Golgi-localized type II membrane protein that lacks classical chitinase activity. Quantitative real-time polymerase chain reaction and β-glucuronidase activity analyses indicated that BC15/OsCTL1 is ubiquitously expressed. Investigation of the global expression profile of wild-type and bc15 plants, using Illumina RNA sequencing, further suggested a possible mechanism by which BC15/OsCTL1 mediates cellulose biosynthesis and cell wall remodeling. Our findings provide genetic evidence of a role for plant chitinases in cellulose biosynthesis in rice, which appears to differ from their roles as revealed by analysis of Arabidopsis (Arabidopsis thaliana).

Dasheng is one of the highest copy number long terminal repeat elements and one of the most recent elements to amplify in the rice (Oryza sativa) genome. However, the absence of any significant coding capacity for retroviral proteins, including gag and pol, suggests that Dasheng is a nonautonomous element. Here, we have exploited the availability of 360 Mb of rice genomic sequence to identify a candidate autonomous element. RIRE2 is a previously described gypsy-like long terminal repeat retrotransposon with significant sequence similarity to Dasheng in the regions where putative cis factors for retrotransposition are thought to be located. Dasheng and RIRE2 elements have similar chromosomal distribution patterns and similar target site sequences, suggesting that they use the same transposition machinery. In addition, the presence of several RIRE2-Dasheng element chimeras in the genome is consistent with the copackaging of element mRNAs in the same virus-like particle. Finally, both families have recently amplified members, suggesting that they could have been co-expressed, a necessary prerequisite for RIRE2 to serve as the source of transposition machinery for Dasheng. Consistent with this hypothesis, transcripts from both elements were found in the same expressed sequence tag library.

BackgroundNon-coding RNA (ncRNA) genes do not encode proteins but produce functional RNA molecules that play crucial roles in many key biological processes. Recent genome-wide transcriptional profiling studies using tiling arrays in organisms such as human and Arabidopsis have revealed a great number of transcripts, a large portion of which have little or no capability to encode proteins. This unexpected finding suggests that the currently known repertoire of ncRNAs may only represent a small fraction of ncRNAs of the organisms. Thus, efficient and effective prediction of ncRNAs has become an important task in bioinformatics in recent years. Among the available computational methods, the comparative genomic approach seems to be the most powerful to detect ncRNAs. The recent completion of the sequencing of several major plant genomes has made the approach possible for plants.ResultsWe have developed a pipeline to predict novel ncRNAs in the Arabidopsis (Arabidopsis thaliana) genome. It starts by comparing the expressed intergenic regions of Arabidopsis as provided in two whole-genome high-density oligo-probe arrays from the literature with the intergenic nucleotide sequences of all completely sequenced plant genomes including rice (Oryza sativa), poplar (Populus trichocarpa), grape (Vitis vinifera), and papaya (Carica papaya). By using multiple sequence alignment, a popular ncRNA prediction program (RNAz), wet-bench experimental validation, protein-coding potential analysis, and stringent screening against various ncRNA databases, the pipeline resulted in 16 families of novel ncRNAs (with a total of 21 ncRNAs).ConclusionIn this paper, we undertake a genome-wide search for novel ncRNAs in the genome of Arabidopsis by a comparative genomics approach. The identified novel ncRNAs are evolutionarily conserved between Arabidopsis and other recently sequenced plants, and may conduct interesting novel biological functions.

The Euonymus lectin (EUL) family is a unique group of carbohydrate-binding proteins that is omnipresent in plants. Sequences encoding EUL-related lectins have been retrieved from all completely sequenced plant genomes. The rice (Oryza sativa) genome contains 5 functional EUL genes referred to as OsEULS2, OsEULS3, OsEULD1a, OsEULD1b, and OsEULD2. In this study we focused on the tissue specific expression, stress inducibility and subcellular localization of the rice EULs. Even though the EUL domain sequence is highly conserved among the rice EULs (at least 80% sequence similarity) different biotic and abiotic stress treatments yielded unique responses for the different EULs. Transcript levels for OsEULs were differentially affected by drought and salt stress, ABA treatment, pathogen infection or insect infestation. Analysis of promoter activity revealed differential expression and tissue specificity for the 5 OsEUL genes, with most expression observed in the vascular system of roots and shoots, as well as in the root tips and seeds. At cell level, all OsEULs are located in the nucleus whereas OsEULD1b and OsEULD2 also locate to the cytoplasm. This paper contributes to the functional characterization of the EULs and provides insight in the biological importance of this family of proteins for rice.

The maize (Zea mays L.) CRINKLY4 (CR4) gene encodes a serine/threonine receptor-like kinase that controls an array of developmental processes in the plant and endosperm. The Arabidopsis thaliana (L.) Heynh. genome encodes an ortholog of CR4, ACR4, and four CRINKLY4-RELATED (CRR) proteins: AtCRRl, AtCRR2, AtCRR3 and AtCRKl. The available genome sequence of rice (Oryza sativa L.) encodes a CR4 ortholog, OsCR4, and four CRR proteins: OsCRRl, OsCRR2, OsCRR3 and OsCRR4, not necessarily orthologous to the Arabidopsis CRRs. A phylogenetic study showed that AtCRRl and AtCRR2 form a

The rice (Oryza sativa) genome harbours three genes encoding CysCysHisCys (CCHC)-type zinc finger-containing glycine-rich RNA-binding proteins, designated OsRZ proteins, but their importance and physiological functions remain largely unknown. Here, the stress-responsive expression patterns of OsRZs were assessed, and the biological and cellular functions of OsRZs were evaluated under low temperature conditions. The expression levels of the three OsRZs were up-regulated by cold stress, whereas drought or high salt stress did not significantly alter its transcript level. OsRZ2 complemented the cold sensitivity of BX04 Escherichia coli cells under low temperatures, and had DNA-melting activity and transcription anti-termination activity, thereby indicating that OsRZ2 possesses an RNA chaperone activity. By contrast, neither OsRZ1 nor OsRZ3 harboured these activities. Ectopic expression of OsRZ2, but not OsRZ3, in cold-sensitive Arabidopsis grp7 knockout plants rescued the grp7 plants from cold and freezing damage, and OsRZ2 complemented the defect in mRNA export from the nucleus to the cytoplasm in grp7 mutant during cold stress. The present findings support the emerging idea that the regulation of mRNA export is one of the adaptive processes in plants under stress conditions, and RNA chaperone functions as a regulator in mRNA export under cold stress conditions.

Sequencing and computational annotation revealed several features, including high gene numbers, unusual composition of the predicted genes and a large number of genes lacking homology to known genes, that distinguish the rice (Oryza sativa) genome from that of other fully sequenced model species. We report here a full-genome transcription analysis of the indica rice subspecies using high-density oligonucleotide tiling microarrays. Our results provided expression data support for the existence of 35,970 (81.9%) annotated gene models and identified 5,464 unique transcribed intergenic regions that share similar compositional properties with the annotated exons and have significant homology to other plant proteins. Elucidating and mapping of all transcribed regions revealed an association between global transcription and cytological chromosome features, and an overall similarity of transcriptional activity between duplicated segments of the genome. Collectively, our results provide the first whole-genome transcription map useful for further understanding the rice genome.

Nucleotide landscapes, which are the way base composition is distributed along a genome, strongly vary among species. The underlying causes of these variations have been much debated. Though mutational bias and selection were initially invoked, GC-biased gene conversion (gBGC), a recombination-associated process favoring the G and C over A and T bases, is increasingly recognized as a major factor. As opposed to vertebrates, evolution of GC content is less well known in plants. Most studies have focused on the GC-poor and homogeneous Arabidopsis thaliana genome and the much more GC-rich and heterogeneous rice (Oryza sativa) genome and have often been generalized as a dicot/monocot dichotomy. This vision is clearly phylogenetically biased and does not allow understanding the mechanisms involved in GC content evolution in plants. To tackle these issues, we used EST data from more than 200 species and provided the most comprehensive description of gene GC content across the seed plant phylogeny so far available. As opposed to the classically assumed dicot/monocot dichotomy, we found continuous variations in GC content from the probably ancestral GC-poor and homogeneous genomes to the more derived GC-rich and highly heterogeneous ones, with several independent enrichment episodes. Our results suggest that gBGC could play a significant role in the evolution of GC content in plant genomes.

Iron-sulphur (Fe-S) cluster assembly occurs in chloroplasts, mitochondria and cytosol, involving dozens of genes in higher plants. In this study, we have identified 41 putative Fe-S cluster assembly genes in rice (Oryza sativa) genome, and the expression of all genes was verified. To investigate the role of Fe-S cluster assembly as a metabolic pathway, we applied abiotic stresses to rice seedlings and analysed Fe-S cluster assembly gene expression by qRT-PCR. Our data showed that genes for Fe-S cluster assembly in chloroplasts of leaves are particularly sensitive to heavy metal treatments, and that Fe-S cluster assembly genes in roots were up-regulated in response to iron toxicity, oxidative stress and some heavy metal assault. The effect of each stress treatment on the Fe-S cluster assembly machinery demonstrated an unexpected tissue or organelle specificity, suggesting that the physiological relevance of the Fe-S cluster assembly is more complex than thought. Furthermore, our results may reveal potential candidate genes for molecular breeding of rice.

The comprehensive annotation of regulatory elements in linear genomes is needed to elucidate the molecular mechanisms underlying chromatin loop formation in plants. Here, we characterized a novel family of conserved noncoding sequences (CNSs) in the rice (Oryza sativa) genome. These sequences, known as AT-rich pincer-like elements (APEs), are composed of 13-bp repeat unit arrays in a reverse-forward configuration. Our findings revealed that there are 611 APE copies across the japonica genome. Deletion of single APEs disrupted the long-range chromatin loops anchoring target-APE regions and moderately remodeled the profile of A/B compartments, topologically associating domains (TADs), and chromatin loops, thereby rewiring the expression of looped gene(s) including those controlling important agronomic traits. Thus, APEs function as hub motifs directly mediating chromatin looping and maintaining 3D genome integrity and stability at the levels of compartments, TADs, and loops. Moreover, neighboring genomic regions harboring numerous paired non-APE (NA) CNSs were more likely to interact with each other. This finding suggests that NA CNS pairs might play a helper role in determining loop frequency in a dose-dependent manner, likely by ensuring the pairing selectivity of anchor sites. Our study highlights the importance of APEs and NA CNSs in maintaining 3D genome structure, thereby providing the framework required to link many noncoding repetitive elements to their molecular functions in plants.

Posttranslational glycosylphosphatidylinositol (GPI) lipid anchoring is common not only for animal and fungal but also for plant proteins. The attachment of the GPI moiety to the carboxyl-terminus after proteolytic cleavage of a C-terminal propeptide is performed by the transamidase complex. Its four known subunits also have obvious full-length orthologs in the Arabidopsis and rice (Oryza sativa) genomes; thus, the mechanism of substrate protein processing appears similar for all eukaryotes. A learning set of plant proteins (substrates for the transamidase complex) has been collected both from the literature and plant sequence databases. We find that the plant GPI lipid anchor motif differs in minor aspects from the animal signal (e.g. the plant hydrophobic tail region can contain a higher fraction of aromatic residues). We have developed the "big-Pi plant" program for prediction of compatibility of query protein C-termini with the plant GPI lipid anchor motif requirements. Validation tests show that the sensitivity for transamidase targets is approximately 94%, and the rate of false positive prediction is about 0.1%. Thus, the big-Pi predictor can be applied as unsupervised genome annotation and target selection tool. The program is also suited for the design of modified protein constructs to test their GPI lipid anchoring capacity. The big-Pi plant predictor Web server and lists of potential plant precursor proteins in Swiss-Prot, SPTrEMBL, Arabidopsis, and rice proteomes are available at http://mendel.imp.univie.ac.at/gpi/plants/gpi_plants.html. Arabidopsis and rice protein hits have been functionally classified. Several GPI lipid-anchored arabinogalactan-related proteins have been identified in rice.