The complete chloroplast genome sequence of Pelargonium x hortorum: organization and evolution of the largest and most highly rearranged chloroplast genome of land plants

The chloroplast genome of Pelargonium x hortorum has been completely sequenced. It maps as a circular molecule of 217,942 bp and is both the largest and most rearranged land plant chloroplast genome yet sequenced. It features 2 copies of a greatly expanded inverted repeat (IR) of 75,741 bp each and, consequently, diminished single-copy regions of 59,710 and 6,750 bp. Despite the increase in size and complexity of the genome, the gene content is similar to that of other angiosperms, with the exceptions of a large number of pseudogenes, the recognition of 2 open reading frames (ORF56 and ORF42) in the trnA intron with similarities to previously identified mitochondrial products (ACRS and pvs-trnA), the losses of accD and trnT-ggu and, in particular, the presence of a highly divergent set of rpoA-like ORFs rather than a single, easily recognized gene for rpoA. The 3-fold expansion of the IR (relative to most angiosperms) accounts for most of the size increase of the genome, but an additional 10% of the size increase is related to the large number of repeats found. The Pelargonium genome contains 35 times as many 31 bp or larger repeats than the unrearranged genome of Spinacia. Most of these repeats occur near the rearrangement hotspots, and 2 different associations of repeats are localized in these regions. These associations are characterized by full or partial duplications of several genes, most of which appear to be nonfunctional copies or pseudogenes. These duplications may also be linked to the disruption of at least 1 but possibly 2 or 3 operons. We propose simple models that account for the major rearrangements with a minimum of 8 IR boundary changes and 12 inversions in addition to several insertions of duplicated sequence.     

There large inversions in the chloroplast genomes and one loss of the chloroplast generps16 suggest an early evolutionary split in the genusAdonis (Ranunculaceae)

Detailed chloroplast DNA restriction site maps for two species in the genusAdonis (Ranunculaceae),A. annua andA. vernalis, were constructed using single and double digests and the sizes of these genomes are 151.3 and 156.5 kilobases, respectively. Three inversions were found inAdonis, relative to the gene order in the majority of land plants. These rearrangements represent two different gene orders and mark an ancient split in the evolutionary history of this genus. Gene probes were used in order to map the endpoints of the inversions and the inverted repeat regions. The inverted repeat is approximately 400 base pairs shorter inA. annua than inA. vernalis. Two inversions, 39 kilobases and 24 kilobases in size, occur inA. annua and one inversion, 42 kilobases in size, is present in the remaining investigated species ofAdonis. The generps16 is absent from the chloroplast genome inAdonis annua. Restriction sites for eleven restriction endonucleases were mapped forA. annua, A. vernalis and four additional species ofAdonis and two species ofTrollius. Eighty-six phylogenetically informative sites were analysed cladistically in order to evaluate the main clades withinAdonis.     

Gene rearrangements in Chlamydomonas chloroplast DNAs are accounted for by inversions and by the expansion/contraction of the inverted repeat

To gain insight into the mutational events responsible for the extensive variation of chloroplast DNA (cpDNA) within the green algal genus Chlamydomonas, we have investigated the chloroplast gene organization of Chlamydomonas pitschmannii, a close relative of the interfertile species C. eugametos and C. moewusii whose cpDNAs have been well characterized. At 187 kb, the circular cpDNA of C. pitschmannii is the smallest Chlamydomonas cpDNA yet reported; it is 56 and 105 kb smaller than those of its C. eugametos and C. moewusii counterparts, respectively. Despite this substantial size difference, the arrangement of 77 genes on the C. pitschmannii cpDNA displays only three noticeable differences from the organization of the corresponding genes on the collinear C. eugametos and C. moewusii cpDNAs. These changes in gene order are accounted for by the expansion/contraction of the inverted repeat and one or two inversions in a single-copy region. In land plant cpDNAs, these kinds of events are also responsible for gene rearrangements. The large size difference between the C. pitschmannii and C. eugametos/C. moewusii cpDNAs is mainly attributed to multiple events of deletions/additions as opposed to the usually observed expansion/contraction of the inverted repeat in land plant cpDNAs. We also found that the mitochondrial genome of C. pitschmannii is a circular DNA molecule of 16.5 kb which is 5.5 and 7.5 kb smaller than its C. moewusii and C. eugametos counterparts, respectively.     

