Neurons That Underlie Drosophila melanogaster Reproductive Behaviors: Detection of a Large Male-Bias in Gene Expression in fruitless-Expressing Neurons

Sex-Biased Expression Is Associated With Chromatin State in Drosophila melanogaster and Drosophila simulans

In Drosophila melanogaster and D. simulans head tissue, 60% of orthologous genes show evidence of sex-biased expression in at least one species. Of these, ∼39% (2,192) are conserved in direction. We hypothesize enrichment of open chromatin in the sex where we see expression bias and closed chromatin in the opposite sex. Male-biased orthologs are significantly enriched for H3K4me3 marks in males of both species (∼89% of male-biased orthologs vs. ∼76% of unbiased orthologs). Similarly, female-biased orthologs are significantly enriched for H3K4me3 marks in females of both species (∼90% of female-biased orthologs vs. ∼73% of unbiased orthologs). The sex-bias ratio in female-biased orthologs was similar in magnitude between the two species, regardless of the closed chromatin (H3K27me2me3) marks in males. However, in male-biased orthologs, the presence of H3K27me2me3 in both species significantly reduced the correlation between D. melanogaster sex-bias ratio and the D. simulans sex-bias ratio. Male-biased orthologs are enriched for evidence of positive selection in the D. melanogaster group. There are more male-biased genes than female-biased genes in both species. For orthologs with gains/losses of sex-bias between the two species, there is an excess of male-bias compared to female-bias, but there is no consistent pattern in the relationship between H3K4me3 or H3K27me2me3 chromatin marks and expression. These data suggest chromatin state is a component of the maintenance of sex-biased expression and divergence of sex-bias between species is reflected in the complexity of the chromatin status.

Single-cell transcriptome profiles of Drosophila fruitless-expressing neurons from both sexes

Drosophila melanogaster reproductive behaviors are orchestrated by fruitless neurons. We performed single-cell RNA-sequencing on pupal neurons that produce sex-specifically spliced fru transcripts, the fru P1-expressing neurons. Uniform Manifold Approximation and Projection (UMAP) with clustering generates an atlas containing 113 clusters. While the male and female neurons overlap in UMAP space, more than half the clusters have sex differences in neuron number, and nearly all clusters display sex-differential expression. Based on an examination of enriched marker genes, we annotate clusters as circadian clock neurons, mushroom body Kenyon cell neurons, neurotransmitter- and/or neuropeptide-producing, and those that express doublesex. Marker gene analyses also show that genes that encode members of the immunoglobulin superfamily of cell adhesion molecules, transcription factors, neuropeptides, neuropeptide receptors, and Wnts have unique patterns of enriched expression across the clusters. In vivo spatial gene expression links to the clusters are examined. A functional analysis of fru P1 circadian neurons shows they have dimorphic roles in activity and period length. Given that most clusters are comprised of male and female neurons indicates that the sexes have fru P1 neurons with common gene expression programs. Sex-specific expression is overlaid on this program, to build the potential for vastly different sex-specific behaviors.

Sex-biased expression is associated with chromatin state in D. melanogaster and D. simulans

We propose a new model for the association of chromatin state and sex-bias in expression. We hypothesize enrichment of open chromatin in the sex where we see expression bias (OS) and closed chromatin in the opposite sex (CO). In this study of D. melanogaster and D. simulans head tissue, sex-bias in expression is associated with H3K4me3 (open mark) in males for male-biased genes and in females for female-biased genes in both species. Sex-bias in expression is also largely conserved in direction and magnitude between the two species on the X and autosomes. In male-biased orthologs, the sex-bias ratio is more divergent between species if both species have H3K27me2me3 marks in females compared to when either or neither species has H3K27me2me3 in females. H3K27me2me3 marks in females are associated with male-bias in expression on the autosomes in both species, but on the X only in D. melanogaster . In female-biased orthologs the relationship between the species for the sex-bias ratio is similar regardless of the H3K27me2me3 marks in males. Female-biased orthologs are more similar in the ratio of sex-bias than male-biased orthologs and there is an excess of male-bias in expression in orthologs that gain/lose sex-bias. There is an excess of male-bias in sex-limited expression in both species suggesting excess male-bias is due to rapid evolution between the species. The X chromosome has an enrichment in male-limited H3K4me3 in both species and an enrichment of sex-bias in expression compared to the autosomes.

Molecular Mechanisms of Sexually Dimorphic Nervous System Patterning in Flies and Worms

