Monitoring bloom-forming Aphanizomenon using environmental DNA metabarcoding: Method development, validation, and field application.

Harmful algal blooms (HABs) have emerged as a critical global environmental and ecological concern. Timely and accurate monitoring of the prevalent bloom-forming genera is crucial for HAB management. Conventional microscope-based methods are time-consuming, labor-intensive, and specialized expertise-dependent, often making them impractical for large-scale surveillance. Molecular methods, such as metabarcoding, provide efficient technical solutions; however, the lack of competent PCR primers and further field validation present obstacles to their wide use. Here, we successfully developed Aphanizomenon-specific primers and validated the application of environmental DNA (eDNA) metabarcoding for field-based monitoring of Aphanizomenon in 37 sites across lentic and lotic freshwater ecosystems in Beijing. The sensitivity and specificity tests of newly developed primers demonstrated high performance - comprehensive recovery of biodiversity in Aphanizomenon communities and high ratios (>95%) of Aphanizomenon sequences in datasets. We observed significant correlations between the sequence abundance derived from eDNA metabarcoding and the total cell density determined through microscopic identification across all the sampling sites, both in the spring (r = 0.8086, p < 0.0001) and summer (r = 0.7902, p < 0.0001), thus validating the utility of eDNA metabarcoding based on the newly developed primers for monitoring in the field. Further, we identified key environmental variables that were primary drivers responsible for the spatiotemporal distribution of Aphanizomenon abundance. These variables included temperature, total nitrogen, and dissolved oxygen in lentic ecosystems, and total phosphorus in lotic ecosystems. The method developed and validated here offers an accurate, efficient, and high-throughput tool for the monitoring of Aphanizomenon blooms in freshwater ecosystems.

Environmental DNA-Based Identification of Non-Native Fish in Beijing: Diversity, Geographical Distribution, and Interactions with Native Taxa.

Rapid urbanization and its associated human activities have facilitated the colonization and spread of non-native species, rendering urban ecosystems, particularly in megacities such as Beijing, highly susceptible to biological invasions. This study employed environmental DNA (eDNA) metabarcoding to evaluate the biodiversity and geographical distribution of non-native fish, as well as their interactions with native fish species, across three river basins in Beijing pertaining to the Daqing River, the North Canal, and the Ji Canal. Across all the 67 sampling sites, we identified 60 fish taxa, representing 11 orders, 23 families, and 40 genera, with an average of 33.0 taxa per site. Of these, 40 taxa were native, accounting for only 47.1% of the historically recorded native fish species. Additionally, we detected 20 non-native fish taxa, spanning 11 orders, 13 families, and 17 genera. Native fish exhibited geographical homogenization across the basins, while non-native taxa displayed varied geographical distributions. Non-metric multidimensional scaling (NMDS) and analysis of similarities (ANOSIM) revealed no significant variation in the non-native communities across the river basins. Although most of the non-native taxa were widespread, some were restricted to specific sites or basins. The North Canal exhibited significantly lower non-native biodiversity compared with the Ji Canal across all alpha diversity indices. Simple linear regression analyses indicated positive correlations between the number of taxa and species richness for both native and non-native taxa. Interestingly, species co-occurrence analyses revealed predominantly positive interactions among both native and non-native species pairs, with only two negative relationships involving one native and two non-native taxa. This study provides insights into the biodiversity and geographical distribution of non-native fish in Beijing and establishes a baseline for future biomonitoring and conservation efforts. The findings underscore the need for further investigation into the mechanisms and dynamics of biological invasions within urban environments in Beijing.

Revealing assembly mechanisms of algal communities in aquatic microniches: Shifts in diversity patterns, microbial interactions and stability along nutrient gradients.

Algal blooms threaten water quality and ecosystem stability in aquatic habitats globally, yet dynamics regulating phytoplankton community assembly, the basis of blooms, remain poorly characterized in small water bodies. Here, we employed high-throughput sequencing to analyze drivers structuring phytoplankton across a trophic gradient of 10 small water bodies over 12 consecutive months. Cyanobacteria and Chlorophyta were identified as potential seed banks priming blooms. Temporal variation in community composition was muted in nutrient-limited waters given Cyanobacteria dominance. Environmental factors and interspecific relationships jointly governed temporal phytoplankton dynamics. Phytoplankton, exhibiting greater sensitivity, responded more rapidly than bacterioplankton to environmental and biological fluctuations. This research provides a robust bench mark characterizing planktonic successional trajectories across small water bodies varying in trophic status. Results reinforce ecological mechanisms underpinning biological control strategies to mitigate algal proliferation and inform water quality management of these ubiquitous aquatic ecosystems.

