Extensive hybridization and past introgression between divergent lineages in a quasi-clonal hermaphroditic fish: Ramifications for species concepts and taxonomy.

Extreme inbreeding is expected to reduce the incidence of hybridization, serving as a prezygotic barrier. Mangrove rivulus is a small killifish that reproduces predominantly by self-fertilization, producing highly homozygous lines throughout its geographic range. The Bahamas and Caribbean are inhabited by two highly diverged phylogeographic lineages of mangrove rivulus, Kryptolebias marmoratus and a 'Central clade' closely related to K. hermaphroditus from Brazil. The two lineages are largely allopatric, but recently were found in syntopy on San Salvador, Bahamas, where a single hybrid was reported. To better characterize the degree of hybridization and the possibility of secondary introgression, here we conducted a detailed genetic analysis of the contact zone on San Salvador. Two mixed populations were identified, one of which contained sexually mature hybrids. The distribution of heterozygosity at diagnostic microsatellite loci in hybrids showed that one of these hybrids was an immediate offspring from the K. marmoratus x Central clade cross, whereas the remaining five hybrids were products of reproduction by self-fertilization for 1-3 generations following the initial cross. Two hybrids had mitochondrial haplotypes of K. marmoratus and the remaining four hybrids had a haplotype of the Central clade, indicating that crosses go in both directions. In hybrids, alleles of parental lineages were represented in equal proportions suggesting lack of recent backcrossing to either of the parental lineages. However, sympatric populations of two lineages were less diverged than allopatric populations, consistent with introgression. Results are discussed in terms of applicability of the biological species concept for isogenic, effectively clonal, organisms.

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

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

Mitochondria, sex and variation in routine metabolic rate.

Variation in the metabolic costs associated with organismal maintenance may play a key role in determining fitness, and thus these differences among individuals are likely to be subject to natural selection. Although the evolvability of maintenance metabolism depends on its underlying genetic architecture, relatively little is known about the nature of genetic variation that underlies this trait. To address this, we measured variation in routine metabolic rate (MO2routine ), an index of maintenance metabolism, within and among three populations of Atlantic killifish, Fundulus heteroclitus, including a population from a region of genetic admixture between two subspecies. Polygenic association tests among individuals from the admixed population identified 54 single nucleotide polymorphisms (SNPs) that were associated with MO2routine , and these SNPs accounted for 43% of interindividual variation in this trait. However, genetic associations with MO2routine involved different SNPs if females and males were analysed separately, and there was a sex-dependent effect of mitochondrial genotype on variation in routine metabolism. These results imply that there are sex-specific genetic mechanisms, and potential mitonuclear interactions, that underlie variation in MO2routine . Additionally, there was evidence for epistatic interactions between 17% of the possible pairs of trait-associated SNPs, suggesting that epistatic effects on MO2routine are common. These data demonstrate not only that phenotypic variation in this ecologically important trait has a polygenic basis with considerable epistasis among loci, but also that these underlying genetic mechanisms, and particularly the role of mitochondrial genotype, may be sex-specific.

Evolutionary Toxicogenomics of the Striped Killifish (Fundulus majalis) in the New Bedford Harbor (Massachusetts, USA).

In this paper, we used a Genotyping-by-Sequencing (GBS) approach to find and genotype more than 4000 genome-wide SNPs (Single Nucleotide Polymorphisms) from striped killifish exposed to a variety of polychlorinated biphenyls (PCBs) and other aromatic pollutants in New Bedford Harbor (NBH, Massachusetts, USA). The aims of this study were to identify the genetic consequences of exposure to aquatic pollutants and detect genes that may be under selection. Low genetic diversity (HE and π) was found in the site exposed to the highest pollution level, but the pattern of genetic diversity did not match the pollution levels. Extensive connectivity was detected among sampling sites, which suggests that balanced gene flow may explain the lack of genetic variation in response to pollution levels. Tests for selection identified 539 candidate outliers, but many of the candidate outliers were not shared among tests. Differences among test results likely reflect different test assumptions and the complex pollutant mixture. Potentially, selectively important loci are associated with 151 SNPs, and enrichment analysis suggests a likely involvement of these genes with pollutants that occur in NBH. This result suggests that selective processes at genes targeted by pollutants may be occurring, even at a small geographical scale, and may allow the local striped killifish to resist the high pollution levels.

Sequence and functional characterization of hypoxia-inducible factors, HIF1α, HIF2αa, and HIF3α, from the estuarine fish, Fundulus heteroclitus.

