Twelve complete chloroplast genomes of wild peanuts: great genetic resources and a better understanding of Arachis phylogeny

The complete chloroplast genome of Arachis lutescens Krapov. & Rigoni (Fabaceae)

Arachis lutescens Krapov. & Rigoni 1958 is an important species due to their potentially extensive applications for cultivated peanut breeding. The whole chloroplast genome of A. lutescens was successfully assembled and annotated for the first time. The complete chloroplast genome of A. lutescens is a typically circular structure of 156,398 bp with a GC content of 36.3%. It comprises a large single-copy (LSC) region of 85,950 bp, a small single-copy (SSC) region of 18,800 bp, and two inverted repeat regions (IRs) of 25,824 bp, each. The plastome of A. lutescens contains a total of 125 genes, including 81 protein-coding genes, 36 tRNAs, and eight rRNAs. The phylogenetic analysis strongly supports the close relationship between A. lutescens and cultivated peanut clades. This study contributes to our understanding of the molecular characteristics and evolutionary relationships of this plant species.

Chloroplast Genome Sequences and Phylogenetic Analysis of Eight Newly Sequenced Caryophyllaceae Species

BACKGROUND

Caryophyllaceae is a big family composed of many economic and medicinal species. However, the phylogeny of the family is insufficient and genome data are lacking for many species.

OBJECTIVE

Using next-generation sequencing (NGS) to acquire the chloroplast (cp) genomes of Eremogone acicularis (F.N.Williams) Ikonn., E. brevipetala (Tsui & L.H.Zhou) Sadeghian & Zarre, E. bryophylla (Fernald) Pusalkar & D.K.Singh, E. kansuensis (Maxim.) Dillenb. & Kadereit, Shivparvatia glanduligera (Edgew.) Pusalkar & D.K.Singh, Silene atsaensis (Marq.) Bocquet, S. caespitella Williams, and S. lhassana (Williams) Majumdar.

METHODS

Bioinformatic software was used to conduct the comparative genome and phylogeny analysis of these cp genomes.

RESULTS

The eight cp genomes were 132 188-151 919 bp in length, containing 130-132 genes. A/T was dominant in simple sequence repeats (SSRs). Forward repeats and palindromic repeats were the most frequent in long terminal repeats (LTRs). Compared with the four species of Eremogone Fenzl, the inverted repeat (IR) boundaries of S. caespitella, S. atsaensis, S. lhassana, and Sh. glanduligera were significantly expanded. Four and one mutational hotspots were identified in the large single copy (LSC) region and small single copy (SSC) region, respectively. The ratio of nonsynonymous substitution to synonymous substitution (Ka/Ks ratio) showed these cp genomes may have undergone strong purifying selection. In the phylogenetic trees, both Silene L. and Eremogone were monophyletic groups. However, Sh. glanduligera was closely related to Amaranthus hypochondriacus.

CONCLUSION

These results have provided new evidence and useful information for species identification, evolution, and genetic research on the Caryophyllaceae.

HIGHLIGHTS

In this study, eight newly sequenced cp genomes of Caryophyllaceae species were reported for the first time.

Deciphering peanut complex genomes paves a way to understand its origin and domestication

Peanut (Arachis) is a key oil and protein crop worldwide with large genome. The genomes of diploid and tetraploid peanuts have been sequenced, which were compared to decipher their genome structures, evolutionary, and life secrets. Genome sequencing efforts showed that different cultivars, although Bt homeologs being more privileged in gene retention and gene expression. This subgenome bias, extended to sequence variation and point mutation, might be related to the long terminal repeat (LTR) explosions after tetraploidization, especially in At subgenomes. Except that, whole-genome sequences revealed many important genes, for example, fatty acids and triacylglycerols pathway, NBS-LRR (nucleotide-binding site-leucine-rich repeats), and seed size decision genes, were enriched after recursive polyploidization. Each ancestral polyploidy, with old ones having occurred hundreds of thousand years ago, has thousands of duplicated genes in extant genomes, contributing to genetic novelty. Notably, although full genome sequences are available, the actual At subgenome ancestor has still been elusive, highlighted with new debate about peanut origin. Although being an orphan crop lagging behind other crops in genomic resources, the genome sequencing achievement has laid a solid foundation for advancing crop enhancement and system biology research of peanut.

Molecular characterizations of genes in chloroplast genomes of the genus Arachis L. (Fabaceae) based on the codon usage divergence

