Understanding how organisms adapt to changing environments is a core focus of research in evolutionary biology. One common mechanism is adaptive introgression, which has received increasing attention as a potential route to rapid adaptation in populations struggling in the face of ecological change, particularly global climate change. However, hybridization can also result in deleterious genetic interactions that may limit the benefits of adaptive introgression. Here, we used a combination of genome-wide quantitative trait locus mapping and differential gene expression analyses between the swordtail fish species Xiphophorus malinche and X. birchmanni to study the consequences of hybridization on thermotolerance. While these two species are adapted to different thermal environments, we document a complicated architecture of thermotolerance in hybrids. We identify a region of the genome that contributes to reduced thermotolerance in individuals heterozygous for X. malinche and X. birchmanni ancestry, as well as widespread misexpression in hybrids of genes that respond to thermal stress in the parental species, particularly in the circadian clock pathway. We also show that a previously mapped hybrid incompatibility between X. malinche and X. birchmanni contributes to reduced thermotolerance in hybrids. Together, our results highlight the challenges of understanding the impact of hybridization on complex ecological traits and its potential impact on adaptive introgression.

The domestication process in long-lived plant perennials differs dramatically from that of annuals, with a huge amount of genetic exchange between crop and wild populations. Though apple is a major fruit crop grown worldwide, the contribution of wild apple species to the genetic makeup of the cultivated apple genome remains a topic of intense study. We used population genomics approaches to investigate the contributions of several wild apple species to European and Chinese rootstock and dessert genomes, with a focus on the extent of wild-crop gene flow. Population genetic structure inferences revealed that the East Asian wild apples, Malus baccata (L.) Borkh and M. hupehensis (Pamp.), form a single panmictic group, and that the European dessert and rootstock apples form a specific gene pool whereas the Chinese dessert and rootstock apples were a mixture of three wild gene pools, suggesting different evolutionary histories of European and Chinese apple varieties. Coalescent-based inferences and gene flow estimates indicated that M. baccata - M. hupehensis contributed to the genome of both European and Chinese cultivated apples through wild-to-crop introgressions, and not as an initial contributor as previously supposed. We also confirmed the contribution through wild-to-crop introgressions of Malus sylvestris Mill. to the cultivated apple genome. Apple tree domestication is therefore one example in woody perennials that involved gene flow from several wild species from multiple geographical areas. This study provides an example of a complex protracted process of domestication in long-lived plant perennials, and is a starting point for apple breeding programmes.

Understanding the ecological rules structuring the organization of species interactions is a prerequisite to predicting how ecosystems respond to environmental changes. While the ecological determinants of single networks have been documented, it remains unclear whether network ecological rules are conserved along spatial and environmental gradients. To address this gap, we reconstructed 48 plant-herbivore interaction networks along six elevation gradients in the Central European Alps in Switzerland, using DNA metabarcoding on orthoptera faeces. We developed hypotheses on the ecological mechanisms expected to structure interaction networks, based on plant phylogeny, plant abundance, leaf toughness, leaf nitrogen content and plant metabolomics. We show that plant phylogenetic relationships and species abundance have the greatest explanatory power regarding the structure of the ecological networks. Moreover, we found that leaf nitrogen content is a key determinant of interactions in warmer environments, while phenolic compounds and tannins are more important in colder environments, suggesting that determinants of species interactions can shift along environmental gradients. With this work, we propose an approach to study the mechanisms that structure the way species interact with each other between bioregions and ecosystems.

