Phylogenetic analysis of 47 chloroplast genomes clarifies the contribution of wild species to the domesticated apple maternal line

Comparative chloroplast genome analyses of cultivated and wild Capsicum species shed light on evolution and phylogeny

The chloroplast (cp.) genome, also known as plastome, plays crucial roles in plant survival, adaptation, and evolution. The stable genetic structure of cp. genomes provides an ideal system for investigating species evolution. We sequenced three complete cp. genome sequences of Capsicum species and analyzed them using sequences of various Capsicum species retrieved from the NCBI database. The cp. genome of Capsicum species maintains a well-preserved quadripartite structure consisting of two inverted repeats (IRs) flanked by a large single copy (LSC) region and a small single copy (SSC) region. The sizes of cp. genome sequences ranged from 156,583 bp (C. lycianthoides) to 157,390 bp (C.pubescens). A total of 127-132 unique genes, including 83-87 protein-coding, 36-37 tRNA, and eight rRNA genes, were predicted. Comparison of cp. genomes of 10 Capsicum species revealed high sequence similarity in genome-wide organization and gene arrangements. Fragments of trnT-UGU/trnL-UAA, ccsA, ndhD, rps12, and ycf1 were identified as variable regions, and nucleotide variability of LSC and SSC was higher than that of IR. Phylogenetic speciation analysis showed that the major domesticated C. annuum species were the most extensively divergent species and closely related to C. tovarii and C. frutescens. Analysis of divergent times suggested that a substantial range of speciation events started occurring ~ 25.79 million years ago (Mya). Overall, comparative analysis of cp. genomes of Capsicum species not only offers new insights into their genetic variation and phylogenetic relationships, but also lays a foundation for evolutionary history, genetic diversity, conservation, and biological breeding of Capsicum species.

Pan-chloroplast genomes for accession-specific marker development in Hibiscus syriacus

Hibiscus syriacus L. is a renowned ornamental plant. We constructed 95 chloroplast genomes of H. syriacus L. cultivars using a short-read sequencing platform (Illumina) and a long-read sequencing platform (Oxford Nanopore Technology). The following genome assembly, we delineate quadripartite structures encompassing large single-copy, small single-copy, and inverted repeat (IRa and IRb) regions, from 160,231 bp to 161,041 bp. Our comprehensive analyses confirmed the presence of 79 protein-coding genes, 30 tRNA genes, and 4 rRNA genes in the pan-chloroplast genome, consistent with prior research on the H. syriacus chloroplast genome. Subsequent pangenome analysis unveiled widespread genome sequence conservation alongside unique cultivar-specific variant patterns consisting of 193 single-nucleotide polymorphisms and 61 insertions or deletions. The region containing intra-species variant patterns, as identified in this study, has the potential to develop accession-specific molecular markers, enhancing precision in cultivar classification. These findings are anticipated to drive advancements in breeding strategies, augment biodiversity, and unlock the agricultural potential inherent in H. syriacus.

A haplotype resolved chromosome-scale assembly of North American wild apple Malus fusca and comparative genomics of the fire blight Mfu10 locus

The Pacific crabapple (Malus fusca) is a wild relative of the commercial apple (Malus × domestica). With a range extending from Alaska to Northern California, M. fusca is extremely hardy and disease resistant. The species represents an untapped genetic resource for the development of new apple cultivars with enhanced stress resistance. However, gene discovery and utilization of M. fusca have been hampered by the lack of genomic resources. Here, we present a high‐quality, haplotype‐resolved, chromosome‐scale genome assembly and annotation for M. fusca. The genome was assembled using high‐fidelity long‐reads and scaffolded using genetic maps and high‐throughput chromatin conformation capture sequencing, resulting in one of the most contiguous apple genomes to date. We annotated the genome using public transcriptomic data from the same species taken from diverse plant structures and developmental stages. Using this assembly, we explored haplotypic structural variation within the genome of M. fusca, identifying thousands of large variants. We further showed high sequence co‐linearity with other domesticated and wild Malus species. Finally, we resolve a known quantitative trait locus associated with resistance to fire blight (Erwinia amylovora). Insights gained from the assembly of a reference‐quality genome of this hardy wild apple relative will be invaluable as a tool to facilitate DNA‐informed introgression breeding.

Comparative Analysis of Tylosema esculentum Mitochondrial DNA Revealed Two Distinct Genome Structures

Tylosema esculentum, commonly known as the marama bean, is an underutilized legume with nutritious seeds, holding potential to enhance food security in southern Africa due to its resilience to prolonged drought and heat. To promote the selection of this agronomically valuable germplasm, this study assembled and compared the mitogenomes of 84 marama individuals, identifying variations in genome structure, single-nucleotide polymorphisms (SNPs), insertions/deletions (indels), heteroplasmy, and horizontal transfer. Two distinct germplasms were identified, and a novel mitogenome structure consisting of three circular molecules and one long linear chromosome was discovered. The structural variation led to an increased copy number of specific genes, nad5, nad9, rrnS, rrn5, trnC, and trnfM. The two mitogenomes also exhibited differences at 230 loci, with only one notable nonsynonymous substitution in the matR gene. Heteroplasmy was concentrated at certain loci on chromosome LS1 (OK638188). Moreover, the marama mitogenome contained an over 9 kb insertion of cpDNA, originating from chloroplast genomes, but had accumulated mutations and lost gene functionality. The evolutionary and comparative genomics analysis indicated that mitogenome divergence in marama might not be solely constrained by geographical factors. Additionally, marama, as a member from the Cercidoideae subfamily, tends to possess a more complete set of mitochondrial genes than Faboideae legumes.

Molecular Structure and Variation Characteristics of the Plastomes from Six Malus baccata (L.) Borkh. Individuals and Comparative Genomic Analysis with Other Malus Species

Malus baccata (L.) Borkh. is an important wild species of Malus. Its rich variation types and population history are not well understood. Chloroplast genome mining plays an active role in germplasm identification and genetic evolution. In this study, by assembly and annotation, six complete cp genome sequences, ranging in size from 160,083 to 160,295 bp, were obtained. The GC content of stable IR regions (42.7%) was significantly higher than that of full length (36.5%) and SC regions (LSC-34.2%, SSC-30.4%). Compared with other Malus species, it was found that there were more sites of polymorphisms and hotspots of variation in LSC and SSC regions, with high variation sites including trnR/UCU-atpA, trnT/UGU-trnL/UAA, ndhF-rpl32 and ccsA-ndhD. The intraspecific and interspecific collinearity was good, and no structural rearrangement was observed. A large number of repeating elements and different boundary expansions may be involved in shaping the cp genome size. Up to 77 or 78 coding genes were annotated in the cp genomes of M. baccata, and high frequency codons such as UUA (Leu), GCU (Ala) and AGA (Arg) were identified by relative synonymous codon usage analysis. Phylogeographic analysis showed that 12 individuals of M. baccata clustered into three different groups with complex structure, whereas variant xiaojinensis (M.H. Cheng & N.G. Jiang) was not closely related to M. baccata evolutionarily. The phylogenetic analysis suggested that two main clades of different M. baccata in the genus Malus were formed and that I and II diverged about 9.7 MYA. In conclusion, through cp genome assembly and comparison, the interspecific relationships and molecular variations of M. baccata were further elucidated, and the results of this study provide valuable information for the phylogenetic evolution and germplasm conservation of M. baccata and Malus.

DNA profiling with the 20K apple SNP array reveals Malus domestica hybridization and admixture in M. sieversii, M. orientalis, and M. sylvestris genebank accessions

The USDA-ARS National Plant Germplasm System (NPGS) apple collection in Geneva, NY, USA maintains accessions of the primary Malus domestica (Suckow) Borkh. progenitor species M. sieversii (Ledeb.) M. Roem., M. orientalis Uglitzk., and M. sylvestris (L.) Mill. Many of these accessions originated from seeds that were collected from wild populations in the species' centers of diversity. Some of these accessions have fruit phenotypes that suggest recent M. domestica hybridization, which if true would represent crop contamination of wild species populations and mislabeled species status of NPGS accessions. Pedigree connections and admixture between M. domestica and its progenitor species can be readily identified with apple SNP array data, despite such arrays not being designed for these purposes. To investigate species purity, most (463 accessions) of the NPGS accessions labeled as these three progenitor species were genotyped using the 20K apple SNP array. DNA profiles obtained were compared with a dataset of more than 5000 unique M. domestica apple cultivars. Only 212 accessions (151 M. sieversii, 26 M. orientalis, and 35 M. sylvestris) were identified as "pure" species representatives because their DNA profiles did not exhibit genotypic signatures of recent hybridization with M. domestica. Twenty-one accessions (17 M. sieversii, 1 M. orientalis, and 3 M. sylvestris) previously labeled as wild species were instead fully M. domestica. Previously unrealized hybridization and admixture between wild species and M. domestica was identified in 230 accessions (215 M. sieversii, 9 M. orientalis, and 6 M. sylvestris). Among these species-mislabeled accessions, 'Alexander', 'Gold Reinette', 'Charlamoff', 'Rosmarina Bianca', and 'King of the Pippins' were the most frequently detected M. domestica parents or grandparents. These results have implications for collection management, including germplasm distribution, and might affect conclusions of previous research focused on these three progenitor species in the NPGS apple collection. Specifically, accessions received from the NPGS for breeding and genomics, genetics, and evolutionary biology research might not be truly representative of their previously assigned species.

