Sex-specific splicing in Drosophila: widespread occurrence, tissue specificity and evolutionary conservation

Tissue-specific metabolomic signatures for a doublesex model of reduced sexual dimorphism

Sex has a major effect on the metabolome. However, we do not yet understand the degree to which these quantitative sex differences in metabolism are associated with anatomical dimorphism and modulated by sex-specific tissues. In the fruit fly, Drosophila melanogaster, knocking out the doublesex (dsx) gene gives rise to adults with intermediate sex characteristics. Here we sought to determine the degree to which this key node in sexual development leads to sex differences in the fly metabolome. We measured 91 metabolites across head, thorax and abdomen in Drosophila, comparing the differences between distinctly sex-dimorphic flies with those of reduced sexual dimorphism: dsx null flies. Notably, in the reduced dimorphism flies, we observed a sex difference in only 1 of 91 metabolites, kynurenate, whereas 51% of metabolites (46/91) were significantly different between wildtype XX and XY flies in at least one tissue, suggesting that dsx plays a major role in sex differences in fly metabolism. Kynurenate was consistently higher in XX flies in both the presence and absence of functioning dsx. We observed tissue-specific consequences of knocking out dsx. Metabolites affected by sex were significantly enriched in branched chain amino acid metabolism and the mTOR pathway. This highlights the importance of considering variation in genes that cause anatomical sexual dimorphism when analyzing sex differences in metabolic profiles and interpreting their biological significance.

Identification of antennal alternative splicing by combining genome and full-length transcriptome analysis in Bactrocera dorsalis

Alternative splicing is an essential post-transcriptional regulatory mechanism that diversifies gene function by generating multiple protein isoforms from a single gene and act as a crucial role in insect environmental adaptation. Olfaction, a key sense for insect adaptation, relies heavily on the antennae, which are the primary olfactory organs expressing most of the olfactory genes. Despite the extensive annotation of olfactory genes within insect antennal tissues facilitated by high-throughput sequencing technology advancements, systematic analyses of alternative splicing are still relatively less. In this study, we focused on the oriental fruit fly (Bactrocera dorsalis), a significant pest of fruit crops. We performed a detailed analysis of alternative splicing in its antennae by utilizing the full-length transcriptome of its antennal tissue and the insect's genome. The results revealed 8600 non-redundant full-length transcripts identified in the oriental fruit fly antennal full-length transcriptome, spanning 4,145 gene loci. Over 40% of these loci exhibited multiple isoforms. Among these, 161 genes showed sex-biased isoform switching, involving seven different types of alternative splicing. Notably, events involving alternative transcription start sites (ATSS) and alternative transcription termination sites (ATTS) were the most common. Of all the genes undergoing ATSS and ATTS alternative splicing between male and female, 32 genes were alternatively spliced in protein coding regions, potentially affecting protein function. These genes were categorized based on the length of the sex-biased isoforms, with the highest difference in isoform fraction (dIF) associated with the ATSS type, including genes such as BdorABCA13, BdorCAT2, and BdorTSN3. Additionally, transcription factor binding sites for doublesex were identified upstream of both BdorABCA13 and BdorCAT2. Besides being expressed in the antennal tissues, BdorABCA13 and BdorCAT2 are also expressed in the mouthparts, legs, and genitalia of both female and male adults, suggesting their functional diversity. This study reveals alternative splicing events in the antennae of Bactrophora dorsalis from two aspects: odorant receptor genes and other types of genes expressed in the antennae. This study not only provides a research foundation for understanding the regulation of gene function by alternative splicing in the oriental fruit fly but also offers new insights for utilizing olfaction-based behavioral manipulation techniques to manage this pest.

Nucleotide-level distance metrics to quantify alternative splicing implemented in TranD

Advances in affordable transcriptome sequencing combined with better exon and gene prediction has motivated many to compare transcription across the tree of life. We develop a mathematical framework to calculate complexity and compare transcript models. Structural features, i.e. intron retention (IR), donor/acceptor site variation, alternative exon cassettes, alternative 5'/3' UTRs, are compared and the distance between transcript models is calculated with nucleotide level precision. All metrics are implemented in a PyPi package, TranD and output can be used to summarize splicing patterns for a transcriptome (1GTF) and between transcriptomes (2GTF). TranD output enables quantitative comparisons between: annotations augmented by empirical RNA-seq data and the original transcript models; transcript model prediction tools for longread RNA-seq (e.g. FLAIR versus Isoseq3); alternate annotations for a species (e.g. RefSeq vs Ensembl); and between closely related species. In C. elegans, Z. mays, D. melanogaster, D. simulans and H. sapiens, alternative exons were observed more frequently in combination with an alternative donor/acceptor than alone. Transcript models in RefSeq and Ensembl are linked and both have unique transcript models with empirical support. D. melanogaster and D. simulans, share many transcript models and long-read RNAseq data suggests that both species are under-annotated. We recommend combined references.

Shifts in isoform usage underlie transcriptional differences in regulatory T cells in type 1 diabetes

Genome-wide association studies have identified numerous loci with allelic associations to Type 1 Diabetes (T1D) risk. Most disease-associated variants are enriched in regulatory sequences active in lymphoid cell types, suggesting that lymphocyte gene expression is altered in T1D. Here we assay gene expression between T1D cases and healthy controls in two autoimmunity-relevant lymphocyte cell types, memory CD4+/CD25+ regulatory T cells (Treg) and memory CD4+/CD25- T cells, using a splicing event-based approach to characterize tissue-specific transcriptomes. Limited differences in isoform usage between T1D cases and controls are observed in memory CD4+/CD25- T-cells. In Tregs, 402 genes demonstrate differences in isoform usage between cases and controls, particularly RNA recognition and splicing factor genes. Many of these genes are regulated by the variable inclusion of exons that can trigger nonsense mediated decay. Our results suggest that dysregulation of gene expression, through shifts in alternative splicing in Tregs, contributes to T1D pathophysiology.

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..

Sex differences in gene expression and alternative splicing in the Chinese horseshoe bat

Sexually dimorphic traits are common in sexually reproducing organisms and can be encoded by differential gene regulation between males and females. Although alternative splicing is common mechanism in generating transcriptional diversity, its role in generating sex differences relative to differential gene expression is less clear. Here, we investigate the relative roles of differential gene expression and alternative splicing between male and female the horseshoe bat species, Rhinolophus sinicus. Horseshoe bats are an excellent model to study acoustic differences between sexes. Using RNA-seq analyses of two somatic tissues (brain and liver) from males and females of the same population, we identified 3,471 and 2,208 differentially expressed genes between the sexes (DEGs) in the brain and liver, respectively. DEGs were enriched with functional categories associated with physiological difference of the sexes (e.g.,gamete generation and energy production for reproduction in females). In addition, we also detected many differentially spliced genes between the sexes (DSGs, 2,231 and 1,027 in the brain and liver, respectively) which were mainly involved in regulation of RNA splicing and mRNA metabolic process. Interestingly, we found a significant enrichment of DEGs on the X chromosome, but not for DSGs. As for the extent of overlap between the two sets of genes, more than expected overlap of DEGs and DSGs was observed in the brain but not in the liver. This suggests that more complex tissues, such as the brain, may require the intricate and simultaneous interplay of both differential gene expression and splicing of genes to govern sex-specific functions. Overall, our results support that variation in gene expression and alternative splicing are important and complementary mechanisms governing sex differences.

