Plastome evolution and phylogenomic insights into the evolution of Lysimachia (Primulaceae: Myrsinoideae).

BACKGROUND: Lysimachia L., the second largest genus within the subfamily Myrsinoideae of Primulaceae, comprises approximately 250 species worldwide. China is the species diversity center of Lysimachia, containing approximately 150 species. Despite advances in the backbone phylogeny of Lysimachia, species-level relationships remain poorly understood due to limited genomic information. This study analyzed 50 complete plastomes for 46 Lysimachia species. We aimed to identify the plastome structure features and hypervariable loci of Lysimachia. Additionally, the phylogenetic relationships and phylogenetic conflict signals in Lysimachia were examined. RESULTS: These fifty plastomes within Lysimachia had the typical quadripartite structure, with lengths varying from 152,691 to 155,784 bp. Plastome size was positively correlated with IR and intron length. Thirteen highly variable regions in Lysimachia plastomes were identified. Additionally, ndhB, petB and ycf2 were found to be under positive selection. Plastid ML trees and species tree strongly supported that L. maritima as sister to subg. Palladia + subg. Lysimachia (Christinae clade), while the nrDNA ML tree clearly placed L. maritima and subg. Palladia as a sister group. CONCLUSIONS: The structures of these plastomes of Lysimachia were generally conserved, but potential plastid markers and signatures of positive selection were detected. These genomic data provided new insights into the interspecific relationships of Lysimachia, including the cytonuclear discordance of the position of L. maritima, which may be the result of ghost introgression in the past. Our findings have established a basis for further exploration of the taxonomy, phylogeny and evolutionary history within Lysimachia.

The phylogeny and global biogeography of Primulaceae based on high-throughput DNA sequence data.

The angiosperm family Primulaceae is morphologically diverse and distributed nearly worldwide. However, phylogenetic uncertainty has obstructed the identification of major morphological and biogeographic transitions within the clade. We used target capture sequencing with the Angiosperms353 probes, taxon-sampling encompassing nearly all genera of the family, tree-based sequence curation, and multiple phylogenetic approaches to investigate the major clades of Primulaceae and their relationship to other Ericales. We generated dated phylogenetic trees and conducted broad-scale biogeographic analyses as well as stochastic character mapping of growth habit. We show that Ardisia, a pantropical genus and the largest in the family, is not monophyletic, with at least 19 smaller genera nested within it. Neotropical members of Ardisia and several smaller genera form a clade, an ancestor of which arrived in the Neotropics and began diversifying about 20 Ma. This Neotropical clade is most closely related to Elingamita and Tapeinosperma, which are most diverse on islands of the Pacific. Both Androsace and Primula are non-monophyletic by the inclusion of smaller genera. Ancestral state reconstructions revealed that there have either been parallel transitions to an herbaceous habit in Primuloideae, Samolus, and at least three lineages of Myrsinoideae, or a common ancestor of nearly all Primulaceae was herbaceous. Our results provide a robust estimate of phylogenetic relationships across Primulaceae and show that a revised classification of Myrsinoideae and several other clades within the family is necessary to render all genera monophyletic.

Pervasive Introgression During Rapid Diversification of the European Mountain Genus Soldanella (L.) (Primulaceae).

Hybridization is a key mechanism involved in lineage diversification and speciation, especially in ecosystems that experienced repeated environmental oscillations. Recently radiated plant groups, which have evolved in mountain ecosystems impacted by historical climate change provide an excellent model system for studying the impact of gene flow on speciation. We combined organellar (whole-plastome) and nuclear genomic data (RAD-seq) with a cytogenetic approach (rDNA FISH) to investigate the effects of hybridization and introgression on evolution and speciation in the genus Soldanella (snowbells, Primulaceae). Pervasive introgression has already occurred among ancestral lineages of snowbells and has persisted throughout the entire evolutionary history of the genus, regardless of the ecology, cytotype, or distribution range size of the affected species. The highest extent of introgression has been detected in the Carpathian species, which is also reflected in their extensive karyotype variation. Introgression occurred even between species with dysploid and euploid cytotypes, which were considered to be reproductively isolated. The magnitude of introgression detected in snowbells is unprecedented in other mountain genera of the European Alpine System investigated hitherto. Our study stresses the prominent evolutionary role of hybridization in facilitating speciation and diversification on the one hand, but also enriching previously isolated genetic pools. [chloroplast capture; diversification; dysploidy; European Alpine system; introgression; nuclear-cytoplasmic discordance; ribosomal DNA.].