Genetic structure and geneflow of Malus across the Korean Peninsula using genotyping-by-sequencing

This study was to understand the genetic structure and diversity of the Korean Malus species. We used genotyping-by-sequencing (GBS) technology to analyze samples of 112 individuals belonging to 18 populations of wild Malus spp. Using GBS, we identified thousands of single nucleotide polymorphisms in the species analyzed. M. baccata and M. toringo, two dominant mainland species of the Korean Peninsula, were distinguishable based on their genetic structure. However, M. toringo collected from Jeju Island exhibited a different genetic profile than that from the mainland. We identified M. cf. micromalus as a hybrid resulting from the Jeju Island M. toringo (pollen donor) and the mainland M. baccata, (pollen recipient). Putative M. mandshurica distributed on the Korean Peninsula showed a high structural and genetic similarity with M. baccata, indicating that it might be an ecotype. Overall, this study contributes to the understanding of the population history and genetic structure of Malus in the Korean Peninsula.

Rearrangement and domestication as drivers of Rosaceae mitogenome plasticity

Background

The mitochondrion is an important cellular component in plants and that functions in producing vital energy for the cell. However, the evolution and structure of mitochondrial genomes (mitogenomes) remain unclear in the Rosaceae family. In this study, we assembled 34 Rosaceae mitogenomes and characterized genome variation, rearrangement rate, and selection signal variation within these mitogenomes.

Results

Comparative analysis of six genera from the Amygdaloideae and five from the Rosoideae subfamilies of Rosaceae revealed that three protein-coding genes were absent from the mitogenomes of five Rosoideae genera. Positive correlations between genome size and repeat content were identified in 38 Rosaceae mitogenomes. Twenty repeats with high recombination frequency (> 50%) provided evidence for predominant substoichiometric conformation of the mitogenomes. Variations in rearrangement rates were identified between eleven genera, and within the Pyrus, Malus, Prunus, and Fragaria genera. Based on population data, phylogenetic inferences from Pyrus mitogenomes supported two distinct maternal lineages of Asian cultivated pears. A Pyrus-specific deletion (DEL-D) in selective sweeps was identified based on the assembled genomes and population data. After the DEL-D sequence fragments originally arose, they may have experienced a subsequent doubling event via homologous recombination and sequence transfer in the Amygdaloideae; afterwards, this variant sequence may have significantly expanded to cultivated groups, thereby improving adaptation during the domestication process.

Conclusions

This study characterizes the variations in gene content, genome size, rearrangement rate, and the impact of domestication in Rosaceae mitogenomes and provides insights into their structural variation patterns and phylogenetic relationships.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01383-3.

Phylogenetic Analysis of Wild Species and the Maternal Origin of Cultivars in the Genus Lilium Using 114 Plastid Genomes

Lilies are one of the most important ornamental flowers worldwide with approximately 100 wild species and numerous cultivars, but the phylogenetic relationships among wild species and their contributions to these cultivars are poorly resolved. We collected the major Lilium species and cultivars and assembled their plastome sequences. Our phylogenetic reconstruction using 114 plastid genomes, including 70 wild species representing all sections and 42 cultivars representing six hybrid divisions and two outgroups, uncovered well-supported genetic relationships within Lilium. The wild species were separated into two distinct groups (groups A and B) associated with geographical distribution, which further diversified into eight different clades that were phylogenetically well supported. Additional support was provided by the distributions of indels and single-nucleotide variants, which were consistent with the topology. The species of sections Archelirion, Sinomartagon III, and Leucolirion 6a and 6b were the maternal donors for Oriental hybrids, Asiatic hybrids, Trumpet hybrids, and Longiflorum hybrids, respectively. The maternal donors of the OT hybrids originated from the two sections Archelirion and Leucolirion 6a, and LA hybrids were derived from the two sections Leucolirion 6b and Sinomartagon. Our study provides an important basis for clarifying the infrageneric classification and the maternal origin of cultivars in Lilium.

Phylogenomic conflict analyses in the apple genus Malus s.l. reveal widespread hybridization and allopolyploidy driving diversification, with insights into the complex biogeographic history in the Northern Hemisphere

Phylogenomic evidence from an increasing number of studies has demonstrated that different data sets and analytical approaches often reconstruct strongly supported but conflicting relationships. In this study, 785 single-copy nuclear genes and 75 complete plastomes were used to infer the phylogenetic relationships and estimate the historical biogeography of the apple genus Malus sensu lato, an economically important lineage disjunctly distributed in the Northern Hemisphere and involved in known and suspected hybridization and allopolyploidy events. The nuclear phylogeny recovered the monophyly of Malus s.l. (including Docynia); however, the genus was supported to be biphyletic in the plastid phylogeny. An ancient chloroplast capture event in the Eocene in western North America best explains the cytonuclear discordance. Our conflict analysis demonstrated that ILS, hybridization, and allopolyploidy could explain the widespread nuclear gene tree discordance. One deep hybridization event (Malus doumeri) and one recent event (Malus coronaria) were detected in Malus s.l. Furthermore, our historical biogeographic analysis integrating living and fossil data supported a widespread East Asian-western North American origin of Malus s.l. in the Eocene, followed by several extinction and dispersal events in the Northern Hemisphere. We also propose a general workflow for assessing phylogenomic discordance and biogeographic analysis using deep genome skimming data sets.

Identification of the Complete Chloroplast Genome of Malus zhaojiaoensis Jiang and Its Comparison and Evolutionary Analysis with Other Malus Species

The genus Malus is rich in species and many of its plastid genomes have been released. However, limited resources and few markers are not conducive to the comparison of differences among species and resource identification and evaluation. In this study, the complete chloroplast genome of Malus zhaojiaoensis was studied by NGS sequencing, with a total length of 159998 bp. It consists of four regions, LSC (88,070 bp), IRB (26,359 bp), SSC (19,210 bp) and IRA (26,359 bp). M. zhaojiaoensis cp genome contained a total of 111 genes made up of three classes: 78 coding sequences, 29 tRNA genes, and four rRNA genes. In addition, a total of 91 SSRs and 43 INEs were found in the M. zhaojiaoensis cp genome, which was slightly different from M. baccata and M. hupehensis in number. The analysis of codon usage and RNA editing showed that high-frequency codons tended to end at A/U bases and RNA editing mainly occurred at the second codon. Comparative genome analysis suggested that the cp genomes of eight Malus species had higher overall similarity, but there were more variation hotspots (rps16_trnK-UUU, trnG-UCC_atpA, atpH_atpF, trnT-GGU_psbD, etc.) in the LSC region. By building evolutionary trees, it can be clearly observed that M. zhaojiaoensis formed a large group with eight species of Malus, but was relatively independent in differentiation. In conclusion, this study provides high-quality chloroplast genome resources of M. zhaojiaoensis and discusses the genetic variation characteristics of Malus genus. The findings of this study will provide a good reference for plastid genome assembly and interspecific comparison in the future.

Complete chloroplast genome studies of different apple varieties indicated the origin of modern cultivated apples from Malus sieversii and Malus sylvestris

Background

Apple is one of the most important temperate deciduous fruit trees worldwide, with a wide range of cultivation. In this study, we assessed the variations and phylogenetic relationships between the complete chloroplast genomes of wild and cultivated apples (Malus spp.).

Method

We obtained the complete chloroplast genomes of 24 apple varieties using next-generation sequencing technology and compared them with genomes of (downloaded from NCBI) the wild species.

Result

The chloroplast genome of Malus is highly conserved, with a genome length of 160,067-160,290 bp, and all have a double-stranded circular tetrad structure. The gene content and sequences of genomes of wild species and cultivated apple were almost the same, but several mutation hotspot regions (psbI-atpA, psbM-psbD, and ndhC-atpE) were detected in these genomes. These regions can provide valuable information for solving specific molecular markers in taxonomic research. Phylogenetic analysis revealed that Malus formed a new clade and four cultivated varieties clustered into a branch with M. sylvestris and M. sieversii, which indicated that M. sylvestris and M. sieversii were the ancestor species of the cultivated apple.

