The evolutionary processes of mitochondrial and chloroplast genomes differ from those of nuclear genomes

Comparison of chloroplast genomes and phylogenetic analysis of four species in Quercus section Cyclobalanopsis

The identification in Quercus L. species was considered to be difficult all the time. The fundamental phylogenies of Quercus have already been discussed by morphological and molecular means. However, the morphological characteristics of some Quercus groups may not be consistent with the molecular results (such as the group Helferiana), which may lead to blurring of species relationships and prevent further evolutionary researches. To understand the interspecific relationships and phylogenetic positions, we sequenced and assembled the CPGs (160,715 bp-160842 bp) of four Quercus section Cyclobalanopsis species by Illumina pair-end sequencing. The genomic structure, GC content, and IR/SC boundaries exhibited significant conservatism. Six highly variable hotspots were detected in comparison analysis, among which rpoC1, clpP and ycf1 could be used as molecular markers. Besides, two genes (petA, ycf2) were detected to be under positive selection pressure. The phylogenetic analysis showed: Trigonobalanus genus and Fagus genus located at the base of the phylogeny tree; The Quercus genus species were distincted to two clades, including five sections. All Compound Trichome Base species clustered into a single branch, which was in accordance with the results of the morphological studies. But neither of group Gilva nor group Helferiana had formed a monophyly. Six Compound Trichome Base species gathered together in pairs to form three branch respectively (Quercus kerrii and Quercus chungii; Quercus austrocochinchinensis with Quercus gilva; Quercus helferiana and Quercus rex). Due to a low support rate (0.338) in the phylogeny tree, the interspecies relationship between the two branches differentiated by this node remained unclear. We believe that Q. helferiana and Q. kerrii can exist as independent species due to their distance in the phylogeny tree. Our study provided genetic information in Quercus genus, which could be applied to further studies in taxonomy and phylogenetics.

Insights into phylogenetic relationships in Pinus inferred from a comparative analysis of complete chloroplast genomes

BACKGROUND

Pinus is the largest genus of Pinaceae and the most primitive group of modern genera. Pines have become the focus of many molecular evolution studies because of their wide use and ecological significance. However, due to the lack of complete chloroplast genome data, the evolutionary relationship and classification of pines are still controversial. With the development of new generation sequencing technology, sequence data of pines are becoming abundant. Here, we systematically analyzed and summarized the chloroplast genomes of 33 published pine species.

RESULTS

Generally, pines chloroplast genome structure showed strong conservation and high similarity. The chloroplast genome length ranged from 114,082 to 121,530 bp with similar positions and arrangements of all genes, while the GC content ranged from 38.45 to 39.00%. Reverse repeats showed a shrinking evolutionary trend, with IRa/IRb length ranging from 267 to 495 bp. A total of 3,205 microsatellite sequences and 5,436 repeats were detected in the studied species chloroplasts. Additionally, two hypervariable regions were assessed, providing potential molecular markers for future phylogenetic studies and population genetics. Through the phylogenetic analysis of complete chloroplast genomes, we offered novel opinions on the genus traditional evolutionary theory and classification.

CONCLUSION

We compared and analyzed the chloroplast genomes of 33 pine species, verified the traditional evolutionary theory and classification, and reclassified some controversial species classification. This study is helpful in analyzing the evolution, genetic structure, and the development of chloroplast DNA markers in Pinus.

Plastid Phylogenomic Analyses Reveal a Cryptic Species of Ligusticopsis (Apiaceae, Angiosperms)

Ligusticopsis litangensis is identified and described as a cryptic species from Sichuan Province, China. Although the distribution of this cryptic species overlaps with that of Ligusticopsis capillacea and Ligusticopsis dielsiana, the morphological boundaries between them are explicit and have obviously distinguishable characters. The main distinguishing features of the cryptic species are as follows: long conical multi-branched roots, very short pedicels in compound umbels, unequal rays, oblong-globose fruits, 1-2 vittae per furrow and 3-4 vittae on the commissure. The above-mentioned features differ somewhat from other species within the genus Ligusticopsis, but generally coincide with the morphological boundaries defined for the genus Ligusticopsis. To determine the taxonomic position of L. litangensis, we sequenced and assembled the plastomes of L. litangensis and compared them with the plastomes of 11 other species of the genus Ligusticopsis. Notably, both phylogenetic analyses based on ITS sequences and the complete chloroplast genome robustly supported that three accessions of L. litangensis are monophyletic clade and then nested in Ligusticopsis genus. Moreover, the plastid genomes of 12 Ligusticopsis species, including the new species, were highly conserved in terms of gene order, gene content, codon bias, IR boundaries and SSR content. Overall, the integration of morphological, comparative genomic and phylogenetic evidence indicates that Ligusticopsis litangensis actually represents a new species.

Chloroplast genomic comparison provides insights into the evolution of seagrasses

BACKGROUND

Seagrasses are a polyphyletic group of monocotyledonous angiosperms that have evolved to live entirely submerged in marine waters. Thus, these species are ideal for studying plant adaptation to marine environments. Herein, we sequenced the chloroplast (cp) genomes of two seagrass species (Zostera muelleri and Halophila ovalis) and performed a comparative analysis of them with 10 previously published seagrasses, resulting in various novel findings.

RESULTS

The cp genomes of the seagrasses ranged in size from 143,877 bp (Zostera marina) to 178,261 bp (Thalassia hemprichii), and also varied in size among different families in the following order: Hydrocharitaceae > Cymodoceaceae > Ruppiaceae > Zosteraceae. The length differences between families were mainly related to the expansion and contraction of the IR region. In addition, we screened out 2,751 simple sequence repeats and 1,757 long repeat sequence types in the cp genome sequences of the 12 seagrass species, ultimately finding seven hot spots in coding regions. Interestingly, we found nine genes with positive selection sites, including two ATP subunit genes (atpA and atpF), three ribosome subunit genes (rps4, rps7, and rpl20), one photosystem subunit gene (psbH), and the ycf2, accD, and rbcL genes. These gene regions may have played critical roles in the adaptation of seagrasses to diverse environments. In addition, phylogenetic analysis strongly supported the division of the 12 seagrass species into four previously recognized major clades. Finally, the divergence time of the seagrasses inferred from the cp genome sequences was generally consistent with previous studies.

CONCLUSIONS

In this study, we compared chloroplast genomes from 12 seagrass species, covering the main phylogenetic clades. Our findings will provide valuable genetic data for research into the taxonomy, phylogeny, and species evolution of seagrasses.

Evidence from Phylogenomics and Morphology Provide Insights into the Phylogeny, Plastome Evolution, and Taxonomy of Kitagawia

