Multitissue H3K27ac profiling of GTEx samples links epigenomic variation to disease.

Genetic variants associated with complex traits are primarily noncoding, and their effects on gene-regulatory activity remain largely uncharacterized. To address this, we profile epigenomic variation of histone mark H3K27ac across 387 brain, heart, muscle and lung samples from Genotype-Tissue Expression (GTEx). We annotate 282 k active regulatory elements (AREs) with tissue-specific activity patterns. We identify 2,436 sex-biased AREs and 5,397 genetically influenced AREs associated with 130 k genetic variants (haQTLs) across tissues. We integrate genetic and epigenomic variation to provide mechanistic insights for disease-associated loci from 55 genome-wide association studies (GWAS), by revealing candidate tissues of action, driver SNPs and impacted AREs. Lastly, we build ARE-gene linking scores based on genetics (gLink scores) and demonstrate their unique ability to prioritize SNP-ARE-gene circuits. Overall, our epigenomic datasets, computational integration and mechanistic predictions provide valuable resources and important insights for understanding the molecular basis of human diseases/traits such as schizophrenia.

Publisher Correction: The immune cell landscape in kidneys of patients with lupus nephritis.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

The immune cell landscape in kidneys of patients with lupus nephritis.

Lupus nephritis is a potentially fatal autoimmune disease for which the current treatment is ineffective and often toxic. To develop mechanistic hypotheses of disease, we analyzed kidney samples from patients with lupus nephritis and from healthy control subjects using single-cell RNA sequencing. Our analysis revealed 21 subsets of leukocytes active in disease, including multiple populations of myeloid cells, T cells, natural killer cells and B cells that demonstrated both pro-inflammatory responses and inflammation-resolving responses. We found evidence of local activation of B cells correlated with an age-associated B-cell signature and evidence of progressive stages of monocyte differentiation within the kidney. A clear interferon response was observed in most cells. Two chemokine receptors, CXCR4 and CX3CR1, were broadly expressed, implying a potentially central role in cell trafficking. Gene expression of immune cells in urine and kidney was highly correlated, which would suggest that urine might serve as a surrogate for kidney biopsies.

Mutational processes shape the landscape of TP53 mutations in human cancer.

Unlike most tumor suppressor genes, the most common genetic alterations in tumor protein p53 (TP53) are missense mutations1,2. Mutant p53 protein is often abundantly expressed in cancers and specific allelic variants exhibit dominant-negative or gain-of-function activities in experimental models3-8. To gain a systematic view of p53 function, we interrogated loss-of-function screens conducted in hundreds of human cancer cell lines and performed TP53 saturation mutagenesis screens in an isogenic pair of TP53 wild-type and null cell lines. We found that loss or dominant-negative inhibition of wild-type p53 function reliably enhanced cellular fitness. By integrating these data with the Catalog of Somatic Mutations in Cancer (COSMIC) mutational signatures database9,10, we developed a statistical model that describes the TP53 mutational spectrum as a function of the baseline probability of acquiring each mutation and the fitness advantage conferred by attenuation of p53 activity. Collectively, these observations show that widely-acting and tissue-specific mutational processes combine with phenotypic selection to dictate the frequencies of recurrent TP53 mutations.

A cloning and expression system to probe T-cell receptor specificity and assess functional avidity to neoantigens.

Recent studies have highlighted the promise of targeting tumor neoantigens to generate potent antitumor immune responses and provide strong motivation for improving our understanding of antigen-T-cell receptor (TCR) interactions. Advances in single-cell sequencing technologies have opened the door for detailed investigation of the TCR repertoire, providing paired information from TCRα and TCRß, which together determine specificity. However, a need remains for efficient methods to assess the specificity of discovered TCRs. We developed a streamlined approach for matching TCR sequences with cognate antigen through on-demand cloning and expression of TCRs and screening against candidate antigens. Here, we first demonstrate the system's capacity to identify viral-antigen-specific TCRs and compare the functional avidity of TCRs specific for a given antigen target. We then apply this system to identify neoantigen-specific TCR sequences from patients with melanoma treated with personalized neoantigen vaccines and characterize functional avidity of neoantigen-specific TCRs. Furthermore, we use a neoantigen-prediction pipeline to show that an insertion-deletion mutation in a putative chronic lymphocytic leukemia (CLL) driver gives rise to an immunogenic neoantigen mut-MGA, and use this approach to identify the mut-MGA-specific TCR sequence. This approach provides a means to identify and express TCRs, and then rapidly assess antigen specificity and functional avidity of a reconstructed TCR, which can be applied for monitoring antigen-specific T-cell responses, and potentially for guiding the design of effective T-cell-based immunotherapies.

Comprehensive analysis of cancer-associated somatic mutations in class I HLA genes.

Detection of somatic mutations in human leukocyte antigen (HLA) genes using whole-exome sequencing (WES) is hampered by the high polymorphism of the HLA loci, which prevents alignment of sequencing reads to the human reference genome. We describe a computational pipeline that enables accurate inference of germline alleles of class I HLA-A, B and C genes and subsequent detection of mutations in these genes using the inferred alleles as a reference. Analysis of WES data from 7,930 pairs of tumor and healthy tissue from the same patient revealed 298 nonsilent HLA mutations in tumors from 266 patients. These 298 mutations are enriched for likely functional mutations, including putative loss-of-function events. Recurrence of mutations suggested that these 'hotspot' sites were positively selected. Cancers with recurrent somatic HLA mutations were associated with upregulation of signatures of cytolytic activity characteristic of tumor infiltration by effector lymphocytes, supporting immune evasion by altered HLA function as a contributory mechanism in cancer.

