Identification of transcriptome SNPs between Xiphophorus lines and species for assessing allele specific gene expression within F₁ interspecies hybrids

Maternal patterns of inheritance alter transcript expression in eggs

BACKGROUND

Modifications to early development can lead to evolutionary diversification. The early stages of development are under maternal control, as mothers produce eggs loaded with nutrients, proteins and mRNAs that direct early embryogenesis. Maternally provided mRNAs are the only expressed genes in initial stages of development and are tightly regulated. Differences in maternal mRNA provisioning could lead to phenotypic changes in embryogenesis and ultimately evolutionary changes in development. However, the extent that maternal mRNA expression in eggs can vary is unknown for most developmental models. Here, we use a species with dimorphic development- where females make eggs and larvae of different sizes and life-history modes-to investigate the extent of variation in maternal mRNA provisioning to the egg.

RESULTS

We find that there is significant variation in gene expression across eggs of different development modes, and that there are both qualitative and quantitative differences in mRNA expression. We separate parental effects from allelic effects, and find that both mechanisms contribute to mRNA expression differences. We also find that offspring of intraspecific crosses differentially provision their eggs based on the parental cross direction (a parental effect), which has not been previously demonstrated in reproductive traits like oogenesis.

CONCLUSION

We find that maternally controlled initiation of development is functionally distinct between eggs of different sizes and maternal genotypes. Both allele-specific effects and parent-of-origin effects contribute to gene expression differences in eggs. The latter indicates an intergenerational effect where a parent's genotype can affect gene expression in an egg made by the next generation.

Unraveling regulatory divergence, heterotic malleability, and allelic imbalance switching in rice due to drought stress

The indica ecotypes, IR64, an elite drought-susceptible variety adapted to irrigated ecosystem, and Apo (IR55423-01 or NSIC RC9), a moderate drought-tolerant upland genotype together with their hybrid (IR64 × Apo) were exposed to non- and water-stress conditions. By sequencing (RNA-seq) these genotypes, we were able to map genes diverging in cis and/or trans factors. Under non-stress condition, cis dominantly explains (11.2%) regulatory differences, followed by trans (8.9%). Further analysis showed that water-limiting condition largely affects trans and cis + trans factors. On the molecular level, cis and/or trans regulatory divergence explains their genotypic differences and differential drought response. Between the two parental genotypes, Apo appears to exhibit more photosynthetic efficiency even under water-limiting condition and is ascribed to trans. Statistical analyses showed that regulatory divergence is significantly influenced by environmental conditions. Likewise, the mode of parental expression inheritance which drives heterosis (HET) is significantly affected by environmental conditions indicating the malleability of heterosis to external factors. Further analysis revealed that the HET class, dominance, was significantly enriched under water-stress condition. We also identified allelic imbalance switching in which several genes prefer IR64- (or Apo-) specific allele under non-stress condition but switched to Apo- (or IR64-) specific allele when exposed to water-stress condition.

Allele-specific expression variation at different ploidy levels in Squalius alburnoides

Allopolyploid plants are long known to be subject to a homoeolog expression bias of varying degree. The same phenomenon was only much later suspected to occur also in animals based on studies of single selected genes in an allopolyploid vertebrate, the Iberian fish Squalius alburnoides. Consequently, this species became a good model for understanding the evolution of gene expression regulation in polyploid vertebrates. Here, we analyzed for the first time genome-wide allele-specific expression data from diploid and triploid hybrids of S. alburnoides and compared homoeolog expression profiles of adult livers and of juveniles. Co-expression of alleles from both parental genomic types was observed for the majority of genes, but with marked homoeolog expression bias, suggesting homoeolog specific reshaping of expression level patterns in hybrids. Complete silencing of one allele was also observed irrespective of ploidy level, but not transcriptome wide as previously speculated. Instead, it was found only in a restricted number of genes, particularly ones with functions related to mitochondria and ribosomes. This leads us to hypothesize that allelic silencing may be a way to overcome intergenomic gene expression interaction conflicts, and that homoeolog expression bias may be an important mechanism in the achievement of sustainable genomic interactions, mandatory to the success of allopolyploid systems, as in S. alburnoides.

