Dimerization of doublesex is mediated by a cryptic ubiquitin-associated domain fold: implications for sex-specific gene regulation

Sexual dimorphic eyestalk transcriptome of kuruma prawn Marsupenaeus japonicus

Kuruma prawn (Marsupenaeus japonicus) is a benthic decapod crustacean that is widely distributed in the Indo-West Pacific region. It is one of the most important fishery resources in Japan, but its annual catches have declined sharply since the 1990s. To increase stocks, various approaches such as seed production and aquaculture were attempted. Since the demand for important fishery species, including kuruma prawn, is expected to increase worldwide in the future, there is a need to develop new technologies that will make aquaculture more efficient. Historically, the eyestalk endocrine organ is known to consist of the X-organ and sinus gland (XO/SG) complex that synthesizes and secrets various neuropeptide hormones that regulate growth, molt, sexual maturation, reproduction, and changes in body color. In the current study, eyestalk-derived neuropeptides were identified in the transcriptome. In addition, most orthologs of sex-determination genes were expressed in eyestalks. We identified two doublesex genes (MjapDsx1 and MjapDsx2) and found that MjapDsx1 showed male-biased expression in the eyestalk ganglion with no sex-specific splicing, unlike insect species. Therefore, this study will provide an opportunity to advance the research of neuropeptides and sex determination in the kuruma prawn.

When does the female bias arise? Insights from the sex determination cascade of a flea beetle with a strongly skewed sex ratio

Reproduction-manipulating bacteria like Wolbachia can shift sex ratios in insects towards females, but skewed sex ratios may also arise from genetic conflicts. The flea beetle Altica lythri harbors three main mtDNA strains that are coupled to three different Wolbachia infections. Depending on the mtDNA types, the females produce either offspring with a balanced sex ratio or exclusively daughters. To obtain markers that can monitor when sex bias arises in the beetle's ontogeny, we elucidated the sex determination cascade of A. lythri. We established a RT-PCR method based on length variants of dsx (doublesex) transcripts to determine the sex of morphologically indistinguishable eggs and larvae. In females of one mtDNA type (HT1/HT1*) known to produce only daughters, male offspring were already missing at the egg stage while for females of another type (HT2), the dsx splice variants revealed a balanced sex ratio among eggs and larvae. Our data suggest that the sex determination cascade in A. lythri is initiated by maternally transmitted female-specific tra (transformer) mRNA as primary signal. This tra mRNA seems to be involved in a positive feedback loop that maintains the production of the female splice variant, as known for female offspring in Tribolium castaneum. The translation of the maternally transmitted female tra mRNA must be inhibited in male offspring, but the underlying primary genetic signal remains to be identified. We discuss which differences between the mtDNA types can influence sex determination and lead to the skewed sex ratio of HT1.

Identification and characterization of the Doublesex gene and its mRNA isoforms in the brine shrimp Artemia franciscana

Doublesex (DSX) proteins are members of the Doublesex/mab-3-related (DMRT) protein family and play crucial roles in sex determination and differentiation among the animal kingdom. In the present study, we identified two Doublesex (Dsx)-like mRNA isoforms in the brine shrimp Artemia franciscana (Kellogg 1906), which are generated by the combination of alternative promoters, alternative splicing and alternative polyadenylation. The two transcripts exhibited sex-biased enrichment, which we termed AfrDsxM and AfrDsxF. They share a common region which encodes an identical N-terminal DNA-binding (DM) domain. RT-qPCR analyses showed that AfrDsxM is dominantly expressed in male Artemia while AfrDsxF is specifically expressed in females. Expression levels of both isoforms increased along with the developmental stages of their respective sexes. RNA interference with dsRNA showed that the knockdown of AfrDsxM in male larvae led to the appearance of female traits including an ovary-like structure in the original male reproductive system and an elevated expression of vitellogenin. However, silencing of AfrDsxF induced no clear phenotypic change in female Artemia. These results indicated that the male AfrDSXM may act as inhibiting regulator upon the default female developmental mode in Artemia. Furthermore, electrophoretic mobility shift assay analyses revealed that the unique DM domain of AfrDSXs can specifically bind to promoter segments of potential downstream target genes like AfrVtg. These data show that AfrDSXs play crucial roles in regulating sexual development in Artemia, and further provide insight into the evolution of sex determination/differentiation in sexual organisms.

The Sex-Specific Splicing of Doublesex in Brine Shrimp Artemia franciscana

The understanding of sex determination and differentiation in animals has recently made remarkable strides through the use of advanced research tools. At the gene level, the Mab-3-related transcription factor (Dmrt) gene family, which encodes for the typical DNA-binding doublesex/Mab-3 (DM) domain in their protein, is known for its contribution to sex determination and differentiation in insects. In this study, DNA-binding DM domain screening has identified eight transcripts from Artemia franciscana transcriptomic that encode proteins containing one conserved DNA-binding DM domain. The genome mapping confirmed that these eight transcripts are transcribed from six different loci on the A. franciscana genome assembly. One of those loci, the Af.dsx-4 locus, is closely related to Doublesex, a gene belonging to the Dmrt gene family. This locus could be transcribed into three alternative transcripts, namely Af.dsx⁴, Af.dsx^F and Af.dsx^M. While Af.dsx⁴ and Af.dsx^F could putatively be translated to form an identical Af.dsx^F protein of 186 aa long, Af.dsx^M translates for an Af.dsx^M protein of 289 aa long but shares a DNA-binding DM domain. Interestingly, Af.dsx^F and Af.dsx^M are confirmed as sex-specific transcripts, Af.dsx^F is only present in females, and Af.dsx^M is only present in male individuals. The results suggest that the sex-specific splicing mechanism of the doublesex described in insects is also present in A. franciscana. Af.dxs-4 locus can be used in further studies to clarify the sex determination pathways in A. fracnciscana.

Identification and Expression Analysis of Dsx and Its Positive Transcriptional Regulation of IAG in Black Tiger Shrimp (Penaeus monodon)

