A distribution model for heritability.

A regression model to predict quantiles of narrow sense individual and family mean heritabilities is developed and used to predict confidence intervals either directly or via a generalized beta distribution model. Extensive simulations of balanced sib analysis trials in randomized complete block designs and normal distributed environmental and additive genetic effects confirmed that heritabilities follow a beta distribution even in cases with up to 10% of the data missing at random. The new model is both more accurate and more precise than commonly used alternatives based on "exact" chi 2 distributions and Satterthwaites approximations to the degrees of freedom. Estimates of the expected heritability and a Taylor approximation of the standard error of the heritability are needed as input to the quantile model. Applications of the presented models for estimating confidence intervals and as an aid in the design of experiments are provided.     

Up-regulation of cell-surface alpha4beta2 neuronal nicotinic receptors by lower temperature and expression of chimeric subunits.

The predominant nicotinic acetylcholine receptor (nAChR) expressed in vertebrate brain is a pentamer containing alpha4 and beta2 subunits. In this study we have examined how temperature and the expression of subunit chimeras can influence the efficiency of cell-surface expression of the rat alpha4beta2 nAChR. Functional recombinant alpha4beta2 nAChRs, showing high affinity binding of nicotinic radioligands (K(d) = 41 +/- 22 pM for [(3)H]epibatidine), are expressed in both stably and transiently transfected mammalian cell lines. Despite this, only very low levels of alpha4beta2 nAChRs can be detected on the cell surface of transfected mammalian cells maintained at 37 degrees C. At 30 degrees C, however, cells expressing alpha4beta2 nAChRs show a 12-fold increase in radioligand binding (with no change in affinity), and a 5-fold up-regulation in cell-surface receptors with no increase in total subunit protein. In contrast to "wild-type" alpha4 and beta2 subunits, chimeric nicotinic/serotonergic subunits ("alpha4chi" and "beta2chi") are expressed very efficiently on the cell surface (at 30 degrees C or 37 degrees C), either as hetero-oligomeric complexes (e.g. alpha4chi+beta2 or alpha4chi+beta2chi) or when expressed alone. Compared with alpha4beta2 nAChRs, expression of complexes containing chimeric subunits typically results in up to 20-fold increase in nicotinic radioligand binding sites (with no change in affinity) and a similar increase in cell-surface receptor, despite a similar level of total chimeric and wild-type protein.     

A note on the conditional approach to interval estimation in the calibration problem.

In the calibration problem, the need to construct a confidence interval to estimate the unknown chi 0 arises when the null hypothesis of zero slope is rejected. Otherwise, the resulting confidence interval will be infinite to reflect the fact that the slope of the regression line may be zero. Under the condition of rejecting the hypothesis of zero slope, we study the properties of the conditional coverage rate of the calibration confidence interval. The conditional coverage rate (P1) is a function of the slope, distance between chi 0 and the mean of the trailing sample means, the sum of squares of chi, and n. When the true slope is close to 0 and chi 0 is away from means, P1 can go down to 0. On the other hand, as the power of testing zero slope reaches 1, with or without chi 0 close to means, P1 will tend to the desired nominal coverage rate. In summary, one should choose a reasonably small alpha in testing zero slope to avoid constructing a confidence interval for chi 0 when the true slope is 0. In addition, it is desirable to have high power in testing zero slope so that the resulting confidence interval will maintain the desired coverage rate when using the conditional approach in the calibration problem.     

Modulation of radiation-induced chromosome aberrations by DMSO, an OH radical scavenger. 1: Dose-response studies in human lymphocytes exposed to 220 kV X-rays

Human G0 lymphocytes were exposed to 220 kV X-radiation in the presence or absence of DMSO, an efficient selective scavenger of OH radicals. Our studies demonstrate that DMSO affects a concentration-dependent modulation of induced asymmetrical aberrations in human lymphocytes exposed to approximately 3.0 Gy, with maximum protectible fractions of approximately 70 percent at DMSO concentrations of greater than or equal to 1 M. The dose dependency for dicentrics in lymphocytes acutely exposed to X-ray doses of 0.51 to 4.98 Gy in the absence of DMSO is adequately described by the linear-quadratic dose-response function Y = alpha D + beta D2. Data from duplicate cultures exposed in the presence of 1 M DMSO produce an excellent fit to the regression function modified as follows: Y(+ DMSO) = alpha(delta D) + beta(delta D)2 where the 'dose modifying' factor delta = 0.501. We interpret these findings as providing evidence that OH radical-mediated lesions in DNA account for approximately 50 percent of the dose dependency for dicentrics resulting from either one-track or two-track events, following exposures of non-cycling cells to moderate-to-high doses of low LET radiation. These data may be used in additional calculations to derive an estimate of approximately 6 x 10(8) s-1 for the rate of reaction of OH radicals with DNA targets involved in aberration formation.     