The hypoxia-inducible factor (HIF) family of transcription factors plays central roles in the development, physiology, pathology, and environmental adaptation of animals. Because many aquatic habitats are characterized by episodes of low dissolved oxygen, fish represent ideal models to study the roles of HIF in the response to aquatic hypoxia. The estuarine fish Fundulus heteroclitus is found in habitats prone to hypoxia. It responds to low oxygen via behavioral, physiological, and molecular changes, and one member of the HIF family, HIF2α, has been previously described. Herein, cDNA sequencing, phylogenetic analyses, and genomic approaches were used to determine other members of the HIFα family from F. heteroclitus and their relationships to HIFα subunits from other vertebrates. In vitro and cellular approaches demonstrated that full-length forms of HIF1α, HIF2α, and HIF3α independently formed complexes with the ß-subunit, aryl hydrocarbon receptor nuclear translocator, to bind to hypoxia response elements and activate reporter gene expression. Quantitative PCR showed that HIFα mRNA abundance varied among organs of normoxic fish in an isoform-specific fashion. Analysis of the F. heteroclitus genome revealed a locus encoding a second HIF2α-HIF2αb-a predicted protein lacking oxygen sensing and transactivation domains. Finally, sequence analyses demonstrated polymorphism in the coding sequence of each F. heteroclitus HIFα subunit, suggesting that genetic variation in these transcription factors may play a role in the variation in hypoxia responses among individuals or populations.

Acclimation and acute temperature effects on population differences in oxidative phosphorylation.

Temperature changes affect metabolism on acute, acclamatory, and evolutionary time scales. To better understand temperature's affect on metabolism at these different time scales, we quantified cardiac oxidative phosphorylation (OxPhos) in three Fundulus taxa acclimated to 12 and 28°C and measured at three acute temperatures (12, 20, and 28°C). The Fundulus taxa (northern Maine and southern Georgia F. heteroclitus, and a sister taxa, F. grandis) were used to identify evolved changes in OxPhos. Cardiac OxPhos metabolism was quantified by measuring six traits: state 3 (ADP and substrate-dependent mitochondrial respiration); E state (uncoupled mitochondrial activity); complex I, II, and IV activities; and LEAK ratio. Acute temperature affected all OxPhos traits. Acclimation only significantly affected state 3 and LEAK ratio. Populations were significantly different for state 3. In addition to direct effects, there were significant interactions between acclimation and population for complex I and between population and acute temperature for state 3. Further analyses suggest that acclimation alters the acute temperature response for state 3, E state, and complexes I and II: at the low acclimation temperature, the acute response was dampened at low assay temperatures, and at the high acclimation temperature, the acute response was dampened at high assay temperatures. Closer examination of the data also suggests that differences in state 3 respiration and complex I activity between populations were greatest between fish acclimated to low temperatures when assayed at high temperatures, suggesting that differences between the populations become more apparent at the edges of their thermal range.

Rivulus berovidesi, a new killifish species (Teleostei: Rivulidae) from western Cuba.

Rivulus berovidesi, a new killifish species, is described from a small stream in Sierra de Cajalbana, northwestern Cuba. It is readily distinguished from Rivulus cylindraceus Poey by the combination of an exclusive color pattern and meristic characters such as a d-type frontal scalation pattern (versus e-type pattern in Rivulus cylindraceus). The current diagnosis of Rivulus berovidesi based on chromatic, morphological and meristic characters is consistent with a recent molecular analysis of this genus in Cuba.

Regulation of pregnane-X-receptor, CYP3A and P-glycoprotein genes in the PCB-resistant killifish (Fundulus heteroclitus) population from New Bedford Harbor.

