Search

We aimed to study the histopathological and immunohistochemistry features in clinically diagnosed cases of nanophthalmos using light microscopy. This was an observational comparative study. We enrolled four eyes of four consecutive patients with nanophthalmos and visually significant cataract, who underwent cataract surgery with prophylactic posterior sclerostomy. Histological analysis of the excised scleral tissue was done and compared with age-matched cadaver controls between January 2021 and October 2021. Hematoxylin and Eosin (H&E) stains were used for histological analysis, and was further supplemented with immunohistochemistry (IHC) and immunofluorescence (IF) analyses using a simple light microscope. The immunostained sections were analyzed using confocal microscope for the fibronectin expression level. The main outcome measure was demonstration of histological changes of sclera in nanophthalmic eyes undergoing cataract surgery. Light microscopic features of nanophthalmos revealed thick fibers with fraying and lightly stained cores, irregular serrated edges, and randomly interspersed fibroblasts compared to regular arrangement of collagen fibers seen in cadaver controls. Immunohistochemistry analysis with anti-fibronectin antibody showed strong positivity in clustered fibers in nanophthalmos, and less intense diffuse staining in cadaver tissue. Histoclinical correlation was observed in one nanophthalmic scleral tissue with axial length less than 17 mm showing severe disorganization with diffuse collagenization, loss of fibrillary architecture compared to another specimen with axial length more than 17 mm. Simple, cost-effective light microscopy using basic stains was effective in identifying the characteristic histopathological features in nanophthalmic eyes, and this was further highlighted by immunohistochemistry and immunofluorescence analyses.

We aimed to study the histopathological and immunohistochemistry features in clinically diagnosed cases of nanophthalmos using light microscopy. This was an observational comparative study. We enrolled four eyes of four consecutive patients with nanophthalmos and visually significant cataract, who underwent cataract surgery with prophylactic posterior sclerostomy. Histological analysis of the excised scleral tissue was done and compared with age-matched cadaver controls between January 2021 and October 2021. Hematoxylin and Eosin (H&E) stains were used for histological analysis, and was further supplemented with immunohistochemistry (IHC) and immunofluorescence (IF) analyses using a simple light microscope. The immunostained sections were analyzed using confocal microscope for the fibronectin expression level. The main outcome measure was demonstration of histological changes of sclera in nanophthalmic eyes undergoing cataract surgery. Light microscopic features of nanophthalmos revealed thick fibers with fraying and lightly stained cores, irregular serrated edges, and randomly interspersed fibroblasts compared to regular arrangement of collagen fibers seen in cadaver controls. Immunohistochemistry analysis with anti-fibronectin antibody showed strong positivity in clustered fibers in nanophthalmos, and less intense diffuse staining in cadaver tissue. Histoclinical correlation was observed in one nanophthalmic scleral tissue with axial length less than 17 mm showing severe disorganization with diffuse collagenization, loss of fibrillary architecture compared to another specimen with axial length more than 17 mm. Simple, cost-effective light microscopy using basic stains was effective in identifying the characteristic histopathological features in nanophthalmic eyes, and this was further highlighted by immunohistochemistry and immunofluorescence analyses. [ABSTRACT FROM AUTHOR]

Abstract Background Genome-wide transcriptional profiling of patient blood samples offers a powerful tool to investigate underlying disease mechanisms and personalized treatment decisions. Most studies are based on analysis of total peripheral blood mononuclear cells (PBMCs), a mixed population. In this case, accuracy is inherently limited since cell subset-specific differential expression of gene signatures will be diluted by RNA from other cells. While using specific PBMC subsets for transcriptional profiling would improve our ability to extract knowledge from these data, it is rarely obvious which cell subset(s) will be the most informative. Results We have developed a computational method (Subset Prediction from Enrichment Correlation, SPEC) to predict the cellular source for a pre-defined list of genes (i.e. a gene signature) using only data from total PBMCs. SPEC does not rely on the occurrence of cell subset-specific genes in the signature, but rather takes advantage of correlations with subset-specific genes across a set of samples. Validation using multiple experimental datasets demonstrates that SPEC can accurately identify the source of a gene signature as myeloid or lymphoid, as well as differentiate between B cells, T cells, NK cells and monocytes. Using SPEC, we predict that myeloid cells are the source of the interferon-therapy response gene signature associated with HCV patients who are non-responsive to standard therapy. Conclusions SPEC is a powerful technique for blood genomic studies. It can help identify specific cell subsets that are important for understanding disease and therapy response. SPEC is widely applicable since only gene expression profiles from total PBMCs are required, and thus it can easily be used to mine the massive amount of existing microarray or RNA-seq data.

Phosphorylation of proteins on tyrosine (Tyr) residues evolved in metazoan organisms as a mechanism of coordinating tissue growth 1 . Multicellular eukaryotes typically have more than 50 distinct protein Tyr kinases that catalyse the phosphorylation of thousands of Tyr residues throughout the proteome 1-3 . How a given Tyr kinase can phosphorylate a specific subset of proteins at unique Tyr sites is only partially understood 4-7 . Here we used combinatorial peptide arrays to profile the substrate sequence specificity of all human Tyr kinases. Globally, the Tyr kinases demonstrate considerable diversity in optimal patterns of residues surrounding the site of phosphorylation, revealing the functional organization of the human Tyr kinome by substrate motif preference. Using this information, Tyr kinases that are most compatible with phosphorylating any Tyr site can be identified. Analysis of mass spectrometry phosphoproteomic datasets using this compendium of kinase specificities accurately identifies specific Tyr kinases that are dysregulated in cells after stimulation with growth factors, treatment with anti-cancer drugs or expression of oncogenic variants. Furthermore, the topology of known Tyr signalling networks naturally emerged from a comparison of the sequence specificities of the Tyr kinases and the SH2 phosphotyrosine (pTyr)-binding domains. Finally we show that the intrinsic substrate specificity of Tyr kinases has remained fundamentally unchanged from worms to humans, suggesting that the fidelity between Tyr kinases and their protein substrate sequences has been maintained across hundreds of millions of years of evolution., (© 2024. The Author(s).)

