Genetic correlates of longevity and selected age-related phenotypes: a genome-wide association study in the Framingham Study

BACKGROUND: Family studies and heritability estimates provide evidence for a genetic contribution to variation in the human life span. METHODS:We conducted a genome wide association study (Affymetrix 100K SNP GeneChip) for longevity-related traits in a community-based sample. We report on 5 longevity and aging traits in up to 1345 Framingham Study participants from 330 families. Multivariable-adjusted residuals were computed using appropriate models (Cox proportional hazards, logistic, or linear regression) and the residuals from these models were used to test for association with qualifying SNPs (70, 987 autosomal SNPs with genotypic call rate [greater than or equal to]80%, minor allele frequency [greater than or equal to]10%, Hardy-Weinberg test p [greater than or equal to] 0.001).RESULTS:In family-based association test (FBAT) models, 8 SNPs in two regions approximately 500 kb apart on chromosome 1 (physical positions 73,091,610 and 73, 527,652) were associated with age at death (p-value < 10-5). The two sets of SNPs were in high linkage disequilibrium (minimum r2 = 0.58). The top 30 SNPs for generalized estimating equation (GEE) tests of association with age at death included rs10507486 (p = 0.0001) and rs4943794 (p = 0.0002), SNPs intronic to FOXO1A, a gene implicated in lifespan extension in animal models. FBAT models identified 7 SNPs and GEE models identified 9 SNPs associated with both age at death and morbidity-free survival at age 65 including rs2374983 near PON1. In the analysis of selected candidate genes, SNP associations (FBAT or GEE p-value < 0.01) were identified for age at death in or near the following genes: FOXO1A, GAPDH, KL, LEPR, PON1, PSEN1, SOD2, and WRN. Top ranked SNP associations in the GEE model for age at natural menopause included rs6910534 (p = 0.00003) near FOXO3a and rs3751591 (p = 0.00006) in CYP19A1. Results of all longevity phenotype-genotype associations for all autosomal SNPs are web posted at http://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?id=phs000007. CONCLUSION: Longevity and aging traits are associated with SNPs on the Affymetrix 100K GeneChip. None of the associations achieved genome-wide significance. These data generate hypotheses and serve as a resource for replication as more genes and biologic pathways are proposed as contributing to longevity and healthy aging.

Characterisation of the microbiota of traditional fermented beverages and screening these and other populations for novel antimicrobial producers and gene clusters

To screen for novel ribosomally synthesised antimicrobials, in-silico genome mining was performed on all publically available fully sequenced bacterial genomes. 49 novel type 1 lantibiotic clusters were identified from a number of species, genera and phyla not usually associated with lantibiotic production, and indicates high prevalence. A crucial step towards the commercialisation of fermented beverages is the characterisation of the microbial content. To achieve this goal, we applied next-generation sequencing techniques to analyse the bacterial and yeast populations of the organic, symbiotically-fermented beverages kefir, water kefir and kombucha. A number of minor components were revealed, many of which had not previously been associated with these beverages. The dominant microorganism in each of the water kefir grains and fermentates was Zymomonas, an ethanol-producing bacterium that had not previously been detected on such a scale. These studies represent the most accurate description of these populations to date, and should aid in future starter design and in determining which species are responsible for specific attributes of the beverages. Finally, high-throughput robotics was applied to screen for the presence of antimicrobial producers associated with these beverages. This revealed a low frequency of bacteriocin production amongst the bacterial isolates, with only lactococcins A, B and LcnN of lactococcin M being identified. However, a proteinaceous antimicrobial produced by the yeast Dekkera bruxellensis, isolated from kombucha, was found to be active against Lactobacillus bulgaricus. This peptide was patially purified.

Hepcidin as a therapeutic tool to limit iron overload and improve anemia in β-thalassemic mice.

Excessive iron absorption is one of the main features of β-thalassemia and can lead to severe morbidity and mortality. Serial analyses of β-thalassemic mice indicate that while hemoglobin levels decrease over time, the concentration of iron in the liver, spleen, and kidneys markedly increases. Iron overload is associated with low levels of hepcidin, a peptide that regulates iron metabolism by triggering degradation of ferroportin, an iron-transport protein localized on absorptive enterocytes as well as hepatocytes and macrophages. Patients with β-thalassemia also have low hepcidin levels. These observations led us to hypothesize that more iron is absorbed in β-thalassemia than is required for erythropoiesis and that increasing the concentration of hepcidin in the body of such patients might be therapeutic, limiting iron overload. Here we demonstrate that a moderate increase in expression of hepcidin in β-thalassemic mice limits iron overload, decreases formation of insoluble membrane-bound globins and reactive oxygen species, and improves anemia. Mice with increased hepcidin expression also demonstrated an increase in the lifespan of their red cells, reversal of ineffective erythropoiesis and splenomegaly, and an increase in total hemoglobin levels. These data led us to suggest that therapeutics that could increase hepcidin levels or act as hepcidin agonists might help treat the abnormal iron absorption in individuals with β-thalassemia and related disorders.

