Constitutive gene expression profile segregates toxicity in locally advanced breast cancer patients treated with high-dose hyperfractionated radical radiotherapy

[EN] Breast cancer patients show a wide variation in normal tissue reactions after radiotherapy. The individual sensitivity to x-rays limits the efficiency of the therapy. Prediction of individual sensitivity to radiotherapy could help to select the radiation protocol and to improve treatment results. The aim of this study was to assess the relationship between gene expression profiles of ex vivo un-irradiated and irradiated lymphocytes and the development of toxicity due to high-dose hyperfractionated radiotherapy in patients with locally advanced breast cancer. Raw data from microarray experiments were uploaded to the Gene Expression Omnibus Database http://www.ncbi.nlm.nih.gov/geo/ (GEO accession GSE15341). We obtained a small group of 81 genes significantly regulated by radiotherapy, lumped in 50 relevant pathways. Using ANOVA and t-test statistical tools we found 20 and 26 constitutive genes (0 Gy) that segregate patients with and without acute and late toxicity, respectively. Non-supervised hierarchical clustering was used for the visualization of results. Six and 9 pathways were significantly regulated respectively. Concerning to irradiated lymphocytes (2 Gy), we founded 29 genes that separate patients with acute toxicity and without it. Those genes were gathered in 4 significant pathways. We could not identify a set of genes that segregates patients with and without late toxicity. In conclusion, we have found an association between the constitutive gene expression profile of peripheral blood lymphocytes and the development of acute and late toxicity in consecutive, unselected patients. These observations suggest the possibility of predicting normal tissue response to irradiation in high-dose non-conventional radiation therapy regimens. Prospective studies with higher number of patients are needed to validate these preliminary results.

New approaches to open problems in gene expression microarray data

In the past decade, the advent of efficient genome sequencing tools and high-throughput experimental biotechnology has lead to enormous progress in the life science. Among the most important innovations is the microarray tecnology. It allows to quantify the expression for thousands of genes simultaneously by measurin the hybridization from a tissue of interest to probes on a small glass or plastic slide. The characteristics of these data include a fair amount of random noise, a predictor dimension in the thousand, and a sample noise in the dozens. One of the most exciting areas to which microarray technology has been applied is the challenge of deciphering complex disease such as cancer. In these studies, samples are taken from two or more groups of individuals with heterogeneous phenotypes, pathologies, or clinical outcomes. these samples are hybridized to microarrays in an effort to find a small number of genes which are strongly correlated with the group of individuals. Eventhough today methods to analyse the data are welle developed and close to reach a standard organization (through the effort of preposed International project like Microarray Gene Expression Data -MGED- Society [1]) it is not unfrequant to stumble in a clinician's question that do not have a compelling statistical method that could permit to answer it.The contribution of this dissertation in deciphering disease regards the development of new approaches aiming at handle open problems posed by clinicians in handle specific experimental designs. In Chapter 1 starting from a biological necessary introduction, we revise the microarray tecnologies and all the important steps that involve an experiment from the production of the array, to the quality controls ending with preprocessing steps that will be used into the data analysis in the rest of the dissertation. While in Chapter 2 a critical review of standard analysis methods are provided stressing most of problems that In Chapter 3 is introduced a method to adress the issue of unbalanced design of miacroarray experiments. In microarray experiments, experimental design is a crucial starting-point for obtaining reasonable results. In a two-class problem, an equal or similar number of samples it should be collected between the two classes. However in some cases, e.g. rare pathologies, the approach to be taken is less evident. We propose to address this issue by applying a modified version of SAM [2]. MultiSAM consists in a reiterated application of a SAM analysis, comparing the less populated class (LPC) with 1,000 random samplings of the same size from the more populated class (MPC) A list of the differentially expressed genes is generated for each SAM application. After 1,000 reiterations, each single probe given a "score" ranging from 0 to 1,000 based on its recurrence in the 1,000 lists as differentially expressed. The performance of MultiSAM was compared to the performance of SAM and LIMMA [3] over two simulated data sets via beta and exponential distribution. The results of all three algorithms over low- noise data sets seems acceptable However, on a real unbalanced two-channel data set reagardin Chronic Lymphocitic Leukemia, LIMMA finds no significant probe, SAM finds 23 significantly changed probes but cannot separate the two classes, while MultiSAM finds 122 probes with score >300 and separates the data into two clusters by hierarchical clustering. We also report extra-assay validation in terms of differentially expressed genes Although standard algorithms perform well over low-noise simulated data sets, multi-SAM seems to be the only one able to reveal subtle differences in gene expression profiles on real unbalanced data. In Chapter 4 a method to adress similarities evaluation in a three-class prblem by means of Relevance Vector Machine [4] is described. In fact, looking at microarray data in a prognostic and diagnostic clinical framework, not only differences could have a crucial role. In some cases similarities can give useful and, sometimes even more, important information. The goal, given three classes, could be to establish, with a certain level of confidence, if the third one is similar to the first or the second one. In this work we show that Relevance Vector Machine (RVM) [2] could be a possible solutions to the limitation of standard supervised classification. In fact, RVM offers many advantages compared, for example, with his well-known precursor (Support Vector Machine - SVM [3]). Among these advantages, the estimate of posterior probability of class membership represents a key feature to address the similarity issue. This is a highly important, but often overlooked, option of any practical pattern recognition system. We focused on Tumor-Grade-three-class problem, so we have 67 samples of grade I (G1), 54 samples of grade 3 (G3) and 100 samples of grade 2 (G2). The goal is to find a model able to separate G1 from G3, then evaluate the third class G2 as test-set to obtain the probability for samples of G2 to be member of class G1 or class G3. The analysis showed that breast cancer samples of grade II have a molecular profile more similar to breast cancer samples of grade I. Looking at the literature this result have been guessed, but no measure of significance was gived before.

