Chloroplast Genome Sequences and Phylogenetic Analysis of Eight Newly Sequenced Caryophyllaceae Species.

BACKGROUND: Caryophyllaceae is a big family composed of many economic and medicinal species. However, the phylogeny of the family is insufficient and genome data are lacking for many species. OBJECTIVE: Using next-generation sequencing (NGS) to acquire the chloroplast (cp) genomes of Eremogone acicularis (F.N.Williams) Ikonn., E. brevipetala (Tsui & L.H.Zhou) Sadeghian & Zarre, E. bryophylla (Fernald) Pusalkar & D.K.Singh, E. kansuensis (Maxim.) Dillenb. & Kadereit, Shivparvatia glanduligera (Edgew.) Pusalkar & D.K.Singh, Silene atsaensis (Marq.) Bocquet, S. caespitella Williams, and S. lhassana (Williams) Majumdar. METHODS: Bioinformatic software was used to conduct the comparative genome and phylogeny analysis of these cp genomes. RESULTS: The eight cp genomes were 132 188-151 919 bp in length, containing 130-132 genes. A/T was dominant in simple sequence repeats (SSRs). Forward repeats and palindromic repeats were the most frequent in long terminal repeats (LTRs). Compared with the four species of Eremogone Fenzl, the inverted repeat (IR) boundaries of S. caespitella, S. atsaensis, S. lhassana, and Sh. glanduligera were significantly expanded. Four and one mutational hotspots were identified in the large single copy (LSC) region and small single copy (SSC) region, respectively. The ratio of nonsynonymous substitution to synonymous substitution (Ka/Ks ratio) showed these cp genomes may have undergone strong purifying selection. In the phylogenetic trees, both Silene L. and Eremogone were monophyletic groups. However, Sh. glanduligera was closely related to Amaranthus hypochondriacus. CONCLUSION: These results have provided new evidence and useful information for species identification, evolution, and genetic research on the Caryophyllaceae. HIGHLIGHTS: In this study, eight newly sequenced cp genomes of Caryophyllaceae species were reported for the first time.

Complete chloroplast genomes of Cerastium alpinum, C. arcticum and C. nigrescens: genome structures, comparative and phylogenetic analysis.

The genus Cerastium includes about 200 species that are mostly found in the temperate climates of the Northern Hemisphere. Here we report the complete chloroplast genomes of Cerastium alpinum, C. arcticum and C. nigrescens. The length of cp genomes ranged from 147,940 to 148,722 bp. Their quadripartite circular structure had the same gene organization and content, containing 79 protein-coding genes, 30 tRNA genes, and four rRNA genes. Repeat sequences varied from 16 to 23 per species, with palindromic repeats being the most frequent. The number of identified SSRs ranged from 20 to 23 per species and they were mainly composed of mononucleotide repeats containing A/T units. Based on Ka/Ks ratio values, most genes were subjected to purifying selection. The newly sequenced chloroplast genomes were characterized by a high frequency of RNA editing, including both C to U and U to C conversion. The phylogenetic relationships within the genus Cerastium and family Caryophyllaceae were reconstructed based on the sequences of 71 protein-coding genes. The topology of the phylogenetic tree was consistent with the systematic position of the studied species. All representatives of the genus Cerastium were gathered in a single clade with C. glomeratum sharing the least similarity with the others.

Evolution of the Sabulina verna group (Caryophyllaceae) in Europe: A deep split, followed by secondary contacts, multiple allopolyploidization and colonization of challenging substrates.