Comparative and phylogenetic analyses of Swertia L. (Gentianaceae) medicinal plants (from Qinghai,China) based on complete chloroplast genomes

Swertia L. is a large genus in Swertiinae (Gentianaceae). In China, many Swertia species are used as traditional Tibetan medicines, known as “Zangyinchen” or “Dida”. However, the phylogenetic relationships among Swertia medicinal plants and their wild relatives have remained unclear. In this study, we sequenced and assembled 16 complete chloroplast (cp) genomes of 10 Swertia species, mainly distributed in Qinghai Province, China. The results showed that these species have typical structures and characteristics of plant cp genomes. The sizes of Swertia cp genomes are ranging from 149,488 bp to 154,097 bp. Most Swertia cp genomes presented 134 genes, including 85 protein coding genes, eight rRNA genes, 37 tRNA genes, and four pseudogenes. Furthermore, the GC contents and boundaries of cp genomes are similar among Swertia species. The phylogenetic analyses indicated that Swertia is a complex polyphyletic group. In addition, positive selection was found in psaI and petL genes, indicating the possible adaptation of Qinghai Swertia species to the light environment of the Qinghai-Tibet plateau. These new cp genome data could be further investigated to develop DNA barcodes for Swertia medicinal plants and for additional systematic studies of Swertia and Swertiinae species.     

The chloroplast genomes of Bryopsis plumosa and Tydemania expeditiones (Bryopsidales,Chlorophyta): compact genomes and genes of bacterial origin

Background

Species of Bryopsidales form ecologically important components of seaweed communities worldwide. These siphonous macroalgae are composed of a single giant tubular cell containing millions of nuclei and chloroplasts, and harbor diverse bacterial communities. Little is known about the diversity of chloroplast genomes (cpDNAs) in this group, and about the possible consequences of intracellular bacteria on genome composition of the host. We present the complete cpDNAs of Bryopsis plumosa and Tydemania expeditiones, as well as a re-annotated cpDNA of B. hypnoides, which was shown to contain a higher number of genes than originally published. Chloroplast genomic data were also used to evaluate phylogenetic hypotheses in the Chlorophyta, such as monophyly of the Ulvophyceae (the class in which the order Bryopsidales is currently classified).

Results

Both DNAs are circular and lack a large inverted repeat. The cpDNA of B. plumosa is 106,859 bp long and contains 115 unique genes. A 13 kb region was identified with several freestanding open reading frames (ORFs) of putative bacterial origin, including a large ORF (>8 kb) closely related to bacterial rhs-family genes. The cpDNA of T. expeditiones is 105,200 bp long and contains 125 unique genes. As in B. plumosa, several regions were identified with ORFs of possible bacterial origin, including genes involved in mobile functions (transposases, integrases, phage/plasmid DNA primases), and ORFs showing close similarity with bacterial DNA methyltransferases. The cpDNA of B. hypnoides differs from that of B. plumosa mainly in the presence of long intergenic spacers, and a large tRNA region. Chloroplast phylogenomic analyses were largely inconclusive with respect to monophyly of the Ulvophyceae, and the relationship of the Bryopsidales within the Chlorophyta.

Conclusions

The cpDNAs of B. plumosa and T. expeditiones are amongst the smallest and most gene dense chloroplast genomes in the core Chlorophyta. The presence of bacterial genes, including genes typically found in mobile elements, suggest that these have been acquired through horizontal gene transfer, which may have been facilitated by the occurrence of obligate intracellular bacteria in these siphonous algae.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1418-3) contains supplementary material, which is available to authorized users.     

Unraveling the nuclear and chloroplast genomes of an agar producing red macroalga, Gracilaria changii (Rhodophyta,Gracilariales)

Agar and agarose have wide applications in food and pharmaceutical industries. Knowledge on the genome of red seaweeds that produce them is still lacking. To fill the gap in genome analyses of these red algae, we have sequenced the nuclear and organellar genomes of an agarophyte, Gracilaria changii. The partial nuclear genome sequence of G. changii has a total length of 35.8 Mb with 10,912 predicted protein coding sequences. Only 39.4% predicted proteins were found to have significant matches to protein sequences in SwissProt. The chloroplast genome of G. changii is 183,855 bp with a total of 201 open reading frames (ORFs), 29 tRNAs and 3 rRNAs predicted. Five genes: ssrA, leuC and leuD CP76_p173 (orf139) and pbsA were absent in the chloroplast genome of G. changii. The genome information is valuable in accelerating functional studies of individual genes and resolving evolutionary relationship of red seaweeds.     