Kitagawia Pimenov is one of the segregate genera of Peucedanum sensu lato within the Apiaceae. The phylogenetic position and morphological delimitation of Kitagawia have been controversial. In this study, we used plastid genome (plastome) and nuclear ribosomal DNA (nrDNA) sequences to reconstruct the phylogeny of Kitagawia, along with comparative plastome and morphological analyses between Kitagawia and related taxa. The phylogenetic results identified that all examined Kitagawia species were divided into Subclade I and Subclade II within the tribe Selineae, and they were all distant from the representative members of Peucedanum sensu stricto. The plastomes of Kitagawia and related taxa showed visible differences in the LSC/IRa junction (JLA) and several hypervariable regions, which separated Subclade I and Subclade II from other taxa. Fruit anatomical and micromorphological characteristics, as well as general morphological characteristics, distinguished the four Kitagawia species within Subclade I from Subclade II and other related genera. This study supported the separation of Kitagawia from Peucedanum sensu lato, confirmed that Kitagawia belongs to Selineae, and two species (K. praeruptora and K. formosana) within Subclade II should be placed in a new genus. We believe that the "core" Kitagawia should be limited to Subclade I, and this genus can be distinguished by the association of a series of morphological characteristics. Overall, our study provides new insights into the phylogeny, plastome evolution, and taxonomy of Kitagawia.

Phylogenomic analysis and development of molecular markers for the determination of twelve plum cultivars (Prunus, Rosaceae)

Background

Plums are one of the most important economic crops of the Rosaceae family and are produced all over the world. China has many local varieties, but the genomic information is limited for genetic studies. Here, we first sequenced, assembled, and analyzed the plastomes of twelve plum cultivars and developed molecular markers to distinguish them.

Results

The twelve plastomes of plum cultivars have a circular structure of 157,863–157,952 bp containing a large single-copy region (LSC) of 86,109–86,287 bp, a small copy region (SSC) of 18,927–19,031 bp, and two inverted repeats (IR) of 26,353–26,387 bp each. The plastomes of plum cultivars encode 131 genes, including 86 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. We detected 50, 54, 54, 53, 53, 50, 54, 54, 54, 49, 50, 54 SSRs in the twelve analyzed varieties, respectively. For repeat sequences, we identified 553 tandem repeats, 204 direct repeats, and 270 palindromic repeats. We also analyzed the expansion/contraction of IR regions. The genes rpl22, rps19, rpl2, ycf1, ndhF, and the trnH span on or near the boundary of IR and single-copy regions. Phylogenetic analysis showed that the twelve cultivars were clustered with the P. salicina and P. domestica. We developed eight markers LZ01 to LZ08 based on whole plastomes and nuclear genes and validated them successfully with six repetitions.

Conclusions

The results obtained here could fill in the blanks of the plastomes of these twelve plum cultivars and provide a wider perspective based on the basis of the plastomes of Prunus to the molecular identification and phylogenetic construction accurately. The analysis from this study provides an important and valuable resource for studying the genetic basis for agronomic and adaptive differentiation of the Prunus species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08965-z.

Infrageneric Plastid Genomes of Cotoneaster (Rosaceae): Implications for the Plastome Evolution and Origin of C. wilsonii on Ulleung Island

Cotoneaster is a taxonomically and ornamentally important genus in the family Rosaceae; however, phylogenetic relationships among its species are complicated owing to insufficient morphological diagnostic characteristics and hybridization associated with polyploidy and apomixis. In this study, we sequenced the complete plastomes of seven Cotoneaster species (C. dielsianus, C. hebephyllus, C. integerrimus, C. mongolicus, C. multiflorus, C. submultiflorus, and C. tenuipes) and included the available complete plastomes in a phylogenetic analysis to determine the origin of C. wilsonii, which is endemic to Ulleung Island, Korea. Furthermore, based on 15 representative lineages within the genus, we carried out the first comparative analysis of Cotoneaster plastid genomes to gain an insight into their molecular evolution. The plastomes were highly conserved, with sizes ranging from 159,595 bp (C. tenuipes) to 160,016 bp (C. hebephyllus), and had a GC content of 36.6%. The frequency of codon usage showed similar patterns among the 15 Cotoneaster species, and 24 of the 35 protein-coding genes were predicted to undergo RNA editing. Eight of the 76 common protein-coding genes, including ccsA, matK, ndhD, ndhF, ndhK, petA, rbcL, and rpl16, were positively selected, implying their potential roles in adaptation and speciation. Of the 35 protein-coding genes, 24 genes (15 photosynthesis-related, seven self-replications, and three others) were found to harbor RNA editing sites. Furthermore, several mutation hotspots were identified, including trnG-UCC/trnR-UCU/atpA and trnT-UGU/trnL-UAA. Maximum likelihood analysis based on 57 representative plastomes of Cotoneaster and two Heteromeles plastomes as outgroups revealed two major lineages within the genus, which roughly correspond to two subgenera, Chaenopetalum and Cotoneaster. The Ulleung Island endemic, C. wilsonii, shared its most recent common ancestor with two species, C. schantungensis and C. zabelii, suggesting its potential origin from geographically close members of the subgenus Cotoneaster, section Integerrimi.

Comparative Chloroplast Genomes of Zosteraceae Species Provide Adaptive Evolution Insights Into Seagrass

Seagrasses are marine flowering plants found in tropical and sub-tropical areas that live in coastal regions between the sea and land. All seagrass species evolved from terrestrial monocotyledons, providing the opportunity to study plant adaptation to sea environments. Here, we sequenced the chloroplast genomes (cpGenomes) of three Zostera species, then analyzed and compared their cpGenome structures and sequence variations. We also performed a phylogenetic analysis using published seagrass chloroplasts and calculated the selection pressure of 17 species within seagrasses and nine terrestrial monocotyledons, as well as estimated the number of shared genes of eight seagrasses. The cpGenomes of Zosteraceae species ranged in size from 143,877 bp (Zostera marina) to 152,726 bp (Phyllospadix iwatensis), which were conserved and displayed similar structures and gene orders. Additionally, we found 17 variable hotspot regions as candidate DNA barcodes for Zosteraceae species, which will be helpful for studying the phylogenetic relationships and interspecies differences between seagrass species. Interestingly, nine genes had positive selection sites, including two ATP subunit genes (atpA and atpF), two ribosome subunit genes (rps4 and rpl20), two DNA-dependent RNA polymerase genes (rpoC1 and rpoC2), as well as accD, clpP, and ycf2. These gene regions may have played key roles in the seagrass adaptation to diverse environments. The Branch model analysis showed that seagrasses had a higher rate of evolution than terrestrial monocotyledons, suggesting that seagrasses experienced greater environmental pressure. Moreover, a branch-site model identified positively selected sites (PSSs) in ccsA, suggesting their involvement in the adaptation to sea environments. These findings are valuable for further investigations on Zosteraceae cpGenomes and will serve as an excellent resource for future studies on seagrass adaptation to sea environments.

