Insights into the evolution of Darwin's finches from comparative analysis of the Geospiza magnirostris genome sequence.

BACKGROUND: A classical example of repeated speciation coupled with ecological diversification is the evolution of 14 closely related species of Darwin's (Galápagos) finches (Thraupidae, Passeriformes). Their adaptive radiation in the Galápagos archipelago took place in the last 2-3 million years and some of the molecular mechanisms that led to their diversification are now being elucidated. Here we report evolutionary analyses of genome of the large ground finch, Geospiza magnirostris. RESULTS: 13,291 protein-coding genes were predicted from a 991.0 Mb G. magnirostris genome assembly. We then defined gene orthology relationships and constructed whole genome alignments between the G. magnirostris and other vertebrate genomes. We estimate that 15% of genomic sequence is functionally constrained between G. magnirostris and zebra finch. Genic evolutionary rate comparisons indicate that similar selective pressures acted along the G. magnirostris and zebra finch lineages suggesting that historical effective population size values have been similar in both lineages. 21 otherwise highly conserved genes were identified that each show evidence for positive selection on amino acid changes in the Darwin's finch lineage. Two of these genes (Igf2r and Pou1f1) have been implicated in beak morphology changes in Darwin's finches. Five of 47 genes showing evidence of positive selection in early passerine evolution have cilia related functions, and may be examples of adaptively evolving reproductive proteins. CONCLUSIONS: These results provide insights into past evolutionary processes that have shaped G. magnirostris genes and its genome, and provide the necessary foundation upon which to build population genomics resources that will shed light on more contemporaneous adaptive and non-adaptive processes that have contributed to the evolution of the Darwin's finches.

Progressive genome-wide introgression in agricultural Campylobacter coli.

Hybridization between distantly related organisms can facilitate rapid adaptation to novel environments, but is potentially constrained by epistatic fitness interactions among cell components. The zoonotic pathogens Campylobacter coli and C. jejuni differ from each other by around 15% at the nucleotide level, corresponding to an average of nearly 40 amino acids per protein-coding gene. Using whole genome sequencing, we show that a single C. coli lineage, which has successfully colonized an agricultural niche, has been progressively accumulating C. jejuni DNA. Members of this lineage belong to two groups, the ST-828 and ST-1150 clonal complexes. The ST-1150 complex is less frequently isolated and has undergone a substantially greater amount of introgression leading to replacement of up to 23% of the C. coli core genome as well as import of novel DNA. By contrast, the more commonly isolated ST-828 complex bacteria have 10-11% introgressed DNA, and C. jejuni and nonagricultural C. coli lineages each have <2%. Thus, the C. coli that colonize agriculture, and consequently cause most human disease, have hybrid origin, but this cross-species exchange has so far not had a substantial impact on the gene pools of either C. jejuni or nonagricultural C. coli. These findings also indicate remarkable interchangeability of basic cellular machinery after a prolonged period of independent evolution.

Genome sequencing and mapping reveal loss of heterozygosity as a mechanism for rapid adaptation in the vegetable pathogen Phytophthora capsici.

The oomycete vegetable pathogen Phytophthora capsici has shown remarkable adaptation to fungicides and new hosts. Like other members of this destructive genus, P. capsici has an explosive epidemiology, rapidly producing massive numbers of asexual spores on infected hosts. In addition, P. capsici can remain dormant for years as sexually recombined oospores, making it difficult to produce crops at infested sites, and allowing outcrossing populations to maintain significant genetic variation. Genome sequencing, development of a high-density genetic map, and integrative genomic or genetic characterization of P. capsici field isolates and intercross progeny revealed significant mitotic loss of heterozygosity (LOH) in diverse isolates. LOH was detected in clonally propagated field isolates and sexual progeny, cumulatively affecting >30% of the genome. LOH altered genotypes for more than 11,000 single-nucleotide variant sites and showed a strong association with changes in mating type and pathogenicity. Overall, it appears that LOH may provide a rapid mechanism for fixing alleles and may be an important component of adaptability for P. capsici.

The genome sequence of the wine yeast VIN7 reveals an allotriploid hybrid genome with Saccharomyces cerevisiae and Saccharomyces kudriavzevii origins.

