TLR7 gain-of-function genetic variation causes human lupus.

Although circumstantial evidence supports enhanced Toll-like receptor 7 (TLR7) signalling as a mechanism of human systemic autoimmune disease1-7, evidence of lupus-causing TLR7 gene variants is lacking. Here we describe human systemic lupus erythematosus caused by a TLR7 gain-of-function variant. TLR7 is a sensor of viral RNA8,9 and binds to guanosine10-12. We identified a de novo, previously undescribed missense TLR7Y264H variant in a child with severe lupus and additional variants in other patients with lupus. The TLR7Y264H variant selectively increased sensing of guanosine and 2',3'-cGMP10-12, and was sufficient to cause lupus when introduced into mice. We show that enhanced TLR7 signalling drives aberrant survival of B cell receptor (BCR)-activated B cells, and in a cell-intrinsic manner, accumulation of CD11c+ age-associated B cells and germinal centre B cells. Follicular and extrafollicular helper T cells were also increased but these phenotypes were cell-extrinsic. Deficiency of MyD88 (an adaptor protein downstream of TLR7) rescued autoimmunity, aberrant B cell survival, and all cellular and serological phenotypes. Despite prominent spontaneous germinal-centre formation in Tlr7Y264H mice, autoimmunity was not ameliorated by germinal-centre deficiency, suggesting an extrafollicular origin of pathogenic B cells. We establish the importance of TLR7 and guanosine-containing self-ligands for human lupus pathogenesis, which paves the way for therapeutic TLR7 or MyD88 inhibition.

Machine Learning Improves Upon Clinicians' Prediction of End Stage Kidney Disease.

Background and Objectives: Chronic kidney disease progression to ESKD is associated with a marked increase in mortality and morbidity. Its progression is highly variable and difficult to predict. Methods: This is an observational, retrospective, single-centre study. The cohort was patients attending hospital and nephrology clinic at The Canberra Hospital from September 1996 to March 2018. Demographic data, vital signs, kidney function test, proteinuria, and serum glucose were extracted. The model was trained on the featurised time series data with XGBoost. Its performance was compared against six nephrologists and the Kidney Failure Risk Equation (KFRE). Results: A total of 12,371 patients were included, with 2,388 were found to have an adequate density (three eGFR data points in the first 2 years) for subsequent analysis. Patients were divided into 80%/20% ratio for training and testing datasets.ML model had superior performance than nephrologist in predicting ESKD within 2 years with 93.9% accuracy, 60% sensitivity, 97.7% specificity, 75% positive predictive value. The ML model was superior in all performance metrics to the KFRE 4- and 8-variable models.eGFR and glucose were found to be highly contributing to the ESKD prediction performance. Conclusions: The computational predictions had higher accuracy, specificity and positive predictive value, which indicates the potential integration into clinical workflows for decision support.

Efficacy of computational predictions of the functional effect of idiosyncratic pharmacogenetic variants.

BACKGROUND: Pharmacogenetic variation is important to drug responses through diverse and complex mechanisms. Predictions of the functional impact of missense pharmacogenetic variants primarily rely on the degree of sequence conservation between species as a primary discriminator. However, idiosyncratic or off-target drug-variant interactions sometimes involve effects that are peripheral or accessory to the central systems in which a gene functions. Given the importance of sequence conservation to functional prediction tools-these idiosyncratic pharmacogenetic variants may violate the assumptions of predictive software commonly used to infer their effect. METHODS: Here we exhaustively assess the effectiveness of eleven missense mutation functional inference tools on all known pharmacogenetic missense variants contained in the Pharmacogenomics Knowledgebase (PharmGKB) repository. We categorize PharmGKB entries into sub-classes to catalog likely off-target interactions, such that we may compare predictions across different variant annotations. RESULTS: As previously demonstrated, functional inference tools perform variably across the complete set of PharmGKB variants, with large numbers of variants incorrectly classified as 'benign'. However, we find substantial differences amongst PharmGKB variant sub-classes, particularly in variants known to cause off-target, type B adverse drug reactions, that are largely unrelated to the main pharmacological action of the drug. Specifically, variants associated with off-target effects (hence referred to as off-target variants) were most often incorrectly classified as 'benign'. These results highlight the importance of understanding the underlying mechanism of pharmacogenetic variants and how variants associated with off-target effects will ultimately require new predictive algorithms. CONCLUSION: In this work we demonstrate that functional inference tools perform poorly on pharmacogenetic variants, particularly on subsets enriched for variants causing off-target, type B adverse drug reactions. We describe how to identify variants associated with off-target effects within PharmGKB in order to generate a training set of variants that is needed to develop new algorithms specifically for this class of variant. Development of such tools will lead to more accurate functional predictions and pave the way for the increased wide-spread adoption of pharmacogenetics in clinical practice.

