Skip to main content

Disorders of sex development expose transcriptional autonomy of genetic sex and androgen-programmed hormonal sex in human blood leukocytes

Abstract

Background

Gender appears to be determined by independent programs controlled by the sex-chromosomes and by androgen-dependent programming during embryonic development. To enable experimental dissection of these components in the human, we performed genome-wide profiling of the transcriptomes of peripheral blood mononuclear cells (PBMC) in patients with rare defined "disorders of sex development" (DSD, e.g., 46, XY-females due to defective androgen biosynthesis) compared to normal 46, XY-males and 46, XX-females.

Results

A discrete set of transcripts was directly correlated with XY or XX genotypes in all individuals independent of male or female phenotype of the external genitalia. However, a significantly larger gene set in the PBMC only reflected the degree of external genital masculinization independent of the sex chromosomes and independent of concurrent post-natal sex steroid hormone levels. Consequently, the architecture of the transcriptional PBMC-"sexes" was either male, female or even "intersex" with a discordant alignment of the DSD individuals' genetic and hormonal sex signatures.

Conclusion

A significant fraction of gene expression differences between males and females in the human appears to have its roots in early embryogenesis and is not only caused by sex chromosomes but also by long-term sex-specific hormonal programming due to presence or absence of androgen during the time of external genital masculinization. Genetic sex and the androgen milieu during embryonic development might therefore independently modulate functional traits, phenotype and diseases associated with male or female gender as well as with DSD conditions.

Background

Biological sex is the product of genetic and hormonal programs activated during embryonic development, with the most striking gender differences occurring in the internal and external genitalia. Biological sex is a core determinant of sex-specific behavior [1–3] and also a major modifier of a variety of maladies including cancer [4, 5] and autoimmune diseases [6]. Delineation of the determinants of biological sex could provide valuable insights into gender-specific behaviors, diseases, treatment responses and outcomes [7]. Sex-specific differences in gene expression have been documented and they are hypothesized to reflect gender-specific development [8, 9]. However, sex-specific gene expression appears to be the product of two distinct pathways: sex chromosome genotype (46,XX versus 46,XY) [8, 9] and androgen-dependent programming during embryogenesis [10, 11]. In contrast to established concepts of androgen programming in rodents [12–14] it has not been possible to tease out the contributions of these two pathways experimentally in normal male and normal female humans since they always work in parallel.

Disorders of sex development (DSD, previously named "intersex") are the product of chromosomal and hormonal abnormalities during development that result in genital ambiguity and incongruent combinations of gender-specific external genitalia, reproductive ducts and gonads in a single individual. Therefore, DSD can serve as an "experiment of nature" to decipher the genetic and hormonal components of sex-specific gene expression in the human. Genetic and hormonal influences on gender extend beyond the reproductive tract, modifying body proportions, hair distribution, gender identity, sex-specific behavior and many other features that also affect DSD patients clinically in multiple ways [1–3, 15–17]. Therefore, improved understanding of the architecture of sex specific gene expression may also contribute to better clinical assessment of DSD-patients. In the present study, we hypothesized that sex chromosome controlled gene expression and androgen action during embryonic development in the first trimester as documented by the clinical degree of genital masculinization after birth could be distinguished by comparing normal males and females to individuals with DSD having well-defined clinical, hormonal and molecular characteristics. Because of simple accessibility, lack of overt organ-specific topography biasing gene expression [10, 18] and documented androgen receptor expression [19] we used peripheral blood mononuclear cells (PBMC) for genome-wide transcriptome analysis. Our study is the first that profiles blood gene expression data in DSD individuals.

Results

Genome-wide transcription profiling

Transcript profiles for 9 normal male and 10 normal female controls were measured using spotted cDNA microarrays with more than 44,000 cDNA elements representing 26,000 unique genes. The Significance Analysis of Microarrays (SAM) procedure [20] was used to identify 157 sex-specific transcripts (121 unique named genes) (false discovery rate < 0.09): 136 transcripts with higher expression in males (86.6%) and 21 (13.4%) with higher expression in females. Gene expression across these 157 sex-dimorphic transcripts was analyzed across all 19 controls and 14 DSD-patients. The DSD patients were comprised of ten 46,XY subjects whose genitalia ranged from normal female to Prader 4, three 46,XX subjects with masculinized genitalia from high prenatal androgenic steroid levels due to congenital adrenal hyperplasia (CAH, 21-hydroxlase deficiency) and one Prader 4 45,X0/46,XY subject. In addition, one normally masculinized 46,XY male with CAH was investigated (Additional file 1 and Figure 1A, B, C, additional files 2, 3, 4, 5).

Figure 1
figure 1

Hierarchical clustering analysis of 34 microarray experiments on PBMC total RNA derived from 9 normal males, 10 normal females, one 46,XY-male with CAH and 14 individuals with DSD [see additional files 2, 3, 4, 5]. (A) Transcript levels of 157 transcripts identified by SAM as differing significantly between PBMC from 9 normal males and 10 normal females. Transcripts are grouped by hierarchical cluster analysis and displayed in rows while experiments are displayed in columns. Expression values per gene are centered by the mean log2 Cy5/Cy3 normalized ratio across the arrays. Increasing blue intensity in the heat map corresponds to higher relative transcript levels compared to the mean expression level across all 34 array experiments. Increasing yellow intensity corresponds to relatively decreased transcript levels compared to the mean. Dark grey corresponds to missing data. Examples taken from these gene clusters are marked by color within the gene tree on the left of the heat map. The red gene tree corresponds to the enlarged gene cluster D, the blue gene tree corresponds to F, the pink gene tree corresponds to G, and the green gene tree corresponds to E, respectively. Gene symbols of the named transcripts are shown on the right. (B) Enlarged cluster dendrogram of the PBMC samples demonstrating the degree of relatedness (Pearson correlation) between the expression patterns of the 157 transcripts. The length of the arms of the dendrogram reflects the degree of correlation between experiments. Samples are color coded to reflect the degree of external genital virilization according to Prader that had been applied to both XX and XY individuals to enable comparability of genital phenotypes independent of the sex chromosomes. The dark grey bar below the experiment cluster represents individuals with an XX karyotype, light grey represents an XY karyotype. The major subdivison of the individuals corresponds strictly to the karyotype. The second level of subdivision in each of the major arms reflects mostly the phenotype of external genital virilization independent of the karyotype. (C) Schematic representation of the Prader stages of external genital virilization. (D) X-chromosome gene cluster (E) Y-chromosome gene cluster (F) Genes with predominantly higher transcript levels in the phenotypic females and the normal females independent of the karyotype (G) Genes with higher transcript levels in the strongly virilized individuals and the normal males independent of the karyotype.

