Genomic basis of ecological niche divergence among cryptic sister species of non-biting midges

Cryptic Species in Ecotoxicology

The advent of genetic methods has led to the discovery of an increasing number of species that previously could not be distinguished from each other on the basis of morphological characteristics. Even though there has been an exponential growth of publications on cryptic species, such species are rarely considered in ecotoxicology. Thus, the particular question of ecological differentiation and the sensitivity of closely related cryptic species is rarely addressed. Tackling this question, however, is of key importance for evolutionary ecology, conservation biology, and, in particular, regulatory ecotoxicology. At the same time, the use of species with (known or unknown) cryptic diversity might be a reason for the lack of reproducibility of ecotoxicological experiments and implies a false extrapolation of the findings. Our critical review includes a database and literature search through which we investigated how many of the species most frequently used in ecotoxicological assessments show evidence of cryptic diversity. We found a high proportion of reports indicating overlooked species diversity, especially in invertebrates. In terrestrial and aquatic realms, at least 67% and 54% of commonly used species, respectively, were identified as cryptic species complexes. The issue is less prominent in vertebrates, in which we found evidence for cryptic species complexes in 27% of aquatic and 6.7% of terrestrial vertebrates. We further exemplified why different evolutionary histories may significantly determine cryptic species' ecology and sensitivity to pollutants. This in turn may have a major impact on the results of ecotoxicological tests and, consequently, the outcome of environmental risk assessments. Finally, we provide a brief guideline on how to deal practically with cryptic diversity in ecotoxicological studies in general and its implementation in risk assessment procedures in particular. Environ Toxicol Chem 2023;42:1889-1914. © 2023 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Two novel, tightly linked, and rapidly evolving genes underlie Aedes aegypti mosquito reproductive resilience during drought

Female Aedes aegypti mosquitoes impose a severe global public health burden as vectors of multiple viral pathogens. Under optimal environmental conditions, Aedes aegypti females have access to human hosts that provide blood proteins for egg development, conspecific males that provide sperm for fertilization, and freshwater that serves as an egg-laying substrate suitable for offspring survival. As global temperatures rise, Aedes aegypti females are faced with climate challenges like intense droughts and intermittent precipitation, which create unpredictable, suboptimal conditions for egg-laying. Here, we show that under drought-like conditions simulated in the laboratory, females retain mature eggs in their ovaries for extended periods, while maintaining the viability of these eggs until they can be laid in freshwater. Using transcriptomic and proteomic profiling of Aedes aegypti ovaries, we identify two previously uncharacterized genes named tweedledee and tweedledum, each encoding a small, secreted protein that both show ovary-enriched, temporally-restricted expression during egg retention. These genes are mosquito-specific, linked within a syntenic locus, and rapidly evolving under positive selection, raising the possibility that they serve an adaptive function. CRISPR-Cas9 deletion of both tweedledee and tweedledum demonstrates that they are specifically required for extended retention of viable eggs. These results highlight an elegant example of taxon-restricted genes at the heart of an important adaptation that equips Aedes aegypti females with 'insurance' to flexibly extend their reproductive schedule without losing reproductive capacity, thus allowing this species to exploit unpredictable habitats in a changing world.

Pathotypes of Fusarium oxysporum f. sp. fragariae express discrete repertoires of accessory genes and induce distinct host transcriptional responses during root infection

Convergent evolution of phytopathogenicity is poorly described, especially among multiple strains of a single microbial species. We investigated this phenomenon with genetically diverse isolates of Fusarium oxysporum f. sp. fragariae (Fof) that cause one of two syndromes: chlorosis and wilting (the 'yellows-fragariae' pathotype), or only wilting (the 'wilt-fragariae' pathotype). We challenged strawberry (Fragaria × ananassa) plants to root infection by five fungal isolates: three yellows-fragariae, one wilt-fragariae and one that is not pathogenic to strawberry. All Fof isolates had chromosome-level assemblies; three were newly generated. The two pathotypes triggered distinct host responses, especially among phytohormone-associated genes; yellows-fragariae isolates strongly induced jasmonic acid-associated genes, whereas the wilt-fragariae isolate primarily induced ethylene biosynthesis and signalling. The differentially expressed genes on fungal accessory chromosomes were almost entirely distinct between pathotypes. We identified an ~150 kbp 'pathogenicity island' that was horizontally transferred between wilt-fragariae strains. This predicted pathogenicity island was enriched with differentially expressed genes whose predicted functions were related to plant infection, and only one of these genes was also upregulated in planta by yellows-fragariae isolates. These results support the conclusion that wilt- and yellows-fragariae cause physiologically distinct syndromes by the expression of discrete repertoires of genes on accessory chromosomes.

