Regulatory and evolutionary impact of DNA methylation in two songbird species and their naturally occurring F1 hybrids

BACKGROUND

Regulation of transcription by DNA methylation in 5'-CpG-3' context is a widespread mechanism allowing differential expression of genetically identical cells to persist throughout development. Consequently, differences in DNA methylation can reinforce variation in gene expression among cells, tissues, populations, and species. Despite a surge in studies on DNA methylation, we know little about the importance of DNA methylation in population differentiation and speciation. Here we investigate the regulatory and evolutionary impact of DNA methylation in five tissues of two Ficedula flycatcher species and their naturally occurring F1 hybrids.

RESULTS

We show that the density of CpG in the promoters of genes determines the strength of the association between DNA methylation and gene expression. The impact of DNA methylation on gene expression varies among tissues with the brain showing unique patterns. Differentially expressed genes between parental species are predicted by genetic and methylation differentiation in CpG-rich promoters. However, both these factors fail to predict hybrid misexpression suggesting that promoter mismethylation is not a main determinant of hybrid misexpression in Ficedula flycatchers. Using allele-specific methylation estimates in hybrids, we also determine the genome-wide contribution of cis- and trans effects in DNA methylation differentiation. These distinct mechanisms are roughly balanced in all tissues except the brain, where trans differences predominate.

CONCLUSIONS

Overall, this study provides insight on the regulatory and evolutionary impact of DNA methylation in songbirds.

Transgressive gene expression and expression plasticity under thermal stress in a stable hybrid zone

Interspecific hybridization can lead to myriad outcomes, including transgressive phenotypes in which the hybrids are more fit than either parent species. Such hybrids may display important traits in the context of climate change, able to respond to novel environmental conditions not previously experienced by the parent populations. While this has been evaluated in an agricultural context, the role of transgressive hybrids under changing conditions in the wild remains largely unexplored; this is especially true regarding transgressive gene expression. Using the blue mussel species complex (genus Mytilus) as a model system, we investigated the effects of hybridization on temperature induced gene expression plasticity by comparing expression profiles in parental species and their hybrids following a 2-week thermal challenge. Hybrid expression plasticity was most often like one parent or the other (50%). However, a large fraction of genes (26%) showed transgressive expression plasticity (i.e. the change in gene expression was either greater or lesser than that of both parent species), while only 2% were intermediately plastic in hybrids. Despite their close phylogenetic relationship, there was limited overlap in the differentially expressed genes responding to temperature, indicating interspecific differences in the responses to high temperature in which responses from hybrids are distinct from both parent species. We also identified differentially expressed long non-coding RNAs (lncRNAs), which we suggest may contribute to species-specific differences in thermal tolerance. Our findings provide important insight into the impact of hybridization on gene expression under warming. We propose transgressive hybrids may play an important role in population persistence under future warming conditions.

Hybridization and gene expression: Beyond differentially expressed genes

Gene expression has a key role in reproductive isolation, and studies of hybrid gene expression have identified mechanisms causing hybrid sterility. Here, we review the evidence for altered gene expression following hybridization and outline the mechanisms shown to contribute to altered gene expression in hybrids. Transgressive gene expression, transcending that of both parental species, is pervasive in early generation sterile hybrids, but also frequently observed in viable, fertile hybrids. We highlight studies showing that hybridization can result in transgressive gene expression, also in established hybrid lineages or species. Such extreme patterns of gene expression in stabilized hybrid taxa suggest that altered hybrid gene expression may result in hybridization-derived evolutionary novelty. We also conclude that while patterns of misexpression in hybrids are well documented, the understanding of the mechanisms causing misexpression is lagging. We argue that jointly assessing differences in cell composition and cell-specific changes in gene expression in hybrids, in addition to assessing changes in chromatin and methylation, will significantly advance our understanding of the basis of altered gene expression. Moreover, uncovering to what extent evolution of gene expression results in altered expression for individual genes, or entire networks of genes, will advance our understanding of how selection moulds gene expression. Finally, we argue that jointly studying the dual roles of altered hybrid gene expression, serving both as a mechanism for reproductive isolation and as a substrate for hybrid ecological adaptation, will lead to significant advances in our understanding of the evolution of gene expression.

Separating phases of allopolyploid evolution with resynthesized and natural Capsella bursa-pastoris

Allopolyploidization is a frequent evolutionary transition in plants that combines whole-genome duplication (WGD) and interspecific hybridization. The genome of an allopolyploid species results from initial interactions between parental genomes and long-term evolution. Distinguishing the contributions of these two phases is essential to understanding the evolutionary trajectory of allopolyploid species. Here, we compared phenotypic and transcriptomic changes in natural and resynthesized Capsella allotetraploids with their diploid parental species. We focused on phenotypic traits associated with the selfing syndrome and on transcription-level phenomena such as expression-level dominance (ELD), transgressive expression (TRE), and homoeolog expression bias (HEB). We found that selfing syndrome, high pollen, and seed quality in natural allotetraploids likely resulted from long-term evolution. Similarly, TRE and most down-regulated ELD were only found in natural allopolyploids. Natural allotetraploids also had more ELD toward the self-fertilizing parental species than resynthesized allotetraploids, mirroring the establishment of the selfing syndrome. However, short-term changes mattered, and 40% of the cases of ELD in natural allotetraploids were already observed in resynthesized allotetraploids. Resynthesized allotetraploids showed striking variation of HEB among chromosomes and individuals. Homoeologous synapsis was its primary source and may still be a source of genetic variation in natural allotetraploids. In conclusion, both short- and long-term mechanisms contributed to transcriptomic and phenotypic changes in natural allotetraploids. However, the initial gene expression changes were largely reshaped during long-term evolution leading to further morphological changes.

Gene regulation and speciation in a migratory divide between songbirds

Behavioral variation abounds in nature. This variation is important for adaptation and speciation, but its molecular basis remains elusive. Here, we use a hybrid zone between two subspecies of songbirds that differ in migration - an ecologically important and taxonomically widespread behavior---to gain insight into this topic. We measure gene expression in five brain regions. Differential expression between migratory states was dominated by circadian genes in all brain regions. The remaining patterns were largely brain-region specific. For example, expression differences between the subspecies that interact with migratory state likely help maintain reproductive isolation in this system and were documented in only three brain regions. Contrary to existing work on regulatory mechanisms underlying species-specific traits, two lines of evidence suggest that trans- (vs. cis) regulatory changes underlie these patterns - no evidence for allele-specific expression in hybrids and minimal associations between genomic differentiation and expression differences. Additional work with hybrids shows expression levels were often distinct (transgressive) from parental forms. Behavioral contrasts and functional enrichment analyses allowed us to connect these patterns to mitonuclear incompatibilities and compensatory responses to stress that could exacerbate selection on hybrids and contribute to speciation.

MarFERReT, an open-source, version-controlled reference library of marine microbial eukaryote functional genes

