Human health risk assessment for (re)emerging protozoan parasites in surface water used for public supply and recreational activities.

Gastrointestinal diseases caused by protozoan parasites remain a major challenge in developing countries and ingestion of contaminated surface water represents one of the main sources by which these diseases are contracted. This study assessed the risk of infection and diseases caused by Cryptosporidium spp. and Giardia sp. due to ingestion of surface water used for public supply and recreational activities, focusing on the southeastern Brazilian Pardo River and applying the USEPA 1623 method to quantify (oo)cyst concentrations. Infection and disease probabilities due to ingestion of drinking water or during recreational activities were estimated using the Quantitative Microbial Risk Assessment (QMRA) approach. Mean concentrations of Cryptosporidium spp. and Giardia sp. in surface water ranged from 0.2 to 0.4 oocysts L-1 and 0.2 to 4.4 cysts L-1, respectively. Considering public water supply, annual infection probabilities were higher for adults than children and exceeded the USEPA limit; also, disease probabilities were higher for adults than children. For recreational activities, annual infection and disease probabilities were higher for children, followed by men and women. The occurrence of both parasites likely reflects raw sewage discharge, effluent from sewage treatment plants, and diffuse sources of pollution, such as runoff from pasture lands and deforested riparian forest corridors. Our results highlight substantial infection risks by both parasite types after conventional treatment of water used for public supply and also call for careful monitoring of water bodies used for recreational purposes.

Shifts in bacterial communities and antibiotic resistance genes in surface water and gut microbiota of guppies (Poecilia reticulata) in the upper Rio Uberabinha, Brazil.

Anthropogenic activities especially water pollution can affect the diversity and composition of microbial communities and promote the spread of antibiotic resistance genes (ARGs). In this study, water samples and guppies (Poecilia reticulata) were sampled from six sampling sites along the Uberabinha River in southeastern Brazil, both microbial communities and ARGs of surface waters and intestinal microbiota of guppies (Poecilia reticulata) were detected. According to the results of 16S rRNA amplicon sequencing, Proteobacteria, Bacteroidetes, Firmicutes and Actinobacteria were dominant phyla in both water and intestinal microbiota, but the abundance of putative pathogens was higher at heavily polluted sites. Up to 83% of bacteria in intestinal microbiota originated from water microbiota; this proportion was relatively higher in less polluted compared to polluted environments. ARGs providing resistance of tetracyclines and quinolones were dominant in both water and gut microbiota. The relative abundances of class I integrons and ARGs were as high as 1.74 × 10-1/16S rRNA copies and 3.61 × 10-1/16S rRNA copies, respectively, at heavily polluted sites. Correlation analysis suggests that integrons and bacteria play key roles in explaining the widespread occurrence of ARGs in the surface, but not in intestinal microbiota. We could rule out the class I integrons a potential intermediary bridge for ARGs between both types of microbiomes. Our results highlight the tight link in microbial communities and ARGs between ambient microbiota of stream ecosystems and intestinal microbiota of fish. Our study could have far-reaching consequences for fisheries and consumer safety and calls for investigations of gut microbiota of target species of both commercial fisheries and recreational (hobby) angling.

Geographical and temporal variation of multiple paternity in invasive mosquitofish (Gambusia holbrooki, Gambusia affinis).

Multiple paternity (MP) increases offspring's genetic variability, which could be linked to invasive species' evolvability in novel distribution ranges. Shifts in MP can be adaptive, with greater MP in harsher/colder environments or towards the end of the reproductive season, but climate could also affect MP indirectly via its effect on reproductive life histories. We tested these hypotheses by genotyping N = 2,903 offspring from N = 306 broods of two closely related livebearing fishes, Gambusia holbrooki and Gambusia affinis. We sampled pregnant females across latitudinal gradients in their invasive ranges in Europe and China, and found more sires per brood and a greater reproductive skew towards northern sampling sites. Moreover, examining monthly sampling from two G. affinis populations, we found MP rates to vary across the reproductive season in a northern Chinese, but not in a southern Chinese population. While our results confirm an increase of MP in harsher/more unpredictable environments, path analysis indicated that, in both cases, the effects of climate are likely to be indirect, mediated by altered life histories. In both species, which rank amongst the 100 most invasive species worldwide, higher MP at the northern edge of their distribution probably increases their invasive potential and favours range expansions, especially in light of the predicted temperature increases due to global climate changes.

