Scan statistic-based analysis of exome sequencing data identifies FAN1 at 15q13.3 as a susceptibility gene for schizophrenia and autism.

We used a family-based cluster detection approach designed to localize significant rare disease-risk variants clusters within a region of interest to systematically search for schizophrenia (SCZ) susceptibility genes within 49 genomic loci previously implicated by de novo copy number variants. Using two independent whole-exome sequencing family datasets and a follow-up autism spectrum disorder (ASD) case/control whole-exome sequencing dataset, we identified variants in one gene, Fanconi-associated nuclease 1 (FAN1), as being associated with both SCZ and ASD. FAN1 is located in a region on chromosome 15q13.3 implicated by a recurrent copy number variant, which predisposes to an array of psychiatric and neurodevelopmental phenotypes. In both SCZ and ASD datasets, rare nonsynonymous risk variants cluster significantly in affected individuals within a 20-kb window that spans several key functional domains of the gene. Our finding suggests that FAN1 is a key driver in the 15q13.3 locus for the associated psychiatric and neurodevelopmental phenotypes. FAN1 encodes a DNA repair enzyme, thus implicating abnormalities in DNA repair in the susceptibility to SCZ or ASD.

Sequence kernel association tests for the combined effect of rare and common variants.

Recent developments in sequencing technologies have made it possible to uncover both rare and common genetic variants. Genome-wide association studies (GWASs) can test for the effect of common variants, whereas sequence-based association studies can evaluate the cumulative effect of both rare and common variants on disease risk. Many groupwise association tests, including burden tests and variance-component tests, have been proposed for this purpose. Although such tests do not exclude common variants from their evaluation, they focus mostly on testing the effect of rare variants by upweighting rare-variant effects and downweighting common-variant effects and can therefore lose substantial power when both rare and common genetic variants in a region influence trait susceptibility. There is increasing evidence that the allelic spectrum of risk variants at a given locus might include novel, rare, low-frequency, and common genetic variants. Here, we introduce several sequence kernel association tests to evaluate the cumulative effect of rare and common variants. The proposed tests are computationally efficient and are applicable to both binary and continuous traits. Furthermore, they can readily combine GWAS and whole-exome-sequencing data on the same individuals, when available, and are also applicable to deep-resequencing data of GWAS loci. We evaluate these tests on data simulated under comprehensive scenarios and show that compared with the most commonly used tests, including the burden and variance-component tests, they can achieve substantial increases in power. We next show applications to sequencing studies for Crohn disease and autism spectrum disorders. The proposed tests have been incorporated into the software package SKAT.

Scan-statistic approach identifies clusters of rare disease variants in LRP2, a gene linked and associated with autism spectrum disorders, in three datasets.

Cluster-detection approaches, commonly used in epidemiology and astronomy, can be applied in the context of genetic sequence data for the identification of genetic regions significantly enriched with rare disease-risk variants (DRVs). Unlike existing association tests for sequence data, the goal of cluster-detection methods is to localize significant disease mutation clusters within a gene or region of interest. Here, we focus on a chromosome 2q replicated linkage region that is associated with autism spectrum disorder (ASD) and that has been sequenced in three independent datasets. We found that variants in one gene, LRP2, residing on 2q are associated with ASD in two datasets (the combined variable-threshold-test p value is 1.2 × 10(-5)). Using a cluster-detection method, we show that in the discovery and replication datasets, variants associated with ASD cluster preponderantly in 25 kb windows (adjusted p values are p(1) = 0.003 and p(2) = 0.002), and the two windows are highly overlapping. Furthermore, for the third dataset, a 25 kb region similar to those in the other two datasets shows significant evidence of enrichment of rare DRVs. The region implicated by all three studies is involved in ligand binding, suggesting that subtle alterations in either LRP2 expression or LRP2 primary sequence modulate the uptake of LRP2 ligands. BMP4 is a ligand of particular interest given its role in forebrain development, and modest changes in BMP4 binding, which binds to LRP2 near the mutation cluster, might subtly affect development and could lead to autism-associated phenotypes.

Finding disease variants in Mendelian disorders by using sequence data: methods and applications.

