The mycorrhizal symbiosis: research frontiers in genomics, ecology, and agricultural application

Mycorrhizal symbioses between plants and fungi are vital for the soil structure, nutrient cycling, plant diversity, and ecosystem sustainability. More than 250 000 plant species are associated with mycorrhizal fungi. Recent advances in genomics and related approaches have revolutionized our understanding of the biology and ecology of mycorrhizal associations. The genomes of 250+ mycorrhizal fungi have been released and hundreds of genes that play pivotal roles in regulating symbiosis development and metabolism have been characterized. rDNA metabarcoding and metatranscriptomics provide novel insights into the ecological cues driving mycorrhizal communities and functions expressed by these associations, linking genes to ecological traits such as nutrient acquisition and soil organic matter decomposition. Here, we review genomic studies that have revealed genes involved in nutrient uptake and symbiosis development, and discuss adaptations that are fundamental to the evolution of mycorrhizal lifestyles. We also evaluated the ecosystem services provided by mycorrhizal networks and discuss how mycorrhizal symbioses hold promise for sustainable agriculture and forestry by enhancing nutrient acquisition and stress tolerance. Overall, unraveling the intricate dynamics of mycorrhizal symbioses is paramount for promoting ecological sustainability and addressing current pressing environmental concerns. This review ends with major frontiers for further research.

Cynipid wasps systematically reprogram host metabolism and restructure cell walls in developing galls

Many insects have evolved the ability to manipulate plant growth to generate extraordinary structures called galls, in which insect larva can develop while being sheltered and feeding on the plant. In particular, cynipid (Hymenoptera: Cynipidae) wasps have evolved to form morphologically complex galls and generate an astonishing array of gall shapes, colors, and sizes. However, the biochemical basis underlying these remarkable cellular and developmental transformations remains poorly understood. A key determinant in plant cellular development is cell wall deposition that dictates the physical form and physiological function of newly developing cells, tissues, and organs. However, it is unclear to what degree cell walls are restructured to initiate and support the formation of new gall tissue. Here, we characterize the molecular alterations underlying gall development using a combination of metabolomic, histological, and biochemical techniques to elucidate how valley oak (Quercus lobata) leaf cells are reprogrammed to form galls. Strikingly, gall development involves an exceptionally coordinated spatial deposition of lignin and xylan to form de novo gall vasculature. Our results highlight how cynipid wasps can radically change the metabolite profile and restructure the cell wall to enable the formation of galls, providing insights into the mechanism of gall induction and the extent to which plants can be entirely reprogrammed to form unique structures and organs.

Centralized Interactive Phenomics Resource: an integrated online phenomics knowledgebase for health data users

OBJECTIVE

Development of clinical phenotypes from electronic health records (EHRs) can be resource intensive. Several phenotype libraries have been created to facilitate reuse of definitions. However, these platforms vary in target audience and utility. We describe the development of the Centralized Interactive Phenomics Resource (CIPHER) knowledgebase, a comprehensive public-facing phenotype library, which aims to facilitate clinical and health services research.

MATERIALS AND METHODS

The platform was designed to collect and catalog EHR-based computable phenotype algorithms from any healthcare system, scale metadata management, facilitate phenotype discovery, and allow for integration of tools and user workflows. Phenomics experts were engaged in the development and testing of the site.

RESULTS

The knowledgebase stores phenotype metadata using the CIPHER standard, and definitions are accessible through complex searching. Phenotypes are contributed to the knowledgebase via webform, allowing metadata validation. Data visualization tools linking to the knowledgebase enhance user interaction with content and accelerate phenotype development.

DISCUSSION

The CIPHER knowledgebase was developed in the largest healthcare system in the United States and piloted with external partners. The design of the CIPHER website supports a variety of front-end tools and features to facilitate phenotype development and reuse. Health data users are encouraged to contribute their algorithms to the knowledgebase for wider dissemination to the research community, and to use the platform as a springboard for phenotyping.

CONCLUSION

CIPHER is a public resource for all health data users available at https://phenomics.va.ornl.gov/ which facilitates phenotype reuse, development, and dissemination of phenotyping knowledge.

Thermal, water, and land cover factors led to contrasting urban and rural vegetation resilience to extreme hot months