Male and female brains display anatomical and functional differences. Such differences are observed in species across the animal kingdom, including humans, but have been particularly well-studied in two classic animal model systems, the fruit fly Drosophila melanogaster and the nematode Caenorhabditis elegans. Here we summarize recent advances in understanding how the worm and fly brain acquire sexually dimorphic features during development. We highlight the advantages of each system, illustrating how the precise anatomical delineation of sexual dimorphisms in worms has enabled recent analysis into how these dimorphisms become specified during development, and how focusing on sexually dimorphic neurons in the fly has enabled an increasingly detailed understanding of sex-specific behaviors.

Sexual antagonism, temporally fluctuating selection, and variable dominance affect a regulatory polymorphism in Drosophila melanogaster

Understanding how genetic variation is maintained within species is a major goal of evolutionary genetics that can shed light on the preservation of biodiversity. Here, we examined the maintenance of a regulatory single-nucleotide polymorphism (SNP) of the X-linked Drosophila melanogaster gene fezzik. The derived variant at this site is at intermediate frequency in many worldwide populations but absent in populations from the ancestral species range in sub-Saharan Africa. We collected and genotyped wild-caught individuals from a single European population biannually over a period of 5 years, which revealed an overall difference in allele frequency between the sexes and a consistent change in allele frequency across seasons in females but not in males. Modeling based on the observed allele and genotype frequencies suggested that both sexually antagonistic and temporally fluctuating selection may help maintain variation at this site. The derived variant is predicted to be female-beneficial and mostly recessive; however, there was uncertainty surrounding our dominance estimates and long-term modeling projections suggest that it is more likely to be dominant. By examining gene expression phenotypes, we found that phenotypic dominance was variable and dependent upon developmental stage and genetic background, suggesting that dominance may be variable at this locus. We further determined that fezzik expression and genotype are associated with starvation resistance in a sex-dependent manner, suggesting a potential phenotypic target of selection. By characterizing the mechanisms of selection acting on this SNP, our results improve our understanding of how selection maintains genetic and phenotypic variation in natural populations.

Dopamine neuron-specific LRRK2 G2019S effects on gene expression revealed by translatome profiling

Leucine-rich repeat kinase 2 (LRRK2) mutations are the most common genetic cause of late-onset Parkinson's disease. The pathogenic G2019S mutation enhances LRRK2 kinase activity and induces neurodegeneration in C. elegans, Drosophila and rodent models through unclear mechanisms. Gene expression profiling has the potential to provide detailed insight into the biological pathways modulated by LRRK2 kinase activity. Prior in vivo studies have surveyed the effects of LRRK2 G2019S on genome-wide mRNA expression in complex brain tissues with high cellular heterogeneity, limiting their power to detect more restricted gene expression changes occurring in a cell type-specific manner. Here, we used translating ribosome affinity purification (TRAP) coupled to RNA-seq to profile dopamine neuron-specific gene expression changes caused by LRRK2 G2019S in the Drosophila CNS. A number of genes were differentially expressed in the presence of mutant LRRK2 that represent a broad range of molecular functions including DNA repair (RfC3), mRNA metabolism and translation (Ddx1 and lin-28), calcium homeostasis (MCU), and other categories (Ugt37c1, disp, l(1)G0196, CG6602, CG1126 and CG11068). Further analysis on a subset of these genes revealed that LRRK2 G2019S did not alter their expression across the whole brain, consistent with dopamine neuron-specific effects uncovered by the TRAP approach that may yield insight into the neurodegenerative process. To our knowledge, this is the first study to profile the effects of LRRK2 G2019S specifically on DA neuron gene expression in vivo. Beyond providing a set of differentially expressed gene candidates relevant to LRRK2, we demonstrate the effective use of TRAP to perform high-resolution assessment of dopamine neuron gene expression for the study of PD.

Analysis of cell-type-specific chromatin modifications and gene expression in Drosophila neurons that direct reproductive behavior

Examining the role of chromatin modifications and gene expression in neurons is critical for understanding how the potential for behaviors are established and maintained. We investigate this question by examining Drosophila melanogaster fru P1 neurons that underlie reproductive behaviors in both sexes. We developed a method to purify cell-type-specific chromatin (Chromatag), using a tagged histone H2B variant that is expressed using the versatile Gal4/UAS gene expression system. Here, we use Chromatag to evaluate five chromatin modifications, at three life stages in both sexes. We find substantial changes in chromatin modification profiles across development and fewer differences between males and females. Additionally, we find chromatin modifications that persist in different sets of genes from pupal to adult stages, which may point to genes important for cell fate determination in fru P1 neurons. We generated cell-type-specific RNA-seq data sets, using translating ribosome affinity purification (TRAP). We identify actively translated genes in fru P1 neurons, revealing novel stage- and sex-differences in gene expression. We also find chromatin modification enrichment patterns that are associated with gene expression. Next, we use the chromatin modification data to identify cell-type-specific super-enhancer-containing genes. We show that genes with super-enhancers in fru P1 neurons differ across development and between the sexes. We validated that a set of genes are expressed in fru P1 neurons, which were chosen based on having a super-enhancer and TRAP-enriched expression in fru P1 neurons.