The Complete Chloroplast Genomes of Primula obconica Provide Insight That Neither Species nor Natural Section Represent Monophyletic Taxa in Primula (Primulaceae).

The genus Primula (Primulaceae) comprises more than 500 species, with 300 species distributed in China. The contradictory results between systematic analyses and morphology-based taxonomy make taxonomy studies difficult. Furthermore, frequent introgression between closely related species of Primula can result in non-monophyletic species. In this study, the complete chloroplast genome of sixteen Primula obconica subsp. obconica individuals were assembled and compared with 84 accessions of 74 species from 21 sections of the 24 sections of the genus in China. The plastome sizes of P. obconica subsp. obconica range from 153,584 bp to 154,028 bp. Genome-wide variations were detected, and 1915 high-quality SNPs and 346 InDels were found. Most SNPs were detected in downstream and upstream gene regions (45.549% and 41.91%). Two cultivated accessions, ZP1 and ZP2, were abundant with SSRs. Moreover, 12 SSRs shared by 9 accessions showed variations that may be used as molecular markers for population genetic studies. The phylogenetic tree showed that P. obconica subsp. obconica cluster into two independent clades. Two subspecies have highly recognizable morphological characteristics, isolated geographical distribution areas, and distinct phylogenetic relationships compared with P. obconica subsp. obconica. We elevate the two subspecies of P. obconica to separate species. Our phylogenetic tree is largely inconsistent with morphology-based taxonomy. Twenty-one sections of Primula were mainly divided into three clades. The monophyly of Sect. Auganthus, Sect. Minutissimae, Sect. Sikkimensis, Sect. Petiolares, and Sect. Ranunculoides are well supported in the phylogenetic tree. The Sect. Obconicolisteri, Sect. Monocarpicae, Sect. Carolinella, Sect. Cortusoides, Sect. Aleuritia, Sect. Denticulata, Sect. Proliferae Pax, and Sect. Crystallophlomis are not a monophyletic group. The possible explanations for non-monophyly may be hybridization, polyploidization, recent introgression, incorrect taxonomy, or chloroplast capture. Multiple genomic data and population genetic studies are therefore needed to reveal the evolutionary history of Primula. Our results provided valuable information for intraspecific variation and phylogenetic relationships within Primula.

Population subdivision promoted by a sea-level-change-driven bottleneck: A glimpse from the evolutionary history of the mangrove plant Aegiceras corniculatum.

Historic climate changes drive geographical populations of coastal plants to contract and recover dynamically, even die out completely. Species suffering from such bottlenecks usually lose intraspecific genetic diversity, but how do these events influence population subdivision patterns of coastal plants? Here, we investigated this question in the typical coastal plant: mangrove species Aegiceras corniculatum. Inhabiting the intertidal zone of the tropical and subtropical coast of the Indo-West Pacific oceans, its populations are deemed to be greatly shaped by historic sea-level fluctuations. Using dual methods of Sanger and Illumina sequencing, we found that the 18 sampled populations were structured into two groups, namely, the "Indo-Malayan" group, comprising three subgroups (the northern South China Sea, Gulf of Bengal, and Bali), and the "Pan-Australasia" group, comprising the subgroups of the southern South China Sea and Australasia. Based on the approximate Bayesian computations and Stairway Plot, we inferred that the southern South China Sea subgroup, which penetrates the interior of the "Indo-Malayan" group, originated from the Australasia subgroup, accompanied by a severe bottleneck event, with a spot of gene flow from both the Australasia and "Indo-Malayan" groups. Geographical barriers such as the Sundaland underlie the genetic break between Indian and Pacific Oceans, but the discontinuity between southern and northern South China Sea was originated from genetic drift in the bottleneck event. Hence, we revealed a case evidencing that the bottleneck event promoted population subdivision. This conclusion may be applicable in other taxa beyond coastal plants.

