Non-coding sequence variation reveals fragility within interleukin 2 feedback circuitry and shapes autoimmune disease risk.

Genetic variants associated with human autoimmune diseases commonly map to non-coding control regions, particularly enhancers that function selectively in immune cells and fine-tune gene expression within a relatively narrow range of values. How such modest, cell-type-selective changes can meaningfully shape organismal disease risk remains unclear. To explore this issue, we experimentally manipulated species-conserved enhancers within the disease-associated IL2RA locus and studied accompanying changes in the progression of autoimmunity. Perturbing distinct enhancers with restricted activity in conventional T cells (Tconvs) or regulatory T cells (Tregs)-two functionally antagonistic T cell subsets-caused only modest, cell-type-selective decreases in IL2ra expression parameters. However, these same perturbations had striking and opposing effects in vivo , completely preventing or severely accelerating disease in a murine model of type 1 diabetes. Quantitative tissue imaging and computational modelling revealed that each enhancer manipulation impinged on distinct IL-2-dependent feedback circuits. These imbalances altered the intracellular signaling and intercellular communication dynamics of activated Tregs and Tconvs, producing opposing spatial domains that amplified or constrained ongoing autoimmune responses. These findings demonstrate how subtle changes in gene regulation stemming from non-coding variation can propagate across biological scales due to non-linearities in intra- and intercellular feedback circuitry, dramatically shaping disease risk at the organismal level.

A functional map of HIV-host interactions in primary human T cells.

Human Immunodeficiency Virus (HIV) relies on host molecular machinery for replication. Systematic attempts to genetically or biochemically define these host factors have yielded hundreds of candidates, but few have been functionally validated in primary cells. Here, we target 426 genes previously implicated in the HIV lifecycle through protein interaction studies for CRISPR-Cas9-mediated knock-out in primary human CD4+ T cells in order to systematically assess their functional roles in HIV replication. We achieve efficient knockout (>50% of alleles) in 364 of the targeted genes and identify 86 candidate host factors that alter HIV infection. 47 of these factors validate by multiplex gene editing in independent donors, including 23 factors with restrictive activity. Both gene editing efficiencies and HIV-1 phenotypes are highly concordant among independent donors. Importantly, over half of these factors have not been previously described to play a functional role in HIV replication, providing numerous novel avenues for understanding HIV biology. These data further suggest that host-pathogen protein-protein interaction datasets offer an enriched source of candidates for functional host factor discovery and provide an improved understanding of the mechanics of HIV replication in primary T cells.

XYZeq: Spatially resolved single-cell RNA sequencing reveals expression heterogeneity in the tumor microenvironment.

Single-cell RNA sequencing (scRNA-seq) of tissues has revealed remarkable heterogeneity of cell types and states but does not provide information on the spatial organization of cells. To better understand how individual cells function within an anatomical space, we developed XYZeq, a workflow that encodes spatial metadata into scRNA-seq libraries. We used XYZeq to profile mouse tumor models to capture spatially barcoded transcriptomes from tens of thousands of cells. Analyses of these data revealed the spatial distribution of distinct cell types and a cell migration-associated transcriptomic program in tumor-associated mesenchymal stem cells (MSCs). Furthermore, we identify localized expression of tumor suppressor genes by MSCs that vary with proximity to the tumor core. We demonstrate that XYZeq can be used to map the transcriptome and spatial localization of individual cells in situ to reveal how cell composition and cell states can be affected by location within complex pathological tissue.

Functional CRISPR dissection of gene networks controlling human regulatory T cell identity.

Human regulatory T (Treg) cells are essential for immune homeostasis. The transcription factor FOXP3 maintains Treg cell identity, yet the complete set of key transcription factors that control Treg cell gene expression remains unknown. Here, we used pooled and arrayed Cas9 ribonucleoprotein screens to identify transcription factors that regulate critical proteins in primary human Treg cells under basal and proinflammatory conditions. We then generated 54,424 single-cell transcriptomes from Treg cells subjected to genetic perturbations and cytokine stimulation, which revealed distinct gene networks individually regulated by FOXP3 and PRDM1, in addition to a network coregulated by FOXO1 and IRF4. We also discovered that HIVEP2, to our knowledge not previously implicated in Treg cell function, coregulates another gene network with SATB1 and is important for Treg cell-mediated immunosuppression. By integrating CRISPR screens and single-cell RNA-sequencing profiling, we have uncovered transcriptional regulators and downstream gene networks in human Treg cells that could be targeted for immunotherapies.