Metropolitan pressures: Significant biodiversity declines and strong filtering of functional traits in fish assemblages.

Accelerating global urbanization is leading to drastic losses and restructuring of biodiversity. Although it is crucial to understand urban impacts on biodiversity to develop mitigation strategies, there is a dearth of knowledge on the functional structure of fish assemblages spanning the entire city-scale spectrum of urbanization intensity. Here, using environmental DNA sampled from 109 water sites in Beijing, we investigated the taxonomic and functional diversity patterns of fish assemblages across the city and uncovered community-, trait-, and species-level responses to various environmental stressors. By ranking sampling sites into three disturbance levels according to water physiochemical and landcover conditions, we found that both native and non-native fish taxonomic and functional α-diversity decreased significantly with elevating disturbance, as strong disturbance led to the disappearance of many species. However, the quantitative taxonomic and functional ß-diversity components of native and non-native fish showed distinct patterns; assemblage turnover dominated native fish ß-diversity and decreased with increasing disturbance, whereas species/trait richness differences dominated non-native fish ß-diversity and increased with disturbance intensity particularly in lotic waters. RLQ and fourth-corner analyses revealed that fish size, fecundity, diet, and reproductive behaviors were significantly correlated with water quality, with pollution-tolerant, larger-sized native and omnivorous non-native fishes being urban winners, which indicates strong trait-dependent environmental filtering. Potential ecological indicator species were identified based on the sensitivity of fish responses to pollution loads; these were mostly small native species, and many have bivalve-dependent reproduction. Our results demonstrate that, along with native fish assemblage simplification and homogenization, urban stressors exert profound impacts on community trait composition, highlighting the need to consider both biodiversity loss and functional reorganization in combating disturbance of aquatic ecosystems under global urbanization. Furthermore, correlations between cropland cover and water nutrient level suggested that the management of agricultural runoff might be critically important for safeguarding urban water quality.

Context-dependent antioxidant defense system (ADS)-based stress memory in response to recurrent environmental challenges in congeneric invasive species.

Marine ecosystems are facing escalating environmental fluctuations owing to climate change and human activities, imposing pressures on marine species. To withstand recurring environmental challenges, marine organisms, especially benthic species lacking behavioral choices to select optimal habitats, have to utilize well-established strategies such as the antioxidant defense system (ADS) to ensure their survival. Therefore, understanding of the mechanisms governing the ADS-based response is essential for gaining insights into adaptive strategies for managing environmental challenges. Here we conducted a comparative analysis of the physiological and transcriptional responses based on the ADS during two rounds of 'hypersalinity-recovery' challenges in two model congeneric invasive ascidians, Ciona robusta and C. savignyi. Our results demonstrated that C. savignyi exhibited higher tolerance and resistance to salinity stresses at the physiological level, while C. robusta demonstrated heightened responses at the transcriptional level. We observed distinct transcriptional responses, particularly in the utilization of two superoxide dismutase (SOD) isoforms. Both Ciona species developed physiological stress memory with elevated total SOD (T-SOD) and glutathione (GSH) responses, while only C. robusta demonstrated transcriptional stress memory. The regulatory distinctions within the Nrf2-Keap1 signalling pathway likely explain the formation disparity of transcriptional stress memory between both Ciona species. These findings support the 'context-dependent stress memory hypothesis', emphasizing the emergence of species-specific stress memory at diverse regulatory levels in response to recurrent environmental challenges. Our results enhance our understanding of the mechanisms of environmental challenge management in marine species, particularly those related to the ADS. Supplementary Information: The online version contains supplementary material available at 10.1007/s42995-024-00228-y.

Spatial distribution of benthic taxonomic and functional diversity in the Yellow River Basin: From ecological processes to associated determinant factors.