Hierarchical cluster analysis sorted the subjects into 2 groups that directly correlated with the sex chromosomes (Figure 1A, B, C). Segregation of the 46,XX and 46,XY subjects into 2 major groups was driven by large differences in expression of 11 of the 157 transcripts corresponding to 8 unique genes (Figure 1B, D, E). High expression of XIST (X (inactive)-specific transcript) and LOC554203 (Xq13.2) characterized genetic females (Figure 1A, D), while DDX3Y (DEAD (Asp-Glu-Ala-Asp) box polypeptide 3, Y-linked) and five additional Y-linked genes characterized genetic males (Figure 1E). RT-PCR confirmed abundant expression of XIST mRNA in two female controls (N-020, N-031) and in two 46,XX virilized CAH-individuals (P-103, P-105) and showed an absence of expression in two 46,XY-controls (N-003, N-015) and two sex-reversed 46,XY DSD-females (P-019, P033). RT-PCR also confirmed that DDX3Y transcript levels were abundant in these same 46,XY individuals, but absent in the 46,XX subjects (data not shown).

Within each of the 46,XX and 46,XY groups were 2 subgroups that largely correlated with the degree of genital masculinization (Figure 1B, F, G). The 46,XX normal female subjects clustered separately from the three 46,XX subjects with genital masculinization and CAH based on differences in expression across 146 of the 157 sex-specific transcripts. Similarly, the 46,XY individuals separated into a group comprised of normal males and highly masculinized individuals with DSD (with one exception) and a group with predominantly female external genitalia, again based on expression differences across the same 146 transcripts. Only two patients did not follow this pattern. Subject P-115, who clustered with phenotypic males, had 46,XY DSD and Prader 1–2 genitalia due to an unclassified defect in androgen biosynthesis. Subject P-088 had mixed gonadal dysgenesis due to 45,X0/46,XY mosaicism and Prader 4 genitalia yet clustered with phenotypic females. Since mosaicism can vary between tissue compartments and since developmental genes may modify the site-specific expression level of the androgen receptor in the external genitalia [21, 22], some discrepancies between genital phenotype and systemic androgenization are not unexpected. Using qRT-PCR, we evaluated transcript levels of FZD6 (Frizzled 6), since levels measured by microarray analysis correlated with genital phenotype (Figure 1G). Similar to the microarray experiments, mean FZD6 transcript levels were significantly higher in normal 46,XY-males and in strongly virilized 46,XX-CAH females compared with normal 46,XX-females and 46,XY DSD-females (Figure 2).

Figure 2
figure 2

Ratio of FZD6 transcript levels in 2 46,XY normal males (N-003; N-015) and 46,XX strongly virilized females (P-103; P-105) compared with two 46,XX normal females (N-020; N-031) and two 46,XY females (P-019; P-033) by semi-quantitative RT-PCR. The y-axis reflects the ratios of expression levels of the normal male – and the strongly virilized DSD individuals, respectively, divided through the phenotypic female individuals (normal females and 46,XY-sex reversed females) as indicated. RT-PCR confirms higher expression of FZD6 in virilized individuals.

Influence of X- and Y-chromosome genes

We excluded the 11 X- and Y-linked transcripts that correlated with the two sex chromosome clusters and reanalyzed the samples by hierarchical clustering. With the remaining 146 sex-specific transcripts, subjects no longer grouped according to karyotype, but separated into 2 groups largely distinguished by external genitalia masculinization (Figure 3, additional files 6, 7, 8, 9). Four of the five highly virilized Prader 4 DSD-individuals and all normal males clustered together as did the normal females and subjects with feminized external genitalia. The mean Prader stage of the DSD-individuals excluding the normal controls was 1.6 in the "female" group and 3.2 in those that clustered with the normal males (p = 0.02) (Figure 3). When we excluded all normal controls from the cluster analysis, DSD-individuals separated into the same 2 groups (data not shown) indicating that the normal controls used for SAM had not significantly influenced the clustering of the DSD-patients.

Figure 3
figure 3

(A) Hierarchical clustering analysis using the 146 transcripts that did not correlate with karyotype. Genes of the two sex chromosome gene clusters (see D and E in Figure 1) were excluded from calculating the dendrograms, but are displayed in the figure to show their distribution across the samples. All color codes are the same as in figure 1. (B) Enlarged dendrogram demonstrating that major subdivision of the individuals mirrors the degree of genital virilization independent of karyotype. (C). Mean Prader stage of genital virilization (± 1 SD) of the DSD individuals the main branches of the dendrogram excluding the normal controls. The mean Prader stage on the left side is Prader "1.6" and "3.2" on the right (p = 0.02) [see also additional files 6, 7, 8, 9].

Influence of concurrent sex hormone levels

At the time of collection of PBMC RNA the 46,XY DSD-subjects with feminized genitalia varied widely in their sex steroid hormone status. Individuals who clustered in the female (left main branch) of Figure 3 included pre- and post-pubertal normal individuals with intact steroid biosynthetic pathways and individuals lacking steroid hormones (subjects P-017, P-032 and P-089 with 46,XY DSD who had been gonadectomized and were not on hormone replacement therapy and subject P-033 who had an inherited defect in, and almost complete lack of, steroid biosynthesis). This cluster also included gonadectomized 46,XY DSD-individuals on estrogen replacement, such as subject P-027 who showed normal adult female plasma estradiol of 492 pmol/L (normal range 73–1065 pmol/L).

Similarly, the phenotypic "male" cluster (right main branch of Figure 3) contained normal postpubertal males (several with documented fertility and presumably adult male testosterone concentrations) and prepubertal individuals (N-028, P-103, P-105, P-110, P-113). The prepubertal subjects P-103, P-105 and P-110 had low serum testosterone levels compared with age – and male sex adjusted normal ranges at the time of blood sampling (P-103: 0.1 nmol/L; P-105: 0 nmol/l; P-110: 0 nmol/l; normal range 0.14–1.32 nmol/l). Subject P-112, who was 46,XX and had a 21-hydroxylase deficiency usually resulting in high levels of androgenic steroids, had been treated prenatally with dexamethasone. Despite having only slightly virilized Prader 1 genitalia, transcript levels of the PBMCs more closely resembled phenotypic males, suggesting residual effects due to elevated adrenal androgens. One subject (P-004), born with Prader 4 external genitalia due to 5α-reductase type II deficiency and gonadectomized in early childhood, was on estrogens and gestagens at the time of blood sampling (plasma estradiol: 822 pmol/L). Therefore, the assignment to either the "female" or "male" clusters based on the 146 sex-specific transcript levels reflected neither pubertal status nor concurrent serum androgen or estrogen levels.