Evolution of the sex-related locus and genomic features shared in microsporidia and fungi.

BACKGROUND: Microsporidia are obligate intracellular, eukaryotic pathogens that infect a wide range of animals from nematodes to humans, and in some cases, protists. The preponderance of evidence as to the origin of the microsporidia reveals a close relationship with the fungi, either within the kingdom or as a sister group to it. Recent phylogenetic studies and gene order analysis suggest that microsporidia share a particularly close evolutionary relationship with the zygomycetes. METHODOLOGY/PRINCIPAL FINDINGS: Here we expanded this analysis and also examined a putative sex-locus for variability between microsporidian populations. Whole genome inspection reveals a unique syntenic gene pair (RPS9-RPL21) present in the vast majority of fungi and the microsporidians but not in other eukaryotic lineages. Two other unique gene fusions (glutamyl-prolyl tRNA synthetase and ubiquitin-ribosomal subunit S30) that are present in metazoans, choanoflagellates, and filasterean opisthokonts are unfused in the fungi and microsporidians. One locus previously found to be conserved in many microsporidian genomes is similar to the sex locus of zygomycetes in gene order and architecture. Both sex-related and sex loci harbor TPT, HMG, and RNA helicase genes forming a syntenic gene cluster. We sequenced and analyzed the sex-related locus in 11 different Encephalitozoon cuniculi isolates and the sibling species E. intestinalis (3 isolates) and E. hellem (1 isolate). There was no evidence for an idiomorphic sex-related locus in this Encephalitozoon species sample. According to sequence-based phylogenetic analyses, the TPT and RNA helicase genes flanking the HMG genes are paralogous rather than orthologous between zygomycetes and microsporidians. CONCLUSION/SIGNIFICANCE: The unique genomic hallmarks between microsporidia and fungi are independent of sequence based phylogenetic comparisons and further contribute to define the borders of the fungal kingdom and support the classification of microsporidia as unusual derived fungi. And the sex/sex-related loci appear to have been subject to frequent gene conversion and translocations in microsporidia and zygomycetes.

Use of 16S ribosomal RNA gene analyses to characterize the bacterial signature associated with poor oral health in West Virginia.

BACKGROUND: West Virginia has the worst oral health in the United States, but the reasons for this are unclear. This pilot study explored the etiology of this disparity using culture-independent analyses to identify bacterial species associated with oral disease. METHODS: Bacteria in subgingival plaque samples from twelve participants in two independent West Virginia dental-related studies were characterized using 16S rRNA gene sequencing and Human Oral Microbe Identification Microarray (HOMIM) analysis. Unifrac analysis was used to characterize phylogenetic differences between bacterial communities obtained from plaque of participants with low or high oral disease, which was further evaluated using clustering and Principal Coordinate Analysis. RESULTS: Statistically different bacterial signatures (P<0.001) were identified in subgingival plaque of individuals with low or high oral disease in West Virginia based on 16S rRNA gene sequencing. Low disease contained a high frequency of Veillonella and Streptococcus, with a moderate number of Capnocytophaga. High disease exhibited substantially increased bacterial diversity and included a large proportion of Clostridiales cluster bacteria (Selenomonas, Eubacterium, Dialister). Phylogenetic trees constructed using 16S rRNA gene sequencing revealed that Clostridiales were repeated colonizers in plaque associated with high oral disease, providing evidence that the oral environment is somehow influencing the bacterial signature linked to disease. CONCLUSIONS: Culture-independent analyses identified an atypical bacterial signature associated with high oral disease in West Virginians and provided evidence that the oral environment influenced this signature. Both findings provide insight into the etiology of the oral disparity in West Virginia.

In vivo ventricular gene delivery of a beta-adrenergic receptor kinase inhibitor to the failing heart reverses cardiac dysfunction.

