Immunoseq: the identification of functionally relevant variants through targeted capture and sequencing of active regulatory regions in human immune cells

$\textbf{BACKGROUND}$: The observation that the genetic variants identified in genome-wide association studies (GWAS) frequently lie in non-coding regions of the genome that contain cis-regulatory elements suggests that altered gene expression underlies the development of many complex traits. In order to efficiently make a comprehensive assessment of the impact of non-coding genetic variation in immune related diseases we emulated the whole-exome sequencing paradigm and developed a custom capture panel for the known DNase I hypersensitive site (DHS) in immune cells - "Immunoseq". $\textbf{RESULTS}$: We performed Immunoseq in 30 healthy individuals where we had existing transcriptome data from T cells. We identified a large number of novel non-coding variants in these samples. Relying on allele specific expression measurements, we also showed that our selected capture regions are enriched for functional variants that have an impact on differential allelic gene expression. The results from a replication set with 180 samples confirmed our observations. $\textbf{CONCLUSIONS}$: We show that Immunoseq is a powerful approach to detect novel rare variants in regulatory regions. We also demonstrate that these novel variants have a potential functional role in immune cells.This work was supported by grants from the Canadian Institute of Health Research (CIHR), the UK Medical Research Council (G1100125), the Swedish Research Council (DO283001) and Knut and Alice Wallenberg Foundation (KAW). We also acknowledge the use of subjects from the Cambridge BioResource and the support of the Cambridge NIHR Biomedical Research Centre. AM was supported by the Fond de Recherche Santé Québec Doctoral training award. TP and CL holds a Canada Research Chair

Genetic variations in A20 DUB domain provide a genetic link to citrullination and neutrophil extracellular traps in systemic lupus erythematosus

Objectives: Genetic variations in TNFAIP3 (A20) de-ubiquitinase (DUB) domain increase the risk of systemic lupus erythematosus (SLE) and rheumatoid arthritis. A20 is a negative regulator of NF-κB but the role of its DUB domain and related genetic variants remain unclear. We aimed to study the functional effects of A20 DUB-domain alterations in immune cells and understand its link to SLE pathogenesis. Methods: CRISPR/Cas9 was used to generate human U937 monocytes with A20 DUB-inactivating C103A knock-in (KI) mutation. Whole genome RNA-sequencing was used to identify differentially expressed genes between WT and C103A KI cells. Functional studies were performed in A20 C103A U937 cells and in immune cells from A20 C103A mice and genotyped healthy individuals with A20 DUB polymorphism rs2230926. Neutrophil extracellular trap (NET) formation was addressed ex vivo in neutrophils from A20 C103A mice and SLE-patients with rs2230926. Results: Genetic disruption of A20 DUB domain in human and murine myeloid cells did not give rise to enhanced NF-κB signalling. Instead, cells with C103A mutation or rs2230926 polymorphism presented an upregulated expression of PADI4, an enzyme regulating protein citrullination and NET formation, two key mechanisms in autoimmune pathology. A20 C103A cells exhibited enhanced protein citrullination and extracellular trap formation, which could be suppressed by selective PAD4 inhibition. Moreover, SLE-patients with rs2230926 showed increased NETs and increased frequency of autoantibodies to citrullinated epitopes. Conclusions: We propose that genetic alterations disrupting the A20 DUB domain mediate increased susceptibility to SLE through the upregulation of PADI4 with resultant protein citrullination and extracellular trap formation

Characterization of functional methylomes by next-generation capture sequencing identifies novel disease-associated variants.

