Positive diversifying selection is a pervasive adaptive force throughout the Drosophila radiation

Gene expression, evolution, and the genetics of electrosensing in the smalltooth sawfish, Pristis pectinata

Sawfishes (Pristidae) are large, highly threatened rays named for their tooth-studded rostrum, which is used for prey sensing and capture. Of all five species, the smalltooth sawfish, Pristis pectinata, has experienced the greatest decline in range, currently found in only ~20% of its historic range. To better understand the genetic underpinnings of these taxonomically and morphologically unique animals, we collected transcriptomic data from several tissue types, mapped them to the recently completed reference genome, and contrasted the patterns observed with comparable data from other elasmobranchs. Evidence of positive selection was detected in 79 genes in P. pectinata, several of which are involved in growth factor/receptor tyrosine kinase signaling and body symmetry and may be related to the unique morphology of sawfishes. Changes in these genes may impact cellular responses to environmental conditions such as temperature, dissolved oxygen, and salinity. Data acquired also allow for examination of the molecular components of P. pectinata electrosensory systems, which are highly developed in sawfishes and have likely been influential in their evolutionary success.

The evolution of intramitochondriality in Midichloria bacteria

Midichloria spp. are intracellular bacterial symbionts of ticks. Representatives of this genus colonise mitochondria in the cells of their hosts. To shed light on this unique interaction we evaluated the presence of an intramitochondrial localization for three Midichloria in the respective tick host species and generated eight high-quality draft genomes and one closed genome, showing that this trait is non-monophyletic, either due to losses or multiple acquisitions. Comparative genomics supports the first hypothesis, as the genomes of non-mitochondrial symbionts are reduced subsets of those capable of colonising the organelles. We detect genomic signatures of mitochondrial tropism, including the differential presence of type IV secretion system and flagellum, which could allow the secretion of unique effectors and/or direct interaction with mitochondria. Other genes, including adhesion molecules, proteins involved in actin polymerisation, cell wall and outer membrane proteins, are only present in mitochondrial symbionts. The bacteria could use these to manipulate host structures, including mitochondrial membranes, to fuse with the organelles or manipulate the mitochondrial network.

Evolutionary dynamics of genome size and content during the adaptive radiation of Heliconiini butterflies

Heliconius butterflies, a speciose genus of Müllerian mimics, represent a classic example of an adaptive radiation that includes a range of derived dietary, life history, physiological and neural traits. However, key lineages within the genus, and across the broader Heliconiini tribe, lack genomic resources, limiting our understanding of how adaptive and neutral processes shaped genome evolution during their radiation. Here, we generate highly contiguous genome assemblies for nine Heliconiini, 29 additional reference-assembled genomes, and improve 10 existing assemblies. Altogether, we provide a dataset of annotated genomes for a total of 63 species, including 58 species within the Heliconiini tribe. We use this extensive dataset to generate a robust and dated heliconiine phylogeny, describe major patterns of introgression, explore the evolution of genome architecture, and the genomic basis of key innovations in this enigmatic group, including an assessment of the evolution of putative regulatory regions at the Heliconius stem. Our work illustrates how the increased resolution provided by such dense genomic sampling improves our power to generate and test gene-phenotype hypotheses, and precisely characterize how genomes evolve.

Evolutionary Divergence and Radula Diversification in Two Ecomorphs from an Adaptive Radiation of Freshwater Snails

(1) Background: Adaptive diversification of complex traits plays a pivotal role in the evolution of organismal diversity. In the freshwater snail genus Tylomelania, adaptive radiations were likely promoted by trophic specialization via diversification of their key foraging organ, the radula. (2) Methods: To investigate the molecular basis of radula diversification and its contribution to lineage divergence, we used tissue-specific transcriptomes of two sympatric Tylomelania sarasinorum ecomorphs. (3) Results: We show that ecomorphs are genetically divergent lineages with habitat-correlated abundances. Sequence divergence and the proportion of highly differentially expressed genes are significantly higher between radula transcriptomes compared to the mantle and foot. However, the same is not true when all differentially expressed genes or only non-synonymous SNPs are considered. Finally, putative homologs of some candidate genes for radula diversification (hh, arx, gbb) were also found to contribute to trophic specialization in cichlids and Darwin’s finches. (4) Conclusions: Our results are in line with diversifying selection on the radula driving Tylomelania ecomorph divergence and indicate that some molecular pathways may be especially prone to adaptive diversification, even across phylogenetically distant animal groups.

Differential gene expression in Drosophila melanogaster and D. nigrosparsa infected with the same Wolbachia strain

