Rounds Of Genome Duplication Research Articles

Comparative Analysis of Transposable Elements in Strawberry Genomes of Different Ploidy Levels.

Transposable elements (TEs) make up a large portion of plant genomes and play a vital role in genome structure, function, and evolution. Cultivated strawberry (Fragaria x ananassa) is one of the most important fruit crops, and its octoploid genome was formed through several rounds of genome duplications from diploid ancestors. Here, we built a pan-genome TE library for the Fragaria genus using ten published strawberry genomes at different ploidy levels, including seven diploids, one tetraploid, and two octoploids, and performed comparative analysis of TE content in these genomes. The TEs comprise 51.83% (F. viridis) to 60.07% (F. nilgerrensis) of the genomes. Long terminal repeat retrotransposons (LTR-RTs) are the predominant TE type in the Fragaria genomes (20.16% to 34.94%), particularly in F. iinumae (34.94%). Estimating TE content and LTR-RT insertion times revealed that species-specific TEs have shaped each strawberry genome. Additionally, the copy number of different LTR-RT families inserted in the last one million years reflects the genetic distance between Fragaria species. Comparing cultivated strawberry subgenomes to extant diploid ancestors showed that F. vesca and F. iinumae are likely the diploid ancestors of the cultivated strawberry, but not F. viridis. These findings provide new insights into the TE variations in the strawberry genomes and their roles in strawberry genome evolution.

Genomic analysis of hypoxia inducible factor alpha in ray-finned fishes reveals missing Ohnologs and evidence of widespread positive selection

As aquatic hypoxia worsens on a global scale, fishes will become increasingly challenged by low oxygen, and understanding the molecular basis of their response to hypoxia may help to better define the capacity of fishes to cope with this challenge. The hypoxia inducible factor (HIF) plays a critical role in the molecular response to hypoxia by activating the transcription of genes that serve to improve oxygen delivery to the tissues or enhance the capacity of tissues to function at low oxygen. The current study examines the molecular evolution of genes encoding the oxygen-dependent HIFα subunit (HIFA) in the ray-finned fishes (Actinopterygii). Genomic analyses demonstrate that several lineages retain four paralogs of HIFA predicted from two rounds of genome duplication at the base of vertebrate evolution, broaden the known distribution of teleost-specific HIFA paralogs, and provide evidence for salmonid-specific HIFA duplicates. Evolution of the HIFA gene family is characterized by widespread episodic positive selection at amino acid sites that potentially mediate protein stability, protein–protein interactions, and transcriptional regulation. HIFA transcript abundance depends upon paralog, tissue, and fish lineage. A phylogenetically-informed gene nomenclature is proposed along with avenues for future research on this critical family of transcription factors.

Knockout of microphthalmia-associated transcription factor (mitf) confers a red and yellow tilapia with few pigmented melanophores

As the most important genetic factor for melanophore differentiation and melanogenesis, mitf has been studied for over 30 years. The phenotypes of mitf mutants have also received continuous attention for nearly a century. Due to the third round of genome duplication, teleosts have two mitf genes. Mutation analysis demonstrates that mitfa plays an important role in zebrafish pigmentation. However, there have been no functional studies on body color regulation by mitfb or studies of mitfa;mitfb double mutants. In the present study, we mutated both mitf genes in tilapia using CRISPR/Cas9. Disruption of mitfa resulted in light yellow body color with weak gray vertical bars due to significantly reduced numbers of melanophores and increased sizes of xanthophores, while disruption of mitfb led to slight hypo-pigmentation. Double mutation of mitfa and mitfb resulted in dramatic hypo-pigmentation in trunk and fins due to loss of additional melanophores and increase in the number of iridophores compared to the single mutants, but had no influence on retinal pigment epithelium (RPE) pigmentation. The mitfa−/−;mitfb−/− mutants were yellow with black spots at early juvenile stage (60 dpf), yellow-reddish at late juvenile stage (120 dpf) and red and yellow (including iris) at adult stage (180 dpf), due to increased erythrophores and enlarged xanthophores. Our results demonstrated that both mitf genes are important for body color formation in tilapia, but mitfa plays a more important role than mitfb. Additionally, by comparing the phenotypes of the mitfa−/−;mitfb−/−, pmela−/−;pmelb−/− and hps4−/− mutants, we found that disruption of genes with different functions in different aspects of melanogenesis could lead to different body colors in tilapia. Mitf is the most important gene for melanophore differentiation and is probably also critical for differentiation of erythrophores, xanthophores and iridophores. To our knowledge, this is the first report showing that both mitf genes are involved in body color formation by loss of function study, and the first report showing that erythrophore numbers and xanthophore sizes were affected by mitf. Our research not only serves as a model for studying mitf function in tilapia and the closely related cichlids, but also provides new strategies for breeding red and yellow tilapia for aquaculture.