With continuing global warming and urbanization, it is increasingly important to understand the resilience of urban vegetation to extreme high temperatures, but few studies have examined urban vegetation at large scale or both concurrent and delayed responses. In this study, we performed an urban-rural comparison using the Enhanced Vegetation Index and months that exceed the historical 90th percentile in mean temperature (referred to as "hot months") across 85 major cities in the contiguous United States. We found that hot months initially enhanced vegetation greenness but could cause a decline afterwards, especially for persistent (≥4 months) and intense (≥+2 °C) episodes in summer. The urban responses were more positive than rural in the western United States or in winter, but more negative during spring-autumn in the eastern United States. The east-west difference can be attributed to the higher optimal growth temperatures and lower water stress levels of the western urban vegetation than the rural. The urban responses also had smaller magnitudes than the rural responses, especially in deciduous forest biomes, and least in evergreen forest biomes. Within each biome, analysis at 1 km pixel level showed that impervious fraction and vegetation cover, local urban heat island intensity, and water stress were the key drivers of urban-rural differences. These findings advance our understanding of how prolonged exposure to warm extremes, particularly within urban environments, affects vegetation greenness and vitality. Urban planners and ecosystem managers should prioritize the long and intense events and the key drivers in fostering urban vegetation resilience to heat waves.

Development and Assessment of SNP Genotyping Arrays for Citrus and Its Close Relatives

Rapid advancements in technologies provide various tools to analyze fruit crop genomes to better understand genetic diversity and relationships and aid in breeding. Genome-wide single nucleotide polymorphism (SNP) genotyping arrays offer highly multiplexed assays at a relatively low cost per data point. We report the development and validation of 1.4M SNP Axiom® Citrus HD Genotyping Array (Citrus 15AX 1 and Citrus 15AX 2) and 58K SNP Axiom® Citrus Genotyping Arrays for Citrus and close relatives. SNPs represented were chosen from a citrus variant discovery panel consisting of 41 diverse whole-genome re-sequenced accessions of Citrus and close relatives, including eight progenitor citrus species. SNPs chosen mainly target putative genic regions of the genome and are accurately called in both Citrus and its closely related genera while providing good coverage of the nuclear and chloroplast genomes. Reproducibility of the arrays was nearly 100%, with a large majority of the SNPs classified as the most stringent class of markers, "PolyHighResolution" (PHR) polymorphisms. Concordance between SNP calls in sequence data and array data average 98%. Phylogenies generated with array data were similar to those with comparable sequence data and little affected by 3 to 5% genotyping error. Both arrays are publicly available.

Phycobilisome protein ApcG interacts with PSII and regulates energy transfer in Synechocystis

Photosynthetic organisms harvest light using pigment-protein complexes. In cyanobacteria, these are water-soluble antennae known as phycobilisomes (PBSs). The light absorbed by PBS is transferred to the photosystems in the thylakoid membrane to drive photosynthesis. The energy transfer between these complexes implies that protein-protein interactions allow the association of PBS with the photosystems. However, the specific proteins involved in the interaction of PBS with the photosystems are not fully characterized. Here, we show in Synechocystis sp. PCC 6803 that the recently discovered PBS linker protein ApcG (sll1873) interacts specifically with PSII through its N-terminal region. Growth of cyanobacteria is impaired in apcG deletion strains under light-limiting conditions. Furthermore, complementation of these strains using a phospho-mimicking version of ApcG causes reduced growth under normal growth conditions. Interestingly, the interaction of ApcG with PSII is affected when a phospho-mimicking version of ApcG is used, targeting the positively charged residues interacting with the thylakoid membrane, suggesting a regulatory role mediated by phosphorylation of ApcG. Low-temperature fluorescence measurements showed decreased PSI fluorescence in apcG deletion and complementation strains. The PSI fluorescence was the lowest in the phospho-mimicking complementation strain, while the pull-down experiment showed no interaction of ApcG with PSI under any tested condition. Our results highlight the importance of ApcG for selectively directing energy harvested by the PBS and imply that the phosphorylation status of ApcG plays a role in regulating energy transfer from PSII to PSI.

Declines in ice cover are accompanied by light limitation responses and community change in freshwater diatoms

The rediscovery of diatom blooms embedded within and beneath the Lake Erie ice cover (2007-2012) ignited interest in psychrophilic adaptations and winter limnology. Subsequent studies determined the vital role ice plays in winter diatom ecophysiology as diatoms partition to the underside of ice, thereby fixing their location within the photic zone. Yet, climate change has led to widespread ice decline across the Great Lakes, with Lake Erie presenting a nearly "ice-free" state in several recent winters. It has been hypothesized that the resultant turbid, isothermal water column induces light limitation amongst winter diatoms and thus serves as a competitive disadvantage. To investigate this hypothesis, we conducted a physiochemical and metatranscriptomic survey that spanned spatial, temporal, and climatic gradients of the winter Lake Erie water column (2019-2020). Our results suggest that ice-free conditions decreased planktonic diatom bloom magnitude and altered diatom community composition. Diatoms increased their expression of various photosynthetic genes and iron transporters, which suggests that the diatoms are attempting to increase their quantity of photosystems and light-harvesting components (a well-defined indicator of light limitation). We identified two gene families which serve to increase diatom fitness in the turbid ice-free water column: proton-pumping rhodopsins (a potential second means of light-driven energy acquisition) and fasciclins (a means to "raft" together to increase buoyancy and co-locate to the surface to optimize light acquisition). With large-scale climatic changes already underway, our observations provide insight into how diatoms respond to the dynamic ice conditions of today and shed light on how they will fare in a climatically altered tomorrow.

