Cubic-spline interpolation to estimate effects of inbreeding on milk yield in first lactation Holstein cows

Milk yield records (305d, 2X, actual milk yield) of 123,639 registered first lactation Holstein cows were used to compare linear regression (y = β(0) + β(1)X + e), quadratic regression, (y = β(0) + β(1)X + β(2)X(2) + e) cubic regression (y = β(0) + β(1)X + β(2)X(2) + β(3)X(3) +e) and fixed factor models, with cubic-spline interpolation models, for estimating the effects of inbreeding on milk yield. Ten animal models, all with herd-year-season of calving as fixed effect, were compared using the Akaike corrected-Information Criterion (AICc). The cubic-spline interpolation model with seven knots had the lowest AICc, whereas for all those labeled as "traditional", AICc was higher than the best model. Results from fitting inbreeding using a cubic-spline with seven knots were compared to results from fitting inbreeding as a linear covariate or as a fixed factor with seven levels. Estimates of inbreeding effects were not significantly different between the cubic-spline model and the fixed factor model, but were significantly different from the linear regression model. Milk yield decreased significantly at inbreeding levels greater than 9%. Variance component estimates were similar for the three models. Ranking of the top 100 sires with daughter records remained unaffected by the model used.     

Fire,landscape change and models of small mammal habitat suitability at multiple spatial scales

Fire is an important process in many ecosystems, but inappropriate fire regimes can adversely affect biodiversity. We identified a naturally flammable heathy woodland ecosystem where the use of planned fire had increased the extent of older vegetation, and quantified the abundance of two small native mammals in this landscape (silky mouse Pseudomys apodemoides and heath rat P. shortridgei). We defined four time‐since‐fire (TSF) categories representing a 2‐ to 55‐year post‐fire sequence and, on the basis of a habitat accommodation model, predicted that both species would select younger age‐classes over older ones. We also predicted that (i) much of the variance in vegetation structure would remain unexplained by TSF and (ii) statistical models of mammal abundance and occupancy including structural variables as predictors would be better than models including TSF. Pseudomys apodemoides selected 17‐ to 23‐year‐old sites, while there was no evidence that P. shortridgei selected a particular TSF category, findings that were inconsistent with our predictions. In line with our predictions, relatively large portions of the variance in vegetation structure remained unexplained by TSF (adjustedr² for four structural variables: 0.24, 0.29, 0.35 and 0.57), and in three of four cases there was strong evidence that statistical models of mammal abundance and occupancy including structural variables were better than those including TSF. At the site scale (hectares), P. shortridgei abundance was positively related to the cover of dead material at the base of Xanthorrhoea plants and at the trap scale (metres), the trapability of both species was significantly related to vegetation volume at 0–20 cm. Our findings suggest that TSF may not be a good proxy for either vegetation structure or species abundance/occupancy.     

Homology and ontogeny: pattern and process in comparative developmental biology

In this article the interface between development and homology is discussed. Development is here interpreted as a sequence of evolutionarily independent stages. Any approach stressing the importance of specific developmental stages is rejected. A homology definition is favoured which includes similarity, and complexity serves as a test for homology. Complexity is seen as the possibility of subdividing a character into evolutionarily independent corresponding substructures. Topology as a test for homology is critically discussed because corresponding positions are not necessarily indicative of homology. Complexity can be used twofold for homology assessments of development: either stages or processes of development are homologised. These two approaches must not be conflated. This distinction leads to the conclusion that there is no ontogenetic homology “criterion”.     

Parametric spatial cure rate models for interval-censored time-to-relapse data

Several recent papers (e.g., Chen, Ibrahim, and Sinha, 1999, Journal of the American Statistical Association 94, 909-919; Ibrahim, Chen, and Sinha, 2001a, Biometrics 57, 383-388) have described statistical methods for use with time-to-event data featuring a surviving fraction (i.e., a proportion of the population that never experiences the event). Such cure rate models and their multivariate generalizations are quite useful in studies of multiple diseases to which an individual may never succumb, or from which an individual may reasonably be expected to recover following treatment (e.g., various types of cancer). In this article we extend these models to allow for spatial correlation (estimable via zip code identifiers for the subjects) as well as interval censoring. Our approach is Bayesian, where posterior summaries are obtained via a hybrid Markov chain Monte Carlo algorithm. We compare across a broad collection of rather high-dimensional hierarchical models using the deviance information criterion, a tool recently developed for just this purpose. We apply our approach to the analysis of a smoking cessation study where the subjects reside in 53 southeastern Minnesota zip codes. In addition to the usual posterior estimates, our approach yields smoothed zip code level maps of model parameters related to the relapse rates over time and the ultimate proportion of quitters (the cure rates).     

Value of heartburn for diagnosing gastroesophageal reflux disease in severely obese patients

Objective : To evaluate the prevalence of gastroesophageal reflux disease (GERD) in severely obese patients and the association between symptoms and objective data of GERD in this population. Research Methods and Procedures : A total of 158 consecutive severely obese patients (BMI ≥ 40 kg/m²) were prospectively evaluated. Symptoms were evaluated by a structured clinical questionnaire. Objective assessment was made by ambulatory 24‐hour esophageal pH monitoring and endoscopy. GERD was defined by the presence of symptoms or complications (esophagitis). The clinical criterion defining GERD was the presence of at least two episodes of heartburn per week. Results : The mean age of the 138 patients subjected to complete study was 42.6 ± 10.2 years, with a BMI of 50.1 ± 6.9 kg/m² (range, 40.6 to 69.4 kg/m²); 78% were women. The prevalence of GERD evaluated by symptoms and/or esophagitis was 33.3% (46/138). Clinical criteria of GERD were present in 31/138 cases (22.5%), and 26 (18.8%) had esophagitis. In 69/138 patients (50%), pHmetry was abnormal. Fifty‐three patients with esophagitis and/or abnormal pHmetry were asymptomatic. The sensitivity of heartburn as a diagnostic criterion of GERD in patients with severe obesity was 29.3%, with a specificity of 85.7%. No significant association was observed between severe obesity grade and the prevalence of symptoms and/or objective data. Discussion : Asymptomatic gastroesophageal reflux (abnormal esophageal acid exposure and/or reflux esophagitis) is more common than symptomatic gastroesophageal reflux in severely obese patients. Increased BMI is not associated with a greater prevalence of GERD in these patients.     