Complete Chloroplast Genome Sequence of Fortunella venosa (Champ. ex Benth.) C.C.Huang (Rutaceae): Comparative Analysis, Phylogenetic Relationships, and Robust Support for Its Status as an Independent Species

Fortunella venosa (Rutaceae) is an endangered species endemic to China and its taxonomic status has been controversial. The genus Fortunella contains a variety of important economic plants with high value in food, medicine, and ornamental. However, the placement of Genus Fortunella into Genus Citrus has led to controversy on its taxonomy and Systematics. In this present research, the Chloroplast genome of F. venosa was sequenced using the second-generation sequencing, and its structure and phylogenetic relationship analyzed. The results showed that the Chloroplast genome size of F. venosa was 160,265 bp, with a typical angiosperm four-part ring structure containing a large single copy region (LSC) (87,597 bp), a small single copy region (SSC) (18,732 bp), and a pair of inverted repeat regions (IRa\IRb) (26,968 bp each). There are 134 predicted genes in Chloroplast genome, including 89 protein-coding genes, 8 rRNAs, and 37 tRNAs. The GC-content of the whole Chloroplast genome was 43%, with the IR regions having a higher GC content than the LSC and the SSC regions. There were no rearrangements present in the Chloroplast genome; however, the IR regions showed obvious contraction and expansion. A total of 108 simple sequence repeats (SSRs) were present in the entire chloroplast genome and the nucleotide polymorphism was high in LSC and SSC. In addition, there is a preference for codon usage with the non-coding regions being more conserved than the coding regions. Phylogenetic analysis showed that species of Fortunella are nested in the genus of Citrus and the independent species status of F. venosa is supported robustly, which is significantly different from F. japonica. These findings will help in the development of DNA barcodes that can be useful in the study of the systematics and evolution of the genus Fortunella and the family Rutaceae.

Integrated analysis of three newly sequenced fern chloroplast genomes: Genome structure and comparative analysis

BACKGROUND

Some ferns have medicinal properties and are used in therapeutic interventions. However, the classification and phylogenetic relationships of ferns remain incompletely reported. Considering that chloroplast genomes provide ideal information for species identification and evolution, in this study, three unpublished and one published ferns were sequenced and compared with other ferns to obtain comprehensive information on their classification and evolution.

MATERIALS AND METHODS

The complete chloroplast genomes of Dryopteris goeringiana (Kunze) Koidz, D. crassirhizoma Nakai, Athyrium brevifrons Nakai ex Kitagawa, and Polystichum tripteron (Kunze) Presl were sequenced using the Illumina HiSeq 4,000 platform. Simple sequence repeats (SSRs), nucleotide diversity analysis, and RNA editing were investigated in all four species. Genome comparison and inverted repeats (IR) boundary expansion and contraction analyses were also performed. The relationships among the ferns were studied by phylogenetic analysis based on the whole chloroplast genomes.

RESULTS

The whole chloroplast genomes ranged from 148,539 to 151,341 bp in size and exhibited typical quadripartite structures. Ten highly variable loci with parsimony informative (Pi) values of > 0.02 were identified. A total of 75-108 SSRs were identified, and only six SSRs were present in all four ferns. The SSRs contained a higher number of A + T than G + C bases. C-to-U conversion was the most common type of RNA editing event. Genome comparison analysis revealed that single-copy regions were more highly conserved than IR regions. IR boundary expansion and contraction varied among the four ferns. Phylogenetic analysis showed that species in the same genus tended to cluster together with and had relatively close relationships.

CONCLUSION

The results provide valuable information on fern chloroplast genomes that will be useful to identify and classify ferns, and study their phylogenetic relationships and evolution.

Phylogeny and evolution of chloroplast tRNAs in Adoxaceae

Chloroplasts are semiautonomous organelles found in photosynthetic plants. The major functions of chloroplasts include photosynthesis and carbon fixation, which are mainly regulated by its circular genomes. In the highly conserved chloroplast genome, the chloroplast transfer RNA genes (cp tRNA) play important roles in protein translation within chloroplasts. However, the evolution of cp tRNAs remains unclear. Thus, in the present study, we investigated the evolutionary characteristics of chloroplast tRNAs in five Adoxaceae species using 185 tRNA gene sequences. In total, 37 tRNAs encoding 28 anticodons are found in the chloroplast genome in Adoxaceae species. Some consensus sequences are found within the Ψ-stem and anticodon loop of the tRNAs. Some putative novel structures were also identified, including a new stem located in the variable region of tRNATyr in a similar manner to the anticodon stem. Furthermore, phylogenetic and evolutionary analyses indicated that synonymous tRNAs may have evolved from multiple ancestors and frequent tRNA duplications during the evolutionary process may have been primarily caused by positive selection and adaptive evolution. The transition and transversion rates are uneven among different tRNA isotypes. For all tRNAs, the transition rate is greater with a transition/transversion bias of 3.13. Phylogenetic analysis of cp tRNA suggested that the type I introns in different taxa (including eukaryote organisms and cyanobacteria) share the conserved sequences "U-U-x2-C" and "U-x-G-x2-T," thereby indicating the diverse cyanobacterial origins of organelles. This detailed study of cp tRNAs in Adoxaceae may facilitate further investigations of the evolution, phylogeny, structure, and related functions of chloroplast tRNAs.

Plastid Genomes of Flowering Plants: Essential Principles

The plastid genome (plastome ) has proved a valuable source of data for evaluating evolutionary relationships among angiosperms. Through basic and applied approaches, plastid transformation technology offers the potential to understand and improve plant productivity, providing food, fiber, energy, and medicines to meet the needs of a burgeoning global population. The growing genomic resources available to both phylogenetic and biotechnological investigations is allowing novel insights and expanding the scope of plastome research to encompass new species. In this chapter, we present an overview of some of the seminal and contemporary research that has contributed to our current understanding of plastome evolution and attempt to highlight the relationship between evolutionary mechanisms and the tools of plastid genetic engineering.

Complete Plastid Genome Sequencing of Eight Species from Hansenia, Haplosphaera and Sinodielsia (Apiaceae): Comparative Analyses and Phylogenetic Implications

