Fibroblast Growth Factor 2 Protein Stability Provides Decreased Dependence on Heparin for Induction of FGFR Signaling and Alters ERK Signaling Dynamics.

Fibroblast growth factor 2 (FGF2) plays important roles in tissue development and repair. Using heparan sulfates (HS)/heparin as a cofactor, FGF2 binds to FGF receptor (FGFR) and induces downstream signaling pathways, such as ERK pathway, that regulate cellular behavior. In most cell lines, FGF2 signaling displays biphasic dose-response profile, reaching maximal response to intermediate concentrations, but weak response to high levels of FGF2. Recent reports demonstrated that the biphasic cellular response results from competition between binding of FGF2 to HS and FGFR that impinge upon ERK signaling dynamics. However, the role of HS/heparin in FGF signaling has been controversial. Several studies suggested that heparin is not required for FGF-FGFR complex formation and that the main role of heparin is to protect FGF from degradation. In this study, we investigated the relationship between FGF2 stability, heparin dependence and ERK signaling dynamics using FGF2 variants with increased thermal stability (FGF2-STABs). FGF2-STABs showed higher efficiency in induction of FGFR-mediated proliferation, lower affinity to heparin and were less dependent on heparin than wild-type FGF2 (FGF2-wt) for induction of FGFR-mediated mitogenic response. Interestingly, in primary mammary fibroblasts, FGF2-wt displayed a sigmoidal dose-response profile, while FGF2-STABs showed a biphasic response. Moreover, at low concentrations, FGF2-STABs induced ERK signaling more potently and displayed a faster dynamics of full ERK activation and higher amplitudes of ERK signaling than FGF2-wt. Our results suggest that FGF2 stability and heparin dependence are important factors in FGF-FGFR signaling complex assembly and ERK signaling dynamics.

A locus at 7p14.3 predisposes to refractory celiac disease progression from celiac disease.

BACKGROUND: Approximately 5% of patients with celiac disease (CeD) do not respond to a gluten-free diet and progress to refractory celiac disease (RCD), a severe progression that is characterized by infiltration of intraepithelial T lymphocytes. Patients with RCD type II (RCDII) show clonal expansions of intraepithelial T lymphocytes that result in a poor prognosis and a high mortality rate through development of aggressive enteropathy-associated T-cell lymphoma. It is not known whether genetic variations play a role in severe progression of CeD to RCDII. PATIENTS AND METHODS: We performed the first genome-wide association study to identify the causal genes for RCDII and the molecular pathways perturbed in RCDII. The genome-wide association study was performed in 38 Dutch patients with RCDII, and the 15 independent top-associated single nucleotide polymorphism (SNP) variants (P<5×10) were replicated in 56 independent French and Dutch patients with RCDII. RESULTS: After replication, SNP rs2041570 on chromosome 7 was significantly associated with progression to RCDII (P=2.37×10, odds ratio=2.36) but not with CeD susceptibility. SNP rs2041570 risk allele A was associated with lower levels of FAM188B expression in blood and small intestinal biopsies. Stratification of RCDII biopsies based on rs2041570 genotype showed differential expression of innate immune and antibacterial genes that are expressed in Paneth cells. CONCLUSION: We have identified a novel SNP associated with the severe progression of CeD to RCDII. Our data suggest that genetic susceptibility to CeD might be distinct from the progression to RCDII and suggest a role for Paneth cells in RCDII progression.

Genome-wide Analysis of STAT3-Mediated Transcription during Early Human Th17 Cell Differentiation.

The development of therapeutic strategies to combat immune-associated diseases requires the molecular mechanisms of human Th17 cell differentiation to be fully identified and understood. To investigate transcriptional control of Th17 cell differentiation, we used primary human CD4+ T cells in small interfering RNA (siRNA)-mediated gene silencing and chromatin immunoprecipitation followed by massive parallel sequencing (ChIP-seq) to identify both the early direct and indirect targets of STAT3. The integrated dataset presented in this study confirms that STAT3 is critical for transcriptional regulation of early human Th17 cell differentiation. Additionally, we found that a number of SNPs from loci associated with immune-mediated disorders were located at sites where STAT3 binds to induce Th17 cell specification. Importantly, introduction of such SNPs alters STAT3 binding in DNA affinity precipitation assays. Overall, our study provides important insights for modulating Th17-mediated pathogenic immune responses in humans.

