The complete chloroplast genome sequence of Mikania micrantha (Asteraceae), a noxious invasive weed to South China.

It is the first report on complete chloroplast genome of Mikania micrantha (Asteraceae), a noxious invasive weed to South China. The genome is a circular molecule of 152,092 bp in length with 37.58% average GC content, and includes a large single copy region (83,793 bp), a small single copy region (18,287 bp), and two inverted repeat regions (25,006 bp). The M. micrantha cp genome encodes 135 unique genes, including 90 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. ML tree based on 16 complete cpDNA from Asteraceae indicated that M. micrantha has a close sister relationship with Ageratina adenophora and Praxelis clematidea. The complete cpDNA of M. micrantha provides useful molecular data for further phylogenetic and evolutionary analysis.

Plastome sequences of Lygodium japonicum and Marsilea crenata reveal the genome organization transformation from basal ferns to core leptosporangiates.

Previous studies have shown that core leptosporangiates, the most species-rich group of extant ferns (monilophytes), have a distinct plastid genome (plastome) organization pattern from basal fern lineages. However, the details of genome structure transformation from ancestral ferns to core leptosporangiates remain unclear because of limited plastome data available. Here, we have determined the complete chloroplast genome sequences of Lygodium japonicum (Lygodiaceae), a member of schizaeoid ferns (Schizaeales), and Marsilea crenata (Marsileaceae), a representative of heterosporous ferns (Salviniales). The two species represent the sister and the basal lineages of core leptosporangiates, respectively, for which the plastome sequences are currently unavailable. Comparative genomic analysis of all sequenced fern plastomes reveals that the gene order of L. japonicum plastome occupies an intermediate position between that of basal ferns and core leptosporangiates. The two exons of the fern ndhB gene have a unique pattern of intragenic copy number variances. Specifically, the substitution rate heterogeneity between the two exons is congruent with their copy number changes, confirming the constraint role that inverted repeats may play on the substitution rate of chloroplast gene sequences.

The complete chloroplast genome sequence of Cephalotaxus oliveri (Cephalotaxaceae): evolutionary comparison of cephalotaxus chloroplast DNAs and insights into the loss of inverted repeat copies in gymnosperms.

We have determined the complete chloroplast (cp) genome sequence of Cephalotaxus oliveri. The genome is 134,337 bp in length, encodes 113 genes, and lacks inverted repeat (IR) regions. Genome-wide mutational dynamics have been investigated through comparative analysis of the cp genomes of C. oliveri and C. wilsoniana. Gene order transformation analyses indicate that when distinct isomers are considered as alternative structures for the ancestral cp genome of cupressophyte and Pinaceae lineages, it is not possible to distinguish between hypotheses favoring retention of the same IR region in cupressophyte and Pinaceae cp genomes from a hypothesis proposing independent loss of IRA and IRB. Furthermore, in cupressophyte cp genomes, the highly reduced IRs are replaced by short repeats that have the potential to mediate homologous recombination, analogous to the situation in Pinaceae. The importance of repeats in the mutational dynamics of cupressophyte cp genomes is also illustrated by the accD reading frame, which has undergone extreme length expansion in cupressophytes. This has been caused by a large insertion comprising multiple repeat sequences. Overall, we find that the distribution of repeats, indels, and substitutions is significantly correlated in Cephalotaxus cp genomes, consistent with a hypothesis that repeats play a role in inducing substitutions and indels in conifer cp genomes.

Microsatellite loci for an old rare species, Pseudotaxus chienii, and transferability in Taxus wallichiana var. mairei (Taxaceae).

PREMISE OF THE STUDY: Microsatellite loci were developed for Pseudotaxus chienii, an old rare species endemic to China, and which provided a useful tool for investigating the patterns of population genetic structure, phylogeography, evolutionary history, and adaptive potential. Transferability was assayed in the related species, Taxus wallichiana var. mairei. • METHODS AND RESULTS: A total of 15 microsatellite loci were targeted in P. chienii using the Fast Isolation by AFLP of Sequences COntaining Repeats (FIASCO) protocol. Polymorphism was evaluated in five populations of P. chienii and five populations of T. wallichiana var. mairei. Of these loci, 13 were polymorphic in P. chienii, whereas 15 were polymorphic in T. wallichiana var. mairei. • CONCLUSIONS: The 15 microsatellite loci developed lay a solid foundation for further studies on population genetic variability and investigations of local adaptation. Additionally, cross-species amplification in T. wallichiana var. mairei showed that these loci may also have potential utility in other genera of Taxaceae.

