Microfluidics-free single-cell genomics with templated emulsification

Current single-cell RNA-sequencing approaches have limitations that stem from the microfluidic devices or fluid handling steps required for sample processing. We develop a method that does not require specialized microfluidic devices, expertise or hardware. Our approach is based on particle-templated emulsification, which allows single-cell encapsulation and barcoding of cDNA in uniform droplet emulsions with only a vortexer. Particle-templated instant partition sequencing (PIP-seq) accommodates a wide range of emulsification formats, including microwell plates and large-volume conical tubes, enabling thousands of samples or millions of cells to be processed in minutes. We demonstrate that PIP-seq produces high-purity transcriptomes in mouse-human mixing studies, is compatible with multiomics measurements and can accurately characterize cell types in human breast tissue compared to a commercial microfluidic platform. Single-cell transcriptional profiling of mixed phenotype acute leukemia using PIP-seq reveals the emergence of heterogeneity within chemotherapy-resistant cell subsets that were hidden by standard immunophenotyping. PIP-seq is a simple, flexible and scalable next-generation workflow that extends single-cell sequencing to new applications.

A dual sgRNA library design to probe genetic modifiers using genome-wide CRISPRi screens

Mapping genetic interactions is essential for determining gene function and defining novel biological pathways. We report a simple to use CRISPR interference (CRISPRi) based platform, compatible with Fluorescence Activated Cell Sorting (FACS)-based reporter screens, to query epistatic relationships at scale. This is enabled by a flexible dual-sgRNA library design that allows for the simultaneous delivery and selection of a fixed sgRNA and a second randomized guide, comprised of a genome-wide library, with a single transduction. We use this approach to identify epistatic relationships for a defined biological pathway, showing both increased sensitivity and specificity than traditional growth screening approaches.

Systematic functional interrogation of SARS-CoV-2 host factors using Perturb-seq

Genomic and proteomic screens have identified numerous host factors of SARS-CoV-2, but efficient delineation of their molecular roles during infection remains a challenge. Here we use Perturb-seq, combining genetic perturbations with a single-cell readout, to investigate how inactivation of host factors changes the course of SARS-CoV-2 infection and the host response in human lung epithelial cells. Our high-dimensional data resolve complex phenotypes such as shifts in the stages of infection and modulations of the interferon response. However, only a small percentage of host factors showed such phenotypes upon perturbation. We further identified the NF-κB inhibitor IκBα (NFKBIA), as well as the translation factors EIF4E2 and EIF4H as strong host dependency factors acting early in infection. Overall, our study provides massively parallel functional characterization of host factors of SARS-CoV-2 and quantitatively defines their roles both in virus-infected and bystander cells.

A dual sgRNA library design to probe genetic modifiers using genome-wide CRISPRi screens

Mapping genetic interactions is essential for determining gene function and defining novel biological pathways. We report a simple to use CRISPR interference (CRISPRi) based platform, compatible with Fluorescence Activated Cell Sorting (FACS)-based reporter screens, to query epistatic relationships at scale. This is enabled by a flexible dual-sgRNA library design that allows for the simultaneous delivery and selection of a fixed sgRNA and a second randomized guide, comprised of a genome-wide library, with a single transduction. We use this approach to identify epistatic relationships for a defined biological pathway, showing both increased sensitivity and specificity than traditional growth screening approaches.

MTCH2 is a mitochondrial outer membrane protein insertase

In the mitochondrial outer membrane, α-helical transmembrane proteins play critical roles in cytoplasmic-mitochondrial communication. Using genome-wide CRISPR screens, we identified MTCH2, and its paralog MTCH1, and showed that it is required for insertion of biophysically diverse tail-anchored (TA), signal-anchored, and multipass proteins, but not outer membrane β-barrel proteins. Purified MTCH2 was sufficient to mediate insertion into reconstituted proteoliposomes. Functional and mutational studies suggested that MTCH2 has evolved from a solute carrier transporter. MTCH2 uses membrane-embedded hydrophilic residues to function as a gatekeeper for the outer membrane, controlling mislocalization of TAs into the endoplasmic reticulum and modulating the sensitivity of leukemia cells to apoptosis. Our identification of MTCH2 as an insertase provided a mechanistic explanation for the diverse phenotypes and disease states associated with MTCH2 dysfunction.