The influence of a short-term gluten-free diet on the human gut microbiome.

BACKGROUND: A gluten-free diet (GFD) is the most commonly adopted special diet worldwide. It is an effective treatment for coeliac disease and is also often followed by individuals to alleviate gastrointestinal complaints. It is known there is an important link between diet and the gut microbiome, but it is largely unknown how a switch to a GFD affects the human gut microbiome. METHODS: We studied changes in the gut microbiomes of 21 healthy volunteers who followed a GFD for four weeks. We collected nine stool samples from each participant: one at baseline, four during the GFD period, and four when they returned to their habitual diet (HD), making a total of 189 samples. We determined microbiome profiles using 16S rRNA sequencing and then processed the samples for taxonomic and imputed functional composition. Additionally, in all 189 samples, six gut health-related biomarkers were measured. RESULTS: Inter-individual variation in the gut microbiota remained stable during this short-term GFD intervention. A number of taxon-specific differences were seen during the GFD: the most striking shift was seen for the family Veillonellaceae (class Clostridia), which was significantly reduced during the intervention (p = 2.81 × 10(-05)). Seven other taxa also showed significant changes; the majority of them are known to play a role in starch metabolism. We saw stronger differences in pathway activities: 21 predicted pathway activity scores showed significant association to the change in diet. We observed strong relations between the predicted activity of pathways and biomarker measurements. CONCLUSIONS: A GFD changes the gut microbiome composition and alters the activity of microbial pathways.

Systematic annotation of celiac disease loci refines pathological pathways and suggests a genetic explanation for increased interferon-gamma levels.

Although genome-wide association studies and fine mapping have identified 39 non-HLA loci associated with celiac disease (CD), it is difficult to pinpoint the functional variants and susceptibility genes in these loci. We applied integrative approaches to annotate and prioritize functional single nucleotide polymorphisms (SNPs), genes and pathways affected in CD. CD-associated SNPs were intersected with regulatory elements categorized by the ENCODE project to prioritize functional variants, while results from cis-expression quantitative trait loci (eQTL) mapping in 1469 blood samples were combined with co-expression analyses to prioritize causative genes. To identify the key cell types involved in CD, we performed pathway analysis on RNA-sequencing data from different immune cell populations and on publicly available expression data on non-immune tissues. We discovered that CD SNPs are significantly enriched in B-cell-specific enhancer regions, suggesting that, besides T-cell processes, B-cell responses play a major role in CD. By combining eQTL and co-expression analyses, we prioritized 43 susceptibility genes in 36 loci. Pathway and tissue-specific expression analyses on these genes suggested enrichment of CD genes in the Th1, Th2 and Th17 pathways, but also predicted a role for four genes in the intestinal barrier function. We also discovered an intricate transcriptional connectivity between CD susceptibility genes and interferon-γ, a key effector in CD, despite the absence of CD-associated SNPs in the IFNG locus. Using systems biology, we prioritized the CD-associated functional SNPs and genes. By highlighting a role for B cells in CD, which classically has been described as a T-cell-driven disease, we offer new insights into the mechanisms and pathways underlying CD.

Expression profiles of long non-coding RNAs located in autoimmune disease-associated regions reveal immune cell-type specificity.

BACKGROUND: Although genome-wide association studies (GWAS) have identified hundreds of variants associated with a risk for autoimmune and immune-related disorders (AID), our understanding of the disease mechanisms is still limited. In particular, more than 90% of the risk variants lie in non-coding regions, and almost 10% of these map to long non-coding RNA transcripts (lncRNAs). lncRNAs are known to show more cell-type specificity than protein-coding genes. METHODS: We aimed to characterize lncRNAs and protein-coding genes located in loci associated with nine AIDs which have been well-defined by Immunochip analysis and by transcriptome analysis across seven populations of peripheral blood leukocytes (granulocytes, monocytes, natural killer (NK) cells, B cells, memory T cells, naive CD4(+) and naive CD8(+) T cells) and four populations of cord blood-derived T-helper cells (precursor, primary, and polarized (Th1, Th2) T-helper cells). RESULTS: We show that lncRNAs mapping to loci shared between AID are significantly enriched in immune cell types compared to lncRNAs from the whole genome (α <0.005). We were not able to prioritize single cell types relevant for specific diseases, but we observed five different cell types enriched (α <0.005) in five AID (NK cells for inflammatory bowel disease, juvenile idiopathic arthritis, primary biliary cirrhosis, and psoriasis; memory T and CD8(+) T cells in juvenile idiopathic arthritis, primary biliary cirrhosis, psoriasis, and rheumatoid arthritis; Th0 and Th2 cells for inflammatory bowel disease, juvenile idiopathic arthritis, primary biliary cirrhosis, psoriasis, and rheumatoid arthritis). Furthermore, we show that co-expression analyses of lncRNAs and protein-coding genes can predict the signaling pathways in which these AID-associated lncRNAs are involved. CONCLUSIONS: The observed enrichment of lncRNA transcripts in AID loci implies lncRNAs play an important role in AID etiology and suggests that lncRNA genes should be studied in more detail to interpret GWAS findings correctly. The co-expression results strongly support a model in which the lncRNA and protein-coding genes function together in the same pathways.

