MIPS PlantsDB: a database framework for comparative plant genome research

Gossypium barbadense and Gossypium hirsutum genomes provide insights into the origin and evolution of allotetraploid cotton

Allotetraploid cotton is an economically important natural-fiber-producing crop worldwide. After polyploidization, Gossypium hirsutum L. evolved to produce a higher fiber yield and to better survive harsh environments than Gossypium barbadense, which produces superior-quality fibers. The global genetic and molecular bases for these interspecies divergences were unknown. Here we report high-quality de novo-assembled genomes for these two cultivated allotetraploid species with pronounced improvement in repetitive-DNA-enriched centromeric regions. Whole-genome comparative analyses revealed that species-specific alterations in gene expression, structural variations and expanded gene families were responsible for speciation and the evolutionary history of these species. These findings help to elucidate the evolution of cotton genomes and their domestication history. The information generated not only should enable breeders to improve fiber quality and resilience to ever-changing environmental conditions but also can be translated to other crops for better understanding of their domestication history and use in improvement.

Exploring miRNAs for developing climate-resilient crops: A perspective review

Climate changes and environmental stresses have significant implications on global crop production and necessitate developing crops that can withstand an array of climate changes and environmental perturbations such as irregular water-supplies leading to drought or water-logging, hyper soil-salinity, extreme and variable temperatures, ultraviolet radiations and metal stress. Plants have intricate molecular mechanisms to cope with these dynamic environmental changes, one of the most common and effective being the reprogramming of expression of stress-responsive genes. Plant microRNAs (miRNAs) have emerged as key post-transcriptional and translational regulators of gene-expression for modulation of stress implications. Recent reports are establishing their key roles in epigenetic regulations of stress/adaptive responses as well as in providing plants genome-stability. Several stress responsive miRNAs are being identified from different crop plants and miRNA-driven RNA-interference (RNAi) is turning into a technology of choice for improving crop traits and providing phenotypic plasticity in challenging environments. Here we presents a perspective review on exploration of miRNAs as potent targets for engineering crops that can withstand multi-stress environments via loss-/gain-of-function approaches. This review also shed a light on potential roles plant miRNAs play in genome-stability and their emergence as potent target for genome-editing. Current knowledge on plant miRNAs, their biogenesis, function, their targets, and latest developments in bioinformatics approaches for plant miRNAs are discussed. Though there are recent reviews discussing primarily the individual miRNAs responsive to single stress factors, however, considering practical limitation of this approach, special emphasis is given in this review on miRNAs involved in responses and adaptation of plants to multi-stress environments including at epigenetic and/or epigenomic levels.

RepetDB: a unified resource for transposable element references

Background

Thanks to their ability to move around and replicate within genomes, transposable elements (TEs) are perhaps the most important contributors to genome plasticity and evolution. Their detection and annotation are considered essential in any genome sequencing project. The number of fully sequenced genomes is rapidly increasing with improvements in high-throughput sequencing technologies. A fully automated de novo annotation process for TEs is therefore required to cope with the deluge of sequence data.

However, all automated procedures are error-prone, and an automated procedure for TE identification and classification would be no exception. It is therefore crucial to provide not only the TE reference sequences, but also evidence justifying their classification, at the scale of the whole genome. A few TE databases already exist, but none provides evidence to justify TE classification. Moreover, biological information about the sequences remains globally poor.

Results

We present here the RepetDB database developed in the framework of GnpIS, a genetic and genomic information system. RepetDB is designed to store and retrieve detected, classified and annotated TEs in a standardized manner. RepetDB is an implementation with extensions of InterMine, an open-source data warehouse framework used here to store, search, browse, analyze and compare all the data recorded for each TE reference sequence. InterMine can display diverse information for each sequence and allows simple to very complex queries. Finally, TE data are displayed via a worldwide data discovery portal. RepetDB is accessible at urgi.versailles.inra.fr/repetdb.

Conclusions

RepetDB is designed to be a TE knowledge base populated with full de novo TE annotations of complete (or near-complete) genome sequences. Indeed, the description and classification of TEs facilitates the exploration of specific TE families, superfamilies or orders across a large range of species. It also makes possible cross-species searches and comparisons of TE family content between genomes.

Domesticating Vigna Stipulacea: A Potential Legume Crop With Broad Resistance to Biotic Stresses

Though crossing wild relatives to modern cultivars is a usual means to introduce alleles of stress tolerance, an alternative is de novo domesticating wild species that are already tolerant to various kinds of stresses. As a test case, we chose Vigna stipulacea Kuntze, which has fast growth, short vegetative stage, and broad resistance to pests and diseases. We developed an ethyl methanesulfonate-mutagenized population and obtained three mutants with reduced seed dormancy and one with reduced pod shattering. We crossed one of the mutants of less seed dormancy to the wild type and confirmed that the phenotype was inherited in a Mendelian manner. De novo assembly of V. stipulacea genome, and the following resequencing of the F2 progenies successfully identified a Single Nucleotide Polymorphism (SNP) associated with seed dormancy. By crossing and pyramiding the mutant phenotypes, we will be able to turn V. stipulacea into a crop which is yet primitive but can be cultivated without pesticides.

Chromosome-based survey sequencing reveals the genome organization of wild wheat progenitor Triticum dicoccoides

Wild emmer wheat (Triticum turgidum ssp. dicoccoides) is the progenitor of wheat. We performed chromosome-based survey sequencing of the 14 chromosomes, examining repetitive sequences, protein-coding genes, miRNA/target pairs and tRNA genes, as well as syntenic relationships with related grasses. We found considerable differences in the content and distribution of repetitive sequences between the A and B subgenomes. The gene contents of individual chromosomes varied widely, not necessarily correlating with chromosome size. We catalogued candidate agronomically important loci, along with new alleles and flanking sequences that can be used to design exome sequencing. Syntenic relationships and virtual gene orders revealed several small-scale evolutionary rearrangements, in addition to providing evidence for the 4AL-5AL-7BS translocation in wild emmer wheat. Chromosome-based sequence assemblies contained five novel miRNA families, among 59 families putatively encoded in the entire genome which provide insight into the domestication of wheat and an overview of the genome content and organization.

The chromosome-level quality genome provides insights into the evolution of the biosynthesis genes for aroma compounds of Osmanthus fragrans

Sweet osmanthus (Osmanthus fragrans) is a very popular ornamental tree species throughout Southeast Asia and USA particularly for its extremely fragrant aroma. We constructed a chromosome-level reference genome of O. fragrans to assist in studies of the evolution, genetic diversity, and molecular mechanism of aroma development. A total of over 118 Gb of polished reads was produced from HiSeq (45.1 Gb) and PacBio Sequel (73.35 Gb), giving 100× depth coverage for long reads. The combination of Illumina-short reads, PacBio-long reads, and Hi-C data produced the final chromosome quality genome of O. fragrans with a genome size of 727 Mb and a heterozygosity of 1.45 %. The genome was annotated using de novo and homology comparison and further refined with transcriptome data. The genome of O. fragrans was predicted to have 45,542 genes, of which 95.68 % were functionally annotated. Genome annotation found 49.35 % as the repetitive sequences, with long terminal repeats (LTR) being the richest (28.94 %). Genome evolution analysis indicated the evidence of whole-genome duplication 15 million years ago, which contributed to the current content of 45,242 genes. Metabolic analysis revealed that linalool, a monoterpene is the main aroma compound. Based on the genome and transcriptome, we further demonstrated the direct connection between terpene synthases (TPSs) and the rich aromatic molecules in O. fragrans. We identified three new flower-specific TPS genes, of which the expression coincided with the production of linalool. Our results suggest that the high number of TPS genes and the flower tissue- and stage-specific TPS genes expressions might drive the strong unique aroma production of O. fragrans.

LncRNAs in polyploid cotton interspecific hybrids are derived from transposon neofunctionalization

BACKGROUND

