Plastomes of Mimosoideae: structural and size variation,sequence divergence,and phylogenetic implication

Plastomes of Fabaceae show both significant structural and size variation; however, most published plastomes are from subfamily Papilionoideae and only a few are from the other two subfamilies. In order to address the plastome structural and size variation of subfamily Mimosoideae, we integrated 11 newly sequenced plastomes from representing genera with three previously published ones. Each mimosoid plastome presented a typical quadripartite structure and contained 111 unique genes. Their inverted repeats (IRs) experienced multiple expansion/contraction; a ca. 13-kb IR expansion into small single copy (SSC) was detected in plastomes of a clade formed by tribe Ingeae and Acacia sensu stricto (s.s.), and a ca. 1.7-kb IR expansion into and a ca. 1.9-kb contraction out of large single copy (LSC) were found in Pithecellobium flexicaule and Acacia dealbata, respectively. Linear regression analysis showed decreased synonymous substitution rates of genes relocating from SSC into IR. A loss of both introns of clpP occurred in A. dealbata and Faidherbia albida, and a duplicated clpP copy was detected in A. dealbata. Furthermore, a 421-bp inversion that containing rps18 was found in A. dealbata. The size of mimosoid plastomes was found significantly affected by a IR-SC boundary shift, and also associated with repeat content. Plastome coding and noncoding regions with variable sequence divergence may supply valuable markers for molecular evolutionary and phylogenetic studies at different taxonomic levels. Plastid phylogenomics well resolved relationships among sampled mimosoid species.     

Comparative plastome data analysis of Dendrosicyos socotranus and Corallocarpus boehmii

Comparison of the Dendrosicyos socotranus and Corallocarpus boehmii (tribe Coniandreae, family Cucurbitaceae) plastome data was of interest. Data on RNA, tRNA, GC%, plastome size, CDS and pseudogene were tabulated for the two species. The total length of 1,57,380 bp and 1,58,744 bp which includes LSC, SSC, IRa, and IRb, while their GC content was 37.1% and 37% respectively. The variation in the length of genes e.g. ndhD, ndhI, rpl22, rpoC2, rps16, rps19, rps8, ycf1and ycf2 noted. Data help to document the genetic differences between usual (climber) with those of tree cucurbits.     

Complete sequencing and comparative analyses of the pepper (Capsicum annuum L.) plastome revealed high frequency of tandem repeats and large insertion/deletions on pepper plastome

Plants in the family Solanaceae are used as model systems in comparative and evolutionary genomics. The complete chloroplast genomes of seven solanaceous species have been sequenced, including tobacco, potato and tomato, but not peppers. We analyzed the complete chloroplast genome sequence of the hot pepper, Capsicum annuum. The pepper chloroplast genome was 156,781 bp in length, including a pair of inverted repeats (IR) of 25,783 bp. The content and the order of 133 genes in the pepper chloroplast genome were identical to those of other solanaceous plastomes. To characterize pepper plastome sequence, we performed comparative analysis using complete plastome sequences of pepper and seven solanaceous plastomes. Frequency and contents of large indels and tandem repeat sequences and distribution pattern of genome-wide sequence variations were investigated. In addition, a phylogenetic analysis using concatenated alignments of coding sequences was performed to determine evolutionary position of pepper in Solanaceae. Our results revealed two distinct features of pepper plastome compared to other solanaceous plastomes. Firstly, large indels, including insertions on accD and rpl20 gene sequences, were predominantly detected in the pepper plastome compared to other solanaceous plastomes. Secondly, tandem repeat sequences were particularly frequent in the pepper plastome. Taken together, our study represents unique features of evolution of pepper plastome among solanaceous plastomes.     

Complete Plastid Genome Sequencing of Trochodendraceae Reveals a Significant Expansion of the Inverted Repeat and Suggests a Paleogene Divergence between the Two Extant Species