BACKGROUND: Genetic manipulation to reverse molecular abnormalities associated with dysfunctional myocardium may provide novel treatment. This study aimed to determine the feasibility and functional consequences of in vivo beta-adrenergic receptor kinase (betaARK1) inhibition in a model of chronic left ventricular (LV) dysfunction after myocardial infarction (MI). METHODS AND RESULTS: Rabbits underwent ligation of the left circumflex (LCx) marginal artery and implantation of sonomicrometric crystals. Baseline cardiac physiology was studied 3 weeks after MI; 5x10(11) viral particles of adenovirus was percutaneously delivered through the LCx. Animals received transgenes encoding a peptide inhibitor of betaARK1 (Adeno-betaARKct) or an empty virus (EV) as control. One week after gene delivery, global LV and regional systolic function were measured again to assess gene treatment. Adeno-betaARKct delivery to the failing heart through the LCx resulted in chamber-specific expression of the betaARKct. Baseline in vivo LV systolic performance was improved in Adeno-betaARKct-treated animals compared with their individual pre-gene delivery values and compared with EV-treated rabbits. Total beta-AR density and betaARK1 levels were unchanged between treatment groups; however, beta-AR-stimulated adenylyl cyclase activity in the LV was significantly higher in Adeno-betaARKct-treated rabbits compared with EV-treated animals. CONCLUSIONS: In vivo delivery of Adeno-betaARKct is feasible in the infarcted/failing heart by coronary catheterization; expression of betaARKct results in marked reversal of ventricular dysfunction. Thus, inhibition of betaARK1 provides a novel treatment strategy for improving the cardiac performance of the post-MI heart.

Genetic variants in IGF-I, IGF-II, IGFBP-3, and adiponectin genes and colon cancer risk in African Americans and Whites.

PURPOSE: Evaluating genetic susceptibility may clarify effects of known environmental factors and also identify individuals at high risk. We evaluated the association of four insulin-related pathway gene polymorphisms in insulin-like growth factor-1 (IGF-I) (CA)( n ) repeat, insulin-like growth factor-2 (IGF-II) (rs680), insulin-like growth factor-binding protein-3 (IGFBP-3) (rs2854744), and adiponectin (APM1 rs1501299) with colon cancer risk, as well as relationships with circulating IGF-I, IGF-II, IGFBP-3, and C-peptide in a population-based study. METHODS: Participants were African Americans (231 cases and 306 controls) and Whites (297 cases, 530 controls). Consenting subjects provided blood specimens and lifestyle/diet information. Genotyping for all genes except IGF-I was performed by the 5'-exonuclease (Taqman) assay. The IGF-I (CA)(n) repeat was assayed by PCR and fragment analysis. Circulating proteins were measured by enzyme immunoassays. Odds ratios (ORs) and 95 % confidence intervals (CIs) were calculated by logistic regression. RESULTS: The IGF-I (CA)( 19 ) repeat was higher in White controls (50 %) than African American controls (31 %). Whites homozygous for the IGF-I (CA)(19) repeat had a nearly twofold increase in risk of colon cancer (OR = 1.77; 95 % CI = 1.15-2.73), but not African Americans (OR = 0.73, 95 % CI 0.50-1.51). We observed an inverse association between the IGF-II Apa1 A-variant and colon cancer risk (OR = 0.49, 95 % CI 0.28-0.88) in Whites only. Carrying the IGFBP-3 variant alleles was associated with lower IGFBP-3 protein levels, a difference most pronounced in Whites (p-trend <0.05). CONCLUSIONS: These results support an association between insulin pathway-related genes and elevated colon cancer risk in Whites but not in African Americans.

Expression of a beta-adrenergic receptor kinase 1 inhibitor prevents the development of myocardial failure in gene-targeted mice.

Heart failure is accompanied by severely impaired beta-adrenergic receptor (betaAR) function, which includes loss of betaAR density and functional uncoupling of remaining receptors. An important mechanism for the rapid desensitization of betaAR function is agonist-stimulated receptor phosphorylation by the betaAR kinase (betaARK1), an enzyme known to be elevated in failing human heart tissue. To investigate whether alterations in betaAR function contribute to the development of myocardial failure, transgenic mice with cardiac-restricted overexpression of either a peptide inhibitor of betaARK1 or the beta2AR were mated into a genetic model of murine heart failure (MLP-/-). In vivo cardiac function was assessed by echocardiography and cardiac catheterization. Both MLP-/- and MLP-/-/beta2AR mice had enlarged left ventricular (LV) chambers with significantly reduced fractional shortening and mean velocity of circumferential fiber shortening. In contrast, MLP-/-/betaARKct mice had normal LV chamber size and function. Basal LV contractility in the MLP-/-/betaARKct mice, as measured by LV dP/dtmax, was increased significantly compared with the MLP-/- mice but less than controls. Importantly, heightened betaAR desensitization in the MLP-/- mice, measured in vivo (responsiveness to isoproterenol) and in vitro (isoproterenol-stimulated membrane adenylyl cyclase activity), was completely reversed with overexpression of the betaARK1 inhibitor. We report here the striking finding that overexpression of this inhibitor prevents the development of cardiomyopathy in this murine model of heart failure. These findings implicate abnormal betaAR-G protein coupling in the pathogenesis of the failing heart and point the way toward development of agents to inhibit betaARK1 as a novel mode of therapy.