One of the fundamental objectives in ecology is to investigate the ecological processes and associated factors governing the abundance and spatial distribution patterns of biodiversity. However, the reaction of biological communities to environmental degradation remains relatively unknown, even for ecologically crucial communities like macroinvertebrates in aquatic ecosystems. Here, we sampled 117 locations to quantify relative contributions of geographical and environmental factors, including water quality, land use, climate, and hydrological factors, to determine the absolute and relative compositions of macroinvertebrate communities and their spatial distribution in the Yellow River Basin (YRB), the sixth-longest river system on Earth. We assessed relative roles of species sorting and dispersal in determining macroinvertebrate community structure along YRB. Our results demonstrated that alpha and beta diversity indices showed an increase from the up- to low-reaches of YRB. The middle and low-reaches exhibited elevated species diversity and both regions exhibited relatively stable community compositions. The biodiversity of macroinvertebrates was influenced by a combination of geographical factors and environmental variables, with environmental factors predominantly serving as the principal determinants. Results of multiple linear regression and variance decomposition showed that the effect of environmental factors was approximately three times greater than that of spatial factors. These findings provide support for the hypothesis that species sorting, driven by environmental gradients, plays a significant role in shaping the community structure of macroinvertebrates in running water ecosystems at the basin scales. Moreover, the factors contributing to substantial shifts in biodiversity across different segments of YRB indicate that distinct river sections have been influenced by varying stressors, with downstream areas being more susceptible to the impacts of water pollution and urbanization resulting from human activities.

Deciphering protein-mediated underwater adhesion in an invasive biofouling ascidian: Discovery, validation, and functional mechanism of an interfacial protein.

Discovering macromolecules and understanding the associated mechanisms involved in underwater adhesion are essential for both studying the fundamental ecology of benthos in aquatic ecosystems and developing biomimetic adhesive materials in industries. Here, we employed quantitative proteomics to assess protein expression variations during the development of the distinct adhesive structure - stolon in the model fouling ascidian, Ciona robusta. We found 16 adhesive protein candidates with increased expression in the stolon, with ascidian adhesive protein 1 (AAP1) being particularly rich in adhesion-related signal peptides, amino acids, and functional domains. Western blot and immunolocalization analyses confirmed the prominent AAP1 signals in the mantle, tunic, stolon, and adhesive footprints, indicating the interfacial role of this protein. Surface coating and atomic force microscopy experiments verified AAP1's adhesion to diverse materials, likely through the specific electrostatic and hydrophobic amino acid interactions with various substrates. In addition, molecular docking calculations indicated the AAP1's potential for cross-linking via hydrogen bonds and salt bridges among Von Willebrand factor type A domains, enhancing its adhesion capability. Altogether, the newly discovered interfacial protein responsible for permanent underwater adhesion, along with the elucidated adhesion mechanisms, are expected to contribute to the development of biomimetic adhesive materials and anti-fouling strategies. STATEMENT OF SIGNIFICANCE: Discovering macromolecules and studying their associated mechanisms involved in underwater adhesion are essential for understanding the fundamental ecology of benthos in aquatic ecosystems and developing innovative bionic adhesive materials in various industries. Using multidisciplinary analytical methods, we identified an interfacial protein - Ascidian Adhesive Protein 1 (AAP1) from the model marine fouling ascidian, Ciona robusta. The interfacial functions of AAP1 are achieved by electrostatic and hydrophobic interactions, and the Von Willebrand factor type A domain-based cross-linking likely enhances AAP1's interfacial adhesion. The identification and validation of the interfacial functions of AAP1, combined with the elucidation of adhesion mechanisms, present a promising target for the development of biomimetic adhesive materials and the formulation of effective anti-fouling strategies.

An overlooked source of false positives in eDNA-based biodiversity assessment and management.

Biodiversity conservation and management in urban aquatic ecosystems is crucial to human welfare, and environmental DNA (eDNA)-based methods have become popular in biodiversity assessment. Here we report a highly overlooked source of significant false positives for eDNA-based biodiversity assessment in urban aquatic ecosystems supplied with treated wastewater - eDNA pollution originating from treated wastewater represents a noteworthy source of false positives. To investigate whether eDNA pollution is specific to a certain treatment or prevalent across methods employed by wastewater treatment plants, we conducted tests on effluent treated using three different secondary processes, both before and after upgrades to tertiary treatment. We metabarcoded eDNA collected from effluent immediately after full treatment and detected diverse native and non-native, commercial and ornamental fishes (48 taxa) across all treatment processes before and after upgrades. Thus, eDNA pollution occurred irrespective of the treatment processes applied. Release of eDNA pollution into natural aquatic ecosystems could translate into false positives for eDNA-based analysis. We discuss and propose technical solutions to minimize these false positives in environmental nucleic acid-based biodiversity assessments and conservation programs.