The early-diverging eudicot order Trochodendrales contains only two monospecific genera, Tetracentron and Trochodendron. Although an extensive fossil record indicates that the clade is perhaps 100 million years old and was widespread throughout the Northern Hemisphere during the Paleogene and Neogene, the two extant genera are both narrowly distributed in eastern Asia. Recent phylogenetic analyses strongly support a clade of Trochodendrales, Buxales, and Gunneridae (core eudicots), but complete plastome analyses do not resolve the relationships among these groups with strong support. However, plastid phylogenomic analyses have not included data for Tetracentron. To better resolve basal eudicot relationships and to clarify when the two extant genera of Trochodendrales diverged, we sequenced the complete plastid genome of Tetracentron sinense using Illumina technology. The Tetracentron and Trochodendron plastomes possess the typical gene content and arrangement that characterize most angiosperm plastid genomes, but both genomes have the same unusual ∼4 kb expansion of the inverted repeat region to include five genes (rpl22, rps3, rpl16, rpl14, and rps8) that are normally found in the large single-copy region. Maximum likelihood analyses of an 83-gene, 88 taxon angiosperm data set yield an identical tree topology as previous plastid-based trees, and moderately support the sister relationship between Buxaceae and Gunneridae. Molecular dating analyses suggest that Tetracentron and Trochodendron diverged between 44-30 million years ago, which is congruent with the fossil record of Trochodendrales and with previous estimates of the divergence time of these two taxa. We also characterize 154 simple sequence repeat loci from the Tetracentron sinense and Trochodendron aralioides plastomes that will be useful in future studies of population genetic structure for these relict species, both of which are of conservation concern.     

The complete plastome of Panax stipuleanatus: Comparative and phylogenetic analyses of the genus Panax (Araliaceae)

Panax stipuleanatus (Araliaceae) is an endangered and medicinally important plant endemic to China. However, phylogenetic relationships within the genus Panax have remained unclear. In this study, we sequenced the complete plastome of P. stipuleanatus and included previously reported Panax plastomes to better understand the relationships between species and plastome evolution within the genus Panax. The plastome of P. stipuleanatus is 156,069 base pairs (bp) in length, consisting of a pair of inverted repeats (IRs, each 25,887 bp) that divide the plastome into a large single copy region (LSC, 86,126 bp) and a small single copy region (SSC, 8169 bp). The plastome contains 114 unigenes (80 protein-coding genes, 30 tRNA genes, and 4 rRNA genes). Comparative analyses indicated that the plastome gene content and order, as well as the expansion/contraction of the IR regions, are all highly conserved within Panax. No significant positive selection in the plastid protein-coding genes was observed across the eight Panax species, suggesting the Panax plastomes may have undergone a strong purifying selection. Our phylogenomic analyses resulted in a phylogeny with high resolution and supports for Panax. Nine proteincoding genes and 10 non-coding regions presented high sequence divergence, which could be useful for identifying different Panax species.     

Do nonasterid holoparasitic flowering plants have plastid genomes?

Past work involving the plastid genome (plastome) of holoparasitic plants has been confined to Scrophulariaceae (or Orobanchaceae) which have truncated plastomes owing to loss of photosynthetic and other genes. Nonasterid holoparasites from Balanophoraceae (Corynaea), Hydnoraceae (Hydnora) and Cytinaceae (Cytinus) were tested for the presence of plastid genes and a plastome. Using PCR, plastid 16S rDNA was successfully amplified and sequenced from the above three holoparasites. The sequence of Cytinus showed 121 single base substitutions relative to Nicotiana (8% of the molecule) whereas higher sequence divergence was observed in Hydnora and Corynaea (287 and 513 changes, respectively). Secondary structural models for these 16S rRNAs show that most changes are compensatory, thus suggesting they are functional. Probes constructed for 16S rDNA and for four plastid-encoded ribosomal protein genes (rps2, rps4, rps7 and rpl16) were used in Southern blots of digested genomic DNA from the three holoparasites. Positive hybridizations were obtained using each of the five probes only for Cytinus. For SmaI digests, all plastid gene probes hybridized to a common fragment ca. 20 kb in length in this species. Taken together, these data provide preliminary evidence suggestive of the retention of highly diverged and truncated plastid genome in Cytinus. The greater sequence divergence for 16S rDNA and the negative hybridization results for Hydnora and Corynaea suggests two possibilities: the loss of typically conserved elements of their plastomes or the complete absence of a plastome.     

