Hematopoietic cell types: prototype for a revised cell ontology

The Cell Ontology (CL) aims for the representation of in vivo and in vitro cell types from all of biology. The CL is a candidate reference ontology of the OBO Foundry and requires extensive revision to bring it up to current standards for biomedical ontologies, both in its structure and its coverage of various subfields of biology. We have now addressed the specific content of one area of the CL, the section of the ontology dealing with hematopoietic cells. This section has been extensively revised to improve its content and eliminate multiple inheritance in the asserted hierarchy, and the groundwork has been laid for structuring the hematopoietic cell type terms as cross-products incorporating logical definitions built from relationships to external ontologies, such as the Protein Ontology and the Gene Ontology. The methods and improvements to the CL in this area represent a paradigm for improvement of the entire ontology over time.

Conserved cryptic recombination signals in Vkappa gene segments are cleaved in small pre-B cells

Background

The cleavage of recombination signals (RS) at the boundaries of immunoglobulin V, D, and J gene segments initiates the somatic generation of the antigen receptor genes expressed by B lymphocytes. RS contain a conserved heptamer and nonamer motif separated by non-conserved spacers of 12 or 23 nucleotides. Under physiologic conditions, V(D)J recombination follows the "12/23 rule" to assemble functional antigen-receptor genes, i.e., cleavage and recombination occur only between RS with dissimilar spacer types. Functional, cryptic RS (cRS) have been identified in V_H gene segments; these V_H cRS were hypothesized to facilitate self-tolerance by mediating V_H → V_HDJ_H replacements. At the Igκ locus, however, secondary, de novo rearrangements can delete autoreactive VκJκ joins. Thus, under the hypothesis that V-embedded cRS are conserved to facilitate self-tolerance by mediating V-replacement rearrangements, there would be little selection for Vκ cRS. Recent studies have demonstrated that V_H cRS cleavage is only modestly more efficient than V(D)J recombination in violation of the 12/23 rule and first occurs in pro-B cells unable to interact with exogenous antigens. These results are inconsistent with a model of cRS cleavage during autoreactivity-induced V_H gene replacement.

Results

To test the hypothesis that cRS are absent from Vκ gene segments, a corollary of the hypothesis that the need for tolerizing V_H replacements is responsible for the selection pressure to maintain V_H cRS, we searched for cRS in mouse Vκ gene segments using a statistical model of RS. Scans of 135 mouse Vκ gene segments revealed highly conserved cRS that were shown to be cleaved in the 103/BCL2 cell line and mouse bone marrow B cells. Analogous to results for V_H cRS, we find that Vκ cRS are conserved at multiple locations in Vκ gene segments and are cleaved in pre-B cells.

Conclusion

Our results, together with those for V_H cRS, support a model of cRS cleavage in which cleavage is independent of BCR-specificity. Our results are inconsistent with the hypothesis that cRS are conserved solely to support receptor editing. The extent to which these sequences are conserved, and their pattern of conservation, suggest that they may serve an as yet unidentified purpose.

An improved ontological representation of dendritic cells as a paradigm for all cell types

BACKGROUND

Recent increases in the volume and diversity of life science data and information and an increasing emphasis on data sharing and interoperability have resulted in the creation of a large number of biological ontologies, including the Cell Ontology (CL), designed to provide a standardized representation of cell types for data annotation. Ontologies have been shown to have significant benefits for computational analyses of large data sets and for automated reasoning applications, leading to organized attempts to improve the structure and formal rigor of ontologies to better support computation. Currently, the CL employs multiple is_a relations, defining cell types in terms of histological, functional, and lineage properties, and the majority of definitions are written with sufficient generality to hold across multiple species. This approach limits the CL's utility for computation and for cross-species data integration.

RESULTS

To enhance the CL's utility for computational analyses, we developed a method for the ontological representation of cells and applied this method to develop a dendritic cell ontology (DC-CL). DC-CL subtypes are delineated on the basis of surface protein expression, systematically including both species-general and species-specific types and optimizing DC-CL for the analysis of flow cytometry data. We avoid multiple uses of is_a by linking DC-CL terms to terms in other ontologies via additional, formally defined relations such as has_function.

CONCLUSION

This approach brings benefits in the form of increased accuracy, support for reasoning, and interoperability with other ontology resources. Accordingly, we propose our method as a general strategy for the ontological representation of cells. DC-CL is available from http://www.obofoundry.org.

Multiple, conserved cryptic recombination signals in VH gene segments: detection of cleavage products only in pro B cells

