Sex-specific embryonic gene expression in species with newly evolved sex chromosomes

Function of doublesex and transformer-2 genes and its respond to environment factors in Dendroctonus armandi Tsai et Li (Coleoptera: Curculionidae: Scolytinae)

The doublesex (dsx) and transformer-2 (tra-2) genes play important roles in insect sex determination. However, the sex determination of Dendroctonus armandi in the native and coniferous forests of the Qinling Mountains is still unclear. In this study, we cloned and identified the full-length sequences of two Dadsx isoform (Dadsx1 and Dadsx2) and two isoform Datra-2 (Datra-2A and Datra-2B) in D. armandi. These four isoforms are expressed in larvae, pupae, as well as in males and females, and their expression levels are highest during the pupae stage. The relative expression levels of Dadsx and Datra-2 genes showed significant differences under different environment treatments (including temperature, nutrient, feeding duration, and terpenoid), with male Dadsx1 relative expression levels significantly higher than females and female Dadsx2 relative expression levels significantly higher than males. The silencing of Datra-2A and Datra-2B isoforms both leads to the downregulation of Dadsx1 isoform. Compared to the control group (average mortality rate: 36.7 %; average deformity rate: 0.2 %; emergence rate: 60.9 %; sex ratio: 1.11), when Dadsx1 (average mortality rate: 69.0 %; average deformity rate: 16.5 %; emergence rate: 29.6 %; sex ratio: 0.28), Dadsx2 (average mortality rate: 64.4 %; average deformity rate: 17.4 %; emergence rate: 25.4 %; sex ratio: 2.41), Datra-2A (average mortality rate: 65.1 %; average deformity rate: 16.1 %; emergence rate: 15.4 %; sex ratio: 0.31) and Datra-2B (average mortality rate: 65.4 %; average deformity rate: 17.2 %; emergence rate: 13.8 %; sex ratio: 0.33) isoform s were silenced, it showed a significant increase in mortality and deformity rates, a significant decrease in emergence rate, and a severe sex imbalance. The results indicate that the relative expression levels of Dadsx and Datra-2 genes are influenced by external factors and play a crucial role in maintaining the sex ratio of D. armandi and ensuring its lifecycle development.

The role of unbalanced segmental duplication in sex chromosome evolution in Australian ridge-tailed goannas

Varanids are known for conserved sex chromosomes, but there are differences in the size of the W chromosome but not in morphology among species representing varying stages of sex chromosome evolution. We tested for homology of the ZW sex chromosome system with size differences in varanids among four species from two lineages in Australia, the Odatria and the Gouldii. We found that while DNA sequences of the sex chromosomes are conserved in the species we tested, we also identified a homologous region on an enlarged autosomal microchromosome that shares sequences with the W chromosome in some isolated populations of V. acanthurus and V. citrinus from the Odatria lineage. The enlarged microchromosome was unpaired in all individuals tested and is likely an unbalanced segmental duplication translocated between chromosome 1, the W, and another microchromosome. This suggests an ancient balanced duplication homologous to the W and the terminal region of the long arm of chromosome 1. The most parsimonious explanation is that the duplicated region likely originated on chromosome 1. We hypothesised in our reconstruction that genes and related DNA sequences associated with the sex-linkage group have likely originated on an autosome. Subsequently, the sequences may have undergone duplication and translocation to the W chromosome, followed by the accumulation of lineage specific repeat elements and amplifications on the W at different rates in various lineages. Lastly, these sequences are likely to have undergone duplication and translocation to another autosomal microchromosome. Given the role of segmental duplications and translocations as important evolutionary drivers of speciation in other taxa, together with the rapid speciation that has occurred in Australian varanids, our findings provide broader insight into the evolutionary pathway leading to rapid chromosomal and genic divergence of species.

Genomic Resources for the Scuttle Fly Megaselia abdita: A Model Organism for Comparative Developmental Studies in Flies