Doublesex (Dsx) is a polymorphic transcription factor of the DMRTs family, which is involved in male sex trait development and controls sexual dimorphism at different developmental stages in arthropods. However, the transcriptional regulation of the Dsx gene is largely unknown in decapods. In this study, we reported the cDNA sequence of PmDsx in Penaeus monodon, which encodes a 257 amino acid polypeptide. It shared many similarities with Dsx homologs and has a close relationship in the phylogeny of different species. We demonstrated that the expression of the male sex differentiation gene Dsx was predominantly expressed in the P. monodon testis, and that PmDsx dsRNA injection significantly decreased the expression of the insulin-like androgenic gland hormone (IAG) and male sex-determining gene while increasing the expression of the female sex-determining gene. We also identified a 5′-flanking region of PmIAG that had two potential cis-regulatory elements (CREs) for the PmDsx transcription. Further, the dual-luciferase reporter analysis and truncated mutagenesis revealed that PmDsx overexpression significantly promoted the transcriptional activity of the PmIAG promoter via a specific CRE. These results suggest that PmDsx is engaged in male reproductive development and positively regulates the transcription of the PmIAG by specifically binding upstream of the promoter of the PmIAG. It provides a theoretical basis for exploring the sexual regulation pathway and evolutionary dynamics of Dmrt family genes in P. monodon.

Transcriptome analysis of different life-history stages and screening of male-biased genes in Daphnia sinensis

Background

In the life history of Daphnia, the reproductive mode of parthenogenesis and sexual reproduction alternate in aquatic ecosystem, which are often affected by environmental and genetic factors. Recently, the sex-biased genes are of great significance for clarifying the origin and evolution of reproductive transformation and the molecular regulation mechanism of sex determination in Daphnia. Although some genes on reproductive transition of Daphnia had been researched, molecular mechanism on the maintenance of sexually dimorphic phenotypes of Daphnia are still not well known, including differentially expressed genes in different life-history stages.

Results

In this study, four life-history stages of Daphnia sinensis, juvenile female (JF), parthenogenetic female (PF), sexual female (SF) and male (M), were performed for transcriptome, and male-biased genes were screened. A total of 110437 transcripts were obtained and assembled into 22996 unigenes. In the four life-history stages (JF, PF, SF and M), the number of unique unigenes is respectively 2863, 445, 437 and 586, and the number of common unigenes is 9708. The differentially expressed genes (DEGs) between male and other three female stages (M vs JF, M vs PF and M vs SF) were 4570, 4358 and 2855, respectively. GO gene enrichment analysis showed that the up-regulated genes in male were mainly enriched in hydrolase activity and peptidase activity. Thirty-six genes in male were significantly higher expression than in the three female stages, including one Doublesex (Dsx) gene, one laminin gene, five trypsin genes and one serine protease genes, and one chitin synthase gene and two chitinase genes.

Conclusions

Our results showed that thirty-six candidate genes may be as the male-biased genes involving in the maintenance of sexually dimorphic phenotypes. This work will provide a reference for further exploring the functional genes related to sex differentiation in Daphnia species. Moreover, according to previous investigations, we thought that the expression level of functional genes may be related to the life-history stages of organisms, and may be also affected by different Daphnia species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08824-x.

Regulation of Doublesex1 Expression for Environmental Sex Determination in the Cladoceran Crustacean Daphnia

The cladoceran crustacean Daphnia produces only females by parthenogenesis in a healthy population. However, in response to environmental declines such as crowding and lack of foods, it produces eggs destined to become males that are genetically identical to females. During the development of the sexually committed eggs, DM domain-containing transcription factor Doublesex1 (Dsx1) orchestrates male trait formation globally both in somatic and gonadal tissues. Recent studies have revealed that Dsx1 expression is tightly controlled at transcriptional, post-transcriptional, and epigenetic levels to avoid sexual ambiguity. In this review, together with basic information on Dsx1 structure and expression, we introduce the multi-layered Dsx1 regulation and discuss how each regulation is interconnected for controlling male development in environmental sex-determining Daphnia.

Comparative phylotranscriptomics reveals putative sex differentiating genes across eight diverse bivalve species

Mollusks, especially bivalves, exhibit a great diversity of sex determining mechanisms, including both genetic and environmental sex determination. Some bivalve species can be gonochoristic (separate sexes), while others are hermaphroditic (sequential or simultaneous). Several models have been proposed for specific bivalve species, utilizing information gained from gene expression data, as well as limited RAD-seq data (e.g., from Crassostrea gigas). However, these mechanisms are not as well studied as those in model organisms (e.g., Mus musculus, Drosophila melanogaster, Caenorhabditis elegans) and many genes involved in sex differentiation are not well characterized. We used phylotranscriptomics to better understand which possible sex differentiating genes are in bivalves and how these genes relate to similar genes in diverse phyla. We collected RNAseq data from eight phylogenetically diverse bivalve species: Argopecten irradians, Ensis directus, Geukensia demissa, Macoma tenta, Mercenaria mercenaria, Mya arenaria, Mytilus edulis, and Solemya velum. Using these data, we assembled representative transcriptomes for each species. We then searched for candidate sex differentiating genes using BLAST and confirmed the identity of nine genes using phylogenetics analyses from nine phyla. To increase the confidence of identification, we included ten bivalve genomes in our analyses. From the analysis of doublesex and mab-3 related transcription factor (DMRT) genes, we confirmed the identify of a Mollusk-specific sex determining DMRT gene: DMRT1L. Based on gene expression data from M. edulis and previous research, DMRT1L and FoxL2 are key genes for male and female development, respectively.

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.

Identification of potential candidate genes involved in the sex determination cascade in an aquatic firefly, Sclerotia aquatilis (Coleoptera, Lampyridae)