The spatial distribution of chloroplast DNA and allozyme polymorphisms within A population of Silene alba (Caryophyllaceae)

Genetic structure within a population of Silene alba was studied using a chloroplast DNA (cpDNA) and six allozyme polymorphisms. A 20 × 65 m area was sampled by constructing a 5 × 5 m grid and determining the genotype of the plant nearest to each grid node. Analysis of the spatial distribution of genotypes by Moran's I and join-counts showed a significant degree of association of like cpDNA haplotypes (I = 0.52, S.N.D. = 2.55) but a random or slightly repulsed distribution of allozyme genotypes. A second sample was taken by collecting all individuals from within a 1-m wide transect established along the perimeter of the grid. Genotype and allele frequencies were calculated by grouping individuals from arbitrary 3-m intervals along the transect, and analyzed by Wright's F statistics. The F_st value calculated for cpDNA (0.875) differed markedly from that based on all allozymes (0.027). Taken together, the results suggest that in this population gene flow in the maternally inherited cpDNA is limited by restricted seed dispersal, whereas gene flow in the nuclear DNA based allozymes is more pervasive owing to the added effects of pollen dispersal. The utility of cpDNA polymorphism for the study of fine-scale gene flow is discussed.     

Complete chloroplast genomes of Saccharum spontaneum,Saccharum officinarum and Miscanthus floridulus (Panicoideae: Andropogoneae) reveal the plastid view on sugarcane origins

Sugarcane (Saccharum hybrid cultivar) ranks among the world's top 10 food crops and annually provides 60–70% of the sugar produced worldwide. Despite its economic importance there has been no large-scale systematics study of genus Saccharum and the existing model of sugarcane origins has remained largely unchallenged for almost 50 years. For the first time, we have assembled the complete plastid genomes of Miscanthus floridulus (first report for this genus), Saccharum spontaneum and Saccharum officinarum allowing us to elucidate the phylogenetic origins of Saccharum s.s. species. We demonstrate that Saccharum s.s. is divided into four species, with S. spontaneum diverging from the remainder of the genus about 1.5 million years ago and S. robustum diverging 750,000 years ago. Two separate lineages, one leading to S. officinarum and the other leading to modern hybrid cultivars diverged from S. robustum 640,000 years ago. These findings overturn all previous hypotheses on sugarcane origins, demonstrating that sugarcane's antecedents could not have arisen by human action. All modern cultivars share a common Polynesian origin, whereas Old World canes, S. barberi and S. sinense, cluster as a distinct S. officinarum lineage. This makes modern cultivars a distinct species of genus Saccharum, and we formally propose the name Saccharum cultum for the ancestor of all lineages currently classified as Saccharum hybrid cultivars.     

珍稀濒危飘带兜兰叶绿体全基因组种内变异研究(附录)

飘带兜兰(Paphiopedilum parishii)分布范围狭窄,仅在中国、缅甸、泰国以及老挝有少量分布。近年来,因生境破坏和人为滥采而导致飘带兜兰野生种群极度缩减。为开发种内多态性的分子标记用于保护生物学研究,该研究对飘带兜兰4个野生个体经测序、组装、注释获得的叶绿体基因组序列,与已公布的飘带兜兰2个个体的叶绿体全基因组序列进行比对,分析飘带兜兰叶绿体基因组的种内差异。结果表明:(1)飘带兜兰叶绿体基因组具有典型被子植物叶绿体基因组环状四分体结构,基因组长度为154 403～154 809 bp,共编码129个基因,包括78个蛋白质编码基因、39个tRNA基因、8个rRNA基因,以及4个假基因。(2)在飘带兜兰6个个体叶绿体基因组中检测到103～107个SSRs(simple sequence repeats)位点,其中21个SSR位点具有多态性。此外,在6个个体叶绿体基因组中还检测到60个长序列重复,包括17～21个正向重复、18～29个反向重复、9～16个回文重复、4～9个互补重复。(3)通过比较6个个体叶绿体基因组序列的核苷酸多样性,共发现70处变异,包括10个SNPs(single nucleotide polymorphism)、60个插入缺失(InDels)。其中,有3个SNP位点发生了非同义替换,导致编码功能基因的氨基酸发生改变; 19个插入缺失多态性较高,具有开发为分子标记的潜力。(4)通过计算核苷酸多样性值(P_i)共发现8个有变异的区域,P_i值为0～0.006 32,其中变异度较大的是rps3-rpl22、trnL-UAC-rpl32、rpoB-trnC-GCA以及ycf4,这些高变区可开发为分子标记用于评估飘带兜兰遗传多样性。(5)系统发生分析结果表明,飘带兜兰6个个体叶绿体基因组序列聚在一起,与长瓣兜兰互为姐妹群。综上表明,飘带兜兰叶绿体基因组的SSRs、长序列重复、SNPs、InDels以及核苷酸序列呈现了足够的种内多样性,可开发成分子标记用于该种的系统演化及保护生物学研究。     

