Transcriptome sequencing reveals population differentiation in gene expression linked to functional traits and environmental gradients in the South African shrub Protea repens

Functional traits shape plant-plant interactions and recruitment in a hotspot of woody plant diversity

Understanding and predicting recruitment in species-rich plant communities requires identifying functional determinants of both density-independent performance and interactions. In a common-garden field experiment with 25 species of the woody plant genus Protea, we varied the initial spatial and taxonomic arrangement of seedlings and followed their survival and growth during recruitment. Neighbourhood models quantified how six key functional traits affect density-independent performance, interaction effects and responses. Trait-based neighbourhood models accurately predicted individual survival and growth from the initial spatial and functional composition of species-rich experimental communities. Functional variation among species caused substantial variation in density-independent survival and growth that was not correlated with interaction effects and responses. Interactions were spatially restricted but had important, predominantly competitive, effects on recruitment. Traits increasing the acquisition of limiting resources (water for survival and soil P for growth) mediated trade-offs between interaction effects and responses. Moreover, resprouting species had higher survival but reduced growth, likely reinforcing the survival-growth trade-off in adult plants. Resource acquisition of juvenile plants shapes Protea community dynamics with acquisitive species with strong competitive effects suffering more from competition. Together with functional determinants of density-independent performance, this makes recruitment remarkably predictable, which is critical for efficient restoration and near-term ecological forecasts of species-rich communities.

Morphological and transcriptomic responses/acclimations of erect-type submerged macrophyte Hydrilla verticillata both at low-light exposure and light recovery phases

Light intensity is a determinant for submerged macrophytes. Little is known about their molecular responses to low-light exposure, despite more informative and responsive than morphological traits. For erect-type submerged macrophytes, the stem is more crucial relative to the leaf in acclimation to low-light stress, but receives less attention. We determined morphological and stem transcriptomic responses/acclimations of Hydrilla verticillata to extremely and mildly low light (7.2 and 36 μmol photons m-2 s-1, respectively), that is, EL and ML, with the radiation intensity of 180 μmol photons m-2 s-1 as the control. Low-light exposure continued for 9 days, followed by a 7-day recovery phase (180 μmol photons m-2 s-1). At the exposure phase, the low-light treatments, in particular the EL, decreased the relative growth ratio, but induced greater height and longer stem internode distance and epidermal cell. Such responses/acclimations continued into the recovery phase, despite more or less changes in the magnitude. Transcriptome showed that the photosynthetic system was inhibited at the exposure phase, but the macrophyte adjusted hormone synthesis relating to cell division and elongation. Moreover, the EL activated cell stress responses such as DNA repair. Following light recovery, the macrophyte exhibited a strong-light response, although energy metabolism enhanced. Especially, the EL enriched the pathways relating to anthocyanin synthesis at such phase, indicating an activation of photoprotective mechanism. Our findings suggest that negative influences of low light occur at both low-light exposure and recovery phases, but submerged macrophytes would acclimate to light environments. Transcriptome can show molecular basis of plant responses/acclimations, including but not limited to morphology. This study establishes a bridge connecting morphological and molecular responses/acclimations.

The genome of the king protea, Protea cynaroides

The king protea (Protea cynaroides), an early-diverging eudicot, is the most iconic species from the Megadiverse Cape Floristic Region, and the national flower of South Africa. Perhaps best known for its iconic flower head, Protea is a key genus for the South African horticulture industry and cut-flower market. Ecologically, the genus and the family Proteaceae are important models for radiation and adaptation, particularly to soils with limited phosphorus bio-availability. Here, we present a high-quality chromosome-scale assembly of the P. cynaroides genome as the first representative of the fynbos biome. We reveal an ancestral whole-genome duplication event that occurred in the Proteaceae around the late Cretaceous that preceded the divergence of all crown groups within the family and its extant diversity in all Southern continents. The relatively stable genome structure of P. cynaroides is invaluable for comparative studies and for unveiling paleopolyploidy in other groups, such as the distantly related sister group Ranunculales. Comparative genomics in sequenced genomes of the Proteales shows loss of key arbuscular mycorrhizal symbiosis genes likely ancestral to the family, and possibly the order. The P. cynaroides genome empowers new research in plant diversification, horticulture and adaptation, particularly to nutrient-poor soils.

Transcriptomics of developing wild sunflower seeds from the extreme ends of a latitudinal gradient differing in seed oil composition