Small-scale phenotypic differentiation along complex stream gradients in a non-native amphipod.

BACKGROUND: Selective landscapes in rivers are made up by an array of selective forces that vary from source to downstream regions or between seasons, and local/temporal variation in fitness maxima can result in gradual spatio-temporal variation of phenotypic traits. This study aimed at establishing freshwater amphipods as future model organisms to study adaptive phenotypic diversification (evolutionary divergence and/or adaptive plasticity) along stream gradients. METHODS: We collected Gammarus roeselii from 16 sampling sites in the Rhine catchment during two consecutive seasons (summer and winter). Altogether, we dissected n = 1648 individuals and quantified key parameters related to morphological and life-history diversification, including naturally selected (e.g., gill surface areas) as well as primarily sexually selected traits (e.g., male antennae). Acknowledging the complexity of selective regimes in streams and the interrelated nature of selection factors, we assessed several abiotic (e.g., temperature, flow velocity) and biotic ecological parameters (e.g., conspecific densities, sex ratios) and condensed them into four principal components (PCs). RESULTS: Generalized least squares models revealed pronounced phenotypic differentiation in most of the traits investigated herein, and components of the stream gradient (PCs) explained parts of the observed differences. Depending on the trait under investigation, phenotypic differentiation could be ascribed to variation in abiotic conditions, anthropogenic disturbance (influx of thermally polluted water), or population parameters. For example, female fecundity showed altitudinal variation and decreased with increasing conspecific densities, while sexual dimorphism in the length of male antennae-used for mate finding and assessment-increased with increasing population densities and towards female-biased sex ratios. CONCLUSIONS: We provide a comprehensive protocol for comparative analyses of intraspecific variation in life history traits in amphipods. Whether the observed phenotypic differentiation over small geographical distances reflects evolutionary divergence or plasticity (or both) remains to be investigated in future studies. Independent of the mechanisms involved, variation in several traits is likely to have consequences for ecosystem functions. For example, leaf-shredding in G. roeselii strongly depends on body size, which varied in dependence of several ecological parameters.

Elevated temperatures translate into reduced dispersal abilities in a natural population of an aquatic insect.

Rising global temperatures force many species to shift their distribution ranges. However, whether or not (and how fast) such range shifts occur depends on species' dispersal capacities. In most ecological studies, dispersal-related traits (such as the wing size or wing loading in insects) are treated as fixed, species-specific characteristics, ignoring the important role of phenotypic plasticity during insect development. We tested the hypothesis that dispersal-related traits themselves vary in dependence of ambient environmental conditions (temperature regimes, discharge patterns and biotic interactions during individual development). We collected data over 8 years from a natural population of the crane fly Tipula maxima in central Germany. Using linear mixed-effect models, we analysed how phenotypic traits, phenological characteristics and population densities are affected by environmental conditions during the preceding 3, 6 and 12 months. We found a moderate (5.6%) increase in wing length per 1°C increase in mean annual temperatures during the previous year. At the same time, body weight increased by as much as 17.8% in females and 26.9% in males per 1°C, likely driven by increased habitat productivity, which resulted in a 16.4% (female) and 19.3% (male) increased wing loading. We further found a shorter, more synchronized emergence period (i.e. a narrower time frame for dispersal) with increasing temperatures. Altogether, our results suggest that dispersal abilities of T. maxima were negatively affected by elevated temperatures, and we discuss how similar patterns might affect the persistence of populations of other aquatic insects, especially stenoecious taxa with narrow distribution ranges. Our study calls for integration of information on temperature-induced phenotypic plasticity of dispersal-related traits into models forecasting range shifts in the face of climate change. Furthermore, the patterns reported here are likely to affect metapopulation dynamics of aquatic insects under climate change conditions and may contribute to the ongoing decline of insect biomass and diversity.

Biological water quality assessment in the degraded Mutara rangelands, northeastern Rwanda.