Transcriptome signature of liver tissue with divergent mutton odour and flavour using RNA deep sequencing

Mutton consumption is less popular in many Asian countries including Indonesia, whose consumers often complain about the unpleasant flavour and odour of the meat. The main causes of mutton odour are the two compounds of branched chain fatty acid (BCFA): methylnonanoic (MNA), phenol, 3-methyl (MP), 4-methylnonanoic (MNA) and 4-ethyloctanoic (EOA) present in all the adipose tissue; and the 3-methylindole (MI) or skatole and indole, which are originated from pastoral diets. It is crucial to understand the genetic mechanism of mutton odour and flavour (MOF) to select sheep for lower BCFA and indole thus reduce the unpleasant flavour of meat. The aim of the present study was to investigate transcriptome profiling in liver tissue with divergent MOF using RNA deep sequencing. Liver tissues from higher (n = 3) and lower (n = 3) MOF sheep were analysed using Illumina HiSeq 2500. The total number of reads produced for each liver sample ranged from 21.37 to 25.37 million. Approximately 103 genes were differentially expressed (DEGs) with significance level of p-adjusted value <0.05. Among them, 60 genes were up-regulated, and 43 were down-regulated (p < 0.01, FC > 1.5) in higher MOF group. Differentially regulated genes in high MOF liver samples were enriched in biological processes such as cellular response to chemical stimulus and endogenous stimulus; cellular components such as such as basement membrane and extracellular matrix; and molecular functions such as haeme binding and oxidoreductase activity. Among the DEGs, metabolic phase I related genes belonging to the cytochrome P450 CYP2A6 were dominantly expressed. Additionally, phase II conjugation genes including UDP glucuronosyltransferases UGT2B18, sulfotransferase SULT1C1, and glutathione S-transferase GSTM1 were identified. The dominant candidate genes for SOF could be cytochrome P450, sodium-channel protein, transmembrane protein, glutathione transferase, UDP glucuronosyltransferases and sulfotransferase. Pathway analysis identified steroid hormone biosynthesis and chemical carcinogenesis by cytochrome P450 pathways which may play important roles in MOF-related molecules metabolism. This work highlighted potential genes and gene-networks that may affect meat off flavour and odour in sheep.

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

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.

A second generation SNP and SSR integrated linkage map and QTL mapping for the Chinese mitten crab Eriocheir sinensis

The Chinese mitten crab Eriocheir sinensis is the most economically important cultivated crab species in China, and its genome has a high number of chromosomes (2n = 146). To obtain sufficient markers for construction of a dense genetic map for this species, we employed the recently developed specific-locus amplified fragment sequencing (SLAF-seq) method for large-scale SNPs screening and genotyping in a F1 full-sib family of 149 individuals. SLAF-seq generated 127,677 polymorphic SNP markers, of which 20,803 valid markers were assigned into five segregation types and were used together with previous SSR markers for linkage map construction. The final integrated genetic map included 17,680 SNP and 629 SSR markers on the 73 linkage groups (LG), and spanned 14,894.9 cM with an average marker interval of 0.81 cM. QTL mapping localized three significant growth-related QTL to a 1.2 cM region in LG53 as well as 146 sex-linked markers in LG48. Genome-wide QTL-association analysis further identified four growth-related QTL genes named LNX2, PAK2, FMRFamide and octopamine receptors. These genes are involved in a variety of different signaling pathways including cell proliferation and growth. The map and SNP markers described here will be a valuable resource for the E. sinensis genome project and selective breeding programs.

Gene copy silencing and DNA methylation in natural and artificially produced allopolyploid fish