Genetic variation in the non-coding genome: Involvement of micro-RNAs and long non-coding RNAs in disease.

It has been found that the majority of disease-associated genetic variants identified by genome-wide association studies are located outside of protein-coding regions, where they seem to affect regions that control transcription (promoters, enhancers) and non-coding RNAs that also can influence gene expression. In this review, we focus on two classes of non-coding RNAs that are currently a major focus of interest: micro-RNAs and long non-coding RNAs. We describe their biogenesis, suggested mechanism of action, and discuss how these non-coding RNAs might be affected by disease-associated genetic alterations. The discovery of these alterations has already contributed to a better understanding of the etiopathology of human diseases and yielded insight into the function of these non-coding RNAs. We also provide an overview of available databases, bioinformatics tools, and high-throughput techniques that can be used to study the mechanism of action of individual non-coding RNAs. This article is part of a Special Issue entitled: From Genome to Function.

Haplotypes of the IL-1 gene cluster are associated with gastroesophageal reflux disease and Barrett's esophagus.

OBJECTIVES: Gastroesophageal reflux (GERD) is a one of the major public health problem that can lead to reflux esophagitis (RE), Barrett's esophagus (BE), and esophageal adenocarcinoma (EAC). The aim of our study was to determine the impact of IL-1 gene polymorphisms on the development of GERD, RE and BE. METHODS: Three hundred and thirty-three Czech patients with gastroesophageal reflux and 165 healthy controls were included in this case-control study. Four polymorphisms in the genes of the IL-1 cluster [IL-1A(-889C/T), IL-1B(-511C/T), IL-1B(+3953C/T), and IL-1RN(VNTR)] were analyzed. RESULTS: Significant differences were found in IL-1RN 1/2 genotype between patients with GERD/RE and controls and in IL-1B+3953 T allele between patients with BE and healthy subjects. In addition, complex analysis revealed differences in IL-1 haplotype frequencies between the groups. Specifically, the haplotype TCCL was significantly more frequent (p = 0.016) in GERD patients than in controls and the haplotype CCCL more frequent (p = 0.008) in RE patients than in controls. However, in patients with BE, frequency of haplotype TCTL was lower (p = 0.05) and haplotypes CTCL and TCCL were higher (p = 0.03 and p = 0.02) in comparison with the controls. CONCLUSIONS: Our results suggest that IL-1 haplotypes may be associated with susceptibility to GERD, RE and BE.

Human disease-associated genetic variation impacts large intergenic non-coding RNA expression.

Recently it has become clear that only a small percentage (7%) of disease-associated single nucleotide polymorphisms (SNPs) are located in protein-coding regions, while the remaining 93% are located in gene regulatory regions or in intergenic regions. Thus, the understanding of how genetic variations control the expression of non-coding RNAs (in a tissue-dependent manner) has far-reaching implications. We tested the association of SNPs with expression levels (eQTLs) of large intergenic non-coding RNAs (lincRNAs), using genome-wide gene expression and genotype data from five different tissues. We identified 112 cis-regulated lincRNAs, of which 45% could be replicated in an independent dataset. We observed that 75% of the SNPs affecting lincRNA expression (lincRNA cis-eQTLs) were specific to lincRNA alone and did not affect the expression of neighboring protein-coding genes. We show that this specific genotype-lincRNA expression correlation is tissue-dependent and that many of these lincRNA cis-eQTL SNPs are also associated with complex traits and diseases.

