Phylogenomic and expression analysis of Colossoma macropomum cyp19a1a and cyp19a1b and their non-classical role in tambaqui sex differentiation

Abstract

The genes coding for Cytochrome P450 aromatase (cyp19a1a and cyp19a1b) and estrogen (E₂) receptors (esr1, esr2a and esr2b) play a conserved role in ovarian differentiation and development among teleosts. Classically, the "gonad form" of aromatase, coded by the cyp19a1a, is responsible for the ovarian differentiation in genetic females via ligation and activation of the Esr, which mediates the endocrine and exocrine signaling to allow or block the establishment of the feminine phenotype. However, in neotropical species, studies on the molecular and endocrine processes involved in gonad differentiation as well as on the effects of sex modulators are recent and scarce. In this study, we combined in silico analysis, real-time quantitative PCR (qPCR) assay and quantification of E₂ plasma levels of differentiating tambaqui (Colossoma macropomum) to unveil the roles of the paralogs cypa19a1a and cyp19a1b during sex differentiation. Although the synteny of each gene is very conserved among characids, the genomic environment displays striking differences in comparison to model teleost species, with many rearrangements in cyp19a1a and cyp19a1b adjacencies and transposable element traces in both regulatory regions. The high dissimilarity (DI) of SF-1 binding motifs in cyp19a1a (DI = 10.06 to 14.90 %) and cyp19a1b (DI = 8.41 to 13.50 %) regulatory region, respectively, may reflect in an alternative pathway in tambaqui. Indeed, while low transcription of cyp19a1a was detected prior to sex differentiation, the expression of cyp19a1b and esr2a presented a large variation at this phase, which could be associated with sex-specific differential expression. Histological analysis revealed that anti-estradiol treatments did not affect gonadal sex ratios, although Fadrozole (50 mg kg^-1 of food) reduced E₂ plasma levels (p < 0,005) as well cyp19a1a transcription; and tamoxifen (200 mg kg^-1 of food) down regulated both cyp19a1a and cyp19a1b but did not influence E₂ levels. Altogether, our results bring into light new insights about the evolutionary fate of cyp19a1 paralogs in neotropical fish, which may have generated uncommon roles for the gonadal and brain forms of cyp19a1 genes and the unexpected lack of effect of endocrine disruptors on tambaqui sexual differentiation.

Keywords: Cyp19a1; Neotropical fish; Promoter Structure; Sex development; Synteny; aromatase.