Multi-tissue RNA-Seq Analysis and Long-read-based Genome Assembly Reveal Complex Sex-specific Gene Regulation and Molecular Evolution in the Manila Clam

The molecular factors and gene regulation involved in sex determination and gonad differentiation in bivalve molluscs are unknown. It has been suggested that doubly uniparental inheritance (DUI) of mitochondria may be involved in these processes in species such as the ubiquitous and commercially relevant Manila clam, Ruditapes philippinarum. We present the first long-read-based de novo genome assembly of a Manila clam, and a RNA-Seq multi-tissue analysis of 15 females and 15 males. The highly contiguous genome assembly was used as reference to investigate gene expression, alternative splicing, sequence evolution, tissue-specific co-expression networks, and sexual contrasting SNPs. Differential expression (DE) and differential splicing (DS) analyses revealed sex-specific transcriptional regulation in gonads, but not in somatic tissues. Co-expression networks revealed complex gene regulation in gonads, and genes in gonad-associated modules showed high tissue specificity. However, male gonad-associated modules showed contrasting patterns of sequence evolution and tissue specificity. One gene set was related to the structural organization of male gametes and presented slow sequence evolution but high pleiotropy, whereas another gene set was enriched in reproduction-related processes and characterized by fast sequence evolution and tissue specificity. Sexual contrasting SNPs were found in genes overrepresented in mitochondrial-related functions, providing new candidates for investigating the relationship between mitochondria and sex in DUI species. Together, these results increase our understanding of the role of DE, DS, and sequence evolution of sex-specific genes in an understudied taxon. We also provide resourceful genomic data for studies regarding sex diagnosis and breeding in bivalves.

Long-read RNA sequencing reveals widespread sex-specific alternative splicing in threespine stickleback fish

Alternate isoforms are important contributors to phenotypic diversity across eukaryotes. Although short-read RNA-sequencing has increased our understanding of isoform diversity, it is challenging to accurately detect full-length transcripts, preventing the identification of many alternate isoforms. Long-read sequencing technologies have made it possible to sequence full-length alternative transcripts, accurately characterizing alternative splicing events, alternate transcription start and end sites, and differences in UTR regions. Here, we use Pacific Biosciences (PacBio) long-read RNA-sequencing (Iso-Seq) to examine the transcriptomes of five organs in threespine stickleback fish (Gasterosteus aculeatus), a widely used genetic model species. The threespine stickleback fish has a refined genome assembly in which gene annotations are based on short-read RNA sequencing and predictions from coding sequence of other species. This suggests some of the existing annotations may be inaccurate or alternative transcripts may not be fully characterized. Using Iso-Seq we detected thousands of novel isoforms, indicating many isoforms are absent in the current Ensembl gene annotations. In addition, we refined many of the existing annotations within the genome. We noted many improperly positioned transcription start sites that were refined with long-read sequencing. The Iso-Seq-predicted transcription start sites were more accurate and verified through ATAC-seq. We also detected many alternative splicing events between sexes and across organs. We found a substantial number of genes in both somatic and gonadal samples that had sex-specific isoforms. Our study highlights the power of long-read sequencing to study the complexity of transcriptomes, greatly improving genomic resources for the threespine stickleback fish.

Gene expression and alternative splicing dynamics are perturbed in female head transcriptomes following heterospecific copulation

BACKGROUND

Despite the growing interest in the female side of copulatory interactions, the roles played by differential expression and alternative splicing mechanisms of pre-RNA on tissues outside of the reproductive tract have remained largely unknown. Here we addressed these questions in the context of con- vs heterospecific matings between Drosophila mojavensis and its sister species, D. arizonae. We analyzed transcriptional responses in female heads using an integrated investigation of genome-wide patterns of gene expression, including differential expression (DE), alternative splicing (AS) and intron retention (IR).

RESULTS

Our results indicated that early transcriptional responses were largely congruent between con- and heterospecific matings but are substantially perturbed over time. Conspecific matings induced functional pathways related to amino acid balance previously associated with the brain's physiology and female postmating behavior. Heterospecific matings often failed to activate regulation of some of these genes and induced expression of additional genes when compared with those of conspecifically-mated females. These mechanisms showed functional specializations with DE genes mostly linked to pathways of proteolysis and nutrient homeostasis, while AS genes were more related to photoreception and muscle assembly pathways. IR seems to play a more general role in DE regulation during the female postmating response.

CONCLUSIONS

We provide evidence showing that AS genes substantially perturbed by heterospecific matings in female heads evolve at slower evolutionary rates than the genome background. However, DE genes evolve at evolutionary rates similar, or even higher, than those of male reproductive genes, which highlights their potential role in sexual selection and the evolution of reproductive barriers.

Sex-Specific Transcriptome Differences in Human Adipose Mesenchymal Stem Cells

In humans, sexual dimorphism can manifest in many ways and it is widely studied in several knowledge fields. It is increasing the evidence that also cells differ according to sex, a correlation still little studied and poorly considered when cells are used in scientific research. Specifically, our interest is on the sex-related dimorphism on the human mesenchymal stem cells (hMSCs) transcriptome. A systematic meta-analysis of hMSC microarrays was performed by using the Transcriptome Mapper (TRAM) software. This bioinformatic tool was used to integrate and normalize datasets from multiple sources and allowed us to highlight chromosomal segments and genes differently expressed in hMSCs derived from adipose tissue (hADSCs) of male and female donors. Chromosomal segments and differentially expressed genes in male and female hADSCs resulted to be related to several processes as inflammation, adipogenic and neurogenic differentiation and cell communication. Obtained results lead us to hypothesize that the donor sex of hADSCs is a variable influencing a wide range of stem cell biologic processes. We believe that it should be considered in biologic research and stem cell therapy.

Modulated Expression of the Protein Kinase GSK3 in Motor and Dopaminergic Neurons Increases Female Lifespan in Drosophila melanogaster

Most eukaryotic genes express multiple transcripts and proteins, and a sophisticated gene expression strategy plays a crucial role in ensuring the cell-specificity of genetic information and the correctness of phenotypes. The Drosophila melanogaster gene shaggy encodes several isoforms of the conserved glycogen synthase kinase 3 (GSK3), which is vitally important for multiple biological processes. To characterize the phenotypic effects of differential shaggy expression, we explored how the multidirectional modulation of the expression of the main GSK3 isoform, Shaggy-PB, in different tissues and cells affects lifespan. To this end, we used lines with transgenic constructs that encode mutant variants of the protein. The effect of shaggy misexpression on lifespan depended on the direction of the presumed change in GSK3 activity and the type of tissue/cell. The modulation of GSK3 activity in motor and dopaminergic neurons improved female lifespan but caused seemingly negative changes in the structural (mitochondrial depletion; neuronal loss) and functional (perturbed locomotion) properties of the nervous system, indicating the importance of analyzing the relationship between lifespan and healthspan in invertebrate models. Our findings provide new insights into the molecular and cellular bases of lifespan extension, demonstrating that the fine-tuning of transcript-specific shaggy expression in individual groups of neurons is sufficient to provide a sex-specific increase in survival and slow aging.

