Telomeres are chromosome protectors that shorten during cell replication and in stressful conditions. Developing individuals are susceptible to telomere erosion when their growth is fast and resources limited. This is critical because the rate of telomere attrition in early life is linked to health and life span of adults. The metabolic telomere attrition hypothesis (MeTA) suggests that telomere dynamics can respond to biochemical signals conveying information about the organism’s energetic state. Among these signals are glucocorticoids (hormones that promote catabolic processes, potentially impairing costly telomere maintenance) and nucleotides, which activate anabolic pathways though the cellular enzyme target of rapamycin (TOR) preventing telomere attrition. During the energetically demanding growth phase, the regulation of telomeres in response to two contrasting signals—one promoting telomere maintenance and the other inducing attrition—provides an ideal experimental setting to test MeTa. We studied nestlings of a rapidly developing free-living passerine, the great tit (Parus major), that either received glucocorticoids (Cort-chicks), nucleotides (Nuc-chicks), or a combination of both (NucCort-chicks) all compared with controls (Cnt-chicks). Contrary to Cort-chicks, which showed telomere attrition, NucCort-chicks, did not. NucCort-chicks was the only group showing increased gene expression of telo2 (proxy for TOR activation), of mitochondrial enzymes linked to ATP production (atp5f1a-atp5f1b-cox6a1-cox4) and a higher efficiency in aerobically producing ATP. NucCort-chicks had also a higher expression of telomere maintenance genes (trf2) and of enzymatic antioxidant genes (gpx4-sod1). The findings show that nucleotides availability is crucial for preventing telomere erosion during fast growth in stressful environments.

Introduced and invasive species make excellent natural experiments for investigating rapid evolution. Here, we describe the effects of genetic drift and rapid genetic adaptation in pink salmon (Oncorhynchus gorbuscha) that were accidentally introduced to the Great Lakes via a single introduction event 31-generations ago. Using whole-genome resequencing for 134 fish spanning five sample groups across the native and introduced range, we estimate that the progenitor population’s effective population size was 146,886 at the time of introduction, whereas the founding population’s effective population size was just 72—a 2040-fold decrease. As expected with a severe bottleneck, we show reductions in genome-wide measures of genetic diversity, specifically a 37.7% reduction in the number of SNPs and an 8.2% reduction in observed heterozygosity. Despite this decline in genetic diversity, we provide evidence for putative selection at 47 loci across multiple chromosomes in the introduced populations, including missense variants in genes associated with circadian rhythm, immunological response, and maturation, which match expected or known phenotypic changes in the Great Lakes. For one of these genes, we use a species-specific agent-based model to rule out genetic drift and conclude that a strong response to selection occurred in a period gene (per2) that plays a predominant role in determining an organism’s daily clock, matching large day length differences experienced by introduced salmon during important phenological periods. Together, these results inform how populations might evolve rapidly to new environments, even with a small pool of standing genetic variation.

Epigenetic modifications, like DNA methylation, generate phenotypic diversity in fish and ultimately lead to adaptive evolutionary processes. Euryhaline marine species that migrate between salinity contrasted habitats have received little attention regarding the role of salinity on whole-genome DNA methylation. Investigation of salinity-induced DNA methylation in fish will help to better understand the potential role of this process in salinity acclimation. Using whole genome bisulfite sequencing, we compared DNA methylation patterns in European sea bass (Dicentrarchus labrax) juveniles in seawater and after freshwater transfer. We targeted the gill as a crucial organ involved in plastic responses to environmental changes. To investigate the function of DNA methylation in gills, we performed RNAseq and assessed DNA methylome-transcriptome correlations. We showed a negative correlation between gene expression levels and DNA methylation levels in promoters, first introns and exons. A significant effect of salinity on DNA methylation dynamics with an overall DNA hypomethylation in freshwater-transferred fish compared to seawater controls was demonstrated. This suggests a role of DNA methylation changes in salinity acclimation. Genes involved in key functions as metabolism, ion transport and transepithelial permeability (junctional complexes) were differentially methylated and expressed between salinity conditions. Expression of genes involved in mitochondrial metabolism was increased as well as the expression of DNA methyltransferases 3a. This study reveals novel aspects on the link of DNA methylation and gene expression patterns.

