The gonadal transcriptome of the unisexual Amazon molly Poecilia formosa in comparison to its sexual ancestors, Poecilia mexicana and Poecilia latipinna

Genetic and phenotypic diversification in a widespread fish, the Sailfin Molly (Poecilia latipinna)

Widespread species often experience significant environmental clines over the area they naturally occupy. We investigated a widespread livebearing fish, the Sailfin molly (Poecilia latipinna) combining genetic, life-history, and environmental data, asking how structured populations are. Sailfin mollies can be found in coastal freshwater and brackish habitats from roughly Tampico, Veracruz in Mexico to Wilmington, North Carolina, in the USA. In addition, they are found inland on the Florida peninsula. Using microsatellite DNA, we genotyped 168 individuals from 18 populations covering most of the natural range of the Sailfin molly. We further determined standard life-history parameters for both males and females for these populations. Finally, we measured biotic and abiotic parameters in the field. We found six distinct genetic clusters based on microsatellite data, with very strong indication of isolation by distance. However, we also found significant numbers of migrants between adjacent populations. Despite genetic structuring we did not find evidence of cryptic speciation. The genetic clusters and the migration patterns do not match paleodrainages. Life histories vary between populations but not in a way that is easy to interpret. We suggest a role of humans in migration in the sailfin molly, for example in the form of a ship channel that connects southern Texas with Louisiana which might be a conduit for fish migration.

Hepatic transcriptome, transcriptional effects and antioxidant responses in Poecilia vivipara exposed to sanitary sewage

Aquatic environments are subject to threats from multiple human activities, particularly through the release of untreated sanitary sewage into the coastal environments. These effluents contain a large group of natural or synthetic compounds referred to as emerging contaminants. Monitoring the types and quantities of toxic substances in the environment, especially complex mixtures, is an exhausting and challenging task. Integrative effect-based tools, such as biomarkers, are recommended for environmental quality monitoring programs. In this study, fish Poecilia vivipara were exposed for 24 and 96 h to raw untreated sewage diluted 33 % (v/v) in order to identify hepatic genes to be used as molecular biomarkers. Through a de novo hepatic transcriptome assembly, using Illumina MiSeq, 54,285 sequences were assembled creating a reference transcriptome for this guppy species. Transcripts involved in biotransformation systems, antioxidant defenses, ABC transporters, nuclear and xenobiotic receptors were identified and evaluated by qPCR. Sanitary sewage induced transcriptional changes in AhR, PXR, CYP2K1, CYP3A30, NQO1, UGT1A1, GSTa3, GSTmu, ST1C1, SOD, ABCC1 and SOX9 genes from liver of fish, particularly after 96 h of exposure. Changes in hepatic enzyme activities were also observed. The enzymes showed differences in fish exposed to both periods, while in the gills there was a prevalence of significant results after 96 h. The observed differences were associated to gender and/or to sewage exposure. The obtained results support the use of P. vivipara as sentinel and model organism for ecotoxicological studies and evidence the importance of understanding the differential responses associated to gender.

Reproduction-associated pathways in females of gibel carp (Carassius gibelio) shed light on the molecular mechanisms of the coexistence of asexual and sexual reproduction

Gibel carp (Carassius gibelio) is a cyprinid fish that originated in eastern Eurasia and is considered as invasive in European freshwater ecosystems. The populations of gibel carp in Europe are mostly composed of asexually reproducing triploid females (i.e., reproducing by gynogenesis) and sexually reproducing diploid females and males. Although some cases of coexisting sexual and asexual reproductive forms are known in vertebrates, the molecular mechanisms maintaining such coexistence are still in question. Both reproduction modes are supposed to exhibit evolutionary and ecological advantages and disadvantages. To better understand the coexistence of these two reproduction strategies, we performed transcriptome profile analysis of gonad tissues (ovaries) and studied the differentially expressed reproduction-associated genes in sexual and asexual females. We used high-throughput RNA sequencing to generate transcriptomic profiles of gonadal tissues of triploid asexual females and males, diploid sexual males and females of gibel carp, as well as diploid individuals from two closely-related species, C. auratus and Cyprinus carpio. Using SNP clustering, we showed the close similarity of C. gibelio and C. auratus with a basal position of C. carpio to both Carassius species. Using transcriptome profile analyses, we showed that many genes and pathways are involved in both gynogenetic and sexual reproduction in C. gibelio; however, we also found that 1500 genes, including 100 genes involved in cell cycle control, meiosis, oogenesis, embryogenesis, fertilization, steroid hormone signaling, and biosynthesis were differently expressed in the ovaries of asexual and sexual females. We suggest that the overall downregulation of reproduction-associated pathways in asexual females, and their maintenance in sexual ones, allows the populations of C. gibelio to combine the evolutionary and ecological advantages of the two reproductive strategies. However, we showed that many sexual-reproduction-related genes are maintained and expressed in asexual females, suggesting that gynogenetic gibel carp retains the genetic toolkits for meiosis and sexual reproduction. These findings shed new light on the evolution of this asexual and sexual complex.

