Effects of Gene Dose, Chromatin, and Network Topology on Expression in Drosophila melanogaster

Change in RhoGAP and RhoGEF availability drives transitions in cortical patterning and excitability in Drosophila

Actin cortex patterning and dynamics are critical for cell shape changes. These dynamics undergo transitions during development, often accompanying changes in collective cell behavior. Although mechanisms have been established for individual cells' dynamic behaviors, the mechanisms and specific molecules that result in developmental transitions in vivo are still poorly understood. Here, we took advantage of two developmental systems in Drosophila melanogaster to identify conditions that altered cortical patterning and dynamics. We identified a Rho guanine nucleotide exchange factor (RhoGEF) and Rho GTPase activating protein (RhoGAP) pair required for actomyosin waves in egg chambers. Specifically, depletion of the RhoGEF, Ect2, or the RhoGAP, RhoGAP15B, disrupted actomyosin wave induction, and both proteins relocalized from the nucleus to the cortex preceding wave formation. Furthermore, we found that overexpression of a different RhoGEF and RhoGAP pair, RhoGEF2 and Cumberland GAP (C-GAP), resulted in actomyosin waves in the early embryo, during which RhoA activation precedes actomyosin assembly by ∼4 s. We found that C-GAP was recruited to actomyosin waves, and disrupting F-actin polymerization altered the spatial organization of both RhoA signaling and the cytoskeleton in waves. In addition, disrupting F-actin dynamics increased wave period and width, consistent with a possible role for F-actin in promoting delayed negative feedback. Overall, we showed a mechanism involved in inducing actomyosin waves that is essential for oocyte development and is general to other cell types, such as epithelial and syncytial cells.

Aneuploidy-induced cellular behaviors: Insights from Drosophila

A balanced gene complement is crucial for proper cell function. Aneuploidy, the condition of having an imbalanced chromosome set, alters the stoichiometry of gene copy numbers and protein complexes and has dramatic consequences at the cellular and organismal levels. In humans, aneuploidy is associated with different pathological conditions including cancer, microcephaly, mental retardation, miscarriages, and aging. Over the last century, Drosophila has provided a valuable system for studying the consequences of systemic aneuploidies. More recently, it has contributed to the identification and molecular dissection of aneuploidy-induced cellular behaviors and their impact at the tissue and organismal levels. In this perspective, we review this active field of research, first by comparing knowledge from yeast, mouse, and human cells, then by highlighting the contributions of Drosophila. The aim of these discussions was to further our understanding of the functional interplay between aneuploidy, cell physiology, and tissue homeostasis in human development and disease.

The sex-specific factor SOA controls dosage compensation in Anopheles mosquitoes

The Anopheles mosquito is one of thousands of species in which sex differences play a central part in their biology, as only females need a blood meal to produce eggs. Sex differentiation is regulated by sex chromosomes, but their presence creates a dosage imbalance between males (XY) and females (XX). Dosage compensation (DC) can re-equilibrate the expression of sex chromosomal genes. However, because DC mechanisms have only been fully characterized in a few model organisms, key questions about its evolutionary diversity and functional necessity remain unresolved1. Here we report the discovery of a previously uncharacterized gene (sex chromosome activation (SOA)) as a master regulator of DC in the malaria mosquito Anopheles gambiae. Sex-specific alternative splicing prevents functional SOA protein expression in females. The male isoform encodes a DNA-binding protein that binds the promoters of active X chromosomal genes. Expressing male SOA is sufficient to induce DC in female cells. Male mosquitoes lacking SOA or female mosquitoes ectopically expressing the male isoform exhibit X chromosome misregulation, which is compatible with viability but causes developmental delay. Thus, our molecular analyses of a DC master regulator in a non-model organism elucidates the evolutionary steps that lead to the establishment of a chromosome-specific fine-tuning mechanism.

Viruses in Laboratory Drosophila and Their Impact on Host Gene Expression

Drosophila melanogaster has one of the best characterized antiviral immune responses among invertebrates. However, relatively few easily transmitted natural virus isolates are available, and so many Drosophila experiments have been performed using artificial infection routes and artificial host-virus combinations. These may not reflect natural infections, especially for subtle phenotypes such as gene expression. Here, to explore the laboratory virus community and to better understand how natural virus infections induce changes in gene expression, we have analysed seven publicly available D. melanogaster transcriptomic sequencing datasets that were originally sequenced for projects unrelated to virus infection. We have found ten known viruses-including five that have not been experimentally isolated-but no previously unknown viruses. Our analysis of host gene expression revealed that numerous genes were differentially expressed in flies that were naturally infected with a virus. For example, flies infected with nora virus showed patterns of gene expression consistent with intestinal vacuolization and possible host repair via the upd3 JAK/STAT pathway. We also found marked sex differences in virus-induced differential gene expression. Our results show that natural virus infection in laboratory Drosophila does indeed induce detectable changes in gene expression, suggesting that this may form an important background condition for experimental studies in the laboratory.

