Signatures of differential selection in chloroplast genome between japonica and indica

Comparative Analysis of Chloroplast Pan-Genomes and Transcriptomics Reveals Cold Adaptation in Medicago sativa

Alfalfa (Medicago sativa) is a perennial forage legume that is widely distributed all over the world; therefore, it has an extremely complex genetic background. Though population structure and phylogenetic studies have been conducted on a large group of alfalfa nuclear genomes, information about the chloroplast genomes is still lacking. Chloroplast genomes are generally considered to be conservative and play an important role in population diversity analysis and species adaptation in plants. Here, 231 complete alfalfa chloroplast genomes were successfully assembled from 359 alfalfa resequencing data, on the basis of which the alfalfa chloroplast pan-genome was constructed. We investigated the genetic variations of the alfalfa chloroplast genome through comparative genomic, genetic diversity, phylogenetic, population genetic structure, and haplotype analysis. Meanwhile, the expression of alfalfa chloroplast genes under cold stress was explored through transcriptome analysis. As a result, chloroplast genomes of 231 alfalfa lack an IR region, and the size of the chloroplast genome ranges from 125,192 bp to 126,105 bp. Using population structure, haplotypes, and construction of a phylogenetic tree, it was found that alfalfa populations could be divided into four groups, and multiple highly variable regions were found in the alfalfa chloroplast genome. Transcriptome analysis showed that tRNA genes were significantly up-regulated in the cold-sensitive varieties, while rps7, rpl32, and ndhB were down-regulated, and the editing efficiency of ycf1, ycf2, and ndhF was decreased in the cold-tolerant varieties, which may be due to the fact that chloroplasts store nutrients through photosynthesis to resist cold. The huge number of genetic variants in this study provide powerful resources for molecular markers.

Chloroplast Pan-Genomes and Comparative Transcriptomics Reveal Genetic Variation and Temperature Adaptation in the Cucumber

Although whole genome sequencing, genetic variation mapping, and pan-genome studies have been done on a large group of cucumber nuclear genomes, organelle genome information is largely unclear. As an important component of the organelle genome, the chloroplast genome is highly conserved, which makes it a useful tool for studying plant phylogeny, crop domestication, and species adaptation. Here, we have constructed the first cucumber chloroplast pan-genome based on 121 cucumber germplasms, and investigated the genetic variations of the cucumber chloroplast genome through comparative genomic, phylogenetic, haplotype, and population genetic structure analysis. Meanwhile, we explored the changes in expression of cucumber chloroplast genes under high- and low-temperature stimulation via transcriptome analysis. As a result, a total of 50 complete chloroplast genomes were successfully assembled from 121 cucumber resequencing data, ranging in size from 156,616-157,641 bp. The 50 cucumber chloroplast genomes have typical quadripartite structures, consisting of a large single copy (LSC, 86,339-86,883 bp), a small single copy (SSC, 18,069-18,363 bp), and two inverted repeats (IRs, 25,166-25,797 bp). Comparative genomic, haplotype, and population genetic structure results showed that there is more genetic variation in Indian ecotype cucumbers compared to other cucumber cultivars, which means that many genetic resources remain to be explored in Indian ecotype cucumbers. Phylogenetic analysis showed that the 50 cucumber germplasms could be classified into 3 types: East Asian, Eurasian + Indian, and Xishuangbanna + Indian. The transcriptomic analysis showed that matK were significantly up-regulated under high- and low-temperature stresses, further demonstrating that cucumber chloroplasts respond to temperature adversity by regulating lipid metabolism and ribosome metabolism. Further, accD has higher editing efficiency under high-temperature stress, which may contribute to the heat tolerance. These studies provide useful insight into genetic variation in the chloroplast genome, and established the foundation for exploring the mechanisms of temperature-stimulated chloroplast adaptation.

Multiple domestication events explain the origin of Gossypium hirsutum landraces in Mexico

Several Mesoamerican crops constitute wild-to-domesticated complexes generated by multiple initial domestication events, and continuous gene flow among crop populations and between these populations and their wild relatives. It has been suggested that the domestication of cotton (Gossypium hirsutum) started in the northwest of the Yucatán Peninsula, from where it spread to other regions inside and outside of Mexico. We tested this hypothesis by assembling chloroplast genomes of 23 wild, landraces, and breeding lines (transgene-introgressed and conventional). The phylogenetic analysis showed that the evolutionary history of cotton in Mexico involves multiple events of introgression and genetic divergence. From this, we conclude that Mexican landraces arose from multiple wild populations. Our results also revealed that their structural and functional chloroplast organizations had been preserved. However, genetic diversity decreases as a consequence of domestication, mainly in transgene-introgressed (TI) individuals (π = 0.00020, 0.00001, 0.00016, 0, and 0, of wild, TI-wild, landraces, TI-landraces, and breeding lines, respectively). We identified homologous regions that differentiate wild from domesticated plants and indicate a relationship among the samples. A decrease in genetic diversity associated with transgene introgression in cotton was identified for the first time, and our outcomes are therefore relevant to both biosecurity and agrobiodiversity conservation.