Achiasmatic meiosis in the unisexual Amazon molly, Poecilia formosa

Unisexual reproduction, which generates clonal offspring, is an alternative strategy to sexual breeding and occurs even in vertebrates. A wide range of non-sexual reproductive modes have been described, and one of the least understood questions is how such pathways emerged and how they mechanistically proceed. The Amazon molly, Poecilia formosa, needs sperm from males of related species to trigger the parthenogenetic development of diploid eggs. However, the mechanism, of how the unreduced female gametes are produced, remains unclear. Cytological analyses revealed that the chromosomes of primary oocytes initiate pachytene but do not proceed to bivalent formation and meiotic crossovers. Comparing ovary transcriptomes of P. formosa and its sexual parental species revealed expression levels of meiosis-specific genes deviating from P. mexicana but not from P. latipinna. Furthermore, several meiosis genes show biased expression towards one of the two alleles from the parental genomes. We infer from our data that in the Amazon molly diploid oocytes are generated by apomixis due to a failure in the synapsis of homologous chromosomes. The fact that this failure is not reflected in the differential expression of known meiosis genes suggests the underlying molecular mechanism may be dysregulation on the protein level or misexpression of a so far unknown meiosis gene, and/or hybrid dysgenesis because of compromised interaction of proteins from diverged genomes.

Patterns of gene expression in ovaries of sexual vs. asexual lineages of a freshwater snail

Why sexual reproduction is so common when asexual reproduction should be much more efficient and less costly remains an open question in evolutionary biology. Comparisons between otherwise similar sexual and asexual taxa allow us to characterize the genetic architecture underlying asexuality, which can, in turn, illuminate how this reproductive mode transition occurred and the mechanisms by which it is maintained or disrupted. Here, we used transcriptome sequencing to compare patterns of ovarian gene expression between actively reproducing obligately sexual and obligately asexual females from multiple lineages of Potamopyrgus antipodarum, a freshwater New Zealand snail characterized by frequent separate transitions to asexuality and coexistence of otherwise similar sexual and asexual lineages. We also used these sequence data to evaluate whether population history accounts for variation in patterns of gene expression. We found that source population was a major source of gene expression variation, and likely more influential than reproductive mode. This outcome for these common garden-raised snails is strikingly similar to earlier results from field-collected snails. While we did not identify a likely set of candidate genes from expression profiles that could plausibly explain how transitions to asexuality occurred, we identified around 1,000 genes with evidence of differential expression between sexual and asexual reproductive modes, and 21 genes that appear to exhibit consistent expression differences between sexuals and asexuals across genetic backgrounds. This second smaller set of genes provides a good starting point for further exploration regarding a potential role in the transition to asexual reproduction. These results mark the first effort to characterize the causes of asexuality in P. antipodarum, demonstrate the apparently high heritability of gene expression patterns in this species, and hint that for P. antipodarum, transitions to asexuality might not necessarily be strongly associated with broad changes in gene expression.

Premeiotic endoreplication is essential for obligate parthenogenesis in geckos

Obligate parthenogenesis evolved in reptiles convergently several times, mainly through interspecific hybridization. The obligate parthenogenetic complexes typically include both diploid and triploid lineages. Offspring of parthenogenetic hybrids are genetic copies of their mother; however, the cellular mechanism enabling the production of unreduced cells is largely unknown. Here, we show that oocytes go through meiosis in three widespread, or even strongly invasive, obligate parthenogenetic complexes of geckos, namely in diploid and triploid Lepidodactylus lugubris, and triploid Hemiphyllodactylus typus and Heteronotia binoei. In all four lineages, the majority of oocytes enter the pachytene at the original ploidy level, but their chromosomes cannot pair properly and instead form univalents, bivalents and multivalents. Unreduced eggs with clonally inherited genomes are formed from germ cells that had undergone premeiotic endoreplication, in which appropriate segregation is ensured by the formation of bivalents made from copies of identical chromosomes. We conclude that the induction of premeiotic endoreplication in reptiles was independently co-opted at least four times as an essential component of parthenogenetic reproduction and that this mechanism enables the emergence of fertile polyploid lineages within parthenogenetic complexes.

