Correlated binary regression with covariates specific to each binary observation

Regression methods are considered for the analysis of correlated binary data when each binary observation may have its own covariates. It is argued that binary response models that condition on some or all binary responses in a given "block" are useful for studying certain types of dependencies, but not for the estimation of marginal response probabilities or pairwise correlations. Fully parametric approaches to these latter problems appear to be unduly complicated except in such special cases as the analysis of paired binary data. Hence, a generalized estimating equation approach is advocated for inference on response probabilities and correlations. Illustrations involving both small and large block sizes are provided.     

Vegetation responses to past climatic variation

Vegetation responses to climatic change can be studied retrospectively by utilizing the Quaternary fossil record. There has been controversy over the extent to which major changes in vegetation patterns at the continental scale lag behind the climatic changes that drive them, and to what extent vegetation can ever be said to be in equilibrium with climate. The equilibrium question has no single answer. The predominant mode of vegetation response to climatic change depends on the space and time frame and resolution of the data set in which the response is observed.Vegetation (as observed on particular space and time scales) can be in dynamic equilibrium with climate if its response time is sufficiently fast in relation to the rate of climatic change to which it is observed to be responding. Several processes can be involved in the response: successional, migrational, edaphic, and evolutionary. Successional response times can be deduced from forest succession models. The other processes are less well understood and different ideas exist concerning their rates. According to one hypothesis, migrational lags caused delays of thousands of years in the postglacial dynamics of forest composition. The alternative hypothesis explains these changes as dynamic equilibrium responses to changes in climatic seasonality and climatic anomaly patterns. Neither hypothesis need be universally true; gradient analysis and forest succession models are among the techniques that can be used in inferential tests of these hypotheses for particular space-time regions.Dynamic equilibrium may often be a reasonable approximation for the responses of the broadest continental-scale forest patterns to orbitally induced climatic changes. But as spatial and temporal frames of observation are diminished and resolution increased, biotic processes must eventually come to dominate. At sufficiently fine scales the main observable phenomena are successional responses to natural disturbance events. The late-Quaternary record of vegetation change allows a choice of observation scales and thus provides a continuum of possibilities for study, ranging from long-term dynamic bioclimatology to more conventional vegetation dynamics.I thank Margaret Davis, Honor Prentice, Jim Ritchie, Al Solomon, Geoff Spaulding and Tom Webb for their reviews of earlier drafts. Research supported by a US Department of Energy, Carbon Dioxide Research Division, grant to Brown University and a Swedish Natural Science Research Council grant to the project SlsSimulation of Natural Forest Dynamics.I thank Margaret Davis, Honor Prentice, Jim Ritchie, Al Solomon, Geoff Spaulding and Tom Webb for their reviews of earlier drafts. Research supported by a US Department of Energy, Carbon Dioxide Research Division, grant to Brown University and a Swedish Natural Science Research Council grant to the project SlsSimulation of Natural Forest Dynamics.     

SB subregion of the human major histocompatibility complex: gene organization, allelic polymorphism and expression in transformed cells

The SB region of the human major histocompatibility complex (MHC) has been cloned from cosmid and lambda phage libraries made from the human B-lymphoblastoid cell line Priess (DR4/4, DC4/4, SB3/4). Two alpha genes and two beta genes are encoded in the 100 kb long SB region in the order SB alpha-SB beta-SX alpha-SX beta. The SB alpha and SB beta genes encode the alpha and beta subunits of the SB subset of class II MHC molecules. Both the SX alpha and the SX beta genes are pseudogenes in the haplotype examined. From the isolated clones, the two haplotypes of the Priess cell line, SB3 and SB4, are distinguished by nucleotide sequencing and blot hybridization analyses. Restriction site polymorphisms between the SB3 and SB4 clones were observed only in relatively small regions of the SB beta and SX beta genes. A mouse macrophage cell line was transfected with one of the cosmid clones containing both SB alpha and SB beta genes. Expression of the alpha and beta genes was detected by fluorescene-activated cell sorting (FACS) and two-dimensional gel electrophoresis using SB-specific monoclonal antibodies.     

