Next Article in Journal
Effects of Fatty-Type and Lean-Type on Growth Performance and Lipid Droplet Metabolism in Pekin Ducks
Previous Article in Journal
Prepubertal Use of Long-Term GnRH Agonists in Dogs: Current Knowledge and Recommendations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genome-Wide Identification and Characterization of the BRD Family in Nile Tilapia (Oreochromis niloticus)

Key Laboratory of Freshwater Fish Reproduction and Development (Ministry of Education), Key Laboratory of Aquatic Science of Chongqing, School of Life Sciences, Southwest University, Chongqing 400715, China
*
Authors to whom correspondence should be addressed.
Animals 2022, 12(17), 2266; https://doi.org/10.3390/ani12172266
Submission received: 14 July 2022 / Revised: 22 August 2022 / Accepted: 31 August 2022 / Published: 1 September 2022
(This article belongs to the Section Animal Genetics and Genomics)

Abstract

:

Simple Summary

Nile tilapia is a good model for genome-wide identification and examination of the expression and role of gene families. In this study, we identified 54 bromodomain genes (BRDs) divided into eight subfamilies in Nile tilapia. Phylogenetic analysis revealed a high conservation of the BRDs family in vertebrates, with BRDs expansion due to fish-specific duplications. Most of the BRDs displayed sexually dimorphic expression in the gonads at 90 and 180 dah (days after hatching), including 21 testis-dominated genes (brdt, brd4a and brd2b, etc.), and 9 ovary-dominated genes (brd3b, brd2a and kat2a, etc.). Male fish treated with JQ1 (BET subfamily inhibitor) displayed abnormal spermatogenesis. The numbers of germ cells were reduced and the expression of steroidogenic enzyme genes was downregulated, while the expression of apoptosis-promoting genes was elevated in the testes of treated fish.

Abstract

The bromodomain (BRD) proteins specifically recognize the N-acetyllysine motifs, which is a key event in the reading process of epigenetic marks. BRDs are evolutionarily highly conserved. Over recent years, BRDs attracted great interest because of their important roles in biological processes. However, the genome-wide identification of this family was not carried out in many animal groups, in particular, in teleosts. Moreover, the expression patterns were not reported for any of the members in this family, and the role of the BRD family was not extensively studied in fish reproduction. In this study, we identified 16 to 120 BRD genes in 24 representative species. BRDs expanded significantly in vertebrates. Phylogenetic analysis showed that the BRD family was divided into eight subfamilies (I–VIII). Transcriptome analysis showed that BRDs in Nile tilapia (Oreochromis niloticus) exhibited different expression patterns in different tissues, suggesting that these genes may play different roles in growth and development. Gonadal transcriptome analysis showed that most of the BRDs display sexually dimorphic expression in the gonads at 90 and 180 dah (days after hatching), including 21 testis-dominated genes (brdt, brd4a and brd2b, etc.), and nine ovary-dominated genes (brd3b, brd2a and kat2a, etc.). Consistent with transcriptomic data, the results of qRT-PCR and fluorescence in situ hybridization showed that brdt expression was higher in the testis than in the ovary, suggesting its critical role in the spermatogenesis of the tilapia. Male fish treated with JQ1 (BET subfamily inhibitor) displayed abnormal spermatogenesis. The numbers of germ cells were reduced, and the expression of steroidogenic enzyme genes was downregulated, while the expression of apoptosis-promoting genes was elevated in the testis tissue of treated fish. Our data provide insights into the evolution and expression of BRD genes, which is helpful for understanding their critical roles in sex differentiation and gonadal development in teleosts.

1. Introduction

Bromodomain (BRD) proteins include a protein domain with approximately 110 amino acids [1,2], comprising four α-helices (αZ, αA, αB, αC) linked by ZA and BC loops. Recognizing the acetylated lysine binding site [3,4], BRD proteins play different roles in biological processes through various mechanisms. For example, these proteins can act as a scaffold for larger protein complexes to facilitate their assembly, function as transcription factors and transcriptional co-regulators [5], and engage in various chromatin modifications [6,7]. Bromodomain was first identified in the fruit fly [8], and later in humans [3,5] and the buffalo [9]. However, to date, no study was systematically performed to analyze the evolution of the BRD family.
BRDs are ubiquitously expressed in the brain, testis, ovary, liver, spleen, lung, kidney and other tissues, with varying expression levels. For example, brd2 is expressed in the zebrafish egg, early embryo, and nervous system [10]. Brdt is specifically expressed in murine testis [11]. Brd1 is widely expressed in the mammalian brain and plays several roles in brain development [12]. Brd7 is expressed in human ovarian cells and involved in epithelial ovarian carcinoma [13]. Brd8 is expressed in the lung and coordinates the transcriptional network in human airway epithelial cells [14]. Despite its functional importance, the existing studies mainly focused on specific members of BRDs in humans, mice, and zebrafish. The expression level of all genes in this family and in other vertebrates, especially in non-model teleosts, was not yet studied. Advances in transcriptome sequencing make it possible to study the expression of all BRDs in different tissues at different developmental stages in vertebrates.
The BET subfamily is characterized by an N-terminal dibromodomain [15] and a unique C-terminal extra terminal (ET) domain [16], comprising Brd2, Brd3, Brd4, and the testis-specific protein, Brdt. This subfamily plays a key role in regulating gene transcription through the interaction between bromodomains and acetylated histones in cell multiplication and differentiation [17,18]. The BET subfamily is required for proper spermatogenesis in vertebrates. BET subfamily members were reported to be commonly expressed in the Sertoli cells of mouse testes [19] and involved in the regulation of testis development and spermatogenesis [20,21,22]. BET subfamily members are highly expressed in mature buffalo Sertoli cells, indicating their critical roles in spermatogenesis [19]. In gilthead seabream and European seabass, brdt had higher expression in gonads and relatively weak expression in other tissues [23]. Moreover, inhibition of brdt can induce contraception in male mice [22]. Disruption of brdt led to infertility in male mice [24]. A genome-wide association study indicated that single nucleotide polymorphisms of brdt were associated with oligospermia or azoospermia in humans [25]. However, to date, these studies were limited to the expression and function of individual genes or several members of the BET subfamily. Recently, treatment of BET inhibitor JQ1 (a thieno-triazolo-1,4-diazapine), which competitively binds to the acetyllysine recognition sites of BET family proteins, resulted in decreased numbers of spermatid cells and changes in chromatin accessibility in the mouse [26]. In vivo studies indicated that administration of JQ1 in male mice can lead to defective spermatogenesis and lower motility of spermatozoa [22]. However, the role of the BET subfamily was not extensively studied in fish reproduction.
Nile tilapia is an important aquaculture fish in the world, with the growth rate of males being faster than that of females. Because of the key role of BRDs in regulating gene expression, it is important to identify BRDs and characterize their expression patterns in different tissues in tilapia. The publication of the whole genome sequence, transcriptomes of eight tissues [27], and transcriptomes of gonads at four developmental stages [28], make it a good model for the genome-wide identification and examination of the expression and role of BRDs. In this study, we performed genome-wide identification of the BRD family and analyzed expression patterns. We injected, intraperitoneally, male tilapia with JQ1 to explore the possible function of the BET subfamily in male reproduction. Our study provides a new perspective on the evolution of BRDs and lays a foundation for elucidating their roles during fish spermatogenesis.

2. Materials and Methods

2.1. Fish and Experimental Conditions

The Nile tilapia used in this experiment were cultivated in the Key Laboratory of Freshwater Fish Resources and Reproductive Development of the Ministry of Education (Beibei, Chongqing, China), Southwest University. Fish were raised in a water circulation system at 26 °C, with the natural photoperiod (usually 12 h). Animal experiments were authorized by the Committee of Laboratory Animal Experimentation at Southwest University, China (No. IACUC-20181015-12).

2.2. Identification of BRD Genes in Representative Animals