Evaluation of SARS-CoV-2 serology assays reveals a range of test performance.

Appropriate use and interpretation of serological tests for assessments of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) exposure, infection and potential immunity require accurate data on assay performance. We conducted a head-to-head evaluation of ten point-of-care-style lateral flow assays (LFAs) and two laboratory-based enzyme-linked immunosorbent assays to detect anti-SARS-CoV-2 IgM and IgG antibodies in 5-d time intervals from symptom onset and studied the specificity of each assay in pre-coronavirus disease 2019 specimens. The percent of seropositive individuals increased with time, peaking in the latest time interval tested (>20 d after symptom onset). Test specificity ranged from 84.3% to 100.0% and was predominantly affected by variability in IgM results. LFA specificity could be increased by considering weak bands as negative, but this decreased detection of antibodies (sensitivity) in a subset of SARS-CoV-2 real-time PCR-positive cases. Our results underline the importance of seropositivity threshold determination and reader training for reliable LFA deployment. Although there was no standout serological assay, four tests achieved more than 80% positivity at later time points tested and more than 95% specificity.

CRISPR screen in regulatory T cells reveals modulators of Foxp3.

Regulatory T (Treg) cells are required to control immune responses and maintain homeostasis, but are a significant barrier to antitumour immunity1. Conversely, Treg instability, characterized by loss of the master transcription factor Foxp3 and acquisition of proinflammatory properties2, can promote autoimmunity and/or facilitate more effective tumour immunity3,4. A comprehensive understanding of the pathways that regulate Foxp3 could lead to more effective Treg therapies for autoimmune disease and cancer. The availability of new functional genetic tools has enabled the possibility of systematic dissection of the gene regulatory programs that modulate Foxp3 expression. Here we developed a CRISPR-based pooled screening platform for phenotypes in primary mouse Treg cells and applied this technology to perform a targeted loss-of-function screen of around 500 nuclear factors to identify gene regulatory programs that promote or disrupt Foxp3 expression. We identified several modulators of Foxp3 expression, including ubiquitin-specific peptidase 22 (Usp22) and ring finger protein 20 (Rnf20). Usp22, a member of the deubiquitination module of the SAGA chromatin-modifying complex, was revealed to be a positive regulator that stabilized Foxp3 expression; whereas the screen suggested that Rnf20, an E3 ubiquitin ligase, can serve as a negative regulator of Foxp3. Treg-specific ablation of Usp22 in mice reduced Foxp3 protein levels and caused defects in their suppressive function that led to spontaneous autoimmunity but protected against tumour growth in multiple cancer models. Foxp3 destabilization in Usp22-deficient Treg cells could be rescued by ablation of Rnf20, revealing a reciprocal ubiquitin switch in Treg cells. These results reveal previously unknown modulators of Foxp3 and demonstrate a screening method that can be broadly applied to discover new targets for Treg immunotherapies for cancer and autoimmune disease.

Test performance evaluation of SARS-CoV-2 serological assays.