Biological function of sex-specific transcripts

PANTHER (Protein ANalysis THrough Evolutionary Relationships) [23] identified 121 named genes from the 157 sex-specific transcripts. Compared with the NCBI Homo sapiens reference list, 28 biological processes were significantly enriched (p < 0.05) and included developmental processes, extracellular matrix, cell communication and several that are involved in proliferation (table 1). GSEA (Gene Set Enrichment Analysis) [24] was performed using all interpretable genes from normal males and females using 1,107 gene sets from the Molecular Signature Database http://www.broad.mit.edu/gsea/msigdb. 18 gene sets were significantly enriched in males, one in females, of which the majority had been associated with proliferation, often in a context of cancer (e.g., Figure 4[25], http://microarray-pubs.stanford.edu/sob). Interestingly, many of the sex-specific genes identified in our study have published roles in human cancer, including TFPI2 (Tissue factor pathway inhibitor 2) (Figure 1G) [26], LUM (Lumican) [27], PDPN (Podoplanin) [28], SDC1 (Syndecan-1) [29], ESCO2 (Establishment of cohesion 1 homolog 2) (Figure 1G) [30], SSX2IP (Synovial sarcoma, X breakpoint 2 interacting protein) [31], CTGF (Connective tissue growth factor) [32], MMP2 (Matrix metallopeptidase 2) [33]. IGF1R (Figure 1F) [34], SERPINA1 (Serpin peptidase inhibitor, clade A) [35] and SEPT6 (septin 6) (Figure 1F) [36]. A subgroup of these genes function in DNA damage repair, such as CHEK1 (CHK1 checkpoint homolog) [37], XRCC1 (X-ray repair complementing defective repair in Chinese hamster cells 1) [38], H2AFX (H2A histone family, member X) [39], ALKBH (alkylation repair homolog 1) [40] and DCC1 (Defective in sister chromatid cohesion homolog 1) [41].

Table 1 Biological processes identified by PANTHER [23] as being significantly over-represented (p < 0.05) in the 157-transcript sex-specific SAM-list compared to the NCBI human reference gene list.
Figure 4
figure 4

Gene Set Enrichment Analysis. Enrichment plot is shown for a set of upregulated genes of a transcriptional profile commonly activated in undifferentiated cancer ("CANCER_UNDIFFERENTIATED_META_UP" as compared to differentiated cancer [22]). The enrichment score (ES, green line) reflects the degree to which the gene set is over-represented at the top or bottom of the ranked list of genes. Black bars illustrate the position of genes belonging to the gene set in the ranked list of genes included in the analysis. The ranked list metric shown in gray measures a gene's correlation with a phenotype. A positive value indicates correlation with "normal female "- phenotype, a negative value with "normal male" – phenotype.

Discussion

Differences in gene expression that correlate with gender can be dissected into a small set of sex chromosome genes and a larger set of genes that appear to be programmed by androgens, most likely during embryonic development. PBMCs from individuals with diverse well-defined etiologies of DSD including gonadal dysgenesis, defects of androgen biosynthesis, and high levels of androgenic steroids due to 21-hydroxylase deficiency showed strikingly consistent patterns of gene expression across 146 transcripts that correlated most strongly with the genital phenotype. Expression across all 157 transcripts identified as differentially expressed between normal males and females allowed separation of PBMCs into four categories of "PBMC-sexes": chromosomal and phenotypical male, chromosomal and phenotypical female and two intersex constellations in which the chromosomal and steroid influenced transcripts were discordant. While it is possible that transcript profiling of a higher number of DSD-patients could lead to demarcation of even diagnosis-specific patient clusters, the consistency of the transcript profiles across the different diagnoses argues strongly that these four categories are relatively robust.

Animal studies on rodents underlined the importance of sex specificity of the pattern of GH secretion and GH signaling involving STAT5b for maintenance of a large part of sexual dimorphism of gene expression in the liver [42]. Since human PBMC express the GH receptor [43], this mechanism could have a potential influence on the genes detected in the present study. However, in contrast to the extensive sex specific differences in rats with males showing pulsatile GH secretion and females showing more continuous GH levels, data on sex specificity of GH secretion in humans is less clear. While age and puberty status but not sex are important determinants for GH secretion during growth and puberty [44], there is other data showing that GH pulse amplitude but not pulse frequency differ between women and men in adulthood [45]. However, looking at the data of the current study, the correlation between gene expression profiles and genital masculinization was independent of the individuals' sex hormone levels at the time of expression profiling. Puberty status, presence of gonads and sex hormone replacement did not influence expression across the sex-specific transcripts. Additional GSEA on transcription factor targets gene sets did not reveal sex specific enrichment of STAT5b target genes in our study (data not shown). Therefore, we have no experimental evidence that sex specific differences of GH secretion or changes of GH secretion due to age or puberty could have influenced the sex specific gene set of the PBMC of our study. However, it cannot be excluded that a considerably larger set of investigated patients would allow the identification of some sex-dimorphic genes influenced by the GH signaling pathway in PBMCs.

Genital masculinization is a direct consequence of early embryonic androgen action between the 7th and 12th week of gestation. Our findings therefore support the existence of lasting programming of sex-specific genes implemented by presence or absence of androgen during the first trimester of embryogenesis. Sex-specific gene expression patterns have been observed in other tissues [8–11, 46–48], yet the relative contribution of the sex chromosomes and androgen programming was unknown. Based on expression patterns readily apparent in PBMCs, cells not usually thought to differ significantly between males and females, our data suggest that many of the sex-specific differences in gene expression are due to the influence of presence or absence of androgens during embryonic development. Given the limited lifespan of PBMCs, the lasting changes in transcript programs due to androgens are likely embedded in progenitor stem cells of fetal leukopoesis. It cannot be excluded that some of the genes detected in our study also played functional roles in androgen programming of lasting PBMC expression patterns during embryogenesis. However, more research is needed to understand the molecular mechanisms of long term hormonal programming, e.g., by studying the effects of androgens on embryonic stem cell differentiation [49]. Whether such programming affects the stem cells of other tissue types is unknown, but highly likely.