The chicken or the egg? Plastome evolution and an independent loss of the inverted repeat in papilionoid legumes

The plastid genome (plastome), while surprisingly constant in gene order and content across most photosynthetic angiosperms, exhibits variability in several unrelated lineages. During the diversification history of the legume family Fabaceae, plastomes have undergone many rearrangements, including inversions, expansion, contraction and loss of the typical inverted repeat (IR), gene loss and repeat accumulation in both shared and independent events. While legume plastomes have been the subject of study for some time, most work has focused on agricultural species in the IR-lacking clade (IRLC) and the plant model Medicago truncatula. The subfamily Papilionoideae, which contains virtually all of the agricultural legume species, also comprises most of the plastome variation detected thus far in the family. In this study three non-papilioniods were included among 34 newly sequenced legume plastomes, along with 33 publicly available sequences, to assess plastome structural evolution in the subfamily. In an effort to examine plastome variation across the subfamily, approximately 20% of the sampling represents the IRLC with the remainder selected to represent the early-branching papilionoid clades. A number of IR-related and repeat-mediated changes were identified and examined in a phylogenetic context. Recombination between direct repeats associated with ycf2 resulted in intraindividual plastome heteroplasmy. Although loss of the IR has not been reported in legumes outside of the IRLC, one genistoid taxon was found to completely lack the typical plastome IR. The role of the IR and non-IR repeats in the progression of plastome change is discussed.     

Characterization of chloroplast DNA microsatellites from Saccharum spp and related species

Microsatellites, or simple sequence repeats (SSRs), and their flanking regions in chloroplast genomes (plastomes) of some species of the family Poaceae were analyzed in silico to look for DNA sequence variations. Comparison of the complete chloroplast DNA sequences (cpDNAs) of sugarcane (Saccharum hybrid cv. SP-80-3280 and S. officinarum cv. NCo310) and related species, Agrostis stolonifera, Brachypodium distachyon, Hordeum vulgare subsp vulgare, Lolium perenne, Oryza nivara, O. sativa subsp indica, O. sativa subsp japonica, Sorghum bicolor, Triticum aestivum, Zea mays, and Z. mays cv. B73, allowed us to examine the organization of chloroplast SSRs (cpSSRs) in genic and intergenic regions. We identified 204 cpSSRs in the sugarcane cpDNA; 22.5% were in genic regions. The ndh, rps, trn, and rpl gene clusters of the chloroplasts had the most repeats. Mononucleotide repeats were the most abundant cpSSRs in these species; however, di-, tri-, tetra-, penta-, and hexanucleotide repeats were also identified. Many base substitutions and deletions/insertions were identified in the cpSSR loci and their flanking regions. Multiple alignments of all cpSSR sequences of Poaceae species made identification of nucleotide variability possible; repeat motifs are not uniformly distributed across the Poaceae plastomes, but are mostly confined to intergenic regions. Phylogeny was determined by maximum parsimony and neighbor-joining inference methods. The cpSSRs of these species were found to be polymorphic. It appears that individual cpSSRs in the Poaceae are stable, at least over short periods of evolutionary time. We conclude that the plastome database can be exploited for phylogenetic analysis and biotechnological development.     

Mutational meltdown or controlled chain reaction: The dynamics of rapid plastome evolution in the hyperdiversity of Poaceae

