An integrated gene-to-outcome multimodal database for metabolic dysfunction-associated steatotic liver disease.

Metabolic dysfunction-associated steatotic liver disease (MASLD) is the commonest cause of chronic liver disease worldwide and represents an unmet precision medicine challenge. We established a retrospective national cohort of 940 histologically defined patients (55.4% men, 44.6% women; median body mass index 31.3; 32% with type 2 diabetes) covering the complete MASLD severity spectrum, and created a secure, searchable, open resource (SteatoSITE). In 668 cases and 39 controls, we generated hepatic bulk RNA sequencing data and performed differential gene expression and pathway analysis, including exploration of gender-specific differences. A web-based gene browser was also developed. We integrated histopathological assessments, transcriptomic data and 5.67 million days of time-stamped longitudinal electronic health record data to define disease-stage-specific gene expression signatures, pathogenic hepatic cell subpopulations and master regulator networks associated with adverse outcomes in MASLD. We constructed a 15-gene transcriptional risk score to predict future hepatic decompensation events (area under the receiver operating characteristic curve 0.86, 0.81 and 0.83 for 1-, 3- and 5-year risk, respectively). Additionally, thyroid hormone receptor beta regulon activity was identified as a critical suppressor of disease progression. SteatoSITE supports rational biomarker and drug development and facilitates precision medicine approaches for patients with MASLD.

An acidic microenvironment in Tuberculosis increases extracellular matrix degradation by regulating macrophage inflammatory responses.

Mycobacterium tuberculosis (M.tb) infection causes marked tissue inflammation leading to lung destruction and morbidity. The inflammatory extracellular microenvironment is acidic, however the effect of this acidosis on the immune response to M.tb is unknown. Using RNA-seq we show that acidosis produces system level transcriptional change in M.tb infected human macrophages regulating almost 4000 genes. Acidosis specifically upregulated extracellular matrix (ECM) degradation pathways with increased expression of Matrix metalloproteinases (MMPs) which mediate lung destruction in Tuberculosis. Macrophage MMP-1 and -3 secretion was increased by acidosis in a cellular model. Acidosis markedly suppresses several cytokines central to control of M.tb infection including TNF-α and IFN-γ. Murine studies demonstrated expression of known acidosis signaling G-protein coupled receptors OGR-1 and TDAG-8 in Tuberculosis which are shown to mediate the immune effects of decreased pH. Receptors were then demonstrated to be expressed in patients with TB lymphadenitis. Collectively, our findings show that an acidic microenvironment modulates immune function to reduce protective inflammatory responses and increase extracellular matrix degradation in Tuberculosis. Acidosis receptors are therefore potential targets for host directed therapy in patients.

Fur removal promotes an earlier expression of involution-related genes in mammary gland of lactating mice.

Peak lactation occurs when milk production is at its highest. The factors limiting peak lactation performance have been subject of intense debate. Milk production at peak lactation appears limited by the capacity of lactating females to dissipate body heat generated as a by-product of processing food and producing milk. As a result, manipulations that enhance capacity to dissipate body heat (such as fur removal) increase peak milk production. We investigated the potential correlates of shaving-induced increases in peak milk production in laboratory mice. By transcriptomic profiling of the mammary gland, we searched for the mechanisms underlying experimentally increased milk production and its consequences for mother-young conflict over weaning, manifested by advanced or delayed involution of mammary gland. We demonstrated that shaving-induced increases in milk production were paradoxically linked to reduced expression of some milk synthesis-related genes. Moreover, the mammary glands of shaved mice had a gene expression profile indicative of earlier involution relative to unshaved mice. Once provided with enhanced capacity to dissipate body heat, shaved mice were likely to rear their young to independence faster than unshaved mothers.

High neural activity accelerates the decline of cognitive plasticity with age in Caenorhabditis elegans.

The ability to learn progressively declines with age. Neural hyperactivity has been implicated in impairing cognitive plasticity with age, but the molecular mechanisms remain elusive. Here, we show that chronic excitation of the Caenorhabditis elegans O2-sensing neurons during ageing causes a rapid decline of experience-dependent plasticity in response to environmental O2 concentration, whereas sustaining lower activity of O2-sensing neurons retains plasticity with age. We demonstrate that neural activity alters the ageing trajectory in the transcriptome of O2-sensing neurons, and our data suggest that high-activity neurons redirect resources from maintaining plasticity to sustaining continuous firing. Sustaining plasticity with age requires the K+-dependent Na+/Ca2+ (NCKX) exchanger, whereas the decline of plasticity with age in high-activity neurons acts through calmodulin and the scaffold protein Kidins220. Our findings demonstrate directly that the activity of neurons alters neuronal homeostasis to govern the age-related decline of neural plasticity and throw light on the mechanisms involved.