The molecular mechanisms of the liver metastasis of colorectal cancer (CRLM) remain poorly understood. Here, we applied machine learning and bioinformatics trajectory inference to analyze a gene expression dataset of CRLM. We studied the co-regulation patterns at the gene level, the potential paths of tumor development, their functional context, and their prognostic relevance. Our analysis confirmed the subtyping of five liver metastasis subtypes (LMS). We provide gene-marker signatures for each LMS, and a comprehensive functional characterization that considers both the hallmarks of cancer and the tumor microenvironment. The ordering of CRLMs along a pseudotime-tree revealed a continuous shift in expression programs, suggesting a developmental relationship between the subtypes. Notably, trajectory inference and personalized analysis discovered a range of epigenetic states that shape and guide metastasis progression. By constructing prognostic maps that divided the expression landscape into regions associated with favorable and unfavorable prognoses, we derived a prognostic expression score. This was associated with critical processes such as epithelial-mesenchymal transition, treatment resistance, and immune evasion. These factors were associated with responses to neoadjuvant treatment and the formation of an immuno-suppressive, mesenchymal state. Our machine learning-based molecular profiling provides an in-depth characterization of CRLM heterogeneity with possible implications for treatment and personalized diagnostics.

Purpose: To study the refractive profile of children after they received intravitreal injection of bevacizumab for retinopathy of prematurity (ROP)., Methods: The study was conducted at a tertiary eye care hospital in South India. ROP patients of more than 1 year of age, presenting to the Pediatric Ophthalmology Clinic and Retina Clinic and having history of treatment for type I ROP with intravitreal bevacizumab (IVB) or intravitreal bevacizumab and laser photocoagulation were included in the study. Cycloplegic refraction was done, and the refractive status was evaluated. The refractive status of age-matched, full-term children with uneventful perinatal and neonatal history was also recorded and compared to the study group., Results: Among 134 eyes of 67 study subjects, the major refractive error was myopia in 93 eyes (69.4%; spherical equivalent [SE] = -2.89 ± 3.1, range = -11.5 to -0.5 D). There were 75 eyes (56%) with low-to-moderate myopia; high myopia was seen in 13.4%, emmetropia in 18.7%, and hypermetropia in 11.9% of eyes. The majority of them (87%) had with-the-rule (WTR) astigmatism. In 134 eyes, the SE was -1.78 ± 3.2 (range = -11.5 to 4 D); the SE of the 75 eyes with low-to-moderate myopia was -1.53 ± 1.2 (range = -0.50 to -5 D). In the control group, the majority had emmetropia (91.8%). There was no significant association between the age at which IVB had been injected and the development of refractive errors (P = 0.078). The prevalence of low-to-moderate myopia was more than high myopia in patients with zone I and zone II ROP before treatment (60.0% and 54.5%, respectively)., Conclusion: Myopia was the major refractive error seen in post-IVB pediatric patients. WTR astigmatism was more commonly seen. The age at which IVB injection had been given had no effect on the development of refractive errors., Competing Interests: None

Transformation to aggressive disease histologies generates formidable clinical challenges across cancers, but biological insights remain few. We modeled the genetic heterogeneity of chronic lymphocytic leukemia (CLL) through multiplexed in vivo CRISPR-Cas9 B-cell editing of recurrent CLL loss-of-function drivers in mice and recapitulated the process of transformation from indolent CLL into large cell lymphoma [i.e., Richter syndrome (RS)]. Evolutionary trajectories of 64 mice carrying diverse combinatorial gene assortments revealed coselection of mutations in Trp53, Mga, and Chd2 and the dual impact of clonal Mga/Chd2 mutations on E2F/MYC and interferon signaling dysregulation. Comparative human and murine RS analyses demonstrated tonic PI3K signaling as a key feature of transformed disease, with constitutive activation of the AKT and S6 kinases, downmodulation of the PTEN phosphatase, and convergent activation of MYC/PI3K transcriptional programs underlying enhanced sensitivity to MYC/mTOR/PI3K inhibition. This robust experimental system presents a unique framework to study lymphoid biology and therapy., Significance: Mouse models reflective of the genetic complexity and heterogeneity of human tumors remain few, including those able to recapitulate transformation to aggressive disease histologies. Herein, we model CLL transformation into RS through multiplexed in vivo gene editing, providing key insight into the pathophysiology and therapeutic vulnerabilities of transformed disease. This article is highlighted in the In This Issue feature, p. 101., (© 2022 American Association for Cancer Research.)

Protein phosphorylation is one of the most widespread post-translational modifications in biology 1,2 . With advances in mass-spectrometry-based phosphoproteomics, 90,000 sites of serine and threonine phosphorylation have so far been identified, and several thousand have been associated with human diseases and biological processes 3,4 . For the vast majority of phosphorylation events, it is not yet known which of the more than 300 protein serine/threonine (Ser/Thr) kinases encoded in the human genome are responsible 3 . Here we used synthetic peptide libraries to profile the substrate sequence specificity of 303 Ser/Thr kinases, comprising more than 84% of those predicted to be active in humans. Viewed in its entirety, the substrate specificity of the kinome was substantially more diverse than expected and was driven extensively by negative selectivity. We used our kinome-wide dataset to computationally annotate and identify the kinases capable of phosphorylating every reported phosphorylation site in the human Ser/Thr phosphoproteome. For the small minority of phosphosites for which the putative protein kinases involved have been previously reported, our predictions were in excellent agreement. When this approach was applied to examine the signalling response of tissues and cell lines to hormones, growth factors, targeted inhibitors and environmental or genetic perturbations, it revealed unexpected insights into pathway complexity and compensation. Overall, these studies reveal the intrinsic substrate specificity of the human Ser/Thr kinome, illuminate cellular signalling responses and provide a resource to link phosphorylation events to biological pathways., (© 2023. The Author(s).)

Chronic lymphocytic leukemia (CLL) is characterized by disordered DNA methylation, suggesting these epigenetic changes might play a critical role in disease onset and progression. The methyltransferase DNMT3A is a key regulator of DNA methylation. Although DNMT3A somatic mutations in CLL are rare, we found that low DNMT3A expression is associated with more aggressive disease. A conditional knockout mouse model showed that homozygous depletion of Dnmt3a from B cells results in the development of CLL with 100% penetrance at a median age of onset of 5.3 months, and heterozygous Dnmt3a depletion yields a disease penetrance of 89% with a median onset at 18.5 months, confirming its role as a haploinsufficient tumor suppressor. B1a cells were confirmed as the cell of origin of disease in this model, and Dnmt3a depletion resulted in focal hypomethylation and activation of Notch and Myc signaling. Amplification of chromosome 15 containing the Myc gene was detected in all CLL mice tested, and infiltration of high- Myc -expressing CLL cells in the spleen was observed. Notably, hyperactivation of Notch and Myc signaling was exclusively observed in the Dnmt3a CLL mice, but not in three other CLL mouse models tested ( Sf3b1-Atm , Ikzf3 , and MDR ), and Dnmt3a -depleted CLL were sensitive to pharmacologic inhibition of Notch signaling in vitro and in vivo . Consistent with these findings, human CLL samples with lower DNMT3A expression were more sensitive to Notch inhibition than those with higher DNMT3A expression. Altogether, these results suggest that Dnmt3a depletion induces CLL that is highly dependent on activation of Notch and Myc signaling. SIGNIFICANCE: Loss of DNMT3A expression is a driving event in CLL and is associated with aggressive disease, activation of Notch and Myc signaling, and enhanced sensitivity to Notch inhibition., (©2021 American Association for Cancer Research.)

