Entropy Based Approach to Finding Interacting Genes Responsible for Complex Human Disease

2000 Mathematics Subject Classification: 62P10, 92D10, 92D30, 94A17, 62L10.

NOTCH1 Gain of Function in Germ Cells Causes Failure of Spermatogenesis in Male Mice

NOTCH1 is a member of the NOTCH receptor family, a group of single-pass trans-membrane receptors. NOTCH signaling is highly conserved in evolution and mediates communication between adjacent cells. NOTCH receptors have been implicated in cell fate determination, as well as maintenance and differentiation of stem cells. In the mammalian testis expression of NOTCH1 in somatic and germ cells has been demonstrated, however its role in spermatogenesis was not clear. To study the significance of NOTCH1 in germ cells, we applied a cre/loxP approach in mice to induce NOTCH1 gain- or loss-of function specifically in male germ cells. Using a Stra8-icretransgene we produced mice with conditional activation of the NOTCH1 intracellular domain (NICD) in germ cells. Spermatogenesis in these mutants was progressively affected with age, resulting in decreased testis weight and sperm count. Analysis of downstream target genes of NOTCH1 signaling showed an increased expression of Hes5, with a reduction of the spermatogonial differentiation marker, Neurog3 expression in the mutant testis. Apoptosis was significantly increased in mouse germ cells with the corresponding elevation of pro-apoptotic Trp53 and Trp63genes' expression. We also showed that the conditional germ cell-specific ablation of Notch1 had no effect on spermatogenesis or male fertility. Our data suggest the importance of NOTCH signaling regulation in male germ cells for their survival and differentiation.

Multiplex PCR for Assessment of Tetracycline Resistance (tet) Genes in Antibiotic-Resistant Soil Bacteria and the Ecological Implications

Antibiotics are becoming increasingly prevalent in bacterial communities due to clinical and agricultural misuse and overuse in their environment. As exposure increases, so does the incidence of microbial resistance. Such is the case with bacterial resistance to tetracyclines, a phenotype often acquired through the horizontal gene transfer of tet genes between bacteria. The objective of this project was to analyze the bacterial diversity of tet resistance genes in soil from Miami-Dade County. Bacterial isolates were Gram-stained and the Kirby-Bauer antibiotic disk diffusion test was performed to determine each bacterium’s degree of resistance. The 16S rRNA gene from antibiotic-resistant isolates was amplified by PCR and sequenced to identify the isolates. All isolates’ tet genes were amplified by multiplex PCR, sequenced, and compared. Among eight isolates, three distinct species were positively identified based on their 16S rRNA sequences and four distinct tet genes were identified, though all tested susceptible to tetracycline via the Kirby-Bauer test. This project clarifies some aspects of the ecology of antibiotic resistance genes, their natural ecological function and the potential for the expansion of intrinsic multi-antibiotic resistance into new ecosystems and/or hosts.

Function of the PDLIM2 protein in oncogenic Wnt signalling pathways

Aberrant regulation of the Wnt signalling pathway is a recurrent theme in cancer biology. Hyper activation due to oncogenic mutations and paracrine activity has been found in both colon cancer and breast cancer, and continues to evolve as a central mechanism in oncogenesis. PDLIM2, a cytoskeletal PDZ protein, is an IGF-1 regulated gene that is highly expressed in cancer cell lines derived from metastatic tumours. Suppression of PDLIM2 inhibits polarized cell migration, reverses the Epithelial to Mesenchymal transition (EMT) phenotype, suppresses the transcription of β-catenin target genes, and regulates gene expression of key transcription factors in EMT. This thesis investigates the mechanism by which PDLIM2 contributes to the maintenance of Wnt signalling in cancer cells. Here we show that PDLIM2 is a critical regulator of the Wnt pathway by regulating β-catenin at the adherens juctions, as also its transcriptional activity by the interaction of PDLIM2 with TCF4 at the nucleus. Evaluation of PDLIM2 in macrophages and co-culture studies with cancer cells and fibroblasts showed the influence exerted on PDLIM2 by paracrine cues. Thus, PDLIM2 integrates cytoskeleton signalling with gene expression by modulating the Wnt signalling pathway and reconciling microenvironmental cues with signals in epithelial cells. Negative correlation of mRNA and protein levels in the triple negative breast cancer cell BT549 suggests that PDLIM2 is part of a more complex mechanism that involves transcription and posttranslational modifications. GST pulldown studies and subsequent mass spectrometry analysis showed that PDLIM2 interacts with 300 proteins, with a high biological function in protein biosynthesis and Ubiquitin/proteasome pathways, including 13 E3 ligases. Overall, these data suggest that PDLIM2 has two distinct functions depending of its location. Located at the cytoplasm mediates cytoskeletal re-arrangements, whereas at the nucleus PDLIM2 acts as a signal transduction adaptor protein mediating transcription and ubiquitination of key transcription factors in cancer development.