Sexual differentiation, dimorphism, and courtship behavior are the downstream developmental programs of the sex determination cascade. The sex determination cascade in arthropods often involves key genes, transformer (tra), doublesex (dsx), transformer-2 (tra2), and fruitless (fru). These genes are conserved among insect taxa; however, they have never been reported in fireflies. In this study, the candidate genes for these key genes were identified for the first time in an aquatic firefly, Sclerotia aquatilis using transcriptome analysis. A comparative protein-protein interaction (PPI) network of sex determination cascade was reconstructed for S. aquatilis based on a network of a model insect, Drosophila melanogaster. Subsequently, a sex determination cascade in S. aquatilis was proposed based on the amino acid sequence structures and expression profiles of these candidates. This study describes the first efforts toward understanding the molecular control of sex determination cascade in fireflies.

Phylogenetic analysis and embryonic expression of panarthropod Dmrt genes

Background

One set of the developmentally important Doublesex and Male-abnormal-3 Related Transcription factors (Dmrt) is subject of intense research, because of their role in sex-determination and sexual differentiation. This likely non-monophyletic group of Dmrt genes is represented by the Drosophila melanogaster gene Doublesex (Dsx), the Caenorhabditis elegans Male-abnormal-3 (Mab-3) gene, and vertebrate Dmrt1 genes. However, other members of the Dmrt family are much less well studied, and in arthropods, including the model organism Drosophila melanogaster, data on these genes are virtually absent with respect to their embryonic expression and function.

Results

Here we investigate the complete set of Dmrt genes in members of all main groups of Arthropoda and a member of Onychophora, extending our data to Panarthropoda as a whole. We confirm the presence of at least four families of Dmrt genes (including Dsx-like genes) in Panarthropoda and study their expression profiles during embryogenesis. Our work shows that the expression patterns of Dmrt11E, Dmrt93B, and Dmrt99B orthologs are highly conserved among panarthropods. Embryonic expression of Dsx-like genes, however, is more derived, likely as a result of neo-functionalization after duplication.

Conclusions

Our data suggest deep homology of most of the panarthropod Dmrt genes with respect to their function that likely dates back to their last common ancestor. The function of Dsx and Dsx-like genes which are critical for sexual differentiation in animals, however, appears to be much less conserved.

Electronic supplementary material

The online version of this article (10.1186/s12983-019-0322-0) contains supplementary material, which is available to authorized users.

Origins and structural properties of novel and de novo protein domains during insect evolution

Over long time scales, protein evolution is characterized by modular rearrangements of protein domains. Such rearrangements are mainly caused by gene duplication, fusion and terminal losses. To better understand domain emergence mechanisms we investigated 32 insect genomes covering a speciation gradient ranging from ~ 2 to ~ 390 mya. We use established domain models and foldable domains delineated by hydrophobic cluster analysis (HCA), which does not require homologous sequences, to also identify domains which have likely arisen de novo, that is, from previously noncoding DNA. Our results indicate that most novel domains emerge terminally as they originate from ORF extensions while fewer arise in middle arrangements, resulting from exonization of intronic or intergenic regions. Many novel domains rapidly migrate between terminal or middle positions and single- and multidomain arrangements. Young domains, such as most HCA-defined domains, are under strong selection pressure as they show signals of purifying selection. De novo domains, linked to ancient domains or defined by HCA, have higher degrees of intrinsic disorder and disorder-to-order transition upon binding than ancient domains. However, the corresponding DNA sequences of the novel domains of de novo origins could only rarely be found in sister genomes. We conclude that novel domains are often recruited by other proteins and undergo important structural modifications shortly after their emergence, but evolve too fast to be characterized by cross-species comparisons alone.

The doublesex gene integrates multi-locus complementary sex determination signals in the Japanese ant, Vollenhovia emeryi

A female diploid, male haploid sex determination system (haplodiploidy) is found in hymenopteran taxa, such as ants, wasps, bees and sawflies. In this system, a single, complementary sex-determination (sl-CSD) locus functions as the primary sex-determination signal. In the taxa that has evolved this system, females and males are heterozygous and hemi/homozygous at the CSD locus, respectively. While the sl-CSD system enables females to alter sex ratios in the nest, it carries a high cost in terms of inbreeding, as individuals that are homozygous at the CSD locus become sterile diploid males. To counter this risk, some of hymenopteran species have evolved a multi-locus CSD (ml-CSD) system, which effectively reduces the proportion of sterile males. However, the mechanism by which these multiple primary signals are integrated and how they affect the terminal sex-differentiation signal of the molecular cascade have not yet been clarified. To resolve these questions, we examined the molecular cascade in the Japanese ant Vollenhovia emeryi, which we previously confirmed has two CSD loci. Here, we showed that the sex-determination gene, doublesex (dsx), which is highly conserved among phylogenetically distant taxa, is responsible for integrating two CSD signals in V. emeryi. After identifying and characterizing dsx, genotypes containing two CSD loci and splicing patterns of dsx were found to correspond to the sexual phenotype, suggesting that two primary signals are integrated into dsx. These findings will facilitate future molecular and functional studies of the sex determination cascade in V. emeryi, and shed light on the evolution and diversification of sex determination systems in insects.

The ubiquity and ancestry of insect doublesex

The doublesex (dsx) gene functions as a molecular switch at the base of the insect sex determination cascade, and triggers male or female somatic sexual differentiation in Drosophila. Having been reported from only seven current insect orders, the exact phylogenetic distribution of dsx within the largest Arthropod sub-phylum, the Hexapoda, is unknown. To understand the evolution of this integral gene relative to other arthropods, we tested for the presence of dsx within public EST and genome sequencing projects representative of all 32 hexapod orders. We find the dsx gene to be ubiquitous, with putative orthologs recovered from 30 orders. Additionally, we recovered both alternatively spliced and putative paralogous dsx transcripts from several orders of hexapods, including basal lineages, indicating the likely presence of these characteristics in the hexapod common ancestor. Of note, other arthropods such as chelicerates and crustaceans express two dsx genes, both of which are shown to lack alternative splicing. Furthermore, we discovered a large degree of length heterogeneity in the common region of dsx coding sequences within and among orders, possibly resulting from lineage-specific selective pressures inherent to each taxon. Our work serves as a valuable resource for understanding the evolution of sex determination in insects.

