Scale morphology and phylogenetic character mapping of scale-surface microstructures in sixteen Aphanius species (Teleostei: Aphaniidae)

Recent studies on the analysis of the evolutionary relationships between the various taxa of the teleost fishes are based on the evaluation of the efficacy of different approaches, verifying the taxonomical congruence between the hypotheses inferred from molecular data and those from morphological data. The first objective of this study was to investigate the macro- and microscopic characterization of the body key scales for the sixteen Aphanius species, and to discuss their taxonomic significance. In addition, to reveal the phylogenetic origins of particular adaptations and to identify a possible cases of convergent evolution in scales of Aphanius species, we used a phylogenetic character mapping procedure that uses an independent data source, namely molecular data that is presumed to be functionally independent of the scale traits. Comparison of the resulted phylogenetic trees and scale data dendrogram implied that neither scale shape/scale surface morphology, nor scale surface microstructure can be used for the species identification and inferring phylogenetic relationships in the genus Aphanius. This conclusion confirms previous studies, indicating that the scale shape and ornamentations varies with habitat, age and size of the fish. We also concluded that the polygonal state in the scales of Aphanius species was a primitive shared character state.

Genotoxicity by long-term exposure to the auxinic herbicides 2,4-dichlorophenoxyacetic acid and dicamba on Cnesterodon decemmaculatus (Pisces: Poeciliidae)

Long-term genotoxic effects of two auxinic herbicide formulations, namely, the 58.4% 2,4-dichlorophenoxyacetic acid (2,4-D)-based DMA^® and the 57.7% dicamba (DIC)-based Banvel^® were evaluated on Cnesterodon decemmaculatus. Primary DNA lesions were analyzed by the single-cell gel electrophoresis methodology. Two sublethal concentrations were tested for each herbicide corresponding to 2.5% and 5% of the LC50_96h values. Accordingly, fish were exposed to 25.2 and 50.4 mg/L or 41 and 82 mg/L for 2,4-D and DIC, respectively. Fish were continuously exposed for 28 days with replacement of test solutions every 3 days. Genotoxicity was evaluated in ten individuals from each experimental point at the beginning of the exposure period (0 day) and at 7, 14, 21 and 28 days thereafter. Results demonstrated for first time that 2,4-D-based formulation DMA^® induced primary DNA strand breaks after 7-28 days exposure on C. decemmaculatus regardless its concentration. On the other hand, DIC-based formulation Banvel^® exerted its genotoxic effect after exposure during 7-14 days and 7 days of 2.5 and 5% LC50_96h, respectively. The present study represents the first evidence of primary DNA lesions induced by two widely employed auxinic herbicides on C. decemmaculatus, namely 2,4-D and DIC, following long-term exposure.

Characterization and expression of Xiphophorus maculatus microsatellite Msb069 full sequence in subgenus Poecilia

Msb069 primer pairs encompassed region is believed to be associated with a quantitative trait loci (QTL) of dorsal fin length in subgenus Poecilia. However, detailed investigation on Msb069 which originated from Xiphophorus on subgenus Poecilia remains unexplored. In this study, full sequence of Msb069 was characterized by sequencing bioinformatics analysis and gene expression. The sequence analysis of Msb069 primer pairs encompassed region on three species of Poecilia revealed higher number of microsatellite tandem repeats in Poecilia latipinna (ATG₁₆) compared to P. sphenops (ATG_13-14). There is no notable pattern of ATGtandem repeats discovered in the hybrids. The full sequence of Msb069 is 734 bp in length and showed a 233 bp conserved region between Xiphophorus and Poecilia. BLAST search performed on this sequence revealed no significant similarities. Nonquantitative RT-PCR exhibited the presence of Msb069 transcripts in three different tissues in subgenus Poecilia. Meanwhile, quantitative RTPCR expression on two different tissues showed relatively higher expression of Msb069 transcript in P. latipinna dorsal fin tissues in both male and female fishes, suggesting a repressive function of this transcript with respect to dorsal fin length. However the exact gene expression event of Msb069 is still unknown and requires further investigation.

Natural and sexual selection drive multivariate phenotypic divergence along climatic gradients in an invasive fish

