Xiao-Fei Ma

Disentangling the Roles of History and Local Selection in Shaping Clinal Variation of Allele Frequencies and Gene Expression in Norway Spruce (Picea abies)

Understanding the genetic basis of local adaptation is challenging due to the subtle balance among conflicting evolutionary forces that are involved in its establishment and maintenance. One system with which to tease apart these difficulties is clines in adaptive characters. Here we analyzed genetic and phenotypic variation in bud set, a highly heritable and adaptive trait, among 18 populations of Norway spruce (Picea abies), arrayed along a latitudinal gradient ranging from 47°N to 68°N. We confirmed that variation in bud set is strongly clinal, using a subset of five populations. Genotypes for 137 single-nucleotide polymorphisms (SNPs) chosen from 18 candidate genes putatively affecting bud set and 308 control SNPs chosen from 264 random genes were analyzed for patterns of genetic structure and correlation to environment. Population genetic structure was low (FST = 0.05), but latitudinal patterns were apparent among Scandinavian populations. Hence, part of the observed clinal variation should be attributable to population demography. Conditional on patterns of genetic structure, there was enrichment of SNPs within candidate genes for correlations with latitude. Twenty-nine SNPs were also outliers with respect to FST. The enrichment for clinal variation at SNPs within candidate genes (i.e., SNPs in PaGI, PaPhyP, PaPhyN, PaPRR7, and PaFTL2) indicated that local selection in the 18 populations, and/or selection in the ancestral populations from which they were recently derived, shaped the observed cline. Validation of these genes using expression studies also revealed that PaFTL2 expression is significantly associated with latitude, thereby confirming the central role played by this gene in the control of phenology in plants.

Genetic differentiation, clinal variation and phenotypic associations with growth cessation across the Populus tremula Photoperiodic Pathway

Perennial plants monitor seasonal changes through changes in environmental conditions such as the quantity and quality of light. To ensure a correct initiation of critical developmental processes, such as the initiation and cessation of growth, plants have adapted to a spatially variable light regime and genes in the photoperiodic pathway have been implicated as likely sources for these adaptations. Here we examine genetic variation in genes from the photoperiodic pathway in Populus tremula (Salicaceae) for signatures diversifying selection in response to varying light regimes across a latitudinal gradient. We fail to identify any loci with unusually high levels of genetic differentiation among populations despite identifying four SNPs that show significant allele frequency clines with latitude. We do, however, observe large covariance in allelic effects across populations for growth cessation, a highly adaptive trait in P. tremula. High covariance in allelic effects is a signature compatible with diversifying selection along an environmental gradient. We also observe significantly higher heterogeneity in genetic differentiation among SNPs from the photoperiod genes than among SNPs from randomly chosen genes. This suggests that spatially variable selection could be affecting genes from the photoperiod pathway even if selection is not strong enough to cause individual loci to be identified as outliers. SNPs from three genes in the photoperiod pathway (PHYB2, LHY1, and LHY2) show significant associations with natural variation in growth cessation. Collectively these SNPs explain 10–15% of the phenotypic variation in growth cessation. Covariances in allelic effects across populations help explain an additional 5–7% of the phenotypic variation in growth cessation.

Intragenomic and interspecific 5S rDNA sequence variation in five Asian pines

Patterns of intragenomic and interspecific variation of 5S rDNA in Pinus (Pinaceae) were studied by cloning and sequencing multiple 5S rDNA repeats from individual trees. Five pines, from both subgenera, Pinus and Strobus, were selected. The 5S rDNA repeat in pines has a conserved 120-base pair (bp) transcribed region and an intergenic spacer region of variable length (382-608 bp). The evolutionary rate in the spacer region is three- to sevenfold higher than in the genic region. We found substantial sequence divergence between the two subgenera. Intragenomic sequence heterogeneity was high for all species, and more than 86% of the clones within each individual were unique. The 5S gene tree revealed that different 5S repeats within individuals are polyphyletic, indicating that their ancestral divergence preceded the speciation events. The degrees of interspecific and intragenomic divergence among diploxylon pines are similar. The observed sequence patterns suggest that concerted evolution has been acting after the diversification of the two subgenera but very weak after the speciation of the four diploxylon pines. Sequence patterns in P. densata are consistent with hybrid origin. It had higher intragenomic diversity and maintained polymorphic copies of the parental types in addition to new and recombinant types unique to the hybrid.

Transcriptomic Analysis of a Tertiary Relict Plant, Extreme Xerophyte Reaumuria soongorica to Identify Genes Related to Drought Adaptation

