Testing the ectomycorrhizal-dominance hypothesis for ecosystem multifunctionality in a subtropical mountain forest.

Mycorrhizal associations are key mutualisms that shape the structure of forest communities and multiple ecosystem functions. However, we lack a framework for predicting the varying dominance of distinct mycorrhizal associations in an integrated proxy of multifunctionality across ecosystems. Here, we used the datasets containing diversity of mycorrhizal associations and 18 ecosystem processes related to supporting, provisioning, and regulating services to examine how the dominance of ectomycorrhiza (EcM) associations affects ecosystem multifunctionality in subtropical mountain forests in Southwest China. Meanwhile, we synthesized the prevalence of EcM-dominant effects on ecosystem functioning in forest biomes. Our results demonstrated that elevation significantly modified the distributions of EcM trees and fungal dominance, which in turn influenced multiple functions simultaneously. Multifunctionality increased with increasing proportion of EcM associations, supporting the ectomycorrhizal-dominance hypothesis. Meanwhile, we observed that the impacts of EcM dominance on individual ecosystem functions exhibited different relationships among forest biomes. Our findings highlight the importance of ectomycorrhizal dominance in regulating multifunctionality in subtropical forests. However, this ectomycorrhizal feedback in shaping ecosystem functions cannot necessarily be generalized across forests. Therefore, we argue that the predictions for ecosystem multifunctionality in response to the shifts of mycorrhizal composition could vary across space and time.

Differential expressions of anthocyanin synthesis genes underlie flower color divergence in a sympatric Rhododendron sanguineum complex.

