A numerical method for calculating moments of coalescent times in finite populations with selection

A numerical method is developed for solving a nonstandard singular system of second-order differential equations arising from a problem in population genetics concerning the coalescent process for a sample from a population undergoing selection. The nonstandard feature of the system is that there are terms in the equations that approach infinity as one approaches the boundary. The numerical recipe is patterned after the LU decomposition for tridiagonal matrices. Although there is no analytic proof that this method leads to the correct solution, various examples are presented that suggest that the method works. This method allows one to calculate the expected number of segregating sites in a random sample of n genes from a population whose evolution is described by a model which is not selectively neutral.     

Staining Senile Plaques using Bodian's Method Modified with Methenamine

A new method is presented for staining various types of senile plaques isolated from the brains of patients with Alzheimer type dementia and related diseases in paraffin embedded sections using a modified Bodian's method with methenamine. This methenamine-Bodian method made it possible to observe diffuse plaques and other amyloid deposits which are barely detected by Bodian's original method. The staining of senile plaques by the method presented here was comparable to that of immunostaining with anti-β-protein. The new method also stained neurofibrillary tangles. Therefore, the methenamine-Bodian method could be widely used for the detection of senile changes in paraffin embedded sections from autopsied human brains.     

Random-effects models for serial observations with binary response

This paper presents a general mixed model for the analysis of serial dichotomous responses provided by a panel of study participants. Each subject's serial responses are assumed to arise from a logistic model, but with regression coefficients that vary between subjects. The logistic regression parameters are assumed to be normally distributed in the population. Inference is based upon maximum likelihood estimation of fixed effects and variance components, and empirical Bayes estimation of random effects. Exact solutions are analytically and computationally infeasible, but an approximation based on the mode of the posterior distribution of the random parameters is proposed, and is implemented by means of the EM algorithm. This approximate method is compared with a simpler two-step method proposed by Korn and Whittemore (1979, Biometrics 35, 795-804), using data from a panel study of asthmatics originally described in that paper. One advantage of the estimation strategy described here is the ability to use all of the data, including that from subjects with insufficient data to permit fitting of a separate logistic regression model, as required by the Korn and Whittemore method. However, the new method is computationally intensive.     

New approaches to the prediction of the folding of membrane proteins with redox function

A new method is elaborated for determining the hydropathy profile of membrane haemoproteins. The method is called membrane propensity for haemoproteins (MPH) and is based on the statistical analysis of the amino acid composition of the predicted transmembrane regions of cytochrome b from the bc1 and the b6f complexes. The accuracy of the MPH method in predicting the ends of the known transmembrane segments of the reaction center of Rhodopseudomonas viridis is higher than that obtained by hydropathy methods based on physico-chemical parameters. The MPH method is able to clearly exclude from the membrane polypeptides that are not consistently predicted to be transmembrane by other methods or techniques, for instance the region corresponding to helix IV of mitochondrial cytochrome b. A correlation has been found between the shape of the hydropathy profile of the transmembrane segments predicted by this new method and the known structure of the membrane-spanning helices of Rhodobacter reaction centers. From the above correlation it is proposed that the haem-coordinating domain of mitochondrial cytochrome b is folded in a novel structure, called "clepsydra domain", which is formed by distorted transmembrane helices packed in a waisted antiparallel bundle.     

Mapping epistatic quantitative trait loci with one-dimensional genome searches

The discovery of epistatically interacting QTL is hampered by the intractability and low power to detect QTL in multidimensional genome searches. We describe a new method that maps epistatic QTL by identifying loci of high QTL by genetic background interaction. This approach allows detection of QTL involved not only in pairwise but also higher-order interaction, and does so with one-dimensional genome searches. The approach requires large populations derived from multiple related inbred-line crosses as is more typically available for plants. Using maximum likelihood, the method contrasts models in which QTL allelic values are either nested within, or fixed over, populations. We apply the method to simulated doubled-haploid populations derived from a diallel among three inbred parents and illustrate the power of the method to detect QTL of different effect size and different levels of QTL by genetic background interaction. Further, we show how the method can be used in conjunction with standard two-locus QTL detection models that use two-dimensional genome searches and find that the method may double the power to detect first-order epistasis.     

Using Gaussian process with test rejection to detect T-cell epitopes in pathogen genomes

