Population-size history inferences from the coho salmon (Oncorhynchus kisutch) genome.

Coho salmon (Oncorhynchus kisutch) are a culturally and economically important species that return from multiyear ocean migrations to spawn in rivers that flow to the Northern Pacific Ocean. Southern stocks of coho salmon in Canada and the United States have significantly declined over the past quarter century, and unfortunately, conservation efforts have not reversed this trend. To assist in stock management and conservation efforts, we generated a chromosome-level genome assembly. We also resequenced the genomes of 83 coho salmon across the North American range to identify nucleotide variants and understand the demographic histories of these salmon by modeling effective population size from genome-wide data. From demographic history modeling, we observed reductions in effective population sizes between 3,750 and 8,000 years ago for several northern sampling sites, which may correspond to bottleneck events during recolonization after glacial retreat.

Pulmonary surfactant protein B carried by HDL predicts incident CVD in patients with type 1 diabetes.

Atherosclerotic CVD is the major cause of death in patients with type 1 diabetes mellitus (T1DM). Alterations in the HDL proteome have been shown to associate with prevalent CVD in T1DM. We therefore sought to determine which proteins carried by HDL might predict incident CVD in patients with T1DM. Using targeted MS/MS, we quantified 50 proteins in HDL from 181 T1DM subjects enrolled in the prospective Coronary Artery Calcification in Type 1 Diabetes study. We used Cox proportional regression analysis and a case-cohort design to test associations of HDL proteins with incident CVD (myocardial infarction, coronary artery bypass grafting, angioplasty, or death from coronary heart disease). We found that only one HDL protein-SFTPB (pulmonary surfactant protein B)-predicted incident CVD in all the models tested. In a fully adjusted model that controlled for lipids and other risk factors, the hazard ratio was 2.17 per SD increase of SFTPB (95% confidence interval, 1.12-4.21, P = 0.022). In addition, plasma fractionation demonstrated that SFTPB is nearly entirely bound to HDL. Although previous studies have shown that high plasma levels of SFTPB associate with prevalent atherosclerosis only in smokers, we found that SFTPB predicted incident CVD in T1DM independently of smoking status and a wide range of confounding factors, including HDL-C, LDL-C, and triglyceride levels. Because SFTPB is almost entirely bound to plasma HDL, our observations support the proposal that SFTPB carried by HDL is a marker-and perhaps mediator-of CVD risk in patients with T1DM.

Whole Genome Linkage Disequilibrium and Effective Population Size in a Coho Salmon (Oncorhynchus kisutch) Breeding Population Using a High-Density SNP Array.

The estimation of linkage disequilibrium between molecular markers within a population is critical when establishing the minimum number of markers required for association studies, genomic selection, and inferring historical events influencing different populations. This work aimed to evaluate the extent and decay of linkage disequilibrium in a coho salmon breeding population using a high-density SNP array. Linkage disequilibrium was estimated between a total of 93,502 SNPs found in 64 individuals (33 dams and 31 sires) from the breeding population. The markers encompass all 30 coho salmon chromosomes and comprise 1,684.62 Mb of the genome. The average density of markers per chromosome ranged from 48.31 to 66 per 1 Mb. The minor allele frequency averaged 0.26 (with a range from 0.22 to 0.27). The overall average linkage disequilibrium among SNPs pairs measured as r 2 was 0.10. The Average r 2 value decreased with increasing physical distance, with values ranging from 0.21 to 0.07 at a distance lower than 1 kb and up to 10 Mb, respectively. An r 2 threshold of 0.2 was reached at distance of approximately 40 Kb. Chromosomes Okis05, Okis15 and Okis28 showed high levels of linkage disequilibrium (>0.20 at distances lower than 1 Mb). Average r 2 values were lower than 0.15 for all chromosomes at distances greater than 4 Mb. An effective population size of 43 was estimated for the population 10 generations ago, and 325, for 139 generations ago. Based on the effective number of chromosome segments, we suggest that at least 74,000 SNPs would be necessary for an association mapping study and genomic predictions. Therefore, the SNP panel used allowed us to capture high-resolution information in the farmed coho salmon population. Furthermore, based on the contemporary N e, a new mate allocation strategy is suggested to increase the effective population size.

Design and characterization of an 87k SNP genotyping array for Arctic charr (Salvelinus alpinus).

