Widespread misregulation of inter-species hybrid transcriptomes due to sex-specific and sex-chromosome regulatory evolution

When gene regulatory networks diverge between species, their dysfunctional expression in inter-species hybrid individuals can create genetic incompatibilities that generate the developmental defects responsible for intrinsic post-zygotic reproductive isolation. Both cis- and trans-acting regulatory divergence can be hastened by directional selection through adaptation, sexual selection, and inter-sexual conflict, in addition to cryptic evolution under stabilizing selection. Dysfunctional sex-biased gene expression, in particular, may provide an important source of sexually-dimorphic genetic incompatibilities. Here, we characterize and compare male and female/hermaphrodite transcriptome profiles for sibling nematode species Caenorhabditis briggsae and C. nigoni, along with allele-specific expression in their F-1 hybrids, to deconvolve features of expression divergence and regulatory dysfunction. Despite evidence of widespread stabilizing selection on gene expression, misexpression of sex-biased genes pervades F-1 hybrids of both sexes. This finding implicates greater fragility of male genetic networks to produce dysfunctional organismal phenotypes. Spermatogenesis genes are especially prone to high divergence in both expression and coding sequences, consistent with a "faster male" model for Haldane's rule and elevated sterility of hybrid males. Moreover, underdominant expression pervades male-biased genes compared to female-biased and sex-neutral genes and an excess of cis-trans compensatory regulatory divergence for X-linked genes underscores a "large-X effect" for hybrid male expression dysfunction. Extensive regulatory divergence in sex determination pathway genes likely contributes to demasculinization of XX hybrids. The evolution of genetic incompatibilities due to regulatory versus coding sequence divergence, however, are expected to arise in an uncorrelated fashion. This study identifies important differences between the sexes in how regulatory networks diverge to contribute to sex-biases in how genetic incompatibilities manifest during the speciation process. Author summary As species diverge, many mutations that affect traits do so by altering gene expression. Such gene regulatory changes also accumulate in the control of static traits, due to compensatory effects of mutation on multiple regulatory elements. Theory predicts many of these changes to cause inter-species hybrids to experience dysfunctional gene expression that leads to reduced fitness, disproportionately affecting genes biased toward expression in one sex and that localize to sex chromosomes. Our analyses of genome-wide gene expression from Caenorhabditis nematode roundworms support these predictions. We find widespread rewiring of gene regulation, despite the extensive morphological stasis and conserved expression profiles that are hallmarks of these animals. Misregulation of expression in inter-species hybrids of both sexes is most severe for genes linked to the X-chromosome, but male organismal phenotypes are most disrupted in hybrids. This fragility of male genetic networks and sex differences in regulatory evolution of local versus distant elements may underlie feminized and sterile phenotypes among hybrids. Our work clarifies how distinct components of regulatory networks evolve and contribute to sex differences in the manifestation of genetic incompatibilities in the speciation process.