Metatranscriptomics generates large volumes of sequence data about transcribed genes in natural environments. Taxonomic annotation of these datasets depends on availability of curated reference sequences. For marine microbial eukaryotes, current reference libraries are limited by gaps in sequenced organism diversity and barriers to updating libraries with new sequence data, resulting in taxonomic annotation of about half of eukaryotic environmental transcripts. Here, we introduce Marine Functional EukaRyotic Reference Taxa (MarFERReT), a marine microbial eukaryotic sequence library designed for use with taxonomic annotation of eukaryotic metatranscriptomes. We gathered 902 publicly accessible marine eukaryote genomes and transcriptomes and assessed their sequence quality and cross-contamination issues, selecting 800 validated entries for inclusion in MarFERReT. Version 1.1 of MarFERReT contains reference sequences from 800 marine eukaryotic genomes and transcriptomes, covering 453 species- and strain-level taxa, totaling nearly 28 million protein sequences with associated NCBI and PR2 Taxonomy identifiers and Pfam functional annotations. The MarFERReT project repository hosts containerized build scripts, documentation on installation and use case examples, and information on new versions of MarFERReT.

Different complex regulatory phenotypes underlie hybrid male sterility in divergent rodent crosses

Hybrid incompatibilities are a critical component of species barriers and may arise due to negative interactions between divergent regulatory elements in parental species. We used a comparative approach to identify common themes in the regulatory phenotypes associated with hybrid male sterility in two divergent rodent crosses, dwarf hamsters and house mice. We investigated three potential characteristic regulatory phenotypes in hybrids including the propensity towards over or underexpression relative to parental species, the influence of developmental stage on the extent of misexpression, and the role of the sex chromosomes on misexpression phenotypes. In contrast to near pervasive overexpression in hybrid house mice, we found that misexpression in hybrid dwarf hamsters was dependent on developmental stage. In both house mouse and dwarf hamster hybrids, however, misexpression increased with the progression of spermatogenesis, although to varying extents and with potentially different consequences. In both systems, we detected sex-chromosome specific overexpression in stages of spermatogenesis where inactivated X chromosome expression was expected, but the hybrid overexpression phenotypes were fundamentally different. Importantly, misexpression phenotypes support the presence of multiple histological blocks to spermatogenesis in dwarf hamster hybrids, including a potential role of meiotic stalling early in spermatogenesis. Collectively, we demonstrate that while there are some similarities in hybrid regulatory phenotypes of house mice and dwarf hamsters, there are also clear differences that point towards unique mechanisms underlying hybrid male sterility in each system. Our results highlight the potential of comparative approaches in helping to understand the importance of disrupted gene regulation in speciation.

Learning and memory in hybrid migratory songbirds: cognition as a reproductive isolating barrier across seasons

Hybrid zones can be used to identify traits that maintain reproductive isolation and contribute to speciation. Cognitive traits may serve as post-mating reproductive isolating barriers, reducing the fitness of hybrids if, for example, misexpression occurs in hybrids and disrupts important neurological mechanisms. We tested this hypothesis in a hybrid zone between two subspecies of Swainson's thrushes (Catharus ustulatus) using two cognitive tests-an associative learning spatial test and neophobia test. We included comparisons across the sexes and seasons (spring migration and winter), testing if hybrid females performed worse than males (as per Haldane's rule) and if birds (regardless of ancestry or sex) performed better during migration, when they are building navigational maps and encountering new environments. We documented reduced cognitive abilities in hybrids, but this result was limited to males and winter. Hybrid females did not perform worse than males in either season. Although season was a significant predictor of performance, contrary to our prediction, all birds learned faster during the winter. The hypothesis that cognitive traits could serve as post-mating isolating barriers is relatively new; this is one of the first tests in a natural hybrid zone and non-food-caching species. We also provide one of the first comparisons of cognitive abilities between seasons. Future neurostructural and neurophysiological work should be used to examine mechanisms underlying our behavioral observations.

Transcriptomic insight into the hybridization mechanism of the Tambacu, a hybrid from Colossoma macropomum (Tambaqui) and Piaractus mesopotamicus (Pacu)

Interspecific hybrids are highly complex organisms, especially considering aspects related to the organization of genetic material. The diversity of possibilities created by the genetic combination between different species makes it difficult to establish a large-scale analysis methodology. An example of this complexity is Tambacu, an interspecific hybrid of Colossoma macropomum (Tambaqui) and Piaractus mesopotamicus (Pacu). Either genotype represents an essential role in South American aquaculture. However, despite this importance, the genetic information for these genotypes is still highly scarce in specialized databases. Using RNA-Seq analysis, we characterized the transcriptome of white muscle from Pacu, Tambaqui, and their interspecific hybrid (Tambacu). The sequencing process allowed us to obtain a significant number of reads (approximately 53 billion short reads). A total of annotated contigs were 37,285, 96,738, and 158,709 for Pacu, Tambaqui, and Tambacu. After that, we performed a comparative analysis of the transcriptome of the three genotypes, where we evaluated the differential expression (Tambacu vs Pacu = 11,156, and Tambacu vs Tambaqui = 876) profile of the transcript and the degree of similarity between the nucleotide sequences between the genotypes. We assessed the intensity and pattern of expression across genotypes using differential expression information. Clusterization analysis showed a closer relationship between Tambaqui and Tambacu. Furthermore, digital differential expression analysis selected some target genes related to essential cellular processes to evaluate and validate the expression through the RT-qPCR. The RT-qPCR analysis demonstrated significantly (p < 0.05) elevated expression of the mafbx, foxo1a, and rgcc genes in the hybrid compared to the parents. Likewise, we can observe genes significantly more expressed in Pacu (mtco1 and mylpfa) and mtco2 in Tambaqui. Our results showed that the phenotype presented by Tambacu might be associated with changes in the gene expression profile and not necessarily with an increase in gene variability. Thus, the molecular mechanisms underlying these "hybrid effects" may be related to additive and, in some cases, dominant regulatory interactions between parental alleles that act directly on gene regulation in the hybrid transcripts.

Higher temperature induces oxidative stress in hybrids but not in parental species: A case study of crested newts

Ectotherms are particularly sensitive to global warming due to their limited capacity to thermoregulate, which can impact their performance and fitness. From a physiological standpoint, higher temperatures often enhance biological processes that can induce the production of reactive oxygen species and result in a state of cellular oxidative stress. Temperature alters interspecific interactions, including species hybridization. Hybridization under different thermal conditions could amplify parental (genetic) incompatibilities, thus affecting a hybrid's development and distribution. Understanding the impact of global warming on the physiology of hybrids and particularly their oxidative status could help in predicting future scenarios in ecosystems and in hybrids. In the present study, we investigated the effect of water temperature on the development, growth and oxidative stress of two crested newt species and their reciprocal hybrids. Larvae of Triturus macedonicus and T. ivanbureschi, and their T. macedonicus-mothered and T. ivanbureschi-mothered hybrids were exposed for 30 days to temperatures of 19°C and 24°C. Under the higher temperature, the hybrids experienced increases in both growth and developmental rates, while parental species exhibited accelerated growth (T. macedonicus) or development (T. ivanbureschi). Warm conditions also had different effects on the oxidative status of hybrid and parental species. Parental species had enhanced antioxidant responses (catalase, glutathione peroxidase, glutathione S-transferase and SH groups), which allowed them to alleviate temperature-induced stress (revealed by the absence of oxidative damage). However, warming induced an antioxidant response in the hybrids, including oxidative damage in the form of lipid peroxidation. These findings point to a greater disruption of redox regulation and metabolic machinery in hybrid newts, which can be interpreted as the cost of hybridization that is likely linked to parental incompatibilities expressed under a higher temperature. Our study aims to improve mechanistic understanding of the resilience and distribution of hybrid species that cope with climate-driven changes.