The order Diptera (true flies) holds promise as a model taxon in evolutionary developmental biology due to the inclusion of the model organism, Drosophila melanogaster, and the ability to cost-effectively rear many species in laboratories. One of them, the scuttle fly Megaselia abdita (Phoridae) has been used in evolutionary developmental biology for 30 years and is an excellent phylogenetic intermediate between fruit flies and mosquitoes but remains underdeveloped in genomic resources. Here, we present a de novo chromosome-level assembly and annotation of M. abdita and transcriptomes of 9 embryonic and 4 postembryonic stages. We also compare 9 stage-matched embryonic transcriptomes between M. abdita and D. melanogaster. Our analysis of these resources reveals extensive chromosomal synteny with D. melanogaster, 28 orphan genes with embryo-specific expression including a novel F-box LRR gene in M. abdita, and conserved and diverged features of gene expression dynamics between M. abdita and D. melanogaster. Collectively, our results provide a new reference for studying the diversification of developmental processes in flies.

Dosage compensation in non-model insects - progress and perspectives

In many multicellular eukaryotes, heteromorphic sex chromosomes are responsible for determining the sexual characteristics and reproductive functions of individuals. Sex chromosomes can cause a dosage imbalance between sexes, which in some species is re-equilibrated by dosage compensation (DC). Recent genomic advances have extended our understanding of DC mechanisms in insects beyond model organisms such as Drosophila melanogaster. We review current knowledge of insect DC, focusing on its conservation and divergence across orders, the evolutionary dynamics of neo-sex chromosomes, and the diversity of molecular mechanisms. We propose a framework to uncover DC regulators in non-model insects that relies on integrating evolutionary, genomic, and functional approaches. This comprehensive approach will facilitate a deeper understanding of the evolution and essentiality of gene regulatory mechanisms.

Evidence for a Novel X Chromosome in Termites

Termites, together with cockroaches, belong to the Blattodea. They possess an XX/XY sex determination system which has evolved from an XX/X0 system present in other Blattodean species, such as cockroaches and wood roaches. Little is currently known about the sex chromosomes of termites, their gene content, or their evolution. We here investigate the X chromosome of multiple termite species and compare them with the X chromosome of cockroaches using genomic and transcriptomic data. We find that the X chromosome of the termite Macrotermes natalensis is large and differentiated showing hall marks of sex chromosome evolution such as dosage compensation, while this does not seem to be the case in the other two termite species investigated here where sex chromosomes may be evolutionary younger. Furthermore, the X chromosome in M. natalensis is different from the X chromosome found in the cockroach Blattella germanica indicating that sex chromosome turn-over events may have happened during termite evolution.

Single-cell RNA-seq of Drosophila miranda testis reveals the evolution and trajectory of germline sex chromosome regulation

Although sex chromosomes have evolved from autosomes, they often have unusual regulatory regimes that are sex- and cell-type-specific such as dosage compensation (DC) and meiotic sex chromosome inactivation (MSCI). The molecular mechanisms and evolutionary forces driving these unique transcriptional programs are critical for genome evolution but have been, in the case of MSCI in Drosophila, subject to continuous debate. Here, we take advantage of the younger sex chromosomes in D. miranda (XR and the neo-X) to infer how former autosomes acquire sex-chromosome-specific regulatory programs using single-cell and bulk RNA sequencing and ribosome profiling, in a comparative evolutionary context. We show that contrary to mammals and worms, the X down-regulation through germline progression is most consistent with the shutdown of DC instead of MSCI, resulting in half gene dosage at the end of meiosis for all 3 X's. Moreover, lowly expressed germline and meiotic genes on the neo-X are ancestrally lowly expressed, instead of acquired suppression after sex linkage. For the young neo-X, DC is incomplete across all tissue and cell types and this dosage imbalance is rescued by contributions from Y-linked gametologs which produce transcripts that are translated to compensate both gene and protein dosage. We find an excess of previously autosomal testis genes becoming Y-specific, showing that the neo-Y and its masculinization likely resolve sexual antagonism. Multicopy neo-sex genes are predominantly expressed during meiotic stages of spermatogenesis, consistent with their amplification being driven to interfere with mendelian segregation. Altogether, this study reveals germline regulation of evolving sex chromosomes and elucidates the consequences these unique regulatory mechanisms have on the evolution of sex chromosome architecture.