Interspecific hybridization and whole genome duplication are driving forces of genomic and organism diversification. But the effect of interspecific hybridization and whole genome duplication on the non-coding portion of the genome in particular remains largely unknown. In this study, we examine the profile of long non-coding RNAs (lncRNAs), comparing them with that of coding genes in allotetraploid cotton (Gossypium hirsutum), its putative diploid ancestors (G. arboreum; G. raimondii), and an F1 hybrid (G. arboreum × G. raimondii, AD).

RESULTS

We find that most lncRNAs (80%) that were allelic expressed in the allotetraploid genome. Moreover, the genome shock of hybridization reprograms the non-coding transcriptome in the F1 hybrid. Interestingly, the activated lncRNAs are predominantly transcribed from demethylated TE regions, especially from long interspersed nuclear elements (LINEs). The DNA methylation dynamics in the interspecies hybridization are predominantly associated with the drastic expression variation of lncRNAs. Similar trends of lncRNA bursting are also observed in the progress of polyploidization. Additionally, we find that a representative novel lncRNA XLOC_409583 activated after polyploidization from a LINE in the A subgenome of allotetraploid cotton was involved in control of cotton seedling height.

CONCLUSION

Our results reveal that the processes of hybridization and polyploidization enable the neofunctionalization of lncRNA transcripts, acting as important sources of increased plasticity for plants.

Genomic Analysis and Prediction within a US Public Collaborative Winter Wheat Regional Testing Nursery

The development of inexpensive, whole-genome profiling enables a transition to allele-based breeding using genomic prediction models. These models consider alleles shared between lines to predict phenotypes and select new lines based on estimated breeding values. This approach can leverage highly unbalanced datasets that are common to breeding programs. The Southern Regional Performance Nursery (SRPN) is a public nursery established by the USDA-ARS in 1931 to characterize performance and quality of near-release wheat ( L.) varieties from breeding programs in the US Central Plains. New entries are submitted annually and can be re-entered only once. The trial is grown at >30 locations each year and lines are evaluated for grain yield, disease resistance, and agronomic traits. Overall genetic gain is measured across years by including common check cultivars for comparison. We have generated whole-genome profiles via genotyping-by-sequencing (GBS) for 939 SPRN entries dating back to 1992 to explore the potential use of the nursery as a genomic selection (GS) training population (TP). The GS prediction models across years (average = 0.33) outperformed year-to-year phenotypic correlation for yield ( = 0.27) for a majority of the years evaluated, suggesting that genomic selection has the potential to outperform low heritability selection on yield in these highly variable environments. We also examined the predictability of programs using both program-specific and whole-set TPs. Generally, the predictability of a program was similar with both approaches. These results suggest that wheat breeding programs can collaboratively leverage the immense datasets that are generated from regional testing networks.

Comparative Analysis of Homologous Sequences of Saccharum officinarum and Saccharum spontaneum Reveals Independent Polyploidization Events

Sugarcane (Saccharum spp. hybrids) is an economically important crop widely grown in tropical and subtropical regions for sugar and ethanol production. However, the large genome size, high ploidy level, interspecific hybridization and aneuploidy make sugarcane one of the most complex genomes and have long hampered genome research in sugarcane. Modern sugarcane cultivars are derived from interspecific hybridization between S. officinarum and S. spontaneum with 80-90% of the genome from S. officinarum and 10-20% of the genome from S. spontaneum. We constructed bacterial artificial chromosome (BAC) libraries of S. officinarum variety LA Purple (2n = 8x = 80) and S. spontaneum haploid clone AP85-441 (2n = 4x = 32), and selected and sequenced 97 BAC clones from the two Saccharum BAC libraries. A total of 5,847,280 bp sequence from S. officinarum and 5,011,570 bp from S. spontaneum were assembled and 749 gene models were annotated in these BACs. A relatively higher gene density and lower repeat content were observed in S. spontaneum BACs than in S. officinarum BACs. Comparative analysis of syntenic regions revealed a high degree of collinearity in genic regions between Saccharum and Sorghum bicolor and between S. officinarum and S. spontaneum. In the syntenic regions, S. spontaneum showed expansion relative to S. officinarum, and both S. officinarum and S. spontaneum showed expansion relative to sorghum. Among the 75 full-length LTR retrotransposons identified in the Saccharum BACs, none of them are older than 2.6 mys and no full-length LTR elements are shared between S. officinarum and S. spontaneum. In addition, divergence time estimated using a LTR junction marker and a syntenic gene shared by 3 S. officinarum and 1 S. spontaneum BACs revealed that the S. spontaneum intergenic region was distant to those from the 3 homologous regions in S. officinarum. Our results suggested that S. officinarum and S. spontaneum experienced at least two rounds of independent polyploidization in each lineage after their divergence from a common ancestor.

Draft genome sequence of wild Prunus yedoensis reveals massive inter-specific hybridization between sympatric flowering cherries

BACKGROUND

Hybridization is an important evolutionary process that results in increased plant diversity. Flowering Prunus includes popular cherry species that are appreciated worldwide for their flowers. The ornamental characteristics were acquired both naturally and through artificially hybridizing species with heterozygous genomes. Therefore, the genome of hybrid flowering Prunus presents important challenges both in plant genomics and evolutionary biology.

RESULTS

We use long reads to sequence and analyze the highly heterozygous genome of wild Prunus yedoensis. The genome assembly covers > 93% of the gene space; annotation identified 41,294 protein-coding genes. Comparative analysis of the genome with 16 accessions of six related taxa shows that 41% of the genes were assigned into the maternal or paternal state. This indicates that wild P. yedoensis is an F1 hybrid originating from a cross between maternal P. pendula f. ascendens and paternal P. jamasakura, and it can be clearly distinguished from its confusing taxon, Yoshino cherry. A focused analysis of the S-locus haplotypes of closely related taxa distributed in a sympatric natural habitat suggests that reduced restriction of inter-specific hybridization due to strong gametophytic self-incompatibility is likely to promote complex hybridization of wild Prunus species and the development of a hybrid swarm.

CONCLUSIONS

We report the draft genome assembly of a natural hybrid Prunus species using long-read sequencing and sequence phasing. Based on a comprehensive comparative genome analysis with related taxa, it appears that cross-species hybridization in sympatric habitats is an ongoing process that facilitates the diversification of flowering Prunus.

Molecular characterization and differential expression suggested diverse functions of P-type II Ca2+ATPases in Triticum aestivum L

Background

Plant P-type II Ca²⁺ATPases are formed by two distinct groups of proteins (ACAs and ECAs) that perform pumping of Ca²⁺ outside the cytoplasm during homeostasis, and play vital functions during development and stress management. In the present study, we have performed identification and characterisation of P-type II Ca²⁺ATPase gene family in an important crop plant Triticum aestivum.

Results

Herein, a total of 33 TaACA and 9 TaECA proteins were identified from the various chromosomes and sub-genomes of Triticum aestivum. Phylogenetic analysis revealed clustering of the homoeologous TaACA and TaECA proteins into 11 and 3 distinct groups that exhibited high sequence homology and comparable structural organization as well. Both TaACA and TaECA group proteins consisted of eight to ten transmembrane regions, and their respective domains and motifs. Prediction of sub-cellular localization was found variable for most of the proteins; moreover, it was consistent with the evolutionarily related proteins from rice and Arabidopsis in certain cases. The occurrence of assorted sets of cis-regulatory elements indicated their diverse functions. The differential expression of various TaACA and TaECA genes during developmental stages suggested their roles in growth and development. The modulated expression during heat, drought, salt and biotic stresses along with the occurrence of various stress specific cis-regulatory elements suggested their association with stress response. Interaction of these genes with numerous development and stress related genes indicated their decisive role in various biological processes and signaling.

Conclusion

T. aestivum genome consisted of a maximum of 42 P-type II Ca²⁺ATPase genes, derived from each A, B and D sub-genome. These genes may play diverse functions during plant growth and development. They may also be involved in signalling during abiotic and biotic stresses. The present study provides a comprehensive insight into the role of P-type II Ca²⁺ATPase genes in T. aestivum. However, the specific function of each gene needs to be established, which could be utilized in future crop improvement programs.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4792-9) contains supplementary material, which is available to authorized users.

A comprehensive review of web-based resources of non-coding RNAs for plant science research