Background Reaumuria soongorica is an extreme xerophyte shrub widely distributed in the desert regions including sand dune, Gobi and marginal loess of central Asia which plays a crucial role to sustain and restore fragile desert ecosystems. However, due to the lacking of the genomic sequences, studies on R. soongorica had mainly limited in physiological responses to drought stress. Here, a deep transcriptomic sequencing of R. soongorica will facilitate molecular functional studies and pave the path to understand drought adaptation for a desert plant. Methodology/Principal Findings A total of 53,193,660 clean paired-end reads was generated from the Illumina HiSeq™ 2000 platform. By assembly with Trinity, we got 173,700 contigs and 77,647 unigenes with mean length of 677 bp and N50 of 1109 bp. Over 55% (43,054) unigenes were successfully annotated based on sequence similarity against public databases as well as Rfam and Pfam database. Local BLAST and Kyoto Encyclopedia of Genes and Genomes (KEGG) maps were used to further exhausting seek for candidate genes related to drought adaptation and a set of 123 putative candidate genes were identified. Moreover, all the C4 photosynthesis genes existed and were active in R. soongorica, which has been regarded as a typical C3 plant. Conclusion/Significance The assembled unigenes in present work provide abundant genomic information for the functional assignments in an extreme xerophyte R. soongorica, and will help us exploit the genetic basis of how desert plants adapt to drought environment in the near future.

Adaptive evolution of the Populus tremula photoperiod pathway

Perennial plants monitor seasonal changes through changes in environmental conditions such as the quantity and quality of light and genes in the photoperiodic pathway are known to be involved in controlling these processes. Here, we examine 25 of genes from the photoperiod pathway in Populus tremula (Salicaceae) for signatures of adaptive evolution. Overall, levels of synonymous polymorphism in the 25 genes are lower than at control loci selected randomly from the genome. This appears primarily to be caused by lower levels of synonymous polymorphism in genes associated with the circadian clock. Natural selection appears to play an important role in shaping protein evolution at several of the genes in the photoperiod pathways, which is highlighted by the fact that approximately 40% of the genes from the photoperiod pathway have estimates of selection on nonsynonymous polymorphisms that are significantly different from zero. A surprising observation we make is that circadian clock‐associated genes appear to be over‐represented among the genes showing elevated rates of protein evolution; seven genes are evolving under positive selection and all but one of these genes are involved in the circadian clock of Populus.

Identification of Differentially Expressed Genes in Leaf of Reaumuria soongorica under PEG-Induced Drought Stress by Digital Gene Expression Profiling

Reaumuria soongorica (Pall.) Maxim., a resurrection semi-shrub, is a typical constructive and dominant species in desert ecosystems in northwestern China. However, the gene expression characteristics of R. soongorica under drought stress have not been elucidated. Digital gene expression analysis was performed using Illumina technique to investigate differentially expressed genes (DEGs) between control and PEG-treated samples of R. soongorica. A total of 212,338 and 211,052 distinct tags were detected in the control and PEG-treated libraries, respectively. A total of 1,325 genes were identified as DEGs, 379 (28.6%) of which were up-regulated and 946 (71.4%) were down-regulated in response to drought stress. Functional annotation analysis identified numerous drought-inducible genes with various functions in response to drought stress. A number of regulatory proteins, functional proteins, and proteins induced by other stress factors in R. soongorica were identified. Alteration in the regulatory proteins (transcription factors and protein kinase) may be involved in signal transduction. Functional proteins, including flavonoid biosynthetic proteins, late embryogenesis abundant (LEA) proteins, small heat shock proteins (sHSP), and aquaporin and proline transporter may play protective roles in response to drought stress. Flavonoids, LEA proteins and sHSP function as reactive oxygen species scavenger or molecular chaperone. Aquaporin and proline transporters regulate the distribution of water and proline throughout the whole plant. The tolerance ability of R. soongorica may be gained through effective signal transduction and enhanced protection of functional proteins to reestablish cellular homeostasis. DEGs obtained in this study may provide useful insights to help further understand the drought-tolerant mechanism of R. soongorica.

A Picea abies linkage map based on SNP markers identifies QTLs for four aspects of resistance to Heterobasidion parviporum infection.

A consensus linkage map of Picea abies, an economically important conifer, was constructed based on the segregation of 686 SNP markers in a F1 progeny population consisting of 247 individuals. The total length of 1889.2 cM covered 96.5% of the estimated genome length and comprised 12 large linkage groups, corresponding to the number of haploid P. abies chromosomes. The sizes of the groups (from 5.9 to 9.9% of the total map length) correlated well with previous estimates of chromosome sizes (from 5.8 to 10.8% of total genome size). Any locus in the genome has a 97% probability to be within 10 cM from a mapped marker, which makes the map suited for QTL mapping. Infecting the progeny trees with the root rot pathogen Heterobasidion parviporum allowed for mapping of four different resistance traits: lesion length at the inoculation site, fungal spread within the sapwood, exclusion of the pathogen from the host after initial infection, and ability to prevent the infection from establishing at all. These four traits were associated with two, four, four and three QTL regions respectively of which none overlapped between the traits. Each QTL explained between 4.6 and 10.1% of the respective traits phenotypic variation. Although the QTL regions contain many more genes than the ones represented by the SNP markers, at least four markers within the confidence intervals originated from genes with known function in conifer defence; a leucoanthocyanidine reductase, which has previously been shown to upregulate during H. parviporum infection, and three intermediates of the lignification process; a hydroxycinnamoyl CoA shikimate/quinate hydroxycinnamoyltransferase, a 4-coumarate CoA ligase, and a R2R3-MYB transcription factor.