Potentiation of beta-adrenergic signaling by adenoviral-mediated gene transfer in adult rabbit ventricular myocytes.

Our laboratory has been testing the hypothesis that genetic modulation of the beta-adrenergic signaling cascade can enhance cardiac function. We have previously shown that transgenic mice with cardiac overexpression of either the human beta2-adrenergic receptor (beta2AR) or an inhibitor of the beta-adrenergic receptor kinase (betaARK), an enzyme that phosphorylates and uncouples agonist-bound receptors, have increased myocardial inotropy. We now have created recombinant adenoviruses encoding either the beta2AR (Adeno-beta2AR) or a peptide betaARK inhibitor (consisting of the carboxyl terminus of betaARK1, Adeno-betaARKct) and tested their ability to potentiate beta-adrenergic signaling in cultured adult rabbit ventricular myocytes. As assessed by radioligand binding, Adeno-beta2AR infection led to approximately 20-fold overexpression of beta-adrenergic receptors. Protein immunoblots demonstrated the presence of the Adeno-betaARKct transgene. Both transgenes significantly increased isoproterenol-stimulated cAMP as compared to myocytes infected with an adenovirus encoding beta-galactosidase (Adeno-betaGal) but did not affect the sarcolemmal adenylyl cyclase response to Forskolin or NaF. beta-Adrenergic agonist-induced desensitization was significantly inhibited in Adeno-betaARKct-infected myocytes (16+/-2%) as compared to Adeno-betaGal-infected myocytes (37+/-1%, P < 0.001). We conclude that recombinant adenoviral gene transfer of the beta2AR or an inhibitor of betaARK-mediated desensitization can potentiate beta-adrenergic signaling.

A functional polymorphism in the 5HTR2C gene associated with stress responses also predicts incident cardiovascular events.

Previously we have shown that a functional nonsynonymous single nucleotide polymorphism (rs6318) of the 5HTR2C gene located on the X-chromosome is associated with hypothalamic-pituitary-adrenal axis response to a stress recall task, and with endophenotypes associated with cardiovascular disease (CVD). These findings suggest that individuals carrying the rs6318 Ser23 C allele will be at higher risk for CVD compared to Cys23 G allele carriers. The present study examined allelic variation in rs6318 as a predictor of coronary artery disease (CAD) severity and a composite endpoint of all-cause mortality or myocardial infarction (MI) among Caucasian participants consecutively recruited through the cardiac catheterization laboratory at Duke University Hospital (Durham, NC) as part of the CATHGEN biorepository. Study population consisted of 6,126 Caucasian participants (4,036 [65.9%] males and 2,090 [34.1%] females). A total of 1,769 events occurred (1,544 deaths and 225 MIs; median follow-up time = 5.3 years, interquartile range = 3.3-8.2). Unadjusted Cox time-to-event regression models showed, compared to Cys23 G carriers, males hemizygous for Ser23 C and females homozygous for Ser23C were at increased risk for the composite endpoint of all-cause death or MI: Hazard Ratio (HR) = 1.47, 95% confidence interval (CI) = 1.17, 1.84, p = .0008. Adjusting for age, rs6318 genotype was not related to body mass index, diabetes, hypertension, dyslipidemia, smoking history, number of diseased coronary arteries, or left ventricular ejection fraction in either males or females. After adjustment for these covariates the estimate for the two Ser23 C groups was modestly attenuated, but remained statistically significant: HR = 1.38, 95% CI = 1.10, 1.73, p = .005. These findings suggest that this functional polymorphism of the 5HTR2C gene is associated with increased risk for CVD mortality and morbidity, but this association is apparently not explained by the association of rs6318 with traditional risk factors or conventional markers of atherosclerotic disease.