Molecular evolution of psbA gene in ferns: unraveling selective pressure and co-evolutionary pattern.

BACKGROUND: The photosynthetic oxygen-evolving photo system II (PS II) produces almost the entire oxygen in the atmosphere. This unique biochemical system comprises a functional core complex that is encoded by psbA and other genes. Unraveling the evolutionary dynamics of this gene is of particular interest owing to its direct role in oxygen production. psbA underwent gene duplication in leptosporangiates, in which both copies have been preserved since. Because gene duplication is often followed by the non-fictionalization of one of the copies and its subsequent erosion, preservation of both psbA copies pinpoint functional or regulatory specialization events. The aim of this study was to investigate the molecular evolution of psbA among fern lineages. RESULTS: We sequenced psbA , which encodes D1 protein in the core complex of PSII, in 20 species representing 8 orders of extant ferns; then we searched for selection and convolution signatures in psbA across the 11 fern orders. Collectively, our results indicate that: (1) selective constraints among D1 protein relaxed after the duplication in 4 leptosporangiate orders; (2) a handful positively selected codons were detected within species of single copy psbA, but none in duplicated ones; (3) a few sites among D1 protein were involved in co-evolution process which may intimate significant functional/structural communications between them. CONCLUSIONS: The strong competition between ferns and angiosperms for light may have been the main cause for a continuous fixation of adaptive amino acid changes in psbA , in particular after its duplication. Alternatively, a single psbA copy may have undergone bursts of adaptive changes at the molecular level to overcome angiosperms competition. The strong signature of positive Darwinian selection in a major part of D1 protein is testament to this. At the same time, species own two psbA copies hardly have positive selection signals among the D1 protein coding sequences. In this study, eleven co-evolving sites have been detected via different molecules, which may be more important than others.

AFLP genome scan to detect genetic structure and candidate loci under selection for local adaptation of the invasive weed Mikania micrantha.

Why some species become successful invaders is an important issue in invasive biology. However, limited genomic resources make it very difficult for identifying candidate genes involved in invasiveness. Mikania micrantha H.B.K. (Asteraceae), one of the world's most invasive weeds, has adapted rapidly in response to novel environments since its introduction to southern China. In its genome, we expect to find outlier loci under selection for local adaptation, critical to dissecting the molecular mechanisms of invasiveness. An explorative amplified fragment length polymorphism (AFLP) genome scan was used to detect candidate loci under selection in 28 M. micrantha populations across its entire introduced range in southern China. We also estimated population genetic parameters, bottleneck signatures, and linkage disequilibrium. In binary characters, such as presence or absence of AFLP bands, if all four character combinations are present, it is referred to as a character incompatibility. Since character incompatibility is deemed to be rare in populations with extensive asexual reproduction, a character incompatibility analysis was also performed in order to infer the predominant mating system in the introduced M. micrantha populations. Out of 483 AFLP loci examined using stringent significance criteria, 14 highly credible outlier loci were identified by Dfdist and Bayescan. Moreover, remarkable genetic variation, multiple introductions, substantial bottlenecks and character compatibility were found to occur in M. micrantha. Thus local adaptation at the genome level indeed exists in M. micrantha, and may represent a major evolutionary mechanism of successful invasion. Interactions between genetic diversity, multiple introductions, and reproductive modes contribute to increase the capacity of adaptive evolution.

Development of microsatellite loci for Cephalotaxus oliveri (Cephalotaxaceae) and cross-amplification in Cephalotaxus.

