Whole-genome Microarray Research Articles

Benign breast tissue composition in breast cancer patients: association with risk factors, clinical variables, and gene expression.

Breast tissue composition (epithelium, non-fatty stroma, and adipose) changes qualitatively and quantitatively throughout the lifespan, and may mediate relationships between risk factors and breast cancer initiation. We sought to identify relationships between tissue composition, risk factors, tumor characteristics, and gene expression. Participants were 146 patients from the Polish Breast Cancer Study, with data on risk factor and clinicopathological characteristics. Benign breast tissue composition was evaluated using digital image analysis of histologic sections. Whole-genome microarrays were performed on the same tissue blocks. Mean epithelial, non-fatty stromal, and adipose proportions were 8.4% (SD = 4.9%), 27.7% (SD = 24.0%), and 64.0% (SD = 24.0%), respectively. Among women <50 years old, stroma proportion decreased and adipose proportion increased with age, with approximately 2% difference per year (P < 0.01). The variation in epithelial proportion with age was modest (0.1% per year). Higher epithelial proportion was associated with obesity (7.6% in nonobese vs. 10.1% in obese; P = 0.02) and with poorly differentiated tumors (7.8% in well/moderate vs. 9.9% in poor; P = 0.05). Gene expression signatures associated with epithelial and stromal proportion were identified and validated. Stroma-associated genes were in metabolism and stem cell maintenance pathways, whereas epithelial genes were enriched for cytokine and immune response pathways. Breast tissue composition was associated with age, body mass index, and tumor grade, with consequences for breast gene expression. Breast tissue morphologic factors may influence breast cancer etiology. Composition and gene expression may act as biomarkers of breast cancer risk and progression.

Activation of the JAK/STAT pathway in Behcet's disease.

Th1/Th17-type T-cell responses are upregulated in Behcet's disease (BD). However, signaling pathways associated with this aberrant immune response are not clarified. Whole-genome microarray profiling was performed with human U133 (Plus 2.0) chips using messenger RNA of isolated CD14(+) monocytes and CD4(+) T cells from peripheral blood mononucleated cell (PBMC) in patients with BD (n = 9) and healthy controls (HCs) (n = 9). Flow cytometric analysis of unstimulated (US) and stimulated (phytohaemagglutinin) signal transducer and activator of transcription (STAT3) and pSTAT3 expressions of PBMCs were also analyzed (BD and HC, both n = 26). Janus family of kinase (JAK1) was observed to be upregulated in both CD14(+) monocytes (1.95-fold) and CD4(+) T lymphocytes (1.40-fold) of BD patients. Using canonical pathway enrichment analysis, JAK/STAT signaling was identified as activated in both CD14(+) monocytes (P = 9.55E-03) and in CD4(+) lymphocytes (P =8.13E-04) in BD. Interferon signaling was also prominent among upregulated genes in CD14(+) monocytes (P = 5.62E-05). Glucocorticoid receptor signaling and interleukin (IL-6) signaling were among the most enriched pathways in differentially expressed genes in CD14+ monocytes (P = 2.45E-09 and 1.00E-06, respectively). Basal US total STAT3 expression was significantly higher in BD (1.2 vs 3.45, P < 0.05). The JAK1/STAT3 signaling pathway is activated in BD, possibly through the activation of Th1/Th17-type cytokines such as IL-2, interferon (IFN-γ), IL-6, IL-17 and IL-23.

Carbon nanotubes induce apoptosis resistance of human lung epithelial cells through FLICE-inhibitory protein.

Chronic exposure to single-walled carbon nanotubes (SWCNT) has been reported to induce apoptosis resistance of human lung epithelial cells. As resistance to apoptosis is a foundation of neoplastic transformation and cancer development, we evaluated the apoptosis resistance characteristic of the exposed lung cells to understand the pathogenesis mechanism. Passage control and SWCNT-transformed human lung epithelial cells were treated with known inducers of apoptosis via the intrinsic (antimycin A and CDDP) or extrinsic (FasL and TNF-α) pathway and analyzed for apoptosis by DNA fragmentation, annexin-V expression, and caspase activation assays. Whole-genome microarray was performed to aid the analysis of apoptotic gene signaling network. The SWCNT-transformed cells exhibited defective death receptor pathway in association with cellular FLICE-inhibitory protein (c-FLIP) overexpression. Knockdown or chemical inhibition of c-FLIP abrogated the apoptosis resistance of SWCNT-transformed cells. Whole-genome expression signature analysis confirmed these findings. This study is the first to demonstrate carbon nanotube-induced defective death receptor pathway and the role of c-FLIP in the process.