Unraveling the ecological threads: How invasive alien plants influence soil carbon dynamics.

Invasive alien plants (IAPs) pose significant threats to native ecosystems and biodiversity worldwide. However, the understanding of their precise impact on soil carbon (C) dynamics in invaded ecosystems remains a crucial area of research. This review comprehensively explores the mechanisms through which IAPs influence soil C pools, fluxes, and C budgets, shedding light on their effects and broader consequences. Key mechanisms identified include changes in litter inputs, rates of organic matter decomposition, alterations in soil microbial communities, and shifts in nutrient cycling, all driving the impact of IAPs on soil C dynamics. These mechanisms affect soil C storage, turnover rates, and ecosystem functioning. Moreover, IAPs tend to increase gross primary productivity and net primary productivity leading to the alterations in fluxes and C budgets. The implications of IAP-induced alterations in soil C dynamics are significant and extend to plant-soil interactions, ecosystem structure, and biodiversity. Additionally, they have profound consequences for C sequestration, potentially impacting climate change mitigation. Restoring native plant communities, promoting soil health, and implementing species-specific management are essential measures to significantly mitigate the impacts of IAPs on soil C dynamics. Overall, understanding and mitigating the effects of IAPs on soil C storage, nutrient cycling, and related processes will contribute to the conservation of native biodiversity and complement global C neutrality efforts.

eDNA-Based Early Detection Illustrates Rapid Spread of the Non-Native Golden Mussel Introduced into Beijing via Water Diversion.

The world's largest water diversion, the South-to-North Water Transfer Project (SNWTP) in China, has created an "invasion highway" to introduce invasive golden mussels (Limnoperna fortunei) from the Yangtze River basin to Beijing. To examine the spread and colonization patterns of this newly introduced invasive species, we conducted comprehensive environmental DNA (eDNA)-based early detection and conventional field surveys across all water bodies in five river basins in Beijing from 2020 to 2023. Our results indicated a rapid spread over the past four years. Among the 130 tested sites, the number of sites with positive signals from eDNA analysis exhibited an annual increase: Commencing with four infested sites identified through field surveys in 2019, eDNA analysis detected an additional 13, 11, and 10 positive sites in 2020, 2021, and 2022, respectively, and a substantial rise comprising an additional 28 sites in 2023. Conventional field surveys detected mussels 1-3 years later than eDNA-based analysis at 16 sites. Across all 16 sites, we detected a low population density ranging from 1 to 30 individuals/m2. These findings collectively indicate that the invasions by golden mussels in Beijing are still in their early stages. To date, golden mussels have successfully colonized four out of the five investigated river basins, including the Jiyun River (22.2% positive sites), North Canal River (59.6% positive sites), Chaobai River (40% positive sites), and Yongding River (63.6% positive sites), with the North Canal River and Yongding River being the most heavily infested. Currently, only the Daqing River basin remains uninfested. Given the significant number of infested sites and the ongoing transport of large new propagules via SNWTP, further rapid spread and colonization are anticipated across aquatic ecosystems in Beijing and beyond. Consequently, we call for the proper implementation of effective management strategies, encompassing early detection, risk assessment, and the use of appropriate control measures to mitigate the potential ecological and economic damages in invaded ecosystems.

Influence of calcium concentration on larval adhesion in a highly invasive fouling ascidian: From morphological changes to molecular mechanisms.