One of the major goals of contemporary evolutionary biology is to elucidate the relative roles of allopatric and ecological differentiation and polyploidy in speciation. In this study, we address the taxonomically intricate Sabulina verna group, which has a disjunct Arctic-alpine postglacial range in Europe and occupies a broad range of ecological niches, including substrates toxic to plants. Using genome-wide ddRAD sequencing combined with morphometric analyses based on extensive sampling of 111 natural populations, we aimed to disentangle internal evolutionary relationships and examine their correspondence with the pronounced edaphic and ploidy diversity within the group. We identified two spatially distinct groups of diploids: a widespread Arctic-alpine group and a spatially restricted yet diverse Balkan group. Most tetraploids exhibited a considerably admixed ancestry derived from both these groups, suggesting their allopolyploid origin. Four genetic clusters in congruence with geography and mostly supported by morphological traits were recognized in the diploid Arctic-alpine group. Tetraploids are split into two distinct and geographically vicariant groups, indicating their repeated polytopic origin. Furthermore, our results also revealed at least five-fold parallel colonization of toxic substrates (serpentine and metalliferous), altogether demonstrating a complex interaction between geography, challenging substrates and polyploidy in the evolution of the group. Finally, we propose a new taxonomic treatment of this complex.

[Cloning and functional verification of PhAEP gene, a key enzyme for biosynthesis of heterophyllin A in Pseudostellaria heterophylla].

This paper aimed to study the role of asparagine endopeptidase(AEP) gene in the biosynthesis mechanism of cyclic peptide compounds in Pseudostellaria heterophylla. The transcriptome database of P. heterophylla was systematically mined and screened, and an AEP gene, tentatively named PhAEP, was successfully cloned. The heterologous function verification by Nicotiana benthamiana showed that the expression of the gene played a role in the biosynthesis of heterophyllin A in P. heterophylla. Bioinformatics analysis showed that the cDNA of PhAEP was 1 488 bp in length, encoding 495 amino acids with a molecular weight of 54.72 kDa. The phylogenetic tree showed that the amino acid sequence encoded by PhAEP was highly similar to that of Butelase-1 in Clitoria ternatea, reaching 80%. The sequence homology and cyclase active site analysis revealed that the PhAEP enzyme may specifically hydrolyse the C-terminal Asn/Asp(Asx) site of the core peptide in the HA linear precursor peptide of P. heterophylla, thereby participating in the ring formation of the linear precursor peptide. The results of real-time quantitative polymerase chain reaction(RT-qPCR) showed that the expression level of PhAEP was the highest in fruits, followed by in roots, and the lowest in leaves. The heterophyllin A of P. heterophylla was detected in N. benthamiana that co-expressed PrePhHA and PhAEP genes instantaneously. In this study, the PhAEP gene, a key enzyme in the biosynthesis of heterophyllin A in P. heterophylla, has been successfully cloned, which lays a foundation for further analysis of the molecular mechanism of PhAEP enzyme in the biosynthesis of heterophyllin A in P. heterophylla and has important significance for the study of synthetic biology of cyclic peptide compounds in P. heterophylla.

Parallel adaptation to lower altitudes is associated with enhanced plasticity in Heliosperma pusillum (Caryophyllaceae).

High levels of phenotypic plasticity are thought to be inherently costly in stable or extreme environments, but enhanced plasticity may evolve as a response to new environments and foster novel phenotypes. Heliosperma pusillum forms glabrous alpine and pubescent montane ecotypes that diverged recurrently and polytopically (parallel evolution) and can serve as evolutionary replicates. The specific alpine and montane localities are characterized by distinct temperature conditions, available moisture, and light. Noteworthy, the ecotypes show a home-site fitness advantage in reciprocal transplantations. To disentangle the relative contribution of constitutive versus plastic gene expression to altitudinal divergence, we analyze the transcriptomic profiles of two parallely evolved ecotype pairs, grown in reciprocal transplantations at native altitudinal sites. In this incipient stage of divergence, only a minor proportion of genes appear constitutively differentially expressed between the ecotypes in both pairs, regardless of the growing environment. Both derived, montane populations bear comparatively higher plasticity of gene expression than the alpine populations. Genes that change expression plastically or constitutively underlie similar ecologically relevant pathways, related to response to drought and trichome formation. Other relevant processes, such as photosynthesis, rely mainly on plastic changes. The enhanced plasticity consistently observed in the montane ecotype likely evolved as a response to the newly colonized, drier, and warmer niche. We report a striking parallelism of directional changes in gene expression plasticity. Thus, plasticity appears to be a key mechanism shaping the initial stages of phenotypic evolution, likely fostering adaptation to novel environments.