An eQTL in the cystathionine beta synthase gene is linked to osteoporosis in laying hens.

BACKGROUND: Skeletal damage is a challenge for laying hens because the physiological adaptations required for egg laying make them susceptible to osteoporosis. Previously, we showed that genetic factors explain 40% of the variation in end of lay bone quality and we detected a quantitative trait locus (QTL) of large effect on chicken chromosome 1. The aim of this study was to combine data from the commercial founder White Leghorn population and the F2 mapping population to fine-map this QTL and understand its function in terms of gene expression and physiology. RESULTS: Several single nucleotide polymorphisms on chromosome 1 between 104 and 110 Mb (galGal6) had highly significant associations with tibial breaking strength. The alternative genotypes of markers of large effect that flanked the region had tibial breaking strengths of 200.4 vs. 218.1 Newton (P < 0.002) and, in a subsequent founder generation, the higher breaking strength genotype was again associated with higher breaking strength. In a subsequent generation, cortical bone density and volume were increased in individuals with the better bone genotype but with significantly reduced medullary bone quality. The effects on cortical bone density were confirmed in a further generation and was accompanied by increased mineral maturity of the cortical bone as measured by infrared spectrometry and there was evidence of better collagen cross-linking in the cortical bone. Comparing the transcriptome of the tibia from individuals with good or poor bone quality genotypes indicated four differentially-expressed genes at the locus, one gene, cystathionine beta synthase (CBS), having a nine-fold higher expression in the genotype for low bone quality. The mechanism was cis-acting and although there was an amino-acid difference in the CBS protein between the genotypes, there was no difference in the activity of the enzyme. Plasma homocysteine concentration, the substrate of CBS, was higher in the poor bone quality genotype. CONCLUSIONS: Validated markers that predict bone strength have been defined for selective breeding and a gene was identified that may suggest alternative ways to improve bone health in addition to genetic selection. The identification of how genetic variants affect different aspects of bone turnover shows potential for translational medicine.

Lack of long-term acclimation in Antarctic encrusting species suggests vulnerability to warming.

Marine encrusting communities play vital roles in benthic ecosystems and have major economic implications with regards to biofouling. However, their ability to persist under projected warming scenarios remains poorly understood and is difficult to study under realistic conditions. Here, using heated settlement panel technologies, we show that after 18 months Antarctic encrusting communities do not acclimate to either +1 °C or +2 °C above ambient temperatures. There is significant up-regulation of the cellular stress response in warmed animals, their upper lethal temperatures decline with increasing ambient temperature and population genetic analyses show little evidence of differential survival of genotypes with treatment. By contrast, biofilm bacterial communities show no significant differences in community structure with temperature. Thus, metazoan and bacterial responses differ dramatically, suggesting that ecosystem responses to future climate change are likely to be far more complex than previously anticipated.

Dissection of Ramularia Leaf Spot Disease by Integrated Analysis of Barley and Ramularia collo-cygni Transcriptome Responses.

Ramularia leaf spot disease (RLS), caused by the ascomycete fungus Ramularia collo-cygni, has emerged as a major economic disease of barley. No substantial resistance has been identified, so far, among barley genotypes and, based on the epidemiology of the disease, a quantitative genetic determinacy of RLS has been suggested. The relative contributions of barley and R. collo-cygni genetics to disease infection and epidemiology are practically unknown. Here, we present an integrated genome-wide analysis of host and pathogen transcriptome landscapes identified in a sensitive barley cultivar following infection by an aggressive R. collo-cygni isolate. We compared transcriptional responses in the infected and noninfected leaf samples in order to identify which molecular events are associated with RLS symptom development. We found a large proportion of R. collo-cygni genes to be expressed in planta and that many were also closely associated with the infection stage. The transition from surface to apoplastic colonization was associated with downregulation of cell wall-degrading genes and upregulation of nutrient uptake and resistance to oxidative stresses. Interestingly, the production of secondary metabolites was dynamically regulated within the fungus, indicating that R. collo-cygni produces a diverse panel of toxic compounds according to the infection stage. A defense response against R. collo-cygni was identified in barley at the early, asymptomatic infection and colonization stages. We found activation of ethylene signaling, jasmonic acid signaling, and phenylpropanoid and flavonoid pathways to be highly induced, indicative of a classical response to necrotrophic pathogens. Disease development was found to be associated with gene expression patterns similar to those found at the onset of leaf senescence, when nutrients, possibly, are used by the infecting fungus. These analyses, combining both barley and R. collo-cygni transcript profiles, demonstrate the activation of complex transcriptional programs in both organisms.

