Search

Publication of the complete genome sequence of Arabidopsis thaliana, the first plant reference genome, in December 2000 heralded the beginning of the plant genome era. Over the past 20 years reference genomes have been generated for hundreds of plant species, spanning non-vascular to flowering plants. Releasing these plant genomes has dramatically advanced studies in all disciplines of plant biology. Importantly, multiple reference-level genomes have been generated for the major crops and their progenitors, enabling the creation of pan-genomes and exploration of domestication history and natural variations that can be adopted by modern crop breeding. We summarize the progress of plant genome sequencing and the challenges of sequencing more complex plant genomes and generating pan-genomes., Competing Interests: Declaration of Interests The authors have no interests to declare., (Copyright © 2021 Elsevier Ltd. All rights reserved.)

As one of the great survivors of the plant kingdom, barnyard grasses (Echinochloa spp.) are the most noxious and common weeds in paddy ecosystems. Meanwhile, at least two Echinochloa species have been domesticated and cultivated as millets. In order to better understand the genomic forces driving the evolution of Echinochloa species toward weed and crop characteristics, we assemble genomes of three Echinochloa species (allohexaploid E. crus-galli and E. colona, and allotetraploid E. oryzicola) and re-sequence 737 accessions of barnyard grasses and millets from 16 rice-producing countries. Phylogenomic and comparative genomic analyses reveal the complex and reticulate evolution in the speciation of Echinochloa polyploids and provide evidence of constrained disease-related gene copy numbers in Echinochloa. A population-level investigation uncovers deep population differentiation for local adaptation, multiple target-site herbicide resistance mutations of barnyard grasses, and limited domestication of barnyard millets. Our results provide genomic insights into the dual roles of Echinochloa species as weeds and crops as well as essential resources for studying plant polyploidization, adaptation, precision weed control and millet improvements., (© 2022. The Author(s).)

Momilactones are bioactive diterpenoids that contribute to plant defense against pathogens and allelopathic interactions between plants. Both cultivated and wild grass species of Oryza and Echinochloa crus-galli (barnyard grass) produce momilactones using a biosynthetic gene cluster (BGC) in their genomes. The bryophyte Calohypnum plumiforme (formerly Hypnum plumaeforme ) also produces momilactones, and the bifunctional diterpene cyclase gene CpDTC1/HpDTC1, which is responsible for the production of the diterpene framework, has been characterized. To understand the molecular architecture of the momilactone biosynthetic genes in the moss genome and their evolutionary relationships with other momilactone-producing plants, we sequenced and annotated the C. plumiforme genome. The data revealed a 150-kb genomic region that contains two cytochrome P450 genes, the CpDTC1 / HpDTC1 gene and the "dehydrogenase momilactone A synthase" gene tandemly arranged and inductively transcribed following stress exposure. The predicted enzymatic functions in yeast and recombinant assay and the successful pathway reconstitution in Nicotiana benthamiana suggest that it is a functional BGC responsible for momilactone production. Furthermore, in a survey of genomic sequences of a broad range of plant species, we found that momilactone BGC is limited to the two grasses ( Oryza and Echinochloa ) and C. plumiforme , with no synteny among these genomes. These results indicate that while the gene cluster in C. plumiforme is functionally similar to that in rice and barnyard grass, it is likely a product of convergent evolution. To the best of our knowledge, this report of a BGC for a specialized plant defense metabolite in bryophytes is unique., Competing Interests: The authors declare no competing interest.

Brassica napus (AACC, 2n = 38) is an important oilseed crop grown worldwide. However, little is known about the population evolution of this species, the genomic difference between its major genetic groups, such as European and Asian rapeseed, and the impacts of historical large-scale introgression events on this young tetraploid. In this study, we reported the de novo assembly of the genome sequences of an Asian rapeseed (B. napus), Ningyou 7, and its four progenitors and compared these genomes with other available genomic data from diverse European and Asian cultivars. Our results showed that Asian rapeseed originally derived from European rapeseed but subsequently significantly diverged, with rapid genome differentiation after hybridization and intensive local selective breeding. The first historical introgression of B. rapa dramatically broadened the allelic pool but decreased the deleterious variations of Asian rapeseed. The second historical introgression of the double-low traits of European rapeseed (canola) has reshaped Asian rapeseed into two groups (double-low and double-high), accompanied by an increase in genetic load in the double-low group. This study demonstrates distinctive genomic footprints and deleterious SNP (single nucleotide polymorphism) variants for local adaptation by recent intra- and interspecies introgression events and provides novel insights for understanding the rapid genome evolution of a young allopolyploid crop., (© 2019 The Authors. Plant Biotechnology Journal published by Society for Experimental Biology and The Association of Applied Biologists and John Wiley & Sons Ltd.)

