Regulatory divergence in Drosophila revealed by mRNA-seq

Evolution of the regulation of developmental gene expression in blind Mexican cavefish

Developmental evolution and diversification of morphology can arise through changes in the regulation of gene expression or protein-coding sequence. To unravel mechanisms underlying early developmental evolution in cavefish of the species Astyanax mexicanus, we compared transcriptomes of surface-dwelling and blind cave-adapted morphs at the end of gastrulation. Twenty percent of the transcriptome was differentially expressed. Allelic expression ratios in cave X surface hybrids showed that cis-regulatory changes are the quasi-exclusive contributors to inter-morph variations in gene expression. Among a list of 108 genes with change at the cis-regulatory level, we explored the control of expression of rx3, which is a master eye gene. We discovered that cellular rx3 levels are cis-regulated in a cell-autonomous manner, whereas rx3 domain size depends on non-autonomous Wnt and Bmp signalling. These results highlight how uncoupled mechanisms and regulatory modules control developmental gene expression and shape morphological changes. Finally, a transcriptome-wide search for fixed coding mutations and differential exon use suggested that variations in coding sequence have a minor contribution. Thus, during early embryogenesis, changes in gene expression regulation are the main drivers of cavefish developmental evolution.

The role of heterochronic gene expression and regulatory architecture in early developmental divergence

New developmental programs can evolve through adaptive changes to gene expression. The annelid Streblospio benedicti has a developmental dimorphism, which provides a unique intraspecific framework for understanding the earliest genetic changes that take place during developmental divergence. Using comparative RNAseq through ontogeny, we find that only a small proportion of genes are differentially expressed at any time, despite major differences in larval development and life history. These genes shift expression profiles across morphs by either turning off any expression in one morph or changing the timing or amount of gene expression. We directly connect the contributions of these mechanisms to differences in developmental processes. We examine F1 offspring - using reciprocal crosses - to determine maternal mRNA inheritance and the regulatory architecture of gene expression. These results highlight the importance of both novel gene expression and heterochronic shifts in developmental evolution, as well as the trans-acting regulatory factors in initiating divergence.

Pervasive tissue-, genetic background-, and allele-specific gene expression effects in Drosophila melanogaster

The pervasiveness of gene expression variation and its contribution to phenotypic variation and evolution is well known. This gene expression variation is context dependent, with differences in regulatory architecture often associated with intrinsic and environmental factors, and is modulated by regulatory elements that can act in cis (linked) or in trans (unlinked) relative to the genes they affect. So far, little is known about how this genetic variation affects the evolution of regulatory architecture among closely related tissues during population divergence. To address this question, we analyzed gene expression in the midgut, hindgut, and Malpighian tubule as well as microbiome composition in the two gut tissues in four Drosophila melanogaster strains and their F1 hybrids from two divergent populations: one from the derived, European range and one from the ancestral, African range. In both the transcriptome and microbiome data, we detected extensive tissue- and genetic background-specific effects, including effects of genetic background on overall tissue specificity. Tissue-specific effects were typically stronger than genetic background-specific effects, although the two gut tissues were not more similar to each other than to the Malpighian tubules. An examination of allele specific expression revealed that, while both cis and trans effects were more tissue-specific in genes expressed differentially between populations than genes with conserved expression, trans effects were more tissue-specific than cis effects. Despite there being highly variable regulatory architecture, this observation was robust across tissues and genetic backgrounds, suggesting that the expression of trans variation can be spatially fine-tuned as well as or better than cis variation during population divergence and yielding new insights into cis and trans regulatory evolution.

Evolution and genetics of accessory gland transcriptome divergence between Drosophila melanogaster and D. simulans

Studies of allele-specific expression in interspecific hybrids have provided important insights into gene-regulatory divergence and hybrid incompatibilities. Many such investigations in Drosophila have used transcriptome data from complex mixtures of many tissues or from gonads, however, regulatory divergence may vary widely among species, sexes, and tissues. Thus, we lack sufficiently broad sampling to be confident about the general biological principles of regulatory divergence. Here, we seek to fill some of these gaps in the literature by characterizing regulatory evolution and hybrid misexpression in a somatic male sex organ, the accessory gland, in F1 hybrids between Drosophila melanogaster and D. simulans. The accessory gland produces seminal fluid proteins, which play an important role in male and female fertility and may be subject to adaptive divergence due to male-male or male-female interactions. We find that trans differences are relatively more abundant than cis, in contrast to most of the interspecific hybrid literature, though large effect-size trans differences are rare. Seminal fluid protein genes have significantly elevated levels of expression divergence and tend to be regulated through both cis and trans divergence. We find limited misexpression (over- or underexpression relative to both parents) in this organ compared to most other Drosophila studies. As in previous studies, male-biased genes are overrepresented among misexpressed genes and are much more likely to be underexpressed. ATAC-Seq data show that chromatin accessibility is correlated with expression differences among species and hybrid allele-specific expression. This work identifies unique regulatory evolution and hybrid misexpression properties of the accessory gland and suggests the importance of tissue-specific allele-specific expression studies.

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.

Human gene regulatory evolution is driven by the divergence of regulatory element function in both cis and trans

Gene regulatory divergence between species can result from cis-acting local changes to regulatory element DNA sequences or global trans-acting changes to the regulatory environment. Understanding how these mechanisms drive regulatory evolution has been limited by challenges in identifying trans-acting changes. We present a comprehensive approach to directly identify cis- and trans-divergent regulatory elements between human and rhesus macaque lymphoblastoid cells using assay for transposase-accessible chromatin coupled to self-transcribing active regulatory region (ATAC-STARR) sequencing. In addition to thousands of cis changes, we discover an unexpected number (∼10,000) of trans changes and show that cis and trans elements exhibit distinct patterns of sequence divergence and function. We further identify differentially expressed transcription factors that underlie ∼37% of trans differences and trace how cis changes can produce cascades of trans changes. Overall, we find that most divergent elements (67%) experienced changes in both cis and trans, revealing a substantial role for trans divergence-alone and together with cis changes-in regulatory differences between species.

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.

Tissue-specific and cis-regulatory changes underlie parallel, adaptive gene expression evolution in house mice

Changes in gene regulation have long been appreciated as a driving force of adaptive evolution, however the relative contributions of cis- and trans-acting changes to gene regulation over short evolutionary timescales remain unclear. Instances of recent, parallel phenotypic evolution provide an opportunity to assess whether parallel patterns are seen at the level of gene expression, and to assess the relative contribution of cis- and trans- changes to gene regulation in the early stages of divergence. Here, we studied gene expression in liver and brown adipose tissue in two wild-derived strains of house mice that independently adapted to cold, northern environments, and we compared them to a strain of house mice from a warm, tropical environment. To investigate gene regulatory evolution, we studied expression in parents and allele-specific expression in F1 hybrids of crosses between warm-adapted and cold-adapted strains. First, we found that the different cold-adapted mice showed both unique and shared changes in expression, but that the proportion of shared changes (i.e. parallelism) was greater than expected by chance. Second, we discovered that expression evolution occurred largely at tissue-specific and cis-regulated genes, and that these genes were over-represented in parallel cases of evolution. Finally, we integrated the expression data with scans for selection in natural populations and found substantial parallelism in the two northern populations for genes under selection. Furthermore, selection outliers were associated with cis-regulated genes more than expected by chance; cis-regulated genes under selection influenced phenotypes such as body size, immune functioning, and activity level. These results demonstrate that parallel patterns of gene expression in mice that have independently adapted to cold environments are driven largely by tissue-specific and cis-regulatory changes, providing insight into the mechanisms of adaptive gene regulatory evolution at the earliest stages of divergence.