Dense genotyping identifies and localizes multiple common and rare variant association signals in celiac disease.

Using variants from the 1000 Genomes Project pilot European CEU dataset and data from additional resequencing studies, we densely genotyped 183 non-HLA risk loci previously associated with immune-mediated diseases in 12,041 individuals with celiac disease (cases) and 12,228 controls. We identified 13 new celiac disease risk loci reaching genome-wide significance, bringing the number of known loci (including the HLA locus) to 40. We found multiple independent association signals at over one-third of these loci, a finding that is attributable to a combination of common, low-frequency and rare genetic variants. Compared to previously available data such as those from HapMap3, our dense genotyping in a large sample collection provided a higher resolution of the pattern of linkage disequilibrium and suggested localization of many signals to finer scale regions. In particular, 29 of the 54 fine-mapped signals seemed to be localized to single genes and, in some instances, to gene regulatory elements. Altogether, we define the complex genetic architecture of the risk regions of and refine the risk signals for celiac disease, providing the next step toward uncovering the causal mechanisms of the disease.

Multiplexed array-based and in-solution genomic enrichment for flexible and cost-effective targeted next-generation sequencing.

The unprecedented increase in the throughput of DNA sequencing driven by next-generation technologies now allows efficient analysis of the complete protein-coding regions of genomes (exomes) for multiple samples in a single sequencing run. However, sample preparation and targeted enrichment of multiple samples has become a rate-limiting and costly step in high-throughput genetic analysis. Here we present an efficient protocol for parallel library preparation and targeted enrichment of pooled multiplexed bar-coded samples. The procedure is compatible with microarray-based and solution-based capture approaches. The high flexibility of this method allows multiplexing of 3-5 samples for whole-exome experiments, 20 samples for targeted footprints of 5 Mb and 96 samples for targeted footprints of 0.4 Mb. From library preparation to post-enrichment amplification, including hybridization time, the protocol takes 5-6 d for array-based enrichment and 3-4 d for solution-based enrichment. Our method provides a cost-effective approach for a broad range of applications, including targeted resequencing of large sample collections (e.g., follow-up genome-wide association studies), and whole-exome or custom mini-genome sequencing projects. This protocol gives details for a single-tube procedure, but scaling to a manual or automated 96-well plate format is possible and discussed.

Relationship between the 17q21 locus and adult asthma in a Czech population.

Several whole-genome association studies have shown a significant link between childhood asthma and the 17q12 chromosome region. We selected tagging single nucleotide polymorphisms (SNPs) in the ORMDL3 gene (17q12) to investigate gene variability in relation to adult allergic asthma and asthma/atopy traits in a Czech Caucasian population of adults. We conducted a case-control association study comprising 668 unrelated subjects (337 asthmatic and 331 control subjects). Four selected SNPs (rs17608925, rs12603332, rs8076131, and rs3169572) were genotyped using the TaqMan SNP Genotyping Assays. The single locus analysis showed only a borderline association between rs3169572 variant and asthma (p = 0.030, p(corr) > 0.05). However, seven different haplotypes were identified; among them, the TTAA haplotype was marginally associated with asthma (p = 0.045, p(corr) > 0.05) and TCAG haplotype was significantly associated with asthma in males (p = 0.009, p(corr) < 0.05, odds ratio = 1.48, 95% confidence interval = 1.10-2.00). In addition, associations between the ORMDL3 genotypes and the total IgE level (p = 0.05, p(corr) > 0.05) and hypersensitivity to the pollen (p = 0.007, p(corr) < 0.05) were established. However, no relationship between ORMDL3 SNPs and the pulmonary functions was found (p > 0.05). These findings suggest that the genetic variability in the 17q21 region may be one of the risk factors also for adult asthma, especially in male individuals.

Haplotype analysis of the interleukin-18 gene in Czech patients with allergic disorders.