Wolbachia are maternally inherited endosymbionts that infect nearly half of all arthropod species. Wolbachia manipulate their hosts to maximize their transmission, but they can also provide benefits such as nutrients and resistance against viruses to their hosts. The Wolbachia strain wMel was recently found to increase locomotor activities and possibly trigger cytoplasmic incompatibility in the transinfected fly Drosophila nigrosparsa. Here, we investigated, in females of both D. melanogaster and D. nigrosparsa, the gene expression between animals uninfected and infected with wMel, using RNA sequencing to see if the two Drosophila species respond to the infection in the same or different ways. A total of 2164 orthologous genes were used. The two fly species responded to the infection in different ways. Significant changes shared by the fly species belong to the expression of genes involved in processes such as oxidation-reduction process, iron-ion binding, and voltage-gated potassium-channel activity. We discuss our findings also in the light of how Wolbachia survive within both the native and the novel host.

Molecular Evolution of Ecological Specialisation: Genomic Insights from the Diversification of Murine Rodents

Adaptive radiations are characterized by the diversification and ecological differentiation of species, and replicated cases of this process provide natural experiments for understanding the repeatability and pace of molecular evolution. During adaptive radiation, genes related to ecological specialization may be subject to recurrent positive directional selection. However, it is not clear to what extent patterns of lineage-specific ecological specialization (including phenotypic convergence) are correlated with shared signatures of molecular evolution. To test this, we sequenced whole exomes from a phylogenetically dispersed sample of 38 murine rodent species, a group characterized by multiple, nested adaptive radiations comprising extensive ecological and phenotypic diversity. We found that genes associated with immunity, reproduction, diet, digestion, and taste have been subject to pervasive positive selection during the diversification of murine rodents. We also found a significant correlation between genome-wide positive selection and dietary specialization, with a higher proportion of positively selected codon sites in derived dietary forms (i.e., carnivores and herbivores) than in ancestral forms (i.e., omnivores). Despite striking convergent evolution of skull morphology and dentition in two distantly related worm-eating specialists, we did not detect more genes with shared signatures of positive or relaxed selection than in a nonconvergent species comparison. Although a small number of the genes we detected can be incidentally linked to craniofacial morphology or diet, protein-coding regions are unlikely to be the primary genetic basis of this complex convergent phenotype. Our results suggest a link between positive selection and derived ecological phenotypes, and highlight specific genes and general functional categories that may have played an integral role in the extensive and rapid diversification of murine rodents.

Genomic Signature of Shifts in Selection in a Subalpine Ant and Its Physiological Adaptations

Understanding how organisms adapt to extreme environments is fundamental and can provide insightful case studies for both evolutionary biology and climate-change biology. Here, we take advantage of the vast diversity of lifestyles in ants to identify genomic signatures of adaptation to extreme habitats such as high altitude. We hypothesized two parallel patterns would occur in a genome adapting to an extreme habitat: 1) strong positive selection on genes related to adaptation and 2) a relaxation of previous purifying selection. We tested this hypothesis by sequencing the high-elevation specialist Tetramorium alpestre and four other phylogenetically related species. In support of our hypothesis, we recorded a strong shift of selective forces in T. alpestre, in particular a stronger magnitude of diversifying and relaxed selection when compared with all other ants. We further disentangled candidate molecular adaptations in both gene expression and protein-coding sequence that were identified by our genome-wide analyses. In particular, we demonstrate that T. alpestre has 1) a higher level of expression for stv and other heat-shock proteins in chill-shock tests and 2) enzymatic enhancement of Hex-T1, a rate-limiting regulatory enzyme that controls the entry of glucose into the glycolytic pathway. Together, our analyses highlight the adaptive molecular changes that support colonization of high-altitude environments.

Gene diversity explains variation in biological features of insect killing fungus, Beauveria bassiana

Beauveria bassiana is a species complex whose isolates show considerable natural genetic variability. However, little is known about how this genetic diversity affects the fungus performance. Herein, we characterized the diversity of genes involved in various mechanisms of the infective cycle of 42 isolates that have different growth rates, thermotolerance and virulence. The analysed genes showed general genetic diversity measured as non-synonymous changes (NSC) and copy number variation (CNV), with most of them being subjected to positive episodic diversifying selection. Correlation analyses between NSC or CNV and the isolate virulence, thermotolerance and growth rate revealed that various genes shaped the biological features of the fungus. Lectin-like, mucin signalling, Biotrophy associated and chitinase genes NSCs correlated with the three biological features of B. bassiana. In addition, other genes (i.e. DNA photolyase and cyclophilin B) that had relatively conserved sequences, had variable CNs across the isolates which were correlated with the variability of either virulence or thermotolerance of B. bassiana isolates. The data obtained is important for a better understanding of population structure, ecological and potential impact when isolates are used as mycoinsecticides and can justify industrialization of new isolates.

Elevated rates of positive selection drive the evolution of pestiferousness in the Colorado potato beetle (Leptinotarsa decemlineata, Say)