Investigation of regulatory divergence between homoeologs in the recently formed allopolyploids, Tragopogon mirus and T. miscellus (Asteraceae)

Polyploidy is an important evolutionary process throughout eukaryotes, particularly in flowering plants. Duplicated gene pairs (homoeologs) in allopolyploids provide additional genetic resources for changes in molecular, biochemical, and physiological mechanisms that result in evolutionary novelty. Therefore, understanding how divergent genomes and their regulatory networks reconcile is vital for unraveling the role of polyploidy in plant evolution. Here, we compared the leaf transcriptomes of recently formed natural allotetraploids (Tragopogon mirus and T. miscellus) and their diploid parents (T. porrifolius X T. dubius and T. pratensis X T. dubius, respectively). Analysis of 35 400 expressed loci showed a significantly higher level of transcriptomic additivity compared to old polyploids; only 22% were non-additively expressed in the polyploids, with 5.9% exhibiting transgressive expression (lower or higher expression in the polyploids than in the diploid parents). Among approximately 7400 common orthologous regions (COREs), most loci in both allopolyploids exhibited expression patterns that were vertically inherited from their diploid parents. However, 18% and 20.3% of the loci showed novel expression bias patterns in T. mirus and T. miscellus, respectively. The expression changes of 1500 COREs were explained by cis-regulatory divergence (the condition in which the two parental subgenomes do not interact) between the diploid parents, whereas only about 423 and 461 of the gene expression changes represent trans-effects (the two parental subgenomes interact) in T. mirus and T. miscellus, respectively. The low degree of both non-additivity and trans-effects on gene expression may present the ongoing evolutionary processes of the newly formed Tragopogon polyploids (~80-90 years).

Evolution of gene expression across brain regions in behaviourally divergent deer mice

The evolution of innate behaviours is ultimately due to genetic variation likely acting in the nervous system. Gene regulation may be particularly important because it can evolve in a modular brain-region specific fashion through the concerted action of cis- and trans-regulatory changes. Here, to investigate transcriptional variation and its regulatory basis across the brain, we perform RNA sequencing (RNA-Seq) on ten brain subregions in two sister species of deer mice (Peromyscus maniculatus and P. polionotus)-which differ in a range of innate behaviours, including their social system-and their F1 hybrids. We find that most of the variation in gene expression distinguishes subregions, followed by species. Interspecific differential expression (DE) is pervasive (52-59% of expressed genes), whereas the number of DE genes between sexes is modest overall (~3%). Interestingly, the identity of DE genes varies considerably across brain regions. Much of this modularity is due to cis-regulatory divergence, and while 43% of genes were consistently assigned to the same gene regulatory class across subregions (e.g. conserved, cis- or trans-regulatory divergence), a similar number were assigned to two or more different gene regulatory classes. Together, these results highlight the modularity of gene expression differences and divergence in the brain, which may be key to explain how the evolution of brain gene expression can contribute to the astonishing diversity of animal behaviours.

Non-Additive and Asymmetric Allelic Expression of p38 mapk in Hybrid Tilapia (Oreochromis niloticus ♀ × O. aureus ♂)

Hybridization is a widely used breeding technique in fish species that enhances desirable traits in cultured species through heterosis. However, the mechanism by which hybrids alter gene expression to form heterosis remains unclear. In this study, a group of hybrid tilapia was used to elucidate heterosis through interspecies crossing. Specifically, p38 was analyzed to describe the regulation of gene expression variation in hybrid tilapia. Transcripts from the Nile tilapia allele were found to be significantly higher than those from the blue tilapia allele in hybrid individuals, indicating that the expression of p38 was dominated by Nile tilapia sub-genomic alleles. The study also found a compensatory interaction of cis- and trans-acting elements of the Nile tilapia and blue tilapia sub-genomes, inducing a non-additive expression of p38 in hybrids. Eight specific SNPs were identified in the p38 promoter regions of Nile tilapia and blue tilapia, and were found to be promoter differences leading to differences in gene expression efficiencies between parental alleles using a dual-luciferase reporter system. This study provides insights into the non-additive expression patterns of key functional genes in fish hybrids related to growth and immunity response.

Non-additive genetic components contribute significantly to population-wide gene expression variation

Gene expression variation, an essential step between genotype and phenotype, is collectively controlled by local (cis) and distant (trans) regulatory changes. Nevertheless, how these regulatory elements differentially influence gene expression variation remains unclear. Here, we bridge this gap by analyzing the transcriptomes of a large diallel panel consisting of 323 unique hybrids originating from genetically divergent Saccharomyces cerevisiae isolates. Our analysis across 5,087 transcript abundance traits showed that non-additive components account for 36% of the gene expression variance on average. By comparing allele-specific read counts in parent-hybrid trios, we found that trans-regulatory changes underlie the majority of gene expression variation in the population. Remarkably, most cis-regulatory variations are also exaggerated or attenuated by additional trans effects. Overall, we showed that the transcriptome is globally buffered at the genetic level mainly due to trans-regulatory variation in the population.

Hybrid Epigenomes Reveal Extensive Local Genetic Changes to Chromatin Accessibility Contribute to Divergence in Embryonic Gene Expression Between Species

Chromatin accessibility plays an important role in shaping gene expression, yet little is known about the genetic and molecular mechanisms that influence the evolution of chromatin configuration. Both local (cis) and distant (trans) genetic influences can in principle influence chromatin accessibility and are based on distinct molecular mechanisms. We, therefore, sought to characterize the role that each of these plays in altering chromatin accessibility in 2 closely related sea urchin species. Using hybrids of Heliocidaris erythrogramma and Heliocidaris tuberculata, and adapting a statistical framework previously developed for the analysis of cis and trans influences on the transcriptome, we examined how these mechanisms shape the regulatory landscape at 3 important developmental stages, and compared our results to similar analyses of the transcriptome. We found extensive cis- and trans-based influences on evolutionary changes in chromatin, with cis effects generally larger in effect. Evolutionary changes in accessibility and gene expression are correlated, especially when expression has a local genetic basis. Maternal influences appear to have more of an effect on chromatin accessibility than on gene expression, persisting well past the maternal-to-zygotic transition. Chromatin accessibility near gene regulatory network genes appears to be distinctly regulated, with trans factors appearing to play an outsized role in the configuration of chromatin near these genes. Together, our results represent the first attempt to quantify cis and trans influences on evolutionary divergence in chromatin configuration in an outbred natural study system and suggest that chromatin regulation is more genetically complex than was previously appreciated.

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.

Environmentally robust cis-regulatory changes underlie rapid climatic adaptation

Changes in gene expression are thought to play a major role in adaptive evolution. While it is known that gene expression is highly sensitive to the environment, very few studies have determined the influence of genetic and environmental effects on adaptive gene expression differences in natural populations. Here, we utilize allele-specific expression to characterize cis and trans gene regulatory divergence in temperate and tropical house mice in two metabolic tissues under two thermal conditions. First, we show that gene expression divergence is pervasive between populations and across thermal conditions, with roughly 5 to 10% of genes exhibiting genotype-by-environment interactions. Second, we found that most expression divergence was due to cis-regulatory changes that were stable across temperatures. In contrast, patterns of expression plasticity were largely attributable to trans-effects, which showed greater sensitivity to temperature. Nonetheless, we found a small subset of temperature-dependent cis-regulatory changes, thereby identifying loci underlying expression plasticity. Finally, we performed scans for selection in wild house mice to identify genomic signatures of rapid adaptation. Genomic outliers were enriched in genes with evidence for cis-regulatory divergence. Notably, these genes were associated with phenotypes that affected body weight and metabolism, suggesting that cis-regulatory changes are a possible mechanism for adaptive body size evolution between populations. Our results show that gene expression plasticity, largely controlled in trans, may facilitate the colonization of new environments, but that evolved changes in gene expression are largely controlled in cis, illustrating the genetic and nongenetic mechanisms underlying the establishment of populations in new environments.

