Cryptic diversity hides host and habitat specialization in a gorgonian‐algal symbiosis

Shallow water anthozoans, the major builders of modern coral reefs, enhance their metabolic and calcification rates with algal symbionts. Controversy exists over whether these anthozoan–algae associations are flexible over the lifetimes of individual hosts, promoting acclimative plasticity, or are closely linked, such that hosts and symbionts co‐evolve across generations. Given the diversity of algal symbionts and the morphological plasticity of many host species, cryptic variation within either partner could potentially confound studies of anthozoan‐algal associations. Here, we used ribosomal, organelle and nuclear sequences, along with microsatellite variation, to study the relationship between lineages of a common Caribbean gorgonian and its algal symbionts. The gorgonian Eunicea flexuosa is a broadcast spawner, composed of two recently diverged, genetically distinct lineages largely segregated by depth. We sampled colonies of the two lineages across depth gradients at three Caribbean locations. We find that each host lineage is associated with a unique Symbiodinium B1/184 phylotype. This relationship between host and symbiont is maintained when host colonies are reciprocally transplanted, although cases of within phylotype switching were also observed. Even when the phylotypes of both partners are present at intermediate depths, the specificity between host and symbiont lineages remained absolute. Unrecognized cryptic diversity may mask host‐symbiont specificity and change the inference of evolutionary processes in mutualistic associations. Symbiotic specificity thus likely contributes to the ecological divergence of the two partners, generating species diversity within coral reefs.     

Genomic signatures of evolution in Nautilus—An endangered living fossil

Living fossils are survivors of previously more diverse lineages that originated millions of years ago and persisted with little morphological change. Therefore, living fossils are model organisms to study both long‐term and ongoing adaptation and speciation processes. However, many aspects of living fossil evolution and their persistence in the modern world remain unclear. Here, we investigate three major aspects of the evolutionary history of living fossils: cryptic speciation, population genetics and effective population sizes, using members of the genera Nautilus and Allonautilus as classic examples of true living fossils. For this, we analysed genomewide ddRAD‐Seq data for all six currently recognized nautiloid species throughout their distribution range. Our analyses identified three major allopatric Nautilus clades: a South Pacific clade, subdivided into three subclades with no signs of admixture between them; a Coral Sea clade, consisting of two genetically distinct populations with significant admixture; and a widespread Indo‐Pacific clade, devoid of significant genetic substructure. Within these major clades, we detected five Nautilus groups, which likely correspond to five distinct species. With the exception of Nautilus macromphalus, all previously described species are at odds with genomewide data, testifying to the prevalence of cryptic species among living fossils. Detailed F_ST analyses further revealed significant genome‐wide and locus‐specific signatures of selection between species and differentiated populations, which is demonstrated here for the first time in a living fossil. Finally, approximate Bayesian computation (ABC) simulations suggest large effective population sizes, which may explain the low levels of population differentiation commonly observed in living fossils.     

Genomewide signatures of selection in Epichloë reveal candidate genes for host specialization

Host specialization is a key process in ecological divergence and speciation of plant‐associated fungi. The underlying determinants of host specialization are generally poorly understood, especially in endophytes, which constitute one of the most abundant components of the plant microbiome. We addressed the genetic basis of host specialization in two sympatric subspecies of grass‐endophytic fungi from the Epichloë typhina complex: subsp. typhina and clarkii. The life cycle of these fungi entails unrestricted dispersal of gametes and sexual reproduction before infection of a new host, implying that the host imposes a selective barrier on viability of the progeny. We aimed to detect genes under divergent selection between subspecies, experiencing restricted gene flow due to adaptation to different hosts. Using pooled whole‐genome sequencing data, we combined F_ST and D_XY population statistics in genome scans and detected 57 outlier genes showing strong differentiation between the two subspecies. Genomewide analyses of nucleotide diversity (π), Tajima's D and dN/dS ratios indicated that these genes have evolved under positive selection. Genes encoding secreted proteins were enriched among the genes showing evidence of positive selection, suggesting that molecular plant–fungus interactions are strong drivers of endophyte divergence. We focused on five genes encoding secreted proteins, which were further sequenced in 28 additional isolates collected across Europe to assess genetic variation in a larger sample size. Signature of positive selection in these isolates and putative identification of pathogenic function supports our findings that these genes represent strong candidates for host specialization determinants in Epichloë endophytes. Our results highlight the role of secreted proteins as key determinants of host specialization.     

Genomic signatures of rapid adaptive evolution in the bluespotted cornetfish,a Mediterranean Lessepsian invader