Polygenic routes lead to parallel altitudinal adaptation in Heliosperma pusillum (Caryophyllaceae).

Understanding how organisms adapt to the environment is a major goal of modern biology. Parallel evolution-the independent evolution of similar phenotypes in different populations-provides a powerful framework to investigate the evolutionary potential of populations, the constraints of evolution, its repeatability and therefore its predictability. Here, we quantified the degree of gene expression and functional parallelism across replicated ecotype formation in Heliosperma pusillum (Caryophyllaceae), and gained insights into the architecture of adaptive traits. Population structure analyses and demographic modelling support a previously formulated hypothesis of parallel polytopic divergence of montane and alpine ecotypes. We detect a large proportion of differentially expressed genes (DEGs) underlying divergence within each replicate ecotype pair, with a strikingly low number of shared DEGs across pairs. Functional enrichment of DEGs reveals that the traits affected by significant expression divergence are largely consistent across ecotype pairs, in strong contrast to the nonshared genetic basis. The remarkable redundancy of differential gene expression indicates a polygenic architecture for the diverged adaptive traits. We conclude that polygenic traits appear key to opening multiple routes for adaptation, widening the adaptive potential of organisms.

Genetic variation of Cerastium alpinum L. from Babia Góra, a critically endangered species in Poland.

Babia Góra massif is the only site of occurrence of the Cerastium alpinum L. in Poland, an arctic-alpine perennial plant with a wide distribution in North America, northwestern Asia, and Europe. To determine whether the isolated Polish populations are genetically distinct, we have performed an evaluation of C. alpinum from Babia Góra with the use of iPBS markers. A total number of 133 individuals of C. alpinum from seven populations representing four localizations of the species were analyzed, i.e., from Babia Góra (Poland), Alps (Switzerland), Nuolja massif (Sweden), and Kaffiøyra (Svalbard, Norway). Genetic analysis of all C. alpinum samples using eight PBS primers identified 262 bands, 79.4% of which were polymorphic. iPBS markers revealed low genetic diversity (average He = 0.085) and high population differentiation (FST = 0.617). AMOVA results confirmed that the majority of the genetic variation (62%) was recorded among populations. The grouping revealed by PCoA showed that C. alpinum from Svalbard is the most diverged population, C. alpinum from Switzerland and Sweden form a pair of similar populations, whereas C. alpinum from the Babia Góra form a heterogeneous group of four populations. Results of isolation by distance analysis suggested that the spatial distance is the most probable cause of the observed differentiation among populations. Although significant traces of a bottleneck effect were noted for all populations of C. alpinum from Babia Góra, the populations still maintain a low but significant level of genetic polymorphism. These results are of great importance for developing conservation strategies for this species in Poland.

Cloning and functional verification of PhAEP gene, a key enzyme for biosynthesis of heterophyllin A in Pseudostellaria heterophylla / 中国中药杂志

This paper aimed to study the role of asparagine endopeptidase(AEP) gene in the biosynthesis mechanism of cyclic peptide compounds in Pseudostellaria heterophylla. The transcriptome database of P. heterophylla was systematically mined and screened, and an AEP gene, tentatively named PhAEP, was successfully cloned. The heterologous function verification by Nicotiana benthamiana showed that the expression of the gene played a role in the biosynthesis of heterophyllin A in P. heterophylla. Bioinformatics analysis showed that the cDNA of PhAEP was 1 488 bp in length, encoding 495 amino acids with a molecular weight of 54.72 kDa. The phylogenetic tree showed that the amino acid sequence encoded by PhAEP was highly similar to that of Butelase-1 in Clitoria ternatea, reaching 80%. The sequence homology and cyclase active site analysis revealed that the PhAEP enzyme may specifically hydrolyse the C-terminal Asn/Asp(Asx) site of the core peptide in the HA linear precursor peptide of P. heterophylla, thereby participating in the ring formation of the linear precursor peptide. The results of real-time quantitative polymerase chain reaction(RT-qPCR) showed that the expression level of PhAEP was the highest in fruits, followed by in roots, and the lowest in leaves. The heterophyllin A of P. heterophylla was detected in N. benthamiana that co-expressed PrePhHA and PhAEP genes instantaneously. In this study, the PhAEP gene, a key enzyme in the biosynthesis of heterophyllin A in P. heterophylla, has been successfully cloned, which lays a foundation for further analysis of the molecular mechanism of PhAEP enzyme in the biosynthesis of heterophyllin A in P. heterophylla and has important significance for the study of synthetic biology of cyclic peptide compounds in P. heterophylla.