Rwanda is a heavily overpopulated country that also suffers from overstocking with livestock, especially following the return of war refuges after the civil war (1991-1995). At present, approximately 20% of the human population in Nyagatare District in northeastern Rwanda has no access to clean drinking water and sanitation. We used a biotic index based on the presence of selected families of aquatic macroinvertebrates, derived from the "Tanzania River Scoring System" (TARISS), to assess water quality at N = 55 sites in the Mutara grasslands in Nyagatare District. Poor water quality became evident across most sampling sites both in the Muvumba (mean ± SE TARISS score 5.25 ± 0.10) and Karangazi Rivers (4.79 ± 0.12). Using a general linear model, we asked whether direct effects of land use forms and input of anthropogenic wastewater have an impact on water quality. Our results found no immediate effects of both forms of disturbance/pollution, probably because overall water quality was already poor. Our study is intended to serve as a starting point for continuous monitoring of water quality in the Mutara rangelands in NE Rwanda. The method applied here is cost-efficient, requires only basic equipment, and training local students to apply this technique can provide a solid basis for its implementation in future surveys related to public health.

CircFUT10 reduces proliferation and facilitates differentiation of myoblasts by sponging miR-133a.

Circular RNAs (circRNAs) have been identified in various tissues and cell types from human, monkey, porcine, and mouse. However, knowledge on circRNAs in bovine muscle development is limited. We downloaded and analyzed the circRNAs sequencing data of bovine skeletal muscle tissue, and further characterized the role of a candidate circRNA (circFUT10) in muscle development. Quantitative real-time PCR (qPCR) and Western blot assays were used to confirm the expression of genes involved in myoblasts differentiation and proliferation. Flow cytometry was performed to assess cell cycle distribution and cell apoptosis. EdU incorporation and CCK-8 assay were performed to demonstrate cell proliferation. We demonstrated that circFUT10 was highly (but differentially) expressed in embryonic and adult skeletal muscle tissue. circFUT10 induced bovine primary myoblasts differentiation and increased the expression of MyoD, MyoG, and MyhC in mRNA and protein levels. circFUT10 increased the number of myoblasts in the G0/G1 phase of the cell cycle, and decreased the proportion of cells in the S-phase. circFUT10 inhibited the proliferation of myoblasts and promoted them apoptosis. Via a luciferase screening assay, circFUT10 is observed to sponge to miR-133a with three potential binding sites. Specifically, we show that circFUT10 regulated myoblasts differentiation and cell survival by directly binding to miR-133a and inhibiting miR-133a activity. Modulation of circFUT10 expression in muscle tissue may emerge as a potential target in breeding strategies attempting to control muscle development in cattle.

Extreme environments and the origins of biodiversity: Adaptation and speciation in sulphide spring fishes.

Organisms adapted to physiochemical stressors provide ideal systems to study evolutionary mechanisms that drive adaptation and speciation. This review study focuses on livebearing fishes of the Poecilia mexicana species complex (Poeciliidae), members of which have repeatedly colonized hydrogen sulphide (H2 S)-rich springs. H2 S is a potent respiratory toxicant that creates extreme environmental conditions in aquatic ecosystems. There is also a rich history of research on H2 S in toxicology and biomedicine, which has facilitated the generation of a priori hypotheses about the proximate mechanisms of adaptation. Testing these hypotheses through the application of high-throughput genomic and transcriptomic analyses has led to the identification of the physiological underpinnings mediating adaptation to H2 S-rich environments. In addition, systematic natural history studies have provided a nuanced understanding of how the presence of a physiochemical stressor interacts with other sources of selection to drive evolutionary change in a variety of organismal traits, including physiology, morphology, behaviour and life history. Adaptation to extreme environments in P. mexicana also coincides with ecological speciation, and evolutionarily independent lineages span almost the full range of the speciation continuum from panmixia to complete reproductive isolation. Multiple mechanisms of reproductive isolation are involved in reducing gene flow between adjacent populations that are adapted to contrasting environmental conditions. Comparative studies among evolutionarily independent lineages within the P. mexicana species complex and, more recently, other members of the family Poeciliidae that have colonized H2 S-rich environments will provide insights into the factors facilitating or impeding convergent evolution, providing tangible links between micro-evolutionary processes and macro-evolutionary patterns.

Consistent individual differences in associative learning speed are not linked to boldness in female Atlantic mollies.