Many sequencing studies are now underway to identify the genetic causes for both Mendelian and complex traits. Via exome-sequencing, genes harboring variants implicated in several Mendelian traits have already been identified. The underlying methodology in these studies is a multistep algorithm based on filtering variants identified in a small number of affected individuals and depends on whether they are novel (not yet seen in public resources such as dbSNP), shared among affected individuals, and other external functional information on the variants. Although intuitive, these filter-based methods are nonoptimal and do not provide any measure of statistical uncertainty. We describe here a formal statistical approach that has several distinct advantages: (1) it provides fast computation of approximate p values for individual genes, (2) it adjusts for the background variation in each gene, (3) it allows for incorporation of functional or linkage-based information, and (4) it accommodates designs based on both affected relative pairs and unrelated affected individuals. We show via simulations that the proposed approach can be used in conjunction with the existing filter-based methods to achieve a substantially better ranking of a gene relevant for disease when compared to currently used filter-based approaches, this is especially so in the presence of disease locus heterogeneity. We revisit recent studies on three Mendelian diseases and show that the proposed approach results in the implicated gene being ranked first in all studies, and approximate p values of 10(-6) for the Miller Syndrome gene, 1.0 × 10(-4) for the Freeman-Sheldon Syndrome gene, and 3.5 × 10(-5) for the Kabuki Syndrome gene.

The MIF promoter SNP rs755622 is associated with immune activation in glioblastoma.

Intratumoral heterogeneity is a defining hallmark of glioblastoma, driving drug resistance and ultimately recurrence. Many somatic drivers of microenvironmental change have been shown to affect this heterogeneity and, ultimately, the treatment response. However, little is known about how germline mutations affect the tumoral microenvironment. Here, we find that the single-nucleotide polymorphism (SNP) rs755622 in the promoter of the cytokine macrophage migration inhibitory factor (MIF) is associated with increased leukocyte infiltration in glioblastoma. Furthermore, we identified an association between rs755622 and lactotransferrin expression, which could also be used as a biomarker for immune-infiltrated tumors. These findings demonstrate that a germline SNP in the promoter region of MIF may affect the immune microenvironment and further reveal a link between lactotransferrin and immune activation.

Spatiotemporal evolution of the clear cell renal cell carcinoma microenvironment links intra-tumoral heterogeneity to immune escape.

BACKGROUND: Intratumoral heterogeneity (ITH) is a hallmark of clear cell renal cell carcinoma (ccRCC) that reflects the trajectory of evolution and influences clinical prognosis. Here, we seek to elucidate how ITH and tumor evolution during immune checkpoint inhibitor (ICI) treatment can lead to therapy resistance. METHODS: Here, we completed a single-arm pilot study to examine the safety and feasibility of neoadjuvant nivolumab in patients with localized RCC. Primary endpoints were safety and feasibility of neoadjuvant nivolumab. Then, we spatiotemporally profiled the genomic and immunophenotypic characteristics of 29 ccRCC patients, including pre- and post-therapy samples from 17 ICI-treated patients. Deep multi-regional whole-exome and transcriptome sequencing were performed on 29 patients at different time points before and after ICI therapy. T cell repertoire was also monitored from tissue and peripheral blood collected from a subset of patients to study T cell clonal expansion during ICI therapy. RESULTS: Angiogenesis, lymphocytic infiltration, and myeloid infiltration varied significantly across regions of the same patient, potentially confounding their utility as biomarkers of ICI response. Elevated ITH associated with a constellation of both genomic features (HLA LOH, CDKN2A/B loss) and microenvironmental features, including elevated myeloid expression, reduced peripheral T cell receptor (TCR) diversity, and putative neoantigen depletion. Hypothesizing that ITH may itself play a role in shaping ICI response, we derived a transcriptomic signature associated with neoantigen depletion that strongly associated with response to ICI and targeted therapy treatment in several independent clinical trial cohorts. CONCLUSIONS: These results argue that genetic and immune heterogeneity jointly co-evolve and influence response to ICI in ccRCC. Our findings have implications for future biomarker development for ICI response across ccRCC and other solid tumors and highlight important features of tumor evolution under ICI treatment. TRIAL REGISTRATION: The study was registered on ClinicalTrial.gov (NCT02595918) on November 4, 2015.

Rare structural variants disrupt multiple genes in neurodevelopmental pathways in schizophrenia.

Schizophrenia is a devastating neurodevelopmental disorder whose genetic influences remain elusive. We hypothesize that individually rare structural variants contribute to the illness. Microdeletions and microduplications >100 kilobases were identified by microarray comparative genomic hybridization of genomic DNA from 150 individuals with schizophrenia and 268 ancestry-matched controls. All variants were validated by high-resolution platforms. Novel deletions and duplications of genes were present in 5% of controls versus 15% of cases and 20% of young-onset cases, both highly significant differences. The association was independently replicated in patients with childhood-onset schizophrenia as compared with their parents. Mutations in cases disrupted genes disproportionately from signaling networks controlling neurodevelopment, including neuregulin and glutamate pathways. These results suggest that multiple, individually rare mutations altering genes in neurodevelopmental pathways contribute to schizophrenia.