Search

Whole-exome sequencing of two unrelated kindreds with systemic autoimmune disease featuring antinuclear antibodies with IgG4 elevation uncovered an identical ultrarare heterozygous TNIP1 Q333P variant segregating with disease. Mice with the orthologous Q346P variant developed antinuclear autoantibodies, salivary gland inflammation, elevated IgG2c, spontaneous germinal centers and expansion of age-associated B cells, plasma cells and follicular and extrafollicular helper T cells. B cell phenotypes were cell-autonomous and rescued by ablation of Toll-like receptor 7 (TLR7) or MyD88. The variant increased interferon-β without altering nuclear factor kappa-light-chain-enhancer of activated B cells signaling, and impaired MyD88 and IRAK1 recruitment to autophagosomes. Additionally, the Q333P variant impaired TNIP1 localization to damaged mitochondria and mitophagosome formation. Damaged mitochondria were abundant in the salivary epithelial cells of Tnip1 Q346P mice. These findings suggest that TNIP1-mediated autoimmunity may be a consequence of increased TLR7 signaling due to impaired recruitment of downstream signaling molecules and damaged mitochondria to autophagosomes and may thus respond to TLR7-targeted therapeutics., (© 2024. The Author(s).)

Although circumstantial evidence supports enhanced Toll-like receptor 7 (TLR7) signalling as a mechanism of human systemic autoimmune disease 1-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 RNA 8 , 9 and binds to guanosine 10 - 12 . We identified a de novo, previously undescribed missense TLR7 Y264H variant in a child with severe lupus and additional variants in other patients with lupus. The TLR7 Y264H variant selectively increased sensing of guanosine and 2',3'-cGMP 10-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 Tlr7 Y264H 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., (© 2022. The Author(s).)

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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Chuah, Walters, Christiadi, Karpe, Kennard, Singer, Talaulikar, Ge, Suominen, Andrews and Jiang.)

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., Competing Interests: The authors declare that they have no competing interests., (© 2021 McConnell et al.)

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 ( Ikzf1 L132P ) 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. IKZF1 L132P 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 IKZF1 L132P 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 Ikzf1 L132P/L132P mice., (Copyright © 2021 by The American Association of Immunologists, Inc.)

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 response 1-3 . The underlying mechanism of PMR, however, is unknown. Here we show that the cell-surface NINJ1 protein 4-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.

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., Competing Interests: Declaration of interests S.C., N.K., B.K.L.S., and S.L.H. are salaried employees of Sanaria, the developer and owner of PfSPZ Vaccine and the investigational new drug (IND) application sponsor of the clinical trials. S.L.H. and B.K.L.S. have a financial interest in Sanaria. All other authors declare no competing interests., (Copyright © 2021 The Author(s). Published by Elsevier Inc. All rights reserved.)

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.

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.

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., (Copyright © 2019 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

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. IKK2 V203I results in enhanced NF-κB signaling, as well as T and B cell functional defects. IKK2 V203 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., (© 2018 Crown copyright. The government of Australia, Canada, or the UK ("the Crown") owns the copyright interests of authors who are government employees. The Crown Copyright is not transferable.)

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.

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.

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.

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.

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., (© The Author 2015. Published by Oxford University Press.)

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.

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.

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.

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.

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.

Druggable proteins required for B lymphocyte survival and immune responses are an emerging source of new treatments for autoimmunity and lymphoid malignancy. In this study, we show that mice with an inactivating mutation in the intramembrane protease signal peptide peptidase-like 2A (SPPL2A) unexpectedly exhibit profound humoral immunodeficiency and lack mature B cell subsets, mirroring deficiency of the cytokine B cell-activating factor (BAFF). Accumulation of Sppl2a-deficient B cells was rescued by overexpression of the BAFF-induced survival protein B cell lymphoma 2 (BCL2) but not BAFF and was distinguished by low surface BAFF receptor and IgM and IgD B cell receptors. CD8-negative dendritic cells were also greatly decreased. SPPL2A deficiency blocked the proteolytic processing of CD74 MHC II invariant chain in both cell types, causing dramatic build-up of the p8 product of Cathepsin S and interfering with earlier steps in CD74 endosomal retention and processing. The findings illuminate an important role for the final step in the CD74-MHC II pathway and a new target for protease inhibitor treatment of B cell diseases.

FoxO transcription factors, inhibited by insulin/insulin-like growth factor signalling (IIS), are crucial players in numerous organismal processes including lifespan. Using genomic tools, we uncover over 700 direct dFOXO targets in adult female Drosophila. dFOXO is directly required for transcription of several IIS components and interacting pathways, such as TOR, in the wild-type fly. The genomic locations occupied by dFOXO in adults are different from those observed in larvae or cultured cells. These locations remain unchanged upon activation by stresses or reduced IIS, but the binding is increased and additional targets activated upon genetic reduction in IIS. We identify the part of the IIS transcriptional response directly controlled by dFOXO and the indirect effects and show that parts of the transcriptional response to IIS reduction do not require dfoxo. Promoter analyses revealed GATA and other forkhead factors as candidate mediators of the indirect and dfoxo-independent effects. We demonstrate genome-wide evolutionary conservation of dFOXO targets between the fly and the worm Caenorhabditis elegans, enriched for a second tier of regulators including the dHR96/daf-12 nuclear hormone receptor.

Copy number variants (CNVs) account for a major proportion of human genetic polymorphism and have been predicted to have an important role in genetic susceptibility to common disease. To address this we undertook a large, direct genome-wide study of association between CNVs and eight common human diseases. Using a purpose-designed array we typed approximately 19,000 individuals into distinct copy-number classes at 3,432 polymorphic CNVs, including an estimated approximately 50% of all common CNVs larger than 500 base pairs. We identified several biological artefacts that lead to false-positive associations, including systematic CNV differences between DNAs derived from blood and cell lines. Association testing and follow-up replication analyses confirmed three loci where CNVs were associated with disease-IRGM for Crohn's disease, HLA for Crohn's disease, rheumatoid arthritis and type 1 diabetes, and TSPAN8 for type 2 diabetes-although in each case the locus had previously been identified in single nucleotide polymorphism (SNP)-based studies, reflecting our observation that most common CNVs that are well-typed on our array are well tagged by SNPs and so have been indirectly explored through SNP studies. We conclude that common CNVs that can be typed on existing platforms are unlikely to contribute greatly to the genetic basis of common human diseases.

Structural variations of DNA greater than 1 kilobase in size account for most bases that vary among human genomes, but are still relatively under-ascertained. Here we use tiling oligonucleotide microarrays, comprising 42 million probes, to generate a comprehensive map of 11,700 copy number variations (CNVs) greater than 443 base pairs, of which most (8,599) have been validated independently. For 4,978 of these CNVs, we generated reference genotypes from 450 individuals of European, African or East Asian ancestry. The predominant mutational mechanisms differ among CNV size classes. Retrotransposition has duplicated and inserted some coding and non-coding DNA segments randomly around the genome. Furthermore, by correlation with known trait-associated single nucleotide polymorphisms (SNPs), we identified 30 loci with CNVs that are candidates for influencing disease susceptibility. Despite this, having assessed the completeness of our map and the patterns of linkage disequilibrium between CNVs and SNPs, we conclude that, for complex traits, the heritability void left by genome-wide association studies will not be accounted for by common CNVs.

Multicellular animals match costly activities, such as growth and reproduction, to the environment through nutrient-sensing pathways. The insulin/IGF signaling (IIS) pathway plays key roles in growth, metabolism, stress resistance, reproduction, and longevity in diverse organisms including mammals. Invertebrate genomes often contain multiple genes encoding insulin-like ligands, including seven Drosophila insulin-like peptides (DILPs). We investigated the evolution, diversification, redundancy, and functions of the DILPs, combining evolutionary analysis, based on the completed genome sequences of 12 Drosophila species, and functional analysis, based on newly-generated knock-out mutations for all 7 dilp genes in D. melanogaster. Diversification of the 7 DILPs preceded diversification of Drosophila species, with stable gene diversification and family membership, suggesting stabilising selection for gene function. Gene knock-outs demonstrated both synergy and compensation of expression between different DILPs, notably with DILP3 required for normal expression of DILPs 2 and 5 in brain neurosecretory cells and expression of DILP6 in the fat body compensating for loss of brain DILPs. Loss of DILP2 increased lifespan and loss of DILP6 reduced growth, while loss of DILP7 did not affect fertility, contrary to its proposed role as a Drosophila relaxin. Importantly, loss of DILPs produced in the brain greatly extended lifespan but only in the presence of the endosymbiontic bacterium Wolbachia, demonstrating a specific interaction between IIS and Wolbachia in lifespan regulation. Furthermore, loss of brain DILPs blocked the responses of lifespan and fecundity to dietary restriction (DR) and the DR response of these mutants suggests that IIS extends lifespan through mechanisms that both overlap with those of DR and through additional mechanisms that are independent of those at work in DR. Evolutionary conservation has thus been accompanied by synergy, redundancy, and functional differentiation between DILPs, and these features may themselves be of evolutionary advantage., Competing Interests: The authors have declared that no competing interests exist.

Background: Large-scale high throughput studies using microarray technology have established that copy number variation (CNV) throughout the genome is more frequent than previously thought. Such variation is known to play an important role in the presence and development of phenotypes such as HIV-1 infection and Alzheimer's disease. However, methods for analyzing the complex data produced and identifying regions of CNV are still being refined., Results: We describe the presence of a genome-wide technical artifact, spatial autocorrelation or 'wave', which occurs in a large dataset used to determine the location of CNV across the genome. By removing this artifact we are able to obtain both a more biologically meaningful clustering of the data and an increase in the number of CNVs identified by current calling methods without a major increase in the number of false positives detected. Moreover, removing this artifact is critical for the development of a novel model-based CNV calling algorithm - CNVmix - that uses cross-sample information to identify regions of the genome where CNVs occur. For regions of CNV that are identified by both CNVmix and current methods, we demonstrate that CNVmix is better able to categorize samples into groups that represent copy number gains or losses., Conclusion: Removing artifactual 'waves' (which appear to be a general feature of array comparative genomic hybridization (aCGH) datasets) and using cross-sample information when identifying CNVs enables more biological information to be extracted from aCGH experiments designed to investigate copy number variation in normal individuals.

Copy number variation (CNV) of DNA sequences is functionally significant but has yet to be fully ascertained. We have constructed a first-generation CNV map of the human genome through the study of 270 individuals from four populations with ancestry in Europe, Africa or Asia (the HapMap collection). DNA from these individuals was screened for CNV using two complementary technologies: single-nucleotide polymorphism (SNP) genotyping arrays, and clone-based comparative genomic hybridization. A total of 1,447 copy number variable regions (CNVRs), which can encompass overlapping or adjacent gains or losses, covering 360 megabases (12% of the genome) were identified in these populations. These CNVRs contained hundreds of genes, disease loci, functional elements and segmental duplications. Notably, the CNVRs encompassed more nucleotide content per genome than SNPs, underscoring the importance of CNV in genetic diversity and evolution. The data obtained delineate linkage disequilibrium patterns for many CNVs, and reveal marked variation in copy number among populations. We also demonstrate the utility of this resource for genetic disease studies.

Recent work has shown that copy number polymorphism is an important class of genetic variation in human genomes. Here we report a new method that uses SNP genotype data from parent-offspring trios to identify polymorphic deletions. We applied this method to data from the International HapMap Project to produce the first high-resolution population surveys of deletion polymorphism. Approximately 100 of these deletions have been experimentally validated using comparative genome hybridization on tiling-resolution oligonucleotide microarrays. Our analysis identifies a total of 586 distinct regions that harbor deletion polymorphisms in one or more of the families. Notably, we estimate that typical individuals are hemizygous for roughly 30-50 deletions larger than 5 kb, totaling around 550-750 kb of euchromatic sequence across their genomes. The detected deletions span a total of 267 known and predicted genes. Overall, however, the deleted regions are relatively gene-poor, consistent with the action of purifying selection against deletions. Deletion polymorphisms may well have an important role in the genetics of complex traits; however, they are not directly observed in most current gene mapping studies. Our new method will permit the identification of deletion polymorphisms in high-density SNP surveys of trio or other family data.

The human X chromosome has a unique biology that was shaped by its evolution as the sex chromosome shared by males and females. We have determined 99.3% of the euchromatic sequence of the X chromosome. Our analysis illustrates the autosomal origin of the mammalian sex chromosomes, the stepwise process that led to the progressive loss of recombination between X and Y, and the extent of subsequent degradation of the Y chromosome. LINE1 repeat elements cover one-third of the X chromosome, with a distribution that is consistent with their proposed role as way stations in the process of X-chromosome inactivation. We found 1,098 genes in the sequence, of which 99 encode proteins expressed in testis and in various tumour types. A disproportionately high number of mendelian diseases are documented for the X chromosome. Of this number, 168 have been explained by mutations in 113 X-linked genes, which in many cases were characterized with the aid of the DNA sequence.

As more genomes are sequenced, there is an increasing need for automated first-pass annotation which allows timely access to important genomic information. The Ensembl gene-building system enables fast automated annotation of eukaryotic genomes. It annotates genes based on evidence derived from known protein, cDNA, and EST sequences. The gene-building system rests on top of the core Ensembl (MySQL) database schema and Perl Application Programming Interface (API), and the data generated are accessible through the Ensembl genome browser (http://www.ensembl.org). To date, the Ensembl predicted gene sets are available for the A. gambiae, C. briggsae, zebrafish, mouse, rat, and human genomes and have been heavily relied upon in the publication of the human, mouse, rat, and A. gambiae genome sequence analysis. Here we describe in detail the gene-building system and the algorithms involved. All code and data are freely available from http://www.ensembl.org.

Ensembl (http://www.ensembl.org/) is a bioinformatics project to organize biological information around the sequences of large genomes. It is a comprehensive source of stable automatic annotation of individual genomes, and of the synteny and orthology relationships between them. It is also a framework for integration of any biological data that can be mapped onto features derived from the genomic sequence. Ensembl is available as an interactive Web site, a set of flat files, and as a complete, portable open source software system for handling genomes. All data are provided without restriction, and code is freely available. Ensembl's aims are to continue to "widen" this biological integration to include other model organisms relevant to understanding human biology as they become available; to "deepen" this integration to provide an ever more seamless linkage between equivalent components in different species; and to provide further classification of functional elements in the genome that have been previously elusive.

The PilE protein is the major component of the Neisseria meningitidis pilus, which is encoded by the pilE/pilS locus that includes an expressed gene and eight homologous silent fragments. The silent gene fragments have been shown to recombine through gene conversion with the expressed gene and thereby provide a means by which novel antigenic variants of the PilE protein can be generated. We have analyzed the evolutionary rate of the pilE gene using the nucleotide sequence of two complete pilE/pilS loci. The very high rate of evolution displayed by the PilE protein appears driven by both recombination and positive selection. Within the semivariable region of the pilE and pilS genes, recombination appears to occur within multiple small sequence blocks that lie between conserved sequence elements. Within the hypervariable region, positive selection was identified from comparison of the silent and expressed genes. The unusual gene conversion mechanism that operates at the pilE/pilS locus is a strategy employed by N. meningitidis to enhance mutation of certain regions of the PilE protein. The silent copies of the gene effectively allow "parallelized" evolution of pilE, thus enabling the encoded protein to rapidly explore a large area of sequence space in an effort to find novel antigenic variants.