Plasmids provide an efficient vehicle for gene sharing among bacterial populations, playing a key role in bacterial evolution. Network approaches are particularly suitable to represent multipartite relationships and are useful tools to characterize plasmid-mediated gene sharing events. The bacterial family Lysobacteraceae includes plant commensal, plant pathogenic and opportunistic human pathogens for which plasmid-mediated adaptation has been reported. We searched for homologues of plasmid gene sequences from this family in the entire diversity of available bacterial genome sequences and built a network of plasmid gene sharing from the results. While plasmid genes are openly shared between the bacteria of the family Lysobacteraceae, taxonomy strongly defined the boundaries of these exchanges, which only barely reached other families. Most inferred plasmid gene sharing events involved a few genes only, and evidence of full plasmid transfers were restricted to taxonomically closely related taxa. We detected multiple plasmid-chromosome gene transfers, including the known sharing of a heavy metal resistance transposon. In the network, bacterial lifestyles shaped substructures of isolates colonizing specific ecological niches and harbouring specific types of resistance genes. Genes associated with pathogenicity or antibiotic and metal resistance were among those that most importantly structured the network, highlighting the imprints of human-mediated selective pressure on pathogenic populations. A massive sequencing effort on environmental Lysobacteraceae is therefore required to refine our understanding of how this reservoir fuels the emergence and the spread of genes among this family and its potential impact on plant, animal and human health.

Ecologically distinct species may respond to landscape changes in different ways. In addition to basic ecological data, the extent of the geographic range has been successfully used as an indicator of species sensitivity to anthropogenic landscapes, with widespread species usually found to be less sensitive compared to range-restricted species. In this study, we investigate connectivity patterns of two closely related but ecologically distinct newt species - the range-restricted, Lissotriton montandoni and the widespread, L. vulgaris - using genomic data, a highly replicated setting (six geographic regions per species), and tools from landscape genetics. Our results show the importance of forest for connectivity in both species, but at the same time suggest differential use of forested habitat, with L. montandoni and L. vulgaris showing the highest connectivity at forest-core and forest-edges, respectively. Anthropogenic landscapes (i.e., higher crop- or urban-cover) increased resistance in both species, but the effect was one to three orders of magnitude stronger in L. montandoni than in L. vulgaris. This result is consistent with a view of L. vulgaris as an ecological generalist. Even so, currently, the negative impact of anthropogenic landscapes is mainly seen in connectivity among L. vulgaris populations, which show significantly stronger isolation and lower effective sizes relative to L. montandoni. Overall, this study emphasizes how habitat destruction is compromising genetic connectivity not only in endemic, range-restricted species of conservation concern but also in widespread generalist species, despite their comparatively lower sensitivity to anthropogenic landscape changes.

High-throughput sequencing has substantially improved our understanding of fungal diversity. However, the short read (<500 bp) length of current second-generation sequencing approaches provides limited taxonomic and phylogenetic resolution for species discrimination. Longer sequences containing more information are highly desired to provide greater taxonomic resolution. Here, we amplified full-length rRNA operons (~5.5 kb) and established a corresponding fungal rRNA operon database for ONT sequences (FRODO), which contains ONT sequences representing eight phyla, 41 classes, 109 orders, 256 families, 524 genera and 1116 species. We also benchmarked the optimal method for sequence classification and determined that the RDP classifier based on our FRODO database was capable of improving the classification of ONT reads, with an average of 98%-99% reads correctly classified at the genus or species level. We investigated the applicability of our approach in three representative mycobiomes, namely, the soil, marine and human gut mycobiomes, and found that the gut contains the largest number of unknown species (over 90%), followed by the marine (42%) and soil (33.8%) mycobiomes. We also observed a distinct difference in the composition of the marine and soil mycobiomes, with the highest richness and diversity detected in soils. Overall, our study provides a systematic approach for mycobiome studies and revealed that the previous methods might have underestimated the diversity of mycobiome species. Future application of this method will lead to a better understanding of the taxonomic and functional diversity of fungi in environmental and health-related mycobiomes.