A Point Mutation in IKAROS ZF1 Causes a B Cell Deficiency in Mice.

IKZF1 (IKAROS) is essential for normal lymphopoiesis in both humans and mice. Previous Ikzf1 mouse models have demonstrated the dual role for IKZF1 in both B and T cell development and have indicated differential requirements of each zinc finger. Furthermore, mutations in IKZF1 are known to cause common variable immunodeficiency in patients characterized by a loss of B cells and reduced Ab production. Through N-ethyl-N-nitrosourea mutagenesis, we have discovered a novel Ikzf1 mutant mouse with a missense mutation (L132P) in zinc finger 1 (ZF1) located in the DNA binding domain. Unlike other previously reported murine Ikzf1 mutations, this L132P point mutation (Ikzf1L132P ) conserves overall protein expression and has a B cell-specific phenotype with no effect on T cell development, indicating that ZF1 is not required for T cells. Mice have reduced Ab responses to immunization and show a progressive loss of serum Igs compared with wild-type littermates. IKZF1L132P overexpressed in NIH3T3 or HEK293T cells failed to localize to pericentromeric heterochromatin and bind target DNA sequences. Coexpression of wild-type and mutant IKZF1, however, allows for localization to pericentromeric heterochromatin and binding to DNA indicating a haploinsufficient mechanism of action for IKZF1L132P Furthermore, Ikzf1+/L132P mice have late onset defective Ig production, similar to what is observed in common variable immunodeficiency patients. RNA sequencing revealed a total loss of Hsf1 expression in follicular B cells, suggesting a possible functional link for the humoral immune response defects observed in Ikzf1L132P/L132P mice.

Detecting Causal Variants in Mendelian Disorders Using Whole-Genome Sequencing.

Increasingly affordable sequencing technologies are revolutionizing the field of genomic medicine. It is now feasible to interrogate all major classes of variation in an individual across the entire genome for less than $1000 USD. While the generation of patient sequence information using these technologies has become routine, the analysis and interpretation of this data remains the greatest obstacle to widespread clinical implementation. This chapter summarizes the steps to identify, annotate, and prioritize variant information required for clinical report generation. We discuss methods to detect each variant class and describe strategies to increase the likelihood of detecting causal variant(s) in Mendelian disease. Lastly, we describe a sample workflow for synthesizing large amount of genetic information into concise clinical reports.

Atypical B cells are part of an alternative lineage of B cells that participates in responses to vaccination and infection in humans.

The diversity of circulating human B cells is unknown. We use single-cell RNA sequencing (RNA-seq) to examine the diversity of both antigen-specific and total B cells in healthy subjects and malaria-exposed individuals. This reveals two B cell lineages: a classical lineage of activated and resting memory B cells and an alternative lineage, which includes previously described atypical B cells. Although atypical B cells have previously been associated with disease states, the alternative lineage is common in healthy controls, as well as malaria-exposed individuals. We further track Plasmodium-specific B cells after malaria vaccination in naive volunteers. We find that alternative lineage cells are primed after the initial immunization and respond to booster doses. However, alternative lineage cells develop an atypical phenotype with repeated boosts. The data highlight that atypical cells are part of a wider alternative lineage of B cells that are a normal component of healthy immune responses.

