The chloroplast genome of Amygdalus L. (Rosaceae) reveals the phylogenetic relationship and divergence time

Background Limited access to genetic information has greatly hindered our understanding of the molecular evolution, phylogeny, and differentiation time of subg. Amygdalus. This study reported complete chloroplast (cp) genome sequences of subg. Amygdalus, which further enriched the available valuable resources of complete cp genomes of higher plants and deepened our understanding of the divergence time and phylogenetic relationships of subg. Amygdalus. Results The results showed that subg. Amygdalus species exhibited a tetrad structure with sizes ranging from 157,736 bp (P. kansuensis) to 158,971 bp (P. davidiana), a pair of inverted repeat regions (IRa/IRb) that ranged from 26,137–26,467 bp, a large single-copy region that ranged from 85,757–86,608 bp, and a small single-copy region that ranged from 19,020–19,133 bp. The average GC content of the complete cp genomes in the 12 species was 36.80%. We found that the structure of the subg. Amygdalus complete cp genomes was highly conserved, and the 12 subg. Amygdalus species had an rps19 pseudogene. There was not rearrangement of the complete cp genome in the 12 subg. Amygdalus species. All 12 subg. Amygdalus species clustered into one clade based on both Bayesian inference and maximum likelihood. The divergence time analyses based on the complete cp genome sequences showed that subg. Amygdalus species diverged approximately 15.65 Mya. Conclusion Our results provide data on the genomic structure of subg. Amygdalus and elucidates their phylogenetic relationships and divergence time.

Phylogenomic Analyses of Hepatica Species and Comparative Analyses Within Tribe Anemoneae (Ranunculaceae)

Hepatica is a small genus of Ranunculaceae with medicinal and horticultural value. We characterized nine complete chloroplast (cp) genomes of Hepatica, which ranged from 159,549 to 161,081 bp in length and had a typical quadripartite structure with a large single-copy region (LSC; 80,270–81,249 bp), a small single-copy region (SSC; 17,029–17,838 bp), and two copies of inverted repeat (IR; 31,008–31,100 bp). The cp genomes of Hepatica possess 76 protein-coding genes (PCGs), 29 tRNAs, and four rRNA genes. Comparative analyses revealed a conserved ca. 5-kb IR expansion in Hepatica and other Anemoneae; moreover, multiple inversion events occurred in Hepatica and its relatives. Analyses of selection pressure (dN/dS) showed that most of the PCGs are highly conserved except for rpl20 and rpl22 in Hepatica falconeri, Hepatica americana, and Hepatica acutiloba. Two genes (rps16 and infA) were identified as pseudogenes in Hepatica. In contrast, rpl32 gene was completely lost. The phylogenetic analyses based on 76 PCGs resolved the phylogeny of Hepatica and its related genera. Non-monophyly of Anemone s.l. indicates that Hepatica should be reclassified as an independent genus. In addition, Hepatica nobilis var. japonica is not closely related to H. nobilis.

Structural rearrangements in plastid genomes of Apiaceae as phylogenetic markers

Apiaceae belong to angiosperm families with frequent plastome structural rearrangements, some of which are generally regarded as synapomorphic for large clades, although typically with limited taxon sampling. Our study aims to improve understanding of the structural rearrangements in plastome within the Tordylieae tribe (ApiaceaeApioideae) with a dense sampling scheme of its species. We showed that presence of psbA pseudogene in inverted repeats near the border with a large single-copy region, which is found in the Tordylieae tribe, may be a clade-specific synapomorphy.

Comparative Chloroplast Genomics of Endangered Euphorbia Species: Insights into Hotspot Divergence, Repetitive Sequence Variation, and Phylogeny

Euphorbia is one of the largest genera in the Euphorbiaceae family, comprising 2000 species possessing commercial, medicinal, and ornamental importance. However, there are very little data available on their molecular phylogeny and genomics, and uncertainties still exist at a taxonomic level. Herein, we sequence the complete chloroplast (cp) genomes of two species, E. larica and E. smithii, of the genus Euphorbia through next-generation sequencing and perform a comparative analysis with nine related genomes in the family. The results revealed that the cp genomes had similar quadripartite structure, gene content, and genome organization with previously reported genomes from the same family. The size of cp genomes ranged from 162,172 to 162,358 bp with 132 and 133 genes, 8 rRNAs, 39 tRNA in E. smithii and E. larica, respectively. The numbers of protein-coding genes were 85 and 86, with each containing 19 introns. The four-junction regions were studied and results reveal that rps19 was present at JLB (large single copy region and inverted repeat b junction) in E. larica where its complete presence was located in the IRb (inverted repeat b) region in E. smithii. The sequence comparison revealed that highly divergent regions in rpoC1, rpocB, ycf3, clpP, petD, ycf1, and ndhF of the cp genomes might provide better understanding of phylogenetic inferences in the Euphorbiaceae and order Malpighiales. Phylogenetic analyses of this study illustrate sister clades of E. smithii with E. tricullii and these species form a monophyletic clade with E. larica. The current study might help us to understand the genome architecture, genetic diversity among populations, and evolutionary depiction in the genera.