Rice (Oryza sativa L.) is one of the most important crops worldwide. Its relatives, including phylogenetically related species of rice and paddy weeds with a similar ecological niche, can provide crucial genetic resources (such as resistance to biotic and abiotic stresses and high photosynthetic efficiency) for rice research. Although many rice genomic databases have been constructed, a database providing large-scale curated genomic data from rice relatives and offering specific gene resources is still lacking. Here, we present RiceRelativesGD, a user-friendly genomic database of rice relatives. RiceRelativesGD integrates large-scale genomic resources from 2 cultivated rice and 11 rice relatives, including 208 321 specific genes and 13 643 genes related to photosynthesis and responsive to external stimuli. Diverse bioinformatics tools are embedded in the database, which allow users to search, visualize and download the information of interest. To our knowledge, this is the first genomic database providing a centralized genetic resource of rice relatives. RiceRelativesGD will serve as a significant and comprehensive knowledgebase for the rice community., (© The Author(s) 2019. Published by Oxford University Press.)

Agronomically critical weeds that have evolved alongside crop species are characterized by rapid adaptation and invasiveness, which can result in an enormous reduction in annual crop yield worldwide. We discuss here recent genome-based research studies on agricultural weeds and crop-weed interactions that reveal several major evolutionary innovations such as de-domestication, interactions mediated by allelochemical secondary metabolites, and parasitic genetic elements that play crucial roles in enhancing weed invasiveness in agricultural settings. We believe that these key studies will guide future research into the evolution of crop-weed interactions, and further the development of practical applications in agricultural weed control and crop breeding., (Copyright © 2018 Elsevier Ltd. All rights reserved.)

Background: Long noncoding RNAs (lncRNAs) are transcripts longer than 200 bp that do not encode proteins but nonetheless have been shown to play important roles in various biological processes in plants. Brassica napus is an important seed oil crop worldwide and the target of many genetic improvement activities. To understand better the function of lncRNAs in regulating plant metabolic activities, we carried out a genome-wide lncRNA identification of lncRNAs in Brassica napus with a focus on lncRNAs involved in lipid metabolism. Twenty ribosomal RNA depleted strand specific RNA-seq (ssRNA-seq) datasets were generatred using RNAs isolated from B. napus seeds at four developmental stages. For comparison we also included 30 publically available RNA-seq datasets generated from poly(A) enriched mRNAs isolated from from various Brassica napus tissues in our analysis., Results: A total of 8905 lncRNA loci were identified, including 7100 long intergenic noncoding RNA (lincRNA) loci and 1805 loci generating long noncoding natural antisense transcript (lncNAT). Many lncRNAs were identified only in the ssRNA-seq and poly(A) RNA-seq dataset, suggesting that B. napus has a large lncRNA repertoire and it is necessary to use libraries prepared from different tissues and developmental stages as well as different library preparation approaches to capture the whole spectrum of lncRNAs. Analysis of coexpression networks revealed that among the regulatory modules are networks containing lncRNAs and protein-coding genes related to oil biosynthesis indicating a possible role of lncRNAs in the control of lipid metabolism. One such example is that several lncRNAs are potential regulators of BnaC08g11970D that encodes oleosin1, a protein found in oil bodies and involved in seed lipid accumulation. We also observed that the expression levels of B. napus lncRNAs is positively correlated with their conservation levels., Conclusions: We demonstrated that the B. napus genome has a large number of lncRNA and that these lncRNAs are expressed broadly across many developmental times and in different tissue types. We also provide evidence indicating that specific lncRNAs appear to be important regulators of lipid biosynthesis forming regulatory networks with transcripts involved in lipid biosynthesis. We also provide evidence that these lncRNAs are conserved in other species of the Brassicaceae family.

Background: Characterization of genetic variations in maize has been challenging, mainly due to deterioration of collinearity between individual genomes in the species. An international consortium of maize research groups combined resources to develop the maize haplotype version 3 (HapMap 3), built from whole-genome sequencing data from 1218 maize lines, covering predomestication and domesticated Zea mays varieties across the world., Results: A new computational pipeline was set up to process more than 12 trillion bp of sequencing data, and a set of population genetics filters was applied to identify more than 83 million variant sites., Conclusions: We identified polymorphisms in regions where collinearity is largely preserved in the maize species. However, the fact that the B73 genome used as the reference only represents a fraction of all haplotypes is still an important limiting factor.

