Untangling the influence of phylogeny, soil and climate on leaf element concentrations in a biodiversity hotspot

Combined full-length transcriptomic and metabolomic analysis reveals the molecular mechanisms underlying nutrients and taste components development in Primulina juliae

BACKGROUND

Primulina juliae has recently emerged as a novel functional vegetable, boasting a significant biomass and high calcium content. Various breeding strategies have been employed to the domestication of P. juliae. However, the absence of genome and transcriptome information has hindered the research of mechanisms governing the taste and nutrients in this plant. In this study, we conducted a comprehensive analysis, combining the full-length transcriptomics and metabolomics, to unveil the molecular mechanisms responsible for the development of nutrients and taste components in P. juliae.

RESULTS

We obtain a high-quality reference transcriptome of P. juliae by combing the PacBio Iso-seq and Illumina sequencing technologies. A total of 58,536 cluster consensus sequences were obtained, including 28,168 complete protein coding transcripts and 8,021 Long Non-coding RNAs. Significant differences were observed in the composition and content of compounds related to nutrients and taste, particularly flavonoids, during the leaf development. Our results showed a decrease in the content of most flavonoids as leaves develop. Malate and succinate accumulated with leaf development, while some sugar metabolites were decreased. Furthermore, we identified the different accumulation of amino acids and fatty acids, which are associated with taste traits. Moreover, our transcriptomic analysis provided a molecular basis for understanding the metabolic variations during leaf development. We identified 4,689 differentially expressed genes in the two developmental stages, and through a comprehensive transcriptome and metabolome analysis, we discovered the key structure genes and transcription factors involved in the pathways.

CONCLUSIONS

This study provides a high-quality reference transcriptome and reveals molecular mechanisms associated with the development of nutrients and taste components in P. juliae. These findings will enhance our understanding of the breeding and utilization of P. juliae as a vegetable.

Characterization of Firmiana danxiaensis plastomes and comparative analysis of Firmiana: insight into its phylogeny and evolution

BACKGROUND

Firmiana danxiaensis is a critically endangered and ecologically important tree currently only found in four locations in Danxia or Karst habitats in northern Guangdong Province, China. The specialized habitat preference makes it an ideal model species for study of adaptive evolution. Meanwhile, the phylogenetic relationships of F. danxiaensis in four locations under two landforms are unclear. Therefore, we sequenced its complete chloroplast (cp.) genomes and conducted comprehensive interspecific and intrageneric plastome studies.

RESULTS

The F. danxiaensis plastomes in four locations showed a typical quadripartite and circular structure that ranged from 160,832 to 161,206 bp in size, with 112 unique genes encoded. Comparative genomics showed that the plastomes of F. danxiaensis were relatively conserved with high similarity of genome organization, gene number, GC content and SSRs. While the genomes revealed higher biased codon preferences in Karst habitat than those in Danxia habitats. Eighteen and 11 divergent hotpots were identified at interspecific and intrageneric levels for species identification and further phylogenetic studies. Seven genes (clpP, accD, ccsA, ndhH, rpl20, rpoC2, and rps4) were under positive selection and may be related to adaptation. Phylogenetic analysis revealed that F. danxiaensis is sister to F. major and F. simplex. However, the interspecific relationships are not consistent with the habitat types.

CONCLUSIONS

The characteristics and interspecific relationship of F. danxiaensis plastomes provide new insights into further integration of geographical factors, environmental factors, and genetic variations on the genomic study of F. danxiaensis. Together, our study will contribute to the study of species identification, population genetics, and conservation biology of F. danxiaensis.

Genome-Wide Analysis of MYB Genes in Primulina eburnea (Hance) and Identification of Members in Response to Drought Stress

Due to periodic water deficiency in karst environments, Primulina eburnea experiences sporadic drought stress in its habitat. Despite being one of the largest gene families and functionally diverse in terms of plant growth and development, MYB transcription factors in P. eburnea have not been studied. Here, a total of 230 MYB genes were identified in P. eburnea, including 67 1R-MYB, 155 R2R3-MYB, six 3R-MYB, and two 4R-MYB genes. The R2R3-type PebMYB genes could be classified into 16 subgroups, while the remaining PebMYB genes (1R-MYB, 3R-MYB, and 4R-MYB genes) were divided into 10 subgroups. Notably, the results of the phylogenetic analysis were further supported by the motif and gene structure analysis, which showed that individuals in the same subgroup had comparable motif and structure organization. Additionally, gene duplication and synteny analyses were performed to better understand the evolution of PebMYB genes, and 291 pairs of segmental duplicated genes were found. Moreover, RNA-seq analysis revealed that the PebMYB genes could be divided into five groups based on their expression characteristics. Furthermore, 11 PebMYB genes that may be involved in drought stress response were identified through comparative analysis with Arabidopsis thaliana. Notably, seven of these genes (PebMYB3, PebMYB13, PebMYB17, PebMYB51, PebMYB142, PebMYB69, and PebMYB95) exhibited significant differences in expression between the control and drought stress treatments, suggesting that they may play important roles in drought stress response. These findings clarified the characteristics of the MYB gene family in P. eburnea, augmenting our comprehension of their potential roles in drought stress adaptation.

An integrative approach to species delimitation sinks three Chinese limestone karst Elatostema (Urticaceae) species

Elatostema is recognized as a taxonomically difficult group due to the reduced nature of the tiny flowers and inflorescences, also the large number of species (ca 650 to 700). Different opinions on morphological species delimitation have resulted in instability, which is problematic in such a speciose group. In this paper, the taxonomic status of three putative species, E.robustipes, E.scaposum, E.conduplicatum and their hypothetical closest relatives, was revised using morphological and molecular observations. Morphological comparison suggested high similarity between E.robustipes & E.retrohirtum, E.scaposum & E.oblongifolium, E.conduplicatum & E.coriaceifolium, respectively. Phylogenetic analyses of four universal DNA barcodes (ITS, trnH-psbA, matK and rbcL) suggested that each species pair represents a single evolutionary lineage. Taking these two findings together, we propose E.robustipes to be a synonym of E.retrohirtum, E.scaposum a synonym of E.oblongifolium, and E.conduplicatum a synonym of E.coriaceifolium. Our results recover the number, shape and size of the bracts and bracteoles to be relatively stable characters, and the disposition of the male inflorescences on modified stems to be an unstable character, unsuitable for species delimitation in Elatostema.

Primulinapingnanensis, a new species of Gesneriaceae from Guangxi, China

A new species of Primulina, P.pingnanensis, from the Guangxi Zhuangzu Autonomous Region, China, is described and illustrated here. It is morphologically similar to P.orthandra but has significant differences in the bracts, corolla tube and lobes shape, as well as in the indumentum of the outer surface of the corolla, the filaments, the staminodes and the anthers. Colorful photographs and essential information of this new taxon are also provided, including detailed taxonomic description, distribution, habitat, the comparison table, and the IUCN conservation status. We also discuss a validation of new combination P.crassifolia and Chiritacrassifolia.