Biological invasions are increasingly creating ecological and economical problems both on land and in aquatic environments. For over a century, the Mediterranean Sea has steadily been invaded by Indian Ocean/Red Sea species (called Lessepsian invaders) via the Suez Canal, with a current estimate of ~450 species. The bluespotted cornetfish, Fistularia commersonii, considered a ‘Lessepsian sprinter’, entered the Mediterranean in 2000 and by 2007 had spread through the entire basin from Israel to Spain. The situation is unique and interesting both because of its unprecedented rapidity and by the fact that it took this species c. 130 years to immigrate into the Mediterranean. Using genome scans, with restriction site‐associated DNA (RAD) sequencing, we evaluated neutral and selected genomic regions for Mediterranean vs. Red Sea cornetfish individuals. We found that few fixed neutral changes were detectable among populations. However, almost half of the genes associated with the 47 outlier loci (potentially under selection) were related to disease resistance and osmoregulation. Due to the short time elapsed from the beginning of the invasion to our sampling, we interpret these changes as signatures of rapid adaptation that may be explained by several mechanisms including preadaptation and strong local selection. Such genomic regions are therefore good candidates to further study their role in invasion success.     

Patterns of ecological specialization among microbial populations in the Red Sea and diverse oligotrophic marine environments

Large swaths of the nutrient‐poor surface ocean are dominated numerically by cyanobacteria (Prochlorococcus), cyanobacterial viruses (cyanophage), and alphaproteobacteria (SAR11). How these groups thrive in the diverse physicochemical environments of different oceanic regions remains poorly understood. Comparative metagenomics can reveal adaptive responses linked to ecosystem‐specific selective pressures. The Red Sea is well‐suited for studying adaptation of pelagic‐microbes, with salinities, temperatures, and light levels at the extreme end for the surface ocean, and low nutrient concentrations, yet no metagenomic studies have been done there. The Red Sea (high salinity, high light, low N and P) compares favorably with the Mediterranean Sea (high salinity, low P), Sargasso Sea (low P), and North Pacific Subtropical Gyre (high light, low N). We quantified the relative abundance of genetic functions among Prochlorococcus, cyanophage, and SAR11 from these four regions. Gene frequencies indicate selection for phosphorus acquisition (Mediterranean/Sargasso), DNA repair and high‐light responses (Red Sea/Pacific Prochlorococcus), and osmolyte C1 oxidation (Red Sea/Mediterranean SAR11). The unexpected connection between salinity‐dependent osmolyte production and SAR11 C1 metabolism represents a potentially major coevolutionary adaptation and biogeochemical flux. Among Prochlorococcus and cyanophage, genes enriched in specific environments had ecotype distributions similar to nonenriched genes, suggesting that inter‐ecotype gene transfer is not a major source of environment‐specific adaptation. Clustering of metagenomes using gene frequencies shows similarities in populations (Red Sea with Pacific, Mediterranean with Sargasso) that belie their geographic distances. Taken together, the genetic functions enriched in specific environments indicate competitive strategies for maintaining carrying capacity in the face of physical stressors and low nutrient availability.     

Genomic and phenotypic signatures of climate adaptation in an Anolis lizard

Integrated knowledge on phenotype, physiology, and genomic adaptations is required to understand the effects of climate on evolution. The functional genomic basis of organismal adaptation to changes in the abiotic environment, its phenotypic consequences, and its possible convergence across vertebrates are still understudied. In this study, we use a comparative approach to verify predicted gene functions for vertebrate thermal adaptation with observed functions underlying repeated genomic adaptations in response to elevation in the lizard Anolis cybotes. We establish a direct link between recurrently evolved phenotypes and functional genomics of altitude‐related climate adaptation in three highland and lowland populations in the Dominican Republic. We show that across vertebrates, genes contained in this interactome are expressed within the brain, the endocrine system, and during development. These results are relevant to elucidate the effect of global climate change across vertebrates and might aid in furthering insight into gene–environment relationships under disturbances to homeostasis.     

Familial social structure and socially driven genetic differentiation in Hawaiian short‐finned pilot whales

