Genetic Drivers of Epigenetic and Transcriptional Variation in Human Immune Cells.

Characterizing the multifaceted contribution of genetic and epigenetic factors to disease phenotypes is a major challenge in human genetics and medicine. We carried out high-resolution genetic, epigenetic, and transcriptomic profiling in three major human immune cell types (CD14+ monocytes, CD16+ neutrophils, and naive CD4+ T cells) from up to 197 individuals. We assess, quantitatively, the relative contribution of cis-genetic and epigenetic factors to transcription and evaluate their impact as potential sources of confounding in epigenome-wide association studies. Further, we characterize highly coordinated genetic effects on gene expression, methylation, and histone variation through quantitative trait locus (QTL) mapping and allele-specific (AS) analyses. Finally, we demonstrate colocalization of molecular trait QTLs at 345 unique immune disease loci. This expansive, high-resolution atlas of multi-omics changes yields insights into cell-type-specific correlation between diverse genomic inputs, more generalizable correlations between these inputs, and defines molecular events that may underpin complex disease risk.

Targeting inflammatory monocytes by immune-modifying nanoparticles prevents acute kidney allograft rejection.

Inflammatory monocytes are a major component of the cellular infiltrate in acutely rejecting human kidney allografts. Since immune-modifying nanoparticles (IMPs) bind to circulating inflammatory monocytes via the specific scavenger receptor MARCO, causing diversion to the spleen and subsequent apoptosis, we investigated the therapeutic potential of negatively charged, 500-nm diameter polystyrene IMPs to prevent kidney allograft rejection. Kidney transplants were performed from BALB/c (H2d) to C57BL/6 (H2b) mice in two groups: controls (allo) and allo mice infused with IMPs. Groups were studied for 14 (acute rejection) or 100 (chronic rejection) days. Allo mice receiving IMPs exhibited superior survival and markedly less acute rejection, with better kidney function, less tubulitis, and diminished inflammatory cell density, cytokine and cytotoxic molecule expression in the allograft and lower titers of donor-specific IgG2c antibody in serum at day 14, as compared to allo mice. Cells isolated from kidneys from allo mice receiving IMPs showed reduced Ly6Chi monocytes, CD11b+ cells and NKT+ cells compared to allo mice. IMPs predominantly bound CD11b+ cells in the bloodstream and CD11b+ and CD11c-B220+ marginal zone B cells in the spleen. In the spleen, IMPs were found predominantly in red pulp, colocalized with MARCO and expression of cleaved caspase-3. At day 100, allo mice receiving IMPs exhibited reduced macrophage M1 responses but were not protected from chronic rejection. IMPs afforded significant protection from acute rejection, inhibiting both innate and adaptive alloimmunity. Thus, our current experimental findings, coupled with our earlier demonstration of IMP-induced protection in kidney ischemia-reperfusion injury, identify IMPs as a potential induction agent in kidney transplantation.

Whole-genome sequencing identifies EN1 as a determinant of bone density and fracture.

The extent to which low-frequency (minor allele frequency (MAF) between 1-5%) and rare (MAF ≤ 1%) variants contribute to complex traits and disease in the general population is mainly unknown. Bone mineral density (BMD) is highly heritable, a major predictor of osteoporotic fractures, and has been previously associated with common genetic variants, as well as rare, population-specific, coding variants. Here we identify novel non-coding genetic variants with large effects on BMD (ntotal = 53,236) and fracture (ntotal = 508,253) in individuals of European ancestry from the general population. Associations for BMD were derived from whole-genome sequencing (n = 2,882 from UK10K (ref. 10); a population-based genome sequencing consortium), whole-exome sequencing (n = 3,549), deep imputation of genotyped samples using a combined UK10K/1000 Genomes reference panel (n = 26,534), and de novo replication genotyping (n = 20,271). We identified a low-frequency non-coding variant near a novel locus, EN1, with an effect size fourfold larger than the mean of previously reported common variants for lumbar spine BMD (rs11692564(T), MAF = 1.6%, replication effect size = +0.20 s.d., Pmeta = 2 × 10(-14)), which was also associated with a decreased risk of fracture (odds ratio = 0.85; P = 2 × 10(-11); ncases = 98,742 and ncontrols = 409,511). Using an En1(cre/flox) mouse model, we observed that conditional loss of En1 results in low bone mass, probably as a consequence of high bone turnover. We also identified a novel low-frequency non-coding variant with large effects on BMD near WNT16 (rs148771817(T), MAF = 1.2%, replication effect size = +0.41 s.d., Pmeta = 1 × 10(-11)). In general, there was an excess of association signals arising from deleterious coding and conserved non-coding variants. These findings provide evidence that low-frequency non-coding variants have large effects on BMD and fracture, thereby providing rationale for whole-genome sequencing and improved imputation reference panels to study the genetic architecture of complex traits and disease in the general population.

