Drivers and dynamics of a massive adaptive radiation in cichlid fishes.

Adaptive radiation is the likely source of much of the ecological and morphological diversity of life1-4. How adaptive radiations proceed and what determines their extent remains unclear in most cases1,4. Here we report the in-depth examination of the spectacular adaptive radiation of cichlid fishes in Lake Tanganyika. On the basis of whole-genome phylogenetic analyses, multivariate morphological measurements of three ecologically relevant trait complexes (body shape, upper oral jaw morphology and lower pharyngeal jaw shape), scoring of pigmentation patterns and approximations of the ecology of nearly all of the approximately 240 cichlid species endemic to Lake Tanganyika, we show that the radiation occurred within the confines of the lake and that morphological diversification proceeded in consecutive trait-specific pulses of rapid morphospace expansion. We provide empirical support for two theoretical predictions of how adaptive radiations proceed, the 'early-burst' scenario1,5 (for body shape) and the stages model1,6,7 (for all traits investigated). Through the analysis of two genomes per species and by taking advantage of the uneven distribution of species in subclades of the radiation, we further show that species richness scales positively with per-individual heterozygosity, but is not correlated with transposable element content, number of gene duplications or genome-wide levels of selection in coding sequences.

Loss of Dnmt3a increases self-renewal and resistance to pegIFN-α in JAK2-V617F-positive myeloproliferative neoplasms.

ABSTRACT: Pegylated interferon alfa (pegIFN-α) can induce molecular remissions in patients with JAK2-V617F-positive myeloproliferative neoplasms (MPNs) by targeting long-term hematopoietic stem cells (LT-HSCs). Additional somatic mutations in genes regulating LT-HSC self-renewal, such as DNMT3A, have been reported to have poorer responses to pegIFN-α. We investigated whether DNMT3A loss leads to alterations in JAK2-V617F LT-HSC functions conferring resistance to pegIFN-α treatment in a mouse model of MPN and in hematopoietic progenitors from patients with MPN. Long-term treatment with pegIFN-α normalized blood parameters and reduced splenomegaly and JAK2-V617F chimerism in single-mutant JAK2-V617F (VF) mice. However, pegIFN-α in VF;Dnmt3aΔ/Δ (VF;DmΔ/Δ) mice worsened splenomegaly and failed to reduce JAK2-V617F chimerism. Furthermore, LT-HSCs from VF;DmΔ/Δ mice compared with VF were less prone to accumulate DNA damage and exit dormancy upon pegIFN-α treatment. RNA sequencing showed that IFN-α induced stronger upregulation of inflammatory pathways in LT-HSCs from VF;DmΔ/Δ than from VF mice, indicating that the resistance of VF;DmΔ/Δ LT-HSC was not due to failure in IFN-α signaling. Transplantations of bone marrow from pegIFN-α-treated VF;DmΔ/Δ mice gave rise to more aggressive disease in secondary and tertiary recipients. Liquid cultures of hematopoietic progenitors from patients with MPN with JAK2-V617F and DNMT3A mutation showed increased percentages of JAK2-V617F-positive colonies upon IFN-α exposure, whereas in patients with JAK2-V617F alone, the percentages of JAK2-V617F-positive colonies decreased or remained unchanged. PegIFN-α combined with 5-azacytidine only partially overcame resistance in VF;DmΔ/Δ mice. However, this combination strongly decreased the JAK2-mutant allele burden in mice carrying VF mutation only, showing potential to inflict substantial damage preferentially to the JAK2-mutant clone.

Gene expression remodelling and immune response during adaptive divergence in an African cichlid fish.