PREMISE OF THE STUDY: Microsatellite loci were developed for Cephalotaxus oliveri, an endemic and endangered conifer in China, which will allow assessment of the levels of genetic diversity and a means to understand the genetic consequences of habitat fragmentation. METHODS AND RESULTS: Using the Fast Isolation by AFLP of Sequences COntaining (FIASCO) Repeats protocol, 19 microsatellite loci were identified in C. oliveri, 13 of which were polymorphic within a sample of 52 individuals representing five natural populations. The actual number of alleles per locus ranged from one to five. Twelve polymorphic loci were also successfully amplified in C. fortunei. CONCLUSIONS: These microsatellite loci will provide a useful tool for further investigating genetic variation in natural populations of C. oliveri, which will inform future conservation and management strategies. Additionally, cross-amplification in C. fortunei suggested the potential utility of these loci in this and other congeneric species.

Molecular evolution of rbcL in three gymnosperm families: identifying adaptive and coevolutionary patterns.

BACKGROUND: The chloroplast-localized ribulose-1, 5-biphosphate carboxylase/oxygenase (Rubisco), the primary enzyme responsible for autotrophy, is instrumental in the continual adaptation of plants to variations in the concentrations of CO2. The large subunit (LSU) of Rubisco is encoded by the chloroplast rbcL gene. Although adaptive processes have been previously identified at this gene, characterizing the relationships between the mutational dynamics at the protein level may yield clues on the biological meaning of such adaptive processes. The role of such coevolutionary dynamics in the continual fine-tuning of RbcL remains obscure. RESULTS: We used the timescale and phylogenetic analyses to investigate and search for processes of adaptive evolution in rbcL gene in three gymnosperm families, namely Podocarpaceae, Taxaceae and Cephalotaxaceae. To understand the relationships between regions identified as having evolved under adaptive evolution, we performed coevolutionary analyses using the software CAPS. Importantly, adaptive processes were identified at amino acid sites located on the contact regions among the Rubisco subunits and on the interface between Rubisco and its activase. Adaptive amino acid replacements at these regions may have optimized the holoenzyme activity. This hypothesis was pinpointed by evidence originated from our analysis of coevolution that supported the correlated evolution between Rubisco and its activase. Interestingly, the correlated adaptive processes between both these proteins have paralleled the geological variation history of the concentration of atmospheric CO2. CONCLUSIONS: The gene rbcL has experienced bursts of adaptations in response to the changing concentration of CO2 in the atmosphere. These adaptations have emerged as a result of a continuous dynamic of mutations, many of which may have involved innovation of functional Rubisco features. Analysis of the protein structure and the functional implications of such mutations put forward the conclusion that this evolutionary scenario has been possible through a complex interplay between adaptive mutations, often structurally destabilizing, and compensatory mutations. Our results unearth patterns of evolution that have likely optimized the Rubisco activity and uncover mutational dynamics useful in the molecular engineering of enzymatic activities. REVIEWERS: This article was reviewed by Prof. Christian Blouin (nominated by Dr W Ford Doolittle), Dr Endre Barta (nominated by Dr Sandor Pongor), and Dr Nicolas Galtier.

Evolution of the rpoB-psbZ region in fern plastid genomes: notable structural rearrangements and highly variable intergenic spacers.

BACKGROUND: The rpoB-psbZ (BZ) region of some fern plastid genomes (plastomes) has been noted to go through considerable genomic changes. Unraveling its evolutionary dynamics across all fern lineages will lead to clarify the fundamental process shaping fern plastome structure and organization. RESULTS: A total of 24 fern BZ sequences were investigated with taxon sampling covering all the extant fern orders. We found that: (i) a tree fern Plagiogyria japonica contained a novel gene order that can be generated from either the ancestral Angiopteris type or the derived Adiantum type via a single inversion; (ii) the trnY-trnE intergenic spacer (IGS) of the filmy fern Vandenboschia radicans was expanded 3-fold due to the tandem 27-bp repeats which showed strong sequence similarity with the anticodon domain of trnY; (iii) the trnY-trnE IGSs of two horsetail ferns Equisetum ramosissimum and E. arvense underwent an unprecedented 5-kb long expansion, more than a quarter of which was consisted of a single type of direct repeats also relevant to the trnY anticodon domain; and (iv) ycf66 has independently lost at least four times in ferns. CONCLUSIONS: Our results provided fresh insights into the evolutionary process of fern BZ regions. The intermediate BZ gene order was not detected, supporting that the Adiantum type was generated by two inversions occurring in pairs. The occurrence of Vandenboschia 27-bp repeats represents the first evidence of partial tRNA gene duplication in fern plastomes. Repeats potentially forming a stem-loop structure play major roles in the expansion of the trnY-trnE IGS.

