Sarracenia pitcher plant-associated microbial communities differ primarily by host species across a longitudinal gradient.

Plant-associated microbial communities can profoundly affect plant health and success, and research is still uncovering factors driving the assembly of these communities. Here, we examine how geography versus host species affects microbial community structure and differential abundances of individual taxa. We use metabarcoding to characterize the bacteria and eukaryotes associated with five, often co-occurring species of Sarracenia pitcher plants (Sarraceniaceae) and three natural hybrids along the longitudinal gradient of the U.S. Gulf Coast, as well as samples from S. purpurea in Massachusetts. To tease apart the effects of geography versus host species, we focus first on sites with co-occurring species and then on species located across different sites. Our analyses show that bacterial and eukaryotic community structures are clearly and consistently influenced by host species identity, with geographic factors also playing a role. Naturally occurring hybrids appear to also host unique communities, which are in some ways intermediate between their parent species. We see significant effects of geography (site and longitude), but these generally explain less of the variation among pitcher communities. Overall, in Sarracenia pitchers, host plant phenotype significantly affects the pitcher microbiomes and other associated organisms.

Interactions between strains govern the eco-evolutionary dynamics of microbial communities.

Genomic data has revealed that genotypic variants of the same species, that is, strains, coexist and are abundant in natural microbial communities. However, it is not clear if strains are ecologically equivalent, and at what characteristic genetic distance they might exhibit distinct interactions and dynamics. Here, we address this problem by tracking 10 taxonomically diverse microbial communities from the pitcher plant Sarracenia purpurea in the laboratory for more than 300 generations. Using metagenomic sequencing, we reconstruct their dynamics over time and across scales, from distant phyla to closely related genotypes. We find that most strains are not ecologically equivalent and exhibit distinct dynamical patterns, often being significantly more correlated with strains from another species than their own. Although even a single mutation can affect laboratory strains, on average, natural strains typically decouple in their dynamics beyond a genetic distance of 100 base pairs. Using mathematical consumer-resource models, we show that these taxonomic patterns emerge naturally from ecological interactions between community members, but only if the interactions are coarse-grained at the level of strains, not species. Finally, by analyzing genomic differences between strains, we identify major functional hubs such as transporters, regulators, and carbohydrate-catabolizing enzymes, which might be the basis for strain-specific interactions. Our work suggests that fine-scale genetic differences in natural communities could be created and stabilized via the rapid diversification of ecological interactions between strains.

Transcriptome-wide shift from photosynthesis and energy metabolism upon endogenous fluid protein depletion in young Nepenthes ampullaria pitchers.

Carnivorous pitcher plants produce specialised pitcher organs containing secretory glands, which secrete acidic fluids with hydrolytic enzymes for prey digestion and nutrient absorption. The content of pitcher fluids has been the focus of many fluid protein profiling studies. These studies suggest an evolutionary convergence of a conserved group of similar enzymes in diverse families of pitcher plants. A recent study showed that endogenous proteins were replenished in the pitcher fluid, which indicates a feedback mechanism in protein secretion. This poses an interesting question on the physiological effect of plant protein loss. However, there is no study to date that describes the pitcher response to endogenous protein depletion. To address this gap of knowledge, we previously performed a comparative RNA-sequencing experiment of newly opened pitchers (D0) against pitchers after 3 days of opening (D3C) and pitchers with filtered endogenous proteins (>10 kDa) upon pitcher opening (D3L). Nepenthes ampullaria was chosen as a model study species due to their abundance and unique feeding behaviour on leaf litters. The analysis of unigenes with top 1% abundance found protein translation and stress response to be overrepresented in D0, compared to cell wall related, transport, and signalling for D3L. Differentially expressed gene (DEG) analysis identified DEGs with functional enrichment in protein regulation, secondary metabolism, intracellular trafficking, secretion, and vesicular transport. The transcriptomic landscape of the pitcher dramatically shifted towards intracellular transport and defence response at the expense of energy metabolism and photosynthesis upon endogenous protein depletion. This is supported by secretome, transportome, and transcription factor analysis with RT-qPCR validation based on independent samples. This study provides the first glimpse into the molecular responses of pitchers to protein loss with implications to future cost/benefit analysis of carnivorous pitcher plant energetics and resource allocation for adaptation in stochastic environments.

DNA metabarcoding reveals changes in the contents of carnivorous plants along an elevation gradient.