Quantifying the Survival of Multiple Salmonella enterica Serovars In Vivo via Massively Parallel Whole-Genome Sequencing To Predict Zoonotic Risk.

Salmonella enterica is an animal and zoonotic pathogen of worldwide importance. Salmonella serovars that differ in their host and tissue tropisms exist. Cattle are an important reservoir of human nontyphoidal salmonellosis, and contaminated bovine peripheral lymph nodes enter the food chain via ground beef. The relative abilities of different serovars to survive within the bovine lymphatic system are poorly understood and constrain the development of control strategies. This problem was addressed by developing a massively parallel whole-genome sequencing method to study mixed-serovar infections in vivoSalmonella serovars differ genetically by naturally occurring single nucleotide polymorphisms (SNPs) in certain genes. It was hypothesized that these SNPs could be used as markers to simultaneously identify serovars in mixed populations and quantify the abundance of each member in a population. The performance of the method was validated in vitro using simulated pools containing up to 11 serovars in various proportions. It was then applied to study serovar survival in vivo in cattle challenged orally with the same 11 serovars. All the serovars successfully colonized the bovine lymphatic system, including the peripheral lymph nodes, and thus pose similar risks of zoonosis. This method enables the fates of multiple genetically unmodified strains to be evaluated simultaneously in a single animal. It could be useful in reducing the number of animals required to study mixed-strain infections and in testing the cross-protective efficacy of vaccines and treatments. It also has the potential to be applied to diverse bacterial species which possess shared but polymorphic alleles.IMPORTANCE While some Salmonella serovars are more frequently isolated from lymph nodes rather than the feces and environment of cattle, the relative abilities of serovars to survive within the lymphatic system of cattle remain ill defined. A sequencing-based method which used available information from sequenced Salmonella genomes to study the dynamics of mixed-serovar infections in vivo was developed. The main advantages of the method include the simultaneous identification and quantification of multiple strains without any genetic modification and minimal animal use. This approach could be used in vaccination trials or in epidemiological surveys where an understanding of the dynamics of closely related strains of a pathogen in mixed populations could inform the prediction of zoonotic risk and the development of intervention strategies.

Maintaining their genetic distance: Little evidence for introgression between widely hybridizing species of Geum with contrasting mating systems.

Within the plant kingdom, many genera contain sister lineages with contrasting outcrossing and inbreeding mating systems that are known to hybridize. The evolutionary fate of these sister lineages is likely to be influenced by the extent to which they exchange genes. We measured gene flow between outcrossing Geum rivale and selfing Geum urbanum, sister species that hybridize in contemporary populations. We generated and used a draft genome of G. urbanum to develop dd-RAD data scorable in both species. Coalescent analysis of RAD data from allopatric populations indicated that the species diverged 2-3 Mya, and that historical gene flow between them was extremely low (1 migrant every 25 generations). Comparison of genetic divergence between species in sympatry and allopatry, together with an analysis of allele frequencies in potential parental and hybrid populations, provided no evidence of contemporary introgression in sympatric populations. Cluster- and species-specific marker analyses revealed that, apart from four early-generation hybrids, individuals in sympatric populations fell into two genetically distinct groups that corresponded exactly to their morphological species classification with maximum individual admixture estimates of only 1-3%. However, we did observe joint segregation of four putatively introgressed SNPs across two scaffolds in the G. urbanum population that was associated with significant morphological variation, interpreted as tentative evidence for rare, recent interspecific gene flow. Overall, our results indicate that despite the presence of hybrids in contemporary populations, genetic exchange between G. rivale and G. urbanum has been extremely limited throughout their evolutionary history.

