The History and Diversity of Rice Domestication as Resolved From 1464 Complete Plastid Genomes

Plant organellar genomes: much done, much more to do

Plastids and mitochondria are the only organelles that possess genomes of endosymbiotic origin. In recent decades, advances in sequencing technologies have contributed to a meteoric rise in the number of published organellar genomes, and have revealed greatly divergent evolutionary trajectories. In this review, we quantify the abundance and distribution of sequenced plant organellar genomes across the plant tree of life. We compare numerous genomic features between the two organellar genomes, with an emphasis on evolutionary trajectories, transfers, the current state of organellar genome editing by transcriptional activator-like effector nucleases (TALENs), transcription activator-like effector (TALE)-mediated deaminase, and clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein (Cas), as well as genetic transformation. Finally, we propose future research to understand these different evolutionary trajectories, and genome-editing strategies to promote functional studies and eventually improve organellar genomes.

The pan-plastome of Prunus mume: insights into Prunus diversity, phylogeny, and domestication history

Backgrounds

Prunus mume in the Rosaceae and commonly referred to as mei or Chinese plum is widely used as a traditional ornamental flowering plant and fruit tree in China. Although some population and genetic analyses have been conducted for this species, no extensive comparisons of genetic variation from plastomes have yet been investigated.

Methods

We de novo assembled a total of 322 complete P. mume plastomes in this study and did a series of comparative analyses to better resolve pan-plastomic patterns of P. mume. To determine the phylogeny and domestication history of this species, we reconstructed the phylogenetic tree of Prunus genus, and resolved the population structure of P. mume. We also examined the nucleotide variation of P. mume to find potential DNA barcodes.

Results

The assembled plastomes exhibited a typical quadripartite structure and ranged from 157,871 bp to 158,213 bp in total size with a GC content ranging from 36.73 to 36.75%. A total of 112 unique genes were identified. Single nucleotide variants (SNVs) were the most common variants found among the plastomes, followed by nucleotide insertions/deletions (InDels), and block substitutions with the intergenic spacer (IGS) regions containing the greatest number of variants. From the pan-plastome data six well-supported genetic clusters were resolved using multiple different population structure analyses. The different cultivars were unevenly distributed among multiple clades. We also reconstructed a phylogeny for multiple species of Prunus to better understand genus level diversity and history from which a complex introgressive relationship between mei and other apricots/plums was resolved.

Conclusion

This study constructed the pan-plastome of P. mume, which indicated the domestication of P. mume involved multiple genetic origins and possible matrilineal introgression from other species. The phylogenetic analysis in Prunus and the population structure of P. mume provide an important maternal history for Prunus and the groundwork for future studies on intergenomic sequence transfers, cytonuclear incompatibility, and conservation genetics.

The pan-plastome of Hemerocallis citrina reveals new insights into the genetic diversity and cultivation history of an economically important food plant

BACKGROUND

Hemerocallis citrina Baroni (Huang hua cai in Chinese) is a perennial herbaceous plant grown for its flower buds that are eaten fresh or dried and is known as the vegetarian three treasures. The nuclear genome of H. citrina has been reported, but the intraspecific variation of the plastome (plastid genome) has not yet been studied. Therefore, the panplastome of this species collected from diverse locations is reported here for the first time.

RESULTS

In this study, 65 H. citrina samples were resequenced, de novo assembled, and aligned with the published plastome of H. citrina to resolve the H. citrina panplastome. The sizes of the 65 newly assembled complete plastomes of H. citrina ranged from 156,048 bp to 156,263 bp, and the total GC content ranged from 37.31 to 37.34%. The structure of the complete plastomes showed a typical tetrameric structure, including a large single copy (LSC), a small single copy (SSC), and a pair of inverted repeat regions (IRA and IRB). Many nucleotide variants were identified between plastomes, among which the variants in the intergenic spacer region were the most abundant, with the highest number of variants concentrated in the LSC region. Based on the phylogenetic tree constructed using the ML method, population structure analysis, and principal component analysis (PCA), the panplastome data were subdivided into five genetic clusters. The C5 genetic cluster was mostly represented by samples from Qidong, Hunan Province, while samples from Shanxi and Shaanxi Provinces were classified into the C4 genetic cluster. The greatest genetic diversity was found in the C1 genetic cluster, and the greatest genetic distance between any two clusters was found between the C4 and C5 clusters.