A major challenge in the development of peptide-based vaccines is finding the right immunogenic element, with efficient and long-lasting immunization effects, from large potential targets encoded by pathogen genomes. Computer models are convenient tools for scanning pathogen genomes to preselect candidate immunogenic peptides for experimental validation. Current methods predict many false positives resulting from a low prevalence of true positives. We develop a test reject method based on the prediction uncertainty estimates determined by Gaussian process regression. This method filters false positives among predicted epitopes from a pathogen genome. The performance of stand-alone Gaussian process regression is compared to other state-of-the-art methods using cross validation on 11 benchmark data sets. The results show that the Gaussian process method has the same accuracy as the top performing algorithms. The combination of Gaussian process regression with the proposed test reject method is used to detect true epitopes from the Vaccinia virus genome. The test rejection increases the prediction accuracy by reducing the number of false positives without sacrificing the method's sensitivity. We show that the Gaussian process in combination with test rejection is an effective method for prediction of T-cell epitopes in large and diverse pathogen genomes, where false positives are of concern.     

Preparation of single-stranded DNA from PCR products with streptavidin magnetic beads

The preparation of single-stranded DNA from double-stranded PCR products is an essential step in the identification of aptamers by Systematic Evolution of Ligands by EXponential enrichment (SELEX). The most widely used method for producing single-stranded DNA is alkaline denaturation of biotinylated PCR products attached to streptavidin-coated magnetic beads. Recently, it has been suggested that this method may be unsuitable due to the release of interfering amounts of streptavidin and biotinylated DNA. In this article, the alkaline method is compared with a thermal method that is known to release significant amounts of streptavidin and biotinylated DNA. Results show that trace amounts of streptavidin and biotinylated DNA are released in the alkaline method, but this can be curtailed by preconditioning the beads in aqueous sodium hydroxide. The main product in the alkaline method is single-stranded DNA, which is produced in high yield.     

Phylogenetic dating with confidence intervals using mean path lengths

The mean path length (MPL) method, a simple method for dating nodes in a phylogenetic tree, is presented. For small trees the age estimates and corresponding confidence intervals, calibrated with fossil data, can be calculated by hand, and for larger trees a computer program gives the results instantaneously (a Pascal program is available upon request). Necessary input data are a rooted phylogenetic tree with edge lengths (internode lengths) approximately corresponding to the number of substitutions between the nodes. Given this, the MPL method produces relative age estimates with confidence intervals for all nodes of the tree. With the age of one or several nodes of the tree being known from reference fossils, the relative age estimates induce absolute age estimates and confidence intervals of the nodes of the tree. The MPL method relies on the assumptions that substitutions occur randomly and independently in different sites in the DNA sequence and that the substitution rates are approximately constant in time, i.e., assuming a molecular clock. A method is presented for identification of the nodes in the tree at which significant deviations from the clock assumption occur, such that dating may be done using different rates in different parts of the tree. The MPL method is illustrated with the Liliales, a group of monocot flowering plants.     

Analysis of failure time data with dependent interval censoring

This article develops a method for the analysis of screening data for which the chance of being screened is dependent on the event of interest (informative censoring). Because not all subjects make all screening visits, the data on the failure of interest is interval censored. We propose a model that will properly adjust for the dependence to obtain an unbiased estimate of the nonparametric failure time function, and we provide an extension for applying the method for estimation of the regression parameters from a (discrete time) proportional hazards regression model. The method is applied on a data set from an observational study of cytomegalovirus shedding in a population of HIV-infected subjects who participated in a trial conducted by the AIDS Clinical Trials Group.     

Comparative evaluation of fibrinogen and fibrin hydrolysis with plasmin

A high-sensitive method is developed for determining the degree of plasmin-catalyzed fibrinogen hydrolysis by the released amino groups stained with trinitrobenzene sulphoacid. The method permits determining 0.02-0.08 casein units of plasmin. The method made it possible to establish that after streptokinase activation plasmin hydrolyzes equally fibrinogen and fibrin in solution and as gel. When a tissue activator is used, fibrin intensifies significantly the plasminogen activation. Inhibition of plasmin by an inhibitor produced from soya is considerably slowed down in fibrin gel.     

Quantitative analysis of sinistrin in serum with high-performance liquid chromatography for renal function testing

Sinistrin, as inulin, is widely used as a marker for renal function testing. A reliable and accurate method with a simple sample preparation for the quantitative determination of sinistrin would be of advantage. We developed a high-performance liquid chromatography (HPLC)-based method with electrochemical detection for the quantitative measurement of sinistrin in serum and plasma. Sample preparation is easy and includes enzymatic removal of glucose, deproteinization, acid hydrolysis of sinistrin, and HPLC separation of fructose. The recovery of sinistrin from serum is the same as that from water and is near 100%. The method presented has a linear range up to 500 mg/L. The results from our method are in agreement with a fully enzymatic quantification (regression coefficient 1.01, coefficient of variation 0.97). Sinistrin concentrations in aqueous solutions can be measured down to 2mg/L with a coefficient of variation of 5.7%. Quantification in serum is primarily limited by its physiological fructose content. The sensitivity of the described method is sufficient for its use in renal function testing. We describe a method for quantification of sinistrin which allows accurate measurements especially at low concentrations and low sample volumes. This laboratory method may be used for obtaining sinistrin pharmacokinetics in renal function testing.     