Investigation of Drosophila fruitless neurons that express Dpr/DIP cell adhesion molecules

Drosophila reproductive behaviors are directed by fruitless neurons. A reanalysis of genomic studies shows that genes encoding dpr and DIP immunoglobulin superfamily (IgSF) members are expressed in fru P1 neurons. We find that each fru P1 and dpr/DIP (fru P1 ∩ dpr/DIP) overlapping expression pattern is similar in both sexes, but there are dimorphisms in neuronal morphology and cell number. Behavioral studies of fru P1 ∩ dpr/DIP perturbation genotypes indicate that the mushroom body functions together with the lateral protocerebral complex to direct courtship behavior. A single-cell RNA-seq analysis of fru P1 neurons shows that many DIPs have high expression in a small set of neurons, whereas the dprs are often expressed in a larger set of neurons at intermediate levels, with a myriad of dpr/DIP expression combinations. Functionally, we find that perturbations of sex hierarchy genes and of DIP-ε change the sex-specific morphologies of fru P1 ∩ DIP-α neurons.

Fruitless decommissions regulatory elements to implement cell-type-specific neuronal masculinization

In the fruit fly Drosophila melanogaster, male-specific splicing and translation of the Fruitless transcription factor (Fru^M) alters the presence, anatomy, and/or connectivity of >60 types of central brain neurons that interconnect to generate male-typical behaviors. While the indispensable function of Fru^M in sex-specific behavior has been understood for decades, the molecular mechanisms underlying its activity remain unknown. Here, we take a genome-wide, brain-wide approach to identifying regulatory elements whose activity depends on the presence of Fru^M. We identify 436 high-confidence genomic regions differentially accessible in male fruitless neurons, validate candidate regions as bona fide, differentially regulated enhancers, and describe the particular cell types in which these enhancers are active. We find that individual enhancers are not activated universally but are dedicated to specific fru⁺ cell types. Aside from fru itself, genes are not dedicated to or common across the fru circuit; rather, Fru^M appears to masculinize each cell type differently, by tweaking expression of the same effector genes used in other circuits. Finally, we find Fru^M motifs enriched among regulatory elements that are open in the female but closed in the male. Together, these results suggest that Fru^M acts cell-type-specifically to decommission regulatory elements in male fruitless neurons.

Courtship behavior in male Drosophila melanogaster is controlled by a well-defined neural circuit that is labeled by the male-specific transcription factor Fruitless (Fru^M). While Fru^M is known to change the number, anatomy and connectivity of neurons which comprise the circuit and has been suggested to repress the expression of a few gene targets, the mechanism of how Fru^M regulates genes across many different kinds of neurons is unknown. Using an approach to identify gene regulatory elements based on their chromatin accessibility states (ATAC-seq), we identified a large set of chromatin accessibility changes downstream of Fruitless. By examining the activity of these regulatory elements in vivo, we found that their activity was 1) sexually dimorphic and 2) specific to a single class of Fru^M neurons, suggesting that Fru^M acts on different chromatin targets in different neuron classes comprising the courtship circuit. Further, we found a known Fru^M-regulated enhancer of the Fru^M-repressed gene Lgr3 to have closed chromatin specifically in Fru^M neurons. Combined with an enrichment of Fru^M motifs in regions which are closed in Fru^M neurons, we present a mechanism where Fru^M directs the decommissioning of sex-shared regulatory elements to masculinize neurons in a cell-type specific manner.

Two Doublesex1 mutants revealed a tunable gene network underlying intersexuality in Daphnia magna

In recent years, the binary definition of sex is being challenged by repetitive reports about individuals with ambiguous sexual identity from various animal groups. This has created an urge to decode the molecular mechanism underlying sexual development. However, sexual ambiguities are extremely uncommon in nature, limiting their experimental value. Here, we report the establishment of a genetically modified clone of Daphnia magna from which intersex daphniids can be readily generated. By mutating the conserved central sex determining factor Doublesex1, body-wide feminization of male daphniid could be achieved. Comparative transcriptomic analysis also revealed a genetic network correlated with Doublesex1 activity which may account for the establishment of sexual identity in D. magna. We found that Dsx1 repressed genes related to growth and promoted genes related to signaling. We infer that different intersex phenotypes are the results of fluctuation in activity of these Dsx1 downstream factors. Our results demonstrated that the D. magna genome is capable of expressing sex in a continuous array, supporting the idea that sex is actually a spectrum.

The neurogenetics of sexually dimorphic behaviors from a postdevelopmental perspective

