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.

Synthesis and Scope of the Role of Postmating Prezygotic Isolation in Speciation.

How barriers to gene flow arise and are maintained are key questions in evolutionary biology. Speciation research has mainly focused on barriers that occur either before mating or after zygote formation. In comparison, postmating prezygotic (PMPZ) isolation-a barrier that acts after gamete release but before zygote formation-is less frequently investigated but may hold a unique role in generating biodiversity. Here we discuss the distinctive features of PMPZ isolation, including the primary drivers and molecular mechanisms underpinning PMPZ isolation. We then present the first comprehensive survey of PMPZ isolation research, revealing that it is a widespread form of prezygotic isolation across eukaryotes. The survey also exposes obstacles in studying PMPZ isolation, in part attributable to the challenges involved in directly measuring PMPZ isolation and uncovering its causal mechanisms. Finally, we identify outstanding knowledge gaps and provide recommendations for improving future research on PMPZ isolation. This will allow us to better understand the nature of this often-neglected reproductive barrier and its contribution to speciation.

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.

Intravenous Bacille Calmette-Guérin vaccination protects simian immunodeficiency virus-infected macaques from tuberculosis.

Tuberculosis, caused by Mycobacterium tuberculosis (Mtb), is the most common cause of death in people living with human immunodeficiency virus (HIV). Intra-dermal Bacille Calmette-Guérin (BCG) delivery is the only licensed vaccine against tuberculosis; however, it offers little protection from pulmonary tuberculosis in adults and is contraindicated in people living with HIV. Intravenous BCG confers protection against Mtb infection in rhesus macaques; we hypothesized that it might prevent tuberculosis in simian immunodeficiency virus (SIV)-infected macaques, a model for HIV infection. Here intravenous BCG-elicited robust airway T cell influx and elevated plasma and airway antibody titres in both SIV-infected and naive animals. Following Mtb challenge, all 7 vaccinated SIV-naive and 9 out of 12 vaccinated SIV-infected animals were protected, without any culturable bacteria detected from tissues. Peripheral blood mononuclear cell responses post-challenge indicated early clearance of Mtb in vaccinated animals, regardless of SIV infection. These data support that intravenous BCG is immunogenic and efficacious in SIV-infected animals.

Quantitative Analyses of Coupling in Hybrid Zones.

In hybrid zones, whether barrier loci experience selection mostly independently or as a unit depends on the ratio of selection to recombination as captured by the coupling coefficient. Theory predicts a sharper transition between an uncoupled and coupled system when more loci affect hybrid fitness. However, the extent of coupling in hybrid zones has rarely been quantified. Here, we use simulations to characterize the relationship between the coupling coefficient and variance in clines across genetic loci. We then reanalyze 25 hybrid zone data sets and find that cline variances and estimated coupling coefficients form a smooth continuum from high variance and weak coupling to low variance and strong coupling. Our results are consistent with low rates of hybridization and a strong genome-wide barrier to gene flow when the coupling coefficient is much greater than 1, but also suggest that this boundary might be approached gradually and at a near constant rate over time.

The Endogenous Retinoic Acid Receptor Pathway Is Exploited by Mycobacterium tuberculosis during Infection, Both In Vitro and In Vivo.

Retinoic acid (RA) is a fundamental vitamin A metabolite involved in regulating immune responses through the nuclear RA receptor (RAR) and retinoid X receptor. While performing experiments using THP-1 cells as a model for Mycobacterium tuberculosis infection, we observed that serum-supplemented cultures displayed high levels of baseline RAR activation in the presence of live, but not heat-killed, bacteria, suggesting that M. tuberculosis robustly induces the endogenous RAR pathway. Using in vitro and in vivo models, we have further explored the role of endogenous RAR activity in M. tuberculosis infection through pharmacological inhibition of RARs. We found that M. tuberculosis induces classical RA response element genes such as CD38 and DHRS3 in both THP-1 cells and human primary CD14+ monocytes via a RAR-dependent pathway. M. tuberculosis-stimulated RAR activation was observed with conditioned media and required nonproteinaceous factor(s) present in FBS. Importantly, RAR blockade by (4-[(E)-2-[5,5-dimethyl-8-(2-phenylethynyl)-6H-naphthalen-2-yl]ethenyl]benzoic acid), a specific pan-RAR inverse agonist, in a low-dose murine model of tuberculosis significantly reduced SIGLEC-F+CD64+CD11c+high alveolar macrophages in the lungs, which correlated with 2× reduction in tissue mycobacterial burden. These results suggest that the endogenous RAR activation axis contributes to M. tuberculosis infection both in vitro and in vivo and reveal an opportunity for further investigation of new antituberculosis therapies.