Stress-mediated convergence of splicing landscapes in male and female rock doves

BACKGROUND

The process of alternative splicing provides a unique mechanism by which eukaryotes are able to produce numerous protein products from the same gene. Heightened variability in the proteome has been thought to potentiate increased behavioral complexity and response flexibility to environmental stimuli, thus contributing to more refined traits on which natural and sexual selection can act. While it has been long known that various forms of environmental stress can negatively affect sexual behavior and reproduction, we know little of how stress can affect the alternative splicing associated with these events, and less still about how splicing may differ between sexes. Using the model of the rock dove (Columba livia), our team previously uncovered sexual dimorphism in the basal and stress-responsive gene transcription of a biological system necessary for facilitating sexual behavior and reproduction, the hypothalamic-pituitary-gonadal (HPG) axis. In this study, we delve further into understanding the mechanistic underpinnings of how changes in the environment can affect reproduction by testing the alternative splicing response of the HPG axis to an external stressor in both sexes.

RESULTS

This study reveals dramatic baseline differences in HPG alternative splicing between males and females. However, after subjecting subjects to a restraint stress paradigm, we found a significant reduction in these differences between the sexes. In both stress and control treatments, we identified a higher incidence of splicing activity in the pituitary in both sexes as compared to other tissues. Of these splicing events, the core exon event is the most abundant form of splicing and more frequently occurs in the coding regions of the gene. Overall, we observed less splicing activity in the 3'UTR (untranslated region) end of transcripts than the 5'UTR or coding regions.

CONCLUSIONS

Our results provide vital new insight into sex-specific aspects of the stress response on the HPG axis at an unprecedented proximate level. Males and females uniquely respond to stress, yet exhibit splicing patterns suggesting a convergent, optimal splicing landscape for stress response. This information has the potential to inform evolutionary theory as well as the development of highly-specific drug targets for stress-induced reproductive dysfunction.

Rapid sex-specific adaptation to high temperature in Drosophila

The pervasive occurrence of sexual dimorphism demonstrates different adaptive strategies of males and females. While different reproductive strategies of the two sexes are well-characterized, very little is known about differential functional requirements of males and females in their natural habitats. Here, we study the impact environmental change on the selection response in both sexes. Exposing replicated Drosophila populations to a novel temperature regime, we demonstrate sex-specific changes in gene expression, metabolic and behavioral phenotypes in less than 100 generations. This indicates not only different functional requirements of both sexes in the new environment but also rapid sex-specific adaptation. Supported by computer simulations we propose that altered sex-biased gene regulation from standing genetic variation, rather than new mutations, is the driver of rapid sex-specific adaptation. Our discovery of environmentally driven divergent functional requirements of males and females has important implications-possibly even for gender aware medical treatments.

Sex biased expression and co-expression networks in development, using the hymenopteran Nasonia vitripennis

Sexual dimorphism requires regulation of gene expression in developing organisms. These developmental differences are caused by differential expression of genes and isoforms. The effect of expressing a gene is also influenced by which other genes are simultaneously expressed (functional interactions). However, few studies have described how these processes change across development. We compare the dynamics of differential expression, isoform switching and functional interactions in the sexual development of the model parasitoid wasp Nasonia vitripennis, a system that permits genome wide analysis of sex bias from early embryos to adults. We find relatively little sex-bias in embryos and larvae at the gene level, but several sub-networks show sex-biased functional interactions in early developmental stages. These networks provide new candidates for hymenopteran sex determination, including histone modification. In contrast, sex-bias in pupae and adults is driven by the differential expression of genes. We observe sex-biased isoform switching consistently across development, but mostly in genes that are already differentially expressed. Finally, we discover that sex-biased networks are enriched by genes specific to the Nasonia clade, and that those genes possess the topological properties of key regulators. These findings suggest that regulators in sex-biased networks evolve more rapidly than regulators of other developmental networks.

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.

Building a new research framework for social evolution: intralocus caste antagonism

The breeding and non‐breeding ‘castes’ of eusocial insects provide a striking example of role‐specific selection, where each caste maximises fitness through different morphological, behavioural and physiological trait values. Typically, queens are long‐lived egg‐layers, while workers are short‐lived, largely sterile foragers. Remarkably, the two castes are nevertheless produced by the same genome. The existence of inter‐caste genetic correlations is a neglected consequence of this shared genome, potentially hindering the evolution of caste dimorphism: alleles that increase the productivity of queens may decrease the productivity of workers and vice versa, such that each caste is prevented from reaching optimal trait values. A likely consequence of this ‘intralocus caste antagonism’ should be the maintenance of genetic variation for fitness and maladaptation within castes (termed ‘caste load’), analogous to the result of intralocus sexual antagonism. The aim of this review is to create a research framework for understanding caste antagonism, drawing in part upon conceptual similarities with sexual antagonism. By reviewing both the social insect and sexual antagonism literature, we highlight the current empirical evidence for caste antagonism, discuss social systems of interest, how antagonism might be resolved, and challenges for future research. We also introduce the idea that sexual and caste antagonism could interact, creating a three‐way antagonism over gene expression. This includes unpacking the implications of haplodiploidy for the outcome of this complex interaction.

The Y Chromosome Modulates Splicing and Sex-Biased Intron Retention Rates in Drosophila

The Drosophila Y chromosome is a 40-Mb segment of mostly repetitive DNA; it harbors a handful of protein-coding genes and a disproportionate amount of satellite repeats, transposable elements, and multicopy DNA arrays. Intron retention (IR) is a type of alternative splicing (AS) event by which one or more introns remain within the mature transcript. IR recently emerged as a deliberate cellular mechanism to modulate gene expression levels and has been implicated in multiple biological processes. However, the extent of sex differences in IR and the contribution of the Y chromosome to the modulation of AS and IR rates has not been addressed. Here we showed pervasive IR in the fruit fly Drosophila melanogaster with thousands of novel IR events, hundreds of which displayed extensive sex bias. The data also revealed an unsuspected role for the Y chromosome in the modulation of AS and IR. The majority of sex-biased IR events introduced premature termination codons and the magnitude of sex bias was associated with gene expression differences between the sexes. Surprisingly, an extra Y chromosome in males (X^YY genotype) or the presence of a Y chromosome in females (X^XY genotype) significantly modulated IR and recapitulated natural differences in IR between the sexes. Our results highlight the significance of sex-biased IR in tuning sex differences and the role of the Y chromosome as a source of variable IR rates between the sexes. Modulation of splicing and IR rates across the genome represent new and unexpected outcomes of the Drosophila Y chromosome.

Mating Changes Sexually Dimorphic Gene Expression in the Seed Beetle Callosobruchus maculatus