Skipping sex: A nonrecombinant genomic assemblage of complementary reproductive modules

The unusual occurrence and developmental diversity of asexual eukaryotes remain a puzzle. De novo formation of a functioning asexual genome requires a unique assembly of sets of genes or gene states to disrupt cellular mechanisms of meiosis and gametogenesis, and to affect discrete components of sexuality and produce clonal or hemiclonal offspring. We highlight two usually overlooked but essential conditions to understand the molecular nature of clonal organisms, that is, a nonrecombinant genomic assemblage retaining modifiers of the sexual program, and a complementation between altered reproductive components. These subtle conditions are the basis for physiologically viable and genetically balanced transitions between generations. Genomic and developmental evidence from asexual animals and plants indicates the lack of complementation of molecular changes in the sexual reproductive program is likely the main cause of asexuals' rarity, and can provide an explanatory frame for the developmental diversity and lability of developmental patterns in some asexuals as well as for the discordant time to extinction estimations.

Preliminary Identification of Candidate Genes Related to Survival of Gynogenetic Rainbow Trout (Oncorhynchus mykiss) Based on Comparative Transcriptome Analysis

In the present research, the eggs from four rainbow trout females were used to provide four groups of gynogenetic doubled haploids (DHs). The quality of the eggs from different clutches was comparable, however, interclutch differences were observed in the gynogenetic variants of the experiment and the survival of DH specimens from different groups varied from 3% to 57% during embryogenesis. Transcriptome analysis of the eggs from different females exhibited inter-individual differences in the maternal genes' expression. Eggs originating from females whose gynogenetic offspring had the highest survival showed an increased expression of 46 genes when compared to the eggs from three other females. Eggs with the highest survival of gynogenetic embryos showed an up-regulation of genes that are associated with cell survival, migration and differentiation (tyrosine-protein kinase receptor TYRO3-like gene), triglyceride metabolism (carnitine O-palmitoyltransferase 1 gene), biosynthesis of polyunsaturated fat (3-oxoacyl-acyl-carrier-protein reductase gene), early embryogenic development (protein argonaute-3 gene, leucine-rich repeat-containing protein 3-like gene), 5S RNA binding (ribosome biogenesis regulatory protein homolog) as well as senescence and aging (telomerase reverse transcriptase, TERT gene), among others. Positive correlation between the genotypic efficiency and egg transcriptome profiles indicated that at least some of the differentially expressed genes should be considered as potential candidate genes for the efficiency of gynogenesis in rainbow trout.

Using asexual vertebrates to study genome evolution and animal physiology: Banded (Fundulus diaphanus) x Common Killifish (F. heteroclitus) hybrid lineages as a model system

Wild, asexual, vertebrate hybrids have many characteristics that make them good model systems for studying how genomes evolve and epigenetic modifications influence animal physiology. In particular, the formation of asexual hybrid lineages is a form of reproductive incompatibility, but we know little about the genetic and genomic mechanisms by which this mode of reproductive isolation proceeds in animals. Asexual lineages also provide researchers with the ability to produce genetically identical individuals, enabling the study of autonomous epigenetic modifications without the confounds of genetic variation. Here, we briefly review the cellular and molecular mechanisms leading to asexual reproduction in vertebrates and the known genetic and epigenetic consequences of the loss of sex. We then specifically discuss what is known about asexual lineages of Fundulus diaphanus x F. heteroclitus to highlight gaps in our knowledge of the biology of these clones. Our preliminary studies of F. diaphanus and F. heteroclitus karyotypes from Porter's Lake (Nova Scotia, Canada) agree with data from other populations, suggesting a conserved interspecific chromosomal arrangement. In addition, genetic analyses suggest that: (a) the same major clonal lineage (Clone A) of F. diaphanus x F. heteroclitus has remained dominant over the past decade, (b) some minor clones have also persisted, (c) new clones may have recently formed, and iv) wild clones still mainly descend from F. diaphanus ♀ x F. heteroclitus ♂ crosses (96% in 2017-2018). These data suggest that clone formation may be a relatively rare, but continuous process, and there are persistent environmental or genetic factors causing a bias in cross direction. We end by describing our current research on the genomic causes and consequences of a transition to asexuality and the potential physiological consequences of epigenetic variation.