Hansenia Turcz., Haplosphaera Hand.-Mazz. and Sinodielsia H.Wolff are three Apiaceae genera endemic to the Hengduan Mountains and the Himalayas, which usually inhabit elevations greater than 2000 m. The phylogenetic relationships between and within the genera were uncertain, especially the placement of Hap. himalayensis and S. microloba. Therefore, we aimed to conduct comparative (simple sequence repeat (SSR) structure, codon usage bias, nucleotide diversity (Pi) and inverted repeat (IR) boundaries) and phylogenetic analyses of Hansenia, Haplosphaera and Sinodielsia (also compared with Chamaesium and Bupleurum) to reduce uncertainties in intergeneric and interspecific relationships. We newly assembled eight plastid genomes from Hansenia, Haplosphaera and Sinodielsia species, and analyzed them with two plastid genomes from GenBank of Hap. phaea,S. yunnanensis. Phylogenetic analyses used these ten genomes and another 22 plastid genome sequences of Apiaceae. We found that the newly assembled eight genomes ranged from 155,435 bp to 157,797 bp in length and all had a typical quadripartite structure. Fifty-five to 75 SSRs were found in Hansenia, Haplosphaera and Sinodielsia species, and the most abundant SSR was mononucleotide, which accounted for 58.47% of Hansenia, 60.21% of Haplosphaera and 48.01% of Sinodielsia. There was no evident divergence of codon usage frequency between the three genera, where codons ranged from 21,134 to 21,254. The Pi analysis showed that trnE(UUC)-trnT(GGU), trnH(GUG)-psbA and trnE(UUC)-trnT(GGU) spacer regions had the highest Pi values in the plastid genomes of Hansenia (0.01889), Haplosphaera (0.04333) and Sinodielsia (0.01222), respectively. The ndhG-ndhI spacer regions were found in all three genera to have higher diversity values (Pi values: 0.01028-0.2), and thus may provide potential DNA barcodes in phylogenetic analysis. IR boundary analysis showed that the length of rps19 and ycf1 genes entering IRs were usually stable in the same genus. Our phylogenetic tree demonstrated that Hap. himalayensis is sister to Han. weberbaueriana; meanwhile, Haplosphaera and Hansenia are nested together in the East Asia clade, and S. microloba is nested within individuals of S. yunnanensis in the Acronema clade. This study will enrich the complete plastid genome dataset of the Apiaceae genera and has provided a new insight into phylogeny reconstruction using complete plastid genomes of Hansenia, Haplosphaera and Sinodielsia.

Plastome evolution and phylogenetic relationships among Malvaceae subfamilies

Malvaceae s.l. is the largest family of Malvales, comprising more than 4225 species. Within Malvaceae, the phylogenetic relationships among subfamilies remain controversial. To resolve these relationships and explore plastome evolution in Malvaceae, we assembled a complete plastome data set of 39 Malvales species, including newly reported sequences for 13 Malvaceae and two Dipterocarpaceae species. All species possessed the typical quadripartite structure of angiosperm plastomes, but significant independent expansions of the Inverted Repeat regions were detected in Abelmoschus esculentus and Durio zibethinus. Nine coding sequences were identified with positively selected sites in Malvaceae. Several highly variable noncoding and coding regions were identified in the plastomes of Malvaceae that may be valuable for phylogenetic reconstruction at lower taxonomic levels. Phylogenetic reconstructions based on 78 protein-coding genes strongly supported nearly all relationships among Malvaceae subfamilies. The diversification of the subfamilies of Malvaceae was dated to the late Cretaceous and early Eocene, during a time of global warmth.

Analyses of Hybrid Viability across a Hybrid Zone between Two Alnus Species Using Microsatellites and cpDNA Markers

Diploid Alnus glutinosa s. str. and autotetraploid A. rohlenae form a narrow hybrid zone in a study area in southern Serbia, which results in triploid hybrid formation. The vast majority of previous studies have been focused on studies of maternal plants, but the offspring resulting from their crossing have not been much studied. Here, we use the variability of microsatellites and chloroplast DNA between these species and their putative hybrids to create an overall picture of the development of the hybrid zone and its predicted type. To elucidate the gene transfer within both species, the origins of individual ploidies and especially the role of triploid hybrids, a germination experiment was carried out linked with a flow cytometry study of the resulting seedlings. The tension zone model seems to offer the most adequate explanation of our observations, with selection against triploid hybrids and the spatial positioning of the hybrid zone. Despite selection against them, the triploid hybrids play an important role in the exchange of genes between the two species and therefore serve as a bridge for introgression. The presence of fertile triploids is essential for enriching the haplotype diversity between these species and for the development of new genetic lineages.

Determination of the evolutionary pressure on Camellia oleifera on Hainan Island using the complete chloroplast genome sequence

Camellia oleifera is one of the four largest woody edible oil plants in the world with high ecological and medicinal values. Due to frequent interspecific hybridization, it was difficult to study its genetics and evolutionary history. This study used C. oleifera that was collected on Hainan Island to conduct our research. The unique island environment makes the quality of tea oil higher than that of other species grown in the mainland. Moreover, a long-term geographic isolation might affect gene structure. In order to better understand the molecular biology of this species, protect excellent germplasm resources, and promote the population genetics and phylogenetic studies of Camellia plants, high-throughput sequencing technology was used to obtain the chloroplast genome sequence of Hainan C. oleifera. The results showed that the whole chloroplast genome of C. oleifera in Hainan was 156,995 bp in length, with a typical quadripartite structure of a large single copy (LSC) region of 86,648 bp, a small single copy (SSC) region of 18,297 bp, and a pair of inverted repeats (IRs) of 26,025 bp. The whole genome encoded a total of 141 genes (115 different genes), including 88 protein-coding genes, 45 tRNA genes, and eight rRNA genes. Among these genes, nine genes contained one intron, two genes contained two introns, and four overlapping genes were also detected. The total GC content of Hainan C. oleifera's chloroplast genome was 37.29%. The chloroplast genome structure characteristics of Hainan C. oleifera were compared with mainland C. oleifera and those of the other eight closely related Theaceae species; it was found that the contractions and expansions of the IR/LSC and IR/SSC regions affected the length of chloroplast genome. The chloroplast genome sequences of these Theaceae species were highly similar. A comparative analysis indicated that the Theaceae species were conserved in structure and evolution. A total of 51 simple sequence repeat (SSR) loci were detected in the chloroplast genome of Hainan C. oleifera, and all Camellia plants did not have pentanucleotide repeats, which could be used as a good marker in phylogenetic studies. We also detected seven long repeats, the base composition of all repeats was biased toward A/T, which was consistent with the codon bias. It was found that Hainan C. oleifera had a similar evolutionary relationship with C. crapnelliana, through the use of codons and phylogenetic analysis. This study can provide an effective genomic resource for the evolutionary history of Theaceae family.

Heteroplasmy and Patterns of Cytonuclear Linkage Disequilibrium in Wild Carrot