Sexually dimorphic phenotypes arise largely from sex-specific gene expression, which has mainly been characterized in sexually naïve adults. However, we expect sexual dimorphism in transcription to be dynamic and dependent on factors such as reproductive status. Mating induces many behavioral and physiological changes distinct to each sex and is therefore expected to activate regulatory changes in many sex-biased genes. Here, we first characterized sexual dimorphism in gene expression in Callosobruchus maculatus seed beetles. We then examined how females and males respond to mating and how it affects sex-biased expression, both in sex-limited (abdomen) and sex-shared (head and thorax) tissues. Mating responses were largely sex-specific and, as expected, females showed more genes responding compared with males (∼2,000 vs. ∼300 genes in the abdomen, ∼500 vs. ∼400 in the head and thorax, respectively). Of the sex-biased genes present in virgins, 16% (1,041 genes) in the abdomen and 17% (243 genes) in the head and thorax altered their relative expression between the sexes as a result of mating. Sex-bias status changed in 2% of the genes in the abdomen and 4% in the head and thorax following mating. Mating responses involved de-feminization of females and, to a lesser extent, de-masculinization of males relative to their virgin state: mating decreased rather than increased dimorphic expression of sex-biased genes. The fact that regulatory changes of both types of sex-biased genes occurred in both sexes suggests that male- and female-specific selection is not restricted to male- and female-biased genes, respectively, as is sometimes assumed.

The Goddard and Saturn Genes Are Essential for Drosophila Male Fertility and May Have Arisen De Novo

New genes arise through a variety of mechanisms, including the duplication of existing genes and the de novo birth of genes from noncoding DNA sequences. While there are numerous examples of duplicated genes with important functional roles, the functions of de novo genes remain largely unexplored. Many newly evolved genes are expressed in the male reproductive tract, suggesting that these evolutionary innovations may provide advantages to males experiencing sexual selection. Using testis-specific RNA interference, we screened 11 putative de novo genes in Drosophila melanogaster for effects on male fertility and identified two, goddard and saturn, that are essential for spermatogenesis and sperm function. Goddard knockdown (KD) males fail to produce mature sperm, while saturn KD males produce few sperm, and these function inefficiently once transferred to females. Consistent with a de novo origin, both genes are identifiable only in Drosophila and are predicted to encode proteins with no sequence similarity to any annotated protein. However, since high levels of divergence prevented the unambiguous identification of the noncoding sequences from which each gene arose, we consider goddard and saturn to be putative de novo genes. Within Drosophila, both genes have been lost in certain lineages, but show conserved, male-specific patterns of expression in the species in which they are found. Goddard is consistently found in single-copy and evolves under purifying selection. In contrast, saturn has diversified through gene duplication and positive selection. These data suggest that de novo genes can acquire essential roles in male reproduction.

Two spliced isoforms of the sex-determination gene fruitless in the Chinese mitten crab Eriocheir sinensis

The fruitless (fru) gene plays an important role in the sex determination pathway and courtship behavior of Drosophila melanogaster. In the present study, two fru isoforms (Esfru1 and Esfru2) were identified from the Chinese mitten crab Eriocheir sinensis. Sequence analysis showed that Esfru1 and Esfru2 were encoded by the same genomic locus and generated by alternative splicing of pre-mRNA. Esfru1 had all introns completely spliced out, while Esfru2 had a longer exon1 with an additional 78bp sequence. They both contained a conserved BTB domain which was also found in D. melanogaster Fru isoforms. Analysis on temporal expression profiles of Esfru1-2 (a common region of Esfru1 and Esfru2) and Esfru2 (a specific region of Esfru2) showed that they expressed in similar patterns during embryonic development but in different patterns during larval development. The expression of Esfru1-2 decreased gradually from zoea III stage, however, Esfru2 increased from zoea IV stage and reached the peak at megalopa stage. Expression distribution in tissues and in situ hybridization analysis revealed that they showed sexually dimorphic expression in gonads, hepatopancreas and brains. Esfru1-2 showed significantly higher expression in female gonads and hepatopancreas than in males, but highly expressed in male brains than that in females. Interestingly, Esfru2 was displayed in a male-specific manner in each tissue, especially in gonads and hepatopancreas. Our results indicate that Esfru1 might be involved in both sexual brain neuronal structure development and related to female-specific character development. Esfru2 may participate in male-specific character development. This is the first report that characterizes two spliced variants of fru in crustaceans and provides basic information for further functional studies of the crab sex-determination mechanism.

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.

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

Male and female reproductive behaviors in Drosophila melanogaster are vastly different, but neurons that express sex-specifically spliced fruitless transcripts (fru P1) underlie these behaviors in both sexes. How this set of neurons can generate such different behaviors between the two sexes is an unresolved question. A particular challenge is that fru P1-expressing neurons comprise only 2-5% of the adult nervous system, and so studies of adult head tissue or whole brain may not reveal crucial differences. Translating Ribosome Affinity Purification (TRAP) identifies the actively translated pool of mRNAs from fru P1-expressing neurons, allowing a sensitive, cell-type-specific assay. We find four times more male-biased than female-biased genes in TRAP mRNAs from fru P1-expressing neurons. This suggests a potential mechanism to generate dimorphism in behavior. The male-biased genes may direct male behaviors by establishing cell fate in a similar context of gene expression observed in females. These results suggest a possible global mechanism for how distinct behaviors can arise from a shared set of neurons.

Buffering of Genetic Regulatory Networks in Drosophila melanogaster

Regulatory variation in gene expression can be described by cis- and trans-genetic components. Here we used RNA-seq data from a population panel of Drosophila melanogaster test crosses to compare allelic imbalance (AI) in female head tissue between mated and virgin flies, an environmental change known to affect transcription. Indeed, 3048 exons (1610 genes) are differentially expressed in this study. A Bayesian model for AI, with an intersection test, controls type I error. There are ∼200 genes with AI exclusively in mated or virgin flies, indicating an environmental component of expression regulation. On average 34% of genes within a cross and 54% of all genes show evidence for genetic regulation of transcription. Nearly all differentially regulated genes are affected in cis, with an average of 63% of expression variation explained by the cis-effects. Trans-effects explain 8% of the variance in AI on average and the interaction between cis and trans explains an average of 11% of the total variance in AI. In both environments cis- and trans-effects are compensatory in their overall effect, with a negative association between cis- and trans-effects in 85% of the exons examined. We hypothesize that the gene expression level perturbed by cis-regulatory mutations is compensated through trans-regulatory mechanisms, e.g., trans and cis by trans-factors buffering cis-mutations. In addition, when AI is detected in both environments, cis-mated, cis-virgin, and trans-mated-trans-virgin estimates are highly concordant with 99% of all exons positively correlated with a median correlation of 0.83 for cis and 0.95 for trans We conclude that the gene regulatory networks (GRNs) are robust and that trans-buffering explains robustness.

Molecular cloning of the sex-related gene PSI in Bemisia tabaci and its alternative splicing properties

