Molecular mechanisms of interleukin-2 gene inducibility: developmental control and combinatorial action of transcription factors

Interleukin-2 is one of the lymphokines secreted by T helper type 1 cells upon activation mediated by T-cell receptor (TCR) and accessory molecules. The ability to express IL-2 is correlated with T-lineage commitment and is regulated during T cell development and differentiation. Understanding the molecular mechanism of how IL-2 gene inducibility is controlled at each transition and each differentiation process of T-cell development is to understand one aspect of T-cell development. In the present study, we first attempted to elucidate the molecular basis for the developmental changes of IL-2 gene inducibility. We showed that IL-2 gene inducibility is acquired early in immature CD4- CD8-TCR- thymocytes prior to TCR gene rearrangement. Similar to mature T cells, a complete set of transcription factors can be induced at this early stage to activate IL-2 gene expression. The progression of these cells to cortical CD4^+CD8^+TCR^(1o) cells is accompanied by the loss of IL-2 gene inducibility. We demonstrated that DNA binding activities of two transcription factors AP-1 and NF-AT are reduced in cells at this stage. Further, the loss of factor binding, especially AP-1, is attributable to the reduced ability to activate expression of three potential components of AP-1 and NF-AT, including c-Fos, FosB, and Fra-2. We next examined the interaction of transcription factors and the IL-2 promoter in vivo by using the EL4 T cell line and two non-T cell lines. We showed an all-or-none phenomenon regarding the factor-DNA interaction, i.e., in activated T cells, the IL-2 promoter is occupied by sequence-specific transcription factors when all the transcription factors are available; in resting T cells or non-T cells, no specific protein-DNA interaction is observed when only a subset of factors are present in the nuclei. Purposefully reducing a particular set of factor binding activities in stimulated T cells using pharmacological agents cyclosporin A or forskolin also abolished all interactions. The results suggest that a combinatorial and coordinated protein-DNA interaction is required for IL-2 gene activation. The thymocyte experiments clearly illustrated that multiple transcription factors are regulated during intrathymic T-cell development, and this regulation in tum controls the inducibility of the lineage-specific IL-2 gene. The in vivo study of protein-DNA interaction stressed the combinatorial action of transcription factors to stably occupy the IL-2 promoter and to initiate its transcription, and provided a molecular mechanism for changes in IL-2 gene inducibility in T cells undergoing integration of multiple environmental signals.

Proteomic analysis of the Cdc48/Ubx network identifies a role for Ubx2 in the regulation of lipid biosynthesis

Cdc48/p97 is an essential, highly abundant hexameric member of the AAA (ATPase associated with various cellular activities) family. It has been linked to a variety of processes throughout the cell but it is best known for its role in the ubiquitin proteasome pathway. In this system it is believed that Cdc48 behaves as a segregase, transducing the chemical energy of ATP hydrolysis into mechanical force to separate ubiquitin-conjugated proteins from their tightly-bound partners.

Current models posit that Cdc48 is linked to its substrates through a variety of adaptor proteins, including a family of seven proteins (13 in humans) that contain a Cdc48-binding UBX domain. As such, due to the complexity of the network of adaptor proteins for which it serves as the hub, Cdc48/p97 has the potential to exert a profound influence on the ubiquitin proteasome pathway. However, the number of known substrates of Cdc48/p97 remains relatively small, and smaller still is the number of substrates that have been linked to a specific UBX domain protein. As such, the goal of this dissertation research has been to discover new substrates and better understand the functions of the Cdc48 network. With this objective in mind, we established a proteomic screen to assemble a catalog of candidate substrate/targets of the Ubx adaptor system.

Here we describe the implementation and optimization of a cutting-edge quantitative mass spectrometry method to measure relative changes in the Saccharomyces cerevisiae proteome. Utilizing this technology, and in order to better understand the breadth of function of Cdc48 and its adaptors, we then performed a global screen to identify accumulating ubiquitin conjugates in cdc48-3 and ubxΔ mutants. In this screen different ubx mutants exhibited reproducible patterns of conjugate accumulation that differed greatly from each other, pointing to various unexpected functional specializations of the individual Ubx proteins.