Organellar genomes are considered to be strictly uniparentally-inherited. Uniparental inheritance allows for cytonuclear coevolution and the development of highly coordinated cytonuclear interactions. Yet, instances of biparental inheritance have been documented across eukaryotes. Biparental inheritance in otherwise uniparentally-inherited organelles is termed leakage (maternal or paternal) and allows for the presence of multiple variants of the same organellar genome within an individual, called heteroplasmy. It is unclear what, if any, evolutionary consequences are placed on nuclear and/or organellar genomes due to heteroplasmy. One way of accessing cytonuclear interactions and potential coevolution is through calculating cytonuclear linkage disequilibrium (cnLD), or the non-random association of alleles between nuclear and organellar genomes. Patterns of cnLD can indicate positive or negative cytonuclear selection, coevolution between the nuclear and organellar genomes, non-traditional organellar inheritance, or instances of ancestral heteroplasmy. In plants, cytonuclear interactions have been shown to play a role in cytoplasmic male sterility which occurs in gynodioecious species and is associated with leakage. We used the gynodioecious species, Daucus carota L. spp. carota, or wild carrot, to investigate cnLD. We genotyped a total of 265 individuals from two regions of the USA at 15 nuclear microsatellites, the mitochondrial genes cox1 and atp9, and an intergenic region between trnS and trnG (StoG) in the plastid genome to calculate nuclear–nuclear LD (nucLD), cnLD, and organellar LD (i.e., within the mtDNA and between mtDNA and ptDNA) within the two regions. We were further able to identify cox1 and StoG heteroplasmy and calculate some of the same LD measures within heteroplasmic and homoplasmic (non-heteroplasmic) datasets. We used a Z-transformation test to demonstrate that heteroplasmic individuals display significantly higher levels of cnLD within both regions. In spite of this, within and between organellar LD is low to moderate. Given these patterns of LD in two regions of the USA in which gene flow has been shown to occur between crop and wild carrot, we suggest that heteroplasmy is an evolutionary mechanism which permits the maintenance of cnLD while also acting to disrupt organellar LD.

Complete chloroplast genome sequence of Dryopteris fragrans (L.) Schott and the repeat structures against the thermal environment

Dryopteris fragrans (L.) Schott is a fern growing on the surface of hot rocks and lava. It is exposed to sunlight directly and bears local hot environment. We sequenced the complete nucleotide sequence of its chloroplast (cp) genome. The cp genome was 151,978 bp in length, consisting of a large single-copy region (85,332 bp), a small single-copy region (31,947 bp) and a pair of inverted repeats (17,314 bp). The cp genome contained 112 genes and 345 RNA editing sites in protein-coding genes. Simple sequence repeats (SSRs) and long repeat structure pairs (30–55 bp) were identified. The number and percent of repeat structures are extremely high in ferns. Thermal denaturation experiments showed its cp genome to have numerous, dispersed and high GC percent repeat structures, which conferred the strongest thermal stability. This repeat-heavy genome may provide the molecular basis of how D. fragrans cp survives its hot environment.

Comparative Chloroplast Genomics of Dipsacales Species: Insights Into Sequence Variation, Adaptive Evolution, and Phylogenetic Relationships

In general, the chloroplast genomes of angiosperms are considered to be highly conserved and affected little by adaptive evolution. In this study, we tested this hypothesis based on sequence differentiation and adaptive variation in the plastid genomes in the order Dipsacales. We sequenced the plastid genomes of one Adoxaceae species and six Caprifoliaceae species, and together with seven previously released Dipsacales chloroplasts, we determined the sequence variations, evolutionary divergence of the plastid genomes, and phylogeny of Dipsacales species. The chloroplast genomes of Adoxaceae species ranged in size from 157,074 bp (Sinadoxa corydalifolia) to 158,305 bp (Sambucus williamsii), and the plastid genomes of Caprifoliaceae varied from 154,732 bp (Lonicera fragrantissima var. lancifolia) to 156,874 bp (Weigela florida). The differences in the number of genes in Caprifoliaceae and Adoxaceae species were largely due to the expansion and contraction of inverted repeat regions. In addition, we found that the number of dispersed repeats (Adoxaceae = 37; Caprifoliaceae = 384) was much higher than that of tandem repeats (Adoxaceae = 34; Caprifoliaceae = 291) in Dipsacales species. Interestingly, we determined 19 genes with positive selection sites, including three genes encoding ATP protein subunits (atpA, atpB, and atpI), four genes for ribosome protein small subunits (rps3, rps7, rps14, and rps15), four genes for photosystem protein subunits (psaA, psaJ, psbC, and pabK), two genes for ribosome protein large subunits (rpl22 and rpl32), and the clpP, infA, matK, rbcL, ycf1, and ycf2 genes. These gene regions may have played key roles in the adaptation of Dipsacales to diverse environments. In addition, phylogenetic analysis based on the plastid genomes strongly supported the division of 14 Dipsacales species into two previously recognized sections. The diversification of Adoxaceae and Caprifoliaceae was dated to the late Cretaceous and Tertiary periods. The availability of these chloroplast genomes provides useful genetic information for studying taxonomy, phylogeny, and species evolution in Dipsacales.

Chimeric mitochondrial peptides from contiguous regular and swinger RNA

Previous mass spectrometry analyses described human mitochondrial peptides entirely translated from swinger RNAs, RNAs where polymerization systematically exchanged nucleotides. Exchanges follow one among 23 bijective transformation rules, nine symmetric exchanges (X ↔ Y, e.g. A ↔ C) and fourteen asymmetric exchanges (X → Y → Z → X, e.g. A → C → G → A), multiplying by 24 DNA's protein coding potential. Abrupt switches from regular to swinger polymerization produce chimeric RNAs. Here, human mitochondrial proteomic analyses assuming abrupt switches between regular and swinger transcriptions, detect chimeric peptides, encoded by part regular, part swinger RNA. Contiguous regular- and swinger-encoded residues within single peptides are stronger evidence for translation of swinger RNA than previously detected, entirely swinger-encoded peptides: regular parts are positive controls matched with contiguous swinger parts, increasing confidence in results. Chimeric peptides are 200 × rarer than swinger peptides (3/100,000 versus 6/1000). Among 186 peptides with > 8 residues for each regular and swinger parts, regular parts of eleven chimeric peptides correspond to six among the thirteen recognized, mitochondrial protein-coding genes. Chimeric peptides matching partly regular proteins are rarer and less expressed than chimeric peptides matching non-coding sequences, suggesting targeted degradation of misfolded proteins. Present results strengthen hypotheses that the short mitogenome encodes far more proteins than hitherto assumed. Entirely swinger-encoded proteins could exist.

•

Chimeric peptides are translated from contiguous regular and swinger RNA

•

They are 200x rarer than mitochondrial swinger peptides

•

Chimeric peptides integrated in regular mitochondrial proteins are downregulated

•

Contiguous regular parts are matched positive controls for swinger parts

•

The last point validates results beyond other statistical tests for robustness

Multiplexed shotgun sequencing reveals congruent three-genome phylogenetic signals for four botanical sections of the flax genus Linum