Chloroplast phylogenomics in Camelina (Brassicaceae) reveals multiple origins of polyploid species and the maternal lineage of C. sativa

The genus Camelina (Brassicaceae) comprises 7-8 diploid, tetraploid, and hexaploid species. Of particular agricultural interest is the biofuel crop, C. sativa (gold-of-pleasure or false flax), an allohexaploid domesticated from the widespread weed, C. microcarpa. Recent cytogenetics and genomics work has uncovered the identity of the parental diploid species involved in ancient polyploidization events in Camelina. However, little is known about the maternal subgenome ancestry of contemporary polyploid species. To determine the diploid maternal contributors of polyploid Camelina lineages, we sequenced and assembled 84 Camelina chloroplast genomes for phylogenetic analysis. Divergence time estimation was used to infer the timing of polyploidization events. Chromosome counts were also determined for 82 individuals to assess ploidy and cytotypic variation. Chloroplast genomes showed minimal divergence across the genus, with no observed gene-loss or structural variation. Phylogenetic analyses revealed C. hispida as a maternal diploid parent to the allotetraploid Camelina rumelica, and C. neglecta as the closest extant diploid contributor to the allohexaploids C. microcarpa and C. sativa. The tetraploid C. rumelica appears to have evolved through multiple independent hybridization events. Divergence times for polyploid lineages closely related to C. sativa were all inferred to be very recent, at only ~65 thousand years ago. Chromosome counts confirm that there are two distinct cytotypes within C. microcarpa (2n = 38 and 2n = 40). Based on these findings and other recent research, we propose a model of Camelina subgenome relationships representing our current understanding of the hybridization and polyploidization history of this recently-diverged genus.

Comparative Chloroplast Genome Analyses of the Winter-Blooming Eastern Asian Endemic Genus Chimonanthus (Calycanthaceae) With Implications For Its Phylogeny and Diversification

Chimonanthus of Calycanthaceae is a small endemic genus in China, with unusual winter-blooming sweet flowers widely cultivated for ornamentals and medicinal uses. The evolution of Chimonanthus plastomes and its phylogenetic relationships remain unresolved due to limited availability of genetic resources. Here, we report fully assembled and annotated chloroplast genomes of five Chimonanthus species. The chloroplast genomes of the genus (size range 153,010 – 153,299 bp) reveal high similarities in gene content, gene order, GC content, codon usage, amino acid frequency, simple sequence repeats, oligonucleotide repeats, synonymous and non-synonymous substitutions, and transition and transversion substitutions. Signatures of positive selection are detected in atpF and rpoB genes in C. campanulatus. The correlations among substitutions, InDels, and oligonucleotide repeats reveal weak to strong correlations in distantly related species at the intergeneric levels, and very weak to weak correlations among closely related Chimonanthus species. Chloroplast genomes are used to reconstruct a well-resolved phylogenetic tree, which supports the monophyly of Chimonanthus. Within Chimonanthus, C. praecox and C. campanulatus form one clade, while C. grammatus, C. salicifolius, C. zhejiangensis, and C. nitens constitute another clade. Chimonanthus nitens appears paraphyletic and is closely related to C. salicifolius and C. zhejiangensis, suggesting the need to reevaluate the species delimitation of C. nitens. Chimonanthus and Calycanthus diverged in mid-Oligocene; the radiation of extant Chimonanthus species was dated to the mid-Miocene, while C. grammatus diverged from other Chimonanthus species in the late Miocene. C. salicifolius, C. nitens(a), and C. zhejiangensis are inferred to have diverged in the Pleistocene of the Quaternary period, suggesting recent speciation of a relict lineage in the subtropical forest regions in eastern China. This study provides important insights into the chloroplast genome features and evolutionary history of Chimonanthus and family Calycanthaceae.

Phylogenetic relationships in Chinese oaks (Fagaceae, Quercus): Evidence from plastid genome using low-coverage whole genome sequencing

China is a second center of oak diversity but with less intensively systematic studies. Here, with 49 species representing all four sections in China, we firstly gave insight into the comprehensive phylogenetic relationships of Chinese oaks based on 54 complete plastid genomes. Our results recovered a robust phylogenetic framework and provided strong support for most nodes. The phylogenetic tree supported Quercus section Ilex as not monophyletic, in which Quercus section Cyclobalanopsis and Quercus section Cerris were nested. Most likely, incomplete lineage sorting and/or introgression among ancestral lineages in these three sections resulted in this complex pattern. The current distribution, diversification and molecular differentiation of Q. sect. Ilex in China are likely consequences of local adaptation to the geographic and paleoclimatic changes, which were driven by the uplift of Tibetan Plateau, the Hengduan Mountains and the Himalayas.

Transcriptional regulation of bark freezing tolerance in apple (Malus domestica Borkh.)

Freezing tolerance is a significant trait in plants that grow in cold environments and survive through the winter. Apple (Malus domestica Borkh.) is a cold-tolerant fruit tree, and the cold tolerance of its bark is important for its survival at low temperatures. However, little is known about the gene activity related to its freezing tolerance. To better understand the gene expression and regulation properties of freezing tolerance in dormant apple trees, we analyzed the transcriptomic divergences in the bark from 1-year-old branches of two apple cultivars, "Golden Delicious" (G) and "Jinhong" (H), which have different levels of cold resistance, under chilling and freezing treatments. "H" can safely overwinter below -30 °C in extremely low-temperature regions, whereas "G" experiences severe freezing damage and death in similar environments. Based on 28 bark transcriptomes (from the epidermis, phloem, and cambium) from 1-year-old branches under seven temperature treatments (from 4 to -29 °C), we identified 4173 and 7734 differentially expressed genes (DEGs) in "G" and "H", respectively, between the chilling and freezing treatments. A gene coexpression network was constructed according to this expression information using weighted gene correlation network analysis (WGCNA), and seven biologically meaningful coexpression modules were identified from the network. The expression profiles of the genes from these modules suggested the gene regulatory pathways that are responsible for the chilling and freezing stress responses of "G" and/or "H." Module 7 was probably related to freezing acclimation and freezing damage in "H" at the lower temperatures. This module contained more interconnected hub transcription factors (TFs) and cold-responsive genes (CORs). Modules 6 and 7 contained C-repeat binding factor (CBF) TFs, and many CBF-dependent homologs were identified as hub genes. We also found that some hub TFs had higher intramodular connectivity (KME) and gene significance (GS) than CBFs. Specifically, most hub TFs in modules 6 and 7 were activated at the beginning of the early freezing stress phase and maintained upregulated expression during the whole freezing stress period in "G" and "H". The upregulation of DEGs related to methionine and carbohydrate biosynthetic processes in "H" under more severe freezing stress supported the maintenance of homeostasis in the cellular membrane. This study improves our understanding of the transcriptional regulation patterns underlying freezing tolerance in the bark of apple branches.

The complete chloroplast genome sequences of five pinnate-leaved Primula species and phylogenetic analyses

The six pinnate-leaved species are a very particular group in the genus Primula. In the present paper, we sequenced, assembled and annotated the chloroplast genomes of five of them (P. cicutarrifolia, P. hubeiensis, P. jiugongshanensis, P. merrilliana, P. ranunculoides). The five chloroplast genomes ranged from ~ 150 to 152 kb, containing 113 genes (four ribosomal RNA genes, 29 tRNA genes and 80 protein-coding genes). The six pinnate-leaved species exhibited synteny of gene order and possessed similar IR boundary regions in chloroplast genomes. The gene accD was pseudogenized in P. filchnerae. In the chloroplast genomes of the six pinnate-leaved Primula species, SSRs, repeating sequences and divergence hotspots were identified; ycf1 and trnH-psbA were the most variable markers among CDSs and noncoding sequences, respectively. Phylogenetic analyses showed that the six Primula species were separated into two distant clades: one was formed by P. filchnerae and P. sinensis and the other clade was consisting of two subclades, one formed by P. hubeiensis and P. ranunculoides, the other by P. merrilliana, P. cicutarrifolia and P. jiugongshanensis. P. hubeiensis was closely related with P. ranunculoides and therefore it should be placed into Sect. Ranunculoides. P. cicutarrifolia did not group first with P. ranunculoides but with P. merrilliana, although the former two were once united in one species, our results supported the separation of P. ranunculoides from P. cicutarrifolia as one distinct species.

Plastome Characterization and Phylogenomics of East Asian Beeches with a Special Emphasis on Fagus multinervis on Ulleung Island, Korea