Sex-specific estimation of cis and trans regulation of gene expression in heads and gonads of Drosophila melanogaster

The regulatory architecture of gene expression is known to differ substantially between sexes in Drosophila, but most studies performed so far used whole-body data and only single crosses, which may have limited their scope to detect patterns that are robust across tissues and biological replicates. Here, we use allele-specific gene expression of parental and reciprocal hybrid crosses between 6 Drosophila melanogaster inbred lines to quantify cis- and trans-regulatory variation in heads and gonads of both sexes separately across 3 replicate crosses. Our results suggest that female and male heads, as well as ovaries, have a similar regulatory architecture. On the other hand, testes display more and substantially different cis-regulatory effects, suggesting that sex differences in the regulatory architecture that have been previously observed may largely derive from testis-specific effects. We also examine the difference in cis-regulatory variation of genes across different levels of sex bias in gonads and heads. Consistent with the idea that intersex correlations constrain expression and can lead to sexual antagonism, we find more cis variation in unbiased and moderately biased genes in heads. In ovaries, reduced cis variation is observed for male-biased genes, suggesting that cis variants acting on these genes in males do not lead to changes in ovary expression. Finally, we examine the dominance patterns of gene expression and find that sex- and tissue-specific patterns of inheritance as well as trans-regulatory variation are highly variable across biological crosses, although these were performed in highly controlled experimental conditions. This highlights the importance of using various genetic backgrounds to infer generalizable patterns.

Non-additive genetic components contribute significantly to population-wide gene expression variation

Gene expression variation, an essential step between genomic variation and phenotypic landscape, is collectively controlled by local (cis) and distant (trans) regulatory changes. Nevertheless, how these regulatory elements differentially influence the heritability of expression traits remains unclear. Here, we bridge this gap by analyzing the transcriptomes of a large diallel panel consisting of 323 unique hybrids originated from genetically divergent yeast isolates. We estimated the broad- and narrow-sense heritability across 5,087 transcript abundance traits and showed that non-additive components account for 36% of the phenotypic variance on average. By comparing allelic expression ratios in the hybrid and the corresponding parental pair, we identified regulatory changes in 25% of all cases, with a majority acting in trans. We further showed that trans-regulation could underlie coordinated expression variation across highly connected genes, resulting in significantly higher non-additive variance and most likely in some of the missing heritability of gene expression traits.

Speciation and development

Understanding general principles about the origin of species remains one of the foundational challenges in evolutionary biology. The genomic divergence between groups of individuals can spawn hybrid inviability and hybrid sterility, which presents a tantalizing developmental problem. Divergent developmental programs may yield either conserved or divergent phenotypes relative to ancestral traits, both of which can be responsible for reproductive isolation during the speciation process. The genetic mechanisms of developmental evolution involve cis- and trans-acting gene regulatory change, protein-protein interactions, genetic network structures, dosage, and epigenetic regulation, all of which also have roots in population genetic and molecular evolutionary processes. Toward the goal of demystifying Darwin's "mystery of mysteries," this review integrates microevolutionary concepts of genetic change with principles of organismal development, establishing explicit links between population genetic process and developmental mechanisms in the production of macroevolutionary pattern. This integration aims to establish a more unified view of speciation that binds process and mechanism.

Contributions of mutation and selection to regulatory variation: lessons from the Saccharomyces cerevisiae TDH3 gene

Heritable variation in gene expression is common within and among species and contributes to phenotypic diversity. Mutations affecting either cis- or trans-regulatory sequences controlling gene expression give rise to variation in gene expression, and natural selection acting on this variation causes some regulatory variants to persist in a population for longer than others. To understand how mutation and selection interact to produce the patterns of regulatory variation we see within and among species, my colleagues and I have been systematically determining the effects of new mutations on expression of the TDH3 gene in Saccharomyces cerevisiae and comparing them to the effects of polymorphisms segregating within this species. We have also investigated the molecular mechanisms by which regulatory variants act. Over the past decade, this work has revealed properties of cis- and trans-regulatory mutations including their relative frequency, effects, dominance, pleiotropy and fitness consequences. Comparing these mutational effects to the effects of polymorphisms in natural populations, we have inferred selection acting on expression level, expression noise and phenotypic plasticity. Here, I summarize this body of work and synthesize its findings to make inferences not readily discernible from the individual studies alone. This article is part of the theme issue 'Interdisciplinary approaches to predicting evolutionary biology'.

Evaluation of MYBL1 as the master regulator for pachytene spermatocyte genes dysregulated in interspecific hybrid dzo

Misregulation of spermatogenesis transcription factors (TF) in hybrids can lead to misexpression, which is a mechanism for hybrid male sterility (HMS). We used dzo (male offspring of Bos taurus ♂ × Bos grunniens ♀) in bovines to investigate the relationship of the key TF with HMS via RNA sequencing and assay for transposase-accessible chromatin with high-throughput sequencing analyses. RNA sequencing showed that the widespread misexpression in dzo was associated with spermatogenesis-related genes and somatic or progenitor genes. The transition from leptotene or zygotene spermatocytes to pachytene spermatocytes may be the key stage for meiosis arrest in dzo. The analysis of TF-binding motif enrichment revealed that the male meiosis-specific master TF MYB proto-oncogene like 1 (MYBL1, known as A-MYB) motif was enriched on the promoters of downregulated pachytene spermatocyte genes in dzo. Assay for transposase-accessible chromatin with high-throughput sequencing revealed that TF-binding sites for MYBL1, nuclear transcription factor Y, and regulatory factor X were enriched in the low-chromatin accessibility region of dzo. The target genes of the MYBL1-binding motif were associated with meiosis-specific genes and significantly downregulated in dzo testis. The transcription factor MYBL1 may be the candidate master regulator for pachytene spermatocyte genes dysregulated in interspecific HMS dzo. This study reported that a few upstream TF regulation changes might exert a cascading effect downstream in a regulatory network as a mechanism for HMS.

Maternal patterns of inheritance alter transcript expression in eggs

BACKGROUND

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

RESULTS

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

CONCLUSION

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

Expression inheritance and constraints on cis- and trans-regulatory mutations underlying lotus color variation

