A pan-grass transcriptome reveals patterns of cellular divergence in crops.

Different plant species within the grasses were parallel targets of domestication, giving rise to crops with distinct evolutionary histories and traits1. Key traits that distinguish these species are mediated by specialized cell types2. Here we compare the transcriptomes of root cells in three grass species-Zea mays, Sorghum bicolor and Setaria viridis. We show that single-cell and single-nucleus RNA sequencing provide complementary readouts of cell identity in dicots and monocots, warranting a combined analysis. Cell types were mapped across species to identify robust, orthologous marker genes. The comparative cellular analysis shows that the transcriptomes of some cell types diverged more rapidly than those of others-driven, in part, by recruitment of gene modules from other cell types. The data also show that a recent whole-genome duplication provides a rich source of new, highly localized gene expression domains that favour fast-evolving cell types. Together, the cell-by-cell comparative analysis shows how fine-scale cellular profiling can extract conserved modules from a pan transcriptome and provide insight on the evolution of cells that mediate key functions in crops.

The strength and pattern of natural selection on gene expression in rice.

Levels of gene expression underpin organismal phenotypes1,2, but the nature of selection that acts on gene expression and its role in adaptive evolution remain unknown1,2. Here we assayed gene expression in rice (Oryza sativa)3, and used phenotypic selection analysis to estimate the type and strength of selection on the levels of more than 15,000 transcripts4,5. Variation in most transcripts appears (nearly) neutral or under very weak stabilizing selection in wet paddy conditions (with median standardized selection differentials near zero), but selection is stronger under drought conditions. Overall, more transcripts are conditionally neutral (2.83%) than are antagonistically pleiotropic6 (0.04%), and transcripts that display lower levels of expression and stochastic noise7-9 and higher levels of plasticity9 are under stronger selection. Selection strength was further weakly negatively associated with levels of cis-regulation and network connectivity9. Our multivariate analysis suggests that selection acts on the expression of photosynthesis genes4,5, but that the efficacy of selection is genetically constrained under drought conditions10. Drought selected for earlier flowering11,12 and a higher expression of OsMADS18 (Os07g0605200), which encodes a MADS-box transcription factor and is a known regulator of early flowering13-marking this gene as a drought-escape gene11,12. The ability to estimate selection strengths provides insights into how selection can shape molecular traits at the core of gene action.

Evolutionary systems biology reveals patterns of rice adaptation to drought-prone agro-ecosystems.

Rice (Oryza sativa) was domesticated around 10,000 years ago and has developed into a staple for half of humanity. The crop evolved and is currently grown in stably wet and intermittently dry agro-ecosystems, but patterns of adaptation to differences in water availability remain poorly understood. While previous field studies have evaluated plant developmental adaptations to water deficit, adaptive variation in functional and hydraulic components, particularly in relation to gene expression, has received less attention. Here, we take an evolutionary systems biology approach to characterize adaptive drought resistance traits across roots and shoots. We find that rice harbors heritable variation in molecular, physiological, and morphological traits that is linked to higher fitness under drought. We identify modules of co-expressed genes that are associated with adaptive drought avoidance and tolerance mechanisms. These expression modules showed evidence of polygenic adaptation in rice subgroups harboring accessions that evolved in drought-prone agro-ecosystems. Fitness-linked expression patterns allowed us to identify the drought-adaptive nature of optimizing photosynthesis and interactions with arbuscular mycorrhizal fungi. Taken together, our study provides an unprecedented, integrative view of rice adaptation to water-limited field conditions.

Genome editing retraces the evolution of toxin resistance in the monarch butterfly.