Dnmt3aa but Not Dnmt3ab Is Required for Maintenance of Gametogenesis in Nile Tilapia (Oreochromis niloticus).

Dnmt3a, a de novo methyltransferase, is essential for mammalian germ line DNA methylation. Only one Dnmt3a is identified in mammals, and homozygous mutants of Dnmt3a are lethal, while two Dnmt3a paralogs, dnmt3aa and dnmt3ab, are identified in teleosts due to the third round of genome duplication, and homozygous mutants of dnmt3aa and dnmt3ab are viable in zebrafish. The expression patterns and roles of dnmt3aa and dnmt3ab in gonadal development remain poorly understood in teleosts. In this study, we elucidated the precise expression patterns of dnmt3aa and dnmt3ab in tilapia gonads. Dnmt3aa was highly expressed in oogonia, phase I and II oocytes and granulosa cells in ovaries and spermatogonia and spermatocytes in testes, while dnmt3ab was mainly expressed in ovarian granulosa cells and testicular spermatocytes. The mutation of dnmt3aa and dnmt3ab was achieved by CRISPR/Cas9 in tilapia. Lower gonadosomatic index (GSI), increased apoptosis of oocytes and spermatocytes and significantly reduced sperm quality were observed in dnmt3aa−/− mutants, while normal gonadal development was observed in dnmt3ab−/− mutants. Consistently, the expression of apoptotic genes was significantly increased in dnmt3aa−/− mutants. In addition, the 5-methylcytosine (5-mC) level in dnmt3aa−/− gonads was decreased significantly, compared with that of dnmt3ab−/− and wild type (WT) gonads. Taken together, our results suggest that dnmt3aa, not dnmt3ab, plays important roles in maintaining gametogenesis in teleosts.

Genomes Vary in Size and Spatial Patterns Within Chimeric Blades of Porphyra spp.

Genome size variation is of crucial biological importance, however variation in genome sizes within a single individual/organism is rarely reported except for some species groups such as algae where polygenomy, endopolyploidy and mixopolyploidy have previously been reported. The red algal genus Porphyra forms part of very profitable marine food products commonly known as ‘Nori’. Farming of these valuable marine crops was revolutionized by the discovery of their life cycle in the 40’s. One of the most remarkable characteristics of these taxa is the formation of chimeric gametophytic thalli. After meiosis, the four meiotic products are not released as individuals spores, but instead develop together into a single leafy thallus through successive mitotic divisions. In this study, we used flow cytometry to estimate genome sizes in 670 vegetative thallus sections from 195 blades from three Porphyra species, to determine if this chimerism could be related to the presence of multiple genome sizes and mixoploidy within thalli. Our results show a wide variety of genome sizes both within and between thalli. We interpreted these results as the presence of two different genome types of different sizes (a and b) with separate rounds of genome duplications within the vegetative thalli. By analyzing several sections per thallus, we were able to show that the different genome types and ploidy levels are not distributed haphazardly through the thallus, but are distributed along the thallus in a sectorial way in mosaics. In some individuals, the 2C genome size can either be interpreted as diploids or alternatively as haploid cells that are arrested at the G2-stage of the mitotic cycle, acting as diploid with two copies of their genome during most of their life-time. We conclude that Porphyra species belong to an aneuploid/euploid system, where genome duplications, mixoploidy, chromosomal dynamics and the presence of different genome types in the chimeric thalli play a role in shaping the genetic diversity of these taxa. Our results may have important implications to understand red algae biology and evolution and raise further questions on concepts of what constitutes an individual.

Specific cyprinid HIF isoforms contribute to cellular mitochondrial regulation