Purpose: Whole-exome (WES) and RNA sequencing (RNA-seq) are key components of cancer immunogenomic analyses. To evaluate the consistency of tumor WES and RNA-seq profiling platforms across different centers, the Cancer Immune Monitoring and Analysis Centers (CIMAC) and the Cancer Immunologic Data Commons (CIDC) conducted a systematic harmonization study., Experimental Design: DNA and RNA were centrally extracted from fresh frozen and formalin-fixed paraffin-embedded non-small cell lung carcinoma tumors and distributed to three centers for WES and RNA-seq profiling. In addition, two 10-plex HapMap cell line pools with known mutations were used to evaluate the accuracy of the WES platforms., Results: The WES platforms achieved high precision (> 0.98) and recall (> 0.87) on the HapMap pools when evaluated on loci using > 50× common coverage. Nonsynonymous mutations clustered by tumor sample, achieving an index of specific agreement above 0.67 among replicates, centers, and sample processing. A DV200 > 24% for RNA, as a putative presequencing RNA quality control (QC) metric, was found to be a reliable threshold for generating consistent expression readouts in RNA-seq and NanoString data. MedTIN > 30 was likewise assessed as a reliable RNA-seq QC metric, above which samples from the same tumor across replicates, centers, and sample processing runs could be robustly clustered and HLA typing, immune infiltration, and immune repertoire inference could be performed., Conclusions: The CIMAC collaborating laboratory platforms effectively generated consistent WES and RNA-seq data and enable robust cross-trial comparisons and meta-analyses of highly complex immuno-oncology biomarker data across the NCI CIMAC-CIDC Network., (©2020 American Association for Cancer Research.)

Background: Immune checkpoint inhibitors have revolutionized cancer treatment, but the benefits in refractory patients with esophageal cancer have been modest. Predictors of response as well as new targets for novel therapeutic combinations are needed. In this phase 2 clinical trial, we tested single-agent pembrolizumab in patients with advanced esophageal cancer, who received at least one prior line of therapy., Methods: Pembrolizumab 200 mg every 3 weeks was tested in 49 patients with refractory esophageal cancer: 39 with adenocarcinoma and 10 with esophageal squamous cell carcinoma. Major endpoints were radiological response by Immune-related Response Evaluation Criteria In Solid Tumors and survival. Tumor samples were evaluated for programmed cell death ligand 1 (PD-L1) expression, tumor mutational burden (TMB), and immune contexture by both NanoString mRNA expression analysis and flow cytometry. Peripheral blood mononuclear cells and a panel of circulating chemokines were also analyzed., Results: The overall response rate (ORR) was 8% (4 of 49 patients; 95% CI 2.3% to 19.6%). Median overall survival (OS) was 5.8 months (95% CI 4.0 to 9.5). ORR and OS were not associated with histology. For PD-L1-positive patients, ORR was 13.3% (95% CI 1.7% to 40.5%) and median OS was 7.9 months (95% CI 4.7 to 15.5). A trend toward improved OS was observed in seven patients with a TMB ≥10 mut/Mb (p=0.086). Tumors with a PD-L1 Combined Positive Score ≥1 showed enrichment of LAG3 (p=0.005) and IDO1 (p=0.04) gene expression. Baseline levels of circulating CXCL10, interleukin 2 (IL2) receptor α (IL2RA) and IL6 were associated with survival: CXCL10 favorably, (HR 0.37, p=0.002 (progression-free survival); HR 0.55, p=0.018 (OS)); IL2RA and IL6 unfavorably (HR 1.57, p=0.020 for IL6 (OS); HR 2.36, p=0.025 for IL2RA (OS))., Conclusions: Pembrolizumab monotherapy was modestly effective in refractory esophageal cancer. Circulating CXCL10 at baseline appeared to be a robust predictor of response. Other T cell exhaustion markers are upregulated in PD-L1-positive patients, suggesting that immunotherapy combinations such as anti-LAG3/programmed cell death protein 1 (PD-1) or anti-IDO1/PD-1 may be of promise in refractory esophageal cancer., Competing Interests: Competing interests: AJB received funding support from Bayer, Merck and Novartis. JMC received research funding to his institution from Abbvie, Merus, Roche, and Bristol Myers Squibb, received research funding from Merck, Astrazeneca, Esperas Pharma, and Tesaro, received consulting fees from Bristol Myers Squibb, and received travel funding from Bristol Myers Squibb. PCE has consulted for and has received honoraria from ALX Oncology, Arcus Bioscience, Astellas, Astra-Zeneca, Blueprint Medicines, Bristol-Myers Squibb, Celgene, Daiichi-Sankyo, Five Prime, Ideaya, Istari, Legend, Lilly, Loxo, Merck, Ono, Taiho, Takeda, Turning Point Therapeutics, Xencor and Zymeworks., (© Author(s) (or their employer(s)) 2021. Re-use permitted under CC BY-NC. No commercial re-use. See rights and permissions. Published by BMJ.)