Beech trees of the genus Fagus (Fagaceae) are monoecious and distributed in the Northern Hemisphere. They represent an important component of mixed broad-leaved evergreen-deciduous forests and are an economically important source of timber. Despite their ecological and economical importance, however, little is known regarding the overall plastome evolution among Fagus species in East Asia. In particular, the taxonomic position and status of F. multinervis, a beech species endemic to Ulleung Island of Korea, remains unclear even today. Therefore, in this study, we characterized four newly completed plastomes of East Asian Fagus species (one accession each of F. crenata and F. multinervis and two accessions of F. japonica). Moreover, we performed phylogenomic analyses comparing these four plastomes with F. sylvatica (European beech) plastome. The four plastomes were highly conserved, and their size ranged from 158,163 to 158,348 base pair (bp). The overall GC content was 37.1%, and the sequence similarity ranged from 99.8% to 99.99%. Codon usage patterns were similar among species, and 7 of 77 common protein-coding genes were under positive selection. Furthermore, we identified five highly variable hotspot regions of the Fagus plastomes (ccsA/ndhD, ndhD/psaC, ndhF/rpl32, trnS-GCU/trnG-UCC, and ycf1). Phylogenetic analysis revealed the monophyly of Fagus as well as early divergence of the subgenus Fagus and monophyletic Engleriana. Finally, phylogenetic results supported the taxonomic distinction of F. multinervis from its close relatives F. engleriana and F. japonica. However, the sister species and geographic origin of F. multinervis on Ulleung Island could not be determined.

Threat to Asian wild apple trees posed by gene flow from domesticated apple trees and their "pestified" pathogens

Secondary contact between crops and their wild relatives poses a threat to wild species, not only through gene flow between plants, but also through the dispersal of crop pathogens and genetic exchanges involving these pathogens, particularly those that have become more virulent by indirect selection on resistant crops, a phenomenon known as "pestification." Joint analyses of wild and domesticated hosts and their pathogens are essential to address this issue, but such analyses remain rare. We used population genetics approaches, demographic inference and pathogenicity tests on host-pathogen pairs of wild or domesticated apple trees from Central Asia and their main fungal pathogen, Venturia inaequalis, which itself has differentiated agricultural and wild-type populations. We confirmed the occurrence of gene flow from cultivated (Malus domestica) to wild (Malus sieversii) apple trees in Asian forests, potentially threatening the persistence of Asian wild apple trees. Pathogenicity tests demonstrated the pestification of V. inaequalis, the agricultural-type population being more virulent on both wild and domesticated trees. Single nucleotide polymorphism (SNP) markers and the demographic modelling of pathogen populations revealed hybridization following secondary contact between agricultural and wild-type fungal populations, and dispersal of the agricultural-type pathogen population in wild forests, increasing the threat of disease in the wild apple species. We detected an SNP potentially involved in pathogen pestification, generating an early stop codon in a gene encoding a small secreted protein in the agricultural-type fungal population. Our findings, based on joint analyses of paired host and pathogen data sets, highlight the threat posed by cultivating a crop near its centre of origin, in terms of pestified pathogen invasions in wild plant populations and introgression in the wild-type pathogen population.

Genetic integrity is still maintained in natural populations of the indigenous wild apple species Malus sylvestris (Mill.) in Saxony as demonstrated with nuclear SSR and chloroplast DNA markers

Malus sylvestris (Mill.) is the only indigenous wild apple species in Central Europe. Agriculture, forestry, and urbanization increasingly endanger Malus sylvestris natural habitats. In addition, the risks of cross-hybridization associated with increase in the cultivation of the domesticated apple Malus × domestica (Borkh.) threaten the genetic integrity of M. sylvestris. The present study investigated the number of hybrids, genetic diversity, and genetic structure of 292 putative M. sylvestris that originate from five different natural M. sylvestris populations in Saxony, Germany. All samples were genetically analyzed using nine nuclear microsatellite markers (ncSSR) and four maternally inherited chloroplast markers (cpDNA) along with 56 apple cultivars commonly cultivated in Saxony. Eighty-seven percent of the wild apple accessions were identified as pure M. sylvestris. The cpDNA analysis showed six private haplotypes for M. sylvestris, whereas three haplotypes were present in M. sylvestris and M. × domestica. The analysis of molecular variance (AMOVA) resulted in a moderate (ncSSR) and great (cpDNA) variation among pure M. sylvestris and M. × domestica individuals indicating a low gene flow between both species. The genetic diversity within the pure M. sylvestris populations was high with a weak genetic structure between the M. sylvestris populations indicating an unrestricted genetic exchange between these M. sylvestris populations. The clear distinguishing of M. sylvestris and M. ×domestica confirms our expectation of the existence of pure M. sylvestris accessions in this area and supports the argument for the implementation of preservation measures to protect the M. sylvestris populations in Saxony.

Comparison of the cytoplastic genomes by resequencing: insights into the genetic diversity and the phylogeny of the agriculturally important genus Brassica

Background

The genus Brassica mainly comprises three diploid and three recently derived allotetraploid species, most of which are highly important vegetable, oil or ornamental crops cultivated worldwide. Despite being extensively studied, the origination of B. napus and certain detailed interspecific relationships within Brassica genus remains undetermined and somewhere confused. In the current high-throughput sequencing era, a systemic comparative genomic study based on a large population is necessary and would be crucial to resolve these questions.

Results

The chloroplast DNA and mitochondrial DNA were synchronously resequenced in a selected set of Brassica materials, which contain 72 accessions and maximally integrated the known Brassica species. The Brassica genomewide cpDNA and mtDNA variations have been identified. Detailed phylogenetic relationships inside and around Brassica genus have been delineated by the cpDNA- and mtDNA- variation derived phylogenies. Different from B. juncea and B. carinata, the natural B. napus contains three major cytoplasmic haplotypes: the cam-type which directly inherited from B. rapa, polima-type which is close to cam-type as a sister, and the mysterious but predominant nap-type. Certain sparse C-genome wild species might have primarily contributed the nap-type cytoplasm and the corresponding C subgenome to B. napus, implied by their con-clustering in both phylogenies. The strictly concurrent inheritance of mtDNA and cpDNA were dramatically disturbed in the B. napus cytoplasmic male sterile lines (e.g., mori and nsa). The genera Raphanus, Sinapis, Eruca, Moricandia show a strong parallel evolutional relationships with Brassica.

Conclusions

The overall variation data and elaborated phylogenetic relationships provide further insights into genetic understanding of Brassica, which can substantially facilitate the development of novel Brassica germplasms.

Characterization and Comparison of Two Complete Plastomes of Rosaceae Species (Potentilla dickinsii var. glabrata and Spiraea insularis) Endemic to Ulleung Island, Korea

Potentilla dickinsii var. glabrata and Spiraea insularis in the family Rosaceae are species endemic to Ulleung Island, Korea, the latter of which is listed as endangered. In this study, we characterized the complete plastomes of these two species and compared these with previously reported plastomes of other Ulleung Island endemic species of Rosaceae (Cotoneaster wilsonii, Prunus takesimensis, Rubus takesimensis, and Sorbus ulleungensis). The highly conserved complete plastomes of P. dickinsii var. glabrata and S. insularis are 158,637 and 155,524 base pairs with GC contents of 37% and 36.9%, respectively. Comparative phylogenomic analysis identified three highly variable intergenic regions (trnT-UGU/trnL-UAA, rpl32/trnL-UAG, and ndhF/rpl32) and one variable genic region (ycf1). Only 14 of the 75 protein-coding genes have been subject to strong purifying selection. Phylogenetic analysis of 23 representative plastomes within the Rosaceae supported the monophyly of Potentilla and the sister relationship between Potentilla and Fragaria and indicated that S. insularis is sister to a clade containing Cotoneaster, Malus, Pyrus, and Sorbus. The plastome resources generated in this study will contribute to elucidating the plastome evolution of insular endemic Rosaceae on Ulleung Island and also in assessing the genetic consequences of anagenetic speciation for various endemic lineages on the island.

Overexpression of chalcone isomerase in apple reduces phloridzin accumulation and increases susceptibility to herbivory by two-spotted mites