[Phylogenetic analysis and detection of selective pressure on the CVNH domains].

Cyanovirin-N (CV-N) is a novel protein with broad-spectrum antiviral activity. Its homologs constitute a protein family known as CVNH (Cyanovirin-N homology), which possess the evolutionarily conserved anti-HIV (Human immunodeficiency virus) domain. In this study, more details about the patchy organism distribution of CVNH domain were explored by reconstructing gene trees. Duplicated CVNH sequences were also identified in a wide range of species including Aspergillus niger, Neosartorya fischeri NRRL 181, Penicillium chrysogenum Wisconsin 54-1255, Neurospora crassa, Cyanothece sp. PCC, and Ceratopteris richardii. Besides these findings, both the mechanistic and mechanistic-empirical combination (MEC) models were used to analyze the adaptive evolution of amino acid sites in the CVNH domain. Our results showed that: (1) neither model reveals significant sites undergoing positive selection; (2) purifying selection has played a dominant role during CVNH evolution; and (3) the MEC model better fits the CVNH data set. Furthermore, the ancestral branch leading to Cyanothece sp. PCC 7822 and 7424 were examined using "branch-specific" and "branch-site" models. Six positively selected sites (34L, 63L, 13H, 76C, 78K, and 80I) were identified on the branch.

Complete chloroplast genome sequence of a tree fern Alsophila spinulosa: insights into evolutionary changes in fern chloroplast genomes.

BACKGROUND: Ferns have generally been neglected in studies of chloroplast genomics. Before this study, only one polypod and two basal ferns had their complete chloroplast (cp) genome reported. Tree ferns represent an ancient fern lineage that first occurred in the Late Triassic. In recent phylogenetic analyses, tree ferns were shown to be the sister group of polypods, the most diverse group of living ferns. Availability of cp genome sequence from a tree fern will facilitate interpretation of the evolutionary changes of fern cp genomes. Here we have sequenced the complete cp genome of a scaly tree fern Alsophila spinulosa (Cyatheaceae). RESULTS: The Alsophila cp genome is 156,661 base pairs (bp) in size, and has a typical quadripartite structure with the large (LSC, 86,308 bp) and small single copy (SSC, 21,623 bp) regions separated by two copies of an inverted repeat (IRs, 24,365 bp each). This genome contains 117 different genes encoding 85 proteins, 4 rRNAs and 28 tRNAs. Pseudogenes of ycf66 and trnT-UGU are also detected in this genome. A unique trnR-UCG gene (derived from trnR-CCG) is found between rbcL and accD. The Alsophila cp genome shares some unusual characteristics with the previously sequenced cp genome of the polypod fern Adiantum capillus-veneris, including the absence of 5 tRNA genes that exist in most other cp genomes. The genome shows a high degree of synteny with that of Adiantum, but differs considerably from two basal ferns (Angiopteris evecta and Psilotum nudum). At one endpoint of an ancient inversion we detected a highly repeated 565-bp-region that is absent from the Adiantum cp genome. An additional minor inversion of the trnD-GUC, which is possibly shared by all ferns, was identified by comparison between the fern and other land plant cp genomes. CONCLUSION: By comparing four fern cp genome sequences it was confirmed that two major rearrangements distinguish higher leptosporangiate ferns from basal fern lineages. The Alsophila cp genome is very similar to that of the polypod fern Adiantum in terms of gene content, gene order and GC content. However, there exist some striking differences between them: the trnR-UCG gene represents a putative molecular apomorphy of tree ferns; and the repeats observed at one inversion endpoint may be a vestige of some unknown rearrangement(s). This work provided fresh insights into the fern cp genome evolution as well as useful data for future phylogenetic studies.

Genetic differentiation of relictual populations of Alsophila spinulosa in southern China inferred from cpDNA trnL-F noncoding sequences.