Hotspot mutation of IKZF3 (IKZF3-L162R) has been identified as a putative driver of chronic lymphocytic leukemia (CLL), but its function remains unknown. Here, we demonstrate its driving role in CLL through a B cell-restricted conditional knockin mouse model. Mutant Ikzf3 alters DNA binding specificity and target selection, leading to hyperactivation of B cell receptor (BCR) signaling, overexpression of nuclear factor κB (NF-κB) target genes, and development of CLL-like disease in elderly mice with a penetrance of ~40%. Human CLL carrying either IKZF3 mutation or high IKZF3 expression was associated with overexpression of BCR/NF-κB pathway members and reduced sensitivity to BCR signaling inhibition by ibrutinib. Our results thus highlight IKZF3 oncogenic function in CLL via transcriptional dysregulation and demonstrate that this pro-survival function can be achieved by either somatic mutation or overexpression of this CLL driver. This emphasizes the need for combinatorial approaches to overcome IKZF3-mediated BCR inhibitor resistance., Competing Interests: Declaration of interests C.J.W. is an equity holder of Biontech, Inc. and receives research funding from Pharmacyclics. D.S.N. has been a consultant for H3 Biomedicine and received research funding from Celgene. J.A.B. reports receiving grant support and advisory board fees from Pharmacyclics, grant support, advisory board fees, and lecture fees from Gilead, advisory board fees from AstraZeneca, and lecture fees and travel support from Janssen. T.J.K. has received research funding and/or has served as an advisor to Ascerta/AstraZeneca, Celgene, Genentech/Roche, Gilead, Janssen, Loxo Oncology, Octernal Therapeutics, Pharmacyclics/AbbVie, TG Therapeutics, VelosBio, and Verastem. Cirmtuzumab was developed by T.J.K. and licensed by the University of California to Oncternal Therapeutics, Inc., which has provided stock/options to the university and T.J.K. All other authors do not have any relevant conflict of interest., (Copyright © 2021 Elsevier Inc. All rights reserved.)

Personal neoantigen vaccines have been envisioned as an effective approach to induce, amplify and diversify antitumor T cell responses. To define the long-term effects of such a vaccine, we evaluated the clinical outcome and circulating immune responses of eight patients with surgically resected stage IIIB/C or IVM1a/b melanoma, at a median of almost 4 years after treatment with NeoVax, a long-peptide vaccine targeting up to 20 personal neoantigens per patient ( NCT01970358 ). All patients were alive and six were without evidence of active disease. We observed long-term persistence of neoantigen-specific T cell responses following vaccination, with ex vivo detection of neoantigen-specific T cells exhibiting a memory phenotype. We also found diversification of neoantigen-specific T cell clones over time, with emergence of multiple T cell receptor clonotypes exhibiting distinct functional avidities. Furthermore, we detected evidence of tumor infiltration by neoantigen-specific T cell clones after vaccination and epitope spreading, suggesting on-target vaccine-induced tumor cell killing. Personal neoantigen peptide vaccines thus induce T cell responses that persist over years and broaden the spectrum of tumor-specific cytotoxicity in patients with melanoma.

Background: Angiosarcoma is an uncommon endothelial malignancy and a highly aggressive soft tissue sarcoma. Due to its infiltrative nature, successful management of localized angiosarcoma is often challenging. Systemic chemotherapy is used in the metastatic setting and occasionally in patients with high-risk localized disease in neoadjuvant or adjuvant settings. However, responses tend to be short-lived and most patients succumb to metastatic disease. Novel therapies are needed for patients with angiosarcomas., Methods: We performed a retrospective analysis of patients with locally advanced or metastatic angiosarcoma, who were treated with checkpoint inhibitors at our institution. We collected their clinical information and outcome measurements. In one patient with achieved complete response, we analyzed circulating and infiltrating T cells within peripheral blood and tumor tissue., Results: We have treated seven angiosarcoma (AS) patients with checkpoint inhibitors either in the context of clinical trials or off label [Pembrolizumab + Axitinib (NCT02636725; n = 1), AGEN1884, a CTLA-4 inhibitor (NCT02694822; n = 2), Pembrolizumab (n = 4)]. Five patients had cutaneous angiosarcoma, one primary breast angiosarcoma and one radiation-associated breast angiosarcoma. At 12 weeks, 5/7 patients (71%) had partial response of their lesions either on imaging and/or clinical exam and two (29%) had progressive disease. 6/7 patients are alive to date and, thus far, 3/7 patients (43%) have progressed (median 3.4 months)- one achieved partial response after pembrolizumab was switched to ongoing Nivolumab/Ipilimumab, one died of progressive disease at 31 weeks (primary breast angiosarcoma) and one was placed on pazopanib. One patient had a complete response (CR) following extended treatment with monotherapy AGEN1884. No patient experienced any ≥ grade 2 toxicities., Conclusions: This case series underscores the value of targeted immunotherapy in treating angiosarcoma. It also identifies genetic heterogeneity of cutaneous angiosarcomas and discusses specific genetic findings that may explain reported benefits from immunotherapy.

Mutation-derived neoantigens distinguish tumor from normal cells. T cells can sense the HLA-presented mutations, recognize tumor cells as non-self and destroy them. Therapeutically, immunotherapy antibodies can increase the virulence of the immune system by increasing T-cell cytotoxicity targeted toward neoantigens. Neoantigen vaccines act through antigen-presenting cells, such as dendritic cells, to activate patient-endogenous T cells that recognize vaccine-encoded mutations. Infusion of mutation-targeting T cells by adoptive cell therapy (ACT) directly increases the number and frequency of cytotoxic T cells recognizing and killing tumor cells. At the same time, publicly-funded consortia have profiled tumor genomes across many indications, identifying mutations in each tumor. For example, we find basal and HER2 positive tumors contain more mutated proteins and more TP53 mutations than luminal A/B breast tumors. HPV negative tumors have more mutated proteins than HPV positive head and neck tumors and in agreement with the hypothesis that HPV activity interferes with p53 activity, only 14% of the HPV positive mutations have TP53 mutations vs. 86% of the HPV negative tumors. Lung adenocarcinomas in smokers have over four times more mutated proteins relative to those in never smokers (median 248 vs. 61, respectively). With an eye toward immunotherapy applications, we review the spectrum of mutations in multiple indications, show variations in indication sub-types, and examine intra- and inter-indication prevalence of re-occurring mutation neoantigens that could be used for warehouse vaccines and ACT.

Human memory B cells and marginal zone (MZ) B cells share common features such as the expression of CD27 and somatic mutations in their IGHV and BCL6 genes, but the relationship between them is controversial. Here, we show phenotypic progression within lymphoid tissues as MZ B cells emerge from the mature naïve B cell pool via a precursor CD27 - CD45RB MEM55+ population distant from memory cells. By imaging mass cytometry, we find that MZ B cells and memory B cells occupy different microanatomical niches in organised gut lymphoid tissues. Both populations disseminate widely between distant lymphoid tissues and blood, and both diversify their IGHV repertoire in gut germinal centres (GC), but nevertheless remain largely clonally separate. MZ B cells are therefore not developmentally contiguous with or analogous to classical memory B cells despite their shared ability to transit through GC, where somatic mutations are acquired.