Gene expression responses to acute and chronic heat stress in the common reef-building coral Pocillopora verrucosa

Global climate change is impacting coral reefs worldwide, with approximately 19% of reefs being permanently degraded, 15% showing symptoms of imminent collapse, and 20% at risk of becoming critically affected in the next few decades. This alarming level of reef degradation is mainly due to an increase in frequency and intensity of natural and anthropogenic disturbances. Recent evidence has called into question whether corals have the capacity to acclimatize or adapt to climate changes and some groups of corals showed inherent physiological tolerance to environmental stressors. The aim of the present study was to evaluate mRNA expression patterns underlying differences in thermal tolerance in specimen of the common reef-building coral Pocillopora verrucosa collected at different locations in Bangka Island waters (North Sulawesi, Indonesia). Part of the experimental work was carried out at the CoralEye Reef Research Outpost (Bangka Island). This includes sampling of corals at selected sites and at different depths (3 and 12 m) as well as their experimental exposure to an increased water temperature under controlled conditions for 3 and 7 days. Levels of mRNAs encoding ATP synthase (ATPs) NADH dehydrogenase (NDH) and a 70kDa Heat Shock Protein (HSP70) were evaluated by quantitative real time PCR. Transcriptional profiles evaluated under field conditions suggested an adaptation to peculiar local environmental conditions in corals collected at different sites and at the low depth. Nevertheless, high–depth collected corals showed a less pronounced site-to-site separation suggesting more homogenous environmental conditions. Exposure to an elevated temperature under controlled conditions pointed out that corals adapted to the high depth are more sensitive to the effects of thermal stress, so that reacted to thermal challenge by significantly over-expressing the selected gene products. Being continuously exposed to fluctuating environmental conditions, low-depth adapted corals are more resilient to the stress stimulus, and indeed showed unaffected or down-regulated mRNA expression profiles. Overall these results highlight that transcriptional profiles of selected genes involved in cellular stress response are modulated by natural seasonal temperature changes in P. verrucosa. Moreover, specimens living in more variable habitats (low-depth) exhibit higher basal HSP70 mRNA levels, possibly enhancing physiological tolerance to environmental stressors.

Changes in CXCL12/CXCR4-chemokine expression during onset of colorectal malignancies