Reference gene selection for quantitative real-time PCR normalization in Reaumuria soongorica.

Despite its superiority for evaluating gene expression, real-time quantitative polymerase chain reaction (qPCR) results can be significantly biased by the use of inappropriate reference genes under different experimental conditions. Reaumuria soongorica is a dominant species of desert ecosystems in arid central Asia. Given the increasing interest in ecological engineering and potential genetic resources for arid agronomy, it is important to analyze gene function. However, systematic evaluation of stable reference genes should be performed prior to such analyses. In this study, the stabilities of 10 candidate reference genes were analyzed under 4 kinds of abiotic stresses (drought, salt, dark, and heat) within 4 accessions (HG010, HG020, XGG030, and XGG040) from 2 different habitats using 3 algorithms (geNorm, NormFinder, and BestKeeper). After validation of the ribulose-1,5-bisphosphate carboxylase/oxygenase large unite (rbcL) expression pattern, our data suggested that histone H2A (H2A) and eukaryotic initiation factor 4A-2 (EIF4A2) were the most stable reference genes, cyclophilin (CYCL) was moderate, and elongation factor 1α (EF1α) was the worst choice. This first systematic analysis for stably expressed genes will facilitate future functional analyses and deep mining of genetic resources in R. soongorica and other species of the Reaumuria genus.

Transcriptomic analysis of a psammophyte food crop, sand rice (Agriophyllum squarrosum) and identification of candidate genes essential for sand dune adaptation

BackgroundSand rice (Agriophyllum squarrosum) is an annual desert plant adapted to mobile sand dunes in arid and semi-arid regions of Central Asia. The sand rice seeds have excellent nutrition value and have been historically consumed by local populations in the desert regions of northwest China. Sand rice is a potential food crop resilient to ongoing climate change; however, partly due to the scarcity of genetic information, this species has undergone only little agronomic modifications through classical breeding during recent years.ResultsWe generated a deep transcriptomic sequencing of sand rice, which uncovers 67,741 unigenes. Phylogenetic analysis based on 221 single-copy genes showed close relationship between sand rice and the recently domesticated crop sugar beet. Transcriptomic comparisons also showed a high level of global sequence conservation between these two species. Conservation of sand rice and sugar beet orthologs assigned to response to salt stress gene ontology term suggests that sand rice is also a potential salt tolerant plant. Furthermore, sand rice is far more tolerant to high temperature. A set of genes likely relevant for resistance to heat stress, was functionally annotated according to expression levels, sequence annotation, and comparisons corresponding transcriptome profiling results in Arabidopsis.ConclusionsThe present work provides abundant genomic information for functional dissection of the important traits in sand rice. Future screening the genetic variation among different ecotypes and constructing a draft genome sequence will further facilitate agronomic trait improvement and final domestication of sand rice.

Population dynamics of Agriophyllum squarrosum, a pioneer annual plant endemic to mobile sand dunes, in response to global climate change

Climate change plays an important role in the transition of ecosystems. Stratigraphic investigations have suggested that the Asian interior experienced frequent transitions between grassland and desert ecosystems as a consequence of global climate change. Using maternally and bi-parentally inherited markers, we investigated the population dynamics of Agriophyllum squarrosum (Chenopodiaceae), an annual pioneer plant endemic to mobile sand dunes. Phylogeographic analysis revealed that A. squarrosum could originate from Gurbantunggut desert since ~1.6 Ma, and subsequently underwent three waves of colonisation into other deserts and sandy lands corresponding to several glaciations. The rapid population expansion and distribution range shifts of A. squarrosum from monsoonal climate zones suggested that the development of the monsoonal climate significantly enhanced the population growth and gene flow of A. squarrosum. These data also suggested that desertification of the fragile grassland ecosystems in the Qinghai-Tibetan Plateau was more ancient than previously suggested and will be aggravated under global warming in the future. This study provides new molecular phylogeographic insights into how pioneer annual plant species in desert ecosystems respond to global climate change, and facilitates evaluation of the ecological potential and genetic resources of future crops for non-arable dry lands to mitigate climate change.