Hypoxia-inducible factor 1 (HIF-1) functions as a master regulator of the cellular response to hypoxic stress. Two HIF-1α paralogs, HIF-1αA and HIF-1αB, were generated in euteleosts by the specific, third round of genome duplication, but one paralog was later lost in most families with the exception of cyprinid fish. How these duplicates function in mitochondrial regulation and whether their preservation contributes to the hypoxia tolerance demonstrated by cyprinid fish in freshwater environments is not clear. Here we demonstrated the divergent function of these two zebrafish Hif-1a paralogs through cellular approaches. The results showed that Hif-1aa played a role in tricarboxylic acid cycle by increasing the expression of Citrate synthase and the activity of mitochondrial complex II, and it also enhanced mitochondrial membrane potential and ROS production by reducing free Ca2+ in the cytosol. Hif-1ab promoted intracellular ATP content by up-regulating the activity of mitochondrial complexes I, III and IV and the expression of related genes. Furthermore, both the two zebrafish Hif-1a paralogs promoted mitochondrial mass and the expression level of mtDNA, contributing to mitochondrial biogenesis. Our study reveals the divergent functions of Hif-1aa and Hif-1ab in cellular mitochondrial regulation.

Discovery of four Noggin genes in lampreys suggests two rounds of ancient genome duplication

The secreted protein Noggin1 was the first discovered natural embryonic inducer produced by cells of the Spemann organizer. Thereafter, it was shown that vertebrates have a whole family of Noggin genes with different expression patterns and functional properties. For example, Noggin1 and Noggin2 inhibit the activity of BMP, Nodal/Activin and Wnt-beta-catenin signalling, while Noggin4 cannot suppress BMP but specifically modulates Wnt signalling. In this work, we described and investigated phylogeny and expression patterns of four Noggin genes in lampreys, which represent the most basally divergent group of extant vertebrates, the cyclostomes, belonging to the superclass Agnatha. Assuming that lampreys have Noggin homologues in all representatives of another superclass of vertebrates, the Gnathostomata, we propose a model for Noggin family evolution in vertebrates. This model is in agreement with the hypotheses suggesting two rounds of genome duplication in the ancestor of vertebrates before the divergence of Agnatha and Gnathostomata.

Evolution of the genes mediating phototransduction in rod and cone photoreceptors.

This paper reviews current knowledge of the evolution of the multiple genes encoding proteins that mediate the process of phototransduction in rod and cone photoreceptors of vertebrates. The approach primarily involves molecular phylogenetic analysis of phototransduction protein sequences, combined with analysis of the syntenic arrangement of the genes. At least 35 of these phototransduction genes appear to reside on no more than five paralogons - paralogous regions that each arose from a common ancestral region. Furthermore, it appears that such paralogs arose through quadruplication during the two rounds of genome duplication (2R WGD) that occurred in a chordate ancestor prior to the vertebrate radiation, probably around 600 millions years ago. For several components of the phototransduction cascade, it is shown that distinct isoforms already existed prior to WGD, with the likely implication that separate classes of scotopic and photopic photoreceptor cells had already evolved by that stage. The subsequent quadruplication of the entire genome then permitted the refinement of multiple distinct protein isoforms in rods and cones. A unified picture of the likely pattern and approximate timing of all the important gene duplications is synthesised, and the implications for our understanding of the evolution of rod and cone phototransduction are presented.

OHNOLOGS v2: a comprehensive resource for the genes retained from whole genome duplication in vertebrates.

All vertebrates including human have evolved from an ancestor that underwent two rounds of whole genome duplication (2R-WGD). In addition, teleost fish underwent an additional third round of genome duplication (3R-WGD). The genes retained from these genome duplications, so-called ohnologs, have been instrumental in the evolution of vertebrate complexity, development and susceptibility to genetic diseases. However, the identification of vertebrate ohnologs has been challenging, due to lineage specific genome rearrangements since 2R- and 3R-WGD. We previously identified vertebrate ohnologs using a novel synteny comparison across multiple genomes. Here, we refine and apply this approach on 27 vertebrate genomes to identify ohnologs from both 2R- and 3R-WGD, while taking into account the phylogenetically biased sampling of available species. We assemble vertebrate ohnolog pairs and families in an expanded OHNOLOGS v2 database. We find that teleost fish have retained more 2R-WGD ohnologs than mammals and sauropsids, and that these 2R-ohnologs have retained significantly more ohnologs from the subsequent 3R-WGD than genes without 2R-ohnologs. Interestingly, species with fewer extant genes, such as sauropsids, have retained similar or higher proportions of ohnologs. OHNOLOGS v2 should allow deeper evolutionary genomic analysis of the impact of WGD on vertebrates and can be freely accessed at http://ohnologs.curie.fr.

Epithelial RNase H2 Maintains Genome Integrity and Prevents Intestinal Tumorigenesis in Mice

