Microbiome and epigenetic variation in wild fish with low genetic diversity.

Non-genetic sources of phenotypic variation, such as the epigenome and the microbiome, could be important contributors to adaptive variation for species with low genetic diversity. However, little is known about the complex interaction between these factors and the genetic diversity of the host, particularly in wild populations. Here, we examine the skin microbiome composition of two closely-related mangrove killifish species with different mating systems (self-fertilising and outcrossing) under sympatric and allopatric conditions. This allows us to partition the influence of the genotype and the environment on their microbiome and (previously described) epigenetic profiles. We find the diversity and community composition of the skin microbiome are strongly shaped by the environment and, to a lesser extent, by species-specific influences. Heterozygosity and microbiome alpha diversity, but not epigenetic variation, are associated with the fluctuating asymmetry of traits related to performance (vision) and behaviour (aggression). Our study identifies that a proportion of the epigenetic diversity and microbiome differentiation is unrelated to genetic variation, and we find evidence for an associative relationship between microbiome and epigenetic diversity in these wild populations. This suggests that both mechanisms could potentially contribute to variation in species with low genetic diversity.

The complete genome sequence of the protogynous hermaphrodite North American annual killifish Millerichthys robustus (Rivulidae, Cyprinodontiformes).

The Mexican Rivulus, Millerichthys robustus, is an enigmatic species of seasonal killifish endemic of the Southeast Mexico that has changed paradigms on the evolution of annualism in killifishes. This species survives in ephemeral environments that experience a period of seasonal drought that causes the death of all adult fish. However, populations persist due to their drought-resistant embryos capable of arresting their development in diapause until the next the rainy season. Millerichthys evolved unique characteristics within annual killifishes as a functional sequential hermaphroditism, in which females change to males (protogynous) under perceived conditions of mate competition. Also, M. robustus express different color phenotypes in both sexes: five-color phenotypes continuously distributed in various perceptual units between yellow and red in males, and different number of ocelli disposed in caudal peduncle in females. The phylogenetic relationships of M. robustus revealed that it is a sister clade to two non-annual species found exclusively in Cuba (Rivulus cylindraceus and Rivulus berovidesi), indicating that the annual life cycle, through the acquisition of embryonic diapause, has evolved independently in this species. Here, we present the complete genome sequences for the North American annual killifish Millerichthys robustus. The raw data and assembled genome are available in GeneBank.

Influence of the mutation load on the genomic composition of hybrids between outcrossing and self-fertilizing species.

Hybridization is a natural process whereby two diverging evolutionary lineages reproduce and create offspring of mixed ancestry. Differences in mating systems (e.g., self-fertilization and outcrossing) are expected to affect the direction and extent of hybridization and introgression in hybrid zones. Among other factors, selfers and outcrossers are expected to differ in their mutation loads. This has been studied both theoretically and empirically; however, conflicting predictions have been made on the effects mutation loads of parental species with different mating systems can have on the genomic composition of hybrids. Here, we develop a multi-locus, selective model to study how the different mutation load built up in selfers and outcrossers as a result of selective interference and homozygosity impact the long-term genetic composition of hybrid populations. Notably, our results emphasize that genes from the parental population with lesser mutation load get rapidly overrepresented in hybrid genomes, regardless of the hybrids own mating system. When recombination tends to be more important than mutation, outcrossers' genomes tend to be of higher quality and prevail. When recombination rates are low, however, selfers' genomes may reach higher quality than outcrossers' genomes and prevail in the hybrids. Taken together, these results provide concrete insights into one of the multiple factors influencing hybrid genome ancestry and introgression patterns in hybrid zones containing species with different mating systems.

The Complete Genome Sequences of three species from the killifish genus Kryptolebias (Rivulidae, Cyprinodontiformes).