Killifish survive and reproduce in the New Bedford Harbor (NBH) in Massachusetts (MA), USA, a site severely contaminated with polychlorinated biphenyls (PCBs) for decades. Levels of 22 different PCB congeners were analyzed in liver from killifish collected in 2008. Concentrations of dioxin-like PCBs in liver of NBH killifish were ∼400 times higher, and the levels of non-dioxin-like PCBs ∼3000 times higher than in killifish from a reference site, Scorton Creek (SC), MA. The NBH killifish are known to be resistant to the toxicity of dioxin-like compounds and to have a reduced aryl hydrocarbon receptor (AhR) signaling response. Little is known about the responses of these fish to non-dioxin-like PCBs, which are at extraordinarily high levels in NBH fish. In mammals, some non-dioxin-like PCB congeners act through nuclear receptor 1I2, the pregnane-X-receptor (PXR). To explore this pathway in killifish, a PXR cDNA was sequenced and its molecular phylogenetic relationship to other vertebrate PXRs was determined. Killifish were also collected in 2009 from NBH and SC, and after four months in the laboratory they were injected with a single dose of either the dioxin-like PCB 126 (an AhR agonist) or the non-dioxin-like PCB 153 (a mammalian PXR agonist). Gills and liver were sampled three days after injection and transcript levels of genes encoding PXR, cytochrome P450 3A (CYP3A), P-glycoprotein (Pgp), AhR2 and cytochrome P450 1A (CYP1A) were measured by quantitative PCR. As expected, there was little effect of PCB exposure on mRNA expression of AhR2 or CYP1A in liver and gills of NBH fish. In NBH fish, but not in SC fish, there was increased mRNA expression of hepatic PXR, CYP3A and Pgp upon exposure to either of the two PCB congeners. However, basal PXR and Pgp mRNA levels in liver of NBH fish were significantly lower than in SC fish. A different pattern was seen in gills, where there were no differences in basal mRNA expression of these genes between the two populations. In SC fish, but not in NBH fish, there was increased mRNA expression of branchial PXR and CYP3A upon exposure to PCB126 and of CYP3A upon exposure to PCB153. The results suggest a difference between the two populations in non-AhR transcription factor signaling in liver and gills, and that this could involve killifish PXR. It also implies possible cross-regulatory interactions between that factor (presumably PXR) and AhR2 in liver of these fish.

Targeted approach to identify genetic loci associated with evolved dioxin tolerance in Atlantic killifish (Fundulus heteroclitus).

BACKGROUND: The most toxic aromatic hydrocarbon pollutants are categorized as dioxin-like compounds (DLCs) to which extreme tolerance has evolved independently and contemporaneously in (at least) four populations of Atlantic killifish (Fundulus heteroclitus). Surprisingly, the magnitude and phenotype of DLC tolerance is similar among these killifish populations that have adapted to varied, but highly aromatic hydrocarbon-contaminated urban/industrialized estuaries of the US Atlantic coast. Multiple tolerant and neighboring sensitive killifish populations were compared with the expectation that genetic loci associated with DLC tolerance would be revealed. RESULTS: Since the aryl hydrocarbon receptor (AHR) pathway partly or fully mediates DLC toxicity in vertebrates, single nucleotide polymorphisms (SNPs) from 42 genes associated with the AHR pathway were identified to serve as targeted markers. Wild fish (N = 36/37) from four highly tolerant killifish populations and four nearby sensitive populations were genotyped using 59 SNP markers. Similar to other killifish population genetic analyses, strong genetic differentiation among populations was detected, consistent with isolation by distance models. When DLC-sensitive populations were pooled and compared to pooled DLC-tolerant populations, multi-locus analyses did not distinguish the two groups. However, pairwise comparisons of nearby tolerant and sensitive populations revealed high differentiation among sensitive and tolerant populations at these specific loci: AHR 1 and 2, cathepsin Z, the cytochrome P450s (CYP1A and 3A30), and the NADH dehydrogenase subunits. In addition, significant shifts in minor allele frequency were observed at AHR2 and CYP1A loci across most sensitive/tolerant pairs, but only AHR2 exhibited shifts in the same direction across all pairs. CONCLUSIONS: The observed differences in allelic composition at the AHR2 and CYP1A SNP loci were identified as significant among paired sensitive/tolerant populations of Atlantic killifish with multiple statistical tests. The genetic patterns reported here lend support to the argument that AHR2 and CYP1A play a role in the adaptive response to extreme DLC contamination. Additional functional assays are required to isolate the exact mechanism of DLC tolerance.

Contrasting phylogeographic histories between broadly sympatric topminnows in the Fundulus notatus species complex.