Genomic Insights into Adaptation to Karst Limestone and Incipient Speciation in East Asian Platycarya spp. (Juglandaceae)

When challenged by similar environmental conditions, phylogenetically distant taxa often independently evolve similar traits (convergent evolution). Meanwhile, adaptation to extreme habitats might lead to divergence between taxa that are otherwise closely related. These processes have long existed in the conceptual sphere, yet molecular evidence, especially for woody perennials, is scarce. The karst endemic Platycarya longipes and its only congeneric species, Platycarya strobilacea, which is widely distributed in the mountains in East Asia, provide an ideal model for examining the molecular basis of both convergent evolution and speciation. Using chromosome-level genome assemblies of both species, and whole-genome resequencing data from 207 individuals spanning their entire distribution range, we demonstrate that P. longipes and P. strobilacea form two species-specific clades, which diverged around 2.09 million years ago. We find an excess of genomic regions exhibiting extreme interspecific differentiation, potentially due to long-term selection in P. longipes, likely contributing to the incipient speciation of the genus Platycarya. Interestingly, our results unveil underlying karst adaptation in both copies of the calcium influx channel gene TPC1 in P. longipes. TPC1 has previously been identified as a selective target in certain karst-endemic herbs, indicating a convergent adaptation to high calcium stress among karst-endemic species. Our study reveals the genic convergence of TPC1 among karst endemics and the driving forces underneath the incipient speciation of the two Platycarya lineages.

Genetic and ecophysiological evidence that hybridization facilitated lineage diversification in yellow Camellia (Theaceae) species: a case study of natural hybridization between C. micrantha and C. flavida

BACKGROUND

Hybridization is generally considered an important creative evolutionary force, yet this evolutionary process is still poorly characterized in karst plants. In this study, we focus on natural hybridization in yellow Camellia species, a group of habitat specialists confined to karst/non-karst habitats in southwestern China.

RESULTS

Based on population genome data obtain from double digest restriction-site associated DNA (ddRAD) sequencing, we found evidence for natural hybridization and introgression between C. micrantha and C. flavida, and specifically confirmed their hybrid population, C. "ptilosperma". Ecophysiological results suggested that extreme hydraulic traits were fixed in C. "ptilosperma", these being consistent with its distinct ecological niche, which lies outside its parental ranges.

CONCLUSION

The identified hybridization event is expected to have played a role in generating novel variation during, in which the hybrid population displays different phenological characteristics and novel ecophysiological traits associated with the colonization of a new niche in limestone karst.

Changes in root chemical diversity along an elevation gradient of Changbai Mountain, China

Root chemical traits play a critical role in plant resource use strategies and ecosystem nutrient cycling; however, the chemical diversity of multiple elements of fine root and community chemical assembly belowground are poorly understood. Here, we measured 13 elements (C, N, K, Ca, Mg, S, P, Al, Fe, Na, Mn, Zn, and Cu) in the fine roots of 204 plant species along elevational transect from 540 to 2357 m of Changbai Mountain, China to explore the variation, diversity, and community assembly of root chemical traits. At the species level, the concentrations of macronutrients (N, K, Ca, Mg, S, and P) decreased, whereas the trace metals (Fe, Mn, and Zn) increased with elevation. Root chemical traits at the community level systematically shifted along elevational gradients showing a pattern similar to that at the species level, which were mainly influenced by climate and soil rather than species diversity. In general, the interactions of climate and soil were the main drivers of root chemical assembly for woody layers, whereas soil factors played significant role for root chemical assembly for herb layer. The chemical assembly of rock-derived element P was mainly driven by soil factors. Meanwhile, root chemical diversities were mainly regulated by species diversity, the interactions of climate and soil, and soil factors in the tree, shrub, and herb layers, respectively. A better understanding of plant root chemical diversity and community chemical assembly will help to reveal the role of chemical traits in ecosystem functioning.

Recent and ancient evolutionary events shaped plant elemental composition of edaphic endemics: a phylogeny-wide analysis of Iberian gypsum plants

The analysis of plant elemental composition and the underlying factors affecting its variation are a current hot topic in ecology. Ecological adaptation to atypical soils may shift plant elemental composition. However, no previous studies have evaluated its relevance against other factors such as phylogeny, climate or individual soil conditions. We evaluated the effect of the phylogeny, environment (climate, soil), and affinity to gypsum soils on the elemental composition of 83 taxa typical of Iberian gypsum ecosystems. We used a new statistical procedure (multiple phylogenetic variance decomposition, MPVD) to decompose total explained variance by different factors across all nodes in the phylogenetic tree of target species (covering 120 million years of Angiosperm evolution). Our results highlight the relevance of phylogeny on the elemental composition of plants both at early (with the development of key preadaptive traits) and recent divergence times (diversification of the Iberian gypsum flora concurrent with Iberian gypsum deposit accumulation). Despite the predominant phylogenetic effect, plant adaptation to gypsum soils had a strong impact on the elemental composition of plants, particularly on sulphur concentrations, while climate and soil effects were smaller. Accordingly, we detected a convergent evolution of gypsum specialists from different lineages on increased sulphur and magnesium foliar concentrations.

Immigration dynamics of tropical and subtropical Southeast Asian limestone karst floras

Ex situ origins and dispersal of taxa have played important roles in the assembly of island-like biodiversity hotspots. Insular limestone karsts in Southeast Asia are hotspots of biodiversity and endemism, but the immigration processes of their unique floras are still poorly known. Here, we used Gesneriaceae as a proxy to investigate the immigration dynamics of tropical and subtropical Southeast Asian karst floras. We present the most comprehensive phylogenetic analysis of the Old World gesneriads to date based on twelve loci. By estimating divergence times and reconstructing ancestral states (habitat, soil type and range), we found that immigration into subtropical Southeast Asian karst floras first occurred in the Early Miocene, with two peaks in the Early-Middle Miocene and the Pliocene-Early Pleistocene, whereas immigration into tropical Southeast Asian karsts initiated in the Late Eocene, with two peaks in the Late Oligocene and the Late Miocene. We also discover that Southeast Asian karst biodiversity comprises immigrant pre-adapted lineages and descendants from local acid soil ancestors, although niche shift from acid soil to karst in tropical Southeast Asian islands was lacking. This study advances our understanding of the historical assembly of Southeast Asian karst floras.

Global patterns of leaf construction traits and their covariation along climate and soil environmental gradients