Monogeneans in intergeneric hybrids of leuciscid fish: Is parasite infection driven by hybrid heterosis, genetic incompatibilities, or host-parasite coevolutionary interactions?

BACKGROUND

Several hypotheses have been proposed to explain parasite infection in parental species and their hybrids. Hybrid heterosis is generally applied to explain the advantage for F1 generations of hybrids exhibiting a lower level of parasite infection when compared to parental species. Post-F1 generations often suffer from genetic incompatibilities potentially reflected in the higher level of parasite infection when compared to parental species. However, the presence of specific parasites in an associated host is also limited by close coevolutionary genetic host-parasite associations. This study focused on monogenean parasites closely associated with two leuciscid fish species-common bream and roach-with the aim of comparing the level of monogenean infection between parental species and hybrids representing two F1 generations with different mtDNA and two backcross generations with different cyto-nuclear compositions.

RESULTS

Monogenean infection in F1 generations of hybrids was lower when compared to parental species, in line with the hybrid heterosis hypothesis. Monogenean infection in backcross generations exhibited similarities with the parental species whose genes contributed more to the backcross genotype. The distribution of monogeneans associated with one or the other parental species showed the same asymmetry with a higher proportion of roach-associated monogeneans in both F1 generations and backcross generation with roach in the paternal position. A higher proportion of common bream-associated monogeneans was found in backcross generation with common bream in the paternal position.

CONCLUSIONS

Our study indicated that cyto-nuclear incompatibilities in hybrids do not induce higher monogenean infection in backcross generations when compared to parental species. However, as backcross hybrids with a higher proportion of the genes of one parental taxon also exhibited high level of this parental taxon-associated parasites, host-parasite coevolutionary interactions seem to play an obvious role in determining the level of infection of host-specific monogeneans in hybrids.

Patterns of Performance Variation Between Animal Hybrids and their Parents: A Meta-analysis

Hybridization is a widespread phenomenon in animals, and hybrid heterosis/breakdown could be key processes determining the evolutionary dynamics of hybrids. Indeed, hybrids are not consistently disadvantaged compared to the parental lineages, as was historically assumed. Multiple processes could lead to performance differences between parental lineages and their hybrids. Despite many studies evaluated the performance of hybrids, a quantitative synthesis is required to assess the general pattern. Here we used meta-analytic and meta-regression approaches to quantify the fitness differences between parental lineages and their hybrids, and to identify possible processes that could lead to these differences. Specifically, we tested biological and methodological parameters that could determine differences in performance between hybrids and parental lineages. Hybrid performance was extremely variable across studies, being often significantly higher or lower compared to the mean performance of their parents. Nevertheless, the averaged hybrid performance was similar to the fitness of parental lineages, with differences across studies related to how performance was assessed. Genetic divergence between parental lineages, and the approach used to identify hybrids were the parameters most strongly related to variation in hybrid performance. Performance was lower for hybrids between distantly related lineages. Furthermore, study settings and the use of imprecise approaches for hybrid identification (e.g. morphology-based) can bias assessments of performance. Studies performed on wild populations and using genetic approaches for hybrid identification detected more often a decreased hybrid performance, compared to laboratory studies. We highlight the importance of appropriate settings for a realistic understanding of the evolutionary impacts of hybridization.

Genomic Signatures of Mitonuclear Coevolution in Mammals

Mitochondrial (mt) and nuclear-encoded proteins are integrated in aerobic respiration, requiring co-functionality among gene products from fundamentally different genomes. Different evolutionary rates, inheritance mechanisms, and selection pressures set the stage for incompatibilities between interacting products of the two genomes. The mitonuclear coevolution hypothesis posits that incompatibilities may be avoided if evolution in one genome selects for complementary changes in interacting genes encoded by the other genome. Nuclear compensation, in which deleterious mtDNA changes are offset by compensatory nuclear changes, is often invoked as the primary mechanism for mitonuclear coevolution. Yet, direct evidence supporting nuclear compensation is rare. Here, we used data from 58 mammalian species representing eight orders to show strong correlations between evolutionary rates of mt and nuclear-encoded mt-targeted (N-mt) proteins, but not between mt and non-mt-targeted nuclear proteins, providing strong support for mitonuclear coevolution across mammals. N-mt genes with direct mt interactions also showed the strongest correlations. Although most N-mt genes had elevated dN/dS ratios compared to mt genes (as predicted under nuclear compensation), N-mt sites in close contact with mt proteins were not overrepresented for signs of positive selection compared to noncontact N-mt sites (contrary to predictions of nuclear compensation). Furthermore, temporal patterns of N-mt and mt amino acid substitutions did not support predictions of nuclear compensation, even in positively selected, functionally important residues with direct mitonuclear contacts. Overall, our results strongly support mitonuclear coevolution across ∼170 million years of mammalian evolution but fail to support nuclear compensation as the major mode of mitonuclear coevolution.

Physiological aspects of sex differences and Haldane's rule in Rumex hastatulus

Haldane's rule (HR, impairment of fertility and/or viability of interracial hybrids) seems to be one of few generalizations in evolutionary biology. The validity of HR has been confirmed in animals, and more recently in some dioecious plants (Silene and Rumex). Dioecious Rumex hastatulus has two races differing in the sex chromosome system: Texas (T) and North Carolina (NC), and T × NC males showed both reduced pollen fertility and rarity-two classical symptoms of Haldane's rule (HR). The reduced fertility of these plants has a simple mechanistic explanation, but the reason for their rarity was not elucidated. Here, we measured selected physiological parameters related to the antioxidant defense system in parental races and reciprocal hybrids of R. hastatulus. We showed that the X-autosome configurations, as well as asymmetries associated with Y chromosomes and cytoplasm, could modulate this system in hybrids. The levels and quantitative patterns of the measured parameters distinguish the T × NC hybrid from the other analyzed forms. Our observations suggest that the rarity of T × NC males is caused postzygotically and most likely related to the higher level of oxidative stress induced by the chromosomal incompatibilities. It is the first report on the physiological aspects of HR in plants.

Extensive transgressive gene expression in testis but not ovary in the homoploid hybrid Italian sparrow

Hybridization can result in novel allelic combinations which can impact the hybrid phenotype through changes in gene expression. While misexpression in F₁ hybrids is well documented, how gene expression evolves in stabilized hybrid taxa remains an open question. As gene expression evolves in a stabilizing manner, break‐up of co‐evolved cis‐ and trans‐regulatory elements could lead to transgressive patterns of gene expression in hybrids. Here, we address to what extent gonad gene expression has evolved in an established and stable homoploid hybrid, the Italian sparrow (Passer italiae). Through comparison of gene expression in gonads from individuals of the two parental species (i.e., house and Spanish sparrow) to that of Italian sparrows, we find evidence for strongly transgressive expression in male Italian sparrows—2530 genes (22% of testis genes tested for inheritance) exhibit expression patterns outside the range of both parent species. In contrast, Italian sparrow ovary expression was similar to that of one of the parent species, the house sparrow (Passer domesticus). Moreover, the Italian sparrow testis transcriptome is 26 times as diverged from those of the parent species as the parental transcriptomes are from each other, despite being genetically intermediate. This highlights the potential for regulation of gene expression to produce novel variation following hybridization. Genes involved in mitochondrial respiratory chain complexes and protein synthesis are enriched in the subset that is over‐dominantly expressed in Italian sparrow testis, suggesting that selection on key functions has moulded the hybrid Italian sparrow transcriptome.