BACKGROUND: Serological tests are crucial tools for assessments of SARS-CoV-2 exposure, infection and potential immunity. Their appropriate use and interpretation require accurate assay performance data. METHOD: We conducted an evaluation of 10 lateral flow assays (LFAs) and two ELISAs to detect anti-SARS-CoV-2 antibodies. The specimen set comprised 128 plasma or serum samples from 79 symptomatic SARS-CoV-2 RT-PCR-positive individuals; 108 pre-COVID-19 negative controls; and 52 recent samples from individuals who underwent respiratory viral testing but were not diagnosed with Coronavirus Disease 2019 (COVID-19). Samples were blinded and LFA results were interpreted by two independent readers, using a standardized intensity scoring system. RESULTS: Among specimens from SARS-CoV-2 RT-PCR-positive individuals, the percent seropositive increased with time interval, peaking at 81.8-100.0% in samples taken >20 days after symptom onset. Test specificity ranged from 84.3-100.0% in pre-COVID-19 specimens. Specificity was higher when weak LFA bands were considered negative, but this decreased sensitivity. IgM detection was more variable than IgG, and detection was highest when IgM and IgG results were combined. Agreement between ELISAs and LFAs ranged from 75.7-94.8%. No consistent cross-reactivity was observed. CONCLUSION: Our evaluation showed heterogeneous assay performance. Reader training is key to reliable LFA performance, and can be tailored for survey goals. Informed use of serology will require evaluations covering the full spectrum of SARS-CoV-2 infections, from asymptomatic and mild infection to severe disease, and later convalescence. Well-designed studies to elucidate the mechanisms and serological correlates of protective immunity will be crucial to guide rational clinical and public health policies.

A large CRISPR-induced bystander mutation causes immune dysregulation.

A persistent concern with CRISPR-Cas9 gene editing has been the potential to generate mutations at off-target genomic sites. While CRISPR-engineering mice to delete a ~360 bp intronic enhancer, here we discovered a founder line that had marked immune dysregulation caused by a 24 kb tandem duplication of the sequence adjacent to the on-target deletion. Our results suggest unintended repair of on-target genomic cuts can cause pathogenic "bystander" mutations that escape detection by routine targeted genotyping assays.

Genome-wide CRISPR Screens in Primary Human T Cells Reveal Key Regulators of Immune Function.

Human T cells are central effectors of immunity and cancer immunotherapy. CRISPR-based functional studies in T cells could prioritize novel targets for drug development and improve the design of genetically reprogrammed cell-based therapies. However, large-scale CRISPR screens have been challenging in primary human cells. We developed a new method, single guide RNA (sgRNA) lentiviral infection with Cas9 protein electroporation (SLICE), to identify regulators of stimulation responses in primary human T cells. Genome-wide loss-of-function screens identified essential T cell receptor signaling components and genes that negatively tune proliferation following stimulation. Targeted ablation of individual candidate genes characterized hits and identified perturbations that enhanced cancer cell killing. SLICE coupled with single-cell RNA sequencing (RNA-seq) revealed signature stimulation-response gene programs altered by key genetic perturbations. SLICE genome-wide screening was also adaptable to identify mediators of immunosuppression, revealing genes controlling responses to adenosine signaling. The SLICE platform enables unbiased discovery and characterization of functional gene targets in primary cells.

Author Correction: Discovery of stimulation-responsive immune enhancers with CRISPR activation.

In this Letter, analysis of steady-state regulatory T (Treg) cell percentages from Il2ra enhancer deletion (EDEL) and wild-type (WT) mice revealed no differences between them (Extended Data Fig. 9d). This analysis included two mice whose genotypes were incorrectly assigned. Even after correction of the genotypes, no significant differences in Treg cell percentages were seen when data across experimental cohorts were averaged (as was done in Extended Data Fig. 9d). However, if we normalize the corrected data to account for variation among experimental cohorts, a subtle decrease in EDEL Treg cell percentages is revealed and, using the corrected and normalized data, we have redrawn Extended Data Fig. 9d in Supplementary Fig. 1. The Supplementary Information to this Amendment contains the corrected and reanalysed Extended Data Fig. 9d. The sentence "This enhancer deletion (EDEL) strain also had no obvious T cell phenotypes at steady state (Extended Data Fig. 9)." should read: "This enhancer deletion (EDEL) strain had a small decrease in the percentage of Treg cells (Extended Data Fig. 9).". This error does not affect any of the main figures in the Letter or the data from mice with the human autoimmune-associated single nucleotide polymorphism (SNP) knocked in or with a 12-base-pair deletion at the site (12DEL). In addition, we stated in the Methods that we observed consistent immunophenotypes of EDEL mice across three founders, but in fact, we observed consistent phenotypes in mice from two founders. This does not change any of our conclusions and the original Letter has not been corrected.