Extensive Copy Number Variation Explains Genome Size Variation in the Unicellular Zygnematophycean Alga, Closterium peracerosum-strigosum-littorale Complex

Genome sizes are known to vary within and among closely related species, but the knowledge about genomic factors contributing to the variation and their impacts on gene functions is limited to only a small number of species. This study identified a more than 2-fold heritable genome size variation among the unicellular Zygnematophycean alga, Closterium peracerosum-strigosum-littorale (C. psl.) complex, based on short-read sequencing analysis of 22 natural strains and F1 segregation analysis. Six de novo assembled genomes revealed that genome size variation is largely attributable to genome-wide copy number variation (CNV) among strains rather than mating type-linked genomic regions or specific repeat sequences such as rDNA. Notably, about 30% of genes showed CNV even between strains that can mate with each other. Transcriptome and gene ontology analysis demonstrated that CNV is distributed nonrandomly in terms of gene functions, such that CNV was more often observed in the gene set with stage-specific expression. Furthermore, in about 30% of these genes with CNV, the expression level does not increase proportionally with the gene copy number, suggesting presence of dosage compensation, which was overrepresented in genes involved in basic biological functions, such as translation. Nonrandom patterns in gene duplications and corresponding expression changes in terms of gene functions may contribute to maintaining the high level of CNV associated with extensive genome size variation in the C. psl. complex, despite its possible detrimental effects.

NetREX-CF integrates incomplete transcription factor data with gene expression to reconstruct gene regulatory networks

The inference of Gene Regulatory Networks (GRNs) is one of the key challenges in systems biology. Leading algorithms utilize, in addition to gene expression, prior knowledge such as Transcription Factor (TF) DNA binding motifs or results of TF binding experiments. However, such prior knowledge is typically incomplete, therefore, integrating it with gene expression to infer GRNs remains difficult. To address this challenge, we introduce NetREX-CF-Regulatory Network Reconstruction using EXpression and Collaborative Filtering-a GRN reconstruction approach that brings together Collaborative Filtering to address the incompleteness of the prior knowledge and a biologically justified model of gene expression (sparse Network Component Analysis based model). We validated the NetREX-CF using Yeast data and then used it to construct the GRN for Drosophila Schneider 2 (S2) cells. To corroborate the GRN, we performed a large-scale RNA-Seq analysis followed by a high-throughput RNAi treatment against all 465 expressed TFs in the cell line. Our knockdown result has not only extensively validated the GRN we built, but also provides a benchmark that our community can use for evaluating GRNs. Finally, we demonstrate that NetREX-CF can infer GRNs using single-cell RNA-Seq, and outperforms other methods, by using previously published human data.

Chromosome-level Assembly, Dosage Compensation and Sex-biased Gene Expression in the Small Brown Planthopper, Laodelphax striatellus

In insects, sex chromosome differentiation often results in unequal gene dosages between sexes. Dosage compensation mechanisms evolve to balance gene expression, but the degree and mechanism of regulation often vary by insect species. In hemipteran species, the small brown planthopper (SBPH), Laodelphax striatellus, is an injurious crop pest, with a sex chromosome type XX in females and XO in males. This species offers the opportunity to study dosage compensation and sex-biased gene expression. In this study, we generated a chromosome-level genome of SBPH using Oxford Nanopore Technologies and high-throughput chromatin conformation capture (Hi-C) technology. We also sequenced RNA-seq data from 16 tissue samples to annotate the genome and analyze gene dosage compensation. We finally obtained a 510.2 megabases (Mb) genome with 99.12% of the scaffolds anchored on 15 chromosomes (14 autosomes and 1 X chromosome) and annotated 16,160 protein-coding genes based on full-length cDNA sequencing data. Furthermore, we found complete dosage compensation in all L. striatellus somatic tissues, but lack of dosage compensation in gonad tissue testis. We also found that female-biased genes were significantly enriched on the X chromosome in all tissues, whereas male-biased genes in gonad tissues were enriched on autosomes. This study not only provides a high-quality genome assembly but also lays a foundation for a better understanding of the sexual regulatory network in hemipteran insects.

