Alex C. King
Armin K. Zenker- University of Applied Sciences and Arts North-Western Switzerland (FHNW), Muttenz, Switzerland
The sex of both humans and Danio rerio has previously been shown to affect the way individuals respond to drug exposure. Genes which allow identification of sex in juvenile zebrafish show potential to reveal these confounding variables between sex in toxicological and preclinical trials but the link between these is so far missing. These sex-specific, early expressed genes where expression is not altered by drug exposure must be carefully selected for this purpose. We aimed to discover genes which can be used in pharmaceutical trials and environmental toxicology studies to uncover sex-related variations in gene expression with drug application using the model organism Danio rerio. Previously published early sex determining genes from King et al. were evaluated as well as additional genes selected from our zebrafish Next-generation sequencing (NGS) data which are known from previously published works not to be susceptible to changes in expression with drug exposure. NGS revealed a further ten female-specific genes (vtg1, cyp17a1, cyp19a1a, igf3, ftz-f1, gdf9, foxl2a, Nr0b1, ipo4, lhcgr) and five male related candidate genes (FKBP5, apobb1, hbaa1, dmrt1, spata6) which are also expressed in juvenile zebrafish, 28 days post fertilisation (dpf). Following this, a literature review was performed to classify which of these early-expressed sex specific genes are already known to be affected by drug exposure in order to determine candidate genes to be used in pharmaceutical trials or environmental toxicology testing studies. Discovery of these early sex-determining genes in Danio rerio will allow identification of sex-related responses to drug testing to improve sex-specific healthcare and the medical treatment of human patients.
1 Introduction
Today, Danio rerio (Zebrafish) are popular model organisms for use in many scientific fields of study. The characteristics of embryonic development in zebrafish including their rapid growth phase, high fecundity and transparency make them ideal as a model organism (Grunwald and Eisen, 2002; Segner, 2009). Zebrafish are frequently used as a pharmacological tool in drug discovery (Goldsmith, 2004; Zon and Peterson, 2005), in environmental risk assessment studies and to aid understanding of human form, function and mechanisms of disease (Nagabhushana and Mishra, 2016; Chen et al., 2017).
Adult, sexually mature zebrafish are used in the 21-day Fish Assay, a short-term screening for oestrogenic and androgenic activity, and aromatase inhibition (OECD, 2009). In these tests, sex can be easily determined both morphologically and genetically. However, during both Fish Embryo Acute Toxicity (FET) tests (OECD, 2013), which measure acute toxicity and multi-effect detecting Fish Early-Life Stage (FELS) toxicity test (OECD, 2010), the sex of the embryos is indistinguishable (Kimmel et al., 1995). The importance of knowing the sex of the individual being tested comes as sex-specific differences to drug exposure are revealed (Kling et al., 2008; Dlugos et al., 2011). Puzzling results from FET and FELS tests are often overlooked as being caused by differing effects of various drugs on sex (Holdcroft, 2007). Males and females react in a different way to exposure but so far FET and FELS tests fail to adequately link sex-related variations in gene expression to drug exposure (Liu et al., 2021; Fisher et al., 2022). Knowledge of sex at juvenile stages during these tests is fundamental to prevent false or unclear results being presented due to confounding sex variables which mask differences to drug application. By revealing sex during these early-stage tests we can discover the effect of the drugs in later life as well as the sensitivity of male and females to each drug. FELS and FET tests have several advantages when compared to the 21-day fish assay, they both use younger, smaller fish which means that fish are more susceptible to the effects of drugs (Vasconcelos et al., 2010; Mohammed, 2011), as well as being cheaper to use than adult fish. The 3Rs, reduce, refine, replace philosophy aims to decrease the use of animals in laboratory tests and minimise their pain by improving their welfare standards; this theory arguable points to use of juvenile fish over adults (Flecknell, 2002; Leconte et al., 2011).
In mammals and fish, there is considerable evidence that sex creates variations in susceptibility and response to exposure of drugs (Gartlehner et al., 2010). Use of a certain pharmaceuticals or exposure to specific chemicals results in changes to hormone levels as well as altering the regulation of specific sex differentiation genes (Kokane and Perrotti, 2020). Pharmaceuticals with low concentrations of active progesterone (P4) or norgestrel (NGT) as well as those containing substances which alter steroidogenic genes change the balance of sex hormones, namely, oestrogens and androgens, in juvenile zebrafish. These are thought to consequently affect the expression level of sex-determining genes (Sox9a, gapdhs, atp1b1a, cyp26b1, gyg1a, rdh10b, pdia, KPNA2, ccnb1, ctsla, Chr4, bmp15 and zbp3) (King et al., 2020); therefore, disrupting sex differentiation in zebrafish.
The action of some drugs (anti-inflammatories) is driven by the oestrogen pathway, thereby making females more susceptible to these drugs than males (Houben et al., 2020; Kokane and Perrotti, 2020). Opposingly oestrogen can have an inhibitory effect on the way a female liver breaks down medication (Osman, 2003; Martín-Millán and Castañeda, 2013). Studies have previously shown that developmental exposure of zebrafish to 17-Alpha-ethinyloestradiol (EE2) increases adult anxiety and shoal cohesion as well as changing gene expression, these outcomes were found to differ significantly with sex (Volkova et al., 2015). Other studies have also shown that application of P4 results in a female sex-biased population; whereas NGT results in a male dominated population (Liang et al., 2015). P4 is mostly found in human and animal excretion and is an essential regulator for the growth and maturation of oocytes. Synthetic progestins are often active ingredients of human pharmaceuticals, including oral and implant contraceptives (Kejuan et al., 2007).
Awareness of the differing responses that drug exposure has on male and females is important and can be life changing. Pharmacokinetics and pharmacodynamics as well as differences in genetic make-up, gene expression and hormone regulation of drug metabolising enzymes and transporters, body weight, height, body surface area and total water body content are significantly different in men and women (Soldin and Mattison, 2009; Plank-Bazinet et al., 2016). These factors create sex-related variations in the absorption, distribution, metabolism, and excretion of drugs, which in turn affects the drug efficacy and creates variations in adverse reactions. This points to evidence suggesting that personalised medicine to cater for differences between males and female would be beneficial for improving reactions to drug response and for understanding sex-based disease susceptibilities (U.S. Government Accountability Office, 2001; Nicolson et al., 2010; Lyubimov et al., 2012; Ridker et al., 2015). Previous trials with mice show that the dosage of sodium amytal for males can be half that which is effective to anesthetise females (Yang and Li, 2012). CYP3A4-substrate drugs, including, cyclosporine, erythromycin, tirilazad, verapamil, nifedipine, diazepam and alfentanil, take longer to clear in women than men, even after taking into consideration sex-specific physiological factors (Franconi et al., 2007; Greenblatt and Von Moltke, 2008). Additionally, volunteers for drug tests are normally male, non-smokers aged between 25–35. Therefore, there is a gross underestimation of female susceptibility to drugs because of this (Kim et al., 2010; McCarthy et al., 2017).
The main purpose of this study is to identify whether sex specific genes which are expressed during juvenile stages of zebrafish development are known to be affected by drug exposure. For this purpose, we used transcript data from NGS of zebrafish to select additional appropriate novel, early expressed, sex-specific genetic markers, thereby, further unravelling the puzzle of sex-differentiating genes which are expressed during juvenile development. The carefully selected juvenile-expressed genes shown in this paper were categorised based on a literature review of results from previously published studies in order to indicate any sex-related differences in expression of these genes due to interference from drug exposure. Our research will aid future zebrafish studies as well as helping to uncover the confounding variable of sex during drug testing. This will allow underlying, masked sex-related differences in response to pharmaceuticals to be highlighted. In turn, an advancement in sex-specific drug development can be achieved by providing knowledge regarding sex variations in which to build and develop appropriate medication including providing safe doses of drugs for both males and females.
2 Experimental procedures
2.1 RNA sequencing using NGS
NGS RNA sequences were produced from two 28-day-old zebrafish (which we classified as juvenile in this paper), two adult male and two adult female zebrafish using whole-body tissue samples. Juveniles were tested to indicate whether sex was genetically visible at this stage of development. RNA quality testing and NGS were performed at the Biocentre under the supervision of Philippe Demougin (Life Sciences Training Facility (LSTF), Basel), according to Illumina’s standard protocol for RNA sequencing (Illumina Inc., San Diego, USA, Cat. # RS-100-0801). The NGS transcript data of the sampled zebrafish were compared to identify possible genes that contribute to sex determination as well as those which are expressed early in juvenile development, 28dpf (King et al., 2020). All data analyses steps were carried out using the Galaxy server (usegalaxy.org) (Jalili et al., 2020). Data from transcript sequencing was mapped to the GRCz11 reference genome (GCA_000002035.4) using HISAT2 (Galaxy Version 2.1.0+galaxy4) (Kim et al., 2019). Transcripts were assembled and merged from the mapped reads with StringTie (Galaxy Version 1.3.4) (Pertea et al., 2015). After determining gene expression with featureCount (Galaxy Version 1.6.3+galaxy2), the NGS data was normalised and differential gene expression analysis of juvenile, male and female zebrafish was carried out using DESeq2 (Galaxy Version 2.11.40.2) (Love et al., 2014; King et al., 2020). Normalisation was based on the ‘median of ratios’ - the counts divided by sample-specific size factors determined by median ratio of gene counts relative to geometric mean per gene, to account for sequencing depth and RNA composition. This allowed us to produce a table of gene expression which highlighted genes with the highest transcript count difference between males and females. From this data, genes expressed at juvenile-stage (28dpf) which were also expressed highly in male or female samples were selected as the best gene candidates for indicating early sex identity of zebrafish. Four male (Sox9a, Gapdhs, atp1b1a and cyp26b1) and nine female genes (gyg1a, rdh10b, pdia, KPNA2, ccnb1, ctsla, Chr4, bmp15 and zbp3) which were previously identified as an early sex-determining marker in zebrafish (King et al., 2020) were used within a literature search to access previous information as to whether their level of expression was known to be influenced by certain drugs. In this study a further ten female and five male genes, which could be used in early sex-determination were selected from NGS data. A literature search was also carried out on these genes to discover if they were known to be affected by drug usage. Based on published literature from ScienceDirect in early 2021, we obtained information previous studies from the past 20 years by searching for the key terms; including the gene name, drug or pharmaceutical exposure, sex, or gender influence. Then, we eliminated irrelevant literature by reading the titles and abstracts and supplemented our literature database by reading the references of the selected studies. We used the selected literature to inform whether genes highlighted from our NGS data as early expressed within juveniles could be candidates in drug trials. Finally, we grouped early expressed genes from NGS data based on the literature search into three categories. If the study was shown to affect the expression of the selected genes, then they were deemed not good to use in pharmaceutical trials. Conversely, if it was not known from previous research that drugs affect the expression of a certain gene, it was included as a potential candidate for revealing how sex response differs with drug exposure and further investigation of its use in this area is necessary in future research. Statistics were not performed due to the small sampling size but gene expression tables were produced to indicate genes with high variation between males and females.
2.2 Statement on experimental protocols
NGS data was collected from 28-day-old and adult zebrafish, which had already been euthanised with Tricaine by the breeder. No experiments were carried out with live adult nor subadult fish before euthanasia. Accordingly, animal welfare guidelines were followed while zebrafish were kept by the breeder.
3 Results
Five male and ten female genes which show great potential as novel early sex-determining gene markers were highlighted from zebrafish NGS data (Figures 1, 2, respectively). The 15 genes were put into three categories based on the results of previous published studies (Table 1). The three categories; category one: gene expression influenced by drugs, category two: drug exposure thought not to influence gene expression (sometimes it is unclear whether gene expression is influenced by drugs and more research is needed to clarify this) and category three: gene expression not influenced by drugs, were used to highlight if drugs have been previously found to affect gene expression of the selected early sex-determining genes.
FIGURE 1. Heat map showing the relative transcript expression from NGS data of four previously identified male genes in bold letters (King et al., 2020) and five novel male genes, which can be used for early sex-determination in zebrafish (juveniles age 28dpf). Underscored no drug influence.
TABLE 1. Summary of early-expressed sex-determining genes selected from a previous study and NGS data in this study and their relative category to indicate if their expression is influenced by drug exposure based on a literature review (category one: gene expression influenced by drugs, category two: drug exposure thought not to influence gene expression and category three: gene expression not influenced by drugs).
The expression of a further four male (Sox9a, Gapdhs, atp1b1a and cyp26b1) and nine female genes (gyg1a, rdh10b, pdia, KPNA2, ccnb1, ctsla, Chr4, bmp15 and zbp3), which were previously identified as early-sex determining genes (King et al., 2020), were also compared to drug exposure and categorised into the three above influence categories (Figures 1, 2 respectively).
3.1 Genes selected for early male sex-determination
Five novel genes were selected from zebrafish NGS data; three of the five genes were among the list of transcripts with the highest differentiation between the males and females. The other two genes were selected from our NGS data because of highlighted sex differences in their expression. The following five sex-related genes were selected as potential markers of male sex-determination in juvenile zebrafish, 28dpf (Figure 1). These identifications are all supported by previous literature.
3.1.1 FKBP prolyl isomerase 5 (fkbp5)
Based on NGS data, fkbp5 was expressed 2.8 times higher in males compared to females. Fkbp5 was also expressed in juveniles, 28dpf. The gene, fkbp5, was highlighted in the list of the top 13 male genes with the highest transcript differences between male and female zebrafish. The gene fkbp5 is involved in FK506 and heat shock protein binding, isomerase, peptidyl-prolyl cis-trans isomerase and peptidyl-prolyl cis-trans isomerase activity (Zfn.org—ZFIN ID: ZDB-GENE-030616-630).
3.1.2 Apolipoprotein Bb, tandem duplicate 1 (apobb1)
Using NGS data, apobb1 was selected from the list of genes which showed the highest difference between male and female transcript data. There was also high expression of apobb1 in juvenile samples. On average adult male zebrafish expressed apobb1 2.1 times more than females based on transcript data. Apobb1 is predicted to have lipid transporter activity and is involved in cellular response to xenobiotic stimulus (Zfn.org—ZFIN ID: ZDB-GENE-030131-9732).
3.1.3 Haemoglobin, alpha adult 1 (hbaa1)
Hbaa1 was selected from the table of gene expression which highlighted genes with the highest transcript count difference between males and females. Here, a high expression of hbaa1 was found in juveniles and adult male samples–transcript data of these was almost identical. Expression of hbaa1 was 1.7-fold higher in males compared to females based on NGS data. Hbaa1 is predicted to contribute to haptoglobin binding and peroxidase activity and is involved in response to hypoxia (Zfn.org—ZFIN ID: ZDB-GENE-980526-79).
3.1.4 Double-sex and mab-3 related transcription factor 1 (dmrt1)
There was on average a four-fold higher expression of dmrt1 found in adult males compared with females, shown in NGS transcript data. Dmrt1 was also expressed weakly in juvenile NGS transcript data. It is predicted to have DNA-binding transcription factor activity; metal ion and sequence-specific DNA binding activity (Zfn.org—ZFIN ID: ZDB-GENE-050511-1).
3.1.5 Spermatogenesis associated 6 (spata6)
NGS transcript data revealed an average of two times higher expression of spata6 in male compared to female samples. Expression of spata6 was also seen in juvenile samples of NGS data. Spata6 is predicted to have myosin light chain binding activity and to be involved in spermatogenesis (Zfn.org—ZFIN ID: ZDB-GENE-050417-351).
3.2 Genes selected for early female sex-determination
The following ten female genes were selected from our NGS data and with support from previous literature were highlighted as potential markers of female sex-determination in juvenile zebrafish, 28dpf (Figure 2).
FIGURE 2. Heat map showing the relative transcript expression from NGS data of nine previously identified female genes in bold letters (King et al., 2020) and ten novel female genes, which can be used for early sex-determination in zebrafish (juveniles age 28dpf). Underscored no drug influence.
3.2.1 Vitellogenin 1 (vtg1)
We found from NGS data that there was a significantly higher female transcript expression of vtg1 compared with male zebrafish with on average 60-fold more transcripts for females compared to males. Based on NGS data, the expression of vtg1 is absent in both males and juveniles. Vtgs are female-specific, hepatically synthesised lipoproteins. They are the main carriers of retinoids and carotenoids (Zfn.org—ZFIN ID: ZDB-GENE-001201-1).
3.2.2 Cytochrome P450, family 17, subfamily A, polypeptide 1 (cyp17a1)
Our data indicates a 36-fold increase in cyp17a in female zebrafish compared to males. Slight expression of cyp17a1 was also found in juvenile samples. Cyp17a1 is expressed in Leydig cells; brain; and gonad. Cyp17a1 exhibits 17-alpha, 20-alpha-dihydroxypregn-4-en-3-one dehydrogenase; 17-alpha-hydroxyprogesterone aldolase; and steroid 17-alpha-monooxygenase activity (Zfn.org—ZFIN ID: ZDB-GENE-040213-2).
3.2.3 Cytochrome P450, family 19, subfamily A, polypeptide 1a (cyp19a1a)
Cyp19a1a was expressed on average 26 times higher in female compared to male samples based on our NGS transcript data. There is a lack of cyp19a1a in both male and juvenile samples based on NGS data. Cyp19a1a is predicted to have iron and haem-binding, as well as oxidoreductase activity. It is involved in several processes, including primary sex determination, germ-line fate; regulation of oestrogen biosynthetic process; and response to oestradiol (Zfn.org—ZFIN ID: ZDB-GENE-990415-43).
3.2.4 Insulin-like growth factor 3 (igf3)
Using NGS data from zebrafish, we found that igf3 expression in females was on average 44 times that of males. Igf3 was weakly expressed in juvenile samples. Igf3 exhibits insulin-like growth factor receptor binding activity. It is involved in oocyte maturation and sperm capacitation and is localised to the extracellular space (Zfn.org—ZFIN ID: ZDB-GENE-080611-1).
3.2.5 Nuclear receptor subfamily 5, group A, member 2 (ftz-f1, nr5a2 or ff1a)
Ftz-f1 gene was on average 1.6-fold higher expressed in females compared to males, based on NGS transcript data. Ftz-f1 was also expressed in juvenile zebrafish samples. The ftz-f1 gene exhibits DNA binding and DNA-binding transcription factor activity and is involved in several processes, including cartilage and digestive system development; and hepatoblast differentiation (Zfn.org—ZFIN ID: ZDB-GENE-990415-79).
3.2.6 Growth differentiation factor 9 (gdf9)
NGS data from this study revealed, 240 times more gdf9 transcripts when adult zebrafish females and males were compared. Gdf9 was only slightly expressed in juvenile zebrafish based on transcript data. Gdf9 is predicted to have cytokine and transforming growth factor beta receptor binding activity (Zfn.org—ZFIN ID: ZDB-GENE-050221-7).
3.2.7 Forkhead box L2a (foxl2a)
Foxl2a is on average 1.8 times higher expressed in female than in male samples. It was also expressed in juveniles, 28dpf. Foxl2a is predicted to have DNA-binding transcription factor and RNA polymerase II-specific activity. It is involved in maintenance of animal organ identity (Zfn.org—ZFIN ID: ZDB-GENE-060512-241).
3.2.8 Nuclear receptor subfamily 0, group B, member 1 (nr0b1 or dax1)
Nr0b1 was on average 1.9 times more expressed in female than male zebrafish based on transcript data from this study. Expression of nr0b1 was also seen in juvenile transcript data. It is predicted to have transcription co-repressor and transcription factor binding activity and is involved in negative regulation of male gonad development and water homeostasis (Zfn.org—ZFIN ID: ZDB-GENE-070130-1).
3.2.9 Importin 4 (ipo4)
Transcript data from NGS results indicated that ipo4 is on average expressed 4.2 times more in adult females when compared with adult males. NGS data in this study also showed expression of ipo4 in juvenile samples. Ipo4 is predicted to have nuclear import signal receptor and nuclear localisation sequence binding activity (Zfn.org—ZFIN ID: ZDB-GENE-041014-307).
3.2.10 Luteinizing hormone/choriogonadotropin receptor (lhcgr)
Lhcgr was 5.1 times more expressed female zebrafish compared to males based on NGS transcript data. Based on our data, lhcgr was also shown to be expressed in juveniles, 28dpf. Gonadotropin receptor lhcgr exhibits luteinising hormone receptor activity and is involved in gamete generation and gonad development (Zfn.org—ZFIN ID: ZDB-GENE-040806-3).
4 Discussion
We assume from the gene expression patterns in NGS data that both juveniles are male (King et al., 2020), indicating that the sex was already genetically visible at this stage of development (28dpf). Although our data shows slight differences in transcript number across sex categories for adult danio, the transcript number across each sex is within the same range. The fact that both the juveniles are males explains why some of the highly expressed selected female genes showed no or very little expression in juveniles, namely, vtg1, cyp19a1a, igf3 and gdf9 (Figure 2). Despite the lack of expression in our NGS data, the expression of these four genes in juvenile and female zebrafish have previously been extensively studied and it is a well-established that these are early expressed female-determining genes. Hence why they were selected as early-sex determining markers for their comparison against drug exposure in this study. Studies referenced in this paper which were used to categorise gene expression change influenced by drug exposure include a range of model organisms. Danio rerio and humans as well as other species have orthologous genes which means these organisms can be used for gene expression comparison (Moyer et al., 2021).
The process of NGS is expensive (Schwarze et al., 2020) but it produces a huge amount of data, thus two samples per group were deemed sufficient in this study.
4.1 Category one: drug exposure influences expression of certain early sex-determining genes in Danio rerio
It is important to determine if the expression of previously discovered (King et al., 2020) and novel sex-determining genes is known to be affected by drug application. Certain drugs, including hormone or steroid-based drugs, are known to commonly affect the expression of female and male-specific genes. In mammals, oestrogen is the major female sex hormone that regulates secondary sex differentiation. Previous studies have shown that pharmacological manipulation of oestrogen production in zebrafish alters sex determination (Guiguen et al., 2010). It is not applicable to use genes which are affected by hormone levels as markers for revealing confounding sex differences in pharmaceutical trials.
4.1.1 Formerly selected male sex-determining genes influenced by drug activity
4.1.1.1 Sox9a
King et al. (2020) Sox9a is an early expressed, male sex-determining gene which plays a role in spermatogenesis and is predominantly expressed in the testis (Groh et al., 2011; Chen et al., 2017; Lin et al., 2017). A high expression of sox9a has been found in both zebrafish in males, juveniles, 28dpf and during earlier development (Jørgensen et al., 2008; King et al., 2020). Despite the fact that sox9a has strong potential as a candidate gene for early sex determination, previous studies involving long-term exposure to natural and synthetic P4 and norgestrel and other drugs affecting steroid levels have been shown to affect sox9a expression, thereby altering sex-determination of zebrafish. Tamoxifen is a selective oestrogen receptor modulator used to prevent and treat breast cancer. Exposure to tamoxifen during development induces abnormal upregulation of sox9 in sertoli cells of female mice which in turn causes testis differentiation (Yu et al., 2014). Therefore, its use as a sex-determining gene in pharmaceutical studies cannot be justified (van boxtel et al., 2010; Liang et al., 2015).
4.1.1.2 Cyp26b1
Cyp26b1 is expressed highly in male and juvenile zebrafish (King et al., 2020). Retinoic acid (RA) is essential for spermatogenesis in male zebrafish. The relative level of RA in the gonads is correlated with sex-differentiation. Cyp26b1 degrades RA and has a higher expression in the testis. Cyp26b1 is upregulated in individuals of male fate. The loss of RA prevents developing gonads from entering into meiosis to form the ovaries, the opposite is true for female fated individuals (Rodríguez-Marí et al., 2013; Pradhan and Olsson, 2015). In mice, male transcript factors, sf1 and sox9 upregulate and activate cyp26b1. In gonads fated to become ovaries, foxl2 is highly expressed and acts an antagonist to cyp26b1 (Kashimada et al., 2011).
Retinoids have been well studied for their chemo-preventive properties. Retinoids are activated by their conversion to RA (Sedjo et al., 2004). In all vertebrates and some invertebrates, RA is a potent teratogen that, at pharmacological concentrations, can induce defects in congenital processes. RA signalling can be blocked by pharmaceuticals, thereby altering cyp26b1 level and sexual fate (Mark et al., 2006).
Anti-convulsants are a broad group of pharmacological agents which stabilise the level of nerve cell impulses and are used to treat epileptic seizures, bipolar disorder as well as for the treatment of neuropathic pain. Retinoids and vitamin A control have a vital role in sex-determination while the embryo is developing. Pharmaceuticals with anti-convulsant agents can have a huge effect on the retinoid metabolism. Due to the effect of cyp26b1 on RA expression during juvenile development, the altered retinoid metabolism by these drugs may significantly affect sex-determination (Nau et al., 1995). These results taken together indicate that cyp26b1 is not suitable to indicate sex variations during pharmaceutical trials.
4.1.1.3 Atp1b1a
Atp1b1a is highly expressed in adult male, juvenile and developing zebrafish (Wang et al., 2008b; Hatzold et al., 2016; King et al., 2020). The oestrogen medication, Ethinylestradiol (EE) is commonly used in combination with progestins in birth control pills. EE has often been used to treat menopausal symptoms, gynaecological disorders, and certain hormone-sensitive cancers. Sewage effluent sometimes contains concentrations of EE2 which alters reproduction in fish. After zebrafish were exposed to EE2, atp1b1a was downregulated (Martyniuk et al., 2007).
In recent times, graphene oxide (GO) has been researched for application of drug delivery in order to discover its properties for multi-pharmaceutical action including those with photothermal therapy. GO has been found to be useful in cancer therapies. In vitro studies have shown that GO is a valuable nanomaterial for stem cell research in medicine as it has the potential to promote stem cell adhesion, growth and differentiation (Duran et al., 2015; Priyadarsini et al., 2018). After zebrafish were exposed to 100 μg/L GO, the expression of atp1b1a was vastly altered. GO exposure to zebrafish caused atp1b1a to be downregulated. With these points considered, atp1a1b is not suitable as a sex-determining marker due to the alteration in expression of this gene with application of oestrogen-based drugs or those containing GO (Zou et al., 2018).
4.1.2 Novel male sex-determining genes influenced by drug activity
4.1.2.1 Fkbp5
Fkbp5 gene is an important target of androgen signalling in the prostate. Androgen-deprivation leads to a decrease in FKBP51 and androgen stimulation resulting in up to 5-fold upregulation of FKBP5 (Amler et al., 2000). Therefore, FKBP5 is termed androgen-responsive, is conserved among species and plays a role in zebrafish male development (Magee et al., 2006). The transcript data shown in this study indicates juvenile expression of fkbp5 and an almost a 3-fold increase of fkbp5 gene expression in adult male zebrafish compared to females. This is in correlation with previous studies; therefore, it can be concluded that fkbp5 can be used as an early indicator of sex in zebrafish.
Fkbp5 and fkbp51 expression has been linked to stress response regulation and can be used in anti-depressant treatment. Fkbp5, the gene expressing the protein fkbp51 (FK506-binding protein 51), is a strong inhibitor of the glucocorticoid receptor (GR). Restoration of GR signalling has been found to be a treatment for depression. In a previous study, patients who responded to anti-depressant treatment had a significant reduction of fkbp5 gene and associated fkbp51 protein expression, whereas an increasing expression level of fkbp5 was observed in non-responders (Ising et al., 2019). Therefore, the use of fkbp5 as a gene for identifying sex related differences in drug trials cannot be justified.
4.1.2.2 Hbaa1
In zebrafish, hbaa1 is expressed 19 organs including, eye; hematopoietic system; presumptive blood; and vasculature; with highest expression level in the spleen and 13th highest expression in the testis of zebrafish at fully formed stage. Zebrafish hbaa1 is orthologous to several human genes including HBZ (haemoglobin subunit zeta) (Zfn.org—ZFIN ID: ZDB-GENE-980526-79). Hbaa1 is expressed during zebrafish early development from zygote-1-cell-stage to adult-stage (Higano et al., 1997). Highest expression of hbaa1 in early development occurs at 35dpf (Tiedke et al., 2011). There are very few studies to our knowledge in which adult male and female zebrafish hbaa1 expression are compared. Our data indicates a higher expression of hbaa1 in adult male zebrafish compared with females as well as a high expression level in juveniles, therefore, hbaa1 should be considered as an early marker for male sex-determination in zebrafish but more research should be done to confirm this.
Again, few studies were found which show the effect of drugs on the expression of hbaa1 and its orthologs. Valproate (VPA) and its valproic acid, sodium valproate, and valproate semi-sodium forms are medications primarily used to treat epilepsy and bipolar disorder and prevent migraine headaches. A study found that VPA blocked the expression of human HBZ (Belrose et al., 2011; Gazon et al., 2016). More studies should now be carried out to investigate this pattern of expression in hbaa1 and its orthologs to discover for certain whether it can be discounted as a gene for uncovering sex-related differences during drug trials.
4.1.3 Formerly selected female sex-determining genes influenced by drug activity
4.1.3.1 Ccnb1
King et al. (2020) Ccnb1 has shown its potential to be an early expressed gene to indicate female sex (Yasuda et al., 2010; Horie and Kotani, 2016; King et al., 2020). The gene is expressed in the unfertilised egg, oocyte and mature ovarian follicle and involved in maturation and female sex determination in zebrafish (Kondo et al., 1997; Kondo et al., 2001; Knoll-Gellida et al., 2006; Wang et al., 2007; Nagahama and Yamashita, 2008; Kotani et al., 2013; Yasuda et al., 2013; Takahashi et al., 2014; Dingare et al., 2018; Takei et al., 2018; Yi et al., 2019).
Ccnb1 is significantly overexpressed in various cancer types including breast and prostate. A study showed that several tested drugs had an influence the expression of ccnb1; Doxorubicin, Cisplatin, Etoposide, Fluorouracil, Irinotecan, Ochratoxin A, Oxaliplatin and Tretinoin could decrease the expression of ccnb1, while Paclitaxel and Tamoxifen could increase the expression of ccnb1 (Ding et al., 2014). Additionally, resveratrol, a plant compound that acts like an antioxidant, has been used as an inhibitor of ccnb1 and can therefore regulate both the proliferation and apoptosis of pituitary tumour cells as well as altering the expression level of various EMT (epithelial-to-mesenchymal transition) markers (Li et al., 2019a). All things considered, the use of ccnb1 as an early sex determining gene in drug trials is not appropriate.
4.1.3.2 Rdh10b
Rdh10b is more highly expressed in female and juveniles than in males (King et al., 2020). The gene is necessary for oxidisation of retinol to retinal and RA (Sandell et al., 2007). Retinoids are found in eggs and oocytes of marine and freshwater fish, teleosts, amphibians and mammals (Costaridis et al., 1996; Irie and Seki, 2002). RA is highly expressed in the gonad of female zebrafish and is necessary for sexual development in adult female but not adult male mice (Koubova et al., 2006; Belyaeva et al., 2020b). A decrease in RA results in a male phenotype (Rodríguez-Marí et al., 2013). In the ovary of zebrafish, genes which are necessary for the synthesis of retinal are upregulated (Levi et al., 2012). Insulin represses rdh10 expression by inactivating FoxO1 (Obrochta et al., 2015). The alteration rdh10b gene expression with insulin shows that it is not suitable as a gene to reveal the confounding variable of sex in FELS tests.
4.1.3.3 Pdia
Pdia is expressed during zebrafish development. It is found highly expressed in the female gonad at juvenile and adult stages and been suggested as an early female sex-determining gene in zebrafish (Zheng et al., 2018; King et al., 2020). It has been observed that the expression of pdia decreases in the presence of tunicamycin, a mixture of homologous nucleoside antibiotics (Nair et al., 2007). Additionally, application of dithiothreitol (DTT), a reducing agent that causes the formation of unfolded proteins, causes the upregulation of pdia2 (Al-Sheikh et al., 2004). Due to the changing expression of pdia with drug influence its use as an early sex-determining gene for the discovery of drug related differences cannot be applied.
4.1.3.4 Ctsla
Ctsla was previous selected as a female sex-determining gene (King et al., 2020). This maternally inherited gene is known to be involved in oogenesis and is highly expressed in juvenile zebrafish and the ovary and oocytes of adult zebrafish (Tingaud-Sequeira and Cerdà, 2007; Tingaud-Sequeira et al., 2011; King et al., 2020).
The mammalian ortholog of ctsla is ctsl. Nifedipine, a calcium antagonist, was shown to suppress ctsl activity and mRNA expression in mice (Nishimura et al., 2002). Previous studies have demonstrated that the expression of ctsl was upregulated in several cancers; therefore, supporting its use as a therapeutic target for cancer (Rudzińska et al., 2019). There is evidence that ctsl may be central in drug resistance development (Zheng et al., 2004). Further research has shown that by inhibiting ctsl, development of resistance to cancer drugs such as doxorubicin, a chemotherapy medication used to treat cancer, can be reversed, and prevented. No matter what level of resistance cancer cells have reached, their reaction to drugs is restored after ctsl inhibition (Zheng et al., 2009). Additionally, some anti-tuberculous and anti-leprotic drugs inhibit the lysosomal cysteine proteinases such as ctsl; therefore, changing the expression of the gene (Kamboj et al., 2003). These results taken together imply that ctsla is not suitable to indicate sex variations during drug trials.
4.1.4 Novel female sex-determining genes influenced by drug activity
4.1.4.1 Vtg1
Vtg1 exhibits antioxidant activity and is involved in cellular response to oestrogen, oestradiol, and xenobiotic stimulus. It is expressed in several structures, including the unfertilised egg and at various stages through early development; in zygote:1-cell-stage at 0hpf to 0.75hpf, at larval-stage, 5dpf and all the way to breeding adult-stage at 90dpf to 730dpf (Zfn.org—ZFIN ID: ZDB-GENE-001201-1). During vitellogenesis, vtg1 is secreted and transported to the ovaries in the plasma. NGS data in this study as well as previous research showed that vtg1 has a higher overall expression in females compared to males (Levi et al., 2012). Increased expression of vtg1 in females is essential for the maintenance of the female sex and the prevention of testis differentiation in zebrafish.
Observations from previous studies indicate that oestrogen, E2, stimulates gene expression of vtg1 (King et al., 2020). When juvenile zebrafish are exposed to the oestrogen 17alpha-ethynyl oestradiol, the sex ratio is significantly female-biased; in one case an entirely female population was produced after exposure (Örn et al., 2003). Therefore, the use of vtg1 would be unsuitable as a marker for identifying female zebrafish in FET with endocrine disrupting substances, since the gene is upregulated by oestrogenic substances (Leet et al., 2011). However, the evidence of juvenile expression of vtg1 from previous literature as well as the significant female transcript expression from NGS data displayed in this study, indicates the usefulness of vtg1 as an early precursor for sex-determination. For these reasons, vtg1 can be used as a molecular marker for the indication of feminisation but not for uncovering sex-related differences in drug discovery (Jin et al., 2009).
4.1.4.2 Cyp17a1
Cyp17a1 is involved in androst-4-ene-3,17-dione biosynthetic and p4 metabolic process. Human orthologs of this gene are implicated in several diseases, including prostate cancer. In the zygote of zebrafish, cyp17a1 is already expressed during the 1-cell stage. Additionally, its expression is seen in the zygote at larval-stage, 5dpf; to breeding adult, 90dpf to 730dpf (Zfn.org—ZFIN ID: ZDB-GENE-040213-2). Elevated expression of cyp17a1 in female zebrafish can be seen in previous research; whereas cyp17a1a deficient zebrafish had male phenotypes (Zhai et al., 2018). A significantly higher expression of cyp17a1 can be seen in adult female zebrafish compared to males based on transcript data in this study and this is supported by previous research (Lin et al., 2017; Yan et al., 2017). These observations confirm that cyp17a1 is essential for female sexual development. It can be concluded from these findings that cyp17a1 could be suitable as a marker of female sex-specific expression in juvenile zebrafish.
Cyp17a1 has been found to be highly expressed in half about of human prostate carcinomas. Intracellular androgen is synthesised by these cancer cells. Abiraterone, a drug used to treat prostate cancer, can block nuclear accumulation of androgen receptors (ARs) which consequently suppresses expression of cyp17a1 (Giatromanolaki et al., 2019). The alteration of cyp17a1 gene expression with drug application means that it should not be selected to reveal sex-related differences with drug usage.
4.1.4.3 Igf3
In zebrafish, Igf3 is expressed in several structures, including gill; gonad; head; kidney; and skeletal muscle; with fifth highest expression level in the female gonad at fully-formed stage (Zfn.org—ZFIN ID: ZDB-GENE-080611-1). Igf3 is expressed in early ovary development (Wang et al., 2008a; McMenamin et al., 2013) and its expression increases throughout the course of ovary differentiation (Li et al., 2011). Igf3 has been found to be expressed from zygote 1-cell-stage, 0.0 h–0.75h, to adult, 90days to 730 days (Zfn.org—ZFIN ID: ZDB-GENE-080611-1). Igf3 is generally thought of as a female biasing gene (Xia et al., 2018). The increased expression of igf3 in adult female zebrafish can be explained by its role in oocyte maturation (Li et al., 2011). Using transcript data from this study, an increased expression of igf3 can be seen in adult female zebrafish compared to males. Since the whole zebrafish was used for the analysis of gene expression and no tissue-specific expression analysis of the ovaries was performed, the eggs present in the females in our study possibly contributed to the increased expression of igf3. Despite this, analysis of igf3 reveals it as a potentially suitable gene for assigning a sex-specific expression profile in developing zebrafish.
Insulin-like growth factor 1 receptor (IGF-1R) has already been used as a therapeutic target (Dolgin, 2020). Studies have shown that the stable overexpression of igfbp3 could suppress the growth of tumour cells by inducing apoptosis and suppressing survival and growth of cells (Dar et al., 2010). Despite its use as an early genetic marker of sex-determination, due to the change in the expression of igf3 with drug application it cannot be used as a gene to reveal sex-related differences during drug trials.
4.1.4.4 Nr0b1
Nr0b1 is predicted to have transcription corepressor activity and transcription factor binding activity. It is involved in negative regulation of male gonad development and water homeostasis. Nr0b1 is predicted to localise to the cytoplasm and nucleus and is expressed in several structures, including central nervous system; digestive system; fin; inter-renal primordium; and reproductive system, with highest levels of expression in the liver (Zhao et al., 2006). The gene is also expressed in earlier zebrafish development from gastrula-50%-epiboly-stage, 5.25hpf to 5.66hpf until adult-stage, 90dpf to 730dpf (Zfn.org—ZFIN ID: ZDB-GENE-070130-1). In zebrafish, nr0b1 mutants developed mostly as males (Chen et al., 2016). Mutations of the gene nr0b1 in zebrafish larvae, 13dpf, caused a decreased expression of the female-specific aromatase-encoding gene, cyp19a1a. This suggests that nr0b1 works during the early bipotential stage to promote female development (Chen et al., 2016).
Long-term oestradiol (E2) administration to the protandrous black porgy fish resulted in high levels of nr0b1 (Wu et al., 2008). Additionally, RA has been widely used to treat photoaged skin cancer and psoriasis (Lalli, 2014). Vitamin A is also crucial in activating Retinoic Acid Receptors (RAR). Both tretinoin and isotretinoin, are two synthetic forms of vitamin A. RA treatment has been found to increase the expression of nr0b1 (Nagl et al., 2009). These previous findings highlight changes in expression of nr0b1 after application of steroid-based and retinoid drugs, therefore, nr0b1 should not be used to indicate confounding sex identity with drug exposure (Li et al., 2016b).
4.1.4.5 Lhcgr
Human orthologs of the gene, lhcgr, are implicated in leydig cell tumour; breast cancer; and gonadal disease. It is expressed in digestive system; gonad; and hypophysis from juvenile-stage; 45dpf to 89dpf through to adult-stage; 90dpf to 730dpf (Zfn.org—ZFIN ID: ZDB-GENE-040806-3). Recent studies have shown that lhcgr is more important in controlling female reproduction compared to its role in male development (Zhang et al., 2015; Song et al., 2020). In female, lhcgr plays a vital role in triggering final oocyte maturation and ovulation. Lhcgr, a follicle cell-specific marker is expressed in follicles and isolated follicle layers, interacts with female-specific, gdf9, and is thought to be essential for ovulation (Wang et al., 2005; Li et al., 2016a). Lhcgr encodes luteinizing hormone subunit beta, lhβ; the transcription levels of the pituitary glycoprotein hormone, lhβ, were upregulated in cyp17a1-deficient fish pituitary (Zhai et al., 2018).
Oestradiol (E2)-based drugs have been found to alter the expression of lhcgr. Use of 17β-oestradiol (E2), but not testosterone, changed the regulated the expression of lhcgr (Liu et al., 2011). Recent studies have shown that lhcgr is expressed at high levels in various gynaecological cancers and so it was predicted be a potentially important biomarker for cancer diagnosis and therapy, in particular for ovarian cancer (Xiong et al., 2019; Zhong et al., 2019). FSH receptor binding inhibitor (FRBI) suppresses the expression of lhcgr mRNA and protein in sheep cumulus-oocyte complexes (COCs) (Wei et al., 2018). The mechanism of FRBI as an antagonist of lhcgr also has potential in aiding ovarian diseases and ovarian and follicular functioning in order to promote fertility in humans and animals (Wei et al., 2017; Lai et al., 2018). Although lhcgr has been identified as a potential drug target for cancer therapy and methods to apply this have been refined, there have so far been very few studies showing how these anti-cancer drugs affect the expression of lhcgr (Cavalli et al., 2016; Gao et al., 2018). Further studies of lhcgr expression levels in relation to drug application are needed to reveal if it can be used as an early expressed gene for unravelling confounding sex variables during anti-cancer drug development.
4.1.4.6 Cyp19a1a
Cyp19a1a is expressed in brain; eye; female organism; the ovary, ovarian granulosa cells; and the reproductive system; with highest expression level in female gonad (Chiang et al., 2001). The expression of the gene begins at zygote-1-cell-stage, through to adulthood (Zfn.org—ZFIN ID: ZDB-GENE-990415-43). Cyp19a1a is expressed at the highest level in the female gonad at juvenile stage, 30 to 44 days. It is also expressed in the immature gonad, at larval stage, 14 to 20dpf; in the gonad and female organism at larval stage, 21 to 44dpf; and in the mature ovarian follicle at fully-formed stage (Santos et al., 2007; Yu et al., 2014; Valdivieso et al., 2019). Cyp19a1a is a female promoting gene; the importance of cyp19a1a in the process of female sex differentiation was demonstrated when the silencing of the gene produced only male offspring (Lau et al., 2016; Yin et al., 2017). Other studies complement these findings whereby cyp19a1a is stated to be a female sex-determining gene (Santos et al., 2007; Jørgensen et al., 2008; Tong et al., 2010a; Rodríguez-Marí et al., 2010; Chen et al., 2017; Crowder et al., 2018b; Valdivieso et al., 2019) The NGS transcript data described here also complies with these findings. Taken together this implies that cyp19a1a could be used as an earlier genetic marker for indicating female sex-determination.
The expression of cyp19a1a in the gonad of labeo rohita fry changed after treated with fadrozole, an aromatase inhibitor. Equally treatment with 17β-oestradiol (E2) had no effect on the expression of cyp19a1a However, other studies have shown that two endocrine disrupting chemicals (EDCs) - an exogenous oestrogen 17β-oestradiol (E2) and the oestrogen mimic 4-n-nonylphenol (NP) cause a significant reduction in the expression of cyp19a1a in ovarian tissues with complete inhibition at the higher concentrations; therefore, implying that drugs with these actions can affect the expression of cyp19a1a (Shanthanagouda et al., 2013). More studies should be carried out to clarify the effect of different drugs on the expression of cyp19a1a, so that its use as an early sex-determining gene for uncovering sex-related difference with drug application can be definitively stated.
4.1.4.7 Foxl2a
Foxl2a is expressed in female zebrafish; in the gonad; ovarian follicle and ovary at larval-stage, 14dpf to 20dpf until breeding adult, 90dpf to 730dpf (Cocquet et al., 2002; Uchida et al., 2002; Tong et al., 2010a). Our data also shows adult female and juvenile, 28dpf, expression of foxl2a in zebrafish. Foxl2a in combination with other genes is essential to regulate zebrafish hormone producing ovarian cells and is necessary for ovary development, growth, maintenance and folliculogenesis; its loss significantly affects female development (Selman et al., 1993; Uchida et al., 2002; Cocquet et al., 2003; Tong et al., 2010b; Herpin and Schartl, 2015; Kossack and Draper, 2016). This indicates the importance of foxl2a in preventing female-to-male sex reversal by repressing testis determining genes, dmrt1 and amh (Uchida et al., 2002; Uhlenhaut et al., 2009; Matson et al., 2011; Schlessinger et al., 2011; Huang et al., 2017; Yang et al., 2017; Crowder et al., 2018b). Based on transcript data in this study, a 1.8-fold increase of foxl2a gene expression was seen in adult female zebrafish compared to males. This provides valid evidence for its use as an early female gene marker in zebrafish.
Mutations in foxl2a have been found to cause several diseases; blepharophimosis syndrome (BPES), a complex eyelid malformation associated with premature ovarian failure, ovarian dysfunction–polycystic ovary syndrome (PCOS) and granulosa cell tumours and keloid scars (Cocquet et al., 2002; Fleming et al., 2010; Verdin and De Baere, 2012). Regulation of expression of foxl2 or its downstream targets with drugs may provide a potential target for therapy of these diseases but these must first be tested for clarification of this (Verdin and De Baere, 2012). Studies have shown that after treatment with oestradiol and 2-methyl-testosterone, foxl2 was upregulated and downregulated, respectively (Jiang et al., 2011). This suggests that foxl2 would be affected by steroid-based drug application and therefore, not suitable for uncovering sex-related variables to steroid-based drug treatment (Wu et al., 2019).
4.2 Category two: Drug exposure thought not to influence expression of certain previously discovered early sex-determining genes in Danio rerio–More research needed
In literature, it was unclear whether certain sex-determining genes are affected by drug usage. In some cases, previous studies could be found which identify drugs which affect sex determining gene expression levels and in other studies the same genes were identified as potential therapeutic targets, but no studies could be found which link these and indicate that drugs directly alter the expression of these gene. Therefore, further investigation into their expression level with exposure to various drugs is needed to clarify their use in pharmaceutical trials to reveal confounding sex-related differences.
4.2.1 Formerly selected male sex-determining genes not influenced by drug activity
4.2.1.1 Gapdhs
King et al. (2020) The gapdhs gene encodes a protein and the testis-specific enzyme, glyceraldehyde-3-phosphate dehydrogenase spermatogenic, gapdhs. They are vital for energy production, sperm motility to ensure male fertility (Welch et al., 1995; Kuravsky et al., 2011). Gapdhs is expressed in adult male zebrafish as well as during development and in juveniles (Liu et al., 2013; King et al., 2020). This potential target is known to belong to a well-known druggable protein family which implies that gapdhs can be used as an early sex determining gene without its expression level being affected by drug application (Gottwald et al., 2006). More research is needed to clarify this finding.
4.2.2 Novel male sex-determining genes not influenced by drug activity
4.2.2.1 Apobb1
Apobb1 is expressed in intestine; liver; and yolk syncytial layer. Expression occurs in early development from gastrula 50% epiboly-stage and continues to adulthood (Otis et al., 2015). Apobb1 is expressed at lower levels with vitellogenin and not expressed in the ovary (Levi et al., 2012). There is a lack of studies regarding expression levels of apobb1 comparing adult male and female zebrafish. Our data shows on average more than a two-fold increase in transcripts of male zebrafish compared to female zebrafish as well as expression in juvenile. These findings give reason for the selection of apobb1 as a potential early male sex-determining gene. No previous literature indicates that drug activity influences apobb1 expression levels, further studies should also be carried out to prove this in order to highlight apobb1 as an early male sex-determining gene for uncovering sex-related variables with drug application.
4.2.2.2 Spata6
In mice and zebrafish, spata6 is predicted to localise to the sperm connecting piece during spermiogenesis. The sperm connecting piece is essential for linking the developing flagellum to the head in late spermiogenesis. When the expression of spata6 is disrupted males are sterile; this is due to the disruption to the formation of the sperm connecting piece which leads to acephalic spermatozoa in both the epididymis and ejaculates (Yuan et al., 2015; Sujit et al., 2020). In zebrafish, spata6 is expressed in neural tube; spinal cord; and ventral mesoderm, spermatids and mature spermatozoa (Yuan et al., 2015). Spata6 is expressed from 20dpf and increases to adulthood (Oh et al., 2003). Spata6 is a known pro-male gene and was upregulated in the ovaries of fish living in high-population density groups, showing stress from environmental factors causes masculinisation of zebrafish populations (Valdivieso et al., 2019). Previous literature taken with our NGS results show spata6 to be an early male sex-determining gene.
Spata6 is thought to play a key role in testicular germ cell tumours (TGCTs). Downregulation of spata6 results in apoptosis of TGCTs (Huo et al., 2015). This indicates that spata6 could be used as a potential gene target for drug therapy. However, more studies are needed to prove or disprove this and therefore indicate whether it can be used when reveal confounding sex variables during drug trials.
4.2.2.3 Dmrt1
Dmrt1 is involved in gamete generation; male gonad development; and male sex. Dmrt1, which acts downstream of male SYR, has previously been shown to be more strongly expressed in males than in females (Jørgensen et al., 2008; Chen et al., 2017). It is known to maintain the male fate and control gonadal differentiation in vertebrate and non-vertebrate species (Matson et al., 2011; King et al., 2020). This is supported by evidence of higher expression of dmrt1 in male when compared with female zebrafish in NGS data from this study. It is expressed in 14 organs with highest expression level in testis (fully formed stage) of zebrafish. The expression of dmrt1 was also found to be expressed in juveniles throughout development; it is highly expressed in blastula stage of zebrafish (3.33h–3.66 h) through to breeding adult stage (90d-730 d) (Jørgensen et al., 2008).
Dmrt1 is thought to have an influence on sex-specific patterns of childhood asthma, and its expression in testis tissue suggests a potential involvement in hormone regulation (Schieck et al., 2016). It may therefore act as a potential drug target gene for asthma suffers and steroid-based genes may influence its expression. If drugs are found which change the expression of dmrt1 for therapeutic reasons, its use in human drug development as an early sex-determining gene for identification of sex-related difference is not appropriate.
4.2.3 Formerly selected female sex-determining genes not influenced by drug activity
4.2.3.1 Bmp15
King et al. (2020) Bmp15 is expressed in the gonads of juvenile and adult female zebrafish as well as in the ovaries of humans (Fitzpatrick et al., 1998; Paulini and Melo, 2011; King et al., 2020). It negatively regulates oocyte development, function and maturation through its own activity as well as by increasing expression of other genes which convert androgens to oestrogens, namely, cyp19a1a or decreasing the expression of the male determining gene, dmrt1 (Otsuka et al., 2000; Clelland et al., 2006; Matson et al., 2011; Yan et al., 2017; Hosseini et al., 2019). Loss of bmp15 in adult zebrafish results in a female-to-male sex reversal (Dranow et al., 2016; Kossack and Draper, 2016; Crowder et al., 2018a). In Humans, loss of bmp15 leads to disruption in ovary function (Galloway et al., 2000; Di Pasquale et al., 2004). Drug application has not been previously shown to affect bmp15 expression and so its use as a female sex-identifying gene in drug discovery could be warranted. More tests are needed to prove that bmp15 expression is not affected by drug application.
4.2.3.2 Chr4
Chr4 expression has previously been linked to female sex determination in zebrafish and is expressed in juvenile stages (Charlesworth et al., 2005; Anderson et al., 2012; Howe et al., 2013; King et al., 2020). After research was carried out to discover changes in expression of chr4 with drug application, it was found that there were no previous studies which indicate that expression of chr4 is affected by drug application. Further studies should be carried out to clarify this result and its potential use for uncovering sex-related variables during drug trials.
4.2.3.3 Zpb3
Zpb3 gene expression is partly influenced by folliculogenesis-specific basic helix-loop-helix, Fig alpha, Fig α, which is needed in female development (Zfn.org—ZFIN ID: ZDB-GENE-031121-1, King et al., 2020). The germ cell-specific transcription factor, Fig α has been found in mammals and zebrafish where it is necessary for ovary differentiation and preservation (Soyal et al., 2000; Kanamori et al., 2006; Jørgensen et al., 2008; Schlessinger et al., 2011). As previously shown zpb3 is more highly expressed in oocytes and can be used as an early marker for female (Onichtchouk et al., 2003). No previous literature indicates that drug activity influences zpb3 expression levels. Further studies are needed to indicate the use of zpb3 to reveal sex differences with drug exposure during pharmaceutical trials.
4.2.3.4 Gyg1a
Gyg1a is highly expressed in developing, juvenile and female adult zebrafish (Maanen et al., 1999; Wen et al., 2005; King et al., 2020). No previous study has provided evidence that drugs can influence the expression of gyg1a, more research is needed to clarify this finding.
4.2.3.5 Kpna2
Kpna2 is highly expressed in developing and adult female zebrafish gonads and its use as an early female sex-determining gene has been previously discussed (Major et al., 2011; Cui et al., 2020; King et al., 2020). Kpna2 overexpression has been linked to various human cancers, including, breast, ovarian endodermal sinus and ovarian primitive germ cell tumours (Huang et al., 2018; Christiansen and Dyrskjøt, 2020). An inhibitor of kpna2 could in theory act as a prospective drug to reduce growth of cancer cells. The genes, IRF1 has been shown to negatively regulate the expression of kpna2. Therefore, drugs which target IRF1 gene could disrupt normal kpna2 expression. Additionally, some studies have shown that kpna2 may be dependent on p53 nuclear translocation to support autophagy in order to allow chemoresistance and metastasis of certain cancerous cells, OSCC (Oral squamous cell carcinoma) (Lin et al., 2018). Although the use of kpna2 as a therapeutic target has been previously identified, studies regarding the relative gene expression levels with drug application are lacking and are required before kpna2 can be used as an early sex-determining gene for uncovering sex differences during drug trials.
4.2.4 Novel female sex-determining genes not influenced by drug activity
4.2.4.1 Ftz-f1
Ftz-f1 is expressed in several structures, including digestive system; endocrine system; head; nervous system; and slow muscle cell. It is expressed from zygote:1-cell-stage, 0.0h–0.75 h until breeding-adult-stage, 90days to 730 days (Tong et al., 2010b) (Zfn.org—ZFIN ID: ZDB-GENE-990415-79). At around 12dpf, when female germ cells start to rapidly grow, the initial sign of sexual dimorphism occurs. At this time somatic cells start to increase expression of ftz-f1 transcription factors (Tong et al., 2010b). Ftz-f1 has been previously associated with female sex-determination in zebrafish (Cao et al., 2012). The zebrafish transcript data which is presented here showed an increase of 1.6-fold in females compared to males. Ftz-f1 can therefore, be considered for identification of female zebrafish early in development.
Overexpression of ftz-f1 enhances activation, implicating ftz-f1 as a critical component of the juvenile hormone response (Dubrovsky et al., 2011). Additionally, amino acid sequencing of ftz-f1 gene revealed that the gene protein is a member of the nuclear hormone receptor superfamily (Lavorgna et al., 1991). This research could therefore imply that the expression of ftz-f1 can be altered by hormone-based drugs. However, no previous studies revealed that ftz-f1 expression was altered by drug application and more research should be carried into this before it can be classified as an earlier sex-determining gene for use in revealing sex-related variables occurring from use of drugs.
4.2.4.2 gdf9
Gdf9 is an oocyte specific factor encoding proteins which are secreted by oocytes (Clelland et al., 2006; Paulini and Melo, 2011). It is involved in juvenile and adult female sex-determination in zebrafish and is essential for ovarian follicle development in mammals (Silva et al., 2005; Clelland et al., 2006; Gilchrist et al., 2008; Oron et al., 2010; Chen et al., 2017; Strauss and Williams, 2019). It is expressed in fertilised egg; gonad; and immature gonad (Fitzpatrick et al., 1998). The expression of gdf9 starts in the gonads before sex differentiation, and its expression dramatically increases in the differentiated ovary, with a similar expression pattern to that of cyp19a1a. Gdf9 significantly suppresses the expression of male amh while increasing that of activin beta subunits, inhbaa and inhbb. Knock-down expression of gdf9 creates a male-biased sex ratio. Gdf9 promotes follicle-stimulating hormone (FSH)-induced P4 production and STAR expression (Leitão et al., 2014; Li et al., 2019b). P4 is an endogenous steroid and progestogen sex hormone involved in the menstrual cycle, pregnancy, and embryogenesis of humans and other species (Chen et al., 2017). Our transcript data along with previous knowledge supplies evidence for gdf9 as an early gene marker in zebrafish sex-determination.
In humans, gdf9 was recently shown to be a regulator of the development and spread of certain cancer cells including those of the breast and kidney. Gdf9 expression is decreased in clear-cell renal cell carcinoma (CCRCC) which is linked to long-term survival of patients. Although this indicates that gdf9 is a potential tumour suppressor in CCRCC which can be used a potential target for tumour drug development, very few studies have looked at drugs which target the gene for therapeutic reasons nor have they indicated the relative level of gdf9 expression (Hanavadi et al., 2007; Du et al., 2012; Du et al., 2014).
4.2.4.3 Ipo4
Ipo4 has a role in protein import into nucleus. It is expressed in 33 zebrafish organs with forth highest expression in the female gonad and 14th highest expression in the mature ovarian follicle [bgee. org—gene: ipo4— ENSDARG00000041493—Danio rerio (zebrafish)]. Ipo4 is expressed during development from 0dpf until 56dpf, peaking in expression level at 30dpf (Major et al., 2011). In this study, adult female zebrafish transcript data for ipo4 is around four times higher in adult females when compared to the adult males indicating it as a potential marker for early female sex-determination.
Ipo4 was found to be downregulated in both SARS-Cov-2 infection and aging (Belyaeva et al., 2020a). It could, therefore, provide a target for drug therapy for these. Ipo4 is thought to be a key gene contributing to the progression and development of gastric and other cancers (Xu et al., 2019; Zhou et al., 2020). The drugs, guanosine triphosphate and mkc-1 both have both been linked to the gene ipo4. Mkc-1 has an antineoplastic agent, meaning it inhibits the growth and spread of tumours or malignant cells (https://www.genecards.org/cgi-bin/carddisp.pl?gene=IPO4). This also highlights ipo4 as a gene target for cancer drug therapy. However, more studies are also needed for clarification of the action of these drugs on the expression of ipo4.
4.3 Category three: Drug exposure does not influence expression of certain previously discovered early sex-determining genes in Danio rerio–Evidence exists already
None of the previously selected early sex-determining genes (King et al., 2020) as well as none of the additional genes highlighted in this study were found to not be influenced by drug usage. This reveals the need for research to be focused in this area. Genes in which previous studies show that their expression is not influenced by drug exposure would have great potential to be used as candidate genes for revealing confounding sex variables in drug discovery. The genes highlighted in this study which show potential to be used for revealing sex variables in pharmaceutics trials (see category two) should now be evaluated further to clarify which of the two categories (category one or category three) they can be placed into.
5 Conclusion
We reveal here several male and female marker genes in early zebrafish development which are not known to be disrupted during drug exposure. These show great potential to unlock the confounding variable of sex in order to reveal sex-related differences during drug trials. These genes should be tested in the near future to guarantee their usefulness within drug development and to show if the effect on the molecular level would lead to a skewed sex ratio in environmental toxicity studies (Zhang et al., 2013). Results from NGS data and the literature review clearly confirm the differences in sensitivity within a species and indicate that females are more susceptible to drug exposure than males (Chen et al., 2021).
We uncover that sex-related genes are very often affected in some way by drugs. From this we deduce that there will be a sex-based difference in susceptibility to drugs, so we emphasise the importance of producing males and female-specific pharmaceuticals at specific sex doses. The model of the healthy male volunteer in pharmaceutical studies should be re-thought to allow safe dosage of drugs to prescribed for females. Drugs interfere with hormones and the gene expression pathway, likewise, hormones and genes influence drug metabolism. So, while it can be said that the sex of an individual is confounding the drug, likewise the drug is also confounding the sex.
Our findings emphasise the importance of promoting sex-specific drugs which would positively affect human pharmaceutical drug discovery. This paper paves the way in connecting the vital ‘missing link’ between juvenile sex-related genes and drug response which will greatly impact the future of pharmaceutical trials.
Data availability statement
The data discussed in this publication have been deposited in NCBI’s Gene Expression Omnibus (Edgar et al., 2002) and are accessible through GEO Series accession number GSE234104 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE234104).
Ethics statement
Ethical review and approval was not required for the animal study because all experiments were conducted with 28-day-old and adult zebrafish, which had already been euthanised with Tricaine by the breeder. No experiments were carried out with live adult nor subadult fish. Accordingly, welfare guidelines were not applicable.
Author contributions
AK: NGS data processing and manuscript writing. AZ: Laboratory and manuscript planning, input, and overview.
Acknowledgments
We would like to thank Boris Kolvenbach for his support with gene alignment in galaxy.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Al-Sheikh, H., Watson, A. J., Lacey, G. A., Punt, P. J., MacKenzie, D. A., Jeenes, D. J., et al. (2004). Endoplasmic reticulum stress leads to the selective transcriptional downregulation of the glucoamylase gene in Aspergillus Niger. Mol. Microbiol. 53, 1731–1742. doi:10.1111/j.1365-2958.2004.04236.x
Amler, L. C., Agus, D. B., LeDuc, C., Sapinoso, M. L., Fox, W. D., Kern, S., et al. (2000). Dysregulated expression of androgen-responsive and nonresponsive genes in the androgen-independent prostate cancer xenograft model CWR22-R. Cancer Res. 60, 6134–6141. doi:10.1038/86980
Anderson, J., Braasch, I., Amores, A., Hohenlohe, P., Rodrı, A., Batzel, P., et al. (2012). Multiple sex-associated regions and a putative sex chromosome in zebrafish revealed by RAD mapping and population genomics. PLoS One 7, 1–14. doi:10.1371/journal.pone.0040701
Belrose, G., Gross, A., Olindo, S., Lézin, A., Dueymes, M., Komla-Soukha, I., et al. (2011). Effects of valproate on Tax and HBZ expression in HTLV-1 and HAM/TSPT lymphocytes. Blood 118, 2483–2491. doi:10.1182/blood-2010-11-321364
Belyaeva, A., Cammarata, L., Radhakrishnan, A., Squires, C., Yang, K. D., Shivashankar, G. V., et al. (2020a). Causal network models of SARS-CoV-2 expression and aging to identify candidates for drug repurposing. Nat. Commun. 1, 1–40. doi:10.1038/s41467-021-21056-z
Belyaeva, O. V., Adams, M. K., Popov, K. M., and Kedishvili, N. Y. (2020b). Generation of retinaldehyde for retinoic acid biosynthesis. Biomolecules 10, 5–17. doi:10.3390/biom10010005
Cao, J., Chen, J., Jiang, Z., Luo, Y., and Gan, X. (2012). Molecular cloning and expression analysis of FTZ-F1 in the GIFT tilapia, Oreochromis niloticus. Zool. Res. 33, 40–46. doi:10.3724/SP.J.1141.2012.E03-04E40
Cavalli, R., Soster, M., and Argenziano, M. (2016). Nanobubbles: A promising efficient tool for therapeutic delivery. Ther. Deliv. 7, 117–138. doi:10.4155/tde.15.92
Charlesworth, D., Charlesworth, B., and Marais, G. (2005). Steps in the evolution of heteromorphic sex chromosomes. Hered. (Edinb). 95, 118–128. doi:10.1038/sj.hdy.6800697
Chen, K., Wu, M., Chen, C., Xu, H., Wu, X., and Qiu, X. (2021). Impacts of chronic exposure to sublethal diazepam on behavioral traits of female and male zebrafish (Danio rerio). Ecotoxicol. Environ. Saf. 208, 111747. doi:10.1016/j.ecoenv.2020.111747
Chen, S., Zhang, H., Wang, F., Zhang, W., and Peng, G. (2016). nr0b1 (DAX1) mutation in zebrafish causes female-to-male sex reversal through abnormal gonadal proliferation and differentiation. Mol. Cell. Endocrinol. 433, 105–116. doi:10.1016/j.mce.2016.06.005
Chen, W., Liu, L., and Ge, W. (2017). Expression analysis of growth differentiation factor 9 (Gdf9/gdf9), anti-müllerian hormone (Amh/amh) and aromatase (Cyp19a1a/cyp19a1a) during gonadal differentiation of the zebrafish, Danio rerio. Biol. Reprod. 96, 401–413. doi:10.1095/biolreprod.116.144964
Chiang, E. F., Yan, Y., Tong, S., Hsiao, P., Guiguen, Y., Postlethwait, J., et al. (2001). Characterization of duplicated zebrafish cyp19 genes. J. Exp. Zool. 290, 709–714. doi:10.1002/jez.1121
Christiansen, A., and Dyrskjøt, L. (2020). The functional role of the novel biomarker karyopherin α 2 (KPNA2) in cancer. Cancer Lett. 331, 18–23. doi:10.1016/j.canlet.2012.12.013
Clelland, E., Kohli, G., Campbell, R. K., Sharma, S., Shimasaki, S., and Peng, C. (2006). Bone morphogenetic protein-15 in the zebrafish ovary: Complementary deoxyribonucleic acid cloning, genomic organization, tissue distribution, and role in oocyte maturation. Endocrinology 147, 201–209. doi:10.1210/en.2005-1017
Cocquet, J., De Baere, E., Gareil, M., Pannetier, M., Xia, X., Fellous, M., et al. (2003). Structure, evolution and expression of the FOXL2 transcription unit. Cytogenet. Genome Res. 101, 206–211. doi:10.1159/000074338
Cocquet, J., Pailhoux, E., Jaubert, F., Servel, N., Xia, X., Pannetier, M., et al. (2002). Evolution and expression of FOXL2. J. Med. Genet. 39, 916–921. doi:10.1136/jmg.39.12.916
Costaridis, P., Horton, C., Zeitlinger, J., Holder, N., and Maden, M. (1996). Endogenous retinoids in the zebrafish embryo and adult. Dev. Dyn. 205, 41–51. doi:10.1002/(SICI)1097-0177(199601)205:1<41:AID-AJA4>3.0.CO;2-5
Crowder, C. M., Lassiter, C., and Gorelick, D. (2018a). Nuclear androgen receptor regulates testes organization and oocyte maturation in zebrafish. Endocrinology 159, 980–993. doi:10.1210/en.2017-00617
Crowder, C. M., Romano, S. N., and Gorelick, D. A. (2018b). G protein-coupled estrogen receptor is not required for sex determination or ovary function in zebrafish. Endocrinology 159, 3515–3523. doi:10.1210/en.2018-00685
Cui, X., Jing, X., Wu, X., Xu, J., Liu, Z., Huo, K., et al. (2020). Analyses of dna methylation involved in the activation of nuclear karyopherin alpha 2 leading to identify the progression and prognostic significance across human breast cancer. Cancer Manag. Res. 12, 6665–6677. doi:10.2147/CMAR.S261290
Dar, A. A., Majid, S., Nosrati, M., De Semir, D., Federman, S., and Kashani-Sabet, M. (2010). Functional modulation of IGF-binding protein-3 expression in melanoma. J. Invest. Dermatol. 130, 2071–2079. Available at:. doi:10.1038/jid.2010.70
Di Pasquale, E., Beck-Peccoz, P., and Persani, L. (2004). Hypergonadotropic ovarian failure associated with an inherited mutation of human bone morphogenetic protein-15 (BMP15) gene. Am. J. Hum. Genet. 75, 106–111. doi:10.1086/422103
Ding, K., Li, W., Zou, Z., Zou, X., and Wang, C. (2014). CCNB1 is a prognostic biomarker for ER+ breast cancer. Med. Hypotheses 83, 359–364. doi:10.1016/j.mehy.2014.06.013
Dingare, C., Niedzwetzki, A., Klemmt, P. A., Godbersen, S., Fuentes, R., Mullins, M. C., et al. (2018). The hippo pathway effector taz is required for cell morphogenesis and fertilization in zebrafish. Dev 145, dev167023–13. doi:10.1242/dev.167023
Dlugos, C. A., Brown, S. J., and Rabin, R. A. (2011). Gender differences in ethanol-induced behavioral sensitivity in zebrafish. Alcohol 45, 11–18. doi:10.1016/j.alcohol.2010.08.018
Dolgin, E. (2020). IGF-1R drugs travel from cancer cradle to Graves. Nat. Biotechnol. 38, 385–388. doi:10.1038/s41587-020-0481-8
Dranow, D. B., Hu, K., Bird, A. M., Lawry, S. T., Adams, M. T., Sanchez, A., et al. (2016). Bmp15 is an oocyte-produced signal required for maintenance of the adult female sexual phenotype in zebrafish. PLoS Genet. 12, e1006323–e1006324. doi:10.1371/journal.pgen.1006323
Du, P., Ye, L., Li, H., Ruge, F., Yang, Y., and Jiang, W. G. (2012). Loss of expression of growth differentiation factor-9 (GDF9) in human kidney cancer and regulation of growth and migration of kidney cancer cells by GDF9. Anticancer Res. 32, 4375–4383.
Du, P., Ye, L., Yang, Y., and Jiang, W. G. (2014). Reduced expression of growth and differentiation factor-9 (GDF9) is associated with aggressive behaviour of human clear-cell renal cell carcinoma and poor patient survival. Anticancer Res. 34, 6515–6520.
Dubrovsky, E. B., Dubrovskaya, V. A., Bernardo, T., Otte, V., DiFilippo, R., and Bryan, H. (2011). The Drosophila FTZ-F1 nuclear receptor mediates juvenile hormone activation of E75A gene expression through an intracellular pathway. J. Biol. Chem. 286, 33689–33700. doi:10.1074/jbc.M111.273458
Duran, N., Martinez, D., Silveira, C., Duran, M., de Moraes, A., Simoes, M., et al. (2015). Graphene oxide: A carrier for pharmaceuticals and a scaffold for cell interactions. Curr. Top. Med. Chem. 15, 309–327. doi:10.2174/1568026615666150108144217
Edgar, R., Domrachev, M., and Lash, A. E. (2002). Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 30 (1), 207–210.
Fisher, J. L., Jones, E. F., Flanary, V. L., Williams, A. S., Ramsey, E. J., and Lasseigne, B. N. (2022). Considerations and challenges for sex - aware drug repurposing. Biol. Sex. Differ. 13, 13–23. doi:10.1186/s13293-022-00420-8
Fitzpatrick, S. L., Sindoni, D. M., Shughrue, P. J., Lane, M. V., Merchenthaler, I. J., and Frail, D. E. (1998). Expression of growth differentiation factor-9 messenger ribonucleic acid in ovarian and nonovarian rodent and human tissues. Endocrinology 139, 2571–2578. doi:10.1210/endo.139.5.6014
Fleming, N. I., Knower, K. C., Lazarus, K. A., Fuller, P. J., Simpson, E. R., and Clyne, C. D. (2010). Aromatase is a direct target of FOXl2: C134W in granulosa cell tumors via a single highly conserved binding site in the ovarian specific promoter. PLoS One 5, 143899–e14410. doi:10.1371/journal.pone.0014389
Franconi, F., Brunelleschi, S., Steardo, L., and Cuomo, V. (2007). Gender differences in drug responses. Pharmacol. Res. 55, 81–95. doi:10.1016/j.phrs.2006.11.001
Galloway, S. M., McNatty, K. P., Cambridge, L. M., Laitinen, M. P. E., Juengel, J. L., Jokiranta, T. S., et al. (2000). Mutations in an oocyte-derived growth factor gene (BMP15) cause increased ovulation rate and infertility in a dosage-sensitive manner. Nat. Genet. 25, 279–283. doi:10.1038/77033
Gao, Y., Chen, S., Li, W., Wang, H., Xiao, K., Wu, L., et al. (2018). An experimental study of ovarian cancer imaging and therapy by paclitaxel-loaded phase-transformation lipid nanoparticles combined with low-intensity focused ultrasound. Biochem. Biophys. Res. Commun. 504, 203–210. doi:10.1016/j.bbrc.2018.08.156
Gartlehner, G., Chapman, A., Strobelberger, M., and Thaler, K. (2010). Differences in efficacy and safety of pharmaceutical treatments between men and women: An umbrella review. PLoS One 5, e11895–e11898. doi:10.1371/journal.pone.0011895
Gazon, H., Belrose, G., Terol, M., Meniane, J. C., Mesnard, J. M., Césaire, R., et al. (2016). Impaired expression of DICER and some microRNAs in HBZ expressing cells from acute adult T-cell leukemia patients. Oncotarget 7, 30258–30275. doi:10.18632/oncotarget.7162
Giatromanolaki, A., Fasoulaki, V., Kalamida, D., Mitrakas, A., Kakouratos, C., Lialiaris, T., et al. (2019). CYP17A1 and androgen-receptor expression in prostate carcinoma tissues and cancer cell lines. Curr. Urol. 13, 157–165. doi:10.1159/000499276
Gilchrist, R. B., Lane, M., and Thompson, J. G. (2008). Oocyte-secreted factors: Regulators of cumulus cell function and oocyte quality. Hum. Reprod. Update 14, 159–177. doi:10.1093/humupd/dmm040
Goldsmith, P. (2004). Zebrafish as a pharmacological tool: The how, why and when. Curr. Opin. Pharmacol. 4, 504–512. doi:10.1016/j.coph.2004.04.005
Gottwald, U., Davies, B., Fritsch, M., and Habenicht, U. F. (2006). New approaches for male fertility control: HE6 as an example of a putative target. Mol. Cell. Endocrinol. 250, 49–57. doi:10.1016/j.mce.2005.12.024
Greenblatt, D. J., and Von Moltke, L. L. (2008). Gender has a small but statistically significant effect on clearance of CYP3A substrate drugs. J. Clin. Pharmacol. 48, 1350–1355. doi:10.1177/0091270008323754
Groh, K. J., Nesatyy, V. J., Segner, H., Eggen, R. I. L., and Suter, M. J. F. (2011). Global proteomics analysis of testis and ovary in adult zebrafish (Danio rerio). Fish. Physiol. Biochem. 37, 619–647. doi:10.1007/s10695-010-9464-x
Grunwald, D. J., and Eisen, J. S. (2002). Headwaters of the zebrafish — Emergence of a new model vertebrate. Nat. Rev. Genet. 3, 717–724. doi:10.1038/nrg892
Guiguen, Y., Fostier, A., Piferrer, F., and Chang, C. F. (2010). Ovarian aromatase and estrogens: A pivotal role for gonadal sex differentiation and sex change in fish. Gen. Comp. Endocrinol. 165, 352–366. Available at:. doi:10.1016/j.ygcen.2009.03.002
Hanavadi, S., Martin, T. A., Watkins, G., Mansel, R. E., and Jiang, W. G. (2007). The role of growth differentiation factor-9 (GDF-9) and its analog, GDF-9b/BMP-15, in human breast cancer. Ann. Surg. Oncol. 14, 2159–2166. doi:10.1245/s10434-007-9397-5
Hatzold, J., Beleggia, F., Herzig, H., Altmüller, J., Nürnberg, P., Bloch, W., et al. (2016). Tumor suppression in basal keratinocytes via dual non-cell-autonomous functions of a Na,K-ATPase beta subunit. Elife 5, 142777–e14330. doi:10.7554/eLife.14277
Herpin, A., and Schartl, M. (2015). Plasticity of gene-regulatory networks controlling sex determination: Of masters, slaves, usual suspects, newcomers, and usurpators. EMBO Rep. 16, 1260–1274. doi:10.15252/embr.201540667
Higano, B. C. S., Chielens, D., Raskind, W., Bryant, E., Flowers, M. E. D., Radich, J., et al. (1997). Characterization of adult α- and β-globin genes in the zebrafish. Blood 89, 688–700. doi:10.1182/blood.v89.2.688
Holdcroft, A. (2007). Gender bias in research: How does it affect evidence based medicine? J. R. Soc. Med. 100, 2–3. doi:10.1177/014107680710000102
Horie, M., and Kotani, T. (2016). Formation of mos RNA granules in the zebrafish oocyte that differ from cyclin B1 RNA granules in distribution, density and regulation. Eur. J. Cell Biol. 95, 563–573. doi:10.1016/j.ejcb.2016.10.001
Hosseini, S., Ha, N. T., Simianer, H., Falker-Gieske, C., Brenig, B., Franke, A., et al. (2019). Genetic mechanism underlying sexual plasticity and its association with colour patterning in zebrafish (Danio rerio). BMC Genomics 20, 341–417. doi:10.1186/s12864-019-5722-1
Houben, T., Bitorina, A. V., Oligschlaeger, Y., Jeurissen, M. L. J., Rensen, S., Köhler, S. E., et al. (2020). Sex-opposed inflammatory effects of 27-hydroxycholesterol are mediated via differences in estrogen signaling. J. Pathol. 251, 429–439. doi:10.1002/path.5477
Howe, K., Clark, M. D., Torroja, C. F., Torrance, J., Berthelot, C., Muffato, M., et al. (2013). The zebrafish reference genome sequence and its relationship to the human genome. Nature 496, 498–503. doi:10.1038/nature12111
Huang, L., Zhou, Y., Cao, X. P., Lin, J. X., Zhang, L., Huang, S. T., et al. (2018). KPNA2 promotes migration and invasion in epithelial ovarian cancer cells by inducing epithelial-mesenchymal transition via Akt/GSK-3β/Snail activation. J. Cancer 9, 157–165. doi:10.7150/jca.20879
Huang, S., Ye, L., and Chen, H. (2017). Sex determination and maintenance: The role of DMRT1 and FOXL2. Asian J. Androl. 19, 619–624. doi:10.4103/1008-682X.194420
Huo, S., Du, W., Shi, P., Si, Y., and Zhao, S. (2015). The role of spermatogenesis-associated protein 6 in testicular germ cell tumors. Int. J. Clin. Exp. Pathol. 8, 9119–9125.
Irie, T., and Seki, T. (2002). Retinoid composition and retinal localization in the eggs of teleost fishes. Comp. Biochem. Physiol. - B Biochem. Mol. Biol. 131, 209–219. doi:10.1016/s1096-4959(01)00496-1
Ising, M., Maccarrone, G., Brückl, T., Scheuer, S., Hennings, J., Holsboer, F., et al. (2019). FKBP5 gene expression predicts antidepressant treatment outcome in depression. Int. J. Mol. Sci. 20, 485–512. doi:10.3390/ijms20030485
Jalili, V., Afgan, E., Gu, Q., Clements, D., Blankenberg, D., Goecks, J., et al. (2020). The galaxy platform for accessible, reproducible and collaborative biomedical analyses: 2020 update. Nucleic Acids Res. 48, W395-W402–W402. doi:10.1093/nar/gkaa434
Jiang, W., Yang, Y., Zhao, D., Liu, X., Duan, J., Xie, S., et al. (2011). Effects of sexual steroids on the expression of foxl2 in Gobiocypris rarus. Comp. Biochem. Physiol. - B Biochem. Mol. Biol. 160, 187–193. doi:10.1016/j.cbpb.2011.08.005
Jin, Y., Chen, R., Sun, L., Wang, W., Zhou, L., Liu, W., et al. (2009). Enantioselective induction of estrogen-responsive gene expression by permethrin enantiomers in embryo-larval zebrafish. Chemosphere 74, 1238–1244. Available at:. doi:10.1016/j.chemosphere.2008.11.015
Jørgensen, A., Morthorst, J. E., Andersen, O., Rasmussen, L. J., and Bjerregaard, P. (2008). Expression profiles for six zebrafish genes during gonadal sex differentiation. Reprod. Biol. Endocrinol. 6, 25–12. doi:10.1186/1477-7827-6-25
Kamboj, R. C., Raghav, N., Mittal, A., Khurana, S., Sadana, R., and Singh, H. (2003). Effects of some antituberculous and anti-leprotic drugs on cathepsins B, H and L. Indian J. Clin. biochem. 18, 39–47. doi:10.1007/BF02867366
Kanamori, A., Yamamura, A., Koshiba, S., Lee, J., Orlando, E. F., Hori, H., et al. (2006). Methyltestosterone efficiently induces male development in the self-fertilizing hermaphrodite fish, kryptolebias marmoratus. Genesis 44, 495–503. doi:10.1002/dvg.20240
Kashimada, K., Svingen, T., Feng, C., Pelosi, E., Bagheri-fam, S., Harley, V. R., et al. (2011). Antagonistic regulation of Cyp26b1 by transcription factors SOX9/SF1 and FOXL2 during gonadal development in mice. FASEB J. 25, 3561–3569. doi:10.1096/fj.11-184333
Kejuan, F., Meirik, O., Yongang, D., Yan, C., Weijin, Z., and Fajans, P. (2007). Once-a-month contraceptive pills in China: A review of available evidence. Contraception 75, 337–343. doi:10.1016/j.contraception.2007.01.007
Kim, A. M., Tingen, C. M., and Woodruff, T. K. (2010). Sex bias in trials and treatment must end. Nature 465, 688–689. doi:10.1038/465688a
Kim, D., Paggi, J. M., Park, C., Bennett, C., and Salzberg, S. L. (2019). Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat. Biotechnol. 37, 907–915. Available at:. doi:10.1038/s41587-019-0201-4
Kimmel, C. B., Ballard, W. W., Kimmel, S. R., Ullmann, B., and Schilling, T. F. (1995). Stages of embryonic development of the zebrafish. Dev. Dyn. 203, 253–310. doi:10.1002/aja.1002030302
King, A. C., Gut, M., and Zenker, A. K. (2020). Shedding new light on early sex determination in zebrafish. Arch. Toxicol. 94, 4143–4158. doi:10.1007/s00204-020-02915-y
Kling, P., Norman, A., Andersson, P. L., Norrgren, L., and Förlin, L. (2008). Gender-specific proteomic responses in zebrafish liver following exposure to a selected mixture of brominated flame retardants. Ecotoxicol. Environ. Saf. 71, 319–327. doi:10.1016/j.ecoenv.2007.12.010
Knoll-Gellida, A., André, M., Gattegno, T., Forgue, J., Admon, A., and Babin, P. J. (2006). Molecular phenotype of zebrafish ovarian follicle by serial analysis of gene expression and proteomic profiling, and comparison wih the transcriptomes of other animals. BMC Genomics 7, 1–28. doi:10.1186/1471-2164-7-46
Kokane, S. S., and Perrotti, L. I. (2020). Sex differences and the role of estradiol in mesolimbic reward circuits and vulnerability to cocaine and opiate addiction. Front. Behav. Neurosci. 14, 74–21. doi:10.3389/fnbeh.2020.00074
Kondo, T., Kotani, T., and Yamashita, M. (2001). Dispersion of cyclin B mRNA aggregation is coupled with translational activation of the mRNA during zebrafish oocyte maturation. Dev. Biol. 229, 421–431. doi:10.1006/dbio.2000.9990
Kondo, T., Yanagawa, T., Yoshida, N., and Yamashita, M. (1997). Introduction of cyclin B induces activation of the maturation-promoting factor and breakdown of germinal vesicle in growing zebrafish oocytes unresponsive to the maturation-inducing hormone. Dev. Biol. 190, 142–152. doi:10.1006/dbio.1997.8673
Kossack, M., and Draper, B. (2016). Genetic regulation of sex determination and maintenance in zebrafish (Danio rerio). Physiol. Behav. 176, 100–106. doi:10.1016/bs.ctdb.2019.02.004
Kotani, T., Yasuda, K., Ota, R., and Yamashita, M. (2013). Cyclin b1 mRNA translation is temporally controlled through formation and disassembly of RNA granules. J. Cell Biol. 202, 1041–1055. doi:10.1083/jcb.201302139
Koubova, J., Menke, D. B., Zhou, Q., Capel, B., Griswold, M. D., and Page, D. C. (2006). Retinoic acid regulates sex-specific timing of meiotic initiation in mice. PNAS 103, 2474–2479. doi:10.1073/pnas.0510813103
Kuravsky, M. L., Aleshin, V. V., Frishman, D., and Muronetz, V. I. (2011). Testis-specific glyceraldehyde-3-phosphate dehydrogenase: Origin and evolution. BMC Evol. Biol. 11, 1–15. Available at: http://www.biomedcentral.com/1471-2148/11/160.
Lai, L., Shen, X., Liang, H., Deng, Y., Gong, Z., Wei, S., et al. (2018). Innate, T-and B-cell responses in acute human zika patients. Biomed. Res. Int. 1, 1–10. doi:10.1093/cid/cix732
Lalli, E. (2014). Role of orphan nuclear receptor DAX-1/NR0B1 in development, physiology, and disease. Adv. Biol. 1, 1–19. doi:10.1155/2014/582749
Lau, E. S., Zhang, Z., Qin, M., and Ge, W. (2016). Knockout of zebrafish ovarian aromatase gene (cyp19a1a) by TALEN and CRISPR/Cas9 leads to all-male offspring due to failed ovarian differentiation. Nat. Publ. Gr. 6, 37357–37414. doi:10.1038/srep37357
Lavorgna, G., Ueda, H., Clos, J., and Wu, C. (1991). FTZ-F1, a steroid hormone receptor-like protein implicated in the activation of fushi tarazu. Sci. (80) 252, 848–851. doi:10.1126/science.1709303
Leconte, I., Bailey, G., Davis-Bruno, K., Hew, K. W., Kim, J., Silva Lima, B., et al. (2011). Value of juvenile animal studies. Birth Defects Res. Part B - Dev. Reprod. Toxicol. 92, 292–303. doi:10.1002/bdrb.20331
Leet, J. K., Gall, E., and Sepúlveda, M. S. (2011). A review of studies on androgen and estrogen exposure in fish early life stages: Effects on gene and hormonal control of sexual differentiation. Appl. Toxicol. 31, 379–398. doi:10.1002/jat.1682
Leitão, C. C. F., Costa, J. J. D. N., Brito, I. R., De Melo Magalhães-Padilha, D., Almeida, A. P., De Figueiredo, J. R., et al. (2014). Effects of GDF-9 and FSH on mrna expression for FSH-R, GDF-9 and BMPs in in vitro cultured goat preantral follicles. Braz. Arch. Biol. Technol. 57, 200–208. doi:10.1590/s1516-89132014000200007
Levi, L., Ziv, T., Admon, A., Levavi-Sivan, B., and Lubzens, E. (2012). Insight into molecular pathways of retinal metabolism, associated with vitellogenesis in zebrafish. Am. J. Physiol. - Endocrinol. Metab. 302, 626–644. doi:10.1152/ajpendo.00310.2011
Li, B., Zhu, H. B., Song, G. D., Cheng, J. H., Li, C. Z., Zhang, Y., et al. (2019a). Regulating the CCNB1 gene can affect cell proliferation and apoptosis in pituitary adenomas and activate epithelial-to-mesenchymal transition. Oncol. Lett. 18, 4651–4658. doi:10.3892/ol.2019.10847
Li, C., Tan, X. F., Lim, T. K., Lin, Q., and Gong, Z. (2016a). Comprehensive and quantitative proteomic analyses of zebrafish plasma reveals conserved protein profiles between genders and between zebrafish and human. Sci. Rep. 6, 24329–24415. doi:10.1038/srep24329
Li, J., Liu, Z., Wang, D., and Cheng, C. H. K. (2011). Insulin-like growth factor 3 is involved in oocyte maturation in zebrafish. Biol. Reprod. 84, 476–486. doi:10.1095/biolreprod.110.086363
Li, J., Luo, W., Huang, T., and Gong, Y. (2019b). Growth differentiation factor 9 promotes follicle-stimulating hormone-induced progesterone production in chicken follicular granulosa cells. Gen. Comp. Endocrinol. 276, 69–76. Available at:. doi:10.1016/j.ygcen.2019.03.005
Li, Y. C., Luo, M. L., Guo, H., Wang, T. T., Lin, S. R., Chen, J. B., et al. (2016b). Identification of NR0B1 as a novel androgen receptor co-repressor in mouse Sertoli cells. Int. J. Mol. Med. 38, 853–860. doi:10.3892/ijmm.2016.2672
Liang, Y. Q., Huang, G. Y., Liu, S. S., Zhao, J. L., Yang, Y. Y., Chen, X. W., et al. (2015). Long-term exposure to environmentally relevant concentrations of progesterone and norgestrel affects sex differentiation in zebrafish (Danio rerio). Aquat. Toxicol. 160, 172–179. doi:10.1016/j.aquatox.2015.01.006
Lin, F., Gao, L., Su, Z., Cao, X., Zhan, Y., Li, Y., et al. (2018). Knockdown of KPNA2 inhibits autophagy in oral squamous cell carcinoma cell lines by blocking p53 nuclear translocation. Oncol. Rep. 40, 179–194. doi:10.3892/or.2018.6451
Lin, Q., Mei, J., Li, Z., Zhang, X., Zhou, L., and Gui, J. (2017). Distinct and cooperative roles of amh and dmrt1 in self-renewal and differentiation of male germ cells in zebrafish. Genetics 207, 1007–1022. doi:10.1534/genetics.117.300274
Liu, J., Cui, J. Y., Lu, Y. F., Christopher Corton, J., and Klaassen, C. D. (2021). Sex-age-and race/ethnicity-dependent variations in drug-processing and NRF2-regulated genes in human livers. Drug Metab. Dispos. 49, 111–119. doi:10.1124/dmd.120.000181
Liu, J., Sun, C. M., Zhang, C. L., Wang, X., and Li, J. Y. (2013). Location and characterization of GAPDS in male reproduction. Urol. Int. 90, 449–454. doi:10.1159/000345629
Liu, K. C., Lin, S. W., and Ge, W. (2011). Differential regulation of gonadotropin receptors (fshr and lhcgr) by estradiol in the zebrafish ovary involves nuclear estrogen receptors that are likely located on the plasma membrane. Endocrinology 152, 4418–4430. doi:10.1210/en.2011-1065
Love, M. I., Huber, W., and Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 550–621. doi:10.1186/s13059-014-0550-8
Lyubimov, A. V., Chang, T. K. H., and Waxman, D. J. (2012). Sex differences in drug metabolism. Encycl. Drug Metab. Interact. 6, 1–19. doi:10.1002/9780470921920.edm004
Maanen, M. Van, Fournier, P. A., Palmer, T. N., and Abraham, L. J. (1999). Characterization of the human glycogenin-1 gene: Identification of a muscle-specific regulatory domain. Gene 234, 217–226. doi:10.1016/s0378-1119(99)00211-5
Magee, J. A., Chang, L. W., Stormo, G. D., and Milbrandt, J. (2006). Direct, androgen receptor-mediated regulation of the FKBP5 gene via a distal enhancer element. Endocrinology 147, 590–598. doi:10.1210/en.2005-1001
Major, A. T., Whiley, P. A. F., and Loveland, K. L. (2011). Expression of nucleocytoplasmic transport machinery: Clues to regulation of spermatogenic development. Biochim. Biophys. Acta - Mol. Cell Res. 1813, 1668–1688. Available at:. doi:10.1016/j.bbamcr.2011.03.008
Mark, M., Ghyselinck, N. B., and Chambon, P. (2006). Function of retinoid nuclear receptors: Lessons from genetic and pharmacological dissections of the retinoic acid signaling pathway during mouse embryogenesis. Annu. Rev. Pharmacol. Toxicol. 46, 451–480. doi:10.1146/annurev.pharmtox.46.120604.141156
Martín-Millán, M., and Castañeda, S. (2013). Estrogens, osteoarthritis and inflammation. Jt. Bone Spine 80, 368–373. doi:10.1016/j.jbspin.2012.11.008
Matson, C. K., Murphy, M. W., Sarver, A. L., Griswold, M. D., Bardwell, V. J., and Zarkower, D. (2011). DMRT1 prevents female reprogramming in the postnatal mammalian testis. Nature 476, 101–104. doi:10.1038/nature10239
McCarthy, M. M., Woolley, C. S., and Arnold, A. P. (2017). Incorporating sex as a biological variable in neuroscience: What do we gain? Nat. Rev. Neurosci. 18, 707–708. Available at:. doi:10.1038/nrn.2017.137
McMenamin, S. K., Minchin, J. E. N., Gordon, T. N., Rawls, J. F., and Parichy, D. M. (2013). Dwarfism and increased adiposity in the gh1 mutant zebrafish vizzini. Endocrinology 154, 1476–1487. doi:10.1210/en.2012-1734
Mohammed, A. (2011). Why are Early Life Stages of Aquatic Organisms more Sensitive to Toxicants than Adults?, New Insights into Toxicity and Drug Testing. University of the West Indies, St. Augustine. Available at: https://www.intechopen.com/books/advanced-biometric-technologies/liveness-detection-in-biometrics.
Moyer, A. J., Gardiner, K., and Reeves, R. H. (2021). All creatures great and small: New approaches for understanding down syndrome genetics. Trends Genet. 37, 444–459. Available at:. doi:10.1016/j.tig.2020.09.017
Nagabhushana, A., and Mishra, R. K. (2016). Finding clues to the riddle of sex determination in zebrafish. J. Biosci. 41, 145–155. doi:10.1007/s12038-016-9593-1
Nagahama, Y., and Yamashita, M. (2008). Regulation of oocyte maturation in fish. Dev. Growth Differ. 50, 195–219. doi:10.1111/j.1440-169X.2008.01019.x
Nagl, F., Schönhofer, K., Seidler, B., Mages, J., Allescher, H. D., Schmid, R. M., et al. (2009). Retinoic acid-induced nNOS expression depends on a novel PI3K/Akt/DAX1 pathway in human TGW-nu-I neuroblastoma cells. Am. J. Physiol. - Cell Physiol. 297, 1146–1156. doi:10.1152/ajpcell.00034.2009
Nair, S., Xu, C., Shen, G., Hebbar, V., Gopalakrishnan, A., Hu, R., et al. (2007). Toxicogenomics of endoplasmic reticulum stress inducer tunicamycin in the small intestine and liver of Nrf2 knockout and C57BL/6J mice. Toxicol. Lett. 168, 21–39. doi:10.1016/j.toxlet.2006.10.012
Nau, H., Tzimas, G., Mondry, M., Plum, C., and Spohr, H. L. (1995). Antiepileptic drugs alter endogenous retinoid concentrations: A possible mechanism of teratogenesis of anticonvulsant therapy. Life Sci. 57, 53–60. doi:10.1016/0024-3205(95)00242-x
Nicolson, T. J., Mellor, H. R., and Roberts, R. R. A. (2010). Gender differences in drug toxicity. Trends Pharmacol. Sci. 31, 108–114. doi:10.1016/j.tips.2009.12.001
Nishimura, F., Naruishi, H., Naruishi, K., Yamada, T., Sasaki, J., Peters, C., et al. (2002). Cathepsin-L, a key molecule in the pathogenesis of drug-induced and I-cell disease-mediated gingival overgrowth: A study with cathepsin-L-deficient mice. Am. J. Pathol. 161, 2047–2052. Available at:. doi:10.1016/S0002-9440(10)64483-5
Obrochta, K. M., Krois, C. R., Campos, B., and Napoli, J. L. (2015). Insulin regulates retinol dehydrogenase expression and all-trans-retinoic acid biosynthesis through FoxO1. J. Biol. Chem. 290, 7259–7268. doi:10.1074/jbc.M114.609313
OECD (2010). Oecd guidelines 210 for the testing of chemicals. OECD 24, 1. doi:10.1787/9789264203785-en
OECD (2013). Oecd guidelines 236 for the testing of chemicals. OECD. Available at: https://www.oecd-ilibrary.org/environment/oecd-guidelines-for-the-testing-of-chemicals-section-2-effects-on-biotic-systems_20745761.
Oh, C., Aho, H., Shamsadin, R., Nayernia, K., Müller, C., Sancken, U., et al. (2003). Characterization, expression pattern and chromosomal localization of the spermatogenesis associated 6 gene (Spata6). Mol. Hum. Reprod. 9, 321–330. doi:10.1093/molehr/gag047
Onichtchouk, D., Aduroja, K., Belting, H. G., Gnügge, L., and Driever, W. (2003). Transgene driving GFP expression from the promoter of the zona pellucida gene zpc is expressed in oocytes and provides an early marker for gonad differentiation in zebrafish. Dev. Dyn. 228, 393–404. doi:10.1002/dvdy.10392
Örn, S., Holbech, H., Madsen, T. H., Norrgren, L., and Petersen, G. I. (2003). Gonad development and vitellogenin production in zebrafish (Danio rerio) exposed to ethinylestradiol and methyltestosterone. Aquat. Toxicol. 65, 397–411. doi:10.1016/s0166-445x(03)00177-2
Oron, G., Fisch, B., Ao, A., Zhang, X. Y., Farhi, J., Haroush, A. Ben, et al. (2010). Expression of growth-differentiating factor 9 and its type 1 receptor in human ovaries. Reprod. Biomed. Online 21, 109–117. doi:10.1016/j.rbmo.2010.03.011
Osman, M. (2003). Therapeutic implications of sex differences in asthma and atopy. Arch. Dis. Child. 88, 587–590. doi:10.1136/adc.88.7.587
Otis, J. P., Zeituni, E. M., Thierer, J. H., Anderson, J. L., Brown, A. C., Boehm, E. D., et al. (2015). Zebrafish as a model for apolipoprotein biology: Comprehensive expression analysis and a role for ApoA-IV in regulating food intake. DMM Dis. Model. Mech. 8, 295–309. doi:10.1242/dmm.018754
Otsuka, F., Yao, Z., Lee, T. H., Yamamoto, S., Erickson, G. F., and Shimasaki, S. (2000). Bone morphogenetic protein-15: Identification of target cells and biological functions. J. Biol. Chem. 275, 39523–39528. doi:10.1074/jbc.M007428200
Paulini, F., and Melo, E. O. (2011). The role of oocyte-secreted factors GDF9 and BMP15 in follicular development and oogenesis. Reprod. Domest. Anim. 46, 354–361. doi:10.1111/j.1439-0531.2010.01739.x
Pertea, M., Pertea, G. M., Antonescu, C. M., Chang, T. C., Mendell, J. T., and Salzberg, S. L. (2015). StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat. Biotechnol. 33, 290–295. doi:10.1038/nbt.3122
Plank-Bazinet, J. L., Sampson, A., Miller, L. R., Fadiran, E. O., Kallgren, D., Agarwal, R. K., et al. (2016). The science of sex and gender in human health: Online courses to create a foundation for sex and gender accountability in biomedical research and treatment. Biol. Sex. Differ. 7, 47–49. Available at:. doi:10.1186/s13293-016-0100-z
Pradhan, A., and Olsson, P. (2015). Inhibition of retinoic acid synthesis disrupts spermatogenesis and fecundity in zebrafish. Gen. Comp. Endocrinol. 217–218, 81–91. doi:10.1016/j.ygcen.2015.02.002
Priyadarsini, S., Mohanty, S., Mukherjee, S., Basu, S., and Mishra, M. (2018). Graphene and graphene oxide as nanomaterials for medicine and biology application. J. Nanostructure Chem. 8, 123–137. Available at:. doi:10.1007/s40097-018-0265-6
Ridker, P. M., Cook, N. R., I-Min, L., Gordon, D., and Buring, J. E. (2015). A randomized trial of low-dose aspirin in the primary prevention of cardiovascular disease in WomenMichael gaziano. N. Engl. J. Med. 352, 687–696. doi:10.1056/NEJMoa050613
Rodríguez-Marí, A., Cañestro, C., BreMiller, R. A., Catchen, J. M., Yan, Y. L., and Postlethwait, J. H. (2013). Retinoic acid metabolic genes, meiosis, and gonadal sex differentiation in zebrafish. PLoS One 8, 739511–e74019. doi:10.1371/journal.pone.0073951
Rodríguez-Marí, A., Cañestro, C., BreMiller, R. A., Nguyen-Johnson, A., Asakawa, K., Kawakami, K., et al. (2010). Sex reversal in zebrafish fancl mutants is caused by Tp53-mediated germ cell apoptosis. PLoS Genet. 6, 10010344–e1001114. doi:10.1371/journal.pgen.1001034
Rudzińska, M., Parodi, A., Soond, S. M., Vinarov, A. Z., Korolev, D. O., Morozov, A. O., et al. (2019). The role of cysteine cathepsins in cancer progression and drug resistance. Int. J. Mol. Sci. 20, 3602–3621. doi:10.3390/ijms20143602
Sandell, L. L., Sanderson, B. W., Moiseyev, G., Johnson, T., Mushegian, A., Young, K., et al. (2007). RDH10 is essential for synthesis of embryonic retinoic acid and is required for limb, craniofacial, and organ development. Genes Dev. 21, 1113–1124. doi:10.1101/gad.1533407
Santos, E. M., Workman, V. L., Paull, G. C., Filby, A. L., Look, K. J. W. Van, Kille, P., et al. (2007). Molecular basis of sex and reproductive status in breeding zebrafish. Physiol. Genomics 30, 111–122. doi:10.1152/physiolgenomics.00284.2006
Schieck, M., Schouten, J. P., Michel, S., Suttner, K., Toncheva, A. A., Gaertner, V. D., et al. (2016). Doublesex and mab-3 related transcription factor 1 (DMRT1) is a sex-specific genetic determinant of childhood-onset asthma and is expressed in testis and macrophages. J. Allergy Clin. Immunol. 138, 421–431. doi:10.1016/j.jaci.2015.12.1305
Schlessinger, D., Forabosco, A., Uda, M., Pelosi, E., and Crisponi, L. (2011). Determination and stability of gonadal sex. J. Androl. 31, 16–25. doi:10.2164/jandrol.109.008201
Schwarze, K., Buchanan, J., Fermont, J. M., Dreau, H., Tilley, M. W., Taylor, J. M., et al. (2020). The complete costs of genome sequencing: A microcosting study in cancer and rare diseases from a single center in the United Kingdom. Genet. Med. 22, 85–94. Available at:. doi:10.1038/s41436-019-0618-7
Sedjo, R. L., Ranger-Moore, J., Foote, J., Craft, N. E., Alberts, D. S., Xu, M. J., et al. (2004). Circulating endogenous retinoic acid concentrations among participants enrolled in a randomized placebo-controlled clinical trial of retinyl palmitate. Cancer Epidemiol. Biomarkers Prev. 13, 1687–1692. doi:10.1158/1055-9965.1687.13.11
Segner, H. (2009). Zebrafish (Danio rerio) as a model organism for investigating endocrine disruption. Comp. Biochem. Physiol. Part C 149, 187–195. doi:10.1016/j.cbpc.2008.10.099
Selman, K., Wallace, R. A., Sarka, A., and Qi, X. (1993). Stages of oocyte development in the zebrafish, Brachydanio rerio. J. Morphol. 218, 203–224. doi:10.1002/jmor.1052180209
Shanthanagouda, A. H., Patil, J. G., and Nugegoda, D. (2013). Effects of exposure to oestrogenic compounds on aromatase gene expression are gender dependent in the rainbowfish, Melanotaenia fluviatilis. Comp. Biochem. Physiol. - C Toxicol. Pharmacol. 157, 162–171. doi:10.1016/j.cbpc.2012.11.004
Silva, J. R. V., Van Den Hurk, R., Van Tol, H. T. A., Roelen, B. A. J., and Figueiredo, J. R. (2005). Expression of growth differentiation factor 9 (GDF9), bone morphogenetic protein 15 (BMP15), and BMP receptors in the ovaries of goats. Mol. Reprod. Dev. 70, 11–19. doi:10.1002/mrd.20127
Soldin, O. P., and Mattison, D. R. (2009). Sex differences in pharmacokinetics and pharmacodynamics. Clin. Pharmacokinet. 48, 143–157. doi:10.2165/00003088-200948030-00001
Song, Y., Chen, J., Tao, B., Luo, D., Zhu, Z., and Hu, W. (2020). Kisspeptin2 regulates hormone expression in female zebrafish (Danio rerio) pituitary. Mol. Cell. Endocrinol. 513, 110858–8. doi:10.1016/j.mce.2020.110858
Soyal, S. M., Amleh, A., and Dean, J. (2000). FIGalpha, a germ cell-specific transcription factor required for ovarian follicle formation. Development 127, 4645–4654. doi:10.1242/dev.127.21.4645
Strauss, J. F., and Williams, C. J. (2019). Ovarian life cycle. Eighth edi. Netherlands: Elsevier Inc.
Sujit, K. M., Singh, V., Trivedi, S., Singh, K., Gupta, G., and Rajender, S. (2020). Increased DNA methylation in the spermatogenesis-associated (SPATA) genes correlates with infertility. Andrology 8, 602–609. doi:10.1111/andr.12742
Takahashi, K., Kotani, T., Katsu, Y., and Yamashita, M. (2014). Possible involvement of insulin-like growth factor 2 mRNA-binding protein 3 in zebrafish oocyte maturation as a novel cyclin B1 mRNA-binding protein that represses the translation in immature oocytes. Biochem. Biophys. Res. Commun. 448, 22–27. doi:10.1016/j.bbrc.2014.04.020
Takei, N., Nakamura, T., Kawamura, S., Takada, Y., Satoh, Y., Kimura, A. P., et al. (2018). High-sensitivity and high-resolution in situ hybridization of coding and long non-coding RNAs in vertebrate ovaries and testes. Biol. Proced. Online 20, 6–14. doi:10.1186/s12575-018-0071-z
Tiedke, J., Gerlach, F., Mitz, S. A., Hankeln, T., and Burmester, T. (2011). Ontogeny of globin expression in zebrafish (Danio rerio). J. Comp. Physiol. B Biochem. Syst. Environ. Physiol. 181, 1011–1021. doi:10.1007/s00360-011-0588-9
Tingaud-Sequeira, A., Carnevali, O., and Cerdà, J. (2011). Cathepsin B differential expression and enzyme processing and activity during Fundulus heteroclitus embryogenesis. Comp. Biochem. Physiol. - A Mol. Integr. Physiol. 158, 221–228. doi:10.1016/j.cbpa.2010.11.002
Tingaud-Sequeira, A., and Cerdà, J. (2007). Phylogenetic relationships and gene expression pattern of three different cathepsin L (Ctsl) isoforms in zebrafish: Ctsla is the putative yolk processing enzyme. Gene 386, 98–106. doi:10.1016/j.gene.2006.08.018
Tong, S. K., Hsu, H. J., and Chung, B. C. (2010a). Zebrafish monosex population reveals female dominance in sex determination and earliest events of gonad differentiation. Dev. Biol. 344, 849–856. Available at:. doi:10.1016/j.ydbio.2010.05.515
Tong, S. K., Hsu, H. J., and Chung, B. C. (2010b). Zebrafish monosex population reveals female dominance in sex determination and earliest events of gonad differentiation. Dev. Biol. 344, 849–856. doi:10.1016/j.ydbio.2010.05.515
Uchida, D., Yamashita, M., Kitano, T., and Iguchi, T. (2002). Oocyte apoptosis during the transition from ovary-like tissue to testes during sex differentiation of juvenile zebrafish. J. Exp. Biol. 205, 711–718. doi:10.1242/jeb.205.6.711
Uhlenhaut, N. H., Jakob, S., Anlag, K., Eisenberger, T., Sekido, R., Kress, J., et al. (2009). Somatic sex reprogramming of adult ovaries to testes by FOXL2 ablation. Cell 139, 1130–1142. doi:10.1016/j.cell.2009.11.021
U.S. Government Accountability Office (2001). Drug safety: Most drugs withdrawn in recent years had greater health risks for women. Drug Saf. 343, 1–8. Available at: http://www.gao.gov/new.items/d01286r.pdf.
Valdivieso, A., Ribas, L., and Piferrer, F. (2019). Ovarian transcriptomic signatures of zebrafish females resistant to different environmental perturbations. J. Exp. Zool. Part B Mol. Dev. Evol. 332, 55–68. doi:10.1002/jez.b.22848
van boxtel, A. L., Pieterse, B., Cenijn, P., Kamstra, J. H., Brouwer, A., van wieringen, W., et al. (2010). Dithiocarbamates induce craniofacial abnormalities and downregulate sox9a during zebrafish development. Toxicol. Sci. 117, 209–217. doi:10.1093/toxsci/kfq169
Vasconcelos, V., Azevedo, J., Silva, M., and Ramos, V. (2010). Effects of marine toxins on the reproduction and early stages development of aquatic organisms. Mar. Drugs 8, 59–79. doi:10.3390/md8010059
Verdin, H., and De Baere, E. (2012). FOXL2 impairment in human disease. Horm. Res. Paediatr. 77, 2–11. doi:10.1159/000335236
Volkova, K., Reyhanian Caspillo, N., Porseryd, T., Hallgren, S., Dinnétz, P., and Porsch-Hällström, I. (2015). Developmental exposure of zebrafish (Danio rerio) to 17α-ethinylestradiol affects non-reproductive behavior and fertility as adults, and increases anxiety in unexposed progeny. Horm. Behav. 73, 30–38. doi:10.1016/j.yhbeh.2015.05.014
Wang, D. S., Jiao, B., Hu, C., Huang, X., Liu, Z., and Cheng, C. H. K. (2008a). Discovery of a gonad-specific IGF subtype in teleost. Biochem. Biophys. Res. Commun. 367, 336–341. doi:10.1016/j.bbrc.2007.12.136
Wang, H., Tan, J. T. T., Emelyanov, A., Korzh, V., and Gong, Z. (2005). Hepatic and extrahepatic expression of vitellogenin genes in the zebrafish, Danio rerio. Gene 356, 91–100. doi:10.1016/j.gene.2005.03.041
Wang, L., Sewell, W. F., Kim, S. D., Shin, J. T., MacRae, C. A., Zon, L. I., et al. (2008b). Eya4 regulation of Na+/K+-ATPase is required for sensory system development in zebrafish. Development 135, 3425–3434. doi:10.1242/dev.012237
Wang, X. G., Bartfai, R., Sleptsova-Freidrich, I., and Orban, L. (2007). The timing and extent of ‘ juvenile ovary ’ phase are highly variable during zebrafish testis differentiation. J. Fish. Biol. 70, 33–44. doi:10.1111/j.1095-8649.2007.01363.x
Wei, S., Lai, L., Yang, J., and Zhuandi, G. (2018). Expression levels of follicle-stimulating hormone receptor and implication in diagnostic and therapeutic strategy of ovarian cancer. Oncol. Res. Treat. 41, 651–654. doi:10.1159/000490810
Wei, S., Shen, X., Gong, Z., Deng, Y., Lai, L., and Liang, H. (2017). FSHR and LHR expression and signaling as well as maturation and apoptosis of cumulus-oocyte complexes following treatment with FSH receptor binding inhibitor in sheep. Cell. Physiol. biochem. 43, 660–669. doi:10.1159/000480650
Welch, J. E., Brown, P. R., Brien, D. A. O., and Eddy, E. M. (1995). Genomic organization of a mouse glyceraldehyde 3-phosphate dehydrogenase gene (Gapd-s) expressed in post-meiotic spermatogenic cells. Dev. Genet. 16, 179–189. doi:10.1002/dvg.1020160210
Wen, C., Zhang, Z., Ma, W., Xu, M., Wen, Z., and Peng, J. (2005). Genome-wide identification of female-enriched genes in zebrafish. Dev. Dyn. 232, 171–179. doi:10.1002/dvdy.20210
Wu, G. C., Tomy, S., and Chang, C. F. (2008). The expression of nr0b1 and nr5a4 during gonad development and sex change in protandrous black porgy fish, Acanthopagrus schlegeli. Biol. Reprod. 78, 200–210. doi:10.1095/biolreprod.107.062612
Wu, J., Miao, C., Lv, X., Zhang, Y., Li, Y., and Wang, D. (2019). Estrogen regulates forkhead transcription factor 2 to promote apoptosis of human ovarian granulosa-like tumor cells. J. Steroid Biochem. Mol. Biol. 194, 105418–105511. doi:10.1016/j.jsbmb.2019.105418
Xia, H., Zhong, C., Wu, X., Chen, J., Tao, B., Xia, X., et al. (2018). Mettl3 mutation disrupts gamete maturation and reduces fertility in zebrafish. Genetics 208, 729–743. doi:10.1534/genetics.117.300574
Xiong, S., Mhawech-Fauceglia, P., Tsao-Wei, D., Roman, L., Gaur, R. K., Epstein, A. L., et al. (2019). Expression of the luteinizing hormone receptor (LHR) in ovarian cancer. BMC Cancer 19, 1114–1118. doi:10.1186/s12885-019-6153-8
Xu, X., Zhang, X., Xing, H., Liu, Z., Jia, J., Jin, C., et al. (2019). Importin-4 functions as a driving force in human primary gastric cancer. J. Cell. Biochem. 120, 12638–12646. doi:10.1002/jcb.28530
Yan, Y. L., Desvignes, T., Bremiller, R., Wilson, C., Dillon, D., High, S., et al. (2017). Gonadal soma controls ovarian follicle proliferation through Gsdf in zebrafish. Dev. Dyn. 246, 925–945. doi:10.1002/dvdy.24579
Yang, L., Li, Y., Hong, H., Chang, C. W., Guo, L. W., Lyn-Cook, B., et al. (2012). Sex differences in the expression of drug-metabolizing and transporter genes in human liver. J. Drug Metab. Toxicol. 3, 1000119–1000120. doi:10.4172/2157-7609.1000119
Yang, Y. J., Wang, Y., Li, Z., Zhou, L., and Gui, J. F. (2017). Sequential, divergent, and cooperative requirements of foxl2a and foxl2b in ovary development and maintenance of zebrafish. Genetics 205, 1551–1572. doi:10.1534/genetics.116.199133
Yasuda, K., Kotani, T., Ota, R., and Yamashita, M. (2010). Transgenic zebrafish reveals novel mechanisms of translational control of cyclin B1 mRNA in oocytes. Dev. Biol. 348, 76–86. Available at:. doi:10.1016/j.ydbio.2010.09.015
Yasuda, K., Kotani, T., and Yamashita, M. (2013). A cis-acting element in the coding region of cyclin B1 mRNA couples subcellular localization to translational timing. Dev. Biol. 382, 517–529. Available at:. doi:10.1016/j.ydbio.2013.05.014
Yi, X., Yu, J., Ma, C., Dong, G., Shi, W., Li, H., et al. (2019). The effector of Hippo signaling, Taz, is required for formation of the micropyle and fertilization in zebrafish. PLoS Genet. 15, e1007408–e1007418. doi:10.1371/journal.pgen.1007408
Yin, Y., Tang, H., Liu, Y., Chen, Y., Li, G., Liu, X., et al. (2017). Targeted disruption of aromatase reveals dual functions of cyp19a1a during sex differentiation in zebrafish. Endocr. Soc. 1, 3030–3041. doi:10.1210/en.2016-1865
Yu, M., Wang, J., Liu, W., Qin, J., Zhou, Q., Wang, Y., et al. (2014). Effects of tamoxifen on the sex determination gene and the activation of sex reversal in the developing gonad of mice. Toxicology 321, 89–95. Available at:. doi:10.1016/j.tox.2014.04.006
Yuan, S., Stratton, C. J., Bao, J., Zheng, H., Bhetwal, B. P., Yanagimachi, R., et al. (2015). SPATA6 is required for normal assembly of the sperm connecting piece and tight head-tail conjunction. PNAS 112, 430–439. doi:10.1073/pnas.1424648112
Zhai, G., Shu, T., Xia, Y., Lu, Y., Shang, G., Jin, X., et al. (2018). Characterization of sexual trait development in cyp17a1 -deficient zebrafish. Endocrinology 159, 3549–3562. doi:10.1210/en.2018-00551
Zhang, X., Gao, L., Yang, K., Tian, H., Wang, W., and Ru, S. (2013). Monocrotophos pesticide modulates the expression of sexual differentiation genes and causes phenotypic feminization in zebrafish (Danio rerio). Comp. Biochem. Physiol. - C Toxicol. Pharmacol. 157, 33–40. doi:10.1016/j.cbpc.2012.09.004
Zhang, Z., Zhu, B., and Ge, W. (2015). Genetic analysis of zebrafish gonadotropin (FSH and LH) functions by TALEN-mediated gene disruption. Mol. Endocrinol. 29, 76–98. doi:10.1210/me.2014-1256
Zhao, Y., Yang, Z., Phelan, J. K., Wheeler, D. A., Lin, S., and McCabe, E. R. B. (2006). Zebrafish dax1 is required for development of the interrenal organ, the adrenal cortex equivalent. Mol. Endocrinol. 20, 2630–2640. doi:10.1210/me.2005-0445
Zheng, M., Lu, J., and Zhao, D. (2018). Toxicity and transcriptome sequencing (RNA-seq) analyses of adult zebrafish in response to exposure carboxymethyl cellulose stabilized iron sulfide nanoparticles. Sci. Rep. 8, 8083–8111. Available at:. doi:10.1038/s41598-018-26499-x
Zheng, X., Chou, P. M., Mirkin, B. L., and Rebbaa, A. (2004). Senescence-initiated reversal of drug resistance: Specific role of cathepsin L. Cancer Res. 64, 1773–1780. doi:10.1158/0008-5472.can-03-0820
Zheng, X., Chu, F., Chou, P. M., Gallati, C., Dier, U., Mirkin, B. L., et al. (2009). Cathepsin L inhibition suppresses drug resistance in vitro and in vivo: A putative mechanism. Am. J. Physiol. - Cell Physiol. 296, 65–74. doi:10.1152/ajpcell.00082.2008
Zhong, Y., Wang, Y., Huang, J., Xu, X., Pan, W., Gao, S., et al. (2019). Association of hCG and LHCGR expression patterns with clinicopathological parameters in ovarian cancer. Pathol. Res. Pract. 215, 748–754. Available at:. doi:10.1016/j.prp.2019.01.001
Zhou, Y., Liu, F., Xu, Q., Yang, B., Li, X., Jiang, S., et al. (2020). Inhibiting Importin 4-mediated nuclear import of CEBPD enhances chemosensitivity by repression of PRKDC-driven DNA damage repair in cervical cancer. Oncogene 39, 5633–5648. doi:10.1038/s41388-020-1384-3
Zon, L. I., and Peterson, R. T. (2005). In vivo drug discovery in the zebrafish. Nat. Rev. Drug Discov. 4, 35–44. doi:10.1038/nrd1606
Keywords: sex-variation, gene expression, zebrafish, drug influence, human health
Citation: King AC and Zenker AK (2023) Sex blind: bridging the gap between drug exposure and sex-related gene expression in Danio rerio using next-generation sequencing (NGS) data and a literature review to find the missing links in pharmaceutical and environmental toxicology studies. Front. Toxicol. 5:1187302. doi: 10.3389/ftox.2023.1187302
Received: 15 March 2023; Accepted: 01 June 2023;
Published: 16 June 2023.
Edited by:
Steven Jan Van Cruchten, University of Antwerp, BelgiumReviewed by:
Lucia Vergauwen, University of Antwerp, BelgiumDavid Hala, Texas A&M University at Galveston, United States
Copyright © 2023 King and Zenker. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Armin K. Zenker, Armin.Zenker@fhnw.ch