Leaf functional traits and their covariation underlie plant ecological adaptations along environmental gradients, but there is limited information on the global covariation patterns of key leaf construction traits. To explore how leaf construction traits co-vary across diverse climate and soil environmental conditions, we compiled a global dataset including cell wall mass per unit leaf mass (CW_mass ), leaf carbon (C) and calcium (Ca) concentrations, and specific leaf area (SLA) for 2348 angiosperm species from 340 sites world-wide. Our results demonstrated negative correlations between leaf C and Ca concentrations and between leaf C and SLA across diverse nongraminoid angiosperms. Leaf C concentration increased with increasing mean annual temperature (MAT) and mean annual precipitation (MAP) and with decreasing soil pH and calcium carbonate (CaCO₃ ) concentration, whereas leaf Ca concentration and SLA exhibited the opposite responses to these environmental variables. The covariations of leaf Ca-C and of leaf SLA-C were stronger in habitats with lower MAT and MAP, and/or higher soil CaCO₃ content. This global-scale analysis demonstrates that the leaf C and Ca concentrations and SLA together govern the C and biomass investment strategies in leaves of nongraminoids. We conclude that environmental conditions strongly shape leaf construction traits and their covariation patterns.

Gene flow, linked selection, and divergent sorting of ancient polymorphism shape genomic divergence landscape in a group of edaphic specialists

Interpreting the formation of genomic variation landscape, especially genomic regions with elevated differentiation (i.e. islands), is fundamental to a better understanding of the genomic consequences of adaptation and speciation. Edaphic islands provide excellent systems for understanding the interplay of gene flow and selection in driving population divergence and speciation. However, discerning the relative contribution of these factors that modify patterns of genomic variation remains difficult. We analysed 132 genomes from five recently divergent species in Primulina genus, with four species distributed in Karst limestone habitats and the fifth one growing in Danxia habitats. We demonstrated that both gene flow and linked selection have contributed to genome-wide variation landscape, where genomic regions with elevated differentiation (i.e., islands) were largely derived by divergent sorting of ancient polymorphism. Specifically, we identified several lineage-specific genomic islands that might have facilitated adaptation of P. suichuanensis to Danxia habitats. Our study is amongst the first cases disentangling evolutionary processes that shape genomic variation of plant specialists, and demonstrates the important role of ancient polymorphism in the formation of genomic islands that potentially mediate adaptation and speciation of endemic plants in special soil habitats.

Plant adaptability in karst regions

Karst ecosystems are formed by dissolution of soluble rocks, usually with conspicuous landscape features, such as sharp peaks, steep slopes and deep valleys. The plants in karst regions develop special adaptability. Here, we reviewed the research progresses on plant adaptability in karst regions, including drought, high temperature and light, high-calcium stresses responses and the strategies of water utilization for plants, soil nutrients impact, human interference and geographical traits on karst plants. Drought, high temperature and light change their physiological and morphological structures to adapt to karst environments. High-calcium and soil nutrients can transfer surplus nutrients to special parts of plants to avoid damage of high nutrient concentration. Therefore, karst plants can make better use of limited water. Human interference also affects geographical distribution of karst plants and their growing environment. All of these aspects may be analyzed to provide guidance and suggestions for related research on plant adaptability mechanisms.

Soil properties mediate ecosystem intrinsic water use efficiency and stomatal conductance via taxonomic diversity and leaf economic spectrum

The interactions between plants and soils lead to complex feedbacks that regulate intrinsic water use efficiency (iWUE) and stomatal conductance (gs) at ecosystem level and reflect water constraints on plant productivity. However, the relationships among soil properties, biodiversity, and leaf functional traits contributing to the variability in ecosystem iWUE and gs remain largely unknown. To elucidate these relationships, we used principal component analysis to reduce soil properties to a fertility spectrum and a limiting-resource spectrum across grassland, and early-, mid- and late-successional forests in a karst catchment. Leaf functional traits at community level were calculated based on leaf biomass, and were reduced to an economic spectrum and a limiting-resource spectrum. Leaf carbon (δ¹³C) and oxygen (δ¹⁸O) stable isotopes at community levels were used as proxies for ecosystem iWUE and g_s. The effects of soil properties, biodiversity (taxonomic, functional and phylogenetic diversity) and leaf traits on δ¹³C and δ¹⁸O were evaluated using structural equation models. Our results showed that variability in ecosystem iWUE and gs was determined overwhelmingly by indirect effects of soil properties via two different pathways: the soil fertility spectrum, determining the number of coexisting species (taxonomic diversity) and turnover of species (leaf economic spectrum), and the soil limiting-resource spectrum, shaping the specific phylogenetic lineages (phylogenic diversity). In addition, δ¹³C and δ¹⁸O were constrained by the interactive effects of leaf economic spectrum, and taxonomic and phylogenic diversity; total effects of biodiversity on δ¹³C and δ¹⁸O were larger than those of leaf economic spectrum. Our study highlighted the critical role of the evaluating interaction relationships between leaf functional traits, biodiversity metrics and soil properties in understanding the mechanisms of ecosystem function responding to environmental change.

Increase in leaf organic acids to enhance adaptability of dominant plant species in karst habitats

Estimation of leaf nutrient composition of dominant plant species from contrasting habitats (i.e., karst and nonkarst forests) provides an opportunity to understand how plants are adapted to karst habitats from the perspective of leaf traits. Here, we measured leaf traits-specific leaf area (SLA), concentrations of total carbon ([TC]), nitrogen ([TN]), phosphorus ([TP]), calcium ([Ca]), magnesium ([Mg]), manganese ([Mn]), minerals ([Min]), soluble sugars, soluble phenolics, lipids, and organic acids ([OA])-and calculated water-use efficiency (WUE), construction costs (CC), and N/P ratios, and searched for correlations between these traits of 18 abundant plant species in karst and nonkarst forests in southwestern China. Variation in leaf traits within and across the abundant species was both divergent and convergent. Leaf [TC], [Ca], [Min], [OA], and CC were habitat-dependent, while the others were not habitat- but species-specific. The correlations among [TN], [TP], SLA, [TC], CC, [Min], WUE, [OA], and CC were habitat-independent, and inherently associated with plant growth and carbon allocation; those between [CC] and [Lip], between [Ca] and [Mg], and between [Mg] and [WUE] were habitat-dependent. Habitat significantly affected leaf [Ca] and thus indirectly affected leaf [OA], [Min], and CC. Our results indicate that plants may regulate leaf [Ca] to moderate levels via adjusting leaf [OA] under both high and low soil Ca availability, and offer new insights into the abundance of common plant species in contrasting habitats.

Relationships Between Leaf Carbon and Macronutrients Across Woody Species and Forest Ecosystems Highlight How Carbon Is Allocated to Leaf Structural Function