Background & AimsRNase H2 is a holoenzyme, composed of 3 subunits (ribonuclease H2 subunits A, B, and C), that cleaves RNA:DNA hybrids and removes mis-incorporated ribonucleotides from genomic DNA through ribonucleotide excision repair. Ribonucleotide incorporation by eukaryotic DNA polymerases occurs during every round of genome duplication and produces the most frequent type of naturally occurring DNA lesion. We investigated whether intestinal epithelial proliferation requires RNase H2 function and whether RNase H2 activity is disrupted during intestinal carcinogenesis.MethodsWe generated mice with epithelial-specific deletion of ribonuclease H2 subunit B (H2bΔIEC) and mice that also had deletion of tumor-suppressor protein p53 (H2b/p53ΔIEC); we compared phenotypes with those of littermate H2bfl/fl or H2b/p53fl/fl (control) mice at young and old ages. Intestinal tissues were collected and analyzed by histology. We isolated epithelial cells, generated intestinal organoids, and performed RNA sequence analyses. Mutation signatures of spontaneous tumors from H2b/p53ΔIEC mice were characterized by exome sequencing. We collected colorectal tumor specimens from 467 patients, measured levels of ribonuclease H2 subunit B, and associated these with patient survival times and transcriptome data.ResultsThe H2bΔIEC mice had DNA damage to intestinal epithelial cells and proliferative exhaustion of the intestinal stem cell compartment compared with controls and H2b/p53ΔIEC mice. However, H2b/p53ΔIEC mice spontaneously developed small intestine and colon carcinomas. DNA from these tumors contained T>G base substitutions at GTG trinucleotides. Analyses of transcriptomes of human colorectal tumors associated lower levels of RNase H2 with shorter survival times.ConclusionsIn analyses of mice with disruption of the ribonuclease H2 subunit B gene and colorectal tumors from patients, we provide evidence that RNase H2 functions as a colorectal tumor suppressor. H2b/p53ΔIEC mice can be used to study the roles of RNase H2 in tissue-specific carcinogenesis.

Duplication and gene expression patterns of β-catenin in Nile tilapia.

β-Catenin, a key transcriptional coactivator of the Wnt pathway, plays an important role in animal embryonic development and organogenesis. In our earlier study, we have reported that two types of β-catenin (β-catenin-1 and β-catenin-2) were ubiquitously expressed in almost all the tissues examined in tilapia. However, the immunolocalization of β-catenin in those tissues, especially in extra-gonadal tissues, remains unclear. In the present study, we further confirm that these two types of β-catenin gene exist only in teleosts and are derived from 3R (third round of genome duplication) by phylogenetic and syntenic analyses. Moreover, the transcriptome analysis conducted in this investigation reveals that two β-catenins exhibited similar expression patterns in seven adult tissues and four key developmental stages of XX and XY gonads. Finally, immunohistochemistry analysis was performed to detect the cell localization of β-catenin. A positive signal of β-catenin was observed in various tissues of tilapia, such as the intestine, liver, kidney, spleen, eye, brain, and gonads. The results of our study indicate that tilapia β-catenin might be involved in the organ development and play some specific functions in biological processes; all these data will provide basic reference for understanding the molecular mechanism of β-catenin in regulating of teleost organogenesis.

MHC and adaptive immunity in teleost fishes.

The adaptive immune system has long been considered a key evolutionary innovation of the vertebrates, the product of two rounds of genome duplication that gave rise to the raw material necessary for the evolution of a highly specific immune response and immune memory. While comparative studies of a small number of model organisms have led to the commonly held view that the adaptive immune system has remained relatively static since its origin, recent studies of non-model organisms are challenging this notion, highlighting the fact that we have only begun to scratch the surface in terms of our understanding of immune system diversity. Some of the most exciting recent results have come from the comparative analysis of teleost fishes, a group that includes more than 40% of vertebrates, and shows remarkable diversity in immune system structure and function. Despite the repeated loss of key components of the adaptive immune machinery in this group, affected species are capable of mounting a robust response to immune challenge, suggesting that they have evolved alternative mechanisms of immune protection. Such deviations from the canonical model of vertebrate immunity create opportunities to explore common paradigms of immune function, and may contribute to new experimental approaches and methods of treatment.

Characterization of the γ-aminobutyric acid signaling system in the zebrafish (Danio rerio Hamilton) central nervous system by reverse transcription-quantitative polymerase chain reaction

