TraitRateProp: a web server for the detection of trait-dependent evolutionary rate shifts in sequence sites

Distinctive plastome evolution in carnivorous angiosperms

BACKGROUND

Independent origins of carnivory in multiple angiosperm families are fabulous examples of convergent evolution using a diverse array of life forms and habitats. Previous studies have indicated that carnivorous plants have distinct evolutionary trajectories of plastid genome (plastome) compared to their non-carnivorous relatives, yet the extent and general characteristics remain elusive.

RESULTS

We compared plastomes from 9 out of 13 carnivorous families and their non-carnivorous relatives to assess carnivory-associated evolutionary patterns. We identified inversions in all sampled Droseraceae species and four species of Utricularia, Pinguicula, Darlingtonia and Triphyophyllum. A few carnivores showed distinct shifts in inverted repeat boundaries and the overall repeat contents. Many ndh genes, along with some other genes, were independently lost in several carnivorous lineages. We detected significant substitution rate variations in most sampled carnivorous lineages. A significant overall substitution rate acceleration characterizes the two largest carnivorous lineages of Droseraceae and Lentibulariaceae. We also observe moderate substitution rates acceleration in many genes of Cephalotus follicularis, Roridula gorgonias, and Drosophyllum lusitanicum. However, only a few genes exhibit significant relaxed selection.

CONCLUSION

Our results indicate that the carnivory of plants have different effects on plastome evolution across carnivorous lineages. The complex mechanism under carnivorous habitats may have resulted in distinctive plastome evolution with conserved plastome in the Brocchinia hechtioides to strongly reconfigured plastomes structures in Droseraceae. Organic carbon obtained from prey and the efficiency of utilizing prey-derived nutrients might constitute possible explanation.

Convergent resistance to GABA receptor neurotoxins through plant-insect coevolution

The molecular mechanisms of coevolution between plants and insects remain elusive. GABA receptors are targets of many neurotoxic terpenoids, which represent the most diverse array of natural products known. Over deep evolutionary time, as plant terpene synthases diversified in plants, so did plant terpenoid defence repertoires. Here we show that herbivorous insects and their predators evolved convergent amino acid changing substitutions in duplicated copies of the Resistance to dieldrin (Rdl) gene that encodes the GABA receptor, and that the evolution of duplicated Rdl and terpenoid-resistant GABA receptors is associated with the diversification of moths and butterflies. These same substitutions also evolved in pests exposed to synthetic insecticides that target the GABA receptor. We used in vivo genome editing in Drosophila melanogaster to evaluate the fitness effects of each putative resistance mutation and found that pleiotropy both facilitates and constrains the evolution of GABA receptor resistance. The same genetic changes that confer resistance to terpenoids across 300 Myr of insect evolution have re-evolved in response to synthetic analogues over one human lifespan.

Mitonuclear conflict in a macaque species exhibiting phylogenomic discordance

Speciation and hybridization are intertwined processes in the study of evolution. Hybridization between sufficiently diverged populations can result in genomic conflict within offspring, causing reduced viability and fertility, thus increasing divergence between populations. Conflicts between mitochondrial and nuclear genes are increasingly found to play a role in this process in various systems. We examine the possibility of this conflict in the bear macaque, Macaca arctoides (Primates: Cercopithecidae), a primate species exhibiting mitonuclear discordance due to extensive hybridization with species in the sinica and fascicularis groups. Here, divergence, introgression and natural selection of mitonuclear genes (N = 160) relative to nuclear control genes (N = 144) were analysed to determine whether there are evolutionary processes involved in resolving the potential conflict caused by mitonuclear discordance. Nucleotide divergence of mitonuclear genes is increased relative to control nuclear genes between M. arctoides and the species sharing its nuclear ancestry (p = 0.007), consistent with genetic conflict. However, measures of introgression and selection do not identify large-scale co-introgression or co-evolution as means to resolve mitonuclear conflict. Nonetheless, mitochondrial tRNA synthetases stand out in analyses using dN/dS and extended branch lengths as potential targets of selection. The methodology implemented provides a framework that can be used to examine the effects of mitonuclear co-introgression and co-evolution on a genomic scale in a variety of systems.

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 fitness¹. Convergent evolution can be used as a guide to identify candidate mutations that underlie adaptive traits^2-4, and new genome editing technology is facilitating functional validation of these mutations in whole organisms^1,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 toxins^6-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 glycosides¹². Here we describe mutational paths involving three repeatedly changing amino acid sites (111, 119 and 122) in ATPα that are associated with cardiac glycoside specialization^13,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 lineage^11,15. We show in vivo, in vitro and in silico that the path conferred resistance and target-site insensitivity to cardiac glycosides¹⁶, 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 predators^17-19. The order in which the substitutions evolved was explained by amelioration of antagonistic pleiotropy through epistasis^13,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 communities^{2,6-11,15,17-19}.