The P-element somatic inhibitor (PSI) is gene known to regulate the transcription of doublesex (dsx) when transformer (tra) is absent in Bombyx mori. In this study, we identified and characterized a PSI homolog in Bemisia tabaci (BtPSI). BtPSI cDNA had a total length of 5700 bp and contained a predicted open reading frame (ORF) of 2208 nucleotides encoding for 735 amino acids. Multiple sequence alignments of the common regions of PSI proteins from B. tabaci and five other insect species revealed a high degree of sequence conservation. BtPSI is expressed in all stages of B. tabaci development, and expression did not significantly differ between female and male adult. A total of 92 BtPSI isoforms (78 in female and 22 in male) were identified, and a marker indicating the female-specific form was found. These results increase the understanding of genes that may determine sex in B. tabaci and provide a foundation for research on the sex determination mechanism in this insect.

Sxl-Dependent, tra/tra2-Independent Alternative Splicing of the Drosophila melanogaster X-Linked Gene found in neurons

Somatic sexual determination and behavior in Drosophila melanogaster are under the control of a genetic cascade initiated by Sex lethal (Sxl). In the female soma, SXL RNA-binding protein regulates the splicing of transformer (tra) transcripts into a female-specific form. The RNA-binding protein TRA and its cofactor TRA2 function in concert in females, whereas SXL, TRA, and TRA2 are thought to not function in males. To better understand sex-specific regulation of gene expression, we analyzed male and female head transcriptome datasets for expression levels and splicing, quantifying sex-biased gene expression via RNA-Seq and qPCR. Our data uncouple the effects of Sxl and tra/tra2 in females in the-sex-biased alternative splicing of head transcripts from the X-linked locus found in neurons (fne), encoding a pan-neuronal RNA-binding protein of the ELAV family. We show that FNE protein levels are downregulated by Sxl in female heads, also independently of tra/tra2. We argue that this regulation may have important sexually dimorphic consequences for the regulation of nervous system development or function.

Genetic Architecture of Sexual Dimorphism in Humans

Males and females differ across a broad spectrum of morphological, physiological, and behavioral characters. In fact, sexually dimorphic traits typically contribute the largest component of phenotypic variance in most taxa that use sex to reproduce. However, we know very little about the mechanisms that maintain these dimorphic states and how these sexually dimorphic traits evolve. Here, we review our current knowledge of the underlying genetic basis of sexual dimorphism in humans. First, we briefly review the etiology of sex differences starting from sex determination's initial switch early in embryogenesis. We then survey recent sex-biased transcriptomic expression literature in order to provide additional insight into the landscape of sex-biased gene expression in both gonadal and non-gonadal tissues: from overall prevalence to tissue specificity to conservation across species. Finally, we discuss implications of sex-biased genetic architecture to human health and disease in light of the National Institute of Health's recently proposed initiative to promote study samples from both sexes.

Transcriptomic response to estrogen exposure in the male Zhikong scallop, Chlamys farreri

The transcriptomes of Zhikong scallop exposed to 17β-estradiol were determined using the Roche/454. A total of 51,997 unigenes, representing 45,030 contigs and 6967 singlets were obtained. And 14,028, 19,798 and 14,981 of these unigenes were annotated from the non-redundant nucleic acid database, non-redundant protein database and Swiss protein database, respectively. A total of 10,699 unigenes were further annotated to biological processes (9080), molecular functions (8692) and cellular components (7829) using the GO, and 8945 unigenes were mapped to biological pathways including the metabolism (2862) and genetic information processing (2263). Most importantly, 16,692 unigenes and 18,686 unigenes in testis, and 10,492 unigenes and 13,186 unigenes in digestive gland were up-regulated significantly after exposure to 50 and 500 ng E2/L; while 10,212 unigenes and 9409 unigenes in testis and 10,629 unigenes and 9463 unigenes in digestive gland were down-regulated. These valuable information provides insights into the mechanisms in invertebrate exposure to EDCs.

Subcellular localisations of the CPTI collection of YFP-tagged proteins in Drosophila embryos

A key challenge in the post-genomic area is to identify the function of the genes discovered, with many still uncharacterised in all metazoans. A first step is transcription pattern characterisation, for which we now have near whole-genome coverage in Drosophila. However, we have much more limited information about the expression and subcellular localisation of the corresponding proteins. The Cambridge Protein Trap Consortium generated, via piggyBac transposition, over 600 novel YFP-trap proteins tagging just under 400 Drosophila loci. Here, we characterise the subcellular localisations and expression patterns of these insertions, called the CPTI lines, in Drosophila embryos. We have systematically analysed subcellular localisations at cellularisation (stage 5) and recorded expression patterns at stage 5, at mid-embryogenesis (stage 11) and at late embryogenesis (stages 15-17). At stage 5, 31% of the nuclear lines (41) and 26% of the cytoplasmic lines (67) show discrete localisations that provide clues on the function of the protein and markers for organelles or regions, including nucleoli, the nuclear envelope, nuclear speckles, centrosomes, mitochondria, the endoplasmic reticulum, Golgi, lysosomes and peroxisomes. We characterised the membranous/cortical lines (102) throughout stage 5 to 10 during epithelial morphogenesis, documenting their apico-basal position and identifying those secreted in the extracellular space. We identified the tricellular vertices as a specialized membrane domain marked by the integral membrane protein Sidekick. Finally, we categorised the localisation of the membranous/cortical proteins during cytokinesis.

Male-specific phosphorylated SR proteins in adult flies of the Mediterranean fruitfly Ceratitis capitata

Alternative splicing is a widely used mechanism of gene regulation in sex determination pathways of Insects. In species from orders as distant as Diptera, Hymenoptera and Coleoptera, female differentiation relies on the activities of conserved splicing regulators, TRA and TRA-2, promoting female-specific expression of the global effector doublesex (dsx). Less understood is to what extent post-translational modifications of splicing regulators plays a role in this pathway. In Drosophila melanogaster phosphorylation of TRA, TRA-2 and the general RBP1 factor by the LAMMER kinase doa (darkener of apricot) is required for proper female sex determination. To explore whether this is a general feature of the pathway we examined sex-specific differences in phosphorylation levels of SR splicing factors in the dipteran species D. melanogaster, Ceratitis capitata (Medfly) and Musca domestica (Housefly). We found a distinct and reproducible pattern of male-specific phosphorylation on protein extracts enriched for SR proteins in C. capitata suggesting that differential phosphorylation may also contribute to the regulation of sex-specific splicing in the Medfly.

Phospho-regulated Drosophila adducin is a determinant of synaptic plasticity in a complex with Dlg and PIP2 at the larval neuromuscular junction

Adducin is a ubiquitously expressed actin- and spectrin-binding protein involved in cytoskeleton organization, and is regulated through phosphorylation of the myristoylated alanine-rich C-terminal kinase (MARCKS)-homology domain by protein kinase C (PKC). We have previously shown that the Drosophila adducin, Hu-li tai shao (Hts), plays a role in larval neuromuscular junction (NMJ) growth. Here, we find that the predominant isoforms of Hts at the NMJ contain the MARCKS-homology domain, which is important for interactions with Discs large (Dlg) and phosphatidylinositol 4,5-bisphosphate (PIP2). Through the use of Proximity Ligation Assay (PLA), we show that the adducin-like Hts isoforms are in complexes with Dlg and PIP2 at the NMJ. We provide evidence that Hts promotes the phosphorylation and delocalization of Dlg at the NMJ through regulation of the transcript distribution of the PAR-1 and CaMKII kinases in the muscle. We also show that Hts interactions with Dlg and PIP2 are impeded through phosphorylation of the MARCKS-homology domain. These results are further evidence that Hts is a signaling-responsive regulator of synaptic plasticity in Drosophila.