Apples (Malus spp.) accumulate significant quantities of the dihydrochalcone glycoside, phloridzin, whilst pears (Pyrus spp.) do not. To explain this difference, we hypothesized that a metabolic bottleneck in the phenylpropanoid pathway might exist in apple. Expression analysis indicated that transcript levels of early phenylpropanoid pathway genes in apple and pear leaves were similar, except for chalcone isomerase (CHI), which was much lower in apple. Apples also showed very low CHI activity compared with pear. To relieve the bottleneck at CHI, transgenic apple plants overexpressing the Arabidopsis AtCHI gene were produced. Unlike other transgenic apples where phenylpropanoid flux was manipulated, AtCHI overexpression (CHIox) plants were phenotypically indistinguishable from wild-type, except for an increase in red pigmentation in expanding leaves. CHIox plants accumulated slightly increased levels of flavanols and flavan-3-ols in the leaves, but the major change was a 2.8- to 19-fold drop in phloridzin concentrations compared with wild-type. The impact of these phytochemical changes on insect preference was studied using a two-choice leaf assay with the polyphagous apple pest, the two-spotted spider mite (Tetranychus urticae Koch). Transgenic CHIox leaves were more susceptible to herbivory, an effect that could be reversed (complemented) by application of phloridzin to transgenic leaves. Taken together, these findings shed new light on phenylpropanoid biosynthesis in apple and suggest a new physiological role for phloridzin as an antifeedant in leaves.

A Novel Test for Absolute Fit of Evolutionary Models Provides a Means to Correctly Identify the Substitution Model and the Model Tree

A novel test is described that visualizes the absolute model-data fit of the substitution and tree components of an evolutionary model. The test utilizes statistics based on counts of character state matches and mismatches in alignments of observed and simulated sequences. This comparison is used to assess model-data fit. In simulations conducted to evaluate the performance of the test, the test estimator was able to identify both the correct tree topology and substitution model under conditions where the Goldman-Cox test-which tests the fit of a substitution model to sequence data and is also based on comparing simulated replicates with observed data-showed high error rates. The novel test was found to identify the correct tree topology within a wide range of DNA substitution model misspecifications, indicating the high discriminatory power of the test. Use of this test provides a practical approach for assessing absolute model-data fit when testing phylogenetic hypotheses.

Intra-individual heteroplasmy in the Gentiana tongolensis plastid genome (Gentianaceae)

Chloroplasts are typically inherited from the female parent and are haploid in most angiosperms, but rare intra-individual heteroplasmy in plastid genomes has been reported in plants. Here, we report an example of plastome heteroplasmy and its characteristics in Gentiana tongolensis (Gentianaceae). The plastid genome of G. tongolensis is 145,757 bp in size and is missing parts of petD gene when compared with other Gentiana species. A total of 112 single nucleotide polymorphisms (SNPs) and 31 indels with frequencies of more than 2% were detected in the plastid genome, and most were located in protein coding regions. Most sites with SNP frequencies of more than 10% were located in six genes in the LSC region. After verification via cloning and Sanger sequencing at three loci, heteroplasmy was identified in different individuals. The cause of heteroplasmy at the nucleotide level in plastome of G. tongolensis is unclear from the present data, although biparental plastid inheritance and transfer of plastid DNA seem to be most likely. This study implies that botanists should reconsider the heredity and evolution of chloroplasts and be cautious with using chloroplasts as genetic markers, especially in Gentiana.

Plastome of Quercus xanthoclada and comparison of genomic diversity amongst selected Quercus species using genome skimming

The genus Quercus L. contains several of the most economically important species for timber production in the Northern Hemisphere. It was one of the first genera described, but genetic diversity at a global scale within and amongst oak species remains unclear, despite numerous regional or species-specific assessments. To evaluate global plastid diversity in oaks, we sequenced the complete chloroplast of Quercus xanthoclada and compared its sequence with those available from other main taxonomic groups in Quercus. We quantify genomic divergence amongst oaks and performed a sliding window analysis to detect the most variable regions amongst members of the various clades, as well as divergent regions occurring in specific pairs of species. We identified private and shared SNPs amongst oaks species and sections and stress the need for a large global assessment of genetic diversity in this economically and ecologically important genus.

Characterization of complete chloroplast genome of Malus sylvestris L

The European wild apple (Malus sylvestris L.) is an important economical fruit crop. In this present study, we characterized the complete chloroplast (cp) genome sequence of Malus sylvestris L. The complete cp genome is 159,926 bp in length with a typical quadripartite structure, containing a large single-copy region (88,064 bp), a small single-copy region (26,353 bp) and a pair of inverted repeat regions (19,157 bp each). A total of 110 unique genes were found in the newly sequenced genome, including 78 protein-coding genes, 28 tRNA genes, and 4 rRNA genes. Of these, 6 protein-coding genes, 7 tRNA genes, and all 4 rRNA genes are duplicated in the inverted regions. A phylogenetic tree was reconstructed using the neighbor-joining method based on the full length of cp genomes within genus Malus. The result showed that M. sylvestris L. was clustered together with the cultivated apple. The complete cp genome could provide valuable information for understanding the phylogenetic relationships within the genus Malus.

Structure and features of the complete chloroplast genome of Melastoma dodecandrum

Melastoma dodecandrum, the only creeping species in the Melastoma genus, serves as a medicinal herb in southeast China. It belongs to the huge family Melastomataceae, which contains over 5000 species worldwide. In this study, we used next-generation sequencing to determine the complete chloroplast genome sequences of M. dodecandrum, which is a circular molecule of 156,611 bp in length. After annotation, we identified 131 putative genes in total, comprised of 85 protein-coding genes, 38 transfer RNA genes and 8 ribosomal RNA genes. Genome structure, GC content, repeat sequences and codon usage were investigated to gain a comprehensive understanding of this genome. Furthermore, we conducted comparative genome analyses between the M. dodecandrum genome and that of four other Melastomataceae species. Additionally, a phylogenetic analysis was performed based on available chloroplast genomes of Melastomataceae species and several Myrtaceae species, revealing the taxonomic relationships between M. dodecandrum and related species. In conclusion, our study represents the first look into the complete chloroplast genome of M. dodecandrum, providing abundant information for further studies such as species identification, taxonomy and phylogenetic resolution of Melastomataceae species.

Wild Origins of Macadamia Domestication Identified Through Intraspecific Chloroplast Genome Sequencing

Identifying the geographic origins of crops is important for the conservation and utilization of novel genetic variation. Even so, the origins of many food crops remain elusive. The tree nut crop macadamia has a remarkable domestication history, from subtropical rain forests in Australia through Hawaii to global cultivation all within the last century. The industry is based primarily on Macadamia integrifolia and M. integrifolia-M. tetraphylla hybrid cultivars with Hawaiian cultivars the main contributors to world production. Sequence data from the chloroplast genome assembled using a genome skimming strategy was used to determine population structure among remnant populations of the main progenitor species, M. integrifolia. Phylogenetic analysis of a 506 bp chloroplast SNP alignment from 64 wild and cultivated accessions identified phylogeographic structure and deep divergences between clades providing evidence for historical barriers to seed dispersal. High levels of variation were detected among wild accessions. Most Hawaiian cultivars, however, shared a single chlorotype that was also present at two wild sites at Mooloo and Mt Bauple from the northernmost distribution of the species in south-east Queensland. Our results provide evidence for a maternal genetic bottleneck during early macadamia domestication, and pinpoint the likely source of seed used to develop the Hawaiian cultivars. The extensive variability and structuring of M. integrifolia chloroplast genomic variation detected in this study suggests much unexploited genetic diversity is available for improvement of this recently domesticated crop.

Genetic analysis of a major international collection of cultivated apple varieties reveals previously unknown historic heteroploid and inbred relationships

Domesticated apple (Malus x domestica Borkh.) is a major global crop and the genetic diversity held within the pool of cultivated varieties is important for the development of future cultivars. The aim of this study was to investigate the diversity held within the domesticated form, through the analysis of a major international germplasm collection of cultivated varieties, the UK National Fruit Collection, consisting of over 2,000 selections of named cultivars and seedling varieties. We utilised Diversity Array Technology (DArT) markers to assess the genetic diversity within the collection. Clustering attempts, using the software STRUCTURE revealed that the accessions formed a complex and historically admixed group for which clear clustering was challenging. Comparison of accessions using the Jaccard similarity coefficient allowed us to identify clonal and duplicate material as well as revealing pairs and groups that appeared more closely related than a standard parent-offspring or full-sibling relations. From further investigation, we were able to propose a number of new pedigrees, which revealed that some historically important cultivars were more closely related than previously documented and that some of them were partially inbred. We were also able to elucidate a number of parent-offspring relationships that had resulted in a number of important polyploid cultivars. This included reuniting polyploid cultivars that in some cases dated as far back as the 18th century, with diploid parents that potentially date back as far as the 13th century.

Phylogeny of Maleae (Rosaceae) Based on Multiple Chloroplast Regions: Implications to Genera Circumscription