Barnyardgrass (Echinochloa crus-galli) is a pernicious weed in agricultural fields worldwide. The molecular mechanisms underlying its success in the absence of human intervention are presently unknown. Here we report a draft genome sequence of the hexaploid species E. crus-galli, i.e., a 1.27 Gb assembly representing 90.7% of the predicted genome size. An extremely large repertoire of genes encoding cytochrome P450 monooxygenases and glutathione S-transferases associated with detoxification are found. Two gene clusters involved in the biosynthesis of an allelochemical 2,4-dihydroxy-7-methoxy-1,4-benzoxazin-3-one (DIMBOA) and a phytoalexin momilactone A are found in the E. crus-galli genome, respectively. The allelochemical DIMBOA gene cluster is activated in response to co-cultivation with rice, while the phytoalexin momilactone A gene cluster specifically to infection by pathogenic Pyricularia oryzae. Our results provide a new understanding of the molecular mechanisms underlying the extreme adaptation of the weed.

Heterosis is a fundamental biological phenomenon characterized by the superior performance of a hybrid compared with its parents. The underlying molecular basis for heterosis, particularly for allopolyploids, remains elusive. In this study we analyzed the transcriptomes of Brassica napus parental lines and their F 1 hybrids at three stages of early flower development. Phenotypically, the F 1 hybrids show remarkable heterosis in silique number and grain yield. Transcriptome analysis revealed that various phytohormone (auxin and salicylic acid) response genes are significantly altered in the F 1 hybrids relative to the parental lines. We also found evidence for decreased expression divergence of the homoeologous gene pairs in the allopolyploid F 1 hybrids and suggest that high-parental expression-level dominance plays an important role in heterosis. Small RNA and methylation studies aimed at examining the epigenetic effect of the changes in gene expression level in the F 1 hybrids showed that the majority of the small interfering RNA (siRNA) clusters had a higher expression level in the F 1 hybrids than in the parents, and that there was an increase in genome-wide DNA methylation in the F 1 hybrid. Transposable elements associated with siRNA clusters had a higher level of methylation and a lower expression level in the F 1 hybrid, implying that the non-additively expressed siRNA clusters resulted in lower activity of the transposable elements through DNA methylation in the hybrid. Our data provide insights into the role that changes in gene expression pattern and epigenetic mechanisms contribute to heterosis during early flower development in allopolyploid B. napus., (© 2017 The Authors The Plant Journal © 2017 John Wiley & Sons Ltd.)

Circular RNAs (circRNAs) have been identified in diverse eukaryotic species and are characterized by RNA backsplicing events. Current available methods for circRNA identification are able to determine the start and end locations of circRNAs in the genome but not their full-length sequences. In this study, we developed a method to assemble the full-length sequences of circRNAs using the backsplicing RNA-Seq reads and their corresponding paired-end reads. By applying the method to an rRNA-depleted/RNase R-treated RNA-Seq dataset, we for the first time identified full-length sequences of nearly 3,000 circRNAs in rice. We further showed that alternative circularization of circRNA is a common feature in rice and, surprisingly, found that the junction sites of a large number of rice circRNAs are flanked by diverse non-GT/AG splicing signals while most human exonic circRNAs are flanked by canonical GT/AG splicing signals. Our study provides a method for genome-wide identification of full-length circRNAs and expands our understanding of splicing signals of circRNAs.

Despite the importance of host-microbe interactions in natural ecosystems, agriculture and medicine, the impact of long-term (especially decades or longer) microbial colonization on the dynamics of host genomes is not well understood. The vegetable crop 'Jiaobai' with enlarged edible stems was domesticated from wild Zizania latifolia (Oryzeae) approximately 2000 years ago as a result of persistent infection by a fungal endophyte, Ustilago esculenta. Asexual propagation via infected rhizomes is the only means of Jiaobai production, and the Z. latifolia-endophyte complex has been maintained continuously for two centuries. Here, genomic analysis revealed that cultivated Z. latifolia has a significantly smaller repertoire of immune receptors compared with wild Z. latifolia. There are widespread gene losses/mutations and expression changes in the plant-pathogen interaction pathway in Jiaobai. These results show that continuous long-standing endophyte association can have a major effect on the evolution of the structural and transcriptomic components of the host genome., (© 2015 The Authors The Plant Journal © 2015 John Wiley & Sons Ltd.)