Vaccination with intravenous BCG protects macaques with pre-existing SIV infection from tuberculosis.

Tuberculosis (TB) is the most common cause of death in people living with HIV. BCG delivered intradermally (ID) is the only licensed vaccine to prevent TB. However, it offers little protection from pulmonary TB in adults. Intravenous (IV) BCG, but not ID BCG, confers striking protection against Mycobacterium tuberculosis (Mtb) infection and disease in rhesus macaques. We investigated whether IV BCG could protect against TB in macaques with a pre-existing SIV infection. There was a robust influx of airway T cells following IV BCG in both SIV-infected and SIV-naïve animals, with elevated antibody titers in plasma and airways. Following Mtb challenge, all 7 SIV-naïve and 9 out of 12 SIV-infected vaccinated animals were completely protected, without any culturable bacilli in their tissues. PBMC responses post-challenge indicated early clearance of Mtb in vaccinated animals regardless of SIV infection. These data support that IV BCG is immunogenic and efficacious in SIV-infected animals.

Host Immunity to Mycobacterium tuberculosis Infection Is Similar in Simian Immunodeficiency Virus (SIV)-Infected, Antiretroviral Therapy-Treated and SIV-Naïve Juvenile Macaques.

Pre-existing HIV infection increases tuberculosis (TB) risk in children. Antiretroviral therapy (ART) reduces, but does not abolish, this risk in children with HIV. The immunologic mechanisms involved in TB progression in both HIV-naive and HIV-infected children have not been explored. Much of our current understanding is based on human studies in adults and adult animal models. In this study, we sought to model childhood HIV/Mycobacterium tuberculosis (Mtb) coinfection in the setting of ART and characterize T cells during TB progression. Macaques equivalent to 4 to 8 year-old children were intravenously infected with SIVmac239M, treated with ART 3 months later, and coinfected with Mtb 3 months after initiating ART. SIV-naive macaques were similarly infected with Mtb alone. TB pathology and total Mtb burden did not differ between SIV-infected, ART-treated and SIV-naive macaques, although lung Mtb burden was lower in SIV-infected, ART-treated macaques. No major differences in frequencies of CD4+ and CD8+ T cells and unconventional T cell subsets (Vγ9+ γδ T cells, MAIT cells, and NKT cells) in airways were observed between SIV-infected, ART-treated and SIV-naive macaques over the course of Mtb infection, with the exception of CCR5+ CD4+ and CD8+ T cells which were slightly lower. CD4+ and CD8+ T cell frequencies did not differ in the lung granulomas. Immune checkpoint marker levels were similar, although ki-67 levels in CD8+ T cells were elevated. Thus, ART treatment of juvenile macaques, 3 months after SIV infection, resulted in similar progression of Mtb and T cell responses compared to Mtb in SIV-naive macaques.

The contribution of sex chromosome conflict to disrupted spermatogenesis in hybrid house mice.

Incompatibilities on the sex chromosomes are important in the evolution of hybrid male sterility, but the evolutionary forces underlying this phenomenon are unclear. House mice (Mus musculus) lineages have provided powerful models for understanding the genetic basis of hybrid male sterility. X chromosome-autosome interactions cause strong incompatibilities in M. musculus F1 hybrids, but variation in sterility phenotypes suggests a more complex genetic basis. In addition, XY chromosome conflict has resulted in rapid expansions of ampliconic genes with dosage-dependent expression that is essential to spermatogenesis. Here, we evaluated the contribution of XY lineage mismatch to male fertility and stage-specific gene expression in hybrid mice. We performed backcrosses between two house mouse subspecies to generate reciprocal Y-introgression strains and used these strains to test the effects of XY mismatch in hybrids. Our transcriptome analyses of sorted spermatid cells revealed widespread overexpression of the X chromosome in sterile F1 hybrids independent of Y chromosome subspecies origin. Thus, postmeiotic overexpression of the X chromosome in sterile F1 mouse hybrids is likely a downstream consequence of disrupted meiotic X-inactivation rather than XY gene copy number imbalance. Y chromosome introgression did result in subfertility phenotypes and disrupted expression of several autosomal genes in mice with an otherwise nonhybrid genomic background, suggesting that Y-linked incompatibilities contribute to reproductive barriers, but likely not as a direct consequence of XY conflict. Collectively, these findings suggest that rapid sex chromosome gene family evolution driven by genomic conflict has not resulted in strong male reproductive barriers between these subspecies of house mice.