Stoichiometry of leaf macronutrients can provide insight into the tradeoffs between leaf structural and metabolic investments. Structural carbon (C) in cell walls is contained in lignin and polysaccharides (cellulose, hemicellulose, and pectins). Much of leaf calcium (Ca) and a fraction of magnesium (Mg) were further bounded with cell wall pectins. The macronutrients phosphorus (P), potassium (K), and nitrogen (N) are primarily involved in cell metabolic functions. There is limited information on the functional interrelations among leaf C and macronutrients, and the functional dimensions characterizing the leaf structural and metabolic tradeoffs are not widely appreciated. We investigated the relationships between leaf C and macronutrient (N, P, K, Ca, Mg) concentrations in two widespread broad-leaved deciduous woody species Quercus wutaishanica (90 individuals) and Betula platyphylla (47 individuals), and further tested the generality of the observed relationships in 222 woody eudicots from 15 forest ecosystems. In a subsample of 20 broad-leaved species, we also analyzed the relationships among C, Ca, lignin, and pectin concentrations in leaf cell walls. We found a significant leaf C-Ca tradeoff operating within and across species and across ecosystems. This basic relationship was explained by variations in the share of cell wall lignin and pectin investments at the cell scale. The C-Ca tradeoffs were mainly driven by soil pH and mean annual temperature and precipitation, suggesting that leaves were more economically built with less C and more Ca as soil pH increased and at lower temperature and lower precipitation. However, we did not detect consistent patterns among C-N, and C-Mg at different levels of biological organization, suggesting substantial plasticity in N and Mg distribution among cell organelles and cell protoplast and cell wall. We observed two major axes of macronutrient differentiation: the cell-wall structural axis consisting of protein-free C and Ca and the protoplasm metabolic axis consisting of P and K, underscoring the decoupling of structural and metabolic elements inherently linked with cell wall from protoplasm investment strategies. We conclude that the tradeoffs between leaf C and Ca highlight how carbon is allocated to leaf structural function and suggest that this might indicate biogeochemical niche differentiation of species.

A Combined Morphological and Molecular Evolutionary Analysis of Karst-Environment Adaptation for the Genus Urophysa (Ranunculaceae)

The karst environment is characterized by low soil water content, periodic water deficiency, and poor nutrient availability, which provides an ideal natural laboratory for studying the adaptive evolution of its inhabitants. However, how species adapt to such a special karst environment remains poorly understood. Here, transcriptome sequences of two Urophysa species (Urophysa rockii and Urophysa henryi), which are Chinese endemics with karst-specific distribution, and allied species in Semiaquilegia and Aquilegia (living in non-karst habitat) were collected. Single-copy genes (SCGs) were extracted to perform the phylogenetic analysis using concatenation and coalescent methods. Positively selected genes (PSGs) and clusters of paralogous genes (Mul_genes) were detected and subsequently used to conduct gene function annotation. We filtered 2,271 SCGs and the coalescent analysis revealed that 1,930 SCGs shared the same tree topology, which was consistent with the topology detected from the concatenated tree. Total of 335 PSGs and 243 Mul_genes were detected, and many were enriched in stress and stimulus resistance, transmembrane transport, cellular ion homeostasis, calcium ion transport, calcium signaling regulation, and water retention. Both molecular and morphological evidences indicated that Urophysa species evolved complex strategies for adapting to hostile karst environments. Our findings will contribute to a new understanding of genetic and phenotypic adaptive mechanisms of karst adaptation in plants.

Leaf trait variation in species-rich tropical Andean forests

Screening species-rich communities for the variation in functional traits along environmental gradients may help understanding the abiotic drivers of plant performance in a mechanistic way. We investigated tree leaf trait variation along an elevation gradient (1000-3000 m) in highly diverse neotropical montane forests to test the hypothesis that elevational trait change reflects a trend toward more conservative resource use strategies at higher elevations, with interspecific trait variation decreasing and trait integration increasing due to environmental filtering. Analysis of trait variance partitioning across the 52 tree species revealed for most traits a dominant influence of phylogeny, except for SLA, leaf thickness and foliar Ca, where elevation was most influential. The community-level means of SLA, foliar N and Ca, and foliar N/P ratio decreased with elevation, while leaf thickness and toughness increased. The contribution of intraspecific variation was substantial at the community level in most traits, yet smaller than the interspecific component. Both within-species and between-species trait variation did not change systematically with elevation. High phylogenetic diversity, together with small-scale edaphic heterogeneity, cause large interspecific leaf trait variation in these hyper-diverse Andean forests. Trait network analysis revealed increasing leaf trait integration with elevation, suggesting stronger environmental filtering at colder and nutrient-poorer sites.

Genomic Mechanisms of Physiological and Morphological Adaptations of Limestone Langurs to Karst Habitats

Knowledge of the physiological and morphological evolution and adaptation of nonhuman primates is critical to understand hominin origins, physiological ecology, morphological evolution, and applications in biomedicine. Particularly, limestone langurs represent a direct example of adaptations to the challenges of exploiting a high calcium and harsh environment. Here, we report a de novo genome assembly (Tfra_2.0) of a male François's langur (Trachypithecus francoisi) with contig N50 of 16.3 Mb and resequencing data of 23 individuals representing five limestone and four forest langur species. Comparative genomics reveals evidence for functional evolution in genes and gene families related to calcium signaling in the limestone langur genome, probably as an adaptation to naturally occurring high calcium levels present in water and plant resources in karst habitats. The genomic and functional analyses suggest that a single point mutation (Lys1905Arg) in the α1c subunit of the L-type voltage-gated calcium channel Cav1.2 (CACNA1C) attenuates the inward calcium current into the cells in vitro. Population genomic analyses and RNA-sequencing indicate that EDNRB is less expressed in white tail hair follicles of the white-headed langur (T. leucocephalus) compared with the black-colored François's langur and hence might be responsible for species-specific differences in body coloration. Our findings contribute to a new understanding of gene-environment interactions and physiomorphological adaptative mechanisms in ecologically specialized primate taxa.

Leaf Multi-Element Network Reveals the Change of Species Dominance Under Nitrogen Deposition

Elements are important functional traits reflecting plant response to climate change. Multiple elements work jointly in plant physiology. Although a large number of studies have focused on the variation and allocation of multiple elements in plants, it remains unclear how these elements co-vary to adapt to environmental change. We proposed a novel concept of the multi-element network including the mutual effects between element concentrations to more effectively explore the alterations in response to long-term nitrogen (N) deposition. Leaf multi-element networks were constructed with 18 elements (i.e., six macronutrients, six micronutrients, and six trace elements) in this study. Multi-element networks were species-specific, being effectively discriminated irrespective of N deposition level. Different sensitive elements and interactions to N addition were found in different species, mainly concentrating on N, Ca, Mg, Mn, Li, Sr, Ba, and their related stoichiometry. Interestingly, high plasticity of multi-element network increased or maintained relative aboveground biomass (species dominance) in community under simulated N deposition, which developed the multi-element network hypothesis. In summary, multi-element networks provide a novel approach for exploring the adaptation strategies of plants and to better predict the change of species dominance under altering nutrient availability or environmental stress associated with future global climate change.