Variation in gene expression contributes to ecological speciation by facilitating population persistence in novel environments. Likewise, immune responses can be of relevance in speciation driven by adaptation to different environments. Previous studies examining gene expression differences between recently diverged ecotypes have often relied on only one pair of populations, targeted the expression of only a subset of genes or used wild-caught individuals. Here, we investigated the contribution of habitat-specific parasites and symbionts and the underlying immunological abilities of ecotype hosts to adaptive divergence in lake-river population pairs of the cichlid fish Astatotilapia burtoni. To shed light on the role of phenotypic plasticity in adaptive divergence, we compared parasite and microbiota communities, immune response, and gene expression patterns of fish from natural habitats and a lake-like pond set-up. In all investigated population pairs, lake fish were more heavily parasitized than river fish, in terms of both parasite taxon composition and infection abundance. The innate immune response in the wild was higher in lake than in river populations and was elevated in a river population exposed to lake parasites in the pond set-up. Environmental differences between lake and river habitat and their distinct parasite communities have shaped differential gene expression, involving genes functioning in osmoregulation and immune response. Most changes in gene expression between lake and river samples in the wild and in the pond set-up were based on a plastic response. Finally, gene expression and bacterial communities of wild-caught individuals and individuals acclimatized to lake-like pond conditions showed shifts underlying adaptive phenotypic plasticity.

Time matters! Developmental shift in gene expression between the head and the trunk region of the cichlid fish Astatotilapia burtoni.

BACKGROUND: Differential gene expression can be translated into differing phenotypic traits. Especially during embryogenesis, specific gene expression networks regulate the development of different body structures. Cichlid fishes, with their impressive phenotypic diversity and propensity to radiate, are an emerging model system in the genomics era. Here we set out to investigate gene expression throughout development in the well-studied cichlid fish Astatotilapia burtoni, native to Lake Tanganyika and its affluent rivers. RESULTS: Combining RNA-sequencing from different developmental time points as well as integrating adult gene expression data, we constructed a new genome annotation for A. burtoni comprising 103,253 transcripts (stemming from 52,584 genomic loci) as well as a new reference transcriptome set. We compared our transcriptome to the available reference genome, redefining transcripts and adding new annotations. We show that about half of these transcripts have coding potential. We also characterize transcripts that are not present in the genome assembly. Next, using our newly constructed comprehensive reference transcriptome, we characterized differential gene expression through time and showed that gene expression is shifted between different body parts. We constructed a gene expression network that identified connected genes responsible for particular phenotypes and made use of it to focus on genes under potential positive selection in A. burtoni, which were implicated in fin development and vision. CONCLUSIONS: We provide new genomic resources for the cichlid fish Astatotilapia burtoni, which will contribute to its further establishment as a model system. Tracing gene expression through time, we identified gene networks underlying particular functions, which will help to understand the genetic basis of phenotypic diversity in cichlids.

Thermoregulation and microhabitat choice in the polymorphic asp viper (Vipera aspis).

In ectotherms, thermoregulation strongly depends on environmental conditions, as well as on intrinsic factors, such as skin colour. Indeed, due to its physical properties, melanin pigments allow melanistic morphs to benefit of a more efficient thermoregulation compared to non-melanistic ones. Despite thermal benefits of melanism have often been highlighted under experimental conditions, such field data remain scarce. In this study, we investigated the influence of colouration on body temperature and microhabitat choice in a montane population of colour polymorphic asp viper (Vipera aspis) characterized by a strong presence of melanism (64%). Results highlighted a difference in internal body temperature, but only within gravid females, with melanistic individuals having higher body temperatures compared to non-melanistic ones. No differences were found when considering both sexes. We also showed that melanistic and non-melanistic vipers were found in different microhabitat types, i.e. melanistic snakes used areas marked by a scarcer sun exposure and by higher vegetation cover compared to non-melanistic ones. This result has important implications. Indeed, besides providing a possible explanation for the lack of difference in body temperature (except for gravid females), it confirms that melanistic individuals can potentially use their efficient thermoregulation in order to inhabit less exposed and thermally unfavourable microhabitats.

Sex-biased gene expression and recent sex chromosome turnover.

Cichlids are well known for their propensity to radiate generating arrays of morphologically and ecologically diverse species in short evolutionary time. Following this rapid evolutionary pace, cichlids show high rates of sex chromosome turnover. We here studied the evolution of sex-biased gene (SBG) expression in 14 recently diverged taxa of the Lake Tanganyika Tropheini cichlids, which show different XY sex chromosomes. Across species, sex chromosome sequence divergence predates divergence in expression between the sexes. Only one sex chromosome, the oldest, showed signs of demasculinization in gene expression and potentially contribution to the resolution of sexual conflict. SBGs in general showed high rates of turnovers and evolved mostly under drift. Sexual selection did not shape the rapid evolutionary changes of SBGs. Male-biased genes evolved faster than female-biased genes, which seem to be under more phylogenetic constraint. We found a relationship between the degree of sex bias and sequence evolution driven by sequence differences among the sexes. Consistent with other species, strong sex bias towards sex-limited expression contributes to resolving sexual conflict in cichlids. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part II)'.