Disentangling the effects of sulfate and other seawater ions on microbial communities and greenhouse gas emissions in a coastal forested wetland

Seawater intrusion into freshwater wetlands causes changes in microbial communities and biogeochemistry, but the exact mechanisms driving these changes remain unclear. Here we use a manipulative laboratory microcosm experiment, combined with DNA sequencing and biogeochemical measurements, to tease apart the effects of sulfate from other seawater ions. We examined changes in microbial taxonomy and function as well as emissions of carbon dioxide, methane, and nitrous oxide in response to changes in ion concentrations. Greenhouse gas emissions and microbial richness and composition were altered by artificial seawater regardless of whether sulfate was present, whereas sulfate alone did not alter emissions or communities. Surprisingly, addition of sulfate alone did not lead to increases in the abundance of sulfate reducing bacteria or sulfur cycling genes. Similarly, genes involved in carbon, nitrogen, and phosphorus cycling responded more strongly to artificial seawater than to sulfate. These results suggest that other ions present in seawater, not sulfate, drive ecological and biogeochemical responses to seawater intrusion and may be drivers of increased methane emissions in soils that received artificial seawater addition. A better understanding of how the different components of salt water alter microbial community composition and function is necessary to forecast the consequences of coastal wetland salinization.

Altered Expression of Two Small Secreted Proteins (ssp4 and ssp6) Affects the Degradation of a Natural Lignocellulosic Substrate by Pleurotus ostreatus

Pleurotus ostreatus is a white-rot fungus that can degrade lignin in a preferential manner using a variety of extracellular enzymes, including manganese and versatile peroxidases (encoded by the vp1-3 and mnp1-6 genes, respectively). This fungus also secretes a family of structurally related small secreted proteins (SSPs) encoded by the ssp1-6 genes. Using RNA sequencing (RNA-seq), we determined that ssp4 and ssp6 are the predominant members of this gene family that were expressed by P. ostreatus during the first three weeks of growth on wheat straw. Downregulation of ssp4 in a strain harboring an ssp RNAi construct (KDssp1) was then confirmed, which, along with an increase in ssp6 transcript levels, coincided with reduced lignin degradation and the downregulation of vp2 and mnp1. In contrast, we observed an increase in the expression of genes related to pectin and side-chain hemicellulose degradation, which was accompanied by an increase in extracellular pectin-degrading capacity. Genome-wide comparisons between the KDssp1 and the wild-type strains demonstrated that ssp silencing conferred accumulated changes in gene expression at the advanced cultivation stages in an adaptive rather than an inductive mode of transcriptional response. Based on co-expression networking, crucial gene modules were identified and linked to the ssp knockdown genotype at different cultivation times. Based on these data, as well as previous studies, we propose that P. ostreatus SSPs have potential roles in modulating the lignocellulolytic and pectinolytic systems, as well as a variety of fundamental biological processes related to fungal growth and development.

Evolutionary Analyses of Gene Expression Divergence in Panicum hallii: Exploring Constitutive and Plastic Responses Using Reciprocal Transplants

The evolution of gene expression is thought to be an important mechanism of local adaptation and ecological speciation. Gene expression divergence occurs through the evolution of cis- polymorphisms and through more widespread effects driven by trans-regulatory factors. Here, we explore expression and sequence divergence in a large sample of Panicum hallii accessions encompassing the species range using a reciprocal transplantation experiment. We observed widespread genotype and transplant site drivers of expression divergence, with a limited number of genes exhibiting genotype-by-site interactions. We used a modified FST-QST outlier approach (QPC analysis) to detect local adaptation. We identified 514 genes with constitutive expression divergence above and beyond the levels expected under neutral processes. However, no plastic expression responses met our multiple testing correction as QPC outliers. Constitutive QPC outlier genes were involved in a number of developmental processes and responses to abiotic environments. Leveraging earlier expression quantitative trait loci results, we found a strong enrichment of expression divergence, including for QPC outliers, in genes previously identified with cis and cis-environment interactions but found no patterns related to trans-factors. Population genetic analyses detected elevated sequence divergence of promoters and coding sequence of constitutive expression outliers but little evidence for positive selection on these proteins. Our results are consistent with a hypothesis of cis-regulatory divergence as a primary driver of expression divergence in P. hallii.