Unprecedented Parallel Photosynthetic Losses in a Heterotrophic Orchid Genus

Heterotrophic plants are evolutionary experiments in genomic, morphological, and physiological change. Yet, genomic sampling gaps exist among independently derived heterotrophic lineages, leaving unanswered questions about the process of genome modification. Here, we have sequenced complete plastid genomes for all species of the leafless orchid genus Hexalectris, including multiple individuals for most, and leafy relatives Basiphyllaea and Bletia. Our objectives are to determine the number of independent losses of photosynthesis and to test hypotheses on the process of genome degradation as a result of relaxed selection. We demonstrate four to five independent losses of photosynthesis in Hexalectris based on degradation of the photosynthetic apparatus, with all but two species displaying evidence of losses, and variation in gene loss extending below the species level. Degradation in the atp complex is advanced in Hexalectris warnockii, whereas only minimal degradation (i.e., physical loss) has occurred among some "housekeeping" genes. We find genomic rearrangements, shifts in Inverted Repeat boundaries including complete loss in one accession of H. arizonica, and correlations among substitutional and genomic attributes. Our unprecedented finding of multiple, independent transitions to a fully mycoheterotrophic lifestyle in a single genus reveals that the number of such transitions among land plants is likely underestimated. This study underscores the importance of dense taxon sampling, which is highly informative for advancing models of genome evolution in heterotrophs. Mycoheterotrophs such as Hexalectris provide forward-genetic opportunities to study the consequences of radical genome evolution beyond what is possible with mutational studies in model organisms alone.

Genomic Signals of Adaptation towards Mutualism and Sociality in Two Ambrosia Beetle Complexes

Mutualistic symbiosis and eusociality have developed through gradual evolutionary processes at different times in specific lineages. Like some species of termites and ants, ambrosia beetles have independently evolved a mutualistic nutritional symbiosis with fungi, which has been associated with the evolution of complex social behaviors in some members of this group. We sequenced the transcriptomes of two ambrosia complexes (Euwallacea sp. near fornicatus⁻Fusarium euwallaceae and Xyleborus glabratus⁻Raffaelea lauricola) to find evolutionary signatures associated with mutualism and behavior evolution. We identified signatures of positive selection in genes related to nutrient homeostasis; regulation of gene expression; development and function of the nervous system, which may be involved in diet specialization; behavioral changes; and social evolution in this lineage. Finally, we found convergent changes in evolutionary rates of proteins across lineages with phylogenetically independent origins of sociality and mutualism, suggesting a constrained evolution of conserved genes in social species, and an evolutionary rate acceleration related to changes in selective pressures in mutualistic lineages.

Investigating the NAD-ME biochemical pathway within C4 grasses using transcript and amino acid variation in C4 photosynthetic genes

Expanding knowledge of the C₄ photosynthetic pathway can provide key information to aid biological improvements to crop photosynthesis and yield. While the C₄ NADP-ME pathway is well characterised, there is increasing agricultural and bioengineering interest in the comparably understudied NAD-ME and PEPCK pathways. Within this study, a systematic identification of key differences across species has allowed us to investigate the evolution of C₄-recruited genes in one C₃ and eleven C₄ grasses (Poaceae) spanning two independent origins of C₄ photosynthesis. We present evidence for C₄-specific paralogs of NAD-malic enzyme 2, MPC1 and MPC2 (mitochondrial pyruvate carriers) via increased transcript abundance and associated rates of evolution, implicating them as genes recruited to perform C₄ photosynthesis within NAD-ME and PEPCK subtypes. We then investigate the localisation of AspAT across subtypes, using novel and published evidence to place AspAT3 in both the cytosol and peroxisome. Finally, these findings are integrated with transcript abundance of previously identified C₄ genes to provide an updated model for C₄ grass NAD-ME and PEPCK photosynthesis. This updated model allows us to develop on the current understanding of NAD-ME and PEPCK photosynthesis in grasses, bolstering our efforts to understand the evolutionary 'path to C₄' and improve C₄ photosynthesis.

Modeling a trait-dependent diversification process coupled with molecular evolution on a random species tree

Understanding the evolution of binary traits, which affects the birth and survival of species and also the rate of molecular evolution, remains challenging. In this work, we present a probabilistic modeling framework for binary trait, random species trees, in which the number of species and their traits are represented by an asymmetric, two-type, continuous time Markov branching process. The model involves a number of different parameters describing both character and molecular evolution on the so-called 'reduced' tree, consisting of only extant species at the time of observation. We expand our model by considering the impact of binary traits on dN/dS, the normalized ratio of nonsynonymous to synonymous substitutions. We also develop mechanisms which enable us to understand the substitution rates on a phylogenetic tree with regards to the observed traits. The properties obtained from the model are illustrated with a phylogeny of outcrossing and selfing plant species, which allows us to investigate not only the branching tree rates, but also the molecular rates and the intensity of selection.