Search

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., Genomes vary from one another in multifarious ways, and the totality of this genetic variation underpins the heritability of human traits. Over the past two years, the human reference sequence [...]

Author(s): Richard Redon [1]; Shumpei Ishikawa [2, 3]; Karen R. Fitch [4]; Lars Feuk [5, 6]; George H. Perry [7]; T. Daniel Andrews [1]; Heike Fiegler [1]; Michael H. Shapero [...]

Author(s): Mark T. Ross (corresponding author) [1]; Darren V. Grafham [1]; Alison J. Coffey [1]; Steven Scherer [2]; Kirsten McLay [1]; Donna Muzny [2]; Matthias Platzer [3]; Gareth R. Howell [...]

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.

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Objective Systemic lupus erythematosus (SLE) is a 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-α (IFNα) 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 neutralizing anti-IFNα 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. [ABSTRACT FROM AUTHOR]

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 ∼19,000 individuals into distinct copy-number classes at 3,432 polymorphic CNVs, including an estimated ∼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. [ABSTRACT FROM AUTHOR]

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 insulinlike 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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Pyroptosis requires the induction of gasdermin D expression by the transcription factor IRF2. IRF2 induces gasdermin D: In response to activation of canonical and noncanonical inflammasomes, a subset of caspases processes the protein gasdermin D (GSDMD) to release N-terminal fragments that oligomerize and form pores in the plasma membrane. Assembly of the GSDMD pore leads to release of the inflammatory cytokine IL-1β and causes cell death by pyroptosis. Kayagaki et al. found that loss of the transcriptional regulator IRF2 reduced GSDMD mRNA and protein abundance in mice and in human cells, resulting in decreased IL-1β secretion and reduced pyroptosis in response to inflammasome activation. Given that loss of GSDMD in mice results in ameliorated disease in models of inflammasome-driven pathologies, these findings suggest that IRF2 might be a therapeutic target for the treatment of sepsis and other inflammasome-mediated diseases. 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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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.)

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.

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.)

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.

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.

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.

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.