Identifying the genetic mechanisms of adaptation requires the elucidation of links between the evolution of DNA sequence, phenotype, and fitness1. Convergent evolution can be used as a guide to identify candidate mutations that underlie adaptive traits2-4, and new genome editing technology is facilitating functional validation of these mutations in whole organisms1,5. We combined these approaches to study a classic case of convergence in insects from six orders, including the monarch butterfly (Danaus plexippus), that have independently evolved to colonize plants that produce cardiac glycoside toxins6-11. Many of these insects evolved parallel amino acid substitutions in the α-subunit (ATPα) of the sodium pump (Na+/K+-ATPase)7-11, the physiological target of cardiac glycosides12. Here we describe mutational paths involving three repeatedly changing amino acid sites (111, 119 and 122) in ATPα that are associated with cardiac glycoside specialization13,14. We then performed CRISPR-Cas9 base editing on the native Atpα gene in Drosophila melanogaster flies and retraced the mutational path taken across the monarch lineage11,15. We show in vivo, in vitro and in silico that the path conferred resistance and target-site insensitivity to cardiac glycosides16, culminating in triple mutant 'monarch flies' that were as insensitive to cardiac glycosides as monarch butterflies. 'Monarch flies' retained small amounts of cardiac glycosides through metamorphosis, a trait that has been optimized in monarch butterflies to deter predators17-19. The order in which the substitutions evolved was explained by amelioration of antagonistic pleiotropy through epistasis13,14,20-22. Our study illuminates how the monarch butterfly evolved resistance to a class of plant toxins, eventually becoming unpalatable, and changing the nature of species interactions within ecological communities2,6-11,15,17-19.

Toxic to the touch: The makings of lethal mantles in pitohui birds and poison dart frogs.

How do chemically defended animals resist their own toxins? This intriguing question on the concept of autotoxicity is at the heart of how species interactions evolve. In this issue of Molecular Ecology (Molecular Ecology, 2024, 33), Bodawatta and colleagues report on how Papua New Guinean birds coopted deadly neurotoxins to create lethal mantles that protect against predators and parasites. Combining chemical screening of the plumage of a diverse collection of passerine birds with genome sequencing, the researchers unlocked a deeper understanding of how some birds sequester deadly batrachotoxin (BTX) from their food without poisoning themselves. They identified that birds impervious to BTX bear amino acid substitutions in the toxin-binding site of the voltage-gated sodium channel Nav1.4, whose function is essential for proper contraction and relaxation of vertebrate muscles. Comparative genetic and molecular docking analyses show that several of the substitutions associated with insensitivity to BTX may have become prevalent among toxic birds through positive selection. Intriguingly, poison dart frogs that also co-opted BTX in their lethal mantles were found to harbour similar toxin insensitivity substitutions in their Nav1.4 channels. Taken together, this sets up a powerful model system for studying the mechanisms behind convergent molecular evolution and how it may drive biological diversity.

Rapid evolutionary changes in gene expression in response to climate fluctuations.

There is now abundant evidence of rapid evolution in natural populations, but the genetic mechanisms of these changes remain unclear. One possible route to rapid evolution is through changes in the expression of genes that influence traits under selection. We examined contemporary evolutionary gene expression changes in plant populations responding to environmental fluctuations. We compared genome-wide gene expression, using RNA-seq, in two populations of Brassica rapa collected over four time points between 1997 and 2014, during which precipitation in southern California fluctuated dramatically and phenotypic and genotypic changes occurred. By combining transcriptome profiling with the resurrection approach, we directly examined evolutionary changes in gene expression over time. For both populations, we found a substantial number of differentially expressed genes between generations, indicating rapid evolution in the expression of many genes. Using existing gene annotations, we found that many changes occurred in genes involved in regulating stress responses and flowering time. These appeared related to the fluctuations in precipitation and were potentially adaptive. However, the evolutionary changes in gene expression differed across generations within and between populations, indicating largely independent evolutionary trajectories across populations and over time. Our study provides strong evidence for rapid evolution in gene expression, and indicates that changes in gene expression can be one mechanism of rapid evolutionary responses to selection episodes. This study also illustrates that combining resurrection studies with transcriptomics is a powerful approach for investigating evolutionary changes at the gene regulatory level, and will provide new insights into the genetic basis of contemporary evolution.

Evolutionary and ecological functional genomics, from lab to the wild.

Plant phenotypes are the result of both genetic and environmental forces that act to modulate trait expression. Over the last few years, numerous approaches in functional genomics and systems biology have led to a greater understanding of plant phenotypic variation and plant responses to the environment. These approaches, and the questions that they can address, have been loosely termed evolutionary and ecological functional genomics (EEFG), and have been providing key insights on how plants adapt and evolve. In particular, by bringing these studies from the laboratory to the field, EEFG studies allow us to gain greater knowledge of how plants function in their natural contexts.