Discovery of stimulation-responsive immune enhancers with CRISPR activation.

The majority of genetic variants associated with common human diseases map to enhancers, non-coding elements that shape cell-type-specific transcriptional programs and responses to extracellular cues. Systematic mapping of functional enhancers and their biological contexts is required to understand the mechanisms by which variation in non-coding genetic sequences contributes to disease. Functional enhancers can be mapped by genomic sequence disruption, but this approach is limited to the subset of enhancers that are necessary in the particular cellular context being studied. We hypothesized that recruitment of a strong transcriptional activator to an enhancer would be sufficient to drive target gene expression, even if that enhancer was not currently active in the assayed cells. Here we describe a discovery platform that can identify stimulus-responsive enhancers for a target gene independent of stimulus exposure. We used tiled CRISPR activation (CRISPRa) to synthetically recruit a transcriptional activator to sites across large genomic regions (more than 100 kilobases) surrounding two key autoimmunity risk loci, CD69 and IL2RA. We identified several CRISPRa-responsive elements with chromatin features of stimulus-responsive enhancers, including an IL2RA enhancer that harbours an autoimmunity risk variant. Using engineered mouse models, we found that sequence perturbation of the disease-associated Il2ra enhancer did not entirely block Il2ra expression, but rather delayed the timing of gene activation in response to specific extracellular signals. Enhancer deletion skewed polarization of naive T cells towards a pro-inflammatory T helper (TH17) cell state and away from a regulatory T cell state. This integrated approach identifies functional enhancers and reveals how non-coding variation associated with human immune dysfunction alters context-specific gene programs.

A Cas9 Ribonucleoprotein Platform for Functional Genetic Studies of HIV-Host Interactions in Primary Human T Cells.

New genetic tools are needed to understand the functional interactions between HIV and human host factors in primary cells. We recently developed a method to edit the genome of primary CD4+ T cells by electroporation of CRISPR/Cas9 ribonucleoproteins (RNPs). Here, we adapted this methodology to a high-throughput platform for the efficient, arrayed editing of candidate host factors. CXCR4 or CCR5 knockout cells generated with this method are resistant to HIV infection in a tropism-dependent manner, whereas knockout of LEDGF or TNPO3 results in a tropism-independent reduction in infection. CRISPR/Cas9 RNPs can furthermore edit multiple genes simultaneously, enabling studies of interactions among multiple host and viral factors. Finally, in an arrayed screen of 45 genes associated with HIV integrase, we identified several candidate dependency/restriction factors, demonstrating the power of this approach as a discovery platform. This technology should accelerate target validation for pharmaceutical and cell-based therapies to cure HIV infection.

A pharmacogenetic candidate gene study of tenofovir-associated Fanconi syndrome.