Sex- and tissue-specific functions of Drosophila doublesex transcription factor target genes

Primary sex-determination "switches" evolve rapidly, but Doublesex (DSX)-related transcription factors (DMRTs) act downstream of these switches to control sexual development in most animal species. Drosophila dsx encodes female- and male-specific isoforms (DSX(F) and DSX(M)), but little is known about how dsx controls sexual development, whether DSX(F) and DSX(M) bind different targets, or how DSX proteins direct different outcomes in diverse tissues. We undertook genome-wide analyses to identify DSX targets using in vivo occupancy, binding site prediction, and evolutionary conservation. We find that DSX(F) and DSX(M) bind thousands of the same targets in multiple tissues in both sexes, yet these targets have sex- and tissue-specific functions. Interestingly, DSX targets show considerable overlap with targets identified for mouse DMRT1. DSX targets include transcription factors and signaling pathway components providing for direct and indirect regulation of sex-biased expression.

Developmental link between sex and nutrition; doublesex regulates sex-specific mandible growth via juvenile hormone signaling in stag beetles

Sexual dimorphisms in trait expression are widespread among animals and are especially pronounced in ornaments and weapons of sexual selection, which can attain exaggerated sizes. Expression of exaggerated traits is usually male-specific and nutrition sensitive. Consequently, the developmental mechanisms generating sexually dimorphic growth and nutrition-dependent phenotypic plasticity are each likely to regulate the expression of extreme structures. Yet we know little about how either of these mechanisms work, much less how they might interact with each other. We investigated the developmental mechanisms of sex-specific mandible growth in the stag beetle Cyclommatus metallifer, focusing on doublesex gene function and its interaction with juvenile hormone (JH) signaling. doublesex genes encode transcription factors that orchestrate male and female specific trait development, and JH acts as a mediator between nutrition and mandible growth. We found that the Cmdsx gene regulates sex differentiation in the stag beetle. Knockdown of Cmdsx by RNA-interference in both males and females produced intersex phenotypes, indicating a role for Cmdsx in sex-specific trait growth. By combining knockdown of Cmdsx with JH treatment, we showed that female-specific splice variants of Cmdsx contribute to the insensitivity of female mandibles to JH: knockdown of Cmdsx reversed this pattern, so that mandibles in knockdown females were stimulated to grow by JH treatment. In contrast, mandibles in knockdown males retained some sensitivity to JH, though mandibles in these individuals did not attain the full sizes of wild type males. We suggest that moderate JH sensitivity of mandibular cells may be the default developmental state for both sexes, with sex-specific Dsx protein decreasing sensitivity in females, and increasing it in males. This study is the first to demonstrate a causal link between the sex determination and JH signaling pathways, which clearly interact to determine the developmental fates and final sizes of nutrition-dependent secondary-sexual characters.

Sexual dimorphisms such as the exaggerated antlers of deer, the enormous clawed chelae of crabs, and the horns and mandibles of beetles, are widespread across animal taxa and have fascinated biologists for centuries. Much recent work has uncovered the importance of the role of the sex-determination pathway in the expression of sexually dimorphic traits. However, critical interactions between this pathway and other growth regulatory mechanisms – for example, the physiological mechanisms involved in nutrition-dependent expression of these traits – are less well understood. In this study, we provide evidence of a developmental link between nutrition-sensitivity and sexual differentiation in the giant mandibles of the sexually dimorphic stag beetle, Cyclommatus metallifer. We examined the regulation and function of a key sex determination gene in animals, doublesex (dsx), and its interaction with juvenile hormone (JH), an important insect hormone known to regulate insect polyphenisms including the regulation of the disproportionate growth of male stag beetle mandibles. We found that Cmdsx changes mandibular responsiveness to JH in a sex-specific pattern. Based on these results, we hypothesize that sex-specific regulation of JH responsiveness is a developmental link between nutrition and sexual differentiation in stag beetles.

Nuclear genes with sex bias in Ruditapes philippinarum (Bivalvia, veneridae): Mitochondrial inheritance and sex determination in DUI species

Mitochondria are inherited maternally in most metazoans, but in bivalves with Doubly Uniparental Inheritance (DUI) a mitochondrial lineage is transmitted through eggs (F-type), and another through sperm (M-type). In DUI species, a sex-ratio distortion of the progeny was observed: some females produce a female-biased offspring (female-biased family), others a male-biased progeny (male-biased family), and others a 50:50 sex-ratio. A peculiar segregation pattern of M-type mitochondria in DUI organisms appears to be correlated with the sex bias of these families. According to a proposed model for the inheritance of M-type mitochondria in DUI, the transmission of sperm mitochondria is controlled by three nuclear genes, named W, X, and Z. An additional S gene with different dosage effect would be involved in sex determination. In this study, we analyzed structure and localization of three transcripts (psa, birc, and anubl1) with specific sex and family biases in the Manila clam Ruditapes philippinarum. In situ hybridization confirmed the localization of these transcripts in gametogenic cells. In other animals, homologs of these genes are involved in reproduction and ubiquitination. We hypothesized that these genes may have a role in sex determination and could also be responsible for the maintenance/degradation of spermatozoon mitochondria during embryo development of the DUI species R. philippinarum, so that we propose them as candidate factors of the W/X/Z/S system.

Molecular cloning of doublesex genes of four cladocera (water flea) species

Background