Chemokines have been proposed to contribute to tumour growth and metastatic spread of several cancer entities. Here, we examined the relative levels of CXCL12/CXCR4 in resection specimens from patients with different malignant and non-malignant colorectal diseases as well as colorectal liver metastases (CRLM). CXCL12/CXCR4 mRNA and protein expression profiles were assessed by quantitative real-time PCR, Western blot analysis, enzyme-linked immunosorbent assay (ELISA) and immunohistochemistry in resection specimens from patients with ulcerative colitis (UC; n = 15), colorectal adenoma (CRA; n = 15), colorectal adenocarcinoma (CRC; n = 47) and CRLM (n = 16). Corresponding non-affected tissues served as control. In contrast to UC tissues, CXCL12 showed a distinct down-regulation in CRA, CRC and CRLM specimens, whereas the corresponding receptor CXCR4 demonstrated a significant up-regulation in CRC and CRLM related to corresponding non-affected tissues (p < 0.05, respectively). Our results strongly suggest an association between CXCL12/CXCR4 expression and the induction of CRA, CRC and the development of CRLM. Therefore, CXCR4 may be a potential target for specific therapeutic interventions.

Structural differences between the bladder dome and trigone revealed by mRNA expression analysis of cold-cut biopsies

To establish the mRNA expression profiles of selected genes involved in bladder contractility and epithelial permeability in the bladder dome and trigone in order to evaluate the use of cold-cut biopsies for comparative quantitative studies into the anatomical differences between these two bladder regions.

Single-cell gene-expression profiling reveals qualitatively distinct CD8 T cells elicited by different gene-based vaccines

CD8 T cells play a key role in mediating protective immunity against selected pathogens after vaccination. Understanding the mechanism of this protection is dependent upon definition of the heterogeneity and complexity of cellular immune responses generated by different vaccines. Here, we identify previously unrecognized subsets of CD8 T cells based upon analysis of gene-expression patterns within single cells and show that they are differentially induced by different vaccines. Three prime-boost vector combinations encoding HIV Env stimulated antigen-specific CD8 T-cell populations of similar magnitude, phenotype, and functionality. Remarkably, however, analysis of single-cell gene-expression profiles enabled discrimination of a majority of central memory (CM) and effector memory (EM) CD8 T cells elicited by the three vaccines. Subsets of T cells could be defined based on their expression of Eomes, Cxcr3, and Ccr7, or Klrk1, Klrg1, and Ccr5 in CM and EM cells, respectively. Of CM cells elicited by DNA prime-recombinant adenoviral (rAd) boost vectors, 67% were Eomes(-) Ccr7(+) Cxcr3(-), in contrast to only 7% and 2% stimulated by rAd5-rAd5 or rAd-LCMV, respectively. Of EM cells elicited by DNA-rAd, 74% were Klrk1(-) Klrg1(-)Ccr5(-) compared with only 26% and 20% for rAd5-rAd5 or rAd5-LCMV. Definition by single-cell gene profiling of specific CM and EM CD8 T-cell subsets that are differentially induced by different gene-based vaccines will facilitate the design and evaluation of vaccines, as well as enable our understanding of mechanisms of protective immunity.

MicroRNA Expression Array Identifies Novel Diagnostic Markers for Conventional and Oncocytic Follicular Thyroid Carcinomas

Objective:The most difficult thyroid tumors to be diagnosed by cytology and histology are conventional follicular carcinomas (cFTCs) and oncocytic follicular carcinomas (oFTCs). Several microRNAs (miRNAs) have been previously found to be consistently deregulated in papillary thyroid carcinomas; however, very limited information is available for cFTC and oFTC. The aim of this study was to explore miRNA deregulation and find candidate miRNA markers for follicular carcinomas that can be used diagnostically.Design:Thirty-eight follicular thyroid carcinomas (21 cFTCs, 17 oFTCs) and 10 normal thyroid tissue samples were studied for expression of 381 miRNAs using human microarray assays. Expression of deregulated miRNAs was confirmed by individual RT-PCR assays in all samples. In addition, 11 follicular adenomas, two hyperplastic nodules (HNs), and 19 fine-needle aspiration samples were studied for expression of novel miRNA markers detected in this study.Results:The unsupervised hierarchical clustering analysis demonstrated individual clusters for cFTC and oFTC, indicating the difference in miRNA expression between these tumor types. Both cFTCs and oFTCs showed an up-regulation of miR-182/-183/-221/-222/-125a-3p and a down-regulation of miR-542-5p/-574-3p/-455/-199a. Novel miRNA (miR-885-5p) was found to be strongly up-regulated (>40-fold) in oFTCs but not in cFTCs, follicular adenomas, and HNs. The classification and regression tree algorithm applied to fine-needle aspiration samples demonstrated that three dysregulated miRNAs (miR-885-5p/-221/-574-3p) allowed distinguishing follicular thyroid carcinomas from benign HNs with high accuracy.Conclusions:In this study we demonstrate that different histopathological types of follicular thyroid carcinomas have distinct miRNA expression profiles. MiR-885-5p is highly up-regulated in oncocytic follicular carcinomas and may serve as a diagnostic marker for these tumors. A small set of deregulated miRNAs allows for an accurate discrimination between follicular carcinomas and hyperplastic nodules and can be used diagnostically in fine-needle aspiration biopsies.