Both cis- and trans-regulatory mutations drive changes in gene expression that underpin plant phenotypic evolution. However, how and why these two major types of regulatory mutations arise in different genes and how gene expression is inherited and associated with these regulatory changes are unclear. Here, by studying allele-specific expression in F1 hybrids of pink-flowered sacred lotus (Nelumbo nucifera) and yellow-flowered American lotus (N. lutea), we reveal the relative contributions of cis- and trans-regulatory changes to interspecific expression rewiring underlying petal color change and how the expression is inherited in hybrids. Although cis-only variants influenced slightly more genes, trans-only variants had a stronger impact on expression differences between species. In F1 hybrids, genes under cis-only and trans-only regulatory effects showed a propensity toward additive and dominant inheritance, respectively, whereas transgressive inheritance was observed in genes carrying both cis- and trans-variants acting in opposite directions. By investigating anthocyanin and carotenoid coexpression networks in petals, we found that the same category of regulatory mutations, particularly trans-variants, tend to rewire hub genes in coexpression modules underpinning flower color differentiation between species; we identified 45 known genes with cis- and trans-regulatory variants significantly correlated with flower coloration, such as ANTHOCYANIN 5-AROMATIC ACYLTRANSFERASE (ACT), GLUTATHIONE S-TRANSFERASE F11 (GSTF11), and LYCOPENE Ε-CYCLASE (LCYE). Notably, the relative abundance of genes in different categories of regulatory divergence was associated with the inferred magnitude of constraints like expression level and breadth. Overall, our study suggests distinct selective constraints and modes of gene expression inheritance among different regulatory mutations underlying lotus petal color divergence.

Human gene regulatory evolution is driven by the divergence of regulatory element function in both cis and trans

Gene regulatory divergence between species can result from cis-acting local changes to regulatory element DNA sequences or global trans-acting changes to the regulatory environment. Understanding how these mechanisms drive regulatory evolution has been limited by challenges in identifying trans-acting changes. We present a comprehensive approach to directly identify cis- and trans-divergent regulatory elements between human and rhesus macaque lymphoblastoid cells using ATAC-STARR-seq. In addition to thousands of cis changes, we discover an unexpected number (~10,000) of trans changes and show that cis and trans elements exhibit distinct patterns of sequence divergence and function. We further identify differentially expressed transcription factors that underlie >50% of trans differences and trace how cis changes can produce cascades of trans changes. Overall, we find that most divergent elements (67%) experienced changes in both cis and trans, revealing a substantial role for trans divergence-alone and together with cis changes-to regulatory differences between species.

New cup out of old coffee: contribution of parental gene expression legacy to phenotypic novelty in coffee beans of the allopolyploid Coffea arabica L

BACKGROUND AND AIMS

Allopolyploidization is a widespread phenomenon known to generate novel phenotypes by merging evolutionarily distinct parental genomes and regulatory networks in a single nucleus. The objective of this study was to investigate the transcriptional regulation associated with phenotypic novelty in coffee beans of the allotetraploid Coffea arabica.

METHODS

A genome-wide comparative transcriptomic analysis was performed in C. arabica and its two diploid progenitors, C. canephora and C. eugenioides. Gene expression patterns and homeologue expression were studied on seeds at five different maturation stages. The involvement of homeologue expression bias (HEB) in specific traits was addressed both by functional enrichment analyses and by the study of gene expression in the caffeine and chlorogenic acid biosynthesis pathways.

KEY RESULTS

Expression-level dominance in C. arabica seed was observed for most of the genes differentially expressed between the species. Approximately a third of the genes analysed showed HEB. This proportion increased during seed maturation but the biases remained equally distributed between the sub-genomes. The relative expression levels of homeologues remained relatively constant during maturation and were correlated with those estimated in leaves of C. arabica and interspecific hybrids between C. canephora and C. eugenioides. Functional enrichment analyses performed on genes exhibiting HEB enabled the identification of processes potentially associated with physiological traits. The expression profiles of the genes involved in caffeine biosynthesis mirror the differences observed in the caffeine content of mature seeds of C. arabica and its parental species.

CONCLUSIONS

Neither of the two sub-genomes is globally preferentially expressed in C. arabica seeds, and homeologues appear to be co-regulated by shared trans-regulatory mechanisms. The observed HEBs are thought to be a legacy of gene expression differences inherited from diploid progenitor species. Pre-existing functional divergences between parental species appear to play an important role in controlling the phenotype of C. arabica seeds.

Hybrid epigenomes reveal extensive local genetic changes to chromatin accessibility contribute to divergence in embryonic gene expression between species

Chromatin accessibility plays an important role in shaping gene expression patterns across development and evolution; however, little is known about the genetic and molecular mechanisms that influence chromatin configuration itself. Because cis and trans influences can both theoretically influence the accessibility of the epigenome, we sought to better characterize the role that both of these mechanisms play in altering chromatin accessibility in two closely related sea urchin species. Using hybrids of the two species, and adapting a statistical framework previously developed for the analysis of cis and trans influences on the transcriptome, we examined how these mechanisms shape the regulatory landscape at three important developmental stages, and compared our results to similar patterns in the transcriptome. We found extensive cis- and trans-based influences on evolutionary changes in chromatin, with cis effects slightly more numerous and larger in effect. Genetic mechanisms influencing gene expression and chromatin configuration are correlated, but differ in several important ways. Maternal influences also appear to have more of an effect on chromatin accessibility than on gene expression, persisting well past the maternal-to-zygotic transition. Furthermore, chromatin accessibility near GRN genes appears to be regulated differently than the rest of the epigenome, and indicates that trans factors may play an outsized role in the configuration of chromatin near these genes. Together, our results represent the first attempt to quantify cis and trans influences on evolutionary divergence in chromatin configuration in an outbred natural study system, and suggest that the regulation of chromatin is more genetically complex than was previously appreciated.

Heterogeneous Genomic Divergence Landscape in Two Commercially Important European Scallop Species

Two commercially important scallop species of the genus Pecten are found in Europe: the north Atlantic Pecten maximus and the Mediterranean Pecten jacobaeus whose distributions abut at the Almeria-Orán front. Whilst previous studies have quantified genetic divergence between these species, the pattern of differentiation along the Pecten genome is unknown. Here, we mapped RADseq data from 235 P. maximus and 27 P. jacobaeus to a chromosome-level reference genome, finding a heterogeneous landscape of genomic differentiation. Highly divergent genomic regions were identified across 14 chromosomes, while the remaining five showed little differentiation. Demographic and comparative genomics analyses suggest that this pattern resulted from an initial extended period of isolation, which promoted divergence, followed by differential gene flow across the genome during secondary contact. Single nucleotide polymorphisms present within highly divergent genomic regions were located in areas of low recombination and contrasting patterns of LD decay were found between the two species, hinting at the presence of chromosomal inversions in P. jacobaeus. Functional annotations revealed that highly differentiated regions were enriched for immune-related processes and mRNA modification. While future work is necessary to characterize structural differences, this study provides new insights into the speciation genomics of P. maximus and P. jacobaeus.

Network Topology Can Explain Differences in Pleiotropy Between Cis- and Trans-regulatory Mutations

A mutation's degree of pleiotropy (i.e., the number of traits it alters) is predicted to impact the probability of the mutation being detrimental to fitness. For mutations that impact gene expression, mutations acting in cis have been hypothesized to generally be less pleiotropic than mutations affecting the same gene's expression in trans, suggesting that cis-regulatory mutations should be less deleterious and more likely to fix over evolutionary time. Here, we use expression and fitness data from Saccharomyces cerevisiae gene deletion strains to test these hypotheses. By treating deletion of each gene as a cis-regulatory mutation affecting its own expression and deletions of other genes affecting expression of this focal gene as trans-regulatory mutations, we find that cis-acting mutations do indeed tend to be less pleiotropic than trans-acting mutations affecting expression of the same gene. This pattern was observed for the vast majority of genes in the data set and could be explained by the topology of the regulatory network controlling gene expression. Comparing the fitness of cis- and trans-acting mutations affecting expression of the same gene also confirmed that trans-acting deletions tend to be more deleterious. These findings provide strong support for pleiotropy playing a role in the preferential fixation of cis-regulatory alleles over evolutionary time.

