RNA-seq reveals transcriptomic shock involving transposable elements reactivation in hybrids of young lake whitefish species.

Identifying the molecular basis of reproductive isolation among diverging lineages represents an essential step toward understanding speciation in natural populations. Postzygotic barriers can lead to hybrid breakdown, a syndrome that has been documented in several systems, potentially involving the reactivation of transposable elements. In northeastern North America, two lake whitefish lineages have repeatedly colonized postglacial lakes ~12,000 years ago, and a dwarf limnetic species has evolved multiple times from the normal benthic species. Reproductive isolation is incomplete between them; viable hybrids can be generated in the laboratory but significant mortality occurs and is associated with a malformed phenotype in backcross embryos, thus revealing a hybrid breakdown syndrome. By means of RNA-seq analyses, the objective of this study was to determine which genes were misregulated in hybrids and rigorously test the hypothesis of transposable element reactivation. We compared the transcriptomic landscape in pure embryos, F1-hybrids, and healthy and malformed backcrosses at the late embryonic stage. Extensive expression differences consistent with previously documented adaptive divergence between pure normal and dwarf embryos were identified for the first time. Pronounced transcriptome-wide deregulation in malformed backcrosses was observed, with over 15% of transcripts differentially expressed in all comparisons, compared with 1.5% between pure parental forms. Convincing evidence of transposable elements and noncoding transcripts reactivation in malformed backcrosses is presented. We propose that hybrid breakdown likely results from extensive genomic incompatibilities, plausibly encompassing transposable elements. Combined with previous studies, these results reveal synergy among many reproductive barriers, thus maintaining divergence between these two young whitefish species.