Isolation and characterization of a large,neurite-associated glycoconjugate from neuroblastoma cells

A high molecular weight glycoconjugate has been isolated from neurite-producing neuronal tumor cells in culture and has been designated as I(0) based on its elution characteristics in gel filtration chromatography. This molecule cannot be found in a variety of nonneuronal cells. I(0) is found in the substratum-attached material or cell fraction of neurite-producing neuroblastoma cells, depending upon culture conditions. It is found in the substratum-bound fraction of B104 rat neuroblastoma cells during serum starvation and in the EGTA-detached cell fraction of B104 cells grown in chemically defined N2 medium. It occurs only in the cell fraction of the human neuroblastoma line Platt. Examination of behavioral variants of the B104 rat line further strengthens the association of I(0) with neurite production; the constitutive neurite-producing E(R)B9 variant contains I(0) while the non-neurite-producing E(R)A11 variant does not. I(0) is large, eluting in the void volume of sepharose-CL2B columns. Radioiodination of intact cells with lactoperoxidase shows I(0) to be a cell surface component. Metabolic radiolabeling studies show that it contains a high proportion of polysaccharide to protein, does not contain mannose, and is unsulfated. Alkaline borohydride reduction release two size classes of large polysaccharide chain. The alkaline reduction results, along with the mannose incorporation studies, show the presence of O-glycosidic linkages and few, if any, N-linkages. Resistance to nitrous acid deamination, insensitivity to glycosaminoglycan lyases, and the absence of sulfation, indicate that I(0) does not contain the glycosaminoglycans hyaluronic acid, chondroitin-, dermatan-, or heparin- sulfates. Affinity column chromatography reveals high binding affinity of I(0) to polyornithine and no binding to gelatin (collagen) or the glycosaminoglycans hyaluronate and heparin. These studies describe a unique high molecular weight glycoconjugate on the surface of neurite-producing neuroblastoma cell lines from two species.     

Serum and Urine Fibrinogen-Fibrin-related Antigen (F.R.-Antigen) Levels in Renal Disease

The concentration of serum fibrinogen-fibrin-related antigen (F.R.-antigen) was measured in a group of 142 patients with various renal disorders, in 38 of whom urine F.R.-antigen was also estimated. Raised serum F.R.-antigen levels were present in 48% of the patients, with no particular preponderance in any diagnostic category apart from acute reversible intrinsic renal failure in which high levels were invariably present. Significantly-raised serum levels were also present in the patients with microangiopathic haemolytic anaemia and in those with the more severe degrees of renal impairment. Urine F.R.-antigen was increased in 34 of the 38 patients. The amount of F.R.-antigen in the urine correlated with the degree of proteinuria but not with the serum F.R.-antigen levels. The evidence relating to intravascular coagulation in renal disease is reviewed, and it is suggested that there is a high incidence of localized fibrinogen or fibrin degradation in the kidney, which is related more to factors such as the presence of uraemia and microangiopathic haemolytic anaemia rather than to the diagnostic category.     

The Distribution of Acidic and Neutral Peptidases in Barley and Wheat Kernels

The activity of barley and wheat peptidases which hydrolyze alpha-N-benzoyl-dl-arginine-p-nitronnilide (BAPA) and alpha-N-benzoyl-l-arginine ethyl ester (BAEE) has been measured in proximal and distal portions of ungerminated grain and in these tissues during 6 and 7 day incubations. The proximal portion of ungerminated barleys contained the major part of both the acidic (BAPA-ase and acidic BAEE-ase) and neutral (neutral BAEE-ase) peptidases. In ungcrminated wheat these acidic and neutral peptidases were nearly evenly distributed between the proximal and dislal portions. Commercial wheat embryo was very high in acidic peptidase but contained no neutral peptidase. During the germination of both wheat and barleys, acidic and neutral peptidase activity in the seedlings increased with time. No such consistent increase was observed for aleurone and starchy endosperm tissue for any of these grains. Aleurone and starchy endosperm tissue incubated in the absence of the proximal portion of the kernel showed reduced peptidase activities.     