Phylogeographic pattern of a cryptoviviparous mangrove, Aegiceras corniculatum, in the Indo-West Pacific, provides insights for conservation actions.

MAIN CONCLUSION: This study identified the historical geoclimatic factors which caused low genetic diversity and strong phylogeographic structure in a cryptoviviparous mangrove. The phylogeographic pattern was used to suggest conservation actions. Phylogeographic studies are used to understand the spatial distribution and evolution of genetic diversity, and have major conservation implications, especially for threatened taxa like the mangroves. This study aimed to assess the phylogeographic pattern of Aegiceras corniculatum, a cryptoviviparous mangrove, across its distribution range in the Indo-West Pacific (IWP) region. We genotyped 398 samples, collected from 37 populations, at four chloroplast DNA (cpDNA) loci, and identified the influence of historical processes on the contemporary population structure of the species. Low genetic diversity at the population level was observed. The evolutionary relationship between 12 cpDNA haplotypes suggested a strong phylogeographic structure, which was further validated by the clustering algorithms and proportioning of maximum variation among hierarchical population groups. The magnitude and direction of historical gene flow indicated that the species attained its wide distribution from its likely ancestral area of the Malay Archipelago. The divergence time estimates of the haplotypes indicated that the geoclimatic changes during the Pleistocene, especially the glacial sea-level changes and emergence of landmasses, hindered genetic exchange and created genetic differentiation between the phylogenetic groups. The species overwintered the last glacial maxima in multiple refugia in the IWP, as identified by the environmental niche modelling. Overall, our findings indicated that ancient glacial vicariance had influenced the present genetic composition of A. corniculatum, which was maintained by the current demographic features of this region. We discussed how these findings can be used to prioritize areas for conservation actions, restore disturbed habitats and prevent further genetic erosion.

Identification, characterization and expression analysis of lineage-specific genes within mangrove species Aegiceras corniculatum.

Lineage-specific genes (LSGs) are the genes that have no recognizable homology to any sequences in other species, which are important drivers for the generation of new functions, phenotypic changes, and facilitating species adaptation to environment. Aegiceras corniculatum is one of major mangrove plant species adapted to waterlogging and saline conditions, and the exploration of aegiceras-specific genes (ASGs) is important to reveal its adaptation to the harsh environment. Here, we performed a systematic analysis on ASGs, focusing on their sequence characterization, origination and expression patterns. Our results reveal that there are 4823 ASGs in the genome, approximately 11.84% of all protein-coding genes. High proportion (45.78%) of ASGs originate from gene duplication, and the time of gene duplication of ASGs is consistent with the timing of two genome-wide replication (WGD) events that occurred in A. corniculatum, and also coincides with a short period of global warming during the Paleocene-Eocene Maximum (PETM, 55.5 million years ago). Gene structure analysis showed that ASGs have shorter protein lengths, fewer exons, and higher isoelectric point. Expression patterns analysis showed that ASGs had low levels of expression and more tissue-specific expression. Weighted gene co-expression network analysis (WGCNA) revealed that 86 ASGs co-expressed gene modules were primarily involved in pathways related to adversity stress, including plant hormone signal transduction, phenylpropanoid biosynthesis, photosynthesis, peroxisome and pentose phosphate pathway. This study provides a comprehensive analysis of the characteristics and potential functions of ASGs and identifies key candidate genes, which will contribute to the subsequent further investigation of the adaptation of A. corniculatum to intertidal coastal wetland habitats.

Genomic insights into molecular adaptation to intertidal environments in the mangrove Aegiceras corniculatum.