Development of a Medium Density Combined-Species SNP Array for Pacific and European Oysters (Crassostrea gigas and Ostrea edulis).

SNP arrays are enabling tools for high-resolution studies of the genetic basis of complex traits in farmed and wild animals. Oysters are of critical importance in many regions from both an ecological and economic perspective, and oyster aquaculture forms a key component of global food security. The aim of our study was to design a combined-species, medium density SNP array for Pacific oyster (Crassostrea gigas) and European flat oyster (Ostrea edulis), and to test the performance of this array on farmed and wild populations from multiple locations, with a focus on European populations. SNP discovery was carried out by whole-genome sequencing (WGS) of pooled genomic DNA samples from eight C. gigas populations, and restriction site-associated DNA sequencing (RAD-Seq) of 11 geographically diverse O. edulis populations. Nearly 12 million candidate SNPs were discovered and filtered based on several criteria, including preference for SNPs segregating in multiple populations and SNPs with monomorphic flanking regions. An Affymetrix Axiom Custom Array was created and tested on a diverse set of samples (n = 219) showing ∼27 K high quality SNPs for C. gigas and ∼11 K high quality SNPs for O. edulis segregating in these populations. A high proportion of SNPs were segregating in each of the populations, and the array was used to detect population structure and levels of linkage disequilibrium (LD). Further testing of the array on three C. gigas nuclear families (n = 165) revealed that the array can be used to clearly distinguish between both families based on identity-by-state (IBS) clustering parental assignment software. This medium density, combined-species array will be publicly available through Affymetrix, and will be applied for genome-wide association and evolutionary genetic studies, and for genomic selection in oyster breeding programs.

Implementation and outcomes of a new safeguarding pathway for vulnerable adults presenting to sexual health services.

Adult safeguarding is the process of protecting vulnerable adults from harm or exploitation. In 2014, our sexual health clinic implemented a new adult safeguarding pathway, including an adult safeguarding proforma, an electronic database and a monthly adult safeguarding meeting. We conducted a retrospective case note review of patients entered onto the safeguarding database and found that greater numbers of adults were identified as vulnerable following the introduction of this pathway. Many required referral for onward support, highlighting the importance of robust safeguarding procedures in a sexual health setting.

Design and development of exome capture sequencing for the domestic pig (Sus scrofa).

BACKGROUND: The domestic pig (Sus scrofa) is both an important livestock species and a model for biomedical research. Exome sequencing has accelerated identification of protein-coding variants underlying phenotypic traits in human and mouse. We aimed to develop and validate a similar resource for the pig. RESULTS: We developed probe sets to capture pig exonic sequences based upon the current Ensembl pig gene annotation supplemented with mapped expressed sequence tags (ESTs) and demonstrated proof-of-principle capture and sequencing of the pig exome in 96 pigs, encompassing 24 capture experiments. For most of the samples at least 10x sequence coverage was achieved for more than 90% of the target bases. Bioinformatic analysis of the data revealed over 236,000 high confidence predicted SNPs and over 28,000 predicted indels. CONCLUSIONS: We have achieved coverage statistics similar to those seen with commercially available human and mouse exome kits. Exome capture in pigs provides a tool to identify coding region variation associated with production traits, including loss of function mutations which may explain embryonic and neonatal losses, and to improve genomic assemblies in the vicinity of protein coding genes in the pig.

Genomic variation in macrophage-cultured European porcine reproductive and respiratory syndrome virus Olot/91 revealed using ultra-deep next generation sequencing.

BACKGROUND: Porcine Reproductive and Respiratory Syndrome (PRRS) is a disease of major economic impact worldwide. The etiologic agent of this disease is the PRRS virus (PRRSV). Increasing evidence suggest that microevolution within a coexisting quasispecies population can give rise to high sequence heterogeneity in PRRSV. FINDINGS: We developed a pipeline based on the ultra-deep next generation sequencing approach to first construct the complete genome of a European PRRSV, strain Olot/9, cultured on macrophages and then capture the rare variants representative of the mixed quasispecies population. Olot/91 differs from the reference Lelystad strain by about 5% and a total of 88 variants, with frequencies as low as 1%, were detected in the mixed population. These variants included 16 non-synonymous variants concentrated in the genes encoding structural and nonstructural proteins; including Glycoprotein 2a and 5. CONCLUSION: Using an ultra-deep sequencing methodology, the complete genome of Olot/91 was constructed without any prior knowledge of the sequence. Rare variants that constitute minor fractions of the heterogeneous PRRSV population could successfully be detected to allow further exploration of microevolutionary events.