Comparative expression profiling of E. coli and S. aureus inoculated primary mammary gland cells sampled from cows with different genetic predispositions for somatic cell score

BACKGROUND: During the past ten years many quantitative trait loci (QTL) affecting mastitis incidence and mastitis related traits like somatic cell score (SCS) were identified in cattle. However, little is known about the molecular architecture of QTL affecting mastitis susceptibility and the underlying physiological mechanisms and genes causing mastitis susceptibility. Here, a genome-wide expression analysis was conducted to analyze molecular mechanisms of mastitis susceptibility that are affected by a specific QTL for SCS on Bos taurus autosome 18 (BTA18). Thereby, some first insights were sought into the genetically determined mechanisms of mammary gland epithelial cells influencing the course of infection. METHODS: Primary bovine mammary gland epithelial cells (pbMEC) were sampled from the udder parenchyma of cows selected for high and low mastitis susceptibility by applying a marker-assisted selection strategy considering QTL and molecular marker information of a confirmed QTL for SCS in the telomeric region of BTA18. The cells were cultured and subsequently inoculated with heat-inactivated mastitis pathogens Escherichia coli and Staphylococcus aureus, respectively. After 1, 6 and 24 h, the cells were harvested and analyzed using the microarray expression chip technology to identify differences in mRNA expression profiles attributed to genetic predisposition, inoculation and cell culture. RESULTS: Comparative analysis of co-expression profiles clearly showed a faster and stronger response after pathogen challenge in pbMEC from less susceptible animals that inherited the favorable QTL allele 'Q' than in pbMEC from more susceptible animals that inherited the unfavorable QTL allele 'q'. Furthermore, the results highlighted RELB as a functional and positional candidate gene and related non-canonical Nf-kappaB signaling as a functional mechanism affected by the QTL. However, in both groups, inoculation resulted in up-regulation of genes associated with the Ingenuity pathways 'dendritic cell maturation' and 'acute phase response signaling', whereas cell culture affected biological processes involved in 'cellular development'. CONCLUSIONS: The results indicate that the complex expression profiling of pathogen challenged pbMEC sampled from cows inheriting alternative QTL alleles is suitable to study genetically determined molecular mechanisms of mastitis susceptibility in mammary epithelial cells in vitro and to highlight the most likely functional pathways and candidate genes underlying the QTL effect.

14-3-3- and II/3-10-gene expression as molecular markers to address viability and growth activity of Echinococcus multilocularis metacestodes

The present study aimed to search for and characterize parasite molecules, whose expression levels correlate with the viability and growth activity of Echinococcus multilocularis metacestodes. We focused on the expression profiles of 2 parasite-derived genes, 14-3-3 and II/3-10, as putative molecular markers for viability and growth activity of the larval parasite. In experiments in vivo, gene expression levels of 14-3-3 and II/3-10 were relatively quantified by real-time reverse transcription-PCR using a housekeeping gene, beta-actin, as a reference reaction. All three reactions were compared with growth activity of the parasite developing in permissive nu/nu and in non-permissive wild type BALB/c mice. At 2 months p.i., the transcription level of 14-3-3 was significantly higher in parasites actively proliferating in nu/nu mice compared to parasites moderately growing in wild type mice. Immunoblotting experiments confirmed at the protein level that 14-3-3 was over-expressed in parasites derived from nu/nu mice at 2 months p.i. In vitro treatment of E. multilocularis with an anti-echinococcal drug nitazoxanide resulted in a significant decrease of both 14-3-3 and II/3-10 transcription levels found after 8 days of treatment, which correlated with the kinetics of a housekeeping gene, beta-actin. The conclusion is that 14-3-3, combined with II/3-10, exhibits good potential as a molecular marker to assess viability and growth activity of the parasite.

Gene expression profiling reveals consistent differences between clinical samples of human leukaemias and their model cell lines