Revealing remote protein homology with sequence similarity and a modularity-based approach

An important task in functional genomics is to cluster homologous proteins, which may share common functions. Annotating proteins of unknown function by transferring annotations from their homologues of known annotations is one of the most efficient ways to predict protein function. In this paper, we use a modularity-based method called CD for grouping together homologous proteins. The method employs a global heuristic search strategy to find the partitioning of the weighted adjacency graph with the largest modularity. The weighted adjacency graph is constructed by the sigmodal transformation of all pairwise sequence similarities between all protein sequences in a given dataset. The method has been extensively tested on several subsets from the superfamily level of the SCOP (Structural Classification of Proteins) database, where some homologous proteins have very low sequence similarity. Compared with a widely used method MCL, we observe that the number of clusters obtained by CD is closer to the number of superfamilies in the dataset, the value of the F-measure given by CD is 10% better than MCL on average, and CD is more tolerant to noise to the sequence similarity. The experiment results indicate that CD is ideally suitable for clustering homologous proteins when sequence similarity is low.     

土壤宏转录组RNA的提取方法评价

【目的】比较手工法和试剂盒法提取土壤RNA对原核微生物多样性的影响,评价两种方法研究土壤宏转录组的技术偏好性。【方法】针对黑龙江海伦砂质黑壤土、江苏滨海黏心夹砂土、江西鹰潭第四纪红黏土不同母质发育形成的三种水稻土,采用手工和试剂盒法分别获得微生物总RNA,利用高通量测序原核微生物16S r RNA多样性,通过紫外分光和琼脂糖凝胶电泳分析RNA质量,比较手工和试剂盒法提取RNA对水稻土原核微生物群落的影响规律。【结果】三种水稻土中试剂盒提取RNA的纯度皆高于手工法,但RNA提取总量不一致。试剂盒提取黑龙江和江苏水稻土的RNA总量高于手工法,而手工法提取的江西水稻土RNA总量高于试剂盒法。高通量测序发现土壤类型而不是RNA提取方法决定了原核微生物多样性指数和群落结构。3个水稻土共检测到27门和409属,两种RNA提取方法偏好性提取了19门和181属,这些门和属占水稻土微生物总丰度平均值分别为40.4%和44.4%。与手工提取RNA相比,试剂盒法发现11个微生物门的丰度显著偏高(P0.05),但仅有Armatimonadetes门在三种水稻土同时存在专一偏好性;手工法提取也发现11个微生物门的丰度显著高于试剂盒法,并且仅有Firmicutes门在三种水稻土中同时存在专一偏好性。在微生物属水平,三种水稻土中均发现试剂盒法偏好性提取了2属,而手工法偏好性提取了5属。进一步针对72个优势微生物属(丰度0.1%且同时存在于三种水稻土),发现其占所有微生物丰度80%强,且48属的变化规律与RNA提取方法无关。例如,手工法提取水稻土好氧甲烷氧化菌的规律为:黑龙江(1.68%)江西(0.90%)江苏(0.59%),而试剂盒法也得到一致结果,黑龙江(0.52%)江西(0.18%)江苏(0.13%)。【结论】两种RNA提取方法本身的特异偏好性较小。三种水稻土27门409属中,仅有2门7属可能存在方法本身的专一偏好性,即在三种水稻土中,这些微生物均被试剂盒法或手工法特异性偏好提取,占所有原核微生物门和属比例仅为7%和1%左右。此外,针对同一水稻土,手工和试剂盒法提取RNA的总量、纯度、微生物相对丰度明显不同,在原核微生物门和属的分类水平,约70%的门和22%的属的丰度具有统计显著性差异,但两种RNA提http://journals.im.ac.cn/actamicro取方法均能反映优势微生物类群在三种水稻土中的变异规律,土壤类型对原核微生物多样性的影响远大于RNA提取方法本身的偏好性。未来原核微生物宏转录组研究中,应优先考虑科学问题及其实验处理可能导致的微生物组及其转录差异,结合微生物RNA提取特点,最大限度地发挥宏转录组技术优势。     

Gamete selection with an inbred tester

Summary A common inbred tester was used to evaluate gametes selected from three complex hybrid populations using two different inbreds as elite lines. The results support Stadler'S contention that gamete selection is an efficient method for extracting superior gene combinations from hybrid populations. The inbred tester increases the resolving power of the method by eliminating extraneous genetic variability and by providing a homogeneous check population for comparative purposes.Technical Article No. TA7973, Texas Agricultural Experiment Station, College Station, Texas. Project No. 1280. Parts of dissertations submitted by the first two authors in partial fulfillment of requirements for Ph. D. degrees.     