Non-coding RNAs (ncRNAs) are transcribed from genome but not translated into proteins. Many ncRNAs are key regulators of plants growth and development, metabolism and stress tolerance. In order to make the web-based ncRNA resources for plant science research be more easily accessible and understandable, we made a comprehensive review for 83 web-based resources of three types, including genome databases containing ncRNA data, microRNA (miRNA) databases and long non-coding RNA (lncRNA) databases. To facilitate effective usage of these resources, we also suggested some preferred resources of miRNAs and lncRNAs for performing meaningful analysis.

De novo genome assembly of the red silk cotton tree (Bombax ceiba)

Background

Bombax ceiba L. (the red silk cotton tree) is a large deciduous tree that is distributed in tropical and sub-tropical Asia as well as northern Australia. It has great economic and ecological importance, with several applications in industry and traditional medicine in many Asian countries. To facilitate further utilization of this plant resource, we present here the draft genome sequence for B. ceiba.

Findings

We assembled a relatively intact genome of B. ceiba by using PacBio single-molecule sequencing and BioNano optical mapping technologies. The final draft genome is approximately 895 Mb long, with contig and scaffold N50 sizes of 1.0 Mb and 2.06 Mb, respectively.

Conclusions

The high-quality draft genome assembly of B. ceiba will be a valuable resource enabling further genetic improvement and more effective use of this tree species.

Dehydration Stress Contributes to the Enhancement of Plant Defense Response and Mite Performance on Barley

Under natural conditions, plants suffer different stresses simultaneously or in a sequential way. At present, the combined effect of biotic and abiotic stressors is one of the most important threats to crop production. Understanding how plants deal with the panoply of potential stresses affecting them is crucial to develop biotechnological tools to protect plants. As well as for drought stress, the economic importance of the spider mite on agriculture is expected to increase due to climate change. Barley is a host of the polyphagous spider mite Tetranychus urticae and drought produces important yield losses. To obtain insights on the combined effect of drought and mite stresses on the defensive response of this cereal, we have analyzed the transcriptomic responses of barley plants subjected to dehydration (water-deficit) treatment, spider mite attack, or to the combined dehydration-spider mite stress. The expression patterns of mite-induced responsive genes included many jasmonic acid responsive genes and were quickly induced. In contrast, genes related to dehydration tolerance were later up-regulated. Besides, a higher up-regulation of mite-induced defenses was showed by the combined dehydration and mite treatment than by the individual mite stress. On the other hand, the performance of the mite in dehydration stressed and well-watered plants was tested. Despite the stronger defensive response in plants that suffer dehydration and mite stresses, the spider mite demonstrates a better performance under dehydration condition than in well-watered plants. These results highlight the complexity of the regulatory events leading to the response to a combination of stresses and emphasize the difficulties to predict their consequences on crop production.

Gene architecture and expression analyses provide insights into the role of glutathione peroxidases (GPXs) in bread wheat (Triticum aestivum L.)

Glutathione peroxidases (GPXs) are redox sensor proteins that maintain a steady-state of H₂O₂ in plant cells. They exhibit distinct sub-cellular localization and have diverse functionality in response to different stimuli. In this study, a total of 14 TaGPX genes and three splice variants were identified in the genome of Triticum aestivum and evaluated for various physicochemical properties. The TaGPX genes were scattered on the various chromosomes of the A, B, and D sub-genomes and clustered into five homeologous groups based on high sequence homology. The majority of genes were derived from the B sub-genome and localized on chromosome 2. The intron-exon organization, motif and domain architecture, and phylogenetic analyses revealed the conserved nature of TaGPXs. The occurrence of both development-related and stress-responsive cis-acting elements in the promoter region, the differential expression of these genes during various developmental stages, and the modulation of expression in the presence of biotic and abiotic stresses suggested their diverse role in T. aestivum. The majority of TaGPX genes showed higher expression in various leaf developmental stages. However, TaGPX1-A1 was upregulated in the presence of each abiotic stress treatment. A co-expression analysis revealed the interaction of TaGPXs with numerous development and stress-related genes, which indicated their vital role in numerous biological processes. Our study revealed the opportunities for further characterization of individual TaGPX proteins, which might be useful in designing future crop improvement strategies.

Overexpression of HvIcy6 in Barley Enhances Resistance against Tetranychus urticae and Entails Partial Transcriptomic Reprogramming

Cystatins have been largely used for pest control against phytophagous species. However, cystatins have not been commonly overexpressed in its cognate plant species to test their pesticide capacity. Since the inhibitory role of barley HvCPI-6 cystatin against the phytophagous mite Tetranychus urticae has been previously demonstrated, the purpose of our study was to determine if barley transgenic lines overexpressing its own HvIcy6 gene were more resistant against this phytophagous infestation. Besides, a transcriptomic analysis was done to find differential expressed genes among wild-type and transformed barley plants. Barley plants overexpressing HvIcy6 cystatin gene remained less susceptible to T. urticae attack when compared to wild-type plants, with a significant lesser foliar damaged area and a lower presence of the mite. Transcriptomic analysis revealed a certain reprogramming of cellular metabolism and a lower expression of several genes related to photosynthetic activity. Therefore, although caution should be taken to discard potential deleterious pleiotropic effects, cystatins may be used as transgenes with impact on agricultural crops by conferring enhanced levels of resistance to phytophagous pests.

Overexpression of sweet sorghum cryptochrome 1a confers hypersensitivity to blue light, abscisic acid and salinity in Arabidopsis

This work provides the bioinformatics, expression pattern and functional analyses of cryptochrome 1a from sweet sorghum (SbCRY1a), together with an exploration of the signaling mechanism mediated by SbCRY1a. Sweet sorghum [Sorghum bicolor (L.) Moench] is considered to be an ideal candidate for biofuel production due to its high efficiency of photosynthesis and the ability to maintain yield under harsh environmental conditions. Blue light receptor cryptochromes regulate multiple aspects of plant growth and development. Here, we reported the function and signal mechanism of sweet sorghum cryptochrome 1a (SbCRY1a) to explore its potential for genetic improvement of sweet sorghum varieties. SbCRY1a transcripts experienced almost 24 h diurnal cycling; however, its protein abundance showed no oscillation. Overexpression of SbCRY1a in Arabidopsis rescued the phenotype of cry1 mutant in a blue light-specific manner and regulated HY5 accumulation under blue light. SbCRY1a protein was present in both nucleus and cytoplasm. The photoexcited SbCRY1a interacted directly with a putative RING E3 ubiquitin ligase constitutive photomorphogenesis 1 (COP1) from sweet sorghum (SbCOP1) instead of SbSPA1 to suppress SbCOP1-SbHY5 interaction responding to blue light. These observations indicate that the function and signaling mechanism of cryptochromes are basically conservative between monocotyledons and dicotyledons. Moreover, SbCRY1a-overexpressed transgenic Arabidopsis showed oversensitive to abscisic acid (ABA) and salinity. The ABA-responsive gene ABI5 was up-regulated evidently in SbCRY1a transgenic lines, suggesting that SbCRY1a might regulate ABA signaling through the HY5-ABI5 regulon.

Selecting Hypomethylated Genomic Regions Using MRE-Seq

Here, we describe a method capable of filtering the hypomethylated part of plant genomes, the so-called hypomethylome. The principle of the method is based on the filtration and sequence analysis of small DNA fragments generated by methylation-sensitive four-cutter restriction endonucleases, possessing ((5me))CpG motifs in their recognition sites. The majority of these fragments represent genes and their flanking regions containing also regulatory elements-the gene space of the genome. Besides the enrichment of the gene space, another advantage of the method is the simultaneous depletion of repetitive elements due to their methylated nature and its easy application on complex and large plant genomes. Additionally to the wet lab procedure, we describe how to analyze the data using bioinformatics methods and how to apply the method to comparative studies.

A comparative synteny analysis tool for target-gene SNP marker discovery: connecting genomics data to breeding in Solanaceae