Synchronous endometrial and ovarian malignancies occur in 5% of women presenting with endometrial cancer and 10% of patients presenting with ovarian malignancy. When a high-grade serous carcinoma concurrently involves both ovary and endometrium, pathological determination of whether they are synchronous primaries or metastatic tumors from one primary site can be challenging. MicroRNAs (miRNA) are 22-nucleotide noncoding RNAs that are aberrantly expressed in cancer cells and may inherit their cellular lineage characteristics. We explored possible differential miRNA signatures that may separate high-grade ovarian serous carcinoma from primary endometrial serous carcinoma. Forty-seven samples of histologically pure high-grade serous carcinoma of both uterine (16 case) and ovarian primaries (31 cases) were included. Expression of 384 mature miRNAs was analyzed using ABI TaqMan Low-Density Arrays technology. A random forest model was used to identify miRNAs that together could differentiate between uterine and ovarian serous carcinomas. Among 150 miRNAs detectable at various levels in the study cases, a panel of 11-miRNA signatures was identified to significantly discriminate between ovarian and uterine serous carcinoma (P < .05). A nested cross-validated convergent forest plot using 6 of the 11 miRNA signature was eventually established to classify the tumors with 91.5% accuracy. In conclusion, we have characterized a miRNA signature panel in this exploratory study that shows significant discriminatory power in separating primary ovarian high-grade serous carcinoma from its endometrial counterpart., (Copyright © 2018 Elsevier Inc. All rights reserved.)

Background: The genes that produce antibodies and the immune receptors expressed on lymphocytes are not germline encoded; rather, they are somatically generated in each developing lymphocyte by a process called V(D)J recombination, which assembles specific, independent gene segments into mature composite genes. The full set of composite genes in an individual at a single point in time is referred to as the immune repertoire. V(D)J recombination is the distinguishing feature of adaptive immunity and enables effective immune responses against an essentially infinite array of antigens. Characterization of immune repertoires is critical in both basic research and clinical contexts. Recent technological advances in repertoire profiling via high-throughput sequencing have resulted in an explosion of research activity in the field. This has been accompanied by a proliferation of software tools for analysis of repertoire sequencing data. Despite the widespread use of immune repertoire profiling and analysis software, there is currently no standardized format for output files from V(D)J analysis. Researchers utilize software such as IgBLAST and IMGT/High V-QUEST to perform V(D)J analysis and infer the structure of germline rearrangements. However, each of these software tools produces results in a different file format, and can annotate the same result using different labels. These differences make it challenging for users to perform additional downstream analyses., Results: To help address this problem, we propose a standardized file format for representing V(D)J analysis results. The proposed format, VDJML, provides a common standardized format for different V(D)J analysis applications to facilitate downstream processing of the results in an application-agnostic manner. The VDJML file format specification is accompanied by a support library, written in C++ and Python, for reading and writing the VDJML file format., Conclusions: The VDJML suite will allow users to streamline their V(D)J analysis and facilitate the sharing of scientific knowledge within the community. The VDJML suite and documentation are available from https://vdjserver.org/vdjml/ . We welcome participation from the community in developing the file format standard, as well as code contributions.

Unlabelled: An effective preventive vaccine is highly sought after in order to stem the current HIV-1 pandemic. Both conservation of contiguous gp41 membrane-proximal external region (MPER) amino acid sequences across HIV-1 clades and the ability of anti-MPER broadly neutralizing antibodies (BNAbs) to block viral hemifusion/fusion establish the MPER as a prime vaccination target. In earlier studies, we described the development of an MPER vaccine formulation that takes advantage of liposomes to array the MPER on a lipid bilayer surface, paralleling its native configuration on the virus membrane while also incorporating molecular adjuvant and CD4 T cell epitope cargo. Here we demonstrate that several immunizations with MPER/liposomes induce high levels of bone marrow long-lived plasma cell (LLPC) antibody production. Single-cell immunoglobulin gene retrieval analysis shows that these plasma cells are derived from a germ line repertoire of B cells with a diverse representation of immunoglobulin genes, exhibiting antigen-driven positive selection. Characterization of LLPC recombinant monoclonal antibodies (rMAbs) indicates that antigen recognition is achieved through convergence on a common epitopic focus by utilizing various complementarity-determining region H3 (CDRH3) lengths. Importantly, the vast majority of rMAbs produced from these cells lack polyreactivity yet manifest antigen specificity in the context of lipids, shaping MPER-specific paratopes through selective pressure. Taken together, these findings demonstrate that the MPER is a vaccine target with minimal risk of generating off-target autoimmunity., Importance: A useful vaccine must generate desired long-term, antigen-specific antibody responses devoid of polyreactivity or autoreactivity. The common polyreactive features of some HIV-1 BNAbs have raised concern about elicitation of anti-MPER antibodies. Utilizing single-LLPC repertoire analysis and biophysical characterization of anti-MPER rMAbs, we show that their fine specificities require a structural fitness of the antibody combining site involving heavy and light chain variable domains shaped by somatic hypermutation and affinity maturation of B cells in the germinal center. Perhaps more importantly, our results demonstrate that the majority of MPER-specific antibodies are not inherently polyspecific and/or autoreactive, suggesting that polyreactivity of MPER-specific antibodies is separable from their antigen specificity., (Copyright © 2016 Donius et al.)

Type 1 diabetes mellitus is caused by the killing of insulin-producing β cells by CD8+T cells. The disease progression, which is chronic, does not follow a course like responses to conventional antigens such as viruses, but accelerates as glucose tolerance deteriorates. To identify the unique features of the autoimmune effectors that may explain this behavior, we analyzed diabetogenic CD8+ T cells that recognize a peptide from the diabetes antigen IGRP (NRP-V7-reactive) in prediabetic NOD mice and compared them to others that shared their phenotype (CD44(+)CD62L(lo)PD-1(+)CXCR3(+)) but negative for diabetes antigen tetramers and to LCMV (lymphocytic choriomeningitis)-reactive CD8+ T cells. There was an increase in the frequency of the NRP-V7-reactive cells coinciding with the time of glucose intolerance. The T cells persisted in hyperglycemic NOD mice maintained with an insulin pellet despite destruction of β cells. We compared gene expression in the three groups of cells compared with the other two subsets of cells, and the NRP-V7-reactive cells exhibited gene expression of memory precursor effector cells. They had reduced cellular proliferation and were less dependent on oxidative phosphorylation. When prediabetic NOD mice were treated with 2-deoxyglucose to block aerobic glycolysis, there was a reduction in the diabetes antigen versus other cells of similar phenotype and loss of lymphoid cells infiltrating the islets. In addition, treatment of NOD mice with 2-deoxyglucose resulted in improved β cell granularity. These findings identify a link between metabolic disturbances and autoreactive T cells that promotes development of autoimmune diabetes., (© 2016 by The American Society for Biochemistry and Molecular Biology, Inc.)

