Vincent Butty

Braveheart, a Long Noncoding RNA Required for Cardiovascular Lineage Commitment

Long noncoding RNAs (lncRNAs) are often expressed in a development-specific manner, yet little is known about their roles in lineage commitment. Here, we identified Braveheart (Bvht), a heart-associated lncRNA in mouse. Using multiple embryonic stem cell (ESC) differentiation strategies, we show that Bvht is required for progression of nascent mesoderm toward a cardiac fate. We find that Bvht is necessary for activation of a core cardiovascular gene network and functions upstream of mesoderm posterior 1 (MesP1), a master regulator of a common multipotent cardiovascular progenitor. We also show that Bvht interacts with SUZ12, a component of polycomb-repressive complex 2 (PRC2), during cardiomyocyte differentiation, suggesting that Bvht mediates epigenetic regulation of cardiac commitment. Finally, we demonstrate a role for Bvht in maintaining cardiac fate in neonatal cardiomyocytes. Together, our work provides evidence for a long noncoding RNA with critical roles in the establishment of the cardiovascular lineage during mammalian development.

Digital genotyping and haplotyping with polymerase colonies

Polymerase colony (polony) technology amplifies multiple individual DNA molecules within a thin acrylamide gel attached to a microscope slide. Each DNA molecule included in the reaction produces an immobilized colony of double-stranded DNA. We genotype these polonies by performing single base extensions with dye-labeled nucleotides, and we demonstrate the accurate quantitation of two allelic variants. We also show that polony technology can determine the phase, or haplotype, of two single- nucleotide polymorphisms (SNPs) by coamplifying distally located targets on a single chromosomal fragment. We correctly determine the genotype and phase of three different pairs of SNPs. In one case, the distance between the two SNPs is 45 kb, the largest distance achieved to date without separating the chromosomes by cloning or somatic cell fusion. The results indicate that polony genotyping and haplotyping may play an important role in understanding the structure of genetic variation.

Altered Natural Killer Cells in Type 1 Diabetic Patients

Evidence from animal models suggests that natural killer (NK) cells can be important players in the development of type 1 diabetes, although data in humans are still sparse. We studied the frequency and activation state of blood NK cells at different stages of human type 1 diabetes, and whether genetic or phenotypic NK cell peculiarities could be associated with an early onset of diabetes. The onset period is marked by a slight reduction in blood NK cells, but these are unusually activated in some patients (γ-interferon expression). This activation status does not correlate, however, with a particularly young age at onset. In contrast, NK cells in patients with long-standing type 1 diabetes had a markedly lower expression of p30/p46 NK-activating receptor molecules compared with those of control subjects. A slightly decreased expression of NKG2D in all type 1 diabetic patients relative to control subjects was observed, independent of the duration of disease, parallel to prior observations in the NOD mouse. Finally, type 1 diabetic patients had an increased frequency of KIR gene haplotypes that include the activating KIR2DS3 gene, with a genetic interaction between the KIR and HLA complexes. The reduced activation of NK cells in individuals with long-standing type 1 diabetes would seem to be a consequence rather than a cause, but other peculiarities may relate to type 1 diabetes pathogenesis.

Transcriptional Analysis of Murine Macrophages Infected with Different Toxoplasma Strains Identifies Novel Regulation of Host Signaling Pathways

Most isolates of Toxoplasma from Europe and North America fall into one of three genetically distinct clonal lineages, the type I, II and III lineages. However, in South America these strains are rarely isolated and instead a great variety of other strains are found. T. gondii strains differ widely in a number of phenotypes in mice, such as virulence, persistence, oral infectivity, migratory capacity, induction of cytokine expression and modulation of host gene expression. The outcome of toxoplasmosis in patients is also variable and we hypothesize that, besides host and environmental factors, the genotype of the parasite strain plays a major role. The molecular basis for these differences in pathogenesis, especially in strains other than the clonal lineages, remains largely unexplored. Macrophages play an essential role in the early immune response against T. gondii and are also the cell type preferentially infected in vivo. To determine if non-canonical Toxoplasma strains have unique interactions with the host cell, we infected murine macrophages with 29 different Toxoplasma strains, representing global diversity, and used RNA-sequencing to determine host and parasite transcriptomes. We identified large differences between strains in the expression level of known parasite effectors and large chromosomal structural variation in some strains. We also identified novel strain-specifically regulated host pathways, including the regulation of the type I interferon response by some atypical strains. IFNβ production by infected cells was associated with parasite killing, independent of interferon gamma activation, and dependent on endosomal Toll-like receptors in macrophages and the cytoplasmic receptor retinoic acid-inducible gene 1 (RIG-I) in fibroblasts.