Sympatrically distributed closely related species provide opportunities for studying evolutionary patterns of diversification. Such studies must account for historical contingencies in interpreting contemporary patterns of variation. Topminnows in the Fundulus notatus species complex are distributed sympatrically across much of the southern and Midwestern United States. Throughout most of their ranges F. olivaceus is often found in headwater stream habitats, and F. notatus is more typically distributed along the margins of larger river habitats. However, in some drainages, ecological associations of the respective species are reversed, with F. notatus populations isolated in headwater streams and F. olivaceus in downstream river habitats. Phylogeographic analyses of AFLP marker and multi-locus sequence data detected historical isolation in F. notatus consistent with pre-Pleistocene drainage patterns. Four F. notatus clades corresponded to (i) the Western Gulf Slope, (ii) the southwestern Ouachita Highlands, (iii) the Mobile Basin, and (iv) central Coastal Plain and Mississippi River Basin. In contrast, a relative lack of range-wide geographic structure in F. olivaceus is consistent with recent range expansion over much of the same geographic area. The southwestern Ouachita Highlands and Mobile Basin F. notatus clades corresponded to regions where ecological associations between the two species are reversed, providing evidence of the independent evolution of variation in contemporary habitat associations. Fundulus olivaceus from several drainages demonstrated introgression of mitochondrial DNA from F. notatus, but none of the sites in this study included individuals with hybrid ancestry in their nuclear genome. Phylogenetic analyses that included only nuclear loci supported the reciprocal monophyly of F. notatus, F. olivaceus and a third narrowly endemic species, Fundulus euryzonus, and supported a sister relationship between F. olivaceus and F. euryzonus.

Genetic evaluation of captive populations of endangered species and merging of populations: Gila topminnows as an example.

To avoid extinction, captive populations of a number of endangered species have been established. While in captivity, these populations have been managed to retain genetic variation although direct evaluation of this strategy using molecular markers is not common. In addition, when the number of founders for a captive population is small, other founders or populations may be added to increase genetic variation. Here we examined refugial populations of the endangered Gila topminnow (Poeciliopsis occidentalis) from 4 locations in the southwestern United States. We found that over 5 years (about 10 generations), genetic variation as measured by 5 microsatellite loci was not lost, presumably because the adult census population size was 500 or greater. In addition, some variation not initially observed was observed later. Some of these variants may have been missed because of sampling but it appears that some may have been contributed by new mutations. In addition, 2 populations of successfully merged ancestry from the 4 source populations were examined. Based on population-specific markers and a quantitative evaluation of ancestry using a likelihood approach, it appears that ancestry from each of the source populations was retained in both populations.

Seasonal diet composition of Fundulus lima (Cyprinodontiformes: Fundulidae) in two oasis systems of Baja California Sur, Mexico

Fundulus lima inhabits river drainage systems and is threatened after the introduction of cichlids in the area. To support conservation programs, the spatial and temporal variation of the diet composition of this endangered killifish, was determined in two oasis systems of Baja California Sur, Mexico (San Ignacio and La Purisima river drainages), during rainy and dry seasons. F. lima was captured by using passive and active capture techniques. A total of 192 stomach contents of F. lima was analyzed. The contribution of each prey item in the diet composition was quantified by means of the indices of occurrence frequency (% OF), numerical (%N) and volume (%V) percentages. The relative importance of each prey item was determined according to the percentage of the Relative Importance Index (%RII). The similarity of the diet was calculated between hydrological basins (populations combined by basin), seasons (rainy versus dry months), sexes and size classes, by using Schoener’s resource overlap index. We used two ecological indices to determine the type of feeding strategy exhibited by the fish: (1) niche breadth of Levins and (2) proportional similarity of Feisinger. Sand was the most abundant item in the stomach content of killifishes from both drainages (39% and 47%, respectively). Diet composition was similar for both drainages (74%) as well as among their respective size classes; however, it was different between sexes. In both drainages, F. lima predated mainly on diatom algae, dipterous and trichopteran larvae, and fish scales during the dry season; while it preferred dipterous larvae, filamentous algae and ostracods in the rainy season. A feeding strategy of opportunist type was exhibited by F. lima during the rainy season, changing to specialist type during the dry season. This information will be the basis for future investigations related to the conservation of this endangered species and its habitat. Rev. Biol. Trop. 59 (4): 1669-1678. ...

Fundulus lima habita sistemas de aguas continentales, y se encuentra amenazado tras la introducción de cíclidos en el área de estudio. Para poder llevar a cabo programas de conservación, fue determinada la dieta de la sardinilla peninsular en peligro de extinción, en dos sistemas de oasis de Baja California Sur, México (cuencas de los rios San Ignacio y La Purísima), durante las épocas secas y de lluvias. Una gran proporción de arena fue encontrada en el contenido estomacal de este pez para ambas cuencas (39% y 47%, respectivamente). La composición de la dieta fue similar entre ambas cuencas (74%), como también dentro de sus respectivas clases de tallas; sin embargo, fue diferente entre sexos. Durante la época de secas, F. lima consumió principalmente diatomeas, larvas de dípteros y tricópteros, además de escamas de peces; mientras que en la época de lluvias tuvo preferencia por las larvas de dípteros, algas filamentosas y ostrácodos. Este pez exhibe una estrategia alimentaria de tipo oportunista durante la época de lluvias, la cual cambia a tipo especialista en la época de secas. Esta información será la base para futuras investigaciones relacionadas con la conservación de esta especie, en peligro de extinción, y su hábitat.