Long-read sequencing is driving a new reality for genome science in which highly contiguous assemblies can be produced efficiently with modest resources. Genome assemblies from long-read sequences are particularly exciting for understanding the evolution of complex genomic regions that are often difficult to assemble. In this study, we utilized long-read sequencing data to generate a high-quality genome assembly for an Antarctic eelpout, Ophthalmolycus amberensis, the first for the globally distributed family Zoarcidae. We used this assembly to understand how O. amberensis has adapted to the harsh Southern Ocean and compared it to another group of Antarctic fishes: the notothenioids. We showed that selection has largely acted on different targets in eelpouts relative to notothenioids. However, we did find some overlap; in both groups, genes involved in membrane structure, thermal tolerance and vision have evidence of positive selection. We found evidence for historical shifts of transposable element activity in O. amberensis and other polar fishes, perhaps reflecting a response to environmental change. We were specifically interested in the evolution of two complex genomic loci known to underlie key adaptations to polar seas: haemoglobin and antifreeze proteins (AFPs). We observed unique evolution of the haemoglobin MN cluster in eelpouts and related fishes in the suborder Zoarcoidei relative to other Perciformes. For AFPs, we identified the first species in the suborder with no evidence of afpIII sequences (Cebidichthys violaceus) in the genomic region where they are found in all other Zoarcoidei, potentially reflecting a lineage-specific loss of this cluster. Beyond polar fishes, our results highlight the power of long-read sequencing to understand genome evolution.

The seafloor contains complex ecosystems where habitat heterogeneity influences biodiversity. Natural biological and geological features including vents, seeps and reefs create habitats that select for distinct populations of micro- and macrofauna. While largely studied for macrobiological diversity, built habitats may also select distinct microbiomes. Built habitat density on the seafloor is increasing with ocean sprawl expanding in the continental shelf and slope, potentially having widespread effects on benthic communities. This study addresses one type of built habitat, shipwrecks, on microbiomes in surrounding sediment. Using deep-sea sediment samples (762 total) from the Gulf of Mexico, we report elevated diversity and a predictable core microbiome around nine shipwrecks. We show the sphere of influence of built habitats extends up to 300 m onto the seafloor. Supervised learning made predictions of sample proximity to structures based on frequency of taxa. Strongest predictions occurred in sediments nearest and furthest from sites for archaea and mid-transect for bacteria. The response of archaea to built habitats was consistent across sites, while bacteria showed greater between site variability. The archaeal core shipwreck microbiome was enriched in taxa (e.g., Bathyarchaeia, Lokiarchaeia, Thermoplasmata) not present in the surrounding seafloor. Shipwrecks shaped microbiomes in expected ways, providing insight on how built habitats impact microbiome biodiversity in the Anthropocene.

Potato cyst nematodes (PCNs), an umbrella term used for two species, Globodera pallida and G. rostochiensis, belong worldwide to the most harmful pathogens of potato. Pathotype-specific host plant resistances are essential for PCN control. However, the poor delineation of G. pallida pathotypes has hampered the efficient use of available host plant resistances. Long-read sequencing technology allowed us to generate a new reference genome of G. pallida population D383 and, as compared to the current reference, the new genome assembly is 42 times less fragmented. For comparison of diversification patterns of six effector families between G. pallida and G. rostochiensis, an additional reference genome was generated for an outgroup, the beet cyst nematode Heterodera schachtii (IRS population). Large evolutionary contrasts in effector family topologies were observed. While VAPs (venom allergen-like proteins) diversified before the split between the three cyst nematode species, the families GLAND5 and GLAND13 only expanded in PCNs after their separation from the genus Heterodera. Although DNA motifs in the promoter regions thought to be involved in the orchestration of effector expression ("DOG boxes") were present in all three cyst nematode species, their presence is not a necessity for dorsal gland-produced effectors. Notably, DOG box dosage was only loosely correlated with the expression level of individual effector variants. Comparison of the G. pallida genome with those of two other cyst nematodes underlined the fundamental differences in evolutionary history between effector families. Resequencing of PCN populations with different virulence characteristics will allow for the linking of these characteristics to the composition of the effector repertoire as well as for the mapping of PCN diversification patterns resulting from extreme anthropogenic range expansion.