Allelic silencing is an important mechanism for coping with gene dosage changes in polyploid organisms that is well known in allopolyploid plants. Only recently, it was shown in the allotriploid fish Squalius alburnoides that this process also occurs in vertebrates. However, it is still unknown whether this silencing mechanism is common to other allopolyploid fish, and which mechanisms might be responsible for allelic silencing. We addressed these questions in a comparative study between Squalius alburnoides and another allopolyploid complex, the Amazon molly (Poecilia formosa). We examined the allelic expression patterns for three target genes in four somatic tissues of natural allo-anorthoploids and laboratory-produced tri-genomic hybrids of S. alburnoides and P. formosa. Also, for both complexes, we evaluated the correlation between total DNA methylation level and the ploidy status and genomic composition of the individuals. We found that allelic silencing also occurs in other allopolyploid organisms besides the single one that was previously known. We found and discuss disparities within and between the two considered complexes concerning the pattern of allele-specific expression and DNA methylation levels. Disparities might be due to intrinsic characteristics of each genome involved in the hybridization process. Our findings also support the idea that long-term evolutionary processes have an effect on the allele expression patterns and possibly also on DNA methylation levels.

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.

Generation and analysis of expressed sequence tags (ESTs) of Camelina sativa to mine drought stress-responsive genes

Camelina sativa is an oil-producing crop belonging to the family of Brassicaceae. Due to exceptionally high content of omega fatty acid, it is commercially grown around the world as edible oil, biofuel, and animal feed. A commonly referred 'false flax' or gold-of-pleasure Camelina sativa has been interested as one of biofuel feedstocks. The species can grow on marginal land due to its superior drought tolerance with low requirement of agricultural inputs. This crop has been unexploited due to very limited transcriptomic and genomic data. Use of gene-specific molecular markers is an important strategy for new cultivar development in breeding program. In this study, Illumina paired-end sequencing technology and bioinformatics tools were used to obtain expression profiling of genes responding to drought stress in Camelina sativa BN14. A total of more than 60,000 loci were assembled, corresponding to approximately 275 K transcripts. When the species was exposed to 10 kPa drought stress, 100 kPa drought stress, and rehydrated conditions, a total of 107, 2,989, and 982 genes, respectively, were up-regulated, while 146, 3,659, and 1189 genes, respectively, were down-regulated compared to control condition. Some unknown genes were found to be highly expressed under drought conditions, together with some already reported gene families such as senescence-associated genes, CAP160, and LEA under 100 kPa soil water condition, cysteine protease, 2OG, Fe(II)-dependent oxygenase, and RAD-like 1 under rehydrated condition. These genes will be further validated and mapped to determine their function and loci. This EST library will be favorably applied to develop gene-specific molecular markers and discover genes responsible for drought tolerance in Camelina species.

The evolution of larval developmental mode: insights from hybrids between species with obligately and facultatively planktotrophic larvae

Life history characteristics play a pervasive role in the ecology and evolution of species. Transitions between feeding and non-feeding larval development have occurred many times in both terrestrial and marine phyla, however we lack a comprehensive understanding of how such shifts occur. The sea biscuits Clypeaster rosaceus and Clypeaster subdepressus employ different life history strategies (facultatively feeding larvae and obligately feeding larvae, respectively) but can hybridize. In this study, we examined the development of hybrid larvae between these two species in order to investigate the inheritance of larval developmental mode. Our results show that both reciprocal hybrid crosses developed via the feeding mode of their maternal species. However, as feeding larvae can obtain both energy and hormones from algal food, we tested how hormones alone affected development by setting up a treatment where we added exogenous thyroid hormone, but no food. In this treatment the offspring of all four crosses (two homospecific and two heterospecific crosses) were able to metamorphose without algal food. Therefore we hypothesize that although hybrid developmental mode was inherited from the maternal species, this result was not solely due to energetic constraints of egg size.

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.

Tempo and mode of regulatory evolution in Drosophila

Genetic changes affecting gene expression contribute to phenotypic divergence; thus, understanding how regulatory networks controlling gene expression change over time is critical for understanding evolution. Prior studies of expression differences within and between species have identified properties of regulatory divergence, but technical and biological differences among these studies make it difficult to assess the generality of these properties or to understand how regulatory changes accumulate with divergence time. Here, we address these issues by comparing gene expression among strains and species of Drosophila with a range of divergence times and use F₁ hybrids to examine inheritance patterns and disentangle cis- and trans-regulatory changes. We find that the fixation of compensatory changes has caused the regulation of gene expression to diverge more rapidly than gene expression itself. Specifically, we observed that the proportion of genes with evidence of cis-regulatory divergence has increased more rapidly with divergence time than the proportion of genes with evidence of expression differences. Surprisingly, the amount of expression divergence explained by cis-regulatory changes did not increase steadily with divergence time, as was previously proposed. Rather, one species (Drosophila sechellia) showed an excess of cis-regulatory divergence that we argue most likely resulted from positive selection in this lineage. Taken together, this work reveals not only the rate at which gene expression evolves, but also the molecular and evolutionary mechanisms responsible for this evolution.