Contextualizing evolutionary history and identifying genomic features of an insect that might contribute to its pest status is important in developing early detection and control tactics. In order to understand the evolution of pestiferousness, which we define as the accumulation of traits that contribute to an insect population's success in an agroecosystem, we tested the importance of known genomic properties associated with rapid adaptation in the Colorado potato beetle (CPB), Leptinotarsa decemlineata Say. Within the leaf beetle genus Leptinotarsa, only CPB, and a few populations therein, has risen to pest status on cultivated nightshades, Solanum. Using whole genomes from ten closely related Leptinotarsa species native to the United States, we reconstructed a high-quality species tree and used this phylogenetic framework to assess evolutionary patterns in four genomic features of rapid adaptation: standing genetic variation, gene family expansion and contraction, transposable element abundance and location, and positive selection at protein-coding genes. Throughout approximately 20 million years of history, Leptinotarsa species show little evidence of gene family turnover and transposable element variation. However, there is a clear pattern of CPB experiencing higher rates of positive selection on protein-coding genes. We determine that these rates are associated with greater standing genetic variation due to larger effective population size, which supports the theory that the demographic history contributes to rates of protein evolution. Furthermore, we identify a suite of coding genes under positive selection that are putatively associated with pestiferousness in the Colorado potato beetle lineage. They are involved in the biological processes of xenobiotic detoxification, chemosensation and hormone function.

Characterization of two MHC II genes (DOB, DRB) in white-tailed deer (Odocoileus virginianus)

BACKGROUND

The major histocompatibility complex (MHC) is responsible for detecting and addressing foreign pathogens inside the body. While the general structure of MHC genes is relatively well conserved among mammalian species, it is notably different among ruminants due to a chromosomal inversion that splits MHC type II genes into two subregions (IIa, IIb). Recombination rates are reportedly high between these subregions, and a lack of linkage has been documented in domestic ruminants. However, no study has yet examined the degree of linkage between these subregions in a wild ruminant. The white-tailed deer (Odocoileus virginianus), a popular ruminant of the Cervidae family, is habitually plagued by pathogens in its natural environment (e.g. Haemonchus contortus, Elaeophora). Due to the association between MHC haplotypes and disease susceptibility, a deeper understanding of MHC polymorphism and linkage between MHC genes can further aid in this species' successful management. We sequenced MHC-DRB exon 2 (IIa) and MHC-DOB exon 2 (IIb) on the MiSeq platform from an enclosed white-tailed deer population located in Alabama.

RESULTS

We identified 12 new MHC-DRB alleles, and resampled 7 alleles, which along with other published alleles brings the total number of documented alleles in white-tailed deer to 30 for MHC-DRB exon 2. The first examination of MHC-DOB in white-tailed deer found significantly less polymorphism (11 alleles), as was expected of a non-classical MHC gene. While MHC-DRB was found to be under positive, diversifying selection, MHC-DOB was found to be under purifying selection for white-tailed deer. We found no significant linkage disequilibrium between MHC-DRB and MHC-DOB, suggesting that these loci are unlikely to be closely linked.

CONCLUSIONS

Overall, this study identified 12 new MHC-DRB exon 2 alleles and characterized a new, non-classical, MHC II gene (MHC-DOB) for white-tailed deer. We also found a lack of significant linkage between these two loci, which supports previous findings of a chromosomal inversion within the MHC type II gene region in ruminants, and suggests that white-tailed deer may have a recombination hotspot between these MHC regions similar to that found for Bos taurus.

Integrated structural and evolutionary analysis reveals common mechanisms underlying adaptive evolution in mammals

Understanding the molecular basis of adaptation to the environment is a central question in evolutionary biology, yet linking detected signatures of positive selection to molecular mechanisms remains challenging. Here we demonstrate that combining sequence-based phylogenetic methods with structural information assists in making such mechanistic interpretations on a genomic scale. Our integrative analysis shows that positively selected sites tend to colocalize on protein structures and that positively selected clusters are found in functionally important regions of proteins, indicating that positive selection can contravene the well-known principle of evolutionary conservation of functionally important regions. This unexpected finding, along with our discovery that positive selection acts on structural clusters, opens previously unexplored strategies for the development of better models of protein evolution. Remarkably, proteins where we detect the strongest evidence of clustering belong to just two functional groups: Components of immune response and metabolic enzymes. This gives a coherent picture of pathogens and xenobiotics as important drivers of adaptive evolution of mammals.

Major range loss predicted from lack of heat adaptability in an alpine Drosophila species