A genome-wide detection of phylogenetic signals by next generation sequencing (NGS) has recently emerged as a promising genomic approach for phylogenetic analysis of non-model organisms. Here we explored the use of a multiplexed shotgun sequencing method to assess the phylogenetic relationships of 18 Linum samples representing 16 species within four botanical sections of the flax genus Linum. The whole genome DNAs of 18 Linum samples were fragmented, tagged, and sequenced using an Illumina MiSeq. Acquired sequencing reads per sample were further separated into chloroplast, mitochondrial and nuclear sequence reads. SNP calls upon genome-specific sequence data sets revealed 6143 chloroplast, 2673 mitochondrial, and 19,562 nuclear SNPs. Phylogenetic analyses based on three-genome SNP data sets with and without missing observations showed congruent three-genome phylogenetic signals for four botanical sections of the Linum genus. Specifically, two major lineages showing a separation of Linum-Dasylinum sections and Linastrum-Syllinum sections were confirmed. The Linum section displayed three major branches representing two major evolutionary stages leading to cultivated flax. Cultivated flax and its immediate progenitor were formed as its own branch, genetically more closely related to L. decumbens and L. grandiflorum with chromosome count of eight, and distantly apart from six other species with chromosome count of nine. Five species of the Linastrum and Syllinum sections were genetically more distant from cultivated flax, but they appeared to be more closely related to each other, even with variable chromosome counts. These findings not only provide the first evidence of congruent three-genome phylogenetic pathways within the Linum genus, but also demonstrate the utility of the multiplexed shotgun sequencing in acquisition of three-genome phylogenetic signals of non-model organisms.

Complete sequencing of five araliaceae chloroplast genomes and the phylogenetic implications

Background

The ginseng family (Araliaceae) includes a number of economically important plant species. Previously phylogenetic studies circumscribed three major clades within the core ginseng plant family, yet the internal relationships of each major group have been poorly resolved perhaps due to rapid radiation of these lineages. Recent studies have shown that phyogenomics based on chloroplast genomes provides a viable way to resolve complex relationships.

Methodology/Principal Findings

We report the complete nucleotide sequences of five Araliaceae chloroplast genomes using next-generation sequencing technology. The five chloroplast genomes are 156,333–156,459 bp in length including a pair of inverted repeats (25,551–26,108 bp) separated by the large single-copy (86,028–86,566 bp) and small single-copy (18,021–19,117 bp) regions. Each chloroplast genome contains the same 114 unique genes consisting of 30 transfer RNA genes, four ribosomal RNA genes, and 80 protein coding genes. Gene size, content, and order, AT content, and IR/SC boundary structure are similar among all Araliaceae chloroplast genomes. A total of 140 repeats were identified in the five chloroplast genomes with palindromic repeat as the most common type. Phylogenomic analyses using parsimony, likelihood, and Bayesian inference based on the complete chloroplast genomes strongly supported the monophyly of the Asian Palmate group and the Aralia-Panax group. Furthermore, the relationships among the sampled taxa within the Asian Palmate group were well resolved. Twenty-six DNA markers with the percentage of variable sites higher than 5% were identified, which may be useful for phylogenetic studies of Araliaceae.

Conclusion

The chloroplast genomes of Araliaceae are highly conserved in all aspects of genome features. The large-scale phylogenomic data based on the complete chloroplast DNA sequences is shown to be effective for the phylogenetic reconstruction of Araliaceae.

Comparative analyses of plastid sequences between native and introduced populations of aquatic weeds Elodea canadensis and E. nuttallii

Non-indigenous species (NIS) are species living outside their historic or native range. Invasive NIS often cause severe environmental impacts, and may have large economical and social consequences. Elodea (Hydrocharitaceae) is a New World genus with at least five submerged aquatic angiosperm species living in fresh water environments. Our aim was to survey the geographical distribution of cpDNA haplotypes within the native and introduced ranges of invasive aquatic weeds Elodea canadensis and E. nuttallii and to reconstruct the spreading histories of these invasive species. In order to reveal informative chloroplast (cp) genome regions for phylogeographic analyses, we compared the plastid sequences of native and introduced individuals of E. canadensis. In total, we found 235 variable sites (186 SNPs, 47 indels and two inversions) between the two plastid sequences consisting of 112,193 bp and developed primers flanking the most variable genomic areas. These 29 primer pairs were used to compare the level and pattern of intraspecific variation within E. canadensis to interspecific variation between E. canadensis and E. nuttallii. Nine potentially informative primer pairs were used to analyze the phylogeographic structure of both Elodea species, based on 70 E. canadensis and 25 E. nuttallii individuals covering native and introduced distributions. On the whole, the level of variation between the two Elodea species was 53% higher than that within E. canadensis. In our phylogeographic analysis, only a single haplotype was found in the introduced range in both species. These haplotypes H1 (E. canadensis) and A (E. nuttallii) were also widespread in the native range, covering the majority of native populations analyzed. Therefore, we were not able to identify either the geographic origin of the introduced populations or test the hypothesis of single versus multiple introductions. The divergence between E. canadensis haplotypes was surprisingly high, and future research may clarify mechanisms that structure native E. canadensis populations.

Microhomology-mediated and nonhomologous repair of a double-strand break in the chloroplast genome of Arabidopsis

Chloroplast DNA (cpDNA) is under great photooxidative stress, yet its evolution is very conservative compared with nuclear or mitochondrial genomes. It can be expected that DNA repair mechanisms play important roles in cpDNA survival and evolution, but they are poorly understood. To gain insight into how the most severe form of DNA damage, a double-strand break (DSB), is repaired, we have developed an inducible system in Arabidopsis that employs a psbA intron endonuclease from Chlamydomonas, I-CreII, that is targeted to the chloroplast using the rbcS1 transit peptide. In Chlamydomonas, an I-CreII-induced DSB in psbA was repaired, in the absence of the intron, by homologous recombination between repeated sequences (20-60 bp) abundant in that genome; Arabidopsis cpDNA is very repeat poor, however. Phenotypically strong and weak transgenic lines were examined and shown to correlate with I-CreII expression levels. Southern blot hybridizations indicated a substantial loss of DNA at the psbA locus, but not cpDNA as a whole, in the strongly expressing line. PCR analysis identified deletions nested around the I-CreII cleavage site indicative of DSB repair using microhomology (6-12 bp perfect repeats, or 10-16 bp with mismatches) and no homology. These results provide evidence of alternative DSB repair pathways in the Arabidopsis chloroplast that resemble the nuclear, microhomology-mediated and nonhomologous end joining pathways, in terms of the homology requirement. Moreover, when taken together with the results from Chlamydomonas, the data suggest an evolutionary relationship may exist between the repeat structure of the genome and the organelle's ability to repair broken chromosomes.

Comparative phylogeography of two closely related Viola species occurring in contrasting habitats in the Japanese archipelago

