Jessy Labbé

A Multifactor Analysis of Fungal and Bacterial Community Structure in the Root Microbiome of Mature Populus deltoides Trees

Bacterial and fungal communities associated with plant roots are central to the host health, survival and growth. However, a robust understanding of the root-microbiome and the factors that drive host associated microbial community structure have remained elusive, especially in mature perennial plants from natural settings. Here, we investigated relationships of bacterial and fungal communities in the rhizosphere and root endosphere of the riparian tree species Populus deltoides, and the influence of soil parameters, environmental properties (host phenotype and aboveground environmental settings), host plant genotype (Simple Sequence Repeat (SSR) markers), season (Spring vs. Fall) and geographic setting (at scales from regional watersheds to local riparian zones) on microbial community structure. Each of the trees sampled displayed unique aspects to its associated community structure with high numbers of Operational Taxonomic Units (OTUs) specific to an individual trees (bacteria >90%, fungi >60%). Over the diverse conditions surveyed only a small number of OTUs were common to all samples within rhizosphere (35 bacterial and 4 fungal) and endosphere (1 bacterial and 1 fungal) microbiomes. As expected, Proteobacteria and Ascomycota were dominant in root communities (>50%) while other higher-level phylogenetic groups (Chytridiomycota, Acidobacteria) displayed greatly reduced abundance in endosphere compared to the rhizosphere. Variance partitioning partially explained differences in microbiome composition between all sampled roots on the basis of seasonal and soil properties (4% to 23%). While most variation remains unattributed, we observed significant differences in the microbiota between watersheds (Tennessee vs. North Carolina) and seasons (Spring vs. Fall). SSR markers clearly delineated two host populations associated with the samples taken in TN vs. NC, but overall host genotypic distances did not have a significant effect on corresponding communities that could be separated from other measured effects.

Extensive gene flow over Europe and possible speciation over Eurasia in the ectomycorrhizal basidiomycete Laccaria amethystina complex

Biogeographical patterns and large‐scale genetic structure have been little studied in ectomycorrhizal (EM) fungi, despite the ecological and economic importance of EM symbioses. We coupled population genetics and phylogenetic approaches to understand spatial structure in fungal populations on a continental scale. Using nine microsatellite markers, we characterized gene flow among 16 populations of the widespread EM basidiomycete Laccaria amethystina over Europe (i.e. over 2900 km). We also widened our scope to two additional populations from Japan (104 km away) and compared them with European populations through microsatellite markers and multilocus phylogenies, using three nuclear genes (NAR, G6PD and ribosomal DNA) and two mitochondrial ribosomal genes. European L. amethystina populations displayed limited differentiation (average FST = 0.041) and very weak isolation by distance (IBD). This panmictic European pattern may result from effective aerial dispersal of spores, high genetic diversity in populations and mutualistic interactions with multiple hosts that all facilitate migration. The multilocus phylogeny based on nuclear genes confirmed that Japanese and European specimens were closely related but clustered on a geographical basis. By using microsatellite markers, we found that Japanese populations were strongly differentiated from the European populations (FST = 0.416), more than expected by extrapolating the European pattern of IBD. Population structure analyses clearly separated the populations into two clusters, i.e. European and Japanese clusters. We discuss the possibility of IBD in a continuous population (considering some evidence for a ring species over the Northern Hemisphere) vs. an allopatric speciation over Eurasia, making L. amethystina a promising model of intercontinental species for future studies.

Extensive gene flow over Europe and possible speciation

Biogeographical patterns and large‐scale genetic structure have been little studied in ectomycorrhizal (EM) fungi, despite the ecological and economic importance of EM symbioses. We coupled population genetics and phylogenetic approaches to understand spatial structure in fungal populations on a continental scale. Using nine microsatellite markers, we characterized gene flow among 16 populations of the widespread EM basidiomycete Laccaria amethystina over Europe (i.e. over 2900 km). We also widened our scope to two additional populations from Japan (104 km away) and compared them with European populations through microsatellite markers and multilocus phylogenies, using three nuclear genes (NAR, G6PD and ribosomal DNA) and two mitochondrial ribosomal genes. European L. amethystina populations displayed limited differentiation (average FST = 0.041) and very weak isolation by distance (IBD). This panmictic European pattern may result from effective aerial dispersal of spores, high genetic diversity in populations and mutualistic interactions with multiple hosts that all facilitate migration. The multilocus phylogeny based on nuclear genes confirmed that Japanese and European specimens were closely related but clustered on a geographical basis. By using microsatellite markers, we found that Japanese populations were strongly differentiated from the European populations (FST = 0.416), more than expected by extrapolating the European pattern of IBD. Population structure analyses clearly separated the populations into two clusters, i.e. European and Japanese clusters. We discuss the possibility of IBD in a continuous population (considering some evidence for a ring species over the Northern Hemisphere) vs. an allopatric speciation over Eurasia, making L. amethystina a promising model of intercontinental species for future studies.

Survey and analysis of simple sequence repeats in the Laccaria bicolor genome, with development of microsatellite markers