Calcium ion (Ca2+) is involved in the protein-mediated larval adhesion of fouling ascidians, yet the effects of environmental Ca2+ on larval adhesion remain largely unexplored. Here, the larvae of fouling ascidian C. robusta were exposed to different concentrations of Ca2+. Exposures to low-concentration (0 mM and 5 mM) and high-concentration (20 mM and 40 mM) Ca2+ significantly decreased the adhesion rate of larvae, which was primarily attributed to the decreases in adhesive structure length and curvature. Changes in the expressions of genes encoding adhesion-, microvilli-, muscle contraction-, and collagen-related proteins provided a molecular-level explanation for adhesion rate reduction. Additionally, larvae likely prioritized their energy towards immunomodulation in response to Ca2+ stresses, ultimately leading to adhesion reduction. These findings advance our understanding of the influencing mechanisms of environmental Ca2+ on larval adhesion, which are expected to provide references for the development of precise antifouling strategies against ascidians and other fouling species.

eDNA-based detection reveals invasion risks of a biofouling bivalve in the world's largest water diversion project.

Environmental DNA (eDNA) has increasingly been used to detect rare species (e.g., newly introduced nonindigenous species) in both terrestrial and aquatic ecosystems, often with distinct advantages over traditional methods. However, whether water eDNA signals can be used to inform invasion risks remains debatable owing to inherent uncertainties associated with the methods used and the varying conditions among study systems. Here, we sampled eDNA from canals of the central route of the South-to-North Water Diversion Project (hereafter SNWDP) in China to investigate eDNA distribution and efficacy to inform invasion risks in a unique lotic system. We first conducted a total of 16 monthly surveys in this system (two sites in the source reservoir and four sites in the main canal) to test if eDNA could be applied to detect an invasive, biofouling bivalve, the golden mussel Limnoperna fortunei. Second, we initiated a one-time survey in a sub-canal of the SNWDP using refined sampling (12 sites in ~22 km canal) and considered a few environmental predictors. We found that detection of target eDNA in the main canal was achieved up to 1100 km from the putative source population but was restricted to the warmer months (May-November). Detection probability exhibited a significant positive relationship with average daily minimum air temperature and with water temperature, consistent with the expected spawning season. eDNA concentration in the main canal generally fluctuated across months and sites and was generally higher in warmer months. Golden mussel eDNA concentration in the sub-canal decreased significantly with distance from the source and with increasing water temperature and became almost undetectable at ~22 km distance. Given the enormity of the SNWDP, golden mussels may eventually expand their distribution in the main canal, with established "bridgehead" populations facilitating further spread. Our findings suggest an elevated invasion risk of golden mussels in the SNWDP in warm months, highlighting the critical period for spread and, possibly, management.

Interplays between cis- and trans-Acting Factors for Alternative Splicing in Response to Environmental Changes during Biological Invasions of Ascidians.

Alternative splicing (AS), a pivotal biological process contributing to phenotypic plasticity, creates a bridge linking genotypes with phenotypes. Despite its importance, the AS mechanisms underlying environmental response and adaptation have not been well studied, and more importantly, the cis- and trans-acting factors influencing AS variation remain unclear. Using the model invasive congeneric ascidians, Ciona robusta, and Ciona savignyi, we compared their AS responses to environmental changes and explored the potential determinants. Our findings unveiled swift and dynamic AS changes in response to environmental challenges, and differentially alternative spliced genes (DASGs) were functionally enriched in transmembrane transport processes. Interestingly, both the prevalence and level of AS in C. robusta were lower than those observed in C. savignyi. Furthermore, these two indices were higher under temperature stresses compared to salinity stresses in C. savignyi. All the observed patterns underscore the species-specific and environmental context-dependent AS responses to environmental challenges. The dissimilarities in genomic structure and exon/intron size distributions between these two species likely contributed to the observed AS variation. Moreover, we identified a total of 11 and 9 serine/arginine-rich splicing factors (SRSFs) with conserved domains and gene structures in the genomes of C. robusta and C. savignyi, respectively. Intriguingly, our analysis revealed that all detected SRSFs did not exhibit prevalent AS regulations. Instead, we observed AS control over a set of genes related to splicing factors and spliceosome components. Altogether, our results elucidate species-specific and environmental challenge-dependent AS response patterns in closely related invasive ascidians. The identified splicing factors and spliceosome components under AS control offer promising candidates for further investigations into AS-mediated rapid responses to environmental challenges complementary to SRSFs.

Dominance of species sorting over dispersal at microgeographical scales in polluted lotic ecosystems.