Stepwise recombination suppression around the mating-type locus in an ascomycete fungus with self-fertile spores

Recombination is often suppressed at sex-determining loci in plants and animals, and at self-incompatibility or mating-type loci in plants and fungi. In fungal ascomycetes, recombination suppression around the mating-type locus is associated with pseudo-homothallism, i.e. the production of self-fertile dikaryotic sexual spores carrying the two opposite mating types. This has been well studied in two species complexes from different families of Sordariales: Podospora anserina and Neurospora tetrasperma. However, it is unclear whether this intriguing association holds in other species. We show here that Schizothecium tetrasporum, a fungus from a third family in the order Sordariales, also produces mostly self-fertile dikaryotic spores carrying the two opposite mating types. This was due to a high frequency of second meiotic division segregation at the mating-type locus, indicating the occurrence of a single and systematic crossing-over event between the mating-type locus and the centromere, as in P. anserina. The mating-type locus has the typical Sordariales organization, plus a MAT1-1-1 pseudogene in the MAT1-2 haplotype. High-quality genome assemblies of opposite mating types and segregation analyses revealed a suppression of recombination in a region of 1.47 Mb around the mating-type locus. We detected three evolutionary strata, indicating a stepwise extension of recombination suppression. The three strata displayed no rearrangement or transposable element accumulation but gene losses and gene disruptions were present, and precisely at the strata margins. Our findings indicate a convergent evolution of self-fertile dikaryotic sexual spores across multiple ascomycete fungi. The particular pattern of meiotic segregation at the mating-type locus was associated with recombination suppression around this locus, that had extended stepwise. This association between pseudo-homothallism and recombination suppression across lineages and the presence of gene disruption at the strata limits are consistent with a recently proposed mechanism of sheltering deleterious alleles to explain stepwise recombination suppression.

Remarkably coherent population structure for a dominant Antarctic Chlorobium species

BACKGROUND

In Antarctica, summer sunlight enables phototrophic microorganisms to drive primary production, thereby "feeding" ecosystems to enable their persistence through the long, dark winter months. In Ace Lake, a stratified marine-derived system in the Vestfold Hills of East Antarctica, a Chlorobium species of green sulphur bacteria (GSB) is the dominant phototroph, although its seasonal abundance changes more than 100-fold. Here, we analysed 413 Gb of Antarctic metagenome data including 59 Chlorobium metagenome-assembled genomes (MAGs) from Ace Lake and nearby stratified marine basins to determine how genome variation and population structure across a 7-year period impacted ecosystem function.

RESULTS

A single species, Candidatus Chlorobium antarcticum (most similar to Chlorobium phaeovibrioides DSM265) prevails in all three aquatic systems and harbours very little genomic variation (≥ 99% average nucleotide identity). A notable feature of variation that did exist related to the genomic capacity to biosynthesize cobalamin. The abundance of phylotypes with this capacity changed seasonally ~ 2-fold, consistent with the population balancing the value of a bolstered photosynthetic capacity in summer against an energetic cost in winter. The very high GSB concentration (> 10⁸ cells ml^-1 in Ace Lake) and seasonal cycle of cell lysis likely make Ca. Chlorobium antarcticum a major provider of cobalamin to the food web. Analysis of Ca. Chlorobium antarcticum viruses revealed the species to be infected by generalist (rather than specialist) viruses with a broad host range (e.g., infecting Gammaproteobacteria) that were present in diverse Antarctic lakes. The marked seasonal decrease in Ca. Chlorobium antarcticum abundance may restrict specialist viruses from establishing effective lifecycles, whereas generalist viruses may augment their proliferation using other hosts.

CONCLUSION

The factors shaping Antarctic microbial communities are gradually being defined. In addition to the cold, the annual variation in sunlight hours dictates which phototrophic species can grow and the extent to which they contribute to ecosystem processes. The Chlorobium population studied was inferred to provide cobalamin, in addition to carbon, nitrogen, hydrogen, and sulphur cycling, as critical ecosystem services. The specific Antarctic environmental factors and major ecosystem benefits afforded by this GSB likely explain why such a coherent population structure has developed in this Chlorobium species. Video abstract.

Three founding ancestral genomes involved in the origin of sugarcane