Unlabelled: Advances in high-throughput sequencing technologies now allow for large-scale characterization of B cell immunoglobulin (Ig) repertoires. The high germline and somatic diversity of the Ig repertoire presents challenges for biologically meaningful analysis, which requires specialized computational methods. We have developed a suite of utilities, Change-O, which provides tools for advanced analyses of large-scale Ig repertoire sequencing data. Change-O includes tools for determining the complete set of Ig variable region gene segment alleles carried by an individual (including novel alleles), partitioning of Ig sequences into clonal populations, creating lineage trees, inferring somatic hypermutation targeting models, measuring repertoire diversity, quantifying selection pressure, and calculating sequence chemical properties. All Change-O tools utilize a common data format, which enables the seamless integration of multiple analyses into a single workflow., Availability and Implementation: Change-O is freely available for non-commercial use and may be downloaded from http://clip.med.yale.edu/changeo., Contact: steven.kleinstein@yale.edu., (© The Author 2015. Published by Oxford University Press. All rights reserved. For Permissions, please e-mail: journals.permissions@oup.com.)

Individual variation in germline and expressed B-cell immunoglobulin (Ig) repertoires has been associated with aging, disease susceptibility, and differential response to infection and vaccination. Repertoire properties can now be studied at large-scale through next-generation sequencing of rearranged Ig genes. Accurate analysis of these repertoire-sequencing (Rep-Seq) data requires identifying the germline variable (V), diversity (D), and joining (J) gene segments used by each Ig sequence. Current V(D)J assignment methods work by aligning sequences to a database of known germline V(D)J segment alleles. However, existing databases are likely to be incomplete and novel polymorphisms are hard to differentiate from the frequent occurrence of somatic hypermutations in Ig sequences. Here we develop a Tool for Ig Genotype Elucidation via Rep-Seq (TIgGER). TIgGER analyzes mutation patterns in Rep-Seq data to identify novel V segment alleles, and also constructs a personalized germline database containing the specific set of alleles carried by a subject. This information is then used to improve the initial V segment assignments from existing tools, like IMGT/HighV-QUEST. The application of TIgGER to Rep-Seq data from seven subjects identified 11 novel V segment alleles, including at least one in every subject examined. These novel alleles constituted 13% of the total number of unique alleles in these subjects, and impacted 3% of V(D)J segment assignments. These results reinforce the highly polymorphic nature of human Ig V genes, and suggest that many novel alleles remain to be discovered. The integration of TIgGER into Rep-Seq processing pipelines will increase the accuracy of V segment assignments, thus improving B-cell repertoire analyses.

Background: Previous studies of immunoglobulin gene sequences in patients with allergic diseases using low-throughput Sanger sequencing have limited the analytic depth for characterization of IgE repertoires., Objectives: We used a high-throughput, next-generation sequencing approach to characterize immunoglobulin heavy-chain gene (IGH) repertoires in patients with seasonal allergic rhinitis (AR) with the aim of better understanding the underlying disease mechanisms., Methods: IGH sequences in matched peripheral blood and nasal biopsy specimens from nonallergic healthy control subjects (n = 3) and patients with grass pollen-related AR taken in season (n = 3) or out of season (n = 4) were amplified and pyrosequenced on the 454 GS FLX+ System., Results: A total of 97,610 IGH (including 8,135 IgE) sequences were analyzed. Use of immunoglobulin heavy-chain variable region gene families 1 (IGHV1) and 5 (IGHV5) was higher in IgE clonotypic repertoires compared with other antibody classes independent of atopic status. IgE repertoires measured inside the grass pollen season were more diverse and more mutated (particularly in the biopsy specimens) and had more evidence of antigen-driven selection compared with those taken outside of the pollen season or from healthy control subjects. Clonal relatedness was observed for IgE between the blood and nasal biopsy specimens. Furthermore in patients with AR, but not healthy control subjects, we found clonal relatedness between IgE and IgG classes., Conclusion: This is the first report that exploits next-generation sequencing to determine local and peripheral blood IGH repertoires in patients with respiratory allergic disease. We demonstrate that natural pollen exposure was associated with changes in IgE repertoires that were suggestive of ongoing germinal center reactions. Furthermore, these changes were more often apparent in nasal biopsy specimens compared with peripheral blood and in patients with AR compared with healthy control subjects., (Copyright © 2014 The Authors. Published by Elsevier Inc. All rights reserved.)

Unlabelled: Driven by dramatic technological improvements, large-scale characterization of lymphocyte receptor repertoires via high-throughput sequencing is now feasible. Although promising, the high germline and somatic diversity, especially of B-cell immunoglobulin repertoires, presents challenges for analysis requiring the development of specialized computational pipelines. We developed the REpertoire Sequencing TOolkit (pRESTO) for processing reads from high-throughput lymphocyte receptor studies. pRESTO processes raw sequences to produce error-corrected, sorted and annotated sequence sets, along with a wealth of metrics at each step. The toolkit supports multiplexed primer pools, single- or paired-end reads and emerging technologies that use single-molecule identifiers. pRESTO has been tested on data generated from Roche and Illumina platforms. It has a built-in capacity to parallelize the work between available processors and is able to efficiently process millions of sequences generated by typical high-throughput projects., Availability and Implementation: pRESTO is freely available for academic use. The software package and detailed tutorials may be downloaded from http://clip.med.yale.edu/presto., (© The Author 2014. Published by Oxford University Press. All rights reserved. For Permissions, please e-mail: journals.permissions@oup.com.)