Recent studies on consistent individual differences in behavioural tendencies (animal personality) raised the question of whether individual differences in cognitive abilities can be linked to certain personality types. We tested female Atlantic mollies (Poecilia mexicana) in two different classical conditioning experiments. For the first time, we provide evidence for highly consistent individual differences in associative learning speed in fish. We characterized the same individuals for boldness in two experimental situations (latency to emerge from shelter and freezing time after a simulated predator attack) and found high behavioural repeatability. When we tested for a potential correlation between associative learning speed and boldness, however, there was no evidence for a link between them. Our study design included several steps to avoid typical pitfalls of disadvantaging shy individuals during learning tests. We caution that other experimental studies may have suffered from erroneous interpretations due to a more cautious coping style of shy individuals in the respective setup used to assess learning.

Evolution in caves: selection from darkness causes spinal deformities in teleost fishes.

Only few fish species have successfully colonized subterranean habitats, but the underlying biological constraints associated with this are still poorly understood. Here, we investigated the influence of permanent darkness on spinal-column development in one species (Midas cichlid, Amphilophus citrinellus) with no known cave form, and one (Atlantic molly, Poecilia mexicana) with two phylogenetically young cave forms. Specifically, fish were reared under a normal light : dark cycle or in permanent darkness (both species). We also surveyed wild-caught cave and surface ecotypes of P. mexicana In both species, permanent darkness was associated with significantly higher rates of spinal deformities (especially in A. citrinellus). This suggests strong developmental (intrinsic) constraints on the successful colonization of subterranean environments in teleost fishes and might help explain the relative paucity of cave-adapted lineages. Our results add depth to our understanding of the aspects of selection driving trait divergence and maintaining reproductive isolation in cave faunas.

Long non-coding RNA ADNCR suppresses adipogenic differentiation by targeting miR-204.

Adipogenesis is a complex and precisely orchestrated process mediated by a network of adipogenic regulatory factors. Several studies have highlighted the relevance of lncRNAs in adipocyte differentiation, but the precise molecular mechanism has largely remained elusive. In the present study, we performed Ribo-Zero RNA-Seq to investigate both the poly(A)+and poly(A)-lncRNAs of in vitro cultured bovine preadipocytes and differentiated adipocytes. A stringent set of 2882 lncRNAs was finally identified. A comparison of the lncRNAs expression profiles revealed that 16 lncRNAs are differentially expressed during adipocyte differentiation. We focused on the most downregulated lncRNA, which we named adipocyte differentiation-associated long noncoding RNA (ADNCR). Mechanistically, ADNCR inhibited adipocyte differentiation by functioning as a competing endogenous RNA (ceRNA) for miR-204, thereby augmenting the expression of the miR-204 target gene, SIRT1, which is known to inhibit adipocyte differentiation and adipogenic gene expression by docking with NCoR and SMART to repress PPARγ activity. Our data not only provide a valuable genomic resource for the identification of lncRNAs with functional roles in adipocyte differentiation but also reveal new insights into understanding the mechanisms of adipogenic differentiation.

Extremophile Poeciliidae: multivariate insights into the complexity of speciation along replicated ecological gradients.

BACKGROUND: Replicate population pairs that diverge in response to similar selective regimes allow for an investigation of (a) whether phenotypic traits diverge in a similar and predictable fashion, (b) whether there is gradual variation in phenotypic divergence reflecting variation in the strength of natural selection among populations, (c) whether the extent of this divergence is correlated between multiple character suites (i.e., concerted evolution), and (d) whether gradual variation in phenotypic divergence predicts the degree of reproductive isolation, pointing towards a role for adaptation as a driver of (ecological) speciation. Here, we use poeciliid fishes of the genera Gambusia and Poecilia that have repeatedly evolved extremophile lineages able to tolerate high and sustained levels of toxic hydrogen sulfide (H2S) to answer these questions. RESULTS: We investigated evolutionary divergence in response to H2S in Gambusia spp. (and to a lesser extent Poecilia spp.) using a multivariate approach considering the interplay of life history, body shape, and population genetics (nuclear miscrosatellites to infer population genetic differentiation as a proxy for reproductive isolation). We uncovered both shared and unique patterns of evolution: most extremophile Gambusia predictably evolved larger heads and offspring size, matching a priori predictions for adaptation to sulfidic waters, while variation in adult life histories was idiosyncratic. When investigating patterns for both genera (Gambusia and Poecilia), we found that divergence in offspring-related life histories and body shape were positively correlated across populations, but evidence for individual-level associations between the two character suites was limited, suggesting that genetic linkage, developmental interdependencies, or pleiotropic effects do not explain patterns of concerted evolution. We further found that phenotypic divergence was positively correlated with both environmental H2S-concentration and neutral genetic differentiation (a proxy for gene flow). CONCLUSIONS: Our results suggest that higher toxicity exerts stronger selection, and that divergent selection appears to constrain gene flow, supporting a scenario of ecological speciation. Nonetheless, progress toward ecological speciation was variable, partially reflecting variation in the strength of divergent selection, highlighting the complexity of selective regimes even in natural systems that are seemingly governed by a single, strong selective agent.