The vast majority of wine fermentations are performed principally by Saccharomyces cerevisiae. However, there are a growing number of instances in which other species of Saccharomyces play a predominant role. Interestingly, the presence of these other yeast species generally occurs via the formation of interspecific hybrids that contain genomic contributions from both S. cerevisiae and non-S. cerevisiae species. However, despite the large number of wine strains that are characterized at the genomic level, there remains limited information regarding the detailed genomic structure of hybrids used in winemaking. To address this, we describe the genome sequence of the thiol-releasing commercial wine yeast hybrid VIN7. VIN7 is shown to be an almost complete allotriploid interspecific hybrid that is comprised of a heterozygous diploid complement of S. cerevisiae chromosomes and a haploid Saccharomyces kudriavzevii genomic contribution. Both parental strains appear to be of European origin, with the S. cerevisiae parent being closely related to, but distinct from, the commercial wine yeasts QA23 and EC1118. In addition, several instances of chromosomal rearrangement between S. cerevisiae and S. kudriavzevii sequences were observed that may mark the early stages of hybrid genome consolidation.

Comparative and functional genomics of Rhodococcus opacus PD630 for biofuels development.

The Actinomycetales bacteria Rhodococcus opacus PD630 and Rhodococcus jostii RHA1 bioconvert a diverse range of organic substrates through lipid biosynthesis into large quantities of energy-rich triacylglycerols (TAGs). To describe the genetic basis of the Rhodococcus oleaginous metabolism, we sequenced and performed comparative analysis of the 9.27 Mb R. opacus PD630 genome. Metabolic-reconstruction assigned 2017 enzymatic reactions to the 8632 R. opacus PD630 genes we identified. Of these, 261 genes were implicated in the R. opacus PD630 TAGs cycle by metabolic reconstruction and gene family analysis. Rhodococcus synthesizes uncommon straight-chain odd-carbon fatty acids in high abundance and stores them as TAGs. We have identified these to be pentadecanoic, heptadecanoic, and cis-heptadecenoic acids. To identify bioconversion pathways, we screened R. opacus PD630, R. jostii RHA1, Ralstonia eutropha H16, and C. glutamicum 13032 for growth on 190 compounds. The results of the catabolic screen, phylogenetic analysis of the TAGs cycle enzymes, and metabolic product characterizations were integrated into a working model of prokaryotic oleaginy.

Whole-genome comparison reveals novel genetic elements that characterize the genome of industrial strains of Saccharomyces cerevisiae.

Human intervention has subjected the yeast Saccharomyces cerevisiae to multiple rounds of independent domestication and thousands of generations of artificial selection. As a result, this species comprises a genetically diverse collection of natural isolates as well as domesticated strains that are used in specific industrial applications. However the scope of genetic diversity that was captured during the domesticated evolution of the industrial representatives of this important organism remains to be determined. To begin to address this, we have produced whole-genome assemblies of six commercial strains of S. cerevisiae (four wine and two brewing strains). These represent the first genome assemblies produced from S. cerevisiae strains in their industrially-used forms and the first high-quality assemblies for S. cerevisiae strains used in brewing. By comparing these sequences to six existing high-coverage S. cerevisiae genome assemblies, clear signatures were found that defined each industrial class of yeast. This genetic variation was comprised of both single nucleotide polymorphisms and large-scale insertions and deletions, with the latter often being associated with ORF heterogeneity between strains. This included the discovery of more than twenty probable genes that had not been identified previously in the S. cerevisiae genome. Comparison of this large number of S. cerevisiae strains also enabled the characterization of a cluster of five ORFs that have integrated into the genomes of the wine and bioethanol strains on multiple occasions and at diverse genomic locations via what appears to involve the resolution of a circular DNA intermediate. This work suggests that, despite the scrutiny that has been directed at the yeast genome, there remains a significant reservoir of ORFs and novel modes of genetic transmission that may have significant phenotypic impact in this important model and industrial species.

The genome of the domesticated apple (Malus × domestica Borkh.).

We report a high-quality draft genome sequence of the domesticated apple (Malus × domestica). We show that a relatively recent (>50 million years ago) genome-wide duplication (GWD) has resulted in the transition from nine ancestral chromosomes to 17 chromosomes in the Pyreae. Traces of older GWDs partly support the monophyly of the ancestral paleohexaploidy of eudicots. Phylogenetic reconstruction of Pyreae and the genus Malus, relative to major Rosaceae taxa, identified the progenitor of the cultivated apple as M. sieversii. Expansion of gene families reported to be involved in fruit development may explain formation of the pome, a Pyreae-specific false fruit that develops by proliferation of the basal part of the sepals, the receptacle. In apple, a subclade of MADS-box genes, normally involved in flower and fruit development, is expanded to include 15 members, as are other gene families involved in Rosaceae-specific metabolism, such as transport and assimilation of sorbitol.