Pemphigus vulgaris (PV) is a potentially fatal blistering disease caused by autoantibodies (autoAbs) against desmoglein 3 (Dsg3). Here, we clone anti-Dsg3 antibodies (Abs) from four PV patients and identify pathogenic VH1-46 autoAbs from all four patients. Unexpectedly, VH1-46 autoAbs had relatively few replacement mutations. We reverted antibody somatic mutations to their germline sequences to determine the requirement of mutations for autoreactivity. Three of five VH1-46 germline-reverted Abs maintain Dsg3 binding, compared with zero of five non-VH1-46 germline-reverted Abs. Site-directed mutagenesis of VH1-46 Abs demonstrates that acidic amino-acid residues introduced by somatic mutation or heavy chain VDJ recombination are necessary and sufficient for Dsg3 binding. Our data suggest that VH1-46 autoantibody gene usage is commonly found in PV because VH1-46 Abs require few to no mutations to acquire Dsg3 autoreactivity, which may favour their early selection. Common VH gene usage indicates common humoral immune responses, even among unrelated patients.

Objective: To explore the association of a functional germline variant in the 3'-UTR of KRAS with endometrial cancer risk, as well as the association of microRNA (miRNA) signatures and the KRAS-variant with clinical characteristics and survival outcomes in two prospective RTOG endometrial cancer trials., Methods/materials: The association of the KRAS-variant with endometrial cancer risk was evaluated by case-control analysis of 467 women with type 1 or 2 endometrial cancer and 582 age-matched controls. miRNA and DNA were isolated for expression profiling and genotyping from tumor specimens of 46 women with type 1 endometrial cancer enrolled in RTOG trials 9708 and 9905. miRNA expression levels and KRAS-variant genotype were correlated with patient and tumor characteristics, and survival outcomes were evaluated by variant allele type., Results: The KRAS-variant was not significantly associated with overall endometrial cancer risk (14% controls and 17% type 1 cancers), although was enriched in type 2 endometrial cancers (24%, p = 0.2). In the combined analysis of RTOG 9708/9905, miRNA expression differed by age, presence of lymphovascular invasion and KRAS-variant status. Overall survival rates at 3 years for patients with the variant and wild-type alleles were 100% and 77% (HR 0.3, p = 0.24), respectively, favoring the variant., Conclusions: The KRAS-variant may be a genetic marker of risk for type 2 endometrial cancers. In addition, tumor miRNA expression appears to be associated with patient age, lymphovascular invasion and the KRAS-variant, supporting the hypothesis that altered tumor biology can be measured by miRNA expression, and that the KRAS-variant likely impacts endometrial tumor biology.

The adaptive immune system confers protection by generating a diverse repertoire of antibody receptors that are rapidly expanded and contracted in response to specific targets. Next-generation DNA sequencing now provides the opportunity to survey this complex and vast repertoire. In the present work, we describe a set of tools for the analysis of antibody repertoires and their application to elucidating the dynamics of the response to viral vaccination in human volunteers. By analyzing data from 38 separate blood samples across 2 y, we found that the use of the germ-line library of V and J segments is conserved between individuals over time. Surprisingly, there appeared to be no correlation between the use level of a particular VJ combination and degree of expansion. We found the antibody RNA repertoire in each volunteer to be highly dynamic, with each individual displaying qualitatively different response dynamics. By using combinatorial phage display, we screened selected VH genes paired with their corresponding VL library for affinity against the vaccine antigens. Altogether, this work presents an additional set of tools for profiling the human antibody repertoire and demonstrates characterization of the fast repertoire dynamics through time in multiple individuals responding to an immune challenge.

Detecting selection in B cell Ig sequences is critical to understanding affinity maturation and can provide insights into Ag-driven selection in normal and pathologic immune responses. The most common sequence-based methods for detecting selection analyze the ratio of replacement and silent mutations using a binomial statistical analysis. However, these approaches have been criticized for low sensitivity. An alternative method is based on the analysis of lineage trees constructed from sets of clonally related Ig sequences. Several tree shape measures have been proposed as indicators of selection that can be statistically compared across cohorts. However, we show that tree shape analysis is confounded by underlying experimental factors that are difficult to control for in practice, including the sequencing depth and number of generations in each clone. Thus, although lineage tree shapes may reflect selection, their analysis alone is an unreliable measure of in vivo selection. To usefully capture the information provided by lineage trees, we propose a new method that applies the binomial statistical framework to mutations identified based on lineage tree structure. This hybrid method is able to detect selection with increased sensitivity in both simulated and experimental data sets. We anticipate that this approach will be especially useful in the analysis of large-scale Ig sequencing data sets generated by high-throughput sequencing technologies.

Analyses of somatic hypermutation (SHM) patterns in B cell immunoglobulin (Ig) sequences contribute to our basic understanding of adaptive immunity, and have broad applications not only for understanding the immune response to pathogens, but also to determining the role of SHM in autoimmunity and B cell cancers. Although stochastic, SHM displays intrinsic biases that can confound statistical analysis, especially when combined with the particular codon usage and base composition in Ig sequences. Analysis of B cell clonal expansion, diversification, and selection processes thus critically depends on an accurate background model for SHM micro-sequence targeting (i.e., hot/cold-spots) and nucleotide substitution. Existing models are based on small numbers of sequences/mutations, in part because they depend on data from non-coding regions or non-functional sequences to remove the confounding influences of selection. Here, we combine high-throughput Ig sequencing with new computational analysis methods to produce improved models of SHM targeting and substitution that are based only on synonymous mutations, and are thus independent of selection. The resulting "S5F" models are based on 806,860 Synonymous mutations in 5-mer motifs from 1,145,182 Functional sequences and account for dependencies on the adjacent four nucleotides (two bases upstream and downstream of the mutation). The estimated profiles can explain almost half of the variance in observed mutation patterns, and clearly show that both mutation targeting and substitution are significantly influenced by neighboring bases. While mutability and substitution profiles were highly conserved across individuals, the variability across motifs was found to be much larger than previously estimated. The model and method source code are made available at http://clip.med.yale.edu/SHM.

High-throughput immunoglobulin sequencing promises new insights into the somatic hypermutation and antigen-driven selection processes that underlie B-cell affinity maturation and adaptive immunity. The ability to estimate positive and negative selection from these sequence data has broad applications not only for understanding the immune response to pathogens, but is also critical to determining the role of somatic hypermutation in autoimmunity and B-cell cancers. Here, we develop a statistical framework for Bayesian estimation of Antigen-driven SELectIoN (BASELINe) based on the analysis of somatic mutation patterns. Our approach represents a fundamental advance over previous methods by shifting the problem from one of simply detecting selection to one of quantifying selection. Along with providing a more intuitive means to assess and visualize selection, our approach allows, for the first time, comparative analysis between groups of sequences derived from different germline V(D)J segments. Application of this approach to next-generation sequencing data demonstrates different selection pressures for memory cells of different isotypes. This framework can easily be adapted to analyze other types of DNA mutation patterns resulting from a mutator that displays hot/cold-spots, substitution preference or other intrinsic biases.

