Janna Saarela

Lymphatic endothelial reprogramming of vascular endothelial cells by the Prox-1 homeobox transcription factor

Lymphatic vessels are essential for fluid homeostasis, immune surveillance and fat adsorption, and also serve as a major route for tumor metastasis in many types of cancer. We found that isolated human primary lymphatic and blood vascular endothelial cells (LECs and BECs, respectively) show interesting differences in gene expression relevant for their distinct functions in vivo. Although these phenotypes are stable in vitro and in vivo, overexpression of the homeobox transcription factor Prox‐1 in the BECs was capable of inducing LEC‐specific gene transcription in the BECs, and, surprisingly, Prox‐1 suppressed the expression of ∼40% of the BEC‐specific genes. Prox‐1 did not have global effects on the expression of LEC‐specific genes in other cell types, except that it up‐regulated cyclin E1 and E2 mRNAs and activated the cyclin e promoter in various cell types. These data suggest that Prox‐1 acts as a cell proliferation inducer and a fate determination factor for the LECs. Furthermore, the data provide insights into the phenotypic diversity of endothelial cells and into the possibility of transcriptional reprogramming of differentiated endothelial cells.

Somatic STAT3 mutations in large granular lymphocytic leukemia.

BACKGROUND T-cell large granular lymphocytic leukemia is a rare lymphoproliferative disorder characterized by the expansion of clonal CD3+CD8+ cytotoxic T lymphocytes (CTLs) and often associated with autoimmune disorders and immune-mediated cytopenias. METHODS We used next-generation exome sequencing to identify somatic mutations in CTLs from an index patient with large granular lymphocytic leukemia. Targeted resequencing was performed in a well-characterized cohort of 76 patients with this disorder, characterized by clonal T-cell-receptor rearrangements and increased numbers of large granular lymphocytes. RESULTS Mutations in the signal transducer and activator of transcription 3 gene (STAT3) were found in 31 of 77 patients (40%) with large granular lymphocytic leukemia. Among these 31 patients, recurrent mutational hot spots included Y640F in 13 (17%), D661V in 7 (9%), D661Y in 7 (9%), and N647I in 3 (4%). All mutations were located in exon 21, encoding the Src homology 2 (SH2) domain, which mediates the dimerization and activation of STAT protein. The amino acid changes resulted in a more hydrophobic protein surface and were associated with phosphorylation of STAT3 and its localization in the nucleus. In vitro functional studies showed that the Y640F and D661V mutations increased the transcriptional activity of STAT3. In the affected patients, downstream target genes of the STAT3 pathway (IFNGR2, BCL2L1, and JAK2) were up-regulated. Patients with STAT3 mutations presented more often with neutropenia and rheumatoid arthritis than did patients without these mutations. CONCLUSIONS The SH2 dimerization and activation domain of STAT3 is frequently mutated in patients with large granular lymphocytic leukemia; these findings suggest that aberrant STAT3 signaling underlies the pathogenesis of this disease. (Funded by the Academy of Finland and others.).

Familial combined hyperlipidemia is associated with upstream transcription factor 1 ( USF1 )

Familial combined hyperlipidemia (FCHL), characterized by elevated levels of serum total cholesterol, triglycerides or both, is observed in about 20% of individuals with premature coronary heart disease. We previously identified a locus linked to FCHL on 1q21–q23 in Finnish families with the disease. This region has also been linked to FCHL in families from other populations as well as to type 2 diabetes mellitus. These clinical entities have several overlapping phenotypic features, raising the possibility that the same gene may underlie the obtained linkage results. Here, we show that the human gene encoding thioredoxin interacting protein (TXNIP) on 1q, which underlies combined hyperlipidemia in mice, is not associated with FCHL. We show that FCHL is linked and associated with the gene encoding upstream transcription factor 1 (USF1) in 60 extended families with FCHL, including 721 genotyped individuals (P = 0.00002), especially in males with high triglycerides (P = 0.0000009). Expression profiles in fat biopsy samples from individuals with FCHL seemed to differ depending on their carrier status for the associated USF1 haplotype. USF1 encodes a transcription factor known to regulate several genes of glucose and lipid metabolism.

A predominant role for the HLA class II region in the association of the MHC region with multiple sclerosis

Genetic susceptibility to multiple sclerosis is associated with genes of the major histocompatibility complex (MHC), particularly HLA-DRB1 and HLA-DQB1 (ref. 1). Both locus and allelic heterogeneity have been reported in this genomic region. To clarify whether HLA-DRB1 itself, nearby genes in the region encoding the MHC or combinations of these loci underlie susceptibility to multiple sclerosis, we genotyped 1,185 Canadian and Finnish families with multiple sclerosis (n = 4,203 individuals) with a high-density SNP panel spanning the genes encoding the MHC and flanking genomic regions. Strong associations in Canadian and Finnish samples were observed with blocks in the HLA class II genomic region (P < 4.9 × 10−13 and P < 2.0 × 10−16, respectively), but the strongest association was with HLA-DRB1 (P < 4.4 × 10−17). Conditioning on either HLA-DRB1 or the most significant HLA class II haplotype block found no additional block or SNP association independent of the HLA class II genomic region. This study therefore indicates that MHC-associated susceptibility to multiple sclerosis is determined by HLA class II alleles, their interactions and closely neighboring variants.