Altering enhancer-promoter linear distance impacts promoter competition in cis and in trans

In Drosophila, pairing of maternal and paternal homologs can permit trans-interactions between enhancers on one homolog and promoters on another, an example of a phenomenon called transvection. When chromosomes are paired, promoters in cis and in trans to an enhancer can compete for the enhancer's activity, but the parameters that govern this competition are as yet poorly understood. To assess how the linear spacing between an enhancer and promoter can influence promoter competition in Drosophila, we employed transgenic constructs wherein the eye-specific enhancer GMR is placed at varying distances from a heterologous hsp70 promoter driving a fluorescent reporter. While GMR activates the reporter to a high degree when the enhancer and promoter are spaced by a few hundred base pairs, activation is strongly attenuated when the enhancer is moved 3 kilobases away. By examining transcription of endogenous genes near the point of transgene insertion, we show that linear spacing of 3 kb between GMR and the hsp70 promoter results in elevated transcription of neighboring promoters, suggesting a loss of specificity between the enhancer and its intended transgenic target promoter. Furthermore, increasing spacing between GMR and hsp70 by just 100 bp can enhance transvection, resulting in increased activation of a promoter on a paired homolog at the expense of a promoter in cis to the enhancer. Finally, cis-/trans-promoter competition assays in which one promoter carries mutations to key core promoter elements show that GMR will skew its activity toward a wild type promoter, suggesting that an enhancer is in a balanced competition between its potential target promoters in cis and in trans.

The good, the bad, and the ugly: Evolutionary and pathological aspects of gene dosage alterations

Diploid organisms contain a maternal and a paternal genome complement that is thought to provide robustness and allow developmental progression despite genetic perturbations that occur in heterozygosity. However, changes affecting gene dosage from the chromosome down to the individual gene level possess a significant pathological potential and can lead to developmental disorders (DDs). This indicates that expression from a balanced gene complement is highly relevant for proper cellular and organismal function in eukaryotes. Paradoxically, gene and whole chromosome duplications are a principal driver of evolution, while heteromorphic sex chromosomes (XY and ZW) are naturally occurring aneuploidies important for sex determination. Here, we provide an overview of the biology of gene dosage at the crossroads between evolutionary benefit and pathogenicity during disease. We describe the buffering mechanisms and cellular responses to alterations, which could provide a common ground for the understanding of DDs caused by copy number alterations.

3D genome organisation in Drosophila

Ever since Thomas Hunt Morgan's discovery of the chromosomal basis of inheritance by using Drosophila melanogaster as a model organism, the fruit fly has remained an essential model system in studies of genome biology, including chromatin organisation. Very much as in vertebrates, in Drosophila, the genome is organised in territories, compartments and topologically associating domains (TADs). However, these domains might be formed through a slightly different mechanism than in vertebrates due to the presence of a large and potentially redundant set of insulator proteins and the minor role of dCTCF in TAD boundary formation. Here, we review the different levels of chromatin organisation in Drosophila and discuss mechanisms and factors that might be involved in TAD formation. The dynamics of TADs and enhancer-promoter interactions in the context of transcription are covered in the light of currently conflicting results. Finally, we illustrate the value of polymer modelling approaches to infer the principles governing the three-dimensional organisation of the Drosophila genome.

The haplolethality paradox of the wupA gene in Drosophila

Haplolethals (HL) are regions of diploid genomes that in one dose are fatal for the organism. Their biological meaning is obscure because heterozygous loss-of-function mutations result in dominant lethality (DL) and, consequently, should be under strong negative selection. We report an in depth study of the HL associated to the gene wings up A (wupA). It encodes 13 transcripts (A-M) that yield 11 protein isoforms (A-K) of Troponin I (TnI). They are functionally diverse in their control of muscle contraction, cell polarity and cell proliferation. Isoform K transfers to the nucleus where it increases transcription of the cell proliferation related genes CDK2, CDK4, Rap and Rab5. The nuclear translocation of isoform K is prevented by the co-expression of A or B isoforms, which illustrates isoform interactions. The corresponding DL mutations are, either DNA rearrangements clustered towards the gene 3’ end, thus affecting the genomic organization of all transcripts, or CRISPR-induced mutations in one of the two ATG sites which eliminate a subset of wupA products. The joint elimination of isoforms C, F, G and H, however, do not cause DL phenotypes. Genetically driven expression of single isoforms rescue neither DL nor any of the mutants known in the gene, suggesting that normal function requires properly regulated expression of specific combinations, rather than single, TnI isoforms. We conclude that the wupA associated HL results from the combined haploinsufficiency of a large set of TnI isoforms. The qualitative and quantitative normal expression of which, requires the chromosomal integrity of the wupA genomic region. Since all fly TnI isoforms are encoded in the same gene, its HL condition becomes unavoidable. These wupA features are comparable to those of dpp, the only other HL studied to some extent, and reveal a scenario of strict dosage dependence with implications for gene expression regulation and splitting.