Invasive species that rapidly spread throughout novel distribution ranges are prime models to investigate climate-driven phenotypic diversification on a contemporary scale. Previous studies on adaptive diversification along latitudinal gradients in fish have mainly considered body size and reported either increased or decreased body size towards higher latitudes (i.e. Bergmann's rule). Our study is the first to investigate phenotypic divergence in multiple traits, including sexually selected traits (size and shape of the male copulatory organ, the gonopodium) of invasive Gambusia affinis in China. We studied body size, life history traits and morphological variation across populations spanning 17 degrees of latitude and 16 degrees of longitude. Even though we found phenotypic variation along climatic gradients to be strongest in naturally selected traits, some sexually selected traits also showed systematic gradual divergence. For example, males from southern populations possessed wider gonopodia with increased armament. Generally, males and females diverged in response to different components of climatic gradients (latitudinal or longitudinal variation) and in different trait suites. We discuss that not only temperature regimes, but also indirect effects of increased resource and mate competition (as a function of different extrinsic overwinter mortality rates) alter the selective landscape along climatic gradients.

Comparison of Xiphophorus and human melanoma transcriptomes reveals conserved pathway interactions

Comparative analysis of human and animal model melanomas can uncover conserved pathways and genetic changes that are relevant for the biology of cancer cells. Spontaneous melanoma in Xiphophorus interspecies backcross hybrid progeny may be informative in identifying genes and functional pathways that are similarly related to melanoma development in all vertebrates, including humans. To assess functional pathways involved in the Xiphophorus melanoma, we performed gene expression profiling of the melanomas produced in interspecies BC1 and successive backcross generations (i.e., BC5 ) of the cross: X. hellerii × [X. maculatus Jp 163 A × X. hellerii]. Using RNA-Seq, we identified genes that are transcriptionally co-expressed with the driver oncogene, xmrk. We determined functional pathways in the fish melanoma that are also present in human melanoma cohorts that may be related to dedifferentiation based on the expression levels of pigmentation genes. Shared pathways between human and Xiphophorus melanomas are related to inflammation, cell migration, cell proliferation, pigmentation, cancer development, and metastasis. Our results suggest xmrk co-expressed genes are associated with dedifferentiation and highlight these signaling pathways as playing important roles in melanomagenesis.

Multigene phylogeny supports diversification of four-eyed fishes and one-sided livebearers (Cyprinodontiformes: Anablepidae) related to major South American geological events

The high diversity of Neotropical fishes has been attributed to major South American palaeogeographic events, such as Andean uplift, rise of the Isthmus of Panama and marine transgressions. However, the unavailability of temporal information about evolution and diversification of some fish groups prevents the establishment of robust hypotheses about correlations between species diversification and proposed palaeogeographical events. One example is the Anablepidae, a family of teleost fishes found mostly in coastal habitats of Central and South America, but also in some inner river basins of South America. Historical aspects of the distribution patterns of the Anablepidae were never analysed and no accurate estimation of time of its origin and diversification is presently available. A multi-gene analysis was performed to estimate Anablepidae phylogenetic position, age and biogeography, comprising seven nuclear genes. The suborder Cyprinodontoidei was recovered in three major clades, one comprising all the Old World Cyprinodontoidei and two comprising New World lineages. Anablepidae was recovered as the sister group of the New World Poeciliidae, with the Amazonian genus Fluviphylax as their sister group. The ages found for the origin and diversification of Cyprinodontiformes were congruent with the pattern recorded for other vertebrate groups, with an origin anterior to the Cretaceous-Paleogene (K-Pg) transition and diversification during the Paleogene. The age estimated for the split between the Atlantic and Pacific lineages of Anableps was congruent with the rise of Panamanian Isthmus. The results suggest Miocene marine transgressions as determinant to the current distribution of Jenynsia.

Response of PXR signaling pathway to simvastatin exposure in mosquitofish (Gambusia affinis) and its histological changes

As a widely used lipid lowering agent, simvastatin recently has been frequently detected in aquatic environment and the potential adverse effects from simvastatin exposure to non-target organisms such as fish is worthy of more attention. The aim of this study was to reveal the responses of detoxification system in fish to simvastatin exposure. In this investigation a ubiquitous small freshwater fish, mosquito fish (Gambusia affinis), was employed as test organism, and the transcriptional expression of nucleus transcriptional factor pregnane X receptor (PXR) and its downstream genes, including P-glycoprotein (P-gp), cytochrome 3A (CYP3A), multidrug resistance protein 2 (MRP2), UDP-glucuronosyl transferase (UGT) in mosquito fish were investigated by qRT-PCR methods under the exposure of concentrations of simvastatin (0.5 μg L-1, 5 μg L-1, 50 μgL-1, 500 μg L-1) for 24 h, 72 h and 168 h. The related enzyme activity (Erythromycin-N-Demethylase, ERND), the protein expression of PXR and the histological changes of liver tissues in fish were also determined via west blotting and transmission electron microscope approaches in the same conditions. Results showed that the mRNA expression of PXR, CYP3A and P-gp showed significantly changes under simvastatin exposure, exhibiting an obvious time/dose-effect relationship with the prolong of exposure time. ERND activity also showed time-effect at 24 h, and western blotting showed PXR protein displaying a dose-effect relationship to some extent. Hepatocyte cellular of mosquito fish exposed to simvastatin (5 μg L-1, 168 h) exhibited obvious histological changes in form of swelling, incomplete fragmentary structure etc. Overall, simvastatin altered the expression of PXR signaling pathway and subsequently bring about changes in high-levels of mosquito fish.

The transcriptional response of skin to fluorescent light exposure in viviparous (Xiphophorus) and oviparous (Danio, Oryzias) fishes

Differences in light sources are common in animal facilities and potentially can impact experimental results. Here, the potential impact of lighting differences on skin transcriptomes has been tested in three aquatic animal models commonly utilized in biomedical research, (Xiphophorus maculatus (platyfish), Oryzias latipes (medaka) and Danio rerio (zebrafish). Analysis of replicate comparative RNA-Seq data showed the transcriptional response to commonly utilized 4100K or "cool white" fluorescent light (FL) is much greater in platyfish and medaka than in zebrafish. FL induces genes associated with inflammatory and immune responses in both medaka and zebrafish; however, the platyfish exhibit suppression of genes involved with immune/inflammation, as well as genes associated with cell cycle progression. Furthermore, comparative analyses of gene expression data from platyfish UVB exposures, with medaka and zebrafish after exposure to 4100K FL, show comparable effects on the same stress pathways. We suggest the response to light is conserved, but that long-term adaptation to species specific environmental niches has resulted in a shifting of the wavelengths required to incite similar "genetic" responses in skin. We forward the hypothesis that the "genetic perception" of light may have evolved differently than ocular perception and suggest that light type (i.e., wavelengths emitted) is an important parameter to consider in experimental design.

Characterization of basal gene expression trends over a diurnal cycle in Xiphophorus maculatus skin, brain and liver

Evolutionarily conserved diurnal circadian mechanisms maintain oscillating patterns of gene expression based on the day-night cycle. Xiphophorus fish have been used to evaluate transcriptional responses after exposure to various light sources and it was determined that each source incites distinct genetic responses in skin tissue. However, basal expression levels of genes that show oscillating expression patterns in day-night cycle, may affect the outcomes of such experiments, since basal gene expression levels at each point in the circadian path may influence the profile of identified light responsive genes. Lack of knowledge regarding diurnal fluctuations in basal gene expression patterns may confound the understanding of genetic responses to external stimuli (e.g., light) since the dynamic nature of gene expression implies animals subjected to stimuli at different times may be at very different stages within the continuum of genetic homeostasis. We assessed basal gene expression changes over a 24-hour period in 200 select Xiphophorus gene targets known to transcriptionally respond to various types of light exposure. We identified 22 genes in skin, 36 genes in brain and 28 genes in liver that exhibit basal oscillation of expression patterns. These genes, including known circadian regulators, produced the expected expression patterns over a 24-hour cycle when compared to circadian regulatory genes identified in other species, especially human and other vertebrate animal models. Our results suggest the regulatory network governing diurnal oscillating gene expression is similar between Xiphophorus and other vertebrates for the three Xiphophorus organs tested. In addition, we were able to categorize light responsive gene sets in Xiphophorus that do, and do not, exhibit circadian based oscillating expression patterns.

Exposure to 4100K fluorescent light elicits sex specific transcriptional responses in Xiphophorus maculatus skin

It has been reported that exposure to artificial light may affect oxygen intake, heart rate, absorption of vitamins and minerals, and behavioral responses in humans. We have reported specific gene expression responses in the skin of Xiphophorus fish after exposure to ultraviolet light (UV), as well as, both broad spectrum and narrow waveband visible light. In regard to fluorescent light (FL), we have shown that male X. maculatus exposed to 4100K FL (i.e. "cool white") rapidly suppress transcription of many genes involved with DNA replication and repair, chromosomal segregation, and cell cycle progression in skin. We have also detailed sex specific transcriptional responses of Xiphophorus skin after exposure to UVB. However, investigation of gender differences in global gene expression response after exposure to 4100K FL has not been reported, despite common use of this FL source for residential, commercial, and animal facility illumination. Here, we compare RNA-Seq results analyzed to assess changes in the global transcription profiles of female and male X. maculatus skin in response to 4100K FL exposure. Our results suggest 4100K FL exposure incites a sex-biased genetic response including up-modulation of inflammation in females and down modulation of DNA repair/replication in males. In addition, we identify clusters of genes that become oppositely modulated in males and females after FL exposure that are principally involved in cell death and cell proliferation.

A High-Quality Reference Genome for the Invasive Mosquitofish Gambusia affinis Using a Chicago Library

The western mosquitofish, Gambusia affinis, is a freshwater poecilid fish native to the southeastern United States but with a global distribution due to widespread human introduction. Gambusia affinis has been used as a model species for a broad range of evolutionary and ecological studies. We sequenced the genome of a male G. affinis to facilitate genetic studies in diverse fields including invasion biology and comparative genetics. We generated Illumina short read data from paired-end libraries and in vitro proximity-ligation libraries. We obtained 54.9× coverage, N50 contig length of 17.6 kb, and N50 scaffold length of 6.65 Mb. Compared to two other species in the Poeciliidae family, G. affinis has slightly fewer genes that have shorter total, exon, and intron length on average. Using a set of universal single-copy orthologs in fish genomes, we found 95.5% of these genes were complete in the G. affinis assembly. The number of transposable elements in the G. affinis assembly is similar to those of closely related species. The high-quality genome sequence and annotations we report will be valuable resources for scientists to map the genetic architecture of traits of interest in this species.

Diclofenac exposure alter the expression of PXR and its downstream target genes in mosquito fish (Gambusia affinis)

As one of widely used drugs, Diclofenac (DCF) recently has been universally detected in aquatic environment and some negative effects derived from DCF exposure to mammals have been also reported. However, studies about its potential deleterious effects on non-target organisms like fish still require more investigation. In this study an ubiquitous small freshwater invader species in Southern of China, mosquito fish (Gambusia affinis), was employed as test organism. We firstly cloned the crucial partial sequences of nucleus transcriptional factor related genes pregnane X receptor (PXR) and its downstream genes, including P-glycoprotein (P-gp), cytochrome 3A (CYP3A), multidrug resistance protein 2 (MRP2), glutathione peroxidase (GPx) and thioredoxin reductase (TXR) in mosquito fish. The phylogenetic trees of PXR, CYP3A and MRP2 were constructed based on their deduced amino acids sequences, respectively. Phylogenetic trees and blast results showed a high similarity between G. affinis and other killifish species, such as Xiphophorus maculatus. The transcriptional expression of these genes mentioned above and partly related enzymes/proteins activities were then measured under the exposure of environmentally relevant concentrations of DCF (from 0.5μgL^-1 to 500μgL^-1) for 24h and 168h. Results showed that the mRNA expression of PXR, CYP3A, P-gp and TXR showed dramatic induction under DCF exposure, exhibiting an obvious time-effect relationship with the extend of exposure time. In terms of enzyme activity and protein content, no dramatic changes as in transcription were observed. Western blotting showed PXR protein increased at 24h but decreased at 168h with the increasing of DCF concentration, displaying a dose-effect relationship to some extent. GPX activity was continuously induced both at 24h and 168h, exhibiting a good consistency with the performance of GPX gene. GSSH/T-GSH increased in all treatments. Overall, DCF had traceable effects on the expression of PXR and its downstream target genes in mosquito fish.

Maternal-by-environment but not genotype-by-environment interactions in a fish without parental care

The impact of environmental conditions on the expression of genetic variance and on maternal effects variance remains an important question in evolutionary quantitative genetics. We investigate here the effects of early environment on variation in seven adult life history, morphological, and secondary sexual traits (including sperm characteristics) in a viviparous poeciliid fish, the mosquitofish Gambusia holbrooki. Specifically, we manipulated food availability during early development and then assessed additive genetic and maternal effects contributions to the overall phenotypic variance in adults. We found higher heritability for female than male traits, but maternal effects variance for traits in both sexes. An interaction between maternal effects variance and rearing environment affected two adult traits (female age at maturity and male size at maturity), but there was no evidence of trade-offs in maternal effects across environments. Our results illustrate (i) the potential for pre-natal maternal effects to interact with offspring environment during development, potentially affecting traits through to adulthood and (ii) that genotype-by-environment interactions might be overestimated if maternal-by-environment interactions are not accounted for, similar to heritability being overestimated if maternal effects are ignored. We also discuss the potential for dominance genetic variance to contribute to the estimate of maternal effects variance.

Physiological responses and gene expression changes in the western mosquitofish (Gambusia affinis) exposed to progesterone at environmentally relevant concentrations

Progesterone (P4) is a natural and synthetic steroid, widely distributed in the aquatic environments. It can lead to adverse effects on the endocrine system in aquatic organisms. This study investigated the toxicological effects of exposure to environmentally relevant concentrations (4, 44, and 410ng/L) of progesterone for 42 d on adult female mosquitofish, Gambusia affinis. We performed morphological and histological analyses on gonads, anal fins, liver, and gills after the exposure of mosquito fish to P4. The expression levels of genes (vtg, er, and ar isoforms) related to fish reproduction and detoxification (cyp1a) in the liver were quantified by quantitative real-time polymerase chain reaction. The results showed that the progesterone exposure induced slight masculinization in female mosquitofish, influenced the oocyte maturation as revealed by histology of the ovaries, and caused severe damages to the liver and gills of adult female mosquitofish. It also suppressed the mRNAs expression of vtg, er, cyp1a, and significantly enhanced the expression of ar mRNA in the liver. This study reveals the molecular and physiological effects of progesterone at environmentally relevant concentrations, which might further be translated to alterations in the reproduction of mosquitofish.

Let's jump in: A phylogenetic study of the great basin springfishes and poolfishes, Crenichthys and Empetrichthys (Cyprinodontiformes: Goodeidae)

North America’s Great Basin has long been of interest to biologists due to its high level of organismal endemicity throughout its endorheic watersheds. One example of such a group is the subfamily Empetricthyinae. In this paper, we analyzed the relationships of the Empetrichtyinae and assessed the validity of the subspecies designations given by Williams and Wilde within the group using concatenated phylogenetic tree estimation and species tree estimation. Samples from 19 populations were included covering the entire distribution of the three extant species of Empetricthyinae–Crenichthys nevadae, Crenichthys baileyi and Empetricthys latos. Three nuclear introns (S8 intron 4, S7 intron 1, and P0 intron 1) and one mitochondrial gene (Cytb) were sequenced for phylogenetic analysis. Using these sequences, we generated two separate hypotheses of the evolutionary relationships of Empetrichtyinae- one based on the mitochondrial data and one based on the nuclear data using Bayesian phylogenetics. Haplotype networks were also generated to look at the relationships of the populations within Empetrichthyinae. After comparing the two phylogenetic hypotheses, species trees were generated using *BEAST with the nuclear data to further test the validity of the subspecies within Empetrichthyinae. The mitochondrial analyses supported four lineages within C. baileyi and 2 within C. nevadae. The concatenated nuclear tree was more conserved, supporting one clade and an unresolved polytomy in both species. The species tree analysis supported the presence of two species within both C. baileyi and C. nevadae. Based on the results of these analyses, the subspecies designations of Williams and Wilde are not valid, rather a conservative approach suggests there are two species within C. nevadae and two species within C. baileyi. No structure was found for E. latos or the populations of Empetricthyinae. This study represents one of many demonstrating the invalidity of subspecies and their detriment to species identification, conservation, and understanding.

Adverse effects of BDE-47 on in vivo developmental parameters, thyroid hormones, and expression of hypothalamus-pituitary-thyroid (HPT) axis genes in larvae of the self-fertilizing fish Kryptolebias marmoratus

2,2',4,4'-tetrabromodiphenylether (BDE-47) is known to have the potential to disrupt the thyroid endocrine system in fishes due to its structural similarity to the thyroid hormones triiodothyronine (T3) and thyroxine (T4). However, the effects of BDE-47 on thyroid function in fishes remain unclear. In this study, abnormal development (e.g. deformity, hemorrhaging) and an imbalance in thyroid hormone (TH) homeostasis was shown in the early developmental stages of the mangrove killifish Kryptolebias marmoratus in response to BDE-47 exposure. To examine the thyroid endocrinal effect of BDE-47 exposure in mangrove killifish K. marmoratus larvae, transcript levels of genes involved in TH homeostasis and hypothalamus-pituitary-thyroid (HPT) axis-related genes were measured. The expression of thyroid hormone metabolism-related genes (e.g. deiodinases, UGT1ab) and HPT axis-related genes was up-regulated and there were significant changes in TH levels (P < 0.05) in response to BDE-47 exposure. This study provides insights into the regulation of TH homeostasis at the transcriptional level and provides a better understanding on the potential impacts of BDE-47 on the thyroid endocrine system of fishes.

The origin and biogeographic diversification of fishes in the family Poeciliidae

The fish subfamily Poeciliinae (sensu Parenti, 1981) is widely distributed across the Western Hemisphere and a dominant component of the fish communities of Central America. Poeciliids have figured prominently in previous studies on the roles of dispersal and vicariance in shaping current geographic distributions. Most recently, Hrbek et al. combined a DNA-based phylogeny of the family with geological models to provide a biogeographic perspective that emphasized the role of both vicariance and dispersal. Here we expand on that effort with a database enlarged in the quantity of sequence represented per species, in the number of species included, and in an enlarged and more balanced representation of the order Cyprinodontiformes. We combine a robust timetree based upon multiple fossil calibrations with enhanced biogeographic analyses that include ancestral area reconstructions to provide a detailed biogeographic history of this clade. Key features of our results are that the family originated in South America, but its major diversification dates to a later colonization of Central America. We also resolve additional colonizations among South, Central and North America and the Caribbean and consider how this reconstruction contributes to our understanding of the mechanisms of dispersal.

Partial Gene Sequencing of CYP1A, Vitellogenin, and Metallothionein in Mosquitofish Gambusia yucatana and Gambusia sexradiata

Ground characteristics in the Yucatan Peninsula make recovery and treatment of wastewater very expensive. This situation has contributed to an increase of pollutants in the aquifer. Unfortunately, studies related to the effects of those pollutants in native organisms are scarce. The aim of this work was to obtain partial sequences of widely known genes used as biomarkers of pollutant effect in Gambusia yucatana and Gambusia sexradiata. The studied genes were: cytochrome P450 1A (CYP1A); vitellogenin (VTG); metallothionein (MT), and two housekeeping genes, 18S and β-actin. From reported sequences of Gambusia affinis, primers were designed and amplification was done in the local Gambusia species exposed for 48 h to gasoline (100 µL/L, stirred for 24 h pre-exposure). Preliminary results revealed partial sequences of all genes with an approximate average length of 200 bp. BLAST analysis of found sequences indicated a minimum of 97% identity with reported sequences for G. affinis or Gambusia holbrooki showing great similarity.

Transcriptomic and physiological changes in Eastern Mosquitofish (Gambusia holbrooki) after exposure to progestins and anti-progestagens

Endocrine active compounds (EACs) remain an important group of chemicals that require additional evaluation to determine their environmental impacts. While estrogens and androgens were previously demonstrated to impact organisms during environmental exposures, progestagens have recently been shown to have strong impacts on aquatic organisms. To gain an understanding of the impacts of these types of chemicals on aquatic species, experiments evaluating the mechanisms of action of progestagen exposure were conducted with the Eastern Mosquitofish (Gambusia holbrooki). The objective of this study was to conduct hepatic microarray analysis of male and female G. holbrooki exposed to progestins and anti-progestagens. In addition, we evaluated the ability of levonorgestrel, a synthetic progesterone (progestin), to induce anal fin elongation and to determine how anal fin growth is modulated during co-exposures with progesterone and androgen receptor antagonists. Gene expression analyses were conducted on male and female G. holbrooki exposed for 48h to the agonist levonorgestrel, the antagonist mifepristone, or a mixture of the two chemicals. Microarray analysis revealed that mifepristone does not act as an anti-progestagen in G. holbrooki in liver tissues, and that levonorgestrel elicits strong effects on the processes of embryo development and lipid transport. Levonorgestrel was also demonstrated to induce male secondary sexual characteristic formation in females, and co-exposure of either an androgen or levonorgestrel in the presence of the anti-androgen flutamide prevented anal fin elongation. These results provide indications as to the potential impacts of progestins, including non-target effects such as secondary sexual characteristic formation, and demonstrate the importance of this class of chemicals on aquatic organisms.

Toxic and genotoxic effects of the 2,4-dichlorophenoxyacetic acid (2,4-D)-based herbicide on the Neotropical fish Cnesterodon decemmaculatus

Acute toxicity and genotoxicity of the 54.8% 2,4-D-based commercial herbicide DMA® were assayed on Cnesterodon decemmaculatus (Pisces, Poeciliidae). Whereas lethal effect was used as the end point for mortality, frequency of micronuclei (MNs), other nuclear abnormalities and primary DNA damage evaluated by the single cell gel electrophoresis (SCGE) assay were employed as end points for genotoxicity. Mortality studies demonstrated an LC50 96 h value of 1008 mg/L (range, 929-1070) of 2,4-D. Behavioral changes, e.g., gathering at the bottom of the aquarium, slowness in motion, slow reaction and abnormal swimming were observed. Exposure to 2,4-D within the 252-756 mg/L range increased the frequency of MNs in fish exposed for both 48 and 96 h. Whereas blebbed nuclei were induced in treatments lasting for 48 and 96 h, notched nuclei were only induced in fish exposed for 96 h. Regardless of both concentration and exposure time, 2,4-D did not induce lobed nuclei and binucleated erythrocytes. In addition, we found that exposure to 2,4-D within the 252-756 mg/L range increased the genetic damage index in treatments lasting for either 48 and 96 h. The results represent the first experimental evidence of the lethal and several sublethal effects, including behavioral alterations and two genotoxic properties namely the induction of MNs and primary DNA strand breaks, exerted by 2,4-D on an endemic organism as C. decemmaculatus.

Fitness consequences of artificial selection on relative male genital size

Male genitalia often show remarkable differences among related species in size, shape and complexity. Across poeciliid fishes, the elongated fin (gonopodium) that males use to inseminate females ranges from 18 to 53% of body length. Relative genital size therefore varies greatly among species. In contrast, there is often tight within-species allometric scaling, which suggests strong selection against genital-body size combinations that deviate from a species' natural line of allometry. We tested this constraint by artificially selecting on the allometric intercept, creating lines of males with relatively longer or shorter gonopodia than occur naturally for a given body size in mosquitofish, Gambusia holbrooki. We show that relative genital length is heritable and diverged 7.6-8.9% between our up-selected and down-selected lines, with correlated changes in body shape. However, deviation from the natural line of allometry does not affect male success in assays of attractiveness, swimming performance and, crucially, reproductive success (paternity).

Conservation of Hox gene clusters in the self-fertilizing fish Kryptolebias marmoratus (Cyprinodontiformes; Rivulidae)

In this study, whole Hox gene clusters in the self-fertilizing mangrove killifish Kryptolebias marmoratus (Cyprinodontiformes; Rivulidae), a unique hermaphroditic vertebrate in which both sex organs are functional at the same time, were identified from whole genome and transcriptome sequences. The aim was to increase the understanding of the evolutionary status of conservation of this Hox gene cluster across fish species.

Responses of the CYP1A biomarker in Jenynsia multidentata and Phalloceros caudimaculatus and evaluation of a CYP1A refractory phenotype

The level of cytochrome P450 1A (CYP1A) in fish is used as a typical environmental biomarker for the presence of organic contaminants. We used RT-qPCR to investigate CYP1A mRNA levels in the liver, gill and gonopodium of guppies Jenynsia multidentata and Phalloceros caudimaculatus in wetlands within the Rio Grande city (RG) which is under the influence of the Patos Lagoon Estuary (RS, Brazil). The CYP1A mRNA levels evaluated in fish liver from two locations that receive non-treated wastewater effluents (S3 and S4) and another locations near an oil refinery (S6) and an industrial complex (S7), were higher than in locations remote from those sites (S1, S2 and S5). The sum of 16 priority PAHs in sediment confirmed high levels in S4 and S6 (3914.0 and 4414.0 ng g(-1) dw, respectively) comparing to S7>S2>S3>S5>S1 (119.3, 66.3, 62.8, 16.4 and 1.7 ng g(-1) dw). J. multidentata from sites S1 to S4 that were transferred to the laboratory exhibited CYP1A induction after 24 h waterborne exposure to 1 µM betanaphtoflavone (BNF) in all organs compared to controls, except in the liver of fish from site S4. This lack of CYP1A induction by BNF indicates a CYP1A refractory phenotype in guppy. Although this characteristic possibly involves the alteration in AHR signaling or control, the mechanism of resistance is unknown. The present study provides information about the use of the use of CYP1A in South American guppies as an useful biomarker tool for environmental contamination studies.

X. couchianus and X. hellerii genome models provide genomic variation insight among Xiphophorus species

BACKGROUND

Xiphophorus fishes are represented by 26 live-bearing species of tropical fish that express many attributes (e.g., viviparity, genetic and phenotypic variation, ecological adaptation, varied sexual developmental mechanisms, ability to produce fertile interspecies hybrids) that have made attractive research models for over 85 years. Use of various interspecies hybrids to investigate the genetics underlying spontaneous and induced tumorigenesis has resulted in the development and maintenance of pedigreed Xiphophorus lines specifically bred for research. The recent availability of the X. maculatus reference genome assembly now provides unprecedented opportunities for novel and exciting comparative research studies among Xiphophorus species.

RESULTS

We present sequencing, assembly and annotation of two new genomes representing Xiphophorus couchianus and Xiphophorus hellerii. The final X. couchianus and X. hellerii assemblies have total sizes of 708 Mb and 734 Mb and correspond to 98 % and 102 % of the X. maculatus Jp 163 A genome size, respectively. The rates of single nucleotide change range from 1 per 52 bp to 1 per 69 bp among the three genomes and the impact of putatively damaging variants are presented. In addition, a survey of transposable elements allowed us to deduce an ancestral TE landscape, uncovered potential active TEs and document a recent burst of TEs during evolution of this genus.

CONCLUSIONS

Two new Xiphophorus genomes and their corresponding transcriptomes were efficiently assembled, the former using a novel guided assembly approach. Three assembled genome sequences within this single vertebrate order of new world live-bearing fishes will accelerate our understanding of relationship between environmental adaptation and genome evolution. In addition, these genome resources provide capability to determine allele specific gene regulation among interspecies hybrids produced by crossing any of the three species that are known to produce progeny predisposed to tumor development.

Transcriptional control analyses of the Xiphophorus melanoma oncogene

Melanoma development in interspecific hybrids of Xiphophorus is induced by the overexpression of the mutationally activated receptor tyrosine kinase Xmrk in pigment cells. Based on the melanocyte specificity of the transcriptional upregulation, a pigment cell-specific promoter region was postulated for xmrk, the activity of which is controlled in healthy purebred fish by the molecularly still unidentified regulator locus R. However, as yet the xmrk promoter region is still poorly characterized. In order to contribute to a better understanding of xmrk expression regulation, we performed a functional analysis of the entire putative gene regulatory region of the oncogene using conventional plasmid-based reporter systems as well as a newly established method employing BAC-derived luciferase reporter constructs in melanoma and non-melanoma cell lines. Using the melanocyte-specific mitfa promoter as control, we could demonstrate that our in vitro system is able to reliably monitor regulation of transcription through cell type-specific regulatory sequences. We found that sequences within 200kb flanking the xmrk oncogene do not lead to any specific transcriptional activation in melanoma compared to control cells. Hence, xmrk reporter constructs fail to faithfully reproduce the endogenous transcriptional regulation of the oncogene. Our data therefore strongly indicate that the melanocyte-specific transcription of xmrk is not the consequence of pigment cell-specific cis-regulatory elements in the promoter region. This hints at additional regulatory mechanisms involved in transcriptional control of the oncogene, thereby suggesting a key role for epigenetic mechanisms in oncogenic xmrk overexpression and thereby in tumor development in Xiphophorus.

Sex-specific molecular genetic response to UVB exposure in Xiphophorus maculatus skin

In both Xiphophorus fishes and humans, males are reported to have a higher incidence of melanoma than females. To better understand sex-specific differences in the molecular genetic response to UVB, we performed RNA-Seq experiments in skin of female and male Xiphophorus maculatus Jp 163 B following UVB doses of 8 or 16kJ/m(2) exposure. Male X. maculatus differentially express a significantly larger number of transcripts following exposure to 16kJ/m(2) UVB (1293 genes) compared to 8kJ/m(2) UVB (324 genes). Female skin showed differential gene expression in a larger number of transcripts following 8kJ/m(2) UVB (765) than did males; however, both females and males showed similar numbers of differentially expressed genes at 16kJ/m(2) UVB (1167 and1293, respectively). Although most modulated transcripts after UVB exposure represented the same dominant pathways in both females and males (e.g., DNA repair, circadian rhythm, and fatty acid biosynthesis), we identified genes in several pathways that exhibited opposite modulation in female vs. male skin (e.g., synaptic development, cell differentiation, wound healing, and glucose metabolism). The oppositely modulated genes appear related through uncoupling protein 3 (UCP3) that is involved with the regulation of fatty acid oxidation and serves to balance glucose and lipid metabolism. Overall, these results identify gender-specific differences in UVB-induced genetic profiles in the skin of females and males and show female and male X. maculatus respond to UVB differently through pathways involved in reactive oxygen species, wound healing, and energy homeostasis.

Exposure to fluorescent light triggers down regulation of genes involved with mitotic progression in Xiphophorus skin

We report RNA-Seq results from skin of X. maculatus Jp 163 B after exposure to various doses of "cool white" fluorescent light (FL). We show that FL exposure incites a genetic transcriptional response in skin nearly as great as observed for UVB exposure; however, the gene sets modulated due to exposure to the two light sources are quite different. Known light responsive genes involved in maintaining circadian cycling (e.g., clock, cry2a, cry1b, per1b, per2, per3, and arntl1a) exhibited expected shifts in transcriptional expression upon FL exposure. Exposure to FL also resulted in down-regulated transcription of many genes involved with cell cycle progression (e.g., cdc20, cdc45, cdca7b, plk1, cdk1, ccnb-3, and cdca7a) and chromosome segregation (e.g., cenpe, cenpf, cenpi, cenpk, cenpo, cenpp, and cenpu; cep70; knstrm, kntc, mcm2, mcm5; smc2). In addition, several DNA replication and recombination repair genes (e.g., pola1, pole, rec52, rad54l, rpa1, and parpbp) exhibit reduced expression in FL exposed X. maculatus skin. Some genes down modulated by FL are known to be associated with DNA repair and human diseases (e.g., atm2, brip1, fanc1, fancl, and xrcc4). The overall suppression of genes involved with mitotic progression in the skin of adult fish is consistent with entry into the light phase of the circadian cycle. Current efforts are aimed at determining specific wavelengths that may lead to differential expression among the many genes affected by fluorescent light exposure.

Molecular genetic response to varied wavelengths of light in Xiphophorus maculatus skin

Xiphophorus fishes represent a model often utilized to study UVB induced tumorigenesis. Recently, varied genetic responses to UVB exposure have been documented in the skin of female and male Xiphophorus, as have differences in UVB response in the skin of different parental species and for interspecies hybrids produced from crossing them. Additionally, it has been shown that exposure to "cool white" fluorescent light induces a shift in the genetic profiles of Xiphophorus skin that is nearly as robust as the UVB response, but involves a fundamentally different set of genes. Given these results and the use of Xiphophorus interspecies hybrids as an experimental model for UVB inducible melanoma, it is of interest to characterize genes that may be transcriptionally modulated in a wavelength specific manner. The global molecular genetic response of skin upon exposure of the intact animal to specific wavelengths of light has not been investigated. Herein, we report results of RNA-Seq experiments from the skin of male Xiphophorus maculatus Jp 163 B following exposure to varied 50nm wavelengths of light ranging from 300-600nm. We identify two specific wavelength regions, 350-400nm (88 genes) and 500-550nm (276 genes), that exhibit transcriptional modulation of a significantly greater number of transcripts than any of the other 50nm regions in the 300-600nm range. Observed functional sets of genes modulated within these two transcriptionally active light regions suggest different mechanisms of gene modulation.

Molecular genetic response of Xiphophorus maculatus-X. couchianus interspecies hybrid skin to UVB exposure

The phenotypic and genetic similarities between Xiphophorus and human melanoma render Xiphophorus a useful animal model for studying the genetic basis of melanoma etiology. In the Xiphophorus model, melanoma has been shown to be inducible by ultraviolet light (UVB) exposure among interspecies hybrids, but not in parental line fish similarly treated. This leads to questions of what genes are responsive to UVB exposure in the skin of the interspecies hybrids, as well as how parental alleles in hybrids may be differentially regulated and the potential roles they may play in induced melanomagenesis. To address these questions, we produced X. maculatus Jp 163 B×X. couchianus (Sp-Couch) F1 hybrid fish, exposed both hybrid and parental fish to UVB, and performed gene expression profiling of the skin using RNA-Seq methodology. We characterized a group of unique UVB-responsive genes in Sp-Couch hybrid including dct, pmela, tyr, tyrp1a, slc2a11b, rab38a, rab27, tspan10, slc45a2, oca2, slc24a5, ptn and mitfa. These genes are associated with melanin production and melanocyte proliferation. They were also up-regulated in Sp-Couch hybrid, indicating that their UVB response is hybridization initiated. In the hybrid, several melanin production and pigmentation related genes, including slc45a2, tspan10, dct, slc2a11b and ptn showed either X. couchianus or X. maculatus allele specific expression. The finding that these genes exhibit allele specific expression regulatory mechanisms in Sp-Couch hybrids, but do not exhibit a corresponding UVB response in either one of the parental fishes, may suggest UVB targets and imply mechanisms regarding the susceptibility of Sp-Couch to induced melanomagenesis.

Effects of norfloxacin on hepatic genes expression of P450 isoforms (CYP1A and CYP3A), GST and P-glycoprotein (P-gp) in Swordtail fish (Xiphophorus Helleri)

The presence of antibiotics including norfloxacin in the aquatic environment may cause adverse effects in non-target organisms. But the toxic mechanisms of fluoroquinolone to fish species are still not completely elucidated. Thus, it is essential to investigate the response of fish to the exposure of fluoroquinolone at molecular or cellular level for better and earlier prediction of these environmental pollutants toxicity. The sub-chronic toxic effects of norfloxacin (NOR) on swordtail fish (Xiphophoru s helleri) were investigated by measuring mRNA expression of cytochrome P450 1A (CYP1A), cytochrome P450 3A (CYP3A), glutathione S-transferase (GST) and P-glycoprotein (P-gp) and their corresponding enzyme activities (including ethoxyresorufin O-deethylase, erythromycin N-demethylase and GST. Results showed that NOR significantly affected the expression of CYP1A, CYP3A, GST and P-gp genes in swordtails. The gene expressions were more responsive to NOR exposure than their corresponding enzyme activities. Moreover, sexual differences were found in gene expression and enzyme activities of swordtails exposed to NOR. Females displayed more dramatic changes than males. The study further demonstrated that the combined biochemical and molecular parameters were considered as useful biomarkers to improve our understanding of potential ecotoxicological risks of NOR exposure to aquatic organisms.

Mixture effects between different azoles and β-naphthoflavone on the CYP1A biomarker in a fish cell line

The cytochrome P450 1A (CYP1A) biomarker response was studied in the Poeciliopsis lucida hepatocellular carcinoma (PLHC-1) cell line, which represents a good model for studies on aryl hydrocarbon receptor (AhR) - CYP1A signaling. The PLHC-1 cells were exposed to the prototypical CYP1A inducer and AhR agonist β-naphthoflavone (BNF) in combination with different azoles. Two imidazoles (clotrimazole and prochloraz) and two benzimidazoles (nocodazole and omeprazole) were used. Exposure to clotrimazole, prochloraz and nocodazole resulted in 2-4 fold induction of the CYP1A-mediated ethoxyresorufin-O-deethylase (EROD) activities at 24 and 48h, whereas exposure to the omeprazole for 48h had no effect on the EROD activity. Clotrimazole, nocodazole and prochloraz also acted as inhibitors of EROD activities in situ in PLHC-1 cells (IC50=1.3-7.7μM), whereas omeprazole had no effect on this activity (IC50=72μM). Exposure to 10μM prochloraz resulted in 3-fold induction of CYP1A mRNA and exposure to 10μM nocodazole resulted in 16-fold induction of CYP1A mRNA levels at 24h compared to controls. In the mixture experiments, more-than-additive mixture effects between BNF and the azoles clotrimazole, prochloraz and nocodazole on EROD activities were evident, with nocodazole showing the strongest mixture effect. The presence of nocodazole increased the response to BNF up to 200-fold on CYP1A mRNA and up to 16-fold on EROD activities and prolonged the effect of BNF exposure on EROD activities by 24h or longer. This suggests that azoles that are inhibitors and/or competing substrates for the CYP1A enzymes can cause increased sensitivity to exposures to chemicals that depend on CYP1A metabolism for their elimination in situations of mixed chemical exposures. The results also suggest that the EROD biomarker response can be significantly affected in azole-contaminated areas. The responsiveness of the EROD biomarker to BNF exposure was studied in PLHC-1 that had been pre-treated with nocodazole for 5 or 24h at concentrations that are known to disassemble microtubules at 24h in these cells. Pre-treatment of PLHC-1 cells with nocodazole for either 5 or 24h had no effect on the responsiveness to BNF exposure, which implies that the EROD activity can be induced in cells with disassembled microtubules.

Beta-naphthoflavone-inducedCYP1A expression in the guppy Jenynsia multidentata: Time-dependent response, anesthetic MS-222 effect and fin analysis

Cytochrome P450 1A (CYP1A) expression in fish is used as a biomarker of exposure to organic contaminants, such PAHs, PCBs and dioxins, in the aquatic environment. South American guppy fish Jenynsia multidentata were exposed to the prototypical aryl hydrocarbon receptor (AHR) agonist beta-naphthoflavone (BNF; 1μM) and the fins were biopsied to characterize different aspects of CYP1A induction. RTq-PCR was used to quantify CYP1A mRNA levels in fish tissues. CYP1A induction in the gill, liver and anal fin (gonopodium) occurred within the first hour of waterborne exposure to BNF and persisted throughout 2, 4, 8, 24, 48 and 96h compared to controls (DMSO vehicle; p<0.05). The organ-specific temporal pattern of induction was marked by mRNA levels consistently augment as duration of exposure increases and tend to a sustained induction from 24h to 96h for gill and liver (∼15-fold and ∼50-fold over control, respectively). In gonopodium, there was a maximum CYP1A mRNA level at 4h (∼34-fold over control). Basal CYP1A mRNA levels and its induction following BNF exposure were not affected by administration of a chemical anesthetic (fish immersion in 100mgl(-1) MS-222 for 2-5min) in the gill, liver, gonopodium, dorsal or tail fin (p<0.05). In an ex vivo assay, in which small pieces of biopsied fins were exposed to BNF for 4h, high CYP1A induction was observed in the tail and gonopodium (∼49-fold and ∼69-fold, respectively) but not in the dorsal fin compared to controls. To our knowledge, this is the first study to show that a 1h waterborne exposure to an AHR agonist is sufficient to cause CYP1A induction in fish organs and fins. The present study added new information to the field regarding the use of MS-222 as an anesthetic on fish and the analysis of biopsied fins as an alternative non-lethalex vivo assay for evaluating the CYP1A biomarker in fish. This observation could be useful for planning fish toxicological bioassays and biomonitoring studies on the aquatic environments in South America.

Coevolution between Contracaecum (Nematoda, Anisakidae) and Austrolebias (Cyprinodontiformes, Rivulidae) host-parasite complex from SW Atlantic coastal basins

In recent years, molecular studies in host-parasite interactions in terms of coevolution have become important. Larvae (L3) of two species of Contracaecum were found parasitizing species of Rivulidae in the Atlantic coastal basins from Uruguay. The aim of this study is to determine the patterns of differentiation of this host-parasite complex in order to clarify possible coevolutionary events in such interaction throughout phylogeographic approach using both nuclear and mitochondrial molecular markers (internal transcribed spacers (ITS) and cytochrome oxidase subunit 1 (cox-1)). Based on both markers, intraspecific variation in Contracaecum species was lower than 2 %, while interspecific variation was greater than 10 %. Both species of Contracaecum constitute monophyletic groups. Contracaecum resulted in a paraphyletic genus when incorporating other Contracaecum species and closely related nematode sequences from GenBank. ITS regions showed that Contracaecum sp. 1 is more closely related to other species of the same genus than with their counterparts from Atlantic coastal basins in Uruguay. Haplotype network for both markers corroborate the existence of two distinct taxa. While ITS pairwise FST comparisons and the indirect estimate of gene flow confirm the existence of two distinct Contracaecum species, mitochondrial gene detected low levels of migrants between some of the populations from both species. Our results suggest that coevolution in this host-parasite complex species is plausible. Parasite cladogenetic events occur almost simultaneously with the separation of the hypothetical ancestors of each species complex of Austrolebias during Pliocene. Additionally, the two lineages of Contracaecum colonize differently the species within each of the Austrolebias complexes.

Phylogeography of Cuban Rivulus: evidence for allopatric speciation and secondary dispersal across a marine barrier

The genus Rivulus is currently comprised of two species, R. cylindraceus and R. insulaepinorum, which are endemic to Cuba. However, the taxonomic status of the latter species remains dubious because of the poor quality of the original description. In addition, a recent barcoding survey suggests that the two species may be conspecific. The aim of this study was to test the hypothesis that the two species represent a single evolutionary clade. To delimit the species and their evolutionary history, we used a combination of molecular phylogenetic analyses, with both mitochondrial and nuclear sequences, tests of phylogeographic hypotheses, combined with morphological measurements and information on known dispersal barriers and species distribution. None of the data sets support R. insulaepinorum and R. cylindraceus as separate taxa. However, a new species, restricted to the northwestern part of the main island, was identified by phylogenetic analyses, body colour pattern and geographical distribution. The evolutionary distance between the two lineages (cytb, d=15%; CAM-4, d=2.5%) indicates a long period of divergence. Phylogeographic analyses shed light on the dispersal history of R. cylindraceus, which probably originated on the Isla de la Juventud. They also suggest that each lineage had contrasting histories; Rivulus sp. is restricted to a relatively small geographic area whereas R. cylindraceus has dispersed considerably and more than once from its centre of origin, probably facilitated by sea level fluctuations. These results strengthen previous findings, i.e. that the diversity of Cuban freshwater fishes is far from well-known and deserves more in-depth studies, and that vicariance and dispersal events have resulted in a complex biogeographical landscape which has had a significant impact on the freshwater fishes of the Caribbean islands.

The past and present of an estuarine-resident fish, the "four-eyed fish" Anableps anableps (Cyprinodontiformes, Anablepidae), revealed by mtDNA sequences

Historical events, such as changes in sea level during the Pleistocene glacial cycles, had a strong impact on coastal habitats, limiting connectivity and promoting the genetic divergence of various species. In this study, we evaluated the influence of climate oscillations and the possibility of estuary function as a barrier to gene flow among populations of the four-eyed fish, Anableps anableps. This species is fully estuarine-resident, has internal fertilization, is viviparous and does not migrate across long distances. These features make the four-eyed fish an excellent model for the study of evolutionary processes related to genetic differentiation of species and populations in estuaries. The evolutionary history of A. anableps was inferred from phylogeographic and population analyses using sequences of the mitochondrial DNA Control Region of 13 populations distributed in the Amazon and Northeast Coast of Brazil from Calcoene (Amapa) to Parnaiba (Piaui). The 83 retrieved haplotypes show a pattern of four distinct mitochondrial lineages, with up to 3.4% nucleotide divergence among them. The evolutionary reconstruction suggests that these lineages diverged recently in the late Pleistocene/early Holocene after the Atlantic Ocean reaching current levels. Analysis of variability, neutrality and the genetic expansion pattern revealed that the lineages have distinct characteristics, which were shaped by the different geomorphological features of coastal regions combined with sea level oscillations over a very long period of time. Only few neighboring populations show a discreet gene flow. This study may also be helpful for designing new experiments to better understand the geomorphological evolutionary history of the estuaries of the Amazon and the Northeast Coast of Brazil using estuarine-resident species as a model.

Expression patterns of metallothionein, cytochrome P450 1A and vitellogenin genes in western mosquitofish (Gambusia affinis) in response to heavy metals

The aim of this study was to evaluate the effects of three metals (Zn, Cd and Pb) on hepatic metallothionein (MT), cytochrome P450 1A (CYP1A) and vitellogenin (Vtg) mRNA expression in the liver of adult female mosquitofish (Gambusia affinis) after 1, 3 or 8d. Both concentration-response and time-course effects of hepatic MT, CYP1A and Vtg at the transcription level were determined by quantitative real-time PCR. The results from this study showed that Zn, Cd and Pb could significantly induced MT, CYP1A and Vtg mRNA expression levels in mosquitofish. In general, this study demonstrated that heavy metals modulate MT, CYP1A and Vtg mRNA expression levels in a metal-, concentration- or time-dependent manner.

Characterization and differential expression of CPD and 6-4 DNA photolyases in Xiphophorus species and interspecies hybrids

Among the many Xiphophorus interspecies hybrid tumor models are those that exhibit ultraviolet light (UVB) induced melanoma. In previous studies, assessment of UVB induced DNA damage and nucleotide excision DNA repair has been performed in parental lines and interspecies hybrids. Species and hybrid specific differences in the levels of DNA damage induced and the dark repair rates for cyclobutane pyrimidine dimers (CPDs) and 6-4 pyrimidine pyrimidine photoproducts (6-4PPs) have been reported. However, UVB induced DNA lesions in Xiphophorus fishes are thought to primarily be repaired via light dependent CPD and 6-4PP specific photolyases. Photolyases are of evolutionary interest since they are ancient and presumably function solely to ameliorate the deleterious effects of UVB exposure. Herein, we report results from detailed studies of CPD and 6-4PP photolyase gene expression within several Xiphophorus tissues. We determined photolyase gene expression patterns before and after exposure to fluorescent light in X. maculatus, X. couchianus, and for F1 interspecies hybrids produced from crossing these two parental lines (X. maculatus Jp 163 B×X. couchianus). We present novel results showing these two photolyase genes exhibit species, tissue, and hybrid-specific differences in basal and light induced gene expression.

Cortisol release in response to UVB exposure in Xiphophorus fish

Xiphophorus fishes are comprised of 26 known species. Interspecies hybridization between select species has been utilized to produce experimental models to study melanoma development. Xiphophorus melanoma induction protocols utilize ultraviolet light (UVB) to induce DNA damage and associated downstream tumorigenesis. However, the impact of induced stress caused by the UVB treatment of the experimental animals undergoing tumor induction protocols has not been assessed. Stress is an adaptive physiological response to excessive or unpredictable environmental stimuli. The stress response in fishes may be measured by an assay of cortisol released into the water. Here, we present results from investigations of stress response during an experimental treatment and UVB exposure in Xiphophorus maculatus Jp 163 B, Xiphophorus couchianus, and F1 interspecies hybrids produced from the mating X. maculatus Jp 163 B×X. couchianus. Overall, cortisol release rates for males and females after UVB exposure showed no statistical differences. At lower UVB doses (8 and 16kJ/m(2)), X. couchianus exhibited 2 fold higher levels of DNA damage then either X. maculatus or the F1 hybrid. However, based on the cortisol release rates, none of the fish types tested induced a primary stress response at the UVB lower doses (8 and 16kJ/m(2)). In contrast, at a very high UVB dose (32kJ/m(2)) both X. maculatus and the F1 hybrid showed a 5 fold increase in the cortisol release rate. To determine the effect of pigmentation on UVB induced stress, wild type and albino Xiphophorus hellerii were exposed to UVB (32kJ/m(2)). Albino X. hellerii exhibited 3.7 fold increase in the cortisol release while wild type X. hellerii did not exhibit a significant cortisol response to UVB. Overall, the data suggest the rather low UVB doses often employed in tumor induction protocols do not induce a primary stress response in Xiphophorus fishes.

Comparative analysis of melanoma deregulated miRNAs in the medaka and Xiphophorus pigment cell cancer models

Malignant melanoma is the most aggressive and deadly form of skin cancer, with an almost 100% development of resistance to current therapeutic approaches at progression stages. The incidence of melanoma is steadily increasing worldwide. Although many details leading to the development of malignant melanoma are known, the complex process of melanomagenesis is poorly understood. MicroRNAs (miRNAs) are a class of small noncoding-RNAs of ~22nt length that regulate gene expression at the post-transcriptional level. It is now well established that deregulated miRNA expression is seen in many cancers including melanoma. To further study the miRNA functions in melanoma formation and progression we use a transgenic melanoma model in Japanese ricefish (medaka; Oryzias latipes) and the natural Xiphophorus melanoma model. In these fishes, dependent on the genetic background various histo- and patho-types of tumors appear, comparable to human melanoma types. We have studied expression profiles of ten known human melanoma-associated miRNAs and their respective target gene expression in the fish melanoma models. We show that miRNAs of the miR-17-92 cluster (miR-20a2, miR-92a1, miR-17 and miR-18a), miR-126, miR-182, miR-210 and miR-214 are upregulated and their respective target genes (RUNX1, HIF1A, TGFBR2, THBS1 and JAK2) are down-regulated in melanoma. MicroRNA-125b is down-regulated and the target genes (ERBB3a and ERBB3b) are upregulated in fish melanomas. Results provide clear evidence that the fish melanoma-associated miRNAs and respective target genes are deregulated generally like in human melanoma. Our results confirm the value of fish; such as medaka and Xiphophorus as good model systems to identify and decipher molecular mechanisms associated with malignant melanoma.

A RAD-tag genetic map for the platyfish (Xiphophorus maculatus) reveals mechanisms of karyotype evolution among teleost fish

Mammalian genomes can vary substantially in haploid chromosome number even within a small taxon (e.g., 3-40 among deer alone); in contrast, teleost fish genomes are stable (24-25 in 58% of teleosts), but we do not yet understand the mechanisms that account for differences in karyotype stability. Among perciform teleosts, platyfish (Xiphophorus maculatus) and medaka (Oryzias latipes) both have 24 chromosome pairs, but threespine stickleback (Gasterosteus aculeatus) and green pufferfish (Tetraodon nigroviridis) have just 21 pairs. To understand the evolution of teleost genomes, we made a platyfish meiotic map containing 16,114 mapped markers scored on 267 backcross fish. We tiled genomic contigs along the map to create chromosome-length genome assemblies. Genome-wide comparisons of conserved synteny showed that platyfish and medaka karyotypes remained remarkably similar with few interchromosomal translocations but with numerous intrachromosomal rearrangements (transpositions and inversions) since their lineages diverged ∼120 million years ago. Comparative genomics with platyfish shows how reduced chromosome numbers in stickleback and green pufferfish arose by fusion of pairs of ancestral chromosomes after their lineages diverged from platyfish ∼195 million years ago. Zebrafish and human genomes provide outgroups to root observed changes. These studies identify likely genome assembly errors, characterize chromosome fusion events, distinguish lineage-independent chromosome fusions, show that the teleost genome duplication does not appear to have accelerated the rate of translocations, and reveal the stability of syntenies and gene orders in teleost chromosomes over hundreds of millions of years.

UVB-induced gene expression in the skin of Xiphophorus maculatus Jp 163 B

Xiphophorus fish and interspecies hybrids represent long-standing models to study the genetics underlying spontaneous and induced tumorigenesis. The recent release of the Xiphophorus maculatus genome sequence will allow global genetic regulation studies of genes involved in the inherited susceptibility to UVB-induced melanoma within select backcross hybrids. As a first step toward this goal, we report results of an RNA-Seq approach to identify genes and pathways showing modulated transcription within the skin of X. maculatus Jp 163 B upon UVB exposure. X. maculatus Jp 163 B were exposed to various doses of UVB followed by RNA-Seq analysis at each dose to investigate overall gene expression in each sample. A total of 357 genes with a minimum expression change of 4-fold (p-adj<0.05) were identified as responsive to UVB. The molecular genetic response of Xiphophorus skin to UVB exposure permitted assessment of; (1) the basal expression level of each transcript for each skin sample, (2) the changes in expression levels for each gene in the transcriptome upon exposure to increasing doses of UVB, and (3) clusters of genes that exhibit similar patterns of change in expression upon UVB exposure. These data provide a foundation for understanding the molecular genetic response of fish skin to UVB exposure.

Antioxidant, phase II and III responses induced by lipoic acid in the fish Jenynsia multidentata (Anablapidae) and its influence on endolsulfan accumulation and toxicity

Antioxidants like lipoic acid (LA) are known to trigger augmented antioxidant and phase II and III responses. This study aimed to evaluate the effect of LA in P-glycoprotein (Pgp) expression, glutathione-S-transferase (GST) activity, total antioxidant competence, levels of lipid peroxides (TBARS) and accumulation of the organochlorine insecticide endosulfan (Endo: α-, β-isomers and sulfate metabolite) in different organs of the fish Jenynsia multidentata. One hundred and twenty females (1.55±0.07 g) were fed during 8 days with (n=60) or without (n=60) a LA enriched ration (6000 mg/kg). Four experimental groups were defined: -LA/-Endo; +LA/-Endo; -LA/+Endo; and +LA/+Endo. Endo groups were exposed during 24 h to 1.4 μg of insecticide/L. Results showed that only LA induced a significant increment in liver Pgp expression. GST activity was augmented in liver after exposure to LA or Endo. TBARS levels were lowered in liver and gills after LA pre-treatment. Total antioxidant capacity was lowered in liver of Endo exposed fish, a result that was reversed by LA pre-treatment. It is concluded that LA induced the expected effects in terms of Pgp expression, GST activity and reduced TBARS levels although favored α-Endo accumulation in brain. However, the Endo metabolism to the more persistent endosulfan sulfate was not facilitated by LA pre-treatment.

Gene flow and maintenance of genetic diversity in invasive mosquitofish (Gambusia holbrooki)

Genetic analyses contribute to studies of biological invasions by mapping the origin and dispersal patterns of invasive species occupying new territories. Using microsatellite loci, we assessed the genetic diversity and spatial population structure of mosquitofish (Gambusia holbrooki) that had invaded Spanish watersheds, along with the American locations close to the suspected potential source populations. Mosquitofish populations from the Spanish streams that were studied had similar levels of genetic diversity to the American samples; therefore, these populations did not appear to have undergone substantial losses of genetic diversity during the invasion process. Population structure analyses indicated that the Spanish populations fell into four main clusters, which were primarily associated with hydrography. Dispersal patterns indicated that local populations were highly connected upstream and downstream through active dispersal, with an average of 21.5% fish from other locations in each population. After initially introducing fish to one location in a given basin, such dispersal potential might contribute to the spread and colonization of suitable habitats throughout the entire river basin. The two-dimension isolation-by-distance pattern here obtained, indicated that the human-mediated translocation of mosquitofish among the three study basins is a regular occurrence. Overall, both phenomena, high natural dispersal and human translocation, favor gene flow among river basins and the retention of high genetic diversity, which might help retain the invasive potential of mosquitofish populations.

Paleoclimatic modeling and phylogeography of least killifish, Heterandria formosa: insights into Pleistocene expansion-contraction dynamics and evolutionary history of North American Coastal Plain freshwater biota

BACKGROUND

Climatic and sea-level fluctuations throughout the last Pleistocene glacial cycle (~130-0 ka) profoundly influenced present-day distributions and genetic diversity of Northern Hemisphere biotas by forcing range contractions in many species during the glacial advance and allowing expansion following glacial retreat ('expansion-contraction' model). Evidence for such range dynamics and refugia in the unglaciated Gulf-Atlantic Coastal Plain stems largely from terrestrial species, and aquatic species Pleistocene responses remain relatively uninvestigated. Heterandria formosa, a wide-ranging regional endemic, presents an ideal system to test the expansion-contraction model within this biota. By integrating ecological niche modeling and phylogeography, we infer the Pleistocene history of this livebearing fish (Poeciliidae) and test for several predicted distributional and genetic effects of the last glaciation.

RESULTS

Paleoclimatic models predicted range contraction to a single southwest Florida peninsula refugium during the Last Glacial Maximum, followed by northward expansion. We inferred spatial-population subdivision into four groups that reflect genetic barriers outside this refuge. Several other features of the genetic data were consistent with predictions derived from an expansion-contraction model: limited intraspecific divergence (e.g. mean mtDNA p-distance = 0.66%); a pattern of mtDNA diversity (mean Hd = 0.934; mean π = 0.007) consistent with rapid, recent population expansion; a lack of mtDNA isolation-by-distance; and clinal variation in allozyme diversity with higher diversity at lower latitudes near the predicted refugium. Statistical tests of mismatch distributions and coalescent simulations of the gene tree lent greater support to a scenario of post-glacial expansion and diversification from a single refugium than to any other model examined (e.g. multiple-refugia scenarios).

CONCLUSIONS

Congruent results from diverse data indicate H. formosa fits the classic Pleistocene expansion-contraction model, even as the genetic data suggest additional ecological influences on population structure. While evidence for Plio-Pleistocene Gulf Coast vicariance is well described for many freshwater species presently codistributed with H. formosa, this species demography and diversification departs notably from this pattern. Species-specific expansion-contraction dynamics may therefore have figured more prominently in shaping Coastal Plain evolutionary history than previously thought. Our findings bolster growing appreciation for the complexity of phylogeographical structuring within North America's southern refugia, including responses of Coastal Plain freshwater biota to Pleistocene climatic fluctuations.

Effects of steroid hormones on reproduction- and detoxification-related gene expression in adult male mosquitofish, Gambusia affinis

The molecular mechanisms that mediate fish reproduction and detoxification in response to steroid hormones were studied by using adult male western mosquitofish (Gambusia affinis) as sentinel species. The expression patterns of three vitellogenins (VtgA, VtgB and VtgC), two estrogen receptors (ERα and ERβ), two androgen receptors (ARα and ARβ), metallothionein (MT) and cytochrome P450 1A (CYP1A) in the liver and testis of adult male mosquitofish were assessed through exposure treatments with progesterone (P), testosterone (T) and 17β-estradiol (E2), alone and in combination for eight days. The results showed that expression patterns of Vtg subtype, ER subtype, AR subtype, MT and CYP1A genes in male mosquitofish varied according to tissue and specific hormone stress. Vtg subtype mRNA expression was induced in the liver in E2-added treatments, and an up-regulation of ERα mRNA expression was also observed. In addition, hormone treatments increased three Vtg subtype mRNA expression levels in the testis, at least to some extent. All hormone treatments significantly inhibited ERα, ERβ and ARβ mRNA expression in the testis. Some of hormone treatments could affect MT and CYP1A gene expression in mosquitofish. In general, multiple hormone treatments showed different effects on target gene expression compared with corresponding hormone alone. The results from the present study provided valuable information on the toxicological effects of steroid hormones in mosquitofish.

Demographic effects on the use of genetic options for the control of mosquitofish, Gambusia holbrooki

This study tests the sensitivity of genetically based pest control options based on sex ratio distortion to intra- and intersexual aggressive interactions that affect male and female survival and fitness. Data on these interactions and their impacts were gathered for the mosquitofish Gambusia holbrooki (Poeciliidae), a promiscuous species with a strongly male-biased operational sex ratio and well-documented male harassment of females. The experimental design consisted of an orthogonal combination of two population densities and three sex ratios, ranging from strongly male-biased to strongly female-biased, and long-term observations of laboratory populations. Contrary to expectations, the number of males in a population had little evident effect on population demographics. Rather, the density of adult females determined population fecundity (as a result of a stock-recruitment relationship involving females, but not males), constrained male densities (apparently as a result of cannibalism or intersexual aggression), and regulated itself (most likely through effects of intrasexual aggression on female recruitment). The principal effect of males was to constrain their own densities via effects of male-male aggression on adult male mortality rates. Through use of a realistically parameterized genetic/demographic model, we show that of three different genetic options applied to control G. holbrooki, one based on recombinant sex ratio distortion (release of Female Lethal carriers) is marginally more efficient than a sterile male release program, and both outperform an option based on chromosomal sex ratio distortion (Trojan W). Nonlinear dependence of reproductive rate on female density reduces the efficacy of all three approaches. The major effect of intra- and intersexual aggression is mediated through females, whose interactions reduce female numbers and increase the efficacy of a control program based on sex ratio. Socially mediated male mortality has a small impact on control programs due to operational sex ratios that are heavily male-biased. The sensitivity of sex ratio-based control options to social factors will depend on the mating system of the targeted pest, but evidence of widespread density-dependent population regulation suggests that, for most species, the effects of elevated adult mortality (due to intra- and intersexual aggression) on control programs are likely to be slight.

Estrogenic and androgenic effects in mosquitofish (Gambusia affinis) from streams contaminated by municipal effluent in Guangzhou, China

The study reported in this paper used mosquitofish to investigate the estrogenic and androgenic effects of municipal wastewater contamination on the stream system in Guangzhou, China. Western mosquitofish collected from a reference site and five study sites in streams forming part of the Pearl River network were dissected and analyzed for their morphological characteristics (anal fin and hemal spine characteristics) and target mRNA expression of genes (VTGα and ERα mRNA expression). Increased VTGα mRNA expression in males and decreased VTGα mRNA expression in females were observed in samples taken from four of the five study sites, with no such observations being made at the reference site. Correlation analysis indicated a significant correlation between the hemal spine morphology index and the gene transcription relative to the reference site. The multiple index also indicated that both male and female mosquitofish in contaminated streams were altered by discharged wastewater, as reflected in their morphological changes and relative induction of mRNA expression of genes in comparison to fish collected from the reference site.

Hormonal effects of tetrabromobisphenol A using a combination of in vitro and in vivo assays

The aim of this study was to investigate the hormonal effects of tetrabromobisphenol A (TBBPA) in vitro on recombinant yeasts and in vivo on mosquitofish (Gambusia affinis). The in vitro bioassays for (anti-)androgenic activities showed that TBBPA had a weak androgenic activity in vitro with recombinant yeast systems carrying human androgen receptor (hAR). In the in vivo bioassays, the gene expression patterns of vitellogenin (Vtg), estrogen receptors (ERα and ERβ), and androgen receptors (ARα and ARβ) in adult males and juveniles after exposure to TBBPA for 60 days were evaluated. Significant up-regulation of Vtg, ERα, and ERβ mRNAs was observed in the liver after exposure to 500 nM of TBBPA. In the testis, the lowest concentration of TBBPA (50 nM) markedly induced Vtg, ERβ, and ARβ mRNA expression, but the same concentration significantly inhibited ARα mRNA expression. In addition, in juveniles, 100 nM of TBBPA significantly up-regulated the expression of Vtg, ERβ, and ARα mRNAs. However, TBPPA did not cause histological alterations in the liver and testis of adult male mosquitofish. The results from this present study suggest that TBBPA could display low but multiple hormonal activities despite its low toxicity to mosquitofish.

An evaluation of the hybrid speciation hypothesis for Xiphophorus clemenciae based on whole genome sequences

Once thought rare in animal taxa, hybridization has been increasingly recognized as an important and common force in animal evolution. In the past decade, a number of studies have suggested that hybridization has driven speciation in some animal groups. We investigate the signature of hybridization in the genome of a putative hybrid species, Xiphophorus clemenciae, through whole genome sequencing of this species and its hypothesized progenitors. Based on analysis of this data, we find that X. clemenciae is unlikely to have been derived from admixture between its proposed parental species. However, we find significant evidence for recent gene flow between Xiphophorus species. Although we detect genetic exchange in two pairs of species analyzed, the proportion of genomic regions that can be attributed to hybrid origin is small, suggesting that strong behavioral premating isolation prevents frequent hybridization in Xiphophorus. The direction of gene flow between species is potentially consistent with a role for sexual selection in mediating hybridization.

Effects of 17β-trenbolone on Eastern and Western mosquitofish (Gambusia holbrooki and G. affinis) anal fin growth and gene expression patterns

The Eastern and Western mosquitofish (Gambusia holbrooki and G. affinis) are potential bioindicator organisms for endocrine disruptors. Male mosquitofish have an elongated anal fin (gonopodium) used for internal fertilization whose formation is driven by androgens. Normal female mosquitofish have a normal, rounded anal fin which undergoes elongation into a gonopodium structure when female mosquitofish are exposed to androgenic chemicals. Significant issues with using mosquitofish as a bioindicator include the lack of knowledge on how anal fin growth in females corresponds to endpoints relevant to biological integrity and the lack of information on the molecular pathways that regulate anal fin growth. The objectives of this study were to understand how androgen-induced anal fin elongation relates to changes in endpoints related to the female reproductive system and to understand how anal fin elongation occurs in androgen-exposed female mosquitofish. To achieve these objectives, adult female G. holbrooki were exposed to a vehicle control or one of three doses of the androgen 17β-trenbolone (TB) at nominal concentrations of 0.1, 1 or 10 μg TB/L. Anal fin measurements were taken and livers were used for quantitative polymerase chain reaction analysis of vitellogenin (vtg) mRNA expression at multiple time points. 10 μg TB/L induced anal fin elongation after 7 days of treatment (one-way ANOVA, p<0.05) as did 0.1 and 1 μg TB/L at later time points (one-way ANOVA, p<0.05). 10 μg TB/L significantly reduced hepatic vtg gene expression at all time points assessed (one-way ANOVA, p<0.05). There was no correlation between anal fin elongation levels and vtg gene expression (Spearman's ρ, p>0.05). In a separate experiment, female G. holbrooki and G. affinis were exposed to the vehicle control or 1 μg TB/L. Anal fins were used for qualitative gene expression analysis of the genes sonic hedgehog (shh), muscle segment homeobox C (msxC), and fibroblast growth factor receptor 1 (fgfr1) by in situ hybridization. Shh was expressed in the distal tip of the gonopodium while msxC and fgfr1 were more widely expressed along the same anal fin rays during androgen exposure. These data provide insight into the molecular pathways involved in anal fin elongation and pave the way for future work toward developing the mosquitofish into a bioindicator organism for endocrine disruptors.