Resource variation along abiotic gradients influences subsequent trophic interactions and these effects can be transmitted through entire food webs. Interactions along abiotic gradients can provide clues as to how organisms will face changing environmental conditions, such as future range shifts. However, it is challenging to find replicated systems to study these effects. Phytotelmata, such as those found in carnivorous plants, are isolated aquatic communities and thus form a good model for the study of replicated food webs. Due to the degraded nature of the prey, molecular techniques provide a useful tool to study these communities. We studied the pitcher plant Sarracenia purpurea L. in allochthonous populations along an elevational gradient in the Alps and Jura. We predicted that invertebrate richness in the contents of the pitcher plants would decrease with increasing elevation, reflecting harsher environmental conditions. Using metabarcoding of the COI gene, we sequenced the invertebrate contents of these pitcher plants. We assigned Molecular Operational Taxonomic Units at ordinal level as well as recovering species-level data. We found small but significant changes in community composition with elevation. These recovered sequences could belong to invertebrate prey, rotifer inquilines, pollinators and other animals possibly living inside the pitchers. However, we found no directional trend or site-based differences in MOTU richness with elevational gradient. Use of molecular techniques for dietary or contents analysis is a powerful way to examine numerous degraded samples, although factors such as DNA persistence and the relationship with species presence still have to be completely determined.

Convergence between the microcosms of Southeast Asian and North American pitcher plants.

The 'pitchers' of carnivorous pitcher plants are exquisite examples of convergent evolution. An open question is whether the living communities housed in pitchers also converge in structure or function. Using samples from more than 330 field-collected pitchers of eight species of Southeast Asian Nepenthes and six species of North American Sarracenia, we demonstrate that the pitcher microcosms, or miniature ecosystems with complex communities, are strikingly similar. Compared to communities from surrounding habitats, pitcher communities house fewer species. While communities associated with the two genera contain different microbial organisms and arthropods, the species are predominantly from the same phylogenetic clades. Microbiomes from both genera are enriched in degradation pathways and have high abundances of key degradation enzymes. Moreover, in a manipulative field experiment, Nepenthes pitchers placed in a North American bog assembled Sarracenia-like communities. An understanding of the convergent interactions in pitcher microcosms facilitates identification of selective pressures shaping the communities.

Metabolite profiling of the carnivorous pitcher plants Darlingtonia and Sarracenia.

Sarraceniaceae is a New World carnivorous plant family comprising three genera: Darlingtonia, Heliamphora, and Sarracenia. The plants occur in nutrient-poor environments and have developed insectivorous capability in order to supplement their nutrient uptake. Sarracenia flava contains the alkaloid coniine, otherwise only found in Conium maculatum, in which its biosynthesis has been studied, and several Aloe species. Its ecological role and biosynthetic origin in S. flava is speculative. The aim of the current research was to investigate the occurrence of coniine in Sarracenia and Darlingtonia and to identify common constituents of both genera, unique compounds for individual variants and floral scent chemicals. In this comprehensive metabolic profiling study, we looked for compound patterns that are associated with the taxonomy of Sarracenia species. In total, 57 different Sarracenia and D. californica accessions were used for metabolite content screening by gas chromatography-mass spectrometry. The resulting high-dimensional data were studied using a data mining approach. The two genera are characterized by a large number of metabolites and huge chemical diversity between different species. By applying feature selection for clustering and by integrating new biochemical data with existing phylogenetic data, we were able to demonstrate that the chemical composition of the species can be explained by their known classification. Although transcriptome analysis did not reveal a candidate gene for coniine biosynthesis, the use of a sensitive selected ion monitoring method enabled the detection of coniine in eight Sarracenia species, showing that it is more widespread in this genus than previously believed.

Phylogenetic inferences of Nepenthes species in Peninsular Malaysia revealed by chloroplast (trnL intron) and nuclear (ITS) DNA sequences.