Receptor editing is believed to play the major role in purging newly formed B cell compartments of autoreactivity by the induction of secondary V(D)J rearrangements. In the process of immunoglobulin heavy (H) chain editing, these secondary rearrangements are mediated by direct V_H-to-J_H joining or cryptic recombination signals (cRSs) within V_H gene segments. Using a statistical model of RS, we have identified potential cRSs within V_H gene segments at conserved sites flanking complementarity-determining regions 1 and 2. These cRSs are active in extrachromosomal recombination assays and cleaved during normal B cell development. Cleavage of multiple V_H cRSs was observed in the bone marrow of C57BL/6 and RAG2:GFP and μMT congenic animals, and we determined that cRS cleavage efficiencies are 30–50-fold lower than a physiological RS. cRS signal ends are abundant in pro–B cells, including those recovered from μMT mice, but undetectable in pre– or immature B cells. Thus, V_H cRS cleavage regularly occurs before the generation of functional preBCR and BCR. Conservation of cRSs distal from the 3′ end of V_H gene segments suggests a function for these cryptic signals other than V_H gene replacement.

Antigen-independent secondary Igκ rearrangement in B-cell development – implications for receptor editing (B39)

Receptor editing by secondary Vκ-to-Jκ rearrangements is thought to be driven by reactivity to self. However, the evidence supporting antigen-driven editing is open to other interpretations; analysis of Igκ excision circles from birds and mice indicates continuing V(D)J recombination in cis after both functional (F) and non-functional (nF) rearrangements. We have determined the extent and nature of secondary Vκ-to-Jκ rearrangements by analyzing their intermediate cleavage products in the 103/Bcl2 cell line and in purified compartments of developing B cells from normal, Jκ−/+, and IgH transgenic mice. In all cases, we find that the ratio of F:nF VκJκ joints replaced by secondary rearrangement is approximately 1:2, the ratio expected for rearrangement without feedback of any kind. Secondary rearrangement is initiated in small pre-B cells and in BrdU pulse-labeling studies, we detected no evidence for immature to pre-B “retrograde” differentiation. Thus, “receptor editing” occurs in B cells that do not express surface Ig. Significantly, in 3H9 IgH transgenic mice, replacement rearrangements exhibited F:nF ratios of 1:2, regardless of whether the initial rearrangements were permissive or non-permissive for autoreactivity. Our results demonstrate continuing Vκ-to-Jκ rearrangement on a single chromosome without feedback and imply that “receptor editing” is an illusion of antigen-driven selection.

V(D)J recombinase-mediated processing of coding junctions at cryptic recombination signal sequences in peripheral T cells during human development

V(D)J recombinase mediates rearrangements at immune loci and cryptic recombination signal sequences (cRSS), resulting in a variety of genomic rearrangements in normal lymphocytes and leukemic cells from children and adults. The frequency at which these rearrangements occur and their potential pathologic consequences are developmentally dependent. To gain insight into V(D)J recombinase-mediated events during human development, we investigated 265 coding junctions associated with cRSS sites at the hypoxanthine-guanine phosphoribosyltransferase (HPRT) locus in peripheral T cells from 111 children during the late stages of fetal development through early adolescence. We observed a number of specific V(D)J recombinase processing features that were both age and gender dependent. In particular, TdT-mediated nucleotide insertions varied depending on age and gender, including percentage of coding junctions containing N-nucleotide inserts, predominance of GC nucleotides, and presence of inverted repeats (Pr-nucleotides) at processed coding ends. In addition, the extent of exonucleolytic processing of coding ends was inversely related to age. We also observed a coding-partner-dependent difference in exonucleolytic processing and an age-specific difference in the subtypes of V(D)J-mediated events. We investigated these age- and gender-specific differences with recombination signal information content analysis of the cRSS sites in the human HPRT locus to gain insight into the mechanisms mediating these developmentally specific V(D)J recombinase-mediated rearrangements in humans.

Reassignment of the murine 3'TRDD1 recombination signal sequence

T cell receptor genes are assembled in developing T lymphocytes from discrete V, D, and J genes by a site-specific somatic rearrangement mechanism. A flanking recombination signal, composed of a conserved heptamer and a semiconserved nonamer separated by 12 or 23 variable nucleotides, targets the activity of the rearrangement machinery to the adjoining V, D, and J genes. Following the rearrangement of V, D, or J genes, their respective recombination signals are ligated together. Although these signal joints are allegedly invariant, created by the head-to-head abuttal of the heptamers, some do exhibit junctional diversity. Recombination signals were initially identified by comparison and alignment of germ-line sequences with the sequence of rearranged genes. However, their overall low level of sequence conservation makes their characterization solely from sequence data difficult. Recently, computational analysis unraveled correlations between nucleotides at several positions scattered within the spacer and recombination activity, so that it is now possible to identify putative recombination signals and determine and predict their recombination efficiency. In this paper, we analyzed the variability introduced in signal joints generated after rearrangement of the TRDD1 and TRDD2 genes in murine thymocytes. The recurrent presence of identical nucleotides inserted in these signal joints led us to reconsider the location and sequence of the TRDD1 recombination signal. By combining molecular characterization and computational analysis, we show that the functional TRDD1 recombination signal is shifted inside the putative coding sequence of the TRDD1 gene and, consequently, that this gene is shorter than indicated in the databases.