Detecting homology of distantly related proteins with consensus sequences

A simple protocol is described that is suitable for the detection of distantly related members of a protein family. In this procedure, similarity to a consensus sequence is used to distinguish chance similarity from similarity due to common ancestry. The consensus sequence is constructed from the sequences of established members of a protein family and it incorporates features characteristic of the protein fold of this family: conserved residues, the pattern of variable and conserved segments, preferred location of gaps etc. The database is searched with the consensus sequence, using the unitary matrix or log odds matrix for scoring the alignments, with variable gap penalty. The advantage of the method is that it weights key residues, ignores sequence similarity in variable segments (thus partially eliminating "background noise" coming from chance similarity), distinguishes gaps disrupting conserved segments from those occurring in positions known to be tolerant of gap events. The utility of the method was demonstrated in the case of the protein family homologous with the internal repeats of complement B as well as the internal repeats identified in fibroblast proteoglycan PG40. The consensus sequence method succeeded in finding some new members of these protein families that could not be detected by earlier methods of sequence comparison.     

Variational data assimilation with epidemic models

Mathematical modelling is playing an increasing role in developing an understanding of the dynamics of communicable disease and assisting the construction and implementation of intervention strategies. The threat of novel emergent pathogens in human and animal hosts implies the requirement for methods that can robustly estimate epidemiological parameters and provide forecasts. Here, a technique called variational data assimilation is introduced as a means of optimally melding dynamic epidemic models with epidemiological observations and data to provide forecasts and parameter estimates. Using data from a simulated epidemic process the method is used to estimate the start time of an epidemic, to provide a forecast of future epidemic behaviour and estimate the basic reproductive ratio. A feature of the method is that it uses a basic continuous-time SIR model, which is often the first point of departure for epidemiological modelling during the early stages of an outbreak. The method is illustrated by application to data gathered during an outbreak of influenza in a school environment.     

Complete Thymectomy in Adult Rats with Non-invasive Endotracheal Intubation

Thymectomy in neonatal rodents is an established and reliable procedure for immunological studies. However, in adult rats, complications of hemorrhage and pneumothorax from pleural disruption can result in a significant mortality rate. This protocol is a simple method of rat thymectomy that utilizes a mini-sternotomy and endotracheal intubation. Intubation is accomplished with a non-invasive and easily reproducible method and allows for positive pressure ventilation to prevent pneumothorax and a controlled airway that allows sufficient time for careful thymus dissection to minimize pleural disruption. A 1.5 cm sternal incision decreases contact with mediastinal vessels and pleura, while still providing full visualization of the thymus. Following exposure of the mediastinum, the thymus is removed by blunt dissection under magnification. The pleural space is then sealed by suture closure of the pre-tracheal muscles followed by the application of surgical glue. The thorax is then closed by suture closure of the sternum, followed by suture closure of the skin. All thymectomies were complete as evidenced by immunohistochemical (IHC) staining of mediastinal tissue, and absence of naïve T-cells by flow cytometry, and the procedure had a 96% survival rate. This method is suitable when complete thymectomy with minimal complications is desired for further immunological studies in athymic adult rats.     

A general method for dealing with misclassification in regression: the misclassification SIMEX

We have developed a new general approach for handling misclassification in discrete covariates or responses in regression models. The simulation and extrapolation (SIMEX) method, which was originally designed for handling additive covariate measurement error, is applied to the case of misclassification. The statistical model for characterizing misclassification is given by the transition matrix Pi from the true to the observed variable. We exploit the relationship between the size of misclassification and bias in estimating the parameters of interest. Assuming that Pi is known or can be estimated from validation data, we simulate data with higher misclassification and extrapolate back to the case of no misclassification. We show that our method is quite general and applicable to models with misclassified response and/or misclassified discrete regressors. In the case of a binary response with misclassification, we compare our method to the approach of Neuhaus, and to the matrix method of Morrissey and Spiegelman in the case of a misclassified binary regressor. We apply our method to a study on caries with a misclassified longitudinal response.     

A simple imputation method for longitudinal studies with non-ignorable non-responses

Missing data are a common problem in longitudinal studies in the health sciences. Motivated by data from the Muscatine Coronary Risk Factor (MCRF) study, a longitudinal study of obesity, we propose a simple imputation method for handling non-ignorable non-responses (i.e., when non-response is related to the specific values that should have been obtained) in longitudinal studies with either discrete or continuous outcomes. In the proposed approach, two regression models are specified; one for the marginal mean of the response, the other for the conditional mean of the response given non-response patterns. Statistical inference for the model parameters is based on the generalized estimating equations (GEE) approach. An appealing feature of the proposed method is that it can be readily implemented using existing, widely-available statistical software. The method is illustrated using longitudinal data on obesity from the MCRF study.