We showed that MTCH2 was both necessary and sufficient for insertion of diverse α-helical proteins into the mitochondrial outer membrane, and was the defining member of a family of insertases that have co-opted the SLC25 transporter fold.

Lineage tracing reveals the phylodynamics, plasticity, and paths of tumor evolution

Tumor evolution is driven by the progressive acquisition of genetic and epigenetic alterations that enable uncontrolled growth and expansion to neighboring and distal tissues. The study of phylogenetic relationships between cancer cells provides key insights into these processes. Here, we introduced an evolving lineage-tracing system with a single-cell RNA-seq readout into a mouse model of Kras;Trp53(KP)-driven lung adenocarcinoma and tracked tumor evolution from single-transformed cells to metastatic tumors at unprecedented resolution. We found that the loss of the initial, stable alveolar-type2-like state was accompanied by a transient increase in plasticity. This was followed by the adoption of distinct transcriptional programs that enable rapid expansion and, ultimately, clonal sweep of stable subclones capable of metastasizing. Finally, tumors develop through stereotypical evolutionary trajectories, and perturbing additional tumor suppressors accelerates progression by creating novel trajectories. Our study elucidates the hierarchical nature of tumor evolution and, more broadly, enables in-depth studies of tumor progression.

Mapping transcriptomic vector fields of single cells

Single-cell (sc)-RNA-seq, together with RNA-velocity and metabolic labeling, reveals cellular states and transitions at unprecedented resolution. Fully exploiting these data, however, requires kinetic models capable of unveiling governing regulatory functions. Here, we introduce an analytical framework dynamo, that infers absolute RNA velocity, reconstructs continuous vector-field functions that predict cell fates, employs differential geometry to extract underlying regulations, and ultimately predicts optimal reprogramming paths and perturbation outcomes. We highlight dynamo’s power to overcome fundamental limitations of conventional splicing-based RNA velocity analyses to enable accurate velocity estimations on a metabolically-labeled human hematopoiesis scRNA-seq dataset. Furthermore, differential geometry analyses reveal mechanisms driving early megakaryocyte appearance and elucidate asymmetrical regulation within the PU.1–GATA1 circuit. Leveraging the Least-Action-Path method, dynamo accurately predicts drivers of numerous hematopoietic transitions. Finally, in silico perturbations predict cell-fate diversions induced by gene perturbations. Dynamo thus represents an important step in advancing quantitative and predictive theories of cell-state transitions.

Maximizing CRISPRi efficacy and accessibility with dual-sgRNA libraries and optimal effectors

CRISPR interference (CRISPRi) enables programmable, reversible, and titratable repression of gene expression (knockdown) in mammalian cells. Initial CRISPRi-mediated genetic screens have showcased the potential to address basic questions in cell biology, genetics, and biotechnology, but wider deployment of CRISPRi screening has been constrained by the large size of single guide RNA (sgRNA) libraries and challenges in generating cell models with consistent CRISPRi-mediated knockdown. Here, we present next-generation CRISPRi sgRNA libraries and effector expression constructs that enable strong and consistent knockdown across mammalian cell models. First, we combine empirical sgRNA selection with a dual-sgRNA library design to generate an ultra-compact (1-3 elements per gene), highly active CRISPRi sgRNA library. Next, we compare CRISPRi effectors to show that the recently published Zim3-dCas9 provides an excellent balance between strong on-target knockdown and minimal non-specific effects on cell growth or the transcriptome. Finally, we engineer a suite of cell lines with stable expression of Zim3-dCas9 and robust on-target knockdown. Our results and publicly available reagents establish best practices for CRISPRi genetic screening.

Efficient C•G-to-G•C base editors developed using CRISPRi screens, target-library analysis, and machine learning