Most species contain two copies of their genetic endowment, each received from each progenitor. If one of the duplicated genes is non-functional, the remaining copy may supply enough product as to cover the requirements for normal function or, alternatively, may reflect the insufficiency through a visible phenotype. In rare occasions, however, having one copy is so deleterious that causes lethality. These so called “haplolethal regions”, exist across species and represent an evolutionary paradox since they should have been subject to intense negative selection. The inherent difficulties to study haplolethals have precluded their study so far. Here, we analyzed the case of one of the five haplolethal regions of Drosophila, the one associated to the Troponin I encoding gene wupA, by measuring the transcriptional effects of mutations and chromosomal rearrangements affecting this gene. The data show that this haplolethality results from the combined insufficiency of a large number of Troponin I isoforms, which are functionally specialized, show interference and require the integrity of the native chromatin structure for their quantitatively regulated expression. These features unveil novel aspects of gene expression and, possibly, on evolutionary gene splitting. Also, haplolethals underscore the biological significance of protein dosage, in particular for functionally related products.

Dynamic sex chromosome expression in Drosophila male germ cells

Given their copy number differences and unique modes of inheritance, the evolved gene content and expression of sex chromosomes is unusual. In many organisms the X and Y chromosomes are inactivated in spermatocytes, possibly as a defense mechanism against insertions into unpaired chromatin. In addition to current sex chromosomes, Drosophila has a small gene-poor X-chromosome relic (4th) that re-acquired autosomal status. Here we use single cell RNA-Seq on fly larvae to demonstrate that the single X and pair of 4th chromosomes are specifically inactivated in primary spermatocytes, based on measuring all genes or a set of broadly expressed genes in testis we identified. In contrast, genes on the single Y chromosome become maximally active in primary spermatocytes. Reduced X transcript levels are due to failed activation of RNA-Polymerase-II by phosphorylation of Serine 2 and 5.

Independence of chromatin conformation and gene regulation during Drosophila dorsoventral patterning

The relationship between chromatin organization and gene regulation remains unclear. While disruption of chromatin domains and domain boundaries can lead to misexpression of developmental genes, acute depletion of regulators of genome organization has a relatively small effect on gene expression. It is therefore uncertain whether gene expression and chromatin state drive chromatin organization or whether changes in chromatin organization facilitate cell-type-specific activation of gene expression. Here, using the dorsoventral patterning of the Drosophila melanogaster embryo as a model system, we provide evidence for the independence of chromatin organization and dorsoventral gene expression. We define tissue-specific enhancers and link them to expression patterns using single-cell RNA-seq. Surprisingly, despite tissue-specific chromatin states and gene expression, chromatin organization is largely maintained across tissues. Our results indicate that tissue-specific chromatin conformation is not necessary for tissue-specific gene expression but rather acts as a scaffold facilitating gene expression when enhancers become active.

3D genome evolution and reorganization in the Drosophila melanogaster species group

Topologically associating domains, or TADs, are functional units that organize chromosomes into 3D structures of interacting chromatin. TADs play an important role in regulating gene expression by constraining enhancer-promoter contacts and there is evidence that deletion of TAD boundaries leads to aberrant expression of neighboring genes. While the mechanisms of TAD formation have been well-studied, current knowledge on the patterns of TAD evolution across species is limited. Due to the integral role TADs play in gene regulation, their structure and organization is expected to be conserved during evolution. However, more recent research suggests that TAD structures diverge relatively rapidly. We use Hi-C chromosome conformation capture to measure evolutionary conservation of whole TADs and TAD boundary elements between D. melanogaster and D. triauraria, two early-branching species from the melanogaster species group which diverged ∼15 million years ago. We find that the majority of TADs have been reorganized since the common ancestor of D. melanogaster and D. triauraria, via a combination of chromosomal rearrangements and gain/loss of TAD boundaries. TAD reorganization between these two species is associated with a localized effect on gene expression, near the site of disruption. By separating TADs into subtypes based on their chromatin state, we find that different subtypes are evolving under different evolutionary forces. TADs enriched for broadly expressed, transcriptionally active genes are evolving rapidly, potentially due to positive selection, whereas TADs enriched for developmentally-regulated genes remain conserved, presumably due to their importance in restricting gene-regulatory element interactions. These results provide novel insight into the evolutionary dynamics of TADs and help to reconcile contradictory reports related to the evolutionary conservation of TADs and whether changes in TAD structure affect gene expression.