Obestatin as a key regulator of metabolism and cardiovascular function with emerging therapeutic potential for diabetes

Obestatin is a 23-amino acid C-terminally amidated gastrointestinal peptide derived from preproghrelin and which forms an alpha helix. Although obestatin has a short biological half-life and is rapidly degraded, it is proposed to exert wide-ranging pathophysiological actions. Whilst the precise nature of many of its effects is unclear, accumulating evidence supports positive actions on both metabolism and cardiovascular function. For example, obestatin has been reported to inhibit food and water intake, body weight gain, and gastrointestinal motility, and to also mediate promotion of cell survival and prevention of apoptosis. Obestatin-induced increases in β-cell mass, enhanced adipogenesis and improved lipid metabolism have been noted along with upregulation of genes associated with β-cell regeneration, insulin production and adipogenesis. Furthermore, human circulating obestatin levels generally demonstrate an inverse association with obesity and diabetes, whilst the peptide has been shown to confer protective metabolic effects in experimental diabetes, suggesting that it may hold therapeutic potential in this setting. Obestatin also appears to be involved in blood pressure regulation and to exert beneficial effects on endothelial function, with experimental studies indicating that it may also promote cardioprotective actions against, for example, ischaemia-reperfusion injury. This review will present a critical appraisal of the expanding obestatin research area and discuss the emerging therapeutic potential of this peptide for both metabolic and cardiovascular complications of diabetes. 

CTCF modulates Estrogen Receptor function through specific chromatin and nuclear matrix interactions

Enhancer regions and transcription start sites of estrogen-target regulated genes are connected by means of Estrogen Receptor long-range chromatin interactions. Yet, the complete molecular mechanisms controlling the transcriptional output of engaged enhancers and subsequent activation of coding genes remain elusive. Here, we report that CTCF binding to enhancer RNAs is enriched when breast cancer cells are stimulated with estrogen. CTCF binding to enhancer regions results in modulation of estrogen-induced gene transcription by preventing Estrogen Receptor chromatin binding and by hindering the formation of additional enhancer-promoter ER looping. Furthermore, the depletion of CTCF facilitates the expression of target genes associated with cell division and increases the rate of breast cancer cell proliferation. We have also uncovered a genomic network connecting loci enriched in cell cycle regulator genes to nuclear lamina that mediates the CTCF function. The nuclear lamina and chromatin interactions are regulated by estrogen-ER. We have observed that the chromatin loops formed when cells are treated with estrogen establish contacts with the nuclear lamina. Once there, the portion of CTCF associated with the nuclear lamina interacts with enhancer regions, limiting the formation of ER loops and the induction of genes present in the loop. Collectively, our results reveal an important, unanticipated interplay between CTCF and nuclear lamina to control the transcription of ER target genes, which has great implications in the rate of growth of breast cancer cells.

A Rare Loss-of-Function Variant in Scavenger Receptor Class B Type I (SCARB1) Raises HDL Cholesterol and Increases Risk of Coronary Disease

Scavenger receptor BI (SR-BI) is the major receptor for high-density lipoprotein (HDL)
cholesterol (HDL-C). In humans, high amounts of HDL-C in plasma are associated with a
lower risk of coronary heart disease (CHD). Mice that have depleted Scarb1 (SR-BI
knockout mice) have markedly elevated HDL-C levels but, paradoxically, increased
atherosclerosis. The impact of SR-BI on HDL metabolism and CHD risk in humans remains
unclear. Through targeted sequencing of coding regions of lipid-modifying genes in 328
individuals with extremely high plasma HDL-C levels, we identified a homozygote for a lossof-function
variant, in which leucine replaces proline 376 (P376L), in SCARB1, the gene
encoding SR-BI. The P376L variant impairs posttranslational processing of SR-BI and
abrogates selective HDL cholesterol uptake in transfected cells, in hepatocyte-like cells
derived from induced pluripotent stem cells from the homozygous subject, and in mice.
Large population-based studies revealed that subjects who are heterozygous carriers of
the P376L variant have significantly increased levels of plasma HDL-C. P376L carriers have
a profound HDL-related phenotype and an increased risk of CHD (odds ratio = 1.79, which is
statistically significant).