Genetic potential for changes in breeding systems: Predicted and observed trait changes during artificial selection for male and female allocation in a gynodioecious species.

PREMISE: Evolution of separate sexes from hermaphroditism often proceeds through gynodioecy, but genetic constraints on this process are poorly understood. Genetic (co-)variances and between-sex genetic correlations were used to predict evolutionary responses of multiple reproductive traits in a sexually dimorphic gynodioecious species, and predictions were compared with observed responses to artificial selection. METHODS: Schiedea (Caryophyllaceae) is an endemic Hawaiian lineage with hermaphroditic, gynodioecious, subdioecious, and dioecious species. We measured genetic parameters of Schiedea salicaria and used them to predict evolutionary responses of 18 traits in hermaphrodites and females in response to artificial selection for increased male (stamen) biomass in hermaphrodites or increased female (carpel, capsule) biomass in females. Observed responses over two generations were compared with predictions in replicate lines of treatments and controls. RESULTS: In only two generations, both stamen biomass in hermaphrodites and female biomass in females responded markedly to direct selection, supporting a key assumption of models for evolution of dioecy. Other biomass traits, pollen and ovule numbers, and inflorescence characters important in wind pollination evolved indirectly in response to selection on sex allocation. Responses generally followed predictions from multivariate selection models, with some responses unexpectedly large due to increased genetic correlations as selection proceeded. CONCLUSIONS: Results illustrate the power of artificial selection and utility of multivariate selection models incorporating sex differences. They further indicate that pollen and ovule numbers and inflorescence architecture could evolve in response to selection on biomass allocation to male versus female function, producing complex changes in plant phenotype as separate sexes evolve.

Identification and characterization of miRNAome and target genes in Pseudostellaria heterophylla.

miRNAs play a crucial role in the development and growth of plants by inhibiting the function of targeted genes at the post-transcription level. However, no miRNAs in Pseudostellaria heterophylla have been reported and their function in the morphogenesis of organs is still unclear. In this study, a total of 159 conserved miRNAs (belonging to 64 families) and 303 level miRNAs were identified from P. heterophylla. Some of them showed specifically up or down-regulated expression in different tissues and numbers of unigenes involved in Plant-pathogen interaction and MAPK signaling pathway-plant were targeted. The significant negative correlation of expression profiles between 30 miRNAs and their target genes (37 unigenes) was observed, respectively. Further, a large number of genes involved with signal transduction of auxin, zeatin, abscisic acid and, jasmonic acid were targeted. Predicated targets of two miRNAs were validated by 5'RLM-RACE, respectively. A large number of mRNAs from four pathogens were targeted by miRNAs from P. heterophylla and some of them were targeted by miR414. In summary, we reported a population of miRNAs from four different vegetative tissues of P. heterophylla by high throughput sequencing, which was analyzed by combining with the constructed transcriptome. These results may help to explain the function of miRNAs in the morphogenesis of organs and defense of pathogens, and may provide theoretical basis for breeding and genetic improvement of P. heterophylla.

Molecular Diversity and Phylogeny Reconstruction of Genus Colobanthus (Caryophyllaceae) Based on Mitochondrial Gene Sequences.