We have generated a high-density, high-throughput genotyping array for characterizing genome-wide variation in Arctic charr (Salvelinus alpinus). Novel single nucleotide polymorphisms (SNPs) were identified in charr from the Fraser, Nauyuk and Tree River aquaculture strains, which originated from northern Canada and fish from Iceland using high coverage sequencing, reduced representation sequencing and RNA-seq datasets. The array was designed to capture genome-wide variation from a diverse suite of Arctic charr populations. Cross validation of SNPs from various sources and comparison with previously published Arctic charr SNP data provided a set of candidate SNPs that generalize across populations. Further candidate SNPs were identified based on minor allele frequency, association with RNA transcripts, even spacing across intergenic regions and association with the sex determining (sdY) gene. The performance of the 86,503 SNP array was assessed by genotyping Fraser, Nauyuk and Tree River strain individuals, as well as wild Icelandic Arctic charr. Overall, 63,060 of the SNPs were polymorphic within at least one group and 36.8% were unique to one of the four groups, suggesting that the array design allows for characterization of both within and across population genetic diversity. The concordance between sdY markers and known phenotypic sex indicated that the array can accurately determine the sex of individuals based on genotype alone. The Salp87k genotyping array provides researchers and breeders the opportunity to analyze genetic variation in Arctic charr at a more detailed level than previously possible.

The Arctic charr (Salvelinus alpinus) genome and transcriptome assembly.

Arctic charr have a circumpolar distribution, persevere under extreme environmental conditions, and reach ages unknown to most other salmonids. The Salvelinus genus is primarily composed of species with genomes that are structured more like the ancestral salmonid genome than most Oncorhynchus and Salmo species of sister genera. It is thought that this aspect of the genome may be important for local adaptation (due to increased recombination) and anadromy (the migration of fish from saltwater to freshwater). In this study, we describe the generation of a new genetic map, the sequencing and assembly of the Arctic charr genome (GenBank accession: GCF_002910315.2) using the newly created genetic map and a previous genetic map, and present several analyses of the Arctic charr genes and genome assembly. The newly generated genetic map consists of 8,574 unique genetic markers and is similar to previous genetic maps with the exception of three major structural differences. The N50, identified BUSCOs, repetitive DNA content, and total size of the Arctic charr assembled genome are all comparable to other assembled salmonid genomes. An analysis to identify orthologous genes revealed that a large number of orthologs could be identified between salmonids and many appear to have highly conserved gene expression profiles between species. Comparing orthologous gene expression profiles may give us a better insight into which genes are more likely to influence species specific phenotypes.

Regulatory processes that control haploid expression of salmon sperm mRNAs.

OBJECTIVE: Various stages of mRNA processing are necessary for functionally important genes required during late-stage sperm differentiation. Protein-RNA complexes form that edit, stabilize, store, deliver, localize and regulate translation of sperm mRNAs. These regulatory processes are often directed by recognition sequence elements and the particular composition of the proteins associated with the mRNAs. Previous work has shown that the cAMP response element modulator (CREM), estrogen receptor-alpha (ERα) and forkhead box L2A (FOXL2A) proteins are present in late-stage salmon sperm. Here we investigate whether these and other regulatory proteins might control processing of mRNAs not expressed until the haploid stage of development. We also examine regulatory processes that prepare and present mRNAs that generate unique products essential for differentiating sperm (i.e. for flagellar assembly and function). RESULTS: We provide evidence for potential sperm-specific recognition elements in 5'-untranslated regions (utrs) that may bind CREM, ERα, FOXL2A, Y-box and other proteins. We show that changes within the 5'-utrs and open reading frames of some sperm genes lead to distinct protein termini that may provide specific interfaces necessary for localization and function within the paternal gamete.

Genomic Predictions and Genome-Wide Association Study of Resistance Against Piscirickettsia salmonis in Coho Salmon (Oncorhynchus kisutch) Using ddRAD Sequencing.