Global rivers, particularly those in populated urban areas, are challenged by multiple stressors such as pollution from anthropogenic activities. Dissecting the relative role of each ecological process involved in structuring biotic communities is fundamental in both theoretical and applied ecology. The ecological niche-based species sorting and stochastic dynamics-based dispersal are two major competing processes in determining community structure. Studies have reached a common realization on the environmental gradient-geographical scale relationship (EGGSR), resulting in species sorting over dispersal in determining community structure at large geographical scales where significant environmental gradients often exist. However, this relationship has been recently challenged at fine geographical scales when significant environmental gradients are formed by local anthropogenic activities. Here, we used three receiving rivers of wastewater treatment plants (WWTPs) as the promising system to test the hypothesis that species sorting plays a dominant role over dispersal in structuring zooplankton communities at microgeographical scales (∼1.2 km). After WWTP effluent discharge, we consistently detected significant environmental changes in all three receiving rivers, leading to significant variation in both community structure and taxonomic co-occurrence networks. Variation partitioning showed that environmental variables explained higher proportions of community variation than spatial ones, supporting that species sorting played a dominant role over dispersal in structuring zooplankton communities. Thus, our findings here reject EGGSR, illustrating that the effect of species sorting has been overlooked in disturbed aquatic ecosystems at fine spatial scales. More importantly, all analyses in multiple rivers here validate the "microscale species sorting" hypothesis. The validation of such hypothesis provides a novel methodology for point source pollution management by assessing environment-community interactions and functional changes of biological communities. The differed variables underlying species sorting among three rivers illustrate that ecological management should be case-specific, with the full consideration of local water quality background and pollutant composition of each point pollution source.

Multidimensional plasticity jointly contributes to rapid acclimation to environmental challenges during biological invasions.

Rapid plastic response to environmental changes, which involves extremely complex underlying mechanisms, is crucial for organismal survival during many ecological and evolutionary processes such as those in global change and biological invasions. Gene expression is among the most studied molecular plasticity, while co- or posttranscriptional mechanisms are still largely unexplored. Using a model invasive ascidian Ciona savignyi, we studied multidimensional short-term plasticity in response to hyper- and hyposalinity stresses, covering the physiological adjustment, gene expression, alternative splicing (AS), and alternative polyadenylation (APA) regulations. Our results demonstrated that rapid plastic response varied with environmental context, timescales, and molecular regulatory levels. Gene expression, AS, and APA regulations independently acted on different gene sets and corresponding biological functions, highlighting their nonredundant roles in rapid environmental adaptation. Stress-induced gene expression changes illustrated the use of a strategy of accumulating free amino acids under high salinity and losing/reducing them during low salinity to maintain the osmotic homoeostasis. Genes with more exons were inclined to use AS regulations, and isoform switches in functional genes such as SLC2a5 and Cyb5r3 resulted in enhanced transporting activities by up-regulating the isoforms with more transmembrane regions. The extensive 3'-untranslated region (3'UTR) shortening through APA was induced by both salinity stresses, and APA regulation predominated transcriptomic changes at some stages of stress response. The findings here provide evidence for complex plastic mechanisms to environmental changes, and thereby highlight the importance of systemically integrating different levels of regulatory mechanisms in studying initial plasticity in evolutionary trajectories.

(Epi)genomic adaptation driven by fine geographical scale environmental heterogeneity after recent biological invasions.

Elucidating processes and mechanisms involved in rapid local adaptation to varied environments is a poorly understood but crucial component in management of invasive species. Recent studies have proposed that genetic and epigenetic variation could both contribute to ecological adaptation, yet it remains unclear on the interplay between these two components underpinning rapid adaptation in wild animal populations. To assess their respective contributions to local adaptation, we explored epigenomic and genomic responses to environmental heterogeneity in eight recently colonized ascidian (Ciona intestinalis) populations at a relatively fine geographical scale. Based on MethylRADseq data, we detected strong patterns of local environment-driven DNA methylation divergence among populations, significant epigenetic isolation by environment (IBE), and a large number of local environment-associated epigenetic loci. Meanwhile, multiple genetic analyses based on single nucleotide polymorphisms (SNPs) showed genomic footprints of divergent selection. In addition, for five genetically similar populations, we detected significant methylation divergence and local environment-driven methylation patterns, indicating the strong effects of local environments on epigenetic variation. From a functional perspective, a majority of functional genes, Gene Ontology (GO) terms, and biological pathways were largely specific to one of these two types of variation, suggesting partial independence between epigenetic and genetic adaptation. The methylation quantitative trait loci (mQTL) analysis showed that the genetic variation explained only 18.67% of methylation variation, further confirming the autonomous relationship between these two types of variation. Altogether, we highlight the complementary interplay of genetic and epigenetic variation involved in local adaptation, which may jointly promote populations' rapid adaptive capacity and successful invasions in different environments. The findings here provide valuable insights into interactions between invaders and local environments to allow invasive species to rapidly spread, thus contributing to better prediction of invasion success and development of management strategies.