Chromosomes are organized inside the three-dimensional nucleus in a hierarchical fashion, from the wrapping of DNA around nucleosomes to the presence of compartments of active and inactive chromatin. Here, we focus on domains of interacting chromatin known as topologically associating domains, or TADs. TADs are dozens to hundreds of kilobases in size and tend to contain genes with similar epigenetic states. While much research has been done on the mechanisms of TAD formation, less is known about the role of TADs in genome evolution. The research that does exist is contradictory, with some studies claiming that TADs are highly conserved and others reporting that TADs diverge rapidly. We compared the TAD locations between two Drosophila species and found that the majority have been reorganized. We also identified TAD subtypes that differ in terms of their evolutionary conservation: TADs that contain developmental genes are less likely to be reorganized while TADs that contain broadly expressed genes evolve rapidly. We found that TAD reorganization has an effect on gene expression, but only for genes located near the site of disruption. These results provide new insight into the evolution of 3D genome organization and help to reconcile previous contradictions in the literature.

Reduced Neurog3 Gene Dosage Shifts Enteroendocrine Progenitor Towards Goblet Cell Lineage in the Mouse Intestine

BACKGROUND & AIMS

Transient expression of Neurog3 commits intestinal secretory progenitors to become enteroendocrine-biased progenitors and hence drive enteroendocrine differentiation. Loss of Neurog3 in mouse resulted in the depletion of intestinal enteroendocrine cells (EECs) and an increase in goblet cells. Earlier studies in developing mouse pancreas identified a role of Neurog3 gene dosage in endocrine and exocrine cell fate allocation. We aimed to determine whether Neurog3 gene dosage controls fate choice of enteroendocrine progenitors.

METHODS

We acquired mutant Neurog3 reporter mice carrying 2, 1, or null Neurog3 alleles to study Neurog3 gene dosage effect by lineage tracing. Cell types arising from Neurog3+ progenitors were determined by immunohistochemistry using antibodies against intestinal lineage-specific markers. RNA sequencing of sorted Neurog3^+/+, Neurog3^+/-, or bulk intestinal cells were performed and differentially expressed genes were analyzed.

RESULTS

We identified 2731 genes enriched in sorted Neurog3^+/+-derived cells in the Neurog3^+/+EYFP mouse intestine when compared with bulk duodenum epithelial cells. In the intestine of Neurog3^+/-EGFP heterozygous mouse, we observed a 63% decrease in EEC numbers. Many Neurog3-derived cells stained for goblet marker Mucin 2. RNA sequencing of sorted Neurog3^+/- cells uncovered enriched expression of genes characteristic for both goblet and enteroendocrine cells, indicating the mixed lineages arose from Neurog3+ progenitors. Consistent with this hypothesis, deletion of both Neurog3 alleles resulted in the total absence of EECs. All Neurog3⁺-derived cells stained for Mucin 2.

CONCLUSIONS

We identified that the fate of Neurog3⁺ enteroendocrine progenitors is dependent on Neurog3 gene dosage. High Neurog3 gene dosage enforces the commitment of secretory progenitors to an EE lineage, while constraining their goblet cell lineage potential. Transcriptome profiling data was deposited to Gene Ontology omnibus, accession number: GSE149203.

Drosophila Heterochromatin Stabilization Requires the Zinc-Finger Protein Small Ovary

Heterochromatin-mediated repression is essential for controlling the expression of transposons and for coordinated cell type-specific gene regulation. The small ovary (sov) locus was identified in a screen for female-sterile mutations in Drosophila melanogaster, and mutants show dramatic ovarian morphogenesis defects. We show that the null sov phenotype is lethal and map the locus to the uncharacterized gene CG14438, which encodes a nuclear zinc-finger protein that colocalizes with the essential Heterochromatin Protein 1 (HP1a). We demonstrate Sov functions to repress inappropriate gene expression in the ovary, silence transposons, and suppress position-effect variegation in the eye, suggesting a central role in heterochromatin stabilization.