Virus Infection of Plants Alters Pollinator Preference: A Payback for Susceptible Hosts?

Plant volatiles play important roles in attraction of certain pollinators and in host location by herbivorous insects. Virus infection induces changes in plant volatile emission profiles, and this can make plants more attractive to insect herbivores, such as aphids, that act as viral vectors. However, it is unknown if virus-induced alterations in volatile production affect plant-pollinator interactions. We found that volatiles emitted by cucumber mosaic virus (CMV)-infected tomato (Solanum lycopersicum) and Arabidopsis thaliana plants altered the foraging behaviour of bumblebees (Bombus terrestris). Virus-induced quantitative and qualitative changes in blends of volatile organic compounds emitted by tomato plants were identified by gas chromatography-coupled mass spectrometry. Experiments with a CMV mutant unable to express the 2b RNA silencing suppressor protein and with Arabidopsis silencing mutants implicate microRNAs in regulating emission of pollinator-perceivable volatiles. In tomato, CMV infection made plants emit volatiles attractive to bumblebees. Bumblebees pollinate tomato by 'buzzing' (sonicating) the flowers, which releases pollen and enhances self-fertilization and seed production as well as pollen export. Without buzz-pollination, CMV infection decreased seed yield, but when flowers of mock-inoculated and CMV-infected plants were buzz-pollinated, the increased seed yield for CMV-infected plants was similar to that for mock-inoculated plants. Increased pollinator preference can potentially increase plant reproductive success in two ways: i) as female parents, by increasing the probability that ovules are fertilized; ii) as male parents, by increasing pollen export. Mathematical modeling suggested that over a wide range of conditions in the wild, these increases to the number of offspring of infected susceptible plants resulting from increased pollinator preference could outweigh underlying strong selection pressures favoring pathogen resistance, allowing genes for disease susceptibility to persist in plant populations. We speculate that enhanced pollinator service for infected individuals in wild plant populations might provide mutual benefits to the virus and its susceptible hosts.

A genomic perspective on the generation and maintenance of genetic diversity in herbivorous insects.

Understanding the processes that generate and maintain genetic variation within populations is a central goal in evolutionary biology. Theory predicts that some of this variation is maintained as a consequence of adapting to variable habitats. Studies in herbivorous insects have played a key role in confirming this prediction. Here, we highlight theoretical and conceptual models for the maintenance of genetic diversity in herbivorous insects, empirical genomic studies testing these models, and pressing questions within the realm of evolutionary and functional genomic studies. To address key gaps, we propose an integrative approach combining population genomic scans for adaptation, genome-wide characterization of targets of selection through experimental manipulations, mapping the genetic architecture of traits influencing fitness, and functional studies. We also stress the importance of studying the maintenance of genetic variation across biological scales-from variation within populations to divergence among populations-to form a comprehensive view of adaptation in herbivorous insects.

Cucumber mosaic virus and its 2b protein alter emission of host volatile organic compounds but not aphid vector settling in tobacco.