NINJ1 mediates plasma membrane rupture during lytic cell death.

Plasma membrane rupture (PMR) is the final cataclysmic event in lytic cell death. PMR releases intracellular molecules known as damage-associated molecular patterns (DAMPs) that propagate the inflammatory response1-3. The underlying mechanism of PMR, however, is unknown. Here we show that the cell-surface NINJ1 protein4-8, which contains two transmembrane regions, has an essential role in the induction of PMR. A forward-genetic screen of randomly mutagenized mice linked NINJ1 to PMR. Ninj1-/- macrophages exhibited impaired PMR in response to diverse inducers of pyroptotic, necrotic and apoptotic cell death, and were unable to release numerous intracellular proteins including HMGB1 (a known DAMP) and LDH (a standard measure of PMR). Ninj1-/- macrophages died, but with a distinctive and persistent ballooned morphology, attributable to defective disintegration of bubble-like herniations. Ninj1-/- mice were more susceptible than wild-type mice to infection with Citrobacter rodentium, which suggests a role for PMR in anti-bacterial host defence. Mechanistically, NINJ1 used an evolutionarily conserved extracellular domain for oligomerization and subsequent PMR. The discovery of NINJ1 as a mediator of PMR overturns the long-held idea that cell death-related PMR is a passive event.

Recurrent miscalling of missense variation from short-read genome sequence data.

BACKGROUND: Short-read resequencing of genomes produces abundant information of the genetic variation of individuals. Due to their numerous nature, these variants are rarely exhaustively validated. Furthermore, low levels of undetected variant miscalling will have a systematic and disproportionate impact on the interpretation of individual genome sequence information, especially should these also be carried through into in reference databases of genomic variation. RESULTS: We find that sequence variation from short-read sequence data is subject to recurrent-yet-intermittent miscalling that occurs in a sequence intrinsic manner and is very sensitive to sequence read length. The miscalls arise from difficulties aligning short reads to redundant genomic regions, where the rate of sequencing error approaches the sequence diversity between redundant regions. We find the resultant miscalled variants to be sensitive to small sequence variations between genomes, and thereby are often intrinsic to an individual, pedigree, strain or human ethnic group. In human exome sequences, we identify 2-300 recurrent false positive variants per individual, almost all of which are present in public databases of human genomic variation. From the exomes of non-reference strains of inbred mice, we identify 3-5000 recurrent false positive variants per mouse - the number of which increasing with greater distance between an individual mouse strain and the reference C57BL6 mouse genome. We show that recurrently miscalled variants may be reproduced for a given genome from repeated simulation rounds of read resampling, realignment and recalling. As such, it is possible to identify more than two-thirds of false positive variation from only ten rounds of simulation. CONCLUSION: Identification and removal of recurrent false positive variants from specific individual variant sets will improve overall data quality. Variant miscalls arising are highly sequence intrinsic and are often specific to an individual, pedigree or ethnicity. Further, read length is a strong determinant of whether given false variants will be called for any given genome - which has profound significance for cohort studies that pool datasets collected and sequenced at different points in time.

IRF2 transcriptionally induces GSDMD expression for pyroptosis.

Gasdermin-D (GSDMD) is cleaved by caspase-1, caspase-4, and caspase-11 in response to canonical and noncanonical inflammasome activation. Upon cleavage, GSDMD oligomerizes and forms plasma membrane pores, resulting in interleukin-1ß (IL-1ß) secretion, pyroptotic cell death, and inflammatory pathologies, including periodic fever syndromes and septic shock-a plague on modern medicine. Here, we showed that IRF2, a member of the interferon regulatory factor (IRF) family of transcription factors, was essential for the transcriptional activation of GSDMD. A forward genetic screen with N-ethyl-N-nitrosourea (ENU)-mutagenized mice linked IRF2 to inflammasome signaling. GSDMD expression was substantially attenuated in IRF2-deficient macrophages, endothelial cells, and multiple tissues, which corresponded with reduced IL-1ß secretion and inhibited pyroptosis. Mechanistically, IRF2 bound to a previously uncharacterized but unique site within the GSDMD promoter to directly drive GSDMD transcription for the execution of pyroptosis. Disruption of this single IRF2-binding site abolished signaling by both the canonical and noncanonical inflammasomes. Together, our data illuminate a key transcriptional mechanism for expression of the gene encoding GSDMD, a critical mediator of inflammatory pathologies.