Relative apparent synapomorphy analysis (RASA). I: The statistical measurement of phylogenetic signal

We have developed a new approach to the measurement of phylogenetic signal in character state matrices called relative apparent synapomorphy analysis (RASA). RASA provides a deterministic, statistical measure of natural cladistic hierarchy (phylogenetic signal) in character state matrices. The method works by determining whether a measure of the rate of increase of cladistic similarity among pairs of taxa as a function of phenetic similarity is greater than a null equiprobable rate of increase. Our investigation of the utility and limitations of RASA using simulated and bacteriophage T7 data sets indicates that the method has numerous advantages over existing measures of signal. A first advantage is computational efficiency. A second advantage is that RASA employs known methods of statistical inference, providing measurable sensitivity and power. The performance of RASA is examined under various conditions of branching evolution as the number of characters, character states per character, and mutations per branch length are varied. RASA appears to provide an unbiased and reliable measure of phylogenetic signal, and the general approach promises to be useful in the development of new techniques that should increase the rigor and reliability of phylogenetic estimates.      

Molecular evolution of voltage-sensitive ion channel genes: on the origins of electrical excitability

We have analyzed nucleic acid and amino acid sequence alignments of a variety of voltage-sensitive ion channels, using several methods for phylogenetic tree reconstruction. Ancient duplications within this family gave rise to three distantly related groups, one consisting of the Na+ and Ca++ channels, another the K+ channels, and a third including the cyclic nucleotide-binding channels. A series of gene duplications produced at least seven mammalian homologues of the Drosophila Shaker K+ channel; clones of only three of these genes are available from all three mammalian species examined (mouse, rat, and human), pointing to specific genes that have yet to be recovered in one or another of these species. The Shaw-related K+ channels and the Na+ channel family have also undergone considerable expansion in mammals, relative to flies. These expansions presumably reflect the needs of the high degree of physiological and neuronal complexity of mammals. Analysis of the separate domains of the four-domain channels (Ca++ and Na+) supports their having evolved by two sequential gene duplications and implies the historical existence of a functional two-domain channel.      

mtDNA diversity in rhesus monkeys reveals overestimates of divergence time and paraphyly with neighboring species

Reconstructions of the human-African great ape phylogeny by using mitochondrial DNA (mtDNA) have been subject to considerable debate. One confounding factor may be the lack of data on intraspecific variation. To test this hypothesis, we examined the effect of intraspecific mtDNA diversity on the phylogenetic reconstruction of another Plio- Pleistocene radiation of higher primates, the fascicularis group of macaque (Macaca) monkey species. Fifteen endonucleases were used to identify 10 haplotypes of 40-47 restriction sites in M. mulatta, which were compared with similar data for the other members of this species group. Interpopulational, intraspecific mtDNA diversity was large (0.5%- 4.5%), and estimates of divergence time and branching order incorporating this variation were substantially different from those based on single representatives of each species. We conclude that intraspecific mtDNA diversity is substantial in at least some primate species. Consequently, without prior information on the extent of genetic diversity within a particular species, intraspecific variation must be assessed and accounted for when reconstructing primate phylogenies. Further, we question the reliability of hominoid mtDNA phylogenies, based as they are on one or a few representatives of each species, in an already depauperate superfamily of primates.      

Genetic differentiation in Hippocrepis emerus (Leguminosae): allozyme and DNA fingerprint variation in disjunct Scandinavian populations

The structure of genetic variation in disjunct Scandinavian populations of Hippocrepis emerus was studied using allozymes and DNA fingerprinting. Variation in the three native regional populations in Scandinavia was compared with that in a recently introduced population in Sweden. In contrast to the recently introduced population, the native Scandinavian isolates of H. emerus showed high levels of allozyme fixation and low levels of DNA diversity. Variation in allozymes and at DNA fingerprint loci showed closely congruent patterns of geographic variation, with pronounced differentiation between the native Norwegian and Swedish isolates of the species. The structure of genetic variation in native Scandinavian H. emerus is interpreted in terms of historical population bottlenecks and founder events during the species' postglacial immigration into Scandinavia.