Most genome-wide methylation studies (EWAS) of multifactorial disease traits use targeted arrays or enrichment methodologies preferentially covering CpG-dense regions, to characterize sufficiently large samples. To overcome this limitation, we present here a new customizable, cost-effective approach, methylC-capture sequencing (MCC-Seq), for sequencing functional methylomes, while simultaneously providing genetic variation information. To illustrate MCC-Seq, we use whole-genome bisulfite sequencing on adipose tissue (AT) samples and public databases to design AT-specific panels. We establish its efficiency for high-density interrogation of methylome variability by systematic comparisons with other approaches and demonstrate its applicability by identifying novel methylation variation within enhancers strongly correlated to plasma triglyceride and HDL-cholesterol, including at CD36. Our more comprehensive AT panel assesses tissue methylation and genotypes in parallel at ∼4 and ∼3 M sites, respectively. Our study demonstrates that MCC-Seq provides comparable accuracy to alternative approaches but enables more efficient cataloguing of functional and disease-relevant epigenetic and genetic variants for large-scale EWAS.This work was supported by a Canadian Institute of Health Research (CIHR) team grant awarded to E.G., A.T., M.C.V. and M.L. (TEC-128093) and the CIHR funded Epigeneome Mapping Centre at McGill University (EP1-120608) awarded to T.P. and M.L., and the Swedish Research Council, Knut and Alice Wallenberg Foundation and the Torsten Söderberg Foundation awarded to L.R. F.A. holds studentship from The Research Institute of the McGill University Health Center (MUHC). F.G. is a recipient of a research fellowship award from the Heart and Stroke Foundation of Canada. A.T. is the director of a Research Chair in Bariatric and Metabolic Surgery. M.C.V. is the recipient of the Canada Research Chair in Genomics Applied to Nutrition and Health (Tier 1). E.G. and T.P. are recipients of a Canada Research Chair Tier 2 award. The MuTHER Study was funded by a programme grant from the Wellcome Trust (081917/Z/07/Z) and core funding for the Wellcome Trust Centre for Human Genetics (090532). TwinsUK was funded by the Wellcome Trust; European Community's Seventh Framework Programme (FP7/2007-2013). The study also receives support from the National Institute for Health Research (NIHR)-funded BioResource, Clinical Research Facility and Biomedical Research Centre based at Guy's and St Thomas' NHS Foundation Trust in partnership with King's College London. T.D.S. is a holder of an ERC Advanced Principal Investigator award. SNP genotyping was performed by The Wellcome Trust Sanger Institute and National Eye Institute via NIH/CIDR. Finally, we thank the NIH Roadmap Epigenomics Consortium and the Mapping Centers (http://nihroadmap.nih.gov/epigenomics/) for the production of publicly available reference epigenomes. Specifically, we thank the mapping centre at MGH/BROAD for generation of human adipose reference epigenomes used in this study.This is the final version. It was first published by NPG at http://www.nature.com/ncomms/2015/150529/ncomms8211/full/ncomms8211.html#abstrac

Global public perceptions of genomic data sharing: what shapes the willingness to donate DNA and health data?

Analyzing genomic data across populations is central to understanding the role of genetic factors in health and disease. Successful data sharing relies on public support, which requires attention to whether people around the world are willing to donate their data that are then subsequently shared with others for research. However, studies of such public perceptions are geographically limited and do not enable comparison. This paper presents results from a very large public survey on attitudes toward genomic data sharing. Data from 36,268 individuals across 22 countries (gathered in 15 languages) are presented. In general, publics across the world do not appear to be aware of, nor familiar with, the concepts of DNA, genetics, and genomics. Willingness to donate one's DNA and health data for research is relatively low, and trust in the process of data's being shared with multiple users (e.g., doctors, researchers, governments) is also low. Participants were most willing to donate DNA or health information for research when the recipient was specified as a medical doctor and least willing to donate when the recipient was a for-profit researcher. Those who were familiar with genetics and who were trusting of the users asking for data were more likely to be willing to donate. However, less than half of participants trusted more than one potential user of data, although this varied across countries. Genetic information was not uniformly seen as different from other forms of health information, but there was an association between seeing genetic information as special in some way compared to other health data and increased willingness to donate. The global perspective provided by our "Your DNA, Your Say" study is valuable for informing the development of international policy and practice for sharing genomic data. It highlights that the research community not only needs to be worthy of trust by the public, but also urgent steps need to be taken to authentically communicate why genomic research is necessary and how data donation, and subsequent sharing, is integral to this

The DAL10 gene from Norway spruce (Picea abies) belongs to a potentially gymnosperm-specific subclass of MADS-box genes and is specifically active in seed cones and pollen cones.

Transcription factors encoded by different members of the MADS-box gene family have evolved central roles in the regulation of reproductive organ development in the flowering plants, the angiosperms. Development of the stamens and carpels, the pollen- and seed-bearing organs, involves the B- and C-organ-identity MADS-box genes. B- and C-type gene orthologs with activities specifically in developing pollen- and seed-bearing organs are also present in the distantly related gymnosperms: the conifers and the gneto-phytes. We now report on the characterization of DAL10, a novel MADS-box gene from the conifer Norway spruce, which unlike the B- and C-type conifer genes shows no distinct orthology relationship to any angiosperm gene or clade in phylogenetic analyses. Like the B- and C-type genes, it is active specifically in developing pollen cones and seed cones. In situ RNA localization experiments show DAL10 to be expressed in the cone axis, which carry the microsporophylls of the young pollen cone. In contrast, in the seed cone it is expressed both in the cone axis and in the bracts, which subtend the ovuliferous scales. Expression data and the phenotype of transgenic Arabidopsis plants expressing DAL10 suggest that the gene may act upstream to or in concert with the B- and C-type genes in the establishment of reproductive identity of developing cones

Novel risk genes for systemic lupus erythematosus predicted by random forest classification

Genome-wide association studies have identified risk loci for SLE, but a large proportion of the genetic contribution to SLE still remains unexplained. To detect novel risk genes, and to predict an individual's SLE risk we designed a random forest classifier using SNP genotype data generated on the "Immunochip" from 1,160 patients with SLE and 2,711 controls. Using gene importance scores defined by the random forest classifier, we identified 15 potential novel risk genes for SLE. Of them 12 are associated with other autoimmune diseases than SLE, whereas three genes (ZNF804A, CDK1, and MANF) have not previously been associated with autoimmunity. Random forest classification also allowed prediction of patients at risk for lupus nephritis with an area under the curve of 0.94. By allele-specific gene expression analysis we detected cis-regulatory SNPs that affect the expression levels of six of the top 40 genes designed by the random forest analysis, indicating a regulatory role for the identified risk variants. The 40 top genes from the prediction were overrepresented for differential expression in B and T cells according to RNA-sequencing of samples from five healthy donors, with more frequent over-expression in B cells compared to T cells