Gain-of-function IKBKB mutation causes human combined immune deficiency.

Genetic mutations account for many devastating early onset immune deficiencies. In contrast, less severe and later onset immune diseases, including in patients with no prior family history, remain poorly understood. Whole exome sequencing in two cohorts of such patients identified a novel heterozygous de novo IKBKB missense mutation (c.607G>A) in two separate kindreds in whom probands presented with immune dysregulation, combined T and B cell deficiency, inflammation, and epithelial defects. IKBKB encodes IKK2, which activates NF-κB signaling. IKK2V203I results in enhanced NF-κB signaling, as well as T and B cell functional defects. IKK2V203 is a highly conserved residue, and to prove causation, we generated an accurate mouse model by introducing the precise orthologous codon change in Ikbkb using CRISPR/Cas9. Mice and humans carrying this missense mutation exhibit remarkably similar cellular and biochemical phenotypes. Accurate mouse models engineered by CRISPR/Cas9 can help characterize novel syndromes arising from de novo germline mutations and yield insight into pathogenesis.

DeepSNVMiner: a sequence analysis tool to detect emergent, rare mutations in subsets of cell populations.

Background. Massively parallel sequencing technology is being used to sequence highly diverse populations of DNA such as that derived from heterogeneous cell mixtures containing both wild-type and disease-related states. At the core of such molecule tagging techniques is the tagging and identification of sequence reads derived from individual input DNA molecules, which must be first computationally disambiguated to generate read groups sharing common sequence tags, with each read group representing a single input DNA molecule. This disambiguation typically generates huge numbers of reads groups, each of which requires additional variant detection analysis steps to be run specific to each read group, thus representing a significant computational challenge. While sequencing technologies for producing these data are approaching maturity, the lack of available computational tools for analysing such heterogeneous sequence data represents an obstacle to the widespread adoption of this technology. Results. Using synthetic data we successfully detect unique variants at dilution levels of 1 in a 1,000,000 molecules, and find DeeepSNVMiner obtains significantly lower false positive and false negative rates compared to popular variant callers GATK, SAMTools, FreeBayes and LoFreq, particularly as the variant concentration levels decrease. In a dilution series with genomic DNA from two cells lines, we find DeepSNVMiner identifies a known somatic variant when present at concentrations of only 1 in 1,000 molecules in the input material, the lowest concentration amongst all variant callers tested. Conclusions. Here we present DeepSNVMiner; a tool to disambiguate tagged sequence groups and robustly identify sequence variants specific to subsets of starting DNA molecules that may indicate the presence of a disease. DeepSNVMiner is an automated workflow of custom sequence analysis utilities and open source tools able to differentiate somatic DNA variants from artefactual sequence variants that likely arose during DNA amplification. The workflow remains flexible such that it may be customised to variants of the data production protocol used, and supports reproducible analysis through detailed logging and reporting of results. DeepSNVMiner is available for academic non-commercial research purposes at https://github.com/mattmattmattmatt/DeepSNVMiner.

Heterogeneity of Human Neutrophil CD177 Expression Results from CD177P1 Pseudogene Conversion.