The study of genomic structural evolution associated with accelerated evolutionary rates that result in avoidance of meltdown and increase biodiversity is becoming ever more possible as the number of available plastomes increases. To more comprehensively analyze rate heterogeneity among monocots and within Poaceae, we sequenced plastomes from four Poaceae species, combined them with publicly available data from ~200 plastomes, and conducted comparative analyses to quantify the pattern of rate heterogeneity between different lineages, functional groups, and periods of evolutionary time. We compared structural differences across the Poaceae to quantify how changes in plastome size correspond to different genomic subunits and the evolution of IR–SC junction boundaries. The substitution rates among ancestral Poaceae were inferred to be exceptionally rapid compared to other monocots but slowed after divergence into extant lineages, which could not be sufficiently explained by positive selection. As such, rapid rates in the ancestral lineage leading to Poaceae might be more closely linked to large-scale structural changes like the loss of ycf1 and ycf2. The total increase in plastome size across Poaceae was positively correlated with the total length of intergenic spacers, tandem repeats, and dispersed repeats as well as large single copy, and inverted repeats (IRs). The continuous evolution of IR–SC junction boundaries was asynchronous with sizes of total genome and subunits across Poaceae. Future work is needed to better understand what factors in ancestral Poaceae evolved to harness such rapid rates of plastome evolution, avoid a mutational meltdown, and escape the stagnation of strong purifying selection as well as if these factors could be utilized to synthetically control rates.     

The complete chloroplast genome sequence of Mahonia bealei (Berberidaceae) reveals a significant expansion of the inverted repeat and phylogenetic relationship with other angiosperms

Mahonia bealei (Berberidaceae) is a frequently-used traditional Chinese medicinal plant with efficient anti-inflammatory ability. This plant is one of the sources of berberine, a new cholesterol-lowering drug with anti-diabetic activity. We have sequenced the complete nucleotide sequence of the chloroplast (cp) genome of M. bealei. The complete cp genome of M. bealei is 164,792 bp in length, and has a typical structure with large (LSC 73,052 bp) and small (SSC 18,591 bp) single-copy regions separated by a pair of inverted repeats (IRs 36,501 bp) of large size. The Mahonia cp genome contains 111 unique genes and 39 genes are duplicated in the IR regions. The gene order and content of M. bealei are almost unarranged which is consistent with the hypothesis that large IRs stabilize cp genome and reduce gene loss-and-gain probabilities during evolutionary process. A large IR expansion of over 12 kb has occurred in M. bealei, 15 genes (rps19, rpl22, rps3, rpl16, rpl14, rps8, infA, rpl36, rps11, petD, petB, psbH, psbN, psbT and psbB) have expanded to have an additional copy in the IRs. The IR expansion rearrangement occurred via a double-strand DNA break and subsequence repair, which is different from the ordinary gene conversion mechanism. Repeat analysis identified 39 direct/inverted repeats 30 bp or longer with a sequence identity ≥ 90%. Analysis also revealed 75 simple sequence repeat (SSR) loci and almost all are composed of A or T, contributing to a distinct bias in base composition. Comparison of protein-coding sequences with ESTs reveals 9 putative RNA edits and 5 of them resulted in non-synonymous modifications in rpoC1, rps2, rps19 and ycf1. Phylogenetic analysis using maximum parsimony (MP) and maximum likelihood (ML) was performed on a dataset composed of 65 protein-coding genes from 25 taxa, which yields an identical tree topology as previous plastid-based trees, and provides strong support for the sister relationship between Ranunculaceae and Berberidaceae. Molecular dating analyses suggest that Ranunculaceae and Berberidaceae diverged between 90 and 84 mya, which is congruent with the fossil records and with recent estimates of the divergence time of these two taxa.     

Phylogenomics and biogeography of Wisteria: Implications on plastome evolution among inverted repeat-lacking clade (IRLC) legumes