Variable selection for logistic regression using a prediction-focused information criterion

In biostatistical practice, it is common to use information criteria as a guide for model selection. We propose new versions of the focused information criterion (FIC) for variable selection in logistic regression. The FIC gives, depending on the quantity to be estimated, possibly different sets of selected variables. The standard version of the FIC measures the mean squared error of the estimator of the quantity of interest in the selected model. In this article, we propose more general versions of the FIC, allowing other risk measures such as the one based on L(p) error. When prediction of an event is important, as is often the case in medical applications, we construct an FIC using the error rate as a natural risk measure. The advantages of using an information criterion which depends on both the quantity of interest and the selected risk measure are illustrated by means of a simulation study and application to a study on diabetic retinopathy.     

Seemingly unrelated measurement error models, with application to nutritional epidemiology

Motivated by an important biomarker study in nutritional epidemiology, we consider the combination of the linear mixed measurement error model and the linear seemingly unrelated regression model, hence Seemingly Unrelated Measurement Error Models. In our context, we have data on protein intake and energy (caloric) intake from both a food frequency questionnaire (FFQ) and a biomarker, and wish to understand the measurement error properties of the FFQ for each nutrient. Our idea is to develop separate marginal mixed measurement error models for each nutrient, and then combine them into a larger multivariate measurement error model: the two measurement error models are seemingly unrelated because they concern different nutrients, but aspects of each model are highly correlated. As in any seemingly unrelated regression context, the hope is to achieve gains in statistical efficiency compared to fitting each model separately. We show that if we employ a "full" model (fully parameterized), the combination of the two measurement error models leads to no gain over considering each model separately. However, there is also a scientifically motivated "reduced" model that sets certain parameters in the "full" model equal to zero, and for which the combination of the two measurement error models leads to considerable gain over considering each model separately, e.g., 40% decrease in standard errors. We use the Akaike information criterion to distinguish between the two possibilities, and show that the resulting estimates achieve major gains in efficiency. We also describe theoretical and serious practical problems with the Bayes information criterion in this context.     

Structure-based phylogeny of polyene macrolide antibiotic glycosyltransferases

Antibiotic glycosyltransferases (AGts) attach unusual deoxy-sugars to aglycons so antibiotics can exert function. It has been reported that polyene macrolide (PEM) AGts have different evolutionary origin when compared with other polyketide AGts, and our previous analysis have suggested that they could be results of horizontal gene transfer (HGT) from eukaryotes. In this paper, we compared the structures of PEM AGts with structures of eukaryotes and other AGts, and then built models of the representative PEM AGts and GT-1 glycosyltransferases. We also constructed the Neighbor-Joining (NJ) trees based on the normalized Root Mean Square (RMS) distance, the Bayesian tree guided by structural alignments, and carried out analysis on several key conserved residues in PEM AGts. The NJ tree showed a close relationship between PEM AGts and eukaryotic glycosyltransferases, and Bayesian tree further supported their affinity with UDP-glucuronosyltransferases (UGTs). Analysis on key conserved residues showed that PEM AGts may have similar interaction mechanism such as in the formation of hydrogen bonds as eukaryotic glycosyltransferases. Using structure-based phylogenetic approaches, this study further supported that PEM AGts were the result of HGT between prokaryotes and eukaryotes.     

Systematic investigation of Amphioxus (Branchiostoma floridae) microRNAs

Ascaris lumbricoides and Ascaris suum are parasitic nematodes living in the small intestine of humans and pigs, and can cause the disease ascariasis. For long, there has been controversy as to whether the two ascaridoid taxa represent the same species due to their significant resemblances in morphology. However, the complete mitochondrial (mt) genome data have been lacking for A. lumbricoides in spite of human and animal health significance and socio-economic impact globally of these parasites. In the present study, we sequenced the complete mt genomes of A. lumbricoides and A. suum (China isolate), which was 14,303 bp and 14,311 bp in size, respectively. The identity of the mt genomes was 98.1% between A. lumbricoides and A. suum (China isolate), and 98.5% between A. suum (China isolate) and A. suum (USA isolate). Both genomes are circular, and consist of 36 genes, including 12 genes for proteins, 2 genes for rRNA and 22 genes for tRNA, which are consistent with that of all other species of ascaridoid studied to date. All genes are transcribed in the same direction and have a nucleotide composition high in A and T (71.7% for A. lumbricoides and 71.8% for A. suum). The AT bias had a significant effect on both the codon usage pattern and amino acid composition of proteins. Phylogenetic analyses of A. lumbricoides and A. suum using concatenated amino acid sequences of 12 protein-coding genes, with three different computational algorithms (Bayesian analysis, maximum likelihood and maximum parsimony) all clustered in a clade with high statistical support, indicating that A. lumbricoides and A. suum was very closely related. These mt genome data and the results provide some additional genetic evidence that A. lumbricoides and A. suum may represent the same species. The mt genome data presented in this study are also useful novel markers for studying the molecular epidemiology and population genetics of Ascaris.