Assessment of insert sizes and adapter content in fastq data from NexteraXT libraries.

The Illumina NexteraXT transposon protocol is a cost effective way to generate paired end libraries. However, the resulting insert size is highly sensitive to the concentration of DNA used, and the variation of insert sizes is often large. One consequence of this is some fragments may have an insert shorter than the length of a single read, particularly where the library is designed to produce overlapping paired end reads in order to produce longer continuous sequences. Such small insert sizes mean fewer longer reads, and also result in the presence of adapter at the end of the read. Here is presented a protocol to use publicly available tools to identify read pairs with small insert sizes and so likely to contain adapter, to check the sequence of the adapter, and remove adapter sequence from the reads. This protocol does not require a reference genome or prior knowledge of the sequence to be trimmed. Whilst the presence of fragments with small insert sizes may be a particular problem for NexteraXT libraries, the principle can be applied to any Illumina dataset in which the presence of such small inserts is suspected.

Development of a high density 600K SNP genotyping array for chicken.

BACKGROUND: High density (HD) SNP genotyping arrays are an important tool for genetic analyses of animals and plants. Although the chicken is one of the most important farm animals, no HD array is yet available for high resolution genetic analysis of this species. RESULTS: We report here the development of a 600 K Affymetrix® Axiom® HD genotyping array designed using SNPs segregating in a wide variety of chicken populations. In order to generate a large catalogue of segregating SNPs, we re-sequenced 243 chickens from 24 chicken lines derived from diverse sources (experimental, commercial broiler and layer lines) by pooling 10-15 samples within each line. About 139 million (M) putative SNPs were detected by mapping sequence reads to the new reference genome (Gallus_gallus_4.0) of which ~78 M appeared to be segregating in different lines. Using criteria such as high SNP-quality score, acceptable design scores predicting high conversion performance in the final array and uniformity of distribution across the genome, we selected ~1.8 M SNPs for validation through genotyping on an independent set of samples (n = 282). About 64% of the SNPs were polymorphic with high call rates (>98%), good cluster separation and stable Mendelian inheritance. Polymorphic SNPs were further analysed for their population characteristics and genomic effects. SNPs with extreme breach of Hardy-Weinberg equilibrium (P < 0.00001) were excluded from the panel. The final array, designed on the basis of these analyses, consists of 580,954 SNPs and includes 21,534 coding variants. SNPs were selected to achieve an essentially uniform distribution based on genetic map distance for both broiler and layer lines. Due to a lower extent of LD in broilers compared to layers, as reported in previous studies, the ratio of broiler and layer SNPs in the array was kept as 3:2. The final panel was shown to genotype a wide range of samples including broilers and layers with over 100 K to 450 K informative SNPs per line. A principal component analysis was used to demonstrate the ability of the array to detect the expected population structure which is an important pre-investigation step for many genome-wide analyses. CONCLUSIONS: This Affymetrix® Axiom® array is the first SNP genotyping array for chicken that has been made commercially available to the public as a product. This array is expected to find widespread usage both in research and commercial application such as in genomic selection, genome-wide association studies, selection signature analyses, fine mapping of QTLs and detection of copy number variants.

A high resolution genome-wide scan for significant selective sweeps: an application to pooled sequence data in laying chickens.

In most studies aimed at localizing footprints of past selection, outliers at tails of the empirical distribution of a given test statistic are assumed to reflect locus-specific selective forces. Significance cutoffs are subjectively determined, rather than being related to a clear set of hypotheses. Here, we define an empirical p-value for the summary statistic by means of a permutation method that uses the observed SNP structure in the real data. To illustrate the methodology, we applied our approach to a panel of 2.9 million autosomal SNPs identified from re-sequencing a pool of 15 individuals from a brown egg layer line. We scanned the genome for local reductions in heterozygosity, suggestive of selective sweeps. We also employed a modified sliding window approach that accounts for gaps in the sequence and increases scanning resolution by moving the overlapping windows by steps of one SNP only, and suggest to call this a "creeping window" strategy. The approach confirmed selective sweeps in the region of previously described candidate genes, i.e. TSHR, PRL, PRLHR, INSR, LEPR, IGF1, and NRAMP1 when used as positive controls. The genome scan revealed 82 distinct regions with strong evidence of selection (genome-wide p-value<0.001), including genes known to be associated with eggshell structure and immune system such as CALB1 and GAL cluster, respectively. A substantial proportion of signals was found in poor gene content regions including the most extreme signal on chromosome 1. The observation of multiple signals in a highly selected layer line of chicken is consistent with the hypothesis that egg production is a complex trait controlled by many genes.