Patterns and Constraints in the Evolution of Sperm Individualization Genes in Insects, with an Emphasis on Beetles

Gene expression profiles can change dramatically between sexes and sex bias may contribute specific macroevolutionary dynamics for sex-biased genes. However, these dynamics are poorly understood at large evolutionary scales due to the paucity of studies that have assessed orthology and functional homology for sex-biased genes and the pleiotropic effects possibly constraining their evolutionary potential. Here, we explore the correlation of sex-biased expression with macroevolutionary processes that are associated with sex-biased genes, including duplications and accelerated evolutionary rates. Specifically, we examined these traits in a group of 44 genes that orchestrate sperm individualization during spermatogenesis, with both unbiased and sex-biased expression. We studied these genes in the broad evolutionary framework of the Insecta, with a particular focus on beetles (order Coleoptera). We studied data mined from 119 insect genomes, including 6 beetle models, and from 19 additional beetle transcriptomes. For the subset of physically and/or genetically interacting proteins, we also analyzed how their network structure may condition the mode of gene evolution. The collection of genes was highly heterogeneous in duplication status, evolutionary rates, and rate stability, but there was statistical evidence for sex bias correlated with faster evolutionary rates, consistent with theoretical predictions. Faster rates were also correlated with clocklike (insect amino acids) and non-clocklike (beetle nucleotides) substitution patterns in these genes. Statistical associations (higher rates for central nodes) or lack thereof (centrality of duplicated genes) were in contrast to some current evolutionary hypotheses, highlighting the need for more research on these topics.

Highly rearranged chromosomes reveal uncoupling between genome topology and gene expression

Chromatin topology is intricately linked to gene expression, yet its functional requirement remains unclear. Here, we comprehensively assessed the interplay between genome topology and gene expression using highly rearranged chromosomes (balancers) spanning ~75% of the Drosophila genome. Using transheterozyte (balancer/wild-type) embryos, we measured allele-specific changes in topology and gene expression in cis, whilst minimizing trans effects. Through genome sequencing, we resolved eight large nested inversions, smaller inversions, duplications, and thousands of deletions. These extensive rearrangements caused many changes to chromatin topology, including long-range loops, TADs and promoter interactions, yet these are not predictive of changes in expression. Gene expression is generally not altered around inversion breakpoints, indicating that mis-appropriate enhancer-promoter activation is a rare event. Similarly, shuffling or fusing TADs, changing intra-TAD connections and disrupting long-range inter-TAD loops, does not alter expression for the majority of genes. Our results suggest that properties other than chromatin topology ensure productive enhancer-promoter interactions.

Heightened condition-dependence of the sexual transcriptome as a function of genetic quality in Drosophila melanogaster head tissue

Theory suggests sexual traits should show heightened condition-dependent expression. This prediction has been tested extensively in experiments where condition has been manipulated through environmental quality. Condition-dependence as a function of genetic quality has, however, only rarely been addressed, despite its central importance in evolutionary theory. To address the effect of genetic quality on expression of sexual and non-sexual traits, we here compare gene expression in Drosophila melanogaster head tissue between flies with intact genomes (high condition) and flies carrying a major deleterious mutation (low condition). We find that sex-biased genes show heightened condition-dependent expression in both sexes, and that expression in low condition males and females regresses towards a more similar expression profile. As predicted, sex-biased expression was more sensitive to condition in males compared to females, but surprisingly female-biased, rather than male-biased, genes show higher sensitivity to condition in both sexes. Our results thus support the fundamental predictions of the theory of condition-dependence when condition is a function of genetic quality.

Coevolution of Genome Architecture and Social Behavior

Although social behavior can have a strong genetic component, it can also result in selection on genome structure and function, thereby influencing the evolution of the genome itself. Here we explore the bidirectional links between social behavior and genome architecture by considering variation in social and/or mating behavior among populations (social polymorphisms) and across closely related species. We propose that social behavior can influence genome architecture via associated demographic changes due to social living. We establish guidelines to exploit emerging whole-genome sequences using analytical approaches that examine genome structure and function at different levels (regulatory vs structural variation) from the perspective of both molecular biology and population genetics in an ecological context.

Global impacts of chromosomal imbalance on gene expression in Arabidopsis and other taxa