Neonate-primed CD8+ memory cells rival adult-primed memory cells in antigen-driven expansion and anti-viral protection

Immunizations early in life, when the host is most susceptible to infection, allow protective immunological memory to develop. Decreasing the dose of Cas-Br-E murine leukemia virus when priming neonatal mice results in adult-like, Type 1 protective responses, but the resulting memory cell populations are smaller than after adult priming. After secondary challenge, virus-specific CD8+ memory cell populations expand twice as much in neonate-primed mice as in adult-primed mice. We found that when equivalent numbers of virus-specific cells were transferred into virus-susceptible mice, protection from disease was similar whether donor, immune mice were primed as neonates or adults, and IL-4 did not alter in vivo virus-specific CD8+ memory cell effector function. Hence, neonate-primed CD8+ cells develop into memory cells that rival adult-primed cells in proliferation and effector function.

SoDA: implementation of a 3D alignment algorithm for inference of antigen receptor recombinations

MOTIVATION

The antigen receptors of adaptive immunity-T-cell receptors and immunoglobulins-are encoded by genes assembled stochastically from combinatorial libraries of gene segments. Immunoglobulin genes then experience further diversification through hypermutation. Analysis of the somatic genetics of the immune response depends explicitly on inference of the details of the recombinatorial process giving rise to each of the participating antigen receptor genes. We have developed a dynamic programming algorithm to perform this reconstruction and have implemented it as web-accessible software called SoDA (Somatic Diversification Analysis).

RESULTS

We tested SoDA against a set of 120 artificial immunoglobulin sequences generated by simulation of recombination and compared the results with two other widely used programs. SoDA inferred the correct gene segments more frequently than the other two programs. We further tested these programs using 30 human immunoglobulin genes from Genbank and here highlight instances where the recombinations inferred by the three programs differ. SoDA appears generally to find more likely recombinations.

Computational tools for understanding sequence variability in recombination signals

The recombination signals (RSs) that guide V(D)J rearrangement are remarkably diverse. In mice, fewer than 16% of RSs carry consensus heptamers and nonamers and none also contain a consensus spacer sequence. It is increasingly clear that this variability regulates recombination: genetic variability in RSs may help enforce allelic exclusion, determine the general nature of antigen receptor repertoires, and mitigate autoreactivity in B lymphocytes. The great diversity of RSs has largely precluded, however, empiric determinations of how RS sequence affects recombination. For example, 4(39) unique 23-RSs are possible or approximately 3 x 10(23) sequences; some 7 x 10(13) unique 23-RSs can be produced just by changes in the spacer. In contrast, the recombination activities of only 100 or so RSs have been measured, and it is unlikely that the activities of even a tiny fraction of extant RSs can be determined. We have addressed the problem of how sequence determines the efficiency of RS templates by generating computational models that describe the correlation structure of mouse RSs. These models successfully predict RS activity and identify functional, cryptic RSs (cRSs). These models permit studies to identify RSs and cRSs for empiric study and constitute a tool useful for understanding RS structure and function.

A functional analysis of the spacer of V(D)J recombination signal sequences

During lymphocyte development, V(D)J recombination assembles antigen receptor genes from component V, D, and J gene segments. These gene segments are flanked by a recombination signal sequence (RSS), which serves as the binding site for the recombination machinery. The murine Jbeta2.6 gene segment is a recombinationally inactive pseudogene, but examination of its RSS reveals no obvious reason for its failure to recombine. Mutagenesis of the Jbeta2.6 RSS demonstrates that the sequences of the heptamer, nonamer, and spacer are all important. Strikingly, changes solely in the spacer sequence can result in dramatic differences in the level of recombination. The subsequent analysis of a library of more than 4,000 spacer variants revealed that spacer residues of particular functional importance are correlated with their degree of conservation. Biochemical assays indicate distinct cooperation between the spacer and heptamer/nonamer along each step of the reaction pathway. The results suggest that the spacer serves not only to ensure the appropriate distance between the heptamer and nonamer but also regulates RSS activity by providing additional RAG:RSS interaction surfaces. We conclude that while RSSs are defined by a "digital" requirement for absolutely conserved nucleotides, the quality of RSS function is determined in an "analog" manner by numerous complex interactions between the RAG proteins and the less-well conserved nucleotides in the heptamer, the nonamer, and, importantly, the spacer. Those modulatory effects are accurately predicted by a new computational algorithm for "RSS information content." The interplay between such binary and multiplicative modes of interactions provides a general model for analyzing protein-DNA interactions in various biological systems.