In the vertebrate brain, inhibition is largely mediated by γ-aminobutyric acid (GABA). This neurotransmitter comprises a signaling machinery of GABAA, GABAB receptors, transporters, glutamate decarboxylases (gads) and 4-aminobutyrate aminotransferase (abat), and associated proteins. Chloride is intimately related to GABAA receptor conductance, GABA uptake, and GADs activity. The response of target neurons to GABA stimuli is shaped by chloride-cation co-transporters (CCCs), which strictly control Cl− gradient across plasma membranes. This research profiled the expression of forty genes involved in GABA signaling in the zebrafish (Danio rerio) brain, grouped brain regions and retinas. Primer pairs were developed for reverse transcription-quantitative polymerase chain reaction (RT-qPCR). The mRNA levels of the zebrafish GABA system share similarities with that of mammals, and confirm previous studies in non-mammalian species. Proposed GABAA receptors are α1β2γ2, α1β2δ, α2bβ3γ2, α2bβ3δ, α4β2γ2, α4β2δ, α6bβ2γ2 and α6bβ2δ. Regional brain differences were documented. Retinal hetero- or homomeric ρ-composed GABAA receptors could exist, accompanying α1βyγ2, α1βyδ, α6aβyγ2, α6aβyδ. Expression patterns of α6a and α6b were opposite, with the former being more abundant in retinas, the latter in brains. Given the stoichiometry α6wβyγz, α6a- or α6b-containing receptors likely have different regulatory mechanisms. Different gene isoforms could originate after the rounds of genome duplication during teleost evolution. This research depicts that one isoform is generally more abundantly expressed than the other. Such observations also apply to GABAB receptors, GABA transporters, GABA-related enzymes, CCCs and GABAA receptor-associated proteins, whose presence further strengthens the proof of a GABA system in zebrafish.

Most Compositae (Asteraceae) are descendants of a paleohexaploid and all share a paleotetraploid ancestor with the Calyceraceae.

Like many other flowering plants, members of the Compositae (Asteraceae) have a polyploid ancestry. Previous analyses found evidence for an ancient duplication or possibly triplication in the early evolutionary history of the family. We sought to better place this paleopolyploidy in the phylogeny and assess its nature. We sequenced new transcriptomes for Barnadesia, the lineage sister to all other Compositae, and four representatives of closely related families. Using a recently developed algorithm, MAPS, we analyzed nuclear gene family phylogenies for evidence of paleopolyploidy. We found that the previously recognized Compositae paleopolyploidy is also in the ancestry of the Calyceraceae. Our phylogenomic analyses uncovered evidence for a successive second round of genome duplication among all sampled Compositae except Barnadesia. Our analyses of new samples with new tools provide a revised view of paleopolyploidy in the Compositae. Together with results from a high density Lactuca linkage map, our results suggest that the Compositae and Calyceraceae have a common paleotetraploid ancestor and that most Compositae are descendants of a paleohexaploid. Although paleohexaploids have been previously identified, this is the first example where the paleotetraploid and paleohexaploid lineages have survived over tens of millions of years. The complex polyploidy in the ancestry of the Compositae and Calyceraceae represents a unique opportunity to study the long-term evolutionary fates and consequences of different ploidal levels.

Editorial: Marching Toward 100% Whole Genome Sequencing.

Editorial: Marching Toward 100% Whole Genome Sequencing.

Analysis of Yeast Sporulation Efficiency, Spore Viability, and Meiotic Recombination on Solid Medium.

Under conditions of nutrient deprivation, yeast cells initiate a differentiation program in which meiosis is induced and spores are formed. During meiosis, one round of genome duplication is followed by two rounds of chromosome segregation (meiosis I and meiosis II) to generate four haploid nuclei. Meiotic recombination occurs during prophase I. During sporogenesis, each nucleus becomes surrounded by an individual spore wall, and all four haploid spores become contained as a tetrad within an ascus. Important insights into the meiotic function(s) of a gene of interest can be gained by observing the effects of gene mutations on spore viability and viability patterns among tetrads. Moreover, recombination frequencies among viable spores can reveal potential involvement of the gene during meiotic exchange between homologous chromosomes. Here, we describe methods for inducing spore formation on solid medium, determining spore viability, and measuring, via tetrad analysis, frequencies of crossing over and gene conversion as indicators of meiotic chromosome exchange.

Analysis of Recombination and Chromosome Structure during Yeast Meiosis.