A draft sequence of the Neandertal genome.

Neandertals, the closest evolutionary relatives of present-day humans, lived in large parts of Europe and western Asia before disappearing 30,000 years ago. We present a draft sequence of the Neandertal genome composed of more than 4 billion nucleotides from three individuals. Comparisons of the Neandertal genome to the genomes of five present-day humans from different parts of the world identify a number of genomic regions that may have been affected by positive selection in ancestral modern humans, including genes involved in metabolism and in cognitive and skeletal development. We show that Neandertals shared more genetic variants with present-day humans in Eurasia than with present-day humans in sub-Saharan Africa, suggesting that gene flow from Neandertals into the ancestors of non-Africans occurred before the divergence of Eurasian groups from each other.

Organismal, genetic, and transcriptional variation in the deeply sequenced gut microbiomes of identical twins.

We deeply sampled the organismal, genetic, and transcriptional diversity in fecal samples collected from a monozygotic (MZ) twin pair and compared the results to 1,095 communities from the gut and other body habitats of related and unrelated individuals. Using a new scheme for noise reduction in pyrosequencing data, we estimated the total diversity of species-level bacterial phylotypes in the 1.2-1.5 million bacterial 16S rRNA reads obtained from each deeply sampled cotwin to be approximately 800 (35.9%, 49.1% detected in both). A combined 1.1 million read 16S rRNA dataset representing 281 shallowly sequenced fecal samples from 54 twin pairs and their mothers contained an estimated 4,018 species-level phylotypes, with each sample having a unique species assemblage (53.4 +/- 0.6% and 50.3 +/- 0.5% overlap with the deeply sampled cotwins). Of the 134 phylotypes with a relative abundance of >0.1% in the combined dataset, only 37 appeared in >50% of the samples, with one phylotype in the Lachnospiraceae family present in 99%. Nongut communities had significantly reduced overlap with the deeply sequenced twins' fecal microbiota (18.3 +/- 0.3%, 15.3 +/- 0.3%). The MZ cotwins' fecal DNA was deeply sequenced (3.8-6.3 Gbp/sample) and assembled reads were assigned to 25 genus-level phylogenetic bins. Only 17% of the genes in these bins were shared between the cotwins. Bins exhibited differences in their degree of sequence variation, gene content including the repertoire of carbohydrate active enzymes present within and between twins (e.g., predicted cellulases, dockerins), and transcriptional activities. These results provide an expanded perspective about features that make each of us unique life forms and directions for future characterization of our gut ecosystems.

Metabolic streamlining in an open-ocean nitrogen-fixing cyanobacterium.

Nitrogen (N(2))-fixing marine cyanobacteria are an important source of fixed inorganic nitrogen that supports oceanic primary productivity and carbon dioxide removal from the atmosphere. A globally distributed, periodically abundant N(2)-fixing marine cyanobacterium, UCYN-A, was recently found to lack the oxygen-producing photosystem II complex of the photosynthetic apparatus, indicating a novel metabolism, but remains uncultivated. Here we show, from metabolic reconstructions inferred from the assembly of the complete UCYN-A genome using massively parallel pyrosequencing of paired-end reads, that UCYN-A has a photofermentative metabolism and is dependent on other organisms for essential compounds. We found that UCYN-A lacks a number of major metabolic pathways including the tricarboxylic acid cycle, but retains sufficient electron transport capacity to generate energy and reducing power from light. Unexpectedly, UCYN-A has a reduced genome (1.44 megabases) that is structurally similar to many chloroplasts and some bacteria, in that it contains inverted repeats of ribosomal RNA operons. The lack of biosynthetic pathways for several amino acids and purines suggests that this organism depends on other organisms, either in close association or in symbiosis, for critical nutrients. However, size fractionation experiments using natural populations have so far not provided evidence of a symbiotic association with another microorganism. The UCYN-A cyanobacterium is a paradox in evolution and adaptation to the marine environment, and is an example of the tight metabolic coupling between microorganisms in oligotrophic oceanic microbial communities.

Cdcs1 a major colitis susceptibility locus in mice; subcongenic analysis reveals genetic complexity.