Prospective estimation of recombination signal efficiency and identification of functional cryptic signals in the genome by statistical modeling

The recombination signals (RS) that guide V(D)J recombination are phylogenetically conserved but retain a surprising degree of sequence variability, especially in the nonamer and spacer. To characterize RS variability, we computed the position-wise information, a measure correlated with sequence conservation, for each nucleotide position in an RS alignment and demonstrate that most position-wise information is present in the RS heptamers and nonamers. We have previously demonstrated significant correlations between RS positions and here show that statistical models of the correlation structure that underlies RS variability efficiently identify physiologic and cryptic RS and accurately predict the recombination efficiencies of natural and synthetic RS. In scans of mouse and human genomes, these models identify a highly conserved family of repetitive DNA as an unexpected source of frequent, cryptic RS that rearrange both in extrachromosomal substrates and in their genomic context.

Identification and utilization of arbitrary correlations in models of recombination signal sequences

BACKGROUND

A significant challenge in bioinformatics is to develop methods for detecting and modeling patterns in variable DNA sequence sites, such as protein-binding sites in regulatory DNA. Current approaches sometimes perform poorly when positions in the site do not independently affect protein binding. We developed a statistical technique for modeling the correlation structure in variable DNA sequence sites. The method places no restrictions on the number of correlated positions or on their spatial relationship within the site. No prior empirical evidence for the correlation structure is necessary.

RESULTS

We applied our method to the recombination signal sequences (RSS) that direct assembly of B-cell and T-cell antigen-receptor genes via V(D)J recombination. The technique is based on model selection by cross-validation and produces models that allow computation of an information score for any signal-length sequence. We also modeled RSS using order zero and order one Markov chains. The scores from all models are highly correlated with measured recombination efficiencies, but the models arising from our technique are better than the Markov models at discriminating RSS from non-RSS.

CONCLUSIONS

Our model-development procedure produces models that estimate well the recombinogenic potential of RSS and are better at RSS recognition than the order zero and order one Markov models. Our models are, therefore, valuable for studying the regulation of both physiologic and aberrant V(D)J recombination. The approach could be equally powerful for the study of promoter and enhancer elements, splice sites, and other DNA regulatory sites that are highly variable at the level of individual nucleotide positions.

The "dispensable" portion of RAG2 is necessary for efficient V-to-DJ rearrangement during B and T cell development

Previous in vitro studies defined the minimal regions of RAG1 and RAG2 essential for V(D)J recombination. In order to characterize the role of the C-terminal "dispensable" portion of RAG2, we generated core-RAG2 knock-in mice. We found that the core-RAG2-containing recombinase complex is selectively defective in catalyzing V-to-DJ rearrangement at the IgH and TCRbeta loci, resulting in partial developmental blocks in B and T lymphopoiesis. Analysis of recombination intermediates showed defects at the cleavage phase of the reaction. We also observed a reduction in overall recombinase activity in core-RAG2-expressing thymocytes, leading us to suggest that the interaction of a defective recombinase with RSS sequences unique to VH and Vbeta gene segments may underlie the specific V-to-DJ rearrangement defect in core-RAG2 mice.

The targeting of somatic hypermutation closely resembles that of meiotic mutation

We have compared the microsequence specificity of mutations introduced during somatic hypermutation (SH) and those introduced meiotically during neutral evolution. We have minimized the effects of selection by studying nonproductive (hence unselected) Ig V region genes for somatic mutations and processed pseudogenes for meiotic mutations. We find that the two sets of patterns are very similar: the mutabilities of nucleotide triplets are positively correlated between the somatic and meiotic sets. The major differences that do exist fall into three distinct categories: 1) The mutability is sharply higher at CG dinucleotides under meiotic but not somatic mutation. 2) The complementary triplets AGC and GCT are much more mutable under somatic than under meiotic mutation. 3) Triplets of the form WAN (W = T or A) are uniformly more mutable under somatic than under meiotic mutation. Nevertheless, the relative mutabilities both within this set and within the SAN (S = G or C) triplets are highly correlated with those under meiotic mutation. We also find that the somatic triplet specificity is strongly symmetric under strand exchange for A/T triplets as well as for G/C triplets in spite of the strong predominance of A over T mutations. Thus, we suggest that somatic mutation has at least two distinct components: one that specifically targets AGC/GCT triplets and another that acts as true catalysis of meiotic mutation.