Computational disease gene identification: a concert of methods prioritizes type 2 diabetes and obesity candidate genes.

Genome-wide experimental methods to identify disease genes, such as linkage analysis and association studies, generate increasingly large candidate gene sets for which comprehensive empirical analysis is impractical. Computational methods employ data from a variety of sources to identify the most likely candidate disease genes from these gene sets. Here, we review seven independent computational disease gene prioritization methods, and then apply them in concert to the analysis of 9556 positional candidate genes for type 2 diabetes (T2D) and the related trait obesity. We generate and analyse a list of nine primary candidate genes for T2D genes and five for obesity. Two genes, LPL and BCKDHA, are common to these two sets. We also present a set of secondary candidates for T2D (94 genes) and for obesity (116 genes) with 58 genes in common to both diseases.

Slug regulates integrin expression and cell proliferation in human epidermal keratinocytes.

The human epidermis is a self-renewing epithelial tissue composed of several layers of keratinocytes. Within the epidermis there exists a complex array of cell adhesion structures, and many of the cellular events within the epidermis (differentiation, proliferation, and migration) require that these adhesion structures be remodeled. The link between cell adhesion, proliferation, and differentiation within the epidermis is well established, and in particular, there is strong evidence to link the process of terminal differentiation to integrin adhesion molecule expression and function. In this paper, we have analyzed the role of a transcriptional repressor called Slug in the regulation of adhesion molecule expression and function in epidermal keratinocytes. We report that activation of Slug, which is expressed predominantly in the basal layer of the epidermis, results in down-regulation of a number of cell adhesion molecules, including E-cadherin, and several integrins, including alpha3, beta1, and beta4. We demonstrate that Slug binds to the alpha3 promoter and that repression of alpha3 transcription by Slug is dependent on an E-box sequence within the promoter. This reduction in integrin expression is reflected in decreased cell adhesion to fibronectin and laminin-5. Despite the reduction in integrin expression and function, we do not observe any increase in differentiation. We do, however, find that activation of Slug results in a significant reduction in keratinocyte proliferation.

POCUS: mining genomic sequence annotation to predict disease genes.

Here we present POCUS (prioritization of candidate genes using statistics), a novel computational approach to prioritize candidate disease genes that is based on over-representation of functional annotation between loci for the same disease. We show that POCUS can provide high (up to 81-fold) enrichment of real disease genes in the candidate-gene shortlists it produces compared with the original large sets of positional candidates. In contrast to existing methods, POCUS can also suggest counterintuitive candidates.

Role and regulation of the fibroblast growth factor axis in human thyroid follicular cells.

Thyroidal levels of fibroblast growth factor-2 (FGF-2) and fibroblast growth factor receptor 1 (FGFR1) are elevated in human thyroid hyperplasia. To understand the significance of this, effects of FGFR1 activation on normal human thyrocyte growth and function in vitro and the regulation of FGF-2 and FGFR1 expression have been examined. FGF-2 stimulated cell growth, as measured by cell counting, and inhibited thyroid function as measured by 125I uptake. Sensitivity to FGF-2 disappeared after 7 days, although FGFR1 expression was maintained. Thyroid-stimulating hormone (TSH, 300 mU/l) increased FGFR1 mRNA expression within 4 h and protein expression by 8 h. Exogenous FGF-2 decreased FGFR1 protein. Endogenous FGF-2 levels were low (approximately 1-2 pg/microg protein), and TSH treatment decreased these by 50%. Protein kinase C (PKC) activation increased FGF-2 mRNA and FGF-2 secretion within 2 h. This effect was enhanced (4.4-fold) when cells were cultured in TSH. We conclude that TSH stimulates FGFR1 but not FGF-2 expression. PKC activation stimulates FGF-2 synthesis and secretion, and TSH synergizes with PKC activators. Increases in FGFR1 or FGF-2 or in both may contribute to goitrogenesis.