Structure of the 113Cd3β domains from Homarus americanus metallothionein-1: hydrogen bonding and solvent accessibility of sulfur atoms

The three-dimensional structures of the isolated Cd(3)beta domains from Homarus americanus metallothionein have been determined by NMR methods in order to establish a set of beta-domain structures for comparative analysis. First, it was determined that the Cd-cysteine connectivities forming the Cd(3)S(9) metal center were identical to those observed for the beta(N) domain in the native holoprotein. Time- and temperature-dependence studies of the (113)Cd and (1)H 1D-NMR spectra indicated that the beta(N) domain undergoes slow conformational changes before reaching an equilibrium structure. In addition to structural information provided by the metal-to-cysteine connectivities, Phi, chi(1) and chi(2) angle constraints, three H(N...)S hydrogen bond interactions were also determined from a long-range optimized (1)H(N)-(113)Cd HMQC experiment. A simulated annealing protocol was applied to the distance and angle constraints obtained from the 2D-NMR experiments to calculate the three-dimensional structure of the synthetic Cd(3)beta(N) domain of lobster metallothionein. Structure-reactivity relationships are proposed for the reactions of Cd(3)beta domains with 5,5'-dithiobis(2-nitrobenzoate), based on comparisons of surface exposure of sulfur atoms of the lobster and rabbit Cd(3)beta domain structures. Finally, the surface exposure of the beta domains of lobster is compared with beta domains from mammalian metallothioneins.     

Side-chain modeling with an optimized scoring function

Modeling side-chain conformations on a fixed protein backbone has a wide application in structure prediction and molecular design. Each effort in this field requires decisions about a rotamer set, scoring function, and search strategy. We have developed a new and simple scoring function, which operates on side-chain rotamers and consists of the following energy terms: contact surface, volume overlap, backbone dependency, electrostatic interactions, and desolvation energy. The weights of these energy terms were optimized to achieve the minimal average root mean square (rms) deviation between the lowest energy rotamer and real side-chain conformation on a training set of high-resolution protein structures. In the course of optimization, for every residue, its side chain was replaced by varying rotamers, whereas conformations for all other residues were kept as they appeared in the crystal structure. We obtained prediction accuracy of 90.4% for chi(1), 78.3% for chi(1 + 2), and 1.18 A overall rms deviation. Furthermore, the derived scoring function combined with a Monte Carlo search algorithm was used to place all side chains onto a protein backbone simultaneously. The average prediction accuracy was 87.9% for chi(1), 73.2% for chi(1 + 2), and 1.34 A rms deviation for 30 protein structures. Our approach was compared with available side-chain construction methods and showed improvement over the best among them: 4.4% for chi(1), 4.7% for chi(1 + 2), and 0.21 A for rms deviation. We hypothesize that the scoring function instead of the search strategy is the main obstacle in side-chain modeling. Additionally, we show that a more detailed rotamer library is expected to increase chi(1 + 2) prediction accuracy but may have little effect on chi(1) prediction accuracy.     

Statistical and energetic analysis of side-chain conformations in oligopeptides

The distributions of side-chain conformations in 258 crystal structures of oligopeptides have been analyzed. The sample contains 321 residues having side chains that extend beyond the C beta atom. Statistically observed preferences of side-chain dihedral angles are summarized and correlated with stereochemical and energetic constraints. The distributions are compared with observed distributions in proteins of known X-ray structures and with computed minimum-energy conformations of amino acid derivatives. The distributions are similar in all three sets of data, and they appear to be governed primarily by intraresidue interactions. In side chains with no beta-branching, the most important interactions that determine chi 1 are those between the C gamma H2 group and atoms of the neighboring peptide groups. As a result, the g- conformation (chi 1 congruent to -60 degrees) occurs most frequently for rotation around the C alpha-C beta bond in oligopeptides, followed by the t conformation (chi 1 congruent to 180 degrees), while the g+ conformation (chi 1 congruent to 60 degrees) is least favored. In residues with beta-branching, steric repulsions between the C gamma H2 or C gamma H3 groups and backbone atoms govern the distribution of chi 1. The extended (t) conformation is highly favored for rotation around the C beta-C gamma and C gamma-C delta bonds in unbranched side chains, because the t conformer has a lower energy than the g+ and g- conformers in hydrocarbon chains. This study of the observed side-chain conformations has led to a refinement of one of the energy parameters used in empirical conformational energy computations.     