Sex-specific splicing occurs genome-wide during early Drosophila embryogenesis

Sex-specific splicing is an essential process that regulates sex determination and drives sexual dimorphism. Yet, how early in development widespread sex-specific transcript diversity occurs was unknown because it had yet to be studied at the genome-wide level. We use the powerful Drosophila model to show that widespread sex-specific transcript diversity occurs early in development, concurrent with zygotic genome activation. We also present a new pipeline called time2Splice to quantify changes in alternative splicing over time. Furthermore, we determine that one of the consequences of losing an essential maternally deposited pioneer factor called CLAMP (chromatin-linked adapter for MSL proteins) is altered sex-specific splicing of genes involved in diverse biological processes that drive development. Overall, we show that sex-specific differences in transcript diversity exist even at the earliest stages of development..

Continuous transcriptome analysis reveals novel patterns of early gene expression in Drosophila embryos

The transformative events during early organismal development lay the foundation for body formation and long-term phenotype. The rapid progression of events and the limited material available present major barriers to studying these earliest stages of development. Herein, we report an operationally simple RNA sequencing approach for high-resolution, time-sensitive transcriptome analysis in early (≤3 h) Drosophila embryos. This method does not require embryo staging but relies on single-embryo RNA sequencing and transcriptome ordering along a developmental trajectory (pseudo-time). The resulting high-resolution, time-sensitive mRNA expression profiles reveal the exact onset of transcription and degradation for thousands of transcripts. Further, using sex-specific transcription signatures, embryos can be sexed directly, eliminating the need for Y chromosome genotyping and revealing patterns of sex-biased transcription from the beginning of zygotic transcription. Our data provide an unparalleled resolution of gene expression during early development and enhance the current understanding of early transcriptional processes.

The early embryonic transcriptome of a Hawaiian Drosophila picture-wing fly shows evidence of altered gene expression and novel gene evolution

A massive adaptive radiation on the Hawaiian archipelago has produced approximately one-quarter of the fly species in the family Drosophilidae. The Hawaiian Drosophila clade has long been recognized as a model system for the study of both the ecology of island endemics and the evolution of developmental mechanisms, but relatively few genomic and transcriptomic datasets are available for this group. We present here a differential expression analysis of the transcriptional profiles of two highly conserved embryonic stages in the Hawaiian picture-wing fly Drosophila grimshawi. When we compared our results to previously published datasets across the family Drosophilidae, we identified cases of both gains and losses of gene representation in D. grimshawi, including an apparent delay in Hox gene activation. We also found a high expression of unannotated genes. Most transcripts of unannotated genes with open reading frames do not have identified homologs in non-Hawaiian Drosophila species, although the vast majority have sequence matches in genomes of other Hawaiian picture-wing flies. Some of these unannotated genes may have arisen from noncoding sequence in the ancestor of Hawaiian flies or during the evolution of the clade. Our results suggest that both the modified use of ancestral genes and the evolution of new ones may occur in rapid radiations.

Establishment of H3K9me3-dependent heterochromatin during embryogenesis in Drosophila miranda