The killifish genus Kryptolebias currently contains seven recognized species found in freshwater and mangrove microhabitats in South and Central America, the Caribbean, and Florida. Kryptolebias species have several unique features. Beyond its amphibious nature, two of the Kryptolebias species (K. marmoratus and K. hermaphroditus sensu Costa 2011) are the only known vertebrates capable of self-fertilization. Although reference genomes for self-fertilizing species are readily available, a complete understanding of the genomic basis leading to the transition in mating systems in the genus requires more genomic resources at a broader taxonomic level. Here, we present the complete genome sequences for three Kryptolebias species from Brazil: the endangered freshwater species Kryptolebias brasiliensis and Kryptolebias gracilis and the androdiecious but obligate outcrossing Kryptolebias ocellatus (sensu Costa 2011). The raw data and assembled genomes are available in GenBank.

Additive and non-additive epigenetic signatures of natural hybridization between fish species with different mating systems.

Hybridization is a major source of evolutionary innovation. In plants, epigenetic mechanisms can help to stabilize hybrid genomes and contribute to reproductive isolation, but the relationship between genetic and epigenetic changes in animal hybrids is unclear. We analysed the relationship between genetic background and methylation patterns in natural hybrids of two genetically divergent fish species with different mating systems, Kryptolebias hermaphroditus (self-fertilizing) and K. ocellatus (outcrossing). Co-existing parental species displayed highly distinct genetic (SNPs) and methylation patterns (37,000 differentially methylated cytosines). Hybrids had predominantly intermediate methylation patterns (88.5% of the sites) suggesting additive effects, as expected from hybridization between genetically distant species. The large number of differentially methylated cytosines between hybrids and parental species (n = 5,800) suggests that hybridization may play a role in increasing genetic and epigenetic variation. Although most of the observed epigenetic variation was additive and had a strong genetic component, we also found a small percentage of non-additive, potentially stochastic, methylation differences that might act as an evolutionary bet-hedging strategy and increase fitness under environmental instability.

Phylogenomics reveals extensive introgression and a case of mito-nuclear discordance in the killifish genus Kryptolebias.

Introgression is a widespread evolutionary process leading to phylogenetic inconsistencies among distinct parts of the genomes, particularly between mitochondrial and nuclear-based phylogenetic reconstructions (e.g., mito-nuclear discordances). Here, we used mtDNA and genome-wide nuclear sites to provide the first phylogenomic-based hypothesis on the evolutionary relationships within the killifish genus Kryptolebias. In addition, we tested for evidence of past introgression in the genus given the multiple reports of undergoing hybridization between its members. Our mtDNA phylogeny generally agreed with the relationships previously proposed for the genus. However, our reconstruction based on nuclear DNA revealed an unknown lineage - Kryptolebias sp. 'ESP' - as the sister group of the self-fertilizing mangrove killifishes, K. marmoratus and K. hermaphroditus. All individuals sequenced of Kryptolebias sp. 'ESP' had the same mtDNA haplotype commonly observed in K. hermaphroditus, demonstrating a clear case of mito-nuclear discordance. Our analysis further confirmed extensive history of introgression between Kryptolebias sp. 'ESP' and K. hermaphroditus. Population genomics analyses indicate no current gene flow between the two lineages, despite their current sympatry and history of introgression. We also confirmed introgression between other species pairs in the genus that have been recently reported to form hybrid zones. Overall, our study provides a phylogenomic reconstruction covering most of the Kryptolebias species, reveals a new lineage hidden in a case of mito-nuclear discordance, and provides evidence of multiple events of ancestral introgression in the genus. These findings underscore the importance of investigating different genomic information in a phylogenetic framework, particularly in taxa where introgression is common as in the sexually diverse mangrove killifishes.

Against the Odds: Hybrid Zones between Mangrove Killifish Species with Different Mating Systems.

Different mating systems are expected to affect the extent and direction of hybridization. Due to the different levels of sexual conflict, the weak inbreeder/strong outbreeder (WISO) hypothesis predicts that gametes from self-incompatible (SI) species should outcompete gametes from self-compatible (SC) ones. However, other factors such as timing of selfing and unilateral incompatibilities may also play a role on the direction of hybridization. In addition, differential mating opportunities provided by different mating systems are also expected to affect the direction of introgression in hybrid zones involving outcrossers and selfers. Here, we explored these hypotheses with a unique case of recent hybridization between two mangrove killifish species with different mating systems, Kryptolebias ocellatus (obligately outcrossing) and K. hermaphroditus (predominantly self-fertilizing) in two hybrid zones in southeast Brazil. Hybridization rates were relatively high (~20%), representing the first example of natural hybridization between species with different mating systems in vertebrates. All F1 individuals were sired by the selfing species. Backcrossing was small, but mostly asymmetrical with the SI parental species, suggesting pattern commonly observed in plant hybrid zones with different mating systems. Our findings shed light on how contrasting mating systems may affect the direction and extent of gene flow between sympatric species, ultimately affecting the evolution and maintenance of hybrid zones.