Slowing virulence adaptation in Asian rice planthoppers through migration-based deployment of resistance genes

Advances in molecular biology have accelerated rice breeding for resistance to Asian planthoppers. However, experience shows that planthoppers quickly adapt to resistance in tropical overwintering areas. With only limited sources available, the large-scale deployment of resistance genes can rapidly reduce the utility of these public goods. Planthoppers that migrate from tropical to temperate Asia carry virulence against many resistance genes, but adapt more slowly to resistant rice in temperate regions. Therefore, by restricting deployment of selected genes to temperate regions, virulence-adaptation rates and the volume of migrants returning to overwintering sites could be reduced. The current open exchange of breeding materials throughout Asia urgently requires an international, multidisciplinary, stakeholder coalition to promote a more sustainable deployment of planthopper-resistant rice.

Development of an inclusive 580K SNP array and its application for genomic selection and genome-wide association studies in rice

Rice is a globally cultivated crop and is primarily a staple food source for more than half of the world's population. Various single-nucleotide polymorphism (SNP) arrays have been developed and utilized as standard genotyping methods for rice breeding research. Considering the importance of SNP arrays with more inclusive genetic information for GWAS and genomic selection, we integrated SNPs from eight different data resources: resequencing data from the Korean World Rice Collection (KRICE) of 475 accessions, 3,000 rice genome project (3 K-RGP) data, 700 K high-density rice array, Affymetrix 44 K SNP array, QTARO, Reactome, and plastid and GMO information. The collected SNPs were filtered and selected based on the breeder's interest, covering all key traits or research areas to develop an integrated array system representing inclusive genomic polymorphisms. A total of 581,006 high-quality SNPs were synthesized with an average distance of 200 bp between adjacent SNPs, generating a 580 K Axiom Rice Genotyping Chip (580 K _ KNU chip). Further validation of this array on 4,720 genotypes revealed robust and highly efficient genotyping. This has also been demonstrated in genome-wide association studies (GWAS) and genomic selection (GS) of three traits: clum length, heading date, and panicle length. Several SNPs significantly associated with cut-off, -log10 p-value >7.0, were detected in GWAS, and the GS predictabilities for the three traits were more than 0.5, in both rrBLUP and convolutional neural network (CNN) models. The Axiom 580 K Genotyping array will provide a cost-effective genotyping platform and accelerate rice GWAS and GS studies.

Selection and Validation of 48 KASP Markers for Variety Identification and Breeding Guidance in Conventional and Hybrid Rice (Oryza sativa L.)

BACKGROUND

Breeding of conventional and hybrid rice (Oryza sativa L.) have solved hunger problems and increased farmers' income in the world. Molecular markers have been widely used in marker-assisted breeding and identification of larger numbers of different bred varieties in the past decades. The recently developed SNP markers are applied for more stable and detectable compared with other markers. But the cost of genotyping lots SNPs is high. So, it is essential to select less representative SNPs and inexpensive detecting methods to lower the cost and accelerate variety identification and breeding process. KASP (Kompetitive Allele-Specific PCR) is a flexible method to detect the SNPs, and large number of KASP markers have been widely used in variety identification and breeding. However, the ability of less KASP markers on massive variety identification and breeding remains unknown.

RESULTS

Here, 48 KASP markers were selected from 378 markers to classify and analyze 518 varieties including conventional and hybrid rice. Through analyzing the population structure, the 48 markers could almost represent the 378 markers. In terms of variety identification, the 48 KASP markers had a 100% discrimination rate in 53 conventional indica varieties and 193 hybrid varieties, while they could distinguish 89.1% conventional japonica rice from different breeding institutes. Two more markers added would increase the ratio from 68.38 to 77.94%. Additionally, the 48 markers could be used for classification of subpopulations in the bred variety. Also, 8 markers had almost completely different genotypes between japonica and indica, and 3 markers were found to be very important for japonica hybrid rice. In hybrid varieties, the heterozygosity of chromosomes 3, 6 and 11 was relatively higher than others.

CONCLUSIONS

Our results showed that 48 KASP markers could be used to identify rice varieties, and the panel we tested could provide a database for breeders to identify new breeding lines. Also, the specific markers we found were useful for marker-assisted breeding in rice, including conventional and hybrid.

Chloroplast Genomic Resources and Genetic Divergence of Endangered Species Bretschneidera sinensis (Bretschneideraceae)