Social structure can have a significant impact on divergence and evolution within species, especially in the marine environment, which has few environmental boundaries to dispersal. On the other hand, genetic structure can affect social structure in many species, through an individual preference towards associating with relatives. One social species, the short‐finned pilot whale (Globicephala macrorhynchus), has been shown to live in stable social groups for periods of at least a decade. Using mitochondrial control sequences from 242 individuals and single nucleotide polymorphisms from 106 individuals, we examine population structure among geographic and social groups of short‐finned pilot whales in the Hawaiian Islands, and test for links between social and genetic structure. Our results show that there are at least two geographic populations in the Hawaiian Islands: a Main Hawaiian Islands (MHI) population and a Northwestern Hawaiian Islands/Pelagic population (F_ST and Φ_ST p < .001), as well as an eastern MHI community and a western MHI community (F_ST p = .009). We find genetically driven social structure, or high relatedness among social units and clusters (p < .001), and a positive relationship between relatedness and association between individuals (p < .0001). Further, socially organized clusters are genetically distinct, indicating that social structure drives genetic divergence within the population, likely through restricted mate selection (F_ST p = .05). This genetic divergence among social groups can make the species less resilient to anthropogenic or ecological disturbance. Conservation of this species therefore depends on understanding links among social structure, genetic structure and ecological variability within the species.     

Comparative genetics of Na+/K+‐ATPase in monarch butterfly populations with varying host plant toxicity

Herbivores have evolved numerous behavioural and physiological adaptations to host plants; however, molecular adaptations are still poorly understood. One well‐studied case comprises the specialist insects that feed on cardenolide‐containing plants. Here, convergent molecular evolution in the Na⁺/K⁺‐ATPase results in a reduced sensitivity to cardenolides across four insect orders. Because different plant species and genotypes differ in toxicity, Na⁺/K⁺‐ATPase may be under differential selection from geographically varying host plants. We examined the α subunit of Na⁺/K⁺‐ATPase in monarch butterflies (Danaus plexippus) from six worldwide populations to test whether differences in their host plant chemistry result in local adaptation at the molecular level. Although our study revealed multiple synonymous changes, we did not find these to be population‐specific, nor did we identify nonsynonymous changes. Additionally, we compared the amino acid sequence of this subunit across 19 species. We identified two novel changes at sites 836 (K836N) and 840 (E840R) in the αM7‐αM8 regions in the genus Danaus. Although previous studies focused on the first two trans‐membrane domains, C‐terminal domains may also interact with cardenolides. These results reveal a lack of molecular evolution of Na⁺/K⁺‐ATPase at the population level, and call for additional attention regarding the C‐terminal regions of this important enzyme.     

Role of host genetics and heat‐tolerant algal symbionts in sustaining populations of the endangered coral Orbicella faveolata in the Florida Keys with ocean warming

Identifying which factors lead to coral bleaching resistance is a priority given the global decline of coral reefs with ocean warming. During the second year of back‐to‐back bleaching events in the Florida Keys in 2014 and 2015, we characterized key environmental and biological factors associated with bleaching resilience in the threatened reef‐building coral Orbicella faveolata. Ten reefs (five inshore, five offshore, 179 corals total) were sampled during bleaching (September 2015) and recovery (May 2016). Corals were genotyped with 2bRAD and profiled for algal symbiont abundance and type. O. faveolata at the inshore sites, despite higher temperatures, demonstrated significantly higher bleaching resistance and better recovery compared to offshore. The thermotolerant Durusdinium trenchii (formerly Symbiondinium trenchii) was the dominant endosymbiont type region‐wide during initial (78.0% of corals sampled) and final (77.2%) sampling; >90% of the nonbleached corals were dominated by D. trenchii. 2bRAD host genotyping found no genetic structure among reefs, but inshore sites showed a high level of clonality. While none of the measured environmental parameters were correlated with bleaching, 71% of variation in bleaching resistance and 73% of variation in the proportion of D. trenchii was attributable to differences between genets, highlighting the leading role of genetics in shaping natural bleaching patterns. Notably, D. trenchii was rarely dominant in O. faveolata from the Florida Keys in previous studies, even during bleaching. The region‐wide high abundance of D. trenchii was likely driven by repeated bleaching associated with the two warmest years on record for the Florida Keys (2014 and 2015). On inshore reefs in the Upper Florida Keys, O. faveolata was most abundant, had the highest bleaching resistance, and contained the most corals dominated by D. trenchii, illustrating a causal link between heat tolerance and ecosystem resilience with global change.     

Landscape genomic signatures indicate reduced gene flow and forest‐associated adaptive divergence in an endangered neotropical turtle