Filling the gaps: phylogeography of the self-fertilizing Kryptolebias species (Cyprinodontiformes: Rivulidae) along South American mangroves.

Mangrove killifishes of the genus Kryptolebias have been historically classified as rare because of their small size and cryptic nature. Major gaps in distribution knowledge across mangrove areas, particularly in South America, challenge the understanding of the taxonomic status, biogeographical patterns and genetic structuring of the lineages composing the self-fertilizing "Kryptolebias marmoratus species complex." In this study, the authors combined a literature survey, fieldwork and molecular data to fill major gaps of information about the distribution of mangrove killifishes across western Atlantic mangroves. They found that selfing mangrove killifishes are ubiquitously distributed across the Caribbean, Central and South American mangroves and report 14 new locations in South America, extending the range of both the "Central clade" and "Southern clade" lineages which overlap in the Amazon. Although substantial genetic differences were found between clades, the authors also found further genetic structuring within clades, with populations in Central America, north and northeast Brazil generally showing higher levels of genetic diversity compared to the clonal ones in southeast Brazil. The authors discuss the taxonomic status and update the geographical distribution of the Central and Southern clades, as well as potential dispersal routes and biogeographical barriers influencing the distribution of the selfing mangrove killifishes in the western Atlantic mangroves.

More than meets the eye: syntopic and morphologically similar mangrove killifish species show different mating systems and patterns of genetic structure along the Brazilian coast.

Different mating systems can strongly affect the extent of genetic diversity and population structure among species. Given the increased effects of genetic drift on reduced population size, theory predicts that species undergoing self-fertilisation should have greater population structure than outcrossed species; however, demographic dynamics may affect this scenario. The mangrove killifish clade is composed of the two only known examples of self-fertilising species among vertebrates (Kryptolebias marmoratus and Kryptolebias hermaphroditus). A third species in this clade, Kryptolebias ocellatus, inhabits mangrove forests in southeast Brazil; however, its mating system and patterns of genetic structure have been rarely explored. Here, we examined the genetic structure and phylogeographic patterns of K. ocellatus along its distribution, using mitochondrial DNA and microsatellites to compare its patterns of genetic structure with the predominantly selfing and often-syntopic, K. hermaphroditus. Our results indicate that K. ocellatus reproduces mainly by outcrossing, with no current evidence of selfing, despite being an androdioecious species. Our results also reveal a stronger population subdivision in K. ocellatus compared to K. hermaphroditus, contrary to the theoretical predictions based on reproductive biology of the two species. Our findings indicate that, although morphologically similar, K. ocellatus and K. hermaphroditus had remarkably different evolutionary histories when colonising the same mangrove areas in southeastern Brazil, with other factors (e.g., time of colonisation, dispersal/establishment capacity) having more profound effects on the current population structuring of those species than differences in mating systems.

Environmental enrichment induces intergenerational behavioural and epigenetic effects on fish.

Parental effects influence offspring phenotypes through pre- and post-natal routes but little is known about their molecular basis, and therefore their adaptive significance. Epigenetic modifications, which control gene expression without changes in the DNA sequence and are influenced by the environment, may contribute to parental effects. We investigated the effects of environmental enrichment on the behaviour, metabolic rate and brain DNA methylation patterns of parents and offspring of the highly inbreed mangrove killifish (Kryptolebias marmoratus). Parental fish reared in enriched environments had lower cortisol levels, lower metabolic rates and were more active and neophobic than those reared in barren environments. They also differed in 1,854 methylated cytosines (DMCs). Offspring activity and neophobia were determined by the parental environment. Among the DMCs of the parents, 98 followed the same methylation patterns in the offspring, three of which were significantly influenced by parental environments irrespective of their own rearing environment. Our results suggest that parental environment influences the behaviour and, to some extent, the brain DNA methylation patterns of the offspring.