Microarray gene expression profiles of fresh clinical samples of chronic myeloid leukaemia in chronic phase, acute promyelocytic leukaemia and acute monocytic leukaemia were compared with profiles from cell lines representing the corresponding types of leukaemia (K562, NB4, HL60). In a hierarchical clustering analysis, all clinical samples clustered separately from the cell lines, regardless of leukaemic subtype. Gene ontology analysis showed that cell lines chiefly overexpressed genes related to macromolecular metabolism, whereas in clinical samples genes related to the immune response were abundantly expressed. These findings must be taken into consideration when conclusions from cell line-based studies are extrapolated to patients.

Gene expression in cortex and hippocampus during acute pneumococcal meningitis

BACKGROUND: Pneumococcal meningitis is associated with high mortality (approximately 30%) and morbidity. Up to 50% of survivors are affected by neurological sequelae due to a wide spectrum of brain injury mainly affecting the cortex and hippocampus. Despite this significant disease burden, the genetic program that regulates the host response leading to brain damage as a consequence of bacterial meningitis is largely unknown.We used an infant rat model of pneumococcal meningitis to assess gene expression profiles in cortex and hippocampus at 22 and 44 hours after infection and in controls at 22 h after mock-infection with saline. To analyze the biological significance of the data generated by Affymetrix DNA microarrays, a bioinformatics pipeline was used combining (i) a literature-profiling algorithm to cluster genes based on the vocabulary of abstracts indexed in MEDLINE (NCBI) and (ii) the self-organizing map (SOM), a clustering technique based on covariance in gene expression kinetics. RESULTS: Among 598 genes differentially regulated (change factor > or = 1.5; p < or = 0.05), 77% were automatically assigned to one of 11 functional groups with 94% accuracy. SOM disclosed six patterns of expression kinetics. Genes associated with growth control/neuroplasticity, signal transduction, cell death/survival, cytoskeleton, and immunity were generally upregulated. In contrast, genes related to neurotransmission and lipid metabolism were transiently downregulated on the whole. The majority of the genes associated with ionic homeostasis, neurotransmission, signal transduction and lipid metabolism were differentially regulated specifically in the hippocampus. Of the cell death/survival genes found to be continuously upregulated only in hippocampus, the majority are pro-apoptotic, while those continuously upregulated only in cortex are anti-apoptotic. CONCLUSION: Temporal and spatial analysis of gene expression in experimental pneumococcal meningitis identified potential targets for therapy.

Stem cell-related "self-renewal" signature and high epidermal growth factor receptor expression associated with resistance to concomitant chemoradiotherapy in glioblastoma

PURPOSE: Glioblastomas are notorious for resistance to therapy, which has been attributed to DNA-repair proficiency, a multitude of deregulated molecular pathways, and, more recently, to the particular biologic behavior of tumor stem-like cells. Here, we aimed to identify molecular profiles specific for treatment resistance to the current standard of care of concomitant chemoradiotherapy with the alkylating agent temozolomide. PATIENTS AND METHODS: Gene expression profiles of 80 glioblastomas were interrogated for associations with resistance to therapy. Patients were treated within clinical trials testing the addition of concomitant and adjuvant temozolomide to radiotherapy. RESULTS: An expression signature dominated by HOX genes, which comprises Prominin-1 (CD133), emerged as a predictor for poor survival in patients treated with concomitant chemoradiotherapy (n = 42; hazard ratio = 2.69; 95% CI, 1.38 to 5.26; P = .004). This association could be validated in an independent data set. Provocatively, the HOX cluster was reminiscent of a "self-renewal" signature (P = .008; Gene Set Enrichment Analysis) recently characterized in a mouse leukemia model. The HOX signature and EGFR expression were independent prognostic factors in multivariate analysis, adjusted for the O-6-methylguanine-DNA methyltransferase (MGMT) methylation status, a known predictive factor for benefit from temozolomide, and age. Better outcome was associated with gene clusters characterizing features of tumor-host interaction including tumor vascularization and cell adhesion, and innate immune response. CONCLUSION: This study provides first clinical evidence for the implication of a "glioma stem cell" or "self-renewal" phenotype in treatment resistance of glioblastoma. Biologic mechanisms identified here to be relevant for resistance will guide future targeted therapies and respective marker development for individualized treatment and patient selection.

Differential expression and localization of lipid transporters in the bovine mammary gland during the pregnancy-lactation cycle