Programmable C•G-to-G•C base editors (CGBEs) have broad scientific and therapeutic potential, but their editing outcomes have proved difficult to predict and their editing efficiency and product purity are often low. We describe a suite of engineered CGBEs paired with machine learning models to enable efficient, high-purity C•G-to-G•C base editing. We performed a CRISPR interference (CRISPRi) screen targeting DNA repair genes to identify factors that affect C•G-to-G•C editing outcomes and used these insights to develop CGBEs with diverse editing profiles. We characterized ten promising CGBEs on a library of 10,638 genomically integrated target sites in mammalian cells and trained machine learning models that accurately predict the purity and yield of editing outcomes (R = 0.90) using these data. These CGBEs enable correction to the wild-type coding sequence of 546 disease-related transversion single-nucleotide variants (SNVs) with >90% precision (mean 96%) and up to 70% efficiency (mean 14%). Computational prediction of optimal CGBE-single-guide RNA pairs enables high-purity transversion base editing at over fourfold more target sites than achieved using any single CGBE variant.

Genome-wide programmable transcriptional memory by CRISPR-based epigenome editing

A general approach for heritably altering gene expression has the potential to enable many discovery and therapeutic efforts. Here, we present CRISPRoff-a programmable epigenetic memory writer consisting of a single dead Cas9 fusion protein that establishes DNA methylation and repressive histone modifications. Transient CRISPRoff expression initiates highly specific DNA methylation and gene repression that is maintained through cell division and differentiation of stem cells to neurons. Pairing CRISPRoff with genome-wide screens and analysis of chromatin marks establishes rules for heritable gene silencing. We identify single guide RNAs (sgRNAs) capable of silencing the large majority of genes including those lacking canonical CpG islands (CGIs) and reveal a wide targeting window extending beyond annotated CGIs. The broad ability of CRISPRoff to initiate heritable gene silencing even outside of CGIs expands the canonical model of methylation-based silencing and enables diverse applications including genome-wide screens, multiplexed cell engineering, enhancer silencing, and mechanistic exploration of epigenetic inheritance.

GIGYF2 and 4EHP Inhibit Translation Initiation of Defective Messenger RNAs to Assist Ribosome-Associated Quality Control

Ribosome-associated quality control (RQC) pathways protect cells from toxicity caused by incomplete protein products resulting from translation of damaged or problematic mRNAs. Extensive work in yeast has identified highly conserved mechanisms that lead to degradation of faulty mRNA and partially synthesized polypeptides. Here we used CRISPR-Cas9-based screening to search for additional RQC strategies in mammals. We found that failed translation leads to specific inhibition of translation initiation on that message. This negative feedback loop is mediated by two translation inhibitors, GIGYF2 and 4EHP. Model substrates and growth-based assays established that inhibition of additional rounds of translation acts in concert with known RQC pathways to prevent buildup of toxic proteins. Inability to block translation of faulty mRNAs and subsequent accumulation of partially synthesized polypeptides could explain the neurodevelopmental and neuropsychiatric disorders observed in mice and humans with compromised GIGYF2 function.

Combinatorial single-cell CRISPR screens by direct guide RNA capture and targeted sequencing

Single-cell CRISPR screens enable the exploration of mammalian gene function and genetic regulatory networks. However, use of this technology has been limited by reliance on indirect indexing of single-guide RNAs (sgRNAs). Here we present direct-capture Perturb-seq, a versatile screening approach in which expressed sgRNAs are sequenced alongside single-cell transcriptomes. Direct-capture Perturb-seq enables detection of multiple distinct sgRNA sequences from individual cells and thus allows pooled single-cell CRISPR screens to be easily paired with combinatorial perturbation libraries that contain dual-guide expression vectors. We demonstrate the utility of this approach for high-throughput investigations of genetic interactions and, leveraging this ability, dissect epistatic interactions between cholesterol biogenesis and DNA repair. Using direct capture Perturb-seq, we also show that targeting individual genes with multiple sgRNAs per cell improves efficacy of CRISPR interference and activation, facilitating the use of compact, highly active CRISPR libraries for single-cell screens. Last, we show that hybridization-based target enrichment permits sensitive, specific sequencing of informative transcripts from single-cell RNA-seq experiments.