Multiplexed DNA repair assays for multiple lesions and multiple doses via transcription inhibition and transcriptional mutagenesis

Significance DNA, the blueprint of the cell, is constantly damaged by chemicals and radiation. Because DNA damage may cause cell death or mutations that may lead to diseases such as cancer, cells are armed with an arsenal of several distinct mechanisms for repairing the many types of DNA damage that occur. DNA repair capacity (DRC) varies among individuals, and reduced DRC is associated with disease risk; however, the available DRC assays are labor intensive and measure only one pathway at a time. Herein, we present powerful new assays that measure DRC in multiple pathways in a single assay. We use the assays to measure interindividual DRC differences and inhibition of DNA repair, and to uncover unexpected error-prone transcriptional bypass of a thymine dimer. The capacity to repair different types of DNA damage varies among individuals, making them more or less susceptible to the detrimental health consequences of damage exposures. Current methods for measuring DNA repair capacity (DRC) are relatively labor intensive, often indirect, and usually limited to a single repair pathway. Here, we describe a fluorescence-based multiplex flow-cytometric host cell reactivation assay (FM-HCR) that measures the ability of human cells to repair plasmid reporters, each bearing a different type of DNA damage or different doses of the same type of DNA damage. FM-HCR simultaneously measures repair capacity in any four of the following pathways: nucleotide excision repair, mismatch repair, base excision repair, nonhomologous end joining, homologous recombination, and methylguanine methyltransferase. We show that FM-HCR can measure interindividual DRC differences in a panel of 24 cell lines derived from genetically diverse, apparently healthy individuals, and we show that FM-HCR may be used to identify inhibitors or enhancers of DRC. We further develop a next-generation sequencing-based HCR assay (HCR-Seq) that detects rare transcriptional mutagenesis events due to lesion bypass by RNA polymerase, providing an added dimension to DRC measurements. FM-HCR and HCR-Seq provide powerful tools for exploring relationships among global DRC, disease susceptibility, and optimal treatment.

Signatures of strong population differentiation shape extended haplotypes across the human CD28, CTLA4, and ICOS costimulatory genes

The three members of the costimulatory receptor family, CD28, CTLA-4, and ICOS, have complementary effects on T cell activation, and their balance controls the overall outcome of immune and autoimmune responses. They are encoded in a short genomic interval, and overall activity may result from interplay between allelic variants at each locus. With multiethnic DNA panels that represent a wide spectrum of human populations, we demonstrate long-range linkage disequilibrium among the three genes. A large fraction of the variation found in the locus can be explained by the presence of extended haplotypes encompassing variants at CD28, CTLA4, and the ICOS promoter. There are unusual differences in the distribution of some variants and haplotypes between geographic regions. The differences may reflect demographic events and/or the adaptation to diverse environmental and microbial challenges encountered in the course of human migrations and will be important to consider when interpreting association to immune/autoimmune responsiveness.

Methods for time series analysis of RNA-seq data with application to human Th17 cell differentiation

Motivation: Gene expression profiling using RNA-seq is a powerful technique for screening RNA species’ landscapes and their dynamics in an unbiased way. While several advanced methods exist for differential expression analysis of RNA-seq data, proper tools to anal.yze RNA-seq time-course have not been proposed. Results: In this study, we use RNA-seq to measure gene expression during the early human T helper 17 (Th17) cell differentiation and T-cell activation (Th0). To quantify Th17-specific gene expression dynamics, we present a novel statistical methodology, DyNB, for analyzing time-course RNA-seq data. We use non-parametric Gaussian processes to model temporal correlation in gene expression and combine that with negative binomial likelihood for the count data. To account for experiment-specific biases in gene expression dynamics, such as differences in cell differentiation efficiencies, we propose a method to rescale the dynamics between replicated measurements. We develop an MCMC sampling method to make inference of differential expression dynamics between conditions. DyNB identifies several known and novel genes involved in Th17 differentiation. Analysis of differentiation efficiencies revealed consistent patterns in gene expression dynamics between different cultures. We use qRT-PCR to validate differential expression and differentiation efficiencies for selected genes. Comparison of the results with those obtained via traditional timepoint-wise analysis shows that time-course analysis together with time rescaling between cultures identifies differentially expressed genes which would not otherwise be detected. Availability: An implementation of the proposed computational methods will be available at http://research.ics.aalto.fi/csb/software/ Contact: tarmo.aijo@aalto.fi or harri.lahdesmaki@aalto.fi Supplementary information: Supplementary data are available at Bioinformatics online.

Genomic mapping of phosphorothioates reveals partial modification of short consensus sequences