Trans-driven variation in expression is common among detoxification genes in the extreme generalist herbivore Tetranychus urticae

The extreme adaptation potential of the generalist herbivore Tetranychus urticae (the two-spotted spider mite) to pesticides as well as diverse host plants has been associated with clade-specific gene expansions in known detoxifying enzyme families, and with extensive and rapid transcriptional responses. However, how this broad transcriptional potential is regulated remains largely unknown. Using a parental/F1 design in which four inbred strains were crossed to a common inbred strain, we assessed the genetic basis and inheritance of gene expression variation in T. urticae. Mirroring known phenotypic variation in the progenitor strains of the inbreds, we confirmed that the inbred strains we created were genetically distinct, varied markedly in pesticide resistance, and also captured variation in host plant fitness as is commonly observed in this species. By examining differences in gene expression between parents and allele-specific expression in F1s, we found that variation in RNA abundance was more often explained in trans as compared to cis, with the former associated with dominance in inheritance. Strikingly, in a gene ontology analysis, detoxification genes of the cytochrome P450 monooxygenase (CYP) family, as well as dioxygenases (DOGs) acquired from horizontal gene transfer from fungi, were specifically enriched at the extremes of trans-driven up- and downregulation. In particular, multiple CYPs and DOGs with broad substrate-specificities for pesticides or plant specialized compounds were exceptionally highly upregulated as a result of trans-regulatory variation, or in some cases synergism of cis and trans, in the most multi-pesticide resistant strains. Collectively, our findings highlight the potential importance of trans-driven expression variation in genes associated with xenobiotic metabolism and host plant use for rapid adaptation in T. urticae, and also suggests modular control of these genes, a regulatory architecture that might ameliorate negative pleiotropic effects.

Trans-regulatory changes underpin the evolution of the Drosophila immune response

When an animal is infected, the expression of a large suite of genes is changed, resulting in an immune response that can defend the host. Despite much evidence that the sequence of proteins in the immune system can evolve rapidly, the evolution of gene expression is comparatively poorly understood. We therefore investigated the transcriptional response to parasitoid wasp infection in Drosophila simulans and D. sechellia. Although these species are closely related, there has been a large scale divergence in the expression of immune-responsive genes in their two main immune tissues, the fat body and hemocytes. Many genes, including those encoding molecules that directly kill pathogens, have cis regulatory changes, frequently resulting in large differences in their expression in the two species. However, these changes in cis regulation overwhelmingly affected gene expression in immune-challenged and uninfected animals alike. Divergence in the response to infection was controlled in trans. We argue that altering trans-regulatory factors, such as signalling pathways or immune modulators, may allow natural selection to alter the expression of large numbers of immune-responsive genes in a coordinated fashion.

Symmetric subgenomes and balanced homoeolog expression stabilize the establishment of allopolyploidy in cyprinid fish

BACKGROUND

Interspecific postzygotic reproduction isolation results from large genetic divergence between the subgenomes of established hybrids. Polyploidization immediately after hybridization may reset patterns of homologous chromosome pairing and ameliorate deleterious genomic incompatibility between the subgenomes of distinct parental species in plants and animals. However, the observation that polyploidy is less common in vertebrates raises the question of which factors restrict its emergence. Here, we perform analyses of the genome, epigenome, and gene expression in the nascent allotetraploid lineage (2.95 Gb) derived from the intergeneric hybridization of female goldfish (Carassius auratus, 1.49 Gb) and male common carp (Cyprinus carpio, 1.42 Gb), to shed light on the changes leading to the stabilization of hybrids.

RESULTS

We firstly identify the two subgenomes derived from the parental lineages of goldfish and common carp. We find variable unequal homoeologous recombination in somatic and germ cells of the intergeneric F1 and allotetraploid (F22 and F24) populations, reflecting high plasticity between the subgenomes, and rapidly varying copy numbers between the homoeolog genes. We also find dynamic changes in transposable elements accompanied by genome merger and duplication in the allotetraploid lineage. Finally, we observe the gradual decreases in cis-regulatory effects and increases in trans-regulatory effects along with the allotetraploidization, which contribute to increases in the symmetrical homoeologous expression in different tissues and developmental stages, especially in early embryogenesis.

CONCLUSIONS

Our results reveal a series of changes in transposable elements, unequal homoeologous recombination, cis- and trans-regulations (e.g. DNA methylation), and homoeologous expression, suggesting their potential roles in mediating adaptive stabilization of regulatory systems of the nascent allotetraploid lineage. The symmetrical subgenomes and homoeologous expression provide a novel way of balancing genetic incompatibilities, providing a new insight into the early stages of allopolyploidization in vertebrate evolution.

The Genetic Basis of Gene Expression Divergence in Antennae of Two Closely Related Moth Species, Helicoverpa armigera and Helicoverpa assulta

The closely related species Helicoverpa armigera (H. armigera) and Helicoverpa assulta (H. assulta) have different host plant ranges and share two principal components of sex pheromones but with reversed ratios. The antennae are the main olfactory organ of insects and play a crucial role in host plant selection and mate seeking. However, the genetic basis for gene expression divergence in the antennae of the two species is unclear. We performed an allele-specific expression (ASE) analysis in the antennal transcriptomes of the two species and their F₁ hybrids, examining the connection between gene expression divergence and phenotypic differences. The results show that the proportion of genes classified as all cis was higher than that of all trans in males and reversed in females. The contribution of regulatory patterns to gene expression divergence in males was less than that in females, which explained the functional differentiation of male and female antennae. Among the five groups of F₁ hybrids, the fertile males from the cross of H. armigera female and H. assulta male had the lowest proportion of misexpressed genes, and the inferred regulatory patterns were more accurate. By using this group of F₁ hybrids, we discovered that cis-related regulations play a crucial role in gene expression divergence of sex pheromone perception-related proteins. These results are helpful for understanding how specific changes in the gene expression of olfactory-related genes can contribute to rapid evolutionary changes in important olfactory traits in closely related moths.

Effects of homoeologous exchange on gene expression and alternative splicing in a newly formed allotetraploid wheat

Homoeologous exchange (HE) is a major mechanism generating post-polyploidization genetic variation with important evolutionary consequences. However, the direct impacts of HE on gene expression and transcript diversity in allopolyploids without the intertwined evolutionary processes remain to be fully understood. Here, we analyzed high-throughput RNA-seq data of young leaves from plant groups of a synthetic allotetraploid wheat (AADD), which contained variable numbers of HEs. We aimed to investigate if and to which extent HE directly impacts gene expression and alternative splicing (AS). We found that HE impacts expression of genes located within HE regions primarily via a cis-acting dosage effect, which led to significant changes in the total expression level of homoeologous gene pairs, especially for homoeologs whose original expression was biased. In parallel, HE also influences expression of a large number of genes residing in non-HE regions by trans-regulation leading to convergent expression of homoeologs. Intriguingly, when taking the original relative homoeolog expression states into account, homoeolog pairs under trans-effect are more prone to manifesting a convergent response to the HEs whereas those under cis-regulation tended to show further exacerbated subgenome-biased expression. Moreover, HE-induced quantitative, largely individual-specific, changes of AS events were detected. Similar to homoeologous expression, homoeo-AS events under trans-effect were more responsive to HE. HE therefore exerts multifaceted immediate effects on gene expression and, to a less extent, on individualized transcript diversity in nascent allopolyploidy.