Mangroves have colonised extreme intertidal environments characterised by high salinity, hypoxia and other abiotic stresses. Aegiceras corniculatum, a pioneer mangrove species that has evolved two specialised adaptive traits (salt secretion and crypto-vivipary) is an attractive ecological model to investigate molecular mechanisms underlying adaptation to intertidal environments. We assembled de novo a high-quality reference genome of A. corniculatum and performed comparative genomic and transcriptomic analyses to investigate molecular mechanisms underlying adaptation to intertidal environments. We provide evidence that A. corniculatum experienced a whole-genome duplication (WGD) event c. 35 Ma. We infer that maintenance of cellular environmental homeostasis is an important adaptive process in A. corniculatum. The 14-3-3 and H+ -ATPase protein-coding genes, essential for the salt homeostasis, were preferentially retained after the recent WGD event. Using comparative transcriptomics, we show that genes upregulated under high-salt conditions are involved in salt transport and ROS scavenging. We also found that all homologues of DELAY OF GERMINATION1 (DOG1) had lost their heme-binding ability in A. corniculatum, and that this may contribute to crypto-vivipary. Our study provides insight into the genomic correlates of phenotypic adaptation to intertidal environments. This could contribute not only within the genomics community, but also to the field of plant evolution.

Chromosome-level assembly of the mangrove plant Aegiceras corniculatum genome generated through Illumina, PacBio and Hi-C sequencing technologies.

Aegiceras corniculatum is a major mangrove plant species adapted to waterlogging and saline conditions, grows in the coastal intertidal zone of tropical and subtropical regions. Here, we present a chromosome-level genome assembly of A. corniculatum by incorporating PacBio long-read sequencing and Hi-C technology. The results showed that the PacBio draft genome size is 906.63 Mb. Hi-C scaffolding anchored 885.06 Mb contigs (97.62% of draft assembly) onto 24 pseudochromosomes. The contig N50 and scaffold N50 were 7.1 Mb and 37.74 Mb, respectively. Out of 40,727 protein-coding genes predicted in the study, 89% have functional annotations in public databases. We also showed that of the 603.93 Mb repetitive sequences predicted in the assembled genome, long terminal repeat retrotransposons constitute 41.52%. The genome evolution analysis showed that the A. corniculatum genome experienced two whole-genome duplication events and shared the ancient γ whole-genome triplication event. A comparative genomic analysis revealed an incidence of expansion in 1,488 gene families associated with essential metabolism and biosynthetic pathways, including photosynthesis, oxidative phosphorylation, phenylalanine, glyoxylate, dicarboxylate metabolism, and DNA replication, which probably constitute adaptation traits that allow the A. corniculatum to survive in the intertidal zone. Also, the systematic characterization of genes associated with flavonoid biosynthesis pathway and the AcNHX gene family conducted in this study will provide insight into the adaptation mechanism of A. corniculatum to intertidal environments. The high-quality genome reported here can provide historical insights into genomic transformations that support the survival of A. corniculatum under harsh intertidal habitats.

Chloroplast Genomes and Comparative Analyses among Thirteen Taxa within Myrsinaceae s.str. Clade (Myrsinoideae, Primulaceae).

The Myrsinaceae s.str. clade is a tropical woody representative in Myrsinoideae of Primulaceae and has ca. 1300 species. The generic limits and alignments of this clade are unclear due to the limited number of genetic markers and/or taxon samplings in previous studies. Here, the chloroplast (cp) genomes of 13 taxa within the Myrsinaceae s.str. clade are sequenced and characterized. These cp genomes are typical quadripartite circle molecules and are highly conserved in size and gene content. Three pseudogenes are identified, of which ycf15 is totally absent from five taxa. Noncoding and large single copy region (LSC) exhibit higher levels of nucleotide diversity (Pi) than other regions. A total of ten hotspot fragments and 796 chloroplast simple sequence repeats (SSR) loci are found across all cp genomes. The results of phylogenetic analysis support the notion that the monophyletic Myrsinaceae s.str. clade has two subclades. Non-synonymous substitution rates (dN) are higher in housekeeping (HK) genes than photosynthetic (PS) genes, but both groups have a nearly identical synonymous substitution rate (dS). The results indicate that the PS genes are under stronger functional constraints compared with the HK genes. Overall, the study provides hypervariable molecular markers for phylogenetic reconstruction and contributes to a better understanding of plastid gene evolution in Myrsinaceae s.str. clade.