Frequent Paternal Mitochondrial Inheritance and Rapid Haplotype Frequency Shifts in Copepod Hybrids

Mitochondria are assumed to be maternally inherited in most animal species, and this foundational concept has fostered advances in phylogenetics, conservation, and population genetics. Like other animals, mitochondria were thought to be solely maternally inherited in the marine copepod Tigriopus californicus, which has served as a useful model for studying mitonuclear interactions, hybrid breakdown, and environmental tolerance. However, we present PCR, Sanger sequencing, and Illumina Nextera sequencing evidence that extensive paternal mitochondrial DNA (mtDNA) transmission is occurring in inter-population hybrids of T. californicus. PCR on four types of crosses between three populations (total sample size of 376 F1 individuals) with 20% genome-wide mitochondrial divergence showed 2% to 59% of F1 hybrids with both paternal and maternal mtDNA, where low and high paternal leakage values were found in different cross directions of the same population pairs. Sequencing methods further verified nucleotide similarities between F1 mtDNA and paternal mtDNA sequences. Interestingly, the paternal mtDNA in F1s from some crosses inherited haplotypes that were uncommon in the paternal population. Compared to some previous research on paternal leakage, we employed more rigorous methods to rule out contamination and false detection of paternal mtDNA due to non-functional nuclear mitochondrial DNA fragments. Our results raise the potential that other animal systems thought to only inherit maternal mitochondria may also have paternal leakage, which would then affect the interpretation of past and future population genetics or phylogenetic studies that rely on mitochondria as uniparental markers.

Seminal fluid gene expression and reproductive fitness in Drosophila melanogaster

Background

The rapid evolution of seminal fluid proteins (SFPs) has been suggested to be driven by adaptations to postcopulatory sexual selection (e.g. sperm competition). However, we have recently shown that most SFPs evolve rapidly under relaxed selective pressures. Given the role of SFPs in competition for fertilization phenotypes, like the ability to transfer and store sperm and the modulation of female receptivity and ovulation, the prevalence of selectively relaxed SFPs appears as a conundrum. One possible explanation is that selection on SFPs might be relaxed in terms of protein amino acid content, but adjustments of expression are essential for post-mating function. Interestingly, there is a general lack of systematic implementation of gene expression perturbation assays to monitor their effect on phenotypes related to sperm competition.

Results

We successfully manipulated the expression of 16 SFP encoding genes using tissue-specific knockdowns (KDs) and determined the effect of these genes’ perturbation on three important post-mating phenotypes: female refractoriness to remating, defensive (P1), and offensive (P2) sperm competitive abilities in Drosophila melanogaster. Our analyses show that KDs of tested SFP genes do not affect female refractoriness to remating and P2, however, most gene KDs significantly decreased P1. Moreover, KDs of SFP genes that are selectively constrained in terms of protein-coding sequence evolution have lower P1 than KDs of genes evolving under relaxed selection.

Conclusions

Our results suggest a more predominant role, than previously acknowledged, of variation in gene expression than coding sequence changes on sperm competitive ability in D. melanogaster.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-01975-1.

Mitonuclear mismatch alters nuclear gene expression in naturally introgressed Rhinolophus bats

Background

Mitochondrial function involves the interplay between mitochondrial and nuclear genomes. Such mitonuclear interactions can be disrupted by the introgression of mitochondrial DNA between taxa or divergent populations. Previous studies of several model systems (e.g. Drosophila) indicate that the disruption of mitonuclear interactions, termed mitonuclear mismatch, can alter nuclear gene expression, yet few studies have focused on natural populations.

Results

Here we study a naturally introgressed population in the secondary contact zone of two subspecies of the intermediate horseshoe bat (Rhinolophus affinis), in which individuals possess either mitonuclear matched or mismatched genotypes. We generated transcriptome data for six tissue types from five mitonuclear matched and five mismatched individuals. Our results revealed strong tissue-specific effects of mitonuclear mismatch on nuclear gene expression with the largest effect seen in pectoral muscle. Moreover, consistent with the hypothesis that genes associated with the response to oxidative stress may be upregulated in mitonuclear mismatched individuals, we identified several such gene candidates, including DNASE1L3, GPx3 and HSPB6 in muscle, and ISG15 and IFI6 in heart.

Conclusion

Our study reveals how mitonuclear mismatch arising from introgression in natural populations is likely to have fitness consequences. Underlying the processes that maintain mitonuclear discordance is a step forward to understand the role of mitonuclear interactions in population divergence and speciation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12983-021-00424-x.

Differences in Thermal Tolerance between Parental Species Could Fuel Thermal Adaptation in Hybrid Wood Ants

AbstractGenetic variability is essential for adaptation and could be acquired via hybridization with a closely related lineage. We use ants to investigate thermal adaptation and the link between temperature and genetic variation arising from hybridization. We test for differences in cold and heat tolerance between Finnish Formica polyctena and Formica aquilonia wood ants and their naturally occurring hybrids. Using workers, we find that the parental individuals differ in both cold and heat tolerances and express thermal limits that reflect their global distributions. Hybrids, however, cannot combine thermal tolerance of parental species as they have the same heat tolerance as F. polyctena but not the same cold tolerance as F. aquilonia. We then focus on a single hybrid population to investigate the relationship between temperature variation and genetic variation across 16 years using reproductive individuals. On the basis of the thermal tolerance results, we expected the frequency of putative F. polyctena alleles to increase in warm years and F. aquilonia alleles to increase in cold years. We find support for this in hybrid males but not in hybrid females. These results contribute to understanding the outcomes of hybridization, which may be sex specific or depend on the environment. Furthermore, genetic variability resulting from hybridization could help hybrid wood ants cope with changing thermal conditions.

Recovery from hybrid breakdown reveals a complex genetic architecture of mitonuclear incompatibilities

Reproductive isolation is often achieved when genes that are neutral or beneficial in their genomic background become functionally incompatible in a foreign genomic background, causing inviability, sterility or other forms of low fitness in hybrids. Recent studies suggest that mitonuclear interactions are among the initial incompatibilities to evolve at early stages of population divergence across taxa. Yet, the genomic architecture of mitonuclear incompatibilities has rarely been elucidated. We employ an experimental evolution approach starting with low-fitness F₂ interpopulation hybrids of the copepod Tigriopus californicus, in which frequencies of compatible and incompatible nuclear alleles change in response to an alternative mitochondrial background. After about nine generations, we observe a generalized increase in population size and in survivorship, suggesting efficiency of selection against maladaptive phenotypes. Whole genome sequencing of evolved populations showed some consistent allele frequency changes across three replicates of each reciprocal cross, but markedly different patterns between mitochondrial backgrounds. In only a few regions (~6.5% of the genome), the same parental allele was overrepresented irrespective of the mitochondrial background. About 33% of the genome showed allele frequency changes consistent with divergent selection, with the location of these genomic regions strongly differing between mitochondrial backgrounds. In 87% and 89% of these genomic regions, the dominant nuclear allele matched the associated mitochondrial background, consistent with mitonuclear co-adaptation. These results suggest that mitonuclear incompatibilities have a complex polygenic architecture that differs between populations, potentially generating genome-wide barriers to gene flow between closely related taxa.

