Allele-specific expression at the androgen receptor alpha gene in a hybrid unisexual fish, the Amazon molly (Poecilia formosa)

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

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

Identification and effective regulation of scarb1 gene involved in pigmentation change in autotetraploid Carassius auratus

The autotetraploid Carassius auratus (4nRR, 4 n=200, RRRR) is derived from whole-genome duplication of Carassius auratus red var. (RCC, 2 n=100, RR). In the current study, we demonstrated that chromatophores and pigment changes directly caused the coloration and variation of 4nRR skin (red in RCC, brownish-yellow in 4nRR). To further explore the molecular mechanisms underlying coloration formation and variation in 4nRR, we performed transcriptome profiling and molecular functional verification in RCC and 4nRR. Results revealed that scarb1, associated with carotenoid metabolism, underwent significant down-regulation in 4nRR. Efficient editing of this candidate pigment gene provided clear evidence of its significant role in RCC coloration. Subsequently, we identified four divergent scarb1 homeologs in 4nRR: two original scarb1 homeologs from RCC and two duplicated ones. Notably, three of these homeologs possessed two highly conserved alleles, exhibiting biased and allele-specific expression in the skin. Remarkably, after precise editing of both the original and duplicated scarb1 homeologs and/or alleles, 4nRR individuals, whether singly or multiply mutated, displayed a transition from brownish-yellow skin to a cyan-gray phenotype. Concurrently, the proportional areas of the cyan-gray regions displayed a gene-dose correlation. These findings illustrate the subfunctionalization of duplicated scarb1, with all scarb1 genes synergistically and equally contributing to the pigmentation of 4nRR. This is the first report concerning the functional differentiation of duplicated homeologs in an autopolyploid fish, substantially enriching our understanding of coloration formation and change within this group of organisms.

Non-Additive and Asymmetric Allelic Expression of p38 mapk in Hybrid Tilapia (Oreochromis niloticus ♀ × O. aureus ♂)

Hybridization is a widely used breeding technique in fish species that enhances desirable traits in cultured species through heterosis. However, the mechanism by which hybrids alter gene expression to form heterosis remains unclear. In this study, a group of hybrid tilapia was used to elucidate heterosis through interspecies crossing. Specifically, p38 was analyzed to describe the regulation of gene expression variation in hybrid tilapia. Transcripts from the Nile tilapia allele were found to be significantly higher than those from the blue tilapia allele in hybrid individuals, indicating that the expression of p38 was dominated by Nile tilapia sub-genomic alleles. The study also found a compensatory interaction of cis- and trans-acting elements of the Nile tilapia and blue tilapia sub-genomes, inducing a non-additive expression of p38 in hybrids. Eight specific SNPs were identified in the p38 promoter regions of Nile tilapia and blue tilapia, and were found to be promoter differences leading to differences in gene expression efficiencies between parental alleles using a dual-luciferase reporter system. This study provides insights into the non-additive expression patterns of key functional genes in fish hybrids related to growth and immunity response.

Nonadditive and Asymmetric Allelic Expression of Growth Hormone in Hybrid Tilapia

Hybridization is a common breeding technique that can improve germplasm through heterosis in aquaculture. However, the regulation of key gene expression, including the details of transcriptional level changes at the beginning of hybridization events, remains largely undefined, especially in teleosts. In this study, by interspecies crossing between two pure lines of Nile tilapia and blue tilapia, we obtained a hybrid tilapia population as a model to elucidate heterosis, and we traced the molecular outcomes of growth hormone (GH) expression and allele-specific expression (ASE) in hybrids. The hybrids display growth vigor compared to their parents in the 120-day growth trial. GH mRNA expression was uniquely expressed in the pituitary. Higher GH expression was found in the hybrid than the midparent value, in both males and females, showing a nonadditive pattern. We identified four single-nucleotide polymorphism sites between Nile tilapia and blue tilapia. Subsequently, by pyrosequencing, we found asymmetric allelic expression in hybrids with higher maternal allelic transcript ratios in both males and females. Fasting significantly increased GH expression in hybrids, but asymmetric allelic expression was not affected by feeding or fasting conditions. Finally, we identified cis and trans effects via overall expression and ASE values in the hybrid, which showed that the cis and trans effects promoted the expression of maternal subgenome in the hybrid, contributing to the expression superiority of GH in hybrid tilapia. Taken together, the results of our study first illustrated the concept of GH expression superiority and its formation mechanism in hybrid fish with growth vigor.

Nonadditive and allele-specific expression of insulin-like growth factor 1 in Nile tilapia (Oreochromis niloticus, ♀) × blue tilapia (O. aureus, ♂) hybrids

Hybrid Nile tilapia (Oreochromis niloticus, ♀) × blue tilapia (O. aureus, ♂) is a widely cultured tilapia variety due to its growth vigor compared to the parent species. As a peptide hormone, insulin-like growth factor 1 (IGF-1) plays a critical role in regulating somatic growth. The present study focuses on the expression characteristics of IGF-1 in hybrid tilapia. The cloned complete open reading frame of IGF-1 in hybrid tilapia is 549 bp in length, encoding a protein of 182 amino acids. The deduced protein is highly similar to that of Nile tilapia and blue tilapia. IGF-1 was found to be primarily expressed in the liver and muscle in the hybrid; lower expression levels were found in other tissues such as the intestine, spleen, and head-kidney. Increased mRNA expression was observed in the liver and muscle of the hybrid compared to Nile tilapia and blue tilapia, indicating a nonadditive expression pattern in the hybrid. An IGF-1 SNP site (397 site: C in Nile tilapia, G in blue tilapia) for differentiating the Nile tilapia or blue tilapia subgenome in hybrids was identified. Pyrosequencing analysis of the liver transcriptome indicated that most of the hybrids (9 of 10 individuals) predominantly expressed the G allele, demonstrating bias of the blue tilapia subgenome. The present study provides novel data indicating, for the first time, overall gene expression of IGF-1 and allele-specific expression in hybrid tilapia.

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

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.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4382-2) contains supplementary material, which is available to authorized users.