It is necessary for molecular breeders to overcome the difficulties in applying abundant genomic information to crop breeding. Candidate orthologs would be discovered more efficiently in less-studied crops if the information gained from studies of related crops were used. We developed a comparative analysis tool and web-based genome viewer to identify orthologous genes based synteny as well as sequence similarity between tomato, pepper and potato. The tool has a step-by-step interface with multiple viewing levels to support the easy and accurate exploration of functional orthologs. Furthermore, it provides access to single nucleotide-polymorphism markers from the massive genetic resource pool in order to accelerate the development of molecular markers for candidate orthologs in the Solanaceae. This tool provides a bridge between genome data and breeding by supporting effective marker development, data utilization and communication.

Database URL: http://tgsol.seeders.co.kr/scomp/

Plastome Sequencing of Ten Nonmodel Crop Species Uncovers a Large Insertion of Mitochondrial DNA in Cashew

In plant evolution, intracellular gene transfer (IGT) is a prevalent, ongoing process. While nuclear and mitochondrial genomes are known to integrate foreign DNA via IGT and horizontal gene transfer (HGT), plastid genomes (plastomes) have resisted foreign DNA incorporation and only recently has IGT been uncovered in the plastomes of a few land plants. In this study, we completed plastome sequences for l0 crop species and describe a number of structural features including variation in gene and intron content, inversions, and expansion and contraction of the inverted repeat (IR). We identified a putative in cinnamon ( J. Presl) and other sequenced Lauraceae and an apparent functional transfer of to the nucleus of quinoa ( Willd.). In the orchard tree cashew ( L.), we report the insertion of an ∼6.7-kb fragment of mitochondrial DNA into the plastome IR. BLASTn analyses returned high identity hits to mitogenome sequences including an intact open reading frame. Using three plastome markers for five species of , we generated a phylogeny to investigate the distribution and timing of the insertion. Four species share the insertion, suggesting that this event occurred <20 million yr ago in a single clade in the genus. Our study extends the observation of mitochondrial to plastome IGT to include long-lived tree species. While previous studies have suggested possible mechanisms facilitating IGT to the plastome, more examples of this phenomenon, along with more complete mitogenome sequences, will be required before a common, or variable, mechanism can be elucidated.

Genome sequencing of the staple food crop white Guinea yam enables the development of a molecular marker for sex determination

Background

Root and tuber crops are a major food source in tropical Africa. Among these crops are several species in the monocotyledonous genus Dioscorea collectively known as yam, a staple tuber crop that contributes enormously to the subsistence and socio-cultural lives of millions of people, principally in West and Central Africa. Yam cultivation is constrained by several factors, and yam can be considered a neglected “orphan” crop that would benefit from crop improvement efforts. However, the lack of genetic and genomic tools has impeded the improvement of this staple crop.

Results

To accelerate marker-assisted breeding of yam, we performed genome analysis of white Guinea yam (Dioscorea rotundata) and assembled a 594-Mb genome, 76.4% of which was distributed among 21 linkage groups. In total, we predicted 26,198 genes. Phylogenetic analyses with 2381 conserved genes revealed that Dioscorea is a unique lineage of monocotyledons distinct from the Poales (rice), Arecales (palm), and Zingiberales (banana). The entire Dioscorea genus is characterized by the occurrence of separate male and female plants (dioecy), a feature that has limited efficient yam breeding. To infer the genetics of sex determination, we performed whole-genome resequencing of bulked segregants (quantitative trait locus sequencing [QTL-seq]) in F1 progeny segregating for male and female plants and identified a genomic region associated with female heterogametic (male = ZZ, female = ZW) sex determination. We further delineated the W locus and used it to develop a molecular marker for sex identification of Guinea yam plants at the seedling stage.

Conclusions

Guinea yam belongs to a unique and highly differentiated clade of monocotyledons. The genome analyses and sex-linked marker development performed in this study should greatly accelerate marker-assisted breeding of Guinea yam. In addition, our QTL-seq approach can be utilized in genetic studies of other outcrossing crops and organisms with highly heterozygous genomes. Genomic analysis of orphan crops such as yam promotes efforts to improve food security and the sustainability of tropical agriculture.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-017-0419-x) contains supplementary material, which is available to authorized users.

RNA Sequencing and Co-expressed Long Non-coding RNA in Modern and Wild Wheats

There is an urgent need for the improvement of drought-tolerant bread and durum wheat. The huge and complex genome of bread wheat (BBAADD genome) stands as a vital obstruction for understanding the molecular mechanism underlying drought tolerance. However, tetraploid wheat (Triticum turgidum ssp., BBAA genome) is an ancestor of modern bread wheat and offers an important model for studying the drought response due to its less complex genome. Additionally, several wild relatives of tetraploid wheat have already shown a significant drought tolerance. We sequenced root transcriptome of three tetraploid wheat varieties with varying stress tolerance profiles, and built differential expression library of their transcripts under control and drought conditions. More than 5,000 differentially expressed transcripts were identified from each genotype. Functional characterization of transcripts specific to drought-tolerant genotype, revealed their association with osmolytes production and secondary metabolite pathways. Comparative analysis of differentially expressed genes and their non-coding RNA partners, long noncoding RNAs and microRNAs, provided valuable insight to gene expression regulation in response to drought stress. LncRNAs as well as coding transcripts share similar structural features in different tetraploid species; yet, lncRNAs slightly differ from coding transcripts. Several miRNA-lncRNA target pairs were detected as differentially expressed in drought stress. Overall, this study suggested an important pool of transcripts where their manipulations confer a better performance of wheat varieties under drought stress.

Aberrant splicing in maize rough endosperm3 reveals a conserved role for U12 splicing in eukaryotic multicellular development

RNA splicing of U12-type introns functions in human cell differentiation, but it is not known whether this class of introns has a similar role in plants. The maize ROUGH ENDOSPERM3 (RGH3) protein is orthologous to the human splicing factor, ZRSR2. ZRSR2 mutations are associated with myelodysplastic syndrome (MDS) and cause U12 splicing defects. Maize rgh3 mutants have aberrant endosperm cell differentiation and proliferation. We found that most U12-type introns are retained or misspliced in rgh3 Genes affected in rgh3 and ZRSR2 mutants identify cell cycle and protein glycosylation as common pathways disrupted. Transcripts with retained U12-type introns can be found in polysomes, suggesting that splicing efficiency can alter protein isoforms. The rgh3 mutant protein disrupts colocalization with a known ZRSR2-interacting protein, U2AF2. These results indicate conserved function for RGH3/ZRSR2 in U12 splicing and a deeply conserved role for the minor spliceosome to promote cell differentiation from stem cells to terminal fates.

Protein Degradation Rate in Arabidopsis thaliana Leaf Growth and Development

We applied ¹⁵N labeling approaches to leaves of the Arabidopsis thaliana rosette to characterize their protein degradation rate and understand its determinants. The progressive labeling of new peptides with ¹⁵N and measuring the decrease in the abundance of >60,000 existing peptides over time allowed us to define the degradation rate of 1228 proteins in vivo. We show that Arabidopsis protein half-lives vary from several hours to several months based on the exponential constant of the decay rate for each protein. This rate was calculated from the relative isotope abundance of each peptide and the fold change in protein abundance during growth. Protein complex membership and specific protein domains were found to be strong predictors of degradation rate, while N-end amino acid, hydrophobicity, or aggregation propensity of proteins were not. We discovered rapidly degrading subunits in a variety of protein complexes in plastids and identified the set of plant proteins whose degradation rate changed in different leaves of the rosette and correlated with leaf growth rate. From this information, we have calculated the protein turnover energy costs in different leaves and their key determinants within the proteome.

A High Resolution Radiation Hybrid Map of Wheat Chromosome 4A