The transport of lipids across mammary gland epithelial cells (MEC) determines milk lipid content and composition. We investigated the expression of lipid transporters and their regulators in comparison to blood metabolites during lactation and dry period (DP) in dairy cows. Repeated mammary gland biopsies and blood samples were taken from 10 animals at 7 stages of the pregnancy-lactation cycle. Expression levels of the specific mRNAs were determined by quantitative reverse transcription-PCR, whereas ABCA1 was localized by immunohistochemistry. Blood serum metabolites were determined by common enzymatic chemistries. Elevated mRNA profiles of ABCA1 and ABCA7 were found during DP as compared with lactation and were inversely associated with blood cholesterol levels. Elevated levels of ABCG2, NPC1, SREBP1, SREBP2, LXR alpha, and PPAR gamma were found postpartum, whereas ABCG1 did not differ between the functional stages of the mammary gland. The ABCA1 protein was localized in MEC and showed differential activity between DP and lactation suggesting a role of ABCA1 in the removal of excess cellular cholesterol from MEC during the DP. The expression profiles of ABCA7 and NPC1 may reflect a role of these transporters in the clearance of apoptotic cells and the intracellular redistribution of cholesterol, respectively. Regulation of lipid transporters in the mammary gland is partially associated with transcription factors that control lipid homeostasis.

Estrogen Modulates Hepatic Gene Expression and Survival of Rainbow Trout Infected with Pathogenic Bacteria Yersinia ruckeri

In the aquatic environment, fish are exposed to various stimuli at once and have developed different response mechanisms to deal with these multiple stimuli. The current study assessed the combined impacts of estrogens and bacterial infection on the physiological status of fish. Juvenile rainbow trout were exposed to two different concentrations of 17 beta-estradiol (E2) (2 or 20 mg/kg feed) and then infected with three concentrations of Yersinia ruckeri, a bacterial pathogen causing massive losses in wild and farmed salmonid populations. Organism-level endpoints to assess the impact of the single and combined treatments included hepatic vitellogenin transcript expression to evaluate the E2 exposure efficiency and survival rate of pathogen-challenged fish. The two E2 doses increased vitellogenin levels within the physiological range. Infection with Y. ruckeri caused mortality of trout, and this effect was significantly enhanced by a simultaneous exposure to high E2 dose. The hormone reduced survival at intermediate and high (10(4) and 10(6) colony forming units, cfu) bacterial concentrations, but not for a low one (10(2) cfu). Analysis of hepatic gene expression profiles by a salmonid 2 k cDNA microarray chip revealed complex regulations of pathways involved in immune responses, stress responses, and detoxicification pathways. E2 markedly reduced the expression of several genes implicated in xenobiotic metabolism. The results suggest that the interaction between pathogen and E2 interfered with the fish's capability of clearing toxic compounds. The findings of the current study add to our understanding of multiple exposure responses in fish.

Analysis of protein expression in zebrafish during gonad differentiation by targeted proteomics

The molecular mechanisms governing sex determination and differentiation in the zebrafish (Danio rerio) are not fully understood. To gain more insights into the function of specific genes in these complex processes, the expression of multiple candidates needs to be assessed, preferably on the protein level. Here, we developed a targeted proteomics method based on selected reaction monitoring (SRM) to study the candidate sex-related proteins in zebrafish which were selected based on a global proteomics analysis of adult gonads and representational difference analysis of male and female DNA, as well as on published information on zebrafish and other vertebrates. We employed the developed SRM protocols to acquire time-resolved protein expression profiles during the gonad differentiation period in vas::EGFP transgenic zebrafish. Evidence on protein expression was obtained for the first time for several candidate genes previously studied only on the mRNA level or suggested by bioinformatic predictions. Tuba1b (tubulin alpha 1b), initially included in the study as one of the potential housekeeping proteins, was found to be preferentially expressed in the adult testis with nearly absent expression in the ovary. The revealed changes in protein expression patterns associated with gonad differentiation suggest that several of the examined proteins, especially Ilf2 and Ilf3 (interleukin enhancer-binding factors 2 and 3), Raldh3 (retinaldehyde dehydrogenase type 3), Zgc:195027 (low density lipoprotein-related receptor protein 3) and Sept5a (septin 5a), may play a specific role in the sexual differentiation in zebrafish.