Heterochromatin is a key architectural feature of eukaryotic genomes crucial for silencing of repetitive elements. During Drosophila embryonic cellularization, heterochromatin rapidly appears over repetitive sequences, but the molecular details of how heterochromatin is established are poorly understood. Here, we map the genome-wide distribution of H3K9me3-dependent heterochromatin in individual embryos of Drosophila miranda at precisely staged developmental time points. We find that canonical H3K9me3 enrichment is established prior to cellularization and matures into stable and broad heterochromatin domains through development. Intriguingly, initial nucleation sites of H3K9me3 enrichment appear as early as embryonic stage 3 over transposable elements (TEs) and progressively broaden, consistent with spreading to neighboring nucleosomes. The earliest nucleation sites are limited to specific regions of a small number of recently active retrotransposon families and often appear over promoter and 5' regions of LTR retrotransposons, while late nucleation sites develop broadly across the entirety of most TEs. Interestingly, early nucleating TEs are strongly associated with abundant maternal piRNAs and show early zygotic transcription. These results support a model of piRNA-associated co-transcriptional silencing while also suggesting additional mechanisms for site-restricted H3K9me3 nucleation at TEs in pre-cellular Drosophila embryos.

Epigenetic conflict on a degenerating Y chromosome increases mutational burden in Drosophila males

Large portions of eukaryotic genomes consist of transposable elements (TEs), and the establishment of transcription-repressing heterochromatin during early development safeguards genome integrity in Drosophila. Repeat-rich Y chromosomes can act as reservoirs for TEs ('toxic' Y effect), and incomplete epigenomic defenses during early development can lead to deleterious TE mobilization. Here, we contrast the dynamics of early TE activation in two Drosophila species with vastly different Y chromosomes of different ages. Zygotic TE expression is elevated in male embryos relative to females in both species, mostly due to expression of Y-linked TEs. Interestingly, male-biased TE expression diminishes across development in D. pseudoobscura, but remains elevated in D. miranda, the species with the younger and larger Y chromosome. The repeat-rich Y of D. miranda still contains many actively transcribed genes, which compromise the formation of silencing heterochromatin. Elevated TE expression results in more de novo insertions of repeats in males compared to females. This lends support to the idea that the 'toxic' Y chromosome can create a mutational burden in males when genome-wide defense mechanisms are compromised, and suggests a previously unappreciated epigenetic conflict on evolving Y chromosomes between transcription of essential genes and silencing of selfish DNA.

Evolved Differences in cis and trans Regulation Between the Maternal and Zygotic mRNA Complements in the Drosophila Embryo

How gene expression can evolve depends on the mechanisms driving gene expression. Gene expression is controlled in different ways in different developmental stages; here we ask whether different developmental stages show different patterns of regulatory evolution. To explore the mode of regulatory evolution, we used the early stages of embryonic development controlled by two different genomes, that of the mother and that of the zygote. During embryogenesis in all animals, initial developmental processes are driven entirely by maternally provided gene products deposited into the oocyte. The zygotic genome is activated later, when developmental control is handed off from maternal gene products to the zygote during the maternal-to-zygotic transition. Using hybrid crosses between sister species of Drosophila (Dsimulans, D. sechellia, and D. mauritiana) and transcriptomics, we find that the regulation of maternal transcript deposition and zygotic transcription evolve through different mechanisms. We find that patterns of transcript level inheritance in hybrids, relative to parental species, differ between maternal and zygotic transcripts, and maternal transcript levels are more likely to be conserved. Changes in transcript levels occur predominantly through differences in trans regulation for maternal genes, while changes in zygotic transcription occur through a combination of both cis and trans regulatory changes. Differences in the underlying regulatory landscape in the mother and the zygote are likely the primary determinants for how maternal and zygotic transcripts evolve.

The conserved regulatory basis of mRNA contributions to the early Drosophila embryo differs between the maternal and zygotic genomes