Acquisition of new genes often results in the emergence of novel functions and is a key step in lineage-specific adaptation. As the only group of sessile crustaceans, barnacles establish permanent attachment through initial cement secretion at the larval phase followed by continuous cement secretion in juveniles and adults. However, the origins and evolution of barnacle larval and adult cement proteins remain poorly understood. By performing microdissection of larval cement glands, transcriptome and shotgun proteomics and immunohistochemistry validation, we identified 30 larval and 27 adult cement proteins of the epibiotic turtle barnacle Chelonibia testudinaria, of which the majority are stage- and barnacle-specific. While only two proteins, SIPC and CP100K, were expressed in both larvae and adults, detection of protease inhibitors and the cross-linking enzyme lysyl oxidase paralogs in larvae and adult cement suggested functional convergence. Other barnacle specific cement proteins such as CP100k and CP52k likely share a common origin dating back at least to the divergent of Rhizocephala and Thoracica. Different CP52k paralogs could be detected in larval and adult cement, suggesting stage-specific cement proteins may arise from duplication followed by changes in expression timing of the duplicates. Interestingly, the biochemical properties of larval- and adult-specific CP52k paralogs exhibited remarkable differences, reflecting the composition of cement in different life stages of turtle barnacle might be chemically different. We conclude that de novo gene formation and duplicate neofunctionalization are pivotal to the evolution of lineage-specific cement toolkits in barnacles, which may explain how barnacles can inhabit diverse marine substrata.

Holocentric organisms, unlike typical monocentric organisms, have kinetochore activity distributed along almost the whole length of the chromosome. Because of this, chromosomal rearrangements through fission and fusion are more likely to become fixed in holocentric species, which may account for their extraordinary rates of chromosome evolution. Genome synteny has been reported to be conserved in animals with holocentric chromosomes despite high rates of chromosome rearrangements. Comparing genomes of Carex species and a genome of a distantly related Cyperaceae we have characterised conserved vs. rearranged genome regions across pairs of species that range in time since divergence between 2 and 50 million years. We have compared a C. scoparia genome with a linkage map of the same species to study rearrangements at a population level and suppression of recombination patterns. We found a surprisingly conserved genome synteny even between very distantly related species and extraordinarily high rates of chromosome evolution in genus Carex. Comparing the distribution of repetitive DNA and gene density between conserved and rearranged genomic regions, we found repetitive DNA to be related to holocentromeres and as well as rearranged regions of the genome. This evidence of extremely conserved synteny in sedges and the massive events of chromosome fission and fusion found across the evolution of genus Carex suggests the presence of common genomic hotspots of chromosome evolution related to repetitive DNA.

Male reproductive competition can select for condition-dependent, conspicuous traits that signal some aspect of fighting ability and facilitate assessment of potential rivals. However, the underlying mechanisms that link the signal to a male’s current condition are difficult to investigate in wild populations, often requiring invasive experimental manipulation. Here, we use digital photographs and chest skin samples to investigate mechanisms of a visual signal used in male competition in a wild primate, the red chest patch in geladas (Theropithecus gelada). We analyzed photographs collected during natural (n=144) and anesthetized conditions (n=38) to understand variability in male and female chest redness, and we used chest skin biopsies (n=38) to explore sex differences in gene expression. Male and female geladas showed similar average redness, but males exhibited a wider within-individual range in redness under natural conditions. These sex differences were reflected at the molecular level, with 10.5% of genes exhibiting significant sex differences in expression. Subadult males exhibited intermediate gene expression patters between adult males and females, pointing to mechanisms underlying the development of the red chest patch. We found that genes more highly expressed in males were associated with blood vessel development and maintenance but not with androgen or estrogen activity. Together, our results suggest male gelada redness variability is driven by increased blood vessel branching in the chest skin, providing a potential link between male chest redness and current condition as increased blood circulation to exposed skin could lead to heat loss in the cold, high-altitude environment of geladas.

Adaptive phenotypes are shaped by a combination of genetic and environmental forces, and while the literature is rich with studies focusing on either genetics or environment contributions, those that consider both are rare. Here we utilize the cichlid oral jaw apparatus to fill this knowledge-gap. First, we employed RNA-seq in bony and ligamentous tissues important for jaw opening to identify differentially expressed genes between species and across foraging environments. Our foraging treatments were designed to force animals to employ either suction or biting/scraping, which broadly mimic pelagic or benthic modes of feeding. We found a large number of differentially expressed genes between species, and while we identified relatively few differences between environments, species differences were far more pronounced when reared in the pelagic versus benthic environment. Further, these data carried the signature of genetic assimilation, and implicated cell cycle regulation in shaping the jaw across species and environments. Next, we repeated the foraging experiment and performed ATAC-seq procedures on nuclei harvested from the same tissues. Cross-referencing results from both analyses revealed subsets of genes that were both differentially expressed and differentially accessible in either the pelagic (n=15) or the benthic environment (n=11), as well as loci where differences were robust to foraging environment (n=13). All in all, these data provide novel insights into the epigenetic, genetic, and cellular bases of species- and environment-specific bone shapes, as well as the evolution of phenotypic plasticity in this iconic model system.