Analysis of mitochondrial and chloroplast genomes in two volvocine algae: Eudorina elegans and Eudorina cylindrica (Volvocaceae,Chlorophyta)

The colonial volvocine algae span the full range of organizational complexity, from four-celled species to multicellular species, and this group of algae is often used for the study of evolution. In recent years, many organelle genomes have been sequenced using the application of next generation sequencing technology; however, only a few organelle genomes have been reported in colonial volvocine algae. In this study, we determined the organelle genomes of Eudorina elegans and Eudorina cylindrica and analysed the organelle genome size, structure and gene content between these volvocine species. This provided useful information to help us understand the composition of colonial volvocine organelle genomes. Based on the chloroplast genome protein-coding genes, we conducted a phylogenomics analysis of the volvocine algae. The result revealed an unexpected phylogenetic relationship, namely, E. elegans is more closely related to Pleodorina starrii than to E. cylindrica. The substitution rate of volvocine algae was then calculated based on organelle genome protein-coding genes; our analysis suggested the possibility that the two Eudorina species may be under similar evolutionary pressure. Lastly, the synteny analysis of the mitochondrial genome showed that gene arrangements and contents are highly conserved in the family Volvocaceae, and the synteny analysis of the chloroplast genome indicated that the genus Eudorina may have experienced genomic changes.     

Complete chloroplast genome sequence of Elodea canadensis and comparative analyses with other monocot plastid genomes

Elodea canadensis is an aquatic angiosperm native to North America. It has attracted great attention due to its invasive nature when transported to new areas in its non-native range. We have determined the complete nucleotide sequence of the chloroplast (cp) genome of Elodea. Taxonomically Elodea is a basal monocot, and only few monocot cp genomes representing early lineages of monocots have been sequenced so far. The genome is a circular double-stranded DNA molecule 156,700bp in length, and has a typical structure with large (LSC 86,194bp) and small (SSC 17,810bp) single-copy regions separated by a pair of inverted repeats (IRs 26,348bp each). The Elodea cp genome contains 113 unique genes and 16 duplicated genes in the IR regions. A comparative analysis showed that the gene order and organization of the Elodea cp genome is almost identical to that of Amborella trichopoda, a basal angiosperm. The structure of IRs in Elodea is unique among monocot species with the whole cp genome sequenced. In Elodea and another monocot Lemna minor the borders between IRs and LSC are located upstream of rps19 gene and downstream of trnH-GUG gene, while in most monocots, IR has extended to include both trnH and rps19 genes. A phylogenetic analysis conducted using Bayesian method, based on the DNA sequences of 81 chloroplast genes from 17 monocot taxa provided support for the placement of Elodea together with Lemna as a basal monocot and the next diverging lineage of monocots after Acorales. In comparison with other monocots, the Elodea cp genome has gone through only few rearrangements or gene losses. IR of Elodea has a unique structure among the monocot species studied so far as its structure is similar to that of a basal angiosperm Amborella. This result together with phylogenetic analyses supports the placement of Elodea as a basal monocot to the next diverging lineage of monocots after Acorales. So far, only few cp genomes representing early lineages of monocots have been sequenced and, therefore, this study provides valuable information about the course of evolution in divergence of monocot lineages.     