Bayesian inference of the number of factors in gene-expression analysis: application to human virus challenge studies.

BACKGROUND: Nonparametric Bayesian techniques have been developed recently to extend the sophistication of factor models, allowing one to infer the number of appropriate factors from the observed data. We consider such techniques for sparse factor analysis, with application to gene-expression data from three virus challenge studies. Particular attention is placed on employing the Beta Process (BP), the Indian Buffet Process (IBP), and related sparseness-promoting techniques to infer a proper number of factors. The posterior density function on the model parameters is computed using Gibbs sampling and variational Bayesian (VB) analysis. RESULTS: Time-evolving gene-expression data are considered for respiratory syncytial virus (RSV), Rhino virus, and influenza, using blood samples from healthy human subjects. These data were acquired in three challenge studies, each executed after receiving institutional review board (IRB) approval from Duke University. Comparisons are made between several alternative means of per-forming nonparametric factor analysis on these data, with comparisons as well to sparse-PCA and Penalized Matrix Decomposition (PMD), closely related non-Bayesian approaches. CONCLUSIONS: Applying the Beta Process to the factor scores, or to the singular values of a pseudo-SVD construction, the proposed algorithms infer the number of factors in gene-expression data. For real data the "true" number of factors is unknown; in our simulations we consider a range of noise variances, and the proposed Bayesian models inferred the number of factors accurately relative to other methods in the literature, such as sparse-PCA and PMD. We have also identified a "pan-viral" factor of importance for each of the three viruses considered in this study. We have identified a set of genes associated with this pan-viral factor, of interest for early detection of such viruses based upon the host response, as quantified via gene-expression data.

Rapid behavioral and genomic responses to social opportunity.

From primates to bees, social status regulates reproduction. In the cichlid fish Astatotilapia (Haplochromis) burtoni, subordinate males have reduced fertility and must become dominant to reproduce. This increase in sexual capacity is orchestrated by neurons in the preoptic area, which enlarge in response to dominance and increase expression of gonadotropin-releasing hormone 1 (GnRH1), a peptide critical for reproduction. Using a novel behavioral paradigm, we show for the first time that subordinate males can become dominant within minutes of an opportunity to do so, displaying dramatic changes in body coloration and behavior. We also found that social opportunity induced expression of the immediate-early gene egr-1 in the anterior preoptic area, peaking in regions with high densities of GnRH1 neurons, and not in brain regions that express the related peptides GnRH2 and GnRH3. This genomic response did not occur in stable subordinate or stable dominant males even though stable dominants, like ascending males, displayed dominance behaviors. Moreover, egr-1 in the optic tectum and the cerebellum was similarly induced in all experimental groups, showing that egr-1 induction in the anterior preoptic area of ascending males was specific to this brain region. Because egr-1 codes for a transcription factor important in neural plasticity, induction of egr-1 in the anterior preoptic area by social opportunity could be an early trigger in the molecular cascade that culminates in enhanced fertility and other long-term physiological changes associated with dominance.

Genetic variants in the TEP1 gene are associated with prostate cancer risk and recurrence.

BACKGROUND: Telomere-related genes play an important role in carcinogenesis and progression of prostate cancer (PCa). It is not fully understood whether genetic variations in telomere-related genes are associated with development and progression in PCa patients. METHODS: Six potentially functional single-nucleotide polymorphisms (SNPs) of three key telomere-related genes were evaluated in 1015 PCa cases and 1052 cancer-free controls, to test their associations with risk of PCa. Among 426 PCa patients who underwent radical prostatectomy (RP), the prognostic significance of the studied SNPs on biochemical recurrence (BCR) was also assessed using the Kaplan-Meier analysis and Cox proportional hazards regression model. The relative telomere lengths (RTLs) were measured in peripheral blood leukocytes using real-time PCR in the RP patients. RESULTS: TEP1 rs1760904 AG/AA genotypes were significantly associated with a decreased risk of PCa (odds ratio (OR): 0.77, 95% confidence interval (CI): 0.64-0.93, P=0.005) compared with the GG genotype. By using median RTL as a cutoff level, RP patients with TEP1 rs1760904 AG/AA genotypes tended to have a longer RTL than those with the GG genotype (OR: 1.55, 95% CI: 1.04-2.30, P=0.031). A significant interaction between TEP1 rs1713418 and age in modifying PCa risk was observed (P=0.005). After adjustment for clinicopathologic risk factors, the presence of heterozygotes or rare homozygotes of TEP1 rs1760904 and TNKS2 rs1539042 were associated with BCR in the RP cohorts (hazard ratio: 0.53, 95% CI: 0.36-0.79, P=0.002 and hazard ratio: 1.67, 95% CI: 1.07-2.48, P=0.017, respectively). CONCLUSIONS: These data suggest that genetic variations in the TEP1 gene may be biomarkers for risk of PCa and BCR after RP.