Horizontal chromosome transfer and independent evolution drive diversification in Fusarium oxysporum f. sp. fragariae

The genes required for host-specific pathogenicity in Fusarium oxysporum can be acquired through horizontal chromosome transfer (HCT). However, it is unknown if HCT commonly contributes to the diversification of pathotypes. Using comparative genomics and pathogenicity phenotyping, we explored the role of HCT in the evolution of F. oxysporum f. sp. fragariae, the cause of Fusarium wilt of strawberry, with isolates from four continents. We observed two distinct syndromes: one included chlorosis ('yellows-fragariae') and the other did not ('wilt-fragariae'). All yellows-fragariae isolates carried a predicted pathogenicity chromosome, 'chrY-frag ', that was horizontally transferred at least four times. chrY-frag was associated with virulence on specific cultivars and encoded predicted effectors that were highly upregulated during infection. chrY-frag was not present in wilt-fragariae; isolates causing this syndrome evolved pathogenicity independently. All origins of F. oxysporum f. sp. fragariae occurred outside of the host's native range. Our data support the conclusion that HCT is widespread in F. oxysporum, but pathogenicity can also evolve independently. The absence of chrY-frag in wilt-fragariae suggests that multiple, distinct pathogenicity chromosomes can confer the same host specificity. The wild progenitors of cultivated strawberry (Fragaria × ananassa) did not co-evolve with this pathogen, yet we discovered several sources of genetic resistance.

Genomic divergence landscape in recurrently hybridizing Chironomus sister taxa suggests stable steady state between mutual gene flow and isolation

Divergence is mostly viewed as a progressive process often initiated by selection targeting individual loci, ultimately resulting in ever increasing genomic isolation due to linkage. However, recent studies show that this process may stall at intermediate stable equilibrium states without achieving complete genomic isolation. We tested the extent of genomic isolation between two recurrently hybridizing nonbiting midge sister taxa, Chironomus riparius and Chironomus piger, by analyzing the divergence landscape. Using a principal component-based method, we estimated that only about 28.44% of the genomes were mutually isolated, whereas the rest was still exchanged. The divergence landscape was fragmented into isolated regions of on average 30 kb, distributed throughout the genome. Selection and divergence time strongly influenced lengths of isolated regions, whereas local recombination rate only had minor impact. Comparison of divergence time distributions obtained from several coalescence-simulated divergence scenarios with the observed divergence time estimates in an approximate Bayesian computation framework favored a short and concluded divergence event in the past. Most divergence happened during a short time span about 4.5 million generations ago, followed by a stable equilibrium between mutual gene flow through ongoing hybridization for the larger part of the genome and isolation in some regions due to rapid purifying selection of introgression, supported by high effective population sizes and recombination rates.

The Globin Gene Family in Arthropods: Evolution and Functional Diversity

Globins are small heme-proteins that reversibly bind oxygen. Their most prominent roles in vertebrates are the transport and storage of O₂ for oxidative energy metabolism, but recent research has suggested alternative, non-respiratory globin functions. In the species-rich and ecologically highly diverse taxon of arthropods, the copper-containing hemocyanin is considered the main respiratory protein. However, recent studies have suggested the presence of globin genes and their proteins in arthropod taxa, including model species like Drosophila. To systematically assess the taxonomic distribution, evolution and diversity of globins in arthropods, we systematically searched transcriptome and genome sequence data and found a conserved, widespread occurrence of three globin classes in arthropods: hemoglobin-like (HbL), globin X (GbX), and globin X-like (GbXL) protein lineages. These globin types were previously identified in protostome and deuterostome animals including vertebrates, suggesting their early ancestry in Metazoa. The HbL genes show multiple, lineage-specific gene duplications in all major arthropod clades. Some HbL genes (e.g., Glob2 and 3 of Drosophila) display particularly fast substitution rates, possibly indicating the evolution of novel functions, e.g., in spermatogenesis. In contrast, arthropod GbX and GbXL globin genes show high evolutionary stability: GbXL is represented by a single-copy gene in all arthropod groups except Brachycera, and representatives of the GbX clade are present in all examined taxa except holometabolan insects. GbX and GbXL both show a brain-specific expression. Most arthropod GbX and GbXL proteins, but also some HbL variants, include sequence motifs indicative of potential N-terminal acylation (i.e., N-myristoylation, 3C-palmitoylation). All arthropods except for the brachyceran Diptera harbor at least one such potentially acylated globin copy, confirming the hypothesis of an essential, conserved globin function associated with the cell membrane. In contrast to other animals, the fourth ancient globin lineage, represented by neuroglobin, appears to be absent in arthropods, and the putative arthropod orthologs of the fifth metazoan globin lineage, androglobin, lack a recognizable globin domain. Thus, the remarkable evolutionary stability of some globin variants is contrasted by occasional dynamic gene multiplication or even loss of otherwise strongly conserved globin lineages in arthropod phylogeny.