Maleae consists of economically and ecologically important plants. However, there are considerable disputes on generic circumscription due to the lack of a reliable phylogeny at generic level. In this study, molecular phylogeny of 35 generally accepted genera in Maleae is established using 15 chloroplast regions. Gillenia is the most basal clade of Maleae, followed by Kageneckia + Lindleya, Vauquelinia, and a typical radiation clade, the core Maleae, suggesting that the proposal of four subtribes is reasonable. In the core Maleae including 31 genera, chloroplast gene data support that the four Malus-related genera should better be merged into one genus and the six Sorbus-related genera would be classified into two genera, whereas all Photinia-related genera should be accepted as distinct genera. Although the phylogenetic relationships among the genera in Maleae are much clearer than before, it is still premature to make a formal taxonomic treatment for these genera.

Re-exploration of U's Triangle Brassica Species Based on Chloroplast Genomes and 45S nrDNA Sequences

The concept of U's triangle, which revealed the importance of polyploidization in plant genome evolution, described natural allopolyploidization events in Brassica using three diploids [B. rapa (A genome), B. nigra (B), and B. oleracea (C)] and derived allotetraploids [B. juncea (AB genome), B. napus (AC), and B. carinata (BC)]. However, comprehensive understanding of Brassica genome evolution has not been fully achieved. Here, we performed low-coverage (2-6×) whole-genome sequencing of 28 accessions of Brassica as well as of Raphanus sativus [R genome] to explore the evolution of six Brassica species based on chloroplast genome and ribosomal DNA variations. Our phylogenomic analyses led to two main conclusions. (1) Intra-species-level chloroplast genome variations are low in the three allotetraploids (2~7 SNPs), but rich and variable in each diploid species (7~193 SNPs). (2) Three allotetraploids maintain two 45SnrDNA types derived from both ancestral species with maternal dominance. Furthermore, this study sheds light on the maternal origin of the AC chloroplast genome. Overall, this study clarifies the genetic relationships of U's triangle species based on a comprehensive genomics approach and provides important genomic resources for correlative and evolutionary studies.

Application of Chloroplast Phylogenomics to Resolve Species Relationships Within the Plant Genus Amaranthus

Amaranthus species are an emerging and promising nutritious traditional vegetable food source. Morphological plasticity and poorly resolved dendrograms have led to the need for well resolved species phylogenies. We hypothesized that whole chloroplast phylogenomics would result in more reliable differentiation between closely related amaranth species. The aims of the study were therefore: to construct a fully assembled, annotated chloroplast genome sequence of Amaranthus tricolor; to characterize Amaranthus accessions phylogenetically by comparing barcoding genes (matK, rbcL, ITS) with whole chloroplast sequencing; and to use whole chloroplast phylogenomics to resolve deeper phylogenetic relationships. We generated a complete A. tricolor chloroplast sequence of 150,027 bp. The three barcoding genes revealed poor inter- and intra-species resolution with low bootstrap support. Whole chloroplast phylogenomics of 59 Amaranthus accessions increased the number of parsimoniously informative sites from 92 to 481 compared to the barcoding genes, allowing improved separation of amaranth species. Our results support previous findings that two geographically independent domestication events of Amaranthus hybridus likely gave rise to several species within the Hybridus complex, namely Amaranthus dubius, Amaranthus quitensis, Amaranthus caudatus, Amaranthus cruentus and Amaranthus hypochondriacus. Poor resolution of species within the Hybridus complex supports the recent and ongoing domestication within the complex, and highlights the limitation of chloroplast data for resolving recent evolution. The weedy Amaranthus retroflexus and Amaranthus powellii was found to share a common ancestor with the Hybridus complex. Leafy amaranth, Amaranthus tricolor, Amaranthus blitum, Amaranthus viridis and Amaranthus graecizans formed a stable sister lineage to the aforementioned species across the phylogenetic trees. This study demonstrates the power of next-generation sequencing data and reference-based assemblies to resolve phylogenies, and also facilitated the identification of unknown Amaranthus accessions from a local genebank. The informative phylogeny of the Amaranthus genus will aid in selecting accessions for breeding advanced genotypes to satisfy global food demand.

Plastid Genome Comparative and Phylogenetic Analyses of the Key Genera in Fagaceae: Highlighting the Effect of Codon Composition Bias in Phylogenetic Inference

Fagaceae is one of the largest and economically important taxa within Fagales. Considering the incongruence among inferences from plastid and nuclear genes in the previous Fagaceae phylogeny studies, we assess the performance of plastid phylogenomics in this complex family. We sequenced and assembled four complete plastid genomes (Fagus engleriana, Quercus spinosa, Quercus aquifolioides, and Quercus glauca) using reference-guided assembly approach. All of the other 12 published plastid genomes in Fagaceae were retrieved for genomic analyses (including repeats, sequence divergence and codon usage) and phylogenetic inference. The genomic analyses reveal that plastid genomes in Fagaceae are conserved. Comparing the phylogenetic relationships of the key genera in Fagaceae inferred from different codon positions and gene function datasets, we found that the first two codon sites dataset recovered nearly all relationships and received high support. Thus, the result suggested that codon composition bias had great influence on Fagaceae phylogenetic inference. Our study not only provides basic understanding of Fagaceae plastid genomes, but also illuminates the effectiveness of plastid phylogenomics in resolving relationships of this intractable family.

The Complete Plastome Sequences of Seven Species in Gentiana sect. Kudoa (Gentianaceae): Insights Into Plastid Gene Loss and Molecular Evolution

The chloroplast (cp) genome is useful in the study of phylogenomics, molecular dating, and molecular evolution. Gentiana sect. Kudoa is a predominantly alpine flowering plant that is valued for its contributions to medicine, ecology, and horticulture. Previous evolutionary studies showed that the plastid gene loss pattern and intra-sectional phylogenetics in sect. Kudoa are still unclear. In this study, we compared 11 Gentiana plastomes, including 7 newly sequenced plastomes from sect. Kudoa, to represent its three serious: ser. Ornatae, ser. Verticillatae, and ser. Monanthae. The cp genome sizes of the seven species ranged from 137,278 to 147,156 bp. The plastome size variation mainly occurred in the small single-copy and long single-copy regions rather than the inverted repeat regions. Compared with sect. Cruciata, the plastomes in ser. Ornatae and ser. Verticillatae had lost approximately 11 kb of sequences containing 11 ndh genes. Conversely, far fewer losses were observed in ser. Monanthae. The phylogenetic tree revealed that sect. Kudoa was not monophyletic and that ser. Monanthae was more closely related to other sections rather than sect. Kudoa. The molecular dating analysis indicated that ser. Monanthae and sect. Kudoa diverged around 8.23 Ma. In ser. Ornatae and ser. Verticillatae, the divergence occurred at around 0.07-1.78 Ma. The nucleotide diversity analysis indicated that the intergenic regions trnH-psbA, trnK-trnQ, ycf3-trnS and rpl32-trnL constituted divergence hotspots in both sect. Kudoa and Gentiana, and would be useful for future phylogenetic and population genetic studies.

Combined Analyses of Chloroplast DNA Haplotypes and Microsatellite Markers Reveal New Insights Into the Origin and Dissemination Route of Cultivated Pears Native to East Asia

Asian pear plays an important role in the world pear industry, accounting for over 70% of world total production volume. Commercial Asian pear production relies on four major pear cultivar groups, Japanese pear (JP), Chinese white pear (CWP), Chinese sand pear (CSP), and Ussurian pear (UP), but their origins remain controversial. We estimated the genetic diversity levels and structures in a large sample of existing local cultivars to investigate the origins of Asian pears using twenty-five genome-covering nuclear microsatellite (simple sequence repeats, nSSR) markers and two non-coding chloroplast DNA (cpDNA) regions (trnL-trnF and accD-psaI). High levels of genetic diversity were detected for both nSSRs (H_E = 0.744) and cpDNAs (Hd = 0.792). The major variation was found within geographic populations of cultivated pear groups, demonstrating a close relationship among cultivar groups. CSPs showed a greater genetic diversity than CWPs and JPs, and lowest levels of genetic differentiation were detected among them. Phylogeographical analyses indicated that the CSP, CWP, and JP were derived from the same progenitor of Pyrus pyrifolia in China. A dissemination route of cultivated P. pyrifolia estimated by approximate Bayesian computation suggested that cultivated P. pyrifolia from the Middle Yangtze River Valley area contributed the major genetic resources to the cultivars, excluding those of southwestern China. Three major genetic groups of cultivated Pyrus pyrifolia were revealed using nSSRs and a Bayesian statistical inference: (a) JPs; (b) cultivars from South-Central China northward to northeastern China, covering the main pear production area in China; (c) cultivars from southwestern China to southeastern China, including Yunnan, Guizhou, Guangdong, Guangxi, and Fujian Provinces. This reflected the synergistic effects of ecogeographical factors and human selection during cultivar spread and improvement. The analyses indicated that UP cultivars might be originated from the interspecific hybridization of wild Pyrus ussuriensis with cultivated Pyrus pyrifolia. The combination of uniparental DNA sequences and nuclear markers give us a better understanding of origins and genetic relationships for Asian pear groups and will be beneficial for the future improvement of Asian pear cultivars.