Bacterial phosphorothioate (PT) DNA modifications are incorporated by Dnd proteins A-E and often function with DndF-H as a restriction-modification (R-M) system, as in Escherichia coli B7A. However, bacteria such as Vibrio cyclitrophicus FF75 lack dndF-H, which points to other PT functions. Here we report two novel, orthogonal technologies to map PTs across the genomes of B7A and FF75 with >90% agreement: single molecule, real-time sequencing and deep sequencing of iodine-induced cleavage at PT (ICDS). In B7A, we detect PT on both strands of GpsAAC/GpsTTC motifs, but with only 12% of 40,701 possible sites modified. In contrast, PT in FF75 occurs as a single-strand modification at CpsCA, again with only 14% of 160,541 sites modified. Single-molecule analysis indicates that modification could be partial at any particular genomic site even with active restriction by DndF-H, with direct interaction of modification proteins with GAAC/GTTC sites demonstrated with oligonucleotides. These results point to highly unusual target selection by PT-modification proteins and rule out known R-M mechanisms.

A Broad Analysis of IL1 Polymorphism and Rheumatoid Arthritis

OBJECTIVE It has been suggested that polymorphisms in IL1 are correlated with severe and/or erosive rheumatoid arthritis (RA), but the implicated alleles have differed among studies. The aim of this study was to perform a broad and well-powered search for association between allelic polymorphism in IL1A and IL1B and the susceptibility to or severity of RA. METHODS Key coding and regulatory regions in IL1A and IL1B were sequenced in 24 patients with RA, revealing 4 novel single-nucleotide polymorphisms (SNPs) in IL1B. These and a comprehensive set of 24 SNPs tagging most of the underlying genetic diversity were genotyped in 3 independent RA case-control sample sets and 1 longitudinal RA cohort, totaling 3,561 patients and 3,062 control subjects. RESULTS No fully significant associations were observed. Analysis of the discovery case-control sample sets indicated a potential association of IL1B promoter region SNPs with susceptibility to RA (for RA3/A, odds ratio [OR] 1.27, P = 0.0021) or with the incidence of radiographic erosions (for RA4/C, OR 1.56, P = 0.036), but these findings were not replicated in independent case-control samples. No association with rheumatoid factor, anti-cyclic citrullinated peptide, or the Disease Activity Score in 28 joints was found. None of the associations previously observed in other studies were replicated here. CONCLUSION In spite of a broad and highly powered study, we observed no robust, reproducible association between IL1A/B variants and the susceptibility to or severity of RA in white individuals of European descent. Our results provide evidence that, in the majority of cases, polymorphism in IL1A and IL1B is not a major contributor to genetic susceptibility to RA.

H2A.Z acidic patch couples chromatin dynamics to regulation of gene expression programs during ESC differentiation.

The histone H2A variant H2A.Z is essential for embryonic development and for proper control of developmental gene expression programs in embryonic stem cells (ESCs). Divergent regions of amino acid sequence of H2A.Z likely determine its functional specialization compared to core histone H2A. For example, H2A.Z contains three divergent residues in the essential C-terminal acidic patch that reside on the surface of the histone octamer as an uninterrupted acidic patch domain; however, we know little about how these residues contribute to chromatin structure and function. Here, we show that the divergent amino acids Gly92, Asp97, and Ser98 in the H2A.Z C-terminal acidic patch (H2A.ZAP3) are critical for lineage commitment during ESC differentiation. H2A.Z is enriched at most H3K4me3 promoters in ESCs including poised, bivalent promoters that harbor both activating and repressive marks, H3K4me3 and H3K27me3 respectively. We found that while H2A.ZAP3 interacted with its deposition complex and displayed a highly similar distribution pattern compared to wild-type H2A.Z, its enrichment levels were reduced at target promoters. Further analysis revealed that H2A.ZAP3 was less tightly associated with chromatin, suggesting that the mutant is more dynamic. Notably, bivalent genes in H2A.ZAP3 ESCs displayed significant changes in expression compared to active genes. Moreover, bivalent genes in H2A.ZAP3 ESCs gained H3.3, a variant associated with higher nucleosome turnover, compared to wild-type H2A.Z. We next performed single cell imaging to measure H2A.Z dynamics. We found that H2A.ZAP3 displayed higher mobility in chromatin compared to wild-type H2A.Z by fluorescent recovery after photobleaching (FRAP). Moreover, ESCs treated with the transcriptional inhibitor flavopiridol resulted in a decrease in the H2A.ZAP3 mobile fraction and an increase in its occupancy at target genes indicating that the mutant can be properly incorporated into chromatin. Collectively, our work suggests that the divergent residues in the H2A.Z acidic patch comprise a unique domain that couples control of chromatin dynamics to the regulation of developmental gene expression patterns during lineage commitment.