Most humans harbor both CD177neg and CD177pos neutrophils but 1-10% of people are CD177null, placing them at risk for formation of anti-neutrophil antibodies that can cause transfusion-related acute lung injury and neonatal alloimmune neutropenia. By deep sequencing the CD177 locus, we catalogued CD177 single nucleotide variants and identified a novel stop codon in CD177null individuals arising from a single base substitution in exon 7. This is not a mutation in CD177 itself, rather the CD177null phenotype arises when exon 7 of CD177 is supplied entirely by the CD177 pseudogene (CD177P1), which appears to have resulted from allelic gene conversion. In CD177 expressing individuals the CD177 locus contains both CD177P1 and CD177 sequences. The proportion of CD177hi neutrophils in the blood is a heritable trait. Abundance of CD177hi neutrophils correlates with homozygosity for CD177 reference allele, while heterozygosity for ectopic CD177P1 gene conversion correlates with increased CD177neg neutrophils, in which both CD177P1 partially incorporated allele and paired intact CD177 allele are transcribed. Human neutrophil heterogeneity for CD177 expression arises by ectopic allelic conversion. Resolution of the genetic basis of CD177null phenotype identifies a method for screening for individuals at risk of CD177 isoimmunisation.

Reliably Detecting Clinically Important Variants Requires Both Combined Variant Calls and Optimized Filtering Strategies.

A diversity of tools is available for identification of variants from genome sequence data. Given the current complexity of incorporating external software into a genome analysis infrastructure, a tendency exists to rely on the results from a single tool alone. The quality of the output variant calls is highly variable however, depending on factors such as sequence library quality as well as the choice of short-read aligner, variant caller, and variant caller filtering strategy. Here we present a two-part study first using the high quality 'genome in a bottle' reference set to demonstrate the significant impact the choice of aligner, variant caller, and variant caller filtering strategy has on overall variant call quality and further how certain variant callers outperform others with increased sample contamination, an important consideration when analyzing sequenced cancer samples. This analysis confirms previous work showing that combining variant calls of multiple tools results in the best quality resultant variant set, for either specificity or sensitivity, depending on whether the intersection or union, of all variant calls is used respectively. Second, we analyze a melanoma cell line derived from a control lymphocyte sample to determine whether software choices affect the detection of clinically important melanoma risk-factor variants finding that only one of the three such variants is unanimously detected under all conditions. Finally, we describe a cogent strategy for implementing a clinical variant detection pipeline; a strategy that requires careful software selection, variant caller filtering optimizing, and combined variant calls in order to effectively minimize false negative variants. While implementing such features represents an increase in complexity and computation the results offer indisputable improvements in data quality.

Comparison of predicted and actual consequences of missense mutations.

Each person's genome sequence has thousands of missense variants. Practical interpretation of their functional significance must rely on computational inferences in the absence of exhaustive experimental measurements. Here we analyzed the efficacy of these inferences in 33 de novo missense mutations revealed by sequencing in first-generation progeny of N-ethyl-N-nitrosourea-treated mice, involving 23 essential immune system genes. PolyPhen2, SIFT, MutationAssessor, Panther, CADD, and Condel were used to predict each mutation's functional importance, whereas the actual effect was measured by breeding and testing homozygotes for the expected in vivo loss-of-function phenotype. Only 20% of mutations predicted to be deleterious by PolyPhen2 (and 15% by CADD) showed a discernible phenotype in individual homozygotes. Half of all possible missense mutations in the same 23 immune genes were predicted to be deleterious, and most of these appear to become subject to purifying selection because few persist between separate mouse substrains, rodents, or primates. Because defects in immune genes could be phenotypically masked in vivo by compensation and environment, we compared inferences by the same tools with the in vitro phenotype of all 2,314 possible missense variants in TP53; 42% of mutations predicted by PolyPhen2 to be deleterious (and 45% by CADD) had little measurable consequence for TP53-promoted transcription. We conclude that for de novo or low-frequency missense mutations found by genome sequencing, half those inferred as deleterious correspond to nearly neutral mutations that have little impact on the clinical phenotype of individual cases but will nevertheless become subject to purifying selection.

Reducing the search space for causal genetic variants with VASP.