RNA-Seq technology and its application in fish transcriptomics

High-throughput sequencing technologies, also known as next-generation sequencing (NGS) technologies, have revolutionized the way that genomic research is advancing. In addition to the static genome, these state-of-art technologies have been recently exploited to analyze the dynamic transcriptome, and the resulting technology is termed RNA sequencing (RNA-seq). RNA-seq is free from many limitations of other transcriptomic approaches, such as microarray and tag-based sequencing method. Although RNA-seq has only been available for a short time, studies using this method have completely changed our perspective of the breadth and depth of eukaryotic transcriptomes. In terms of the transcriptomics of teleost fishes, both model and non-model species have benefited from the RNA-seq approach and have undergone tremendous advances in the past several years. RNA-seq has helped not only in mapping and annotating fish transcriptome but also in our understanding of many biological processes in fish, such as development, adaptive evolution, host immune response, and stress response. In this review, we first provide an overview of each step of RNA-seq from library construction to the bioinformatic analysis of the data. We then summarize and discuss the recent biological insights obtained from the RNA-seq studies in a variety of fish species.

Identification of the novel candidate genes and variants in boar liver tissues with divergent skatole levels using RNA deep sequencing

Boar taint is the unpleasant odour of meat derived from non-castrated male pigs, caused by the accumulation of androstenone and skatole in fat. Skatole is a tryptophan metabolite produced by intestinal bacteria in gut and catabolised in liver. Since boar taint affects consumer's preference, the aim of this study was to perform transcriptome profiling in liver of boars with divergent skatole levels in backfat by using RNA-Seq. The total number of reads produced for each liver sample ranged from 11.8 to 39.0 million. Approximately 448 genes were differentially regulated (p-adjusted <0.05). Among them, 383 genes were up-regulated in higher skatole group and 65 were down-regulated (p<0.01, FC>1.5). Differentially regulated genes in the high skatole liver samples were enriched in metabolic processes such as small molecule biochemistry, protein synthesis, lipid and amino acid metabolism. Pathway analysis identified the remodeling of epithelial adherens junction and TCA cycle as the most dominant pathways which may play important roles in skatole metabolism. Differential gene expression analysis identified candidate genes in ATP synthesis, cytochrome P450, keratin, phosphoglucomutase, isocitrate dehydrogenase and solute carrier family. Additionally, polymorphism and association analysis revealed that mutations in ATP5B, KRT8, PGM1, SLC22A7 and IDH1 genes could be potential markers for skatole levels in boars. Furthermore, expression analysis of exon usage of three genes (ATP5B, KRT8 and PGM1) revealed significant differential expression of exons of these genes in different skatole levels. These polymorphisms and exon expression differences may have impacts on the gene activity ultimately leading to skatole variation and could be used as genetic marker for boar taint related traits. However, further validation is required to confirm the effect of these genetic markers in other pig populations in order to be used in genomic selection against boar taint in pig breeding programs.

RNaseIII and T4 polynucleotide Kinase sequence biases and solutions during RNA-seq library construction

Background

RNA-seq is a next generation sequencing method with a wide range of applications including single nucleotide polymorphism (SNP) detection, splice junction identification, and gene expression level measurement. However, the RNA-seq sequence data can be biased during library constructions resulting in incorrect data for SNP, splice junction, and gene expression studies. Here, we developed new library preparation methods to limit such biases.

Results