BACKGROUND: Tenofovir disoproxil fumarate (TDF) is a widely used antiretroviral agent with favorable efficacy, safety, and tolerability profiles. However, renal adverse events, including the rare Fanconi syndrome (FS), may occur in a small subset of patients treated for HIV infections. OBJECTIVES: The aim of this study was to identify genetic variants that may be associated with TDF-associated FS (TDF-FS). METHODS: DNA samples collected from 19 cases with TDF-FS and 36 matched controls were sequenced, and genetic association studies were conducted on eight candidate genes: ATP-binding cassette (ABC) transporters ABCC2 (MRP2) and ABCC4 (MRP4), solute carrier family members SLC22A6 (OAT1) and SLC22A8 (OAT3), adenylate kinases 2 (AK2) and 4 (AK4), chloride transporter CIC-5 CLCN5, and Lowe syndrome protein OCRL. The functional effects of a single nucleotide polymorphism (SNP) predicted to alter the transport of tenofovir were then investigated in cells expressing an identified variant of ABCC4. RESULTS: The case group showed a trend toward a higher proportion of rare alleles. Six SNPs in ABCC2 (three SNPs), ABCC4 (one SNP), and OCRL (two SNPs) were associated with TDF-FS case status; however, this association did not remain significant after correction for multiple testing. Six SNPs, present in OCRL (four SNPs) and ABCC2 (two SNPs), were significantly associated with increased serum creatinine levels in the cases, and this association remained significant after multiple test correction (P < 2 × 10). One synonymous SNP in ABCC2 (rs8187707, P = 2.10 × 10, ß = -73.3 ml/min/1.73 m(2)) was also significantly associated with the decreased estimated glomerular filtration rate of creatinine among cases. However, these results were driven by rare SNPs present in a small number of severely affected cases. Finally, a previously uncharacterized, nonsynonymous SNP, rs11568694, that was predicted to alter MRP4 function had no significant effect on tenofovir cellular accumulation in vitro. CONCLUSION: Although no single predictive genetic marker for the development of TDF-FS was identified, the findings from our study suggest that rare variants in multiple genes involved in the renal handling of tenofovir, and/or renal cell homeostasis, may be associated with increased susceptibility to TDF-FS.

Genomewide linkage analysis in Costa Rican families implicates chromosome 15q14 as a candidate region for OCD.

Obsessive compulsive disorder (OCD) has a complex etiology that encompasses both genetic and environmental factors. However, to date, despite the identification of several promising candidate genes and linkage regions, the genetic causes of OCD are largely unknown. The objective of this study was to conduct linkage studies of childhood-onset OCD, which is thought to have the strongest genetic etiology, in several OCD-affected families from the genetically isolated population of the Central Valley of Costa Rica (CVCR). The authors used parametric and non-parametric approaches to conduct genome-wide linkage analyses using 5,786 single nucleotide repeat polymorphisms (SNPs) in three CVCR families with multiple childhood-onset OCD-affected individuals. We identified areas of suggestive linkage (LOD score ≥ 2) on chromosomes 1p21, 15q14, 16q24, and 17p12. The strongest evidence for linkage was on chromosome 15q14 (LOD = 3.13), identified using parametric linkage analysis with a recessive model, and overlapping a region identified in a prior linkage study using a Caucasian population. Each CVCR family had a haplotype that co-segregated with OCD across a ~7 Mbp interval within this region, which contains 18 identified brain expressed genes, several of which are potentially relevant to OCD. Exonic sequencing of the strongest candidate gene in this region, the ryanodine receptor 3 (RYR3), identified several genetic variants of potential interest, although none co-segregated with OCD in all three families. These findings provide evidence that chromosome 15q14 is linked to OCD in families from the CVCR, and supports previous findings to suggest that this region may contain one or more OCD susceptibility loci.

Genetic profiling reveals cross-contamination and misidentification of 6 adenoid cystic carcinoma cell lines: ACC2, ACC3, ACCM, ACCNS, ACCS and CAC2.

Adenoid cystic carcinoma (ACC) is the second most common malignant neoplasm of the salivary glands. Most patients survive more than 5 years after surgery and postoperative radiation therapy. The 10 year survival rate, however, drops to 40%, due to locoregional recurrences and distant metastases. Improving long-term survival in ACC requires the development of more effective systemic therapies based on a better understanding of the biologic behavior of ACC. Much preclinical research in this field involves the use of cultured cells and, to date, several ACC cell lines have been established. Authentication of these cell lines, however, has not been reported. We performed DNA fingerprint analysis on six ACC cell lines using short tandem repeat (STR) examinations and found that all six cell lines had been contaminated with other cells. ACC2, ACC3, and ACCM were determined to be cervical cancer cells (HeLa cells), whereas the ACCS cell line was composed of T24 urinary bladder cancer cells. ACCNS and CAC2 cells were contaminated with cells derived from non-human mammalian species: the cells labeled ACCNS were mouse cells and the CAC2 cells were rat cells. These observations suggest that future studies using ACC cell lines should include cell line authentication to avoid the use of contaminated or non-human cells.