Are you more than the sum of your parents' genes? Phenotypic plasticity in a clonal vertebrate and F1 hybrids of its parental species

All known vertebrate clones have originated from hybridization events and some have produced distinct evolutionary lineages via hybrid speciation. Amazon mollies (Poecilia formosa) present an excellent study system to investigate how clonal species have adapted to heterogeneous environments because they are the product of a single hybridization event between male sailfin mollies (Poecilia latipinna) and female Atlantic mollies (Poecilia mexicana). Here, we ask whether the hybrid species differs from the combination of its parental species' genes in its plastic response to different environments. Using a three-way factorial design, we exposed neonates produced by Amazon mollies and reciprocal F1 hybrid crosses to different thermal (24°C and 29°C) and salinity (0/2, 12, and 20 ppt) regimes. We measured various ontogenetic and life history characteristics across the life span of females. Our major results were as follows: (1) Reaction norms of growth and maturation to temperature and salinity are quite similar between the two hybrid crosses; (2) Amazon molly reaction norms were qualitatively different than the P. latipinna male and P. mexicana female (L×M) hybrids for the ontogenetic variables; (3) Amazon molly reaction norms in reproductive traits were also quite different from L×M hybrids; and (4) The reaction norms of net fertility were very different between Amazon mollies and L×M hybrids. We conclude that best locale for Amazon mollies is not the best locale for hybrids, which suggests that Amazon mollies are not just an unmodified mix of parental genes but instead have adapted to the variable environments in which they are found. Hybridization resulting in asexuality may represent an underappreciated mechanism of speciation because the unlikely events required to produce such hybrids rarely occur and is dependent upon the genetic distance between parental species.

Placing the hybrid origin of the asexual Amazon molly (Poecilia formosa) based on historical climate data

Asexual hybrids are important model organisms for addressing questions in evolution and ecology, especially for understanding the role of hybridization in speciation. They are rare in nature and several hypotheses have been suggested to explain this. We use an asexual fish, the Amazon molly (Poecilia formosa), to establish the area in which it was formed via hybridization 125 000 years ago. Using species distribution models and climate models for the Last Interglacial (LIG) we found that model projections to the LIG show a similar map to the present climate model and parental species potentially overlapped in a relatively small area near Tampico, Mexico. This makes P. formosa one of a few hybrid species for which we know the parental species, the time of hybridization, and likely the place of hybridization. Based on the small area of overlap, our data is in agreement with the idea that asexual hybrids may be rare not because they are evolutionary dead ends but are formed rarely.

The Legacy of Sexual Ancestors in Phenotypic Variability, Gene Expression, and Homoeolog Regulation of Asexual Hybrids and Polyploids

Hybridization and polyploidization are important evolutionary processes whose impacts range from the alteration of gene expression and phenotypic variation to the triggering of asexual reproduction. We investigated fishes of the Cobitis taenia-elongatoides hybrid complex, which allowed us to disentangle the direct effects of both processes, due to the co-occurrence of parental species with their diploid and triploid hybrids. Employing morphological, ecological, and RNAseq approaches, we investigated the molecular determinants of hybrid and polyploid forms. In contrast with other studies, hybridization and polyploidy induced relatively very little transgressivity. Instead, Cobitis hybrids appeared intermediate with a clear effect of genomic dosing when triploids expressed higher similarity to the parent contributing two genome sets. This dosage effect was symmetric in the germline (oocyte gene expression), interestingly though, we observed an overall bias toward C. taenia in somatic tissues and traits. At the level of individual genes, expression-level dominance vastly prevailed over additivity or transgressivity. Also, trans-regulation of gene expression was less efficient in diploid hybrids than in triploids, where the expression modulation of homoeologs derived from the "haploid" parent was stronger than those derived from the "diploid" parent. Our findings suggest that the apparent intermediacy of hybrid phenotypes results from the combination of individual genes with dominant expression rather than from simple additivity. The efficiency of cross-talk between trans-regulatory elements further appears dosage dependent. Important effects of polyploidization may thus stem from changes in relative concentrations of trans-regulatory elements and their binding sites between hybridizing genomes. Links between gene regulation and asexuality are discussed.