CONCLUSION

The resolution of the panplastome and the analysis of the population structure of H. citrina plastomes provide important data for future breeding projects and germplasm preservation.

A derived weedy rice × ancestral cultivar cross identifies evolutionarily relevant weediness QTLs

Weedy rice (Oryza spp.) is a weedy relative of the cultivated rice that competes with the crop and causes significant production loss. The BHA (blackhull awned) US weedy rice group has evolved from aus cultivated rice and differs from its ancestors in several important weediness traits, including flowering time, plant height and seed shattering. Prior attempts to determine the genetic basis of weediness traits in plants using linkage mapping approaches have not often considered weed origins. However, the timing of divergence between crossed parents can affect the detection of quantitative trait loci (QTL) relevant to the evolution of weediness. Here, we used a QTL-seq approach that combines bulked segregant analysis and high-throughput whole genome resequencing to map the three important weediness traits in an F2 population derived from a cross between BHA weedy rice with an ancestral aus cultivar. We compared these QTLs with those previously detected in a cross of BHA with a more distantly related crop, indica. We identified multiple QTLs that overlapped with regions under selection during the evolution of weedy BHA rice and some candidate genes possibly underlying the evolution weediness traits in BHA. We showed that QTLs detected with ancestor-descendant crosses are more likely to be involved in the evolution of weediness traits than those detected from crosses of more diverged taxa.

Multiple configurations of the plastid and mitochondrial genomes of Caragana spinosa

MAIN CONCLUSION

In this study, we assembled the complete plastome and mitogenome of Caragana spinosa and explored the multiple configurations of the organelle genomes. Caragana spinosa belongs to the Papilionoidea subfamily and has significant pharmaceutical value. To explore the possible interaction between the organelle genomes, we assembled and analyzed the plastome and mitogenome of C. spinosa using the Illumina and Nanopore DNA sequencing data. The plastome of C. spinosa was 129,995 bp belonging to the inverted repeat lacking clade (IRLC), which contained 77 protein-coding genes, 29 tRNA genes, and four rRNA genes. The mitogenome was 378,373 bp long and encoded 54 unique genes, including 33 protein-coding, three ribosomal RNA (rRNA), and 18 transfer RNA (tRNA) genes. In addition to the single circular conformation, alternative conformations mediated by one and four repetitive sequences in the plastome and mitogenome were identified and validated, respectively. The inverted repeat (PDR12, the 12th dispersed repeat sequence in C. spinosa plastome) of plastome mediating recombinant was conserved in the genus Caragana. Furthermore, we identified 14 homologous fragments by comparing the sequences of mitogenome and plastome, including eight complete tRNA genes. A phylogenetic analysis of protein-coding genes extracted from the plastid and mitochondrial genomes revealed congruent topologies. Analyses of sequence divergence found one intergenic region, trnN-GUU-ycf1, exhibiting a high degree of variation, which can be used to develop novel molecular markers to distinguish the nine Caragana species accurately. This plastome and mitogenome of C. spinosa could provide critical information for the molecular breeding of C. spinosa and be used as a reference genome for other species of Caragana. In this study, we assembled the complete plastome and mitogenome of Caragana spinosa and explored the multiple configurations of the organelle genomes.

The pan-plastome of tartary buckwheat (fagopyrum tataricum): key insights into genetic diversity and the history of lineage divergence

BACKGROUND

Tartary buckwheat (Fagopyrum tataricum) is an important food and medicine crop plant, which has been cultivated for 4000 years. A nuclear genome has been generated for this species, while an intraspecific pan-plastome has yet to be produced. As such a detailed understanding of the maternal genealogy of Tartary buckwheat has not been thoroughly investigated.

RESULTS

