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.

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.

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.

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.

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.

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.

Infection of Arabidopsis by cucumber mosaic virus triggers jasmonate-dependent resistance to aphids that relies partly on the pattern-triggered immunity factor BAK1.

Many aphid-vectored viruses are transmitted nonpersistently via transient attachment of virus particles to aphid mouthparts and are most effectively acquired or transmitted during brief stylet punctures of epidermal cells. In Arabidopsis thaliana, the aphid-transmitted virus cucumber mosaic virus (CMV) induces feeding deterrence against the polyphagous aphid Myzus persicae. This form of resistance inhibits prolonged phloem feeding but promotes virus acquisition by aphids because it encourages probing of plant epidermal cells. When aphids are confined on CMV-infected plants, feeding deterrence reduces their growth and reproduction. We found that CMV-induced inhibition of growth as well as CMV-induced inhibition of reproduction of M. persicae are dependent upon jasmonate-mediated signalling. BRASSINOSTEROID INSENSITIVE1-ASSOCIATED KINASE1 (BAK1) is a co-receptor enabling detection of microbe-associated molecular patterns and induction of pattern-triggered immunity (PTI). In plants carrying the mutant bak1-5 allele, CMV induced inhibition of M. persicae reproduction but not inhibition of aphid growth. We conclude that in wildtype plants CMV induces two mechanisms that diminish performance of M. persicae: a jasmonate-dependent and PTI-dependent mechanism that inhibits aphid growth, and a jasmonate-dependent, PTI-independent mechanism that inhibits reproduction. The growth of two crucifer specialist aphids, Lipaphis erysimi and Brevicoryne brassicae, was not affected when confined on CMV-infected A. thaliana. However, B. brassicae reproduction was inhibited on CMV-infected plants. This suggests that in A. thaliana CMV-induced resistance to aphids, which is thought to incentivize virus vectoring, has greater effects on polyphagous than on crucifer specialist aphids.

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.

Characterization of the Leaf Microbiome from Whole-Genome Sequencing Data of the 3000 Rice Genomes Project.

BACKGROUND: The crop microbial communities are shaped by interactions between the host, microbes and the environment, however, their relative contribution is beginning to be understood. Here, we explore these interactions in the leaf bacterial community across 3024 rice accessions. FINDINGS: By using unmapped DNA sequencing reads as microbial reads, we characterized the structure of the rice bacterial microbiome. We identified central bacteria taxa that emerge as microbial "hubs" and may have an influence on the network of host-microbe interactions. We found regions in the rice genome that might control the assembly of these microbial hubs. To our knowledge this is one of the first studies that uses raw data from plant genome sequencing projects to characterize the leaf bacterial communities. CONCLUSION: We showed, that the structure of the rice leaf microbiome is modulated by multiple interactions among host, microbes, and environment. Our data provide insight into the factors influencing microbial assemblage in the rice leaf and also opens the door for future initiatives to modulate rice consortia for crop improvement efforts.

Genomic history and ecology of the geographic spread of rice.

Rice (Oryza sativa) is one of the world's most important food crops, and is comprised largely of japonica and indica subspecies. Here, we reconstruct the history of rice dispersal in Asia using whole-genome sequences of more than 1,400 landraces, coupled with geographic, environmental, archaeobotanical and paleoclimate data. Originating around 9,000 yr ago in the Yangtze Valley, rice diversified into temperate and tropical japonica rice during a global cooling event about 4,200 yr ago. Soon after, tropical japonica rice reached Southeast Asia, where it rapidly diversified, starting about 2,500 yr BP. The history of indica rice dispersal appears more complicated, moving into China around 2,000 yr BP. We also identify extrinsic factors that influence genome diversity, with temperature being a leading abiotic factor. Reconstructing the dispersal history of rice and its climatic correlates may help identify genetic adaptations associated with the spread of a key domesticated species.

Nanopore sequencing-based genome assembly and evolutionary genomics of circum-basmati rice.