During the Late Pleistocene, major parts of North America were periodically covered by ice sheets. However, there are still open questions about whether ice-free refugia were present in the Alexander Archipelago along the Southeast (SE) Alaska coast during the Last Glacial Maximum (LGM). Numerous subfossils have been recovered from caves in SE Alaska, including American black (Ursus americanus) and brown (U. arctos) bears, which today are found in the Alexander Archipelago but are genetically distinct from mainland bear populations. Hence, these bear species offer an ideal system to investigate long-term occupation, potential refugial survival, and lineage turnover. Here we present genetic analyses based on 99 new complete mitochondrial genomes from ancient and modern brown and black bears spanning the last ~45,000 years. Black bears form two SE Alaskan subclades that diverged >100,00 years ago, one preglacial and one postglacial. All postglacial ancient brown bears are closely related to modern bears in the archipelago, while a single preglacial brown bear is found in a distantly related clade. A hiatus in the bear subfossil record around the LGM and the deep split of their pre- and post-glacial subclades fail to support a hypothesis of continuous occupancy in SE Alaska throughout the LGM for either species. Our results are consistent with an absence of refugia along the SE Alaska coast, but indicate that vegetation quickly expanded after deglaciation, allowing bears to recolonize the area after a short-lived LGM peak.

The processes governing soil bacteria biogeography are still not fully understood. It remains unknown how the importance of environmental filtering and dispersal differs between bacterial taxonomic and functional biogeography, and whether their importance is scale-dependent. We sampled soils at 195 plots across the Tibet plateau, with distances among plots ranging from 20 m to 1 550 km. Taxonomic composition of bacterial community was characterized by 16S amplicon sequencing, and functional community composition by qPCR targeting 9 functional groups involved in N dynamics. Twelve climatic and soil characteristics were also measured. Both taxonomic and functional dissimilarities were more related to environmental dissimilarity than geographic distance. Taxonomic dissimilarity was mostly explained by soil pH and organic matter, while functional dissimilarity was mostly linked to moisture, temperature and N, P and C availabilities. The roles of environmental filtering and dispersal were, however, scale-dependent and varied between taxonomic and functional dissimilarities, with distance affecting taxonomic dissimilarity over short distances (<~300 km) and functional dissimilarity over long distances (>~600 km). The importance of different environmental predictors varied across scales more for functional than taxonomic dissimilarity. Our results demonstrate how biodiversity dimension (taxonomic versus functional) and spatial scale strongly influence the conclusions derived from bacterial biogeography studies.

In free-living ecological communities, organismal life histories shape interactions with their environment, which ultimately forms the basis of ecological succession. Individual animals in natural populations tend to host diverse parasite species concurrently over their lifetimes. However, the structure and dynamics of mammalian parasite communities have not been contextualized in terms of primary ecological succession, in part because few datasets track occupancy and abundance of multiple parasites in wild hosts starting at birth. Here, we studied community dynamics of twelve subtypes of protozoan microparasites (Theileria spp.) in a herd of African buffalo. We show that Theileria communities followed predictable patterns of succession underpinned by four different parasite life-history strategies. In contrast to many free-living communities, network complexity decreased with host age. Examining parasite communities through the lens of succession may better inform the effect of complex within host eco-evolutionary dynamics on infection outcomes, including parasite co-existence through the lifetime of the host.