In this study, we de novo assembled 513 complete plastomes of Fagopyrum and compared with 8 complete plastomes of Fagopyrum downloaded from the NCBI database to construct a pan-plastome for F. tartaricum and resolve genomic variation. The complete plastomes of the 513 newly assembled Fagopyrum plastome sizes ranged from 159,253 bp to 159,576 bp with total GC contents ranged from 37.76 to 37.97%. These plastomes all maintained the typical quadripartite structure, consisting of a pair of inverted repeat regions (IRA and IRB) separated by a large single copy region (LSC) and a small single copy region (SSC). Although the structure and gene content of the Fagopyrum plastomes are conserved, numerous nucleotide variations were detected from which population structure could be resolved. The nucleotide variants were most abundant in the non-coding regions of the genome and of those the intergenic regions had the most. Mutational hotspots were primarily found in the LSC regions. The complete 521 Fagopyrum plastomes were divided into five genetic clusters, among which 509 Tartary buckwheat plastomes were divided into three genetic clusters (Ft-I/Ft-II/Ft-III). The genetic diversity in the Tartary buckwheat genetic clusters was the greatest in Ft-III, and the genetic distance between Ft-I and Ft-II was the largest. Based on the results of population structure and genetic diversity analysis, Ft-III was further subdivided into three subgroups Ft-IIIa, Ft-IIIb, and Ft-IIIc. Divergence time estimation indicated that the genera Fagopyrum and Rheum (rhubarb) shared a common ancestor about 48 million years ago (mya) and that intraspecies divergence in Tartary buckwheat began around 0.42 mya.

CONCLUSIONS

The resolution of pan-plastome diversity in Tartary buckwheat provides an important resource for future projects such as marker-assisted breeding and germplasm preservation.

Harnessing intra-varietal variation for agro-morphological and nutritional traits in a popular rice landrace for sustainable food security in tropical islands

Introduction

Rice crop meets the calorie and nutritional requirements of a larger segment of the global population. Here, we report the occurrence of intra-varietal variation in a popular rice landrace C14-8 traditionally grown under the geographical isolation of the Andaman Islands.

Methods

Based on grain husk color, four groups were formed, wherein the extent of intra-varietal variation was studied by employing 22 agro-morphological and biochemical traits.

Results

Among the traits studied, flavonoid and anthocyanin contents and grain yield exhibited a wider spectrum of variability due to more coefficients of variation (>25%). The first five principal components (PCs) of principal components analysis explained a significant proportion of the variation (91%) and the first two PCs explained 63.3% of the total variation, with PC1 and PC2 explaining 35.44 and 27.91%, respectively. A total of 50 highly variable SSR (HvSSR) markers spanning over 12 chromosomes produced 314 alleles, which ranged from 1 to 15 alleles per marker, with an average of 6.28. Of the 314 alleles, 64 alleles were found to be rare among the C14-8 selections. While 62% of HvSSR markers exhibited polymorphism among the C14-8 population, chromosomes 2, 7, 9, and 11 harbored the most polymorphic loci. The group clustering of the selections through HvSSR markers conformed to the grouping based on grain husk coloration.

Discussion

Our studies on the existence and pertinence of intra-varietal variations are expected to be of significance in the realms of evolutionary biology and sustainable food and nutritional security under the changing climate.

Long non-coding RNAs as the regulatory hubs in rice response to salt stress

Salinity seriously constrains growth and fertility of rice worldwide. Long non-coding RNAs (lncRNAs) play crucial roles in plant abiotic stress response. However, salt responsive lncRNAs are poorly understood in rice. Herein, salt responsive lncRNAs (DE-lncRNAs) were identified in FL478 (salt tolerant) compared to its susceptible parent (IR29) using RNA-seq in root tissues at seedling stage. In FL478 and IR29, 8724 and 9235 transcripts with length of > 200 bp were nominated as potential lncRNAs, respectively. Rigorous filtering left four (in FL478) and nine (in IR29) DE-lncRNAs with only 2 DE-lncRNAs in common. ATAC-seq data showed that the genomic regions of all four lncRNAs in FL478 and 6/9 in IR29 are significantly accessible for transcription. Weighted correlation network analysis (WGCNA) revealed that lncRNA.2-FL was highly correlated with 173 mRNAs as trans-targets and a gene encoding pentatricopeptide repeat (PPR) protein was predicted as cis-target of lncRNA.2-FL. In silico mutagenesis analysis proposed the same transcription factor binding sites (TFBSs) in vicinity of the trans- and cis-regulatory target genes of lncRNA.2-FL, which significantly affect their transcription start site (TSS). This study provides new insights into involvement of the DE-lncRNAs in rice response to salt stress. Among them, lncRNA.2-FL may play a significant regulatory role in the salt stress tolerance of FL478.

Genetic Diversity of Transcription Factor Genes in Triticum and Mining for Promising Haplotypes for Beneficial Agronomic Traits