BACKGROUND: The chloroplastic trnL intron and the nuclear internal transcribed spacer (ITS) region were sequenced for 11 Nepenthes species recorded in Peninsular Malaysia to examine their phylogenetic relationship and to evaluate the usage of trnL intron and ITS sequences for phylogenetic reconstruction of this genus. RESULTS: Phylogeny reconstruction was carried out using neighbor-joining, maximum parsimony and Bayesian analyses. All the trees revealed two major clusters, a lowland group consisting of N. ampullaria, N. mirabilis, N. gracilis and N. rafflesiana, and another containing both intermediately distributed species (N. albomarginata and N. benstonei) and four highland species (N. sanguinea, N. macfarlanei, N. ramispina and N. alba). CONCLUSIONS: The trnL intron and ITS sequences proved to provide phylogenetic informative characters for deriving a phylogeny of Nepenthes species in Peninsular Malaysia. To our knowledge, this is the first molecular phylogenetic study of Nepenthes species occurring along an altitudinal gradient in Peninsular Malaysia.

Phylogeographic concordance factors quantify phylogeographic congruence among co-distributed species in the Sarracenia alata pitcher plant system.

Comparative phylogeographic investigations have identified congruent phylogeographic breaks in co-distributed species in nearly every region of the world. The qualitative assessments of phylogeographic patterns traditionally used to identify such breaks, however, are limited because they rely on identifying monophyletic groups across species and do not account for coalescent stochasticity. Only long-standing phylogeographic breaks are likely to be obvious; many species could have had a concerted response to more recent landscape events, yet possess subtle signs of phylogeographic congruence because ancestral polymorphism has not completely sorted. Here, we introduce Phylogeographic Concordance Factors (PCFs), a novel method for quantifying phylogeographic congruence across species. We apply this method to the Sarracenia alata pitcher plant system, a carnivorous plant with a diverse array of commensal organisms. We explore whether a group of ecologically associated arthropods have co-diversified with the host pitcher plant, and identify if there is a positive correlation between ecological interaction and PCFs. Results demonstrate that multiple arthropods share congruent phylogeographic breaks with S. alata, and provide evidence that the level of ecological association can be used to predict the degree of similarity in the phylogeographic pattern. This study outlines an approach for quantifying phylogeographic congruence, a central concept in biogeographic research.

Development of Transcriptomic Markers for Population Analysis Using Restriction Site Associated RNA Sequencing (RARseq).

We describe restriction site associated RNA sequencing (RARseq), an RNAseq-based genotype by sequencing (GBS) method. It includes the construction of RNAseq libraries from double stranded cDNA digested with selected restriction enzymes. To test this, we constructed six single- and six-dual-digested RARseq libraries from six F2 pitcher plant individuals and sequenced them on a half of a Miseq run. On average, the de novo approach of population genome analysis detected 544 and 570 RNA SNPs, whereas the reference transcriptome-based approach revealed an average of 1907 and 1876 RNA SNPs per individual, from single- and dual-digested RARseq data, respectively. The average numbers of RNA SNPs and alleles per loci are 1.89 and 2.17, respectively. Our results suggest that the RARseq protocol allows good depth of coverage per loci for detecting RNA SNPs and polymorphic loci for population genomics and mapping analyses. In non-model systems where complete genomes sequences are not always available, RARseq data can be analyzed in reference to the transcriptome. In addition to enriching for functional markers, this method may prove particularly useful in organisms where the genomes are not favorable for DNA GBS.

Structural and functional characteristics of S-like ribonucleases from carnivorous plants.

Although the S-like ribonucleases (RNases) share sequence homology with the S-RNases involved in the self-incompatibility mechanism in plants, they are not associated with this mechanism. They usually function in stress responses in non-carnivorous plants and in carnivory in carnivorous plants. In this study, we clarified the structures of the S-like RNases of Aldrovanda vesiculosa, Nepenthes bicalcarata and Sarracenia leucophylla, and compared them with those of other plants. At ten positions, amino acid residues are conserved or almost conserved only for carnivorous plants (six in total). In contrast, two positions are specific to non-carnivorous plants. A phylogenetic analysis revealed that the S-like RNases of the carnivorous plants form a group beyond the phylogenetic relationships of the plants. We also prepared and characterized recombinant S-like RNases of Dionaea muscipula, Cephalotus follicularis, A. vesiculosa, N. bicalcarata and S. leucophylla, and RNS1 of Arabidopsis thaliana. The recombinant carnivorous plant enzymes showed optimum activities at about pH 4.0. Generally, poly(C) was digested less efficiently than poly(A), poly(I) and poly(U). The kinetic parameters of the recombinant D. muscipula enzyme (DM-I) and A. thaliana enzyme RNS1 were similar. The k cat/K m of recombinant RNS1 was the highest among the enzymes, followed closely by that of recombinant DM-I. On the other hand, the k cat/K m of the recombinant S. leucophylla enzyme was the lowest, and was ~1/30 of that for recombinant RNS1. The magnitudes of the k cat/K m values or k cat values for carnivorous plant S-like RNases seem to correlate negatively with the dependency on symbionts for prey digestion.