In a previous genome-wide assessment of transcript profiles, we have identified transcripts differentially expressed in cultured genital fibroblasts from normal males and 46,XY females with androgen insensitivity syndrome due to inactivating androgen receptor mutations [11]. Similar to the current study we identified a set of transcripts with reproducible differences in gene expression that correlated with the degree of genital masculinization. Since subjects in those studies all had a 46,XY karyotype, the differentially expressed transcripts were almost certainly programmed during development and remained apparent despite serial passage of the fibroblasts in culture. Interestingly, the overlap was restricted to 7 genes. In both datasets FZD6 (Frizzled 6) (Figure 1G), SSX2IP, SDC1, ALKBH, SPRED1 (Sprouty-related, EVH1 domain containing 1), PYCR2 (Pyrroline-5-carboxylate reductase family, member 2) and TOMM40 (Translocase of outer mitochondrial membrane 40) were differentially expressed between the subjects with male and female genitalia. Recent work by Yang and coworkers [9] in which 334 male and female mice were profiled confirms our findings that a significant number of genes differ between males and females. More importantly, as we observed in genital fibroblasts and PBMCs, the degree of overlap of the sex-specific genes between tissues was limited in this study. Our work extends on that of Yang et al by demonstrating that sex steroid induced programming accounts for the majority of the differences in gene expression between phenotypic males and females. Taken together these data suggest that the sex chromosomes to a small degree and androgen signaling during development to a larger degree account for a large part of the differences in gene expression in different tissue types between males and females. Furthermore, the influences of each of these differ significantly between tissue types as evidenced by the limited overlap in the sex-specific transcripts between tissue types. Additional work will be necessary to test whether these differences in gene expression between analogous tissues between males and females underlie developmental, anatomical, and functional differences between the sexes.

Our data also provide a framework for investigating differences in diseases between the sexes and could have future implications for risk assessment and disease monitoring [7]. PBMC expression profiling has been suggested as a diagnostic tool in that germline mutations in disease-causing genes could produce unique transcriptional signatures detectable in blood cells [50–52]. Disease-linked expression profiles are likely to be modified by biological sex [8, 9] and might need to be accounted for when developing diagnostic tools. In addition, our data raise the possibility that lasting transcriptional programs that differ between males and females could underlie differential susceptibility to a variety of diseases. For instance, we observed a striking over-representation of genes involved in growth control and cancer in PBMCs between males and females. Whether sex-specific gene expression could affect a tissue's transcriptional background and oncogenic potential or response to anti-cancer therapy is unknown, but merits consideration in light of differences in cancer incidence between the sexes and our new findings.

Our data have relevant implications for understanding the biology of DSD as well as potential future applications in its management. In Western societies, every schoolchild learns that "sex chromosomes" decide "who's a boy and who's a girl". Our data broaden the view of what determines sex and gender by revealing distinct male and female transcript profiles that correlate with either the karyotype or the androgen milieu present during embryonic development. Therefore, our data confirm the concept of long term androgen programming as established in rodent models [12–14] for the first time on the level of the human PBMC-transcriptome. Our data also reveal 2 categories of gene expression in PBMCs in individuals with DSD in which the karyotypic and androgen-programmed transcript profiles are discordant. Based on our findings, it is possible that other transcript profiles in PBMC could be identified that correlate with sex or gender specific traits outside the genitalia. Were that the case, it might be possible to specify transcript signatures that reflect sex-specific differences in tissues such as the brain and in turn mirror sex-specific behavior, sexual orientation or gender identity. If so, subsets of the sex specific genes of our study could serve as future transcriptional biomarkers in DSD outcome studies in order to develop novel diagnostic tools contributing to decisions on gender assignment in DSD children.

Conclusion

The sex chromosomes and the long term effects of the early embryonic androgen milieu as evidenced by the degree of genital masculinization in defined patients with DSD contribute independently to sex specific gene expression in PBMCs. This gives rise to 4 transcriptional categories of PBMC sexes: normal male, normal female and two "intersex" conditions with discordant alignments of the DSD individuals' genetic and hormonal sex signatures. This previously unknown architecture of sex specific gene expression in the human provides new perspectives for understanding the biology of sexual dimorphism in normal anatomy, physiology and behavior, its particular influences on diseases and its specific alterations in DSD.

Methods

The study was approved by the ethical committee of the University of Lübeck, Germany.

Normal controls and DSD patients

Blood samples were obtained from 34 individuals including 9 normal 46,XY males, 10 normal 46,XX females, one 46,XY male having CAH due to 21-hydroxylase deficiency and 14 individuals with DSD (Additional file 1). Mean age of normal males and normal females was 33 years, and 26.2 years, respectively (not different by t-test). Informed consent was obtained from control subjects, affected individuals, or their parents.

RNA-isolation, RNA-amplification and aRNA-labeling

Blood samples were taken at 1:00 pm +/- 1 h to control for circadian influences on gene expression. White blood cell counts did not differ between normal males and females (by t-test). PBMC were isolated by Ficoll (Biochrom, Berlin, Germany) gradient centrifugation and total RNA was extracted by TRIZOL (Invitrogen, Karlsruhe, Germany). 3 μg of total RNA and 1 μg of T7-oligo(dT)15 were used for linear RNA-amplification. Reverse transcription was done by SuperScript II (400 U) (Invitrogen, Karlsruhe, Germany) at 42°C for 90 min. Second-strand synthesis buffer containing RNaseH (2 U), DNA-Polymerase-I (40 U), and dNTPs (1 mM) (Roche, Mannheim, Germany) was mixed with first reaction mixture and incubated for 2 h at 15°C. In-vitro-transcription-mixture containing T7-RNA-Polymerase (60 U) and NTPs (1 mM) (Roche) was added and incubated for 5 h at 37°C. 3 μg PBMC aRNA, 7.5 μg random-hexamer-primers (Roche), and Cy5-dUTP (0.1 mM) (Amersham, Freiburg, Germany) were used for Cy5-labeling, and 50 μg of human universal reference-RNA (Stratagene; Amsterdam; The Netherlands), 4.5 μg T20VN-primer (TIB MolBiol, Berlin, Germany), and Cy3-dUTP (0,1 mM) (Amersham) for Cy3-labeling.

Microarrays and data analysis