OPUS-Rota: a fast and accurate method for side-chain modeling

In this paper, we introduce a fast and accurate side-chain modeling method, named OPUS-Rota. In a benchmark comparison with the methods SCWRL, NCN, LGA, SPRUCE, Rosetta, and SCAP, OPUS-Rota is shown to be much faster than all the methods except SCWRL, which is comparably fast. In terms of overall chi (1) and chi (1+2) accuracies, however, OPUS-Rota is 5.4 and 8.8 percentage points better, respectively, than SCWRL. Compared with NCN, which has the best accuracy in the literature, OPUS-Rota is 1.6 percentage points better for overall chi (1+2) but 0.3 percentage points weaker for overall chi (1). Hence, our algorithm is much more accurate than SCWRL with similar execution speed, and it has accuracy comparable to or better than the most accurate methods in the literature, but with a runtime that is one or two orders of magnitude shorter. In addition, OPUS-Rota consistently outperforms SCWRL on the Wallner and Elofsson homology-modeling benchmark set when the sequence identity is greater than 40%. We hope that OPUS-Rota will contribute to high-accuracy structure refinement, and the computer program is freely available for academic users.     

Designed protein G core variants fold to native-like structures: sequence selection by ORBIT tolerates variation in backbone specification

The solution structures of two computationally designed core variants of the beta 1 domain of streptococcal protein G (G beta 1) were solved by (1)H NMR methods to assess the robustness of amino acid sequence selection by the ORBIT protein design package under changes in protein backbone specification. One variant has mutations at three of 10 core positions and corresponds to minimal perturbations of the native G beta 1 backbone. The other, with mutations at six of 10 positions, was calculated for a backbone in which the separation between G beta 1's alpha-helix and beta-sheet was increased by 15% relative to native G beta 1. Exchange broadening of some resonances and the complete absence of others in spectra of the sixfold mutant bespeak conformational heterogeneity in this protein. The NMR data were sufficiently abundant, however, to generate structures of similar, moderately high quality for both variants. Both proteins adopt backbone structures similar to their target folds. Moreover, the sequence selection algorithm successfully predicted all core chi(1) angles in both variants, five of six chi(2) angles in the threefold mutant and four of seven chi(2) angles in the sixfold mutant. We conclude that ORBIT calculates sequences that fold specifically to a geometry close to the template, even when the template is moderately perturbed relative to a naturally occurring structure. There are apparently limits to the size of acceptable perturbations: In this study, the larger perturbation led to undesired dynamic behavior.     

TNFalpha induces chromosomal abnormalities independent of ROS through IKK, JNK, p38 and caspase pathways

A role for pro-inflammatory cytokines in inflammation-related cancers has been suggested, but mechanisms are not defined. Here, we demonstrate that treatment of HeLa cells with TNFalpha increases chromosomal aberration. In contrast, IL-1beta did not increase, but rather decreased chromosomal aberration. TNFalpha and IL-1beta increased the production of H2O2 to similar levels in cells, suggesting that increased production of reactive oxygen species might not be the premier factor involved. Reducing H2O2 through overexpression of catalase or treatment of cells with NAC or BHA did not have an effect on TNF-induced chromosomal aberration. TNFalpha-induced NO production has been implicated in DNA damage. Inhibiting NO did not reduce TNF-induced chromosomal aberration. Inhibiting IKK, JNK, and p38 kinase as well as caspases decreased TNF-induced chromosomal aberration, and a correlation between TNF-induced apoptosis and CA generation was not found. Single-strand DNA breaks give rise to double-strand breaks, which then results in chromosomal breaks, when replication forks reach the single-strand breaks during S-phase. In cells progressing through S-phase, TNFalpha activation of IKK, JNK, and p38 is significantly reduced. However, these kinases were activated by IL-1beta in S-phase. The possibility that these pathways, in a TNF-specific manner, may regulate either the generation of single- and double-strand breaks or their repair, thereby resulting in increased chromosomal aberration, is discussed.     

Synthesis, crystal structure, and molecular conformation of peptide N-Boc-L-Pro-dehydro-Phe-L-Gly-OH