The gene products that drive early development are critical for setting up developmental trajectories in all animals. The earliest stages of development are fueled by maternally provided mRNAs until the zygote can take over transcription of its own genome. In early development, both maternally deposited and zygotically transcribed gene products have been well characterized in model systems. Previously, we demonstrated that across the genus Drosophila, maternal and zygotic mRNAs are largely conserved but also showed a surprising amount of change across species, with more differences evolving at the zygotic stage than the maternal stage. In this study, we use comparative methods to elucidate the regulatory mechanisms underlying maternal deposition and zygotic transcription across species. Through motif analysis, we discovered considerable conservation of regulatory mechanisms associated with maternal transcription, as compared to zygotic transcription. We also found that the regulatory mechanisms active in the maternal and zygotic genomes are quite different. For maternally deposited genes, we uncovered many signals that are consistent with transcriptional regulation at the level of chromatin state through factors enriched in the ovary, rather than precisely controlled gene-specific factors. For genes expressed only by the zygotic genome, we found evidence for previously identified regulators such as Zelda and GAGA-factor, with multiple analyses pointing toward gene-specific regulation. The observed mechanisms of regulation are consistent with what is known about regulation in these two genomes: during oogenesis, the maternal genome is optimized to quickly produce a large volume of transcripts to provide to the oocyte; after zygotic genome activation, mechanisms are employed to activate transcription of specific genes in a spatiotemporally precise manner. Thus the genetic architecture of the maternal and zygotic genomes, and the specific requirements for the transcripts present at each stage of embryogenesis, determine the regulatory mechanisms responsible for transcripts present at these stages.

The Origin and Evolution of Maternal Genes

Proteins and RNA molecules are deposited into the developing egg by the mother. These gene products will drive the first stages of development and are coded by maternal genes. Maternal genes are essential, yet, despite their importance, their evolutionary dynamics is largely unknown. Here I review the current knowledge of maternal gene evolution. The evolutionary origin of maternal genes tends to be more recent than that of zygotic genes. Some studies support the theoretical prediction that maternal genes evolve faster than zygotic genes. However, most studies were done on a limited set of species and genes. I also discuss the way forward to understand the evolution of maternal genes by combining high-throughput genomics and theoretical evolutionary approaches.

Evolution of maternal and zygotic mRNA complements in the early Drosophila embryo

The earliest stages of animal development are controlled by maternally deposited mRNA transcripts and proteins. Once the zygote is able to transcribe its own genome, maternal transcripts are degraded, in a tightly regulated process known as the maternal to zygotic transition (MZT). While this process has been well-studied within model species, we have little knowledge of how the pools of maternal and zygotic transcripts evolve. To characterize the evolutionary dynamics and functional constraints on early embryonic expression, we created a transcriptomic dataset for 14 Drosophila species spanning over 50 million years of evolution, at developmental stages before and after the MZT, and compared our results with a previously published Aedes aegypti developmental time course. We found deep conservation over 250 million years of a core set of genes transcribed only by the zygote. This select group is highly enriched in transcription factors that play critical roles in early development. However, we also identify a surprisingly high level of change in the transcripts represented at both stages over the phylogeny. While mRNA levels of genes with maternally deposited transcripts are more highly conserved than zygotic genes, those maternal transcripts that are completely degraded at the MZT vary dramatically between species. We also show that hundreds of genes have different isoform usage between the maternal and zygotic genomes. Our work suggests that maternal transcript deposition and early zygotic transcription are remarkably dynamic over evolutionary time, despite the widespread conservation of early developmental processes.

Effect of Larval Nutrition on Maternal mRNA Contribution to the Drosophila Egg

Embryonic development begins under the control of maternal gene products, mRNAs and proteins that the mother deposits into the egg; the zygotic genome is activated some time later. Maternal control of early development is conserved across metazoans. Gene products contributed by mothers are critical to many early developmental processes, and set up trajectories for the rest of development. Maternal deposition of these factors is an often-overlooked aspect of parental investment. If the mother experiences challenging environmental conditions, such as poor nutrition, previous studies in Drosophila melanogaster have demonstrated a plastic response wherein these mothers may produce larger eggs to buffer the offspring against the same difficult environment. This additional investment can produce offspring that are more fit in the challenging environment. With this study, we ask whether D. melanogaster mothers who experience poor nutrition during their own development change their gene product contribution to the egg. We perform mRNA-Seq on eggs at a stage where all mRNAs are maternally derived, from mothers with different degrees of nutritional limitation. We find that nutritional limitation produces similar transcript changes at all degrees of limitation tested. Genes that have lower transcript abundance in nutritionally limited mothers are those involved in translation, which is likely one of the most energetically costly processes occurring in the early embryo. We find an increase in transcripts for transport and localization of macromolecules, and for the electron transport chain. The eggs produced by nutrition-limited mothers show a plastic response in mRNA deposition, which may better prepare the future embryo for development in a nutrition-limited environment.