The Evolution of Widespread Recombination Suppression on the Dwarf Hamster (Phodopus) X Chromosome.

The X chromosome of therian mammals shows strong conservation among distantly related species, limiting insights into the distinct selective processes that have shaped sex chromosome evolution. We constructed a chromosome-scale de novo genome assembly for the Siberian dwarf hamster (Phodopus sungorus), a species reported to show extensive recombination suppression across an entire arm of the X chromosome. Combining a physical genome assembly based on shotgun and long-range proximity ligation sequencing with a dense genetic map, we detected widespread suppression of female recombination across ∼65% of the Phodopus X chromosome. This region of suppressed recombination likely corresponds to the Xp arm, which has previously been shown to be highly heterochromatic. Using additional sequencing data from two closely related species (P. campbelli and P. roborovskii), we show that recombination suppression on Xp appears to be independent of major structural rearrangements. The suppressed Xp arm was enriched for several transposable element families and de-enriched for genes primarily expressed in placenta, but otherwise showed similar gene densities, expression patterns, and rates of molecular evolution when compared to the recombinant Xq arm. Phodopus Xp gene content and order was also broadly conserved relative to the more distantly related rat X chromosome. These data suggest that widespread suppression of recombination has likely evolved through the transient induction of facultative heterochromatin on the Phodopus Xp arm without major changes in chromosome structure or genetic content. Thus, substantial changes in the recombination landscape have so far had relatively subtle influences on patterns of X-linked molecular evolution in these species.

Genetic divergence among threespine stickleback that differ in nuptial coloration.

Sexual signals are shaped by their intended and unintended receivers as well as the signalling environment. This interplay between sexual and natural selection can lead to divergence in signals in heterogeneous environments. Yet, the extent to which gene flow is restricted when signalling phenotypes vary across environments and over what spatial scales remains an outstanding question. In this study, we quantify gene flow between two colour morphs, red and black, of freshwater threespine stickleback fish (Gasterosteus aculeatus). We capitalize on the very recent divergence of signalling phenotypes in this system to characterize within-species and among-morph genetic variation and to test for levels of gene flow between colour morphs in Oregon and Washington. Despite limited evidence for assortative mating between allopatric red and black populations, we found that black populations are genetically distinct from nearby red populations and that the black morph appears to have evolved independently at least twice in Oregon and Washington. Surprisingly, we uncovered a group of stickleback in one small coastal stream, Connor Creek, which is genetically and morphologically distinct from the red and black colour morphs and from marine stickleback. Historically, both colour morphs have coexisted in this location and sometimes hybridized, raising new questions about the origins and history of these fish, which were first described as anadromous-black hybrids >50 years ago. Understanding how genetic variation is currently partitioned within and among populations and colour morphs in this system should prompt future studies to assess the relative roles of habitat, ecological and pre- and post-reproductive barriers in the genetic divergence and phenotypic patterns we observe in nature.

Spontaneous Control of SIV Replication Does Not Prevent T Cell Dysregulation and Bacterial Dissemination in Animals Co-Infected with M. tuberculosis.

Individuals co-infected with HIV and Mycobacterium tuberculosis (Mtb) are more likely to develop severe tuberculosis (TB) disease than HIV-naive individuals. To understand how a chronic pre-existing Simian immunodeficiency virus (SIV) infection impairs the early immune response to Mtb, we used the Mauritian cynomolgus macaque (MCM) model of SIV/Mtb co-infection. We examined the relationship between peripheral viral control and Mtb burden, Mtb dissemination, and T cell function between SIV+ spontaneous controllers, SIV+ non-controllers, and SIV-naive MCM who were challenged with a barcoded Mtb Erdman strain 6 months post-SIV infection and necropsied 6 weeks post-Mtb infection. Mycobacterial burden was highest in the SIV+ non-controllers in all assessed tissues. In lung granulomas, the frequency of TNF-α-producing CD4+ T cells was reduced in all SIV+ MCM, but IFNγ-producing CD4+ T cells were only lower in the SIV+ non-controllers. Further, while all SIV+ MCM had more PD1+ and TIGIT+ T cells in the lung granulomas relative to SIV-naive MCM, SIV+ controllers exhibited the highest frequency of cells expressing these markers. To measure the effect of SIV infection on within-host bacterial dissemination, we sequenced the molecular barcodes of Mtb present in each tissue and characterized the Mtb population complexity. While Mtb population complexity was not associated with SIV infection group, lymph nodes had increased complexity when compared with lung granulomas across all groups. These results provide evidence that SIV+ animals, independent of viral control, exhibit a dysregulated T cell immune response and enhanced dissemination of Mtb, likely contributing to the poor TB disease course across all SIV/Mtb co-infected animals. IMPORTANCE HIV and TB remain significant global health issues, despite the availability of treatments. Individuals with HIV, including those who are virally suppressed, are at an increased risk to develop and succumb to severe TB disease when compared with HIV-naive individuals. Our study aims to understand the relationship between the extent of SIV replication, mycobacterial growth, and T cell function in the tissues of co-infected Mauritian cynomolgus macaques during the first 6 weeks of Mtb infection. Here we demonstrate that increased viral replication is associated with increased bacterial burden in the tissues and impaired T cell responses, and that the immunological damage attributed to virus infection is not fully eliminated when animals spontaneously control virus replication.