The peptide N-Boc-L-Pro-dehydro-Phe-L-Gly-OH was synthesized by the usual workup procedure and finally coupling the N-Boc-L-Pro-dehydro-Phe to glycine. The peptide crystallizes in monoclinic space group P2(1) with a = 8.951(4) A, b = 5.677(6) A, c = 21.192(11) A, beta = 96.97(4) degrees, V = 1069(1) A3, Z = 2, dm = 1.295(5) Mgm-3, and dc = 1.297(4) Mgm-3. The structure was determined by direct methods using SHELXS86. The structure was refined by the block-diagonal least-squares procedure to an R value of 0.074 for 1002 observed reflections. The C alpha 2-C beta 2 distance of 1.33(2) A is an appropriate double bond length. The angle C alpha 2-C beta 2-C gamma 2 is 133(1) degrees. The peptide backbone torsion angles are theta 1 = -167(1) degrees, omega 0 = 179(1) degrees, phi 1 = -48(1) degrees, psi 1 = 137(1) degrees, omega 1 = 175(1) degrees, phi 2 = 65(2) degrees, psi 2 = 15(2) degrees, omega 2 = -179(1) degrees, and phi 3 = -166(1) degrees. These values show that the Boc group has a trans-trans conformation while the peptide backbone adopts a beta-turn II conformation, which is stabilized by an intramolecular hydrogen bond of length of 3.05(1) A. The structures of dehydro-Phe containing peptides suggest that the dehydro-Phe promotes the beta-turn II conformation. The five-membered pyrrolidine ring of the Pro residue adopts an ideal C gamma-exo conformation with torsion angles chi 1(1) = -24(1) degrees, chi 2(1) = 34(1) degrees, chi 3(1) = -30(1) degrees, chi 4(1) = 15(1) degrees, and theta 0(1) = 6(1) degrees. The side-chain torsion angles in dehydro-Phe are chi 1(2) = -1(2) degrees, chi 2,1(2) = -176(1) degrees, and chi 2,2(2) = 8(2) degrees. The plane of C alpha 2-C beta 2-C gamma 2 is rotated with respect to the plane of the phenyl ring at 7(1) degrees, which indicates that the atoms of the side chain of dehydro-Phe are essentially coplanar. The molecules form a 2(1) screw axis related hydrogen-bonded rows along the b axis.     

Relation of branching angles to optimality for four cost principles

The literature has suggested that branching angles depend on some principle of optimality. Most often cited are the minimization of lumen surface, volume, power and drag. The predicted angles depend on the principle applied, chi and alpha. Assuming flow o r chi, chi can be determined from r chi 0 = r chi 1 + r chi 2 when the radii of the parent (r0) major (r1) and minor (r2) daughters are known. The term alpha = r2/r1. Using different values for chi and alpha, we present graphs for the major and minor branching angles theta 1 and theta 2 and psi = theta 1 + theta 2 for each of the four optimization principles. Because psi is almost independent of alpha for values of chi and alpha found in 198 junctions taken from a human pulmonary artery, we are able to produce a plot of psi versus chi for each of the four principles on one graph. A junction can be provisionally classified as optimizing for a given principle if, knowing chi, the psi obs - psi pred is least for that principle. We find that this nomographic classification agrees almost perfectly with a previous classification based on a more exacting measure, the percent cost index I, where I = observed cost/minimum cost. We explain why this is to be expected in most but not all cases. First we generate a contoured percent cost surface of c = I - 100 around the optimally located junction, J, and superimpose a surface of equal angular deviations a = psi pred-psi obs. We find that c increases and a usually increases with distance from J as the actual junction moves along a straight line away from J. We then produce a plot of c versus a for two competing principles. A comparison of the principles demonstrates that, for most cases, a is smaller for the principle which has the smaller c value.     

Catalytic properties of human liver alcohol dehydrogenase isoenzymes

Human liver alcohol dehydrogenase (ADH) exists in multiple molecular forms which arise from the association of eight different types of subunits, alpha, beta 1, beta 2, beta 3, gamma 1, gamma 2, pi, and chi, into active dimeric molecules. A genetic model accounts for this multiplicity as products of five gene loci, ADH1 through ADH5. Polymorphism occurs at two loci, ADH2 and ADH3, which encode the beta and gamma subunits. All of the known homodimeric and heterodimeric isoenzymes have been isolated and purified to homogeneity. Analysis of the steady-state kinetic properties and substrate and inhibitor specificities has shown substantial differences in the catalytic properties of the isoenzymes. For example, the Km values for NAD+ and ethanol vary as much as 1,000-fold among the isoenzymes. Some of the differences in catalytic properties can be related to specific amino acid substitutions in the ADH isoenzymes.     