Sex-biased gene expression and sexual conflict throughout development

Sex-biased gene expression is likely to account for most sexually dimorphic traits because males and females share much of their genome. When fitness optima differ between sexes for a shared trait, sexual dimorphism can allow each sex to express their optimum trait phenotype, and in this way, the evolution of sex-biased gene expression is one mechanism that could help to resolve intralocus sexual conflict. Genome-wide patterns of sex-biased gene expression have been identified in a number of studies, which we review here. However, very little is known about how sex-biased gene expression relates to sex-specific fitness and about how sex-biased gene expression and conflict vary throughout development or across different genotypes, populations, and environments. We discuss the importance of these neglected areas of research and use data from a small-scale experiment on sex-specific expression of genes throughout development to highlight potentially interesting avenues for future research.

Diversity and dynamics of the Drosophila transcriptome

Animal transcriptomes are dynamic, with each cell type, tissue and organ system expressing an ensemble of transcript isoforms that give rise to substantial diversity. Here we have identified new genes, transcripts and proteins using poly(A)+ RNA sequencing from Drosophila melanogaster in cultured cell lines, dissected organ systems and under environmental perturbations. We found that a small set of mostly neural-specific genes has the potential to encode thousands of transcripts each through extensive alternative promoter usage and RNA splicing. The magnitudes of splicing changes are larger between tissues than between developmental stages, and most sex-specific splicing is gonad-specific. Gonads express hundreds of previously unknown coding and long non-coding RNAs (lncRNAs), some of which are antisense to protein-coding genes and produce short regulatory RNAs. Furthermore, previously identified pervasive intergenic transcription occurs primarily within newly identified introns. The fly transcriptome is substantially more complex than previously recognized, with this complexity arising from combinatorial usage of promoters, splice sites and polyadenylation sites.

Transcriptomic analysis of sexual differentiation in somatic tissues of Drosophila melanogaster: successes and caveats

The advent of high-throughput technologies to analyze RNA expression levels and transcript structure has brought renewed attention to the age-old question of what differentiates males from females. In Drosophila, the characterized somatic sex determination cascade includes proteins implicated in the regulation of pre-mRNA splicing as well as at least 2 transcription factors at its base. Both DNA microarrays and RNA-Seq have been applied in a number of studies to determine the identities, expression levels and structure of transcripts expressed differentially in the 2 sexes, with remarkably divergent results in the number, structure and identity of affected transcripts. We briefly summarize these reports and discuss the reasons for the apparent discrepancies based upon the different conditions used for sample preparation and data analysis.

Caste load and the evolution of reproductive skew

Reproductive skew theory seeks to explain how reproduction is divided among group members in animal societies. Existing theory is framed almost entirely in terms of selection, though nonadaptive processes must also play some role in the evolution of reproductive skew. Here I propose that a genetic correlation between helper fecundity and breeder fecundity may frequently constrain the evolution of reproductive skew. This constraint is part of a wider phenomenon that I term "caste load," which is defined as the decline in mean fitness caused by caste-specific selection pressures, that is, differential selection on breeding and nonbreeding individuals. I elaborate the caste load hypothesis using quantitative and population genetic arguments and individual-based simulations. Although selection can sometimes erode genetic correlations and resolve caste load, this may be constrained when mutations have similar pleiotropic effects on breeder and helper traits. I document evidence for caste load, identify putative genomic adaptations to it, and suggest future research directions. The models highlight the value of considering adaptation within the boundaries imposed by genetic architecture and incidentally reaffirm that monogamy promotes the evolutionary transition to eusociality.

Identification of male gametogenesis expressed genes from the scallop Nodipecten subnodosus by suppressive subtraction hybridization and pyrosequencing

Despite the great advances in sequencing technologies, genomic and transcriptomic information for marine non-model species with ecological, evolutionary, and economical interest is still scarce. In this work we aimed to identify genes expressed during spermatogenesis in the functional hermaphrodite scallop Nodipecten subnodosus (Mollusca: Bivalvia: Pectinidae), with the purpose of obtaining a panel of genes that would allow for the study of differentially transcribed genes between diploid and triploid scallops in the context of meiotic arrest and reproductive sterility. Because our aim was to isolate genes involved in meiosis and other testis maturation-related processes, we generated suppressive subtractive hybridization libraries of testis vs. inactive gonad. We obtained 352 and 177 ESTs by clone sequencing, and using pyrosequencing (454-Roche) we maximized the identified ESTs to 34,276 reads. A total of 1,153 genes from the testis library had a blastx hit and GO annotation, including genes specific for meiosis, spermatogenesis, sex-differentiation, and transposable elements. Some of the identified meiosis genes function in chromosome pairing (scp2, scp3), recombination and DNA repair (dmc1, rad51, ccnb1ip1/hei10), and meiotic checkpoints (rad1, hormad1, dtl/cdt2). Gene expression analyses in different gametogenic stages in both sexual regions of the gonad of meiosis genes confirmed that the expression was specific or increased towards the maturing testis. Spermatogenesis genes included known testis-specific ones (kelch-10, shippo1, adad1), with some of these known to be associated to sterility. Sex differentiation genes included one of the most conserved genes at the bottom of the sex-determination cascade (dmrt1). Transcript from transposable elements, reverse transcriptase, and transposases in this library evidenced that transposition is an active process during spermatogenesis in N. subnodosus. In relation to the inactive library, we identified 833 transcripts with functional annotation related to activation of the transcription and translation machinery, as well as to germline control and maintenance.

New levels of transcriptome complexity at upper thermal limits in wild Drosophila revealed by exon expression analysis

While the cellular heat-shock response has been a paradigm for studying the impact of thermal stress on RNA metabolism and gene expression, the genome-wide response to thermal stress and its connection to physiological stress resistance remain largely unexplored. Here, we address this issue using an array-based exon expression analysis to interrogate the transcriptome in recently established Drosophila melanogaster stocks during severe thermal stress and recovery. We first demonstrated the efficacy of exon-level analyses to reveal a level of thermally induced transcriptome complexity extending well beyond gene-level analyses. Next, we showed that the upper range of both the cellular and physiological thermal stress response profoundly affected message expression and processing in D. melanogaster, limiting expression to a small subset of transcripts, many that share features of known rapidly responding stress genes. As predicted from cellular heat-shock research, constitutive splicing was blocked in a set of novel genes; we did not detect changes to alternative splicing during heat stress, but rather induction of intronless isoforms of known heat-responsive genes. We observed transcriptome plasticity in the form of differential isoform expression during recovery from heat shock, mediated by multiple mechanisms including alternative transcription and alternative splicing. This affected genes involved in DNA regulation, immune response, and thermotolerance. These patterns highlight the complex nature of innate transcriptome responses under stress and potential for adaptive shifts through plasticity and evolved genetic responses at different hierarchical levels.