The gene doublesex (dsx) is known as a key factor regulating genetic sex determination in many organisms. We previously identified two dsx genes (DapmaDsx1 and DapmaDsx2) from a freshwater branchiopod crustacean, Daphnia magna, which are expressed in males but not in females. D. magna produces males by parthenogenesis in response to environmental cues (environmental sex determination) and we showed that DapmaDsx1 expression during embryonic stages is responsible for the male trait development. The D. magna dsx genes are thought to have arisen by a cladoceran-specific duplication; therefore, to investigate evolutionary conservation of sex specific expression of dsx genes and to further assess their functions in the environmental sex determination, we searched for dsx homologs in four closely related cladoceran species.

Results

We identified homologs of both dsx genes from, D. pulex, D. galeata, and Ceriodaphnia dubia, yet only a single dsx gene was found from Moina macrocopa. The deduced amino acid sequences of all 9 dsx homologs contained the DM and oligomerization domains, which are characteristic for all arthropod DSX family members. Molecular phylogenetic analysis suggested that the dsx gene duplication likely occurred prior to the divergence of these cladoceran species, because that of the giant tiger prawn Penaeus monodon is rooted ancestrally to both DSX1 and DSX2 of cladocerans. Therefore, this result also suggested that M. macrocopa lost dsx2 gene secondarily. Furthermore, all dsx genes identified in this study showed male-biased expression levels, yet only half of the putative 5’ upstream regulatory elements are preserved in D. magna and D. pulex.

Conclusions

The all dsx genes of five cladoceran species examined had similar amino acid structure containing highly conserved DM and oligomerization domains, and exhibited sexually dimorphic expression patterns, suggesting that these genes may have similar functions for environmental sex determination in cladocerans.

Positive and purifying selection influence the evolution of doublesex in the Anastrepha fraterculus species group

The gene doublesex (dsx) is considered to be under strong selective constraint along its evolutionary history because of its central role in somatic sex differentiation in insects. However, previous studies of dsx used global estimates of evolutionary rates to investigate its molecular evolution, which potentially miss signals of adaptive changes in generally conserved genes. In this work, we investigated the molecular evolution of dsx in the Anastrepha fraterculus species group (Diptera, Tephritidae), and test the hypothesis that this gene evolved solely by purifying selection using divergence-based and population-based methods. In the first approach, we compared sequences from Anastrepha and other Tephritidae with other Muscomorpha species, analyzed variation in nonsynonymous to synonymous rate ratios (dN/dS) in the Tephritidae, and investigated radical and conservative changes in amino acid physicochemical properties. We show a general selective constraint on dsx, but with signs of positive selection mainly in the common region. Such changes were localized in alpha-helices previously reported to be involved in dimer formation in the OD2 domain and near the C-terminal of the OD1 domain. In the population-based approach, we amplified a region of 540 bp that spanned almost all of the region common to both sexes from 32 different sites in Brazil. We investigated patterns of selection using neutrality tests based on the frequency spectrum and locations of synonymous and nonsynonymous mutations in a haplotype network. As in the divergence-based approach, these analyses showed that dsx has evolved under an overall selective constraint, but with some events of positive selection. In contrast to previous studies, our analyses indicate that even though dsx has indeed evolved as a conserved gene, the common region of dsx has also experienced bouts of positive selection, perhaps driven by sexual selection, during its evolution.

Structural analysis of the UBA domain of X-linked inhibitor of apoptosis protein reveals different surfaces for ubiquitin-binding and self-association

Background

Inhibitor of apoptosis proteins (IAPs) belong to a pivotal antiapoptotic protein family that plays a crucial role in tumorigenesis, cancer progression, chemoresistance and poor patient-survival. X-linked inhibitor of apoptosis protein (XIAP) is a prominent member of IAPs attracting intense research because it has been demonstrated to be a physiological inhibitor of caspases and apoptosis. Recently, an evolutionarily conserved ubiquitin-associated (UBA) domain was identified in XIAP and a number of RING domain-bearing IAPs. This has placed the IAPs in the group of ubiquitin binding proteins. Here, we explore the three-dimensional structure of the XIAP UBA domain (XIAP-UBA) and how it interacts with mono-ubiquitin and diubiquitin conjugates.

Principal Findings

The solution structure of the XIAP-UBA domain was determined by NMR spectroscopy. XIAP-UBA adopts a typical UBA domain fold of three tightly packed α-helices but with an additional N-terminal 3₁₀ helix. The XIAP-UBA binds mono-ubiquitin as well as Lys48-linked and linear-linked diubiquitins at low-micromolar affinities. NMR analysis of the XIAP-UBA–ubiquitin interaction reveals that it involves the classical hydrophobic patches surrounding Ile44 of ubiquitin and the conserved MGF/LV motif surfaces on XIAP-UBA. Furthermore, dimerization of XIAP-UBA was observed. Mapping of the self-association surface of XIAP-UBA reveals that the dimerization interface is formed by residues in the N-terminal 3₁₀ helix, helix α1 and helix α2, separate from the ubiquitin-binding surface.

Conclusion

Our results provide the first structural information of XIAP-UBA and map its interaction with mono-ubiquitin, Lys48-linked and linear-linked diubiquitins. The notion that XIAP-UBA uses different surfaces for ubiquitin-binding and self-association provides a plausible model to explain the reported selectivity of XIAP in binding polyubiquitin chains with different linkages.

Improving protein structure prediction using multiple sequence-based contact predictions

Although residue-residue contact maps dictate the topology of proteins, sequence-based ab initio contact predictions have been found little use in actual structure prediction due to the low accuracy. We developed a composite set of nine SVM-based contact predictors that are used in I-TASSER simulation in combination with sparse template contact restraints. When testing the strategy on 273 nonhomologous targets, remarkable improvements of I-TASSER models were observed for both easy and hard targets, with p value by Student's t test<0.00001 and 0.001, respectively. In several cases, template modeling score increases by >30%, which essentially converts "nonfoldable" targets into "foldable" ones. In CASP9, I-TASSER employed ab initio contact predictions, and generated models for 26 FM targets with a GDT-score 16% and 44% higher than the second and third best servers from other groups, respectively. These findings demonstrate a new avenue to improve the accuracy of protein structure prediction especially for free-modeling targets.