Piscirickettsia salmonis is one of the main infectious diseases affecting coho salmon (Oncorhynchus kisutch) farming, and current treatments have been ineffective for the control of this disease. Genetic improvement for P. salmonis resistance has been proposed as a feasible alternative for the control of this infectious disease in farmed fish. Genotyping by sequencing (GBS) strategies allow genotyping of hundreds of individuals with thousands of single nucleotide polymorphisms (SNPs), which can be used to perform genome wide association studies (GWAS) and predict genetic values using genome-wide information. We used double-digest restriction-site associated DNA (ddRAD) sequencing to dissect the genetic architecture of resistance against P. salmonis in a farmed coho salmon population and to identify molecular markers associated with the trait. We also evaluated genomic selection (GS) models in order to determine the potential to accelerate the genetic improvement of this trait by means of using genome-wide molecular information. A total of 764 individuals from 33 full-sib families (17 highly resistant and 16 highly susceptible) were experimentally challenged against P. salmonis and their genotypes were assayed using ddRAD sequencing. A total of 9,389 SNPs markers were identified in the population. These markers were used to test genomic selection models and compare different GWAS methodologies for resistance measured as day of death (DD) and binary survival (BIN). Genomic selection models showed higher accuracies than the traditional pedigree-based best linear unbiased prediction (PBLUP) method, for both DD and BIN. The models showed an improvement of up to 95% and 155% respectively over PBLUP. One SNP related with B-cell development was identified as a potential functional candidate associated with resistance to P. salmonis defined as DD.

Subcellular localization and characterization of estrogenic pathway regulators and mediators in Atlantic salmon spermatozoal cells.

Much progress has been made regarding our understanding of aromatase regulation, estrogen synthesis partitioning and communication between the germinal and somatic compartments of the differentiating gonad. We now know that most of the enzymatic and signaling apparatus required for steroidogenesis is endogenously expressed within germ cells. However, less is known about the expression and localization of steroidogenic components within mature spermatozoa. We have assembled a sperm library presenting 197,015 putative transcripts. Co-expression clustering analysis revealed that 6687 genes were present at higher levels in sperm in comparison to fifteen other salmon tissue libraries. The sperm transcriptome is highly complex containing the highest proportion of unannotated genes (45%) of the tissues analyzed. Our analysis of highly expressed genes in late-stage sperm revealed dedication to tasks involving chromatin remodeling, flagellogenesis and proteolysis. In addition, using various different embedding and microscopic techniques, we examined the morphology of salmon spermatozoa and characterized expression and localization of several estrogenic regulatory and signaling proteins by immunohistochemistry. We provide evidence for the endogenous synthesis and localization of aromatase (CYP19A and CYP19B1) and potential mediators of estrogen [i.e., ER-alpha and soluble adenylyl cyclase (sAC)] or phosphate (i.e., CREB and FOXL2A) signaling. Partitioning of select transcripts that encode AR-beta, FSH and the LH receptor, but not AR-alpha, LH or the FSH receptor, further points to localized specificity of function in the steroidogenic circuitry of the sperm cell. These results open new avenues of investigation to further our understanding of the intra- and intercellular regulatory processes that guide sperm development and biology.

Functional Annotation of All Salmonid Genomes (FAASG): an international initiative supporting future salmonid research, conservation and aquaculture.

We describe an emerging initiative - the 'Functional Annotation of All Salmonid Genomes' (FAASG), which will leverage the extensive trait diversity that has evolved since a whole genome duplication event in the salmonid ancestor, to develop an integrative understanding of the functional genomic basis of phenotypic variation. The outcomes of FAASG will have diverse applications, ranging from improved understanding of genome evolution, to improving the efficiency and sustainability of aquaculture production, supporting the future of fundamental and applied research in an iconic fish lineage of major societal importance.

Autopolyploidy genome duplication preserves other ancient genome duplications in Atlantic salmon (Salmo salar).

Salmonids (e.g. Atlantic salmon, Pacific salmon, and trouts) have a long legacy of genome duplication. In addition to three ancient genome duplications that all teleosts are thought to share, salmonids have had one additional genome duplication. We explored a methodology for untangling these duplications from each other to better understand them in Atlantic salmon. In this methodology, homeologous regions (paralogous/duplicated genomic regions originating from a whole genome duplication) from the most recent genome duplication were assumed to have duplicated genes at greater density and have greater sequence similarity. This assumption was used to differentiate duplicated gene pairs in Atlantic salmon that are either from the most recent genome duplication or from earlier duplications. From a comparison with multiple vertebrate species, it is clear that Atlantic salmon have retained more duplicated genes from ancient genome duplications than other vertebrates--often at higher density in the genome and containing fewer synonymous mutations. It may be that polysomic inheritance is the mechanism responsible for maintaining ancient gene duplicates in salmonids. Polysomic inheritance (when multiple chromosomes pair during meiosis) is thought to be relatively common in salmonids compared to other vertebrate species. These findings illuminate how genome duplications may not only increase the number of duplicated genes, but may also be involved in the maintenance of them from previous genome duplications as well.

