Models of DNA sequence evolution

Models of DNA sequence evolution and methods for estimating evolutionary distances are needed for studying the rate and pattern of molecular evolution and for inferring the evolutionary relationships of organisms or genes. In this dissertation, several new models and methods are developed.^ The rate variation among nucleotide sites: To obtain unbiased estimates of evolutionary distances, the rate heterogeneity among nucleotide sites of a gene should be considered. Commonly, it is assumed that the substitution rate varies among sites according to a gamma distribution (gamma model) or, more generally, an invariant+gamma model which includes some invariable sites. A maximum likelihood (ML) approach was developed for estimating the shape parameter of the gamma distribution $(\alpha)$ and/or the proportion of invariable sites $(\theta).$ Computer simulation showed that (1) under the gamma model, $\alpha$ can be well estimated from 3 or 4 sequences if the sequence length is long; and (2) the distance estimate is unbiased and robust against violations of the assumptions of the invariant+gamma model.^ However, this ML method requires a huge amount of computational time and is useful only for less than 6 sequences. Therefore, I developed a fast method for estimating $\alpha,$ which is easy to implement and requires no knowledge of tree. A computer program was developed for estimating $\alpha$ and evolutionary distances, which can handle the number of sequences as large as 30.^ Evolutionary distances under the stationary, time-reversible (SR) model: The SR model is a general model of nucleotide substitution, which assumes (i) stationary nucleotide frequencies and (ii) time-reversibility. It can be extended to SRV model which allows rate variation among sites. I developed a method for estimating the distance under the SR or SRV model, as well as the variance-covariance matrix of distances. Computer simulation showed that the SR method is better than a simpler method when the sequence length $L>1,000$ bp and is robust against deviations from time-reversibility. As expected, when the rate varies among sites, the SRV method is much better than the SR method.^ The evolutionary distances under nonstationary nucleotide frequencies: The statistical properties of the paralinear and LogDet distances under nonstationary nucleotide frequencies were studied. First, I developed formulas for correcting the estimation biases of the paralinear and LogDet distances. The performances of these formulas and the formulas for sampling variances were examined by computer simulation. Second, I developed a method for estimating the variance-covariance matrix of the paralinear distance, so that statistical tests of phylogenies can be conducted when the nucleotide frequencies are nonstationary. Third, a new method for testing the molecular clock hypothesis was developed in the nonstationary case. ^

Biomarker circadian rhythm profiles in critically ill mechanically ventilated patients

Objective: To explore the natural trajectory of core body temperature (CBT) and cortisol (CORT) circadian rhythms in mechanically ventilated intensive care unit (MV ICU) patients. ^ Design: Prospective, observational, time-series pilot study. ^ Setting: Medical-surgical and pulmonary ICUs in a tertiary care hospital. ^ Sample: Nine (F = 3, M = 6) adults who were mechanically ventilated within 12 hrs of ICU admission with mean ± SD age of 65.2 ± 14 years old. ^ Measurements: Core body temperature and environmental measures of light, sound, temperature, and relative humidity were logged in 1-min intervals. Hourly urine specimens and 2-hr interval blood specimens were collected for up to 7 consecutive days for CORT assay. Mechanical ventilation days, ICU length of stay, and ICU mortality were documented. Acute Physiology and Chronic Health Evaluation (APACHE) II scores were computed for each study day. The data of each biologic and environmental variable were analyzed using single cosinor analysis of 24-hr serial segments. One patient did not complete the study because mortality occurred within 8 hrs of enrollment. Nine ICU patients completed the study in 1.6 to 7.0 days. ^ Results: No normal circadian rhythm pattern was found when the cosinor-derived parameters of amplitude (one-half the peak-trough variability) and acrophase (peak time) were compared with cosinor-derived parameter reference ranges of healthy, diurnally active humans, although 83% of patient-day CBT segments showed statistically significant (p ≤ .05) and biologically meaningful (R2≥ 0.30) 24-hr rhythms with abnormal cosinor parameters. Cosinor parameters of the environmental temporal profiles showed 27% of light, 76% of ambient temperature, and 78% of relative humidity serial segments had a significant and meaningful 24-hr diurnal pattern. Average daily light intensity varied from 34 to 187 lx with a maximum light exposure of 1877 lx. No sound measurement segment had a statistically significant cosine pattern, and numerous 1-minute interval peaks ≥ 60 dB occurred around the clock. Average daily ambient temperature and relative humidity varied from 19 to 24°C and from 25% to 61%, respectively. There was no statistically significant association between CBT or clinical outcomes and cosinor-derived parameters of the environmental variables. Circadian rhythms of urine and plasma CORT were deferred for later analysis. ^ Conclusions: The natural trajectory of the CBT circadian rhythm in MV ICU patients demonstrated persistent cosinor parameter alteration, even when a significant and meaningful 24-hr rhythm was present. The ICU environmental measures showed erratic light and sound exposures. Room temperature and relative humidity data produced the highest rate of significant and meaningful diurnal 24-hr patterns. Additional research is needed to clarify relations among the CBT biomarker of the circadian clock and environmental variables of MV ICU patients. ^

MOLECULAR AND GENOMIC BASED INSIGHTS INTO THE EVOLUTION OF ENTEROCOCCUS FAECIUM FROM COMMENSAL TO HOSPITAL-ADAPTED PATHOGEN

The basis for the recent transition of Enterococcus faecium from a primarily commensal organism to one of the leading causes of hospital-acquired infections in the United States is not yet understood. To address this, the first part of my project assessed isolates from early outbreaks in the USA and South America using sequence analysis, colony hybridizations, and minimal inhibitory concentrations (MICs) which showed clinical isolates possess virulence and antibiotic resistance determinants that are less abundant or lacking in community isolates. I also revealed that the level of ampicillin resistance increased over time in clinical strains. By sequencing the pbp5 gene, I demonstrated an ~5% difference in the pbp5 gene between strains with MICs <4ug/ml and those with MICs >4µg/ml, but no specific sequence changes correlated with increases in MICs within the latter group. A 3-10% nucleotide difference was also seen in three other genes analyzed, which suggested the existence of two distinct subpopulations of E. faecium. This led to the second part of my project analyzing concatenated core gene sequences, SNPs, the 16S rRNA, and phylogenetics of 21 E. faecium genomes confirming two distinct clades; a community-associated (CA) clade and hospital-associated (HA) clade. Molecular clock calculations indicate that these two clades likely diverged ~ 300,000 to > 1 million years ago, long before the modern antibiotic era. Genomic analysis also showed that, in addition to core genomic differences, HA E. faecium harbor specific accessory genetic elements that may confer selection advantages over CA E. faecium. The third part of my project discovered 6 E. faecium genes with the newly identified “WxL” domain. My analyses, using RT-PCR, western blots, patient sera, whole-cell ELISA, and immunogold electron microscopy, indicated that E. faecium WxL genes exist in operons, encode bacterial cell surface localized proteins, that WxL proteins are antigenic in humans, and are more exposed on the surface of clinical isolates versus community isolates (even though they are ubiquitous in both clades). ELISAs and BIAcore analyses also showed that proteins encoded by these operons bind several different host extracellular matrix proteins, as well as to each other, suggesting a novel cell-surface complex. In summary, my studies provide new insights into the evolution of E. faecium by showing that there are two distantly related clades; one being more successful in the hospital setting. My studies also identified operons encoding WxL proteins whose characteristics could also contribute to colonization and virulence within this species.