Local parasite pressures and host genotype modulate epigenetic diversity in a mixed-mating fish.

Parasite-mediated selection is one of the main drivers of genetic variation in natural populations. The persistence of long-term self-fertilization, however, challenges the notion that low genetic variation and inbreeding compromise the host's ability to respond to pathogens. DNA methylation represents a potential mechanism for generating additional adaptive variation under low genetic diversity. We compared genetic diversity (microsatellites and AFLPs), variation in DNA methylation (MS-AFLPs), and parasite loads in three populations of Kryptolebias hermaphroditus, a predomintanly self-fertilizing fish, to analyze the potential adaptive value of DNA methylation in relation to genetic diversity and parasite loads. We found strong genetic population structuring, as well as differences in parasite loads and methylation levels among sampling sites and selfing lineages. Globally, the interaction between parasites and inbreeding with selfing lineages influenced DNA methylation, but parasites seemed more important in determining methylation levels at the local scale.

Contrasting DNA methylation responses of inbred fish lines to different rearing environments.

Epigenetic mechanisms generate plastic phenotypes that can become locally adapted across environments. Disentangling genomic from epigenomic variation is challenging in sexual species due to genetic variation among individuals, but it is easier in self-fertilizing species. We analysed DNA methylation patterns of two highly inbred strains of a naturally self-fertilizing fish reared in two contrasting environments to investigate the obligatory (genotype-dependent), facilitated (partially depend on the genotype) or pure (genotype-independent) nature of the epigenetic variation. We found higher methylation differentiation between genotypes than between environments. Most methylation differences between environments common to both strains followed a pattern where the two genotypes (inbred lines) responded to the same environmental context with contrasting DNA methylation levels (facilitated epialleles). Our findings suggest that, at least in part, DNA methylation could depend on the dynamic interaction between the genotype and the environment, which could explain the plasticity of epigenetically mediated phenotypes.

Updated checklist and DNA barcode-based species delimitations reveal taxonomic uncertainties among freshwater fishes from the mid-north-eastern Caatinga ecoregion, north-eastern Brazil.

The mid-north-eastern Caatinga is a semiarid freshwater ecoregion in North-eastern Brazil that is dominated by temporary rivers and is currently classified as one of the least ichthyologically-known ecoregions in the world. The present study aimed to provide an updated checklist of mid-north-eastern Caatinga ecoregion (MNCE) freshwater fish species and evaluate their taxonomic identity using morphology, DNA barcoding and multiple species delimitation approaches. After reviewing published studies and ichthyological collections, 119 species were identified. Among these were 94 putatively valid native and 14 non-native species, five undescribed native species, four new records for the MNCE, 11 potential cases of misidentification and 14 species listed as inquirenda. Additionally, 252 individuals from 49 species were barcoded, revealing three potential taxonomic synonyms. The combined molecular approaches estimated a total of 91 native species, although a finalized species list for the MNCE awaits additional taxonomic revisions and field surveys. This study provides the most up-to-date species checklist for the MNCE and a molecular reference database for identifying MNCE fishes with DNA barcodes. Results highlight the need to integrate traditional taxonomy with molecular approaches to correctly identify species, especially in taxonomically problematic ecoregions such as the MNCE.

Ultrasensitive environmental assessment of xeno-estrogens in water samples using label-free graphene immunosensors.