Incomplete reproductive isolation and strong transcriptomic response to hybridization between sympatric sister species of salmon

Global change is altering ecosystems at an unprecedented rate. The resulting shifts in species ranges and reproductive timing are opening the potential for hybridization between closely related species which could dramatically alter the genetic diversity, adaptive capacity and evolutionary trajectory of interbreeding taxa. Here, we used behavioural breeding experiments, in vitro fertilization experiments, and whole-transcriptome gene expression data to assess the potential for and consequences of hybridization between Chinook and Coho salmon. We show that behavioural and gametic prezygotic barriers between socio-economically valuable Chinook and Coho salmon are incomplete. Postzygotically, we demonstrate a clear transcriptomic response to hybridization among F1 Chinook-Coho offspring. Genes transgressively expressed within hybrids were significantly enriched with genes encoded in the nucleus but localized to the mitochondrion, suggesting a potential role for mito-nuclear incompatibilities as a postzygotic mechanism of hybrid breakdown. Chinook and Coho salmon are expected to continue to respond to climate change with shifts in migration timing and habitat use, potentiating hybridization between these species. The downstream consequences of hybridization on the future of these threatened salmon, and the ecosystems they inhabit, is unknown.

Genetics of autoimmunity in plants: an evolutionary genetics perspective

Autoimmunity in plants has been found in numerous hybrids as a form of hybrid necrosis and mutant panels. Uncontrolled cell death is a main cellular outcome of autoimmunity, which negatively impacts growth. Its occurrence highlights the vulnerable nature of the plant immune system. Genetic investigation of autoimmunity in hybrid plants revealed that extreme variation in the immune receptor repertoire is a major contributor, reflecting an evolutionary conundrum that plants face in nature. In this review, we discuss natural variation in the plant immune system and its contribution to fitness. The value of autoimmunity genetics lies in its ability to identify combinations of a natural immune receptor and its partner that are predisposed to triggering autoimmunity. The network of immune components for autoimmunity becomes instrumental in revealing mechanistic details of how immune receptors recognize cellular invasion and activate signaling. The list of autoimmunity-risk variants also allows us to infer evolutionary processes contributing to their maintenance in the natural population. Our approach to autoimmunity, which integrates mechanistic understanding and evolutionary genetics, has the potential to serve as a prognosis tool to optimize immunity in crops.

Hybrid Incompatibilities and Transgressive Gene Expression Between Two Closely Related Subspecies of Drosophila

Genome-wide assays of expression between species and their hybrids have identified genes that become either over- or underexpressed relative to the parental species (i.e., transgressive). Transgressive expression in hybrids is of interest because it highlights possible changes in gene regulation linked to hybrid dysfunction. Previous studies in Drosophila that used long-diverged species pairs with complete or nearly complete isolation (i.e., full sterility and partial inviability of hybrids) and high-levels of genome misregulation have found correlations between expression and coding sequence divergence. The work highlighted the possible effects of directional selection driving sequence divergence and transgressive expression. Whether the same is true for taxa at early stages of divergence that have only achieved partial isolation remains untested. Here, we reanalyze previously published genome expression data and available genome sequence reads from a pair of partially isolated subspecies of Drosophila to compare expression and sequence divergence. We find a significant correlation in rates of expression and sequence evolution, but no support for directional selection driving transgressive expression in hybrids. We find that most transgressive genes in hybrids show no differential expression between parental subspecies and used SNP data to explore the role of stabilizing selection through compensatory mutations. We also examine possible misregulation through cascade effects that could be driven by interacting gene networks or co-option of off-target cis-regulatory elements.

Allele-Specific Expression and Evolution of Gene Regulation Underlying Acute Heat Stress Response and Local Adaptation in the Copepod Tigriopus californicus

Geographic variation in environmental temperature can select for local adaptation among conspecific populations. Divergence in gene expression across the transcriptome is a key mechanism for evolution of local thermal adaptation in many systems, yet the genetic mechanisms underlying this regulatory evolution remain poorly understood. Here we examine gene expression in 2 locally adapted Tigriopus californicus populations (heat tolerant San Diego, SD, and less tolerant Santa Cruz, SC) and their F1 hybrids during acute heat stress response. Allele-specific expression (ASE) in F1 hybrids was used to determine cis-regulatory divergence. We found that the number of genes showing significant allelic imbalance increased under heat stress compared to unstressed controls. This suggests that there is significant population divergence in cis-regulatory elements underlying heat stress response. Specifically, the number of genes showing an excess of transcripts from the more thermal tolerant (SD) population increased with heat stress while that number of genes with an SC excess was similar in both treatments. Inheritance patterns of gene expression also revealed that genes displaying SD-dominant expression phenotypes increase in number in response to heat stress; that is, across loci, gene expression in F1's following heat stress showed more similarity to SD than SC, a pattern that was absent in the control treatment. The observed patterns of ASE and inheritance of gene expression provide insight into the complex processes underlying local adaptation and thermal stress response.

Haldane's rule in the placenta: Sex-biased misregulation of the Kcnq1 imprinting cluster in hybrid mice

Hybrid phenotypes that contribute to postzygotic reproductive isolation often exhibit pronounced asymmetry, both between reciprocal crosses and between the sexes in accordance with Haldane's rule. Inviability in mammalian hybrids is associated with parent-of-origin placental growth abnormalities for which misregulation of imprinted gene (IGs) is the leading candidate mechanism. However, direct evidence for the involvement of IGs in hybrid growth dysplasia is limited. We used transcriptome and reduced representation bisulfite sequencing to conduct the first genome-scale assessment of the contribution of IGs to parent-of-origin placental growth dysplasia in the cross between the house mouse (Mus musculus domesticus) and the Algerian mouse (Mus spretus). IGs with transgressive expression and methylation were concentrated in the Kcnq1 cluster, which contains causal genes for prenatal growth abnormalities in mice and humans. Hypermethylation of the cluster's imprinting control region, and consequent misexpression of the genes Phlda2 and Ascl2, is a strong candidate mechanism for transgressive placental undergrowth. Transgressive placental and gene regulatory phenotypes, including expression and methylation in the Kcnq1 cluster, were more extreme in hybrid males. Although consistent with Haldane's rule, male-biased defects are unexpected in rodent placenta because the X-chromosome is effectively hemizygous in both sexes. In search of an explanation, we found evidence of leaky imprinted (paternal) X-chromosome inactivation in hybrid female placenta, an epigenetic disturbance that may buffer females from the effects of X-linked incompatibilities to which males are fully exposed. Sex differences in chromatin structure on the X and sex-biased maternal effects are nonmutually exclusive alternative explanations for adherence to Haldane's rule in hybrid placenta. The results of this study contribute to understanding the genetic basis of hybrid inviability in mammals, and the role of IGs in speciation.