BACKGROUND AND AIMS

Modern sugarcane cultivars (Saccharum spp.) are high polyploids, aneuploids (2n = ~12x = ~120) derived from interspecific hybridizations between the domesticated sweet species Saccharum officinarum and the wild species S. spontaneum.

METHODS

To analyse the architecture and origin of such a complex genome, we analysed the sequences of all 12 hom(oe)ologous haplotypes (BAC clones) from two distinct genomic regions of a typical modern cultivar, as well as the corresponding sequence in Miscanthus sinense and Sorghum bicolor, and monitored their distribution among representatives of the Saccharum genus.

KEY RESULTS

The diversity observed among haplotypes suggested the existence of three founding genomes (A, B, C) in modern cultivars, which diverged between 0.8 and 1.3 Mya. Two genomes (A, B) were contributed by S. officinarum; these were also found in its wild presumed ancestor S. robustum, and one genome (C) was contributed by S. spontaneum. These results suggest that S. officinarum and S. robustum are derived from interspecific hybridization between two unknown ancestors (A and B genomes). The A genome contributed most haplotypes (nine or ten) while the B and C genomes contributed one or two haplotypes in the regions analysed of this typical modern cultivar. Interspecific hybridizations likely involved accessions or gametes with distinct ploidy levels and/or were followed by a series of backcrosses with the A genome. The three founding genomes were found in all S. barberi, S. sinense and modern cultivars analysed. None of the analysed accessions contained only the A genome or the B genome, suggesting that representatives of these founding genomes remain to be discovered.

CONCLUSIONS

This evolutionary model, which combines interspecificity and high polyploidy, can explain the variable chromosome pairing affinity observed in Saccharum. It represents a major revision of the understanding of Saccharum diversity.

Phylogenomics of Endogonaceae and evolution of mycorrhizas within Mucoromycota

Endogonales (Mucoromycotina), composed of Endogonaceae and Densosporaceae, is the only known non-Dikarya order with ectomycorrhizal members. They also form mycorrhizal-like association with some nonspermatophyte plants. It has been recently proposed that Endogonales were among the earliest mycorrhizal partners with land plants. It remains unknown whether Endogonales possess genomes with mycorrhizal-lifestyle signatures and whether Endogonales originated around the same time as land plants did. We sampled sporocarp tissue from four Endogonaceae collections and performed shotgun genome sequencing. After binning the metagenome data, we assembled and annotated the Endogonaceae genomes. We performed comparative analysis on plant-cell-wall-degrading enzymes (PCWDEs) and small secreted proteins (SSPs). We inferred phylogenetic placement of Endogonaceae and estimated the ages of Endogonaceae and Endogonales with expanded taxon sampling. Endogonaceae have large genomes with high repeat content, low diversity of PCWDEs, but without elevated SSP/secretome ratios. Dating analysis estimated that Endogonaceae originated in the Permian-Triassic boundary and Endogonales originated in the mid-late Silurian. Mycoplasma-related endobacterium sequences were identified in three Endogonaceae genomes. Endogonaceae genomes possess typical signatures of mycorrhizal lifestyle. The early origin of Endogonales suggests that the mycorrhizal association between Endogonales and plants might have played an important role during the colonization of land by plants.

Population genomics of picophytoplankton unveils novel chromosome hypervariability

Tiny photosynthetic microorganisms that form the picoplankton (between 0.3 and 3 μm in diameter) are at the base of the food web in many marine ecosystems, and their adaptability to environmental change hinges on standing genetic variation. Although the genomic and phenotypic diversity of the bacterial component of the oceans has been intensively studied, little is known about the genomic and phenotypic diversity within each of the diverse eukaryotic species present. We report the level of genomic diversity in a natural population of Ostreococcus tauri (Chlorophyta, Mamiellophyceae), the smallest photosynthetic eukaryote. Contrary to the expectations of clonal evolution or cryptic species, the spectrum of genomic polymorphism observed suggests a large panmictic population (an effective population size of 1.2 × 10⁷) with pervasive evidence of sexual reproduction. De novo assemblies of low-coverage chromosomes reveal two large candidate mating-type loci with suppressed recombination, whose origin may pre-date the speciation events in the class Mamiellophyceae. This high genetic diversity is associated with large phenotypic differences between strains. Strikingly, resistance of isolates to large double-stranded DNA viruses, which abound in their natural environment, is positively correlated with the size of a single hypervariable chromosome, which contains 44 to 156 kb of strain-specific sequences. Our findings highlight the role of viruses in shaping genome diversity in marine picoeukaryotes.