Assessing Loss of Regulatory Divergence, Genome-Transcriptome Incongruence, and Preferential Expression Switching in Abaca × Banana Backcrosses

The Musa textilis var. Abuab has high fiber quality (FQ) but is susceptible to abaca bunchy top virus (AbBTV); the Musa balbisiana var. Pacol has low FQ but is resistant against AbBTV. Their backcrosses (BC2 and BC3) possess both desirable traits. Analysis using RNA-seq showed that the regulatory divergence of Abuab and Pacol is largely explained by cis differences with 27.4% and 22.3% if we are to assess it using BC2 and BC3, respectively. Cis differences between the two genotypes are significantly reduced from BC2 to BC3 due to changes in genomic constitution. Trans, on the other hand, is robust to changes in allelic composition. All these are attributed to the loss of heterozygosity in BC3 relative to BC2. Further analysis showed that both backcrosses exhibited genome-wide preferential expression of Pacol- over Abuab-specific alleles, despite the wider genetic presence of the latter in the hybrids. The ratio of the two genotype-specific expressed transcripts and the ratio of their corresponding genetic make-up are significantly disproportionate, a phenomenon that we refer to here as "genome-transcriptome incongruence". We also observed preferential expression switching in which several genes prefer the Abuab- (or Pacol-) specific allele in BC2 but switched to the Pacol- (or Abuab-) specific allele in the BC3 genome.

Proteotoxicity caused by perturbed protein complexes underlies hybrid incompatibility in yeast

Dobzhansky–Muller incompatibilities represent a major driver of reproductive isolation between species. They are caused when interacting components encoded by alleles from different species cannot function properly when mixed. At incipient stages of speciation, complex incompatibilities involving multiple genetic loci with weak effects are frequently observed, but the underlying mechanisms remain elusive. Here we show perturbed proteostasis leading to compromised mitosis and meiosis in Saccharomyces cerevisiae hybrid lines carrying one or two chromosomes from Saccharomyces bayanus var. uvarum. Levels of proteotoxicity are correlated with the number of protein complexes on replaced chromosomes. Proteomic approaches reveal that multi-protein complexes with subunits encoded by replaced chromosomes tend to be unstable. Furthermore, hybrid defects can be alleviated or aggravated, respectively, by up- or down-regulating the ubiquitin-proteasomal degradation machinery, suggesting that destabilized complex subunits overburden the proteostasis machinery and compromise hybrid fitness. Our findings reveal the general role of impaired protein complex assembly in complex incompatibilities.

Hybrid incompatibility can be an important element of reproductive isolation and speciation. Using chromosome replacement lines of yeast, the authors show that perturbed proteostasis caused by destabilized hybrid protein complexes may represent a general mechanism of hybrid incompatibility.

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.

Gene expression variation explains maize seed germination heterosis

BACKGROUND

Heterosis has been extensively utilized in plant breeding, however, the underlying molecular mechanism remains largely elusive. Maize (Zea mays), which exhibits strong heterosis, is an ideal material for studying heterosis.

RESULTS

In this study, there is faster imbibition and development in reciprocal crossing Zhengdan958 hybrids than in their parent lines during seed germination. To investigate the mechanism of heterosis of maize germination, comparative transcriptomic analyses were conducted. The gene expression patterns showed that 1324 (47.27%) and 1592 (66.44%) of the differential expression genes between hybrids and either parental line display parental dominance up or higher levels in the reciprocal cross of Zhengdan958, respectively. Notably, these genes were mainly enriched in metabolic pathways, including carbon metabolism, glycolysis/gluconeogenesis, protein processing in endoplasmic reticulum, etc. CONCLUSION: Our results provide evidence for the higher expression level genes in hybrid involved in metabolic pathways acting as main contributors to maize seed germinating heterosis. These findings provide new insights into the gene expression variation of maize embryos and improve the understanding of maize seed germination heterosis.

The dynamic effect of genetic variation on the in vivo ER stress transcriptional response in different tissues

The genetic regulation of gene expression varies greatly across tissue-type and individuals and can be strongly influenced by the environment. Many variants, under healthy control conditions, may be silent or even have the opposite effect under diseased stress conditions. This study uses an in vivo mouse model to investigate how the effect of genetic variation changes with cellular stress across different tissues. Endoplasmic reticulum stress occurs when misfolded proteins accumulate in the endoplasmic reticulum. This triggers the unfolded protein response, a large transcriptional response which attempts to restore homeostasis. This transcriptional response, despite being a conserved, basic cellular process, is highly variable across different genetic backgrounds, making it an ideal system to study the dynamic effects of genetic variation. In this study, we sought to better understand how genetic variation alters expression across tissues, in the presence and absence of endoplasmic reticulum stress. The use of different mouse strains and their F1s allow us to also identify context-specific cis- and trans- regulatory variation underlying variable transcriptional responses. We found hundreds of genes that respond to endoplasmic reticulum stress in a tissue- and/or genotype-dependent manner. The majority of the regulatory effects we identified were acting in cis-, which in turn, contribute to the variable endoplasmic reticulum stress- and tissue-specific transcriptional response. This study demonstrates the need for incorporating environmental stressors across multiple different tissues in future studies to better elucidate the effect of any particular genetic factor in basic biological pathways, like the endoplasmic reticulum stress response.

RNA-Seq Analysis Reveals Expression Regulatory Divergence of W-Linked Genes between Two Contrasting Chicken Breeds

Simple Summary

Understanding the mode of gene expression and regulation is essential for understanding the evolutionary process. Many previous studies tried to explain regulatory changes at the autosomal level, but little research has extended these explorations to the field of sex chromosomes due to their complex sex-limit features. Here, we first adopted an innovative method of identifying regulatory divergence of W-linked genes. Compared with cis-regulatory divergence, trans acting genes were more extensive in the W chromosome. We also found that divergent sex specific selection cannot strongly affect the expression evolution of the W chromosome. This insensitivity to selection may be one of the reasons why regulatory divergence is so small between autosomal and sex chromosomes.

Abstract

The regulation of gene expression is a complex process involving organism function and phenotypic diversity, and is caused by cis- and trans- regulation. While prior studies identified the regulatory pattern of the autosome rewiring in hybrids, the role of gene regulation in W sex chromosomes is not clear due to their degradation and sex-limit expression. Here, we developed reciprocal crosses of two chicken breeds, White Leghorn and Cornish Game, which exhibited broad differences in gender-related traits, and assessed the expression of the genes on the W chromosome to disentangle the contribution of cis- and trans-factors to expression divergence. We found that female-specific selection does not have a significant effect on W chromosome gene-expression patterns. For different tissues, there were most parental divergence expression genes in muscle, and also more heterosis compared with two other tissues. Notably, a broader pattern of trans regulation in the W chromosome was observed, which is consistent with autosomes. Taken together, this work describes the regulatory divergence of W-linked genes between two contrasting breeds and indicates sex chromosomes have a unique regulation and expression mechanism.

Insights Into Persistent HIV-1 Infection and Functional Cure: Novel Capabilities and Strategies

