SummaryGenome size variation is of fundamental biological importance and has been a longstanding puzzle in evolutionary biology. Several hypotheses for genome size evolution including neutral, maladaptive, and adaptive models have been proposed, but the relative importance of these models remains controversial.Primulina is a genus that is highly diversified in the Karst region of southern China, where genome size variation and the underlying evolutionary mechanisms are poorly understood. We reconstructed the phylogeny of Primulina using DNA sequences for 104 species and determined the genome sizes of 101 species. We examined the phylogenetic signal in genome size variation, and tested the fit to different evolutionary models and for correlations with variation in latitude and specific leaf area (SLA).The results showed that genome size, SLA and latitudinal variation all displayed strong phylogenetic signals, but were best explained by different evolutionary models. Furthermore, significant positive relationships were detected between genome size and SLA and between genome size and latitude.Our study is the first to investigate genome size evolution on such a comprehensive scale and in the Karst region flora. We conclude that genome size in Primulina is phylogenetically conserved but its variation among species is a combined outcome of both neutral and adaptive evolution.

The genus Primulina is an emerging model system in studying the drivers and mechanisms of species diversification, for its high species richness and endemism, together with high degree of habitat specialization. In this study, we sequenced transcriptomes for eleven Primulina species across the phylogeny of the genus using the Illumina HiSeq 2000 platform. A total of 336 million clean reads were processed into 355 573 unigenes with a mean length of 1336 bp and an N50 value of 2191 bp after pooling and reassembling twelve individual pre-assembled unigene sets. Of these unigenes, 249 973 (70%) were successfully annotated and 256 601 (72%) were identified as coding sequences (CDSs). We identified a total of 38 279 simple sequence repeats (SSRs) and 367 123 single nucleotide polymorphisms (SNPs). Marker validation assay revealed that 354 (27.3%) of the 1296 SSR and 795 (39.6%) of the 2008 SNP loci showed successful genotyping performance and exhibited expected polymorphism profiles. We screened 834 putative single-copy nuclear genes and proved their high effectiveness in phylogeny construction and estimation of ancestral population parameters. We identified a total of 85 candidate orthologs under positive selection for 46 of the 66 species pairs. This study provided an efficient application of RNA-seq in development of genomic resources for a group of 'stone plants' from south China Karst regions, a biodiversity hot spot of the World. The assembled unigenes with annotations and the massive gene-associated molecular markers would help guide further molecular systematic, population genetic and ecological genomics studies in Primulina and its relatives.

The widespread ascorbic acid (AsA) plays a vital role in plant development and abiotic stress tolerance, but AsA concentration varies greatly among different plants. GDP-D-mannose epimerase (GME), which catalyzes GDP-D-mannose to GDP-L-galactose or GDP-L-gulose, is a key enzyme in plant AsA biosynthesis pathway. Functions and expression patterns of GME have been well studied in previous works, however, little information is known about the evolutionary patterns of the gene. In this study, GME gene structure, corresponding conserved protein motifs and evolutionary relationships were systematically analyzed. A total of 111 GME gene sequences were retrieved from 59 plant genomes, which representing almost all the major lineages of Viridiplantae: dicotyledons, monocotyledons, gymnosperms, pteridophytes, bryophytes, and chlorophytes. Results showed that homologs of GME were widely present in Viridiplantae. GME gene structures were conservative in higher plants, while varied greatly in the basal subgroups of the phylogeny including lycophytes, bryophytes, and chlorophytes, suggesting GME gene structure might have undergone severe differentiation at lower plant and then gradually fixed as plant evolution. The basic motifs of GME were strongly conserved throughout Viridiplantae, suggesting the conserved function of the protein. Molecular evolution analysis showed that strong purifying selection was the predominant force in the evolution of GME. A few branches and sites under episodic diversifying selection were identified and most of the branches located in the subgroup of chlorphytes, indicating episodic diversifying selection at a few branches and sites may play a role in the evolution of GME and diversifying selection may have occurred at the early stage of Viridiplantae. Our results provide novel insights into functional conservation and the evolution of GME.