As validation of our mass spectrometry findings, we then examined in detail the endoplasmic-reticulum bound transcription factor Spt23, which we identified as a putative Ubx2 substrate. In these studies ubx2Δ cells were deficient in processing of Spt23 to its active p90 form, and in localizing p90 to the nucleus. Additionally, consistent with reduced processing of Spt23, ubx2Δ cells demonstrated a defect in expression of their target gene OLE1, a fatty acid desaturase. Overall, this work demonstrates the power of proteomics as a tool to identify new targets of various pathways and reveals Ubx2 as a key regulator lipid membrane biosynthesis.

Functional genomic studies of the structure and regulation of eukaryotic transcriptomes

The main focus of this thesis is the use of high-throughput sequencing technologies in functional genomics (in particular in the form of ChIP-seq, chromatin immunoprecipitation coupled with sequencing, and RNA-seq) and the study of the structure and regulation of transcriptomes. Some parts of it are of a more methodological nature while others describe the application of these functional genomic tools to address various biological problems. A significant part of the research presented here was conducted as part of the ENCODE (ENCyclopedia Of DNA Elements) Project.

The first part of the thesis focuses on the structure and diversity of the human transcriptome. Chapter 1 contains an analysis of the diversity of the human polyadenylated transcriptome based on RNA-seq data generated for the ENCODE Project. Chapter 2 presents a simulation-based examination of the performance of some of the most popular computational tools used to assemble and quantify transcriptomes. Chapter 3 includes a study of variation in gene expression, alternative splicing and allelic expression bias on the single-cell level and on a genome-wide scale in human lymphoblastoid cells; it also brings forward a number of critical to the practice of single-cell RNA-seq measurements methodological considerations.

The second part presents several studies applying functional genomic tools to the study of the regulatory biology of organellar genomes, primarily in mammals but also in plants. Chapter 5 contains an analysis of the occupancy of the human mitochondrial genome by TFAM, an important structural and regulatory protein in mitochondria, using ChIP-seq. In Chapter 6, the mitochondrial DNA occupancy of the TFB2M transcriptional regulator, the MTERF termination factor, and the mitochondrial RNA and DNA polymerases is characterized. Chapter 7 consists of an investigation into the curious phenomenon of the physical association of nuclear transcription factors with mitochondrial DNA, based on the diverse collections of transcription factor ChIP-seq datasets generated by the ENCODE, mouseENCODE and modENCODE consortia. In Chapter 8 this line of research is further extended to existing publicly available ChIP-seq datasets in plants and their mitochondrial and plastid genomes.

The third part is dedicated to the analytical and experimental practice of ChIP-seq. As part of the ENCODE Project, a set of metrics for assessing the quality of ChIP-seq experiments was developed, and the results of this activity are presented in Chapter 9. These metrics were later used to carry out a global analysis of ChIP-seq quality in the published literature (Chapter 10). In Chapter 11, the development and initial application of an automated robotic ChIP-seq (in which these metrics also played a major role) is presented.

The fourth part presents the results of some additional projects the author has been involved in, including the study of the role of the Piwi protein in the transcriptional regulation of transposon expression in Drosophila (Chapter 12), and the use of single-cell RNA-seq to characterize the heterogeneity of gene expression during cellular reprogramming (Chapter 13).

The last part of the thesis provides a review of the results of the ENCODE Project and the interpretation of the complexity of the biochemical activity exhibited by mammalian genomes that they have revealed (Chapters 15 and 16), an overview of the expected in the near future technical developments and their impact on the field of functional genomics (Chapter 14), and a discussion of some so far insufficiently explored research areas, the future study of which will, in the opinion of the author, provide deep insights into many fundamental but not yet completely answered questions about the transcriptional biology of eukaryotes and its regulation.

Cell-targeted regulation of gene expression through synthetic RNA devices

The ability to interface with and program cellular function remains a challenging research frontier in biotechnology. Although the emerging field of synthetic biology has recently generated a variety of gene-regulatory strategies based on synthetic RNA molecules, few strategies exist through which to control such regulatory effects in response to specific exogenous or endogenous molecular signals. Here, we present the development of an engineered RNA-based device platform to detect and act on endogenous protein signals, linking these signals to the regulation of genes and thus cellular function.