Validation of microsatellite markers for use in genotyping polyclonal Plasmodium falciparum infections.

Genotyping methods for Plasmodium falciparum drug efficacy trials have not been standardized and may fail to accurately distinguish recrudescence from new infection, especially in high transmission areas where polyclonal infections are common. We developed a simple method for genotyping using previously identified microsatellites and capillary electrophoresis, validated this method using mixtures of laboratory clones, and applied the method to field samples. Two microsatellite markers produced accurate results for single-clone but not polyclonal samples. Four other microsatellite markers were as sensitive as, and more specific than, commonly used genotyping techniques based on merozoite surface proteins 1 and 2. When applied to samples from 15 patients in Burkina Faso with recurrent parasitemia after treatment with sulphadoxine-pyrimethamine, the addition of these four microsatellite markers to msp1 and msp2 genotyping resulted in a reclassification of outcomes that strengthened the association between dhfr 59R, an anti-folate resistance mutation, and recrudescence (P = 0.31 versus P = 0.03). Four microsatellite markers performed well on polyclonal samples and may provide a valuable addition to genotyping for clinical drug efficacy studies in high transmission areas.

Analysis of four DLX homeobox genes in autistic probands.

BACKGROUND: Linkage studies in autism have identified susceptibility loci on chromosomes 2q and 7q, regions containing the DLX1/2 and DLX5/6 bigene clusters. The DLX genes encode homeodomain transcription factors that control craniofacial patterning and differentiation and survival of forebrain inhibitory neurons. We investigated the role that sequence variants in DLX genes play in autism by in-depth resequencing of these genes in 161 autism probands from the AGRE collection. RESULTS: Sequencing of exons, exon/intron boundaries and known enhancers of DLX1, 2, 5 and 6 identified several nonsynonymous variants in DLX2 and DLX5 and a variant in a DLX5/6 intragenic enhancer. The nonsynonymous variants were detected in 4 of 95 families from which samples were sequenced. Two of these four SNPs were not observed in 378 undiagnosed samples from North American populations, while the remaining 2 were seen in one sample each. CONCLUSION: Segregation of these variants in pedigrees did not generally support a contribution to autism susceptibility by these genes, although functional analyses may provide insight into the biological understanding of these important proteins.

Effect of the MDR1 C3435T variant and P-glycoprotein induction on dicloxacillin pharmacokinetics.

This study investigated 2 hypotheses about genotype-phenotype relationships for the efflux transporter, P-glycoprotein: (1) the presence of a synonymous C3435T variant in exon 26 of the MDR1 gene correlates to higher plasma concentrations of a P-glycoprotein substrate, dicloxacillin, and (2) the effects of genotypic differences decrease under conditions of P-glycoprotein induction by rifampin. Eighteen healthy volunteers received two 1-g doses of dicloxacillin, one on the 1st study day and the other on the 11th day of rifampin dosing (600 mg daily). Dicloxacillin and its 5-hydroxymethyl metabolite were analyzed using liquid chromatography/tandem mass spectrometry. Mean dicloxacillin C(max) measurements were 30.5 +/- 13.5, 33.3 +/- 4.7, and 31.1 +/- 12.8 mug/mL in individuals with the CC, CT, and TT genotype at position 3435 in exon 26 of the MDR1 gene. Following rifampin dosing, the mean dicloxacillin C(max) across genotypes decreased from 31.4 +/- 10.8 to 22.9 +/- 7.0 microg/mL (P < .05), whereas the mean oral clearance increased from 235 +/- 82 to 297 +/- 71 mL/min (P < .001), and the mean absorption time increased from 0.71 +/- 0.55 to 1.34 +/- 0.77 h (P < .05). Rifampin treatment increased the formation clearance, C(max), and AUC of the 5-hydroxymethyl metabolite by 135%, 119%, and 59%, respectively. The C3435T variant had no effect on dicloxacillin pharmacokinetics. The data suggested that rifampin induced intestinal P-glycoprotein and increased dicloxacillin metabolism.