Investigating the Function of a Small Secreted Protein Family in Physcomitrella Patens

The correct development of multicellular organisms depends upon the perception of signals secreted by cells in order to co-ordinate cell differentiation. The Physcomitrella patens genome encodes many components of potential signaling systems, including putative receptor proteins and putative secreted protein ligands, yet at present little characterization of these proteins has been carried out. We are currently attempting to characterize the expression pattern and function of a family of 6 secreted proteins exhibiting homology to PrsS, the ligand that controls self-incompatibility (SI) in Papaver rhoeas (field poppy). In poppy, PrsS interacts a receptor on the surface of pollen tubes, PrpS causing SI by programmed cell death. Homologues of this protein (SPH – S-Protein Homologues) exist in dicotyledonous plants and bryophytes but not in other plant taxa. We aim to determine spatiotemporal expression differences between these proteins via reporter gene analysis and qPCR of cDNA. In addition we are in the process of creating targeted gene knockouts for all 6 of the genes in P. patens. We are also searching for receptors of PrpS in Physcomitrella using a bioinformatic strategy alongside phage display. In accomplishing this we hope to determine the function of a small novel secreted protein family in Physcomitrella but in addition we also hope to elucidate the function of SPH proteins in Arabidopsis.

Identification of Genes Modulating Mitochondrial Biogenesis in Caenorhabditis elegans

Since Altmann recognized ubiquitously distributed "bioblasts" in 1890, understanding of mitochondria has evolved from "elementary organisms" living inside cells and carrying out vital functions, over the Harman's "free radical theory" in 1956, to one of the driving forces of aging and cause of multiple associated diseases impacting society today. While a tremendous amount of work has contributed to the understanding of mitochondrial biology in different model organisms, the precise molecular mechanisms of basic mitochondrial function have yet to be deciphered. By employing an RNA interference mediated screen in Caenorhabditis elegans, we identified two transcription factors: SPTF-3, a member of Sp1 family, and an uncharacterized, nematode specific W04D2.4. We propose that both proteins modulate expression of many genes with regard to mitochondrial function including mitochondrial single-stranded binding protein encoded by mtss-1, whose promoter was used as transcriptional reporter in the screen. Further, RNA sequencing data indicate that W04D2.4 indirectly regulates expression of mitochondrial DNA via control of genes functionally related to mitochondrial replication and translation machineries. We also demonstrate that from all interventions targeting cytosolic translation, MTSS-1 levels are elevated only upon knockdown of genes encoding cytosolic ribosomal proteins. Reduction of ribosomes leads to increased sptf-3 translation, most likely in an internal ribosome entry side (IRES) mediated manner, eventually inducing mtss-1 expression. Moreover, we identify a novel role for SPTF-3 in the regulation of mitochondrial unfolded stress response (UPRmt) activation, but not endoplasmatic reticulum or oxidative stress responses. Taken together, this study identifies two transcription factors previously not associated with mitochondrial biogenesis and UPRmt in C. elegans, establishing a basis for further investigation of mito-nuclear interactions.

Antioxidant and detoxify genes polymorphisms in colorectal cancer predisposition