Although HIV-1 replication can be efficiently suppressed to undetectable levels in peripheral blood by combination antiretroviral therapy (cART), lifelong medication is still required in people living with HIV (PLWH). Life expectancies have been extended by cART, but age-related comorbidities have increased which are associated with heavy physiological and economic burdens on PLWH. The obstacle to a functional HIV cure can be ascribed to the formation of latent reservoir establishment at the time of acute infection that persists during cART. Recent studies suggest that some HIV reservoirs are established in the early acute stages of HIV infection within multiple immune cells that are gradually shaped by various host and viral mechanisms and may undergo clonal expansion. Early cART initiation has been shown to reduce the reservoir size in HIV-infected individuals. Memory CD4+ T cell subsets are regarded as the predominant cellular compartment of the HIV reservoir, but monocytes and derivative macrophages or dendritic cells also play a role in the persistent virus infection. HIV latency is regulated at multiple molecular levels in transcriptional and post-transcriptional processes. Epigenetic regulation of the proviral promoter can profoundly regulate the viral transcription. In addition, transcriptional elongation, RNA splicing, and nuclear export pathways are also involved in maintaining HIV latency. Although most proviruses contain large internal deletions, some defective proviruses may induce immune activation by expressing viral proteins or producing replication-defective viral-like particles. In this review article, we discuss the state of the art on mechanisms of virus persistence in the periphery and tissue and summarize interdisciplinary approaches toward a functional HIV cure, including novel capabilities and strategies to measure and eliminate the infected reservoirs and induce immune control.

Nonadditive gene expression is correlated with nonadditive phenotypic expression in interspecific triploid hybrids of willow (Salix spp.)

Many studies have highlighted the complex and diverse basis for heterosis in inbred crops. Despite the lack of a consensus model, it is vital that we turn our attention to understanding heterosis in undomesticated, heterozygous, and polyploid species, such as willow (Salix spp.). Shrub willow is a dedicated energy crop bred to be fast-growing and high yielding on marginal land without competing with food crops. A trend in willow breeding is the consistent pattern of heterosis in triploids produced from crosses between diploid and tetraploid species. Here, we test whether differentially expressed genes are associated with heterosis in triploid families derived from diploid Salix purpurea, diploid Salix viminalis, and tetraploid Salix miyabeana parents. Three biological replicates of shoot tips from all family progeny and parents were collected after 12 weeks in the greenhouse and RNA extracted for RNA-Seq analysis. This study provides evidence that nonadditive patterns of gene expression are correlated with nonadditive phenotypic expression in interspecific triploid hybrids of willow. Expression-level dominance was most correlated with heterosis for biomass yield traits and was highly enriched for processes involved in starch and sucrose metabolism. In addition, there was a global dosage effect of parent alleles in triploid hybrids, with expression proportional to copy number variation. Importantly, differentially expressed genes between family parents were most predictive of heterosis for both field and greenhouse collected traits. Altogether, these data will be used to progress models of heterosis to complement the growing genomic resources available for the improvement of heterozygous perennial bioenergy crops.

cis-Regulatory changes in locomotor genes are associated with the evolution of burrowing behavior

How evolution modifies complex, innate behaviors is largely unknown. Divergence in many morphological traits, and some behaviors, is linked to cis-regulatory changes in gene expression. Given this, we compare brain gene expression of two interfertile sister species of Peromyscus mice that show large and heritable differences in burrowing behavior. Species-level differential expression and allele-specific expression in F₁ hybrids indicate a preponderance of cis-regulatory divergence, including many genes whose cis-regulation is affected by burrowing behavior. Genes related to locomotor coordination show the strongest signals of lineage-specific selection on burrowing-induced cis-regulatory changes. Furthermore, genetic markers closest to these candidate genes associate with variation in burrow shape in a genetic cross, suggesting an enrichment for loci affecting burrowing behavior near these candidate locomotor genes. Our results provide insight into how cis-regulated gene expression can depend on behavioral context and how this dynamic regulatory divergence between species may contribute to behavioral evolution.

Gene age shapes the transcriptional landscape of sexual morphogenesis in mushroom forming fungi (Agaricomycetes)

Multicellularity has been one of the most important innovations in the history of life. The role of gene regulatory changes in driving transitions to multicellularity is being increasingly recognized; however, factors influencing gene expression patterns are poorly known in many clades. Here, we compared the developmental transcriptomes of complex multicellular fruiting bodies of eight Agaricomycetes and Cryptococcus neoformans, a closely related human pathogen with a simple morphology. In-depth analysis in Pleurotus ostreatus revealed that allele-specific expression, natural antisense transcripts, and developmental gene expression, but not RNA editing or a ‘developmental hourglass,’ act in concert to shape its transcriptome during fruiting body development. We found that transcriptional patterns of genes strongly depend on their evolutionary ages. Young genes showed more developmental and allele-specific expression variation, possibly because of weaker evolutionary constraint, suggestive of nonadaptive expression variance in fruiting bodies. These results prompted us to define a set of conserved genes specifically regulated only during complex morphogenesis by excluding young genes and accounting for deeply conserved ones shared with species showing simple sexual development. Analysis of the resulting gene set revealed evolutionary and functional associations with complex multicellularity, which allowed us to speculate they are involved in complex multicellular morphogenesis of mushroom fruiting bodies.

Transcriptomic Variations and Network Hubs Controlling Seed Size and Weight During Maize Seed Development

To elucidate the mechanisms underlying seed development in maize, comprehensive RNA-seq analyses were conducted on Zhengdan1002 (ZD1002), Zhengdan958 (ZD958), and their parental lines during seven seed developmental stages. We found that gene expression levels were largely nonadditive in hybrids and that cis-only or trans × cis pattern played a large role in hybrid gene regulation during seed developmental stage. Weighted gene co-expression network (WGCNA) analysis showed that 36 modules were highly correlated (r = -0.90-0.92, p < 0.05) with kernel weight, length, and width during seed development. Forty-five transcription factors and 38 ribosomal protein genes were identified as major hub genes determining seed size/weight. We also described a network hub, Auxin Response Factor 12 of maize (ZmARF12), a member of a family of transcription factor that mediate gene expression in response to auxin, potentially links auxin signal pathways, cell division, and the size of the seeds. The ZmARF12 mutant exhibited larger seed size and higher grain weight. ZmARF12 transcription was negatively associated with cell division during seed development, which was confirmed by evaluating the yield of protoplasts that isolated from the kernels of the mutant and other inbred lines. Transient knock-down of ZmARF12 in maize plants facilitated cell expansion and division, whereas transient silencing of its potential interactor ZmIAA8 impaired cell division. ZmIAA8 expression was repressed in the ZmARF12 over-expressed protoplasts. The mutant phenotype and the genetics studies presented here illustrated evidence that ZmARF12 is a cell division repressor, and potentially determines the final seed size.

The transcriptomic signature of obligate parthenogenesis

Investigating the origin of parthenogenesis through interspecific hybridization can provide insight into how meiosis may be altered by genetic incompatibilities, which is fundamental for our understanding of the formation of reproductive barriers. Yet the genetic mechanisms giving rise to obligate parthenogenesis in eukaryotes remain understudied. In the microcrustacean Daphnia pulex species complex, obligately parthenogenetic (OP) isolates emerged as backcrosses of two cyclically parthenogenetic (CP) parental species, D. pulex and D. pulicaria, two closely related but ecologically distinct species. We examine the genome-wide expression in OP females at the early resting egg production stage, a life-history stage distinguishing OP and CP reproductive strategies, in comparison to CP females of the same stage from the two parental species. Our analyses of the expression data reveal that underdominant and overdominant genes are abundant in OP isolates, suggesting widespread regulatory incompatibilities between the parental species. More importantly, underdominant genes (i.e., genes with expression lower than both parentals) in the OP isolates are enriched in meiosis and cell-cycle pathways, indicating an important role of underdominance in the origin of obligate parthenogenesis. Furthermore, metabolic and biosynthesis pathways enriched with overdominant genes (i.e., expression higher than both parentals) are another genomic signature of OP isolates.