Gut Microbial Metabolites Induce Donor-Specific Tolerance of Kidney Allografts through Induction of T Regulatory Cells by Short-Chain Fatty Acids.

BACKGROUND: Short-chain fatty acids derived from gut microbial fermentation of dietary fiber have been shown to suppress autoimmunity through mechanisms that include enhanced regulation by T regulatory cells (Tregs). METHODS: Using a murine kidney transplantation model, we examined the effects on alloimmunity of a high-fiber diet or supplementation with the short-chain fatty acid acetate. Kidney transplants were performed from BALB/c(H2d) to B6(H2b) mice as allografts in wild-type and recipient mice lacking the G protein-coupled receptor GPR43 (the metabolite-sensing receptor of acetate). Allograft mice received normal chow, a high-fiber diet, or normal chow supplemented with sodium acetate. We assessed rejection at days 14 (acute) and 100 (chronic), and used 16S rRNA sequencing to determine gut microbiota composition pretransplantation and post-transplantation. RESULTS: Wild-type mice fed normal chow exhibited dysbiosis after receiving a kidney allograft but not an isograft, despite the avoidance of antibiotics and immunosuppression for the latter. A high-fiber diet prevented dysbiosis in allograft recipients, who demonstrated prolonged survival and reduced evidence of rejection compared with mice fed normal chow. Allograft mice receiving supplemental sodium acetate exhibited similar protection from rejection, and subsequently demonstrated donor-specific tolerance. Depletion of CD25+ Tregs or absence of the short-chain fatty acid receptor GPR43 abolished this survival advantage. CONCLUSIONS: Manipulation of the microbiome by a high-fiber diet or supplementation with sodium acetate modified alloimmunity in a kidney transplant model, generating tolerance dependent on Tregs and GPR43. Diet-based therapy to induce changes in the gut microbiome can alter systemic alloimmunity in mice, in part through the production of short-chain fatty acids leading to Treg cell development, and merits study as a potential clinical strategy to facilitate transplant acceptance.

Dietary Fiber Protects against Diabetic Nephropathy through Short-Chain Fatty Acid-Mediated Activation of G Protein-Coupled Receptors GPR43 and GPR109A.

BACKGROUND: Studies have reported "dysbiotic" changes to gut microbiota, such as depletion of gut bacteria that produce short-chain fatty acids (SCFAs) through gut fermentation of fiber, in CKD and diabetes. Dietary fiber is associated with decreased inflammation and mortality in CKD, and SCFAs have been proposed to mediate this effect. METHODS: To explore dietary fiber's effect on development of experimental diabetic nephropathy, we used streptozotocin to induce diabetes in wild-type C57BL/6 and knockout mice lacking the genes encoding G protein-coupled receptors GPR43 or GPR109A. Diabetic mice were randomized to high-fiber, normal chow, or zero-fiber diets, or SCFAs in drinking water. We used proton nuclear magnetic resonance spectroscopy for metabolic profiling and 16S ribosomal RNA sequencing to assess the gut microbiome. RESULTS: Diabetic mice fed a high-fiber diet were significantly less likely to develop diabetic nephropathy, exhibiting less albuminuria, glomerular hypertrophy, podocyte injury, and interstitial fibrosis compared with diabetic controls fed normal chow or a zero-fiber diet. Fiber beneficially reshaped gut microbial ecology and improved dysbiosis, promoting expansion of SCFA-producing bacteria of the genera Prevotella and Bifidobacterium, which increased fecal and systemic SCFA concentrations. Fiber reduced expression of genes encoding inflammatory cytokines, chemokines, and fibrosis-promoting proteins in diabetic kidneys. SCFA-treated diabetic mice were protected from nephropathy, but not in the absence of GPR43 or GPR109A. In vitro, SCFAs modulated inflammation in renal tubular cells and podocytes under hyperglycemic conditions. CONCLUSIONS: Dietary fiber protects against diabetic nephropathy through modulation of the gut microbiota, enrichment of SCFA-producing bacteria, and increased SCFA production. GPR43 and GPR109A are critical to SCFA-mediated protection against this condition. Interventions targeting the gut microbiota warrant further investigation as a novel renoprotective therapy in diabetic nephropathy.