Similar distribution ranges shared by closely related plant species may have been shaped through different migration histories if those species have differing habitat preference. To test this hypothesis, phylogeographical patterns and population genetic structures were compared between two sister Viola species: Viola eizanensis preferring woodland and V. chaerophylloides var. sieboldiana preferring grassland, both being native to the Japanese Archipelago. Amplified fragment length polymorphism (AFLP) was used for phylogenetic reconstruction, together with Bayesian ancestry analysis, AMOVA, analysis of genetic diversity statistics, and analysis of the relative contribution of each population to total diversity. The results indicated that V. eizanensis had two distinct lineages occurring in the western and eastern part of Japan, but such lineages are not distinct in V. chaerophylloides var. sieboldiana. Both species exhibited the low genetic diversity and high between-population differentiation typical of selfing plants. In V. chaerophylloides var. sieboldiana, one particular population made a significantly higher contribution to the total heterozygosity (H (T)), whereas in V. eizanensis, no population was identified as making a particularly higher contribution to H (T). These findings suggest that V. eizanensis had been isolated in two large glacial refugia, whereas populations of V. chaerophylloides var. sieboldiana were restricted to a single small refuge. Different light requirements between these two closely related species probably caused these differing responses to climatic change during the last ice age.

Do organellar genomes function as long-term redox damage sensors?

A small group of proteins that form core components of electron transfer complexes are consistently encoded by organellar genomes in multicellular organisms, suggesting functional constraint. These genomes are costly to maintain and vulnerable to mutation. We propose that they provide cell lineages with sensors of long-term redox damage, and of bioenergetic and genomic competence. This proposed adaptive function sets tonic retrograde signalling to the nucleus and anterograde responses influencing protective and cell death pathways. The nature of the proposed gain-of-function signalling mechanisms is unclear but could involve defective complex assembly. Organellar proteomes therefore provide cumulative feedback on bioenergetic and genomic status within cell lineages, selection of the energetically 'fittest' cells and a means of removing cells that compromise survival of the organism.

Doubly uniparental inheritance (DUI) of mitochondrial DNA in Donax trunculus (Bivalvia: Donacidae) and the problem of its sporadic detection in Bivalvia

Mitochondrial DNA is transmitted maternally in metazoan species. This rule does not hold in several species of bivalves that have two mtDNA types, one that is transmitted maternally and the other paternally. This system of mitochondrial DNA transmission is known as doubly uniparental inheritance (DUI). Here we present evidence of DUI in the clam Donax trunculus making Donacidae the sixth bivalve family in which the phenomenon has been found. In addition, we present the taxonomic affiliation of all species in which DUI is currently known to occur and construct a phylogeny of the maternal and paternal genomes of these species. We use this information to address the question of a single or multiple origins of DUI and to discuss whether failed attempts to demonstrate the presence of DUI in several bivalve species might be due to problems of detection or to genuine absence of the phenomenon.

Traces of ancient range shifts in a mountain plant group (Androsace halleri complex, Primulaceae)

Phylogeographical studies frequently detect range shifts, both expansions (including long-distance dispersal) and contractions (including vicariance), in the studied taxa. These processes are usually inferred from the patterns and distribution of genetic variation, with the potential pitfall that different historical processes may result in similar genetic patterns. Using a combination of DNA sequence data from the plastid genome, AFLP fingerprinting, and rigorous phylogenetic and coalescence-based hypothesis testing, we show that Androsace halleri (currently distributed disjunctly in the northwestern Iberian Cordillera Cantábrica, the eastern Pyrenees, and the French Massif Central and Vosges), or its ancestor, was once more widely distributed in the Pyrenees. While there, it hybridized with Androsace laggeri and Androsace pyrenaica, both of which are currently allopatric with A. halleri. The common ancestor of A. halleri and the north Iberian local endemic Androsace rioxana probably existed in the north Iberian mountain ranges with subsequent range expansion (to the French mountain ranges of the Massif Central and the Vosges) and allopatric speciation (A. rioxana, A. halleri in the eastern Pyrenees, A. halleri elsewhere). We have thus been able to use the reticulate evolution in this species group to help elucidate its phylogeographical history, including evidence of range contraction.

Concordant evolution of mitochondrial and nuclear genomes in the wheat pathogen Phaeosphaeria nodorum

We compared patterns of mitochondrial restriction fragment length polymorphism (RFLP) diversity with patterns of nuclear RFLP diversity to investigate the effects of selection, gene flow, and sexual reproduction on the population genetic structure and evolutionary history of the wheat pathogen Phaeosphaeria nodorum. A total of 315 fungal isolates from Texas, Oregon, and Switzerland were analyzed using seven nuclear RFLP probes that hybridized to discrete loci and purified mitochondrial DNA that hybridized to the entire mtDNA genome. Forty-two different mitochondrial haplotypes and 298 different nuclear haplotypes were detected. The two most frequent mtDNA haplotypes were present in every population and represented 32% of all isolates. High levels of gene flow, low levels of population subdivision, no evidence for either host specificity or cyto-nuclear disequilibrium were inferred from the analysis of both genomes. The concordance in estimates of these population genetic parameters from both genomes suggests that the two genomes experienced similar degrees of migration, genetic drift and selection.

Inherited maternal effects on the drought tolerance of a natural hybrid aquatic plant, Potamogeton anguillanus

We tested whether maternal effects have led to the adaptive divergence of strains of the natural hybrid Potamogeton anguillanus, whose putative parents show contrastingly divergent ecologies. To examine the correlation between phenotypic characters and maternal types, we conducted drought experiments and DNA typing using nuclear and chloroplast genes. In the field, we investigated the distribution of the maternal type along the depth and the inshore-offshore gradient. Hybrids of P. malaianus mothers (M-hybrids) and those of P. perfoliatus mothers (P-hybrids) could not be distinguished morphologically under submerged conditions, but differed in drought tolerance. M-hybrids and P. malaianus formed more terrestrial shoots and exhibited higher survival than P-hybrids and P. perfoliatus in drought experiments. The distribution survey clarified that M-hybrids were dominant in shallow and inshore areas, whereas they were almost absent in deeper and offshore areas. These results indicate that the natural hybrid P. anguillanus differs in adaptive values depending on the maternal type. Bidirectional hybridization and heritable maternal effects may have played important roles in its phenotypic adaptation to local environmental conditions.

A genomic perspective on the shortcomings of mitochondrial DNA for "barcoding" identification

Approximately 600-bp sequences of mitochondrial DNA (mtDNA) have been designated as "DNA barcodes" and have become one of the most contentious and animated issues in the application of genetic information to global biodiversity assessment and species identification. Advocates of DNA barcodes have received extensive attention and promotion in many popular and refereed scientific publications. However, we suggest that the utility of barcodes is suspect and vulnerable to technical challenges that are particularly pertinent to mtDNA. We review the natural history of mtDNA and discuss problems for barcoding which are particularly associated with mtDNA and inheritance, including reduced effective population size, maternal inheritance, recombination, inconsistent mutation rate, heteroplasmy, and compounding evolutionary processes. The aforementioned could significantly limit the application and utility of mtDNA barcoding efforts. Furthermore, global use of barcodes will require application and acceptance of a barcode-based species concept that has not been evaluated in the context of the extensive literature concerning species designation. Implementation of mtDNA barcodes in spite of technical and practical shortcomings we discuss may degrade the longstanding synthesis of genetic and organism-based research and will not advance studies ranging from genomic evolution to biodiversity assessment.