Exploring genetic interaction manifolds constructed from rich single-cell phenotypes

How cellular and organismal complexity emerges from combinatorial expression of genes is a central question in biology. High-content phenotyping approaches such as Perturb-seq (single-cell RNA-sequencing pooled CRISPR screens) present an opportunity for exploring such genetic interactions (GIs) at scale. Here, we present an analytical framework for interpreting high-dimensional landscapes of cell states (manifolds) constructed from transcriptional phenotypes. We applied this approach to Perturb-seq profiling of strong GIs mined from a growth-based, gain-of-function GI map. Exploration of this manifold enabled ordering of regulatory pathways, principled classification of GIs (e.g., identifying suppressors), and mechanistic elucidation of synergistic interactions, including an unexpected synergy between CBL and CNN1 driving erythroid differentiation. Finally, we applied recommender system machine learning to predict interactions, facilitating exploration of vastly larger GI manifolds.

A human microglia-like cellular model for assessing the effects of neurodegenerative disease gene variants

Microglia are emerging as a key cell type in neurodegenerative diseases, yet human microglia are challenging to study in vitro. We developed an in vitro cell model system composed of human monocyte-derived microglia-like (MDMi) cells that recapitulated key aspects of microglia phenotype and function. We then used this model system to perform an expression quantitative trait locus (eQTL) study examining 94 genes from loci associated with Alzheimer's disease, Parkinson's disease, and multiple sclerosis. We found six loci (CD33, PILRB, NUP160, LRRK2, RGS1, and METTL21B) in which the risk haplotype drives the association with both disease susceptibility and altered expression of a nearby gene (cis-eQTL). In the PILRB and LRRK2 loci, the cis-eQTL was found in the MDMi cells but not in human peripheral blood monocytes, suggesting that differentiation of monocytes into microglia-like cells led to the acquisition of a cellular state that could reveal the functional consequences of certain genetic variants. We further validated the effect of risk haplotypes at the protein level for PILRB and CD33, and we confirmed that the CD33 risk haplotype altered phagocytosis by the MDMi cells. We propose that increased LRRK2 gene expression by MDMi cells could be a functional outcome of rs76904798, a single-nucleotide polymorphism in the LRKK2 locus that is associated with Parkinson's disease.

Genetic architecture of age-related cognitive decline in African Americans

OBJECTIVE

To identify genetic risk factors associated with susceptibility to age-related cognitive decline in African Americans (AAs).

METHODS

We performed a genome-wide association study (GWAS) and an admixture-mapping scan in 3,964 older AAs from 5 longitudinal cohorts; for each participant, we calculated a slope of an individual's global cognitive change from neuropsychological evaluations. We also performed a pathway-based analysis of the age-related cognitive decline GWAS.

RESULTS

We found no evidence to support the existence of a genomic region which has a strongly different contribution to age-related cognitive decline in African and European genomes. Known Alzheimer disease (AD) susceptibility variants in the ABCA7 and MS4A loci do influence this trait in AAs. Of interest, our pathway-based analyses returned statistically significant results highlighting a shared risk from lipid/metabolism and protein tyrosine signaling pathways between cognitive decline and AD, but the role of inflammatory pathways is polarized, being limited to AD susceptibility.

CONCLUSIONS

The genetic architecture of aging-related cognitive in AA individuals is largely similar to that of individuals of European descent. In both populations, we note a surprising lack of enrichment for immune pathways in the genetic risk for cognitive decline, despite strong enrichment of these pathways among genetic risk factors for AD.

Trans-pQTL study identifies immune crosstalk between Parkinson and Alzheimer loci

OBJECTIVE

Given evidence from genetic studies, we hypothesized that there may be a shared component to the role of myeloid function in Parkinson and Alzheimer disease (PD and AD) and assessed whether PD susceptibility variants influenced protein expression of well-established AD-associated myeloid genes in human monocytes.

METHODS