Eosinophil microRNAs Play a Regulatory Role in Allergic Diseases Included in the Atopic March.

(1) Background: The atopic march is defined by the increased prevalence of allergic diseases after atopic dermatitis onset. In fact, atopic dermatitis is believed to play an important role in allergen sensitization via the damaged skin barrier, leading to allergic diseases such as allergic asthma and allergic rhinitis. The eosinophil, a pro-inflammatory cell that contributes to epithelial damage, is one of the various cells recruited in the inflammatory reactions characterizing these diseases. Few studies were conducted on the transcriptome of this cell type and even less on their specific microRNA (miRNA) profile, which could modulate pathogenesis of allergic diseases and clinical manifestations post-transcriptionally. Actually, their implication in allergic diseases is not fully understood, but they are believed to play a role in inflammation-related patterns and epithelial cell proliferation. (2) Methods: Next-generation sequencing was performed on RNA samples from eosinophils of individuals with atopic dermatitis, atopy, allergic rhinitis and asthma to obtain differential counts of primary miRNA (pri-miRNA); these were also analyzed for asthma-related phenotypes such as forced expiratory volume in one second (FEV1), immunoglobulin E (IgE) and provocative concentration of methacholine inducing a 20% fall in forced expiratory volume in 1 s (PC20) levels, as well as FEV1 to forced vital capacity (FEV1/FVC) ratio. (3) Results: Eighteen miRNAs from eosinophils were identified to be significantly different between affected individuals and unaffected ones. Based on counts from these miRNAs, individuals were then clustered into groups using Ward's method on Euclidian distances. Groups were found to be explained by asthma diagnosis, familial history of respiratory diseases and allergic rhinitis as well as neutrophil counts. (4) Conclusions: The 18 differential miRNA counts for the studying phenotypes allow a better understanding of the epigenetic mechanisms underlying the development of the allergic diseases included in the atopic march.

Conserved expression of transposon-derived non-coding transcripts in primate stem cells.

BACKGROUND: A significant portion of expressed non-coding RNAs in human cells is derived from transposable elements (TEs). Moreover, it has been shown that various long non-coding RNAs (lncRNAs), which come from the human endogenous retrovirus subfamily H (HERVH), are not only expressed but required for pluripotency in human embryonic stem cells (hESCs). RESULTS: To identify additional TE-derived functional non-coding transcripts, we generated RNA-seq data from induced pluripotent stem cells (iPSCs) of four primate species (human, chimpanzee, gorilla, and rhesus) and searched for transcripts whose expression was conserved. We observed that about 30% of TE instances expressed in human iPSCs had orthologous TE instances that were also expressed in chimpanzee and gorilla. Notably, our analysis revealed a number of repeat families with highly conserved expression profiles including HERVH but also MER53, which is known to be the source of a placental-specific family of microRNAs (miRNAs). We also identified a number of repeat families from all classes of TEs, including MLT1-type and Tigger families, that contributed a significant amount of sequence to primate lncRNAs whose expression was conserved. CONCLUSIONS: Together, these results describe TE families and TE-derived lncRNAs whose conserved expression patterns can be used to identify what are likely functional TE-derived non-coding transcripts in primate iPSCs.

Allelic expression mapping across cellular lineages to establish impact of non-coding SNPs.

Most complex disease-associated genetic variants are located in non-coding regions and are therefore thought to be regulatory in nature. Association mapping of differential allelic expression (AE) is a powerful method to identify SNPs with direct cis-regulatory impact (cis-rSNPs). We used AE mapping to identify cis-rSNPs regulating gene expression in 55 and 63 HapMap lymphoblastoid cell lines from a Caucasian and an African population, respectively, 70 fibroblast cell lines, and 188 purified monocyte samples and found 40-60% of these cis-rSNPs to be shared across cell types. We uncover a new class of cis-rSNPs, which disrupt footprint-derived de novo motifs that are predominantly bound by repressive factors and are implicated in disease susceptibility through overlaps with GWAS SNPs. Finally, we provide the proof-of-principle for a new approach for genome-wide functional validation of transcription factor-SNP interactions. By perturbing NFκB action in lymphoblasts, we identified 489 cis-regulated transcripts with altered AE after NFκB perturbation. Altogether, we perform a comprehensive analysis of cis-variation in four cell populations and provide new tools for the identification of functional variants associated to complex diseases.