The scope and strength of sex-specific selection in genome evolution

Males and females share the vast majority of their genomes and yet are often subject to different, even conflicting, selection. Genomic and transcriptomic developments have made it possible to assess sex-specific selection at the molecular level, and it is clear that sex-specific selection shapes the evolutionary properties of several genomic characteristics, including transcription, post-transcriptional regulation, imprinting, genome structure and gene sequence. Sex-specific selection is strongly influenced by mating system, which also causes neutral evolutionary changes that affect different regions of the genome in different ways. Here, we synthesize theoretical and molecular work in order to provide a cohesive view of the role of sex-specific selection and mating system in genome evolution. We also highlight the need for a combined approach, incorporating both genomic data and experimental phenotypic studies, in order to understand precisely how sex-specific selection drives evolutionary change across the genome.

The shared genome is a pervasive constraint on the evolution of sex-biased gene expression

Males and females share most of their genomes, and differences between the sexes can therefore not evolve through sequence divergence in protein coding genes. Sexual dimorphism is instead restricted to occur through sex-specific expression and splicing of gene products. Evolution of sexual dimorphism through these mechanisms should, however, also be constrained when the sexes share the genetic architecture for regulation of gene expression. Despite these obstacles, sexual dimorphism is prevalent in the animal kingdom and commonly evolves rapidly. Here, we ask whether the genetic architecture of gene expression is plastic and easily molded by sex-specific selection, or if sexual dimorphism evolves rapidly despite pervasive genetic constraint. To address this question, we explore the relationship between the intersexual genetic correlation for gene expression (rMF), which captures how independently genes are regulated in the sexes, and the evolution of sex-biased gene expression. Using transcriptome data from Drosophila melanogaster, we find that most genes have a high rMF and that genes currently exposed to sexually antagonistic selection have a higher average rMF than other genes. We further show that genes with a high rMF have less pronounced sex-biased gene expression than genes with a low rMF within D. melanogaster and that the strength of the rMF in D. melanogaster predicts the degree to which the sex bias of a gene's expression has changed between D. melanogaster and six other species in the Drosophila genus. In sum, our results show that a shared genome constrains both short- and long-term evolution of sexual dimorphism.

Genomics of sex determination in Drosophila

Drosophilists have identified many, or perhaps most, of the key regulatory genes determining sex using classical genetics, however, regulatory genes must ultimately result in the deployment of the genome in a quantitative manner, replete with complex interactions with other regulatory pathways. In the last decade, genomics has provided a rich picture of the transcriptional profile of the sexes that underlies sexual dimorphism. The current challenge is linking transcriptional profiles with the regulatory genes. This will be a complex synthesis, but the prospects for progress are outstanding.

Alternative splicing and trans-splicing events revealed by analysis of the Bombyx mori transcriptome

Alternative splicing and trans-splicing events have not been systematically studied in the silkworm Bombyx mori. Here, the silkworm transcriptome was analyzed by RNA-seq. We identified 320 novel genes, modified 1140 gene models, and found thousands of alternative splicing and 58 trans-splicing events. Studies of three SR proteins show that both their alternative splicing patterns and mRNA products are conserved from insect to human, and one isoform of Srsf6 with a retained intron is expressed sex-specifically in silkworm gonads. Trans-splicing of mod(mdg4) in silkworm was experimentally confirmed. We identified integrations from a common 5'-gene with 46 newly identified alternative 3'-exons that are located on both DNA strands over a 500-kb region. Other trans-splicing events in B. mori were predicted by bioinformatic analysis, in which 12 events were confirmed by RT-PCR, six events were further validated by chimeric SNPs, and two events were confirmed by allele-specific RT-PCR in F(1) hybrids from distinct silkworm lines of JS and L10, indicating that trans-splicing is more widespread in insects than previously thought. Analysis of the B. mori transcriptome by RNA-seq provides valuable information of regulatory alternative splicing events. The conservation of splicing events across species and newly identified trans-splicing events suggest that B. mori is a good model for future studies.

The evolution of gene expression and the transcriptome-phenotype relationship

Changes in gene expression underlie the adaptive evolution in many complex phenotypes, and the recent increase in the availability of multi-species comparative transcriptome data has made it possible to scan whole transcriptomes for loci that have experienced adaptive changes in expression. However, despite the increase in data availability, current models of gene expression evolution often do not account for the complexities and inherent noise associated with transcriptome data. Additionally, in contrast to current models of gene sequence evolution, models of transcriptome evolution often lack the sophistication to effectively determine whether transcriptional differences between species or within a clade are the result of neutral or adaptive processes. In this review, we discuss the tools, methods and models that define our current understanding of the relationship between gene expression and complex phenotype evolution. Our goal is to summarize what we know about the evolution of global gene expression patterns underlying complex traits, as well to identify some of the questions that remain to be answered.

Partitioning transcript variation in Drosophila: abundance, isoforms, and alleles

Multilevel analysis of transcription is facilitated by a new array design that includes modules for assessment of differential expression, isoform usage, and allelic imbalance in Drosophila. The ∼2.5 million feature chip incorporates a large number of controls, and it contains 18,769 3' expression probe sets and 61,919 exon probe sets with probe sequences from Drosophila melanogaster and 60,118 SNP probe sets focused on Drosophila simulans. An experiment in D. simulans identified genes differentially expressed between males and females (34% in the 3' expression module; 32% in the exon module). These proportions are consistent with previous reports, and there was good agreement (κ = 0.63) between the modules. Alternative isoform usage between the sexes was identified for 164 genes. The SNP module was verified with resequencing data. Concordance between resequencing and the chip design was greater than 99%. The design also proved apt in separating alleles based upon hybridization intensity. Concordance between the highest hybridization signals and the expected alleles in the genotype was greater than 96%. Intriguingly, allelic imbalance was detected for 37% of 6579 probe sets examined that contained heterozygous SNP loci. The large number of probes and multiple probe sets per gene in the 3' expression and exon modules allows the array to be used in D. melanogaster and in closely related species. The SNP module can be used for allele specific expression and genotyping of D. simulans.

Gene duplication and the genome distribution of sex-biased genes

In species that have two sexes, a single genome encodes two morphs, as each sex can be thought of as a distinct morph. This means that the same set of genes are differentially expressed in the different sexes. Many questions emanate from this statement. What proportion of genes contributes to sexual dimorphism? How do they contribute to sexual dimorphism? How is sex-biased expression achieved? Which sex and what tissues contribute the most to sex-biased expression? Do sex-biased genes have the same evolutionary patterns as nonbiased genes? We review the current data on sex-biased expression in species with heteromorphic sex chromosomes and comment on the most important hypotheses suggested to explain the origin, evolution, and distribution patterns of sex-biased genes. In this perspective we emphasize how gene duplication serves as an important molecular mechanism to resolve genomic clashes and genetic conflicts by generating sex-biased genes, often sex-specific genes, and contributes greatly to the underlying genetic basis of sexual dimorphism.

Somatic sex-specific transcriptome differences in Drosophila revealed by whole transcriptome sequencing

Background

