Well-Behavior, Well-Posedness and Nonsmooth Analysis

AMS subject classification: 90C30, 90C33.

Mixture models of nucleotide sequence evolution that account for heterogeneity in the substitution process across sites and across lineages

Molecular phylogenetic studies of homologous sequences of nucleotides often assume that the underlying evolutionary process was globally stationary, reversible, and homogeneous (SRH), and that a model of evolution with one or more site-specific and time-reversible rate matrices (e.g., the GTR rate matrix) is enough to accurately model the evolution of data over the whole tree. However, an increasing body of data suggests that evolution under these conditions is an exception, rather than the norm. To address this issue, several non-SRH models of molecular evolution have been proposed, but they either ignore heterogeneity in the substitution process across sites (HAS) or assume it can be modeled accurately using the distribution. As an alternative to these models of evolution, we introduce a family of mixture models that approximate HAS without the assumption of an underlying predefined statistical distribution. This family of mixture models is combined with non-SRH models of evolution that account for heterogeneity in the substitution process across lineages (HAL). We also present two algorithms for searching model space and identifying an optimal model of evolution that is less likely to over- or underparameterize the data. The performance of the two new algorithms was evaluated using alignments of nucleotides with 10 000 sites simulated under complex non-SRH conditions on a 25-tipped tree. The algorithms were found to be very successful, identifying the correct HAL model with a 75% success rate (the average success rate for assigning rate matrices to the tree's 48 edges was 99.25%) and, for the correct HAL model, identifying the correct HAS model with a 98% success rate. Finally, parameter estimates obtained under the correct HAL-HAS model were found to be accurate and precise. The merits of our new algorithms were illustrated with an analysis of 42 337 second codon sites extracted from a concatenation of 106 alignments of orthologous genes encoded by the nuclear genomes of Saccharomyces cerevisiae, S. paradoxus, S. mikatae, S. kudriavzevii, S. castellii, S. kluyveri, S. bayanus, and Candida albicans. Our results show that second codon sites in the ancestral genome of these species contained 49.1% invariable sites, 39.6% variable sites belonging to one rate category (V1), and 11.3% variable sites belonging to a second rate category (V2). The ancestral nucleotide content was found to differ markedly across these three sets of sites, and the evolutionary processes operating at the variable sites were found to be non-SRH and best modeled by a combination of eight edge-specific rate matrices (four for V1 and four for V2). The number of substitutions per site at the variable sites also differed markedly, with sites belonging to V1 evolving slower than those belonging to V2 along the lineages separating the seven species of Saccharomyces. Finally, sites belonging to V1 appeared to have ceased evolving along the lineages separating S. cerevisiae, S. paradoxus, S. mikatae, S. kudriavzevii, and S. bayanus, implying that they might have become so selectively constrained that they could be considered invariable sites in these species.

On Games on Non-Wellfounded Sets and Stationary Sets

In this thesis we study a few games related to non-wellfounded and stationary sets. Games have turned out to be an important tool in mathematical logic ranging from semantic games defining the truth of a sentence in a given logic to for example games on real numbers whose determinacies have important effects on the consistency of certain large cardinal assumptions. The equality of non-wellfounded sets can be determined by a so called bisimulation game already used to identify processes in theoretical computer science and possible world models for modal logic. Here we present a game to classify non-wellfounded sets according to their branching structure. We also study games on stationary sets moving back to classical wellfounded set theory. We also describe a way to approximate non-wellfounded sets with hereditarily finite wellfounded sets. The framework used to do this is domain theory. In the Banach-Mazur game, also called the ideal game, the players play a descending sequence of stationary sets and the second player tries to keep their intersection stationary. The game is connected to precipitousness of the corresponding ideal. In the pressing down game first player plays regressive functions defined on stationary sets and the second player responds with a stationary set where the function is constant trying to keep the intersection stationary. This game has applications in model theory to the determinacy of the Ehrenfeucht-Fraisse game. We show that it is consistent that these games are not equivalent.

Insight on how fishing bats discern prey and adjust their mechanic and sensorial features during the attack sequence

Several insectivorous bats have included fish in their diet, yet little is known about the processes underlying this trophic shift. We performed three field experiments with wild fishing bats to address how they manage to discern fish from insects and adapt their hunting technique to capture fish. We show that bats react only to targets protruding above the water and discern fish from insects based on prey disappearance patterns. Stationary fish trigger short and shallow dips and a terminal echolocation pattern with an important component of the narrowband and low frequency calls. When the fish disappears during the attack process, bats regulate their attack increasing the number of broadband and high frequency calls in the last phase of the echolocation as well as by lengthening and deepening their dips. These adjustments may allow bats to obtain more valuable sensorial information and to perform dips adjusted to the level of uncertainty on the location of the submerged prey. The observed ultrafast regulation may be essential for enabling fishing to become cost-effective in bats, and demonstrates the ability of bats to rapidly modify and synchronise their sensorial and motor features as a response to last minute stimulus variations.

Physiological changes in Campylobacter jejuni on entry into stationary phase

Campylobacter jejuni NCTC 11168 does not exhibit the general increase in cellular stress resistance on entry into stationary phase that is seen in most other bacteria. This is consistent with the lack of global stationary phase regulatory elements in this organism. deduced from an analysis of its genome sequence. We now show that C. jejuni NCTC 11168 does undergo certain changes in stationary phase, of a pattern not previously described. As cells entered stationary phase there was a change in membrane fatty acid composition, principally a decrease in the proportion of unsaturated fatty acids and an increase in the content of cyclopropane and short-chain fatty acids. These changes in membrane composition were accompanied by an increase in the resilience of the cell membrane towards loss of integrity caused by pressure and an increase in cellular pressure resistance. By contrast. there were no major changes in resistance to acid or heat treatment. A similar pattern of changes in stress resistance on entry, into stationary phase was seen in C. jejuni NCTC 11351, the type strain. These changes appear to represent a restricted Physiological response to the conditions existing in stationary phase cultures, in an organism having limited capacity for genetic regulation and adaptation to environment. © 2004 Elsevier B.V. All rights reserved.

CspC and CspD are essential for Caulobacter crescentus stationary phase survival

The cold shock response in bacteria involves the expression of low-molecular weight cold shock proteins (CSPs) containing a nucleic acid-binding cold shock domain (CSD), which are known to destabilize secondary structures on mRNAs, facilitating translation at low temperatures. Caulobacter crescentus cspA and cspB are induced upon cold shock, while cspC and cspD are induced during stationary phase. In this work, we determined a new coding sequence for the cspC gene, revealing that it encodes a protein containing two CSDs. The phenotypes of C. crescentus csp mutants were analyzed, and we found that cspC is important for cells to maintain viability during extended periods in stationary phase. Also, cspC and cspCD strains presented altered morphology, with frequent non-viable filamentous cells, and cspCD also showed a pronounced cell death at late stationary phase. In contrast, the cspAB mutant presented increased viability in this phase, which is accompanied by an altered expression of both cspC and cspD, but the triple cspABD mutant loses this characteristic. Taken together, our results suggest that there is a hierarchy of importance among the csp genes regarding stationary phase viability, which is probably achieved by a fine tune balance of the levels of these proteins.

SpdR, a Response Regulator Required for Stationary-Phase Induction of Caulobacter crescentus cspD

The cold shock protein (CSP) family includes small polypeptides that are induced upon temperature downshift and stationary phase. The genome of the alphaproteobacterium Caulobacter crescentus encodes four CSPs, with two being induced by cold shock and two at the onset of stationary phase. In order to identify the environmental signals and cell factors that are involved in cspD expression at stationary phase, we have analyzed cspD transcription during growth under several nutrient conditions. The results showed that expression of cspD was affected by the medium composition and was inversely proportional to the growth rate. The maximum levels of expression were decreased in a spoT mutant, indicating that ppGpp may be involved in the signalization for carbon starvation induction of cspD. A Tn5 mutant library was screened for mutants with reduced cspD expression, and 10 clones that showed at least a 50% reduction in expression were identified. Among these, a strain with a transposon insertion into a response regulator of a two-component system showed no induction of cspD at stationary phase. This protein (SpdR) was able to acquire a phosphate group from its cognate histidine kinase, and gel mobility shift assay and DNase I footprinting experiments showed that it binds to an inverted repeat sequence of the cspD regulatory region. A mutated SpdR with a substitution of the conserved aspartyl residue that is the probable phosphorylation site is unable to bind to the cspD regulatory region and to complement the spdR mutant phenotype.

Henon-like maps with arbitrary stationary combinatorics

We extend the renormalization operator introduced in [A. de Carvalho, M. Martens and M. Lyubich. Renormalization in the Henon family, I: universality but non-rigidity. J. Stat. Phys. 121(5/6) (2005), 611-669] from period-doubling Henon-like maps to Henon-like maps with arbitrary stationary combinatorics. We show that the renonnalization picture also holds in this case if the maps are taken to be strongly dissipative. We study infinitely renormalizable maps F and show that they have an invariant Cantor set O on which F acts like a p-adic adding machine for some p > 1. We then show, as for the period-doubling case in the work of de Carvalho, Martens and Lyubich [Renormalization in the Henon family, I: universality but non-rigidity. J. Stat. Phys. 121(5/6) (2005), 611-669], that the sequence of renormalizations has a universal form, but that the invariant Cantor set O is non-rigid. We also show that O cannot possess a continuous invariant line field.

The putative Leishmania telomerase RNA (leishTER) undergoes trans-splicing and contains a conserved template sequence

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

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. ^

Fine-mapping and mutation analysis of TRPM1: a candidate gene for leopard complex (LP) spotting and congenital stationary night blindness in horses.

Leopard Complex spotting occurs in several breeds of horses and is caused by an incompletely dominant allele (LP). Homozygosity for LP is also associated with congenital stationary night blindness (CSNB) in Appaloosa horses. Previously, LP was mapped to a 6 cm region on ECA1 containing the candidate gene TRPM1 (Transient Receptor Potential Cation Channel, Subfamily M, Member 1) and decreased expression of this gene, measured by qRT-PCR, was identified as the likely cause of both spotting and ocular phenotypes. This study describes investigations for a mutation causing or associated with the Leopard Complex and CSNB phenotype in horses. Re-sequencing of the gene and associated splice sites within the 105 624 bp genomic region of TRPM1 led to the discovery of 18 SNPs. Most of the SNPs did not have a predictive value for the presence of LP. However, one SNP (ECA1:108,249,293 C>T) found within intron 11 had a strong (P < 0.0005), but not complete, association with LP and CSNB and thus is a good marker but unlikely to be causative. To further localize the association, 70 SNPs spanning over two Mb including the TRPM1 gene were genotyped in 192 horses from three different breeds segregating for LP. A single 173 kb haplotype associated with LP and CSNB (ECA1: 108,197,355- 108,370,150) was identified. Illumina sequencing of 300 kb surrounding this haplotype revealed 57 SNP variants. Based on their localization within expressed sequences or regions of high sequence conservation across mammals, six of these SNPs were considered to be the most likely candidate mutations. While the precise function of TRPM1 remains to be elucidated, this work solidifies its functional role in both pigmentation and night vision. Further, this work has identified several potential regulatory elements of the TRPM1 gene that should be investigated further in this and other species.

Estimation of evolutionary distances under stationary and nonstationary models of nucleotide substitution

Estimation of evolutionary distances has always been a major issue in the study of molecular evolution because evolutionary distances are required for estimating the rate of evolution in a gene, the divergence dates between genes or organisms, and the relationships among genes or organisms. Other closely related issues are the estimation of the pattern of nucleotide substitution, the estimation of the degree of rate variation among sites in a DNA sequence, and statistical testing of the molecular clock hypothesis. Mathematical treatments of these problems are considerably simplified by the assumption of a stationary process in which the nucleotide compositions of the sequences under study have remained approximately constant over time, and there now exist fairly extensive studies of stationary models of nucleotide substitution, although some problems remain to be solved. Nonstationary models are much more complex, but significant progress has been recently made by the development of the paralinear and LogDet distances. This paper reviews recent studies on the above issues and reports results on correcting the estimation bias of evolutionary distances, the estimation of the pattern of nucleotide substitution, and the estimation of rate variation among the sites in a sequence.

Cortical activity flips among quasi-stationary states.

Parallel recordings of spike trains of several single cortical neurons in behaving monkeys were analyzed as a hidden Markov process. The parallel spike trains were considered as a multivariate Poisson process whose vector firing rates change with time. As a consequence of this approach, the complete recording can be segmented into a sequence of a few statistically discriminated hidden states, whose dynamics are modeled as a first-order Markov chain. The biological validity and benefits of this approach were examined in several independent ways: (i) the statistical consistency of the segmentation and its correspondence to the behavior of the animals; (ii) direct measurement of the collective flips of activity, obtained by the model; and (iii) the relation between the segmentation and the pair-wise short-term cross-correlations between the recorded spike trains. Comparison with surrogate data was also carried out for each of the above examinations to assure their significance. Our results indicated the existence of well-separated states of activity, within which the firing rates were approximately stationary. With our present data we could reliably discriminate six to eight such states. The transitions between states were fast and were associated with concomitant changes of firing rates of several neurons. Different behavioral modes and stimuli were consistently reflected by different states of neural activity. Moreover, the pair-wise correlations between neurons varied considerably between the different states, supporting the hypothesis that these distinct states were brought about by the cooperative action of many neurons.

Functional Transfer Theorems for Maxima of Stationary Processes

2000 Mathematics Subject Classification: 60G70, 60F12, 60G10.