Main Conclusion: Whole-genome re-sequencing of weedy rice from southern China reveals that weedy rice can originate from hybridization of domesticated indica and japonica rice. Weedy rice (Oryza sativa f. spontanea Rosh.), which harbors phenotypes of both wild and domesticated rice, has become one of the most notorious weeds in rice fields worldwide. While its formation is poorly understood, massive amounts of rice genomic data may provide new insights into this issue. In this study, we determined genomes of three weedy rice samples from the lower Yangtze region, China, and investigated their phylogenetics, population structure and chromosomal admixture patterns. The phylogenetic tree and principle component analysis based on 46,005 SNPs with 126 other Oryza accessions suggested that the three weedy rice accessions were intermediate between japonica and indica rice. An ancestry inference study further demonstrated that weedy rice had two dominant genomic components (temperate japonica and indica). This strongly suggests that weedy rice originated from indica-japonica hybridization. Furthermore, 22,443 novel fixed single nucleotide polymorphisms were detected in the weedy genomes and could have been generated after indica-japonica hybridization for environmental adaptation.

Semi-wild soybean is a unique type of soybean that retains both wild and domesticated characteristics, which provides an important intermediate type for understanding the evolution of the subgenus Soja population in the Glycine genus. In this study, a semi-wild soybean line (Maliaodou) and a wild line (Lanxi 1) collected from the lower Yangtze regions were deeply sequenced while nine other semi-wild lines were sequenced to a 3-fold genome coverage. Sequence analysis revealed that (1) no independent phylogenetic branch covering all 10 semi-wild lines was observed in the Soja phylogenetic tree; (2) besides two distinct subpopulations of wild and cultivated soybean in the Soja population structure, all semi-wild lines were mixed with some wild lines into a subpopulation rather than an independent one or an intermediate transition type of soybean domestication; (3) high heterozygous rates (0.19-0.49) were observed in several semi-wild lines; and (4) over 100 putative selective regions were identified by selective sweep analysis, including those related to the development of seed size. Our results suggested a hybridization origin for the semi-wild soybean, which makes a complex Soja population structure.

The lack of a MIRNA set and genome sequence of wild rice (Oryza rufipogon) has prevented us from determining the role of MIRNA genes in rice domestication. In this study, a genome, three small RNA populations and a degradome of O. rufipogon were sequenced by Illumina platform and the expression levels of microRNAs (miRNAs) were investigated by miRNA chips. A de novo O. rufipogon genome was assembled using c. 55× coverage of raw sequencing data and a total of 387 MIRNAs were identified in the O. rufipogon genome based on c. 5.2 million unique small RNA reads from three different tissues of O. rufipogon. Of these, O. rufipogon MIRNAs, 259 were not found in the cultivated rice, suggesting a loss of these MIRNAs in the cultivated rice. We also found that 48 MIRNAs were novel in the cultivated rice, suggesting that they were potential targets of domestication selection. Some miRNAs showed significant expression differences between wild and cultivated rice, suggesting that expression of miRNA could also be a target of domestication, as demonstrated for the miR164 family. Our results illustrated that MIRNA genes, like protein-coding genes, might have been significantly shaped during rice domestication and could be one of the driving forces that contributed to rice domestication., (© 2012 The Authors. New Phytologist © 2012 New Phytologist Trust.)

Modern crops have usually experienced domestication selection and subsequent genetic improvement (post-domestication selection). Chinese waxy maize, which originated from non-glutinous domesticated maize (Zea mays ssp. mays), provides a unique model for investigating the post-domestication selection of maize. In this study, the genetic diversity of six key genes in the starch pathway was investigated in a glutinous population that included 55 Chinese waxy accessions, and a selective bottleneck that resulted in apparent reductions in diversity in Chinese waxy maize was observed. Significant positive selection in waxy (wx) but not amylose extender1 (ae1) was detected in the glutinous population, in complete contrast to the findings in non-glutinous maize, which indicated a shift in the selection target from ae1 to wx during the improvement of Chinese waxy maize. Our results suggest that an agronomic trait can be quickly improved into a target trait with changes in the selection target among genes in a crop pathway.