BACKGROUND: Aphids, including the generalist herbivore Myzus persicae, transmit cucumber mosaic virus (CMV). CMV (strain Fny) infection affects M. persicae feeding behavior and performance on tobacco (Nicotiana tabacum), Arabidopsis thaliana and cucurbits in varying ways. In Arabidopsis and cucurbits, CMV decreases host quality and inhibits prolonged feeding by aphids, which may enhance virus transmission rates. CMV-infected cucurbits also emit deceptive, aphid-attracting volatiles, which may favor virus acquisition. In contrast, aphids on CMV-infected tobacco (cv. Xanthi) exhibit increased survival and reproduction. This may not increase transmission but might increase virus and vector persistence within plant communities. The CMV 2b counter-defense protein diminishes resistance to aphid infestation in CMV-infected tobacco plants. We hypothesised that in tobacco CMV and its 2b protein might also alter the emission of volatile organic compounds that would influence aphid behavior. RESULTS: Analysis of headspace volatiles emitted from tobacco plants showed that CMV infection both increased the total quantity and altered the blend produced. Furthermore, experiments with a CMV 2b gene deletion mutant (CMV∆2b) showed that the 2b counter-defense protein influences volatile emission. Free choice bioassays were conducted where wingless M. persicae could choose to settle on infected or mock-inoculated plants under a normal day/night regime or in continual darkness. Settling was recorded at 15 min, 1 h and 24 h post-release. Statistical analysis indicated that aphids showed no marked preference to settle on mock-inoculated versus infected plants, except for a marginally greater settlement of aphids on mock-inoculated over CMV-infected plants under normal illumination. CONCLUSIONS: CMV infection of tobacco plants induced quantitative and qualitative changes in host volatile emission and these changes depended in part on the activity of the 2b counter-defense protein. However, CMV-induced alterations in tobacco plant volatile emission did not have marked effects on the settling of aphids on infected versus mock-inoculated plants even though CMV-infected plants are higher quality hosts for M. persicae.

Pathogen-triggered ethylene signaling mediates systemic-induced susceptibility to herbivory in Arabidopsis.

Multicellular eukaryotic organisms are attacked by numerous parasites from diverse phyla, often simultaneously or sequentially. An outstanding question in these interactions is how hosts integrate signals induced by the attack of different parasites. We used a model system comprised of the plant host Arabidopsis thaliana, the hemibiotrophic bacterial phytopathogen Pseudomonas syringae, and herbivorous larvae of the moth Trichoplusia ni (cabbage looper) to characterize mechanisms involved in systemic-induced susceptibility (SIS) to T. ni herbivory caused by prior infection by virulent P. syringae. We uncovered a complex multilayered induction mechanism for SIS to herbivory. In this mechanism, antiherbivore defenses that depend on signaling via (1) the jasmonic acid-isoleucine conjugate (JA-Ile) and (2) other octadecanoids are suppressed by microbe-associated molecular pattern-triggered salicylic acid (SA) signaling and infection-triggered ethylene signaling, respectively. SIS to herbivory is, in turn, counteracted by a combination of the bacterial JA-Ile mimic coronatine and type III virulence-associated effectors. Our results show that SIS to herbivory involves more than antagonistic signaling between SA and JA-Ile and provide insight into the unexpectedly complex mechanisms behind a seemingly simple trade-off in plant defense against multiple enemies.

The evolution of ethylene signaling in plant chemical ecology.

Ethylene is a key hormone in plant development, mediating plant responses to abiotic environmental stress, and interactions with attackers and mutualists. Here, we provide a synthesis of the role of ethylene in the context of plant ecology and evolution, and a prospectus for future research in this area. We focus on the regulatory function of ethylene in multi-organismal interactions. In general, plant interactions with different types of organisms lead to reduced or enhanced levels of ethylene. This in turn affects not only the plant's response to the interacting organism at hand, but also to other organisms in the community. These community-level effects become observable as enhanced or diminished relationships with future commensals, and systemic resistance or susceptibility to secondary attackers. Ongoing comparative genomic and phenotypic analyses continue to shed light on these interactions. These studies have revealed that plants and interacting organisms from separate kingdoms of life have independently evolved the ability to produce, perceive, and respond to ethylene. This signature of convergent evolution of ethylene signaling at the phenotypic level highlights the central role ethylene metabolism and signaling plays in plant interactions with microbes and animals.

Systems genomics of salinity stress response in rice.

Populations can adapt to stressful environments through changes in gene expression. However, the role of gene regulation in mediating stress response and adaptation remains largely unexplored. Here, we use an integrative field dataset obtained from 780 plants of Oryza sativa ssp. indica (rice) grown in a field experiment under normal or moderate salt stress conditions to examine selection and evolution of gene expression variation under salinity stress conditions. We find that salinity stress induces increased selective pressure on gene expression. Further, we show that trans-eQTLs rather than cis-eQTLs are primarily associated with rice's gene expression under salinity stress, potentially via a few master-regulators. Importantly, and contrary to the expectations, we find that cis-trans reinforcement is more common than cis-trans compensation which may be reflective of rice diversification subsequent to domestication. We further identify genetic fixation as the likely mechanism underlying this compensation/reinforcement. Additionally, we show that cis- and trans-eQTLs are under different selection regimes, giving us insights into the evolutionary dynamics of gene expression variation. By examining genomic, transcriptomic, and phenotypic variation across a rice population, we gain insights into the molecular and genetic landscape underlying adaptive salinity stress responses, which is relevant for other crops and other stresses.