We describe efforts to develop an RNA-based device framework for regulating endogenous genes in human cells. Previously developed RNA control devices have demonstrated programmable ligand-responsive genetic regulation in diverse cell types, and we attempted to adapt this class of cis-acting control elements to function in trans. We divided the device into two strands that reconstitute activity upon hybridization. Device function was optimized using an in vivo model system, and we found that device sequence is not as flexible as previously reported. After verifying the in vitro activity of our optimized design, we attempted to establish gene regulation in a human cell line using additional elements to direct device stability, structure, and localization. The significant limitations of our platform prevented endogenous gene regulation.

We next describe the development of a protein-responsive RNA-based regulatory platform. Employing various design strategies, we demonstrated functional devices that both up- and downregulate gene expression in response to a heterologous protein in a human cell line. The activity of our platform exceeded that of a similar, small-molecule-responsive platform. We demonstrated the ability of our devices to respond to both cytoplasmic- and nuclear-localized protein, providing insight into the mechanism of action and distinguishing our platform from previously described devices with more restrictive ligand localization requirements. Finally, we demonstrated the versatility of our device platform by developing a regulatory device that responds to an endogenous signaling protein.

The foundational tool we present here possesses unique advantages over previously described RNA-based gene-regulatory platforms. This genetically encoded technology may find future applications in the development of more effective diagnostic tools and targeted molecular therapy strategies.

Characterizing the regulation of mitochondrial nucleoids

Mitochondria contain a 16.6 kb circular genome encoding 13 proteins as well as mitochondrial tRNAs and rRNAs. Copies of the genome are organized into nucleoids containing both DNA and proteins, including the machinery required for mtDNA replication and transcription. Although mtDNA integrity is essential for cellular and organismal viability, regulation of proliferation of the mitochondrial genome is poorly understood. To elucidate the mechanisms behind this, we chose to study the interplay between mtDNA copy number and the proteins involved in mitochondrial fusion, another required function in cells. Strikingly, we found that mouse embryonic fibroblasts lacking fusion also had a mtDNA copy number deficit. To understand this phenomenon further, we analyzed the binding of mitochondrial transcription factor A, whose role in transcription, replication, and packaging of the genome is well-established and crucial for cellular maintenance. Using ChIP-seq, we were able to detect largely uniform, non-specific binding across the genome, with no occupancy in the known specific binding sites in the regulatory region. We did detect a single binding site directly upstream of a known origin of replication, suggesting that TFAM may play a direct role in replication. Finally, although TFAM has been previously shown to localize to the nuclear genome, we found no evidence for such binding sites in our system.

To further understand the regulation of mtDNA by other proteins, we analyzed publicly available ChIP-seq datasets from ENCODE, modENCODE, and mouseENCODE for evidence of nuclear transcription factor binding to the mitochondrial genome. We identified eight human transcription factors and three mouse transcription factors that demonstrated binding events with the classical strand asymmetrical morphology of classical binding sites. ChIP-seq is a powerful tool for understanding the interactions between proteins and the mitochondrial genome, and future studies promise to further the understanding of how mtDNA is regulated within the nucleoid.

Sequence-function relationships in E. coli transcriptional regulation

Understanding how transcriptional regulatory sequence maps to regulatory function remains a difficult problem in regulatory biology. Given a particular DNA sequence for a bacterial promoter region, we would like to be able to say which transcription factors bind there, how strongly they bind, and whether they interact with each other and/or RNA polymerase, with the ultimate objective of integrating knowledge of these parameters into a prediction of gene expression levels. The theoretical framework of statistical thermodynamics provides a useful framework for doing so, enabling us to predict how gene expression levels depend on transcription factor binding energies and concentrations. We used thermodynamic models, coupled with models of the sequence-dependent binding energies of transcription factors and RNAP, to construct a genotype to phenotype map for the level of repression exhibited by the lac promoter, and tested it experimentally using a set of promoter variants from E. coli strains isolated from different natural environments. For this work, we sought to ``reverse engineer'' naturally occurring promoter sequences to understand how variations in promoter sequence affects gene expression. The natural inverse of this approach is to ``forward engineer'' promoter sequences to obtain targeted levels of gene expression. We used a high precision model of RNAP-DNA sequence dependent binding energy, coupled with a thermodynamic model relating binding energy to gene expression, to predictively design and verify a suite of synthetic E. coli promoters whose expression varied over nearly three orders of magnitude.