Current transcriptomics in pig immunity research.

Swine performance in the face of disease challenge is becoming progressively more important. To improve the pig’s robustness and resilience against pathogens through selection, a better understanding of the genetic and epigenetic factors in the immune response is required. This review highlights results from the most recent transcriptome research, and the meta-analyses performed, in the context of pig immunity. A technological overview is given including wholegenome microarrays, immune-specific arrays, small-scale high-throughput expression methods, high-density tiling arrays, and next generation sequencing (NGS). Although whole genome microarray techniques will remain complementary to NGS for some time in domestic species, research will transition to sequencing-based methods due to cost-effectiveness and the extra information that such methods provide. Furthermore, upcoming high-throughput epigenomic studies, which will add greatly to our knowledge concerning the impact of epigenetic modifications on pig immune response, are listed in this review. With emphasis on the insights obtained from transcriptomic analyses for porcine immunity, we also discuss the experimental design in pig immunity research and the value of the newly published porcine genome assembly in using the pig as a model for human immune response. We conclude by discussing the importance of establishing community standards to maximize the possibility of integrative computational analyses, such as was clearly beneficial for the human ENCODE project.Electronic supplementary materialThe online version of this article (doi:10.1007/s00335-014-9549-4) contains supplementary material, which is available to authorized users.

Transcriptome response of Lactobacillus sakei to meat protein environment.

Lactobacillus sakei is a heterofermentative species of lactic acid bacteria that is used in industrial meat fermentation. To investigate adaptation in a meat environment, whole-genome DNA microarrays were used to analyze the gene expression related to growth and survival of L. sakei strain La22 when grown in sarcoplasmic (S-) or myofibrillar (M-) protein-supplemented chemically defined medium (CDM). Differential expression was detected in 551 genes. Genes encoding enzymes involved in peptide hydrolysis were differentially upregulated in M-CDM or/and S-CDM, and only oppB and oppC, involved in the amino acid and peptide transport system, were upregulated. Most genes related to metabolism of peptides, amino acids and related molecules were over-expressed in M-CDM and S-CDM, except for glnA and metK. Expression of certain genes was according to the differential substrate environment. The expression of genes involved in the stress response was not induced by growth in M-CDM.

Assessing copy number from exome sequencing and exome array CGH based on CNV spectrum in a large clinical cohort.

Detection of copy-number variation (CNV) is important for investigating many genetic disorders. Testing a large clinical cohort by array comparative genomic hybridization provides a deep perspective on the spectrum of pathogenic CNV. In this context, we describe a bioinformatics approach to extract CNV information from whole-exome sequencing and demonstrate its utility in clinical testing. Exon-focused arrays and whole-genome chromosomal microarray analysis were used to test 14,228 and 14,000 individuals, respectively. Based on these results, we developed an algorithm to detect deletions/duplications in whole-exome sequencing data and a novel whole-exome array. In the exon array cohort, we observed a positive detection rate of 2.4% (25 duplications, 318 deletions), of which 39% involved one or two exons. Chromosomal microarray analysis identified 3,345 CNVs affecting single genes (18%). We demonstrate that our whole-exome sequencing algorithm resolves CNVs of three or more exons. These results demonstrate the clinical utility of single-exon resolution in CNV assays. Our whole-exome sequencing algorithm approaches this resolution but is complemented by a whole-exome array to unambiguously identify intragenic CNVs and single-exon changes. These data illustrate the next advancements in CNV analysis through whole-exome sequencing and whole-exome array.Genet Med 17 8, 623-629.

Microduplication of chromosome Xq25 encompassing STAG2 gene in a boy with intellectual disability

Whole-genome microarray analysis is proven to be useful in the identification of submicroscopic copy number imbalances in families with intellectual disabilities. The first case of Xq25 duplication was identified using genome-wide array comparative genomic hybridization (array-CGH) in a 24-year-old patient with a syndromic intellectual disability. We report a 4-year-old boy with a de novo 591 kb duplication at Xq25. The duplication was first detected by a CytoScan HD array platform (Affymetrix, USA) and was confirmed by real-time quantitative PCR (qPCR) of the STAG2 gene, and by fluorescence in situ hybridization (FISH). The patient had clinical features partially consistent with published cases, including an intellectual disability and speech delay. The identification of this additional patient and a detailed analysis of duplications identified in other patient cohorts and absent in normal individuals support the existence of a rare pathological microduplication at Xq25 that encompasses STAG2.