Bretschneidera sinensis is an endangered woody species found in East and South China. Comprehensive intraspecies chloroplast genome studies have demonstrated novel genetic resources to assess the genetic variation and diversity of this species. Using genome skimming method, we assembled the whole chloroplast genome of 12 genotypes of B. sinensis from different geographical locations, covering most wild populations. The B. sinensis chloroplast genome size ranged from 158,959 to 159,045 base pairs (bp) and displayed a typical circular quadripartite structure. Comparative analyses of 12 B. sinensis chloroplast genome revealed 33 polymorphic simple sequence repeats (SSRs), 105 polymorphic single nucleotide polymorphisms (SNPs), and 55 indels. Phylogenetic analysis showed that the 12 genotypes were grouped into 2 branches, which is consistent with the geographical distribution (Eastern clade and Western clade). Divergence time estimates showed that the two clades were divergent from 0.6 Ma in the late Pleistocene. Ex situ conservation is essential for this species. In this study, we identified SNPs, indels, and microsatellites of B. sinensis by comparative analyses of chloroplast genomes and determined genetic variation between populations using these genomic markers. Chloroplast genomic resources are also important for further domestication, population genetic, and phylogenetic analysis, possibly in combination with molecular markers of mitochondrial and/or nuclear genomes.

A large-scale population based organelle pan-genomes construction and phylogeny analysis reveal the genetic diversity and the evolutionary origins of chloroplast and mitochondrion in Brassica napus L

BACKGROUND

Allotetraploid oilseed rape (Brassica napus L.) is an important worldwide oil-producing crop. The origin of rapeseed is still undetermined due to the lack of wild resources. Despite certain genetic architecture and phylogenetic studies have been done focus on large group of Brassica nuclear genomes, the organelle genomes information under global pattern is largely unknown, which provide unique material for phylogenetic studies of B. napus. Here, based on de novo assemblies of 1,579 B. napus accessions collected globally, we constructed the chloroplast and mitochondrial pan-genomes of B. napus, and investigated the genetic diversity, phylogenetic relationships of B. napus, B. rapa and B. oleracea.

RESULTS

Based on mitotype-specific markers and mitotype-variant ORFs, four main cytoplasmic haplotypes were identified in our groups corresponding the nap, pol, ole, and cam mitotypes, among which the structure of chloroplast genomes was more conserved without any rearrangement than mitochondrial genomes. A total of 2,092 variants were detected in chloroplast genomes, whereas only 326 in mitochondrial genomes, indicating that chloroplast genomes exhibited a higher level of single-base polymorphism than mitochondrial genomes. Based on whole-genome variants diversity analysis, eleven genetic difference regions among different cytoplasmic haplotypes were identified on chloroplast genomes. The phylogenetic tree incorporating accessions of the B. rapa, B. oleracea, natural and synthetic populations of B. napus revealed multiple origins of B. napus cytoplasm. The cam-type and pol-type were both derived from B. rapa, while the ole-type was originated from B. oleracea. Notably, the nap-type cytoplasm was identified in both the B. rapa population and the synthetic B. napus, suggesting that B. rapa might be the maternal ancestor of nap-type B. napus.

CONCLUSIONS

The phylogenetic results provide novel insights into the organelle genomic evolution of Brassica species. The natural rapeseeds contained at least four cytoplastic haplotypes, of which the predominant nap-type might be originated from B. rapa. Besides, the organelle pan-genomes and the overall variation data offered useful resources for analysis of cytoplasmic inheritance related agronomical important traits of rapeseed, which can substantially facilitate the cultivation and improvement of rapeseed varieties.

Exploring the evolutionary characteristics between cultivated tea and its wild relatives using complete chloroplast genomes

Background

Cultivated tea is one of the most important economic and ecological trees distributed worldwide. Cultivated tea suffer from long-term targeted selection of traits and overexploitation of habitats by human beings, which may have changed its genetic structure. The chloroplast is an organelle with a conserved cyclic genomic structure, and it can help us better understand the evolutionary relationship of Camellia plants.

Results

We conducted comparative and evolutionary analyses on cultivated tea and wild tea, and we detected the evolutionary characteristics of cultivated tea. The chloroplast genome sizes of cultivated tea were slightly different, ranging from 157,025 to 157,100 bp. In addition, the cultivated species were more conserved than the wild species, in terms of the genome length, gene number, gene arrangement and GC content. However, comparing Camellia sinensis var. sinensis and Camellia sinensis var. assamica with their cultivars, the IR length variation was approximately 20 bp and 30 bp, respectively. The nucleotide diversity of 14 sequences in cultivated tea was higher than that in wild tea. Detailed analysis on the genomic variation and evolution of Camellia sinensis var. sinensis cultivars revealed 67 single nucleotide polymorphisms (SNPs), 46 insertions/deletions (indels), and 16 protein coding genes with nucleotide substitutions, while Camellia sinensis var. assamica cultivars revealed 4 indels. In cultivated tea, the most variable gene was ycf1. The largest number of nucleotide substitutions, five amino acids exhibited site-specific selection, and a 9 bp sequence insertion were found in the Camellia sinensis var. sinensis cultivars. In addition, phylogenetic relationship in the ycf1 tree suggested that the ycf1 gene has diverged in cultivated tea. Because C. sinensis var. sinensis and its cultivated species were not tightly clustered.