Hybridization leads to interspecific gene flow in Sarracenia (Sarraceniaceae).

PREMISE OF THE STUDY: The role of hybridization in plant evolution remains a source of intense debate. Potential consequences range from genetic dead-ends to species fusion or hybrid speciation. While much has been learned from model systems such as Populus, Iris, and Helianthus, many questions remain. Consisting of 11 species that are all capable of hybridizing, Sarracenia presents an excellent system in which to study hybridization. • METHODS: Using microsatellites, we examined a single field site consisting of three species: S. leucophylla, S. alata, and S. rubra subsp. wherryi. We determined the level of genetic admixture and compared it with allopatric sites of the same taxa. • KEY RESULTS: In contrast to the well-defined clusters formed when we examined the allopatric sites, the sympatric field site exhibited a wide range of admixture. Additionally, when the relative genetic makeup of "pure" species at the site was compared with the makeup of hybrids, we found that Sarracenia alata contributed disproportionately to the hybrid genomes. • CONCLUSIONS: Our study provides further evidence that hybridization is contributing to interspecific gene flow in the genus and that all species do not contribute equally to hybridization. Implications for conservation are discussed.

Genetic diversity in three endangered pitcher plant species (Sarracenia; Sarraceniaceae) is lower than widespread congeners.

PREMISE OF THE STUDY: Narrow-ranging, rare species often exhibit levels of genetic diversity lower than more common or widespread congeners. These taxa are at increased risk of extinction due to threats associated with natural as well as anthropogenic events. We assessed genetic variation in three federally endangered Sarracenia species. We discuss maintenance of genetic diversity and evolutionary implications of rarity. • METHODS: We analyzed three noncoding chloroplast regions and nine microsatellite loci in populations spanning the geographic ranges of S. oreophila, S. alabamensis, and S. jonesii. The same microsatellite loci were used to examine a single field site of three more widespread species (S. alata, S. leucophylla, and S. rubra subsp. wherryi). • KEY RESULTS: All three endangered species have experienced reductions in population size and numbers. All show considerably less variation than more widespread members of the genus. Sarracenia alabamensis maintains the greatest microsatellite variation but has the fewest remaining populations and may be under the greatest threat. More widespread S. oreophila maintains surprising chloroplast diversity, yet exhibits little microsatellite diversity. Sarracenia jonesii lacks chloroplast diversity, yet maintains greater microsatellite diversity than S. oreophila. • CONCLUSIONS: The three endangered species differ in levels and structure of diversity, yet not in predictable ways, emphasizing that unique demographic and ecological histories, rather than current distribution and population size, best explain present patterns of genetic variation. Maintenance of remaining genetic variation is important, but preventing further habitat loss and degradation is critical.

S-like ribonuclease gene expression in carnivorous plants.

Functions of S-like ribonucleases (RNases) differ considerably from those of S-RNases that function in self-incompatibility. Expression of S-like RNases is usually induced by low nutrition, vermin damage or senescence. However, interestingly, an Australian carnivorous plant Drosera adelae (a sundew), which traps prey with a sticky digestive liquid, abundantly secretes an S-like RNase DA-I in the digestive liquid even in ordinary states. Here, using D. adelae, Dionaea muscipula (Venus flytrap) and Cephalotus follicularis (Australian pitcher plant), we show that carnivorous plants use S-like RNases for carnivory: the gene da-I encoding DA-I and its ortholog cf-I of C. follicularis are highly expressed and constitutively active in each trap/digestion organ, while the ortholog dm-I of D. muscipula becomes highly active after trapping insects. The da-I promoter is unmethylated only in its trap/digestion organ, glandular tentacles (which comprise a small percentage of the weight of the whole plant), but methylated in other organs, which explains the glandular tentacles-specific expression of the gene and indicates a very rare gene regulation system. In contrast, the promoters of dm-I, which shows induced expression, and cf-I, which has constitutive expression, were not methylated in any organs examined. Thus, it seems that the regulatory mechanisms of the da-I, dm-I and cf-I genes differ from each other and do not correlate with the phylogenetic relationship. The current study suggests that under environmental pressure in specific habitats carnivorous plants have managed to evolve their S-like RNase genes to function in carnivory.