The genome of a cave plant, Primulina huaijiensis, provides insights into adaptation to limestone karst habitats

Although whole genome duplication (WGD) has been suggested to facilitate adaptive evolution and diversification, the role of specific WGD events in promoting diversification and adaptation in angiosperms remains poorly understood. Primulina, a species-rich genus with > 180 species associated with limestone karst habitat, constitutes an ideal system for studying the impact of WGD events on speciation and evolutionary adaptation. We sequenced and assembled a chromosome-level genome of the cave-dwelling species P. huaijiensis to study gene family expansion and gene retention following WGDs. We provide evidence that P. huaijiensis has undergone two WGDs since the γ triplication event shared by all eudicots. In addition to a WGD shared by almost all Lamiales (L event), we identified a lineage-specific WGD (D event) that occurred in the early Miocene around 20.6-24.2 Myr ago and that is shared by almost the entire subtribe Didymocarpinae. We found that gene retentions following the D event led to gene family proliferation (e.g. WRKYs) that probably facilitated adaptation to the high salinity and drought stress in limestone karst. Our study highlights the role of lineage-specific WGD in species diversification and adaptation of plants from special habitats.

Primulina hochiensis var. ochroleuca (Gesneriaceae), a new variety from a limestone area of Guangxi, China, and errata on five new species of Primulina

Primulina hochiensis var. ochroleuca, a new variety from a limestone hill of karst areas, Guangxi, China is described with color photographs. It resembles P. hochiensis var. hochiensis, P. hochiensis var. ovata and P. hochiensis var. rosulata, but can be easily distinguished by a combination of characteristics, especially by its corolla color. We found only one population with approx. 3000 mature individuals at the type locality. This variety is provisionally assessed as vulnerable [VU C1] using IUCN criteria.

A Preliminary Evaluation of The Karst Flora of Brazil Using Collections Data

Karst is defined as landscapes that are underlain by soluble rock in which there is appreciable water movement arising from a combination of high rock solubility and well-developed secondary (fracture) porosity. Karsts occupy approximately 20% of the planet’s dry ice-free land and are of great socioeconomic importance, as they supply water to up to 25% of the world’s population and represent landscapes of cultural and touristic importance. In Southeast Asia karst is associated with high species-richness and endemism in plants and seen as priority areas for the conservation of biodiversity. There has been little research into the floras associated with karst in South America, most of which occurs in Brazil. We therefore sought to evaluate the importance of Brazilian karst with respect to its species-richness and endemism. We sought to do so using curated plant specimen data in the Botanical Information and Ecology Network (BIEN) dataset. We show that, except for Amazonia, the BIEN dataset is representative of the Brazilian flora with respect to the total number of species and overall patterns of species richness. We found that karst is under-sampled, as is the case for much of Brazil. We also found that whilst karst represent an important source of plant diversity for Brazil, including populations of approximately 1/3 of the Brazilian flora, it is not significantly more species-rich or richer in small-range and endemic species than surrounding landscapes. Similarly, whilst important for conservation, comprising populations of 26.5–37.4% of all Brazilian species evaluated as of conservation concern by International Union for Nature Conservation (IUCN), karst is no more so than the surrounding areas. Whilst experimental error, including map resolution and the precision and accuracy of point data may have under-estimated the species-richness of Brazilian karst, it likely represents an important biodiversity resource for Brazil and one that can play a valuable role in conservation. Our findings are in sharp contrast to those for Southeast Asia where karst represents a more important source of species-richness and endemism. We also show that although BIEN represents a comprehensive and curated source of point data, discrepancies in the application of names compared to current more comprehensive taxonomic backbones, can have profound impacts on estimates of species-richness, distribution ranges and estimates of endemism.

Primulina serrulata (Gesneriaceae), a new species from southeastern Guizhou, China

Primulina serrulata R.B.Zhang & F. Wen, a new species from a limestone area in southeastern Guizhou, China, is described and illustrated here. The new species is morphologically related to P. fimbrisepala (Hand.-Mazz.) Y.Z.Wang. We examined the morphological differences between these congeners and provide illustrations and photographs of this new species in this paper.

The complete chloroplast genome of Primula tsiangii W. W. Smith (Primulaceae): a karst endemic primrose in Southwest China

Primula tsiangii W. W. Smith is a karst endemic species in Southwest China. Here, we report the complete chloroplast genome of P. tsiangii. The chloroplast (cp) genome was determined to be 153,281 bp and the GC contents was 37.2%. The sequence includes a large single copy (LSC) region of 84,524 bp, a small single copy (SSC) region of 17,357 bp, and two separated inverted regions of 25,700 bp each. It contains 127 unique genes, including 82 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. This is the first report of cp genomes for karst endemic primrose, and will be useful for the genetic diversity of karst plants.

The rediscovery and delimitation of Elatostemasetulosum W.T.Wang (Urticaceae)

Of the 280 species of Elatostema documented in China, 189 are known only from a single collection. Elatostemasetulosum is one such species, having been known only from the type collection for nearly half a century, until recent field investigations in Guangxi. Due to its morphological similarity to E.huanjiangense and E.tetracephalum, we undertook a critical review of all three species using morphological and molecular evidence. Our results suggest that all three names refer to the same species, which based on priority should be known as Elatostemasetulosum. We recognize E.huanjiangense and E.tetracephalum as synonyms. A distribution map of E.setulosum and the extinction risk according to the IUCN criteria is provided. After recircumscription, the taxon must be considered as Least Concern (LC).

Family-level leaf nitrogen and phosphorus stoichiometry of global terrestrial plants