Bread wheat has a large and complex allohexaploid genome with low recombination level at chromosome centromeric and peri-centromeric regions. This significantly hampers ordering of markers, contigs of physical maps and sequence scaffolds and impedes obtaining of high-quality reference genome sequence. Here we report on the construction of high-density and high-resolution radiation hybrid (RH) map of chromosome 4A supported by high-density chromosome deletion map. A total of 119 endosperm-based RH lines of two RH panels and 15 chromosome deletion bin lines were genotyped with 90K iSelect single nucleotide polymorphism (SNP) array. A total of 2316 and 2695 markers were successfully mapped to the 4A RH and deletion maps, respectively. The chromosome deletion map was ordered in 19 bins and allowed precise identification of centromeric region and verification of the RH panel reliability. The 4A-specific RH map comprises 1080 mapping bins and spans 6550.9 cR with a resolution of 0.13 Mb/cR. Significantly higher mapping resolution in the centromeric region was observed as compared to recombination maps. Relatively even distribution of deletion frequency along the chromosome in the RH panel was observed and putative functional centromere was delimited within a region characterized by two SNP markers.

Large Differences in Gene Expression Responses to Drought and Heat Stress between Elite Barley Cultivar Scarlett and a Spanish Landrace

Drought causes important losses in crop production every season. Improvement for drought tolerance could take advantage of the diversity held in germplasm collections, much of which has not been incorporated yet into modern breeding. Spanish landraces constitute a promising resource for barley breeding, as they were widely grown until last century and still show good yielding ability under stress. Here, we study the transcriptome expression landscape in two genotypes, an outstanding Spanish landrace-derived inbred line (SBCC073) and a modern cultivar (Scarlett). Gene expression of adult plants after prolonged stresses, either drought or drought combined with heat, was monitored. Transcriptome of mature leaves presented little changes under severe drought, whereas abundant gene expression changes were observed under combined mild drought and heat. Developing inflorescences of SBCC073 exhibited mostly unaltered gene expression, whereas numerous changes were found in the same tissues for Scarlett. Genotypic differences in physiological traits and gene expression patterns confirmed the different behavior of landrace SBCC073 and cultivar Scarlett under abiotic stress, suggesting that they responded to stress following different strategies. A comparison with related studies in barley, addressing gene expression responses to drought, revealed common biological processes, but moderate agreement regarding individual differentially expressed transcripts. Special emphasis was put in the search of co-expressed genes and underlying common regulatory motifs. Overall, 11 transcription factors were identified, and one of them matched cis-regulatory motifs discovered upstream of co-expressed genes involved in those responses.

HopBase: a unified resource for Humulus genomics

Hop (Humulus lupulus L. var lupulus) is a dioecious plant of worldwide significance, used primarily for bittering and flavoring in brewing beer. Studies on the medicinal properties of several unique compounds produced by hop have led to additional interest from pharmacy and healthcare industries as well as livestock production as a natural antibiotic. Genomic research in hop has resulted a published draft genome and transcriptome assemblies. As research into the genomics of hop has gained interest, there is a critical need for centralized online genomic resources. To support the growing research community, we report the development of an online resource "HopBase.org." In addition to providing a gene annotation to the existing Shinsuwase draft genome, HopBase makes available genome assemblies and annotations for both the cultivar "Teamaker" and male hop accession number USDA 21422M. These genome assemblies, gene annotations, along with other common data, coupled with a genome browser and BLAST database enable the hop community to enter the genomic age. The HopBase genomic resource is accessible at http://hopbase.org and http://hopbase.cgrb.oregonstate.edu.

HvPap-1 C1A Protease Participates Differentially in the Barley Response to a Pathogen and an Herbivore

Co-evolutionary processes in plant-pathogen/herbivore systems indicate that protease inhibitors have a particular value in biotic interactions. However, little is known about the defensive role of their targets, the plant proteases. C1A cysteine proteases are the most abundant enzymes responsible for the proteolytic activity during different processes like germination, development and senescence in plants. To identify and characterize C1A cysteine proteases of barley with a potential role in defense, mRNA and protein expression patterns were analyzed in response to biotics stresses. A barley cysteine protease, HvPap-1, previously related to abiotic stresses and grain germination, was particularly induced by flagellin or chitosan elicitation, and biotic stresses such as the phytopathogenic fungus Magnaporthe oryzae or the phytophagous mite Tetranychus urticae. To elucidate the in vivo participation of this enzyme in defense, transformed barley plants overexpressing or silencing HvPap-1 encoding gene were subjected to M. oryzae infection or T. urticae infestation. Whereas overexpressing plants were less susceptible to the fungus than silencing plants, the opposite behavior occurred to the mite. This unexpected result highlights the complexity of the regulatory events leading to the response to a particular biotic stress.

Ensembl Plants: Integrating Tools for Visualizing, Mining, and Analyzing Plant Genomic Data