BACKGROUND: The cytokine-deficiency-induced colitis susceptibility (Cdcs)1 locus is a major modifier of murine inflammatory bowel disease (IBD) and was originally identified in experimental crosses of interleukin-10-deficient (Il10(-/-)) mice. Congenic mice, in which this locus was reciprocally transferred between IBD-susceptible C3H/HeJBir-Il10(-/-) and resistant C57BL/6J-Il10(-/-) mice, revealed that this locus likely acts by inducing innate hypo- and adaptive hyperresponsiveness, associated with impaired NF-kappaB responses of macrophages. The aim of the present study was to dissect the complexity of Cdcs1 by further development and characterization of reciprocal Cdcs1 congenic strains and to identify potential candidate genes in the congenic interval. METHODS: In total, 15 reciprocal congenic strains were generated from Il10(-/-) mice of either C3H/HeJBir or C57BL/6J genetic backgrounds by 10 cycles of backcrossing. Colitis activity was monitored by histological grading. Candidate genes were identified by fine mapping of congenic intervals, sequencing, microarray analysis, and a high-throughput real-time reverse-transcription polymerase chain reaction (RT-PCR) approach using bone marrow-derived macrophages. RESULTS: Within the originally identified Cdcs1-interval, 3 independent regions were detected that likely contain susceptibility-determining genetic factors (Cdcs1.1, Cdcs1.2, and Cdcs1.3). Combining results of candidate gene approaches revealed Fcgr1, Cnn3, Larp7, and Alpk1 as highly attractive candidate genes with polymorphisms in coding or regulatory regions and expression differences between susceptible and resistant mouse strains. CONCLUSIONS: Subcongenic analysis of the major susceptibility locus Cdcs1 on mouse chromosome 3 revealed a complex genetic structure. Candidate gene approaches revealed attractive genes within the identified regions.

Mutation discovery in the mouse using genetically guided array capture and resequencing.

Forward genetics (phenotype-driven approaches) remain the primary source for allelic variants in the mouse. Unfortunately, the gap between observable phenotype and causative genotype limits the widespread use of spontaneous and induced mouse mutants. As alternatives to traditional positional cloning and mutation detection approaches, sequence capture and next-generation sequencing technologies can be used to rapidly sequence subsets of the genome. Application of these technologies to mutation detection efforts in the mouse has the potential to significantly reduce the time and resources required for mutation identification by abrogating the need for high-resolution genetic mapping, long-range PCR, and sequencing of individual PCR amplimers. As proof of principle, we used array-based sequence capture and pyrosequencing to sequence an allelic series from the classically defined Kit locus (approximately 200 kb) from each of five noncomplementing Kit mutants (one known allele and four unknown alleles) and have successfully identified and validated a nonsynonymous coding mutation for each allele. These data represent the first documentation and validation that these new technologies can be used to efficiently discover causative mutations. Importantly, these data also provide a specific methodological foundation for the development of large-scale mutation detection efforts in the laboratory mouse.

A core gut microbiome in obese and lean twins.

The human distal gut harbours a vast ensemble of microbes (the microbiota) that provide important metabolic capabilities, including the ability to extract energy from otherwise indigestible dietary polysaccharides. Studies of a few unrelated, healthy adults have revealed substantial diversity in their gut communities, as measured by sequencing 16S rRNA genes, yet how this diversity relates to function and to the rest of the genes in the collective genomes of the microbiota (the gut microbiome) remains obscure. Studies of lean and obese mice suggest that the gut microbiota affects energy balance by influencing the efficiency of calorie harvest from the diet, and how this harvested energy is used and stored. Here we characterize the faecal microbial communities of adult female monozygotic and dizygotic twin pairs concordant for leanness or obesity, and their mothers, to address how host genotype, environmental exposure and host adiposity influence the gut microbiome. Analysis of 154 individuals yielded 9,920 near full-length and 1,937,461 partial bacterial 16S rRNA sequences, plus 2.14 gigabases from their microbiomes. The results reveal that the human gut microbiome is shared among family members, but that each person's gut microbial community varies in the specific bacterial lineages present, with a comparable degree of co-variation between adult monozygotic and dizygotic twin pairs. However, there was a wide array of shared microbial genes among sampled individuals, comprising an extensive, identifiable 'core microbiome' at the gene, rather than at the organismal lineage, level. Obesity is associated with phylum-level changes in the microbiota, reduced bacterial diversity and altered representation of bacterial genes and metabolic pathways. These results demonstrate that a diversity of organismal assemblages can nonetheless yield a core microbiome at a functional level, and that deviations from this core are associated with different physiological states (obese compared with lean).