Evolutionary and functional analyses of cytochrome P4501A promoter polymorphisms in natural populations.

The functional importance of variable, transcriptional regulatory sequences within and among natural populations is largely unexplored. We analysed the cytochrome P4501A (CYP1A) promoter in three populations of the minnow, Fundulus heteroclitus, because two SNPs in the promoter and first intron of CYP1A are under selection in populations adapted to pollutants. To define the importance of these SNPs, 1630 bp of the CYP1A promoter and first intron and exon were sequenced in eight individuals from three populations: a population from a polluted environment resistant to some aromatic pollutants and two flanking reference populations. CYP1A is induced by many aromatic pollutants, but in populations adapted to pollutants, CYP1A has been shown to be refractory to induction. We were interested in understanding whether variation in the CYP1A promoter explains mechanism(s) of adaptation to these aromatic pollutants. The CYP1A promoter was extremely variable (an average of 9.3% of the promoter nucleotides varied among all populations) and exhibited no fixed differences between populations. As CYP1A is poorly inducible in adapted fish, we hypothesized that CYP1A promoter regions might vary functionally between populations. Unexpectedly, in vitro analysis showed significantly greater transcription from CYP1A promoters found in the population from the polluted environment relative to promoters found in both reference populations. Thus, despite extensive variation among populations and lack of fixed differences between populations, individuals from a polluted environment have significantly enhanced promoter activity. These data demonstrate that intraspecific variation, which provides the raw material for natural selection to act on, can occur while maintaining promoter function.

The effects of mitochondrial genotype on hypoxic survival and gene expression in a hybrid population of the killifish, Fundulus heteroclitus.

The physiological link between oxygen availability and mitochondrial function is well established. However, whether or not fitness variation is associated with mitochondrial genotypes in the field remains a contested topic in evolutionary biology. In this study, we draw on a population of the teleost fish, Fundulus heteroclitus, where functionally distinct subspecies hybridize, likely as a result of past glacial events. We had two specific aims: (i) to determine the effect of mtDNA genotype on survivorship of male and female fish under hypoxic stress and (ii) to determine the effect of hypoxic stress, sex and mtDNA genotype on gene expression. We found an unexpected and highly significant effect of sex on survivorship under hypoxic conditions, but no significant effect of mtDNA genotype. Gene expression analyses revealed hundreds of transcripts differentially regulated by sex and hypoxia. Mitochondrial transcripts and other predicted pathways were among those influenced by hypoxic stress, and a transcript corresponding to the mtDNA control region was the most highly suppressed transcript under the conditions of hypoxia. An RT-PCR experiment on the control region was consistent with microarray results. Effects of mtDNA sequence variation on genome expression were limited; however, a potentially important epistasis between mtDNA sequence and expression of a nuclear-encoded mitochondrial translation protein was discovered. Overall, these results confirm that mitochondrial regulation is a major component of hypoxia tolerance and further suggest that purifying selection has been the predominant selective force on mitochondrial genomes in these two subspecies.

Ecologically and evolutionarily important SNPs identified in natural populations.

Evolution by natural selection acts on natural populations amidst migration, gene-by-environmental interactions, constraints, and tradeoffs, which affect the rate and frequency of adaptive change. We asked how many and how rapidly loci change in populations subject to severe, recent environmental changes. To address these questions, we used genomic approaches to identify randomly selected single nucleotide polymorphisms (SNPs) with evolutionarily significant patterns in three natural populations of Fundulus heteroclitus that inhabit and have adapted to highly polluted Superfund sites. Three statistical tests identified 1.4-2.5% of SNPs that were significantly different from the neutral model in each polluted population. These nonneutral patterns in populations adapted to highly polluted environments suggest that these loci or closely linked loci are evolving by natural selection. One SNP identified in all polluted populations using all tests is in the gene for the xenobiotic metabolizing enzyme, cytochrome P4501A (CYP1A), which has been identified previously as being refractory to induction in the three highly polluted populations. Extrapolating across the genome, these data suggest that rapid evolutionary change in natural populations can involve hundreds of loci, a few of which will be shared in independent events.