Changes in dosage of part of the genome (aneuploidy) have long been known to produce much more severe phenotypic consequences than changes in the number of whole genomes (ploidy). To examine the basis of these differences, global gene expression in mature leaf tissue for all five trisomies and in diploids, triploids, and tetraploids of Arabidopsis thaliana was studied. The trisomies displayed a greater spread of expression modulation than the ploidy series. In general, expression of genes on the varied chromosome ranged from compensation to dosage effect, whereas genes from the remainder of the genome ranged from no effect to reduced expression approaching the inverse level of chromosomal imbalance (2/3). Genome-wide DNA methylation was examined in each genotype and found to shift most prominently with trisomy 4 but otherwise exhibited little change, indicating that genetic imbalance is generally mechanistically unrelated to DNA methylation. Independent analysis of gene functional classes demonstrated that ribosomal, proteasomal, and gene body methylated genes were less modulated compared with all classes of genes, whereas transcription factors, signal transduction components, and organelle-targeted protein genes were more tightly inversely affected. Comparing transcription factors and their targets in the trisomies and in expression networks revealed considerable discordance, illustrating that altered regulatory stoichiometry is a major contributor to genetic imbalance. Reanalysis of published data on gene expression in disomic yeast and trisomic mouse cells detected similar stoichiometric effects across broad phylogenetic taxa, and indicated that these effects reflect normal gene regulatory processes.

Non-canonical Drosophila X chromosome dosage compensation and repressive topologically associated domains

Background

In animals with XY sex chromosomes, X-linked genes from a single X chromosome in males are imbalanced relative to autosomal genes. To minimize the impact of genic imbalance in male Drosophila, there is a dosage compensation complex (MSL) that equilibrates X-linked gene expression with the autosomes. There are other potential contributions to dosage compensation. Hemizygous autosomal genes located in repressive chromatin domains are often derepressed. If this homolog-dependent repression occurs on the X, which has no pairing partner, then derepression could contribute to male dosage compensation.

Results

We asked whether different chromatin states or topological associations correlate with X chromosome dosage compensation, especially in regions with little MSL occupancy. Our analyses demonstrated that male X chromosome genes that are located in repressive chromatin states are depleted of MSL occupancy; however, they show dosage compensation. The genes in these repressive regions were also less sensitive to knockdown of MSL components.

Conclusions

Our results suggest that this non-canonical dosage compensation is due to the same transacting derepression that occurs on autosomes. This mechanism would facilitate immediate compensation during the evolution of sex chromosomes from autosomes. This mechanism is similar to that of C. elegans, where enhanced recruitment of X chromosomes to the nuclear lamina dampens X chromosome expression as part of the dosage compensation response in XX individuals.

Electronic supplementary material

The online version of this article (10.1186/s13072-018-0232-y) contains supplementary material, which is available to authorized users.

Reprogramming of regulatory network using expression uncovers sex-specific gene regulation in Drosophila

Gene regulatory networks (GRNs) describe regulatory relationships between transcription factors (TFs) and their target genes. Computational methods to infer GRNs typically combine evidence across different conditions to infer context-agnostic networks. We develop a method, Network Reprogramming using EXpression (NetREX), that constructs a context-specific GRN given context-specific expression data and a context-agnostic prior network. NetREX remodels the prior network to obtain the topology that provides the best explanation for expression data. Because NetREX utilizes prior network topology, we also develop PriorBoost, a method that evaluates a prior network in terms of its consistency with the expression data. We validate NetREX and PriorBoost using the “gold standard” E. coli GRN from the DREAM5 network inference challenge and apply them to construct sex-specific Drosophila GRNs. NetREX constructed sex-specific Drosophila GRNs that, on all applied measures, outperform networks obtained from other methods indicating that NetREX is an important milestone toward building more accurate GRNs.

For many applications knowledge of context-specific gene regulatory networks (GRNs) is desirable, but their inference remains a challenge. Here, the authors introduce a method for construction of context-specific GRNs, and apply it to construct sex-specific Drosophila GRNs.

Facultative dosage compensation of developmental genes on autosomes in Drosophila and mouse embryonic stem cells