Populus tremula (European aspen) shows no evidence of sexual dimorphism.

BackgroundEvolutionary theory suggests that males and females may evolve sexually dimorphic phenotypic and biochemical traits concordant with each sex having different optimal strategies of resource investment to maximise reproductive success and fitness. Such sexual dimorphism would result in sex biased gene expression patterns in non-floral organs for autosomal genes associated with the control and development of such phenotypic traits.ResultsWe examined morphological, biochemical and herbivory traits to test for sexually dimorphic resource allocation strategies within collections of sexually mature and immature Populus tremula (European aspen) trees. In addition we profiled gene expression in mature leaves of sexually mature wild trees using whole-genome oligonucleotide microarrays and RNA-Sequencing.ConclusionsWe found no evidence of sexual dimorphism or differential resource investment strategies between males and females in either sexually immature or mature trees. Similarly, single-gene differential expression and machine learning approaches revealed no evidence of large-scale sex biased gene expression. However, two significantly differentially expressed genes were identified from the RNA-Seq data, one of which is a robust diagnostic marker of sex in P. tremula.Electronic supplementary materialThe online version of this article (doi:10.1186/s12870-014-0276-5) contains supplementary material, which is available to authorized users.

Expression profile of the Schistosoma japonicum degradome reveals differential protease expression patterns and potential anti-schistosomal intervention targets.

Blood fluke proteases play pivotal roles in the processes of invasion, nutrition acquisition, immune evasion, and other host-parasite interactions. Hundreds of genes encoding putative proteases have been identified in the recently published schistosome genomes. However, the expression profiles of these proteases in Schistosoma species have not yet been systematically analyzed. We retrieved and culled the redundant protease sequences of Schistosoma japonicum, Schistosoma mansoni, Echinococcus multilocularis, and Clonorchis sinensis from public databases utilizing bioinformatic approaches. The degradomes of the four parasitic organisms and Homo sapiens were then comparatively analyzed. A total of 262 S. japonicum protease sequences were obtained and the expression profiles generated using whole-genome microarray. Four main clusters of protease genes with different expression patterns were identified: proteases up-regulated in hepatic schistosomula and adult worms, egg-specific or predominantly expressed proteases, cercaria-specific or predominantly expressed proteases, and constantly expressed proteases. A subset of protease genes with different expression patterns were further validated using real-time quantitative PCR. The present study represents the most comprehensive analysis of a degradome in Schistosoma species to date. These results provide a firm foundation for future research on the specific function(s) of individual proteases and may help to refine anti-proteolytic strategies in blood flukes.

Claudin-18 deficiency results in alveolar barrier dysfunction and impaired alveologenesis in mice.

Claudins are a family of transmembrane proteins that are required for tight junction formation. Claudin (CLDN)-18.1, the only known lung-specific tight junction protein, is the most abundant claudin in alveolar epithelial type (AT) 1 cells, and is regulated by lung maturational agonists and inflammatory mediators. To determine the function of CLDN18 in the alveolar epithelium, CLDN18 knockout (KO) mice were generated and studied by histological, biochemical, and physiological approaches, in addition to whole-genome microarray. Alveolar epithelial barrier function was assessed after knockdown of CLDN18 in isolated lung cells. CLDN18 levels were measured by quantitative PCR in lung samples from fetal and postnatal human infants. We found that CLDN18 deficiency impaired alveolar epithelial barrier function in vivo and in vitro, with evidence of increased paracellular permeability and architectural distortion at AT1-AT1 cell junctions. Although CLDN18 KO mice were born without evidence of a lung abnormality, histological and gene expression analysis at Postnatal Day 3 and Week 4 identified impaired alveolarization. CLDN18 KO mice also had evidence of postnatal lung injury, including acquired AT1 cell damage. Human fetal lungs at 23-24 weeks gestational age, the highest-risk period for developing bronchopulmonary dysplasia, a disease of impaired alveolarization, had significantly lower CLDN18 expression relative to postnatal lungs. Thus, CLDN18 deficiency results in epithelial barrier dysfunction, injury, and impaired alveolarization in mice. Low expression of CLDN18 in human fetal lungs supports further investigation into a role for this tight junction protein in bronchopulmonary dysplasia.