Gene copy number variations in adaptive evolution: The genomic distribution of gene copy number variations revealed by genetic mapping and their adaptive role in an undomesticated species,white spruce (Picea glauca)

Gene copy number variation (CNV) has been associated with phenotypic variability in animals and plants, but a genomewide understanding of their impacts on phenotypes is largely restricted to human and agricultural systems. As such, CNVs have rarely been considered in investigations of the genomic architecture of adaptation in wild species. Here, we report on the genetic mapping of gene CNVs in white spruce, which lacks a contiguous assembly of its large genome (~20 Gb), and their relationships with adaptive phenotypic variation. We detected 3,911 gene CNVs including de novo structural variations using comparative genome hybridization on arrays (aCGH) in a large progeny set. We inferred the heterozygosity at CNV loci within parents by comparing haploid and diploid tissues and genetically mapped 82 gene CNVs. Our analysis showed that CNVs were distributed over 10 linkage groups and identified four CNV hotspots that we predict to occur in other species of the Pinaceae. Significant relationships were found between 29 of the gene CNVs and adaptive traits based on regression analyses with timings of bud set and bud flush, and height growth, suggesting a role for CNVs in climate adaptation. The importance of CNVs in adaptive evolution of white spruce was also indicated by functional gene annotations and the clustering of 31% of the mapped adaptive gene CNVs in CNV hotspots. Taken together, these results illustrate the feasibility of studying CNVs in undomesticated species and represent a major step towards a better understanding of the roles of CNVs in adaptive evolution.     

The complete chloroplast genome sequence of Taxus chinensis var. mairei (Taxaceae): loss of an inverted repeat region and comparative analysis with related species

Taxus chinensis var. mairei (Taxaceae) is a domestic variety of yew species in local China. This plant is one of the sources for paclitaxel, which is a promising antineoplastic chemotherapy drugs during the last decade. We have sequenced the complete nucleotide sequence of the chloroplast (cp) genome of T. chinensis var. mairei. The T. chinensis var. mairei cp genome is 129,513 bp in length, with 113 single copy genes and two duplicated genes (trnI-CAU, trnQ-UUG). Among the 113 single copy genes, 9 are intron-containing. Compared to other land plant cp genomes, the T. chinensis var. mairei cp genome has lost one of the large inverted repeats (IRs) found in angiosperms, fern, liverwort, and gymnosperm such as Cycas revoluta and Ginkgo biloba L. Compared to related species, the gene order of T. chinensis var. mairei has a large inversion of ~ 110 kb including 91 genes (from rps18 to accD) with gene contents unarranged. Repeat analysis identified 48 direct and 2 inverted repeats 30 bp long or longer with a sequence identity greater than 90%. Repeated short segments were found in genes rps18, rps19 and clpP. Analysis also revealed 22 simple sequence repeat (SSR) loci and almost all are composed of A or T.     

Polymorphism of the (CTG)n Repeat of the Myotonin Protein Kinase Gene in Populations of the Sakha Republic (Yakutia) and Central Asia

The allele frequency distribution of the (CTG)_n repeat located in the 3′-terminal region of the myotonin protein kinase gene (DMPK) was compared for populations of Yakutia (three ethnogeographic groups of Yakuts, Evenks, Evens, Yukaghirs, and Dolgans) and Central Asia (Kazakhs, Uzbeks, and Uighurs) and other ethnic groups. The populations of the two regions proved to considerably differ from each other: features characteristic of Asian Mongoloids were more distinct in the populations of Yakutia, while the Central Asian populations were closer to European populations. The (CTG)_n allele spectrum of Yakuts was considered in connection with the high incidence of myotonic dystrophy in Yakutia. The results support the hypothesis of the founder effect for the spread of myotonic dystrophy in Yakuts. Data on the (CTG)_n polymorphism were used to estimate the phylogenetic relationships of the populations under study.__________Translated from Molekulyarnaya Biologiya, Vol. 39, No. 3, 2005, pp. 385–393.Original Russian Text Copyright © 2005 by Fedorova, Khusainova, Kutuev, Sukhomyasova, Nikolaeva, Kulichkin, Akhmetova, Salimova, Svyatova, Berezina, Platonov, Khusnutdinova.