De Novo assembly of the Manila clam Ruditapes philippinarum transcriptome provides new insights into expression bias, mitochondrial doubly uniparental inheritance and sex determination

Males and females share the same genome, thus, phenotypic divergence requires differential gene expression and sex-specific regulation. Accordingly, the analysis of expression patterns is pivotal to the understanding of sex determination mechanisms. Many bivalves are stable gonochoric species, but the mechanism of gonad sexualization and the genes involved are still unknown. Moreover, during the period of sexual rest, a gonad is not present and sex cannot be determined. A mechanism associated with germ line differentiation in some bivalves, including the Manila clam Ruditapes philippinarum, is the doubly uniparental inheritance (DUI) of mitochondria, a variation of strict maternal inheritance. Two mitochondrial lineages are present, one transmitted through eggs and the other through sperm, as well as a mother-dependent sex bias of the progeny. We produced a de novo annotation of 17,186 transcripts from R. philippinarum and compared the transcriptomes of males and females and identified 1,575 genes with strong sex-specific expression and 166 sex-specific single nucleotide polymorphisms, obtaining preliminary information about genes that could be involved in sex determination. Then we compared the transcriptomes between a family producing predominantly females and a family producing predominantly males to identify candidate genes involved in regulation of sex-specific aspects of DUI system, finding a relationship between sex bias and differential expression of several ubiquitination genes. In mammalian embryos, sperm mitochondria are degraded by ubiquitination. A modification of this mechanism is hypothesized to be responsible for the retention of sperm mitochondria in male embryos of DUI species. Ubiquitination can additionally regulate gene expression, playing a role in sex determination of several animals. These data enable us to develop a model that incorporates both the DUI literature and our new findings.

Environmental sex determination in the branchiopod crustacean Daphnia magna: deep conservation of a Doublesex gene in the sex-determining pathway

Sex-determining mechanisms are diverse among animal lineages and can be broadly divided into two major categories: genetic and environmental. In contrast to genetic sex determination (GSD), little is known about the molecular mechanisms underlying environmental sex determination (ESD). The Doublesex (Dsx) genes play an important role in controlling sexual dimorphism in genetic sex-determining organisms such as nematodes, insects, and vertebrates. Here we report the identification of two Dsx genes from Daphnia magna, a freshwater branchiopod crustacean that parthenogenetically produces males in response to environmental cues. One of these genes, designated DapmaDsx1, is responsible for the male trait development when expressed during environmental sex determination. The domain organization of DapmaDsx1 was similar to that of Dsx from insects, which are thought to be the sister group of branchiopod crustaceans. Intriguingly, the molecular basis for sexually dimorphic expression of DapmaDsx1 is different from that of insects. Rather than being regulated sex-specifically at the level of pre–mRNA splicing in the coding region, DapmaDsx1 exhibits sexually dimorphic differences in the abundance of its transcripts. During embryogenesis, expression of DapmaDsx1 was increased only in males and its transcripts were primarily detected in male-specific structures. Knock-down of DapmaDsx1 in male embryos resulted in the production of female traits including ovarian maturation, whereas ectopic expression of DapmaDsx1 in female embryos resulted in the development of male-like phenotypes. Expression patterns of another D. magna Dsx gene, DapmaDsx2, were similar to those of DapmaDsx1, but silencing and overexpression of this gene did not induce any clear phenotypic changes. These results establish DapmaDsx1 as a key regulator of the male phenotype. Our findings reveal how ESD is implemented by selective expression of a fundamental genetic component that is functionally conserved in animals using GSD. We infer that there is an ancient, previously unidentified link between genetic and environmental sex determination.

Sex determination is a fundamental biological process that can be broadly divided into two major categories. In genetic sex determination (GSD), sex-specific differentiation results from intrinsic genetic differences between males and females, whereas environmental sex determination (ESD) relies on environmental signals to induce male or female sex determination. In contrast to model organisms that utilize GSD system, environmental sex-determining organisms are poor genetic models. Therefore, although candidate genes involved in ESD have been found in vertebrates, their functions have remained largely unknown, impairing our understanding of ESD and making the comparison of sex-determining genes between both systems difficult. Here, we report the identification of a gene responsible for the production of males during environmental sex determination in the crustacean Daphnia. This gene is homologous to the Doublesex gene that is functionally conserved in animals that use GSD. Expression of Doublesex was increased primarily in male-specific structures. Gain- and loss-of-function analyses established that Daphnia Doublesex gene is a major effector that regulates the male phenotype in Daphnia. We infer that there is an ancient, previously unidentified link between genetic and environmental sex determination.

Comparative and functional genomics of lipases in holometabolous insects

Lipases have key roles in insect lipid acquisition, storage and mobilisation and are also fundamental to many physiological processes underpinning insect reproduction, development, defence from pathogens and oxidative stress, and pheromone signalling. We have screened the recently sequenced genomes of five species from four orders of holometabolous insects, the dipterans Drosophila melanogaster and Anopheles gambiae, the hymenopteran Apis mellifera, the moth Bombyx mori and the beetle Tribolium castaneum, for the six major lipase families that are also found in other organisms. The two most numerous families in the insects, the neutral and acid lipases, are also the main families in mammals, albeit not in Caenorhabditis elegans, plants or microbes. Total numbers of the lipases vary two-fold across the five insect species, from numbers similar to those in mammals up to numbers comparable to those seen in C. elegans. Whilst there is a high degree of orthology with mammalian lipases in the other four families, the great majority of the insect neutral and acid lipases have arisen since the insect orders themselves diverged. Intriguingly, about 10% of the insect neutral and acid lipases have lost motifs critical for catalytic function. Examination of the length of lid and loop regions of the neutral lipase sequences suggest that most of the insect lipases lack triacylglycerol (TAG) hydrolysis activity, although the acid lipases all have intact cap domains required for TAG hydrolysis. We have also reviewed the sequence databases and scientific literature for insights into the expression profiles and functions of the insect neutral and acid lipases and the orthologues of the mammalian adipose triglyceride lipase which has a pivotal role in lipid mobilisation. These data suggest that some of the acid and neutral lipase diversity may be due to a requirement for rapid accumulation of dietary lipids. The different roles required of lipases at the four discrete life stages of holometabolous insects may also contribute to the diversity of lipases required by insects. In addition, insects use lipases to perform roles for which there are no correlates in mammals, including as yolk and male accessory gland proteins.