Seasonal diet composition of Fundulus lima (Cyprinodontiformes: Fundulidae) in two oasis systems of Baja California Sur, Mexico.

Fundulus lima inhabits river drainage systems and is threatened after the introduction of cichlids in the area. To support conservation programs, the spatial and temporal variation of the diet composition of this endangered killifish, was determined in two oasis systems of Baja California Sur, Mexico (San Ignacio and La Purisima river drainages), during rainy and dry seasons. F. lima was captured by using passive and active capture techniques. A total of 192 stomach contents of F. lima was analyzed. The contribution of each prey item in the diet composition was quantified by means of the indices of occurrence frequency (% OF), numerical (% N) and volume (% V) percentages. The relative importance of each prey item was determined according to the percentage of the Relative Importance Index (% RII). The similarity of the diet was calculated between hydrological basins (populations combined by basin), seasons (rainy versus dry months), sexes and size classes, by using Schoener's resource overlap index. We used two ecological indices to determine the type of feeding strategy exhibited by the fish: (1) niche breadth of Levins and (2) proportional similarity of Feisinger. Sand was the most abundant item in the stomach content of killifishes from both drainages (39% and 47%, respectively). Diet composition was similar for both drainages (74%) as well as among their respective size classes; however, it was different between sexes. In both drainages, F. lima predated mainly on diatom algae, dipterous and trichopteran larvae, and fish scales during the dry season; while it preferred dipterous larvae, filamentous algae and ostracods in the rainy season. A feeding strategy of opportunist type was exhibited by F. lima during the rainy season, changing to specialist type during the dry season. This information will be the basis for future investigations related to the conservation of this endangered species and its habitat.

Molecular cloning and brain distribution of three types of gonadotropin-releasing hormone from mummichog Fundulus heteroclitus.

Complementary DNAs encoding gonadotropin-releasing hormone (GnRH) precursors were cloned from the mummichog Fundulus heteroclitus brain, showing that this species has three GnRH forms, i.e. medaka Oryzias latipes GnRH (mdGnRH), chicken GnRH-II (cGnRH-II) and Atlantic salmon Salmo salar GnRH (sGnRH). The F. heteroclitus prepro GnRHs have common structural architectures of vertebrate GnRHs, consisting of the signal peptide, 10 amino acids of mature peptide, GKR sequence and GnRH-associated peptide (GAP). Phylogenetic analysis of fish prepro GnRHs showed that F. heteroclitus mdGnRH is a homologue of sbGnRHs and mdGnRHs of other acanthopterygian. Quantitative real-time PCR revealed that mdGnRH was abundantly expressed in the olfactory bulb and in olfactory lobe areas and is expressed in the pituitary. The cGnRH-II was mainly expressed in the midbrain and interbrain areas, and the sGnRH was expressed not only in the olfactory bulb but also in other regions of the brain. These results suggest that the mdGnRH is involved in the stimulation of gonadotrophs in the pituitary, whereas cGnRH-II and sGnRH are involved in neurotransmission and neuromodulation.

Molecular identification of Na(+)-H(+) exchanger isoforms (NHE2) in the gills of the euryhaline teleost Fundulus heteroclitus.

In the current study, reverse-transcription polymerase chain reaction (RT-PCR) and rapid amplification of cDNA ends (RACE) PCR were used to clone full-length putative Na(+)-H(+) exchanger isoforms (NHE2a) cDNA from the gills of Fundulus heteroclitus. The 2480 bp cDNA includes a coding region for a protein that shows a 57% amino acid homology to rabbit NHE2. These sequences allowed data mining of available fish genome data, which revealed at least three NHE2 subtypes in some teleost species.

Long-term retention of self-fertilization in a fish clade.

Among vertebrate animals, only the mangrove rivulus (Kryptolebias marmoratus) was known to self-fertilize. Here, we use microsatellite analyses to document a high selfing rate (97%) in a related nominal species, Kryptolebias ocellatus, which likewise is androdioecious (populations consist of males and hermaphrodites). In contrast, we find no evidence of self-fertilization in Kryptolebias caudomarginatus (an androdioecious species closely related to the marmoratus-ocellatus clade) or in Kryptolebias brasiliensis (a dioecious outgroup). These findings indicate that the initiation of self-fertilization predated the origin of the marmoratus-ocellatus clade. From mitochondrial DNA sequences and microsatellite data, we document a substantial genetic distance between Kryptolebias marmoratus and K. ocellatus, implying that the selfing capacity has persisted in these fishes for at least several hundred thousand years.