The genus Wisteria (Fabaceae) is disjunctly distributed in eastern Asian and eastern North American temperate deciduous forests, and it is widely cultivated around the world as spectacular garden plants. It is a member of inverted repeat-lacking clade (IRLC). The IRLC Species are characterized by the loss of an IR region in their plastomes, which has long been of great interest. In this research, we report whole plastome sequences from all four Wisteria species and a Wisteriopsis japonica, combining these with existing data to explore phylogenetic relationships and biogeography of Wisteria, as well as plastome evolution of IRLC species. Phylogenetic analyses recognized a clade containing Glycyrrhiza–Wisteriopsis–Wisteria as sister to the remaining genera of IRLC. North American Wisteria frutescens and the three Asian species formed reciprocal clades, and Wisteria brachybotrys was sister to Wisteria floribunda and Wisteria sinensis. Wisteria may have originated in Japan near the boundary of the Oligocene and Miocene. The disappearance of Bering Land Bridge in the late Miocene might lead to the Eastern Asian–Eastern North American disjunction of Wisteria. Allopatric speciation of Wisteria between the Japanese archipelago and the Asian continent in the Quaternary increased the species richness of eastern Asia in comparison with eastern North America. Synonymous substitution rates (d_S) of protein-coding genes in the IRLC species were around 2-fold (SC genes) or 11-fold (IR genes) higher than those of non-IRLC species. For both SC and IR genes, herbaceous legumes have around 3-fold higher d_S than woody ones. Both loss of one IR region and herbaceous habit elevated substitution rates of the plastomes.     

Passiflora plastome sequencing reveals widespread genomic rearrangements

Although past studies have included Passiflora among angiosperm lineages with highly rearranged plastid genomes (plastomes), knowledge about plastome organization in the genus is limited. So far only one draft and one complete plastome have been published. Expanded sampling of Passiflora plastomes is needed to understand the extent of the genomic rearrangement in the genus, which is also unusual in having biparental plastid inheritance and plastome‐genome incompatibility. We sequenced 15 Passiflora plastomes using either Illumina paired‐end or shotgun cloning and Sanger sequencing approaches. Assembled plastomes were annotated using Dual Organellar GenoMe Annotator (DOGMA) and tRNAscan‐SE. The Populus trichocarpa plastome was used as a reference to estimate genomic rearrangements in Passiflora by performing whole genome alignment in progressiveMauve. The phylogenetic distribution of rearrangements was plotted on the maximum likelihood tree generated from 64 plastid encoded protein genes. Inverted repeat (IR) expansion/contraction and loss of the two largest hypothetical open reading frames, ycf1 and ycf2, account for most plastome size variation, which ranges from 139 262 base pairs (bp) in P. biflora to 161 494 bp in P. pittieri. Passiflora plastomes have experienced numerous inversions, gene and intron losses along with multiple independent IR expansions and contractions resulting in a distinct organization in each of the three subgenera examined. Each Passiflora subgenus has a unique plastome structure in terms of gene content, order and size. The phylogenetic distribution of rearrangements shows that Passiflora has experienced widespread genomic changes, suggesting that such events may not be reliable phylogenetic markers.     

Complete chloroplast genome sequence of Magnolia kwangsiensis (Magnoliaceae): implication for DNA barcoding and population genetics

Here, we report a completely sequenced plastome using Illumina/Solexa sequencing-by-synthesis (SBS) technology. The plastome of Magnolia kwangsiensis Figlar & Noot. is 159?667 bp in length with a typical quadripartite structure: 88?030 bp large single-copy (LSC) and 18?669 bp small single-copy (SSC) regions, separated by two 26?484 bp inverted repeat (IR) regions. The overall predicted gene number is 129, among which 17 genes are duplicated in IR regions. The plastome of M. kwangsiensis is identical in its gene order to previously published plastomes of magnoliids. Furthermore, the C-to-U type RNA editing frequency of 114 seed plants is positively correlated with plastome GC content and plastome length, whereas plastome length is not correlated with GC content. A total of 16 potential putative barcoding or low taxonomic level phylogenetic study markers in Magnoliaceae were detected by comparing the coding and noncoding regions of the plastome of M. kwangsiensis with that of Liriodendron tulipifera L. At least eight markers might be applied not only to Magnoliaceae but also to other taxa. The 86 mononucleotide cpSSRs that distributed in single-copy noncoding regions are highly valuable to study population genetics and conservation genetics of this endangered rare species.     

Total duplication of the small single copy region in the angiosperm plastome: Rearrangement and inverted repeat instability in Asarum

Premise of the Study