It is becoming clear that simple sequence repeats (SSRs) play a significant role in fungal genome organization, and they are a large source of genetic markers for population genetics and meiotic maps. We identified SSRs in the Laccaria bicolor genome by in silico survey and analyzed their distribution in the different genomic regions. We also compared the abundance and distribution of SSRs in L. bicolor with those of the following fungal genomes: Phanerochaete chrysosporium, Coprinopsis cinerea, Ustilago maydis, Cryptococcus neoformans, Aspergillus nidulans, Magnaporthe grisea, Neurospora crassa and Saccharomyces cerevisiae. Using the MISA computer program, we detected 277,062 SSRs in the L. bicolor genome representing 8% of the assembled genomic sequence. Among the analyzed basidiomycetes, L. bicolor exhibited the highest SSR density although no correlation between relative abundance and the genome sizes was observed. In most genomes the short motifs (mono- to trinucleotides) were more abundant than the longer repeated SSRs. Generally, in each organism, the occurrence, relative abundance, and relative density of SSRs decreased as the repeat unit increased. Furthermore, each organism had its own common and longest SSRs. In the L. bicolor genome, most of the SSRs were located in intergenic regions (73.3%) and the highest SSR density was observed in transposable elements (TEs; 6,706 SSRs/Mb). However, 81% of the protein-coding genes contained SSRs in their exons, suggesting that SSR polymorphism may alter gene phenotypes. Within a L. bicolor offspring, sequence polymorphism of 78 SSRs was mainly detected in non-TE intergenic regions. Unlike previously developed microsatellite markers, these new ones are spread throughout the genome; these markers could have immediate applications in population genetics.

Newly identified helper bacteria stimulate ectomycorrhizal formation in Populus

Mycorrhiza helper bacteria (MHB) are known to increase host root colonization by mycorrhizal fungi but the molecular mechanisms and potential tripartite interactions are poorly understood. Through an effort to study Populus microbiome, we isolated 21 Pseudomonas strains from native Populus deltoides roots. These bacterial isolates were characterized and screened for MHB effectiveness on the Populus-Laccaria system. Two additional Pseudomonas strains (i.e., Pf-5 and BBc6R8) from existing collections were included for comparative purposes. We analyzed the effect of co-cultivation of these 23 individual Pseudomonas strains on Laccaria bicolor “S238N” growth rate, mycelial architecture and transcriptional changes. Nineteen of the 23 Pseudomonas strains tested had positive effects on L. bicolor S238N growth, as well as on mycelial architecture, with strains GM41 and GM18 having the most significant effect. Four of seven L. bicolor reporter genes, Tra1, Tectonin2, Gcn5, and Cipc1, thought to be regulated during the interaction with MHB strain BBc6R8, were induced or repressed, while interacting with Pseudomonas strains GM17, GM33, GM41, GM48, Pf-5, and BBc6R8. Strain GM41 promoted the highest roots colonization across three Populus species but most notably in P. deltoides, which is otherwise poorly colonized by L. bicolor. Here we report novel MHB strains isolated from native Populus that improve L. bicolor root colonization on Populus. This tripartite relationship could be exploited for Populus species/genotypes nursery production as a means of improving establishment and survival in marginal lands.

The obscure events contributing to the evolution of an incipient sex chromosome in Populus: a retrospective working hypothesis

Genetic determination of gender is a fundamental developmental and evolutionary process in plants. Although it appears that dioecy in Populus is genetically controlled, the precise gender-determining systems remain unclear. The recently released second draft assembly and annotated gene set of the Populus genome provided an opportunity to revisit this topic. We hypothesized that over evolutionary time, selective pressure has reformed the genome structure and gene composition in the peritelomeric region of the chromosome XIX, which has resulted in a distinctive genome structure and cluster of genes contributing to gender determination in Populus trichocarpa. Multiple lines of evidence support this working hypothesis. First, the peritelomeric region of the chromosome XIX contains significantly fewer single nucleotide polymorphisms than the rest of Populus genome and has a distinct evolutionary history. Second, the peritelomeric end of chromosome XIX contains the largest cluster of the nucleotide-binding site–leucine-rich repeat (NBS–LRR) class of disease resistance genes in the entire Populus genome. Third, there is a high occurrence of small microRNAs on chromosome XIX, which is coincident to the region containing the putative gender-determining locus and the major cluster of NBS–LRR genes. Further, by analyzing the metabolomic profiles of floral bud in male and female Populus trees using a gas chromatography-mass spectrometry, we found that there are gender-specific accumulations of phenolic glycosides. Taken together, these findings led to the hypothesis that resistance to and regulation of a floral pathogen and gender determination coevolved, and that these events triggered the emergence of a nascent sex chromosome. Further studies of chromosome XIX will provide new insights into the genetic control of gender determination in Populus.

Specialized Microbiome of a Halophyte and its Role in Helping Non-Host Plants to Withstand Salinity