Dosage Compensation in Drosophila-a Model for the Coordinate Regulation of Transcription

The sex chromosomes have special significance in the history of genetics. The chromosomal basis of inheritance was firmly established when Calvin Bridges demonstrated that exceptions to Mendel's laws of segregation were accompanied at the cytological level by exceptional sex chromosome segregation. The morphological differences between X and Y exploited in Bridges' experiments arose as a consequence of the evolution of the sex chromosomes. Originally a homologous chromosome pair, the degeneration of the Y chromosome has been accompanied by a requirement for increased expression of the single X chromosome in males. Drosophila has been a model for the study of this dosage compensation and has brought key strengths, including classical genetics, the exceptional cytology of polytene chromosomes, and more recently, comprehensive genomics. The impact of these studies goes beyond sex chromosome regulation, providing valuable insights into mechanisms for the establishment and maintenance of chromatin domains, and for the coordinate regulation of transcription.

Alternative Splicing within and between Drosophila Species, Sexes, Tissues, and Developmental Stages

Alternative pre-mRNA splicing ("AS") greatly expands proteome diversity, but little is known about the evolutionary landscape of AS in Drosophila and how it differs between embryonic and adult stages or males and females. Here we study the transcriptomes from several tissues and developmental stages in males and females from four species across the Drosophila genus. We find that 20-37% of multi-exon genes are alternatively spliced. While males generally express a larger number of genes, AS is more prevalent in females, suggesting that the sexes adopt different expression strategies for their specialized function. While the number of total genes expressed increases during early embryonic development, the proportion of expressed genes that are alternatively spliced is highest in the very early embryo, before the onset of zygotic transcription. This indicates that females deposit a diversity of isoforms into the egg, consistent with abundant AS found in ovary. Cluster analysis by gene expression ("GE") levels shows mostly stage-specific clustering in embryonic samples, and tissue-specific clustering in adult tissues. Clustering embryonic stages and adult tissues based on AS profiles results in stronger species-specific clustering, suggesting that diversification of splicing contributes to lineage-specific evolution in Drosophila. Most sex-biased AS found in flies is due to AS in gonads, with little sex-specific splicing in somatic tissues.

High Throughput Sequencing Identifies Misregulated Genes in the Drosophila Polypyrimidine Tract-Binding Protein (hephaestus) Mutant Defective in Spermatogenesis

The Drosophila polypyrimidine tract-binding protein (dmPTB or hephaestus) plays an important role during spermatogenesis. The heph² mutation in this gene results in a specific defect in spermatogenesis, causing aberrant spermatid individualization and male sterility. However, the array of molecular defects in the mutant remains uncharacterized. Using an unbiased high throughput sequencing approach, we have identified transcripts that are misregulated in this mutant. Aberrant transcripts show altered expression levels, exon skipping, and alternative 5’ ends. We independently verified these findings by reverse-transcription and polymerase chain reaction (RT-PCR) analysis. Our analysis shows misregulation of transcripts that have been connected to spermatogenesis, including components of the actomyosin cytoskeletal apparatus. We show, for example, that the Myosin light chain 1 (Mlc1) transcript is aberrantly spliced. Furthermore, bioinformatics analysis reveals that Mlc1 contains a high affinity binding site(s) for dmPTB and that the site is conserved in many Drosophila species. We discuss that Mlc1 and other components of the actomyosin cytoskeletal apparatus offer important molecular links between the loss of dmPTB function and the observed developmental defect in spermatogenesis. This study provides the first comprehensive list of genes misregulated in vivo in the heph² mutant in Drosophila and offers insight into the role of dmPTB during spermatogenesis.