Phylogenic study of Lemnoideae (duckweeds) through complete chloroplast genomes for eight accessions

Background

Phylogenetic relationship within different genera of Lemnoideae, a kind of small aquatic monocotyledonous plants, was not well resolved, using either morphological characters or traditional markers. Given that rich genetic information in chloroplast genome makes them particularly useful for phylogenetic studies, we used chloroplast genomes to clarify the phylogeny within Lemnoideae.

Methods

DNAs were sequenced with next-generation sequencing. The duckweeds chloroplast genomes were indirectly filtered from the total DNA data, or directly obtained from chloroplast DNA data. To test the reliability of assembling the chloroplast genome based on the filtration of the total DNA, two methods were used to assemble the chloroplast genome of Landoltia punctata strain ZH0202. A phylogenetic tree was built on the basis of the whole chloroplast genome sequences using MrBayes v.3.2.6 and PhyML 3.0.

Results

Eight complete duckweeds chloroplast genomes were assembled, with lengths ranging from 165,775 bp to 171,152 bp, and each contains 80 protein-coding sequences, four rRNAs, 30 tRNAs and two pseudogenes. The identity of L. punctata strain ZH0202 chloroplast genomes assembled through two methods was 100%, and their sequences and lengths were completely identical. The chloroplast genome comparison demonstrated that the differences in chloroplast genome sizes among the Lemnoideae primarily resulted from variation in non-coding regions, especially from repeat sequence variation. The phylogenetic analysis demonstrated that the different genera of Lemnoideae are derived from each other in the following order: Spirodela, Landoltia, Lemna, Wolffiella, and Wolffia.

Discussion

This study demonstrates potential of whole chloroplast genome DNA as an effective option for phylogenetic studies of Lemnoideae. It also showed the possibility of using chloroplast DNA data to elucidate those phylogenies which were not yet solved well by traditional methods even in plants other than duckweeds.

Chloroplast genomes of Lilium lancifolium, L. amabile, L. callosum, and L. philadelphicum: Molecular characterization and their use in phylogenetic analysis in the genus Lilium and other allied genera in the order Liliales

Chloroplast (cp) genomes of Lilium amabile, L. callosum, L. lancifolium, and L. philadelphicum were fully sequenced. Using these four novel cp genome sequences and five other previously sequenced cp genomes, features of the cp genomes were characterized in detail among species in the genus Lilium and other related genera in the order Liliales. The lengths and nucleotide composition showed little variation. No structural variation was found among the cp genomes in Liliales. Gene contents were conserved among four newly sequenced cp genome in Lilium species, the only differences being in two pseudogenes. We identified 112 genes in 13 functional categories, 18 of which carried introns that were conserved among the species in Liliales. There were 16-21 SSR loci (>12 bp, >3 repeats) in the cp genomes in Lilium and the genomic locations of these loci were highly variable among the species. Average mutations were 15 SNPs per 1kb and 5 indels per 1kb, respectively, in the cp genomes of the newly sequenced four Lilium species. Phylogenetic classifications revealed some discrepancies between trees based on the cp genomes and previous classifications based on the morphology and geographic distributions.

A Phylogenetic Analysis of Chloroplast Genomes Elucidates the Relationships of the Six Economically Important Brassica Species Comprising the Triangle of U

The Brassica genus comprises many economically important worldwide cultivated crops. The well-established model of the Brassica genus, U's triangle, consists of three basic diploid plant species (Brassica rapa, Brassica oleracea, and Brassica nigra) and three amphidiploid species (Brassica napus, Brassica juncea, and Brassica carinata) that arose through interspecific hybridizations. Despite being extensively studied because of its commercial relevance, several aspects of the origin of the Brassica species and the relationships within and among these six species still remain open questions. Here, we successfully de novo assembled 60 complete chloroplast genomes of Brassica genotypes of all six species. A complete map of the single nucleotide variants and insertions and deletions in the chloroplast genomes of different Brassica species was produced. The chloroplast genome consists of a Large and a Small Single Copy (LSC and SSC) region between two inverted repeats, and while these regions of chloroplast genomes have very different molecular evolutionary rates, phylogenetic analyses of different regions yielded no contradicting topologies and separated the Brassica genus into four clades. B. carinata and B. juncea share their chloroplast genome with one of their hybridization donors B. nigra and B. rapa, respectively, which fits the U model. B. rapa, surprisingly, shows evidence of two types of chloroplast genomes, with one type specific to some Italian broccoletto accessions. B. napus clearly has evidence for two independent hybridization events, as it contains either B. rapa chloroplast genomes. The divergence estimation suggests that B. nigra and B. carinata diverged from the main Brassica clade 13.7 million years ago (Mya), while B. rapa and B. oleracea diverged at 2.18 Mya. The use of the complete chloroplast DNA sequence not only provides insights into comparative genome analysis but also paves the way for a better understanding of the phylogenetic relationships within the Brassica genus.

Rice Chloroplast Genome Variation Architecture and Phylogenetic Dissection in Diverse Oryza Species Assessed by Whole-Genome Resequencing

BACKGROUND

Chloroplast genome variations have been detected, despite its overall conserved structure, which has been valuable for plant population genetics and evolutionary studies. Here, we described chloroplast variation architecture of 383 rice accessions from diverse regions and different ecotypes, in order to mine the rice chloroplast genome variation architecture and phylogenetic.

RESULTS

A total of 3677 variations across the chloroplast genome were identified with an average density of 27.33 per kb, in which wild rice showing a higher variation density than cultivated groups. Chloroplast genome nucleotide diversity investigation indicated a high degree of diversity in wild rice than in cultivated rice. Genetic distance estimation revealed that African rice showed a low level of breeding and connectivity with the Asian rice, suggesting the big distinction of them. Population structure and principal component analysis revealed the existence of clear clustering of African and Asian rice, as well as the indica and japonica in Asian cultivated rice. Phylogenetic analysis based on maximum likelihood and Bayesian inference methods and the population splits test suggested and supported the independent origins of indica and japonica within Asian cultivated rice. In addition, the African cultivated rice was thought to be domesticated differently from Asian cultivated rice.

CONCLUSIONS

The chloroplast genome variation architecture in Asian and African rice are different, as well as within Asian or African rice. Wild rice and cultivated rice also have distinct nucleotide diversity or genetic distance. In chloroplast level, the independent origins of indica and japonica within Asian cultivated rice were suggested and the African cultivated rice was thought to be domesticated differently from Asian cultivated rice. These results will provide more candidate evidence for the further rice chloroplast genomic and evolution studies.

The complete chloroplast genome sequence of Gentiana lawrencei var. farreri (Gentianaceae) and comparative analysis with its congeneric species

Background

The chloroplast (cp) genome is useful in plant systematics, genetic diversity analysis, molecular identification and divergence dating. The genus Gentiana contains 362 species, but there are only two valuable complete cp genomes. The purpose of this study is to report the characterization of complete cp genome of G. lawrencei var. farreri, which is endemic to the Qinghai-Tibetan Plateau (QTP).

Methods

Using high throughput sequencing technology, we got the complete nucleotide sequence of the G. lawrencei var. farreri cp genome. The comparison analysis including genome difference and gene divergence was performed with its congeneric species G. straminea. The simple sequence repeats (SSRs) and phylogenetics were studied as well.

Results

The cp genome of G. lawrencei var. farreri is a circular molecule of 138,750 bp, containing a pair of 24,653 bp inverted repeats which are separated by small and large single-copy regions of 11,365 and 78,082 bp, respectively. The cp genome contains 130 known genes, including 85 protein coding genes (PCGs), eight ribosomal RNA genes and 37 tRNA genes. Comparative analyses indicated that G. lawrencei var. farreri is 10,241 bp shorter than its congeneric species G. straminea. Four large gaps were detected that are responsible for 85% of the total sequence loss. Further detailed analyses revealed that 10 PCGs were included in the four gaps that encode nine NADH dehydrogenase subunits. The cp gene content, order and orientation are similar to those of its congeneric species, but with some variation among the PCGs. Three genes, ndhB, ndhF and clpP, have high nonsynonymous to synonymous values. There are 34 SSRs in the G. lawrencei var. farreri cp genome, of which 25 are mononucleotide repeats: no dinucleotide repeats were detected. Comparison with the G. straminea cp genome indicated that five SSRs have length polymorphisms and 23 SSRs are species-specific. The phylogenetic analysis of 48 PCGs from 12 Gentianales taxa cp genomes clearly identified three clades, which indicated the potential of cp genomes in phylogenetics.