Comparison of solution-based exome capture methods for next generation sequencing

BackgroundTechniques enabling targeted re-sequencing of the protein coding sequences of the human genome on next generation sequencing instruments are of great interest. We conducted a systematic comparison of the solution-based exome capture kits provided by Agilent and Roche NimbleGen. A control DNA sample was captured with all four capture methods and prepared for Illumina GAII sequencing. Sequence data from additional samples prepared with the same protocols were also used in the comparison.ResultsWe developed a bioinformatics pipeline for quality control, short read alignment, variant identification and annotation of the sequence data. In our analysis, a larger percentage of the high quality reads from the NimbleGen captures than from the Agilent captures aligned to the capture target regions. High GC content of the target sequence was associated with poor capture success in all exome enrichment methods. Comparison of mean allele balances for heterozygous variants indicated a tendency to have more reference bases than variant bases in the heterozygous variant positions within the target regions in all methods. There was virtually no difference in the genotype concordance compared to genotypes derived from SNP arrays. A minimum of 11× coverage was required to make a heterozygote genotype call with 99% accuracy when compared to common SNPs on genome-wide association arrays.ConclusionsLibraries captured with NimbleGen kits aligned more accurately to the target regions. The updated NimbleGen kit most efficiently covered the exome with a minimum coverage of 20×, yet none of the kits captured all the Consensus Coding Sequence annotated exons.

Gender Differences in Genetic Risk Profiles for Cardiovascular Disease

Background Cardiovascular disease (CVD) incidence, complications and burden differ markedly between women and men. Although there is variation in the distribution of lifestyle factors between the genders, they do not fully explain the differences in CVD incidence and suggest the existence of gender-specific genetic risk factors. We aimed to estimate whether the genetic risk profiles of coronary heart disease (CHD), ischemic stroke and the composite end-point of CVD differ between the genders. Methodology/Principal Findings We studied in two Finnish population cohorts, using the case-cohort design the association between common variation in 46 candidate genes and CHD, ischemic stroke, CVD, and CVD-related quantitative risk factors. We analyzed men and women jointly and also conducted genotype-gender interaction analysis. Several allelic variants conferred disease risk for men and women jointly, including rs1801020 in coagulation factor XII (HR = 1.31 (1.08–1.60) for CVD, uncorrected p = 0.006 multiplicative model). Variant rs11673407 in the fucosyltransferase 3 gene was strongly associated with waist/hip ratio (uncorrected p = 0.00005) in joint analysis. In interaction analysis we found statistical evidence of variant-gender interaction conferring risk of CHD and CVD: rs3742264 in the carboxypeptidase B2 gene, p(interaction) = 0.009 for CHD, and rs2774279 in the upstream stimulatory factor 1 gene, p(interaction) = 0.007 for CHD and CVD, showed strong association in women but not in men, while rs2069840 in interleukin 6 gene, p(interaction) = 0.004 for CVD, showed strong association in men but not in women (uncorrected p-values). Also, two variants in the selenoprotein S gene conferred risk for ischemic stroke in women, p(interaction) = 0.003 and 0.007. Importantly, we identified a larger number of gender-specific effects for women than for men. Conclusions/Significance A false discovery rate analysis suggests that we may expect half of the reported findings for combined gender analysis to be true positives, while at least third of the reported genotype-gender interaction results are true positives. The asymmetry in positive findings between the genders could imply that genetic risk loci for CVD are more readily detectable in women, while for men they are more confounded by environmental/lifestyle risk factors. The possible differences in genetic risk profiles between the genders should be addressed in more detail in genetic studies of CVD, and more focus on female CVD risk is also warranted in genome-wide association studies.

Interferon regulatory factor 5 (IRF5) gene variants are associated with multiple sclerosis in three distinct populations

Background: IRF5 is a transcription factor involved both in the type I interferon and the toll-like receptor signalling pathways. Previously, IRF5 has been found to be associated with systemic lupus erythematosus, rheumatoid arthritis and inflammatory bowel diseases. Here we investigated whether polymorphisms in the IRF5 gene would be associated with yet another disease with features of autoimmunity, multiple sclerosis (MS). Methods: We genotyped nine single nucleotide polymorphisms and one insertion-deletion polymorphism in the IRF5 gene in a collection of 2337 patients with MS and 2813 controls from three populations: two case–control cohorts from Spain and Sweden, and a set of MS trio families from Finland. Results: Two single nucleotide polymorphism (SNPs) (rs4728142, rs3807306), and a 5 bp insertion-deletion polymorphism located in the promoter and first intron of the IRF5 gene, showed association signals with values of p<0.001 when the data from all cohorts were combined. The predisposing alleles were present on the same common haplotype in all populations. Using electrophoretic mobility shift assays we observed allele specific differences in protein binding for the SNP rs4728142 and the 5 bp indel, and by a proximity ligation assay we demonstrated increased binding of the transcription factor SP1 to the risk allele of the 5 bp indel. Conclusion: These findings add IRF5 to the short list of genes shown to be associated with MS in more than one population. Our study adds to the evidence that there might be genes or pathways that are common in multiple autoimmune diseases, and that the type I interferon system is likely to be involved in the development of these diseases.