Does personality affect premating isolation between locally-adapted populations?

BACKGROUND: One aspect of premating isolation between diverging, locally-adapted population pairs is female mate choice for resident over alien male phenotypes. Mating preferences often show considerable individual variation, and whether or not certain individuals are more likely to contribute to population interbreeding remains to be studied. In the Poecilia mexicana-species complex different ecotypes have adapted to hydrogen sulfide (H2S)-toxic springs, and females from adjacent non-sulfidic habitats prefer resident over sulfide-adapted males. We asked if consistent individual differences in behavioral tendencies (animal personality) predict the strength and direction of the mate choice component of premating isolation in this system. RESULTS: We characterized focal females for their personality and found behavioral measures of 'novel object exploration', 'boldness' and 'activity in an unknown area' to be highly repeatable. Furthermore, the interaction term between our measures of exploration and boldness affected focal females' strength of preference (SOP) for the resident male phenotype in dichotomous association preference tests. High exploration tendencies were coupled with stronger SOPs for resident over alien mating partners in bold, but not shy, females. Shy and/or little explorative females had an increased likelihood of preferring the non-resident phenotype and thus, are more likely to contribute to rare population hybridization. When we offered large vs. small conspecific stimulus males instead, less explorative females showed stronger preferences for large male body size. However, this effect disappeared when the size difference between the stimulus males was small. CONCLUSIONS: Our results suggest that personality affects female mate choice in a very nuanced fashion. Hence, population differences in the distribution of personality types could be facilitating or impeding reproductive isolation between diverging populations depending on the study system and the male trait(s) upon which females base their mating decisions, respectively.

Phylogenetic analyses of gazelles reveal repeated transitions of key ecological traits and provide novel insights into the origin of the genus Gazella.

African bovids are a famous example of a taxonomic group in which the correlated evolution of body size, feeding mode, gregariousness, and social organization in relation to the preferred habitat type has been investigated. A continuum has been described ranging from small-bodied, sedentary, solitary or socially monogamous, forest- or bush-dwelling, browsing species that seek shelter from predation in dense vegetation, to large-bodied, migratory, highly gregarious, grazing taxa inhabiting open savannahs and relying on flight or group-defense behaviors when facing predators. Here, we examined a geographically widespread clade within the Bovidae (the genus Gazella) that shows minimal interspecific variation in body size and asked if we could still uncover correlated changes of key ecological and behavioral traits during repeated transitions from open-land to mountain-dwelling. Our study used a multi-locus phylogeny (based on sequence variation of Cytb and six nuclear intron markers) of all extant members of the genus Gazella to infer evolutionary patterns of key ecological and behavioral traits and to estimate ancestral character states using Bayesian inference. At the base of the Gazella-phylogeny, open plains were inferred as the most likely habitat type, and three independent transitions toward mountain-dwelling were uncovered. Those shifts coincided with shifts from migratory to sedentary lifestyles. Character estimation for group size was largely congruent with movement patterns in that species forming large groups tended to be migratory, while small group size was correlated with a sedentary lifestyle. Evolutionary patterns of two other conspicuous traits (twinning ability vs. exclusive singleton births and hornless vs. horned females) did not follow this trend in the Gazella-phylogeny. Furthermore, we inferred the genus Gazella to have emerged in the Late Miocene to Pliocene (10-3Mya), and estimating ancestral ranges based on a Dispersal-Extinction-Cladogenesis model found the Middle East to be the most likely origin of the genus.