A whole transcriptome library prepared for the SOLiD system displayed numerous read duplications (pile-ups) and gaps in known exons. The pile-ups and gaps of the whole transcriptome library caused a loss of SNP and splice junction information and reduced the quality of gene expression results. Further, we found clear sequence biases for both 5' and 3' end reads in the whole transcriptome library. To remove this bias, RNaseIII fragmentation was replaced with heat fragmentation. For adaptor ligation, T4 Polynucleotide Kinase (T4PNK) was used following heat fragmentation. However, its kinase and phosphatase activities introduced additional sequence biases. To minimize them, we used OptiKinase before T4PNK. Our study further revealed the specific target sequences of RNaseIII and T4PNK.

Conclusions

Our results suggest that the heat fragmentation removed the RNaseIII sequence bias and significantly reduced the pile-ups and gaps. OptiKinase minimized the T4PNK sequence biases and removed most of the remaining pile-ups and gaps, thus maximizing the quality of RNA-seq data.

Reviewers

This article was reviewed by Dr. A. Kolodziejczyk (nominated by Dr. Sarah Teichmann), Dr. Eugene Koonin, and Dr. Christoph Adami. For the full reviews, see the Reviewers' Comments section.

Allelic imbalance metre (Allim), a new tool for measuring allele-specific gene expression with RNA-seq data

Estimating differences in gene expression among alleles is of high interest for many areas in biology and medicine. Here, we present a user-friendly software tool, Allim, to estimate allele-specific gene expression. Because mapping bias is a major problem for reliable estimates of allele-specific gene expression using RNA-seq, Allim combines two different strategies to account for the mapping biases. In order to reduce the mapping bias, Allim first generates a polymorphism-aware reference genome that accounts for the sequence variation between the alleles. Then, a sequence-specific simulation tool estimates the residual mapping bias. Statistical tests for allelic imbalance are provided that can be used with the bias corrected RNA-seq data.

Identification of transcriptome SNPs for assessing allele-specific gene expression in a super-hybrid rice Xieyou9308

Hybridization, a common process in nature, can give rise to a vast reservoir of allelic variants. Combination of these allelic variants may result in novel patterns of gene action and is thought to contribute to heterosis. In this study, we analyzed genome-wide allele-specific gene expression (ASGE) in the super-hybrid rice variety Xieyou9308 using RNA sequencing technology (RNA-Seq). We identified 9325 reliable single nucleotide polymorphisms (SNPs) distributed throughout the genome. Nearly 68% of the identified polymorphisms were CT and GA SNPs between R9308 and Xieqingzao B, suggesting the existence of DNA methylation, a heritable epigenetic mark, in the parents and their F1 hybrid. Of 2793 identified transcripts with consistent allelic biases, only 480 (17%) showed significant allelic biases during tillering and/or heading stages, implying that trans effects may mediate most transcriptional differences in hybrid offspring. Approximately 67% and 62% of the 480 transcripts showed R9308 allelic expression biases at tillering and heading stages, respectively. Transcripts with higher levels of gene expression in R9308 also exhibited R9308 allelic biases in the hybrid. In addition, 125 transcripts were identified with significant allelic expression biases at both stages, of which 74% showed R9308 allelic expression biases. R9308 alleles may tend to preserve their characteristic states of activity in the hybrid and may play important roles in hybrid vigor at both stages. The allelic expression of 355 transcripts was highly stage-specific, with divergent allelic expression patterns observed at different developmental stages. Many transcripts associated with stress resistance were differently regulated in the F1 hybrid. The results of this study may provide valuable insights into molecular mechanisms of heterosis.

The evolutionary history of Xiphophorus fish and their sexually selected sword: a genome-wide approach using restriction site-associated DNA sequencing