Divergence of X-linked trans regulatory proteins and the misexpression of gene targets in sterile Drosophila pseudoobscura hybrids

BACKGROUND

The genetic basis of hybrid incompatibilities is characterized by pervasive cases of gene interactions. Sex chromosomes play a major role in speciation and X-linked hybrid male sterility (HMS) genes have been identified. Interestingly, some of these genes code for proteins with DNA binding domains, suggesting a capability to act as trans-regulatory elements and disturb the expression of a large number of gene targets. To understand how interactions between trans- and cis-regulatory elements contribute to speciation, we aimed to map putative X-linked trans-regulatory elements and to identify gene targets with disrupted gene expression in sterile hybrids between the subspecies Drosophila pseudoobscura pseudoobscura and D. p. bogotana.

RESULTS

We find six putative trans-regulatory proteins within previously mapped X chromosome HMS loci with sequence changes that differentiate the two subspecies. Among them, the previously characterized HMS gene Overdrive (Ovd) had the largest number of amino acid changes between subspecies, with some substitutions localized within the protein's DNA binding domain. Using an introgression approach, we detected transcriptional responses associated with a sterility/fertility Ovd allele swap. We found a network of 52 targets of Ovd and identified cis-regulatory effects among target genes with disrupted expression in sterile hybrids. However, a combined analysis of polymorphism and divergence in non-coding sequences immediately upstream of target genes found no evidence of changes in candidate regulatory proximal cis-elements. Finally, peptidases were over-represented among target genes.

CONCLUSIONS

We provide evidence of divergence between subspecies within the DNA binding domain of the HMS protein Ovd and identify trans effects on the expression of 52 gene targets. Our results identify a network of trans-cis interactions with possible effects on HMS. This network provides molecular evidence of gene × gene incompatibilities as contributors to hybrid dysfunction.

Population biology of accessory gland-expressed de novo genes in Drosophila melanogaster

Early work on de novo gene discovery in Drosophila was consistent with the idea that many such genes have male-biased patterns of expression, including a large number expressed in the testis. However, there has been little formal analysis of variation in the abundance and properties of de novo genes expressed in different tissues. Here, we investigate the population biology of recently evolved de novo genes expressed in the Drosophila melanogaster accessory gland, a somatic male tissue that plays an important role in male and female fertility and the post mating response of females, using the same collection of inbred lines used previously to identify testis-expressed de novo genes, thus allowing for direct cross tissue comparisons of these genes in two tissues of male reproduction. Using RNA-seq data, we identify candidate de novo genes located in annotated intergenic and intronic sequence and determine the properties of these genes including chromosomal location, expression, abundance, and coding capacity. Generally, we find major differences between the tissues in terms of gene abundance and expression, though other properties such as transcript length and chromosomal distribution are more similar. We also explore differences between regulatory mechanisms of de novo genes in the two tissues and how such differences may interact with selection to produce differences in D. melanogaster de novo genes expressed in the two tissues.

Genome-wide association mapping of transcriptome variation in Mimulus guttatus indicates differing patterns of selection on cis- versus trans-acting mutations

We measured the floral bud transcriptome of 151 fully sequenced lines of Mimulus guttatus from one natural population. Thousands of single nucleotide polymorphisms (SNPs) are implicated as transcription regulators, but there is a striking difference in the allele frequency spectrum of cis-acting and trans-acting mutations. Cis-SNPs have intermediate frequencies (consistent with balancing selection) while trans-SNPs exhibit a rare-alleles model (consistent with purifying selection). This pattern only becomes clear when transcript variation is normalized on a gene-to-gene basis. If a global normalization is applied, as is typically in RNAseq experiments, asymmetric transcript distributions combined with "rarity disequilibrium" produce a superabundance of false positives for trans-acting SNPs. To explore the cause of purifying selection on trans-acting mutations, we identified gene expression modules as sets of coexpressed genes. The extent to which trans-acting mutations influence modules is a strong predictor of allele frequency. Mutations altering expression of genes with high "connectedness" (those that are highly predictive of the representative module expression value) have the lowest allele frequency. The expression modules can also predict whole-plant traits such as flower size. We find that a substantial portion of the genetic (co)variance among traits can be described as an emergent property of genetic effects on expression modules.

Spatial and Temporal Transcriptomic Heredity and Asymmetry in an Artificially Constructed Allotetraploid Wheat (AADD)

Polyploidy, or whole-genome duplication (WGD), often induces dramatic changes in gene expression due to "transcriptome shock. " However, questions remain about how allopolyploidy (the merging of multiple nuclear genomes in the same nucleus) affects gene expression within and across multiple tissues and developmental stages during the initial foundation of allopolyploid plants. Here, we systematically investigated the immediate effect of allopolyploidy on gene expression variation in an artificial allopolyploidy system consisting of a constructed allotetraploid wheat (AADD genome, accession AT2) and its diploid progenitors Triticum urartu and Aegilops tauschii. We performed comprehensive RNA sequencing of 81 samples from different genotypes, tissues, and developmental stages. First, we found that intrinsic interspecific differences between the diploid parents played a major role in establishing the expression architecture of the allopolyploid. Nonetheless, allopolyploidy per se also induced dramatic and asymmetric patterns of differential gene expression between the subgenomes, and genes from the D subgenome exhibited a more drastic response. Second, analysis of homoeolog expression bias (HEB) revealed that the D subgenome exhibited significant expression bias and that de novo-generated HEB was attributed mainly to asymmetrical differential gene expression. Homoeolog-specific expression (HSE) analyses showed that the cis-only regulatory pattern was predominant in AT2, reflecting significant divergence between the parents. Co-expression network analysis revealed that homoeolog expression connectivity (HEC) was significantly correlated with sequence divergence in cis elements between subgenomes. Interestingly, allopolyploidy-induced reconstruction of network modules was also associated with different HSE patterns. Finally, a transcriptome atlas of spike development demonstrated that the phenotypic similarity of AT2 to T. urartu may be attributed to the combination of relatively stable expression of A-subgenome genes and drastic downregulation of their D-subgenome homoeologs. These findings provide a broad, multidimensional characterization of allopolyploidy-induced transcriptomic responses and suggest that allopolyploidy can have immediate and complex regulatory effects on the expression of nuclear genes.

OUP accepted manuscript

Interspecific hybridization is often seen as a genomic stress that may lead to new gene expression patterns and deregulation of transposable elements (TEs). The understanding of expression changes in hybrids compared with parental species is essential to disentangle their putative role in speciation processes. However, to date we ignore the detailed mechanisms involved in genomic deregulation in hybrids. We studied the ovarian transcriptome and epigenome of the Drosophila buzzatii and Drosophila koepferae species together with their F₁ hybrid females. We found a trend toward underexpression of genes and TE families in hybrids. The epigenome in hybrids was highly similar to the parental epigenomes and showed intermediate histone enrichments between parental species in most cases. Differential gene expression in hybrids was often associated only with changes in H3K4me3 enrichments, whereas differential TE family expression in hybrids may be associated with changes in H3K4me3, H3K9me3, or H3K27me3 enrichments. We identified specific genes and TE families, which their differential expression in comparison with the parental species was explained by their differential chromatin mark combination enrichment. Finally, cis–trans compensatory regulation could also contribute in some way to the hybrid deregulation. This work provides the first study of histone content in Drosophila interspecific hybrids and their effect on gene and TE expression deregulation.