DNA microarrays were obtained from Stanford Functional Genomics Facility http://www.microarray.org/sfgf/, hybridized according to established procedures http://cmgm.stanford.edu/pbrown/, scanned with GenePix4000B scanner, analyzed with GenePixPro 4.1 (Axon Instruments, Inc., Union City, CA) and loaded into Stanford Microarray Database http://genome-www5.stanford.edu/[53]. Gene filtering, log2 transformation, microarray normalization and centering of fluorescence ratios were performed as described previously [11]. Genes with significant differences in transcription between normal males and females were identified by Significance Analysis of Microarrays (SAM) [20] and further analyzed by hierarchical clustering analyses [54]. The one 46,XY CAH individual was not part of the group of 9 normal males considered for SAM.

Functional categorization of sex-specific genes

PANTHER was used to identify enrichment of biological processes with p-values < 0.05 considered as significant (http://www.pantherdb.org/[23]). GSEA was carried out with the Molecular Signature Database (http://www.broad.mit.edu/gsea[24]) using gene sets with at least 15 and maximal 500 genes. A p-value < 0.01 and a false discovery rate < 0.05 were considered significant.

RT-PCR

Quantitative RT-PCR of XIST, DDX3Y and FZD6 was performed on two 46,XY male controls (N-003, N-015), two 46,XX female controls (N-020, N-031), two 46,XY DSD-females (P-019, P-033) and two virilized 46,XX CAH-individuals (P-103, P-105). Probes were FAM (F) labeled with dabcyl (DB) quencher. Exon-spanning primers and probes (TIB Molbiol, Berlin, Germany) had the following sequences: XIST-forward: ACCACCACACgTCAAgCTCTT; XIST-reverse: CTATCCTCAAgTgCTAgAgTgCCAg, XIST labeled probe: 6FAM-TTCCTACAAgCAgTgCAgAgAgCTgAgT-DB; DDX3Y-forward: gAAAAACAgAgTggAggAgCAAgTA; DDX3Y-reverse: gTCCACgATCATCAAATCTTCCC, DDX3Y labeled probe: 6FAM-ATAAAgACAgTTCAggTTggAgTTgCAgCA-DB; FZD6-forward: CgTCAgTACCATATCCCATgTCCTTA, FZD6-reverse: gAAATgACCTTCAgCTTgTgTgAAC, FZD6 labeled probe: F-AAAgCAAAAgCTCgACCAgAATTggCT-DB; TATA-binding protein (TBP) was used to normalize expression values (TBP-forward: CACgAACCACggCACTg-ATT, TBP-reverse: TTTTCTTgCTgCCAgTCTggAC; TBP-labeled probe: 6FAM-TgTgCAC-AggAgCCAAgAgTgAAgA-DB). 2 μg of total RNA, 0.5 μM dNTPs, 5 μM random decamers and 200 U of MMLV-RT polymerase were used per reverse transcription reaction. Heat denaturation for 3 min at 85°C was followed by reverse transcription for 1 hr at 44°C and subsequent inactivation for 10 min at 92°C. C-DNA corresponding to 12 ng of initial total RNA served as template. Primers (0.4 μM), probe (0.17 μM) and TaqMan Universal PCR mastermix (Applied Biosystems, Foster City, CA, USA) including Taq polymerase were added (final volume 25 μl). Initial denaturation at 95°C for 10 min was followed by 1 min cycling intervals at 60°C using RotorGene RG-3000 cycler (Corbett-Research, Sydney, Australia). Differential transcription was calculated by the ΔΔCT method [55].

Abbreviations

DSD:

disorders of sex development

CAH:

congenital adrenal hyperplasia

SMD:

Stanford Microarray Database

SAM:

Significance Analysis of Microarrays

GSEA:

Gene Set Enrichment Analysis

PANTHER:

Protein ANalysis THrough Evolutionary Relationships.

References

  1. Berenbaum SA: Effects of early androgens on sex-typed activities and interests in adolescents with congenital adrenal hyperplasia. Horm Behav. 1999, 35: 102-110. 10.1006/hbeh.1998.1503.

    Article  CAS  PubMed  Google Scholar 

  2. Hines M, Ahmed SF, Hughes IA: Psychological outcomes and gender-related development in complete androgen insensitivity syndrome. Arch Sex Behav. 2003, 32: 93-101. 10.1023/A:1022492106974.

    Article  PubMed  Google Scholar 

  3. Jürgensen M, Hiort O, Holterhus PM, Thyen U: Gender role behavior in children with XY karyotype and disorders of sex development. Horm Behav. 2007, 51: 443-453. 10.1016/j.yhbeh.2007.01.001.

    Article  PubMed  Google Scholar 

  4. Decaudin D, Salanoubat C, Carde P: Is mantle cell lymphoma a sex-related disease?. Leuk Lymphoma. 2000, 37: 181-184.

    Article  CAS  PubMed  Google Scholar 

  5. Boerma EG, van Imhoff GW, Appel IM, Veeger NJ, Kluin PM, Kluin-Nelemans JC: Gender and age-related differences in Burkitt lymphoma – epidemiological and clinical data from The Netherlands. Eur J Cancer. 2004, 40: 2781-2787. 10.1016/j.ejca.2004.09.004.

    Article  CAS  PubMed  Google Scholar 

  6. Soto ME, Vallejo M, Guillén F, Simón JA, Arena E, Reyes PA: Gender impact in systemic lupus erythematosus. Clin Exp Rheumatol. 2004, 22: 713-721.

    CAS  PubMed  Google Scholar 

  7. Ober C, Loisel DA, Gilad Y: Sex-specific genetic architecture of human disease. Nature Reviews Genetics. 2008, 9: 911-922. 10.1038/nrg2415.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  8. Whitney AR, Diehn M, Popper SJ, Alizadeh AA, Boldrick JC, Relman DA, Brown PO: Individuality and variation in gene expression patterns in human blood. Proc Natl Acad Sci USA. 2003, 100: 1896-1901. 10.1073/pnas.252784499.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  9. Yang X, Schadt EE, Wang S, Wang H, Arnold AP, Ingram-Drake L, Drake TA, Lusis AJ: Tissue-specific expression and regulation of sexually dimorphic genes in mice. Genome Res. 2006, 16: 995-1004. 10.1101/gr.5217506.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  10. Holterhus PM, Hiort O, Demeter J, Brown PO, Brooks JD: Differential gene-expression patterns in genital fibroblasts of normal males and 46,XY females with androgen insensitivity syndrome: evidence for early programming involving the androgen receptor. Genome Biol. 2003, 4: R37-10.1186/gb-2003-4-6-r37.

    Article  PubMed Central  PubMed  Google Scholar 

  11. Holterhus PM, Deppe U, Werner R, Richter-Unruh A, Bebermeier JH, Wünsch L, Krege S, Schweikert HU, Demeter J, Riepe F, Hiort O, Brooks JD: Intrinsic androgen-dependent gene expression patterns revealed by comparison of genital fibroblasts from normal males and individuals with complete and partial androgen insensitivity syndrome. BMC Genomics. 2007, 8: 376-10.1186/1471-2164-8-376.

    Article  PubMed Central  PubMed  Google Scholar 

  12. Döhler KD, Coquelin A, Davis F, Hines M, Shryne JE, Gorski RA: Differentiation of the sexually dimorphic nucleus in the preoptic area of the rat brain is determined by the perinatal hormone environment. Neurosci Lett. 1982, 33: 295-298. 10.1016/0304-3940(82)90388-3.

    Article  PubMed  Google Scholar 

  13. Roy KA, Chatterjee B: Sexual dimorphism in the liver. Ann Rev Physiol. 1983, 45: 37-50. 10.1146/annurev.ph.45.030183.000345.

    Article  CAS  Google Scholar 

  14. Singh J, Handelsman DJ: Imprinting by neonatal sex steroids on the structure and function of the mature mouse prostate. Biol Reprod. 1999, 61: 200-208. 10.1095/biolreprod61.1.200.

    Article  CAS  PubMed  Google Scholar 

  15. Hrabovszky Z, Hutson JM: Androgen imprinting of the brain in animal models and humans with intersex disorders: review and recommendations. J Urol. 2002, 168: 2142-2148. 10.1016/S0022-5347(05)64338-8.

    Article  CAS  PubMed  Google Scholar 

  16. Hines M: Sex steroids and human behavior: prenatal androgen exposure and sex-typical play behavior in children. Ann N Y Acad Sci. 2003, 1007: 272-282. 10.1196/annals.1286.026.

    Article  CAS  PubMed  Google Scholar 

  17. Cohen-Kettenis PT: Gender change in 46,XY persons with 5alpha-reductase-2 deficiency and 17beta-hydroxysteroid dehydrogenase-3 deficiency. Arch Sex Behav. 2005, 34: 399-410. 10.1007/s10508-005-4339-4.

    Article  PubMed  Google Scholar 

  18. Rinn JL, Bondre C, Gladstone HB, Brown PO, Chang HY: Anatomic demarcation by positional variation in fibroblast gene expression programs. PLoS Genet. 2006, 2 (7): e119-10.1371/journal.pgen.0020119.

    Article  PubMed Central  PubMed  Google Scholar 

  19. Hoppe U, Holterhus PM, Wünsch L, Jocham D, Drechsler T, Thiele S, Marschke C, Hiort O: Tissue-specific transcription profiles of sex steroid biosynthesis enzymes and the androgen receptor. J Mol Med. 2006, 84: 651-659. 10.1007/s00109-006-0049-1.

    Article  CAS  PubMed  Google Scholar 

  20. Tusher VG, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci USA. 2001, 98: 5116-5121. 10.1073/pnas.091062498.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  21. Morgan E, Nguyen S, Scott V, Stadler H: Loss of Bmp7 and Fgf8 signaling in Hoxa13-mutant mice causes hypospadia. Development. 2003, 130: 3095-3109. 10.1242/dev.00530.

    Article  CAS  PubMed  Google Scholar 

  22. Köhler B, Delezoide AL, Boizet-Bonhoure B, McPhaul MJ, Sultan C, Lumbroso SJ: Coexpression of Wilms' tumor suppressor 1 (WT1) and androgen receptor (AR) in the genital tract of human male embryos and regulation of AR promoter activity by WT1. J Mol Endocrinol. 2007, 38: 547-554. 10.1677/JME-06-0020.

    Article  PubMed  Google Scholar 

  23. Thomas PD, Kejariwal A, Campbell MJ, Mi H, Diemer K, Guo N, Ladunga I, Ulitsky-Lazareva B, Muruganujan A, Rabkin S, Vandergriff JA, Doremieux O: PANTHER: a browsable database of gene products organized by biological function, using curated protein family and subfamily classification. Nucleic Acids Res. 2003, 31: 334-341. 10.1093/nar/gkg115.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  24. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP: Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005, 102: 15545-15550. 10.1073/pnas.0506580102.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  25. Rhodes DR, Yu J, Shanker K, Deshpande N, Varambally R, Ghosh D, Barrette T, Pandey A, Chinnaiyan AM: Large-scale meta-analysis of cancer microarray data identifies common transcriptional profiles of neoplastic transformation and progression. Proc Natl Acad Sci USA. 2004, 101: 9309-9314. 10.1073/pnas.0401994101.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  26. Becker J, Volland S, Noskova I, Schramm A, Schweigerer LL, Wilting J: Keratoepithelin reverts the suppression of tissue factor pathway inhibitor 2 by MYCN in human neuroblastoma: a mechanism to inhibit invasion. Int J Oncol. 2008, 32: 235-240.

    CAS  PubMed  Google Scholar 

  27. Nikitovic D, Berdiaki A, Zafiropoulos A, Katonis P, Tsatsakis A, Karamanos NK, Tzanakakis GN: Lumican expression is positively correlated with the differentiation and negatively with the growth of human osteosarcoma cells. FEBS J. 2008, 275: 350-361.

    Article  CAS  PubMed  Google Scholar 

  28. Kono T, Shimoda M, Takahashi M, Matsumoto K, Yoshimoto T, Mizutani M, Tabata C, Okoshi K, Wada H, Kubo H: Immunohistochemical detection of the lymphatic marker podoplanin in diverse types of human cancer cells using a novel antibody. Int J Oncol. 2007, 31: 501-508.

    CAS  PubMed  Google Scholar 

  29. Shah L, Walter KL, Borczuk AC, Kawut SM, Sonett JR, Gorenstein LA, Ginsburg ME, Steinglass KM, Powell CA: Expression of syndecan-1 and expression of epidermal growth factor receptor are associated with survival in patients with nonsmall cell lung carcinoma. Cancer. 2004, 101: 1632-1638. 10.1002/cncr.20542.

    Article  CAS  PubMed  Google Scholar 

  30. Ryu B, Kim DS, Deluca AM, Alani RM: Comprehensive expression profiling of tumor cell lines identifies molecular signatures of melanoma progression. PLoS ONE. 2007, 2: e594-10.1371/journal.pone.0000594.

    Article  PubMed Central  PubMed  Google Scholar 

  31. Breslin A, Denniss FA, Guinn BA: SSX2IP: an emerging role in cancer. Biochem Biophys Res Commun. 2007, 363: 462-465. 10.1016/j.bbrc.2007.09.052.

    Article  CAS  PubMed  Google Scholar 

  32. Croci S, Landuzzi L, Nicoletti G, Palladini A, Antognoli A, De Giovanni C, Nanni P, Lollini PL: Expression of Connective Tissue Growth Factor (CTGF/CCN2) in a Mouse Model of Rhabdomyosarcomagenesis. Pathol Oncol Res. 2007, 13: 336-339. 10.1007/BF02940313.

    Article  CAS  PubMed  Google Scholar 

  33. Kondratiev S, Gnepp DR, Yakirevich E, Sabo E, Annino DJ, Rebeiz E, Laver NV: Expression and prognostic role of MMP2, MMP9, MMP13, and MMP14 matrix metalloproteinases in sinonasal and oral malignant melanomas. Hum Pathol. 2008, 39: 337-343. 10.1016/j.humpath.2007.07.003.

    Article  CAS  PubMed  Google Scholar 

  34. Tao Y, Pinzi V, Bourhis J, Deutsch E: Mechanisms of disease: signaling of the insulin-like growth factor 1 receptor pathway – therapeutic perspectives in cancer. Nat Clin Pract Oncol. 2007, 4: 591-602. 10.1038/ncponc0934.

    Article  CAS  PubMed  Google Scholar 

  35. de Sá SV, Corrêa-Giannella ML, Machado MC, Krogh K, de Almeida MQ, Albergaria Pereira MA, Coelho Siqueira SA, Patzina RA, Ibuki FS, Sogayar MC, Machado MC, Giannella-Neto D: Serpin peptidase inhibitor clade A member 1 as a potential marker for malignancy in insulinomas. Clin Cancer Res. 2007, 13: 5322-5330. 10.1158/1078-0432.CCR-06-1477.

    Article  PubMed  Google Scholar 

  36. Ono R, Taki T, Taketani T, Kawaguchi H, Taniwaki M, Okamura T, Kawa K, Hanada R, Kobayashi M, Hayashi Y: SEPTIN6, a human homologue to mouse Septin6, is fused to MLL in infant acute myeloid leukemia with complex chromosomal abnormalities involving 11q23 and Xq24. Cancer Res. 2002, 62: 333-337.

    CAS  PubMed  Google Scholar 

  37. Wang XQ, Stanbridge EJ, Lao X, Cai Q, Fan ST, Redpath JL: p53-dependent Chk1 phosphorylation is required for maintenance of prolonged G2 Arrest. Radiat Res. 2007, 168: 706-715. 10.1667/RR1075.1.

    Article  CAS  PubMed  Google Scholar 

  38. Pakakasama S, Sirirat T, Kanchanachumpol S, Udomsubpayakul U, Mahasirimongkol S, Kitpoka P, Thithapandha A, Hongeng S: Genetic polymorphisms and haplotypes of DNA repair genes in childhood acute lymphoblastic leukemia. Pediatr Blood Cancer. 2007, 48: 16-20. 10.1002/pbc.20742.

    Article  PubMed  Google Scholar 

  39. Leduc F, Maquennehan V, Nkoma GB, Boissonneault G: DNA damage response during chromatin remodeling in elongating spermatids of mice. Biol Reprod. 2008, 78: 324-332. 10.1095/biolreprod.107.064162.

    Article  CAS  PubMed  Google Scholar 

  40. Sedgwick B: Repairing DNA-methylation damage. Nat Rev Mol Cell Biol. 2004, 5: 148-157. 10.1038/nrm1312.

    Article  CAS  PubMed  Google Scholar 

  41. Game JC, Williamson MS, Baccari C: X-ray survival characteristics and genetic analysis for nine Saccharomyces deletion mutants that show altered radiation sensitivity. Genetics. 2005, 169: 51-63. 10.1534/genetics.104.028613.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  42. Clodfelter KH, Holloway MG, Hodor P, Park SH, Ray WJ, Waxman DJ: Sex-dependent liver gene expression is extensive and largely dependent upon signal transducer and activator of transcription 5b (STAT5b): STAT5b-dependent activation of male genes and repression of female genes revealed by microarray analysis. Mol Endocrinol. 2006, 20: 1333-1351. 10.1210/me.2005-0489.

    Article  CAS  PubMed  Google Scholar 

  43. Valerio G, Bond HM, Badolato R, Petrella A, Di Maio S, Salerno M, Waters MJ, Venuta S, Tenore A: Expression of growth hormone receptor by peripheral blood lymphocytes in children: evaluation in clinical conditions of impaired growth. Clin Endocrinol. 1997, 47: 329-335. 10.1046/j.1365-2265.1997.2571066.x.

    Article  CAS  Google Scholar 

  44. Zadik Z, Chalew SA, McCarter RJ, Meistas M, Kowarski AA: The influence of age on the 24-hour integrated concentration of growth hormone in normal individuals. J Clin Endocrinol Metab. 1985, 60: 513-516.

    Article  CAS  PubMed  Google Scholar 

  45. Berg van den G, Veldhuis JD, Frölich M, Roelfsema F: An amplitude-specific divergence in the pulsatile mode of growth hormone (GH) secretion underlies the gender difference in mean GH concentrations in men and premenopausal women. J Clin Endocrinol Metab. 1996, 81: 2460-2467. 10.1210/jc.81.7.2460.

    PubMed  Google Scholar 

  46. Takemoto N, Iizuka N, Yamada-Okabe H, Hamada K, Tamesa T, Okada T, Hashimoto K, Sakamoto K, Takashima M, Miyamoto T, Uchimura S, Hamamoto Y, Oka M: Sex-based molecular profiling of hepatitis C virus-related hepatocellular carcinoma. Int J Oncol. 2005, 26: 673-678.

    CAS  PubMed  Google Scholar 

  47. Rinn JL, Rozowsky JS, Laurenzi IJ, Petersen PH, Zou K, Zhong W, Gerstein M, Snyder M: Major molecular differences between mammalian sexes are involved in drug metabolism and renal function. Dev Cell. 2004, 6: 791-800. 10.1016/j.devcel.2004.05.005.

    Article  CAS  PubMed  Google Scholar 

  48. Dewing P, Shi T, Horvath S, Vilain E: Sexually dimorphic gene expression in mouse brain precedes gonadal differentiation. Brain Res Mol Brain Res. 2003, 118: 82-90. 10.1016/S0169-328X(03)00339-5.

    Article  CAS  PubMed  Google Scholar 

  49. Ray R, Novotny NM, Crisostomo PR, Lahm T, Abarbanell A, Meldrum DR: Sex steroids and stem cell function. Mol Med. 2008, 14: 493-501. 10.2119/2008-00004.Ray.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  50. Cario G, Stanulla M, Fine BM, Teuffel O, Neuhoff NV, Schrauder A, Flohr T, Schäfer BW, Bartram CR, Welte K, Schlegelberger B, Schrappe M: Distinct gene expression profiles determine molecular treatment response in childhood acute lymphoblastic leukemia. Blood. 2005, 105: 821-826. 10.1182/blood-2004-04-1552.

    Article  CAS  PubMed  Google Scholar 

  51. Batliwalla FM, Baechler EC, Xiao X, Li W, Balasubramanian S, Khalili H, Damle A, Ortmann WA, Perrone A, Kantor AB, Gulko PS, Kern M, Furie R, Behrens TW, Gregersen PK: Peripheral blood gene expression profiling in rheumatoid arthritis. Genes Immun. 2005, 6: 388-397. 10.1038/sj.gene.6364209.

    Article  CAS  PubMed  Google Scholar 

  52. Sharp FR, Xu H, Lit L, Walker W, Apperson M, Gilbert DL, Glauser TA, Wong B, Hershey A, Liu DZ, Pinter J, Zhan X, Liu X, Ran R: The future of genomic profiling of neurological diseases using blood. Arch Neurol. 2006, 63: 1529-1536. 10.1001/archneur.63.11.1529.

    Article  PubMed  Google Scholar 

  53. Sherlock G, Hernandez-Boussard T, Kasarskis A, Binkley G, Matese JC, Dwight SS, Kaloper M, Weng S, Jin H, Ball CA, Eisen MB, Spellman PT, Brown PO, Botstein D, Cherry JM: The Stanford Microarray Database. Nucleic Acids Res. 2001, 29: 152-155. 10.1093/nar/29.1.152.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  54. Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci USA. 1998, 95: 14863-14868. 10.1073/pnas.95.25.14863.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  55. Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-delta delta C(T)) method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