Possible association of interleukin 1 gene locus polymorphisms with low back pain

&NA; Based on a hypothesis that interleukin 1 (IL‐1) activity is associated with low back pain (LBP), we investigated relationships between previously described functional IL‐1 gene polymorphisms and LBP. The subjects were a subgroup of a Finnish study cohort. The IL‐1&agr;(C889–T), IL‐1&bgr;(C3954–T) and IL‐1 receptor antagonist (IL‐1RN)(G1812–A, G1887–C and T11100–C) polymorphisms were genotyped in 131 middle‐aged men from three occupational groups (machine drivers, carpenters and office workers). A questionnaire inquired about individual and lifestyle characteristics and the occurrence of LBP, the number of days with pain and days with limitation of daily activities because of pain, and pain intensity, during the past 12 months. Lumbar disc degeneration was determined with magnetic resonance imaging. Carriers of the IL‐1RNA1812 allele had an increased risk of LBP (OR 2.5, 95% CI 1.0–6.0) and carriers of this allele in combination with the IL‐1&agr;T889 or IL‐1&bgr;T3954 allele had a higher risk of and more days with LBP than non‐carriers. Pain intensity was associated with the simultaneous carriage of the IL‐1&agr;T889 and IL‐1RNA1812 alleles (OR 3.7, 95% CI 1.2–11.9). Multiple regression analyses allowing for occupation and disc degeneration showed that carriage of the IL‐1RNA1812 allele was associated with the occurrence of pain, the number of days with pain and days with limitations of daily activities. Carriage of the IL‐1&bgr;T3954 allele was associated with the number of days with pain. The results suggest a possible contribution of the IL‐1 gene locus polymorphisms to the pathogenesis of LBP. The possibility of chance findings cannot be excluded due to the small sample size.

Interleukin 1 polymorphisms and intervertebral disc degeneration.

Background: Enzymatic breakdown of the extracellular matrix, and possibly local inflammation, contributes to intervertebral disc degeneration. We investigated whether polymorphisms within the IL-1 gene locus are associated with lumbar disc degeneration and whether the effect of occupational physical load on disc degeneration is modified by the polymorphisms. Methods: Genotypes were determined from 133 middle-aged men who underwent magnetic resonance imaging of the lumbar spine. The participants represented 3 occupations: 40 were machine drivers, 42 carpenters, and 51 office workers. We evaluated decreased signal intensity of the nucleus pulposus, disc bulges, and decreased disc height as signs of degeneration in the L2/L3–L5/S1 discs. Results: The odds ratio for disc bulges was 2.4 (95% confidence interval = 1.2– 4.8) and 1.9 (1.0–3.7), in carriers of the IL-1αT889 or IL-1βT3954 alleles, respectively. The TT genotype of the IL-1α gene carried more than 3-fold risk of disc bulges as compared with the CC genotype. Conclusions: IL-1 gene cluster polymorphisms could affect the risk of disc degeneration. The effect of physical workload seems to be modified by the IL-1 gene polymorphisms.

Meta-analysis of genome-wide association studies confirms a susceptibility locus for knee osteoarthritis on chromosome 7q22

Objectives Osteoarthritis (OA) is the most prevalent form of arthritis and accounts for substantial morbidity and disability, particularly in older people. It is characterised by changes in joint structure, including degeneration of the articular cartilage, and its aetiology is multifactorial with a strong postulated genetic component. Methods A meta-analysis was performed of four genome-wide association (GWA) studies of 2371 cases of knee OA and 35 909 controls in Caucasian populations. Replication of the top hits was attempted with data from 10 additional replication datasets. Results With a cumulative sample size of 6709 cases and 44 439 controls, one genome-wide significant locus was identified on chromosome 7q22 for knee OA (rs4730250, p=9.2×10−9), thereby confirming its role as a susceptibility locus for OA. Conclusion The associated signal is located within a large (500 kb) linkage disequilibrium block that contains six genes: PRKAR2B (protein kinase, cAMP-dependent, regulatory, type II, β), HPB1 (HMG-box transcription factor 1), COG5 (component of oligomeric golgi complex 5), GPR22 (G protein-coupled receptor 22), DUS4L (dihydrouridine synthase 4-like) and BCAP29 (B cell receptor-associated protein 29). Gene expression analyses of the (six) genes in primary cells derived from different joint tissues confirmed expression of all the genes in the joint environment.