Meiosis is a diploid-specific differentiation program that consists of a single round of genome duplication followed by two rounds of chromosome segregation. These events result in halving of the genetic complement, which is a requirement for formation of haploid reproductive cells (i.e., spores in yeast and gametes in animals and plants). During meiosis I, homologous maternal and paternal chromosomes (homologs) pair and separate, whereas sister chromatids remain connected at the centromeres and separate during the second meiotic division. In most organisms, accurate homolog disjunction requires crossovers, which are formed as products of meiotic recombination. For the past two decades, studies of yeast meiosis have provided invaluable insights into evolutionarily conserved mechanisms of meiosis.

Analysis of Meiotic Recombination and Homolog Interaction during Yeast Meiosis.

During meiosis, one round of genome duplication is followed by two rounds of chromosome segregation, resulting in the halving of the genetic complement and the formation of haploid reproductive cells. In most organisms, intimate juxtaposition of homologous chromosomes and homologous recombination during meiotic prophase are required for meiotic success. Here we present a general protocol for visualizing chromosomal proteins and homolog interaction on surface-spread nuclei and a widely used protocol for analyzing meiotic recombination based on an engineered hotspot referred to as HIS4::LEU2.

Marsupials and monotremes possess a novel family of MHC class I genes that is lost from the eutherian lineage

BackgroundMajor histocompatibility complex (MHC) class I genes are found in the genomes of all jawed vertebrates. The evolution of this gene family is closely tied to the evolution of the vertebrate genome. Family members are frequently found in four paralogous regions, which were formed in two rounds of genome duplication in the early vertebrates, but in some species class Is have been subject to additional duplication or translocation, creating additional clusters. The gene family is traditionally grouped into two subtypes: classical MHC class I genes that are usually MHC-linked, highly polymorphic, expressed in a broad range of tissues and present endogenously-derived peptides to cytotoxic T-cells; and non-classical MHC class I genes generally have lower polymorphism, may have tissue-specific expression and have evolved to perform immune-related or non-immune functions. As immune genes can evolve rapidly and are subject to different selection pressure, we hypothesised that there may be divergent, as yet unannotated or uncharacterised class I genes.ResultsApplication of a novel method of sensitive genome searching of available vertebrate genome sequences revealed a new, extensive sub-family of divergent MHC class I genes, denoted as UT, which has not previously been characterized. These class I genes are found in both American and Australian marsupials, and in monotremes, at an evolutionary chromosomal breakpoint, but are not present in non-mammalian genomes and have been lost from the eutherian lineage. We show that UT family members are expressed in the thymus of the gray short-tailed opossum and in other immune tissues of several Australian marsupials. Structural homology modelling shows that the proteins encoded by this family are predicted to have an open, though short, antigen-binding groove.ConclusionsWe have identified a novel sub-family of putatively non-classical MHC class I genes that are specific to marsupials and monotremes. This family was present in the ancestral mammal and is found in extant marsupials and monotremes, but has been lost from the eutherian lineage. The function of this family is as yet unknown, however, their predicted structure may be consistent with presentation of antigens to T-cells.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-015-1745-4) contains supplementary material, which is available to authorized users.

Genome-wide identification, evolution and expression analysis of nuclear receptor superfamily in Nile tilapia, Oreochromis niloticus

The nuclear receptor (NR) superfamily, which is divided into 7 subfamilies, constitutes one of the largest classes of transcription factors. In this study, through comprehensive database search, we identified all NRs (including 4 novel members) from the tilapia (75), common carp (137), zebrafish (73), fugu (73), tetraodon (72), stickleback (70), medaka (69), coelacanth (55), spotted gar (51) and elephant shark (50). For 21 NRs, two duplicates were found in teleosts, while only one in tetrapods. These duplicates, except those of DAX1, SHP and GCNF found in the elephant shark, were derived from 3R (third round of genome duplication). The linkage duplication of 5 syntenic blocks (comprising 14 duplicated NR couples) in teleosts further supported their 3R origin. Based on transcriptome data from adult tilapia, 53 NRs were found to be expressed in more than one tissue (brain, head kidney, heart, liver, kidney, muscle, ovary and testis), and 4 were tissue-specific, indicating their essential roles in the corresponding tissue. Based on the XX and XY gonadal transcriptome data from four developmental stages, 65 NRs were detected in gonads, with 21, 31, 11 and 29 expressed sexual dimorphically at 5, 30, 90 and 180days after hatching, respectively. The expression of four selected genes was examined by in situ hybridization (ISH) and quantitative PCR (qPCR) to validate the spatial and temporal expression profiles of NRs. Comparative analyses of the expression profiles of duplicated NRs revealed divergence in gene expression as well as gene function. Our results demonstrated that NRs may play important roles in sex determination and gonadal development in teleosts.