Transcription factor (TF) is a class of the sequence-specific DNA-binding proteins that modulate the transcription of target genes, and thus regulate their expressions. Variations in TF are the crucial determinants for phenotypic traits. Although much progress has been made in the functions of TF genes in wheat, one of the most important staple crops globally, the diversity of TF genes in wheat and its progenitors are not well understood, especially the agronomically promising haplotypes have not yet been characterized. Here, we identified a total of 6,023 TF genes from hexaploid wheat through a genome-search method and classified them into 59 gene families based on the conserved domain. The characteristics and dN/dS values of these genes showed evidently selective effects. Based on re-sequencing data, we found a strong genetic bottleneck among these TF genes on A and D subgenomes while no found in B subgenome during wheat domestication. Combined with selective signals and known QTLs on the whole genome, 21 TF genes were preliminarily found to be associated with yield-related traits. The haplotype frequency of these TF genes was further investigated in bread wheat and its progenitors and 13 major haplotypes were the casual loci related to key traits. Finally, the tissue-specific TF genes were also identified using RNA-seq analysis. This study provided insights into the diversity and evolution of TF genes and the identified TF genes and excellent haplotypes associating with traits will contribute to wheat genetic improvement.

Interspecific Hybridization Is an Important Driving Force for Origin and Diversification of Asian Cultivated Rice Oryza sativa L

As one of the most important crops, Asian cultivated rice has evolved into a complex group including several subgroups adapting various eco-climate-systems around the globe. Here, we pictured a comprehensive view of its original domestication, divergences, and the origin of different subgroups by integrating agriculture, archeology, genetics, nuclear, and cytoplasm genome results. Then, it was highlighted that interspecific hybridization-introgression has played important role in improving the genetic diversity and adaptation of Oryza sativa during its evolution process. Natural hybridization-introgression led to the origin of indica, aus, and basmatic subgroups, which adapted to changing cultivated environments, and produced feral weedy rice coexisting and competing with cultivars under production management. Artificial interspecific hybridization-introgression gained several breakthroughs in rice breeding, such as developing three-line hybrid rice, new rice for Africa (NERICA), and some important pest and disease resistance genes in rice genetic improvement, contributing to the stable increase of rice production to meet the expanding human population. We proposed a series to exploit the virtues of hybridization-introgression in the genetic improvement of Asian cultivated rice. But some key issues such as reproductive barriers especially hybrid sterility should be investigated further, which are conducive to gene exchange between cultivated rice and its relatives, and even is beneficial to exploiting interspecific hybrid vigor. New technologies help introduce favorable genes from distant wild species to Asian cultivated rice, such as transgenic and genome editing systems. Rising introgression lines in a wider range with multi-donor benefits allele mining, understanding genetic network of rice growth and development, yield formation, and environmental adaptation. Then, integration of new tools and interspecific hybridization can be a future direction to develop more usable breeding populations which can make Asian cultivated rice more resilient to the changing climate and world.

Organization, Phylogenetic Marker Exploitation, and Gene Evolution in the Plastome of Thalictrum (Ranunculaceae)

Thalictrum is a phylogenetically and economically important genus in the family Ranunculaceae, but is also regarded as one of the most challengingly difficult in plants for resolving the taxonomical and phylogenetical relationships of constituent taxa within this genus. Here, we sequenced the complete plastid genomes of two Thalictrum species using Illumina sequencing technology via de novo assembly. The two Thalictrum plastomes exhibited circular and typical quadripartite structure that was rather conserved in overall structure and the synteny of gene order. By updating the previously reported plastome annotation of other nine Thalictrum species, we found that the expansion or contraction of the inverted repeat region affect the boundary of the single-copy regions in Thalictrum plastome. We identified eight highly variable noncoding regions-infA-rps8, ccsA-ndhD, trnSUGA-psbZ, trnHGUG-psbA, rpl16-rps3, ndhG-ndhI, ndhD-psaC, and ndhJ-ndhK-that can be further used for molecular identification, phylogenetic, and phylogeographic in different species. Selective pressure and codon usage bias of all the plastid coding genes were also analyzed for the 11 species. Phylogenetic relationships showed Thalictrum is monophyly and divided into two major clades based on 11 Thalictrum plastomes. The availability of these plastomes offers valuable genetic information for accurate identification of species and taxonomy, phylogenetic resolution, and evolutionary studies of Thalictrum, and should assist with exploration and utilization of Thalictrum plants.