Colorectal cancer (CRC) results from histologic and gene alterations can lead to a massive cellular proliferation. Most of the authors assume multifactorial causes to CRC genesis. Low physical activity, a fat diet poor in fibers and smoking habits seems to have an important role in CRC. However, there are also genetic causes associated with CRC risk. It has been described that oxidative stress levels could influence CRC development. Thus, cellular balance reactive species and defense enzymes involved in oxidative stress are crucial to maintain a good tissue function and avoid neoplasic process. Therefore, genome variations on these defense enzymes, such as MNSOD, SOD3, GSTP1, GSTT1 and GSTM1, could be important biomarkers to colorectal adenocarcinomas. We intend to determine frequencies distribution of most common polymorphisms involved on oxidative stress regulation (MNSOD, SOD3, GSTP1, GSTT1 and GSTM1) in patients with sporadic colorectal adenocarcinoma (SCA) and in healthy controls, evaluation their possible correlation with SCA risk. Samples common polymorphisms of antioxidant and detoxify genes (MNSOD T175C, SOD3 R213G, GSTP1 A105G, GSTP1 C114T, GSTT1del and GSTM1del) analysis was done by PCR-SSP techniques. In this study we found a higher prevalence of MNSOD 175CC (55% vs 2%; p<0.0001; OR: 58.5; CI 13.3 to 256.7), SOD3 213GG (31% vs 2%; p<0.0001; OR: 21.89; CI 4.93 to 97.29), GSTP1 105GG (46% vs 12%; p<0.0001; OR: 6.14; CI 2.85 to 13.26), GSTP1 114TT (38% vs 0%; p<0.0001; OR: Infinity) and GSTT1 null (75% vs 28%; p<0.0001; OR: 7.71; CI 3.83 to 15.56) mutated genotypes among SCA patients, while the normal genotypes were associated with SCA absence. Furthermore, we found GSTP1 114TT mutated genotype (52% vs 27%; p=0.003; OR: 2.88; CI: 1.41 to 5.89) and GSTT1 null genotype (87% vs 65%; p=0.003; OR: 3.66; CI 1.51 to 8.84) associated with colon samples. These findings suggest a positive association between most of common polymorphisms involved on oxidative stress regulation and SCA prevalence. Dysregulation of MNSOD, SOD3, GSTP1, GSTT1 and GSTM1 genes could be associated with an increase of ROS in colon and rectum tissue and p53 pathway deregulation, induced by oxidative stress on colonic and rectal cells. The present study also provides preliminary evidence that MNSOD 175C, SOD3 213G, GSTP1 105G, GSTP1 114T and GSTT1 null polymorphisms, may be involved in SCA risk and could be useful to clarify this multifactorial disorder.

Identification and functional analysis of novel genes expressed in the Anterior Visceral Endoderm

During early vertebrate development, the correct establishment of the body axes is critical. The anterior pole of the mouse embryo is established when Distal Visceral Endoderm (DVE) cells migrate to form the Anterior Visceral Endoderm (AVE). Symmetrical expression of Lefty1, Cer1 and Dkk1 determines the direction of DVE migration and the future anterior side. In addition to the establishment of the Anterior-Posterior axis, the AVE has also been implicated in anterior neural specification. To better understand the role of the AVE in these processes, we have performed a differential screening using Affymetrix GeneChip technology with AVE cells isolated from cer1P-EGFP transgenic mouse embryos. We found 175 genes which were upregulated in the AVE and 36 genes in the Proximal-posterior sample. Using DAVID software, we characterized the AVE cell population regarding cellular component, molecular function and biological processes. Among the genes that were found to be upregulated in the AVE, several novel genes were identified. Four of these transcripts displaying high-fold change in the AVE were further characterized by in situ hybridization in early stages of development in order to validate the screening. From those four selected genes, one, denominated Adtk1, was chosen to be functionally characterized by targeted inactivation in ES cells. Adtk1 encodes for a serine/threonine kinase. Adtk1 null mutants are smaller and present short limbs due to decreased mineralization, suggesting a potential role in chondrogenesis during limb development. Taken together, these data point to the importance of reporting novel genes present in the AVE.

Association between common lipid metabolism genes polymorphisms and sporadic colorectal adenocarcinoma risk