Evolution of chemosensory and detoxification gene families across herbivorous Drosophilidae.

Herbivorous insects are exceptionally diverse, accounting for a quarter of all known eukaryotic species, but the genomic basis of adaptations that enabled this dietary transition remains poorly understood. Many studies have suggested that expansions and contractions of chemosensory and detoxification gene families-genes directly mediating interactions with plant chemical defenses-underlie successful plant colonization. However, this hypothesis has been challenging to test because the origins of herbivory in many insect lineages are ancient (>150 million years ago (mya)), obscuring genomic evolutionary patterns. Here, we characterized chemosensory and detoxification gene family evolution across Scaptomyza, a genus nested within Drosophila that includes a recently derived (<15 mya) herbivore lineage of mustard (Brassicales) specialists and carnation (Caryophyllaceae) specialists, and several nonherbivorous species. Comparative genomic analyses revealed that herbivorous Scaptomyza has among the smallest chemosensory and detoxification gene repertoires across 12 drosophilid species surveyed. Rates of gene turnover averaged across the herbivore clade were significantly higher than background rates in over half of the surveyed gene families. However, gene turnover was more limited along the ancestral herbivore branch, with only gustatory receptors and odorant-binding proteins experiencing strong losses. The genes most significantly impacted by gene loss, duplication, or changes in selective constraint were those involved in detecting compounds associated with feeding on living plants (bitter or electrophilic phytotoxins) or their ancestral diet (fermenting plant volatiles). These results provide insight into the molecular and evolutionary mechanisms of plant-feeding adaptations and highlight gene candidates that have also been linked to other dietary transitions in Drosophila.

Evolution of chemosensory and detoxification gene families across herbivorous Drosophilidae.

Herbivorous insects are exceptionally diverse, accounting for a quarter of all known eukaryotic species, but the genetic basis of adaptations that enabled this dietary transition remains poorly understood. Many studies have suggested that expansions and contractions of chemosensory and detoxification gene families - genes directly mediating interactions with plant chemical defenses - underlie successful plant colonization. However, this hypothesis has been challenging to test because the origins of herbivory in many lineages are ancient (>150 million years ago [mya]), obscuring genomic evolutionary patterns. Here, we characterized chemosensory and detoxification gene family evolution across Scaptomyza, a genus nested within Drosophila that includes a recently derived (<15 mya) herbivore lineage of mustard (Brassicales) specialists and carnation (Caryophyllaceae) specialists, and several non-herbivorous species. Comparative genomic analyses revealed that herbivorous Scaptomyza have among the smallest chemosensory and detoxification gene repertoires across 12 drosophilid species surveyed. Rates of gene turnover averaged across the herbivore clade were significantly higher than background rates in over half of the surveyed gene families. However, gene turnover was more limited along the ancestral herbivore branch, with only gustatory receptors and odorant binding proteins experiencing strong losses. The genes most significantly impacted by gene loss, duplication, or changes in selective constraint were those involved in detecting compounds associated with feeding on plants (bitter or electrophilic phytotoxins) or their ancestral diet (yeast and fruit volatiles). These results provide insight into the molecular and evolutionary mechanisms of plant-feeding adaptations and highlight strong gene candidates that have also been linked to other dietary transitions in Drosophila .

The influence of genetic architecture on responses to selection under drought in rice.