However, although thermodynamic models enable predictions of mean levels of gene expression, it has become evident that cell-to-cell variability or ``noise'' in gene expression can also play a biologically important role. In order to address this aspect of gene regulation, we developed models based on the chemical master equation framework and used them to explore the noise properties of a number of common E. coli regulatory motifs; these properties included the dependence of the noise on parameters such as transcription factor binding strength and copy number. We then performed experiments in which these parameters were systematically varied and measured the level of variability using mRNA FISH. The results showed a clear dependence of the noise on these parameters, in accord with model predictions.

Finally, one shortcoming of the preceding modeling frameworks is that their applicability is largely limited to systems that are already well-characterized, such as the lac promoter. Motivated by this fact, we used a high throughput promoter mutagenesis assay called Sort-Seq to explore the completely uncharacterized transcriptional regulatory DNA of the E. coli mechanosensitive channel of large conductance (MscL). We identified several candidate transcription factor binding sites, and work is continuing to identify the associated proteins.

Transcription factor occupancy in differentiating skeletal muscle

With recent advances in high-throughput sequencing, mapping of genome-wide transcription factor occupancy has become feasible. To advance the understanding of skeletal muscle differentiation specifically and transcriptional regulation in general, I determined the genome-wide occupancy map for myogenin in differentiating C2C12 myocyte cells. I then analyzed the myogenin map for underlying sequence content and the association between occupied elements and expression trajectories of adjacent genes. Having determined that myogenin primarily associates with expressed genes, I performed a similar analysis on occupancy maps of other transcription factors active during skeletal muscle differentiation, including an extensive analysis of co-occupancy. This analysis provided strong motif evidence for protein-protein interactions as the primary driving force in the formation of Myogenin / Mef2 and MyoD / AP-1 complexes at jointly-occupied sites. Finally, factor occupancy analysis was extended to include bHLH transcription factors in tissues other than skeletal muscle. The cross-tissue analysis led to the emergence of a motif structure used by bHLH TFs to encode either tissue-specific or "general" (public) access in a variety of lineages.

Regulation of the fitness of Atlantic salmon (Salmo salar) by intra-specific competition amongst the juveniles