RNA sequencing reveals the role of splicing polymorphisms in regulating human gene expression.

Expression levels of many human genes are under the genetic control of expression quantitative trait loci (eQTLs). Despite technological advances, the precise molecular mechanisms underlying most eQTLs remain elusive. Here, we use deep mRNA sequencing of two CEU individuals to investigate those mechanisms, with particular focus on the role of splicing control loci (sQTLs). We identify a large number of genes that are differentially spliced between the two samples and associate many of those differences with nearby single nucleotide polymorphisms (SNPs). Subsequently, we investigate the potential effect of splicing SNPs on eQTL control in general. We find a significant enrichment of alternative splicing (AS) events within a set of highly confident eQTL targets discovered in previous studies, suggesting a role of AS in regulating overall gene expression levels. Next, we demonstrate high correlation between the levels of mature (exonic) and unprocessed (intronic) RNA, implying that ∼75% of eQTL target variance can be explained by control at the level of transcription, but that the remaining 25% may be regulated co- or post-transcriptionally. We focus on eQTL targets with discordant mRNA and pre-mRNA expression patterns and use four examples: USMG5, MMAB, MRPL43, and OAS1, to dissect the exact downstream effects of the associated genetic variants.

Macrophages in renal transplantation: Roles and therapeutic implications.

The presence of macrophages within transplanted renal allografts has been appreciated for some time, whereby macrophages were viewed primarily as participants in the process of cell-mediated allograft rejection. Recent insights into macrophage biology have greatly expanded our conceptual understanding of the multiple roles of macrophages within the allograft. Distinct macrophage subsets are present within the kidney and these sub-serve discrete functions in promoting and attenuating inflammation, immune modulation and tissue repair. Unraveling the complex roles macrophages play in transplantation will allow identification of potential therapeutic targets to prevent and treat allograft rejection and maximize graft longevity.

Global analysis of the impact of environmental perturbation on cis-regulation of gene expression.

Genetic variants altering cis-regulation of normal gene expression (cis-eQTLs) have been extensively mapped in human cells and tissues, but the extent by which controlled, environmental perturbation influences cis-eQTLs is unclear. We carried out large-scale induction experiments using primary human bone cells derived from unrelated donors of Swedish origin treated with 18 different stimuli (7 treatments and 2 controls, each assessed at 2 time points). The treatments with the largest impact on the transcriptome, verified on two independent expression arrays, included BMP-2 (t = 2h), dexamethasone (DEX) (t = 24 h), and PGE2 (t = 24 h). Using these treatments and control, we performed expression profiling for 18,144 RefSeq transcripts on biological replicates of the complete study cohort of 113 individuals (n(total) = 782) and combined it with genome-wide SNP-genotyping data in order to map treatment-specific cis-eQTLs (defined as SNPs located within the gene ± 250 kb). We found that 93% of cis-eQTLs at 1% FDR were observed in at least one additional treatment, and in fact, on average, only 1.4% of the cis-eQTLs were considered as treatment-specific at high confidence. The relative invariability of cis-regulation following perturbation was reiterated independently by genome-wide allelic expression tests where only a small proportion of variance could be attributed to treatment. Treatment-specific cis-regulatory effects were, however, 2- to 6-fold more abundant among differently expressed genes upon treatment. We further followed-up and validated the DEX-specific cis-regulation of the MYO6 and TNC loci and found top cis-regulatory variants located 180 kb and 250 kb upstream of the transcription start sites, respectively. Our results suggest that, as opposed to tissue-specificity of cis-eQTLs, the interactions between cellular environment and cis-variants are relatively rare (∼1.5%), but that detection of such specific interactions can be achieved by a combination of functional genomic approaches as described here.

Interleukin-18 binding protein therapy is protective in adriamycin nephropathy.