Studies on the molecular characteristics of chloroplast genome are generally important for clarifying the evolutionary processes of plant species. The base composition, the effective number of codons, the relative synonymous codon usage, the codon bias index, and their correlation coefficients of a total of 41 genes in 21 chloroplast genomes of the genus Arachis were investigated to further perform the correspondence and clustering analyses, revealing significantly higher variations in genomes of wild species than those of the cultivated taxa. The codon usage patterns of all 41 genes in the genus Arachis were AT-rich, suggesting that the natural selection was the main factor affecting the evolutionary history of these genomes. Five genes (i.e., ndhC, petD, atpF, rpl14, and rps11) and five genes (i.e., atpE, psbD, psaB, ycf2, and rps12) showed higher and lower base usage divergences, respectively. This study provided novel insights into our understanding of the molecular evolution of chloroplast genomes in the genus Arachis.

Comparison of chloroplast genomes and phylogenomics in the Ficus sarmentosa complex (Moraceae)

Due to maternal inheritance and minimal rearrangement, the chloroplast genome is an important genetic resource for evolutionary studies. However, the evolutionary dynamics and phylogenetic performance of chloroplast genomes in closely related species are poorly characterized, particularly in taxonomically complex and species-rich groups. The taxonomically unresolved Ficus sarmentosa species complex (Moraceae) comprises approximately 20 taxa with unclear genetic background. In this study, we explored the evolutionary dynamics, hotspot loci, and phylogenetic performance of thirteen chloroplast genomes (including eleven newly obtained and two downloaded from NCBI) representing the F. sarmentosa complex. Their sequence lengths, IR boundaries, repeat sequences, and codon usage were compared. Both sequence length and IR boundaries were found to be highly conserved. All four categories of long repeat sequences were found across all 13 chloroplast genomes, with palindromic and forward sequences being the most common. The number of simple sequence repeat (SSR) loci varied from 175 (F. dinganensis and F. howii) to 190 (F. polynervis), with the dinucleotide motif appearing the most frequently. Relative synonymous codon usage (RSCU) analysis indicated that codons ending with A/T were prior to those ending with C/T. The majority of coding sequence regions were found to have undergone negative selection with the exception of ten genes (accD, clpP, ndhK, rbcL, rpl20, rpl22, rpl23, rpoC1, rps15, and rps4) which exhibited potential positive selective signatures. Five hypervariable genic regions (rps15, ycf1, rpoA, ndhF, and rpl22) and five hypervariable intergenic regions (trnH-GUG-psbA, rpl32-trnL-UAG, psbZ-trnG-GCC, trnK-UUU-rps16 and ndhF-rpl32) were identified. Overall, phylogenomic analysis based on 123 Ficus chloroplast genomes showed promise for studying the evolutionary relationships in Ficus, despite cyto-nuclear discordance. Furthermore, based on the phylogenetic performance of the F. sarmentosa complex and F. auriculata complex, the chloroplast genome also exhibited a promising phylogenetic resolution in closely related species.

Complete chloroplast genome sequence of Lens ervoides and comparison to Lens culinaris

Lens is a member of the Papilionoideae subfamily of Fabaceae and is generally used as a source of vegetable protein as part of human diets in many regions worldwide. Chloroplast (cp) genomes are highly active genetic components of plants and can be utilized as molecular markers for various purposes. As one of the wild lentil species, the Lens ervoides cp genome has been sequenced for the first time in this study using next-generation sequencing. The de novo assembly of the cp genome resulted in a single 122,722 bp sequence as two separate coexisting structural haplotypes with similar lengths. Results indicated that the cp genome of L. ervoides belongs to the inverted repeat lacking clade. Several noteworthy divergences within the coding regions were observed in ndhB, ndhF, rbcL, rpoC2, and ycf2 genes. Analysis of relative synonymous codon usage showed that certain genes, psbN, psaI, psbI, psbE, psbK, petD, and ndhC, preferred using biased codons more often and therefore might have elevated expression and translation efficiencies. Overall, this study exhibited the divergence level between the wild-type and cultured lentil cp genomes and pointed to certain regions that can be utilized as distinction markers for various goals.

Comparative repeatome analysis reveals new evidence on genome evolution in wild diploid Arachis (Fabaceae) species