MOTIVATION: Increasingly, cost-effective high-throughput DNA sequencing technologies are being utilized to sequence human pedigrees to elucidate the genetic cause of a wide variety of human diseases. While numerous tools exist for variant prioritization within a single genome, the ability to concurrently analyze variants within pedigrees remains a challenge, especially should there be no prior indication of the underlying genetic cause of the disease. Here, we present a tool, variant analysis of sequenced pedigrees (VASP), a flexible data integration environment capable of producing a summary of pedigree variation, providing relevant information such as compound heterozygosity, genome phasing and disease inheritance patterns. Designed to aggregate data across a sequenced pedigree, VASP allows both powerful filtering and custom prioritization of both single nucleotide variants (SNVs) and small indels. Hence, clinical and research users with prior knowledge of a disease are able to dramatically reduce the variant search space based on a wide variety of custom prioritization criteria. AVAILABILITY AND IMPLEMENTATION: Source code available for academic non-commercial research purposes at https://github.com/mattmattmattmatt/VASP.

Interplay of dFOXO and two ETS-family transcription factors determines lifespan in Drosophila melanogaster.

Forkhead box O (FoxO) transcription factors (TFs) are key drivers of complex transcriptional programmes that determine animal lifespan. FoxOs regulate a number of other TFs, but how these TFs in turn might mediate the anti-ageing programmes orchestrated by FoxOs in vivo is unclear. Here, we identify an E-twenty six (ETS)-family transcriptional repressor, Anterior open (Aop), as regulated by the single Drosophila melanogaster FoxO (dFOXO) in the adult gut. AOP, the functional orthologue of the human Etv6/Tel protein, binds numerous genomic sites also occupied by dFOXO and counteracts the activity of an ETS activator, Pointed (Pnt), to prevent the lifespan-shortening effects of co-activation of dFOXO and PNT. This detrimental synergistic effect of dFOXO and PNT appears to stem from a mis-regulation of lipid metabolism. At the same time, AOP activity in another fly organ, the fat body, has further beneficial roles, regulating genes in common with dfoxo, such as the secreted, non-sensory, odorant binding protein (Obp99b), and robustly extending lifespan. Our study reveals a complex interplay between evolutionarily conserved ETS factors and dFOXO, the functional significance of which may extend well beyond animal lifespan.

Identification of a pathogenic variant in TREX1 in early-onset cerebral systemic lupus erythematosus by Whole-exome sequencing.

Objective. Systemic lupus erythematosus (SLE) isa chronic and heterogeneous autoimmune disease. Both twin and sibling studies indicate a strong genetic contribution to lupus, but in the majority of cases the pathogenic variant remains to be identified. The genetic contribution to disease is likely to be greatest in cases with early onset and severe phenotypes. Whole-exome sequencing now offers the possibility of identifying rare alleles responsible for disease in such cases. This study was undertaken to identify genetic causes of SLE using whole-exome sequencing.Methods. We performed whole-exome sequencing in a 4-year-old girl with early-onset SLE and conducted biochemical analysis of the putative defect.Results. Whole-exome sequencing in a 4-year-old girl with cerebral lupus identified a rare, homozygous mutation in the three prime repair exonuclease 1 gene(TREX1) that was predicted to be highly deleterious.The TREX1 R97H mutant protein had a 20-fold reduction in exonuclease activity and was associated with an elevated interferon-alpha signature in the patient.The discovery and characterization of a pathogenic TREX1 variant in our proband has therapeutic implications.The patient is now a candidate for therapy. Conclusion. Our study is the first to demonstrate that whole-exome sequencing can be used to identify rare or novel deleterious variants as genetic causes of SLE and, through a personalized approach, improve therapeutic options.

Zinc-finger protein ZFP318 is essential for expression of IgD, the alternatively spliced Igh product made by mature B lymphocytes.