The Atlantic salmon genome provides insights into rediploidization.

The whole-genome duplication 80 million years ago of the common ancestor of salmonids (salmonid-specific fourth vertebrate whole-genome duplication, Ss4R) provides unique opportunities to learn about the evolutionary fate of a duplicated vertebrate genome in 70 extant lineages. Here we present a high-quality genome assembly for Atlantic salmon (Salmo salar), and show that large genomic reorganizations, coinciding with bursts of transposon-mediated repeat expansions, were crucial for the post-Ss4R rediploidization process. Comparisons of duplicate gene expression patterns across a wide range of tissues with orthologous genes from a pre-Ss4R outgroup unexpectedly demonstrate far more instances of neofunctionalization than subfunctionalization. Surprisingly, we find that genes that were retained as duplicates after the teleost-specific whole-genome duplication 320 million years ago were not more likely to be retained after the Ss4R, and that the duplicate retention was not influenced to a great extent by the nature of the predicted protein interactions of the gene products. Finally, we demonstrate that the Atlantic salmon assembly can serve as a reference sequence for the study of other salmonids for a range of purposes.

Functional Divergence in Teleost Cardiac Troponin Paralogs Guides Variation in the Interaction of TnI Switch Region with TnC.

Gene duplication results in extra copies of genes that must coevolve with their interacting partners in multimeric protein complexes. The cardiac troponin (Tn) complex, containing TnC, TnI, and TnT, forms a distinct functional unit critical for the regulation of cardiac muscle contraction. In teleost fish, the function of the Tn complex is modified by the consequences of differential expression of paralogs in response to environmental thermal challenges. In this article, we focus on the interaction between TnI and TnC, coded for by genes that have independent evolutionary origins, but the co-operation of their protein products has necessitated coevolution. In this study, we characterize functional divergence of TnC and TnI paralogs, specifically the interrelated roles of regulatory subfunctionalization and structural subfunctionalization. We determined that differential paralog transcript expression in response to temperature acclimation results in three combinations of TnC and TnI in the zebrafish heart: TnC1a/TnI1.1, TnC1b/TnI1.1, and TnC1a/TnI1.5. Phylogenetic analysis of these highly conserved proteins identified functionally divergent residues in TnI and TnC. The structural and functional effect of these Tn combinations was modeled with molecular dynamics simulation to link divergent sites to changes in interaction strength. Functional divergence in TnI and TnC were not limited to the residues involved with TnC/TnI switch interaction, which emphasizes the complex nature of Tn function. Patterns in domain-specific divergent selection and interaction energies suggest that substitutions in the TnI switch region are crucial to modifying TnI/TnC function to maintain cardiac contraction with temperature changes. This integrative approach introduces Tn as a model of functional divergence that guides the coevolution of interacting proteins.

Genome-wide association analysis reveals loci associated with resistance against Piscirickettsia salmonis in two Atlantic salmon (Salmo salar L.) chromosomes.

BACKGROUND: Pisciricketssia salmonis is the causal agent of Salmon Rickettsial Syndrome (SRS), which affects salmon species and causes severe economic losses. Selective breeding for disease resistance represents one approach for controlling SRS in farmed Atlantic salmon. Knowledge concerning the architecture of the resistance trait is needed before deciding on the most appropriate approach to enhance artificial selection for P. salmonis resistance in Atlantic salmon. The purpose of the study was to dissect the genetic variation in the resistance to this pathogen in Atlantic salmon. METHODS: 2,601 Atlantic salmon smolts were experimentally challenged against P. salmonis by means of intra-peritoneal injection. These smolts were the progeny of 40 sires and 118 dams from a Chilean breeding population. Mortalities were recorded daily and the experiment ended at day 40 post-inoculation. Fish were genotyped using a 50K Affymetrix® Axiom® myDesignTM Single Nucleotide Polymorphism (SNP) Genotyping Array. A Genome Wide Association Analysis was performed on data from the challenged fish. Linear regression and logistic regression models were tested. RESULTS: Genome Wide Association Analysis indicated that resistance to P. salmonis is a moderately polygenic trait. There were five SNPs in chromosomes Ssa01 and Ssa17 significantly associated with the traits analysed. The proportion of the phenotypic variance explained by each marker is small, ranging from 0.007 to 0.045. Candidate genes including interleukin receptors and fucosyltransferase have been found to be physically linked with these genetic markers and may play an important role in the differential immune response against this pathogen. CONCLUSIONS: Due to the small amount of variance explained by each significant marker we conclude that genetic resistance to this pathogen can be more efficiently improved with the implementation of genetic evaluations incorporating genotype information from a dense SNP array.