Root microbiota is a crucial determinant of plant productivity and stress tolerance. Here, we hypothesize that the superior halo-tolerance of seepweed Suaeda salsa is tightly linked to a specialized belowground microbiome. To test this hypothesis, we performed a phylogenetic trait-based framework analysis based on bacterial 16S rRNA gene and fungal nuclear rRNA internal transcribed spacer profiling. Data showed that the dominant α-proteobacteria and γ-proteobacteria communities in bulk soil and root endosphere tend to be phylogenetically clustered and at the same time exhibit phylogenetic over-dispersion in rhizosphere. Likewise, the dominant fungal genera occurred at high phylogenetic redundancy. Interestingly, we found the genomes of rhizospheric and endophytic bacteria associated with S. salsa to be enriched in genes contributing to salt stress acclimatization, nutrient solubilization and competitive root colonization. A wide diversity of rhizobacteria with similarity to known halotolerant taxa further supported this interpretation. These findings suggest that an ecological patterned root-microbial interaction strategy has been adopted in S. salsa system to confront soil salinity. We also demonstrated that the potential core microbiome members improve non-host plants growth and salt tolerance. This work provides a platform to improve plant fitness with halophytes-microbial associates and novel insights into the functions of plant microbiome under salinity.

Highly Efficient Isolation of Populus Mesophyll Protoplasts and Its Application in Transient Expression Assays

Background Populus is a model woody plant and a promising feedstock for lignocellulosic biofuel production. However, its lengthy life cycle impedes rapid characterization of gene function. Methodology/Principal Findings We optimized a Populus leaf mesophyll protoplast isolation protocol and established a Populus protoplast transient expression system. We demonstrated that Populus protoplasts are able to respond to hormonal stimuli and that a series of organelle markers are correctly localized in the Populus protoplasts. Furthermore, we showed that the Populus protoplast transient expression system is suitable for studying protein-protein interaction, gene activation, and cellular signaling events. Conclusions/Significance This study established a method for efficient isolation of protoplasts from Populus leaf and demonstrated the efficacy of using Populus protoplast transient expression assays as an in vivo system to characterize genes and pathways.

Identification of quantitative trait loci affecting ectomycorrhizal symbiosis in an interspecific F1 poplar cross and differential expression of genes in ectomycorrhizas of the two parents: Populus deltoides and Populus trichocarpa

A Populus deltoides × Populus trichocarpa F1 pedigree was analyzed for quantitative trait loci (QTLs) affecting ectomycorrhizal development and for microarray characterization of gene networks involved in this symbiosis. A 300 genotype progeny set was evaluated for its ability to form ectomycorrhiza with the basidiomycete Laccaria bicolor. The percentage of mycorrhizal root tips was determined on the root systems of all 300 progeny and their two parents. QTL analysis identified four significant QTLs, one on the P. deltoides and three on the P. trichocarpa genetic maps. These QTLs were aligned to the P. trichocarpa genome and each contained several megabases and encompass numerous genes. NimbleGen whole-genome microarray, using cDNA from RNA extracts of ectomycorrhizal root tips from the parental genotypes P. trichocarpa and P. deltoides, was used to narrow the candidate gene list. Among the 1,543 differentially expressed genes (p value ≤ 0.05; ≥5.0-fold change in transcript level) having different transcript levels in mycorrhiza of the two parents, 41 transcripts were located in the QTL intervals: 20 in Myc_d1, 14 in Myc_t1, and seven in Myc_t2, while no significant differences among transcripts were found in Myc_t3. Among these 41 transcripts, 25 were overrepresented in P. deltoides relative to P. trichocarpa; 16 were overrepresented in P. trichocarpa. The transcript showing the highest overrepresentation in P. trichocarpa mycorrhiza libraries compared to P. deltoides mycorrhiza codes for an ethylene-sensitive EREBP-4 protein which may repress defense mechanisms in P. trichocarpa while the highest overrepresented transcripts in P. deltoides code for proteins/genes typically associated with pathogen resistance.

Characterization of Transposable Elements in the Ectomycorrhizal Fungus Laccaria bicolor

Background The publicly available Laccaria bicolor genome sequence has provided a considerable genomic resource allowing systematic identification of transposable elements (TEs) in this symbiotic ectomycorrhizal fungus. Using a TE-specific annotation pipeline we have characterized and analyzed TEs in the L. bicolor S238N-H82 genome. Methodology/Principal Findings TEs occupy 24% of the 60 Mb L. bicolor genome and represent 25,787 full-length and partial copy elements distributed within 171 families. The most abundant elements were the Copia-like. TEs are not randomly distributed across the genome, but are tightly nested or clustered. The majority of TEs exhibits signs of ancient transposition except some intact copies of terminal inverted repeats (TIRS), long terminal repeats (LTRs) and a large retrotransposon derivative (LARD) element. There were three main periods of TE expansion in L. bicolor: the first from 57 to 10 Mya, the second from 5 to 1 Mya and the most recent from 0.5 Mya ago until now. LTR retrotransposons are closely related to retrotransposons found in another basidiomycete, Coprinopsis cinerea. Conclusions This analysis 1) represents an initial characterization of TEs in the L. bicolor genome, 2) contributes to improve genome annotation and a greater understanding of the role TEs played in genome organization and evolution and 3) provides a valuable resource for future research on the genome evolution within the Laccaria genus.