Mutations causing medullary cystic kidney disease type 1 lie in a large VNTR in MUC1 missed by massively parallel sequencing.

Although genetic lesions responsible for some mendelian disorders can be rapidly discovered through massively parallel sequencing of whole genomes or exomes, not all diseases readily yield to such efforts. We describe the illustrative case of the simple mendelian disorder medullary cystic kidney disease type 1 (MCKD1), mapped more than a decade ago to a 2-Mb region on chromosome 1. Ultimately, only by cloning, capillary sequencing and de novo assembly did we find that each of six families with MCKD1 harbors an equivalent but apparently independently arising mutation in sequence markedly under-represented in massively parallel sequencing data: the insertion of a single cytosine in one copy (but a different copy in each family) of the repeat unit comprising the extremely long (∼1.5-5 kb), GC-rich (>80%) coding variable-number tandem repeat (VNTR) sequence in the MUC1 gene encoding mucin 1. These results provide a cautionary tale about the challenges in identifying the genes responsible for mendelian, let alone more complex, disorders through massively parallel sequencing.

Genomic epidemiology of the Escherichia coli O104:H4 outbreaks in Europe, 2011.

The degree to which molecular epidemiology reveals information about the sources and transmission patterns of an outbreak depends on the resolution of the technology used and the samples studied. Isolates of Escherichia coli O104:H4 from the outbreak centered in Germany in May-July 2011, and the much smaller outbreak in southwest France in June 2011, were indistinguishable by standard tests. We report a molecular epidemiological analysis using multiplatform whole-genome sequencing and analysis of multiple isolates from the German and French outbreaks. Isolates from the German outbreak showed remarkably little diversity, with only two single nucleotide polymorphisms (SNPs) found in isolates from four individuals. Surprisingly, we found much greater diversity (19 SNPs) in isolates from seven individuals infected in the French outbreak. The German isolates form a clade within the more diverse French outbreak strains. Moreover, five isolates derived from a single infected individual from the French outbreak had extremely limited diversity. The striking difference in diversity between the German and French outbreak samples is consistent with several hypotheses, including a bottleneck that purged diversity in the German isolates, variation in mutation rates in the two E. coli outbreak populations, or uneven distribution of diversity in the seed populations that led to each outbreak.

A scalable, fully automated process for construction of sequence-ready barcoded libraries for 454.

We present an automated, high throughput library construction process for 454 technology. Sample handling errors and cross-contamination are minimized via end-to-end barcoding of plasticware, along with molecular DNA barcoding of constructs. Automation-friendly magnetic bead-based size selection and cleanup steps have been devised, eliminating major bottlenecks and significant sources of error. Using this methodology, one technician can create 96 sequence-ready 454 libraries in 2 days, a dramatic improvement over the standard method.

Conversion analysis for mutation detection in MLH1 and MSH2 in patients with colorectal cancer.

CONTEXT: The accurate identification and interpretation of germline mutations in mismatch repair genes in colorectal cancer cases is critical for clinical management. Current data suggest that mismatch repair mutations are highly heterogeneous and that many mutations are not detected when conventional DNA sequencing alone is used. OBJECTIVE: To evaluate the potential of conversion analysis compared with DNA sequencing alone to detect heterogeneous germline mutations in MLH1, MSH2, and MSH6 in colorectal cancer patients. DESIGN, SETTING, AND PARTICIPANTS: Multicenter study with patients who participate in the Colon Cancer Family Registry. Mutation analyses were performed in participant samples determined to have a high probability of carrying mismatch repair germline mutations. Samples from a total of 64 hereditary nonpolyposis colorectal cancer cases, 8 hereditary nonpolyposis colorectal cancer-like cases, and 17 cases diagnosed prior to age 50 years were analyzed from June 2002 to June 2003. MAIN OUTCOME MEASURES: Classification of family members as carriers or noncarriers of germline mutations in MLH1, MSH2, or MSH6; mutation data from conversion analysis compared with genomic DNA sequencing. RESULTS: Genomic DNA sequencing identified 28 likely deleterious exon mutations, 4 in-frame deletion mutations, 16 missense changes, and 22 putative splice site mutations. Conversion analysis identified all mutations detected by genomic DNA sequencing--plus an additional exon mutation, 12 large genomic deletions, and 1 exon duplication mutation--yielding an increase of 33% (14/42) in diagnostic yield of deleterious mutations. Conversion analysis also showed that 4 of 16 missense changes resulted in exon skipping in transcripts and that 17 of 22 putative splice site mutations affected splicing or mRNA transcript stability. Conversion analysis provided an increase of 56% (35/63) in the diagnostic yield of genetic testing compared with genomic DNA sequencing alone. CONCLUSIONS: The data confirm the heterogeneity of mismatch repair mutations and reveal that many mutations in colorectal cancer cases would be missed using conventional genomic DNA sequencing alone. Conversion analysis substantially increases the diagnostic yield of genetic testing for mismatch repair mutations in patients diagnosed as having colorectal cancer.