We repurposed data in which AD-related myeloid proteins CD33, TREM1, TREM2, TREML2, TYROBP, and PTK2B were measured by flow cytometry in monocytes from 176 participants of the PhenoGenetic Project (PGP) and Harvard Aging Brain Study. Linear regression was used to identify associations between 24 PD risk variants and protein expression. The 2 cohorts were meta-analyzed in a discovery analysis, and the 4 most strongly suggestive results were validated in an independent cohort of 50 PGP participants.

RESULTS

We discovered and validated an association between the PD risk allele rs12456492(G) in the RIT2 locus and increased CD33 expression (p joint = 3.50 × 10(-5)) and found strongly suggestive evidence that rs11060180(A) in the CCDC62/HIP1R locus decreased PTK2B expression (p joint = 1.12 × 10(-4)). Furthermore, in older individuals, increased CD33 expression on peripheral monocytes was associated with a greater burden of parkinsonism (p = 0.047), particularly bradykinesia (p = 6.64 × 10(-3)).

CONCLUSIONS

We find that the rs12456492 PD risk variant affects expression of AD-associated protein CD33 in peripheral monocytes, which suggests that genetic factors for these 2 diseases may converge to influence overlapping innate immune-mediated mechanisms that contribute to neurodegeneration. Furthermore, the effect of the rs12456492(G) PD risk allele on increased CD33 suggests that the inhibition of certain myeloid functions may contribute to PD susceptibility, as is the case for AD.

Nucleosomes impede Cas9 access to DNA in vivo and in vitro

The prokaryotic CRISPR (clustered regularly interspaced palindromic repeats)-associated protein, Cas9, has been widely adopted as a tool for editing, imaging, and regulating eukaryotic genomes. However, our understanding of how to select single-guide RNAs (sgRNAs) that mediate efficient Cas9 activity is incomplete, as we lack insight into how chromatin impacts Cas9 targeting. To address this gap, we analyzed large-scale genetic screens performed in human cell lines using either nuclease-active or nuclease-dead Cas9 (dCas9). We observed that highly active sgRNAs for Cas9 and dCas9 were found almost exclusively in regions of low nucleosome occupancy. In vitro experiments demonstrated that nucleosomes in fact directly impede Cas9 binding and cleavage, while chromatin remodeling can restore Cas9 access. Our results reveal a critical role of eukaryotic chromatin in dictating the targeting specificity of this transplanted bacterial enzyme, and provide rules for selecting Cas9 target sites distinct from and complementary to those based on sequence properties.

DOI:http://dx.doi.org/10.7554/eLife.12677.001

Many bacteria have a type of immune system known as CRISPR that can target and cut foreign DNA to protect it against viruses. Recently, the CRISPR system was adapted to allow scientists to easily manipulate the genome of humans and many other organisms. However, unlike the loosely organized DNA found in bacteria, the DNA that makes up the human genome is tightly packed and wrapped around complexes of proteins to form structures called nucleosomes. It was not clear whether the CRISPR system was able to effectively target the stretches of DNA in a nucleosome.

In 2013, researchers developed a modified version of CRISPR, known as CRISPR interference, to block gene activity and in 2014 used it to systematically repress many of the genes in the human genome. Now, Horlbeck, Witkowsky et al. – who include several of the researchers from the 2014 work – have analyzed existing data for a specific type of human cell grown in the laboratory and found that CRISPR interference activity was strongest in certain areas around the start of each gene. However, CRISPR interference was much weaker in other areas of genes that coincided well with stretches of DNA that are known to often be bound by nucleosomes. Nucleosomes also appeared to block CRISPR editing, although the effects were less pronounced.

Horlbeck, Witkowsky et al. then directly tested whether nucleosomes could prevent the CRISPR system from binding or modifying the DNA. When the individual components were mixed in test tubes, the CRISPR system could readily target “naked” DNA. However, it could not access nucleosome-bound DNA, unless an enzyme that can move nucleosomes along the DNA in the human genome was also added to the mix. These findings suggest one way that CRISPR can manipulate much of the human genome despite the widespread presence of nucleosomes. Future work will now aim to develop computational methods that take the positions of nucleosomes into account when picking DNA sites to target with CRISPR.