Lipids can modulate the risk of developing sporadic colorectal adenocarcinoma (SCA), since alterations into lipid metabolism and transport pathways influence directly cholesterol and lipids absorption by colonic cells and indirectly reactive oxygen species (ROS) synthesis in rectum cells due to lipid accumulation. Lipid metabolism is regulated by several proteins APOA1, APOB, APOC3, APOE, CETP, NPY, PON1 and PPARG that could influence both metabolism and transport processes. Is been reported that several common single-nucleotide polymorphisms (SNPs) in these genes could influence their function and/or expression, changing lipid metabolism balance. Thus, genetic changes in those genes can influence SCA development, once the majority of them were never studied in this disease. Furthermore, there are contradictory results between some studied polymorphisms and SCA risk. Thus, the aim of this study was to explore and describe lipid metabolism-associated genes common polymorphisms (APOA1 -75 G>A; APOB R3500Q; APOC3 C3175G, APOC3 T3206G; APOE Cys112/158Arg; CETP G279A, CETP R451Q; NPY Leu7Pro; PON1 Q192R; PPARG Pro12Ala) status among SCA, and their relationship with SCA risk. Genotyping of common lipid metabolism genes polymorphisms (APOA1 75 G>A; APOB R3500Q; APOC3 C3175G, APOC3 T3206G; APOE Cys112/158Arg; CETP G279A, CETP R451Q; NPY Leu7Pro; PON1 Q192R; PPARG Pro12Ala) were done by PCR-SSP techniques, from formalin-fixed and paraffin-embedded biopsies of 100 healthy individuals and 68 SCA subjects. Mutant genotypes of APOA1 -75AA (32% vs 12%; p=0.001; OR=3.51; 95% CI 1.59-7.72); APOB 3500AA (7% vs 0%; p=0.01); APOC3 3175GG (19% vs 2%; p=0.0002; OR=11.58; 95% CI 2.52-53.22), APOC3 3206GG (19% vs 0%; p<0.0001); CETP 279AA (12% vs 1%; p=0.003; OR=13.20; 95% CI 1.61-108.17), CETP 451AA (16% vs 0%; p<0.0001); NPY 7CC (15% vs 0%; p<0.0001); PPARG 12GG (10% vs 0%; p=0.001); and heterozygote genotype PON1 192AG (56% vs 22%; p<0.0001; OR=4.49; 95% CI 2.298.80) were found associated with SCA prevalence. While, APOE E4/E4 (0% vs 8%; p=0.02) mutant haplotype seemed to have a protective effect on SCA. Moreover, it also been founded differences between APOB 3500GA, APOC3 3206TG, CETP 279AA genotypes and PPARG 12Ala allele prevalence and tissue localization (colon vs rectum). These findings suggest a positive association between most of common lipid metabolism genes polymorphisms studied and SCA prevalence. Dysregulation of APOA1, APOB, APOC3, CETP, NPY, PON1 and PPARG genes could be associated with lower cholesterol plasma levels and increase ROS among colon and rectum mucosa. Furthermore, these results also support the hypothesis that CRC is related with intestinal lipid absorption decrease and secondary bile acids production increase. Moreover, the polymorphisms studied may play an important role as biomarkers to SCA susceptibility.

Association between common lipid metabolism genes polymorphisms and sporadic colorectal adenocarcinoma risk

Lipids can modulate the risk of developing sporadic colorectal adenocarcinoma (SCA), since alterations into lipid metabolism and transport pathways influence directly cholesterol and lipids absorption by colonic cells and indirectly reactive oxygen species (ROS) synthesis in rectum cells due to lipid accumulation. Lipid metabolism is regulated by several proteins APOA1, APOB, APOC3, APOE, CETP, NPY, PON1 and PPARG that could influence both metabolism and transport processes. Is been reported that several common single-nucleotide polymorphisms (SNPs) in these genes could influence their function and/or expression, changing lipid metabolism balance. Thus, genetic changes in those genes can influence SCA development, once the majority of them were never studied in this disease. Furthermore, there are contradictory results between some studied polymorphisms and SCA risk. Thus, the aim of this study was to explore and describe lipid metabolism-associated genes common polymorphisms (APOA1 -75 G>A; APOB R3500Q; APOC3 C3175G, APOC3 T3206G; APOE Cys112/158Arg; CETP G279A, CETP R451Q; NPY Leu7Pro; PON1 Q192R; PPARG Pro12Ala) status among SCA, and their relationship with SCA risk. Genotyping of common lipid metabolism genes polymorphisms (APOA1 75 G>A; APOB R3500Q; APOC3 C3175G, APOC3 T3206G; APOE Cys112/158Arg; CETP G279A, CETP R451Q; NPY Leu7Pro; PON1 Q192R; PPARG Pro12Ala) were done by PCR-SSP techniques, from formalin-fixed and paraffin-embedded biopsies of 100 healthy individuals and 68 SCA subjects. Mutant genotypes of APOA1 -75AA (32% vs 12%; p=0.001; OR=3.51; 95% CI 1.59-7.72); APOB 3500AA (7% vs 0%; p=0.01); APOC3 3175GG (19% vs 2%; p=0.0002; OR=11.58; 95% CI 2.52-53.22), APOC3 3206GG (19% vs 0%; p<0.0001); CETP 279AA (12% vs 1%; p=0.003; OR=13.20; 95% CI 1.61-108.17), CETP 451AA (16% vs 0%; p<0.0001); NPY 7CC (15% vs 0%; p<0.0001); PPARG 12GG (10% vs 0%; p=0.001); and heterozygote genotype PON1 192AG (56% vs 22%; p<0.0001; OR=4.49; 95% CI 2.298.80) were found associated with SCA prevalence. While, APOE E4/E4 (0% vs 8%; p=0.02) mutant haplotype seemed to have a protective effect on SCA. Moreover, it also been founded differences between APOB 3500GA, APOC3 3206TG, CETP 279AA genotypes and PPARG 12Ala allele prevalence and tissue localization (colon vs rectum). These findings suggest a positive association between most of common lipid metabolism genes polymorphisms studied and SCA prevalence. Dysregulation of APOA1, APOB, APOC3, CETP, NPY, PON1 and PPARG genes could be associated with lower cholesterol plasma levels and increase ROS among colon and rectum mucosa. Furthermore, these results also support the hypothesis that CRC is related with intestinal lipid absorption decrease and secondary bile acids production increase. Moreover, the polymorphisms studied may play an important role as biomarkers to SCA susceptibility.