The study was funded by the Deutsche Forschungsgemeinschaft (DFG) (KFO111/1-D, Ho2073/6-1/2 to PMH) and by EuroDSD (in the European Community's Seventh Framework Programme FP7/2007–2013 under grant agreement n° 201444 (WP04 to PMH and RS). JDB was supported by NIH CA111782. We are grateful to the scientists and staff of the Stanford Functional Genomics Facility and the Stanford Microarray Database without whom the project could not have been realized.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Paul-Martin Holterhus.

Additional information

Authors' contributions

PMH was principal investigator, designed the study, performed clinical evaluation of DSD patients, microarray data analysis, and he was the leading writer. JHB was responsible for technical design of experiments, RNA-preparations, microarray hybridizations, database work, data analysis, and he was a writer. RW did genomic DNA sequence analysis in 46,XY DSD patients, RT-PCR, and was also a writer. JD is the leading database manager, performed data analysis and contributed to writing. ARU evaluated DSD patients, performed genomic DNA sequence analysis of CYP21A2, data interpretation and she was a co-writer. GC undertook the gene set enrichment analyses, data interpretation related to growth control and proliferation and wrote respective passages. MA contributed to microarray works and data analysis. RS performed data analysis and data interpretation. FR performed clinical evaluations of DSD patients and provided expert advice concerning genotype-phenotype correlation of 21-hydroxylase deficiency and other defects of steroid biosynthesis causing DSD. JDB is one of the senior writers and the senior advisor of microarray data interpretation, OH is senior writer, performed clinical evaluation of DSD patients, genomic DNA-sequence analysis of 46,XY DSD patients. He is senior advisor correlating clinical data, hormonal data and gene expression data. All authors read and approved the final manuscript.

Electronic supplementary material

12864_2009_2176_MOESM1_ESM.pdf

Additional file 1:Phenotypic, clinical and molecular background information on the control individuals and the DSD patients corresponding to the PBMC sample numbers as indicated in the first column.(PDF 28 KB)

12864_2009_2176_MOESM2_ESM.txt

Additional file 2: Text file to visualize data for Figure 1 together with additional files 3, 4 and 5 in TreeView [54]. (TXT 52 KB)

12864_2009_2176_MOESM3_ESM.zip

Additional file 3: Zip file containing cdt file. Cluster file to visualize data for Figure 1 together with additional files 2, 4 and 5 in TreeView [54]. (ZIP 24 KB)

12864_2009_2176_MOESM4_ESM.zip

Additional file 4: Zip file containing atr file. Array tree file to visualize data for Figure 1 together with additional files 2, 3 and 5 in TreeView[54] (ZIP 539 bytes)

12864_2009_2176_MOESM5_ESM.zip

Additional file 5: Zip file containing gtr file. Gene tree file to visualize data for Figure 1 together with additional files 2, 3 and 4 in TreeView [54]. (ZIP 2 KB)

12864_2009_2176_MOESM6_ESM.txt

Additional file 6: Text file to visualize data for Figure 3 together with additional files 7, 8 and 9 in TreeView [54]. (TXT 58 KB)

12864_2009_2176_MOESM7_ESM.zip

Additional file 7: Zip file containing cdt file. Cluster file to visualize data for Figure 3 together with additional files 6, 8 and 9 in TreeView [54]. (ZIP 24 KB)

12864_2009_2176_MOESM8_ESM.zip

Additional file 8: Zip file containing atr file. Array tree file to visualize data for Figure 3 together with additional files 6, 7 and 9 in TreeView [54]. (ZIP 549 bytes)

12864_2009_2176_MOESM9_ESM.zip

Additional file 9: Zip file containing a gtr file. Gene tree file to visualize data for Figure 3 together with additional files 6, 7 and 8 in TreeView [54]. (ZIP 2 KB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Holterhus, PM., Bebermeier, JH., Werner, R. et al. Disorders of sex development expose transcriptional autonomy of genetic sex and androgen-programmed hormonal sex in human blood leukocytes. BMC Genomics 10, 292 (2009). https://doi.org/10.1186/1471-2164-10-292

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-10-292

Keywords