Ensembl Plants ( http://plants.ensembl.org ) is an integrative resource presenting genome-scale information for 39 sequenced plant species. Available data includes genome sequence, gene models, functional annotation, and polymorphic loci; for the latter, additional information including population structure, individual genotypes, linkage, and phenotype data is available for some species. Comparative data is also available, including genomic alignments and "gene trees," which show the inferred evolutionary history of each gene family represented in the resource. Access to the data is provided through a genome browser, which incorporates many specialist interfaces for different data types, through a variety of programmatic interfaces, and via a specialist data mining tool supporting rapid filtering and retrieval of bulk data. Genomic data from many non-plant species, including those of plant pathogens, pests, and pollinators, is also available via the same interfaces through other divisions of Ensembl.Ensembl Plants is updated 4-6 times a year and is developed in collaboration with our international partners in the Gramene ( http://www.gramene.org ) and transPLANT projects ( http://www.transplantdb.eu ).

Genome-metabolite associations revealed low heritability, high genetic complexity, and causal relations for leaf metabolites in winter wheat (Triticum aestivum)

We investigated associations between the metabolic phenotype, consisting of quantitative data of 76 metabolites from 135 contrasting winter wheat (Triticum aestivum) lines, and 17 372 single nucleotide polymorphism (SNP) markers. Metabolite profiles were generated from flag leaves of plants from three different environments, with average repeatabilities of 0.5-0.6. The average heritability of 0.25 was unaffected by the heading date. Correlations among metabolites reflected their functional grouping, highlighting the strict coordination of various routes of the citric acid cycle. Genome-wide association studies identified significant associations for six metabolic traits, namely oxalic acid, ornithine, L-arginine, pentose alcohol III, L-tyrosine, and a sugar oligomer (oligo II), with between one and 17 associated SNPs. Notable associations with genes regulating transcription or translation explained between 2.8% and 32.5% of the genotypic variance (pG). Further candidate genes comprised metabolite carriers (pG 32.5-38.1%), regulatory proteins (pG 0.3-11.1%), and metabolic enzymes (pG 2.5-32.5%). The combinatorial use of genomic and metabolic data to construct partially directed networks revealed causal inferences in the correlated metabolite traits and associated SNPs. The evaluated causal relationships will provide a basis for predicting the effects of genetic interferences on groups of correlated metabolic traits, and thus on specific metabolic phenotypes.

Ca2+/Cation Antiporters (CaCA): Identification, Characterization and Expression Profiling in Bread Wheat (Triticum aestivum L.)

The Ca2+/cation antiporters (CaCA) superfamily proteins play vital function in Ca2+ ion homeostasis, which is an important event during development and defense response. Molecular characterization of these proteins has been performed in certain plants, but they are still not characterized in Triticum aestivum (bread wheat). Herein, we identified 34 TaCaCA superfamily proteins, which were classified into TaCAX, TaCCX, TaNCL, and TaMHX protein families based on their structural organization and evolutionary relation with earlier reported proteins. Since the T. aestivum comprises an allohexaploid genome, TaCaCA genes were derived from each A, B, and D subgenome and homeologous chromosome (HC), except chromosome-group 1. Majority of genes were derived from more than one HCs in each family that were considered as homeologous genes (HGs) due to their high similarity with each other. These HGs showed comparable gene and protein structures in terms of exon/intron organization and domain architecture. Majority of TaCaCA proteins comprised two Na_Ca_ex domains. However, TaNCLs consisted of an additional EF-hand domain with calcium binding motifs. Each TaCaCA protein family consisted of about 10 transmembrane and two α-repeat regions with specifically conserved signature motifs except TaNCL, which had single α-repeat. Variable expression of most of the TaCaCA genes during various developmental stages suggested their specified role in development. However, constitutively high expression of a few genes like TaCAX1-A and TaNCL1-B indicated their role throughout the plant growth and development. The modulated expression of certain genes during biotic (fungal infections) and abiotic stresses (heat, drought, salt) suggested their role in stress response. Majority of TaCCX and TaNCL family genes were found highly affected during various abiotic stresses. However, the role of individual gene needs to be established. The present study unfolded the opportunity for detail functional characterization of TaCaCA proteins and their utilization in future crop improvement programs.

Sequencing and comparative analyses of Aegilops tauschii chromosome arm 3DS reveal rapid evolution of Triticeae genomes

Bread wheat (Triticum aestivum, AABBDD) is an allohexaploid species derived from two rounds of interspecific hybridizations. A high-quality genome sequence assembly of diploid Aegilops tauschii, the donor of the wheat D genome, will provide a useful platform to study polyploid wheat evolution. A combined approach of BAC pooling and next-generation sequencing technology was employed to sequence the minimum tiling path (MTP) of 3176 BAC clones from the short arm of Ae. tauschii chromosome 3 (At3DS). The final assembly of 135 super-scaffolds with an N50 of 4.2 Mb was used to build a 247-Mb pseudomolecule with a total of 2222 predicted protein-coding genes. Compared with the orthologous regions of rice, Brachypodium, and sorghum, At3DS contains 38.67% more genes. In comparison to At3DS, the short arm sequence of wheat chromosome 3B (Ta3BS) is 95-Mb large in size, which is primarily due to the expansion of the non-centromeric region, suggesting that transposable element (TE) bursts in Ta3B likely occurred there. Also, the size increase is accompanied by a proportional increase in gene number in Ta3BS. We found that in the sequence of short arm of wheat chromosome 3D (Ta3DS), there was only less than 0.27% gene loss compared to At3DS. Our study reveals divergent evolution of grass genomes and provides new insights into sequence changes in the polyploid wheat genome.

Bioinformatics Approach in Plant Genomic Research

The advance in genomics technology leads to the dramatic change in plant biology research. Plant biologists now easily access to enormous genomic data to deeply study plant high-density genetic variation at molecular level. Therefore, fully understanding and well manipulating bioinformatics tools to manage and analyze these data are essential in current plant genome research. Many plant genome databases have been established and continued expanding recently. Meanwhile, analytical methods based on bioinformatics are also well developed in many aspects of plant genomic research including comparative genomic analysis, phylogenomics and evolutionary analysis, and genome-wide association study. However, constantly upgrading in computational infrastructures, such as high capacity data storage and high performing analysis software, is the real challenge for plant genome research. This review paper focuses on challenges and opportunities which knowledge and skills in bioinformatics can bring to plant scientists in present plant genomics era as well as future aspects in critical need for effective tools to facilitate the translation of knowledge from new sequencing data to enhancement of plant productivity.

Mapping-by-sequencing in complex polyploid genomes using genic sequence capture: a case study to map yellow rust resistance in hexaploid wheat

Previously we extended the utility of mapping-by-sequencing by combining it with sequence capture and mapping sequence data to pseudo-chromosomes that were organized using wheat-Brachypodium synteny. This, with a bespoke haplotyping algorithm, enabled us to map the flowering time locus in the diploid wheat Triticum monococcum L. identifying a set of deleted genes (Gardiner et al., 2014). Here, we develop this combination of gene enrichment and sliding window mapping-by-synteny analysis to map the Yr6 locus for yellow stripe rust resistance in hexaploid wheat. A 110 MB NimbleGen capture probe set was used to enrich and sequence a doubled haploid mapping population of hexaploid wheat derived from an Avalon and Cadenza cross. The Yr6 locus was identified by mapping to the POPSEQ chromosomal pseudomolecules using a bespoke pipeline and algorithm (Chapman et al., 2015). Furthermore the same locus was identified using newly developed pseudo-chromosome sequences as a mapping reference that are based on the genic sequence used for sequence enrichment. The pseudo-chromosomes allow us to demonstrate the application of mapping-by-sequencing to even poorly defined polyploidy genomes where chromosomes are incomplete and sub-genome assemblies are collapsed. This analysis uniquely enabled us to: compare wheat genome annotations; identify the Yr6 locus - defining a smaller genic region than was previously possible; associate the interval with one wheat sub-genome and increase the density of SNP markers associated. Finally, we built the pipeline in iPlant, making it a user-friendly community resource for phenotype mapping.

Long non-coding RNAs are major contributors to transcriptome changes in sunflower meiocytes with different recombination rates

BACKGROUND

Meiosis is a form of specialized cell division that marks the transition from diploid meiocyte to haploid gamete, and provides an opportunity for genetic reassortment through recombination. Experimental data indicates that, relative to their wild ancestors, cultivated sunflower varieties show a higher recombination rate during meiosis. To better understand the molecular basis for this difference, we compared gene expression in male sunflower meiocytes in prophase I isolated from a domesticated line, a wild relative, and a F1 hybrid of the two.

RESULTS

Of the genes that showed differential expression between the wild and domesticated genotypes, 63.62 % could not be identified as protein-coding genes, and of these genes, 70.98 % passed stringent filters to be classified as long non-coding RNAs (lncRNAs). Compared to the sunflower somatic transcriptome, meiocytes express a higher proportion of lncRNAs, and the majority of genes with exclusive expression in meiocytes were lncRNAs. Around 40 % of the lncRNAs showed sequence similarity with small RNAs (sRNA), while 1.53 % were predicted to be sunflower natural antisense transcripts (NATs), and 9.18 % contained transposable elements (TE). We identified 6895 lncRNAs that are exclusively expressed in meiocytes, these lncRNAs appear to have higher conservation, a greater degree of differential expression, a higher proportion of sRNA similarity, and higher TE content relative to lncRNAs that are also expressed in the somatic transcriptome.

CONCLUSIONS

lncRNAs play important roles in plant meiosis and may participate in chromatin modification processes, although other regulatory functions cannot be excluded. lncRNAs could also be related to the different recombination rates seen for domesticated and wild sunflowers.

Elucidating the triplicated ancestral genome structure of radish based on chromosome-level comparison with the Brassica genomes

This study presents a chromosome-scale draft genome sequence of radish that is assembled into nine chromosomal pseudomolecules. A comprehensive comparative genome analysis with the Brassica genomes provides genomic evidences on the evolution of the mesohexaploid radish genome. Radish (Raphanus sativus L.) is an agronomically important root vegetable crop and its origin and phylogenetic position in the tribe Brassiceae is controversial. Here we present a comprehensive analysis of the radish genome based on the chromosome sequences of R. sativus cv. WK10039. The radish genome was sequenced and assembled into 426.2 Mb spanning >98 % of the gene space, of which 344.0 Mb were integrated into nine chromosome pseudomolecules. Approximately 36 % of the genome was repetitive sequences and 46,514 protein-coding genes were predicted and annotated. Comparative mapping of the tPCK-like ancestral genome revealed that the radish genome has intermediate characteristics between the Brassica A/C and B genomes in the triplicated segments, suggesting an internal origin from the genus Brassica. The evolutionary characteristics shared between radish and other Brassica species provided genomic evidences that the current form of nine chromosomes in radish was rearranged from the chromosomes of hexaploid progenitor. Overall, this study provides a chromosome-scale draft genome sequence of radish as well as novel insight into evolution of the mesohexaploid genomes in the tribe Brassiceae.

SorghumFDB: sorghum functional genomics database with multidimensional network analysis

Sorghum (Sorghum bicolor [L.] Moench) has excellent agronomic traits and biological properties, such as heat and drought-tolerance. It is a C₄ grass and potential bioenergy-producing plant, which makes it an important crop worldwide. With the sorghum genome sequence released, it is essential to establish a sorghum functional genomics data mining platform. We collected genomic data and some functional annotations to construct a sorghum functional genomics database (SorghumFDB). SorghumFDB integrated knowledge of sorghum gene family classifications (transcription regulators/factors, carbohydrate-active enzymes, protein kinases, ubiquitins, cytochrome P450, monolignol biosynthesis related enzymes, R-genes and organelle-genes), detailed gene annotations, miRNA and target gene information, orthologous pairs in the model plants Arabidopsis, rice and maize, gene loci conversions and a genome browser. We further constructed a dynamic network of multidimensional biological relationships, comprised of the co-expression data, protein–protein interactions and miRNA-target pairs. We took effective measures to combine the network, gene set enrichment and motif analyses to determine the key regulators that participate in related metabolic pathways, such as the lignin pathway, which is a major biological process in bioenergy-producing plants.

Database URL:http://structuralbiology.cau.edu.cn/sorghum/index.html.

cDNA Library Enrichment of Full Length Transcripts for SMRT Long Read Sequencing

The utility of genome assemblies does not only rely on the quality of the assembled genome sequence, but also on the quality of the gene annotations. The Pacific Biosciences Iso-Seq technology is a powerful support for accurate eukaryotic gene model annotation as it allows for direct readout of full-length cDNA sequences without the need for noisy short read-based transcript assembly. We propose the implementation of the TeloPrime Full Length cDNA Amplification kit to the Pacific Biosciences Iso-Seq technology in order to enrich for genuine full-length transcripts in the cDNA libraries. We provide evidence that TeloPrime outperforms the commonly used SMARTer PCR cDNA Synthesis Kit in identifying transcription start and end sites in Arabidopsis thaliana. Furthermore, we show that TeloPrime-based Pacific Biosciences Iso-Seq can be successfully applied to the polyploid genome of bread wheat (Triticum aestivum) not only to efficiently annotate gene models, but also to identify novel transcription sites, gene homeologs, splicing isoforms and previously unidentified gene loci.

Methods for accurate quantification of LTR-retrotransposon copy number using short-read sequence data: a case study in Sorghum

Transposable elements (TEs) are ubiquitous in eukaryotic genomes and their mobility impacts genome structure and function in myriad ways. Because of their abundance, activity, and repetitive nature, the characterization and analysis of TEs remain challenging, particularly from short-read sequencing projects. To overcome this difficulty, we have developed a method that estimates TE copy number from short-read sequences. To test the accuracy of our method, we first performed an in silico analysis of the reference Sorghum bicolor genome, using both reference-based and de novo approaches. The resulting TE copy number estimates were strikingly similar to the annotated numbers. We then tested our method on real short-read data by estimating TE copy numbers in several accessions of S. bicolor and its close relative S. propinquum. Both methods effectively identify and rank similar TE families from highest to lowest abundance. We found that de novo characterization was effective at capturing qualitative variation, but underestimated the abundance of some TE families, specifically families of more ancient origin. Also, interspecific reference-based mapping of S. propinquum reads to the S. bicolor database failed to fully describe TE content in S. propinquum, indicative of recent TE activity leading to changes in the respective repetitive landscapes over very short evolutionary timescales. We conclude that reference-based analyses are best suited for within-species comparisons, while de novo approaches are more reliable for evolutionarily distant comparisons.

Comparative structural analysis of Bru1 region homeologs in Saccharum spontaneum and S. officinarum

BACKGROUND

Sugarcane is a major sugar and biofuel crop, but genomic research and molecular breeding have lagged behind other major crops due to the complexity of auto-allopolyploid genomes. Sugarcane cultivars are frequently aneuploid with chromosome number ranging from 100 to 130, consisting of 70-80 % S. officinarum, 10-20 % S. spontaneum, and 10 % recombinants between these two species. Analysis of a genomic region in the progenitor autoploid genomes of sugarcane hybrid cultivars will reveal the nature and divergence of homologous chromosomes.

RESULTS

To investigate the origin and evolution of haplotypes in the Bru1 genomic regions in sugarcane cultivars, we identified two BAC clones from S. spontaneum and four from S. officinarum and compared to seven haplotype sequences from sugarcane hybrid R570. The results clarified the origin of seven homologous haplotypes in R570, four haplotypes originated from S. officinarum, two from S. spontaneum and one recombinant.. Retrotransposon insertions and sequences variations among the homologous haplotypes sequence divergence ranged from 18.2 % to 60.5 % with an average of 33.7 %. Gene content and gene structure were relatively well conserved among the homologous haplotypes. Exon splitting occurred in haplotypes of the hybrid genome but not in its progenitor genomes. Tajima's D analysis revealed that S. spontaneum hapotypes in the Bru1 genomic regions were under strong directional selection. Numerous inversions, deletions, insertions and translocations were found between haplotypes within each genome.

CONCLUSIONS

This is the first comparison among haplotypes of a modern sugarcane hybrid and its two progenitors. Tajima's D results emphasized the crucial role of this fungal disease resistance gene for enhancing the fitness of this species and indicating that the brown rust resistance gene in R570 is from S. spontaneum. Species-specific InDel, sequences similarity and phylogenetic analysis of homologous genes can be used for identifying the origin of S. spontaneum and S. officinarum haplotype in Saccharum hybrids. Comparison of exon splitting among the homologous haplotypes suggested that the genome rearrangements in Saccharum hybrids after hybridization. The combined minimum difference at 19.5 % among homologous chromosomes in S. officinarum would be sufficient for proper genome assembly of this autopolyploid genome. Retrotransposon insertions and sequences variations among the homologous haplotypes sequence divergence may allow sequencing and assembling the autopolyploid Saccharum genomes and the auto-allopolyploid hybrid genomes using whole genome shotgun sequencing.

PGSB/MIPS Plant Genome Information Resources and Concepts for the Analysis of Complex Grass Genomes

PGSB (Plant Genome and Systems Biology; formerly MIPS-Munich Institute for Protein Sequences) has been involved in developing, implementing and maintaining plant genome databases for more than a decade. Genome databases and analysis resources have focused on individual genomes and aim to provide flexible and maintainable datasets for model plant genomes as a backbone against which experimental data, e.g., from high-throughput functional genomics, can be organized and analyzed. In addition, genomes from both model and crop plants form a scaffold for comparative genomics, assisted by specialized tools such as the CrowsNest viewer to explore conserved gene order (synteny) between related species on macro- and micro-levels.The genomes of many economically important Triticeae plants such as wheat, barley, and rye present a great challenge for sequence assembly and bioinformatic analysis due to their enormous complexity and large genome size. Novel concepts and strategies have been developed to deal with these difficulties and have been applied to the genomes of wheat, barley, rye, and other cereals. This includes the GenomeZipper concept, reference-guided exome assembly, and "chromosome genomics" based on flow cytometry sorted chromosomes.

Structural organization of fatty acid desaturase loci in linseed lines with contrasting linolenic acid contents

Flax (Linum usitatissimum L.), the richest crop source of omega-3 fatty acids (FAs), is a diploid plant with an estimated genome size of ~370 Mb and is well suited for studying genomic organization of agronomically important traits. In this study, 12 bacterial artificial chromosome clones harbouring the six FA desaturase loci sad1, sad2, fad2a, fad2b, fad3a and fad3b from the conventional variety CDC Bethune and the high linolenic acid line M5791 were sequenced, analysed and compared to determine the structural organization of these loci and to gain insights into the genetic mechanisms underlying FA composition in flax. With one gene every 3.2-4.6 kb, the desaturase loci have a higher gene density than the genome's average of one gene per 7.8-8.2 kb. The gene order and orientation across the two genotypes were generally conserved with the exception of the sad1 locus that was predicted to have additional genes in CDC Bethune. High sequence conservation in both genic and intergenic regions of the sad and fad2b loci contrasted with the significant level of variation of the fad2a and fad3 loci, with SNPs being the most frequently observed mutation type. The fad2a locus had 297 SNPs and 36 indels over ~95 kb contrasting with the fad2b locus that had a mere seven SNPs and four indels in ~110 kb. Annotation of the gene-rich loci revealed other genes of known role in lipid or carbohydrate metabolic/catabolic pathways. The organization of the fad2b locus was particularly complex with seven copies of the fad2b gene in both genotypes. The presence of Gypsy, Copia, MITE, Mutator, hAT and other novel repeat elements at the desaturase loci was similar to that of the whole genome. This structural genomic analysis provided some insights into the genomic organization and composition of the main desaturase loci of linseed and of their complex evolution through both tandem and whole genome duplications.

Reactive Oxygen Species Play a Role in the Infection of the Necrotrophic Fungi, Rhizoctonia solani in Wheat

Rhizoctonia solani is a nectrotrophic fungal pathogen that causes billions of dollars of damage to agriculture worldwide and infects a broad host range including wheat, rice, potato and legumes. In this study we identify wheat genes that are differentially expressed in response to the R. solani isolate, AG8, using microarray technology. A significant number of wheat genes identified in this screen were involved in reactive oxygen species (ROS) production and redox regulation. Levels of ROS species were increased in wheat root tissue following R. solani infection as determined by Nitro Blue Tetrazolium (NBT), 3,3'-diaminobenzidine (DAB) and titanium sulphate measurements. Pathogen/ROS related genes from R. solani were also tested for expression patterns upon wheat infection. TmpL, a R. solani gene homologous to a gene associated with ROS regulation in Alternaria brassicicola, and OAH, a R. solani gene homologous to oxaloacetate acetylhydrolase which has been shown to produce oxalic acid in Sclerotinia sclerotiorum, were highly induced in R. solani when infecting wheat. We speculate that the interplay between the wheat and R. solani ROS generating proteins may be important for determining the outcome of the wheat/R. solani interaction.

Wheat miRNA ancestors: evident by transcriptome analysis of A, B, and D genome donors

MicroRNAs are critical players of post-transcriptional gene regulation with profound effects on the fundamental processes of cellular life. Their identification and characterization, together with their targets, hold great significance in exploring and exploiting their roles on a functional context, providing valuable clues into the regulation of important biological processes, such as stress tolerance or environmental adaptation. Wheat is a hardy crop, extensively harvested in temperate regions, and is a major component of the human diet. With the advent of the next generation sequencing technologies considerably decreasing sequencing costs per base-pair, genomic, and transcriptomic data from several wheat species, including the progenitors and wild relatives have become available. In this study, we performed in silico identification and comparative analysis of microRNA repertoires of bread wheat (Triticum aestivum L.) and its diploid progenitors and relatives, Aegilops sharonensis, Aegilops speltoides, Aegilops tauschii, Triticum monococcum, and Triticum urartu through the utilization of publicly available transcriptomic data. Over 200 miRNA families were identified, majority of which have not previously been reported. Ancestral relationships expanded our understanding of wheat miRNA evolution, while T. monococcum miRNAs delivered important clues on the effects of domestication on miRNA expression. Comparative analyses on wild Ae. sharonensis accessions highlighted candidate miRNAs that can be linked to stress tolerance. The miRNA repertoires of bread wheat and its diploid progenitors and relatives provide important insight into the diversification and distribution of miRNA genes, which should contribute to the elucidation of miRNA evolution of Poaceae family. A thorough understanding of the convergent and divergent expression profiles of miRNAs in different genetic backgrounds can provide unique opportunities to modulation of gene regulation for better crop performance.

CG Methylation Covaries with Differential Gene Expression between Leaf and Floral Bud Tissues of Brachypodium distachyon

DNA methylation has the potential to influence plant growth and development through its influence on gene expression. To date, however, the evidence from plant systems is mixed as to whether patterns of DNA methylation vary significantly among tissues and, if so, whether these differences affect tissue-specific gene expression. To address these questions, we analyzed both bisulfite sequence (BSseq) and transcriptomic sequence data from three biological replicates of two tissues (leaf and floral bud) from the model grass species Brachypodium distachyon. Our first goal was to determine whether tissues were more differentiated in DNA methylation than explained by variation among biological replicates. Tissues were more differentiated than biological replicates, but the analysis of replicated data revealed high (>50%) false positive rates for the inference of differentially methylated sites (DMSs) and differentially methylated regions (DMRs). Comparing methylation to gene expression, we found that differential CG methylation consistently covaried negatively with gene expression, regardless as to whether methylation was within genes, within their promoters or even within their closest transposable element. The relationship between gene expression and either CHG or CHH methylation was less consistent. In total, CG methylation in promoters explained 9% of the variation in tissue-specific expression across genes, suggesting that CG methylation is a minor but appreciable factor in tissue differentiation.

VGSC: A Web-Based Vector Graph Toolkit of Genome Synteny and Collinearity

BACKGROUND

In order to understand the colocalization of genetic loci amongst species, synteny and collinearity analysis is a frequent task in comparative genomics research. However many analysis software packages are not effective in visualizing results. Problems include lack of graphic visualization, simple representation, or inextensible format of outputs. Moreover, higher throughput sequencing technology requires higher resolution image output.

IMPLEMENTATION

To fill this gap, this paper publishes VGSC, the Vector Graph toolkit of genome Synteny and Collinearity, and its online service, to visualize the synteny and collinearity in the common graphical format, including both raster (JPEG, Bitmap, and PNG) and vector graphic (SVG, EPS, and PDF).

RESULT

Users can upload sequence alignments from blast and collinearity relationship from the synteny analysis tools. The website can generate the vector or raster graphical results automatically. We also provide a java-based bytecode binary to enable the command-line execution.

Finding and Characterizing Repeats in Plant Genomes

Plant genomes contain a particularly high proportion of repeated structures of various types. This chapter proposes a guided tour of available software that can help biologists to look for these repeats and check some hypothetical models intended to characterize their structures. Since transposable elements are a major source of repeats in plants, many methods have been used or developed for this large class of sequences. They are representative of the range of tools available for other classes of repeats and we have provided a whole section on this topic as well as a selection of the main existing software. In order to better understand how they work and how repeats may be efficiently found in genomes, it is necessary to look at the technical issues involved in the large-scale search of these structures. Indeed, it may be hard to keep up with the profusion of proposals in this dynamic field and the rest of the chapter is devoted to the foundations of the search for repeats and more complex patterns. The second section introduces the key concepts that are useful for understanding the current state of the art in playing with words, applied to genomic sequences. This can be seen as the first stage of a very general approach called linguistic analysis that is interested in the analysis of natural or artificial texts. Words, the lexical level, correspond to simple repeated entities in texts or strings. In fact, biologists need to represent more complex entities where a repeat family is built on more abstract structures, including direct or inverted small repeats, motifs, composition constraints as well as ordering and distance constraints between these elementary blocks. In terms of linguistics, this corresponds to the syntactic level of a language. The last section introduces concepts and practical tools that can be used to reach this syntactic level in biological sequence analysis.

A Comprehensive Prescription for Plant miRNA Identification

microRNAs (miRNAs) are tiny ribo-regulatory molecules involved in various essential pathways for persistence of cellular life, such as development, environmental adaptation, and stress response. In recent years, miRNAs have become a major focus in molecular biology because of their functional and diagnostic importance. This interest in miRNA research has resulted in the development of many specific software and pipelines for the identification of miRNAs and their specific targets, which is the key for the elucidation of miRNA-modulated gene expression. While the well-recognized importance of miRNAs in clinical research pushed the emergence of many useful computational identification approaches in animals, available software and pipelines are fewer for plants. Additionally, existing approaches suffers from mis-identification and annotation of plant miRNAs since the miRNA mining process for plants is highly prone to false-positives, particularly in cereals which have a highly repetitive genome. Our group developed a homology-based in silico miRNA identification approach for plants, which utilizes two Perl scripts "SUmirFind" and "SUmirFold" and since then, this method helped identify many miRNAs particularly from crop species such as Triticum or Aegliops. Herein, we describe a comprehensive updated guideline by the implementation of two new scripts, "SUmirPredictor" and "SUmirLocator," and refinements to our previous method in order to identify genuine miRNAs with increased sensitivity in consideration of miRNA identification problems in plants. Recent updates enable our method to provide more reliable and precise results in an automated fashion in addition to solutions for elimination of most false-positive predictions, miRNA naming and miRNA mis-annotation. It also provides a comprehensive view to genome/transcriptome-wide location of miRNA precursors as well as their association with transposable elements. The "SUmirPredictor" and "SUmirLocator" scripts are freely available together with a reference high-confidence plant miRNA list.