Surprising spatiotemporal stability of a multi-peak fitness landscape revealed by independent field experiments measuring hybrid fitness

The effect of the environment on fitness in natural populations is a fundamental question in evolutionary biology. However, experimental manipulations of both environment and phenotype at the same time are rare. Thus, the relative importance of the competitive environment versus intrinsic organismal performance in shaping the location, height, and fluidity of fitness peaks and valleys remains largely unknown. Here, we experimentally tested the effect of competitor frequency on the complex fitness landscape driving adaptive radiation of a generalist and two trophic specialist pupfishes, a scale-eater and molluscivore, endemic to hypersaline lakes on San Salvador Island (SSI), Bahamas. We manipulated phenotypes, by generating 3407 F4/F5 lab-reared hybrids, and competitive environment, by altering the frequency of rare transgressive hybrids between field enclosures in two independent lake populations. We then tracked hybrid survival and growth rates across these four field enclosures for 3-11 months. In contrast to competitive speciation theory, we found no evidence that the frequency of hybrid phenotypes affected their survival. Instead, we observed a strikingly similar fitness landscape to a previous independent field experiment, each supporting multiple fitness peaks for generalist and molluscivore phenotypes and a large fitness valley isolating the divergent scale-eater phenotype. These features of the fitness landscape were stable across manipulated competitive environments, multivariate trait axes, and spatiotemporal heterogeneity. We suggest that absolute performance constraints and divergent gene regulatory networks shape macroevolutionary (interspecific) fitness landscapes in addition to microevolutionary (intraspecific) competitive dynamics. This interplay between organism and environment underlies static and dynamic features of the adaptive landscape.

Speciation and changes in male gene expression in Drosophila

It has long been acknowledged that changes in the regulation of gene expression may account for major organismal differences. However, we still do not fully understand how changes in gene expression evolve and how do such changes influence organisms' differences. We are even less aware of the impact such changes might have in restricting gene flow between species. Here, we focus on studies of gene expression and speciation in the Drosophila model. We review studies that have identified gene interactions in post-mating reproductive isolation and speciation, particularly those that modulate male gene expression. We also address studies that have experimentally manipulated changes in gene expression to test their effect in post-mating reproductive isolation. We highlight the need for a more in-depth analysis of the role of selection causing disrupted gene expression of such candidate genes in sterile/inviable hybrids. Moreover, we discuss the relevance to incorporate more routinely assays that simultaneously evaluate the potential effects of environmental factors and genetic background in modulating plastic responses in male genes and their potential role in speciation.

Ecological divergence in sympatry causes gene misexpression in hybrids

Ecological speciation occurs when reproductive isolation evolves as a byproduct of adaptive divergence between populations. Selection favouring gene regulatory divergence between species could result in transgressive levels of gene expression in F1 hybrids that may lower hybrid fitness. We combined 58 resequenced genomes with 124 transcriptomes to identify patterns of hybrid gene misexpression that may be driven by adaptive regulatory divergence within a young radiation of Cyprinodon pupfishes, which consists of a dietary generalist and two trophic specialists-a molluscivore and a scale-eater. We found more differential gene expression between closely related sympatric specialists than between allopatric generalist populations separated by 1,000 km. Intriguingly, 9.6% of genes that were differentially expressed between sympatric species were also misexpressed in F1 hybrids. A subset of these genes were in highly differentiated genomic regions and enriched for functions important for trophic specialization, including head, muscle and brain development. These regions also included genes that showed evidence of hard selective sweeps and were significantly associated with oral jaw length-the most rapidly diversifying skeletal trait in this radiation. Our results indicate that divergent ecological selection in sympatry can contribute to hybrid gene misexpression which may act as a reproductive barrier between nascent species.

Mitonuclear mismatch alters performance and reproductive success in naturally introgressed populations of a montane leaf beetle

Coordination between nuclear and mitochondrial genomes is critical to metabolic processes underlying animals' ability to adapt to local environments, yet consequences of mitonuclear interactions have rarely been investigated in populations where individuals with divergent mitochondrial and nuclear genomes naturally interbreed. Genetic variation in the leaf beetle Chrysomela aeneicollis was assessed along a latitudinal thermal gradient in California's Sierra Nevada. Variation at mitochondrial cytochrome oxidase II (COII) and the nuclear gene phosphoglucose isomerase (PGI) shows concordance and was significantly greater along a 65 km transect than 10 other loci. STRUCTURE analyses using neutral loci identified a southern and northern subpopulation, which interbreed in the central drainage Bishop Creek. COII and PGI were used as indicators of mitochondrial and nuclear genetic variation in field and laboratory experiments conducted on beetles from this admixed population. Fecundity, larval development rate, running speed and male mating frequency were higher for beetles with geographically "matched" than "mismatched" mitonuclear genotypes. Effects of mitonuclear mismatch were largest for individuals with northern nuclear genotypes possessing southern mitochondria and were most pronounced after heat treatment or at high elevation. These findings suggest that mitonuclear incompatibility diminishes performance and reproductive success in nature, effects that could intensify at environmental extremes.

Co-evolution in the Jungle: From Leafcutter Ant Colonies to Chromosomal Ends

Biological entities are multicomponent systems where each part is directly or indirectly dependent on the others. In effect, a change in a single component might have a consequence on the functioning of its partners, thus affecting the fitness of the entire system. In this article, we provide a few examples of such complex biological systems, ranging from ant colonies to a population of amino acids within a single-polypeptide chain. Based on these examples, we discuss one of the central and still challenging questions in biology: how do such multicomponent consortia co-evolve? More specifically, we ask how telomeres, nucleo-protein complexes protecting the integrity of linear DNA chromosomes, originated from the ancestral organisms having circular genomes and thus not dealing with end-replication and end-protection problems. Using the examples of rapidly evolving topologies of mitochondrial genomes in eukaryotic microorganisms, we show what means of co-evolution were employed to accommodate various types of telomere-maintenance mechanisms in mitochondria. We also describe an unprecedented runaway evolution of telomeric repeats in nuclei of ascomycetous yeasts accompanied by co-evolution of telomere-associated proteins. We propose several scenarios derived from research on telomeres and supported by other studies from various fields of biology, while emphasizing that the relevant answers are still not in sight. It is this uncertainty and a lack of a detailed roadmap that makes the journey through the jungle of biological systems still exciting and worth undertaking.

Effects of oxidative stress on sex-specific gene expression in the copepod Tigriopus californicus revealed by single individual RNA-seq