HIV-1 provirus transcription and translation in macrophages differs from pre-integrated cDNA complexes and requires E2F transcriptional programs.

HIV-1 cDNA pre-integration complexes persist for weeks in macrophages and remain transcriptionally active. While previous work has focused on the transcription of HIV-1 genes; our understanding of the cellular milieu that accompanies viral production is incomplete. We have used an in vitro system to model HIV-1 infection of macrophages, and single-cell RNA sequencing (scRNA-seq) to compare the transcriptomes of uninfected cells, cells harboring pre-integration complexes (PIC), and those containing integrated provirus and making late HIV proteins. scRNA-seq can distinguish between provirus and PIC cells because their background transcriptomes vary dramatically. PIC cell transcriptomes are characterized by NFkB and AP-1 promoted transcription, while transcriptomes of cells transcribing from provirus are characterized by E2F family transcription products. We also find that the transcriptomes of PIC cells and Bystander cells (defined as cells not producing any HIV transcript and thus presumably not infected) are indistinguishable except for the presence of HIV-1 transcripts. Furthermore, the presence of pathogen alters the transcriptome of the uninfected Bystander cells, so that they are distinguishable from true control cells (cells not exposed to any pathogen). Therefore, a single cell comparison of transcriptomes from provirus and PIC cells provides a new understanding of the transcriptional changes that accompany HIV-1 integration.

Molecular Evolution across Mouse Spermatogenesis.

Genes involved in spermatogenesis tend to evolve rapidly, but we lack a clear understanding of how protein sequences and patterns of gene expression evolve across this complex developmental process. We used fluorescence-activated cell sorting (FACS) to generate expression data for early (meiotic) and late (postmeiotic) cell types across 13 inbred strains of mice (Mus) spanning ∼7 My of evolution. We used these comparative developmental data to investigate the evolution of lineage-specific expression, protein-coding sequences, and expression levels. We found increased lineage specificity and more rapid protein-coding and expression divergence during late spermatogenesis, suggesting that signatures of rapid testis molecular evolution are punctuated across sperm development. Despite strong overall developmental parallels in these components of molecular evolution, protein and expression divergences were only weakly correlated across genes. We detected more rapid protein evolution on the X chromosome relative to the autosomes, whereas X-linked gene expression tended to be relatively more conserved likely reflecting chromosome-specific regulatory constraints. Using allele-specific FACS expression data from crosses between four strains, we found that the relative contributions of different regulatory mechanisms also differed between cell types. Genes showing cis-regulatory changes were more common late in spermatogenesis, and tended to be associated with larger differences in expression levels and greater expression divergence between species. In contrast, genes with trans-acting changes were more common early and tended to be more conserved across species. Our findings advance understanding of gene evolution across spermatogenesis and underscore the fundamental importance of developmental context in molecular evolutionary studies.

Stage-specific disruption of X chromosome expression during spermatogenesis in sterile house mouse hybrids.

Hybrid sterility is a complex phenotype that can result from the breakdown of spermatogenesis at multiple developmental stages. Here, we disentangle two proposed hybrid male sterility mechanisms in the house mice, Mus musculus domesticus and M. m. musculus, by comparing patterns of gene expression in sterile F1 hybrids from a reciprocal cross. We found that hybrid males from both cross directions showed disrupted X chromosome expression during prophase of meiosis I consistent with a loss of meiotic sex chromosome inactivation (MSCI) and Prdm9-associated sterility, but that the degree of disruption was greater in mice with an M. m. musculus X chromosome consistent with previous studies. During postmeiotic development, gene expression on the X chromosome was only disrupted in one cross direction, suggesting that misexpression at this later stage was genotype-specific and not a simple downstream consequence of MSCI disruption which was observed in both reciprocal crosses. Instead, disrupted postmeiotic expression may depend on the magnitude of earlier disrupted MSCI, or the disruption of particular X-linked genes or gene networks. Alternatively, only hybrids with a potential deficit of Sly copies, a Y-linked ampliconic gene family, showed overexpression in postmeiotic cells, consistent with a previously proposed model of antagonistic coevolution between the X- and Y-linked ampliconic genes contributing to disrupted expression late in spermatogenesis. The relative contributions of these two regulatory mechanisms and their impact on sterility phenotypes await further study. Our results further support the hypothesis that X-linked hybrid sterility in house mice has a variable genetic basis, and that genotype-specific disruption of gene regulation contributes to overexpression of the X chromosome at different stages of development. Overall, these findings underscore the critical role of epigenetic regulation of the X chromosome during spermatogenesis and suggest that these processes are prone to disruption in hybrids.