Leaf nitrogen (N) and phosphorus (P) concentrations are critical for photosynthesis, growth, reproduction and other ecological processes of plants. Previous studies on large-scale biogeographic patterns of leaf N and P stoichiometric relationships were mostly conducted using data pooled across taxa, while family/genus-level analyses are rarely reported. Here, we examined global patterns of family-specific leaf N and P stoichiometry using a global data set of 12,716 paired leaf N and P records which includes 204 families, 1,305 genera, and 3,420 species. After determining the minimum size of samples (i.e., 35 records), we analyzed leaf N and P concentrations, N:P ratios and N∼P scaling relationships of plants for 62 families with 11,440 records. The numeric values of leaf N and P stoichiometry varied significantly across families and showed diverse trends along gradients of mean annual temperature (MAT) and mean annual precipitation (MAP). The leaf N and P concentrations and N:P ratios of 62 families ranged from 6.11 to 30.30 mg g^-1, 0.27 to 2.17 mg g^-1, and 10.20 to 35.40, respectively. Approximately 1/3-1/2 of the families (22-35 of 62) showed a decrease in leaf N and P concentrations and N:P ratios with increasing MAT or MAP, while the remainder either did not show a significant trend or presented the opposite pattern. Family-specific leaf N∼P scaling exponents did not converge to a certain empirical value, with a range of 0.307-0.991 for 54 out of 62 families which indicated a significant N∼P scaling relationship. Our results for the first time revealed large variation in the family-level leaf N and P stoichiometry of global terrestrial plants and that the stoichiometric relationships for at least one-third of the families were not consistent with the global trends reported previously. The numeric values of the family-specific leaf N and P stoichiometry documented in the current study provide critical synthetic parameters for biogeographic modeling and for further studies on the physiological and ecological mechanisms underlying the nutrient use strategies of plants from different phylogenetic taxa.

Effects of nitrogen and phosphorus supply on stoichiometry of six elements in leaves of Arabidopsis thaliana

BACKGROUND AND AIMS

Plant elemental composition is of fundamental importance for plant growth and metabolic functions. However, knowledge of how multi-elemental stoichiometry responds to varying nitrogen (N) and phosphorus (P) availabilities remains limited.

METHODS

We conducted experimental manipulations with nine repeat experiments to investigate the effects of N and P supply on the concentrations and variability of six elements, carbon (C), N, P, potassium (K), calcium (Ca) and magnesium (Mg), in leaves of Arabidopsis thaliana.

KEY RESULTS

N supply increased the concentrations of N, K and Mg, decreased the concentration of P, but exerted little influence on the concentrations of C and Ca in green leaves. P supply increased the concentrations of P and Ca, decreased the concentration of C, initially increased and then decreased the concentration of K, but showed little influence on the concentrations of N and Mg in green leaves. Multivariate patterns among the concentrations of these six elements in green leaves was influenced by the type of nutrient supply (i.e. N or P). Elemental variability decreased with increasing elemental concentrations in green leaves at the intraspecific level, supporting the Stability of Limiting Elements Hypothesis that was originally proposed from a meta-analysis of pooled data across species or communities. Compared with green leaves, the senesced leaves showed greater variability in C, N, P, K and Mg concentrations but lower variability in Ca concentration.

CONCLUSIONS

N and P supplies exerted differential influences on the concentrations of C, N, P, K, Ca and Mg in green leaves. The specific C content should be considered when assessing C cycling under global nutrient changes. Stage-dependent patterns of leaf stoichiometric homeostasis differed among elements with various chemical characteristics. These findings can help to improve our understanding of plant eco-physiological responses and acclimation under global nutrient changes from the stoichiometric perspective of multiple elements.

Genetic architecture of quantitative flower and leaf traits in a pair of sympatric sister species of Primulina

Flowers and leaves each represent suites of functionally interrelated traits that are often involved in species divergence and local adaptation. However, a major unresolved issue is how the individual component traits that make up a complex trait such as a flower evolve in a coordinated fashion to retain a high degree of functionality. We use a quantitative trait loci (QTL) approach to elucidate the genetic architecture of divergence in flower and leaf traits between the sister species Primulina depressa and Primulina danxiaensis, which grow sympatrically but in contrasting microhabitats. We found that flower traits were controlled by multiple QTL of small effect, while leaf physiological and morphological traits tended to be controlled by QTL of larger effect. The observed floral integration, manifested by a high degree overlap in both individual trait QTL and QTL for principal component scores (PCA QTL), may have been critical for evolutionary divergence of floral morphology in relation to their pollinators. This overlap suggests that direct selection on only one or a few of the component traits could have caused substantial divergence in other floral traits due to genetic correlations, while the low QTL overlap between floral and vegetative traits suggests that these trait suites are genetically unlinked and can evolve independently in response to different selective pressures corresponding to their distinct functions.

Variation in the angiosperm ionome

The ionome is defined as the elemental composition of a subcellular structure, cell, tissue, organ or organism. The subset of the ionome comprising mineral nutrients is termed the functional ionome. A 'standard functional ionome' of leaves of an 'average' angiosperm, defined as the nutrient composition of leaves when growth is not limited by mineral nutrients, is presented and can be used to compare the effects of environment and genetics on plant nutrition. The leaf ionome of a plant is influenced by interactions between its environment and genetics. Examples of the effects of the environment on the leaf ionome are presented and the consequences of nutrient deficiencies on the leaf ionome are described. The physiological reasons for (1) allometric relationships between leaf nitrogen and phosphorus concentrations and (2) linear relationships between leaf calcium and magnesium concentrations are explained. It is noted that strong phylogenetic effects on the mineral composition of leaves of angiosperm species are observed even when sampled from diverse environments. The evolutionary origins of traits including (1) the small calcium concentrations of Poales leaves, (2) the large magnesium concentrations of Caryophyllales leaves and (3) the large sulphur concentrations of Brassicales leaves are traced using phylogenetic relationships among angiosperm orders, families and genera. The rare evolution of hyperaccumulation of toxic elements in leaves of angiosperms is also described. Consequences of variation in the leaf ionome for ecology, mineral cycling in the environment, strategies for phytoremediation of contaminated land, sustainable agriculture and the nutrition of livestock and humans are discussed.

Primulina malipoensis (Gesneriaceae), a new species from Sino-Vietnamese border area

Primulina malipoensis, a new species from limestone areas around the Sino-Vietnamese border, is described and illustrated. This new species is morphologically similar to P. maguanensis and P. lungzhouensis, but obviously differs from the latter two species by its pale greenish-yellow flowers (vs. purple, with different colour patterns). The phylogenetic affinity, illustration and photographs of this new species are provided in this paper.

Genomic Signature of Adaptive Divergence despite Strong Nonadaptive Forces on Edaphic Islands: A Case Study of Primulina juliae

Both genetic drift and divergent selection are expected to be strong evolutionary forces driving population differentiation on edaphic habitat islands. However, the relative contribution of genetic drift and divergent selection to population divergence has rarely been tested simultaneously. In this study, restriction-site associated DNA-based population genomic analyses were applied to assess the relative importance of drift and divergent selection on population divergence of Primulina juliae, an edaphic specialist from southern China. All populations were found with low standing genetic variation, small effective population size (NE), and signatures of bottlenecks. Populations with the lowest genetic variation were most genetically differentiated from other populations and the extent of genetic drift increased with geographic distance from other populations. Together with evidence of isolation by distance, these results support neutral drift as a critical evolutionary driver. Nonetheless, redundancy analysis revealed that genomic variation is significantly associated with both edaphic habitats and climatic factors independently of spatial effects. Moreover, more genomic variation was explained by environmental factors than by geographic variables, suggesting that local adaptation might have played an important role in driving population divergence. Finally, outlier tests and environment association analyses identified 31 single-nucleotide polymorphisms as candidates for adaptive divergence. Among these candidates, 26 single-nucleotide polymorphisms occur in/near genes that potentially play a role in adaptation to edaphic specialization. This study has important implications that improve our understanding of the joint roles of genetic drift and adaptation in generating population divergence and diversity of edaphic specialists.