Oxidative stress reflects the imbalance of pro-oxidants and antioxidants. Prolonged oxidative stress can induce cellular damage, diseases and aging, and the effects may be sex-specific. Tigriopus californicus has recently been proposed as an alternative model system for sex-specific studies due to the absence of sex chromosomes. In this study, we used comparative transcriptomic analyses to assess sex-specific transcriptional responses to oxidative stress. Male and female individuals were maintained separately in one of three treatments: 1) control conditions with an algae diet, 2) pro-oxidant (H₂O₂) conditions with an algae diet or 3) decreased antioxidant conditions (reduced carotenoids due to a yeast diet). Single individual RNA-seq was then conducted for twenty-four libraries using Ligation Mediated RNA sequencing (LM-Seq). Variance in gene expression was partitioned into 62.3% between sexes, 26.85% among individuals and 10.85% among treatments. Within each of the three treatments, expression was biased toward females. However, compared to the control treatment, males in both pro-oxidant and decreased antioxidant treatments differentially expressed more genes while females differentially expressed fewer genes but with a greater magnitude of fold change. As the first study of copepods to apply single individual RNA-seq, the findings will contribute to a better understanding of transcriptomic variation among individuals as well as sex-specific response mechanisms to oxidative stress in the absence of sex chromosomes.

Assessing the fitness consequences of mitonuclear interactions in natural populations

Metazoans exist only with a continuous and rich supply of chemical energy from oxidative phosphorylation in mitochondria. The oxidative phosphorylation machinery that mediates energy conservation is encoded by both mitochondrial and nuclear genes, and hence the products of these two genomes must interact closely to achieve coordinated function of core respiratory processes. It follows that selection for efficient respiration will lead to selection for compatible combinations of mitochondrial and nuclear genotypes, and this should facilitate coadaptation between mitochondrial and nuclear genomes (mitonuclear coadaptation). Herein, we outline the modes by which mitochondrial and nuclear genomes may coevolve within natural populations, and we discuss the implications of mitonuclear coadaptation for diverse fields of study in the biological sciences. We identify five themes in the study of mitonuclear interactions that provide a roadmap for both ecological and biomedical studies seeking to measure the contribution of intergenomic coadaptation to the evolution of natural populations. We also explore the wider implications of the fitness consequences of mitonuclear interactions, focusing on central debates within the fields of ecology and biomedicine.

Comparative transcriptomics sheds light on differential adaptation and species diversification between two Melastoma species and their F1 hybrid

Variation in gene expression has been shown to promote adaptive divergence, and can lead to speciation. The plant genus Melastoma, thought to have diversified through adaptive radiation, provides an excellent model for the study of gene expressional changes during adaptive differentiation and following interspecific hybridization. In this study, we performed RNA-seq on M. candidum, M. sanguineum and their F1 hybrid, to investigate the role of gene expression in species diversification within the genus. Reference transcriptomes were assembled using combined data from both parental species, resulting in 50 519 and 48 120 transcripts for the leaf and flower petal, after removing redundancy. Differential expression analysis uncovered 3793 and 2116 differentially expressed (DE) transcripts, most of which are between M. candidum and M. sanguineum. Differential expression was observed for genes related to light responses, as well as genes that regulate the development of leaf trichomes, a trait that among others is thought to protect plants against sunlight, suggesting the differential adaptation of the species to sunlight intensity. The analysis of positively selected genes between the two species also revealed possible differential adaptation to other abiotic stresses such as drought and temperature. In the hybrid, almost all possible modes of expression were observed at the DE transcripts, although at most transcripts, the expression levels were similar to that of either parent instead of being intermediate. A small number of transgressively expressed transcripts that matched genes known to promote plant growth and adaptation to stresses in new environments were also found, possibly explaining the vigour observed in the hybrid. The findings in this study provided insights into the role of gene expression in the diversification of Melastoma, which we believe is an important example for more cross-taxa comparisons in the future.

Individual Culturing of Tigriopus Copepods and Quantitative Analysis of Their Mate-guarding Behavior

Copepods of the genus Tigriopus, which are common zooplankton in rocky tide pools, show precopulatory mate-guarding behavior where a male clasps a potential mate to form a pair. While this phenomenon has attracted interest of researchers, methods for its analysis have not been well described. Here we describe procedures for: 1) individual culturing and staging of Tigriopus juveniles and adults, and 2) video-based analysis of their mate-guarding behavior. The culturing method enables experimental control of paring experience of animals as well as the ability to track their development before behavioral tests. The analysis method allows quantitative evaluation of several aspects of the mate-guarding behavior, including capturing attempts by males and swimming trajectory of mate-guarding pairs. Although these methods were originally established for ethological studies on Tigriopus, with proper modifications they can also be applied to studies of other zooplankton in different research fields, such as physiology, toxicology, and ecological genetics.

Mitochondrial Recombination and Introgression during Speciation by Hybridization

Genome recombination is a major source of genotypic diversity and contributes to adaptation and speciation following interspecies hybridization. The contribution of recombination in these processes has been thought to be largely limited to the nuclear genome because organelles are mostly uniparentally inherited in animals and plants, which prevents recombination. Unicellular eukaryotes such as budding yeasts do, however, transmit mitochondria biparentally, suggesting that during hybridization, both parents could provide alleles that contribute to mitochondrial functions such as respiration and metabolism in hybrid populations or hybrid species. We examined the dynamics of mitochondrial genome transmission and evolution during speciation by hybridization in the natural budding yeast Saccharomyces paradoxus. Using population-scale mitochondrial genome sequencing in two endemic North American incipient species SpB and SpC and their hybrid species SpC*, we found that both parental species contributed to the hybrid mitochondrial genome through recombination. We support our findings by showing that mitochondrial recombination between parental types is frequent in experimental crosses that recreate the early step of this speciation event. In these artificial hybrids, we observed that mitochondrial genome recombination enhances phenotypic variation among diploid hybrids, suggesting that it could play a role in the phenotypic differentiation of hybrid species. Like the nuclear genome, the mitochondrial genome can, therefore, also play a role in hybrid speciation.

The impact of seasonality on niche breadth, distribution range and species richness: a theoretical exploration of Janzen's hypothesis

Being invoked as one of the candidate mechanisms for the latitudinal patterns in biodiversity, Janzen's hypothesis states that the limited seasonal temperature variation in the tropics generates greater temperature stratification across elevations, which makes tropical species adapted to narrower ranges of temperatures and have lower effective dispersal across elevations than species in temperate regions. Numerous empirical studies have documented latitudinal patterns in species elevational ranges and thermal niche breadths that are consistent with the hypothesis, but the theoretical underpinnings remain unclear. This study presents the first mathematical model to examine the evolutionary processes that could back up Janzen's hypothesis and assess the effectiveness of limited seasonal temperature variation to promote speciation along elevation in the tropics. Results suggest that trade-offs in thermal tolerances provide a mechanism for Janzen's hypothesis. Limited seasonal temperature variation promotes gradient speciation not due to the reduction in gene flow that is associated with narrow thermal niche, but due to the pleiotropic effects of more stable divergent selection of thermal tolerance on the evolution of reproductive incompatibility. The proposed modelling approach also provides a potential way to test a speciation model against genetic data.

Widespread hybridization within mound-building wood ants in Southern Finland results in cytonuclear mismatches and potential for sex-specific hybrid breakdown