Abstract 5561: Clinical drug response can be predicted using baseline gene expression levels and in vitro drug sensitivity in cell lines

Abstract Robust prediction of in vivo chemotherapeutic response, using baseline gene expression and drug sensitivity data gathered on cancer cell lines, has been a profoundly important, long standing and controversial problem in pharmacogenomics. Here, we present for the first time, a solution to this problem. Currently, personalizing cancer chemotherapy relies on pathology and more recently molecular biomarker-based approaches (e.g. ERBB2 amplification in breast cancer). However, as the driving biology are normally not fully understood, the majority of existing biomarkers do not capture a substantial proportion of variability in drug response. This partly explains the commonly observed lack of reproducibility of findings (e.g. from many conventional gene expression signatures) when these markers are applied to new datasets. In this study, we developed an approach to predict in vivo drug sensitivity that leverages whole-genome gene expression microarray data and allows the expression of every gene to influence the prediction by a small amount. The method works by fitting a ridge regression model of baseline genome-wide gene expression levels against in vitro drug sensitivity in a very large panel of approximately 700 cancer cell lines. Then, after a (crucial) data homogenization step, these models are applied to baseline expression levels from primary tumor biopsies. Our method successfully predicted patient response to different chemotherapeutic agents in three (of four total suitable) independent, publicly available clinical trials, each investigating different drugs and different types of cancer. In each of these cases, we predicted drug response at least as accurately as previously published models that had been derived from the clinical data itself. Interestingly, our approach could also predict clinical response in the absence of any known drug sensitivity biomarker. We effectively enriched for drug responders in breast, myeloma and lung cancers, treated with docetaxel, bortezomib and erlotinib respectively, thus identifying responders to both cytotoxic and targeted agents. Many previous clinical trials and in vitro assays have attempted to discover biomarkers of drug sensitivity, but found that the genes/aberrations which they had identified, performed poorly as predictors, once applied to out-of-batch sets of samples. Our models, on the other hand, are trained on an independent set of cancer cell lines and performed well on three completely separate and independent clinical trial datasets (all assessed using different microarray platforms). These results have far-reaching implications for personalized medicine and drug development (e.g. for the development of companion diagnostics). All datasets and bioinformatics tools to reproduce our results are publicly available. Citation Format: Paul Geeleher, Nancy Cox, R. Stephanie Huang. Clinical drug response can be predicted using baseline gene expression levels and in vitro drug sensitivity in cell lines. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 5561. doi:10.1158/1538-7445.AM2014-5561

ASCL1 is a lineage oncogene providing therapeutic targets for high-grade neuroendocrine lung cancers.

Aggressive neuroendocrine lung cancers, including small cell lung cancer (SCLC) and non-small cell lung cancer (NSCLC), represent an understudied tumor subset that accounts for approximately 40,000 new lung cancer cases per year in the United States. No targeted therapy exists for these tumors. We determined that achaete-scute homolog 1 (ASCL1), a transcription factor required for proper development of pulmonary neuroendocrine cells, is essential for the survival of a majority of lung cancers (both SCLC and NSCLC) with neuroendocrine features. By combining whole-genome microarray expression analysis performed on lung cancer cell lines with ChIP-Seq data designed to identify conserved transcriptional targets of ASCL1, we discovered an ASCL1 target 72-gene expression signature that (i) identifies neuroendocrine differentiation in NSCLC cell lines, (ii) is predictive of poor prognosis in resected NSCLC specimens from three datasets, and (iii) represents novel "druggable" targets. Among these druggable targets is B-cell CLL/lymphoma 2, which when pharmacologically inhibited stops ASCL1-dependent tumor growth in vitro and in vivo and represents a proof-of-principle ASCL1 downstream target gene. Analysis of downstream targets of ASCL1 represents an important advance in the development of targeted therapy for the neuroendocrine class of lung cancers, providing a significant step forward in the understanding and therapeutic targeting of the molecular vulnerabilities of neuroendocrine lung cancer.

ATP-binding cassette transporter enhances tolerance to DDT in Tetrahymena.