Climate warming is threatening biodiversity worldwide. Climate specialists such as alpine species are especially likely to be vulnerable. Adaptation by rapid evolution is the only long-term option for survival of many species, but the adaptive evolutionary potential of heat resistance has not been assessed in an alpine invertebrate. Here, we show that the alpine fly Drosophila nigrosparsa cannot readily adapt to heat stress. Heat-exposed flies from a regime with increased ambient temperature and a regime with increased temperature plus artificial selection for heat tolerance were less heat tolerant than the control group. Increased ambient temperature affected negatively both fitness and competitiveness. Ecological niche models predicted the loss of three quarters of the climatically habitable areas of this fly by the end of this century. Our findings suggest that, alongside with other climate specialists, species from mountainous regions are highly vulnerable to climate warming and unlikely to adapt through evolutionary genetic changes.

Large-Scale Comparative Analysis of Codon Models Accounting for Protein and Nucleotide Selection

There are numerous sources of variation in the rate of synonymous substitutions inside genes, such as direct selection on the nucleotide sequence, or mutation rate variation. Yet scans for positive selection rely on codon models which incorporate an assumption of effectively neutral synonymous substitution rate, constant between sites of each gene. Here we perform a large-scale comparison of approaches which incorporate codon substitution rate variation and propose our own simple yet effective modification of existing models. We find strong effects of substitution rate variation on positive selection inference. More than 70% of the genes detected by the classical branch-site model are presumably false positives caused by the incorrect assumption of uniform synonymous substitution rate. We propose a new model which is strongly favored by the data while remaining computationally tractable. With the new model we can capture signatures of nucleotide level selection acting on translation initiation and on splicing sites within the coding region. Finally, we show that rate variation is highest in the highly recombining regions, and we propose that recombination and mutation rate variation, such as high CpG mutation rate, are the two main sources of nucleotide rate variation. Although we detect fewer genes under positive selection in Drosophila than without rate variation, the genes which we detect contain a stronger signal of adaptation of dynein, which could be associated with Wolbachia infection. We provide software to perform positive selection analysis using the new model.

Hawaiian Drosophila as an Evolutionary Model Clade: Days of Future Past

The Hawaiian Drosophila have been a model system for evolutionary, ecological, and ethological studies since the inception of the Hawaiian Drosophila Project in the 1960s. Here we review the past and present research on this incredible lineage and provide a prospectus for future directions on genomics and microbial interactions. While the number of publications on this group has waxed and waned over the years, we assert that recent systematic, biogeographic, and ecological studies have reinvigorated Hawaiian Drosophila as an evolutionary model system. The characteristics that distinguish good model clades from good model organisms (e.g., Drosophila melanogaster) are somewhat different so we first define what constitutes a good evolutionary model. We argue that the Hawaiian Drosophila possess many desired aspects of a good evolutionary model, describe how this group of geographically isolated flies have been used in the past, and propose some exciting avenues for future evolutionary research on this diverse, dynamic clade of Drosophila.

Life-history traits and physiological limits of the alpine fly Drosophila nigrosparsa (Diptera: Drosophilidae): A comparative study

Interspecific variation in life-history traits and physiological limits can be linked to the environmental conditions species experience, including climatic conditions. As alpine environments are particularly vulnerable under climate change, we focus on the montane-alpine fly Drosophila nigrosparsa. Here, we characterized some of its life-history traits and physiological limits and compared these with those of other drosophilids, namely Drosophila hydei, Drosophila melanogaster, and Drosophila obscura. We assayed oviposition rate, longevity, productivity, development time, larval competitiveness, starvation resistance, and heat and cold tolerance. Compared with the other species assayed, D. nigrosparsa is less fecund, relatively long-living, starvation susceptible, cold adapted, and surprisingly well heat adapted. These life-history characteristics provide insights into invertebrate adaptations to alpine conditions which may evolve under ongoing climate change.

Identifying the minimum number of microsatellite loci needed to assess population genetic structure: A case study in fly culturing

Small, isolated populations are constantly threatened by loss of genetic diversity due to drift. Such situations are found, for instance, in laboratory culturing. In guarding against diversity loss, monitoring of potential changes in population structure is paramount; this monitoring is most often achieved using microsatellite markers, which can be costly in terms of time and money when many loci are scored in large numbers of individuals. Here, we present a case study reducing the number of microsatellites to the minimum necessary to correctly detect the population structure of two Drosophila nigrosparsa populations. The number of loci was gradually reduced from 11 to 1, using the Allelic Richness (AR) and Private Allelic Richness (PAR) as criteria for locus removal. The effect of each reduction step was evaluated by the number of genetic clusters detectable from the data and by the allocation of individuals to the clusters; in the latter, excluding ambiguous individuals was tested to reduce the rate of incorrect assignments. We demonstrate that more than 95% of the individuals can still be correctly assigned when using eight loci and that the major population structure is still visible when using two highly polymorphic loci. The differences between sorting the loci by AR and PAR were negligible. The method presented here will most efficiently reduce genotyping costs when small sets of loci (“core sets”) for long-time use in large-scale population screenings are compiled.