What explains high plant richness in East Asia? Time and diversification in the tribe Lysimachieae (Primulaceae).

What causes the disparity in biodiversity among regions is a fundamental question in biogeography, ecology, and evolutionary biology. Evolutionary and biogeographic processes (speciation, extinction, dispersal) directly determine species richness patterns, and can be studied using integrative phylogenetic approaches. However, the strikingly high richness of East Asia relative to other Northern Hemisphere regions remains poorly understood from this perspective. Here, for the first time, we test two general hypotheses (older colonization time, faster diversification rate) to explain this pattern, using the plant tribe Lysimachieae (Primulaceae) as a model system. We generated a new time-calibrated phylogeny for Lysimachieae (13 genes, 126 species), to estimate colonization times and diversification rates for each region and to test the relative importance of these two factors for explaining regional richness patterns. We find that neither time nor diversification rates alone explain richness patterns among regions in Lysimachieae. Instead, a new index that combines both factors explains global richness patterns in the group and their high East Asian biodiversity. Based on our results from Lysimachieae, we suggest that the high richness of plants in East Asia may be explained by a combination of older colonization times and faster diversification rates in this region.

Molecular thresholds of ITS2 and their implications for molecular evolution and species identification in seed plants.

Although molecular data have revealed huge amounts of plant diversity, interpreting genetic diversity into entities corresponding to species is still challenging. Taxonomic ranking based on genetic distance has been used extensively, but the results have been open to dispute, while the application of the strategy to plants has been restricted to a small number of cases. Here, levels of internal transcribed spacer 2 (ITS2) sequence variation were examined from 17,203 sequences, representing 5,439 species in 113 genera of seed plants, to ascertain the association between species status and their molecular divergence. Our results showed that, although the average genetic distances of sister species (AGDS) varied among angiosperms, the mean value was 3.98% and seemed not to be influenced by higher-level hierarchical classification or life history. AGDS was also stable within the major lineages of the gymnosperms but at approximately half the value of angiosperms, except for the Gnetidae, where the AGDS almost equaled that of angiosperms. We found that these AGDS discrepancies, associated with the rates of molecular evolution, cannot simply be attributed to generation-time differences, and highlight the complex life histories of plants. Our results provide general ITS2 thresholds in seed plants, and suggest their use in species identification.

Understanding microRNA Regulation Involved in the Metamorphosis of the Veined Rapa Whelk (Rapana venosa).

The veined rapa whelk (Rapana venosa) is widely consumed in China. Nevertheless, it preys on oceanic bivalves, thereby reducing this resource worldwide. Its larval metamorphosis comprises a transition from pelagic to benthic form, which involves considerable physiological and structural changes and has vital roles in its natural populations and commercial breeding. Thus, understanding the endogenous microRNAs (miRNAs) that drive metamorphosis is of great interest. This is the first study to use high-throughput sequencing to examine the alterations in miRNA expression that occur during metamorphosis in a marine gastropod. A total of 195 differentially expressed miRNAs were obtained. Sixty-five of these were expressed during the transition from precompetent to competent larvae. Thirty-three of these were upregulated and the others were downregulated. Another 123 miRNAs were expressed during the transition from competent to postlarvae. Ninety-six of these were upregulated and the remaining 27 were downregulated. The expression of miR-276-y, miR-100-x, miR-183-x, and miR-263-x showed a >100-fold change during development, while the miR-242-x and novel-m0052-3p expression levels changed over 3000-fold. Putative target gene coexpression, gene ontology, and pathway analyses suggest that these miRNAs play important parts in cell proliferation, migration, apoptosis, metabolic regulation, and energy absorption. Twenty miRNAs and their target genes involved in ingestion, digestion, cytoskeleton, cell adhesion, and apoptosis were identified. Nine of them were analyzed with real-time polymerase chain reaction (PCR), which showed an inverse correlation between the miRNAs and their relative expression levels. Our data elucidate the role of miRNAs in R. venosa metamorphic transition and serve as a solid basis for further investigations into regulatory mechanisms of gastropod metamorphosis.