Host-parasite interactions can cause strong demographic fluctuations accompanied by selective sweeps of resistance/infectivity alleles. Both demographic bottlenecks and frequent sweeps are expected to reduce the amount of segregating genetic variation and therefore might constrain adaption during coevolution. Recent studies, however, suggest that the interaction of demographic and selective processes is a key component of coevolutionary dynamics and may rather positively affect levels of genetic diversity available for adaptation. Here, we provide direct experimental testing of this hypothesis by disentangling the effect of demography, selection, and of their interaction in an experimental host-parasite system. We grew 12 populations of unicellular algae (Chlorella variabilis) that experienced either growth followed by constant population sizes (3 populations), demographic fluctuations (3 populations), selection induced by exposure to a virus (3 populations), or demographic fluctuations together with virus-induced selection (3 populations). After 50 days, we conducted whole-genome sequencing of each algal population. We observed more genetic diversity in populations that jointly experienced selection and demographic fluctuations than in populations where these processes were experimentally separated. In addition, in those 3 populations that jointly experienced selection and demographic fluctuations, experimentally measured diversity exceeds expected values of diversity that account for the cultures’ population sizes. Our results suggest that eco-evolutionary feedbacks can positively affect genetic diversity and provide the necessary empirical measures to guide further improvements of theoretical models of adaptation during host-parasite coevolution.

Traditional agrosystems, where humans, crops and microbes have coevolved over long periods, can serve as models to understand the eco-evolutionary determinants of disease dynamics and help the engineering of durably resistant agrosystems. Here, we investigated the genetic and phenotypic relationship between rice (Oryza sativa) landraces and their rice blast pathogen (Pyricularia oryzae) in the traditional Yuanyang terraces of flooded rice paddies in China, where rice landraces have been grown and bred over centuries without significant disease outbreaks. Analyses of genetic subdivision revealed that indica rice plants clustered according to landrace names. Three new diverse lineages of rice blast specific to the Yuanyang terraces coexisted with lineages previously detected at the worldwide scale. Population subdivision in the pathogen population did not mirror pattern of population subdivision in the host. Measuring the pathogenicity of rice blast isolates on landraces revealed generalist life histories. Our results suggest that the implementation of disease control strategies based on the emergence or maintenance of a generalist lifestyle in pathogens may sustainably reduce the burden of disease in crops.

The blacklegged tick, Ixodes scapularis, is a vector of Borrelia burgdorferi sensu stricto (s.s.), the causative agent of Lyme disease, part of a slow-moving epidemic of Lyme borreliosis spreading across the northern hemisphere. There are well-known geographic differences in the vectorial capacity of these ticks associated with genetic variation. Despite the need for detailed genetic information in this disease system, previous phylogeographic studies of these ticks have been restricted to relatively few populations or genetic loci. Here we present the most comprehensive phylogeographic study of I. scapularis conducted by using 3RAD and surveying 353 ticks from 33 counties throughout the range of I. scapularis. We found limited genetic variation among populations from the Northeast and Upper Midwest, where Lyme disease is most common, and higher genetic variation among populations from the South. We identify four genetic clusters of I. scapularis that are consistent with four major geographic regions, plus a distinct Central Florida group. In regions where Lyme disease is increasing in frequency, the I. scapularis populations genetically group with ticks from historically highly Lyme-endemic regions. Finally, we identify ten variable DNA sites that contribute the most to population differentiation. These variable sites cluster on one of the chromosome-scale scaffolds for I. scapularis and are within identified genes. Our findings illuminate the need for additional research to identify loci causing variation in the vectorial capacity of I. scapularis and where additional tick sampling would be most valuable to further understand disease trends caused by pathogens transmitted by I. scapularis.

Anthropogenic changes have altered the historical distributions of many North American taxa. As environments shift, ecological and evolutionary processes can combine in complex ways to either stimulate or inhibit range expansion. Here we examine the role of evolution in a rapid range expansion whose ecological context has been well-documented, Anna’s Hummingbird (Calypte anna). Previous work suggests that the C. anna range expansion is the result of an ecological release facilitated by human-mediated environmental changes, where access to new food sources have allowed further filling of the abiotic niche. We examine the role of gene flow and adaptation during range expansion from their native California north into Canada and east into New Mexico and Texas, USA. Using low coverage whole genome sequencing we found high genetic diversity, low divergence, and little evidence of selection on the northern and eastern expansion fronts. Additionally, there are few (if any) limits to gene flow across the native and expanded range. The lack of selective signals between core and expanded ranges could reflect i) an absence of novel selection pressure in the extended range (supporting the ecological release hypothesis), ii) swamping of adaptive variation due to high gene flow, or iii) limitations of genome scans for detecting small shifts in allele frequencies across many loci. Nevertheless, our results provide an example where strong selection is not apparent during a rapid, contemporary range shift.