Importance of latrine communication in European rabbits shifts along a rural-to-urban gradient.

BACKGROUND: Information transfer in mammalian communication networks is often based on the deposition of excreta in latrines. Depending on the intended receiver(s), latrines are either formed at territorial boundaries (between-group communication) or in core areas of home ranges (within-group communication). The relative importance of both types of marking behavior should depend, amongst other factors, on population densities and social group sizes, which tend to differ between urban and rural wildlife populations. Our study is the first to assess (direct and indirect) anthropogenic influences on mammalian latrine-based communication networks along a rural-to-urban gradient in European rabbits (Oryctolagus cuniculus) living in urban, suburban and rural areas in and around Frankfurt am Main (Germany). RESULTS: The proportion of latrines located in close proximity to the burrow was higher at rural study sites compared to urban and suburban ones. At rural sites, we found the largest latrines and highest latrine densities close to the burrow, suggesting that core marking prevailed. By contrast, latrine dimensions and densities increased with increasing distance from the burrow in urban and suburban populations, suggesting a higher importance of peripheral marking. CONCLUSIONS: Increased population densities, but smaller social group sizes in urban rabbit populations may lead to an increased importance of between-group communication and thus, favor peripheral over core marking. Our study provides novel insights into the manifold ways by which man-made habitat alterations along a rural-to-urban gradient directly and indirectly affect wildlife populations, including latrine-based communication networks.

Cooling water of power plant creates "hot spots" for tropical fishes and parasites.

Thermally altered water bodies can function as "hot spots" where non-native species are establishing self-sustaining populations beyond their tropical and subtropical native regions. Whereas many tropical fish species have been found in these habitats, the introduction of non-native parasites often remains undetected. Here, n = 77 convict cichlids (Amatitlania nigrofasciata) were sampled by electro-fishing at two sites from a thermally altered stream in Germany and examined for parasite fauna and feeding ecology. Stomach content analysis suggests an opportunistic feeding strategy of A. nigrofasciata: while plant material dominated the diet at the warm water inlet (∼30 °C), relative contributions of insects, plants, and crustaceans were balanced 3 km downstream (∼27 °C). The most abundant non-native parasite species was the tropical nematode Camallanus cotti with P = 11.90 % and P = 80.00 % at the inlet and further downstream, respectively. Additionally, nematode larvae of Anguillicoloides crassus and one specimen of the subtropical species Bothriocephalus acheilognathi were isolated. A. nigrofasciata was also highly infected with the native parasite Acanthocephalus anguillae, which could be linked to high numbers of the parasite's intermediate host Asellus aquaticus. The aim of this study was to highlight the risk and consequences of the release and establishment of ornamental fish species for the introduction and spread of non-indigenous metazoan parasites using the convict cichlid as a model species. Furthermore, the spread of non-native parasites into adjacent fish communities needs to be addressed in the future as first evidence of Camallanus cotti in native fish species was also found.

miR-30-5p Regulates Muscle Differentiation and Alternative Splicing of Muscle-Related Genes by Targeting MBNL.

MicroRNAs (miRNAs), a class of single stranded, small (~22 nucleotides), non-coding RNAs, play an important role in muscle development. We focused on the role of the miR-30-5p family during bovine muscle development from previous high-throughput sequencing results and analyzed their expression profiles. MHC and MyoG mRNAs expression as well as their proteins were suppressed in differentiated C2C12 cells, suggesting the importance of miR-30-5p in muscle development. MBNL, the candidate target of miR-30-5p, is an alternative splicing regulation factor. MBNL1 and MBNL3 have opposite effects on muscle differentiation. Our results confirmed that miR-30a-5p and miR-30e-5p repress the expression of MBNL1, MBNL2 and MBNL3, whereas miR-30b-5p inhibits MBNL1 and MBNL2 expression. This provides direct evidence that MBNL expression can be flexibly regulated by miR-30-5p. Previous studies showed that MBNL1 promotes exon inclusion of two muscle-related genes (Trim55 and INSR). Through RNA splicing studies, we found that miR-30-5p had an effect on their alternative splicing, which means miR-30-5p via MBNL1 could be integrated into muscle signaling pathways in which INSR or Trim55 are located. In conclusion, miR-30-5p could inhibit muscle cell differentiation and regulate the alternative splicing of Trim55 and INSR by targeting MBNL. These results promote the understanding of the function of miRNAs in muscle development.