Globally distributed uncultivated oceanic N2-fixing cyanobacteria lack oxygenic photosystem II.

Biological nitrogen (N2) fixation is important in controlling biological productivity and carbon flux in the oceans. Unicellular N2-fixing cyanobacteria have only recently been discovered and are widely distributed in tropical and subtropical seas. Metagenomic analysis of flow cytometry-sorted cells shows that unicellular N2-fixing cyanobacteria in "group A" (UCYN-A) lack genes for the oxygen-evolving photosystem II and for carbon fixation, which has implications for oceanic carbon and nitrogen cycling and raises questions regarding the evolution of photosynthesis and N2 fixation on Earth.

Analysis of copy number variants and segmental duplications in the human genome: Evidence for a change in the process of formation in recent evolutionary history.

Segmental duplications (SDs) are operationally defined as >1 kb stretches of duplicated DNA with high sequence identity. They arise from copy number variants (CNVs) fixed in the population. To investigate the formation of SDs and CNVs, we examine their large-scale patterns of co-occurrence with different repeats. Alu elements, a major class of genomic repeats, had previously been identified as prime drivers of SD formation. We also observe this association; however, we find that it sharply decreases for younger SDs. Continuing this trend, we find only weak associations of CNVs with Alus. Similarly, we find an association of SDs with processed pseudogenes, which is decreasing for younger SDs and absent entirely for CNVs. Next, we find that SDs are significantly co-localized with each other, resulting in a highly skewed "power-law" distribution and chromosomal hotspots. We also observe a significant association of CNVs with SDs, but find that an SD-mediated mechanism only accounts for some CNVs (<28%). Overall, our results imply that a shift in predominant formation mechanism occurred in recent history: approximately 40 million years ago, during the "Alu burst" in retrotransposition activity, non-allelic homologous recombination, first mediated by Alus and then the by newly formed CNVs themselves, was the main driver of genome rearrangements; however, its relative importance has decreased markedly since then, with proportionally more events now stemming from other repeats and from non-homologous end-joining. In addition to a coarse-grained analysis, we performed targeted sequencing of 67 CNVs and then analyzed a combined set of 270 CNVs (540 breakpoints) to verify our conclusions.

Applying gene expression, proteomics and single-nucleotide polymorphism analysis for complex trait gene identification.

Previous quantitative trait locus (QTL) analysis of an intercross involving the inbred mouse strains NZB/BlNJ and SM/J revealed QTL for a variety of complex traits. Many QTL have large intervals containing hundreds of genes, and methods are needed to rapidly sort through these genes for probable candidates. We chose nine QTL: the three most significant for high-density lipoprotein (HDL) cholesterol, gallstone formation, and obesity. We searched for candidate genes using three different approaches: mRNA microarray gene expression technology to assess >45,000 transcripts, publicly available SNPs to locate genes that are not identical by descent and that contain nonsynonymous coding differences, and a mass-spectrometry-based proteomics technology to interrogate nearly 1000 proteins for differential expression in the liver of the two parental inbred strains. This systematic approach reduced the number of candidate genes within each QTL from hundreds to a manageable list. Each of the three approaches selected candidates that the other two approaches missed. For example, candidate genes such as Apoa2 and Acads had differential protein levels although the mRNA levels were similar. We conclude that all three approaches are important and that focusing on a single approach such as mRNA expression may fail to identify a QTL gene.

Oocyte regulation of metabolic cooperativity between mouse cumulus cells and oocytes: BMP15 and GDF9 control cholesterol biosynthesis in cumulus cells.