Human‐induced transformations of ecosystems usually result in fragmented populations subject to increased extinction risk. Fragmentation is also often associated with novel environmental heterogeneity, which in combination with restricted gene flow may increase the opportunity for local adaptation. To manage at‐risk populations in these landscapes, it is important to understand how gene flow is changing, and how populations respond to habitat loss. We conducted a landscape genomics analysis using Restriction‐site Associated DNA sequencing to investigate the evolutionary response of the critically endangered Dahl's Toad‐headed turtle (Mesoclemmys dahli) to severe habitat modification. The species has lost almost all of its natural habitat in the southwestern part of its range and about 70% in the northeast. Based on least cost path analysis across different resistance surfaces for 3,211 SNPs, we found that the landscape matrix is restricting gene flow, causing the fragmentation of the species into at least six populations. Genome scans and allele‐environment association analyses indicate that the population fragments in the deforested grasslands of the southwest are adaptively different from those in the more forested northeast. Populations in areas with no forest had low levels of adaptive genetic diversity and the fixation of ancestrally‐polymorphic SNPs, consistent with directional selection in this novel environment. Our results suggest that this forest‐stream specialist is adapting to pond‐grassland conditions, but it is also suffering from negative consequences of habitat loss, including genetic erosion, isolation, small effective population sizes, and inbreeding. We recommend gene flow restoration via genetic rescue to counteract these threats, and provide guidance for this strategy.     

Low host specificity of root‐associated fungi at an Arctic site

In High Arctic ecosystems, plant growth and reproduction are limited by low soil moisture and nutrient availability, low soil and air temperatures, and a short growing season. Mycorrhizal associations facilitate plant nutrient acquisition and water uptake and may therefore be particularly ecologically important in nutrition‐poor and dry environments, such as parts of the Arctic. Similarly, endophytic root associates are thought to play a protective role, increasing plants' stress tolerance, and likely have an important ecosystem function. Despite the importance of these root‐associated fungi, little is known about their host specificity in the Arctic. We investigated the host specificity of root‐associated fungi in the common, widely distributed arctic plant species Bistorta vivipara, Salix polaris and Dryas octopetala in the High Arctic archipelago Svalbard. High‐throughput sequencing of the internal transcribed spacer 1 (ITS1) amplified from whole root systems generated no evidence of host specificity and no spatial autocorrelation within two 3 m × 3 m sample plots. The lack of spatial structure at small spatial scales indicates that Common Mycelial Networks (CMNs) are rare in marginal arctic environments. Moreover, no significant differences in fungal OTU richness were observed across the three plant species, although their root system characteristics (size, biomass) differed considerably. Reasons for lack of host specificity could be that association with generalist fungi may allow arctic plants to more rapidly and easily colonize newly available habitats, and it may be favourable to establish symbiotic relationships with fungi possessing different physiological attributes.     

Whole‐genome resequencing uncovers molecular signatures of natural and sexual selection in wild bighorn sheep

The identification of genes influencing fitness is central to our understanding of the genetic basis of adaptation and how it shapes phenotypic variation in wild populations. Here, we used whole‐genome resequencing of wild Rocky Mountain bighorn sheep (Ovis canadensis) to >50‐fold coverage to identify 2.8 million single nucleotide polymorphisms (SNPs) and genomic regions bearing signatures of directional selection (i.e. selective sweeps). A comparison of SNP diversity between the X chromosome and the autosomes indicated that bighorn males had a dramatically reduced long‐term effective population size compared to females. This probably reflects a long history of intense sexual selection mediated by male–male competition for mates. Selective sweep scans based on heterozygosity and nucleotide diversity revealed evidence for a selective sweep shared across multiple populations at RXFP2, a gene that strongly affects horn size in domestic ungulates. The massive horns carried by bighorn rams appear to have evolved in part via strong positive selection at RXFP2. We identified evidence for selection within individual populations at genes affecting early body growth and cellular response to hypoxia; however, these must be interpreted more cautiously as genetic drift is strong within local populations and may have caused false positives. These results represent a rare example of strong genomic signatures of selection identified at genes with known function in wild populations of a nonmodel species. Our results also showcase the value of reference genome assemblies from agricultural or model species for studies of the genomic basis of adaptation in closely related wild taxa.     

Strong dispersal in a parasitoid wasp overwhelms habitat fragmentation and host population dynamics