The genetic differentiation and phylogeographical pattern of 11 relictual populations of Alsophila spinulosa distributed across Hainan, Guangdong, and Guangxi in southern China were inferred from sequence variations of trnL-F noncoding regions of chloroplast DNA (cpDNA). The length of trnL-F noncoding sequences varied from 863 to 940 bp. The A + T content was 62.23-63.36%. Sequences were neutral in terms of evolution (Tajima's criterion D=-0.62417, P>0.10 and Fu and Li's test D*=-1.45455, P>0.10; F*=-1.32798, P>0.10). Thirty-four haplotypes were identified based on nucleotide variation. Relatively high levels of haplotype diversity (h=0.929) and nucleotide diversity (Dij=0.022263) were detected in A. spinulosa, probably associated with its long evolutionary history which allowed the accumulation of genetic variation within lineages. Both the minimum spanning network and the strict consensus tree of the most parsimonious trees generated for haplotypes demonstrated that the investigated populations of A. spinulosa were subdivided into two geographical groups: Hainan and Guangdong-Guangxi. An analysis of molecular variance (AMOVA) indicated that most of the genetic variation (87.48%, P<0.001) was partitioned among regions. Spatial structure measurements revealed that population genetic structure was not related to geographical distance. This research suggests that blocked gene flow by Qiongzhou strait and an inbreeding system might result in the geographical subdivision between Hainan and Guangdong-Guangxi (F(ST)=0.92, Nm=0.09). Within each region, the "star like" pattern of phylogeography of haplotypes implied a population expansion process during evolutionary history. Gene genealogies together with coalescent theory were useful tools for uncovering the phylogeography of A. spinulosa.

Genetic structure and variation in the relict populations of Alsophila spinulosa from southern China based on RAPD markers and cpDNA atpB-rbcL sequence data.

RAPD markers and sequences of chloroplast DNA (cpDNA) atpB-rbcL intergenic spacers were used to characterize the pattern of genetic variation and the phylogenetic relationships of the relict populations of Alsophila spinulosa located in Jian Feng Ling (JFL) and Diao Luo Shan (DLS), Hainan, and Tang Lang Shan (TLS), Ding Hu Shan (DHS), and Da Xi Shan (DXS), Guangdong, of southern China. 28 random primers generated 118 bands, out of which 26 (22.03%) were polymorphic loci, distinguishing 17 different RAPD phenotypes. Percentage of polymorphic loci, Shannon phenotypic diversity and Nei's gene diversity comprehensively indicated that JFL possessed the highest diversity, TLS and DHS in intermediate and DLS or DXS the least; the corresponding values of the population appeared correlated with the population size. Differentiation was detected among populations of A. spinulosa (1-Hpop/Hsp=0.7453, GST=0.7763, and phist=0.8145). AMOVA showed that 47.44% of the variance was partitioned among regions (Hainan and Guangdong), 34.01% attributed among populations within regions, whereas only 18.55% occurring within populations. Low level of intra-specific diversity was maintained in A. spinulosa with Shannon diversity and gene diversity merely 0.0560 and 0.0590, repectively. Sequence length of atpB-rbcL intergenic spacer varied from 724 bp to 730 bp. Base composition was with A+T content between 63.17% and 63.70%. 13 haplotypes of atpB-rbcL noncoding spacers were identified. UPGMA dendrogram of RAPD phenotypes, principal components analysis based on RAPD patterns, minimum spanning network and neighbour-joining (NJ) tree established on atpB-rbcL haplotypes consistently suggested the geographical subdivision of populations of A. spinulosa between Hainan and Guangdong. Breeding system and conservation strategy of A. spinulosa was discussed based on the information of population genetic structure and variation.

[A study of the wild and captive green peafowl(Pavo muticus) by random-amplified polymorphic DNA].

Random-amplified polymorphic DNA (RAPD) was used to investigate the genetic diversity of the population of 14 wild green peafowl and 18 captive green peafowl (Pavo muticus). Total of 161 and 166 bands were obtained respectively, and 23 random primers were used to amplify the genomic DNA of the wild and captive green peafowls. The average relative hereditary distance of the wild and captive green peafowls is 0.0555 and 0.1355 respectively; and the Shannon diversity index is 0.4348 and 1.0163 respectively. There is a prominent difference between the two populations by T-Test of HO. All the analyses above show that the genetic diversity is very low in wild green peafowl. It tells us that the two populations come from two families by using UPGMA, which can be useful in the breeding management in the future.