Most sexually reproducing animal species are characterized by two morphologically and behaviorally distinct sexes. The genetic, molecular and cellular processes that produce sexual dimorphisms are phylogenetically diverse, though in most cases they are thought to occur early in development. In some species, however, sexual dimorphisms are manifested after development is complete, suggesting the intriguing hypothesis that sex, more generally, might be considered a continuous trait that is influenced by both developmental and postdevelopmental processes. Here, we explore how biological sex is defined at the genetic, neuronal and behavioral levels, its effects on neuronal development and function, and how it might lead to sexually dimorphic behavioral traits in health and disease. We also propose a unifying framework for understanding neuronal and behavioral sexual dimorphisms in the context of both developmental and postdevelopmental, physiological timescales. Together, these two temporally separate processes might drive sex-specific neuronal functions in sexually mature adults, particularly as it pertains to behavior in health and disease.

Sex and the Single Fly: A Perspective on the Career of Bruce S. Baker

Bruce Baker, a preeminent Drosophila geneticist who made fundamental contributions to our understanding of the molecular genetic basis of sex differences, passed away July 1, 2018 at the age of 72. Members of Bruce's laboratory remember him as an intensely dedicated, rigorous, creative, deep-thinking, and fearless scientist. His trainees also remember his strong commitment to teaching students at every level. Bruce's career studying sex differences had three major epochs, where the laboratory was focused on: (1) sex determination and dosage compensation, (2) the development of sex-specific structures, and (3) the molecular genetic basis for sex differences in behavior. Several members of the Baker laboratory have come together to honor Bruce by highlighting some of the laboratory's major scientific contributions in these areas.

Dynamic changes in gene expression and alternative splicing mediate the response to acute alcohol exposure in Drosophila melanogaster

Environmental changes typically cause rapid gene expression responses in the exposed organisms, including changes in the representation of gene isoforms with different functions or properties. Identifying the genes that respond to environmental change, including in genotype-specific ways, is an important step in treating the undesirable physiological effects of stress, such as exposure to toxins or ethanol. Ethanol is a unique environmental stress in that chronic exposure results in permanent physiological changes and the development of alcohol use disorders. Drosophila is a classic model for deciphering the mechanisms of the response to alcohol exposure, as it meets the criteria for the development of alcohol use disorders, and has similar physiological underpinnings with vertebrates. Because many studies on the response to ethanol have relied on a priori candidate genes, broad surveys of gene expression and splicing are required and have been investigated here. Further, we expose Drosophila to ethanol in an environment that is genetically, socially, and ecologically relevant. Both expression and splicing differences, inasmuch as they can be decomposed, contribute to the response to ethanol in Drosophila melanogaster. However, we find that while D. melanogaster responds to ethanol, there is very little genetic variation in how it responds to ethanol. In addition, the response to alcohol over time is dynamic, suggesting that incorporating time into studies on the response to the environment is important.

Event Analysis: Using Transcript Events To Improve Estimates of Abundance in RNA-seq Data

Alternative splicing leverages genomic content by allowing the synthesis of multiple transcripts and, by implication, protein isoforms, from a single gene. However, estimating the abundance of transcripts produced in a given tissue from short sequencing reads is difficult and can result in both the construction of transcripts that do not exist, and the failure to identify true transcripts. An alternative approach is to catalog the events that make up isoforms (splice junctions and exons). We present here the Event Analysis (EA) approach, where we project transcripts onto the genome and identify overlapping/unique regions and junctions. In addition, all possible logical junctions are assembled into a catalog. Transcripts are filtered before quantitation based on simple measures: the proportion of the events detected, and the coverage. We find that mapping to a junction catalog is more efficient at detecting novel junctions than mapping in a splice aware manner. We identify 99.8% of true transcripts while iReckon identifies 82% of the true transcripts and creates more transcripts not included in the simulation than were initially used in the simulation. Using PacBio Iso-seq data from a mouse neural progenitor cell model, EA detects 60% of the novel junctions that are combinations of existing exons while only 43% are detected by STAR. EA further detects ∼5,000 annotated junctions missed by STAR. Filtering transcripts based on the proportion of the transcript detected and the number of reads on average supporting that transcript captures 95% of the PacBio transcriptome. Filtering the reference transcriptome before quantitation, results in is a more stable estimate of isoform abundance, with improved correlation between replicates. This was particularly evident when EA is applied to an RNA-seq study of type 1 diabetes (T1D), where the coefficient of variation among subjects (n = 81) in the transcript abundance estimates was substantially reduced compared to the estimation using the full reference. EA focuses on individual transcriptional events. These events can be quantitate and analyzed directly or used to identify the probable set of expressed transcripts. Simple rules based on detected events and coverage used in filtering result in a dramatic improvement in isoform estimation without the use of ancillary data (e.g., ChIP, long reads) that may not be available for many studies.