The population dynamics of a parasite depend on species traits, host dynamics and the environment. Those dynamics are reflected in the genetic structure of the population. Habitat fragmentation has a greater impact on parasites than on their hosts because resource distribution is increasingly fragmented for species at higher trophic levels. This could lead to either more or less genetic structure than the host, depending on the relative dispersal rates of species. We examined the spatial genetic structure of the parasitoid wasp Hyposoter horticola, and how it was influenced by dispersal, host population dynamics and habitat fragmentation. The host, the Glanville fritillary butterfly, lives as a metapopulation in a fragmented landscape in the Åland Islands, Finland. We collected wasps throughout the 50 by 70 km archipelago and determined the genetic diversity, spatial population structure and genetic differentiation using 14 neutral DNA microsatellite loci. We compared the genetic structure of the wasp with that of the host butterfly using published genetic data collected over the shared landscape. Using maternity assignment, we also identified full‐siblings among the sampled parasitoids to estimate the dispersal range of individual females. We found that because the parasitoid is dispersive, it has low genetic structure, is not very sensitive to habitat fragmentation and has less spatial genetic structure than its butterfly host. The wasp is sensitive to regional rather than local host dynamics, and there is a geographic mosaic landscape for antagonistic co‐evolution of host resistance and parasite virulence.     

Cryptic ecology among host generalist Campylobacter jejuni in domestic animals

Homologous recombination between bacterial strains is theoretically capable of preventing the separation of daughter clusters, and producing cohesive clouds of genotypes in sequence space. However, numerous barriers to recombination are known. Barriers may be essential such as adaptive incompatibility, or ecological, which is associated with the opportunities for recombination in the natural habitat. Campylobacter jejuni is a gut colonizer of numerous animal species and a major human enteric pathogen. We demonstrate that the two major generalist lineages of C. jejuni do not show evidence of recombination with each other in nature, despite having a high degree of host niche overlap and recombining extensively with specialist lineages. However, transformation experiments show that the generalist lineages readily recombine with one another in vitro. This suggests ecological rather than essential barriers to recombination, caused by a cryptic niche structure within the hosts.     

Genetic costructure in a meta‐community under threat of habitat fragmentation

Habitat fragmentation increasingly threatens the services provided by natural communities and ecosystem worldwide. An understanding of the eco‐evolutionary processes underlying fragmentation‐compromised communities in natural settings is lacking, yet critical to realistic and sustainable conservation. Through integrating the multivariate genetic, biotic and abiotic facets of a natural community module experiencing various degrees of habitat fragmentation, we provide unique insights into the processes underlying community functioning in real, natural conditions. The focal community module comprises a parasitic butterfly of conservation concern and its two obligatory host species, a plant and an ant. We show that both historical dispersal and ongoing habitat fragmentation shape population genetic diversity of the butterfly Phengaris alcon and its most limited host species (the plant Gentiana pneumonanthe). Genetic structure of each species was strongly driven by geographical structure, altitude and landscape connectivity. Strikingly, however, was the strong degree of genetic costructure among the three species that could not be explained by the spatial variables under study. This finding suggests that factors other than spatial configuration, including co‐evolutionary dynamics and shared dispersal pathways, cause parallel genetic structure among interacting species. While the exact contribution of co‐evolution and shared dispersal routes on the genetic variation within and among communities deserves further attention, our findings demonstrate a considerable degree of genetic parallelism in natural meta‐communities. The significant effect of landscape connectivity on the genetic diversity and structure of the butterfly also suggests that habitat fragmentation may threaten the functioning of the community module on the long run.     

Patterns of reproductive isolation within and between two Lygaeus species characterized by sexual conflicts over mating

Theory suggests that, under some circumstances, sexual conflict over mating can lead to divergent sexually antagonistic coevolution among populations for traits associated with mating, and that this can promote reproductive isolation and hence speciation. However, sexual conflict over mating may also select for traits (e.g. male willingness to mate) that enhance gene flow between populations, limiting population divergence. In the present study, we compare pre‐ and post‐mating isolation within and between two species characterized by male–female conflict over mating rate. We quantify sexual isolation among five populations of the seed bug Lygaeus equestris collected from Italy and Sweden, and two replicates of a population of the sister‐species Lygaeus simulans, also collected from Italy. We find no evidence of reproductive isolation amongst populations of L. equestris, suggesting that sexual conflict over mating has not led to population divergence in relevant mating traits in L. equestris. However, there was strong asymmetric pre‐mating isolation between L. equestris and L. simulans: male L. simulans were able to mate successfully with female L. equestris, whereas male L. equestris were largely unable to mate with female L. simulans. We found little evidence for strong post‐mating isolation between the two species, however, with hybrid F₂ offspring being produced. Our results suggest that sexual conflict over mating has not led to population divergence, and indeed perhaps supports the contrary theoretical prediction that male willingness to mate may retard speciation by promoting gene flow.