Avid interactions underlie the Lys63-linked polyubiquitin binding specificities observed for UBA domains

Ubiquitin (denoted Ub) receptor proteins as a group must contain a diverse set of binding specificities to distinguish the many forms of polyubiquitin (polyUb) signals. Previous studies suggested that the large class of ubiquitin-associated (UBA) domains contains members with intrinsic specificity for Lys63-linked polyUb or Lys48-linked polyUb, thus explaining how UBA-containing proteins can mediate diverse signaling events. Here we show that previously observed Lys63-polyUb selectivity in UBA domains is the result of an artifact in which the dimeric fusion partner, glutathione S-transferase (GST), positions two UBAs for higher affinity, avid interactions with Lys63-polyUb, but not with Lys48-polyUb. Freed from GST, these UBAs are either nonselective or prefer Lys48-polyUb. Accordingly, NMR experiments reveal no Lys63-polyUb-specific binding epitopes for these UBAs. We reexamine previous conclusions based on GST-UBAs and present an alternative model for how UBAs achieve a diverse range of linkage specificities.

Identification and characterization of the doublesex gene of Nasonia

The doublesex (dsx) gene of the parasitic wasp Nasonia vitripennis is described and characterized. Differential splicing of dsx transcripts has been shown to induce somatic sexual differentiation in Diptera and Lepidoptera, but not yet in other insect orders. Two spliceforms of Nasonia dsx mRNA are differentially expressed in males and females. In addition, in a gynandromorphic line that produces haploids (normally males) with full female phenotypes, these individuals show the female spliceform, providing the first demonstration of a direct association of dsx with somatic sex differentiation in Hymenoptera. Finally, the DNA binding (DM) domain of Nasonia dsx clusters phylogenetically with dsx from other insects, and Nasonia dsx shows microsynteny with dsx of Apis, further supporting identification of the dsx orthologue in Nasonia.

Doublesex and the regulation of sexual dimorphism in Drosophila melanogaster: structure, function, and mutagenesis of a female-specific domain

The DSX (Doublesex) transcription factor regulates somatic sexual differentiation in Drosophila. Female and male isoforms (DSX F and DSX M) are formed due to sex-specific RNA splicing. DNA recognition, mediated by a shared N-terminal zinc module (the DM domain), is enhanced by a C-terminal dimerization element. Sex-specific extension of this element in DSX F and DSX M leads to assembly of distinct transcriptional preinitiation complexes. Here, we describe the structure of the extended C-terminal dimerization domain of DSX F as determined by multidimensional NMR spectroscopy. The core dimerization element is well ordered, giving rise to a dense network of interresidue nuclear Overhauser enhancements. The structure contains dimer-related UBA folds similar to those defined by x-ray crystallographic studies of a truncated domain. Whereas the proximal portion of the female tail extends helix 3 of the UBA fold, the distal tail is disordered. Ala substitutions in the proximal tail disrupt the sex-specific binding of IX (Intersex), an obligatory partner protein and putative transcriptional coactivator; IX-DSX F interaction is, by contrast, not disrupted by truncation of the distal tail. Mutagenesis of the UBA-like dimer of DSX F highlights the importance of steric and electrostatic complementarity across the interface. Two temperature-sensitive mutations at this interface have been characterized in yeast model systems. One weakens a network of solvated salt bridges, whereas the other perturbs the underlying nonpolar interface. These mutations confer graded gene-regulatory activity in yeast within a physiological temperature range and so may provide novel probes for genetic analysis of a sex-specific transcriptional program in Drosophila development.

Structural Basis for UBA-mediated Dimerization of c-Cbl Ubiquitin Ligase

Ligand-induced down-regulation by the ubiquitin-protein ligases, c-Cbl and Cbl-b, controls signaling downstream from many receptor-tyrosine kinases (RTK). Cbl proteins bind to phosphotyrosine residues on activated RTKs to affect ligand-dependent ubiquitylation of these receptors targeting them for degradation in the lysosome. Both c-Cbl and Cbl-b contain a ubiquitin-associated (UBA) domain, which is important for Cbl dimerization and tyrosine phosphorylation; however, the mechanism of UBA-mediated dimerization and its requirement for Cbl biological activity is unclear. Here, we report the crystal structure of the UBA domain of c-Cbl refined to 2.1-A resolution. The structure reveals the protein is a symmetric dimer tightly packed along a large hydrophobic surface formed by helices 2 and 3. NMR chemical shift mapping reveals heterodimerization can occur with the related Cbl-b UBA domain via the same surface employed for homodimerization. Disruption of c-Cbl dimerization by site-directed mutagenesis impairs c-Cbl phosphorylation following activation of the Met/hepatocyte growth factor RTK and c-Cbl-dependent ubiquitination of Met. This provides direct evidence for a role of Cbl dimerization in terminating signaling following activation of RTKs.

Conformational instability of the MARK3 UBA domain compromises ubiquitin recognition and promotes interaction with the adjacent kinase domain