The complete chloroplast genome of Primulina and two novel strategies for development of high polymorphic loci for population genetic and phylogenetic studies

Background

Primulina Hance is an emerging model for studying evolutionary divergence, adaptation and speciation of the karst flora. However, phylogenetic relationships within the genus have not been resolved due to low variation detected in the cpDNA regions. Chloroplast genomes can provide important information for phylogenetic and population genetic studies. Recent advances in next-generation sequencing (NGS) techniques greatly facilitate sequencing whole chloroplast genomes for multiple individuals. Consequently, novel strategies for development of highly polymorphic loci for population genetic and phylogenetic studies based on NGS data are needed.

Methods

For development of high polymorphic loci for population genetic and phylogenetic studies, two novel strategies are proposed here. The first protocol develops lineage-specific highly variable markers from the true high variation regions (Con_Seas) across whole cp genomes, instead of traditional noncoding regions. The pipeline has been integrated into a single perl script, and named "Con_Sea_Identification_and_PIC_Calculation". The second method assembles chloroplast fragments (poTs) and sub-super-marker (CpContigs) through our "SACRing" pipeline. This approach can fundamentally alter the strategies used in phylogenetic and population genetic studies based on cp markers, facilitating a transition from traditional Sanger sequencing to RAD-Seq. Both of these scripts are available at https://github.com/scbgfengchao/.

Results

Three complete Primulina chloroplast genomes were assembled from genome survey data, and then two novel strategies were developed to yield highly polymorphic markers. For experimental evaluation of the first protocol, a set of Primulina species were used for PCR amplification. The results showed that these newly developed markers are more variable than traditional ones, and seem to be a better choice for phylogenetic and population studies in Primulina. The second method was also successfully applied in population genetic studies of 21 individuals from three natural populations of Primulina.

Conclusions

These two novel strategies may provide a pathway for similar research in other non-model species. The newly developed high polymorphic loci in this study will promote further the phylogenetic and population genetic studies in Primulina and other genera of the family Gesneriaceae.

Electronic supplementary material

The online version of this article (10.1186/s12862-017-1067-z) contains supplementary material, which is available to authorized users.

Both temperature fluctuations and East Asian monsoons have driven plant diversification in the karst ecosystems from southern China

Karst ecosystems in southern China are species-rich and have high levels of endemism, yet little is known regarding the evolutionary processes responsible for the origin and diversification of karst biodiversity. The genus Primulina (Gesneriaceae) comprises ca. 170 species endemic to southern China with high levels of ecological (edaphic) specialization, providing an exceptional model to study the plant diversification in karsts. We used molecular data from nine chloroplast and 11 nuclear regions and macroevolutionary analyses to assess the origin and cause of species diversification due to palaeoenvironmental changes and edaphic specialization in Primulina. We found that speciation was positively associated with changes in past temperatures and East Asian monsoons through the evolutionary history of Primulina. Climatic change around the mid-Miocene triggered an early burst followed by a slowdown of diversification rate towards the present with the climate cooling. We detected different speciation rates among edaphic types, and transitions among soil types were infrequently and did not impact the overall speciation rate. Our findings suggest that both global temperature changes and East Asian monsoons have played crucial roles in floristic diversification within the karst ecosystems in southern China, such that speciation was higher when climate was warmer and wetter. This is the first study to directly demonstrate that past monsoon activity is positively correlated with speciation rate in East Asia. This case study could motivate further investigations to assess the impacts of past environmental changes on the origin and diversification of biodiversity in global karst ecosystems, most of which are under threat.

Radionuclide biokinetics in the Russian sturgeon and phylogenetic consistencies with cartilaginous and bony marine fishes

The biokinetics of eight radionuclides (²⁴¹Am, ¹⁰⁹Cd, ¹³⁴Cs, ⁷⁵Se, ⁵⁴Mn, ^110mAg, ⁶⁵Zn, ⁶⁰Co) absorbed from the aquatic medium by juvenile Russian sturgeon (Acipenser gueldenstaedtii) were experimentally determined in fresh (0.42‰) and brackish (9.0‰) waters, of a similar salinity range to the Caspian Sea, and in conjunction with chemical speciation modelling. Uptake and loss rate constants were determined for each radionuclide for a 14 day exposure at each salinity and during 28 days of exposure to radionuclide-free conditions. Whole body (wet): water concentration factors (CF) achieved over 14 days for these eight radionuclides were used in a comparison with the same radionuclide CFs previously determined experimentally for six species of marine teleosts and chondrichthyans, to further test a phylogeny-based model of multi-nuclide bioaccumulation based on marine chordates. Multivariate analyses (multidimensional scaling and hierarchical clustering) identified the relative affinities among these taxa and also those radionuclides which distinguished most between them, in their differing CFs. They consistently showed that sturgeon aggregated as a group, which was also slightly differentiated with salinity. Sturgeon were distinguished from all teleosts and chondrichthyans but were more dissimilar from chondrichthyans than teleosts, in accordance with sturgeon's different periods of divergence from them in evolutionary time. Variable salinity among experiments may also cause changes in radionuclide bioaccumulation due to variations in (i) bioavailability (ii) osmolarity, and (iii) competitive inhibition of a radionuclide's bioaccumulation by its stable analogue or metabolic model. Their potentially confounding effects on these patterns of radionuclide CFs among taxa were critically evaluated for those radionuclides which discriminated most between sturgeon and teleosts or chondrichthyans. Bioavailability, osmolarity and competitive inhibition effects were identified among salinity treatments, however they were not appreciable enough to override the phylogeny-based signal. The results of this study are thus consistent with a phylogeny-based model of radionuclide bioaccumulation by marine chordates being valid for a fish species living in lower salinity regimes.

Speciation history of a species complex of Primulina eburnea (Gesneriaceae) from limestone karsts of southern China, a biodiversity hot spot