Primary Sex Determination in Drosophila melanogaster Does Not Rely on the Male-Specific Lethal Complex

It has been proposed that the Male Specific Lethal (MSL) complex is active in Drosophila melanogaster embryos of both sexes prior to the maternal-to-zygotic transition. Elevated gene expression from the two X chromosomes of female embryos is proposed to facilitate the stable establishment of Sex-lethal (Sxl) expression, which determines sex and represses further activity of the MSL complex, leaving it active only in males. Important supporting data included female-lethal genetic interactions between the seven msl genes and either Sxl or scute and sisterlessA, two of the X-signal elements (XSE) that regulate early Sxl expression. Here I report contrary findings that there are no female-lethal genetic interactions between the msl genes and Sxl or its XSE regulators. Fly stocks containing the msl3(1) allele were found to exhibit a maternal-effect interaction with Sxl, scute, and sisterlessA mutations, but genetic complementation experiments showed that msl3 is neither necessary nor sufficient for the female-lethal interactions, which appear to be due to an unidentified maternal regulator of Sxl. Published data cited as evidence for an early function of the MSL complex in females, including a maternal effect of msl2, have been reevaluated and found not to support a maternal, or other effect, of the MSL complex in sex determination. These findings suggest that the MSL complex is not involved in primary sex determination or in X chromosome dosage compensation prior to the maternal-to-zygotic transition.

Sex Bias and Maternal Contribution to Gene Expression Divergence in Drosophila Blastoderm Embryos

Early embryogenesis is a unique developmental stage where genetic control of development is handed off from mother to zygote. Yet the contribution of this transition to the evolution of gene expression is poorly understood. Here we study two aspects of gene expression specific to early embryogenesis in Drosophila: sex-biased gene expression prior to the onset of canonical X chromosomal dosage compensation, and the contribution of maternally supplied mRNAs. We sequenced mRNAs from individual unfertilized eggs and precisely staged and sexed blastoderm embryos, and compared levels between D. melanogaster, D. yakuba, D. pseudoobscura and D. virilis. First, we find that mRNA content is highly conserved for a given stage and that studies relying on pooled embryos likely systematically overstate the degree of gene expression divergence. Unlike studies done on larvae and adults where most species show a larger proportion of genes with male-biased expression, we find that transcripts in Drosophila embryos are largely female-biased in all species, likely due to incomplete dosage compensation prior to the activation of the canonical dosage compensation mechanism. The divergence of sex-biased gene expression across species is observed to be often due to lineage-specific decrease of expression; the most drastic example of which is the overall reduction of male expression from the neo-X chromosome in D. pseudoobscura, leading to a pervasive female-bias on this chromosome. We see no evidence for a faster evolution of expression on the X chromosome in embryos (no “faster-X” effect), unlike in adults, and contrary to a previous study on pooled non-sexed embryos. Finally, we find that most genes are conserved in regard to their maternal or zygotic origin of transcription, and present evidence that differences in maternal contribution to the blastoderm transcript pool may be due to species-specific divergence of transcript degradation rates.

The early embryonic stage of development is unique as this stage begins under the genetic control of the mother and transitions to genetic control of the zygote. In order to investigate how these developmental specificities have shaped gene expression evolution, we studied the contribution of maternal gene expression to the pool of transcripts in the early embryo and how gene expression differs in early female and male embryos. We find extensive female-bias in gene expression, something not observed in adults. This is largely due to the delay in the establishment of a mechanism that balances gene expression from the differing number of X chromosomes in females and males. The maternal contribution to the transcript pool in the embryo is largely conserved over species, with some exceptions. However, there may also be some species-specific differences in the degradation rate of maternal transcripts, an often-overlooked contributor to variation in transcript level.