Adriamycin nephropathy (AN) is an experimental model of focal segmental glomerulosclerosis (FSGS) in which macrophages are considered to play a pathogenic role. We hypothesize that interleukin-18 (IL-18), largely derived from macrophages, is a key contributor to kidney injury in AN and a potential therapeutic target. In this study, BALB/c mice received adriamycin (9.6 mg/kg) via tail vein injection and subsequently were treated with either neutralizing IL-18 binding protein (IL-18BP; 250 µg) or normal saline (control). At 5 wk, IL-18 was upregulated in AN, and IL-18BP therapy afforded significant protection against the development of AN, resulting in less proteinuria (P < 0.01), kidney dysfunction (P < 0.01), glomerulosclerosis (P < 0.001), and interstitial accumulation of macrophages and T cells (P < 0.001). Gene expression of IL-18 downstream inflammatory molecules, including inducible nitric oxide synthase (P < 0.001), TNF-α (P < 0.001), and IFN-γ (P < 0.01); IL-17 (P < 0.001) and the chemokines CCL2 (P < 0.01) and CCL5 (P < 0.005), was reduced. We demonstrate that IL-18 plays a significant role in the pathogenesis of AN. The protective effect of IL-18BP therapy illustrates the importance of immune mediators in chronic proteinuric kidney disease and highlights the potential of IL-18BP therapy.

Absence of MyD88 signaling induces donor-specific kidney allograft tolerance.

Toll-like receptors (TLRs) play a fundamental role in innate immunity and provide a link between innate and adaptive responses to an allograft; however, whether the development of acute and chronic allograft rejection requires TLR signaling is unknown. Here, we studied TLR signaling in a fully MHC-mismatched, life-sustaining murine model of kidney allograft rejection. Mice deficient in the TLR adaptor protein MyD88 developed donor antigen-specific tolerance, which protected them from both acute and chronic allograft rejection and increased their survival after transplantation compared with wild-type controls. Administration of an anti-CD25 antibody to MyD88-deficient recipients depleted CD4(+)CD25(+)FoxP3(+) cells and broke tolerance. In addition, defective development of Th17 immune responses to alloantigen both in vitro and in vivo occurred, resulting in an increased ratio of Tregs to Th17 effectors. Thus, MyD88 deficiency was associated with an altered balance of Tregs over Th17 cells, promoting tolerance instead of rejection. This study provides evidence that targeting innate immunity may be a clinically relevant strategy to facilitate transplantation tolerance.

Capturing sex-specific and hypofertility-linked effects of assisted reproductive technologies on the cord blood DNA methylome.

BACKGROUND: Children conceived through assisted reproduction are at an increased risk for growth and genomic imprinting disorders, often linked to DNA methylation defects. It has been suggested that assisted reproductive technology (ART) and underlying parental infertility can induce epigenetic instability, specifically interfering with DNA methylation reprogramming events during germ cell and preimplantation development. To date, human studies exploring the association between ART and DNA methylation defects have reported inconsistent or inconclusive results, likely due to population heterogeneity and the use of technologies with limited coverage of the epigenome. In our study, we explored the epigenetic risk of ART by comprehensively profiling the DNA methylome of 73 human cord blood samples of singleton pregnancies (n = 36 control group, n = 37 ART/hypofertile group) from a human prospective longitudinal birth cohort, the 3D (Design, Develop, Discover) Study, using a high-resolution sequencing-based custom capture panel that examines over 2.4 million autosomal CpGs in the genome. RESULTS: We identified evidence of sex-specific effects of ART/hypofertility on cord blood DNA methylation patterns. Our genome-wide analyses identified ~ 46% more CpGs affected by ART/hypofertility in female than in male infant cord blood. We performed a detailed analysis of three imprinted genes which have been associated with altered DNA methylation following ART (KCNQ1OT1, H19/IGF2 and GNAS) and found that female infant cord blood was associated with DNA hypomethylation. When compared to less invasive procedures such as intrauterine insemination, more invasive ARTs (in vitro fertilization, intracytoplasmic sperm injection, embryo culture) resulted in more marked and distinct effects on the cord blood DNA methylome. In the in vitro group, we found a close to fourfold higher proportion of significantly enriched Gene Ontology terms involved in development than in the in vivo group. CONCLUSIONS: Our study highlights the ability of a sensitive, targeted, sequencing-based approach to uncover DNA methylation perturbations in cord blood associated with hypofertility and ART and influenced by offspring sex and ART technique invasiveness.