Haploinsufficiency and aneuploidy are two phenomena, where gene dosage alterations cause severe defects ultimately resulting in developmental failures and disease. One remarkable exception is the X chromosome, where copy number differences between sexes are buffered by dosage compensation systems. In Drosophila, the Male-Specific Lethal complex (MSLc) mediates upregulation of the single male X chromosome. The evolutionary origin and conservation of this process orchestrated by MSL2, the only male-specific protein within the fly MSLc, have remained unclear. Here, we report that MSL2, in addition to regulating the X chromosome, targets autosomal genes involved in patterning and morphogenesis. Precise regulation of these genes by MSL2 is required for proper development. This set of dosage-sensitive genes maintains such regulation during evolution, as MSL2 binds and similarly regulates mouse orthologues via Histone H4 lysine 16 acetylation. We propose that this gene-by-gene dosage compensation mechanism was co-opted during evolution for chromosome-wide regulation of the Drosophila male X.

In Drosophila the Male-Specific Lethal complex (MSLc) mediates upregulation of the single male X chromosome. Here the authors provide evidence that MSL2 also targets autosomal genes required for proper development and that MSL2 binds and similarly regulates mouse orthologues.

Dosage-Dependent Expression Variation Suppressed on the Drosophila Male X Chromosome

DNA copy number variation is associated with many high phenotypic heterogeneity disorders. We systematically examined the impact of Drosophila melanogaster deletions on gene expression profiles to ask whether increased expression variability owing to reduced gene dose might underlie this phenotypic heterogeneity. Indeed, we found that one-dose genes have higher gene expression variability relative to two-dose genes. We then asked whether this increase in variability could be explained by intrinsic noise within cells due to stochastic biochemical events, or whether expression variability is due to extrinsic noise arising from more complex interactions. Our modeling showed that intrinsic gene expression noise averages at the organism level and thus cannot explain increased variation in one-dose gene expression. Interestingly, expression variability was related to the magnitude of expression compensation, suggesting that regulation, induced by gene dose reduction, is noisy. In a remarkable exception to this rule, the single X chromosome of males showed reduced expression variability, even compared with two-dose genes. Analysis of sex-transformed flies indicates that X expression variability is independent of the male differentiation program. Instead, we uncovered a correlation between occupancy of the chromatin-modifying protein encoded by males absent on the first (mof) and expression variability, linking noise suppression to the specialized X chromosome dosage compensation system. MOF occupancy on autosomes in both sexes also lowered transcriptional noise. Our results demonstrate that gene dose reduction can lead to heterogeneous responses, which are often noisy. This has implications for understanding gene network regulatory interactions and phenotypic heterogeneity. Additionally, chromatin modification appears to play a role in dampening transcriptional noise.

Engineering the Drosophila Genome for Developmental Biology

The recent development of transposon and CRISPR-Cas9-based tools for manipulating the fly genome in vivo promises tremendous progress in our ability to study developmental processes. Tools for introducing tags into genes at their endogenous genomic loci facilitate imaging or biochemistry approaches at the cellular or subcellular levels. Similarly, the ability to make specific alterations to the genome sequence allows much more precise genetic control to address questions of gene function.

Drosophila melanogaster positive transcriptional elongation factors regulate metabolic and sex-biased expression in adults

Background

Transcriptional elongation is a generic function, but is also regulated to allow rapid transcription responses. Following relatively long initiation and promoter clearance, RNA polymerase II can pause and then rapidly elongate following recruitment of positive elongation factors. Multiple elongation complexes exist, but the role of specific components in adult Drosophila is underexplored.

Results

We conducted RNA-seq experiments to analyze the effect of RNAi knockdown of Suppressor of Triplolethal and lilliputian. We similarly analyzed the effect of expressing a dominant negative Cyclin-dependent kinase 9 allele. We observed that almost half of the genes expressed in adults showed reduced expression, supporting a broad role for the three tested genes in steady-state transcript abundance. Expression profiles following lilliputian and Suppressor of Triplolethal RNAi were nearly identical raising the possibility that they are obligatory co-factors. Genes showing reduced expression due to these RNAi treatments were short and enriched for genes encoding metabolic or enzymatic functions. The dominant-negative Cyclin-dependent kinase 9 profiles showed both overlapping and specific differential expression, suggesting involvement in multiple complexes. We also observed hundreds of genes with sex-biased differential expression following treatment.

Conclusion

Transcriptional profiles suggest that Lilliputian and Suppressor of Triplolethal are obligatory cofactors in the adult and that they can also function with Cyclin-dependent kinase 9 at a subset of loci. Our results suggest that transcriptional elongation control is especially important for rapidly expressed genes to support digestion and metabolism, many of which have sex-biased function.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3755-x) contains supplementary material, which is available to authorized users.