Mitochondrial genomes have become an interesting object of evolutionary and systematic study both for animals and plants, including angiosperms. Although the framework of the angiosperm phylogeny was built on the information derived from chloroplast and nuclear genes, mitochondrial sequences also revealed their usefulness in solving the phylogenetic issues at different levels of plant systematics. Here, we report for the first time the complete sequences of 26 protein-coding genes of eight Colobanthus species (Caryophyllaceae). Of these, 23 of them represented core mitochondrial genes, which are directly associated with the primary function of that organelle, and the remaining three genes represented a facultative set of mitochondrial genes. Comparative analysis of the identified genes revealed a generally high degree of sequence conservation. The Ka/Ks ratio was <1 for most of the genes, which indicated purifying selection. Only for rps12 was Ka/Ks > 1 in all studied species, suggesting positive selection. We identified 146−165 potential RNA editing sites in genes of the studied species, which is lower than in most angiosperms. The reconstructed phylogeny based on mitochondrial genes was consistent with the taxonomic position of the studied species, showing the separate character of the family Caryophyllaceae and close relationships between all studied Colobanthus species, with C. lycopodioides sharing less similarity.

Resistance to ALS inhibitors conferred by non-target-site resistance mechanisms in Myosoton aquaticum L.

Myosoton aquaticum L. is a competitive broadleaf weed commonly found in wheat fields in China and has become challenging due to its evolving herbicide resistance. In this study, one subpopulation, RF1 (derived from the tribenuron-methyl-resistant population HN10), with none of the known acetolactate synthase (ALS) resistance mutations was confirmed to exhibit resistance to tribenuron-methyl (SU), pyrithiobac­sodium (PTB), florasulam (TP), flucarbazone-Na (SCT), and diflufenican (PDS). In vitro ALS activity assays showed that the total ALS activity of RF1 was lower than that of the susceptible (S) population. However, there was no difference in ALS gene expression induced by tribenuron-methyl between the two populations. The combination of the cytochrome P450 monooxygenase (P450) inhibitor malathion and tribenuron-methyl resulted in the RF1 population behaving like the S population. The rapid P450-mediated tribenuron-methyl metabolism in RF1 plants was also confirmed by liquid chromatography tandem mass spectrometry (LC-MS/MS) analysis. In addition, approximately equal glutathione S-transferase (GST) activity was observed in RF1 and S plants of untreated and tribenuron-methyl treated groups. This study reported one M. aquaticum L. population without ALS resistance mutations exhibiting resistance to ALS inhibitors and the PDS inhibitor diflufenican, and the non-target-site resistance mechanism played a vital role in herbicide resistance.

Reproductive isolation among lineages of Silene nutans (Caryophyllaceae): A potential involvement of plastid-nuclear incompatibilities.

Early stages of speciation in plants might involve genetic incompatibilities between plastid and nuclear genomes, leading to inter-lineage hybrid breakdown due to the disruption between co-adapted plastid and nuclear genes encoding subunits of the same plastid protein complexes. We tested this hypothesis in Silene nutans, a gynodioecious Caryophyllaceae, where four distinct genetic lineages exhibited strong reproductive isolation among each other, resulting in chlorotic or variegated hybrids. By sequencing the whole gene content of the four plastomes through gene capture, and a large part of the nuclear genes encoding plastid subunits from RNAseq data, we searched for non-synonymous substitutions fixed in each lineage on both genomes. Lineages of S. nutans exhibited a high level of dN/dS ratios for plastid and nuclear genes encoding most plastid complexes, with a strong pattern of coevolution for genes encoding the subunits of ribosome and cytochrome b6/f that could explain the chlorosis of hybrids. Overall, relaxation of selection due to past bottlenecks and positive selection have driven the diversity pattern observed in S. nutans plastid complexes, leading to plastid-nuclear incompatibilities. We discuss the possible role of gynodioecy in the evolutionary dynamics of the plastomes through linked selection.

Root-Associated Bacteria Community Characteristics of Antarctic Plants: Deschampsia antarctica and Colobanthus quitensis-a Comparison.