BACKGROUND: The Rhododendron sanguineum complex is endemic to alpine mountains of northwest Yunnan and southeast Tibet of China. Varieties in this complex exhibit distinct flower colors even at the bud stage. However, the underlying molecular regulations for the flower color variation have not been well characterized. Here, we investigated this via measuring flower reflectance profiles and comparative transcriptome analyses on three coexisting varieties of the R. sanguineum complex, with yellow flush pink, bright crimson, and deep blackish crimson flowers respectively. We compared the expression levels of differentially-expressed-genes (DEGs) of the anthocyanin / flavonoid biosynthesis pathway using RNA-seq and qRT-PCR data. We performed clustering analysis based on transcriptome-derived Single Nucleotide Polymorphisms (SNPs) data, and finally analyzed the promoter architecture of DEGs. RESULTS: Reflectance spectra of the three color morphs varied distinctively in the range between 400 and 700 nm, with distinct differences in saturation, brightness, hue, and saturation/hue ratio, an indirect measurement of anthocyanin content. We identified 15,164 orthogroups that were shared among the three varieties. The SNP clustering analysis indicated that the varieties were not monophyletic. A total of 40 paralogous genes encoding 12 enzymes contributed to the flower color polymorphism. These anthocyanin biosynthesis-related genes were associated with synthesis, modification and transportation properties (RsCHS, RsCHI, RsF3H, RsF3'H, RsFLS, RsANS, RsAT, RsOMT, RsGST), as well as genes involved in catabolism and degradation (RsBGLU, RsPER, RsCAD). Variations in sequence and cis-acting elements of these genes might correlate with the anthocyanin accumulation, thus may contribute to the divergence of flower color in the R. sanguineum complex. CONCLUSIONS: Our results suggested that the varieties are very closely related and flower color variations in the R. sanguineum complex correlate tightly with the differential expression levels of genes involved in the anabolic and catabolic synthesis network of anthocyanin. Our study provides a scenario involving intricate relationships between genetic mechanisms for floral coloration accompanied by gene flow among the varieties that may represent an early case of pollinator-mediated incipient sympatric speciation.

Repeated intercontinental migrations and recurring hybridizations characterise the evolutionary history of yew (Taxus L.).

The genus Taxus (Taxaceae) consists of 16 genetically well-defined lineages that are predominantly distributed across the Northern hemisphere. We investigated its biogeographic origin and evolutionary history by sampling 13 chloroplast gene sequences, the nuclear internal transcribed spacers (ITS) and NEEDLY sequences for all 16 lineages. We applied Maximum Parsimony and Bayesian Inference analyses to infer their phylogenetic relationships, time-calibrated phylogenies using BEAST and inferred the ancestral area of occupancy with BioGeoBEARS. We found strong evidence for the hybrid origin of three lineages and dated these events to a rather narrow time window of 6.8-4.9 million years ago (Mya). The dated phylogenies inferred an Upper Cretaceous origin of the genus, with the extant lineages diversifying in North America much later during the Oligocene/early Miocene. Repeated migrations via the Bering land bridge to Eurasia and back were further inferred, with the return to North America as a possible result of vicariance. The diversification in Eurasia (from ~8 Mya onwards) coincided with the orogeny of the Hengduan Mountains, the intensification of the East Asian summer monsoon and the occupancy of ecological niches by lineages that experienced secondary contacts and hybridizations in the Hengduan Mountains and Qinling Mountain, especially around the Sichuan basin. We provide a hypothesis for the evolution of extant lineages of Taxus, a genus with an old and complex evolutionary history. The study highlights that the history of complex species can be unravelled with a careful dissection of phylogenetic signals.

Development of polymorphic microsatellite markers for tree peony Paeonia delavayi (Paeoniaceae) using ddRAD-seq data.

Microsatellite markers were developed for the tree peony Paeonia delavayi to investigate fine scale population genetics of this species. Using ddRAD-seq data from twenty individuals of P. delavayi, we identified 529 polymorphic microsatellite loci, of which 195 were suitable for designing microsatellite primers. Of the 120 microsatellite loci selected for validation, 20 were successfully amplified with clear peaks and displayed polymorphism. Three populations were genotyped using the 20 polymorphic microsatellites. The number of alleles per locus ranged from two to thirteen. Observed and expected heterozygosity ranged from 0 to 0.941 and 0 to 0.834 respectively. The cross-species amplification test using five individuals from a population of P. ludlowii showed that 15 of the 20 polymorphic loci were successfully amplified, and four loci showed polymorphism. Among the 22 alleles occurring in P. ludlowii across fifteen loci, eight alleles across five loci were exclusive to P. ludlowii. The results demonstrate that ddRAD-seq is an efficient method for the development of microsatellite markers for non-model organisms with large genomes. The newly developed markers will be valuable tools to investigate the genetic diversity, genetic structure, and gene flow of P. delavayi from local to regional spatial scales.

Niche convergence and biogeographic history shape elevational tree community assembly in a subtropical mountain forest.

Niche convergence or conservatism have been proposed as essential mechanisms underlying elevational plant community assembly in tropical mountain ecosystems. Subtropical mountains, compared to tropical mountains, are likely to be shaped by a mixing of different geographic affinities of species and remain somehow unclear. Here, we used 31 0.1-ha permanent plots distributed in subtropical forests on the eastern and western aspects of the Gaoligong Mountains, southwest China between 1498 m and 3204 m a.sl. to evaluate how niche-based and biogeographic processes shape tree community assembly along elevational gradients. We analyzed the elevational patterns of taxonomic, phylogenetic and functional diversity, as well as of individual traits, and assessed the relative importance of environmental effects on these diversity measures. We then classified tree species as being either tropical affiliated or temperate affiliated and estimated their contribution to the composition of biogeographic affinities. Species richness decreased with elevation, and species composition showed apparent turnover across the aspects and elevations. Most traits exhibited convergent patterns across the entire elevational gradient. Phylogenetic and functional diversity showed opposing patterns, with phylogenetic diversity increasing and functional diversity decreasing with elevation. Soil nutrients, especially phosphorus and nitrogen, appeared to be the main abiotic variables driving the elevational diversity patterns. Communities at lower elevations were occupied by tropical genera, while highlands contained species of tropical and temperate biogeographic affinities. Moreover, the high phylogenetic diversity at high elevations were likely due to differences in evolutionary history between temperate and tropical species. Our results highlight the importance of niche convergence of tropical species and the legacy of biogeographic history on the composition and structure of subtropical mountain forests. Furthermore, limited soil phosphorus caused traits divergence and the partitioning for different forms of phosphorus may explain the high biodiversity found in phosphorus-limited subtropical forests.

Variation in gene expression along an elevation gradient of Rhododendron sanguineum var. haemaleum assessed in a comparative transcriptomic analysis.

Selection along environmental gradients may play a vital role in driving adaptive evolution. Nevertheless, genomic variation and genetic adaptation along environmental clines remains largely unknown in plants in alpine ecosystems. To close this knowledge gap, we assayed transcriptomic profiles of late flower bud and early leaf bud of Rhododendron sanguineum var. haemaleum from four different elevational belts between 3,000 m and 3,800 m in the Gaoligong Mountains. By comparing differences in gene expression of these samples, a gene co-expression network (WGCNA) was constructed to identify candidate genes related to elevation. We found that the overall gene expression patterns are organ-specific for the flower and leaf. Differentially expressed unigenes were identified in these organs. In flowers, these were mainly related to terpenoid metabolism (RsHMGR, RsTPS), while in leaves mainly related to anthocyanin biosynthesis (RsCHS, RsF3'5'H). Terpenoids are the main components of flower scent (fragrance) likely attracting insects for pollination. In response to fewer pollinators at higher elevation zone, it seems relatively less scent is produced in flower organs to reduce energy consumption. Secondary metabolites in leaves such as anthocyanins determine the plants' alternative adaptive strategy to extreme environments, such as selective pressures of insect herbivory from environmental changes and substrate competition in biosynthesis pathways at high elevations. Our findings indicated that the gene expression profiles generated from flower and leaf organs showed parallel expression shifts but with different functionality, suggesting the existence of flexibility in response strategies of plants exposed to heterogeneous environments across elevational gradients. The genes identified here are likely to be involved in the adaptation of the plants to these varying mountainous environments. This study thus contributes to our understanding of the molecular mechanisms of adaptation in response to environmental change.

A Multiplex PCR System of Novel Microsatellite Loci for Population Genetic Application in Walnuts.

Multiplex polymerase chain reaction (PCR) of microsatellite loci allows for simultaneous amplification of two or more pairs of primers in a single PCR reaction; hence, it is cost and time effective. However, very few attempts have been reported in non-model species. In this study, by combining a genome-based de novo development and cross-species application approach, a multiplex PCR system comprising 5 PCR reactions of 33 microsatellites consisting of 26 novel genomic and 7 literature-sourced loci was tested for polymorphisms, cross-species transferability, and the ability to assess genetic diversity and population structure of three walnut species (Juglans spp.). We found that the genome-based approach is more efficient than other methods. An allelic ladder was developed for each locus to enhance consistent genotyping among laboratories. The population genetic analysis results showed that all 33 loci were successfully transferred across the three species, showing high polymorphism and a strong genetic structure. Hence, the multiplex PCR system is highly applicable in walnut species. Furthermore, we propose an efficient pipeline to characterize and genotype polymorphic microsatellite loci. The novel toolbox developed here will aid future ecology and evolution studies in walnut and could serve as a model for other plant species.

Multitrophic diversity and biotic associations influence subalpine forest ecosystem multifunctionality.

Biodiversity across multiple trophic levels is required to maintain multiple ecosystem functions. Yet it remains unclear how multitrophic diversity and species interactions regulate ecosystem multifunctionality. Here, combining data from 9 different trophic groups (including trees, shrubs, herbs, leaf mites, small mammals, bacteria, pathogenic fungi, saprophytic fungi, and symbiotic fungi) and 13 ecosystem functions related to supporting, provisioning, and regulating services, we used a multitrophic perspective to evaluate the effects of elevation, diversity, and network complexity on scale-dependent subalpine forest multifunctionality. Our results demonstrated that elevation and soil pH significantly modified species composition and richness across multitrophic groups and influenced multiple functions simultaneously. We present evidence that species richness across multiple trophic groups had stronger effects on multifunctionality than species richness at any single trophic level. Moreover, biotic associations, indicating the complexity of trophic networks, were positively associated with multifunctionality. The relative effects of diversity on multifunctionality increased at the scale of the larger community compared to a scale accounting for neighboring interactions. Our results highlight the paramount importance of scale- and context-dependent multitrophic diversity and interactions for a better understanding of mountain ecosystem multifunctionality in a changing world.

Testing genome skimming for species discrimination in the large and taxonomically difficult genus Rhododendron.

Standard plant DNA barcodes based on 2-3 plastid regions, and nrDNA ITS show variable levels of resolution, and fail to discriminate among species in many plant groups. Genome skimming to recover complete plastid genome sequences and nrDNA arrays has been proposed as a solution to address these resolution limitations. However, few studies have empirically tested what gains are achieved in practice. Of particular interest is whether adding substantially more plastid and nrDNA characters will lead to an increase in discriminatory power, or whether the resolution limitations of standard plant barcodes are fundamentally due to plastid genomes and nrDNA not tracking species boundaries. To address this, we used genome skimming to recover near-complete plastid genomes and nuclear ribosomal DNA from Rhododendron species and compared discrimination success with standard plant barcodes. We sampled 218 individuals representing 145 species of this species-rich and taxonomically difficult genus, focusing on the global biodiversity hotspots of the Himalaya-Hengduan Mountains. Only 33% of species were distinguished using ITS+matK+rbcL+trnH-psbA. In contrast, 55% of species were distinguished using plastid genome and nrDNA sequences. The vast majority of this increase is due to the additional plastid characters. Thus, despite previous studies showing an asymptote in discrimination success beyond 3-4 plastid regions, these results show that a demonstrable increase in discriminatory power is possible with extensive plastid genome data. However, despite these gains, many species remain unresolved, and these results also reinforce the need to access multiple unlinked nuclear loci to obtain transformative gains in species discrimination in plants.

Genetic analysis of walnut cultivars from southwest China: Implications for germplasm improvement.

Walnuts are highly valued for their rich nutritional profile and wide medicinal applications. This demand has led to the intensification of breeding activities in major walnut production areas such as southwest China, in order to develop more superior cultivars. With the increasing number of cultivars, accurate identification becomes fundamental to selecting the right cultivar for grafting, industrial processing or development of new cultivars. To ensure proper identification of cultivars and understand the genetic structure of wild and cultivated material, we genotyped 362 cultivated and wild individuals of walnut trees from southwest China (with two additional populations from Xinjiang, plus three cultivars from Canada, France and Belgium) using 36 polymorphic microsatellite loci. We found relatively low indices of genetic diversity (H O = 0.570, H E = 0.404, N A = 2.345) as well as a high level of clonality (>85% of cultivars), indicating reliance on genetically narrow sources of parental material for breeding. Our STRUCTURE and PCoA analyses generally delineated the two species, though considerable levels of introgression were also evident. More significantly, we detected a distinct genetic group of cultivated Juglans sigillata, which mainly comprised individuals of the popular 'Yangbidapao' landrace. Finally, a core set of 18 SSR loci was selected, which was capable of identifying 32 cultivars. In a nutshell, our results call for more utilization of genetically disparate material, including wild walnut trees, as parental sources to breed for more cultivars. The data reported herein will significantly contribute towards the genetic improvement and conservation of the walnut germplasm in southwest China.

Development of 32 novel microsatellite loci in Juglans sigillata using genomic data.

PREMISE: A novel set of microsatellite markers was developed for Juglans sigillata (Juglandaceae), an endemic walnut species in southwestern China, to facilitate cultivar identification and future investigations into the genetic structure and domestication history of this species and its close relatives. METHODS AND RESULTS: We developed 32 microsatellite loci for J. sigillata using genomic data and used them to examine 60 individuals from three natural populations. A high level of polymorphism was detected by these primers, with up to eight alleles observed per locus, and an average of four alleles across populations. The levels of observed and expected heterozygosity ranged from 0.000-1.000 and 0.000-0.785, respectively. All but two of the loci were also successfully amplified in three closely related Eurasian Juglans species (J. regia, J. cathayensis, and J. mandshurica). CONCLUSIONS: The microsatellite loci identified here provide a powerful resource for examining the genetic structure and domestication history of Juglans, as well as identification of its cultivars.

Prevalence of isomeric plastomes and effectiveness of plastome super-barcodes in yews (Taxus) worldwide.

Taxus (yew) is both the most species-rich and taxonomically difficult genus in Taxaceae. To date, no study has elucidated the complexities of the plastid genome (plastome) or examined the possibility of whole plastomes as super-barcodes across yew species worldwide. In this study, we sequenced plastomes from two to three individuals for each of the 16 recognized yew species (including three potential cryptics) and Pseudotaxus chienii. Our comparative analyses uncovered several gene loss events that independently occurred in yews, resulting in a lower plastid gene number than other Taxaceous genera. In Pseudotaxus and Taxus, we found two isomeric arrangements that differ by the orientation of a 35 kb fragment flanked by "trnQ-IRs". These two arrangements exist in different ratios within each sampled individual, and intraspecific shifts in major isomeric arrangements are first reported here in Taxus. Moreover, we demonstrate that entire plastomes can be used to successfully discriminate all Taxus species with 100% support, suggesting that they are useful as super-barcodes for species identification. We also propose that accD and rrn16-rrn23 are promising special barcodes to discriminate yew species. Our newly developed Taxus plastomic sequences provide a resource for super-barcodes and conservation genetics of several endangered yews and serve as comprehensive data to improve models of plastome complexity in Taxaceae as a whole and authenticate Taxus species.

Integrating a comprehensive DNA barcode reference library with a global map of yews (Taxus L.) for forensic identification.

Rapid and accurate identification of endangered species is a critical component of biosurveillance and conservation management, and potentially policing illegal trades. However, this is often not possible using traditional taxonomy, especially where only small or preprocessed parts of plants are available. Reliable identification can be achieved via a comprehensive DNA barcode reference library, accompanied by precise distribution data. However, these require extensive sampling at spatial and taxonomic scales, which has rarely been achieved for cosmopolitan taxa. Here, we construct a comprehensive DNA barcode reference library and generate distribution maps using species distribution modelling (SDM), for all 15 Taxus species worldwide. We find that trnL-trnF is the ideal barcode for Taxus: It can distinguish all Taxus species and in combination with ITS identify hybrids. Among five analysis methods tested, NJ was the most effective. Among 4,151 individuals screened for trnL-trnF, 73 haplotypes were detected, all species-specific and some population private. Taxonomical, geographical and genetic dimensions of sampling strategy were all found to affect the comprehensiveness of the resulting DNA barcode library. Maps from SDM showed that most species had allopatric distributions, except T. mairei in the Sino-Himalayan region. Using the barcode library and distribution map data, two unknown forensic samples were identified to species (and in one case, population) level and another was determined as a putative interspecific hybrid. This integrated species identification system for Taxus can be used for biosurveillance, conservation management and to monitor and prosecute illegal trade. Similar identification systems are recommended for other IUCN- and CITES-listed taxa.

Development of microsatellite markers for Sargentodoxa cuneata (Lardizabalaceae) using next-generation sequencing technology.

PREMISE OF THE STUDY: Microsatellite loci were developed for a woody deciduous liana, Sargentodoxa cuneata (Lardizabalaceae), to help infer the evolutionary histories of ancient monotypic genera in subtropical China. METHODS AND RESULTS: Using next-generation sequencing (Illumina MiSeq) technology, 21 polymorphic primer sets were identified in three wild populations. The number of alleles per locus ranged from one to seven. The expected and observed heterozygosities varied from 0 to 0.788 and 0 to 0.917, respectively. Transferability analyses were performed in Stauntonia chinensis, Akebia trifoliata, and A. quinata. Eighteen (85.7%), 18 (85.7%), and 17 (81.0%) markers were successfully amplified, respectively. CONCLUSIONS: The newly developed markers will facilitate further studies on genetic diversity and phylogeographic patterns throughout the distributional range of S. cuneata. This set of microsatellite primers represents the second report on molecular markers in Lardizabalaceae.

Development of microsatellite loci for Cyclocarya paliurus (Juglandaceae), a monotypic species in subtropical China.

PREMISE OF THE STUDY: Microsatellite loci were developed for a monotypic species endemic to subtropical China, Cyclocarya paliurus, to help infer the evolutionary histories of ancient monotypic genera in subtropical China. • METHODS AND RESULTS: Using the Fast Isolation by AFLP of Sequences COntaining repeats (FIASCO) protocol, 28 primer sets were identified in two wild populations. All loci were polymorphic, with the number of alleles per locus ranging from two to eight. The expected and observed heterozygosities ranged from 0.153 to 0.802 and from 0 to 0.750, respectively. The transferability of the 28 primer pairs was tested on Juglans regia, Pterocarya stenoptera, and Platycarya strobilacea. Twenty-one (75.0%), 22 (78.6%), and 15 (53.6%) markers were successfully amplified in J. regia, P. stenoptera, and P. strobilacea, respectively. • CONCLUSIONS: These loci will be useful for in-depth analysis of genetic diversity and phylogeographical variation throughout the distribution range of C. paliurus.