De novo assembly of the transcriptome of Aegiceras corniculatum, a mangrove species in the Indo-West Pacific region.

Aegiceras corniculatum (L.) Blanco is one of the most salt tolerant mangrove species and can thrive in 3% salinity at the seaward edge of mangrove forests. Here we sequenced the transcriptome of A. corniculatum used Illumina GA platform to develop its genomic resources for ecological and evolutionary studies. We obtained about 50 million high-quality paired-end reads with 75bp in length. Using the short read assembler Velvet, we yielded 49,437 contigs with the average length of 625bp. A total of 32,744 (66.23%) contigs showed significant similarity to the GenBank non-redundant (NR) protein database. 30,911 and 18,004 of these sequences were assigned to Gene Ontology and eukaryotic orthologous groups of proteins (KOG). A total of 4942 transcripts from our assemblies had significant similarity with KEGG Orthologs and were involved in 144 KEGG pathways, while 9899 unigenes had enzyme commission (EC) numbers. In addition, 9792 transcriptome-derived SSRs were identified from 7342 sequences. With our strict criteria, 4165 candidate SNPs were also identified from 2058 contigs. Some of these SNPs were further validated by Sanger sequencing. Genomic resources generated in this study should be valuable in ecological, evolutionary, and functional genomics studies for this mangrove species.

Complete chloroplast genome sequence of Lysimachia coreana (Primulaceae).

The complete chloroplast genome of Lysimachia coreana, an endemic plant in Korea, has been studied in this article. The genome size is 155,386 bp in length with a typical quadripartite structure, including the large single-copy region (85,325 bp), small single-copy region (17,951 bp) and two inverted repeated regions (26,055 bp). The genome consisted of 132 coding genes, including 87 protein-coding genes, 37 tRNA genes and 8 rRNA genes.

[A Cellular Automata Model for a Community Comprising Two Plant Species of Different Growth Forms].

A cellular automata computer model for the interactions between two plant species of different growth forms--the lime hairgrass Deschampsia caespitosa (L.) P. Beauv., a sod cereal, and the moneywort Lysimachia nummularia L., a ground creeping perennial herb--is considered. Computer experiments on the self-maintenance of the populations of each species against the background of a gradual increase in the share of randomly eliminated individuals, coexistence of the populations of two species, and the effect of the phytogenous field have been conducted. As has been shown, all the studied factors determine the number of individuals and self-sustainability of the simulated populations by the degree of their impact. The limits of action have been determined for individual factors; within these limits, the specific features in plant reproduction and dispersal provide sustainable coexistence of the simulated populations. It has been demonstrated that the constructed model allows for studying the long-term developmental dynamics of the plants belonging to the selected growth forms.

Automation and Evaluation of the SOWH Test with SOWHAT.

The Swofford-Olsen-Waddell-Hillis (SOWH) test evaluates statistical support for incongruent phylogenetic topologies. It is commonly applied to determine if the maximum likelihood tree in a phylogenetic analysis is significantly different than an alternative hypothesis. The SOWH test compares the observed difference in log-likelihood between two topologies to a null distribution of differences in log-likelihood generated by parametric resampling. The test is a well-established phylogenetic method for topology testing, but it is sensitive to model misspecification, it is computationally burdensome to perform, and its implementation requires the investigator to make several decisions that each have the potential to affect the outcome of the test. We analyzed the effects of multiple factors using seven data sets to which the SOWH test was previously applied. These factors include a number of sample replicates, likelihood software, the introduction of gaps to simulated data, the use of distinct models of evolution for data simulation and likelihood inference, and a suggested test correction wherein an unresolved "zero-constrained" tree is used to simulate sequence data. To facilitate these analyses and future applications of the SOWH test, we wrote SOWHAT, a program that automates the SOWH test. We find that inadequate bootstrap sampling can change the outcome of the SOWH test. The results also show that using a zero-constrained tree for data simulation can result in a wider null distribution and higher p-values, but does not change the outcome of the SOWH test for most of the data sets tested here. These results will help others implement and evaluate the SOWH test and allow us to provide recommendations for future applications of the SOWH test. SOWHAT is available for download from https://github.com/josephryan/SOWHAT.