Accurately predicting responses to selection is a major goal in biology and important for successful crop breeding in changing environments. However, evolutionary responses to selection can be constrained by such factors as genetic and cross-environment correlations, linkage, and pleiotropy, and our understanding of the extent and impact of such constraints is still developing. Here, we conducted a field experiment to investigate potential constraints to selection for drought resistance in rice (Oryza sativa) using phenotypic selection analysis and quantitative genetics. We found that traits related to drought response were heritable, and some were under selection, including selection for earlier flowering, which could allow drought escape. However, patterns of selection generally were not opposite under wet and dry conditions, and we did not find individual or closely linked genes that influenced multiple traits, indicating a lack of evidence that antagonistic pleiotropy, linkage, or cross-environment correlations would constrain selection for drought resistance. In most cases, genetic correlations had little influence on responses to selection, with direct and indirect selection largely congruent. The exception to this was seed mass under drought, which was predicted to evolve in the opposite direction of direct selection due to correlations. Because of this indirect effect on selection on seed mass, selection for drought resistance was not accompanied by a decrease in seed mass, and yield increased with fecundity. Furthermore, breeding lines with high fitness and yield under drought also had high fitness and yield under wet conditions, indicating that there was no evidence for a yield penalty on drought resistance. We found multiple genes in which expression influenced both water use efficiency (WUE) and days to first flowering, supporting a genetic basis for the trade-off between drought escape and avoidance strategies. Together, these results can provide helpful guidance for understanding and managing evolutionary constraints and breeding stress-resistant crops.

Mining the plant-herbivore interface with a leafmining Drosophila of Arabidopsis.

Experimental infections of Arabidopsis thaliana (Arabidopsis) with genomically characterized plant pathogens such as Pseudomonas syringae have facilitated the dissection of canonical eukaryotic defence pathways and parasite virulence factors. Plants are also attacked by herbivorous insects, and the development of an ecologically relevant genetic model herbivore that feeds on Arabidopsis will enable the parallel dissection of host defence and reciprocal resistance pathways such as those involved in xenobiotic metabolism. An ideal candidate is Scaptomyza flava, a drosophilid fly whose leafmining larvae are true herbivores that can be found in nature feeding on Arabidopsis and other crucifers. Here, we describe the life cycle of S. flava on Arabidopsis and use multiple approaches to characterize the response of Arabidopsis to S. flava attack. Oviposition choice tests and growth performance assays on different Arabidopsis ecotypes, defence-related mutants, and hormone and chitin-treated plants revealed significant differences in host preference and variation in larval performance across Arabidopsis accessions. The jasmonate and glucosinolate pathways in Arabidopsis are important in mediating quantitative resistance against S. flava, and priming with jasmonate or chitin resulted in increased resistance. Expression of xenobiotic detoxification genes was reduced in S. flava larvae reared on Arabidopsis jasmonate signalling mutants and increased in plants pretreated with chitin. These results and future research directions are discussed in the context of developing a genetic model system to analyse insect-plant interactions.

Convergent evolution of cardiac-glycoside resistance in predators and parasites of milkweed herbivores.

The community of plant-feeding insects (herbivores) that specialize on milkweeds (Apocynaceae) form a remarkable example of convergent evolution across levels of biological organization1. In response to toxic cardiac glycosides produced by these plants, the monarch butterfly (Danaus plexippus) and other specialist herbivores have evolved parallel substitutions in the alpha subunit (ATPA) of the Na+/K+-ATPase. These substitutions render the pump insensitive to cardiac glycosides2,3, allowing the monarch and other specialists, from aphids to beetles, to sequester cardiac glycosides, which in turn provide defense against attacks by enemies from the third trophic level4. The evolution of 'target-site-insensitivity' substitutions in these herbivores poses a fundamental biological question: have predators and parasitoids that feed on cardiac-glycoside-sequestering insects also evolved Na+/K+-ATPases that are similarly insensitive to cardiac glycosides (as predicted by Whiteman and Mooney)5? In other words, can plant toxins cause evolutionary cascades that reach the third trophic level? Here we show that at least four enemies of the monarch and other milkweed herbivores have indeed evolved amino-acid substitutions associated with target-site insensitivity to cardiac glycosides. These attackers represent four major animal clades, implicating cardiac glycosides as keystone molecules6 and establishing ATPalpha, which encodes ATPA, as a keystone gene with effects that reverberate within ecological communities7.