Opposing changes in the abundance of satellite DNA and long terminal repeat (LTR) retroelements are the main contributors to the variation in genome size and heterochromatin amount in Arachis diploids. The South American genus Arachis (Fabaceae) comprises 83 species organized in nine taxonomic sections. Among them, section Arachis is characterized by species with a wide genome and karyotype diversity. Such diversity is determined mainly by the amount and composition of repetitive DNA. Here we performed computational analysis on low coverage genome sequencing to infer the dynamics of changes in major repeat families that led to the differentiation of genomes in diploid species (x = 10) of genus Arachis, focusing on section Arachis. Estimated repeat content ranged from 62.50 to 71.68% of the genomes. Species with different genome composition tended to have different landscapes of repeated sequences. Athila family retrotransposons were the most abundant and variable lineage among Arachis repeatomes, with peaks of transpositional activity inferred at different times in the evolution of the species. Satellite DNAs (satDNAs) were less abundant, but differentially represented among species. High rates of evolution of an AT-rich superfamily of satDNAs led to the differential accumulation of heterochromatin in Arachis genomes. The relationship between genome size variation and the repetitive content is complex. However, largest genomes presented a higher accumulation of LTR elements and lower contents of satDNAs. In contrast, species with lowest genome sizes tended to accumulate satDNAs in detriment of LTR elements. Phylogenetic analysis based on repetitive DNA supported the genome arrangement of section Arachis. Altogether, our results provide the most comprehensive picture on the repeatome dynamics that led to the genome differentiation of Arachis species.

Comparative Chloroplast Genomes of Nicotiana Species (Solanaceae): Insights Into the Genetic Variation, Phylogenetic Relationship, and Polyploid Speciation

Nicotiana L. is a genus rich in polyploidy, which represents an ideal natural system for investigating speciation, biodiversity, and phytogeography. Despite a wealth of phylogenetic work on this genus, a robust evolutionary framework with a dated molecular phylogeny for the genus is still lacking. In this study, the 19 complete chloroplast genomes of Nicotiana species were assembled, and five published chloroplast genomes of Nicotiana were retrieved for comparative analyses. The results showed that the 24 chloroplast genomes of Nicotiana, ranging from 155,327 bp (N. paniculata) to 156,142 bp (N. heterantha) in size, exhibited typical quadripartite structure. The chloroplast genomes were rather conserved in genome structure, GC content, RNA editing sites, and gene content and order. The higher GC content observed in the IR regions could be a result of the presence of abundant rRNA and tRNA genes, which contained a relatively higher GC content. A total of seven hypervariable regions, as new molecular markers for phylogenetic analysis, were uncovered. Based on 78 protein-coding genes, we constructed a well-supported phylogenetic tree, which was largely in agreement with previous studies, except for a slight conflict in several sections. Chloroplast phylogenetic results indicated that the progenitors of diploid N. sylvestris, N. knightiana, and the common ancestor of N. sylvestris and N. glauca might have donated the maternal genomes of allopolyploid N. tabacum, N. rustica, and section Repandae, respectively. Meanwhile, the diploid section Noctiflorae lineages (N. glauca) acted as the most likely maternal progenitor of section Suaveolentes. Molecular dating results show that the polyploid events range considerably in ~0.12 million (section Nicotiana) to ~5.77 million (section Repandae) years ago. The younger polyploids (N. tabacum and N. rustica) were estimated to have arisen ~0.120 and ~0.186 Mya, respectively. The older polyploids (section Repandae and Suaveolentes) were considered to have originated from a single polyploid event at ~5.77 and ~4.49 Mya, respectively. In summary, the comparative analysis of chloroplast genomes of Nicotiana species has not only revealed a series of new insights into the genetic variation and phylogenetic relationships in Nicotiana but also provided rich genetic resources for speciation and biodiversity research in the future.

Chloroplast Phylogenomic Analyses Reveal a Maternal Hybridization Event Leading to the Formation of Cultivated Peanuts

Peanuts (Arachis hypogaea L.) offer numerous healthy benefits, and the production of peanuts has a prominent role in global food security. As a result, it is in the interest of society to improve the productivity and quality of peanuts with transgenic means. However, the lack of a robust phylogeny of cultivated and wild peanut species has limited the utilization of genetic resources in peanut molecular breeding. In this study, a total of 33 complete peanut plastomes were sequenced, analyzed and used for phylogenetic analyses. Our results suggest that sect. Arachis can be subdivided into two lineages. All the cultivated species are contained in Lineage I with AABB and AA are the two predominant genome types present, while species in Lineage II possess diverse genome types, including BB, KK, GG, etc. Phylogenetic studies also indicate that all allotetraploid cultivated peanut species have been derived from a possible maternal hybridization event with one of the diploid Arachis duranensis accessions being a potential AA sub-genome ancestor. In addition, Arachis monticola, a tetraploid wild species, is placed in the same group with all the cultivated peanuts, and it may represent a transitional species, which has been through the recent hybridization event. This research could facilitate a better understanding of the taxonomic status of various Arachis species/accessions and the evolutionary relationship among them, and assists in the correct and efficient use of germplasm resources in breeding efforts to improve peanuts for the benefit of human beings.