Oocyte-derived bone morphogenetic protein 15 (BMP15) and growth differentiation factor 9 (GDF9) are key regulators of follicular development. Here we show that these factors control cumulus cell metabolism, particularly glycolysis and cholesterol biosynthesis before the preovulatory surge of luteinizing hormone. Transcripts encoding enzymes for cholesterol biosynthesis were downregulated in both Bmp15(-/-) and Bmp15(-/-) Gdf9(+/-) double mutant cumulus cells, and in wild-type cumulus cells after removal of oocytes from cumulus-cell-oocyte complexes. Similarly, cholesterol synthesized de novo was reduced in these cumulus cells. This indicates that oocytes regulate cumulus cell cholesterol biosynthesis by promoting the expression of relevant transcripts. Furthermore, in wild-type mice, Mvk, Pmvk, Fdps, Sqle, Cyp51, Sc4mol and Ebp, which encode enzymes required for cholesterol synthesis, were highly expressed in cumulus cells compared with oocytes; and oocytes, in the absence of the surrounding cumulus cells, synthesized barely detectable levels of cholesterol. Furthermore, coincident with reduced cholesterol synthesis in double mutant cumulus cells, lower levels were also detected in cumulus-cell-enclosed double mutant oocytes compared with wild-type oocytes. Levels of cholesterol synthesis in double mutant cumulus cells and oocytes were partially restored by co-culturing with wild-type oocytes. Together, these results indicate that mouse oocytes are deficient in synthesizing cholesterol and require cumulus cells to provide products of the cholesterol biosynthetic pathway. Therefore, oocyte-derived paracrine factors, particularly, BMP15 and GDF9, promote cholesterol biosynthesis in cumulus cells, probably as compensation for oocyte deficiencies in cholesterol production.

A mouse model of conditional lipodystrophy.

Lipodystrophies are syndromes of adipose tissue degeneration associated with severe defects in lipid and glucose homeostasis. We report here the generation and analysis of Pparg(ldi), a targeted allele that confers conditional dominant lipodystrophy in mice. The Pparg(ldi) allele was generated by insertion of the Tet activator (tTA) and a tTA-regulated Flag-Pparg1 transgene into the Pparg gene. Unexpectedly, tTA elicits mild lipodystrophy, insulin resistance, and dyslipidemia, and the Flag-PPARgamma1 transgene surprisingly exacerbates these traits. Doxycycline can both completely prevent and reverse these phenotypes, providing a mouse model of inducible lipodystrophy. Embryonic fibroblasts from either Pparg(ldi/+) or the phenotypically similar aP2-nSrebp1c (Sr) transgenic mice undergo robust adipogenesis, suggesting that neither strain develops lipodystrophy because of defective adipocyte differentiation. In addition, Pparg(ldi/+) adipose tissue shares extensive gene expression aberrations with that of Sr mice, authenticating the phenotype at the molecular level and revealing a common expression signature of lipodystrophic fat. Thus, the Pparg(ldi/+) mouse provides a conditional animal model for studying lipodystrophy and its associated physiology and gene expression.

Paired-end mapping reveals extensive structural variation in the human genome.

Structural variation of the genome involves kilobase- to megabase-sized deletions, duplications, insertions, inversions, and complex combinations of rearrangements. We introduce high-throughput and massive paired-end mapping (PEM), a large-scale genome-sequencing method to identify structural variants (SVs) approximately 3 kilobases (kb) or larger that combines the rescue and capture of paired ends of 3-kb fragments, massive 454 sequencing, and a computational approach to map DNA reads onto a reference genome. PEM was used to map SVs in an African and in a putatively European individual and identified shared and divergent SVs relative to the reference genome. Overall, we fine-mapped more than 1000 SVs and documented that the number of SVs among humans is much larger than initially hypothesized; many of the SVs potentially affect gene function. The breakpoint junction sequences of more than 200 SVs were determined with a novel pooling strategy and computational analysis. Our analysis provided insights into the mechanisms of SV formation in humans.

Transcriptional profile of mouse pre-granulosa and Sertoli cells isolated from early-differentiated fetal gonads.

Gonadal sex determining (GSD) genes that initiate fetal ovarian and testicular development and differentiation are expressed in the cells of the urogenital ridge that differentiate as somatic support cells (SSCs), i.e., granulosa cells of the ovary and Sertoli cells of the testis. To identify potential new mammalian GSD genes, we used the Mouse Genome 430v2.0 GeneChip to analyze gene expression differences between XX and XY SSCs cells isolated from the gonads of embryonic day (E) 13 C57BL/6J fetuses carrying an EGFP reporter transgene expressed specifically in SSCs. In addition, genome wide expression differences between XX and XY E13 whole gonads were examined. Analysis revealed that XX and XY E13 SSCs differentially express 647 transcripts (False Discovery Rate cutoff 1%), including transcripts not previously reported to exhibit a sexually dimorphic expression pattern in this unique cell population. Enrichment for genes controlling cell proliferation was noted in XY SSCs, whereas enrichment for genes controlling cell morphology and metabolic status was identified in XX SSCs. Among the newly identified differentially expressed transcripts are potential GSD genes involved in unexplained human sex reversal cases.