A High-Quality Genome Assembly from Short and Long Reads for the Non-biting Midge Chironomus riparius (Diptera)

Chironomus riparius is of great importance as a study species in various fields like ecotoxicology, molecular genetics, developmental biology and ecology. However, only a fragmented draft genome exists to date, hindering the recent rush of population genomic studies in this species. Making use of 50 NGS datasets, we present a hybrid genome assembly from short and long sequence reads that make C. riparius’ genome one of the most contiguous Dipteran genomes published, the first complete mitochondrial genome of the species, and the respective recombination rate among the first insect recombination rates at all. The genome assembly and associated resources will be highly valuable to the broad community working with dipterans in general and chironomids in particular. The estimated recombination rate will help evolutionary biologists gaining a better understanding of commonalities and differences of genomic patterns in insects.

The genomic footprint of climate adaptation in Chironomus riparius

The gradual heterogeneity of climatic factors poses varying selection pressures across geographic distances that leave signatures of clinal variation in the genome. Separating signatures of clinal adaptation from signatures of other evolutionary forces, such as demographic processes, genetic drift and adaptation, to nonclinal conditions of the immediate local environment is a major challenge. Here, we examine climate adaptation in five natural populations of the harlequin fly Chironomus riparius sampled along a climatic gradient across Europe. Our study integrates experimental data, individual genome resequencing, Pool-Seq data and population genetic modelling. Common-garden experiments revealed significantly different population growth rates at test temperatures corresponding to the population origin along the climate gradient, suggesting thermal adaptation on the phenotypic level. Based on a population genomic analysis, we derived empirical estimates of historical demography and migration. We used an F_ST outlier approach to infer positive selection across the climate gradient, in combination with an environmental association analysis. In total, we identified 162 candidate genes as genomic basis of climate adaptation. Enriched functions among these candidate genes involved the apoptotic process and molecular response to heat, as well as functions identified in studies of climate adaptation in other insects. Our results show that local climate conditions impose strong selection pressures and lead to genomic adaptation despite strong gene flow. Moreover, these results imply that selection to different climatic conditions seems to converge on a functional level, at least between different insect species.

Population genetic structure and hybridization patterns in the cryptic sister species Chironomus riparius and Chironomus piger across differentially polluted freshwater systems

Chironomids are an integral and functionally important part of many freshwater ecosystems. Yet, to date, there is limited understanding of their microevolutionary processes under chemically polluted natural environments. In this study, we investigated the genetic variation within populations of the ecotoxicological model species Chironomus riparius and its cryptic sister species Chironomus piger at 18 metal-contaminated and reference sites in northwestern Portugal. Microsatellite analysis was conducted on 909 samples to answer if metal contamination affects genetic variation in natural chironomid populations as previously suggested from controlled laboratory experiments. Similarly high levels of genetic diversity and significant but weak genetic substructuring were found across all sites and temporal replicates, with no effects of metal contamination on the genetic variation or species' abundance, although C. piger tended to be less frequent at highly contaminated sites. Our results indicate that high levels of gene flow and population dynamic processes may overlay potential pollutant effects. At least for our study species, we conclude that the "genetic erosion hypothesis", which suggests that chemical pollution will reduce genome-wide genetic variability in affected populations, does not hold under natural conditions. Interestingly, our study provides evidence of successful hybridization between the two sister species under natural conditions.

Support for the evolutionary speed hypothesis from intraspecific population genetic data in the non-biting midge Chironomus riparius