In order to identify BRD genes in animals, we examined the genomes of 24 representative animals, including Drosophila melanogaster (GCA_000001215.4), Ciona intestinalis (GCA_000224145.2), Branchiostoma belcheri (GCA_001625305.1), Petromyzon marinus (GCA_010993605.1), Callorhinchus milii (GCA_018977255.1), Latimeria chalumnae (GCA_000225785.1), Python bivittatus (GCA_000186305.2), Xenopus tropicalis (GCA_000004195.4), Gallus gallus (GCA_016699485.1), Mus musculus (GCA_000001635.9), Homo sapiens GCA_000001405.29), Lepisosteus oculatus (GCA_000242695.1), Takifugu rubripes (GCA_901000725.2), Micropterus salmoides (GCA_014851395.1), Oreochromis niloticus (GCA_001858045.3), Oryzias latipes (GCA_002234675.1), Poecilia Formosa (GCA_000485575.1), Xiphophorus maculatus (GCA_002775205.2), Ictalurus punctatus (GCA_001660625.2), Pangasianodon hypophthalmus (GCA_009078355.1), Tachysurus fulvidraco (GCA_022655615.1), Danio rerio (GCA_000002035.4), Astyanax mexicanus (GCA_023375975.1), and Cyprinus carpio (GCA_018340385.1). All BRDs were identified by reciprocal best-BLAST against a protein database of these species, using human BRD-protein sequences as the query sequences [29]. The accession numbers of all BRD family-member gene sequences used in this study are listed in Table S1.
The species tree shown in Figure 1 was constructed using the NCBI database Taxonomy Browser (https://www.ncbi.nlm.nih.gov/Taxonomy/CommonTree/wwwcmt.cgi (accessed on 5 October 2021)) [30], and then opened in Evolview (https://www.evolgenius.info//evolview/#login (accessed on 5 October 2021)) for editing.

2.3. Phylogenetic Analysis, Synteny Analysis, and Functional Domains Analysis of BRD Genes

The amino acid sequences of the BRD proteins of Nile tilapia, zebrafish and humans were aligned with Bioedit [31]. The phylogenetic tree was constructed by MEGA (version 7.0, Mega Limited, Auckland, New Zealand) using the neighbor-joining (NJ) method with a bootstrap of 1000 replicates to assess the confidence in the phylogeny [32], then edited with Adobe Illustrator (version 4.0, Adobe Systems Inc., San Jose, CA, USA). Syntenic analysis was conducted by comparing genomic regions that harbor baz2b, crebbp, brd3, brd8, brd2, brd4, brpf3, ep300, brd,1 and kmt2b (representative genes underwent 3R-WGD) in tilapia with those in representative species. The genomic regions surrounding these genes were determined according to the Ensembl, NCBI, and Genomics databases [33] and edited with Adobe Illustrator CS4 software. NCBI Batch CD-Search (https://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi (accessed on 22 November 2021)) was used to search domains on BRD proteins. Gene structure was displayed by TBtools [34] and edited with Adobe Illustrator CS4 software.

2.4. Spatiotemporal Expression Pattern of BRDs in Tilapia

The transcriptomic data from eight tissues (ovary, testis, skeletal muscle, liver, kidney, head kidney, heart, and brain,) of adult tilapia were obtained from the NCBI database (Accession number: PRJNA78915 and SRR1916191) [27]. In our previous study, eight samples from the gonads of female (XX) and male (XY) fish at 5, 30, 90, and 180 dah (days after hatching) were sequenced using Illumina 2000 HiSeq technology [28]. A normalized measure of FPKM (fragments per kilobase million) was used to normalize the expression of BRD genes. BRDs with FPKM ≤ 1 in each sample were defined as background expression, according to a previous study [35]. The sexually dimorphic expression gene in the gonads was identified by fold-change ≥ 2 at cutoffs of p ≤ 0.01 using edgeR.

2.5. Validation of Spatiotemporal Expression Pattern of BRD Genes by qRT-PCR and Fluorescence In Situ Hybridization

The expression of brdt and brd2b was detected by quantitative real time PCR (qRT-PCR). In total, 12 tissues were isolated from Nile tilapia at 180 dah. Total RNA was extracted and a 1.0 µg total RNA template was reversely transcribed into cDNA following the method described previously [36]. The qRT-PCR was performed following the SYBR Green I Master Mix (Takara, Dalian, China) protocol. The relative abundance of mRNA transcripts was calculated by the formula R = 2 Δ Δ Ct [37]. gapdh was used as an internal control. Data (expressed as the mean ± SD) from all experiments were analyzed using GraphPad Prism 5.0. A one-way ANOVA was performed, followed by the Tukey test to determine the significance (* p < 0.05; ** p < 0.01; *** p < 0.001; ns, no significant changes observed).
Fluorescence in situ hybridization (FISH) was used to detect the expressed cell types of brdt in tilapia gonads at 180 dah. The cds of brdt was amplified with specific primers (Table S2), and the amplified fragments were ligated into pGEM-T Easy Vector. The gonads of male and female tilapia were sampled at 180 dah and fixed in 4% paraformaldehyde in phosphate-buffered saline (PBS), then embedded, and serial paraffin sections were cut at 5 µm thickness. The sections were dewaxed, rehydrated, and digested with proteinase K (4 µg/mL; Roche, Mannheim, Germany) for 15 min (37 °C), and then hybridized with DIG-labeled RNA probes for at least 16 h. After hybridization, the sections were washed with SSC for 10 min, DIG I (0.1 mol/L maleic acid, 0.15 mol/L NaCl, ph7.5) for 5 min, blocked with DIG II (1% bovine serum albumin in DIG I) for 30 min at room temperature, and incubated with Anti-DIG-POD (Roche, Mannheim, Germany) labeled with alkaline phosphatase at 1:2000 dilution for 30 min at room temperature. The sections were then washed with DIG1 buffer. After washing, the TSA Plus Fluorescein System (PerkinElmer, Boston, MA, USA) was used for the amplification of the hybridization signals. The nuclei were stained by DAPI (Invitrogen, Carlsbad, CA, USA) staining. Fluorescence signals were captured by confocal microscopy (Olympus FV3000) (Olympus, Tokyo, Japan). The probes and qPCR primer sequences [38,39] used in this study are listed in Table S2.

2.6. JQ1 Treatment

JQ1 was dissolved in DMSO. DMSO, PEG300, Tween-80, and saline were added in sequence at ratios of 5%, 40%, 5%, and 50%, respectively; the stock solution, which was diluted to 1.25 mg/mL, was prepared for freezing at −80 °C. Fish were treated intraperitoneally with 1.25 mg/mL JQ1 at 10, 20, and 30 dah, respectively. The control group was injected with the same dose of DMSO.
In order to investigate the role of BETs in gonadal development, the histology of tilapia testes treated with JQ1 and DMSO at 60, 90, and 120 dah was examined, respectively. Gonads were dissected from fish anesthetized with MS-222 (Sigma Aldrich, St. Louis, MO, USA) and fixed in Bouin solution for 24 h at room temperature. The fixed samples were processed as follows: gonads were dehydrated in 70%, 80%, 90%, and 95% ethanol for 15 min, and 100% ethanol, 3 times, for 5 min each time. They were then cleared in xylene and ethanol mixture (1:1) for 15 min, and continuously cleared twice in xylene for 30 min each time. Next, the gonads were placed in melted paraffin. Samples were sliced at a thickness of 5 µm using a Leica microtome (Leica Microsystems, Wetzlar, Germany). Finally, sections were stained with hematoxylin and eosin (H&E) [40]. The photos were taken under an Olympus BX51 optical microscope (Olympus, Tokyo, Japan). The qRT-PCR was performed to determine the expression levels of apoptosis-related genes, germ-cell marker genes and Sertoli-cell marker genes in JQ1-treated fish and control fish.

3. Results

3.1. Identification of BRDs from Tilapia and Representative Species

Given the fact that the whole genome duplication (WGD) can drive the expansion of gene families, we isolated BRD genes from the genome of 24 representative species by a homology-based search (Figure 1 and Table S1). The contig N50, scaffold N50, annotated gene numbers, and BUSCO information of 24 representative species are listed in Table S4. In the fruit fly, 16 BRDs were identified. In protozoa (vase tunicate and lancelet), 17–20 BRD genes were isolated. In jawless vertebrates (lamprey), 26 members were identified. In species that experienced the 2R WGD event, there was a significant expansion of the BRD family. In total, 38, 40, 39, 39, 38, and 37 BRDs were identified in the elephant shark, coelacanth, spotted gar, frog, python, and chicken, respectively. In total, 43 BRDs were isolated in humans, and all members in humans, except taf1l, were identified in the mouse. Further expansion of BRDs (54–62) was observed in teleosts experiencing the 3R event. Among those teleosts, 54 BRDs were isolated in tilapia. In common carp, a species that underwent the 4R event, 120 BRDs were isolated. Interestingly, taf1l was also lost during the evolution of teleosts. Furthermore, four speckled protein (SP) genes, sp140, sp140l, sp110, and sp100, were in mammals and channel catfish, but some were lost in other species. For example, sp100 was lost in medaka, while sp100 and sp140l were lost in zebrafish. All the four SP genes were absent in the Mexican tetra, python, and chicken. In addition, two to four copies of brd2, brd3, brd4, baz2b, ep300, crebbp, brpf3, kmt2b, brd1, and brd8 were observed in teleosts, including largemouth bass, tilapia, medaka, Amazon molly, Southern platyfish and common carp, while only one copy was present in tetrapods, suggesting these genes were duplicated in the 3R-WGD. However, some of those genes were lost during the evolution of teleosts. For example, only one copy of brd2 was found in torafugu, channel catfish, striped catfish, and yellow catfish.

3.2. Phylogenetic and Syntenic Analyses of BRD Members

Phylogenetic analysis was performed using the conserved domain of all BRDs from tilapia, humans, and zebrafish. The BRD family was divided into eight subfamilies (Figure 2). Subfamily I included four members (cecr2, bptf, kat2a, and kat2b). Subfamily II (the BET subfamily) included brd2, brd3, brd4, and brdt. Subfamily III included eight members (baz1a, baz1b, brwd3, brwd1, phip, crebbp, ep300, and brd8). Subfamily IV comprised the ATPase proteins (atad2 and atad2b), brd1, brpf1, brpf3, brd7, and brd9. Subfamily V included the speckled proteins (sp100, sp110, sp140, and sp140l), tripartite-motif-containing (TRIM) proteins (trim24, trim28, trim33 and trim66), and baz2. Subfamily VI comprised kmt2a (mll) and kmt2b. Subfamily VII included four members (zmynd8, zmynd8l, taf1, and zmynd11). Finally, subfamily VIII included six members (ash1l, pbrm1, pbrm1l, smarca2, smarca4, and smarca4l). To understand the evolutionary history of BRDs after 3R, we performed syntenic analysis of brd2, brd3, brd4, baz2b, ep300, crebbp, brpf3, kmt2b, brd,1 and brd8 from representative tetrapods and teleosts (Figure S1). In most of the vertebrates investigated, the following genes, tbr1, rbms1, acd302, pkp4, tanci, dap1, and wdsub1, were found in the neighboring regions of baz2b, of which tanci existed in all tetrapods and duplicated in 3R species. Some genes, adjacent to brd2, brd3, brd4, ep300, crebbp, brpf3, and kmt2b, such as tap1, wrdr5, akap8l, xpnpep3, slx4, mapk13, and igflr1, existed with only one copy in the 3R species. Conserved synteny of brd2, brd3, brd4, baz2b, ep300, crebbp, brpf3, kmt2b, brd1, and brd8, and their upstream and downstream genes, was observed in tilapia, cave fish, zebrafish, and medaka. However, the neighboring region of these genes did not share conserved synteny with the fruit fly, vase tunicate or lancelet (data not shown).

3.3. Chromosomal Localization of BRD Genes and BRD Protein Domains in Tilapia

The 54 BRD genes in tilapia were unevenly distributed in its genome, covering 19 of the 22 linkage groups (LGs) (Figure S2). In total, seven BRDs were on LG4, and five BRDs were on LG6, LG20, and LG22, respectively. Only one BRD gene was on LG2, LG8, LG16, LG18, LG19, and LG23. Interestingly, paralogs of brd4 and crebbp were found on LG4 and LG6, and paralogs of brd2 and kmt2b were found on LG11 and LG22. This uneven distribution pattern was also found in zebrafish and humans (Table S5). To further characterize BRDs, we predicted the domain of these proteins in tilapia. As expected, all BRD proteins had at least one bromodomain (Figure S3). The BET subfamily members displayed two bromodomains. Members of subfamily Ⅵ had only one bromodomain, while Pbrml and Pbrm had five and six bromodomains, respectively. In addition to Kat2a, Kat2b, Cecr2, Brd8a, Brd8b, Brd7, and Brd9, BRD proteins had other domains besides the bromodomain, such as the ET domain, PHD domain, BAH, and HMG domain.

3.4. Expression and Distribution of BRDs in Eight Tilapia Tissues

Based on previously published transcriptomic data [28], the expression patterns of BRD genes in eight tissues from adult tilapia were analyzed, including ovary, testis, skeletal muscle, liver, kidney, head kidney, heart, and brain. Most BRD genes were expressed in at least one tissue (Figure 3). In total, 53 out of 54 genes were expressed in testis tissue, while kmt2a showed background expression level. Among these genes, the expression of pbrm1 was the highest in testis tissue, while brd2a showed the highest expression in ovary tissue. Notably, brd2a was also the highest-expressed gene in muscle and heart tissues. Brd8a, brpf3b, and smarca4 were the highest-expressed genes in liver, kidney and brain tissues, respectively.

3.5. Expression of BRDs in Tilapia Gonads at Different Stages

The expression patterns of BRDs in XX and XY tilapia gonads at 5, 30, 90, and 180 dah (days after hatching) indicated all 54 BRD genes were expressed in at least one of the gonad samples (Figure 4). Interestingly, four genes (zmynd11, sp100.1, phip, and brd3a) were highly expressed and peaked at 5 dah (key stage of molecular sex determination), and then gradually decreased in later stages (30, 90, and 180 dah). The expression of other genes was at a relatively low level at early developmental stages (5 and 30 dah), with no significant differences between XX and XY gonads. Their expression gradually increased at 90 and 180 dah, and more than 55% (30/54) of those genes showed sexually dimorphic expression (Table S3) at 90 or 180 dah. Among them, 21 BRD genes (cecr2, brd4a, brdt, phip, ep300a, brd8a, brpf3a, brd1a, atad2l, trim66, baz2b, kmt2a, kmt2bb, zmynd11, zmynd8, pbrm1, smarca2, ash1l, brd9, atad2, and sp100.1) were highly expressed in testis tissue, while 9 BRD members (bptfl, kat2a, kat2b, brd2a, baz1b, ep300b, brpf3b, baz2l and zmynd8l) were dominantly expressed in ovary tissue. Interestingly, some paralogous genes displayed different expression patterns. For example, ep300a was highly expressed in testis tissue, while ep300b was highly expressed in ovary tissue.

3.6. Expression Validation and Cellular Location of BRDs

To validate the transcriptome data, brdt and brd2b were selected for qRT-PCR validation. The expression of brdt and brd2b was consistently higher in testis tissue than in ovary tissue with qRT-PCR (Figure 5A,B). FISH was also performed for brdt to detect its cellular location in gonad tissue. Strong signals of brdt were observed in spermatogonia and spermatocytes of the testis, and signals were also detected in ovarian I, II and III oocytes (Figure 5C,D).
To explore the role of the BET subfamily in tilapia spermatogenesis, we injected Nile tilapia (10, 20, and 30 dah) intraperitoneally with the BET inhibitor JQ1. Morphological observations indicated that JQ1 treatment led to shrinkage of the seminiferous layer, with a reduced number of spermatogonia, spermatocytes, and spermatids, compared with the control fish (Figure 6A–E). The expression of the apoptosis-promoting genes baxa, caspase3a, caspase3b, caspase8, and caspase9 was upregulated, while the expression of the apoptosis-inhibiting gene bcl2 was downregulated in the testis of treated fish (Figure 6F). The qRT-PCR test revealed that the expression of germ cell-related genes oct4 (spermatogonia), rec8a (spermatocytes), prm and spata18 (spermatids) was significantly decreased in the testes of treated fish compared with that of control fish (Figure 6G). Compared with control fish, the expression of Sertoli-cell marker genes (amh and gsdf) in the testes of treated fish was significantly decreased (Figure 6H). Moreover, the expression of steroidogenic enzyme genes cyp11a1, cyp11c1, cyp17a1, and star1 was also decreased in the testes of treated fish (Figure 6I). However, ovaries were not significantly different after JQ1 treatment compared to controls (Figure S4).

4. Discussion

BRDs, a family of evolutionarily conserved proteins, play a vital role in gene regulation by ordered recruitment, anchoring, and regulation of various chromatin-modifying factors. To date, identification and expression of BRDs were undertaken in only a few species. Previously, 43 different BRDs were identified and divided into eight subfamilies according to the sequence similarity in humans [3,5]. However, 22 BRDs were identified and classified into six groups on the basis of the conserved sequences in buffalo [9], which was incomplete and inconclusive. Therefore, the number and evolution of the BRD family require clarification. In this study, we identified BRDs in 24 representative animals using bioinformatic analysis. WGD usually generated duplicate homologous genes during evolution [41]. It is generally believed that three rounds of WGD occurred during vertebrate evolution. The first two rounds of duplication events (1R and 2R) occurred early in the vertebrate lineage, while the third duplication event (3R) occurred only in teleosts [42,43,44,45]. Some teleosts, such as common carp, have even undergone a fourth round of genome duplication (4R) [28]. This study indicates an expansion of the BRD family members via 2R, 3R, and 4R WGD. Although teleosts experienced the 3R-WGD [46], only 10 genes (brd2, brd3, brd4, baz2b, ep300, crebbp, brpf3, kmt2b, brd1, and brd8) had paralogs, and most duplicated genes were lost secondarily under selection pressure, as in the case of many other gene families [44,45]. The exact number of BRD genes might change slightly due to incomplete genome sequencing, assembly, and annotation. However, these will not change the expansion pattern of this family in jawed vertebrates due to whole genome duplication. Interestingly, the duplicated paralogs of brd2 and kmt2b were on LG11 and LG22, respectively, similar to cbx3, rps19, and rps27-1 [47,48], suggesting possible inter-chromosomal duplication in tilapia. Moreover, uneven chromosomal distribution of BRDs was observed in tilapia, zebrafish, and humans.
Bromodomains can recognize histone acetylation sites, and different numbers of bromodomains decide different binding affinity to histones among BRDs [49]. Thus, Pbrml and Pbrm might have a higher binding affinity because they contain five and six bromodomains, respectively. Moreover, the histone-binding sites of BRD proteins are selective. According to previous reports, Kat2a recognizes H3K27Ac [50], and Brd2 can bind H4K12ac [51] and H4K5acK12ac [52]. Brpf1 acetylates histones H2A, H2B, H3 and H4 [53], while Bptf binds specifically to H3K4me3 [54]. Pbrm1 interacts broadly with H3K36ac, H3K14ac, H3K115ac, H4K12ac, H2AK15ac, H2BK15ac, H3K122ac, H2BK24ac, and H2BK116ac [49]. Additional domains (e.g., ET, HMG, and PHD) of different BRDs contribute to specific histone-binding site recognition.
Gene expression analysis is essential for understanding protein functions in biological processes. Previous studies investigated the expression pattern of BRDs and indicated that these genes are widely expressed in different tissues at different stages. For example, brd2, brd3, and brd4 are expressed in testis, ovary, brain, muscle and kidney tissue [3,55,56,57]. BRDs were consistently found to be ubiquitously expressed in different tilapia tissue based on transcriptomic data. Notably, these genes were not uniformly expressed in different tissues. Brd8 was reported to play an important role in the development and progression of the human liver [58]. In the present study, brd8b was the highest expressed gene in the tilapia liver, indicating its possible role in liver development. Brpf3 was shown to be critical in renal clear cell carcinoma (KIRC) [59]. Interestingly, brpf3b was the most highly expressed gene in the tilapia kidney, suggesting it may affect kidney development. Smarca4, a highly expressed gene in tilapia brain tissue, was shown to play an important role in the development of child brain stems [60,61]. Zmynd8, a gene preferentially expressed in the tilapia brain, plays an important role in Xenopus neural differentiation [62]. Brd7, highly expressed in the testes of tilapia, demonstrated its potential role in spermatogenesis and male fertility in the mouse [63]. Among the genes expressed in the tilapia ovary, brd2a exhibited the highest expression. A previous study showed that brd2 was located in oocytes and involved in mitotic and possibly meiotic cell-cycle regulation in the mouse ovary [64]. Interestingly, different expression profiles of brd2a and brd2b paralogs were found in tilapia tissues, as is the case in zebrafish [10], supporting functional diversity of duplicated genes based on differences in structural domain configurations. Among the genes expressed in testis tissue, pbrm1 showed the highest expression level, which was reported to play an important role in gonad primordium formation during embryogenesis of C. elegans [65]. However, no study was undertaken to clarify its function in fish testis tissue.
Most BRDs displayed sexual dimorphic expression during gonadal development in tilapia. In total, 21 BRD genes (cecr2, brd4a, brdt, phip, ep300a, brd8a, brpf3a, brd1a, atad2l, trim66, baz2b, kmt2a, kmt2bb, zmynd11, zmynd8, pbrm1, smarca2, ash1l, brd9, atad2, and sp100.1) were expressed higher in testes than in ovaries, suggesting they were potentially involved in testis development and the maintenance of spermatogenesis. For example, cecr2 was highly expressed in adult spermatogonia, and disruption of cecr2 led to compromised ability to fertilize oocytes in the male mouse [66]. Ep300a was also reported to be expressed higher in the testis tissue of the medaka and Mexican tetra [67]. On the contrary, nine BRD members (bptfl, kat2a, kat2b, brd2a, baz1b, ep300b, brpf3b, baz2l and zmynd8l) were dominantly expressed in ovarian tissue, suggesting their potential roles in ovary development. However, there is no study into the function of these genes in the ovary, which is worthy of further study. Brdt was reported to be testis specific [19] and later study found that it was also expressed at a relatively low level in mouse oocytes [68]. Brdt was proven to be essential for the normal progression of spermatogenesis [69], and loss of brdt produced an arrest of spermatogenesis [24,70], with severely impaired chromatin organization and absence of post-meiotic cells in mice [71]. In gilthead seabream, brdt was specifically expressed in secondary spermatocyte and spermatids, whereas in the ovary, brdt was expressed in late previtellogenic oocytes and early vitellogenic oocytes. In European seabass, brdt was expressed in late previtellogenic oocytes and spermatids [23]. In the present study, strong signals of brdt were observed in spermatogonia and spermatocytes of the testis, and moderate signals were also detected in I, II and III oocytes in tilapia. The cell location of brdt in tilapia testes was different from the previous studies in mouse, gilthead seabream, and European seabass, supporting the theory that there might be some interspecific differences in the functions of brdt [72].
The BET subfamily was reported to be involved in spermatogenesis [21]. In mice, brd2 was expressed in spermatogonia, spermatocytes, and round spermatids [19]. Brd3 was expressed in round spermatids, while brd4 was expressed in spermatogonia [19]. Previous studies demonstrated that the BET bromodomain inhibitor JQ1 could penetrate the blood–testis barrier, and JQ1 treatment would result in reduced testis size, germ-cell defects, and reduced sperm motility and spermatid cells in mice [22,26]. In the present study, JQ1 treatment resulted in abnormal seminiferous layers and a reduced number of spermatogonia, spermatocytes, and spermatids, which closely resembles the findings in the characterization of brdt hypomorphic mice [24] and the potent inhibition of brdt in mice with JQ1 treatment [22]. However, we cannot rule out that JQ1 play its role through inhibition of other BET subfamily members which are also highly expressed in tilapia testes (e.g., brd2b, brd4b and brd3a). Previous studies showed that JQ1 could induce apoptosis [73], which was consistent with the upregulation of pro-apoptosis genes in JQ1-treated tilapia. Inhibition of brdt by JQ1 did not affect testosterone or gonadotropin levels in mice [22]. In the present study, the expression of cyp11a1, cyp17a1, and cyp11c1 was downregulated in fish treated with JQ1, suggesting that the JQ1 treatment possibly blocked androgen synthesis in tilapia. However, the actual mechanism of JQ1 acting on gene expression in tilapia spermatogenesis remains to be elucidated. Because the BET subfamily may function as transcription factors, transcriptional co-regulators, and chromatin remodelers, we speculate that JQ1 treatment changed gene expression by altering chromatin conformation and, ultimately, led to abnormal spermatogenesis in tilapia, as in mice [26].

5. Conclusions

In this study, we completely identified the BRD family in 24 representative species by bioinformatic analyses. Differences in structures of the BRD members may be related to the diversity of their functions. Most BRD genes were found to be sexually dimorphic at different stages, with 21 genes being dominantly expressed in the testis and 9 genes dominantly expressed in the ovary, suggesting their important roles in gonad development. JQ1 treatment resulted in a reduced number of germ cells and abnormal spermatogenesis in tilapia. This study provided a new perspective on the evolution of the BRDs and laid a solid foundation for elucidating the role of the BRDs in vertebrates, especially in spermatogenesis.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ani12172266/s1, Figure S1: Synteny of BRDs in representative species. Figure S2: Chromosomal maps of the BRDs in tilapia. Figure S3: Gene structure of BRDs in Nile tilapia. Figure S4: Morphological and histological analyses of the ovaries of treated fish. Table S1: The accession numbers of BRDs isolated. Table S2: Primers used in the present study. Table S3: The sex-dimorphic genes expressed in XX and XY gonads. Table S4: The contig N50, scaffold N50, annotated gene numbers, and BUSCO of the genomes of 24 representative species. Table S5: Chromosome distribution of BRDs in different species.

Author Contributions

D.W., W.T. and C.X. conceived and designed the experiments; C.X. and M.Y. analyzed the data; C.X., Q.Z., Z.M., K.D. and H.Y. performed the experiments; W.T., C.X., M.Y., Q.Z. and J.W. contributed reagents/materials/analysis tools. W.T. and C.X. wrote the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (31861123001, 31972778, 32172953, 31972776 and 32172969), the National Key Research and Development Program of China (2018YFD0900202) and Chongqing Science and Technology Bureau (cstc2021ycjh-bgzxm0024 and CQYC201903173).

Institutional Review Board Statement

Animal experiments were conducted in accordance with the regulations of the Guide for Care and Use of Laboratory Animals and were approved by the Committee of Laboratory Animal Experimentation at Southwest University, China. (No. IACUC-20181015-12).

Informed Consent Statement

Not applicable.

Data Availability Statement

In this paper, all methods using fish material were carried out in accordance with relevant guidelines and regulations. The data used and/or analyzed during the current study were obtained from the NCBI website (https://www.ncbi.nlm.nih.gov/ (accessed on 5 October 2021)); the Ensembl Genomes database (https://ensemblgenomes.org/ (accessed on 22 November 2021)); and Genomicus (https://www.genomicus.bio.ens.psl.eu/genomicus-100.01/cgi-bin/search.pl (accessed on 22 November 2021)).

Conflicts of Interest

The authors declare that they have no competing interests.

References

  1. Haynes, S.R.; Dollard, C.; Winston, F.; Beck, S.; Trowsdale, J.; Dawid, I.B. The bromodomain: A conserved sequence found in human, Drosophila and yeast proteins. Nucleic Acids Res. 1992, 20, 2603. [Google Scholar] [CrossRef] [PubMed]
  2. Jeanmougin, F.; Wurtz, J.M.; Le Douarin, B.; Chambon, P.; Losson, R. The bromodomain revisited. Trends Biochem. Sci. 1997, 22, 151–153. [Google Scholar] [CrossRef]
  3. Filippakopoulos, P.; Picaud, S.; Mangos, M.; Keates, T.; Lambert, J.P.; Barsyte-Lovejoy, D.; Felletar, I.; Volkmer, R.; Müller, S.; Pawson, T.; et al. Histone recognition and large-scale structural analysis of the human bromodomain family. Cell 2012, 149, 214–231. [Google Scholar] [CrossRef]
  4. Dhalluin, C.; Carlson, J.E.; Zeng, L.; He, C.; Aggarwal, A.K.; Zhou, M.M. Structure and ligand of a histone acetyltransferase bromodomain. Nature 1999, 399, 491–496. [Google Scholar] [CrossRef] [PubMed]
  5. Fujisawa, T.; Filippakopoulos, P. Functions of bromodomain-containing proteins and their roles in homeostasis and cancer. Nat. Rev. Mol. Cell Biol. 2017, 18, 246–262. [Google Scholar] [CrossRef]
  6. Muller, S.; Filippakopoulos, P.; Knapp, S. Bromodomains as therapeutic targets. Expert Rev. Mol. Med. 2011, 13, e29. [Google Scholar] [CrossRef]
  7. Lee, S.Y.; Kim, J.J.; Miller, K.M. Bromodomain proteins: Protectors against endogenous DNA damage and facilitators of genome integrity. Exp. Mol. Med. 2021, 53, 1268–1277. [Google Scholar] [CrossRef]
  8. Tamkun, J.W.; Deuring, R.; Scott, M.P.; Kissinger, M.; Pattatucci, A.M.; Kaufman, T.C.; Kennison, J.A. brahma: A regulator of Drosophila homeotic genes structurally related to the yeast transcriptional activator SNF2/SWI2. Cell 1992, 68, 561–572. [Google Scholar] [CrossRef]
  9. Zhang, J.; Huang, L.; Zhang, P.; Huang, X.; Yang, W.; Liu, R.; Sun, Q.; Lu, Y.; Zhang, M.; Fu, Q. Genomic Identification, Evolution, and Expression Analysis of Bromodomain Genes Family in Buffalo. Genes 2022, 13, 103. [Google Scholar] [CrossRef]
  10. DiBenedetto, A.J.; Guinto, J.B.; Ebert, T.D.; Bee, K.J.; Schmidt, M.M.; Jackman, T.R. Zebrafish brd2a and brd2b are paralogous members of the bromodomain-ET (BET) family of transcriptional coregulators that show structural and expression divergence. BMC Dev. Biol. 2008, 8, 39. [Google Scholar] [CrossRef] [Green Version]
  11. Pivot-Pajot, C.; Caron, C.; Govin, J.; Vion, A.; Rousseaux, S.; Khochbin, S. Acetylation-dependent chromatin reorganization by BRDT, a testis-specific bromodomain-containing protein. Mol. Cell Biol. 2003, 23, 5354–5365. [Google Scholar] [CrossRef] [PubMed]
  12. Severinsen, J.; Bjarkam, C.R.; Kiar-Larsen, S.; Olsen, I.M.; Nielsen, M.M.; Blechingberg, J.; Nielsen, A.L.; Holm, I.E.; Foldagr, L.; Young, B.D. Evidence implicating BRD1 with brain development and susceptibility to both schizophrenia and bipolar affective disorder. Mol. Psychiatry 2006, 11, 1126–1138. [Google Scholar] [CrossRef] [PubMed]
  13. Park, Y.A.; Lee, J.W.; Kim, H.S.; Lee, Y.Y.; Kim, T.J.; Choi, C.H.; Choi, J.J.; Jeon, H.K.; Cho, Y.J.; Ryu, J.Y. Tumor suppressive effects of bromodomain-containing protein 7 (BRD7) in epithelial ovarian carcinoma. Clin. Cancer Res. 2014, 20, 565–575. [Google Scholar] [CrossRef] [PubMed]
  14. Browne, J.A.; NandyMazumdar, M.; Paranjapye, A.; Leir, S.H.; Harris, A. The bromodomain containing 8 (BRD8) transcriptional network in human lung epithelial cells. Mol. Cell Endocrinol. 2021, 524, 111169. [Google Scholar] [CrossRef]
  15. Zaware, N.; Zhou, M.M. Bromodomain biology and drug discovery. Nat. Struct. Mol. Biol. 2019, 26, 870–879. [Google Scholar] [CrossRef] [PubMed]
  16. Russell, J.O.; Ko, S.; Saggi, H.S.; Singh, S.; Poddar, M.; Shin, D.; Monga, S.P. Bromodomain and Extraterminal (BET) Proteins Regulate Hepatocyte Proliferation in Hepatocyte-Driven Liver Regeneration. Am. J. Pathol. 2018, 188, 1389–1405. [Google Scholar] [CrossRef] [PubMed]
  17. Andricovich, J.; Perkail, S.; Kai, Y.; Casasanta, N.; Peng, W.; Tzatsos, A. Loss of KDM6A Activates Super-Enhancers to Induce Gender-Specific Squamous-like Pancreatic Cancer and Confers Sensitivity to BET Inhibitors. Cancer Cell 2018, 33, 512–526. [Google Scholar] [CrossRef]
  18. Piunti, A.; Hashizume, R.; Morgan, M.A.; Bartom, E.T.; Horbinski, C.M.; Marshall, S.A.; Rendleman, E.J.; Ma, Q.; Takahashi, Y.H.; Woodfin, A.R.; et al. Therapeutic targeting of polycomb and BET bromodomain proteins in diffuse intrinsic pontine gliomas. Nat. Med. 2017, 23, 493–500. [Google Scholar] [CrossRef]
  19. Shang, E.; Salazar, G.; Crowley, T.E.; Wang, X.; Lopez, R.A.; Wang, X.; Wolgemuth, D.J. Identification of unique, differentiation stage-specific patterns of expression of the bromodomain-containing genes Brd2, Brd3, Brd4, and Brdt in the mouse testis. Gene Expr. Patterns 2004, 4, 513–519. [Google Scholar] [CrossRef]
  20. Bryant, J.M.; Berger, S.L. Low-hanging fruit: Targeting Brdt in the testes. EMBO J. 2012, 31, 3788–3789. [Google Scholar] [CrossRef] [Green Version]
  21. Berkovits, B.D.; Wolgemuth, D.J. The Role of the Double Bromodomain-Containing BET Genes During Mammalian Spermatogenesis. Curr. Top. Dev. Biol. 2013, 102, 293–326. [Google Scholar]
  22. Matzuk, M.M.; McKeown, M.R.; Filippakopoulos, P.; Li, Q.; Ma, L.; Agno, J.E.; Lemieux, M.E.; Picaud, S.; Richard, N.Y.; Qi, J. Small-molecule inhibition of BRDT for male contraception. Cell 2012, 150, 673–684. [Google Scholar] [CrossRef] [PubMed]
  23. Úbeda-Manzanaro, M.; Ortiz-Delgado, J.B.; Sarasquete, C. The Bromodomain testis-specific gene (Brdt) characterization and expression in gilthead seabream, Sparus aurata, and European seabass, Dicentrarchus labrax. Eur. J. Histochem. 2016, 60, 2638. [Google Scholar] [CrossRef] [PubMed]
  24. Shang, E.; Nickerson, H.D.; Wen, D.; Wang, X.; Wolgemuth, D.J. The first bromodomain of Brdt, a testis-specific member of the BET sub-family of double-bromodomain-containing proteins, is essential for male germ cell differentiation. Development 2007, 134, 3507–3515. [Google Scholar] [CrossRef]
  25. Aston, K.I.; Krausz, C.; Laface, I.; Ruiz-Castane, E.; Carrell, D.T. Evaluation of 172 candidate polymorphisms for association with oligozoospermia or azoospermia in a large cohort of men of European descent. Hum. Reprod. 2010, 25, 1383–1397. [Google Scholar] [CrossRef] [PubMed]
  26. Wang, X.; Sang, M.; Gong, S.; Chen, Z.; Zhao, X.; Wang, G.; Li, Z.; Huang, Y.; Chen, S.; Xie, G. BET bromodomain inhibitor JQ1 regulates spermatid development by changing chromatin conformation in mouse spermatogenesis. Genes Dis. 2021, 9, 1062–1073. [Google Scholar] [CrossRef]
  27. Brawand, D.; Wagner, C.E.; Li, Y.I.; Malinsky, M.; Keller, I.; Fan, S.; Simakov, O.; Ng, A.Y.; Lim, Z.W.; Bezault, E. The genomic substrate for adaptive radiation in African cichlid fish. Nature 2014, 513, 375–381. [Google Scholar] [CrossRef]
  28. Tao, W.; Yuan, J.; Zhou, L.; Sun, L.; Sun, Y.; Yang, S.; Li, M.; Zeng, S.; Huang, B.; Wang, D. Characterization of gonadal transcriptomes from Nile tilapia (Oreochromis niloticus) reveals differentially expressed genes. PLoS ONE 2013, 8, 63604. [Google Scholar] [CrossRef]
  29. Emms, D.M.; Kelly, S. OrthoFinder: Solving fundamental biases in whole genome comparisons dramatically improves orthogroup inference accuracy. Genome Biol. 2015, 16, 157. [Google Scholar] [CrossRef]
  30. Tavares, D.L.M.; Canete, S.C.; Henkin, R.; Estrela, P.C.; Freitas, T.R.O.; Galante1, R.; Freitas, C.M.D.S. Taxonomy Browser: A biodiversity data management system. J. Comput. Interdiscip. Sci. 2011, 2, 37–46. [Google Scholar]
  31. Hall, T. BioEdit: A user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp. Ser. 1999, 41, 95–98. [Google Scholar]
  32. Kumar, S.; Stecher, G.; Tamura, K. MEGA7. Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [PubMed]
  33. Muffato, M.; Louis, A.; Poisnel, C.E.; Roest Crollius, H. Genomicus: A database and a browser to study gene synteny in modern and ancestral genomes. Bioinformatics 2010, 26, 1119–1121. [Google Scholar] [CrossRef]
  34. Chen, C.; Chen, H.; Zhang, Y.; Thomas, H.R.; Frank, M.H.; He, Y.; Xia, R. TBtools: An Integrative Toolkit Developed for Interactive Analyses of Big Biological Data. Mol. Plant 2020, 13, 1194–1202. [Google Scholar] [CrossRef]
  35. Tao, W.; Chen, J.; Tan, D.; Yang, J.; Sun, L.; Wei, J.; Conte, M.A.; Kocher, T.D.; Wang, D. Transcriptome display during tilapia sex determination and differentiation as revealed by RNA-Seq analysis. BMC Genom. 2018, 19, 363. [Google Scholar] [CrossRef] [PubMed]
  36. Sun, L.N.; Jiang, X.L.; Xie, Q.P.; Yuan, J.; Huang, B.F.; Tao, W.J.; Zhou, L.Y.; Nagahama, Y.; Wang, D.S. Transdifferentiation of differentiated ovary into functional testis by long-term treatment of aromatase inhibitor in Nile tilapia. Endocrinology 2014, 155, 1476–1488. [Google Scholar] [CrossRef]
  37. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
  38. Wang, F.; Qin, Z.; Li, Z.; Yang, S.; Gao, T.; Sun, L.; Wang, D. Dnmt3aa but Not Dnmt3ab is Required for Maintenance of Gametogenesis in Nile Tilapia (Oreochromis niloticus). Int J. Mol. Sci. 2021, 22, 10170. [Google Scholar] [CrossRef]
  39. Yang, L.; Li, Y.; Wu, Y.; Sun, S.; Song, Q.; Wei, J.; Sun, L.; Li, M.; Wang, D.; Zhou, L. Rln3a is a prerequisite for spermatogenesis and fertility in male fish. J. Steroid Biochem. Mol. Biol. 2020, 197, 105517. [Google Scholar] [CrossRef]
  40. Wang, D.S.; Kobayashi, T.; Zhou, L.Y.; Paul-Prasanth, B.; Ijiri, S.; Sakai, F.; Okubo, K.; Morohashi, K.; Nagahama, Y. Foxl2 up-regulates aromatase gene transcription in a female-specific manner by binding to the promoter as well as interacting with ad4 binding protein/steroidogenic factor 1. Mol. Endocrinol. 2007, 21, 712–725. [Google Scholar] [CrossRef]
  41. Magadum, S.; Banerjee, U.; Murugan, P.; Gangapur, D.; Ravikesavan, R. Gene duplication as a major force in evolution. J. Genet. 2013, 92, 155–161. [Google Scholar] [CrossRef] [PubMed]
  42. Hoegg, S.; Brinkmann, H.; Taylor, J.S.; Meyer, A. Phylogenetic timing of the fish-specific genome duplication correlates with the diversification of teleost fish. J. Mol. Evol. 2004, 59, 190–203. [Google Scholar] [CrossRef] [PubMed]
  43. Escriva, H.; Manzon, L.; Youson, J.; Laudet, V. Analysis of lamprey and hagfish genes reveals a complex history of gene duplications during early vertebrate evolution. Mol. Biol. Evol. 2002, 19, 1440–1450. [Google Scholar] [CrossRef] [Green Version]
  44. Chiu, C.H.; Dewar, K.; Wagner, G.P.; Takahashi, K.; Ruddle, F.; Ledje, C.; Bartsch, P.; Scemama, J.L.; Stellwag, E.; Fried, C. Bichir HoxA cluster sequence reveals surprising trends in ray-finned fish genomic evolution. Genome Res. 2004, 14, 11–17. [Google Scholar] [CrossRef]
  45. Kuraku, S.; Meyer, A.; Kuratani, S. Timing of genome duplications relative to the origin of the vertebrates: Did cyclostomes diverge before or after? Mol. Biol. Evol. 2009, 26, 47–59. [Google Scholar] [CrossRef] [PubMed]
  46. Gillis, W.Q.; St John, J.; Bowerman, B.; Schneider, S.Q. Whole genome duplications and expansion of the vertebrate GATA transcription factor gene family. BMC Evol. Biol. 2009, 9, 207. [Google Scholar] [CrossRef]
  47. Kuang, G.; Tao, W.; Zheng, S.; Wang, X.; Wang, D. Genome-Wide Identification, Evolution and Expression of the Complete Set of Cytoplasmic Ribosomal Protein Genes in Nile Tilapia. Int. J. Mol. Sci. 2020, 21, 1230. [Google Scholar] [CrossRef]
  48. Liu, X.Y.; Zhang, X.-B.; Li, M.-H.; Zheng, S.Q.; Liu, Z.L.; Cheng, Y.Y.; Wang, D.S. Genome-wide identification, evolution of chromobox family genes and their expression in Nile tilapia. Comp. Biochem. Physiol. B Biochem. Mol. Biol. 2017, 203, 25–34. [Google Scholar] [CrossRef]
  49. Charlop-Powers, Z.; Zeng, L.; Zhang, Q.; Zhou, M.M. Structural insights into selective histone H3 recognition by the human Polybromo bromodomain 2. Cell Res. 2010, 20, 529–538. [Google Scholar] [CrossRef]
  50. Park, S.; Stanfield, R.L.; Martinez-Yamout, M.A.; Dyson, H.J.; Wilson, I.A.; Wright, P.E. Role of the CBP catalytic core in intramolecular SUMOylation and control of histone H3 acetylation. Proc. Natl. Acad. Sci. USA 2017, 114, E5335–E5342. [Google Scholar] [CrossRef]
  51. Umehara, T.; Nakamura, Y.; Jang, M.K.; Nakano, K.; Tanaka, A.; Ozato, K.; Padmanabhan, B.; Yokoyama, S. Structural basis for acetylated histone H4 recognition by the human BRD2 bromodomain. J. Biol. Chem. 2010, 285, 7610–7618. [Google Scholar] [CrossRef] [PubMed]
  52. Huang, H.; Zhang, J.; Shen, W.; Wang, X.; Wu, J.; Wu, J.; Shi, Y. Solution structure of the second bromodomain of Brd2 and its specific interaction with acetylated histone tails. BMC Struct. Biol. 2007, 7, 57. [Google Scholar] [CrossRef] [PubMed]
  53. Lloyd, J.T.; Glass, K.C. Biological function and histone recognition of family IV bromodomain-containing proteins. J. Cell Physiol. 2018, 233, 1877–1886. [Google Scholar] [CrossRef]
  54. Li, H.; Ilin, S.; Wang, W.; Duncan, E.M.; Wysocka, J.; Allis, C.D.; Patel, D.J. Molecular basis for site-specific read-out of histone H3K4me3 by the BPTF PHD finger of NURF. Nature 2006, 442, 91–95. [Google Scholar] [CrossRef] [PubMed]
  55. Belkina, A.C.; Denis, G.V. BET domain co-regulators in obesity, inflammation and cancer. Nat. Rev. Cancer 2012, 12, 465–477. [Google Scholar] [CrossRef]
  56. Wang, F.; Liu, H.; Blanton, W.P.; Belkina, A.; Lebrasseur, N.K.; Denis, G.V. Brd2 disruption in mice causes severe obesity without Type 2 diabetes. Biochem. J. 2009, 425, 71–83. [Google Scholar] [CrossRef]
  57. Shang, E.; Wang, X.; Wen, D.; Greenberg, D.A.; Wolgemuth, D.J. Double bromodomain-containing gene Brd2 is essential for embryonic development in mouse. Dev. Dyn. 2009, 238, 908–917. [Google Scholar] [CrossRef]
  58. Yu, Z.; Chen, T.; Mo, H.; Guo, C.; Liu, Q. BRD8, which is negatively regulated by miR-876–3p, promotes the proliferation and apoptosis resistance of hepatocellular carcinoma cells via KAT5. Arch. Biochem. Biophys. 2020, 693, 108550. [Google Scholar] [CrossRef]
  59. Lu, J.; Qian, C.; Ji, Y.; Bao, Q.; Lu, B. Gene Signature Associated with Bromodomain Genes Predicts the Prognosis of Kidney Renal Clear Cell Carcinoma. Front. Genet. 2021, 12, 643935. [Google Scholar] [CrossRef]
  60. Jones, D.T.; Jager, N.; Kool, M.; Zichner, T.; Hutter, B.; Sultan, M.; Cho, Y.J.; Pugh, T.J.; Hovestadt, V.; Stutz, A.M.; et al. Dissecting the genomic complexity underlying medulloblastoma. Nature 2012, 488, 100–105. [Google Scholar] [CrossRef]
  61. Hasselblatt, M.; Thomas, C.; Federico, A.; Nemes, K.; Johann, P.D.; Bison, B.; Bens, S.; Dahlum, S.; Kordes, U.; Redlich, A.; et al. SMARCB1-deficient and SMARCA4-deficient Malignant Brain Tumors with Complex Copy Number Alterations and TP53 Mutations May Represent the First Clinical Manifestation of Li-Fraumeni Syndrome. Am. J. Surg. Pathol. 2022, 46, 1277–1283. [Google Scholar] [CrossRef] [PubMed]
  62. Zeng, W.; Kong, Q.; Li, C.; Mao, B. Xenopus RCOR2 (REST corepressor 2) interacts with ZMYND8, which is involved in neural differentiation. Biochem. Biophys. Res. Commun. 2010, 394, 1024–1029. [Google Scholar] [CrossRef] [PubMed]
  63. Wang, H.; Zhao, R.; Guo, C.; Jiang, S.; Yang, J.; Xu, Y.; Liu, Y.; Fan, L.; Xiong, W.; Ma, J.; et al. Knockout of BRD7 results in impaired spermatogenesis and male infertility. Sci. Rep. 2016, 6, 21776. [Google Scholar] [CrossRef] [Green Version]
  64. Trousdale, R.K.; Wolgemuth, D.J. Bromodomain containing 2 (Brd2) is expressed in distinct patterns during ovarian folliculogenesis independent of FSH or GDF9 action. Mol. Reprod. Dev. 2004, 68, 261–268. [Google Scholar] [CrossRef]
  65. Shibata, Y.; Uchida, M.; Takeshita, H.; Nishiwaki, K.; Sawa, H. Multiple functions of PBRM-1/Polybromo- and LET-526/Osa-containing chromatin remodeling complexes in C. elegans development. Dev. Biol. 2012, 361, 349–357. [Google Scholar] [CrossRef] [PubMed]
  66. Thompson, P.J.; Norton, K.A.; Niri, F.H.; Dawe, C.E.; McDermid, H.E. CECR2 is involved in spermatogenesis and forms a complex with SNF2H in the testis. J. Mol. Biol. 2012, 415, 793–806. [Google Scholar] [CrossRef] [PubMed]
  67. Wang, X.; Yan, J. Directional divergence of Ep300 duplicates in teleosts and its implications. BMC Evol. Biol. 2020, 20, 140. [Google Scholar] [CrossRef]
  68. Paillisson, A.; Levasseur, A.; Gouret, P.; Callebaut, I.; Bontoux, M.; Pontarotti, P.; Monget, P. Bromodomain testis-specific protein is expressed in mouse oocyte and evolves faster than its ubiquitously expressed paralogs BRD2, -3, and -4. Genomics 2007, 89, 215–223. [Google Scholar] [CrossRef]
  69. Wang, P.J.; Manterola, M.; Brown, T.M.; Oh, M.Y.; Garyn, C.; Gonzalez, B.J.; Wolgemuth, D.J. BRDT is an essential epigenetic regulator for proper chromatin organization, silencing of sex chromosomes and crossover formation in male meiosis. PLOS Genet. 2018, 14, e1007209. [Google Scholar]
  70. Gaucher, J.; Boussouar, F.; Montellier, E.; Curtet, S.; Buchou, T.; Bertrand, S.; Hery, P.; Jounier, S.; Depaux, A.; Vitte, A.L. Bromodomain-dependent stage-specific male genome programming by Brdt. EMBO J. 2012, 31, 3809–3820. [Google Scholar] [CrossRef]
  71. Berkovits, B.D.; Wolgemuth, D.J. The first bromodomain of the testis-specific double bromodomain protein Brdt is required for chromocenter organization that is modulated by genetic background. Dev. Biol. 2011, 360, 358–368. [Google Scholar] [CrossRef] [PubMed]
  72. Barda, S.; Paz, G.; Yogev, L.; Yavetz, H.; Lehavi, O.; Hauser, R.; Botchan, A.; Breitbart, H.; Kleiman, S.E. Expression of BET genes in testis of men with different spermatogenic impairments. Fertil. Steril. 2012, 97, 46–52. [Google Scholar] [CrossRef] [PubMed]
  73. Jostes, S.; Nettersheim, D.; Fellermeyer, M.; Schneider, S.; Hafezi, F.; Honecker, F.; Schumacher, V.; Geyer, M.; Kristiansen, G.; Schorle, H. The bromodomain inhibitor JQ1 triggers growth arrest and apoptosis in testicular germ cell tumours in vitro and in vivo. J. Cell Mol. Med. 2017, 21, 1300–1314. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Phylogenetic relationship and numbers of BRD family members identified in 24 representative animals. The four rounds of whole genome duplication (WGD) are represented by 1R, 2R, 3R, and 4R, respectively. The fruit fly, Drosophila melanogaster; vase tunicate, Ciona intestinalis; lancelet, Branchiostoma belcheri; lamprey, Petromyzon marinus; elephant shark, Callorhinchus milii; coelacanth, Latimeria chalumnae; python, Python bivittatus; frog, Xenopus tropicalis; chicken, Gallus gallus; mouse, Mus musculus; human, Homo sapiens; spotted gar, Lepisosteus oculatus; torafugu, Takifugu rubripes; largemouth bass, Micropterus salmoides; Nile tilapia, Oreochromis niloticus; medaka, Oryzias latipes; Amazon molly, Poecilia formosa; Southern platyfish, Xiphophorus maculatus; Mexican tetra, Astyanax mexicanus; channel catfish, Ictalurus punctatus; striped catfish, Pangasianodon hypophthalmus; yellow catfish, Tachysurus fulvidraco; zebrafish, Danio rerio; and common carp, Cyprinus carpio.
Figure 1. Phylogenetic relationship and numbers of BRD family members identified in 24 representative animals. The four rounds of whole genome duplication (WGD) are represented by 1R, 2R, 3R, and 4R, respectively. The fruit fly, Drosophila melanogaster; vase tunicate, Ciona intestinalis; lancelet, Branchiostoma belcheri; lamprey, Petromyzon marinus; elephant shark, Callorhinchus milii; coelacanth, Latimeria chalumnae; python, Python bivittatus; frog, Xenopus tropicalis; chicken, Gallus gallus; mouse, Mus musculus; human, Homo sapiens; spotted gar, Lepisosteus oculatus; torafugu, Takifugu rubripes; largemouth bass, Micropterus salmoides; Nile tilapia, Oreochromis niloticus; medaka, Oryzias latipes; Amazon molly, Poecilia formosa; Southern platyfish, Xiphophorus maculatus; Mexican tetra, Astyanax mexicanus; channel catfish, Ictalurus punctatus; striped catfish, Pangasianodon hypophthalmus; yellow catfish, Tachysurus fulvidraco; zebrafish, Danio rerio; and common carp, Cyprinus carpio.
Animals 12 02266 g001
Figure 2. Phylogenetic tree of BRD family members of tilapia, humans and zebrafish. BRD family was divided into eight subfamilies (I–VIII). Numbers represent the bootstrap value. BRD protein sequences were aligned using Bioedit. The tree was constructed using NJ method in Mega 7.0. GenBank accession numbers are listed in Table S1.
Figure 2. Phylogenetic tree of BRD family members of tilapia, humans and zebrafish. BRD family was divided into eight subfamilies (I–VIII). Numbers represent the bootstrap value. BRD protein sequences were aligned using Bioedit. The tree was constructed using NJ method in Mega 7.0. GenBank accession numbers are listed in Table S1.
Animals 12 02266 g002
Figure 3. Transcriptome data of BRDs in eight tilapia tissues. Red shows high expression level and blue shows low expression level. Each column represents a different gene, and each row represents a different tissue.
Figure 3. Transcriptome data of BRDs in eight tilapia tissues. Red shows high expression level and blue shows low expression level. Each column represents a different gene, and each row represents a different tissue.
Animals 12 02266 g003
Figure 4. Transcriptomic data of BRDs at four developmental stages (5, 30, 90, and 180 dah) of Nile tilapia. A heatmap showing expression of BRDs in tilapia gonads at different developmental stages. Red shows high expression level and blue shows low expression level. dah, days after hatching.
Figure 4. Transcriptomic data of BRDs at four developmental stages (5, 30, 90, and 180 dah) of Nile tilapia. A heatmap showing expression of BRDs in tilapia gonads at different developmental stages. Red shows high expression level and blue shows low expression level. dah, days after hatching.
Animals 12 02266 g004
Figure 5. Expression validation and cellular location of BRDs. Expression validation of brdt and brd2b by quantitative real time PCR (A,B), and cellular location of brdt by fluorescence in situ hybridization in testis (C) and ovary (D) at 180 dah. Nuclei were counterstained with DAPI (a,d). Red fluorescence represents the signal of brdt (b,e). Green fluorescence represents the Vasa signals (c,f). gapdh was used as an internal control. Results are presented as mean ± SD. Asterisk above the error bar indicates statistical differences at p < 0.05 as determined by one-way ANOVA followed by Tukey test (*** p < 0.001). B, brain; P, pituitary; G, gill; H, heart; S, spleen; L, liver; I, intestinal; O, ovary; T, testis; K, kidney; HK, head kidney; M, muscle; I to III, phase I to III oocytes; SG, spermatogonia; PSC, primary spermatocytes; SSC, secondary spermatocytes. The white arrow indicates co-localization of Vasa and brdt signals.
Figure 5. Expression validation and cellular location of BRDs. Expression validation of brdt and brd2b by quantitative real time PCR (A,B), and cellular location of brdt by fluorescence in situ hybridization in testis (C) and ovary (D) at 180 dah. Nuclei were counterstained with DAPI (a,d). Red fluorescence represents the signal of brdt (b,e). Green fluorescence represents the Vasa signals (c,f). gapdh was used as an internal control. Results are presented as mean ± SD. Asterisk above the error bar indicates statistical differences at p < 0.05 as determined by one-way ANOVA followed by Tukey test (*** p < 0.001). B, brain; P, pituitary; G, gill; H, heart; S, spleen; L, liver; I, intestinal; O, ovary; T, testis; K, kidney; HK, head kidney; M, muscle; I to III, phase I to III oocytes; SG, spermatogonia; PSC, primary spermatocytes; SSC, secondary spermatocytes. The white arrow indicates co-localization of Vasa and brdt signals.
Animals 12 02266 g005
Figure 6. Inhibition of BET subfamily by JQ1 led to abnormal spermatogenesis. HE staining of Nile tilapia testis sections at 90 dah (AD). The number of spermatogonia, spermatocytes, and spermatid in the testes of treated fish was significantly reduced, compared with that in the testes of control fish (E). The expression of apoptosis-related genes (F). Significantly decreased expression of germ-cell marker genes in the testes of treated fish, compared with that of control male fish (G). Significantly decreased expression of Sertoli-cell marker genes in the testes of treated fish, compared with that of control male fish (H). Significantly decreased expression of Leydig-cell marker genes in the testes of treated fish, compared with that of control male fish (I). gapdh was used as an internal control. Results are presented as mean ± SD. Asterisk above the error bar indicates statistical differences at p < 0.05 as determined by one-way ANOVA followed by Tukey test (n = 3–10; * p < 0.05, ** p < 0.01, *** p < 0.001). SG, spermatogonia; SC, spermatocytes; ST, spermatid. Scale bars 50 μm (A,C), 10 μm (B,D).
Figure 6. Inhibition of BET subfamily by JQ1 led to abnormal spermatogenesis. HE staining of Nile tilapia testis sections at 90 dah (AD). The number of spermatogonia, spermatocytes, and spermatid in the testes of treated fish was significantly reduced, compared with that in the testes of control fish (E). The expression of apoptosis-related genes (F). Significantly decreased expression of germ-cell marker genes in the testes of treated fish, compared with that of control male fish (G). Significantly decreased expression of Sertoli-cell marker genes in the testes of treated fish, compared with that of control male fish (H). Significantly decreased expression of Leydig-cell marker genes in the testes of treated fish, compared with that of control male fish (I). gapdh was used as an internal control. Results are presented as mean ± SD. Asterisk above the error bar indicates statistical differences at p < 0.05 as determined by one-way ANOVA followed by Tukey test (n = 3–10; * p < 0.05, ** p < 0.01, *** p < 0.001). SG, spermatogonia; SC, spermatocytes; ST, spermatid. Scale bars 50 μm (A,C), 10 μm (B,D).
Animals 12 02266 g006
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Xu, C.; Yu, M.; Zhang, Q.; Ma, Z.; Du, K.; You, H.; Wei, J.; Wang, D.; Tao, W. Genome-Wide Identification and Characterization of the BRD Family in Nile Tilapia (Oreochromis niloticus). Animals 2022, 12, 2266. https://doi.org/10.3390/ani12172266

AMA Style

Xu C, Yu M, Zhang Q, Ma Z, Du K, You H, Wei J, Wang D, Tao W. Genome-Wide Identification and Characterization of the BRD Family in Nile Tilapia (Oreochromis niloticus). Animals. 2022; 12(17):2266. https://doi.org/10.3390/ani12172266

Chicago/Turabian Style

Xu, Chunmei, Miao Yu, Qingqing Zhang, Zhisheng Ma, Kang Du, Huiqin You, Jing Wei, Deshou Wang, and Wenjing Tao. 2022. "Genome-Wide Identification and Characterization of the BRD Family in Nile Tilapia (Oreochromis niloticus)" Animals 12, no. 17: 2266. https://doi.org/10.3390/ani12172266

APA Style

Xu, C., Yu, M., Zhang, Q., Ma, Z., Du, K., You, H., Wei, J., Wang, D., & Tao, W. (2022). Genome-Wide Identification and Characterization of the BRD Family in Nile Tilapia (Oreochromis niloticus). Animals, 12(17), 2266. https://doi.org/10.3390/ani12172266

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop