Rectangular tolerance regions and multivariate normal reference regions in laboratory medicine

Reference intervals are widely used in the interpretation of results of biochemical and physiological tests of patients. When there are multiple biochemical analytes measured from each subject, a multivariate reference region is needed. Because of their greater specificity against false positives, such reference regions are more desirable than separate univariate reference intervals that disregard the cross-correlations between variables. Traditionally, under multivariate normality, reference regions have been constructed as ellipsoidal regions. This approach suffers from a major drawback: it cannot detect component-wise extreme observations. In the present work, procedures are developed to construct rectangular reference regions in the multivariate normal setup. The construction is based on the criteria for tolerance intervals. The problems addressed include the computation of a rectangular tolerance region and simultaneous tolerance intervals. Also addressed is the computation of mixed reference intervals that include both two-sided and one-sided limits, simultaneously. A parametric bootstrap approach is used in the computations, and the accuracy of the proposed methodology is assessed using estimated coverage probabilities. The problem of sample size determination is also addressed, and the results are illustrated using examples that call for the computation of reference regions.     

Plant adaptability in karst regions

Journal of Plant Research - Karst ecosystems are formed by dissolution of soluble rocks, usually with conspicuous landscape features, such as sharp peaks, steep slopes and deep valleys. The plants...     

Plant diversity in mediterranean-climate regions

The high plant diversity of mediterranean-climate regions has attracted much attention over the past few years. This review discusses patterns and determinants of local, differential and regional plant diversity in all five regions. Local diversity shows great variation within and between regions and explanations for these patterns invoke a wide range of hypotheses. Patterns of regional diversity are the result of differential speciation and extinction rates during the Quaternary. These rates have been influenced more by the incidence of fire and the severity of climate change than by environmental heterogeneity. All regions have a high number of rare and locally endemic taxa that survive as small populations, many of which are threatened by habitat transformation.     

Crystalline regions in collagen fibrils

A new image processing technique, content-dependent anisotropic spatial frequency filtering, has been developed to visualize the location and orientation of crystalline regions in collagen fibril cross-sections. The results show that most crystalline regions are oriented with their approximately 4 nm periodicity directed radially from the fibril centre. This periodicity corresponds to the separation between rows of molecular ends in the quasi-hexagonal molecular packing scheme. The extent of crystallinity increases with radius and frequently the lattice is either continuously distorted or interrupted by sharp discontinuities.     

Molecular evolution of duplicated amylase gene regions in Drosophila melanogaster: evidence of positive selection in the coding regions and selective constraints in the cis-regulatory regions

In this study, we randomly sampled Drosophila melanogaster from Japanese and Kenyan natural populations. We sequenced duplicated (proximal and distal) Amy gene regions to test whether the patterns of polymorphism were consistent with neutral molecular evolution. F(st) between the two geographically distant populations, estimated from Amy gene regions, was 0.084, smaller than reported values for other loci, comparing African and Asian populations. Furthermore, little genetic differentiation was found at a microsatellite locus (DROYANETSB) in these samples (G'st = -0.018). The results of several tests (Tajima's, Fu and Li's, and Wall's tests) were not significantly different from neutrality. However, a significantly higher level of fixed replacement substitutions was detected by a modified McDonald and Kreitman test for both populations. This indicates that positive selection occurred during or immediately after the speciation of D. melanogaster. Sliding-window analysis showed that the proximal region 1, a part of the proximal 5' flanking region, was conserved between D. melanogaster and its sibling species, D. simulans. An HKA test was significant when the proximal region 1 was compared with the 5' flanking region of Alcohol dehydrogenase (Adh), indicating a severe selective constraint on the Amy proximal region 1. These results suggest that natural selection has played an important role in the molecular evolution of Amy gene regions in D. melanogaster.     

Small-sample confidence regions in exponential families

This article presents an algorithm for small-sample conditional confidence regions for two or more parameters for any discrete regression model in the generalized linear interactive model family. Regions are constructed by careful inversion of conditional hypothesis tests. This method presupposes the use of approximate or exact techniques for enumerating the sample space for some components of the vector of sufficient statistics conditional on other components. Such enumeration may be performed exactly or by exact or approximate Monte Carlo, including the algorithms of Kolassa and Tanner (1994, Journal of the American Statistical Association 89, 697-702; 1999, Biometrics 55, 246-251). This method also assumes that one can compute certain conditional probabilities for a fixed value of the parameter vector. Because of a property of exponential families, one can use this set of conditional probabilities to directly compute the conditional probabilities associated with any other value of the vector of the parameters of interest. This observation dramatically reduces the computational effort required to invert the hypothesis test to obtain the confidence region. To construct a region with confidence level 1 - alpha, the algorithm begins with a grid of values for the parameters of interest. For each parameter vector on the grid (corresponding to the current null hypothesis), one transforms the initial set of conditional probabilities using exponential tilting and then calculates the p value for this current null hypothesis. The confidence region is the set of parameter values for which the p value is at least alpha.     

SQUARE--determining reliable regions in sequence alignments

The Server for Quick Alignment Reliability Evaluation (SQUARE) is a Web-based version of the method we developed to predict regions of reliably aligned residues in sequence alignments. Given an alignment between a query sequence and a sequence of known structure, SQUARE is able to predict which residues are reliably aligned. The server accesses a database of profiles of sequences of known three-dimensional structures in order to calculate the scores for each residue in the alignment. SQUARE produces a graphical output of the residue profile-derived alignment scores along with an indication of the reliability of the alignment. In addition, the scores can be compared against template secondary structure, conserved residues and important sites.     

Modeling heterochromatin regions in transgenic mice

Transgenic mice carrying bovine satellite DNA IV were obtained. The size of the transgene integrated into the mouse genome was approximately 390 kb (about 100 transgene copies) as determined by a semiquantitative PCR. Restriction analysis with isoschizomeric restrictases HpaII and MspI, showed that the alien DNA was methylated. In the genome of a transgenic founder male, two integration sites for satellite DNA IV were revealed by in situ hybridization and in situ PCR. These sites are situated on two different chromosomes: in pericentromeric heterochromatin and within a chromosomal arm. In transgenic mice, de novo formation of heterochromatin regions (C-block and the CMA3 disk within the centromeric heterochromatin of another chromosome) was revealed by C-banding and staining with chromomycin A3. This formation is not characteristic of mice, because their chromosomes normally contain no interstitial C-blocks or sequences intensely stained by chromomycin A3.     

Helical regions in transfer ribonucleic acid

Summary By following the thermal denaturation of transfer RNA (t-RNA) over a wide wavelength range it can be demonstrated that t-RNA dissociates in at least two steps. The lower melting step seems to melt cooperatively, while this can not be stated with confidence about the higher melting region. It is shown that the use of difference spectra of two-stranded transitions gives similar, but different results than the use of difference spectra for three stranded-transitions. It is thus shown that difference spectra can not be used (at least at present) to distinguish strandedness in RNA. Possible applications of the method of difference spectra are discussed.     

Potential aggregation prone regions in biotherapeutics

Aggregation of a biotherapeutic is of significant concern and judicious process and formulation development is required to minimize aggregate levels in the final product. Aggregation of a protein in solution is driven by intrinsic and extrinsic factors. In this work we have focused on aggregation as an intrinsic property of the molecule. We have studied the sequences and Fab structures of commercial and non-commercial antibody sequences for their vulnerability towards aggregation by using sequence based computational tools to identify potential aggregation-prone motifs or regions. The mAbs in our dataset contain 2 to 8 aggregation-prone motifs per heavy and light chain pair. Some of these motifs are located in variable domains, primarily in CDRs. Most aggregation-prone motifs are rich in β branched aliphatic and aromatic residues. Hydroxyl-containing Ser/Thr residues are also found in several aggregation-prone motifs while charged residues are rare. The motifs found in light chain CDR3 are glutamine (Q)/asparagine (N) rich. These motifs are similar to the reported aggregation promoting regions found in prion and amyloidogenic proteins that are also rich in Q/N, aliphatic and aromatic residues. The implication is that one possible mechanism for aggregation of mAbs may be through formation of cross-β structures and fibrils. Mapping on the available Fab – receptor/antigen complex structures reveals that these motifs in CDRs might also contribute significantly towards receptor/antigen binding. Our analysis identifies the opportunity and tools for simultaneous optimization of the therapeutic protein sequence for potency and specificity while reducing vulnerability towards aggregation.     

Organizer regions in marine colonial hydrozoans