Colobanthus quitensis (Kunth) Bartl. and Deschampsia antarctica Desv. are the only Magnoliophyta to naturally colonize the Antarctic region. The reason for their sole presence in Antarctica is still debated as there is no definitive consensus on how only two unrelated flowering plants managed to establish breeding populations in this part of the world. In this study, we have explored and compared the rhizosphere and root-endosphere dwelling microbial community of C. quitensis and D. antarctica specimens sampled in maritime Antarctica from sites displaying contrasting edaphic characteristics. Bacterial phylogenetic diversity (high-throughput 16S rRNA gene fragment targeted sequencing) and microbial metabolic activity (Biolog EcoPlates) with a geochemical soil background were assessed. Gathered data showed that the microbiome of C. quitensis root system was mostly site-dependent, displaying different characteristics in each of the examined locations. This plant tolerated an active bacterial community only in severe conditions (salt stress and nutrient deprivation), while in other more favorable circumstances, it restricted microbial activity, with a possibility of microbivory-based nutrient acquisition. The microbial communities of D. antarctica showed a high degree of similarity between samples within a particular rhizocompartment. The grass' endosphere was significantly enriched in plant beneficial taxa of the family Rhizobiaceae, which displayed obligatory endophyte characteristics, suggesting that at least part of this community is transmitted vertically. Ultimately, the ecological success of C. quitensis and D. antarctica in Antarctica might be largely attributed to their associations and management of root-associated microbiota.

De Novo Biosynthesis of the Oleanane-Type Triterpenoids of Tunicosaponins in Yeast.

Tunicosaponins are natural products extracted from Psammosilene tunicoides, which is an important ingredient of Yunnan Baiyao Powder, an ancient and famous Asian herbal medicine. The representative aglycones of tunicosaponins are the oleanane-type triterpenoids of gypsogenin and quillaic acid, which were found to manipulate a broad range of virus-host fusion via wrapping the heptad repeat-2 (HR2) domain prevalent in viral envelopes. However, the unknown biosynthetic pathway and difficulty in chemical synthesis hinder the therapeutic use of tunicosaponins. Here, two novel cytochrome P450-dependent monooxygenases that take part in the biosynthesis of tunicosaponins, CYP716A262 (CYP091) and CYP72A567 (CYP099), were identified from P. tunicoides. In addition, the whole biosynthesis pathway of the tunicosaponin aglycones was reconstituted in yeast by transforming the platform strain BY-bAS with the CYP716A262 and CYP716A567 genes, the resulting strain could produce 146.84 and 314.01 mg/L of gypsogenin and quillaic acid, respectively. This synthetic biology platform for complicated metabolic pathways elucidation and microbial cell factories construction can provide alternative sources of important natural products, helping conserve natural plant resources.

Impact of river dynamics on the genetic variation of Gypsophila repens (Caryophyllaceae): a comparison of heath forest and more dynamic gravel bank populations along an alpine river.

Alpine rivers are, despite anthropogenic water flow regulation, still often highly dynamic ecosystems. Plant species occurring along these rivers are subject to ecological disturbance, mainly caused by seasonal flooding. Gypsophila repens typically grows at higher altitudes in the Alps, but also occurs at lower altitudes on gravel banks directly along the river and in heath forests at larger distances from the river. Populations on gravel banks are considered non-permanent and it is assumed that new individuals originate from seed periodically washed down from higher altitudes. Populations in heath forests are, in contrast, permanent and not regularly provided with seeds from higher altitudes through flooding. If the genetic structure of this plant species is strongly affected by gene flow via seed dispersal, then higher levels of genetic diversity in populations but less differentiation among populations on gravel banks than in heath forests can be expected. In this study, we analysed genetic diversity within and differentiation among 15 populations of G. repens from gravel banks and heath forests along the alpine River Isar using amplified fragment length polymorphisms (AFLP). Genetic diversity was, as assumed, slightly higher in gravel bank than in heath forest populations, but genetic differentiation was, in contrast to our expectations, comparable among populations in both habitat types. Our study provides evidence for increased genetic diversity under conditions of higher ecological disturbance and increased seed dispersal on gravel banks. Similar levels of genetic differentiation among populations in both habitat types can be attributed to the species' long lifetime, a permanent soil seed bank and gene flow by pollinators among different habitats/locations.

Evolutionary dynamics of the chloroplast genome sequences of six Colobanthus species.