Limestone karsts in southern China are characterized by high edaphic and topographic heterogeneity and host high levels of species richness and endemism. However, the evolutionary mechanisms for generating such biodiversity remain poorly understood. Here, we performed species delimitation, population genetic analyses, simulations of gene flow and analyses of floral morphological traits to infer the geographic history of speciation in a species complex of Primulina eburnea from limestone karsts of southern China. Using Bayesian species delimitation, we determined that there are seven distinct species that correspond well to the putative morphological species. Species tree reconstruction, Structure and Neighbor‐Net analyses all recovered four lineages in agreement with currently species geographic boundaries. High levels of genetic differentiation were observed both within and among species. Isolation–migration coalescent analysis provides evidence for significant but low gene flow among species. Approximate Bayesian computation (ABC) analysis supports a scenario of historical gene flow rather than recent contemporary gene flow for most species divergences. Finally, we found no evidence of divergent selection contributing to population differentiation of a suite of flower traits. These results support the prevalence of allopatric speciation and highlight the role of geographic isolation in the diversification process. At small geographic scales, limited hybridization occurred in the past between proximate populations but did not eliminate species boundaries. We conclude that limited gene flow might have been the predominant evolutionary force in promoting population differentiation and speciation.

Population Genetic Structure and Phylogeography of Camellia flavida (Theaceae) Based on Chloroplast and Nuclear DNA Sequences

Camellia flavida is an endangered species of yellow camellia growing in limestone mountains in southwest China. The current classification of C. flavida into two varieties, var. flavida and var. patens, is controversial. We conducted a genetic analysis of C. flavida to determine its taxonomic structure. A total of 188 individual plants from 20 populations across the entire distribution range in southwest China were analyzed using two DNA fragments: a chloroplast DNA fragment from the small single copy region and a single-copy nuclear gene called phenylalanine ammonia-lyase (PAL). Sequences from both chloroplast and nuclear DNA were highly diverse; with high levels of genetic differentiation and restricted gene flow. This result can be attributed to the high habitat heterogeneity in limestone karst, which isolates C. flavida populations from each other. Our nuclear DNA results demonstrate that there are three differentiated groups within C. flavida: var. flavida 1, var. flavida 2, and var. patens. These genetic groupings are consistent with the morphological characteristics of the plants. We suggest that the samples included in this study constitute three taxa and the var. flavida 2 group is the genuine C. flavida. The three groups should be recognized as three management units for conservation concerns.

Adaptive molecular evolution of the two-pore channel 1 gene TPC1 in the karst-adapted genus Primulina (Gesneriaceae)

BACKGROUND AND AIMS

Limestone karst areas possess high floral diversity and endemism. The genus Primulina, which contributes to the unique calcicole flora, has high species richness and exhibit specific soil-based habitat associations that are mainly distributed on calcareous karst soils. The adaptive molecular evolutionary mechanism of the genus to karst calcium-rich environments is still not well understood. The Ca²⁺-permeable channel TPC1 was used in this study to test whether its gene is involved in the local adaptation of Primulina to karst high-calcium soil environments.

METHODS

Specific amplification and sequencing primers were designed and used to amplify the full-length coding sequences of TPC1 from cDNA of 76 Primulina species. The sequence alignment without recombination and the corresponding reconstructed phylogeny tree were used in molecular evolutionary analyses at the nucleic acid level and amino acid level, respectively. Finally, the identified sites under positive selection were labelled on the predicted secondary structure of TPC1.

KEY RESULTS

Seventy-six full-length coding sequences of Primulina TPC1 were obtained. The length of the sequences varied between 2220 and 2286 bp and the insertion/deletion was located at the 5' end of the sequences. No signal of substitution saturation was detected in the sequences, while significant recombination breakpoints were detected. The molecular evolutionary analyses showed that TPC1 was dominated by purifying selection and the selective pressures were not significantly different among species lineages. However, significant signals of positive selection were detected at both TPC1 codon level and amino acid level, and five sites under positive selective pressure were identified by at least three different methods.

CONCLUSIONS

The Ca²⁺-permeable channel TPC1 may be involved in the local adaptation of Primulina to karst Ca²⁺-rich environments. Different species lineages suffered similar selective pressure associated with calcium in karst environments, and episodic diversifying selection at a few sites may play a major role in the molecular evolution of Primulina TPC1.

Significant Phylogenetic Signal and Climate-Related Trends in Leaf Caloric Value from Tropical to Cold-Temperate Forests

Leaf caloric value (LCV) is a useful index to represent the conversion efficiency of leaves for solar energy. We investigated the spatial pattern of LCV and explored the factors (phylogeny, climate, and soil) that influence them at a large scale by determining LCV standardized by leaf area in 920 plant species from nine forest communities along the 3700 km North-South Transect of Eastern China. LCV ranged from 0.024 to 1.056 kJ cm-2 with an average of 0.151 kJ cm-2. LCV declined linearly with increasing latitude along the transect. Altogether, 57.29% of the total variation in LCV was explained by phylogenetic group (44.03% of variation), climate (1.27%), soil (0.02%) and their interacting effects. Significant phylogenetic signals in LCV were observed not only within forest communities but also across the whole transect. This phylogenetic signal was higher at higher latitudes, reflecting latitudinal change in the species composition of forest communities from complex to simple. We inferred that climate influences the spatial pattern of LCV through directly regulating the species composition of plant communities, since most plant species might tolerate only a limited temperature range. Our findings provide new insights into the adaptive mechanisms in plant traits in future studies.

Invariant allometric scaling of nitrogen and phosphorus in leaves, stems, and fine roots of woody plants along an altitudinal gradient

Nitrogen (N) to phosphorus (P) allocation in plant organs is of particular interest, as both elements are important to regulate plant growth. We analyzed the scaling relationship of N and P in leaves, stems and fine roots of 224 plant species along an altitudinal transect (500-2,300 m) on the northern slope of Changbai Mountain, China. We tested whether the scaling relationships of N and P were conserved in response to environmental variations. We found that the N and P concentrations of the leaves, stems and fine roots decreased, whereas the N:P ratios increased with increasing altitude. Allometric scaling relationships of N and P were found in the leaves, stems and fine roots, with allometric exponents of 0.78, 0.71 and 0.87, respectively. An invariant allometric scaling of N and P in the leaves, stems and fine roots was detected for woody plants along the altitudinal gradient. These results may advance our understanding of plant responses to climate change, and provide a basis for practical implication of various ecological models.

The first genetic linkage map of Primulina eburnea (Gesneriaceae) based on EST-derived SNP markers

Primulina eburnea is a promising candidate for domestication and floriculture, since it is easy to culture and has beautiful flowers. An F₂ population of 189 individuals was established for the construction of first-generation linkage maps based on expressed sequence tags-derived single-nucleotide polymorphism markers using the massARRAY genotyping platform. Of the 232 screened markers, 215 were assigned to 18 LG according to the haploid number of chromosomes in the species. The linkage map spanned a total of 3774.7 cM with an average distance of 17.6 cM between adjacent markers. This linkage map provides a framework for identification of important genes in breeding programmes.