Hybridization and gene flow between diverging lineages are increasingly recognized as common evolutionary processes, and their consequences can vary from hybrid breakdown to adaptive introgression. We have previously found a population of wood ant hybrids between Formica aquilonia and F. polyctena that shows antagonistic effects of hybridization: females with introgressed alleles show hybrid vigour, whereas males with the same alleles show hybrid breakdown. Here, we investigate whether hybridization is a general phenomenon in this species pair and analyse 647 worker samples from 16 localities in Finland using microsatellite markers and a 1200-bp mitochondrial sequence. Our results show that 27 sampled nests contained parental-like gene pools (six putative F. polyctena and 21 putative F. aquilonia) and all remaining nests (69), from nine localities, contained hybrids of varying degrees. Patterns of genetic variation suggest these hybrids arise from several hybridization events or, instead, have backcrossed to the parental gene pools to varying extents. In contrast to expectations, the mitochondrial haplotypes of the parental species were not randomly distributed among the hybrids. Instead, nests that were closer to parental-like F. aquilonia for nuclear markers preferentially had F. polyctena's mitochondria and vice versa. This systematic pattern suggests there may be underlying selection favouring cytonuclear mismatch and hybridization. We also found a new hybrid locality with strong genetic differences between the sexes similar to those predicted under antagonistic selection on male and female hybrids. Further studies are needed to determine the selective forces that act on male and female genomes in these newly discovered hybrids.

Signals of selection in conditionally expressed genes in the diversification of three horned beetle species

Species radiations may be facilitated by phenotypic differences already present within populations, such as those arising through sex-specific development or developmental processes biased towards particular reproductive or trophic morphs. We sought to test this hypothesis by utilizing a comparative transcriptomic approach to contrast among- and within-species differentiation using three horned beetle species in the genus Onthophagus. These three species exhibit differences along three phenotypic axes reflective of much of the interspecific diversity present within the genus: horn location, polarity of sexual dimorphism and degree of nutritional sensitivity. Our approach combined de novo transcript assembly, assessment of amino acid substitutions (dN/dS) across orthologous gene pairs and integration of gene function and conditional gene expression data. We identified 17 genes across the three species pairs related to axis patterning, development and metabolism with dN/dS > 1 and detected elevated dN/dS in genes related to metabolism and biosynthesis in the most closely related species pair, which is characterized by a loss of nutritional polyphenism and a reversal of sexual dimorphism. Further, we found that genes that are conditionally expressed (i.e. as a function of sex, nutrition or body region) within one of our focal species also showed significantly stronger signals of positive or relaxed purifying selection between species divergent along the same morphological axis (i.e. polarity of sexual dimorphism, degree of nutritional sensitivity or location of horns). Our findings thus reveal a positive relationship between intraspecific differentiation due to condition-specific development and genetic divergences among species.

A Little Bug with a Big Bite: Impact of Hemlock Woolly Adelgid Infestations on Forest Ecosystems in the Eastern USA and Potential Control Strategies

Hemlock woolly adelgid (Adelges tsugae Annand, HWA) remains the single greatest threat to the health and sustainability of hemlock in the eastern USA. The loss of hemlock trees leads to further negative impacts on the diversity and stability of ecosystems in the eastern part of North America. It is, therefore, urgent to develop effective control measures to reduce HWA populations and promote overall hemlock health. Currently available individual and integrated approaches should continue to be evaluated in the laboratory and in the field along with the development of other new and innovative methods.

The Composite Regulatory Basis of the Large X-Effect in Mouse Speciation

The disruption of meiotic sex chromosome inactivation (MSCI) has been proposed to be a major developmental mechanism underlying the rapid evolution of hybrid male sterility. We tested this idea by analyzing cell-specific gene expression across spermatogenesis in two lineages of house mice and their sterile and fertile reciprocal hybrids. We found pervasive disruption of sex chromosome gene expression in sterile hybrids at every stage of spermatogenesis. Failure of MSCI was developmentally preceded by increased silencing of autosomal genes, supporting the hypothesis that divergence at the hybrid incompatibility gene, Prdm9, results in increased rates of autosomal asynapsis which in turn triggers widespread silencing of unsynapsed chromatin. We also detected opposite patterns of postmeiotic overexpression or hyper-repression of the sex chromosomes in reciprocal hybrids, supporting the hypothesis that genomic conflict has driven functional divergence that leads to deleterious X-Y dosage imbalances in hybrids. Our developmental timeline also exposed more subtle patterns of mitotic misregulation on the X chromosome, a previously undocumented stage of spermatogenic disruption in this cross. These results indicate that multiple hybrid incompatibilities have converged on a common regulatory phenotype, the disrupted expression of the sex chromosomes during spermatogenesis. Collectively, these data reveal a composite regulatory basis to hybrid male sterility in mice that helps resolve the mechanistic underpinnings of the well-documented large X-effect in mice speciation. We propose that the inherent sensitivity of spermatogenesis to X-linked regulatory disruption has the potential to be a major driver of reproductive isolation in species with chromosomal sex determination.

Genomic imprinting, disrupted placental expression, and speciation

The importance of regulatory incompatibilities to the early stages of speciation remains unclear. Hybrid mammals often show extreme parent-of-origin growth effects that are thought to be a consequence of disrupted genetic imprinting (parent-specific epigenetic gene silencing) during early development. Here, we test the long-standing hypothesis that abnormal hybrid growth reflects disrupted gene expression due to loss of imprinting (LOI) in hybrid placentas, resulting in dosage imbalances between paternal growth factors and maternal growth repressors. We analyzed placental gene expression in reciprocal dwarf hamster hybrids that show extreme parent-of-origin growth effects relative to their parental species. In massively enlarged hybrid placentas, we observed both extensive transgressive expression of growth-related genes and biallelic expression of many genes that were paternally silenced in normal sized hybrids. However, the apparent widespread disruption of paternal silencing was coupled with reduced gene expression levels overall. These patterns are contrary to the predictions of the LOI model and indicate that hybrid misexpression of dosage-sensitive genes is caused by other regulatory mechanisms in this system. Collectively, our results support a central role for disrupted gene expression and imprinting in the evolution of mammalian hybrid inviability, but call into question the generality of the widely invoked LOI model.

Chromosome-Wide Impacts on the Expression of Incompatibilities in Hybrids of Tigriopus californicus

Chromosome rearrangements such as inversions have been recognized previously as contributing to reproductive isolation by maintaining alleles together that jointly contribute to deleterious genetic interactions and postzygotic reproductive isolation. In this study, an impact of potential incompatibilities merely residing on the same chromosome was found in crosses of populations of the copepod Tigriopus californicus When genetically divergent populations of this copepod are crossed, hybrids show reduced fitness, and deviations from expected genotypic ratios can be used to determine regions of the genome involved in deleterious interactions. In this study, a set of markers was genotyped for a cross of two populations of T. californicus, and these markers show widespread deviations from Mendelian expectations, with entire chromosomes showing marked skew. Despite the importance of mtDNA/nuclear interactions in incompatibilities in this system in previous studies, in these crosses the expected patterns stemming from these interactions are not widely apparent. Females lack recombination in this species, and a striking difference is observed between male and female backcrosses. This suggests that the maintenance of multiple loci on individual chromosomes can enable some incompatibilities, perhaps playing a similar role in the initial rounds of hybridization to chromosomal rearrangements in preserving sets of alleles together that contribute to incompatibilities. Finally, it was observed that candidate pairs of incompatibility regions are not consistently interacting across replicates or subsets of these crosses, despite the repeatability of the deviations at many of the single loci themselves, suggesting that more complicated models of Dobzhansky-Muller incompatibilities may need to be considered.