Deep phylogeographic structure and environmental differentiation in the carnivorous plant Sarracenia alata.

We collected ~29 kb of sequence data using Roche 454 pyrosequencing in order to estimate the timing and pattern of diversification in the carnivorous pitcher plant Sarracenia alata. Utilizing modified protocols for reduced representation library construction, we generated sequence data from 86 individuals across 10 populations from throughout the range of the species. We identified 76 high-quality and high-coverage loci (containing over 500 SNPs) using the bioinformatics pipeline PRGmatic. Results from a Bayesian clustering analysis indicate that populations are highly structured, and are similar in pattern to the topology of a population tree estimated using *BEAST. The pattern of diversification within Sarracenia alata implies that riverine barriers are the primary factor promoting population diversification, with divergence across the Mississippi River occurring more than 60,000 generations before present. Further, significant patterns of niche divergence and the identification of several outlier loci suggest that selection may contribute to population divergence. Our results demonstrate the feasibility of using next-generation sequencing to investigate intraspecific genetic variation in nonmodel species.

Transcriptome analysis of sarracenia, an insectivorous plant.

Sarracenia species (pitcher plants) are carnivorous plants which obtain a portion of their nutrients from insects captured in the pitchers. To investigate these plants, we sequenced the transcriptome of two species, Sarracenia psittacina and Sarracenia purpurea, using Roche 454 pyrosequencing technology. We obtained 46 275 and 36 681 contigs by de novo assembly methods for S. psittacina and S. purpurea, respectively, and further identified 16 163 orthologous contigs between them. Estimation of synonymous substitution rates between orthologous and paralogous contigs indicates the events of genome duplication and speciation within the Sarracenia genus both occurred ∼2 million years ago. The ratios of synonymous and non-synonymous substitution rates indicated that 491 contigs have been under positive selection (K(a)/K(s) > 1). Significant proportions of these contigs were involved in functions related to binding activity. We also found that the greatest sequence similarity for both of these species was to Vitis vinifera, which is most consistent with a non-current classification of the order Ericales as an asterid. This study has provided new insights into pitcher plants and will contribute greatly to future research on this genus and its distinctive ecological adaptations.

The microbial Phyllogeography of the carnivorous plant Sarracenia alata.

Carnivorous pitcher plants host diverse microbial communities. This plant-microbe association provides a unique opportunity to investigate the evolutionary processes that influence the spatial diversity of microbial communities. Using next-generation sequencing of environmental samples, we surveyed microbial communities from 29 pitcher plants (Sarracenia alata) and compare community composition with plant genetic diversity in order to explore the influence of historical processes on the population structure of each lineage. Analyses reveal that there is a core S. alata microbiome, and that it is similar in composition to animal gut microfaunas. The spatial structure of community composition in S. alata (phyllogeography) is congruent at the deepest level with the dominant features of the landscape, including the Mississippi river and the discrete habitat boundaries that the plants occupy. Intriguingly, the microbial community structure reflects the phylogeographic structure of the host plant, suggesting that the phylogenetic structure of bacterial communities and population genetic structure of their host plant are influenced by similar historical processes.

Proteome analysis of pitcher fluid of the carnivorous plant Nepenthes alata.

The genus Nepenthes comprises carnivorous plants that digest insects in pitcher fluid to supplement their nitrogen uptake. In a recent study, two acid proteinases (nepenthesins I and II) were purified from the pitcher fluid. However, no other enzymes involved in prey digestion have been identified, although several enzyme activities have been reported. To identify all the proteins involved, we performed a proteomic analysis of Nepenthes pitcher fluid. The secreted proteins in pitcher fluid were separated by sodium dodecyl sulfate-polyacrylamide gel electrophoresis, and several protein bands were detected by silver staining. The proteins were identified by in-gel tryptic digestion, de novo peptide sequencing, and homology searches against public databases. The proteins included homologues of beta-D-xylosidase, beta-1,3-glucanase, chitinase, and thaumatin-like protein, most of which are designated "pathogenesis-related proteins". These proteins presumably inhibit bacterial growth in the pitcher fluid to ensure sufficient nutrients for Nepenthes growth.