The interleukin-18 (IL-18) gene on chromosome 11q22 has been suggested as a susceptibility factor for allergies. To test for a possible role of IL-18 polymorphisms in Czech population, case-control study including 958 subjects (633 allergic patients and 325 healthy controls) was performed. An allele-specific polymerase chain reaction was used to analyze variants at positions -607 C/A (rs1946518) and -137 G/C (rs187238) in the promoter region together with the polymerase chain reaction-restriction fragment length polymorphism method for the detection of polymorphism at position -140 C/G (previously -133 C/G, rs360721) in intron 1 of the IL-18 gene. The -1297 C/T (rs360719) polymorphism was genotyped by real-time-polymerase chain reaction, using a predevelopment TaqMan allele discrimination assay. There were no significant differences in distribution of alleles or genotypes in any of four single nucleotide polymorphisms in the IL-18 gene between controls and patients. However, subsequent analysis revealed a significant difference in haplotype frequencies between the allergic patients and healthy subjects (p < 0.01). Haplotype formed by -1297 C/-607 A/-137 C/-140 C alleles occurred significantly more frequently in patients than controls (0.0433 vs 0.0129; p < 0.0003; p(corr)< 0.01, OR = 3.37; 95% CI = 1.59-7.14). In contrast, there was no relationship among the IL-18 variants and total serum IgE level. Our results indicate that promoter polymorphisms in the IL-18 gene act in interaction and could play a role in allergic disorders.

Association of Toll-like receptor 9 haplotypes with chronic periodontitis in Czech population.

AIM: Toll-like receptors (TLRs) belong to the pattern recognition receptors family of signal molecules that recognize conserved microbial structures. The aim of this study was to analyse polymorphisms in the TLR genes and their association with chronic periodontitis (CP). MATERIAL AND METHODS: Two polymorphisms (2408G/A, i.e. Arg753Gln and -16934A/T) in TLR-2 and three variants (-1486C/T, -1237C/T and+2848A/G) in the TLR-9 genes were studied in 222 patients with CP and 259 unrelated controls. All polymorphisms were detected using the polymerase chain reaction-restriction fragment length polymorphism methods. Subgingival bacterial colonization was investigated by the VariOr Dento test. RESULTS: No significant differences were found in allele and genotype frequencies of all polymorphisms between patients and controls. Nevertheless, complex analysis revealed differences in TLR9 haplotype frequencies between both groups (p=0.001). Specifically, the haplotype T(-1486)/T(-1237)/A(2848) was significantly more frequent (9.6%versus 2.8%, p<0.000001) and the haplotype T(-1486)/T(-1237)/G (2848) of the TLR9 gene was less frequent (35.9%versus 43.3%, p=0.01) in patients than in controls. However, no significant relationships between periodontal pathogens, TLR polymorphisms and CP were found. CONCLUSIONS: In conclusion, although no significant role of the TLR2 gene in periodontitis was found, our results indicate that TLR9 haplotypes may be associated with susceptibility to CP.

Association of single nucleotide polymorphisms in the eosinophil peroxidase gene with allergic rhinitis in the Czech population.

BACKGROUND: Eosinophil peroxidase (EPO) gene codes a cationic protein released from the specific granules of activated eosinophils. Eosinophil granulocytes play a central role in the protection of organisms against parasites. They are also regarded as key effector cells in allergic inflammation. We attempted to determine the polymorphisms in the EPO gene typical for the Czech population and to analyze their associations with allergic rhinitis and its intermediary phenotypes. METHODS: We sequenced all 12 exons of the EPO gene, and selected variants were subsequently analyzed by polymerase chain reaction and restriction fragment length polymorphism methods in a case-control study comprising a total of 613 subjects (319 controls and 294 patients with rhinitis). RESULTS: In total, 5 polymorphisms (-1710T/C and -1710T/CTCC, 2649T/C, 3097A/G and 3979A/G) were found in the EPO gene. Polymorphisms 2649T/C and 3097A/G were in complete linkage disequilibrium (D' = 1 for both groups), and both of them were in a strong disequilibrium with the 3979A/G variant (D' = 0.801 for controls, D' = 0.848 for rhinitics). Consequently, these 3 polymorphisms were studied in association with the allergic phenotype. In a single locus analysis, only 3979A/G single nucleotide polymorphism was marginally significantly associated with rhinitis (p = 0.030, p(corr )> 0.05). This polymorphism also showed a marginal association with total serum IgE levels (log(e) IgE, mean +/- SD: genotypes GG = 2.60 +/- 1.20; GA = 2.47 +/- 1.88; AA = 2.38 +/- 1.49; p < 0.05). Furthermore, significant differences in haplotype frequencies between patients and healthy subjects were observed (p < 0.05). CONCLUSIONS: Our study supports the hypothesis that genetic variability in the EPO gene may contribute to the susceptibility to allergic rhinitis (or related phenotypes) in the Czech population.