1H-NMR studies on nucleotide binding to the sarcoplasmic reticulum Ca2+ ATPase. Determination of the conformations of bound nucleotides by the measurement of proton-proton transferred nuclear Overhauser enhancements

The glycosidic bond torsion angles and the conformations of the ribose of Mg2+ATP, Mg2+ADP and Mg2+AdoPP[NH]P (magnesium adenosine 5'-[beta, gamma-imido]triphosphate) bound to Ca2+ATPase, both native and modified with fluorescein isothiocyanate (FITC), in intact sarcoplasmic reticulum have been determined by the measurement of proton-proton transferred nuclear Overhauser enhancements by 1H-NMR spectroscopy. This method shows clearly the existence of a low-affinity ATP binding site after modification of the high-affinity site with FITC. For all three nucleotides bound to both the high-affinity (catalytic) site and the low-affinity site, we find that the conformation about the glycosidic bond is anti, the conformation of the ribose 3'-endo of the N type and the conformation about the ribose C4'-C5' bond either gauche-trans or trans-gauche. The values for the glycosidic bond torsion angles chi (O4'-C1'-N9-C4) for Mg2+ATP, Mg2+ADP and Mg2+AdoPP[NH]P bound to the low-affinity site of FITC-modified Ca2+ATPase are approximately equal to 270 degrees, approximately equal to 260 degrees and approximately equal to 240 degrees respectively. In the case of the nucleotides bound to the high-affinity (catalytic) site of native Ca2+ATPase, chi lies in the range 240-280 degrees.     

Conformations of disulfide bridges in proteins

The conformational characteristics of disulfide bridges in proteins have been analyzed using a dataset of 22 protein structures, available at a resolution of less than or equal to 2.0 A, containing a total of 72 disulfide crosslinks. The parameters used in the analysis include (phi, psi) values at Cys residues, bridge dihedral angles chi ss, chi i1, chi j1, chi i2, and chi j2, the distances C alpha i-C alpha j and C beta i-C beta j between the C alpha and C beta atoms of Cys(i) and Cys(j). Eight families of bridge conformations with three or more occurrences have been identified on the basis of these stereochemical parameters. The most populated family corresponds to the "left handed spiral" identified earlier by Richardson [1981) Adv. Protein Chem. 34, 167-330). Disulfide bridging across antiparallel extended strands is observed in alpha-lytic protease, crambin, and beta-trypsin and this structure is shown to be very similar to those obtained in small cystine peptides. Solvent accessible surface area calculations show that the overwhelming majority of disulfide bridges are inaccessible to solvent.     

The effect of oxygen concentration on the X-ray induction of chromosome aberrations in human lymphocytes.

In vitro dose--response curves, for unstable chromosome aberration induction in human lymphocytes under conditions of full oxygenation or of anoxia, have been obtained using 250 kVp X-rays. Dicentric yields have been fitted to the quadratic function Y = alpha D + beta D2. An Oxygen Enhancement Ratio (OER) ranging from 7.2 to 2.7 was calculated from the coefficients of these curves possibly indicating that the data do not fit to the dose-modifying model of the oxygen effect although the differences are not statistically significant. A similar analysis of total aberration data yields an OER of 3.6 to 2.7 fitting much more satisfactorily to the dose-modifying model. Variations in dicentric yield induced by 3.0 Gy and 0.75 Gy of X-rays with increasing oxygen concentration were plotted and for each dose a constant dicentric yield was observed at oxygen levels of 2 and 250 ppm. Above 250 ppm yields increased steeply up to about 1% oxygen and then more gradually to a maximum at 100% oxygen.     

~(60)Coγ线引起上海真厉螨染色体畸变的研究

本文首次报道用~(6O)Coγ线照射一种革螨——上海真厉螨引起的染色体畸变的研究。用~(6O)Coγ线(剂量1—50Krad)照射雌性革螨,引起的染色体畸变类型有:染色体裂隙、断片、微小体、环形染色体、粉碎化和多倍体,染色体断片是最常见的畸变类型,并观察到微核的形成。染色体畸变率随照射剂量增加而增高,辐射剂量与畸变率之间存在密切相关(相关系数为0.85,P<0.025),配得曲线回归方程为Y=3.27+14.49lg(X+1)。     

Further Tests of a Recombination Model in Which χ Removes the Recd Subunit from the Recbcd Enzyme of Escherichia Coli