Understanding animal development and physiology at a molecular-biological level has been advanced by the ability to determine at high resolution the repertoire of mRNA molecules by whole transcriptome resequencing. This includes the ability to detect and quantify rare abundance transcripts and isoform-specific mRNA variants produced from a gene.

The sex hierarchy consists of a pre-mRNA splicing cascade that directs the production of sex-specific transcription factors that specify nearly all sexual dimorphism. We have used deep RNA sequencing to gain insight into how the Drosophila sex hierarchy generates somatic sex differences, by examining gene and transcript isoform expression differences between the sexes in adult head tissues.

Results

Here we find 1,381 genes that differ in overall expression levels and 1,370 isoform-specific transcripts that differ between males and females. Additionally, we find 512 genes not regulated downstream of transformer that are significantly more highly expressed in males than females. These 512 genes are enriched on the × chromosome and reside adjacent to dosage compensation complex entry sites, which taken together suggests that their residence on the × chromosome might be sufficient to confer male-biased expression. There are no transcription unit structural features, from a set of features, that are robustly significantly different in the genes with significant sex differences in the ratio of isoform-specific transcripts, as compared to random isoform-specific transcripts, suggesting that there is no single molecular mechanism that generates isoform-specific transcript differences between the sexes, even though the sex hierarchy is known to include three pre-mRNA splicing factors.

Conclusions

We identify thousands of genes that show sex-specific differences in overall gene expression levels, and identify hundreds of additional genes that have differences in the abundance of isoform-specific transcripts. No transcription unit structural feature was robustly enriched in the sex-differentially expressed transcript isoforms. Additionally, we found that many genes with male-biased expression were enriched on the × chromosome and reside adjacent to dosage compensation entry sites, suggesting that differences in sex chromosome composition contributes to dimorphism in gene expression. Taken together, this study provides new insight into the molecular underpinnings of sexual differentiation.

Microarrays, deep sequencing and the true measure of the transcriptome

Microarrays first made the analysis of the transcriptome possible, and have produced much important information. Today, however, researchers are increasingly turning to direct high-throughput sequencing - RNA-Seq - which has considerable advantages for examining transcriptome fine structure - for example in the detection of allele-specific expression and splice junctions. In this article, we discuss the relative merits of the two techniques, the inherent biases in each, and whether all of the vast body of array work needs to be revisited using the newer technology. We conclude that microarrays remain useful and accurate tools for measuring expression levels, and RNA-Seq complements and extends microarray measurements.

Hts, the Drosophila homologue of Adducin, physically interacts with the transmembrane receptor Golden goal to guide photoreceptor axons

Neurons steer their axons towards their proper targets during development. Molecularly, a number of guidance receptors have been identified. The transmembrane protein Golden goal (Gogo) was reported previously to guide photoreceptor (R) axons in the Drosophila visual system. Here, we show that Hts, the Drosophila homologue of Adducin, physically interacts with Gogo's cytoplasmic domain via its head-neck domain. hts null mutants show similar defects in R axon guidance as do gogo mutants. Rescue experiments suggest that the C-terminal tail but not the MARCKS homology domain of Hts is required. Overexpression of either gogo or hts causes abnormally thick swellings of R8 axons in the medulla, but if both are co-overexpressed, R8 axons appear normal and the amount of excessive Hts is reduced. Our results fit with a model where Gogo both positively and negatively regulates Hts that affects the Actin-Spectrin cytoskeleton in growth cone filopodia, thereby guiding R axons.

Distinct regulatory programs establish widespread sex-specific alternative splicing in Drosophila melanogaster

In Drosophila melanogaster, female-specific expression of Sex-lethal (SXL) and Transformer (TRA) proteins controls sex-specific alternative splicing and/or translation of a handful of regulatory genes responsible for sexual differentiation and behavior. Recent findings in 2009 by Telonis-Scott et al. document widespread sex-biased alternative splicing in fruitflies, including instances of tissue-restricted sex-specific splicing. Here we report results arguing that some of these novel sex-specific splicing events are regulated by mechanisms distinct from those established by female-specific expression of SXL and TRA. Bioinformatic analysis of SXL/TRA binding sites, experimental analysis of sex-specific splicing in S2 and Kc cells lines and of the effects of SXL knockdown in Kc cells indicate that SXL-dependent and SXL-independent regulatory mechanisms coexist within the same cell. Additional determinants of sex-specific splicing can be provided by sex-specific differences in the expression of RNA binding proteins, including Hrp40/Squid. We report that sex-specific alternative splicing of the gene hrp40/squid leads to sex-specific differences in the levels of this hnRNP protein. The significant overlap between sex-regulated alternative splicing changes and those induced by knockdown of hrp40/squid and the presence of related sequence motifs enriched near subsets of Hrp40/Squid-regulated and sex-regulated splice sites indicate that this protein contributes to sex-specific splicing regulation. A significant fraction of sex-specific splicing differences are absent in germline-less tudor mutant flies. Intriguingly, these include alternative splicing events that are differentially spliced in tissues distant from the germline. Collectively, our results reveal that distinct genetic programs control widespread sex-specific splicing in Drosophila melanogaster.

Gender aspects in skin diseases

Gender differences in medicine have been recognized in anatomy, physiology, as well as in epidemiology and manifestations of various diseases. With respect to skin disorders, males are generally more commonly afflicted with infectious diseases while women are more susceptible to psychosomatic problems, pigmentary disorders, certain hair diseases, and particularly autoimmune as well as allergic diseases. Significantly, more female sex-associated dermatoses can be identified than the male sex-associated dermatoses. Dermatoses in the genital area differ between men and women. Gender differences also exist in the occurrence and prognosis of certain skin malignancies. The mechanisms underlying gender differences in skin diseases remain largely unknown. Differences in the skin structure and physiology, effect of sex hormones, ethnic background, sociocultural behaviour and environmental factors may interact to exert the influences. A better understanding of gender differences in human health and diseases will allow the development of novel concepts for prevention, diagnosis and therapy of skin diseases.

Composite effects of polymorphisms near multiple regulatory elements create a major-effect QTL

Many agriculturally, evolutionarily, and medically important characters vary in a quantitative fashion. Unfortunately, the genes and sequence variants accounting for this variation remain largely unknown due to a variety of biological and technical challenges. Drosophila melanogaster contains high levels of sequence variation and low linkage disequilibrium, allowing us to dissect the effects of many causative variants within a single locus. Here, we take advantage of these features to identify and characterize the sequence polymorphisms that comprise major effect QTL alleles segregating at the bric-a-brac locus. We show that natural bric-a-brac alleles with large effects on cuticular pigmentation reflect a cumulative impact of polymorphisms that affect three functional regions: a promoter, a tissue-specific enhancer, and a Polycomb response element. Analysis of allele-specific expression at the bric-a-brac locus confirms that these polymorphisms modulate transcription at the cis-regulatory level. Our results establish that a single QTL can act through a confluence of multiple molecular mechanisms and that sequence variation in regions flanking experimentally validated functional elements can have significant quantitative effects on transcriptional activity and phenotype. These findings have important design and conceptual implications for basic and medical genomics.