As more plastomes are assembled, it is evident that rearrangements, losses, intergenic spacer expansion and contraction, and syntenic breaks within otherwise functioning plastids are more common than was thought previously, and such changes have developed independently in disparate lineages. However, to date, the magnoliids remain characterized by their highly conserved plastid genomes (plastomes).

Methods

Illumina HiSeq and MiSeq platforms were used to sequence the plastomes of Saruma henryi and those of representative species from each of the six taxonomic sections of Asarum. Sequenced plastomes were compared in a phylogenetic context provided by maximum likelihood and parsimony inferences made using an additional 18 publicly available plastomes from early‐diverging angiosperm lineages.

Key Results

In contrast to previously published magnoliid plastomes and the newly sequenced Saruma henryi plastome published here, Asarum plastomes have undergone extensive disruption and contain extremely lengthy AT‐repeat regions. The entirety of the small single copy region (SSC) of A. canadense and A. sieboldii var. sieboldii has been incorporated into the inverted repeat regions (IR), and the SSC of A. delavayi is only 14 bp long. All sampled Asarum plastomes share an inversion of a large portion of the large single copy region (LSC) such that trnE‐UUC is adjacent to the LSC‐IR boundary.

Conclusions

Plastome divergence in Asarum appears to be consistent with trends seen in highly rearranged plastomes of the monocots and eudicots. We propose that plastome instability in Asarum is due to repetitive motifs that serve as recombinatory substrates and reduce genome stability.     

Evidence for replication slippage in the evolution of Oenothera chloroplast DNA.

Evolutionary relationships of four plastid genomes (plastomes) from different Oenothera species have been assessed by sequence comparisons of two intergenic regions that separate the ribosomal protein genes rpl16, rpl14, and rps8. Sequence changes include base substitutions, the occurrence of a 29-base tandem duplication, and variation in the length of two poly-A stretches. Additions/deletions in chloroplast DNA may not be useful for evolutionary comparisons more distant than these, particularly if the sequences undergo divergence after the initial event, but the length mutations reported here allow a finer resolution of the phylogeny of the closely related Oenothera plastomes than would have been possible if only base substitutions had been considered. Comparisons with the orthogous sequence from tobacco chloroplast DNA indicate the direction of change at most of the sites. The results suggest that plastomes I and II are closely related to each other, as are plastomes III and IV. Replication slippage is proposed as a mechanism to explain the length mutations.     

Comparative analysis of inverted repeats of polypod fern (Polypodiales) plastomes reveals two hypervariable regions

Background

Ferns are large and underexplored group of vascular plants (~ 11 thousands species). The genomic data available by now include low coverage nuclear genomes sequences and partial sequences of mitochondrial genomes for six species and several plastid genomes.

Results

We characterized plastid genomes of three species of Dryopteris, which is one of the largest fern genera, using sequencing of chloroplast DNA enriched samples and performed comparative analysis with available plastomes of Polypodiales, the most species-rich group of ferns. We also sequenced the plastome of Adianthum hispidulum (Pteridaceae). Unexpectedly, we found high variability in the IR region, including duplication of rrn16 in D. blanfordii, complete loss of trnI-GAU in D. filix-mas, its pseudogenization due to the loss of an exon in D. blanfordii. Analysis of previously reported plastomes of Polypodiales demonstrated that Woodwardia unigemmata and Lepisorus clathratus have unusual insertions in the IR region. The sequence of these inserted regions has high similarity to several LSC fragments of ferns outside of Polypodiales and to spacer between tRNA-CGA and tRNA-TTT genes of mitochondrial genome of Asplenium nidus. We suggest that this reflects the ancient DNA transfer from mitochondrial to plastid genome occurred in a common ancestor of ferns. We determined the marked conservation of gene content and relative evolution rate of genes and intergenic spacers in the IRs of Polypodiales. Faster evolution of the four intergenic regions had been demonstrated (trnA- orf42, rrn16-rps12, rps7-psbA and ycf2-trnN).

Conclusions