Expression analysis of sex-determining pathway genes during development in male and female Atlantic salmon (Salmo salar).

We studied the expression of 28 genes that are involved in vertebrate sex-determination or sex-differentiation pathways, in male and female Atlantic salmon (Salmo salar) in 11 stages of development from fertilization to after first feeding. Gene expression was measured in half-sibs that shared the same dam. The sire of family 1 was a sex-reversed female (i.e., genetically female but phenotypically male), and so the progeny of this family are all female. The sire of family 2 was a true male, and so the offspring were 50% male and 50% female. Gene expression levels were compared among three groups: 20 female offspring of the cross between a regular female and the sex-reversed female (family 1, first group), ∼ 10 females from the cross between a regular female and a regular male (family 2, second group) and ∼ 10 males from this same family (family 2, third group). Statistically significant differences in expression levels between males and the two groups of females were observed for two genes, gsdf and amh/mis, in the last four developmental stages examined. SdY, the sex-determining gene in rainbow trout, appeared to be expressed in males from 58 days postfertilization (dpf). Starting at 83 dpf, ovarian aromatase, cyp19a, expression appeared to be greater in both groups of females compared with males, but this difference was not statistically significant. The time course of expression suggests that sdY may be involved in the upregulation of gsdf and amh/mis and the subsequent repression of cyp19a in males via the effect of amh/mis.

Genomic Instability of the Sex-Determining Locus in Atlantic Salmon (Salmo salar).

Atlantic salmon and rainbow trout, like other members of the subfamily Salmoninae, are gonochoristic with male heterogamety. The finding that sex-linked genetic markers varied between species suggested that the sex-determining gene differs among salmonid species, or that there is one sex-determining gene that has the capacity to move around the genome. The discovery of sdY, the sex-determining gene in rainbow trout, and its presence in many male salmonids gave support to the latter. Additional evidence for a salmonid-specific, sex-determining jumping gene came from the mapping of the sex-determining locus to three different chromosomes in Tasmanian male Atlantic salmon lineages. To characterize the sex-determining region, we isolated three sdY containing BACs from an Atlantic salmon male library. Sequencing of these BACs yielded two contigs, one of which contained the sdY gene. Sequence analysis of the borders of male-specific and female/male common regions revealed highly repetitive sequences associated with mobile elements, which may allow an sdY cassette to jump around the genome. FISH analysis using a BAC or a plasmid containing the sdY gene showed that the sdY gene did indeed localize to the chromosomes where SEX had been mapped in different Tasmanian Atlantic salmon families. Moreover, the plasmid sdY gene probe hybridized primarily to one of the sex chromosomes as would be expected of a male-specific gene. Our results suggest that a common salmonid sex-determining gene (sdY) can move between three specific loci on chromosomes 2, 3, and 6, giving the impression that there are multiple SEX loci both within and between salmonid species.

Epithelial Cadherin Determines Resistance to Infectious Pancreatic Necrosis Virus in Atlantic Salmon.

Infectious pancreatic necrosis virus (IPNV) is the cause of one of the most prevalent diseases in farmed Atlantic salmon (Salmo salar). A quantitative trait locus (QTL) has been found to be responsible for most of the genetic variation in resistance to the virus. Here we describe how a linkage disequilibrium-based test for deducing the QTL allele was developed, and how it was used to produce IPN-resistant salmon, leading to a 75% decrease in the number of IPN outbreaks in the salmon farming industry. Furthermore, we describe how whole-genome sequencing of individuals with deduced QTL genotypes was used to map the QTL down to a region containing an epithelial cadherin (cdh1) gene. In a coimmunoprecipitation assay, the Cdh1 protein was found to bind to IPNV virions, strongly indicating that the protein is part of the machinery used by the virus for internalization. Immunofluorescence revealed that the virus colocalizes with IPNV in the endosomes of homozygous susceptible individuals but not in the endosomes of homozygous resistant individuals. A putative causal single nucleotide polymorphism was found within the full-length cdh1 gene, in phase with the QTL in all observed haplotypes except one; the absence of a single, all-explaining DNA polymorphism indicates that an additional causative polymorphism may contribute to the observed QTL genotype patterns. Cdh1 has earlier been shown to be necessary for the internalization of certain bacteria and fungi, but this is the first time the protein is implicated in internalization of a virus.