Transposable elements are associated with the variable response to influenza infection.

Influenza A virus (IAV) infections are frequent every year and result in a range of disease severity. Here, we wanted to explore the potential contribution of transposable elements (TEs) to the variable human immune response. Transcriptome profiling in monocyte-derived macrophages from 39 individuals following IAV infection revealed significant inter-individual variation in viral load post-infection. Using transposase-accessible chromatin using sequencing (ATAC-seq), we identified a set of TE families with either enhanced or reduced accessibility upon infection. Of the enhanced families, 15 showed high variability between individuals and had distinct epigenetic profiles. Motif analysis showed an association with known immune regulators (e.g., BATFs, FOSs/JUNs, IRFs, STATs, NFkBs, NFYs, and RELs) in stably enriched families and with other factors in variable families, including KRAB-ZNFs. We showed that TEs and host factors regulating TEs were predictive of viral load post-infection. Our findings shed light on the role TEs and KRAB-ZNFs may play in inter-individual variation in immunity.

Interaction between genetic and epigenetic variation defines gene expression patterns at the asthma-associated locus 17q12-q21 in lymphoblastoid cell lines.

Phenotypic variation results from variation in gene expression, which is modulated by genetic and/or epigenetic factors. To understand the molecular basis of human disease, interaction between genetic and epigenetic factors needs to be taken into account. The asthma-associated region 17q12-q21 harbors three genes, the zona pellucida binding protein 2 (ZPBP2), gasdermin B (GSDMB) and ORM1-like 3 (ORMDL3), that show allele-specific differences in expression levels in lymphoblastoid cell lines (LCLs) and CD4+ T cells. Here, we report a molecular dissection of allele-specific transcriptional regulation of the genes within the chromosomal region 17q12-q21 combining in vitro transfection, formaldehyde-assisted isolation of regulatory elements, chromatin immunoprecipitation and DNA methylation assays in LCLs. We found that a single nucleotide polymorphism rs4795397 influences the activity of ZPBP2 promoter in vitro in an allele-dependent fashion, and also leads to nucleosome repositioning on the asthma-associated allele. However, variable methylation of exon 1 of ZPBP2 masks the strong genetic effect on ZPBP2 promoter activity in LCLs. In contrast, the ORMDL3 promoter is fully unmethylated, which allows detection of genetic effects on its transcription. We conclude that the cis-regulatory effects on 17q12-q21 gene expression result from interaction between several regulatory polymorphisms and epigenetic factors within the cis-regulatory haplotype region.

Allele-specific chromatin remodeling in the ZPBP2/GSDMB/ORMDL3 locus associated with the risk of asthma and autoimmune disease.

Common SNPs in the chromosome 17q12-q21 region alter the risk for asthma, type 1 diabetes, primary biliary cirrhosis, and Crohn disease. Previous reports by us and others have linked the disease-associated genetic variants with changes in expression of GSDMB and ORMDL3 transcripts in human lymphoblastoid cell lines (LCLs). The variants also alter regulation of other transcripts, and this domain-wide cis-regulatory effect suggests a mechanism involving long-range chromatin interactions. Here, we further dissect the disease-linked haplotype and identify putative causal DNA variants via a combination of genetic and functional analyses. First, high-throughput resequencing of the region and genotyping of potential candidate variants were performed. Next, additional mapping of allelic expression differences in Yoruba HapMap LCLs allowed us to fine-map the basis of the cis-regulatory differences to a handful of candidate functional variants. Functional assays identified allele-specific differences in nucleosome distribution, an allele-specific association with the insulator protein CTCF, as well as a weak promoter activity for rs12936231. Overall, this study shows a common disease allele linked to changes in CTCF binding and nucleosome occupancy leading to altered domain-wide cis-regulation. Finally, a strong association between asthma and cis-regulatory haplotypes was observed in three independent family-based cohorts (p = 1.78 x 10(-8)). This study demonstrates the requirement of multiple parallel allele-specific tools for the investigation of noncoding disease variants and functional fine-mapping of human disease-associated haplotypes.

A probabilistic approach for SNP discovery in high-throughput human resequencing data.