There is a growing interest in the possible environmental health impact posed by endocrine-disrupting chemicals (EDCs). A challenge to the field of endocrine disruption is that these substances are diverse and may not appear to share any structural similarity other than usually being low molecular mass (<1000 Da) compounds. Here we demonstrate the effectiveness of sensor device for the detection of low molecular weight, poorly water soluble, estrogenic compounds E1, E2 and EE2, fabricated by electropolymerization over graphene screen printed electrode (SPE). The PANI/Gr-SPE-devices displayed linear responses to estrogenic substances, in EIS assays, from 0.0975 ng/L to 200 ng/L in water samples, with a detection limit of 0.043 pg/L for E1, 0.19 ng/L for E2 and 0.070 pg/L for EE2 which is lower than other current biosensing techniques. This portable, disposable immunosensor offers a solution for immediate measurement at sample collection sites, due to its excellent sensitivity and selectivity when testing water samples obtained directly from rivers and waste water treatment facilities. The simple screen printing production method will enable the low cost, high volume production required for this type of environmental analysis.

Headwater Capture Evidenced by Paleo-Rivers Reconstruction and Population Genetic Structure of the Armored Catfish (Pareiorhaphis garbei) in the Serra do Mar Mountains of Southeastern Brazil.

Paleo-drainage connections and headwater stream-captures are two main historical processes shaping the distribution of strictly freshwater fishes. Recently, bathymetric-based methods of paleo-drainage reconstruction have opened new possibilities to investigate how these processes have shaped the genetic structure of freshwater organisms. In this context, the present study used paleo-drainage reconstructions and single-locus cluster delimitation analyses to examine genetic structure on the whole distribution of Pareiorhaphis garbei, a 'near threatened' armored catfish from the Fluminense freshwater ecoregion in Southeastern Brazil. Sequences of two mitochondrial genes (cytochrome b and cytochrome c oxidase subunit 1) were obtained from five sampling sites in four coastal drainages: Macaé (KAE), São João (SJO), Guapi-Macacu [sub-basins Guapiaçu (GAC) and Guapimirim (GMI)], and Santo Aleixo (SAL). Pronounced genetic structure was found, involving 10 haplotypes for cytB and 6 for coi, with no haplotypes shared between localities. Coalescent-based delineation methods as well as distance-based methods revealed genetic clusters corresponding to each sample site. Paleo-drainage reconstructions showed two putative paleo-rivers: an eastern one connecting KAE and SJO; and a western one merging in the Guanabara Bay (GAC, GMI, and SAL). A disagreement was uncovered between the inferred past riverine connections and current population genetic structure. Although KAE and SJO belong to the same paleo-river, the latter is more closely related to specimens from the Guanabara paleo-river. This discordance between paleo-drainage connections and phylogenetic structure may indicate an ancient stream-capture event in headwaters of this region. Furthermore, all analyses showed high divergence between KAE and the other lineages, suggesting at least one cryptic species in the latter, and that the nominal species should be restricted to the Macaé river basin, its type locality. In this drainage, impacts such as the invasive species and habitat loss can be especially threatening for such species with a narrow range. Our results also suggest that freshwater fishes from headwaters in the Serra do Mar mountains might have different biogeographical patterns than those from the lowlands, indicating a complex and dynamic climatic and geomorphological history.

First record of a male of Kryptolebias hermaphroditus Costa, 2011 (Cyprinodontiformes: Cynolebiidae)

During an ichthyological survey in September 2015 at the Ceará-Mirim River estuary, Rio Grande do Norte State, northeastern Brazil, we collected a male of Kryptolebias hermaphroditus , a cynolebiid species that had been previously described as containing exclusively self-fertilizing hermaphrodites. This is the first record of a male in this species, over 140 years after the discovery of the mangrove rivulid species from Brazil. Our discovery reinforces the need for more studies in K. hermaphroditus , as well as the potential of this species as a model for evolutionary studies due to its unique mating system.(AU)

Durante uma amostragem ictiológica em setembro de 2015 no estuário do rio Ceará-Mirim, Estado do Rio Grande do Norte, Nordeste do Brasil, nós coletamos um macho de Kryptolebias hermaphroditus , uma espécie de rivulídeo que foi descrita como contendo apenas hermafroditas auto-fertilizantes. Este é o primeiro registro de um macho dessa espécie, mais de 140 anos depois da descoberta das espécies de rivulídeos de manguezais do Brasil. Nossa descoberta reforça a necessidade de mais estudos em K . hermaphroditus , assim como, o potencial desta espécie como um modelo para estudos evolutivos devido ao seu sistema reprodutivo único.(AU)