Complementary genomic and epigenomic adaptation to environmental heterogeneity.

While adaptation is commonly thought to result from selection on DNA sequence-based variation, recent studies have highlighted an analogous epigenetic component as well. However, the relative roles of these mechanisms in facilitating population persistence under environmental heterogeneity remain unclear. To address the underlying genetic and epigenetic mechanisms and their relationship during environmental adaptation, we screened the genomes and epigenomes of nine global populations of a predominately sessile marine invasive tunicate, Botryllus schlosseri. We detected clear population differentiation at the genetic and epigenetic levels. Patterns of genetic and epigenetic structure were significantly influenced by local environmental variables. Among these variables, minimum annual sea surface temperature was identified as the top explanatory variable for both genetic and epigenetic variation. However, patterns of population structure driven by genetic and epigenetic variation were somewhat distinct, suggesting possible autonomy of epigenetic variation. We found both shared and specific genes and biological pathways among genetic and epigenetic loci associated with environmental factors, consistent with complementary and independent contributions of genetic and epigenetic variation to environmental adaptation in this system. Collectively, these mechanisms may facilitate population persistence under environmental change and sustain successful invasions across novel environments.

Neighbours matter: Effects of genomic organization on gene expression plasticity in response to environmental stresses during biological invasions.

Gene expression regulation has been widely recognized as an important molecular mechanism underlying phenotypic plasticity in environmental adaptation. However, it remains largely unexplored on the effects of genomic organization on gene expression plasticity under environmental stresses during biological invasions. Here, we use an invasive model ascidian, Ciona robusta, to investigate how genomic organization affects gene expression in response to salinity stresses during range expansions. Our study showed that neighboring genes were co-expressed and approximately 30% of stress responsive genes were physically clustered on chromosomes. Such coordinated expression was substantially affected by the physical distance and orientation of genes. Interestingly, the overall expression correlation of neighboring genes was significantly decreased under high salinity stresses, illustrating that the co-expression regulation could be disrupted by salinity challenges. Furthermore, the clustering of genes was associated with their function constraints and expression patterns - operon genes enriched in gene expression machinery had the highest transcriptional activity and expression stability. Notably, our analyses showed that the tail-to-tail genes, mainly involved in biological functions related to phosphorylation, homeostatic process, and ion transport, exhibited higher intrinsic expression variability and greater response to salinity challenges. Altogether, the results obtained here provide new insights into the effects of gene organization on gene expression plasticity under environmental challenges, hence improving our knowledge on mechanisms of rapid environmental adaptation during biological invasions.

Environmental DNA captures native and non-native fish community variations across the lentic and lotic systems of a megacity.

Globally, urbanization poses a major threat to terrestrial biodiversity, yet its impact on fish diversity is poorly understood, mainly because of surveying difficulties. In this study, environmental DNA metabarcoding was used to survey fish communities at 109 lentic and lotic sites across Beijing, and how environmental variables affect fish biodiversity at fine urban spatial scales was investigated. We identified 52 native and 23 non-native taxa, with lentic and lotic waters harboring both common and habitat-specific species. Water quality strongly affected native fish diversity, especially in lentic systems, but had little influence on non-native diversity. Fish diversity showed little response to urban land cover variation, but the relative sequence abundance of non-natives in lotic waters increased linearly with distance from the city center. Our findings illustrate the complex effects of urbanization on native versus non-native fishes in different aquatic habitats and highlight the distinctive considerations needed to conserve urban aquatic biodiversity.