Next-generation sequencing (NGS) techniques are now key tools in the detection of population genomic and gene expression differences in a large array of organisms. However, so far few studies have utilized such data for phylogenetic estimations. Here, we use NGS data obtained from genome-wide restriction site-associated DNA (RAD) (∼66000 SNPs) to estimate the phylogenetic relationships among all 26 species of swordtail and platyfish (genus Xiphophorus) from Central America. Past studies, both sequence and morphology-based, have differed in their inferences of the evolutionary relationships within this genus, particularly at the species-level and among monophyletic groupings. We show that using a large number of markers throughout the genome, we are able to infer the phylogenetic relationships with unparalleled resolution for this genus. The relationships among all three major clades and species within each of them are highly resolved and consistent under maximum likelihood, Bayesian inference and maximum parsimony. However, we also highlight the current cautions with this data type and analyses. This genus exhibits a particularly interesting evolutionary history where at least two species may have arisen through hybridization events. Here, we are able to infer the paternal lineages of these putative hybrid species. Using the RAD-marker-based tree we reconstruct the evolutionary history of the sexually selected sword trait and show that it may have been present in the common ancestor of the genus. Together our results highlight the outstanding capacity that RAD sequencing data has for resolving previously problematic phylogenetic relationships, particularly among relatively closely related species.

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.

Linking transcriptomic and genomic variation to growth in brook charr hybrids (Salvelinus fontinalis, Mitchill)

Hybridization can lead to phenotypic differences arising from changes in gene expression patterns or new allele combinations. Variation in gene expression is thought to be controlled by differences in transcription regulation of parental alleles, either through cis- or trans-regulatory elements. A previous study among brook charr hybrids from different populations (Rupert, Laval, and domestic) showing distinct length at age during early life stages also revealed different patterns in transcription regulation inheritance of transcript abundance. In the present study, transcript abundance using RNA-sequencing and quantitative real-time PCR, single-nucleotide polymorphism (SNP) genotypes and allelic imbalance were assessed in order to understand the molecular mechanisms underlying the observed transcriptomic and differences in length at age among domestic × Rupert hybrids and Laval × domestic hybrids. We found 198 differentially expressed genes between the two hybrid crosses, and allelic imbalance could be analyzed for 69 of them. Among these 69 genes, 36 genes exhibited cis-acting regulatory effects in both of the two crosses, thus confirming the prevalent role of cis-acting regulatory elements in the regulation of differentially expressed genes among intraspecific hybrids. In addition, we detected a significant association between SNP genotypes of three genes and length at age. Our study is thus one of the few that have highlighted some of the molecular mechanisms potentially involved in the differential phenotypic expression in intraspecific hybrids for nonmodel species.

Genomic resources for a model in adaptation and speciation research: characterization of the Poecilia mexicana transcriptome

Background

Elucidating the genomic basis of adaptation and speciation is a major challenge in natural systems with large quantities of environmental and phenotypic data, mostly because of the scarcity of genomic resources for non-model organisms. The Atlantic molly (Poecilia mexicana, Poeciliidae) is a small livebearing fish that has been extensively studied for evolutionary ecology research, particularly because this species has repeatedly colonized extreme environments in the form of caves and toxic hydrogen sulfide containing springs. In such extreme environments, populations show strong patterns of adaptive trait divergence and the emergence of reproductive isolation. Here, we used RNA-sequencing to assemble and annotate the first transcriptome of P. mexicana to facilitate ecological genomics studies in the future and aid the identification of genes underlying adaptation and speciation in the system.

Description

We provide the first annotated reference transcriptome of P. mexicana. Our transcriptome shows high congruence with other published fish transcriptomes, including that of the guppy, medaka, zebrafish, and stickleback. Transcriptome annotation uncovered the presence of candidate genes relevant in the study of adaptation to extreme environments. We describe general and oxidative stress response genes as well as genes involved in pathways induced by hypoxia or involved in sulfide metabolism. To facilitate future comparative analyses, we also conducted quantitative comparisons between P. mexicana from different river drainages. 106,524 single nucleotide polymorphisms were detected in our dataset, including potential markers that are putatively fixed across drainages. Furthermore, specimens from different drainages exhibited some consistent differences in gene regulation.

Conclusions

Our study provides a valuable genomic resource to study the molecular underpinnings of adaptation to extreme environments in replicated sulfide spring and cave environments. In addition, this study adds to the increasing number of genomic resources in the family Poeciliidae, which are widely used in comparative analyses of behavior, ecology, evolution, and medical genetics.

Alternative strategies for development of a reference transcriptome for quantification of allele specific expression in organisms having sparse genomic resources

