Open Access
Refine
Year of publication
Document Type
- Article (19)
- Postprint (19)
- Doctoral Thesis (1)
Language
- English (39)
Is part of the Bibliography
- yes (39)
Keywords
- conservation genetics (5)
- Microsatellites (3)
- Phylogeography (3)
- Acoustic-Signals (2)
- Allozymes (2)
- Approximate Bayesian Computation (2)
- Beneficial mutations (2)
- Bombina bombina (2)
- Convergent Evolution (2)
- Differential gene expression (2)
Institute
Genetic divergence is impacted by many factors, including phylogenetic history, gene flow, genetic drift, and divergent selection. Rotifers are an important component of aquatic ecosystems, and genetic variation is essential to their ongoing adaptive diversification and local adaptation. In addition to coding sequence divergence, variation in gene expression may relate to variable heat tolerance, and can impose ecological barriers within species. Temperature plays a significant role in aquatic ecosystems by affecting species abundance, spatio-temporal distribution, and habitat colonization. Recently described (formerly cryptic) species of the Brachionus calyciflorus complex exhibit different temperature tolerance both in natural and in laboratory studies, and show that B. calyciflorus sensu stricto (s.s.) is a thermotolerant species. Even within B. calyciflorus s.s., there is a tendency for further temperature specializations. Comparison of expressed genes allows us to assess the impact of stressors on both expression and sequence divergence among disparate populations within a single species. Here, we have used RNA-seq to explore expressed genetic diversity in B. calyciflorus s.s. in two mitochondrial DNA lineages with different phylogenetic histories and differences in thermotolerance. We identify a suite of candidate genes that may underlie local adaptation, with a particular focus on the response to sustained high or low temperatures. We do not find adaptive divergence in established candidate genes for thermal adaptation. Rather, we detect divergent selection among our two lineages in genes related to metabolism (lipid metabolism, metabolism of xenobiotics).
Genetic divergence is impacted by many factors, including phylogenetic history, gene flow, genetic drift, and divergent selection. Rotifers are an important component of aquatic ecosystems, and genetic variation is essential to their ongoing adaptive diversification and local adaptation. In addition to coding sequence divergence, variation in gene expression may relate to variable heat tolerance, and can impose ecological barriers within species. Temperature plays a significant role in aquatic ecosystems by affecting species abundance, spatio-temporal distribution, and habitat colonization. Recently described (formerly cryptic) species of the Brachionus calyciflorus complex exhibit different temperature tolerance both in natural and in laboratory studies, and show that B. calyciflorus sensu stricto (s.s.) is a thermotolerant species. Even within B. calyciflorus s.s., there is a tendency for further temperature specializations. Comparison of expressed genes allows us to assess the impact of stressors on both expression and sequence divergence among disparate populations within a single species. Here, we have used RNA-seq to explore expressed genetic diversity in B. calyciflorus s.s. in two mitochondrial DNA lineages with different phylogenetic histories and differences in thermotolerance. We identify a suite of candidate genes that may underlie local adaptation, with a particular focus on the response to sustained high or low temperatures. We do not find adaptive divergence in established candidate genes for thermal adaptation. Rather, we detect divergent selection among our two lineages in genes related to metabolism (lipid metabolism, metabolism of xenobiotics).
Background
Organisms are expected to respond to changing environmental conditions through local adaptation, range shift or local extinction. The process of local adaptation can occur by genetic changes or phenotypic plasticity, and becomes especially relevant when dispersal abilities or possibilities are somehow constrained. For genetic changes to occur, mutations are the ultimate source of variation and the mutation rate in terms of a mutator locus can be subject to evolutionary change. Recent findings suggest that the evolution of the mutation rate in a sexual species can advance invasion speed and promote adaptation to novel environmental conditions. Following this idea, this work uses an individual-based model approach to investigate if the mutation rate can also evolve in a sexual species experiencing different conditions of directional climate change, under different scenarios of colored stochastic environmental noise, probability of recombination and of beneficial mutations. The color of the noise mimicked investigating the evolutionary dynamics of the mutation rate in different habitats.
Results
The results suggest that the mutation rate in a sexual species experiencing directional climate change scenarios can evolve and reach relatively high values mainly under conditions of complete linkage of the mutator locus and the adaptation locus. In contrast, when they are unlinked, the mutation rate can slightly increase only under scenarios where at least 50% of arising mutations are beneficial and the rate of environmental change is relatively fast. This result is robust under different scenarios of stochastic environmental noise, which supports the observation of no systematic variation in the mutation rate among organisms experiencing different habitats.
Conclusions
Given that 50% beneficial mutations may be an unrealistic assumption, and that recombination is ubiquitous in sexual species, the evolution of an elevated mutation rate in a sexual species experiencing directional climate change might be rather unlikely. Furthermore, when the percentage of beneficial mutations and the population size are small, sexual species (especially multicellular ones) producing few offspring may be expected to react to changing environments not by adaptive genetic change, but mainly through plasticity. Without the ability for a plastic response, such species may become – at least locally – extinct.
Background
Organisms are expected to respond to changing environmental conditions through local adaptation, range shift or local extinction. The process of local adaptation can occur by genetic changes or phenotypic plasticity, and becomes especially relevant when dispersal abilities or possibilities are somehow constrained. For genetic changes to occur, mutations are the ultimate source of variation and the mutation rate in terms of a mutator locus can be subject to evolutionary change. Recent findings suggest that the evolution of the mutation rate in a sexual species can advance invasion speed and promote adaptation to novel environmental conditions. Following this idea, this work uses an individual-based model approach to investigate if the mutation rate can also evolve in a sexual species experiencing different conditions of directional climate change, under different scenarios of colored stochastic environmental noise, probability of recombination and of beneficial mutations. The color of the noise mimicked investigating the evolutionary dynamics of the mutation rate in different habitats.
Results
The results suggest that the mutation rate in a sexual species experiencing directional climate change scenarios can evolve and reach relatively high values mainly under conditions of complete linkage of the mutator locus and the adaptation locus. In contrast, when they are unlinked, the mutation rate can slightly increase only under scenarios where at least 50% of arising mutations are beneficial and the rate of environmental change is relatively fast. This result is robust under different scenarios of stochastic environmental noise, which supports the observation of no systematic variation in the mutation rate among organisms experiencing different habitats.
Conclusions
Given that 50% beneficial mutations may be an unrealistic assumption, and that recombination is ubiquitous in sexual species, the evolution of an elevated mutation rate in a sexual species experiencing directional climate change might be rather unlikely. Furthermore, when the percentage of beneficial mutations and the population size are small, sexual species (especially multicellular ones) producing few offspring may be expected to react to changing environments not by adaptive genetic change, but mainly through plasticity. Without the ability for a plastic response, such species may become – at least locally – extinct.
Ecological and physiological factors lead to different contamination patterns in individual marine mammals. The objective of the present study was to assess whether variations in contamination profiles are indicative of social structures of young male sperm whales as they might reflect a variation in feeding preferences and/or in utilized feeding grounds. We used a total of 61 variables associated with organic compounds and trace element concentrations measured in muscle, liver, kidney and blubber gained from 24 sperm whales that stranded in the North Sea in January and February 2016. Combining contaminant and genetic data, there is evidence for at least two cohorts with different origin among these stranded sperm whales; one from the Canary Island region and one from the northern part of the Atlantic. While genetic data unravel relatedness and kinship, contamination data integrate over areas, where animals occured during their lifetime. Especially in long-lived animals with a large migratory potential, as sperm whales, contamination data may carry highly relevant information about aggregation through time and space.
Effective communication among sympatric species is often instrumental for behavioural isolation, where the failure to successfully discriminate between potential mates could lead to less fit hybrid offspring. Discrimination between con- and heterospecifics tends to occur more often in the sex that invests more in offspring production, i.e. females, but males may also mediate reproductive isolation. In this study, we show that among two Campylomormyrus Africanweakly electric fish species, males preferentially associate with conspecific females during choice tests using live fish as stimuli, i.e. when all sensory modalities potentially used for communication were present. We then conducted playback experiments to determine whether the species-specific electric organ discharge (EOD) used for electrocommunication serves as the cue for this conspecific association preference. Interestingly, only C. compressirostris males associated significantly more with the conspecific EOD waveform when playback stimuli were provided, while no such association preference was observed in C. tamandua males. Given our results, the EOD appears to serve, in part, as a male-mediated pre-zygotic isolation mechanism among sympatric species. However, the failure of C. tamandua males to discriminate between con- and heterospecific playback discharges suggests that multiple modalities may be necessary for species recognition in some African weakly electric fish species.
The Strange-tailed Tyrant Alectrurus risora (Aves: Tyrannidae) is an endemic species of southern South American grasslands that suffered a 90% reduction of its original distribution due to habitat transformation. This has led the species to be classified as globally Vulnerable. By the beginning of the last century, populations were partially migratory and moved south during the breeding season. Currently, the main breeding population inhabits the Ibera wetlands in the province of Corrientes, north-east Argentina, where it is resident all year round. There are two remaining small populations in the province of Formosa, north-east Argentina, and in southern Paraguay, which are separated from the main population by the Parana-Paraguay River and its continuous riverine forest habitat. The populations of Corrientes and Formosa are separated by 300 km and the grasslands between populations are non-continuous due to habitat transformation. We used mtDNA sequences and eight microsatellite loci to test if there were evidences of genetic isolation between Argentinean populations. We found no evidence of genetic structure between populations (Phi(ST) = 0.004, P = 0.32; Fst = 0.01, P = 0.06), which can be explained by either retained ancestral polymorphism or by dispersal between populations. We found no evidence for a recent demographic bottleneck in nuclear loci. Our results indicate that these populations could be managed as a single conservation unit on a regional scale. Conservation actions should be focused on preserving the remaining network of areas with natural grasslands to guarantee reproduction, dispersal and prevent further decline of populations.
The complete mitochondrial genome of a European fire-bellied toad (Bombina bombina) from Germany
(2019)
The European fire-bellied toad, Bombina bombina, is a small aquatic toad belonging to the family Bombinatoridae. The species is native to the lowlands of Central and Eastern Europe, where population numbers have been in decline in recent past decades. Here, we present the first complete mitochondrial genome of the endangered European fire-bellied toad from Northern Germany recovered using iterative mapping. Phylogenetic analyses including other representatives of the Bombinatoridae placed our German specimen as sister to a Polish B. bombina sequence with high support. This finding is congruent with the postulated Pleistocene history of the species. Our complete mitochondrial genome represents an important resource for further population analysis of the European fire-bellied toad, especially those found within Germany.
The complete mitochondrial genome of a European fire-bellied toad (Bombina bombina) from Germany
(2019)
The European fire-bellied toad, Bombina bombina, is a small aquatic toad belonging to the family Bombinatoridae. The species is native to the lowlands of Central and Eastern Europe, where population numbers have been in decline in recent past decades. Here, we present the first complete mitochondrial genome of the endangered European fire-bellied toad from Northern Germany recovered using iterative mapping. Phylogenetic analyses including other representatives of the Bombinatoridae placed our German specimen as sister to a Polish B. bombina sequence with high support. This finding is congruent with the postulated Pleistocene history of the species. Our complete mitochondrial genome represents an important resource for further population analysis of the European fire-bellied toad, especially those found within Germany.
Communication is key to a wide variety of animal behaviours and multiple modalities are often involved in this exchange of information from sender to receiver. The communication of African weakly electric fish, however, is thought to be predominantly unimodal and is mediated by their electric sense, in which species-specific electric organ discharges (EODs) are generated in a context-dependent and thus variable sequence of pulse intervals (SPI). While the primary function of the electric sense is considered to be electrolocation, both of its components likely carry information regarding identity of the sender. However, a clear understanding of their contribution to species recognition is incomplete. We therefore analysed these two electrocommunication components (EOD waveform and SPI statistics) in two sympatric mormyrid Campylomormyrus species. In a set of five playback conditions, we further investigated which components may drive interspecific recognition and discrimination. While we found that both electrocommunication components are species-specific, the cues necessary for species recognition differ between the two species studied. While the EOD waveform and SPI were both necessary and sufficient for species recognition in C. compressirostris males, C. tamandua males apparently utilize other, non-electric modalities. Mapped onto a recent phylogeny, our results suggest that discrimination by electric cues alone may be an apomorphic trait evolved during a recent radiation in this taxon.
Communication is key to a wide variety of animal behaviours and multiple modalities are often involved in this exchange of information from sender to receiver. The communication of African weakly electric fish, however, is thought to be predominantly unimodal and is mediated by their electric sense, in which species-specific electric organ discharges (EODs) are generated in a context-dependent and thus variable sequence of pulse intervals (SPI). While the primary function of the electric sense is considered to be electrolocation, both of its components likely carry information regarding identity of the sender. However, a clear understanding of their contribution to species recognition is incomplete. We therefore analysed these two electrocommunication components (EOD waveform and SPI statistics) in two sympatric mormyrid Campylomormyrus species. In a set of five playback conditions, we further investigated which components may drive interspecific recognition and discrimination. While we found that both electrocommunication components are species-specific, the cues necessary for species recognition differ between the two species studied. While the EOD waveform and SPI were both necessary and sufficient for species recognition in C. compressirostris males, C. tamandua males apparently utilize other, non-electric modalities. Mapped onto a recent phylogeny, our results suggest that discrimination by electric cues alone may be an apomorphic trait evolved during a recent radiation in this taxon.
Objective
We investigated the potential role of indirect benefits for female mate preferences in a highly promiscuous species of live-bearing fishes, the sailfin molly Poecilia latipinna using an integrative approach that combines methods from animal behavior, life-history evolution, and genetics. Males of this species solely contribute sperm for reproduction, and consequently females do not receive any direct benefits. Despite this, females typically show clear mate preferences. It has been suggested that females can increase their reproductive success through indirect benefits from choosing males of higher quality.
Results
Although preferences for large body size have been recorded as an honest signal for genetic quality, this particular study resulted in female preference being unaffected by male body size. Nonetheless, larger males did sire more offspring, but with no effect on offspring quality. This study presents a methodical innovation by combining preference testing with life history measurements—such as the determination of the dry weight of fish embryos—and paternity analyses on single fish embryos.
Objective
We investigated the potential role of indirect benefits for female mate preferences in a highly promiscuous species of live-bearing fishes, the sailfin molly Poecilia latipinna using an integrative approach that combines methods from animal behavior, life-history evolution, and genetics. Males of this species solely contribute sperm for reproduction, and consequently females do not receive any direct benefits. Despite this, females typically show clear mate preferences. It has been suggested that females can increase their reproductive success through indirect benefits from choosing males of higher quality.
Results
Although preferences for large body size have been recorded as an honest signal for genetic quality, this particular study resulted in female preference being unaffected by male body size. Nonetheless, larger males did sire more offspring, but with no effect on offspring quality. This study presents a methodical innovation by combining preference testing with life history measurements—such as the determination of the dry weight of fish embryos—and paternity analyses on single fish embryos.
Relatedness strongly influences social behaviors in a wide variety of species. For most species, the highest typical degree of relatedness is between full siblings with 50% shared genes. However, this is poorly understood in species with unusually high relatedness between individuals: clonal organisms. Although there has been some investigation into clonal invertebrates and yeast, nothing is known about kin selection in clonal vertebrates. We show that a clonal fish, the Amazon molly (Poecilia formosa), can distinguish between different clonal lineages, associating with genetically identical, sister clones, and use multiple sensory modalities. Also, they scale their aggressive behaviors according to the relatedness to other females: they are more aggressive to non-related clones. Our results demonstrate that even in species with very small genetic differences between individuals, kin recognition can be adaptive. Their discriminatory abilities and regulation of costly behaviors provides a powerful example of natural selection in species with limited genetic diversity.
Background
The unisexual Amazon molly (Poecilia formosa) originated from a hybridization between two sexual species, the sailfin molly (Poecilia latipinna) and the Atlantic molly (Poecilia mexicana). The Amazon molly reproduces clonally via sperm-dependent parthenogenesis (gynogenesis), in which the sperm of closely related species triggers embryogenesis of the apomictic oocytes, but typically does not contribute genetic material to the next generation. We compare for the first time the gonadal transcriptome of the Amazon molly to those of both ancestral species, P. mexicana and P. latipinna.
Results
We sequenced the gonadal transcriptomes of the P. formosa and its parental species P. mexicana and P. latipinna using Illumina RNA-sequencing techniques (paired-end, 100 bp). De novo assembly of about 50 million raw read pairs for each species was performed using Trinity, yielding 106,922 transcripts for P. formosa, 115,175 for P. latipinna, and 133,025 for P. mexicana after eliminating contaminations. On the basis of sequence similarity comparisons to other teleost species and the UniProt databases, functional annotation, and differential expression analysis, we demonstrate the similarity of the transcriptomes among the three species. More than 40% of the transcripts for each species were functionally annotated and about 70% were assigned to orthologous genes of a closely related species. Differential expression analysis between the sexual and unisexual species uncovered 2035 up-regulated and 564 down-regulated genes in P. formosa. This was exemplary validated for six genes by qRT-PCR.
Conclusions
We identified more than 130 genes related to meiosis and reproduction within the apomictically reproducing P. formosa. Overall expression of these genes seems to be down-regulated in the P. formosa transcriptome compared to both ancestral species (i.e., 106 genes down-regulated, 29 up-regulated). A further 35 meiosis and reproduction related genes were not found in the P. formosa transcriptome, but were only expressed in the sexual species. Our data support the hypothesis of general down-regulation of meiosis-related genes in the apomictic Amazon molly. Furthermore, the obtained dataset and identified gene catalog will serve as a resource for future research on the molecular mechanisms behind the reproductive mode of this unisexual species.
Abstract
Background
The unisexual Amazon molly (Poecilia formosa) originated from a hybridization between two sexual species, the sailfin molly (Poecilia latipinna) and the Atlantic molly (Poecilia mexicana). The Amazon molly reproduces clonally via sperm-dependent parthenogenesis (gynogenesis), in which the sperm of closely related species triggers embryogenesis of the apomictic oocytes, but typically does not contribute genetic material to the next generation. We compare for the first time the gonadal transcriptome of the Amazon molly to those of both ancestral species, P. mexicana and P. latipinna.
Results
We sequenced the gonadal transcriptomes of the P. formosa and its parental species P. mexicana and P. latipinna using Illumina RNA-sequencing techniques (paired-end, 100 bp). De novo assembly of about 50 million raw read pairs for each species was performed using Trinity, yielding 106,922 transcripts for P. formosa, 115,175 for P. latipinna, and 133,025 for P. mexicana after eliminating contaminations. On the basis of sequence similarity comparisons to other teleost species and the UniProt databases, functional annotation, and differential expression analysis, we demonstrate the similarity of the transcriptomes among the three species. More than 40% of the transcripts for each species were functionally annotated and about 70% were assigned to orthologous genes of a closely related species. Differential expression analysis between the sexual and unisexual species uncovered 2035 up-regulated and 564 down-regulated genes in P. formosa. This was exemplary validated for six genes by qRT-PCR.
Conclusions
We identified more than 130 genes related to meiosis and reproduction within the apomictically reproducing P. formosa. Overall expression of these genes seems to be down-regulated in the P. formosa transcriptome compared to both ancestral species (i.e., 106 genes down-regulated, 29 up-regulated). A further 35 meiosis and reproduction related genes were not found in the P. formosa transcriptome, but were only expressed in the sexual species. Our data support the hypothesis of general down-regulation of meiosis-related genes in the apomictic Amazon molly. Furthermore, the obtained dataset and identified gene catalog will serve as a resource for future research on the molecular mechanisms behind the reproductive mode of this unisexual species.
Background
The Amazon molly, Poecilia formosa (Teleostei: Poeciliinae) is an unisexual, all-female species. It evolved through the hybridisation of two closely related sexual species and exhibits clonal reproduction by sperm dependent parthenogenesis (or gynogenesis) where the sperm of a parental species is only used to activate embryogenesis of the apomictic, diploid eggs but does not contribute genetic material to the offspring.
Here we provide and describe the first de novo assembled transcriptome of the Amazon molly in comparison with its maternal ancestor, the Atlantic molly Poecilia mexicana. The transcriptome data were produced through sequencing of single end libraries (100 bp) with the Illumina sequencing technique.
Results
83,504,382 reads for the Amazon molly and 81,625,840 for the Atlantic molly were assembled into 127,283 and 78,961 contigs for the Amazon molly and the Atlantic molly, respectively. 63% resp. 57% of the contigs could be annotated with gene ontology terms after sequence similarity comparisons. Furthermore, we were able to identify genes normally involved in reproduction and especially in meiosis also in the transcriptome dataset of the apomictic reproducing Amazon molly.
Conclusions
We assembled and annotated the transcriptome of a non-model organism, the Amazon molly, without a reference genome (de novo). The obtained dataset is a fundamental resource for future research in functional and expression analysis. Also, the presence of 30 meiosis-specific genes within a species where no meiosis is known to take place is remarkable and raises new questions for future research.
Background
Natural accessions of Arabidopsis thaliana are a well-known system to measure levels of intraspecific genetic variation. Leaf starch content correlates negatively with biomass. Starch is synthesized by the coordinated action of many (iso)enzymes. Quantitatively dominant is the repetitive transfer of glucosyl residues to the non-reducing ends of α-glucans as mediated by starch synthases. In the genome of A. thaliana, there are five classes of starch synthases, designated as soluble starch synthases (SSI, SSII, SSIII, and SSIV) and granule-bound synthase (GBSS). Each class is represented by a single gene. The five genes are homologous in functional domains due to their common origin, but have evolved individual features as well. Here, we analyze the extent of genetic variation in these fundamental protein classes as well as possible functional implications on transcript and protein levels.
Findings
Intraspecific sequence variation of the five starch synthases was determined by sequencing the entire loci including promoter regions from 30 worldwide distributed accessions of A. thaliana. In all genes, a considerable number of nucleotide polymorphisms was observed, both in non-coding and coding regions, and several amino acid substitutions were identified in functional domains. Furthermore, promoters possess numerous polymorphisms in potentially regulatory cis-acting regions. By realtime experiments performed with selected accessions, we demonstrate that DNA sequence divergence correlates with significant differences in transcript levels.
Conclusions
Except for AtSSII, all starch synthase classes clustered into two or three groups of haplotypes, respectively. Significant difference in transcript levels among haplotype clusters in AtSSIV provides evidence for cis-regulation. By contrast, no such correlation was found for AtSSI, AtSSII, AtSSIII, and AtGBSS, suggesting trans-regulation. The expression data presented here point to a regulation by common trans-regulatory transcription factors which ensures a coordinated action of the products of these four genes during starch granule biosynthesis. The apparent cis-regulation of AtSSIV might be related to its role in the initiation of de novo biosynthesis of granules.
Introduction: We examined patterns of genetic divergence in 26 Mediterranean populations of the semi-terrestrial beachflea Orchestia montagui using mitochondrial (cytochrome oxidase subunit I), microsatellite (eight loci) and allozymic data. The species typically forms large populations within heaps of dead seagrass leaves stranded on beaches at the waterfront. We adopted a hierarchical geographic sampling to unravel population structure in a species living at the sea-land transition and, hence, likely subjected to dramatically contrasting forces.
Results: Mitochondrial DNA showed historical phylogeographic breaks among Adriatic, Ionian and the remaining basins (Tyrrhenian, Western and Eastern Mediterranean Sea) likely caused by the geological and climatic changes of the Pleistocene. Microsatellites (and to a lesser extent allozymes) detected a further subdivision between and within the Western Mediterranean and the Tyrrhenian Sea due to present-day processes. A pattern of isolation by distance was not detected in any of the analyzed data set.
Conclusions: We conclude that the population structure of O. montagui is the result of the interplay of two contrasting forces that act on the species population genetic structure. On one hand, the species semi-terrestrial life style would tend to determine the onset of local differences. On the other hand, these differences are partially counter-balanced by passive movements of migrants via rafting on heaps of dead seagrass leaves across sites by sea surface currents. Approximate Bayesian Computations support dispersal at sea as prevalent over terrestrial regionalism.
The all-female Amazon molly (Poecilia formosa) is the result of a hybridization of the Atlantic molly (P. mexicana) and the sailfin molly (P. latipinna) approximately 120,000 years ago. As a gynogenetic species, P. formosa needs to copulate with heterospecific males including males from one of its bisexual ancestral species. However, the sperm only triggers embryogenesis of the diploid eggs. The genetic information of the sperm donor typically will not contribute to the next generation of P. formosa. Hence, P. formosa possesses generally one allele from each of its ancestral species at any genetic locus. This raises the question whether both ancestral alleles are equally expressed in P. formosa. Allele-specific expression (ASE) has been previously assessed in various organisms, e.g., human and fish, and ASE was found to be important in the context of phenotypic variability and disease. In this study, we utilized Real-Time PCR techniques to estimate ASE of the androgen receptor alpha (arα) gene in several distinct tissues of Amazon mollies. We found an allelic bias favoring the maternal ancestor (P. mexicana) allele in ovarian tissue. This allelic bias was not observed in the gill or the brain tissue. Sequencing of the promoter regions of both alleles revealed an association between an Indel in a known CpG island and differential expression. Future studies may reveal whether our observed cis-regulatory divergence is caused by an ovary-specific trans-regulatory element, preferentially activating the allele of the maternal ancestor.