Unique evolutionary trajectories in repeated adaptation to hydrogen sulphide-toxic habitats of a neotropical fish (Poecilia mexicana).

Replicated ecological gradients are prime systems to study processes of molecular evolution underlying ecological divergence. Here, we investigated the repeated adaptation of the neotropical fish Poecilia mexicana to habitats containing toxic hydrogen sulphide (H2 S) and compared two population pairs of sulphide-adapted and ancestral fish by sequencing population pools of >200 individuals (Pool-Seq). We inferred the evolutionary processes shaping divergence and tested the hypothesis of increase of parallelism from SNPs to molecular pathways. Coalescence analyses showed that the divergence occurred in the face of substantial bidirectional gene flow. Population divergence involved many short, widely dispersed regions across the genome. Analyses of allele frequency spectra suggest that differentiation at most loci was driven by divergent selection, followed by a selection-mediated reduction of gene flow. Reconstructing allelic state changes suggested that selection acted mainly upon de novo mutations in the sulphide-adapted populations. Using a corrected Jaccard index to quantify parallel evolution, we found a negligible proportion of statistically significant parallel evolution of Jcorr  = 0.0032 at the level of SNPs, divergent genome regions (Jcorr  = 0.0061) and genes therein (Jcorr  = 0.0091). At the level of metabolic pathways, the overlap was Jcorr  = 0.2545, indicating increasing parallelism with increasing level of biological integration. The majority of pathways contained positively selected genes in both sulphide populations. Hence, adaptation to sulphidic habitats necessitated adjustments throughout the genome. The largely unique evolutionary trajectories may be explained by a high proportion of de novo mutations driving the divergence. Our findings favour Gould's view that evolution is often the unrepeatable result of stochastic events with highly contingent effects.

Associations between allelic polymorphism of the BMP Binding Endothelial Regulator and phenotypic variation of cattle.

The BMP Binding Endothelial Regulator (BMPER) is an inhibitor of bone morphogenetic proteins (BMPs), which play fundamental roles in adipocyte differentiation, fat development and energy balance. The objectives of this study were to detect polymorphisms of BMPER gene in four indigenous Chinese cattle populations and to investigate their effects on body size traits. Initially, three SNPs, namely G100597A (SNP1), C105331A (SNP2), and G105521A (SNP3) and eight distinct haplotypes were identified. In a total of 12 SNP-SNP combinations, SNP2-SNP3 had a strong linkage in Qinchuan cattle. These four cattle populations belong to intermediate genetic diversity at three SNP loci except Shuxuan cattle population in SNP3. At SNP1, genotype AA was associated with an increased body size. For SNP2, the heterozygous genotype individuals had a greater rump length than those of two other homozygotic genotypes. At SNP3, individuals with GG genotype had smaller rump length and hip width. A total of seven haplotype combinations were detected in Qinchuan cattle population and association analysis results showed individuals with Haplotype combination 4/2 (AAA/CAA) had greater rump length than those with Hap3/1 and Hap3/3 (P < 0.05). These results strongly suggest that bovine BMPER gene may be used as a genetic marker for cattle breeding.

Colonisation of toxic environments drives predictable life-history evolution in livebearing fishes (Poeciliidae).

New World livebearing fishes (family Poeciliidae) have repeatedly colonised toxic, hydrogen sulphide-rich waters across their natural distribution. Physiological considerations and life-history theory predict that these adverse conditions should favour the evolution of larger offspring. Here, we examined nine poeciliid species that independently colonised toxic environments, and show that these fishes have indeed repeatedly evolved much larger offspring size at birth in sulphidic waters, thus uncovering a widespread pattern of predictable evolution. However, a second pattern, only indirectly predicted by theory, proved additionally common: a reduction in the number of offspring carried per clutch (i.e. lower fecundity). Our analyses reveal that this secondary pattern represents a mere consequence of a classic life-history trade-off combined with strong selection on offspring size alone. With such strong natural selection in extreme environments, extremophile organisms may commonly exhibit multivariate phenotypic shifts even though not all diverging traits necessarily represent adaptations to the extreme conditions.