Factors affecting the fitness of juvenile salmon are discussed. Although fitness from the genetic point of view is defined as the relative capacity of carriers of a given genotype to transmit their genes to the gene pool of the following generations, growth and survival of individuals are also components of fitness, and are influenced by responses to competition, which is the major topic of this article including implications for management. In order to better understand the relationships of density-dependent survival in Newfoundland, egg depositions were manipulated experimentally in the Freshwater River. Figures demonstrate the relationship between stock (number of eggs per 100 m2 of river) and recruitment (number of smolts per l00 m2 of Atlantic salmon, and also the percentage survival from egg to smolt stage related to potential egg depositions.

The E2F2 Transcription Factor Sustains Hepatic Glycerophospholipid Homeostasis in Mice

Increasing evidence links metabolic signals to cell proliferation, but the molecular wiring that connects the two core machineries remains largely unknown. E2Fs are master regulators of cellular proliferation. We have recently shown that E2F2 activity facilitates the completion of liver regeneration after partial hepatectomy (PH) by regulating the expression of genes required for S-phase entry. Our study also revealed that E2F2 determines the duration of hepatectomy-induced hepatic steatosis. A transcriptomic analysis of normal adult liver identified "lipid metabolism regulation" as a major E2F2 functional target, suggesting that E2F2 has a role in lipid homeostasis. Here we use wild-type (E2F2(+/+)) and E2F2 deficient (E2F2(-/-)) mice to investigate the in vivo role of E2F2 in the composition of liver lipids and fatty acids in two metabolically different contexts: quiescence and 48-h post-PH, when cellular proliferation and anabolic demands are maximal. We show that liver regeneration is accompanied by large triglyceride and protein increases without changes in total phospholipids both in E2F2(+/+) and E2F2(-/-) mice. Remarkably, we found that the phenotype of quiescent liver tissue from E2F2(-/-) mice resembles the phenotype of proliferating E2F2(+/+) liver tissue, characterized by a decreased phosphatidylcholine to phosphatidylethanolamine ratio and a reprogramming of genes involved in generation of choline and ethanolamine derivatives. The diversity of fatty acids in total lipid, triglycerides and phospholipids was essentially preserved on E2F2 loss both in proliferating and non-proliferating liver tissue, although notable exceptions in inflammation-related fatty acids of defined phospholipid classes were detected. Overall, our results indicate that E2F2 activity sustains the hepatic homeostasis of major membrane glycerolipid components while it is dispensable for storage glycerolipid balance.

Exploring Genetic Factors Involved in Huntington Disease Age of Onset: E2F2 as a New Potential Modifier Gene

Age of onset (AO) of Huntington disease (HD) is mainly determined by the length of the CAG repeat expansion (CAGexp) in exon 1 of the HTT gene. Additional genetic variation has been suggested to contribute to AO, although the mechanism by which it could affect AO is presently unknown. The aim of this study is to explore the contribution of candidate genetic factors to HD AO in order to gain insight into the pathogenic mechanisms underlying this disorder. For that purpose, two AO definitions were used: the earliest age with unequivocal signs of HD (earliest AO or eAO), and the first motor symptoms age (motor AO or mAO). Multiple linear regression analyses were performed between genetic variation within 20 candidate genes and eAO or mAO, using DNA and clinical information of 253 HD patients from REGISTRY project. Gene expression analyses were carried out by RT-qPCR with an independent sample of 35 HD patients from Basque Country Hospitals. We found suggestive association signals between HD eAO and/or mAO and genetic variation within the E2F2, ATF7IP, GRIN2A, GRIN2B, LINC01559, HIP1 and GRIK2 genes. Among them, the most significant was the association between eAO and rs2742976, mapping to the promoter region of E2F2 transcription factor. Furthermore, rs2742976 T allele patient carriers exhibited significantly lower lymphocyte E2F2 gene expression, suggesting a possible implication of E2F2-dependent transcriptional activity in HD pathogenesis. Thus, E2F2 emerges as a new potential HD AO modifier factor.

Global analysis of mRNA half-lives and de novo transcription in a dinoflagellate, Karenia brevis

Dinoflagellates possess many physiological processes that appear to be under post-transcriptional control. However, the extent to which their genes are regulated post-transcriptionally remains unresolved. To gain insight into the roles of differential mRNA stability and de novo transcription in dinoflagellates, we biosynthetically labeled RNA with 4-thiouracil to isolate newly transcribed and pre-existing RNA pools in Karenia brevis. These isolated fractions were then used for analysis of global mRNA stability and de novo transcription by hybridization to a K. brevis microarray. Global K. brevis mRNA half-lives were calculated from the ratio of newly transcribed to pre-existing RNA for 7086 array features using the online software HALO (Half-life Organizer). Overall, mRNA half-lives were substantially longer than reported in other organisms studied at the global level, ranging from 42 minutes to greater than 144 h, with a median of 33 hours. Consistent with well-documented trends observed in other organisms, housekeeping processes, including energy metabolism and transport, were significantly enriched in the most highly stable messages. Shorter-lived transcripts included a higher proportion of transcriptional regulation, stress response, and other response/regulatory processes. One such family of proteins involved in post-transcriptional regulation in chloroplasts and mitochondria, the pentatricopeptide repeat (PPR) proteins, had dramatically shorter half-lives when compared to the arrayed transcriptome. As transcript abundances for PPR proteins were previously observed to rapidly increase in response to nutrient addition, we queried the newly synthesized RNA pools at 1 and 4 h following nitrate addition to N-depleted cultures. Transcriptome-wide there was little evidence of increases in the rate of de novo transcription during the first 4 h, relative to that in N-depleted cells, and no evidence for increased PPR protein transcription. These results lend support to the growing consensus of post-transcriptional control of gene expression in dinoflagellates.

A de novo originated gene depresses budding yeast mating pathway and is repressed by the protein encoded by its antisense strand

Recent transcription profiling studies have revealed an unexpectedly large proportion of antisense transcripts in eukaryotic genomes. These antisense genes seem to regulate gene expression by interacting with sense genes. Previous studies have focused on the non-coding antisense genes, but the possible regulatory role of the antisense protein is poorly understood. In this study, we found that a protein encoded by the antisense gene ADF1 acts as a transcription suppressor, regulating the expression of sense gene MDF1 in Saccharomyces cerevisiae. Based on the evolutionary, genetic, cytological and biochemical evidence, we show that the protein-coding sense gene MDF1 most likely originated de novo from a previously non-coding sequence and can significantly suppress the mating efficiency of baker's yeast in rich medium by binding MAT alpha 2 and thus promote vegetative growth. These results shed new light on several important issues, including a new sense-antisense interaction mechanism, the de novo origination of a functional gene, and the regulation of yeast mating pathway.

Transcription rates of yeast genes are influenced by the distribution of introns

A comparative analysis on the intron sequence oligonucleotide usages in two sets of yeast genes with higher and lower transcription frequencies, respectively, has shown that the intron sequence structures of the two sets of genes are different. There are more potential binding sites for transcription factors in the introns of the genes with high transcription frequencies. So it is speculated that introns regulate the transcription of genes. But more evidences are needed to favor this speculation. The detailed comparative analyses on the distribution ( length and position) of introns and exons in the two sets of gene sequences also show that there is an obvious boundary between the lengths of the two sets of introns. There is no boundary between the lengths of the two sets of exons, although the means of their lengths are of discrepancy. The situation of the gene lengths ( length of intron and exon) is similar to exon lengths. As far as the relative position, the introns in two sets of genes all have a bias toward the 5' ends of genes. But as the actual position is considered, more introns in high transcription genes have a tendency to be located toward the 5' ends of genes, some even located at 5'-UTR. These results suggest that the gene transcription rates are related to the length of intron, but not to the lengths of exons and genes sequences. The positions of introns may also influence the transcription rates. The transcriptional regulation of introns may be correlative with the transcriptional regulation of the upstream of genes, or be its continuous action.

Statistical analysis of sequence features of introns with positive transcriptional regulation in yeast genes

A great deal of experimental studies have shown that many introns of eukaryotic genes function as regulators of transcription. However, comprehensive studies of this problem have not yet been conducted. After checking the transcription frequencies of some Saccharomyces cerevisiae (yeast), genes and their introns, a remarkable phenomenon was discovered that generally the introns of the genes with higher transcription frequencies are longer, and the introns of the genes with lower transcription frequencies are shorter. This suggests that the longer introns of genes with higher transcription frequencies may contain some characteristic sequence structures, which could enhance the transcription of genes. Therefore, two sets of introns of yeast genes were chosen for further study. The transcription frequencies of the first set of genes are higher (>30), and those of the second set of genes are lower (less than or equal to10). Some oligonucleotides are detected by statistically comparative analyses of the occurrence frequencies of oligonucleotides (mainly tetranucleotides and pentanucleotides), whose occurrence frequencies in the first set of introns; are significantly higher than those in the second set of introns, and are also significantly higher than those in the exons flanking the introns of the first set. Some of these extracted oligonucleotides are the same as the regulatory elements of transcription revealed by experimental analyses. Besides, the distributions of these extracted oligonucleotides in the two sets of introns and the exons show that the sequence structures of the first set of introns are favorable for transcription of genes.

Waterborne exposure to fluorotelomer alcohol 6:2 FTOH alters plasma sex hormone and gene transcription in the hypothalamic-pituitary-gonadal (HPG) axis of zebrafish

Fluorotelomer alcohols (FTOHs) have shown estrogenic activity in vitro and in vivo, but the mechanism of this activity is not known. In this study, 18-week-old zebrafish (Danio rerio) were exposed to 0, 0.03, 0.3 and 3.0 mg/l 1H, 1H, 2H, 2H-perfluorooctan-1-ol (6:2 ETCH) for 7 days, and the effects on plasma sex hormone levels were measured followed by use of real-time PCR to examine selected gene expression in hypothalamic-pituitary-gonadal (HPG) axis and liver. Exposure to 6:2 FTOH significantly increased plasma estradiol (E2) and testosterone (T) levels in both males and females. Furthermore, the ratio of T/E2 was reduced in females while increased in males. In females, the increase of E2 was accompanied by up-regulated hepatic estrogenic receptor alpha (ER alpha) and vitellogenin (VTG1 and VTG3) expression. In males, the elevation of the T level is consistent with the up-regulation of cytochrome P450 c17 alpha-hydroxylase, 17, 20-lase (CYP17) and the down-regulation of cytochrome P450 aromatase A (CYP19A). The present study demonstrated that waterborne exposure to 6:2 FTOH alter plasma sex hormone levels and the ratio of T/E2, as well as the transcriptional profiles of some genes in the HPG axis and liver. The results suggested that FTOHs may disturb fish reproduction through endocrine disrupted activity. (C) 2009 Elsevier B.V. All rights reserved.