Evolution, ecology and systematics of Soldanella (Primulaceae) in the southern Apennines (Italy).

BACKGROUND: The populations of Soldanella (Primulaceae) of the southern Apennines (Italy) are unique within the genus for their distribution and ecology. Their highly fragmented distribution range, with three main metapopulations on some of the highest mountains (Gelbison, Sila and Aspromonte massifs) of the area, poses intriguing questions about their evolutionary history and biogeography, and about the possibility of local endemisms. AIMS AND METHODS: In order to clarify the phylogeny and biogeography of the three metapopulations of Soldanella in the southern Apennines, attributed to S. calabrella to date, and to identify possible local endemisms, a comparative approach based on the study of molecular, morphological and ecological characteristics of the populations was employed. Specifically, one nuclear (total ITS) and two plastid (rbcL and trnL) markers were used for the phylogenetic analyses, performed through both maximum likelihood and Bayesian techniques. Among the morphological features, the glandular hair and leaf biometric traits were analysed, and the environment in which the populations grew was characterised for altitude, forest canopy composition and soil pH, C, N and organic matter. RESULTS AND CONCLUSIONS: Our findings demonstrate that the lineage of Soldanella of southern Italy diverged from the Carpathians lineage during the Middle Pleistocene, and underwent an evolutionary radiation during the Late Pleistocene. The populations of the Sila and Aspromonte massifs diverged from the populations of the Gelbison massif around 380000 years ago and are probably undergoing a progressive differentiation due to their isolation. The populations on the Gelbison massif, moreover, have different morphological features from those of the Sila and Aspromonte massifs and a different ecological niche. The molecular, morphological and ecological data clearly demonstrate that the metapopulation of Soldanella on the Gelbison massif belongs to a new taxonomic unit at the species level, which we name Soldanella sacra A. & L. Bellino from the name of the massif on which it was discovered, the "Holy Mountain".

Lysimachia huangsangensis (Primulaceae), a New Species from Hunan, China.

A new species, Lysimachia huangsangensis (Primulaceae), from Hunan, China is described and illustrated. The new species is closely related to L. carinata because of the crested calyx, but differs in the leaf blades that are ovate to elliptic and (3-)4.5-9 × 2-3.4 cm, 2-5-flowered racemes, and the calyx lobes that are ovate-lanceolate and 5-6 × 3-4 mm. The systematic placement and conservation status are also discussed.

Cloning of the Aegiceras corniculatum class I chitinase gene (AcCHI I) and the response of AcCHI I mRNA expression to cadmium stress.

Chitinases in terrestrial plants have been reported these are involved in heavy metal tolerance/detoxification. This is the first attempt to reveal chitinase gene (AcCHI I) and its function on metal detoxification in mangroves Aegiceras corniculatum. RT-PCR and RACE techniques were used to clone AcCHI I, while real-time quantitative PCR was employed to assess AcCHI I mRNA expressions in response to Cadmium (Cd). The deduced AcCHI I protein consists of 316 amino acids, including a signal peptide region, a chitin-binding domain (CBD) and a catalytic domain. Protein homology modeling was performed to identify potential features in AcCHI I. The CBD structure of AcCHI I might be critical for metal tolerance/homeostasis of the plant. Clear tissue-specific differences in AcCHI I expression were detected, with higher transcript levels detected in leaves. Results demonstrated that a short duration of Cd exposure (e.g., 3 days) promoted AcCHI I expression in roots. Upregulated expression was also detected in leaves under 10 mg/kg Cd concentration stress. The present study demonstrates that AcCHI I may play an important role in Cd tolerance/homeostasis in the plant. Further studies of the AcCHI I protein, gene overexpression, the promoter and upstream regulation will be necessary for clarifying the functions of AcCHI I.