Dynamics of sex chromosome evolution in a rapid radiation of cichlid fishes.

Sex is a fundamental trait determined by environmental and/or genetic factors, including sex chromosomes. Sex chromosomes are studied in species scattered across the tree of life, yet little is known about tempo and mode of sex chromosome evolution among closely related species. Here, we examine sex chromosome evolution in the adaptive radiation of cichlid fishes in Lake Tanganyika. Through the analysis of male and female genomes from 244 cichlid taxa (189 described species with 5 represented with two local variants/populations; 50 undescribed species) and of 396 multitissue transcriptomes from 66 taxa, we identify signatures of sex chromosomes in 79 taxa, involving 12 linkage groups. We find that Tanganyikan cichlids have the highest rates of sex chromosome turnover and heterogamety transitions known to date. We show that sex chromosome recruitment is not at random. Moreover convergently emerged sex chromosomes in cichlids support the "limited options" hypothesis of sex chromosome evolution.

Gene expression dynamics during rapid organismal diversification in African cichlid fishes.

Changes in gene expression play a fundamental role in phenotypic evolution. Transcriptome evolutionary dynamics have so far mainly been compared among distantly related species and remain largely unexplored during rapid organismal diversification, in which gene regulatory changes have been suggested as particularly effective drivers of phenotypic divergence. Here we studied gene expression evolution in a model system of adaptive radiation, the cichlid fishes of African Lake Tanganyika. By comparing gene expression profiles of 6 different organs in 74 cichlid species representing all subclades of this radiation, we demonstrate that the rate of gene expression evolution varies among organs, transcriptome parts and the subclades of the radiation, indicating different strengths of selection. We found that the noncoding part of the transcriptome evolved more rapidly than the coding part, and that the gonadal transcriptomes evolved more rapidly than the somatic ones, with the exception of liver. We further show that the rate of gene expression change was not constant over the course of the radiation but accelerated at its later phase. Finally, we show that-at the per-gene level-the evolution of expression patterns is dominated by stabilizing selection.

The Genome of the Toluene-Degrading Pseudomonas veronii Strain 1YdBTEX2 and Its Differential Gene Expression in Contaminated Sand.

The natural restoration of soils polluted by aromatic hydrocarbons such as benzene, toluene, ethylbenzene and m- and p-xylene (BTEX) may be accelerated by inoculation of specific biodegraders (bioaugmentation). Bioaugmentation mainly involves introducing bacteria that deploy their metabolic properties and adaptation potential to survive and propagate in the contaminated environment by degrading the pollutant. In order to better understand the adaptive response of cells during a transition to contaminated material, we analyzed here the genome and short-term (1 h) changes in genome-wide gene expression of the BTEX-degrading bacterium Pseudomonas veronii 1YdBTEX2 in non-sterile soil and liquid medium, both in presence or absence of toluene. We obtained a gapless genome sequence of P. veronii 1YdBTEX2 covering three individual replicons with a total size of 8 Mb, two of which are largely unrelated to current known bacterial replicons. One-hour exposure to toluene, both in soil and liquid, triggered massive transcription (up to 208-fold induction) of multiple gene clusters, such as toluene degradation pathway(s), chemotaxis and toluene efflux pumps. This clearly underlines their key role in the adaptive response to toluene. In comparison to liquid medium, cells in soil drastically changed expression of genes involved in membrane functioning (e.g., lipid composition, lipid metabolism, cell fatty acid synthesis), osmotic stress response (e.g., polyamine or trehalose synthesis, uptake of potassium) and putrescine metabolism, highlighting the immediate response mechanisms of P. veronii 1YdBTEX2 for successful establishment in polluted soil.