Chloroplast Genomes and Comparative Analyses among Thirteen Taxa within Myrsinaceae s.str. Clade (Myrsinoideae, Primulaceae)

The Myrsinaceae s.str. clade is a tropical woody representative in Myrsinoideae of Primulaceae and has ca. 1300 species. The generic limits and alignments of this clade are unclear due to the limited number of genetic markers and/or taxon samplings in previous studies. Here, the chloroplast (cp) genomes of 13 taxa within the Myrsinaceae s.str. clade are sequenced and characterized. These cp genomes are typical quadripartite circle molecules and are highly conserved in size and gene content. Three pseudogenes are identified, of which ycf15 is totally absent from five taxa. Noncoding and large single copy region (LSC) exhibit higher levels of nucleotide diversity (Pi) than other regions. A total of ten hotspot fragments and 796 chloroplast simple sequence repeats (SSR) loci are found across all cp genomes. The results of phylogenetic analysis support the notion that the monophyletic Myrsinaceae s.str. clade has two subclades. Non-synonymous substitution rates (dN) are higher in housekeeping (HK) genes than photosynthetic (PS) genes, but both groups have a nearly identical synonymous substitution rate (dS). The results indicate that the PS genes are under stronger functional constraints compared with the HK genes. Overall, the study provides hypervariable molecular markers for phylogenetic reconstruction and contributes to a better understanding of plastid gene evolution in Myrsinaceae s.str. clade.

Stable and widespread structural heteroplasmy in chloroplast genomes revealed by a new long-read quantification method

The chloroplast genome usually has a quadripartite structure consisting of a large single copy region and a small single copy region separated by two long inverted repeats. It has been known for some time that a single cell may contain at least two structural haplotypes of this structure, which differ in the relative orientation of the single copy regions. However, the methods required to detect and measure the abundance of the structural haplotypes are labour-intensive, and this phenomenon remains understudied. Here we develop a new method, Cp-hap, to detect all possible structural haplotypes of chloroplast genomes of quadripartite structure using long-read sequencing data. We use this method to conduct a systematic analysis and quantification of chloroplast structural haplotypes in 61 land plant species across 19 orders of Angiosperms, Gymnosperms and Pteridophytes. Our results show that there are two chloroplast structural haplotypes which occur with equal frequency in most land plant individuals. Nevertheless, species whose chloroplast genomes lack inverted repeats or have short inverted repeats have just a single structural haplotype. We also show that the relative abundance of the two structural haplotypes remains constant across multiple samples from a single individual plant, suggesting that the process which maintains equal frequency of the two haplotypes operates rapidly, consistent with the hypothesis that flip-flop recombination mediates chloroplast structural heteroplasmy. Our results suggest that previous claims of differences in chloroplast genome structure between species may need to be revisited.Significance StatementChloroplast genome consists of a large single copy region, a small single copy region, and two inverted repeats. Some decades ago, a discovery showed that there are two types chloroplast genome in some plants, which differ the way that the four regions are put together. However, this phenomenon has been shown in just a small number of species, and many open questions remain. Here, we develop a fast method to measure the chloroplast genome structures, based on long-reads. We show that almost all plants we analysed contain two possible genome structures, while a few plants contain only one structure. Our findings hint at the causes of the phenomenon, and provide a convenient new method with which to make rapid progress.

The contraction in inverted repeat regions in the complete plastome sequence ofCressa creticaL

Plastome studies have been the focus of research in plant molecular evolution and systematics.C. creticaL. (Convolvulaceae) is a halophyte, habitat in the ecologically challenged area with high salinity and drought. The complete physical map of plastome revealed that it is 141,419bp long, circular molecule. It contains typical quadripartite structure of large single copy region (LSC 94,808bp), small single copy region (SSC 32,527bp) separated by a pair of inverted repeat regions (IRs 7042bp). This plastome is compared with the complete plastomes of other members of Convolvulaceae showed notable distinctions. An exceptional shift in IRs to SC regions is experienced inC. creticaled to many genes shift in both SC regions and contraction in IRs. The size of IRs reduced to 2 to 4 times as compared to those of the Convolvulaceae members studied. The shifted IRs regions showed remarkable variation in nucleotides patterns. Further, the shift was from the IR boundaries and in between the IR regions led to segment IRs. It is concluded that the shift in IRs may be the strategic move for adaptation in the harsh environment.