IRs of Polypodiales plastomes are dynamic, driven by such events as gene loss, duplication and putative lateral transfer from mitochondria.

     

Structural features of a wheat plastome as revealed by complete sequencing of chloroplast DNA

Structural features of the wheat plastome were clarified by comparison of the complete sequence of wheat chloroplast DNA with those of rice and maize chloroplast genomes. The wheat plastome consists of a 134,545-bp circular molecule with 20,703-bp inverted repeats and the same gene content as the rice and maize plastomes. However, some structural divergence was found even in the coding regions of genes. These alterations are due to illegitimate recombination between two short direct repeats and/or replication slippage. Overall comparison of chloroplast DNAs among the three cereals indicated the presence of some hot-spot regions for length mutations. Whereas the region with clustered tRNA genes and that downstream of rbcL showed divergence in a species-specific manner, the deletion patterns of ORFs in the inverted-repeat regions and the borders between the inverted repeats and the small single-copy region support the notion that wheat and rice are related more closely to each other than to maize.     

Plastid Genome Assembly Using Long-read data

Although plastid genome (plastome) structure is highly conserved across most seed plants, investigations during the past two decades have revealed several disparately related lineages that experienced substantial rearrangements. Most plastomes contain a large inverted repeat and two single-copy regions, and a few dispersed repeats; however, the plastomes of some taxa harbour long repeat sequences (>300 bp). These long repeats make it challenging to assemble complete plastomes using short-read data, leading to misassemblies and consensus sequences with spurious rearrangements. Single-molecule, long-read sequencing has the potential to overcome these challenges, yet there is no consensus on the most effective method for accurately assembling plastomes using long-read data. We generated a pipeline, plastid Genome Assembly Using Long-read data (ptGAUL), to address the problem of plastome assembly using long-read data from Oxford Nanopore Technologies (ONT) or Pacific Biosciences platforms. We demonstrated the efficacy of the ptGAUL pipeline using 16 published long-read data sets. We showed that ptGAUL quickly produces accurate and unbiased assemblies using only ~50× coverage of plastome data. Additionally, we deployed ptGAUL to assemble four new Juncus (Juncaceae) plastomes using ONT long reads. Our results revealed many long repeats and rearrangements in Juncus plastomes compared with basal lineages of Poales. The ptGAUL pipeline is available on GitHub: https://github.com/Bean061/ptgaul .     

A set of 100 chloroplast DNA primer pairs to study population genetics and phylogeny in monocotyledons

Chloroplast DNA sequences are of great interest for population genetics and phylogenetic studies. However, only a small set of markers are commonly used. Most of them have been designed for amplification in a large range of Angiosperms and are located in the Large Single Copy (LSC). Here we developed a new set of 100 primer pairs optimized for amplification in Monocotyledons. Primer pairs amplify coding (exon) and non-coding regions (intron and intergenic spacer). They span the different chloroplast regions: 72 are located in the LSC, 13 in the Small Single Copy (SSC) and 15 in the Inverted Repeat region (IR). Amplification and sequencing were tested in 13 species of Monocotyledons: Dioscorea abyssinica, D. praehensilis, D. rotundata, D. dumetorum, D. bulbifera, Trichopus sempervirens (Dioscoreaceae), Phoenix canariensis, P. dactylifera, Astrocaryum scopatum, A. murumuru, Ceroxylon echinulatum (Arecaceae), Digitaria excilis and Pennisetum glaucum (Poaceae). The diversity found in Dioscorea, Digitaria and Pennisetum mainly corresponded to Single Nucleotide Polymorphism (SNP) while the diversity found in Arecaceae also comprises Variable Number Tandem Repeat (VNTR). We observed that the most variable loci (rps15-ycf1, rpl32-ccsA, ndhF-rpl32, ndhG-ndhI and ccsA) are located in the SSC. Through the analysis of the genetic structure of a wild-cultivated species complex in Dioscorea, we demonstrated that this new set of primers is of great interest for population genetics and we anticipate that it will also be useful for phylogeny and bar-coding studies.