There are known differences in the properties of hair cells along the tonotopic axis of the avian auditory epithelium, the basilar papilla (BP). To determine the genetic basis of these differences, we compared gene expression between the high- (HF), middle-, and low-frequency (LF) thirds of 0-day-old chick auditory epithelia. RNA amplified from each sample was hybridized to whole-genome chicken arrays and GeneSpring software was used to identify differentially expressed genes. Two thousand six hundred sixty-three genes were found to be differentially expressed between the HF and LF segments, using a fold-change cutoff of 2 and a p value of 0.05. Many ion channel genes were differentially expressed between the HF and LF regions of the BP, an expression pattern that was previously established for some but not all of these genes. Quantitative PCR was used to verify tonotopic expression of 15 genes, including KCNMA1 (Slo) and its alternatively spliced STREX exon. Gene set enrichment analyses (GSEA) were performed on the microarray data and revealed many microRNA gene sets significantly enriched in the HF relative to the LF end, suggesting a tonotopic activity gradient. GSEA also suggested differential activity of the kinases protein kinase C and protein kinase A at the HF and LF ends, an interesting corollary to the observation that there is tonotopic expression of the STREX exon that confers on Slo sensitivity to the activity of kinases. Taken together, these results suggest mechanisms of induction and maintenance of tonotopicity and enhance our understanding of the complex nature of proximal-distal gene expression gradients in the chicken BP.

Background: Genome-wide transcriptional profiling of patient blood samples offers a powerful tool to investigate underlying disease mechanisms and personalized treatment decisions. Most studies are based on analysis of total peripheral blood mononuclear cells (PBMCs), a mixed population. In this case, accuracy is inherently limited since cell subset-specific differential expression of gene signatures will be diluted by RNA from other cells. While using specific PBMC subsets for transcriptional profiling would improve our ability to extract knowledge from these data, it is rarely obvious which cell subset(s) will be the most informative., Results: We have developed a computational method (Subset Prediction from Enrichment Correlation, SPEC) to predict the cellular source for a pre-defined list of genes (i.e. a gene signature) using only data from total PBMCs. SPEC does not rely on the occurrence of cell subset-specific genes in the signature, but rather takes advantage of correlations with subset-specific genes across a set of samples. Validation using multiple experimental datasets demonstrates that SPEC can accurately identify the source of a gene signature as myeloid or lymphoid, as well as differentiate between B cells, T cells, NK cells and monocytes. Using SPEC, we predict that myeloid cells are the source of the interferon-therapy response gene signature associated with HCV patients who are non-responsive to standard therapy., Conclusions: SPEC is a powerful technique for blood genomic studies. It can help identify specific cell subsets that are important for understanding disease and therapy response. SPEC is widely applicable since only gene expression profiles from total PBMCs are required, and thus it can easily be used to mine the massive amount of existing microarray or RNA-seq data.

Background: Auditory hair cells spontaneously regenerate following injury in birds but not mammals. A better understanding of the molecular events underlying hair cell regeneration in birds may allow for identification and eventually manipulation of relevant pathways in mammals to stimulate regeneration and restore hearing in deaf patients., Methodology/principal Findings: Gene expression was profiled in forskolin treated (i.e., proliferating) and quiescent control auditory epithelia of post-hatch chicks using an Affymetrix whole-genome chicken array after 24 (n = 6), 48 (n = 6), and 72 (n = 12) hours in culture. In the forskolin-treated epithelia there was significant (p<0.05; >two-fold change) upregulation of many genes thought to be relevant to cell cycle control and inner ear development. Gene set enrichment analysis was performed on the data and identified myriad microRNAs that are likely to be upregulated in the regenerating tissue, including microRNA181a (miR181a), which is known to mediate proliferation in other systems. Functional experiments showed that miR181a overexpression is sufficient to stimulate proliferation within the basilar papilla, as assayed by BrdU incorporation. Further, some of the newly produced cells express the early hair cell marker myosin VI, suggesting that miR181a transfection can result in the production of new hair cells., Conclusions/significance: These studies have identified a single microRNA, miR181a, that can cause proliferation in the chicken auditory epithelium with production of new hair cells.

B lymphocytes producing high-affinity Abs are critical for protection from extracellular pathogens, such as bacteria and parasites. The process by which high-affinity B cells are selected during the immune response has never been elucidated. Although it has been shown that high-affinity cells directly outcompete low-affinity cells in the germinal center (GC), whether there are also intrinsic differences between these cells has not been addressed. It could be that higher affinity cells proliferate more rapidly or are more likely to enter cell cycle, thereby outgrowing lower affinity cells. Alternatively, higher affinity cells could be relatively more resistant to cell death in the GC. By comparing high- and low-affinity B cells for the same Ag, we show here that low-affinity cells have an intrinsically higher death rate than do cells of higher affinity, even in the absence of competition. This suggests that selection in the GC reaction is due at least in part to the control of survival of higher affinity B cells and not by a proliferative advantage conferred upon these cells compared with lower affinity B cells. Control over survival rather than proliferation of low- and high-affinity B cells in the GC allows greater diversity not only in the primary response but also in the memory response.

Statistical methods based on the relative frequency of replacement mutations in B lymphocyte Ig V region sequences have been widely used to detect the forces of selection that shape the B cell repertoire. However, current methods produce an unexpectedly high frequency of false positives and are sensitive to intrinsic biases of somatic hypermutation that can give the appearance of selection. The new statistical test proposed here provides a better trade-off between sensitivity and specificity compared with previous approaches. The low specificity of existing methods was shown in silico to result from an interaction between the effects of positive and negative selection. False detection of positive selection was confirmed in vivo through a re-analysis of published sequence data from diffuse large B cell lymphomas, highlighting the need for re-analysis of some existing studies. The sensitivity of the proposed method to detect selection was validated using new Ig transgenic mouse models in which positive selection was expected to be a significant force, as well as with a simulation-based approach. Previous concerns that intrinsic biases of somatic hypermutation could give the appearance of selection were addressed by extending the current mutation models to more fully account for the impact of microsequence on relative mutability and to include transition bias. High specificity was confirmed using a large set of non-productively rearranged Ig sequences. These results show that selection can be detected in vivo with high specificity using the new method proposed here, allowing greater insight into the existence and direction of antigen-driven selection.