The nucleotide-replacement spectrum under somatic hypermutation exhibits microsequence dependence that is strand-symmetric and distinct from that under germline mutation

Somatic mutation is a fundamental component of acquired immunity. Although its molecular basis remains undetermined, the sequence specificity with which mutations are introduced has provided clues to the mechanism. We have analyzed data representing over 1700 unselected mutations in V gene introns and nonproductively rearranged V genes to identify the sequence specificity of the mutation spectrum-the distribution of resultant nucleotides. In other words, we sought to determine what effects the neighboring bases have on what a given base mutates "to." We find that both neighboring bases have a significant effect on the mutation spectrum. Their influences are complicated, but much of the effect can be characterized as enhancing homogeneity of the mutated DNA sequence. In contrast to what has been reported for the sequence specificity of the "targeting" mechanism, that of the spectrum is notably symmetric under complementation, indicating little if any strand bias. We compared the spectrum to that found previously for germline mutations as revealed by analyzing pseudogene sequences. We find that the influences of nearest neighbors are quite different in the two datasets. Altogether, our findings suggest that the mechanism of somatic hypermutation is complex, involving two or more stages: introduction of mis-pairs and their subsequent resolution, each with distinct sequence specificity and strand bias.

Enhanced evolvability in immunoglobulin V genes under somatic hypermutation

Darwinian theory requires that mutations be produced in a nonanticipatory manner; it is nonetheless consistent to suggest that mutations that have repeatedly led to nonviable phenotypes would be introduced less frequently than others-if under appropriate genetic control. Immunoglobulins produced during infection acquire point mutations that are subsequently selected for improved binding to the eliciting antigen. We and others have speculated that an enhancement of mutability in the complementarity-determining regions (CDR; where mutations have a greater chance of being advantageous) and/or decrement of mutability in the framework regions (FR; where mutations are more likely to be lethal) may be accomplished by differential codon usage in concert with the known sequence specificity of the hypermutation mechanism. We have examined 115 nonproductively rearranged human Ig sequences. The mutation patterns in these unexpressed genes are unselected and therefore directly reflect inherent mutation biases. Using a chi2 test, we have shown that the number of mutations in the CDRs is significantly higher than the number of mutations found in the FRs, providing direct evidence for the hypothesis that mutations are preferentially targeted into the CDRs.

The distribution of variation in regulatory gene segments, as present in MHC class II promoters

Diversity in the antigen-binding receptors of the immune system has long been a primary interest of biologists. Recently it has been suggested that polymorphism in regulatory (noncoding) gene segments is of substantial importance as well. Here, we survey the level of variation in MHC class II gene promoters in man and mouse using extensive collections of published sequences together with unpublished sequences recently deposited by us in the EMBL gene bank using the Shannon entropy to quantify diversity. For comparison, we also apply our analysis to distantly related MHC class II promoters, as well as to class I promoters and to class II coding regions. We observe a high level of intraspecies variability, which in mouse but not in man is localized to a significant extent near the binding sites of transcription factors-sites that are conserved over longer evolutionary distances. This localization may both indicate and enhance heterozygote advantage, as the presence of two functionally different promoters would be expected to confer flexibility in the immune response.

Density-dependent prenatal androgen exposure as an endogenous mechanism for the generation of cycles in small mammal populations

Small mammal populations exhibit cyclic fluctuations in their population densities. Several hypotheses regarding the mechanisms underlying these population cycles have been advanced, but none has yet gained general approval. We propose here an endogenous mechanism based on the masculinization of female offspring in response to increased population levels. High population levels trigger the non-specific stress response resulting in high levels of circulating androgens in individuals of the population, including pregnant females. These androgens masculinize female offspring in utero, thereby reducing the reproductive capacity of the next generation and subsequently the population size. We have developed and analysed a mathematical model to investigate the possible role of prenatal androgen exposure in the generation of limit cycles. We find the locus of Hopf bifurcations for this model and show that limit cycles depend on three parameters: (1) the delay between birth and sexual maturation; (2) the slope of the function that relates average prenatal androgen exposure to total population density; and (3) the difference between the maximum birth rates of the low- and high-androgen exposed females. We derive the analytical form relating these parameters at the Hopf-bifurcation locus and discuss its biological ramifications. In brief, in each of these three parameters is sufficiently large, population cycles will results from the endogenous mechanism proposed.Copyright 1998 Academic Press Limited Copyright 1998 Academic Press Limited