Seed oil composition, an important agronomic trait in cultivated sunflower, varies latitudinally across the native range of its wild progenitor. This pattern is thought to be driven by selection for a higher proportion of saturated fatty acids in southern populations compared with northern populations, likely due to the different temperatures experienced during seed germination. To investigate whether these differences in fatty acid composition between northern and southern populations correspond to transcriptional variation in the expression of genes involved in fatty acid metabolism, we sequenced RNA from developing seeds of sunflowers from Texas, USA, and Saskatchewan, Canada (the extreme ends of sunflower's latitudinal range) grown in a common garden. We found 4,741 genes to be differentially expressed between Texas and Canada, including several genes involved in lipid metabolism. Several differentially expressed lipid metabolism genes also colocalized with known oil quantitative trait loci (QTL). The genes producing stearoyl-ACP-desaturases (SAD) were of particular interest because of their known role in the conversion of fully saturated into unsaturated fatty acids. Two SAD genes were more highly expressed in seeds from Canadian populations, consistent with the observation of increased levels of unsaturated fatty acids in seeds from that region. We also constructed a gene co-expression network to investigate regional variation in network modules. The results of this analysis revealed regional differentiation for eight of 12 modules but no clear relationship with oil biosynthesis. Overall, the differential expression of SAD genes offers a partial explanation for the observed differences in seed oil composition between Texas and Canada, while the expression patterns of other metabolic genes suggest complex regulation of fatty acid production and usage across latitudes.

Adaptive divergence generates distinct plastic responses in two closely related Senecio species

The evolution of plastic responses to external cues allows species to maintain fitness in response to the environmental variations they regularly experience. However, it remains unclear how plasticity evolves during adaptation. To test whether distinct patterns of plasticity are associated with adaptive divergence, we quantified plasticity for two closely related but ecologically divergent Sicilian daisy species (Senecio, Asteraceae). We sampled 40 representative genotypes of each species from their native range on Mt. Etna and then reciprocally transplanted multiple clones of each genotype into four field sites along an elevational gradient that included the native elevational range of each species, and two intermediate elevations. At each elevation, we quantified survival and measured leaf traits that included investment (specific leaf area), morphology, chlorophyll fluorescence, pigment content, and gene expression. Traits and differentially expressed genes that changed with elevation in one species often showed little changes in the other species, or changed in the opposite direction. As evidence of adaptive divergence, both species performed better at their native site and better than the species from the other habitat. Adaptive divergence is, therefore, associated with the evolution of distinct plastic responses to environmental variation, despite these two species sharing a recent common ancestor.

Climate explains population divergence in drought-induced plasticity of functional traits and gene expression in a South African Protea

Long-term environmental variation often drives local adaptation and leads to trait differentiation across populations. Additionally, when traits change in an environment-dependent way through phenotypic plasticity, the genetic variation underlying plasticity will also be under selection. These processes could create a landscape of differentiation across populations in traits and their plasticity. Here, we performed a dry-down experiment under controlled conditions to measure responses in seedlings of a shrub species from the Cape Floristic Region, the common sugarbush (Protea repens). We measured morphological and physiological traits, and sequenced whole transcriptomes of leaf tissues from eight populations that represent both the climatic and the geographical distribution of this species. We found that there is substantial variation in how populations respond to drought, but we also observed common patterns such as reduced leaf size and leaf thickness, and up-regulation of stress-related and down-regulation of growth-related gene groups. Both high environmental heterogeneity and milder source site climates were associated with higher plasticity in various traits and co-expression gene networks. Associations between traits, trait plasticity, gene networks and the source site climate suggest that temperature may play a greater role in shaping these patterns when compared to precipitation, in line with recent changes in the region due to climate change. We also found that traits respond to climatic variation in an environment-dependent manner: some associations between traits and climate were apparent only under certain growing conditions. Together, our results uncover common responses of P. repens populations to drought, and climatic drivers of population differentiation in functional traits, gene expression and their plasticity.

Chromosome-Scale Assembly and Annotation of the Macadamia Genome (Macadamia integrifolia HAES 741)

Macadamia integrifolia is a representative of the large basal eudicot family Proteaceae and the main progenitor species of the Australian native nut crop macadamia. Since its commercialisation in Hawaii fewer than 100 years ago, global production has expanded rapidly. However, genomic resources are limited in comparison to other horticultural crops. The first draft assembly of M. integrifolia had good coverage of the functional gene space but its high fragmentation has restricted its use in comparative genomics and association studies. Here we have generated an improved assembly of cultivar HAES 741 (4,094 scaffolds, 745 Mb, N50 413 kb) using a combination of Illumina paired and PacBio long read sequences. Scaffolds were anchored to 14 pseudo-chromosomes using seven genetic linkage maps. This assembly has improved contiguity and coverage, with >120 Gb of additional sequence. Following annotation, 34,274 protein-coding genes were predicted, representing 90% of the expected gene content. Our results indicate that the macadamia genome is repetitive and heterozygous. The total repeat content was 55% and genome-wide heterozygosity, estimated by read mapping, was 0.98% or an average of one SNP per 102 bp. This is the first chromosome-scale genome assembly for macadamia and the Proteaceae. It is expected to be a valuable resource for breeding, gene discovery, conservation and evolutionary genomics.

Comparative transcriptomics of an arctic foundation species, tussock cottongrass (Eriophorum vaginatum), during an extreme heat event

Tussock cottongrass (Eriophorum vaginatum) is a foundation species for much of the arctic moist acidic tundra, which is currently experiencing extreme effects of climate change. The Arctic is facing higher summer temperatures and extreme weather events are becoming more common. We used Illumina RNA-Seq to analyse cDNA libraries for differential expression of genes from leaves of ecologically well-characterized ecotypes of tussock cottongrass found along a latitudinal gradient in the Alaskan Arctic and transplanted into a common garden. Plant sampling was performed on a typical summer day and during an extreme heat event. We obtained a de novo assembly that contained 423,353 unigenes. There were 363 unigenes up-regulated and 1,117 down-regulated among all ecotypes examined during the extreme heat event. Of these, 26 HSP unigenes had >log2-fold up-regulation. Several TFs associated with heat stress in previous studies were identified that had >log2-fold up- or down-regulation during the extreme heat event (e.g., DREB, NAC). There was consistent variation in DEGs among ecotypes, but not specifically related to whether plants originated from taiga or tundra ecosystems. As the climate changes it is essential to determine ecotypic diversity at the genomic level, especially for widespread species that impact ecosystem function.

Transcriptome sequencing of the apricot (Prunus armeniaca L.) and identification of differentially expressed genes involved in drought stress

Apricot (Prunus armeniaca L.) is an important fruit crop that is widely planted throughout the world. But drought affects both yield and quality of apricot. In order to study the effects of long-term drought on the molecular and physiological mechanisms of apricot, we used transcriptome sequencing and measured physiological indices. First, 322 million high-quality clean reads were obtained, and 74,892 unigenes were generated for the transcriptome. Among the assembled unigenes, 18,671 simple sequence repeats (SSRs) and 5581 differentially expressed genes (DEGs) were identified. The Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of the DEGs revealed that starch and sucrose metabolism, plant-pathogen interaction and plant hormone signal transduction pathways are enriched. Additionally, we used quantitative real-time PCR (qRT-PCR) to confirm the RNA-seq results with 11 drought-related DEGs. Second, through the physiological analysis of apricot leaves under constant drought stress, and the results show the internal microstructure of apricot leaves changed to withstand drought stress. At the same time, plants exposed to long-term drought stress showed higher degree of membrane damage, which reduced photosynthesis in the damaged leaves. Our findings enrich the genome resources for apricot and refine our understanding of the molecular and physiological mechanisms of drought response in this fruit crop, providing insights into drought adaptation of the apricot.

A population genomics appraisal suggests independent dispersals for bitter and sweet manioc in Brazilian Amazonia

Amazonia is a major world centre of plant domestication, but the genetics of domestication remains unclear for most Amazonian crops. Manioc (Manihot esculenta) is the most important staple food crop that originated in this region. Although manioc is relatively well-studied, little is known about the diversification of bitter and sweet landraces and how they were dispersed across Amazonia. We evaluated single nucleotide polymorphisms (SNPs) in wild and cultivated manioc to identify outlier SNPs putatively under selection and to assess the neutral genetic structure of landraces to make inferences about the evolution of the crop in Amazonia. Some outlier SNPs were in putative manioc genes possibly related to plant architecture, transcriptional regulation and responses to stress. The neutral SNPs revealed contrasting genetic structuring for bitter and sweet landraces. The outlier SNPs may be signatures of the genomic changes resulting from domestication, while the neutral genetic structure suggests independent dispersals for sweet and bitter manioc, possibly related to the earlier domestication and diversification of the former. Our results highlight the role of ancient peoples and current smallholders in the management and conservation of manioc genetic diversity, including putative genes and specific genetic resources with adaptive potential in the context of climate change.

Detailed Phenolic Characterization of Protea Pure and Hybrid Cultivars by Liquid Chromatography-Ion Mobility-High Resolution Mass Spectrometry (LC-IM-HR-MS)

In this study we report a detailed investigation of the polyphenol composition of Protea pure (P. cynaroides and P. neriifolia) and hybrid cultivars (Black beauty and Limelight). Aqueous methanol extracts of leaf and bract tissues were analyzed by ultrahigh pressure liquid chromatography hyphenated to photodiode array and ion mobility-high resolution mass spectrometric (UHPLC-PDA-IM-HR-MS) detection. A total of 67 metabolites were characterized based on their relative reversed phase (RP) retention, UV-vis spectra, low and high collision energy HR-MS data, and collisional cross section (CCS) values. These metabolites included 41 phenolic acid esters and 25 flavonoid derivatives, including 5 anthocyanins. In addition, an undescribed hydroxycinnamic acid-polygalatol ester, caffeoyl-O-polygalatol (1,5-anhydro-[6-O-caffeoyl]-sorbitol(glucitol)) was isolated and characterized by 1D and 2D NMR for the first time. This compound and its isomer are shown to be potential chemo-taxonomic markers.

Molecular responses to freshwater limitation in the mangrove tree Avicennia germinans (Acanthaceae)

Environmental variation along the geographical space can shape populations by natural selection. In the context of global warming and changing precipitation regimes, it is crucial to understand the role of environmental heterogeneity in tropical trees adaptation, given their disproportional contribution to water and carbon biogeochemical cycles. Here, we investigated how heterogeneity in freshwater availability along tropical wetlands has influenced molecular variations of the black mangrove (Avicennia germinans). A total of 57 trees were sampled at seven sites differing markedly in precipitation regime and riverine freshwater inputs. Using 2,297 genome-wide single nucleotide polymorphic markers, we found signatures of natural selection by the association between variations in allele frequencies and environmental variables, including the precipitation of the warmest quarter and the annual precipitation. Additionally, we found candidate loci for selection based on statistical deviations from neutral expectations of interpopulation differentiation. Most candidate loci within transcribed sequences were functionally associated with central aspects of drought tolerance or plant response to drought. Moreover, our results suggest the occurrence of the rapid evolution of a population, probably in response to sudden and persistent limitations in plant access to soil water, following a road construction in 1974. Observations supporting rapid evolution included the reduction in tree size and changes in allele frequencies and in transcript expression associated with increased drought tolerance through the accumulation of osmoprotectants and antioxidants, biosynthesis of cuticles, protection against protein degradation, stomatal closure, photorespiration and photosynthesis. We describe a major role of spatial heterogeneity in freshwater availability in the specialization of this typically tropical tree.

Seedling response to water stress in valley oak (Quercus lobata) is shaped by different gene networks across populations

Drought is a major stress for plants, creating a strong selection pressure for traits that enable plant growth and survival in dry environments. Many drought responses are conserved species-wide responses, while others vary among populations distributed across heterogeneous environments. We tested how six populations of the widely distributed California valley oak (Quercus lobata) sampled from contrasting climates would differ in their response to soil drying relative to well-watered controls in a common environment by measuring ecophysiological traits in 93 individuals and gene expression (RNA-seq) in 42 individuals. Populations did not differ in their adjustment of turgor loss point during soil drying, suggesting a generalized species-wide response. Differential expression analysis identified 689 genes with a common response to treatment across populations and 470 genes with population-specific responses. Weighted gene co-expression network analysis (WGCNA) identified groups of genes with similar expression patterns that may be regulated together (gene modules). Several gene modules responded differently to water stress among populations, suggesting regional differences in gene network regulation. Populations from sites with a high mean annual temperature responded to the imposed water stress with significantly greater changes in gene module expression, indicating that these populations may be locally adapted to respond to drought. We propose that this variation among valley oak populations provides a mechanism for differential tolerance to the increasingly frequent and severe droughts in California.

Unraveling the Developmental and Genetic Mechanisms Underpinning Floral Architecture in Proteaceae

Proteaceae are a basal eudicot family with a highly conserved floral groundplan but which displays considerable variation in other aspects of floral and inflorescence morphology. Their morphological diversity and phylogenetic position make them good candidates for understanding the evolution of floral architecture, in particular the question of the homology of the undifferentiated perianth with the differentiated perianth of core eudicots, and the mechanisms underlying the repeated evolution of zygomorphy. In this paper, we combine a morphological approach to explore floral ontogenesis and a transcriptomic approach to access the genes involved in floral organ identity and development, focusing on Grevillea juniperina, a species from subfamily Grevilleoideae. We present developmental data for Grevillea juniperina and three additional species that differ in their floral symmetry using stereomicroscopy, SEM and High Resolution X-Ray Computed Tomography. We find that the adnation of stamens to tepals takes place at early developmental stages, and that the establishment of bilateral symmetry coincides with the asymmetrical growth of the single carpel. To set a framework for understanding the genetic basis of floral development in Proteaceae, we generated and annotated de novo a reference leaf/flower transcriptome from Grevillea juniperina. We found Grevillea homologs of all lineages of MADS-box genes involved in floral organ identity. Using Arabidopsis thaliana gene expression data as a reference, we found homologs of other genes involved in floral development in the transcriptome of G. juniperina. We also found at least 21 class I and class II TCP genes, a gene family involved in the regulation of growth processes, including floral symmetry. The expression patterns of a set of floral genes obtained from the transcriptome were characterized during floral development to assess their organ specificity and asymmetry of expression.

Evolutionary Transcriptomics and Proteomics: Insight into Plant Adaptation

Comparative transcriptomics and proteomics (T&P) have brought biological insight into development, gene function, and physiological stress responses. However, RNA-seq and high-throughput proteomics remain underutilised in studies of plant adaptation. These methodologies have created discovery tools with the potential to significantly advance our understanding of adaptive diversification. We outline experimental recommendations for their application. We discuss analysis models and approaches that accelerate the identification of adaptive gene sets and integrate transcriptome, proteome, phenotypic, and environmental data. Finally, we encourage widespread uptake and future developments in T&P that will advance our understanding of evolution and adaptation.

Isolation by distance and isolation by environment contribute to population differentiation in Protea repens (Proteaceae L.), a widespread South African species

PREMISE OF THE STUDY

The Cape Floristic Region (CFR) of South Africa is renowned for its botanical diversity, but the evolutionary origins of this diversity remain controversial. Both neutral and adaptive processes have been implicated in driving diversification, but population-level studies of plants in the CFR are rare. Here, we investigate the limits to gene flow and potential environmental drivers of selection in Protea repens L. (Proteaceae L.), a widespread CFR species.

METHODS

We sampled 19 populations across the range of P. repens and used genotyping by sequencing to identify 2066 polymorphic loci in 663 individuals. We used a Bayesian F_ST outlier analysis to identify single-nucleotide polymorphisms (SNPs) marking genomic regions that may be under selection; we used those SNPs to identify potential drivers of selection and excluded them from analyses of gene flow and genetic structure.

RESULTS

A pattern of isolation by distance suggested limited gene flow between nearby populations. The populations of P. repens fell naturally into two or three groupings, which corresponded to an east-west split. Differences in rainfall seasonality contributed to diversification in highly divergent loci, as do barriers to gene flow that have been identified in other species.

CONCLUSIONS

The strong pattern of isolation by distance is in contrast to the findings in the only other widespread species in the CFR that has been similarly studied, while the effects of rainfall seasonality are consistent with well-known patterns. Assessing the generality of these results will require investigations of other CFR species.

Anchored phylogenomics improves the resolution of evolutionary relationships in the rapid radiation of Protea L

PREMISE OF THE STUDY

Estimating phylogenetic relationships in relatively recent evolutionary radiations is challenging, especially if short branches associated with recent divergence result in multiple gene tree histories. We combine anchored enrichment next-generation sequencing with species tree analyses to produce a robust estimate of phylogenetic relationships in the genus Protea (Proteaceae), an iconic radiation in South Africa.

METHODS

We sampled multiple individuals within 59 out of 112 species of Protea and 6 outgroup species for a total of 163 individuals, and obtained sequences for 498 low-copy, orthologous nuclear loci using anchored phylogenomics. We compare several approaches for building species trees, and explore gene tree-species tree discrepancies to determine whether poor phylogenetic resolution reflects a lack of informative sites, incomplete lineage sorting, or hybridization.

KEY RESULTS

Phylogenetic estimates from species tree approaches are similar to one another and recover previously well-supported clades within Protea, in addition to providing well-supported phylogenetic hypotheses for previously poorly resolved intrageneric relationships. Individual gene trees are markedly different from one another and from species trees. Nonetheless, analyses indicate that differences among gene trees occur primarily concerning clades supported by short branches.

CONCLUSIONS

Species tree methods using hundreds of nuclear loci provided strong support for many previously unresolved relationships in the radiation of the genus Protea. In cases where support for particular relationships remains low, these appear to arise from few informative sites and lack of information rather than strongly supported disagreement among gene trees.