The reuse of dichlorodiphenyltrichloroethane (DDT) as an indoor residual spray was permitted by the World Health Organization in 2007, and approximately 14 countries still use DDT to control disease vectors. The extensive exposure of insects to DDT has resulted in the emergence of DDT resistance, especially in mosquitoes, and the mechanism for this resistance in mosquitoes has been widely reported. Spraying can also introduce DDT directly into surface water, and DDT can subsequently accumulate in microorganisms, but the mechanism for the resistance to DDT degradation in microorganisms is unclear. Using whole-genome microarray analysis, we detected an abcb15 gene that was up-regulated in a specific manner by DDT treatment in T. thermophile. The deduced ABCB15 peptide sequence had two transmembrane domains (TMDs) and two nucleotide-binding domains (NBDs) to form the structure TMD-NBD-TMD-NBD, and each NBD contained three conserved motifs: Walker-A, C-loop, and Walker-B, which indicated the T. thermophila abcb15 was a typical ABC transporter gene. The expression of ABCB15 fused with a C-terminal green fluorescent protein was found to be on the periphery of the cell, suggesting that ABCB15 was a membrane pump protein. In addition, cells with abcb15 partially knocked down (abcb15-KD) grew slower than wild-type cells in the presence of 256 mg L(-1) DDT, indicating the tolerance of abcb15-KD strain to DDT exposure was decreased. Thus, we suggest that in Tetrahymena, the membrane pump protein encoded by ABCT gene abcb15 can enhance the tolerance to DDT and protect cells from this exogenous toxin by efficiently pumping it to the extracellular space.

Neurotoxic action of microcystin-LR is reflected in the transcriptional stress response of Caenorhabditis elegans

Cyanobacterial blooms in aquatic environments are frequently characterized by elevated levels of microcystins, a potent hepatotoxin. Here we exposed the nematode Caenorhabditis elegans with environmentally realistic concentrations of MC-LR to explore its non-hepatic toxicity. Lifespan, reproduction and growth assays confirmed the toxic potential of 100μg/L MC-LR even in this liver-lacking invertebrate. Whole-genome microarray analysis revealed that a neuromodulating action was the dominant response in nematodes challenged with 100μg/L MC-LR. Indeed, most of the 201 differentially expressed genes were associated with neurobehavior, neurogenesis, and signaling associated pathways. In addition, a whole-genome miRNA-microarray highlighted that, in particular, members of the let-7 family were differentially regulated. These miRNAs are involved in the developmental timing of cell fates, including neurons, and are probably also part of the stress response system. To conclude, neurological modulation is the main transcriptional stress response in C. elegans exposed to MC-LR.

Effect of desipramine on gene expression in the mouse frontal cortex – Microarray study

BackgroundThese studies aimed to identify the genes differentially expressed in the frontal cortex of mice treated repeatedly with either saline or desipramine (DMI). MethodsDifferences in gene expression in the mouse frontal cortex were studied using a whole-genome microarray approach. ResultsThe analyses revealed a group of 88 transcripts (18 genes) that were differentially expressed between the mice treated with saline and those treated with DMI. These genes include Spnb2, Mef2c, Ncam1, Hsp90ab1, Kif1b, Ddx6 and Gsk3b, which were connected in the gene relationship network. ConclusionsIt appears that one week of DMI administration measurably altered the expression of a small number of genes, including genes connected with neuroplasticity and cytoskeletal changes, the regulation of calcium levels in the cell or translation processes.

Iron- and hemin-dependent gene expression of Porphyromonas gingivalis.

Although iron under anaerobic conditions is more accessible and highly reactive because of its reduced form, iron-dependent regulation is not well known in anaerobic bacteria. Here, we investigated iron- and hemin-dependent gene regulation in Porphyromonas gingivalis, an established periodontopathogen that primarily inhabits anaerobic pockets. Whole-genome microarrays of P.gingivalis genes were used to compare the levels of gene expression under iron-replete and iron-depleted conditions as well as under hemin-replete and hemin-depleted conditions. Under iron-depleted conditions, the expression of genes encoding proteins that participate in iron uptake and adhesion/invasion of host cells was increased, while that of genes encoding proteins involved in iron storage, energy metabolism, and electron transport was decreased. Interestingly, many of the genes with altered expression had no known function. Limiting the amount of hemin also resulted in a reduced expression of the genes encoding proteins involved in energy metabolism and electron transport. However, hemin also had a significant effect on many other biological processes such as oxidative stress protection and lipopolysaccharide synthesis. Overall, comparison of the data from iron-depleted conditions to those from hemin-depleted ones showed that although some regulation is through the iron derived from hemin, there also is significant distinct regulation through hemin only. Furthermore, our data showed that the molecular mechanisms of iron-dependent regulation are novel as the deletion of the putative Fur protein had no effect on the expression of iron-regulated genes. Finally, our functional studies demonstrated greater survivability of host cells in the presence of the iron-stressed bacterium than the iron-replete P.gingivalis cells. The major iron-regulated proteins encoded by PG1019-20 may play a role in this process as deletion of these sequences also resulted in reduced survival of the bacterium when grown with eukaryotic cells. Taken together, the results of this study demonstrated the utility of whole-genome microarray analysis for the identification of genes with altered expression profiles during varying growth conditions and provided a framework for the detailed analysis of the molecular mechanisms of iron and hemin acquisition, metabolism and virulence of P.gingivalis.