Conclusions

The cultivated species were more conserved than the wild species in terms of architecture and linear sequence order. The variation of the chloroplast genome in cultivated tea was mainly manifested in the nucleotide polymorphisms and sequence insertions. These results provided evidence regarding the influence of human activities on tea.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01800-1.

The complete chloroplast genome sequence of Oryza sativa Temperate japonica

In this study, the complete chloroplast genome of O. sativa Temperate Japonica YunJing-24 was assembled using Illumina sequencing data. The complete chloroplast (cp) genome is 134,556 bp in length, including a pair of invert repeats (IRA and IRB) regions of 20,797 bp, large single-copy (LSC) region of 80,615 bp, and small single-copy (SSC) region of 12,347 bp. A total of 129 genes were predicted in the genome, including 87 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. Phylogenetic analysis confirmed the phylogenetic relationship between O. sativa Temperate Japonica and other representative species.

Two divergent chloroplast genome sequence clades captured in the domesticated rice gene pool may have significance for rice production

BACKGROUND

The whole chloroplast genomes of 3018 rice genotypes were assembled from available sequence data by alignment with a reference rice chloroplast genome sequence, providing high quality chloroplast genomes for analysis of diversity on a much larger scale than in any previous plant study.

RESULTS

Updated annotation of the chloroplast genome identified 13 more tRNA genes and 30 more introns and defined the function of more of the genes. Domesticated rice had chloroplast genomes that were distinct from those in wild relatives. Analysis confirms an Australian chloroplast clade as a sister to the domesticated clade. All domesticated rice genotypes could be assigned to one of two main clades suggesting the domestication of two distinct maternal genome clades that diverged long before domestication. These clades were very distinct having 4 polymorphisms between all 1486 accession in clade A and all 1532 accessions in clade B. These would result in expression of 3 proteins with altered amino acid sequences and a tRNA with an altered sequence and may be associated with adaptive evolution of the two chloroplast types. Diversity within these pools may have been captured during domestication with subclades enriched in specific groups such as basmati, tropical japonica and temperate japonica. However the phylogenies of the chloroplast and nuclear genomes differed possibly due to modern rice breeding and reticulate evolution prior to domestication. Indica and aus genotypes were common in both chloroplast clades while japonica genotypes were more likely to be found in the same clade (cladeB).

CONCLUSIONS

The different evolutionary paths of the cytoplasmic and nuclear genomes of rice have resulted in the presence of apparently functional chloroplast genome diversity and the implications for rice crop performance require further investigation.

Development of nuclear SSR and chloroplast genome markers in diverse Liriodendron chinense germplasm based on low-coverage whole genome sequencing

Background

Liriodendron chinense ranges widely in subtropical China and northern Vietnam; however, it inhabits several small, isolated populations and is now an endangered species due to its limited seed production. The objective of this study was to develop a set of nuclear SSR (simple sequence repeats) and multiple chloroplast genome markers for genetic studies in L. chinense and their characterization in diverse germplasm.

Results

We performed low-coverage whole genome sequencing of the L. chinense from four genotypes, assembled the chloroplast genome and identified nuclear SSR loci by searching in contigs for SSR motifs. Comparative analysis of the four chloroplast genomes of L. chinense revealed 45 SNPs, 17 indels, 49 polymorphic SSR loci, and five small inversions. Most chloroplast intraspecific polymorphisms were located in the interspaces of single-copy regions. In total, 6147 SSR markers were isolated from low-coverage whole genome sequences. The most common SSR motifs were dinucleotide (70.09%), followed by trinucleotide motifs (23.10%). The motif AG/TC (33.51%) was the most abundant, followed by TC/AG (25.53%). A set of 13 SSR primer combinations were tested for amplification and their ability to detect polymorphisms in a set of 109 L. chinense individuals, representing distinct varieties or germplasm. The number of alleles per locus ranged from 8 to 28 with an average of 21 alleles. The expected heterozygosity (H_e) varied from 0.19 to 0.93 and the observed heterozygosity (H_o) ranged from 0.11 to 0.79.

Conclusions

The genetic resources characterized and tested in this study provide a valuable tool to detect polymorphisms in L. chinense for future genetic studies and breeding programs.