Identification of superior reference genes for data normalisation of expression studies via quantitative PCR in hybrid roses (Rosa hybrida)

Background: Gene expression studies are a prerequisite for understanding the biological function of genes. Because of its high sensitivity and easy use, quantitative PCR (qPCR) has become the gold standard for gene expression quantification. To normalise qPCR measurements between samples, the most prominent technique is the use of stably expressed endogenous control genes, the so called reference genes. However, recent studies show there is no universal reference gene for all biological questions. Roses are important ornamental plants for which there has been no evaluation of useful reference genes for gene expression studies. Results: We used three different algorithms (BestKeeper, geNorm and NormFinder) to validate the expression stability of nine candidate reference genes in different rose tissues from three different genotypes of Rosa hybrida and in leaves treated with various stress factors. The candidate genes comprised the classical "housekeeping genes" (Actin, EF-1α, GAPDH, Tubulin and Ubiquitin), and genes showing stable expression in studies in Arabidopsis (PP2A, SAND, TIP and UBC). The programs identified no single gene that showed stable expression under all of the conditions tested, and the individual rankings of the genes differed between the algorithms. Nevertheless the new candidate genes, specifically, PP2A and UBC, were ranked higher as compared to the other traditional reference genes. In general, Tubulin showed the most variable expression and should be avoided as a reference gene. Conclusions: Reference genes evaluated as suitable in experiments with Arabidopsis thaliana were stably expressed in roses under various experimental conditions. In most cases, these genes outperformed conventional reference genes, such as EF1-α and Tubulin. We identified PP2A, SAND and UBC as suitable reference genes, which in different combinations may be used for normalisation in expression analyses via qPCR for different rose tissues and stress treatments. However, the vast genetic variation found within the genus Rosa, including differences in ploidy levels, might also influence expression stability of reference genes, so that future research should also consider different genotypes and ploidy levels.

Identification and expression analysis of two Wnt4 genes in the spotted scat ( Scatophagus argus )

Background: WNT4 is a protein that plays a crucial role in ovarian differentiation and development in mammals, with a relatively well understood function in mammalian gonadal differentiation. The role of Wnt4 in teleost fish; however, remains unclear. In the present study, cDNAs of Wnt4a and Wnt4b were cloned and characterized in the spotted scat. The expression patterns of two Wnt4 genes in the gonads at different stages of development and in fish after treatment with 17α-methyltestosterone (MT) were investigated. Results: The tissue distribution showed that Wnt4a was expressed in various tissues, including the gonads, gills, spleen, brain, and fin. Interestingly, Wnt4b not only was expressed in the gills, brain, and spleen, but also was obviously expressed in the ovary. During gonad development, Wnt4a was highly expressed in the testis at stage I and Wnt4b was mainly expressed in the ovary at stages II–III. After MT treatment, the mRNA expression of Wnt4a increased significantly up to 40 d, and the transcript level of Wnt4b decreased at 20 d. Conclusions: These results suggest that Wnt4a may be involved in gonad development and plays a role in the process of spermatogonial proliferation. Our results also demonstrate that Wnt4b is not only expressed in the nervous system, but also in the ovary and it may be involved in ovary development of the spotted scat.