The evolutionary speed hypothesis (ESH) proposes a causal mechanism for the latitudinal diversity gradient. The central idea of the ESH is that warmer temperatures lead to shorter generation times and increased mutation rates. On an absolute time scale, both should lead to an acceleration of selection and drift. Based on the ESH, we developed predictions regarding the distribution of intraspecific genetic diversity: populations of ectothermic species with more generations per year owing to warmer ambient temperatures should be more differentiated from each other, accumulate more mutations and show evidence for increased mutation rates compared with populations in colder regions. We used the multivoltine insect species Chironomus riparius to test these predictions with cytochrome oxidase I (COI) sequence data and found that populations from warmer regions are indeed significantly more differentiated and have significantly more derived haplotypes than populations from colder regions. We also found a significant correlation of the annual mean temperature with the population mutation parameter θ that serves as a proxy for the per generation mutation rate under certain assumptions. This pattern could be corroborated with two nuclear loci. Overall, our results support the ESH and indicate that the thermal regime experienced may be crucially driving the evolution of ectotherms and may thus ultimately govern their speciation rate.

Ecological genomics in Xanthomonas: the nature of genetic adaptation with homologous recombination and host shifts

BACKGROUND

Comparative genomics provides insights into the diversification of bacterial species. Bacterial speciation usually takes place with lasting homologous recombination, which not only acts as a cohering force between diverging lineages but brings advantageous alleles favored by natural selection, and results in ecologically distinct species, e.g., frequent host shift in Xanthomonas pathogenic to various plants.

RESULTS

Using whole-genome sequences, we examined the genetic divergence in Xanthomonas campestris that infected Brassicaceae, and X. citri, pathogenic to a wider host range. Genetic differentiation between two incipient races of X. citri pv. mangiferaeindicae was attributable to a DNA fragment introduced by phages. In contrast to most portions of the genome that had nearly equivalent levels of genetic divergence between subspecies as a result of the accumulation of point mutations, 10% of the core genome involving with homologous recombination contributed to the diversification in Xanthomonas, as revealed by the correlation between homologous recombination and genomic divergence. Interestingly, 179 genes were under positive selection; 98 (54.7%) of these genes were involved in homologous recombination, indicating that foreign genetic fragments may have caused the adaptive diversification, especially in lineages with nutritional transitions. Homologous recombination may have provided genetic materials for the natural selection, and host shifts likely triggered ecological adaptation in Xanthomonas. To a certain extent, we observed positive selection nevertheless contributed to ecological divergence beyond host shifting.

CONCLUSION

Altogether, mediated with lasting gene flow, species formation in Xanthomonas was likely governed by natural selection that played a key role in helping the deviating populations to explore novel niches (hosts) or respond to environmental cues, subsequently triggering species diversification.

The effect of temperature gradients on endocrine signaling and antioxidant gene expression during Chironomus riparius development

Temperature is one of the most important environmental factors affecting the biological processes of aquatic species. To investigate the potential effects of temperature on the developmental processes of aquatic invertebrates, we analyzed biological and molecular transcriptional responses during Chironomus riparius development, including five stages spanning from embryo to adult stages. We assessed the temperature change-induced reduction of survival rate, changes in biological development including the male:female ratio in emerged adults, the success rates of pupation and emergence, and the developmental timing of pupation and emergence. The increased temperature induced expression of endocrine signaling genes, such as the ecdysone receptor, ultraspiracle (ortholog of the RXR), and the estrogen-related receptor in the fourth-instar larval and pupal stages of C. riparius development. Altered temperature also affected the activity of antioxidant genes, including catalase, peroxidase, glutathione peroxidase, and superoxide dismutase during the fourth-instar larval to adult stages of C. riparius development, as a result of altered development. Increased temperature during the fourth-instar larval stage increased oxidative stress in pupae and adults. Responses of antioxidant genes to increased temperature occurred in a developmental stage-dependent manner. However, reduced temperature did not induce the expression of antioxidant genes in a developmental stage-dependent manner, although it did induce oxidative stress during C. riparius development. Increased temperature also caused greater toxicity of di-ethylhexyl phthalate (DEHP) in fourth-instar larvae. Our findings suggest that altered temperatures may disturb the invertebrate hormone system and developmental processes by inducing oxidative stress in aquatic environments.