DOI:http://dx.doi.org/10.7554/eLife.12677.002

Rheumatoid arthritis-associated RBPJ polymorphism alters memory CD4+ T cells

Notch signaling has recently emerged as an important regulator of immune responses in autoimmune diseases. The recombination signal-binding protein for immunoglobulin kappa J region (RBPJ) is a transcriptional repressor, but converts into a transcriptional activator upon activation of the canonical Notch pathway. Genome-wide association studies of rheumatoid arthritis (RA) identified a susceptibility locus, rs874040(CC), which implicated the RBPJ gene. Here, chromatin state mapping generated using the chromHMM algorithm reveals strong enhancer regions containing DNase I hypersensitive sites overlapping the rs874040 linkage disequilibrium block in human memory, but not in naïve CD4(+) T cells. The rs874040 overlapping this chromatin state was associated with increased RBPJ expression in stimulated memory CD4(+) T cells from healthy subjects homozygous for the risk allele (CC) compared with memory CD4(+) T cells bearing the protective allele (GG). Transcriptomic analysis of rs874040(CC) memory T cells showed a repression of canonical Notch target genes IL (interleukin)-9, IL-17 and interferon (IFN)γ in the basal state. Interestingly, activation of the Notch pathway using soluble Notch ligand, Jagged2-Fc, induced IL-9 and IL-17A while delta-like 4Fc, another Notch ligand, induced higher IFNγ expression in the rs874040(CC) memory CD4(+) T cells compared with their rs874040(GG) counterparts. In RA, RBPJ expression is elevated in memory T cells from RA patients compared with control subjects, and this was associated with induced inflammatory cytokines IL-9, IL-17A and IFNγ in response to Notch ligation in vitro. These findings demonstrate that the rs874040(CC) allele skews memory T cells toward a pro-inflammatory phenotype involving Notch signaling, thus increasing the susceptibility to develop RA.

CD33 modulates TREM2: convergence of Alzheimer loci

Here, we report results from a protein quantitative trait analysis in monocytes from 226 individuals to evaluate cross-talk between Alzheimer loci. We find that the NME8 locus influences PTK2B and that the CD33 risk allele leads to greater TREM2 expression. Further, we observe (1) a decreased TREM1/TREM2 ratio with a TREM1 risk allele, (2) decreased TREM2 expression with CD33 suppression, and (3) elevated cortical TREM2 mRNA expression with amyloid pathology.

A pharmacogenetic study implicates SLC9a9 in multiple sclerosis disease activity

OBJECTIVE

A proportion of multiple sclerosis (MS) patients experience disease activity despite treatment. The early identification of the most effective drug is critical to impact long-term outcome and to move toward a personalized approach. The aim of the present study is to identify biomarkers for further clinical development and to yield insights into the pathophysiology of disease activity.

METHODS

We performed a genome-wide association study in interferon-β (IFNβ)-treated MS patients followed by validation in 3 independent cohorts. The role of the validated variant was examined in several RNA data sets, and the function of the presumed target gene was explored using an RNA interference approach in primary T cells in vitro.

RESULTS

We found an association between rs9828519(G) and nonresponse to IFNβ (pdiscovery = 4.43 × 10(-8)) and confirmed it in a meta-analysis across 3 replication data sets (preplication = 7.78 × 10(-4)). Only 1 gene is found in the linkage disequilibrium block containing rs9828519: SLC9A9. Exploring the function of this gene, we see that SLC9A9 mRNA expression is diminished in MS subjects who are more likely to have relapses. Moreover, SLC9A9 knockdown in T cells in vitro leads an increase in expression of IFNγ, which is a proinflammatory molecule.

INTERPRETATION

This study identifies and validates the role of rs9828519, an intronic variant in SLC9A9, in IFNβ-treated subjects, demonstrating a successful pharmacogenetic screen in MS. Functional characterization suggests that SLC9A9, an Na(+) -H(+) exchanger found in endosomes, appears to influence the differentiation of T cells to a proinflammatory fate and may have a broader role in MS disease activity, outside of IFNβ treatment.