Discussion

The “missing” sequence of G. lawrencei var. farreri mainly consistent of ndh genes which could be dispensable under chilling-stressed conditions in the QTP. The complete cp genome sequence of G. lawrencei var. farreri provides intragenic information that will contribute to genetic and phylogenetic research in the Gentianaceae.

Computational analysis of atpB gene promoter from different Pakistani apple varieties

Apple is the fourth most important fruit crop grown in temperate areas of the world belongs to the family Rosaceae. In the present study, the promoter (∼1000bp) region of atpB gene was used to evaluate the genetic diversity and phylogeny of six local apple varieties. atpB gene is one of the large chloroplastic region which encodes β-subunit of ATP synthase and previously it had been used largely in phylogenetic studies. During the present study, atpB promoter was amplified, sequenced and analyzed using various bioinformatics tools including Place Signal Scan, MEGA6 and BLASTn. During the phylogenetic analysis, obtained phylogram divided the studied varieties into two clusters revealing the monophyletic origin of studied apple varieties. Pairwise distance revealed moderate genetic diversity that ranges from 0.047-0.170 with an average of 0.101. While identifying different cis-acting elements present in the atpB promoter region, results exhibited the occurrence of 56 common and 20 unique cis-regulatory elements among studied varieties. The identified cis-acting regulatory elements were mapped as well. It was observed that Kala Kulu has the highest unique features with reference to the availability of cis-acting elements. Moreover, the possible functions of all regulatory elements present on the promoter sequence of atpB gene were predicted based on already reported information regarding their in vivo role.

Genomics of crop wild relatives: expanding the gene pool for crop improvement

Plant breeders require access to new genetic diversity to satisfy the demands of a growing human population for more food that can be produced in a variable or changing climate and to deliver the high-quality food with nutritional and health benefits demanded by consumers. The close relatives of domesticated plants, crop wild relatives (CWRs), represent a practical gene pool for use by plant breeders. Genomics of CWR generates data that support the use of CWR to expand the genetic diversity of crop plants. Advances in DNA sequencing technology are enabling the efficient sequencing of CWR and their increased use in crop improvement. As the sequencing of genomes of major crop species is completed, attention has shifted to analysis of the wider gene pool of major crops including CWR. A combination of de novo sequencing and resequencing is required to efficiently explore useful genetic variation in CWR. Analysis of the nuclear genome, transcriptome and maternal (chloroplast and mitochondrial) genome of CWR is facilitating their use in crop improvement. Genome analysis results in discovery of useful alleles in CWR and identification of regions of the genome in which diversity has been lost in domestication bottlenecks. Targeting of high priority CWR for sequencing will maximize the contribution of genome sequencing of CWR. Coordination of global efforts to apply genomics has the potential to accelerate access to and conservation of the biodiversity essential to the sustainability of agriculture and food production.

High-throughput multiplex cpDNA resequencing clarifies the genetic diversity and genetic relationships among Brassica napus, Brassica rapa and Brassica oleracea

Brassica napus (rapeseed) is a recent allotetraploid plant and the second most important oilseed crop worldwide. The origin of B. napus and the genetic relationships with its diploid ancestor species remain largely unresolved. Here, chloroplast DNA (cpDNA) from 488 B. napus accessions of global origin, 139 B. rapa accessions and 49 B. oleracea accessions were populationally resequenced using Illumina Solexa sequencing technologies. The intraspecific cpDNA variants and their allelic frequencies were called genomewide and further validated via EcoTILLING analyses of the rpo region. The cpDNA of the current global B. napus population comprises more than 400 variants (SNPs and short InDels) and maintains one predominant haplotype (Bncp1). Whole-genome resequencing of the cpDNA of Bncp1 haplotype eliminated its direct inheritance from any accession of the B. rapa or B. oleracea species. The distribution of the polymorphism information content (PIC) values for each variant demonstrated that B. napus has much lower cpDNA diversity than B. rapa; however, a vast majority of the wild and cultivated B. oleracea specimens appeared to share one same distinct cpDNA haplotype, in contrast to its wild C-genome relatives. This finding suggests that the cpDNA of the three Brassica species is well differentiated. The predominant B. napus cpDNA haplotype may have originated from uninvestigated relatives or from interactions between cpDNA mutations and natural/artificial selection during speciation and evolution. These exhaustive data on variation in cpDNA would provide fundamental data for research on cpDNA and chloroplasts.

Complete chloroplast and ribosomal sequences for 30 accessions elucidate evolution of Oryza AA genome species

Cytoplasmic chloroplast (cp) genomes and nuclear ribosomal DNA (nR) are the primary sequences used to understand plant diversity and evolution. We introduce a high-throughput method to simultaneously obtain complete cp and nR sequences using Illumina platform whole-genome sequence. We applied the method to 30 rice specimens belonging to nine Oryza species. Concurrent phylogenomic analysis using cp and nR of several of specimens of the same Oryza AA genome species provides insight into the evolution and domestication of cultivated rice, clarifying three ambiguous but important issues in the evolution of wild Oryza species. First, cp-based trees clearly classify each lineage but can be biased by inter-subspecies cross-hybridization events during speciation. Second, O. glumaepatula, a South American wild rice, includes two cytoplasm types, one of which is derived from a recent interspecies hybridization with O. longistminata. Third, the Australian O. rufipogan-type rice is a perennial form of O. meridionalis.

Comprehensive Survey of Genetic Diversity in Chloroplast Genomes and 45S nrDNAs within Panax ginseng Species

We report complete sequences of chloroplast (cp) genome and 45S nuclear ribosomal DNA (45S nrDNA) for 11 Panax ginseng cultivars. We have obtained complete sequences of cp and 45S nrDNA, the representative barcoding target sequences for cytoplasm and nuclear genome, respectively, based on low coverage NGS sequence of each cultivar. The cp genomes sizes ranged from 156,241 to 156,425 bp and the major size variation was derived from differences in copy number of tandem repeats in the ycf1 gene and in the intergenic regions of rps16-trnUUG and rpl32-trnUAG. The complete 45S nrDNA unit sequences were 11,091 bp, representing a consensus single transcriptional unit with an intergenic spacer region. Comparative analysis of these sequences as well as those previously reported for three Chinese accessions identified very rare but unique polymorphism in the cp genome within P. ginseng cultivars. There were 12 intra-species polymorphisms (six SNPs and six InDels) among 14 cultivars. We also identified five SNPs from 45S nrDNA of 11 Korean ginseng cultivars. From the 17 unique informative polymorphic sites, we developed six reliable markers for analysis of ginseng diversity and cultivar authentication.

DNA barcodes from four loci provide poor resolution of taxonomic groups in the genus Crataegus

DNA barcodes can facilitate identification of organisms especially when morphological characters are limited or unobservable. To what extent this potential is realized in specific groups of plants remains to be determined. Libraries of barcode sequences from well-studied authoritatively identified plants represented by herbarium voucher specimens are needed in order for DNA barcodes to serve their intended purpose, where this is possible, and to understand the reasons behind their failure to do so, when this occurs. We evaluated four loci, widely regarded as universal DNA barcodes for plants, for their utility in hawthorn species identification. Three plastid regions, matK, rbcLa and psbA-trnH, and the internal transcribed spacer 2 (ITS2) of nuclear ribosomal DNA discriminate only some of the species of Crataegus that can be recognized on the basis of their morphology etc. This is, in part, because in Rosaceae tribe Maleae most individual plastid loci yield relatively little taxonomic resolution and, in part, because the effects of allopolyploidization have not been eliminated by concerted evolution of the ITS regions. Although individual plastid markers provided generally poor resolution of taxonomic groups in Crataegus, a few species were notable exceptions. In contrast, analyses of concatenated sequences of the 3 plastid barcode loci plus 11 additional plastid loci gave a well-resolved maternal phylogeny. In the ITS2 tree, different individuals of some species formed groups with taxonomically unrelated species. This is a sign of lineage sorting due to incomplete concerted evolution in ITS2. Incongruence between the ITS2 and plastid trees is best explained by hybridization between different lineages within the genus. In aggregate, limited between-species variation in plastid loci, hybridization and a lack of concerted evolution in ITS2 all combine to limit the utility of standard barcoding markers in Crataegus. These results have implications for authentication of hawthorn materials in natural health products.