The complete plastome sequences of six species were sequenced to better understand the evolutionary relationships and mutation patterns in the chloroplast genome of the genus Colobanthus. The length of the chloroplast genome sequences of C. acicularis, C. affinis, C. lycopodioides, C. nivicola, C. pulvinatus and C. subulatus ranged from 151,050 to 151,462 bp. The quadripartite circular structure of these genome sequences has the same overall organization and gene content with 73 protein-coding genes, 30 tRNA genes, four rRNA genes and five conserved chloroplast open reading frames. A total of 153 repeat sequences were revealed. Forward repeats were dominant, whereas complementary repeats were found only in C. pulvinatus. The mononucleotide SSRs composed of A/T units were most common, and hexanucleotide SSRs were detected least often. Eleven highly variable regions which could be utilized as potential markers for phylogeny reconstruction, species identification or phylogeography were identified within Colobanthus chloroplast genomes. Seventy-three protein-coding genes were used in phylogenetic analyses. Reconstructed phylogeny was consistent with the systematic position of the studied species, and the representatives of the same genus were grouped in one clade. All studied Colobanthus species formed a single group and C. lycopodioides was least similar to the remaining species.

Identification of High Molecular Variation Loci in Complete Chloroplast Genomes of Mammillaria (Cactaceae, Caryophyllales).

In plants, partial DNA sequences of chloroplasts have been widely used in evolutionary studies. However, the Cactaceae family (1500-1800 species) lacks molecular markers that allow a phylogenetic resolution between species and genera. In order to identify sequences with high variation levels, we compared previously reported complete chloroplast genomes of seven species of Mammillaria. We identified repeated sequences (RSs) and two types of DNA variation: short sequence repeats (SSRs) and divergent homologous loci. The species with the highest number of RSs was M. solisioides (256), whereas M. pectinifera contained the highest amount of SSRs (84). In contrast, M. zephyranthoides contained the lowest number (35) of both RSs and SSRs. In addition, five of the SSRs were found in the seven species, but only three of them showed variation. A total of 180 homologous loci were identified among the seven species. Out of these, 20 loci showed a molecular variation of 5% to 31%, and 12 had a length within the range of 150 to 1000 bp. We conclude that the high levels of variation at the reported loci represent valuable knowledge that may help to resolve phylogenetic relationships and that may potentially be convenient as molecular markers for population genetics and phylogeographic studies.

Using species distribution models to locate the potential cradles of the allopolyploid Gypsophila bermejoi G. López (Caryophyllaceae).

Polyploidy has been an influential force in plant evolution, playing a crucial role in diversification. Differences in polyploid and diploid distributions have been long noted, with polyploid taxa especially abundant in harsh environments. These plants have higher photosynthetic rates and/or higher tolerance to water deficits. Moreover, there is data pointing to an increase in the rate of unreduced gamete formation by plants under conditions of stress. Accordingly, a higher frequency of polyploid individuals would be expected in populations living under extreme environments, a phenomenon that may be relevant when considering the origin of allopolyploid species. Hybridization between distinct autopolyploids is known to produce allopolyploids and hence, a high frequency of compatible autopolyploids in an area could enhance the formation of stable populations of the corresponding allopolyploid hybrid. Here we consider the allopolyploid species Gypsophila bermejoi G. López and its parental taxa G. struthium L. subsp. struthium and G. tomentosa L. We have used Species Distribution Models to locate areas with low bioclimatic suitability for both parental taxa during the Last Glacial Maximum (LGM), hypothesizing that the rate of tetraploid hybrid formation would be higher than expected where low suitability areas of both parental species overlap. We selected those areas taking into account the strict gypsophyllic nature of these taxa. There is data pointing to a post-glacial origin of the current G. bermejoi populations and according to our hypothesis, such locations could be centers for hybrid tetraploid formation or potential cradles of this species. Indeed, potential Mid-Holocene cradles were also identified in this manner. The evolution of bioclimatic suitability in both LGM and Mid-Holocene cradles was studied to assess the possible survival of the hybrids, and the current distribution of G. bermejoi proved to be consistent with our hypothesis.