A TREM1 variant alters the accumulation of Alzheimer-related amyloid pathology

OBJECTIVE

Genome-wide association studies have linked variants in TREM2 (triggering receptor expressed on myeloid cells 2) and TREML2 with Alzheimer disease (AD) and AD endophenotypes. Here, we pursue a targeted analysis of the TREM locus in relation to cognitive decline and pathological features of AD.

METHODS

Clinical, cognitive, and neuropathological phenotypes were collected in 3 prospective cohorts on aging (n = 3,421 subjects). Our primary analysis was an association with neuritic plaque pathology. To functionally characterize the associated variants, we used flow cytometry to measure TREM1 expression on monocytes.

RESULTS

We provide evidence that an intronic variant, rs6910730(G) , in TREM1, is associated with an increased burden of neuritic plaques (p = 3.7 × 10(-4) ), diffuse plaques (p = 4.1 × 10(-3) ), and Aβ density (p = 2.6 × 10(-3) ) as well as an increased rate of cognitive decline (p = 5.3 × 10(-3) ). A variant upstream of TREM2, rs7759295(C) , is independently associated with an increased tau tangle density (p = 4.9 × 10(-4) ), an increased burden of neurofibrillary tangles (p = 9.1 × 10(-3) ), and an increased rate of cognitive decline (p = 2.3 × 10(-3) ). Finally, a cytometric analysis shows that the TREM1 rs6910730(G) allele is associated with decreased TREM1 expression on the surface of myeloid cells (p = 1.7 × 10(-3) ).

INTERPRETATION

We provide evidence that 2 common variants within the TREM locus are associated with pathological features of AD and aging-related cognitive decline. Our evidence suggests that these variants are likely to be independent of known AD variants and that they may work through an alteration of myeloid cell function.

Polarization of the effects of autoimmune and neurodegenerative risk alleles in leukocytes

To extend our understanding of the genetic basis of human immune function and dysfunction, we performed an expression quantitative trait locus (eQTL) study of purified CD4(+) T cells and monocytes, representing adaptive and innate immunity, in a multi-ethnic cohort of 461 healthy individuals. Context-specific cis- and trans-eQTLs were identified, and cross-population mapping allowed, in some cases, putative functional assignment of candidate causal regulatory variants for disease-associated loci. We note an over-representation of T cell-specific eQTLs among susceptibility alleles for autoimmune diseases and of monocyte-specific eQTLs among Alzheimer's and Parkinson's disease variants. This polarization implicates specific immune cell types in these diseases and points to the need to identify the cell-autonomous effects of disease susceptibility variants.

CD33: increased inclusion of exon 2 implicates the Ig V-set domain in Alzheimer's disease susceptibility

We previously demonstrated that the Alzheimer's disease (AD) associated risk allele, rs3865444(C), results in a higher surface density of CD33 on monocytes. Here, we find alternative splicing of exon 2 to be the primary mechanism of the genetically driven differential expression of CD33 protein. We report that the risk allele, rs3865444(C), is associated with greater cell surface expression of CD33 in both subjects of European and African-American ancestry and that there is a single haplotype influencing CD33 surface expression. A meta-analysis of the two populations narrowed the number of significant SNPs in high linkage disequilibrium (LD) (r(2) > 0.8) with rs3865444 to just five putative causal variants associated with increased protein expression. Using gene expression data from flow-sorted CD14(+)CD16(-) monocytes from 398 healthy subjects of three populations, we show that the rs3865444(C) risk allele is strongly associated with greater expression of CD33 exon 2 (pMETA = 2.36 × 10(-60)). Western blotting confirms increased protein expression of the full-length CD33 isoform containing exon 2 relative to the rs3865444(C) allele (P < 0.0001). Of the variants in strong LD with rs3865444, rs12459419, which is located in a putative SRSF2 splice site of exon 2, is the most likely candidate to mediate the altered alternative splicing of CD33's Immunoglobulin V-set domain 2 and ultimately influence AD susceptibility.