IgD and IgM are produced by alternative splicing of long primary RNA transcripts from the Ig heavy chain (Igh) locus and serve as the receptors for antigen on naïve mature B lymphocytes. IgM is made selectively in immature B cells, whereas IgD is coexpressed with IgM when the cells mature into follicular or marginal zone B cells, but the transacting factors responsible for this regulated change in splicing have remained elusive. Here, we use a genetic screen in mice to identify ZFP318, a nuclear protein with two U1-type zinc fingers found in RNA-binding proteins and no known role in the immune system, as a critical factor for IgD expression. A point mutation in an evolutionarily conserved lysine-rich domain encoded by the alternatively spliced Zfp318 exon 10 abolished IgD expression on marginal zone B cells, decreased IgD on follicular B cells, and increased IgM, but only slightly decreased the percentage of B cells and did not decrease expression of other maturation markers CD21, CD23, or CD62L. A targeted Zfp318 null allele extinguished IgD expression on mature B cells and increased IgM. Zfp318 mRNA is developmentally regulated in parallel with IgD, with little in pro-B cells, moderate amounts in immature B cells, and high levels selectively in mature follicular B cells. These findings identify ZFP318 as a crucial factor regulating the expression of the two major antibody isotypes on the surface of most mature B cells.

The RNA-binding protein hnRNPLL induces a T cell alternative splicing program delineated by differential intron retention in polyadenylated RNA.

BACKGROUND: Retention of a subset of introns in spliced polyadenylated mRNA is emerging as a frequent, unexplained finding from RNA deep sequencing in mammalian cells. RESULTS: Here we analyze intron retention in T lymphocytes by deep sequencing polyadenylated RNA. We show a developmentally regulated RNA-binding protein, hnRNPLL, induces retention of specific introns by sequencing RNA from T cells with an inactivating Hnrpll mutation and from B lymphocytes that physiologically downregulate Hnrpll during their differentiation. In Ptprc mRNA encoding the tyrosine phosphatase CD45, hnRNPLL induces selective retention of introns flanking exons 4 to 6; these correspond to the cassette exons containing hnRNPLL binding sites that are skipped in cells with normal, but not mutant or low, hnRNPLL. We identify similar patterns of hnRNPLL-induced differential intron retention flanking alternative exons in 14 other genes, representing novel elements of the hnRNPLL-induced splicing program in T cells. Retroviral expression of a normally spliced cDNA for one of these targets, Senp2, partially corrects the survival defect of Hnrpll-mutant T cells. We find that integrating a number of computational methods to detect genes with differentially retained introns provides a strategy to enrich for alternatively spliced exons in mammalian RNA-seq data, when complemented by RNA-seq analysis of purified cells with experimentally perturbed RNA-binding proteins. CONCLUSIONS: Our findings demonstrate that intron retention in mRNA is induced by specific RNA-binding proteins and suggest a biological significance for this process in marking exons that are poised for alternative splicing.

Rasgrp1 mutation increases naive T-cell CD44 expression and drives mTOR-dependent accumulation of Helios⁺ T cells and autoantibodies.

Missense variants are a major source of human genetic variation. Here we analyze a new mouse missense variant, Rasgrp1(Anaef), with an ENU-mutated EF hand in the Rasgrp1 Ras guanine nucleotide exchange factor. Rasgrp1(Anaef) mice exhibit anti-nuclear autoantibodies and gradually accumulate a CD44(hi) Helios(+) PD-1(+) CD4(+) T cell population that is dependent on B cells. Despite reduced Rasgrp1-Ras-ERK activation in vitro, thymocyte selection in Rasgrp1(Anaef) is mostly normal in vivo, although CD44 is overexpressed on naïve thymocytes and T cells in a T-cell-autonomous manner. We identify CD44 expression as a sensitive reporter of tonic mTOR-S6 kinase signaling through a novel mouse strain, chino, with a reduction-of-function mutation in Mtor. Elevated tonic mTOR-S6 signaling occurs in Rasgrp1(Anaef) naïve CD4(+) T cells. CD44 expression, CD4(+) T cell subset ratios and serum autoantibodies all returned to normal in Rasgrp1(Anaef)Mtor(chino) double-mutant mice, demonstrating that increased mTOR activity is essential for the Rasgrp1(Anaef) T cell dysregulation. DOI: http://dx.doi.org/10.7554/eLife.01020.001.