Longevity and the evolution of the mitochondrial DNA-coded proteins in mammals

The amino acids sequences of the mitochondrial DNA-coded peptides of placental mammals evolved at different rates in different branches of the mammalian phylogenetic tree. Adaptive selection was suggested to account for the faster evolution of some mitochondrial DNA-coded proteins in several branches of the mammalian tree, but the driving force(s) for the accelerated evolution has not been elucidated. Mitochondria generate reactive oxygen species (ROS) that appear to constrain the life span of many species. Therefore, I tested the hypothesis that the evolution of mammalian longevity drives the accelerated evolution of mitochondrial DNA-coded peptides. Using rodents as an outgroup for a clad that included most placental mammals (excluding rodents and hedgehogs) the computed rates of amino acid substitution per site were positively correlated with genus longevity (maximal observed averaged life span) for most of the mitochondrial DNA-coded peptides. The substitution per site of ATP6, the proton conducting subunit of ATPsynthase, CYTB, the core subunit of ubiquinone oxidoreductase that participate in both electron and proton transport, and ND3, a subunit of NADH dehydrogenase, showed the strongest correlations with longevity. Additional confirmation for the hypothesis was obtained by the observation that the genetic distances between placental mammals species that belong to different orders are positively correlated with the sum of longevities of the species pairs. The substitutions per site for the entire amino acid sequence coded by the heavy strand mtDNA were also positively correlated with the average longevities of the placental mammals orders. These results support the hypothesis that the evolution of longevity in mammals drove the accelerated evolution of mtDNA-coded peptide. It is suggested that, in mammals, adaptive selection of mutations that decrease the rate of production of reactive oxygen species, directly or indirectly (e.g. by increasing proton leak), increases longevity.

Variable populations within variable populations: quantifying mitochondrial heteroplasmy in natural populations of the gynodioecious plant Silene vulgaris

Populations of mitochondria reside within individuals. Among angiosperms, these populations are rarely considered as genetically variable entities and typically are not found to be heteroplasmic in nature, leading to the widespread assumption that plant mitochondrial populations are homoplasmic. However, empirical studies of mitochondrial variation in angiosperms are relatively uncommon due to a paucity of sequence variation. Recent greenhouse studies of Silene vulgaris suggested that heteroplasmy might occur in this species at a level that it is biologically relevant. Here, we use established qualitative methods and a novel quantitative PCR method to study the intraindividual population genetics of mitochondria across two generations in natural populations of S. vulgaris. We show incidences of heteroplasmy for mitochondrial atpA and patterns of inheritance that are suggestive of more widespread heteroplasmy at both atpA and cox1. Further, our results demonstrate that quantitative levels of mitochondrial variation within individuals are high, constituting 26% of the total in one population. These findings are most consistent with a biparental model of mitochondrial inheritance. However, selection within individuals may be instrumental in the maintenance of variation because S. vulgaris is gynodioecious. Male sterility is, in part, regulated by the mitochondrial genome, and strong selection pressures appear to influence the frequency of females in these populations.

Extensive gene flow blurs phylogeographic but not phylogenetic signal in Olea europaea L

Genetic structure and evolutionary patterns of the wild olive tree (Olea europaea L.) were investigated with AFLP fingerprinting data at three geographic levels: (a) phylogenetic relationships of the six currently recognized subspecies in Eurasia and Africa; (b) lineage identification in subsp. europaea of the Mediterranean basin; and (c) phylogeography in the western Mediterranean. Two statistical approaches (Bayesian inference and analysis of molecular variance) were used to analyse the AFLP fingerprints. To determine the congruency and transferability of results across studies previous RAPD and ISSR data were analysed in a similar manner. Comparisons proved that qualitative results were mostly congruent but quantitative values differed, depending on the method of analysis. Neighbour-Joining analysis of AFLP phenotypes supported current classification of subspecies. At a Mediterranean scale no clear cut phylogeographic pattern was recovered, likely due to extensive gene flow between populations of subsp. europaea. Gene flow estimates calculated with conventional F-statistics showed that reproductive barriers separated neither populations nor lineages of O. europaea. Genetic divergence between eastern and western parts of the Mediterranean basin was observed only when geographical and population information were incorporated into the analyses through hierarchical analysis of molecular variance (AMOVA). Within the western Mediterranean, the highest genetic diversity was found in two regions: on both sides of the Strait of Gibraltar and in the Balearic archipelago. Additionally, long-lasting isolation of the northern-most populations of the Iberian Peninsula appeared to be responsible for a significant divergence.

Studies on Cultivated Ephedra Plants in Inner Mongolia Autonomous Region and Ningxia Hui Autonomous Region

Progression of the desertification in northern China has been causing damage to wild Ephedra plants on which we depend for most of supply of the traditional herbal medicine, "Ma huang." The Chinese government encourages the cultivation of Ephedra plants, and Ephedra fields have been reclaimed in the original Ephedra habitats in recent years. We surveyed 7 Ephedra fields that have been recently developed in the Inner Mongolia Autonomous Region and Ningxia Hui Autonomous Region to collect information on Ephedra plant cultivation, especially pertaining to crop species. Specimens taken from those Ephedra fields were genetically and morphologically analyzed, and their ephedrine alkaloid content was examined. DNA analyses of Ephedra specimens, including DNA sequencing of ITS (internal transcribing sequence of nuclear ribosomal DNA) and trn L/F (intron of trnL and intergenic spacer between the trnL and trnF of chloroplast DNA) region and species-specific amplification of trn L/F were conducted to identify Ephedra species. Based on the results of DNA sequencing and morphological determination, the crops grown in 6 fields ware identified as Ephedra sinica, while co-planting of E. sinica and E. intermedia was found in one field where a higher appearance rate of plants with varied morphology from wild Ephedra plants was observed. Furthermore, direct sequencing of the PCR product of the trn L/F region of some specimens from the field and their species-specific PCR showed ambivalent result. Cloning and sequencing of the PCR product of the trn L/F region of those specimens DNA suggested their heteroplasmy, containing both E. sinica- and E. intermedia-type chloroplasts. On the other hand, the profile of the ephedrine alkaloid content was clearly correlated with the result of direct sequencing of the trn L/F region; the specimens showing the E. sinica-type sequence contained more ephedrine than pseudoephedrine, and the specimens of the E. intermedia-type more pseudoephedrine.