When one of two infecting lambda phage types in a replication-blocked cross is chi + and DNA packaging is divorced from the RecBCD-chi interaction, complementary chi-stimulated recombinants are recovered equally in mass lysates only if the chi + parent is in excess in the infecting parental mixture. Otherwise, the chi 0 recombinant is recovered in excess. This observation implies that, along with the chi 0 chromosome, two chi + parent chromosomes are involved in the formation of each chi + recombinant. The trimolecular nature of chi +-stimulated recombination is manifest in recombination between lambda and a plasmid. When lambda recombines with a plasmid via the RecBCD pathway, the resulting chromosome has an enhanced probability of undergoing lambda x lambda recombination in the interval into which the plasmid was incorporated. These two observations support a model in which DNA is degraded by Exo V from cos, the sequence that determines the end of packaged lambda DNA and acts as point of entry for RecBCD enzyme, to chi, the DNA sequence that stimulates the RecBCD enzyme to effect recombination. The model supposes that chi acts by ejecting the RecD subunit from the RecBCD enzyme with two consequences. (1) ExoV activity is blocked leaving a highly recombinagenic, frayed duplex end near chi, and (2) as the enzyme stripped of the RecD subunit travels beyond chi it is competent to catalyze reciprocal recombination.     

Orientations of the tryptophan 9 and 11 side chains of the gramicidin channel based on deuterium nuclear magnetic resonance spectroscopy.

Deuterium nuclear magnetic resonance spectroscopy was used to investigate the orientations of the indole rings of Trp9 and Trp11 in specific indole-d5-labeled samples of gramicidin A incorporated into dimyristoyl phosphatidylcholine bilayers in the beta 6.3 channel conformation. The magnitudes and signs of the deuterium quadrupolar splittings were fit to the rings and assigned to specific ring bonds, using a full rotation search of the chi 1 and chi 2 angles of each Trp and a least-squares method. Unique assignments were obtained. The data and assignments are in close agreement with four sets of (chi 1, chi 2) angles for each Trp in which the indole N-H is oriented toward the membrane's exterior surface. (Four additional sets of (chi 1, chi 2) angles with the N-H's pointing toward the membrane interior are inconsistent with previous observations.) One of the sets of (chi 1, chi 2) angles for each Trp is consistent with the corresponding Trp orientation found by Arsen'ev et al. (1986. Biol. Membr. 3:1077-1104) for gramicidin in sodium dodecyl sulfate micelles. Together, the 1H and 2H nuclear magnetic resonance methods suggest that the Trp9 and Trp11 side chain orientations could be very similar in dimyristoyl phosphatidylcholine membranes and in sodium dodecyl sulfate micelles. The data for Trp11 could be fit using a static quadrupolar coupling constant of 180 kHz under the assumption that the ring is essentially immobile. By contrast, Trp9 could be fit only under the assumption that the quadrupolar splittings for ring 9 are reduced by approximately 14% due to motional averaging. Such a difference in motional averaging between rings 11 and 9 is also consistent with the 15N data of Hu et al. (1993. Biochemistry. 32:7035-7047).     

Allotypes of the constant region of the rabbit T cell receptor beta-2 chain

Our laboratory previously reported that there was restriction fragment length polymorphism of TCR C beta genes in rabbits. EcoRI digests of DNA from different rabbits gave fragments of 9 and 6 kb (C beta a) or 14 and 6 kb (C beta b) that hybridized to a C beta cDNA probe. We also reported that the 9- and 14-kb types segregated as Mendelian traits and that there were allotypic differences in the first exon of the C beta 1 genes of C beta a and C beta b animals. Here we report the DNA sequence of the C beta 2 gene present in the 6-kb EcoRI fragment from a C beta b animal and compare the exon sequences with that of a cDNA from a C beta a animal. We find replacement changes in the first and third exons that probably represent allotypic forms of the rabbit C beta 2 gene. The genomic DNA 5' of exon 1 of both beta 1 and beta 2 contain alternating purine/pyrimidine repeat sequences. The genomic C beta 2 has an open reading frame of 69 amino acids in frame with exon 1 similar to a longer one previously found 5' of exon 1 of C beta 1. Further 5' of this region, rabbit C beta 1 and C beta 2 DNA sequences are only about 66% similar. Both the C beta 1 and C beta 2 sequences have two chi sequences; one in exon 1 with a perfect match and one in the intron downstream of exon 1 with one mismatch. Alternating purine/pyrimidine repeats and chi sequences found in rabbit C beta 1 and C beta 2 genes may have contributed to process(es) of gene duplication and/or conversion.