Background: MicroRNAs (miRNAs) posttranscriptionally down-regulate gene expression by binding target mRNAs. Analysis of the evolution of miRNA binding sites is helpful in understanding the co-evolution between miRNAs and their targets. To understand this process in plants a comparative analysis of miRNA-targeted duplicated gene pairs derived from a well-documented whole genome duplication (WGD) event in combination with a population genetics study of six experimentally validated miRNA binding sites in rice (O. sativa) was carried out., Results: Of the 1,331 pairs of duplicate genes from the WGD, 41 genes (29 pairs) were computationally predicted to be miRNA targets. Sequence substitution analysis indicated that the synonymous substitution rate was significantly lower in the miRNA binding sites than their 5' and 3' flanking regions. Of the 29 duplicated gene pairs, 17 have only one paralog been targeted by a miRNA. This could be due to either gain of a miRNA binding site after the WGD or because one of the duplicated genes has escaped from being a miRNA target after the WGD (loss of miRNA binding site). These possibilities were distinguished by separating miRNAs conserved in both dicots and monocot plants from rice-specific miRNAs and by phylogenetic analysis of miRNA target gene families. The gain/loss rate of miRNA binding sites was estimated to be 3.0 x 10(-9) gain/loss per year. Most (70.6%) of the gains/losses were due to nucleotide mutation. By analysis of cultivated (O. sativa; n = 30) and wild (O. rufipogon; n = 15) rice populations, no segregating site was observed in six miRNA binding sites whereas 0.12-0.20 SNPs per 21-nt or 1.53-1.80 x 10(-3) of the average pairwise nucleotide diversity (pi) were found in their flanking regions., Conclusion: Both molecular evolution and population genetics support the hypothesis that conservation of miRNA binding sites is maintained by purifying selection through elimination of deleterious alleles. Nucleotide mutations play a major role in the gain/loss of miRNA binding sites during evolution.

Rice (Oryza sativa) is one of three predominant grain crops, and its nuclear and organelle genomes have been sequenced. Following genome analysis revealed many exchanges of DNA sequences between the nuclear and organelle genomes. In this study, a total of 45 chloroplast DNA insertions more than 2 kb in length were detected in rice nuclear genome. A homologous recombination mechanism is expected for those chloroplast insertions with high similarity between their flanking sequences. Only five chloroplast insertions with high sequence similarity between two flanking sequences from an insertion were found in the 45 insertions, suggesting that rice might follow the non-homologous end-joining (NHEJ) repair of double-stranded breaks mechanism, which is suggested to be common to all eukaryotes. Our studies indicate that the most chloroplast insertions occurred at a nuclear region characterized by a sharp change of repetitive sequence density. One potential explanation is that regions such as this might be susceptible target sites or "hotspots" of DNA damage. Our results also suggest that the insertion of retrotransposon elements or non-chloroplast DNA into chloroplast DNA insertions may contribute significantly to their fragmentation process. Moreover, based on chloroplast insertions in nuclear genomes of two subspecies (indica and japonica) of cultivated rice, our results strongly suggest that they diverged during 0.06-0.22 million years ago.

Moso bamboo (Phyllostachys pubescens) is one of the world's most important bamboo species. It has the largest area of all planted bamboo--over two-thirds of the total bamboo forest area--and the highest economic value in China. Moso bamboo is a tetraploid (4x=48) and a special member of the grasses family. Although several genomes have been sequenced or are being sequenced in the grasses family, we know little about the genome of the bambusoids (bamboos). In this study, the moso bamboo genome size was estimated to be about 2034 Mb by flow cytometry (FCM), using maize (cv. B73) and rice (cv. Nipponbare) as internal references. The rice genome has been sequenced and the maize genome is being sequenced. We found that the size of the moso bamboo genome was similar to that of maize but significantly larger than that of rice. To determine whether the bamboo genome had a high proportion of repeat elements, similar to that of the maize genome, approximately 1000 genome survey sequences (GSS) were generated. Sequence analysis showed that the proportion of repeat elements was 23.3% for the bamboo genome, which is significantly lower than that of the maize genome (65.7%). The bamboo repeat elements were mainly Gypsy/DIRS1 and Ty1/Copia LTR retrotransposons (14.7%), with a few DNA transposons. However, more genomic sequences are needed to confirm the above results due to several factors, such as the limitation of our GSS data. This study is the first to investigate sequence composition of the bamboo genome. Our results are valuable for future genome research of moso and other bamboos.