Unraveling patterns of disrupted gene expression across a complex tissue.

Whole tissue RNASeq is the standard approach for studying gene expression divergence in evolutionary biology and provides a snapshot of the comprehensive transcriptome for a given tissue. However, whole tissues consist of diverse cell types differing in expression profiles, and the cellular composition of these tissues can evolve across species. Here, we investigate the effects of different cellular composition on whole tissue expression profiles. We compared gene expression from whole testes and enriched spermatogenesis populations in two species of house mice, Mus musculus musculus and M. m. domesticus, and their sterile and fertile F1 hybrids, which differ in both cellular composition and regulatory dynamics. We found that cellular composition differences skewed expression profiles and differential gene expression in whole testes samples. Importantly, both approaches were able to detect large-scale patterns such as disrupted X chromosome expression, although whole testes sampling resulted in decreased power to detect differentially expressed genes. We encourage researchers to account for histology in RNASeq and consider methods that reduce sample complexity whenever feasible. Ultimately, we show that differences in cellular composition between tissues can modify expression profiles, potentially altering inferred gene ontological processes, insights into gene network evolution, and processes governing gene expression evolution.

Pre-existing Simian Immunodeficiency Virus Infection Increases Expression of T Cell Markers Associated with Activation during Early Mycobacterium tuberculosis Coinfection and Impairs TNF Responses in Granulomas.

Tuberculosis (TB) is the leading infectious cause of death among people living with HIV. People living with HIV are more susceptible to contracting Mycobacterium tuberculosis and often have worsened TB disease. Understanding the immunologic defects caused by HIV and the consequences it has on M. tuberculosis coinfection is critical in combating this global health epidemic. We previously showed in a model of SIV and M. tuberculosis coinfection in Mauritian cynomolgus macaques (MCM) that SIV/M. tuberculosis-coinfected MCM had rapidly progressive TB. We hypothesized that pre-existing SIV infection impairs early T cell responses to M. tuberculosis infection. We infected MCM with SIVmac239, followed by coinfection with M. tuberculosis Erdman 6 mo later. Although similar, TB progression was observed in both SIV+ and SIV-naive animals at 6 wk post-M. tuberculosis infection; longitudinal sampling of the blood (PBMC) and airways (bronchoalveolar lavage) revealed a significant reduction in circulating CD4+ T cells and an influx of CD8+ T cells in airways of SIV+ animals. At sites of M. tuberculosis infection (i.e., granulomas), SIV/M. tuberculosis-coinfected animals had a higher proportion of CD4+ and CD8+ T cells expressing PD-1 and TIGIT. In addition, there were fewer TNF-producing CD4+ T cells in granulomas of SIV/M. tuberculosis-coinfected animals. Taken together, we show that concurrent SIV infection alters T cell phenotypes in granulomas during the early stages of TB disease. As it is critical to establish control of M. tuberculosis replication soon postinfection, these phenotypic changes may distinguish the immune dysfunction that arises from pre-existing SIV infection, which promotes TB progression.

Homage to Felsenstein 1981, or why are there so few/many species?

If there are no constraints on the process of speciation, then the number of species might be expected to match the number of available niches and this number might be indefinitely large. One possible constraint is the opportunity for allopatric divergence. In 1981, Felsenstein used a simple and elegant model to ask if there might also be genetic constraints. He showed that progress towards speciation could be described by the build-up of linkage disequilibrium among divergently selected loci and between these loci and those contributing to other forms of reproductive isolation. Therefore, speciation is opposed by recombination, because it tends to break down linkage disequilibria. Felsenstein then introduced a crucial distinction between "two-allele" models, which are subject to this effect, and "one-allele" models, which are free from the recombination constraint. These fundamentally important insights have been the foundation for both empirical and theoretical studies of speciation ever since.