The Par-1/MARK protein kinases play a pivotal role in establishing cellular polarity. This family of kinases contains a unique domain architecture, in which a ubiquitin-associated (UBA) domain is located C-terminal to the kinase domain. We have used a combination of x-ray crystallography and NMR dynamics experiments to understand the interaction of the human (h) MARK3 UBA domain with the adjacent kinase domain as compared with ubiquitin. The x-ray crystal structure of the linked hMARK3 kinase and UBA domains establishes that the UBA domain forms a stable intramolecular interaction with the N-terminal lobe of the kinase domain. However, solution-state NMR studies of the isolated UBA domain indicate that it is highly dynamic, undergoing conformational transitions that can be explained by a folding-unfolding equilibrium. NMR titration experiments indicated that the hMARK3 UBA domain has a detectable but extremely weak affinity for mono ubiquitin, which suggests that conformational instability of the isolated hMARK3 UBA domain attenuates binding to ubiquitin despite the presence of residues typically involved in ubiquitin recognition. Our data identify a molecular mechanism through which the hMARK3 UBA domain has evolved to bind the kinase domain, in a fashion that stabilizes an open conformation of the N- and C-terminal lobes, at the expense of its capacity to engage ubiquitin. These results may be relevant more generally to the 30% of UBA domains that lack significant ubiquitin-binding activity, and they suggest a unique mechanism by which interaction domains may evolve new binding properties.

Structural basis for ubiquitin-mediated dimerization and activation of the ubiquitin protein ligase Cbl-b

Cbl proteins are E3 ubiquitin ligases that are negative regulators of many receptor tyrosine kinases. Cbl-b and c-Cbl contain a ubiquitin-associated (UBA) domain, which is present in a variety of proteins involved in ubiquitin-mediated processes. Despite high sequence identity, Cbl UBA domains display remarkably different ubiquitin-binding properties. Here, we report the crystal structure of the UBA domain of Cbl-b in complex with ubiquitin at 1.9 A resolution. The structure reveals an atypical mechanism of ubiquitin recognition by the first helix of the UBA. Helices 2 and 3 of the UBA domain form a second binding surface, which mediates UBA dimerization in the crystal and in solution. Site-directed mutagenesis demonstrates that Cbl-b dimerization is regulated by ubiquitin binding and required for tyrosine phosphorylation of Cbl-b and ubiquitination of Cbl-b substrates. These studies demonstrate a role for ubiquitin in regulating biological activity by promoting protein dimerization.

Multiple alternative splicing of mouse Dmrt1 during gonadal differentiation

Evolutionarily conserved Dmrt1 encodes a transcriptional regulator that is expressed exclusively in the gonads and is required for testis differentiation. Here we report that four transcripts of the mouse Dmrt1 were generated in developing gonads and adult testis by alternative splicing. Dmrt1 a encodes the known protein with 374 amino acids. Dmrt1 b, Dmrt1 c, and Dmrt1 d encode predicted proteins with 212, 257, and 194 amino acids, respectively. Dmrt1 a2 and Dmrt1 a3 have anterior alternative polyadenylation signals in 3'UTR than the known Dmrt1 transcript (Dmrt1 a1). Dmrt1 a2 lacks 18 nucleic acids in 3'UTR. Interestingly, Dmrt1 b lacks exon5, Dmrt1 c lacks exon3, and Dmrt1 d lacks both exon1 and exon2 which encode DM domain. RT-PCR showed that all of Dmrt1 transcripts were only detectable in adult testis. However, during gonadogenesis, the multiple alternatively spliced transcripts all had gonad-specific and sexually dimorphic expression profiles. Northern blot and real time fluorescent quantitative RT-PCR further indicated that the expression of Dmrt1 a was dominantly higher than those of Dmrt1 b, c, and d, although they showed a similar up and down pattern of expression during embryo development. The expression of Dmrt1 transcripts climbed up to a climax at 13.5dpc in both male and female gonad. Afterwards, expression in male gonad decreased to a low level and maintained, whereas female one reduced rapidly until undetectable in adult ovary. These results provide new insight into roles of regulation at level of splicing of the Dmrt1 in governing sex differentiation.

Male-biased expression of X-chromosomal DM domain-less Dmrt8 genes in the mouse

The vertebrate DMRT gene family encodes putative transcription factors related to the sexual regulators Doublesex (Drosophila melanogaster) and MAB-3 (Caenorhabditis elegans). They share a highly conserved DNA binding motif, the DM domain. In human and mouse seven DMRT genes (DMRT1-DMRT7) have been analyzed. DMRT8, a gene related to DMRT7, is located on the X chromosome in placental mammals. While DMRT8 is single copy in most mammals, three copies are present in mouse, rat, and rabbit. Despite the loss of the DM domain, DMRT8 genes have been maintained in the mammalian lineage, suggesting a DM domain-independent function. In adult mouse, two Dmrt8 genes are expressed exclusively in testis. Dmrt8.1 mRNA was detected in Sertoli cells by in situ hybridization. In embryos, Dmrt8.2 shows a dynamic expression restricted to male and female gonads and might therefore be involved in sexual development in the mouse.

Regulation of sexual dimorphism: mutational and chemogenetic analysis of the doublesex DM domain

Doublesex (dsx) is a transcription factor in Drosophila that regulates somatic sexual differentiation. Male- and female-specific splicing isoforms of DSX share a novel DNA-binding domain, designated the DM motif. Broadly conserved among metazoan sex-determining factors, the DM domain contains a nonclassical zinc module and binds in the DNA minor groove. Here, we characterize the DM motif by site-directed and random mutagenesis using a yeast one-hybrid (Y1H) system and extend this analysis by chemogenetic complementation in vitro. The Y1H system is based on a sex-specific Drosophila enhancer element and validated through studies of intersexual dsx mutations. We demonstrate that the eight motif-specific histidines and cysteines engaged in zinc coordination are each critical and cannot be interchanged; folding also requires conserved aliphatic side chains in the hydrophobic core. Mutations that impair DNA binding tend to occur at conserved positions, whereas neutral substitutions occur at nonconserved sites. Evidence for a specific salt bridge between a conserved lysine and the DNA backbone is obtained through the synthesis of nonstandard protein and DNA analogs. Together, these results provide molecular links between the structure of the DM domain and its function in the regulation of sexual dimorphism.