Study of gene expression in thyrotropin-stimulated thyroid cells by cDNA expression array: ID3 transcription modulating factor as an early response protein and tumor marker in thyroid carcinomas.

Induction of cell proliferation by mitogen or growth factor stimulation leads to the specific induction or repression of a large number of genes. To identify genes differentially regulated by the cAMP-dependent transduction pathway, which is poorly characterized so far, we used the cDNA expression array technology. Hybridizations of Atlas human cDNA expression arrays with (32)P-labeled cDNA probes derived from control or thyrotropin (TSH)-stimulated dog thyrocytes in primary culture generated expression profiles of hundreds of genes simultaneously. Among the genes that displayed modified expression, we selected the transcription factor ID3, whose expression was increased by a cAMP-dependent stimulus. ID3 overexpression after TSH stimulation was first verified by Northern blotting analysis, and its mRNA regulation was then investigated in response to a variety of agents acting on thyrocyte proliferation and/or differentiation. We show that: (1) ID3 mRNA induction was stronger after stimulation of the cAMP cascade, but was not restricted to this signaling pathway, as phorbol myristate ester (TPA) and insulin also stimulated mRNA accumulation; (2) in contrast, powerful mitogens for thyroid cells, epidermal growth factor and hepatocyte growth factor, did not significantly modify ID3 mRNA levels; (3) ID3 protein levels closely parallelled mRNA levels, as revealed by immunofluorescence experiments showing a nuclear signal regulated by TSH; (4) in papillary thyroid carcinomas, ID3 mRNA was downregulated. Our results suggest that ID3 expression might be more related to the differentiating process induced by TSH than to the proliferative action of this hormone.

Immediate early gene expression in dog thyrocytes in response to growth, proliferation, and differentiation stimuli.

In dog thyroid cells, insulin or IGF-1 induces cell growth and is required for the mitogenic action of TSH through cyclic AMP, of EGF, and of phorbol esters. HGF per se stimulates cell proliferation and is thus the only full mitogenic agent. TSH and cAMP enhance, whereas EGF phorbol esters and HGF repress differentiation expression. In this study, we have investigated for each factor and regulatory cascade of the intermediate step of immediate early gene induction, that is, c-myc, c-jun, jun D, jun B, c-fos, fos B, fra-1, fra-2, and egr1; fra-1 and fra-2 expressions were very low. TSH or forskolin increased the levels of c-myc, jun B, jun D, c-fos, and fos B while decreasing those of c-jun and egr1. Phorbol myristate ester stimulated the expression of all the genes. EGF and HGF stimulated the expression of all the genes except jun D and for EGF fos B. All these effects were obtained in the presence and in the absence of insulin, which shows that insulin is not necessary for the effects of the mitogens on immediate early gene expression. The definition of the repertoire of early immediate genes inductible by the various growth cascades provides a framework for the analysis of gene expression in tumors. (1) Insulin was able to induce all the protooncogenes investigated except fos B. This suggests that fos B could be the factor missing for insulin to induce mitogenesis. (2) No characteristic pattern of immediate early gene expression has been observed for insulin, which induces cell hypertrophy and is permissive for the action of the other growth factors. These effects are therefore not accounted for by a specific immediate early gene expression. On the other hand, insulin clearly enhances the effects of TSH, phorbol ester, and EGF on c-myc, junB, and c-fos expression. This suggests that the effect of insulin on mitogenesis might result from quantitative differences in the transcription complexes formed. (3) c-myc, c-fos, and jun B mRNA induction by all stimulating agents, whether inducing cell hypertrophy, or growth and dedifferentiation, or growth and differentiation, suggests that, although these expressions are not sufficient, they may be necessary for the various growth responses of thyroid cells. (4) The inhibition of c-jun and egr1 mRNA expression, and the marked induction of jun D mRNA appear to be specific features of the TSH cAMP pathway. They might be related to its differentiating action. (5) fos B, which is induced by TSH, forskolin, phorbol ester, and HGF but not by insulin, could be involved in the mitogenic action of the former factors.