Toxicogenomic effects of nano- and bulk-TiO2 particles in the soil nematode Caenorhabditis elegans

The toxicity and toxicogenomics of selected anatase and rutile nanoparticles (NP) and bulk titanium dioxide (TiO2) particles were evaluated in the soil nematode Caenorhabditis elegans. Results indicated that bulk or nano-TiO2 particles were slightly toxic to soil nematode C. elegans, as measured by reproduction EC50 values ranging from 4 to 32 mg/L. Whole-genome microarray results indicated that the regulation of glutathione-S-transferase gst-3, cytochrome P450 cypp33-c11, stress resistance regulator scl-1, oxidoreductase wah-1 and embryonic development pod-2 genes were significantly affected by nano-sized and bulk-TiO2 particles. More specifically, it was determined that anatase particles exerted a greater effect on metabolic pathways, whereas rutile particles had a greater effect on developmental processes. The up-regulation of the pod-2 gene corroborated the phenotypic effect observed in the reproduction test. Our results demonstrated that C. elegans is a good genomic model for nano-TiO2 toxicity assessment.

Genome-wide profiling of changes in gene expression in response to infection of the japonica rice variety Yunyin by Magnaporthe oryzae

Rice blast, caused by the fungus Magnaporthe oryzae, is one of the diseases most responsible for significantly decreasing the yield and quality of commercially grown rice. This study investigated the blast resistance response of japonica rice variety Yunyin (YY) treated with the Sichuan-43 isolate of M. oryzae. Seedlings of YY exhibited stronger resistance to blast than those of Lijiangxintuanheigu (LTH). The gene expression profile of YY treated with Sichuan-43 was determined using whole-genome microarray technology. Bioinformatics was used to identify putative resistance-related genes from the large number of genes assayed. Five candidate genes were further characterized by gene ontology classification analysis and pathway enrichment analysis, and were then integrated into various types of gene network regulation diagrams. Systematic bioinformatic analysis of the microarray provided a transcriptome map of YY, which helped to elucidate the mechanisms involved in blast resistance. Our results enhance the current understanding of the effects of rice blast resistance genes at the transcription level, and will facilitate further exploration of the molecular mechanism of blast resistance in YY.

Gene expression profiling of three tissues of chicken after heat stress treatment by microarray technique

In the present study, tissue samples were collected from the cerebrum, liver, and leg muscle of 8-day-old dwarf chicks that were exposed to a 3 h treatment of 28℃± 1℃ (control group) or 40℃± 1℃ (treatment group). Differentially expressed (DE) genes in these samples were detected using whole-genome microarray chips, and their functions were.

The Sex-Determining Factors SRY and SOX9 Regulate Similar Target Genes and Promote Testis Cord Formation during Testicular Differentiation

Male sex determination is mediated sequentially by sex-determining region Y (SRY) and related SRY-box 9 (SOX9) transcription factors. To understand the gene regulatory hierarchy for SRY and SOX9, a series of chromatin immunoprecipitation and whole-genome promoter tiling microarray (ChIP-Chip) experiments were conducted with mouse gonadal cells at the time of sex determination. SRY and SOX9 bind to the promoters of many common targets involved in testis differentiation and regulate their expression in Sertoli cells. SRY binds to various ovarian differentiation genes and represses their activation through WNT/β-catenin signaling. Sertoli cell-Sertoli cell junction signaling, important for testis cord formation, is the top canonical pathway among the SRY and SOX9 targets. Hence, SRY determines Sertoli cell fate by repressing ovarian and activating testicular differentiation genes, promotes early Sertoli cells to form testis cord, and then passes on its functions to SOX9, which regulates common targets and activates its own gene regulatory program, beyond SRY actions, in sex determination.