Characterization of the complete chloroplast genome of Arachis pintoi Krapov. & W.C.Greg., a perennial leguminous forage

Arachis pintoi Krapov. & W.C.Greg. is an important leguminous forage grass species that have extremely wide ranges of distribution in the tropical and sub-tropical regions. It has high feeding value and horticultural value. In this study, we sequenced and assembled the complete chloroplast genome of A. pintoi. The chloroplast genome is 156,185 bp in length, containing a pair of inverted repeated (IR) regions of 25,820 bp that are separated by a large single copy (LSC) region of 85,637 bp, and a small single copy (SSC) region of 18,908 bp. The complete chloroplast genome contains 112 unique genes, including 80 protein-coding genes, 28 transfer RNA genes (tRNAs), and four ribosomal RNA genes (rRNAs). The overall GC content was 36.4%. The phylogenetic analysis demonstrated that A. pintoi formed a single branch among genus Arachis. The whole chloroplast genome of A. pintoi will be a useful resource for future studies on phylogeny and conservation in Arachis.

New Insight into the Phylogeny and Taxonomy of Cultivated and Related Species of Crataegus in China, Based on Complete Chloroplast Genome Sequencing

Hawthorns (Crataegus L.) are one of the most important processing and table fruits in China, due to their medicinal properties and health benefits. However, the interspecific relationships and evolution history of cultivated Crataegus in China remain unclear. Our previously published data showed C. bretschneideri may be derived from the hybridization of C. pinnatifida with C. maximowiczii, and that introgression occurs between C. hupehensis, C. pinnatifida, and C. pinnatifida var. major. In the present study, chloroplast sequences were used to further elucidate the phylogenetic relationships of cultivated Crataegus native to China. The chloroplast genomes of three cultivated species and one related species of Crataegus were sequenced for comparative and phylogenetic analyses. The four chloroplast genomes of Crataegus exhibited typical quadripartite structures and ranged from 159,607 bp (C. bretschneideri) to 159,875 bp (C. maximowiczii) in length. The plastomes of the four species contained 113 genes consisting of 79 protein-coding genes, 30 tRNA genes, and 4 rRNA genes. Six hypervariable regions (ndhC-trnV(UAC)-trnM(CAU), ndhA, atpH-atpI, ndhF, trnR(UCU)-atpA, and ndhF-rpl32), 196 repeats, and a total of 386 simple sequence repeats were detected as potential variability makers for species identification and population genetic studies. In the phylogenomic analyses, we also compared the entire chloroplast genomes of three published Crataegus species: C. hupehensis (MW201730.1), C. pinnatifida (MN102356.1), and C. marshallii (MK920293.1). Our phylogenetic analyses grouped the seven Crataegus taxa into two main clusters. One cluster included C. bretschneideri, C. maximowiczii, and C. marshallii, whereas the other included C. hupehensis, C. pinnatifida, and C. pinnatifida var. major. Taken together, our findings indicate that C. maximowiczii is the maternal origin of C. bretschneideri. This work provides further evidence of introgression between C. hupehensis, C. pinnatifida, and C. pinnatifida var. major, and suggests that C. pinnatifida var. major might have been artificially selected and domesticated from hybrid populations, rather than evolved from C. pinnatifida.

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.

Interspecific Hybrids Between Pelargonium × hortorum and Species From P. Section Ciconium Reveal Biparental Plastid Inheritance and Multi-Locus Cyto-Nuclear Incompatibility

The genetics underlying Cyto-Nuclear Incompatibility (CNI) was studied in Pelargonium interspecific hybrids. We created hybrids of 12 closely related crop wild relatives (CWR) with the ornamental P. × hortorum. Ten of the resulting 12 (F₁) interspecific hybrids segregate for chlorosis suggesting biparental plastid inheritance. The segregation ratios of the interspecific F₂ populations show nuclear interactions of one, two, or three nuclear genes regulating plastid function dependent on the parents. We further validated that biparental inheritance of plastids is common in section Ciconium, using diagnostic PCR primers. Our results pave the way for using the diverse species from section Ciconium, each with its own set of characteristics, as novel sources of desired breeding traits for P. × hortorum cultivars.