BACKGROUND: The circum-basmati group of cultivated Asian rice (Oryza sativa) contains many iconic varieties and is widespread in the Indian subcontinent. Despite its economic and cultural importance, a high-quality reference genome is currently lacking, and the group's evolutionary history is not fully resolved. To address these gaps, we use long-read nanopore sequencing and assemble the genomes of two circum-basmati rice varieties. RESULTS: We generate two high-quality, chromosome-level reference genomes that represent the 12 chromosomes of Oryza. The assemblies show a contig N50 of 6.32 Mb and 10.53 Mb for Basmati 334 and Dom Sufid, respectively. Using our highly contiguous assemblies, we characterize structural variations segregating across circum-basmati genomes. We discover repeat expansions not observed in japonica-the rice group most closely related to circum-basmati-as well as the presence and absence variants of over 20 Mb, one of which is a circum-basmati-specific deletion of a gene regulating awn length. We further detect strong evidence of admixture between the circum-basmati and circum-aus groups. This gene flow has its greatest effect on chromosome 10, causing both structural variation and single-nucleotide polymorphism to deviate from genome-wide history. Lastly, population genomic analysis of 78 circum-basmati varieties shows three major geographically structured genetic groups: Bhutan/Nepal, India/Bangladesh/Myanmar, and Iran/Pakistan. CONCLUSION: The availability of high-quality reference genomes allows functional and evolutionary genomic analyses providing genome-wide evidence for gene flow between circum-aus and circum-basmati, describes the nature of circum-basmati structural variation, and reveals the presence/absence variation in this important and iconic rice variety group.

An inferred fitness consequence map of the rice genome.

The extent to which sequence variation impacts plant fitness is poorly understood. High-resolution maps detailing the constraint acting on the genome, especially in regulatory sites, would be beneficial as functional annotation of noncoding sequences remains sparse. Here, we present a fitness consequence (fitCons) map for rice (Oryza sativa). We inferred fitCons scores (ρ) for 246 inferred genome classes derived from nine functional genomic and epigenomic datasets, including chromatin accessibility, messenger RNA/small RNA transcription, DNA methylation, histone modifications and engaged RNA polymerase activity. These were integrated with genome-wide polymorphism and divergence data from 1,477 rice accessions and 11 reference genome sequences in the Oryzeae. We found ρ to be multimodal, with ~9% of the rice genome falling into classes where more than half of the bases would probably have a fitness consequence if mutated. Around 2% of the rice genome showed evidence of weak negative selection, frequently at candidate regulatory sites, including a novel set of 1,000 potentially active enhancer elements. This fitCons map provides perspective on the evolutionary forces associated with genome diversity, aids in genome annotation and can guide crop breeding programs.

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.

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.

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.

Engineering resistance to virus transmission.

Engineering plants for resistance to virus transmission by invertebrate vectors has lagged behind other forms of plant protection. Vectors typically transmit more than one virus. Thus, vector resistance could provide a wider range of protection than defenses directed solely against one virus or virus group. We discuss current knowledge of vector-host-virus interactions, the roles of viral gene products in host and vector manipulation, and the effects of semiochemicals on host-vector interactions, and how this knowledge could be employed to disrupt transmission dynamics. We also discuss how resistance to vectors could be generated through genetic engineering or gene editing or indirectly through use of biocontrol using plant-resident viruses that infect vectors.

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.

Multidrug transporters and organic anion transporting polypeptides protect insects against the toxic effects of cardenolides.

In the struggle against dietary toxins, insects are known to employ target site insensitivity, metabolic detoxification, and transporters that shunt away toxins. Specialized insects across six taxonomic orders feeding on cardenolide-containing plants have convergently evolved target site insensitivity via specific amino acid substitutions in the Na/K-ATPase. Nonetheless, in vitro pharmacological experiments have suggested a role for multidrug transporters (Mdrs) and organic anion transporting polypeptides (Oatps), which may provide a basal level of protection in both specialized and non-adapted insects. Because the genes coding for these proteins are evolutionarily conserved and in vivo genetic evidence in support of this hypothesis is lacking, here we used wildtype and mutant Drosophila melanogaster (Drosophila) in capillary feeder (CAFE) assays to quantify toxicity of three chemically diverse, medically relevant cardenolides. We examined multiple components of fitness, including mortality, longevity, and LD50, and found that, while the three cardenolides each stimulated feeding (i.e., no deterrence to the toxin), all decreased lifespan, with the most apolar cardenolide having the lowest LD50 value. Flies showed a clear non-monotonic dose response and experienced high levels of toxicity at the cardenolide concentration found in plants. At this concentration, both Mdr and Oatp knockout mutant flies died more rapidly than wildtype flies, and the mutants also experienced more adverse neurological effects on high-cardenolide-level diets. Our study further establishes Drosophila as a model for the study of cardenolide pharmacology and solidifies support for the hypothesis that multidrug and organic anion transporters are key players in insect protection against dietary cardenolides.