In recent years RNA-Seq technology has been used not only to quantify differences in gene expression but also to understand the underlying mechanisms that lead to these differences. Nucleotide sequence variation arising through evolution may differentially affect the expression profiles of divergent species. RNA-Seq technology, combined with techniques to differentiate parental alleles and quantify their abundance, have recently become popular methods for allele specific gene expression (ASGE) analyses. However, analysis of gene expression within interspecies hybrids may be difficult when one of the two parental genomes represented in the hybrid does not have robust genomic resources or available transcriptome data. Herein, we compare two strategies for analyzing allele specific expression within interspecies hybrids produced from crossing two Xiphophorus fish species. The first strategy relies upon a robust reference transcriptome assembly from one species followed by identification of SNPs and creation of an in silico reference transcriptome for the second species. The second strategy employs de novo assembly of reference transcriptomes for both parental species followed by identification of homologous transcripts prior to mapping hybrid reads to a combined hybrid reference. Our results show that, although both methods are able to achieve balanced allelic distribution upon read mapping of F(1) hybrid fish transcriptomes, the second "de novo" assembly approach is superior for ASGE analyses and leads to results more consistent with those found from quantitative real time PCR assessment of gene expression. In addition, our analysis indicates that indels between the two parental alleles are the major cause of the differences in results observed when employing these two methods.

A phylogeographic investigation of the hybrid origin of a species of swordtail fish from Mexico

Hybrid speciation may contribute significantly to generating biodiversity, but only a few well-documented examples for it exist so far that do not involve polyploidization as a mechanism. The swordtail fish, Xiphophorus clemenciae, shows common hallmarks of a hybrid origin and still overlaps in its current geographic distribution with its putative ancestral species (Xiphophorus hellerii and Xiphophorus maculatus). Xiphophorus clemenciae provides an ideal system for investigating the possible continued genetic interactions between a hybrid and its parental species. Here, we use microsatellite and mitochondrial markers to investigate the population structure of these species of swordtails and search for signs of recent hybridization. Individuals were sampled from 21 localities across the known range of X. clemenciae- the Isthmus of Tehuantepec (IT) Mexico, and several environmental parameters that might represent barriers to dispersal were recorded. The hybridization event that gave rise to X. clemenciae appears to be rather ancient, and a single origin is likely. We find negligible evidence for ongoing hybridization and introgression between the putative ancestral species, because they now occupy distinct ecological niches, and a common haplotype is shared by most populations of X. clemenciae. The population structure within these species shows an isolation-by-distance (IBD) pattern and genetic differentiation between most populations is significant and high. We infer that tectonic evolution in the Isthmus has greatly restricted gene flow between the southern and central IT populations of X. clemenciae and X. helleriii and provide preliminary information to aid in conservation management of this geographically restricted hybrid species, X. clemenciae.

Production of F₁ offspring with vitrified sperm from a live-bearing fish, the green swordtail Xiphophorus hellerii

This study reports the first production of offspring with vitrified sperm from a live-bearing fish Xiphophorus hellerii. The overall goal of this study was to develop streamlined protocols for integration into a standardized approach for vitrification of aquatic species germplasm. The objectives were to (1) estimate acute toxicity of cryoprotectants, (2) evaluate vitrification solutions, (3) compare different thawing methods, (4) evaluate membrane integrity of post-thaw sperm vitrified in different cryoprotectants, and (5) evaluate the fertility of vitrified sperm. Nine cryoprotectants and two commercial vitrification additives were tested for acute toxicity and glass forming ability, alone and in combination. Two vitrification solutions, 40% glycerol (Gly) and 20% Gly+20% ethylene glycol (EG) in 500 mOsmol/kg Hanks' balanced salt solution (HBSS), were selected for vitrification of 10 μL sperm samples using inoculating loops plunged into liquid nitrogen. Samples were thawed at 24°C (one loop in 5 μL of HBSS or three loops in 500 μL of HBSS). Samples thawed in 500 μL were concentrated by centrifugation (1000 g for 5 min at 4°C) into 5 μL for artificial insemination. Offspring were produced from virgin females inseminated with sperm vitrified with 20% Gly+20% EG and concentrated by centrifugation.