New high-throughput sequencing technologies are generating large amounts of sequence data, allowing the development of targeted large-scale resequencing studies. For these studies, accurate identification of polymorphic sites is crucial. Heterozygous sites are particularly difficult to identify, especially in regions of low coverage. We present a new strategy for identifying heterozygous sites in a single individual by using a machine learning approach that generates a heterozygosity score for each chromosomal position. Our approach also facilitates the identification of regions with unequal representation of two alleles and other poorly sequenced regions. The availability of confidence scores allows for a principled combination of sequencing results from multiple samples. We evaluate our method on a gold standard data genotype set from HapMap. We are able to classify sites in this data set as heterozygous or homozygous with 98.5% accuracy. In de novo data our probabilistic heterozygote detection ("ProbHD") is able to identify 93% of heterozygous sites at a <5% false call rate (FCR) as estimated based on independent genotyping results. In direct comparison of ProbHD with high-coverage 1000 Genomes sequencing available for a subset of our data, we observe >99.9% overall agreement for genotype calls and close to 90% agreement for heterozygote calls. Overall, our data indicate that high-throughput resequencing of human genomic regions requires careful attention to systematic biases in sample preparation as well as sequence contexts, and that their impact can be alleviated by machine learning-based sequence analyses allowing more accurate extraction of true DNA variants.

Population genomics in a disease targeted primary cell model.

The common genetic variants associated with complex traits typically lie in noncoding DNA and may alter gene regulation in a cell type-specific manner. Consequently, the choice of tissue or cell model in the dissection of disease associations is important. We carried out an expression quantitative trait loci (eQTL) study of primary human osteoblasts (HOb) derived from 95 unrelated donors of Swedish origin, each represented by two independently derived primary lines to provide biological replication. We combined our data with publicly available information from a genome-wide association study (GWAS) of bone mineral density (BMD). The top 2000 BMD-associated SNPs (P < approximately 10(-3)) were tested for cis-association of gene expression in HObs and in lymphoblastoid cell lines (LCLs) using publicly available data and showed that HObs have a significantly greater enrichment (threefold) of converging cis-eQTLs as compared to LCLs. The top 10 BMD loci with SNPs showing strong cis-effects on gene expression in HObs (P = 6 x 10(-10) - 7 x 10(-16)) were selected for further validation using a staged design in two cohorts of Caucasian male subjects. All 10 variants were tested in the Swedish MrOS Cohort (n = 3014), providing evidence for two novel BMD loci (SRR and MSH3). These variants were then tested in the Rotterdam Study (n = 2090), yielding converging evidence for BMD association at the 17p13.3 SRR locus (P(combined) = 5.6 x 10(-5)). The cis-regulatory effect was further fine-mapped to the proximal promoter of the SRR gene (rs3744270, r(2) = 0.5, P = 2.6 x 10(-15)). Our results suggest that primary cells relevant to disease phenotypes complement traditional approaches for prioritization and validation of GWAS hits for follow-up studies.

Inherited disorders affecting mitochondrial function are associated with glutathione deficiency and hypocitrullinemia.

Disorders affecting mitochondria, including those that directly affect the respiratory chain function or result from abnormalities in branched amino acid metabolism (organic acidemias), have been shown to be associated with impaired redox balance. Almost all of the evidence underlying this conclusion has been obtained from studies on patient biopsies or animal models. Since the glutathione (iGSH) system provides the main protection against oxidative damage, we hypothesized that untreated oxidative stress in individuals with mitochondrial dysfunction would result in chronic iGSH deficiency. We confirm this hypothesis here in studies using high-dimensional flow cytometry (Hi-D FACS) and biochemical analysis of freshly obtained blood samples from patients with mitochondrial disorders or organic acidemias. T lymphocyte subsets, monocytes and neutrophils from organic acidemia and mitochondrial patients who were not on antioxidant supplements showed low iGSH levels, whereas similar subjects on antioxidant supplements showed normal iGSH. Measures of iROS levels in blood were insufficient to reveal the chronic oxidative stress in untreated patients. Patients with organic acidemias showed elevated plasma protein carbonyls, while plasma samples from all patients tested showed hypocitrullinemia. These findings indicate that measurements of iGSH in leukocytes may be a particularly useful biomarker to detect redox imbalance in mitochondrial disorders and organic acidemias, thus providing a relatively non-invasive means to monitor disease status and response to therapies. Furthermore, studies here suggest that antioxidant therapy may be useful for relieving the chronic oxidative stress that otherwise occurs in patients with mitochondrial dysfunction.