How changes in biodiversity affect disease, particularly in the face of large-scale land-use change, is a contentious topic in disease ecology that has implications for public health and conservation policy. The ‘dilution effect’ hypothesis argues that declines in biodiversity are associated with increased disease risk, but this can be challenging to demonstrate because many pathogens have complex life cycles such that changes to the species composition and abundance of hosts can influence the density and infection prevalence of vectors via multiple mechanisms. Key to addressing this debate is a quantification of interactions between hosts, vectors, and pathogens. In their recent study published in Molecular Ecology, Kocher et al. (2022) captured thousands of sandflies, some species of which are vectors for the Leishmania protozoan that causes Leishmaniasis, across a human footprint gradient in French Guiana (Fig. 1). By implementing DNA metabarcoding of vectors combined with an innovative modeling approach, they effectively quantified the nuanced relationships between changes in land-use, mammalian host diversity, vector abundance, and parasite prevalence. In support of the dilution effect hypothesis, Kocher et al. found that sites with higher mammal diversity were associated with lower relative abundance of reservoir hosts and higher Leishmania infection prevalence in sandflies. However, while infection prevalence was lower when mammal diversity was high, the density of sandfly vectors was higher, which resulted in a weak overall effect of mammal diversity on the density of infected vectors, the most important indicator of Leishmania transmission risk.

Adaptation enables natural populations to survive in a changing environment. Understanding the mechanics of adaptation is therefore crucial for learning about the evolution and ecology of natural populations, and for better conservation and management of natural resources such as fish stocks. In this review we focus on the impact of random sweepstakes on selection in highly fecund populations. In random sweepstakes the distribution of individual recruitment success is highly skewed, resulting in a huge variance in the number of offspring contributed by the individuals present in any given generation. We also describe selective sweepstakes which are well approximated by recurrent selective sweeps of strongly beneficial allelic types arising by mutation. We demonstrate that both types of sweepstakes reproduction may facilitate rapid adaptation. Finally, we review an important case study in which a model of recurrent selective sweeps is shown to essentially explain population genomic data of the highly fecund Atlantic cod, with broad implications for studying the evolution and ecology of highly fecund populations across domains of life.

An exhaustive assessment of biodiversity is a major challenge of ecological research, and molecular approaches such as the metabarcoding of environmental DNA are boosting our ability to perform biodiversity inventories. Are we actually able to assess the whole community, to unravel the intricate interactions between organisms and the impacts of global changes on the different trophic levels? The majority of metabarcoding papers published in the last years used just one or two markers and analyzed a limited number of taxonomic groups. Nevertheless, approaches are emerging that might allow “all-taxa biological inventories”. Exhaustive biodiversity assessments can be attempted by combining a large number of specific primers, by exploiting the power of universal primers, or by combining specific and universal primers to obtain good information on key taxa while limiting the overlooked biodiversity. Multiplexes of primers and shotgun sequencing may provide a better coverage of biodiversity compared to standard metabarcoding, but still require major methodological advances. We identify the strengths and limitations of different approaches, and suggest new development lines that might improve broad scale biodiversity analyses in the near future.

Wolbachia are among the most prevalent and widespread endosymbiotic bacteria on earth. Wolbachia’s success in infecting an enormous number of arthropod species is attributed to two features: the range of phenotypes they induce in their hosts, and their ability to switch between host species. Whilst much progress has been made in elucidating their induced phenotypes, our understanding of Wolbachia host shifting is still very limited: we lack answers to even fundamental questions concerning Wolbachia’s routes of transfer and the importance of factors influencing host shifts. Here, we investigate the diversity and host-shift patterns of Wolbachia in scale insects, a group of arthropods with intimate associations with other insects that make them well-suited to studying host shifts. Using Illumina multi-target amplicon sequencing of Wolbachia-infected scale insects and their direct associates we determined the identity of all Wolbachia strains. We then fitted a Generalised Additive Mixed Model (GAMM) to our data to estimate the influence of host phylogeny and the geographic distribution on Wolbachia strain sharing among scale insect species. The model predicts no significant contribution of host geography but strong effects of host phylogeny, with high rates of Wolbachia sharing among closely related species and a sudden drop-off in sharing with increasing phylogenetic distance. We also detected the same Wolbachia strain in scale insects and several intimately associated species (ants, wasps, beetles, and flies). This indicates putative host shifts and potential routes of transfers via these associates and highlights the importance of ecological connectivity in Wolbachia host-shifting.