Genome-wide association study (GWAS) for growth rate and age at sexual maturation in Atlantic salmon (Salmo salar).

Early sexual maturation is considered a serious drawback for Atlantic salmon aquaculture as it retards growth, increases production times and affects flesh quality. Although both growth and sexual maturation are thought to be complex processes controlled by several genetic and environmental factors, selection for these traits has been continuously accomplished since the beginning of Atlantic salmon selective breeding programs. In this genome-wide association study (GWAS) we used a 6.5K single-nucleotide polymorphism (SNP) array to genotype ∼ 480 individuals from the Cermaq Canada broodstock program and search for SNPs associated with growth and age at sexual maturation. Using a mixed model approach we identified markers showing a significant association with growth, grilsing (early sexual maturation) and late sexual maturation. The most significant associations were found for grilsing, with markers located in Ssa10, Ssa02, Ssa13, Ssa25 and Ssa12, and for late maturation with markers located in Ssa28, Ssa01 and Ssa21. A lower level of association was detected with growth on Ssa13. Candidate genes, which were linked to these genetic markers, were identified and some of them show a direct relationship with developmental processes, especially for those in association with sexual maturation. However, the relatively low power to detect genetic markers associated with growth (days to 5 kg) in this GWAS indicates the need to use a higher density SNP array in order to overcome the low levels of linkage disequilibrium observed in Atlantic salmon before the information can be incorporated into a selective breeding program.

Genomic organization and evolution of the trace amine-associated receptor (TAAR) repertoire in Atlantic salmon (Salmo salar).

There is strong evidence that olfaction plays a key role in the homing of salmonids to their natal spawning grounds, particularly in the freshwater phase. However, the physiological and genetic mechanisms behind this biological phenomenon are largely unknown. It has been shown that Pacific salmon respond to dissolved free amino acids from their natal streams. This indicates that amino acids comprise part of the olfcatory cues for imprinting and homing in salmonids. As trace amine-associated receptors (TAARs), a class of olfactory receptors that are close relatives of the G protein-coupled aminergic neurotransmitter receptors, recognize amino acid metabolites, we hypothesize that TAARs play an important role in salmon homing by recognizing olfactory cues. Therefore, to better understand homing in Atlantic salmon, we set out to characterize the TAAR genes in this species. We searched the first assembly of the Atlantic salmon genome for sequences resembling TAARs previously characterized in other teleosts. We identified 27 putatively functional TAAR genes and 25 putative TAAR pseudogenes, which cluster primarily on chromosome 21 (Ssa21). Phylogenetic analysis of TAAR amino acid sequences from 15 vertebrate species revealed the TAAR gene family arose after the divergence of jawed and jawless vertebrates. The TAARs group into three classes with salmon possessing class I and class III TAARs. Within each class, evolution is characterized by species-specific gene expansions, which is in contrast to what is observed in other olfactory receptor families (e.g., OlfCs and oras).

Sex-specific expression and localization of aromatase and its regulators during embryonic and larval development of Atlantic salmon.

The products of dax1, foxl2a and mis have each been shown to have proliferative and/or differentiative activities during mammalian organogenesis. These factors also play a role in regulating the biosynthesis of estrogen, particularly by modulating the activity of aromatase cyp19a. We demonstrate the transcription and translation of these genes during salmon embryogenesis. We were able to track sex-specific differences in these processes through accurate determination of the sex of each embryo and larva examined from genotyped microsatellites. We detected sex- and stage-specific immunolabeling of the embryonic gut, kidney, gonads, neural cord and skeletal muscle by DAX-1, FOXL2A and MIS. These results indicate the potential of these factors to mediate proliferation and/or differentiation programs during development of these tissues. As well, immunolabeling of skeletal muscle by CYP19B1 throughout the study reveals probable neurogenic activity associated with peripheral radial glial cells and the growing embryonic musculature.