Organizers are specific tissue regions of developing organisms that provide accuracy and robustness to the body plan formation. Hydrozoan cnidarians (both solitary and colonial) require organizer regions for maintaining the regular body patterning during continuous tissue dynamics during asexual reproduction and growth. While the hypostomal organizer of the solitary Hydra has been studied relatively well, our knowledge of organizers in colonial hydrozoans remains fragmentary and incomplete. As colonial hydrozoans demonstrate an amazing diversity of morphological and life history traits, it is of special interest to investigate the organizers specific for particular ontogenetic stages and particular types of colonies. In the present study we aimed to assess the inductive capacities of several candidate organizer regions in hydroids with different colony organization. We carried out grafting experiments on colonial hydrozoans belonging to Leptothecata and Anthoathecata. We confirmed that the hypostome tip is an organizer in the colonial Anthoathecata as it is in the solitary polyp Hydra. We also found that the posterior tip of the larva is an organizer in hydroids regardless of the peculiarities of their metamorphosis mode and colony structure. We show for the first time that the shoot growing tip, which can be considered a key evolutionary novelty of Leptothecata, is an organizer region. Taken together, our data demonstrate that organizers function throughout the larval and polypoid stages in colonial hydroids.     

Scaffold attachment regions in centromere-associated DNA

Due to indications that kinetochore proteins are an integral part of the protein scaffold component of the chromosome (Earnshaw et al. 1984), we chose to map the distribution of scaffold attachment regions (SARs) at centromeres. Using the SAR mapping assay of Mirkovitch et al., Southern blots were prepared and probed with ³²P-labeled fragments from the human 1.9 kb centromeric α-satellite repeat unit of chromosome 1 or the 1.7 kb centromeric α-satellite repeat unit of chromosome 16. Our results demonstrated the presence of one SAR site per 1.9 kb repeat unit in chromosome 1, and every 1.7 kb repeat unit in chromosome 16, separated by regions of small DNA loops over the length of the α-satellite regions. We also identified several in vitro vertebrate topoisomerase II and cenP-B consensus sequences throughout the chromosome 1 α-satellite region using computer and base ratio analysis, to address the question as to why some α-satellite regions are SAR related and others are not. To provide in situ indications of SAR localization in the human genome, SAR DNA and non-SAR DNA were prepared following lithium 3,5-di-iodosalicylate extraction. Sequences protected from DNAse I digestion by SAR proteins, as compared with unprotected DNA that was digested by the enzyme, was labeled with biotin-UTP, hybridized to chromosomal DNA in situ, and then detected with fluorescein-avidin-DCS. Both SAR and non-SAR DNA selectively labeled virtually all centromeric regions of the human metaphase karyotype. Chromosomal arms were less strongly bound by SAR DNA, with a pattern that followed the chromosomal axis. In the more condensed chromosomes an R-banding pattern was evident. In general, labeling patterns produced by both SAR and non-SAR fractions were similar, as expected from the indications that SAR DNAs are heterogenous in sequence and do not form a specific class of sequences. We conclude that centromeric regions of several, possibly all, human metaphase chromosomes are also regions where the chromosomal axis contains loops, smaller in size than in the arms and where attachment sites are concentrated. This clustering of SARs may be responsible in part for the tight chromatin packing associated with the primary constriction of the centromeric region. Received: 10 October 1995; in revised form: 10 May 1996 / Accepted: 13 May 1996     

Inferring selection in partially sequenced regions

A common approach for identifying loci influenced by positive selection involves scanning large portions of the genome for regions that are inconsistent with the neutral equilibrium model or represent outliers relative to the empirical distribution of some aspect of the data. Once identified, partial sequence is generated spanning this more localized region in order to quantify the site-frequency spectrum and evaluate the data with tests of neutrality and selection. This method is widely used as partial sequencing is less expensive with regard to both time and money. Here, we demonstrate that this approach can lead to biased maximum likelihood estimates of selection parameters and reduced rejection rates, with some parameter combinations resulting in clearly misleading results. Most significantly, for a commonly used sample size in Drosophila population genetics (i.e., n = 12), the estimate of the target of selection has a large mean square error and the strength of selection is severely under estimated when the true selected site has not been sampled. We propose sequencing approaches that are much more likely to accurately localize the target and estimate the strength of selection. Additionally, we examine the performance of a commonly used test of selection under a variety of recurrent and single sweep models.