Testing the relationship between morphological and molecular rates of change along phylogenies

Molecular evolution has been considered to be essentially a stochastic process, little influenced by the pace of phenotypic change. This assumption was challenged by a study that demonstrated an association between rates of morphological and molecular change estimated for total-evidence phylogenies, a finding that led some researchers to challenge molecular date estimates of major evolutionary radiations. Here we show that Omland's (1997) result is probably due to methodological bias, particularly phylogenetic nonindependence, rather than being indicative of an underlying evolutionary phenomenon. We apply three new methods specifically designed to overcome phylogenetic bias to 13 published phylogenetic datasets for vertebrate taxa, each of which includes both morphological characters and DNA sequence data. We find no evidence of an association between rates of molecular and morphological rates of change.

A generalized mechanistic codon model.

Models of codon evolution have attracted particular interest because of their unique capabilities to detect selection forces and their high fit when applied to sequence evolution. We described here a novel approach for modeling codon evolution, which is based on Kronecker product of matrices. The 61 × 61 codon substitution rate matrix is created using Kronecker product of three 4 × 4 nucleotide substitution matrices, the equilibrium frequency of codons, and the selection rate parameter. The entities of the nucleotide substitution matrices and selection rate are considered as parameters of the model, which are optimized by maximum likelihood. Our fully mechanistic model allows the instantaneous substitution matrix between codons to be fully estimated with only 19 parameters instead of 3,721, by using the biological interdependence existing between positions within codons. We illustrate the properties of our models using computer simulations and assessed its relevance by comparing the AICc measures of our model and other models of codon evolution on simulations and a large range of empirical data sets. We show that our model fits most biological data better compared with the current codon models. Furthermore, the parameters in our model can be interpreted in a similar way as the exchangeability rates found in empirical codon models.

Evaluation of Fly Ash Concrete Durability Containing Class II Durability Aggregates, July 1986

Fly ash was used in this evaluation study to replace 15% of the cement in Class C-3 concrete paving mixes. One Class "c" ash from Iowa approved sources was examined in each mix. Substitution rate was based on 1 to 1 basis, for each pound of cement removed 1.0 pound of ash was added. The freeze/thaw durability of the concrete studied was not adversely affected by the presence of fly ash. This study reveals that the durability of the concrete test specimens made with Class II durability aggregates was slightly increased in all cases by the substitution of cement with 15% Class "c" fly ash. In all cases durability factors either remained the same or slightly improved except for one case where the durability factor decreased from 36 to 34. The expansion decreased in all cases.

Structural and magnetic characterization of maghemites prepared from Al-substituted magnetites

Synthetic aluminum-substituted maghemites were characterized by total chemical analysis, powder X-ray diffraction (XRD), Mössbauer spectroscopy (ME), and vibrating sample magnetometry (VSM). The aim was to determine the structural, magnetic, and hyperfine properties of γ-Fe2-xAl xO3 as the Al concentration is varied. The XRD results of the synthetic products were indexed exclusively as maghemite. Increasing Al for Fe substitution decreased the mean crystalline dimension and shifted all diffraction peaks to higher º2θ angles. The a0 dimension of the cubic unit cell decreased with increasing Al according to the equation a o = 0.8385 - 3.63 x 10-5 Al (R²= 0.94). Most Mössbauer spectra were composed of one sextet, but at the highest substitution rate of 142.5 mmol mol-1 Al, both a doublet and sextet were obtained at 300 K. All hyperfine parameters from the sub-spectra were consistent with high-spin Fe3+ (0.2 a 0.7 mms-1) and suggested a strong superparamagnetic component associated with the doublet. The magnetic hyperfine field of the sextets decreased with the amount of Al-substitution [Bhf (T) = 49.751 - 0.1202Al; R² = 0.94] while the linewidth increased linearly. The saturation magnetization also decreased with increasing isomorphous substitution.

Evaluation of Fly Ash Concrete Durability Containing Class II Durability Aggregates, July 1986

Fly ash was used in this evaluation study to replace 15% of the cement in Class C-3 concrete paving mixes. One Class "c" ash from Iowa approved sources was examined in each mix. Substitution rate was based on 1 to 1 basis, for each pound of cement removed 1.0 pound of ash was added. The freeze/thaw durability of the concrete studied was not adversely affected by the presence of fly ash. This study reveals that the durability of the concrete test specimens made with Class II durability aggregates was slightly increased in all cases by the substitution of cement with 15% Class "c" fly ash. In all cases durability factors either remained the same or slightly improved except for one case where the durability factor decreased from 36 to 34. The expansion decreased in all cases.

MODELING HETEROTACHY IN PHYLOGENETICS

Il a été démontré que l’hétérotachie, variation du taux de substitutions au cours du temps et entre les sites, est un phénomène fréquent au sein de données réelles. Échouer à modéliser l’hétérotachie peut potentiellement causer des artéfacts phylogénétiques. Actuellement, plusieurs modèles traitent l’hétérotachie : le modèle à mélange des longueurs de branche (MLB) ainsi que diverses formes du modèle covarion. Dans ce projet, notre but est de trouver un modèle qui prenne efficacement en compte les signaux hétérotaches présents dans les données, et ainsi améliorer l’inférence phylogénétique. Pour parvenir à nos fins, deux études ont été réalisées. Dans la première, nous comparons le modèle MLB avec le modèle covarion et le modèle homogène grâce aux test AIC et BIC, ainsi que par validation croisée. A partir de nos résultats, nous pouvons conclure que le modèle MLB n’est pas nécessaire pour les sites dont les longueurs de branche diffèrent sur l’ensemble de l’arbre, car, dans les données réelles, le signaux hétérotaches qui interfèrent avec l’inférence phylogénétique sont généralement concentrés dans une zone limitée de l’arbre. Dans la seconde étude, nous relaxons l’hypothèse que le modèle covarion est homogène entre les sites, et développons un modèle à mélanges basé sur un processus de Dirichlet. Afin d’évaluer différents modèles hétérogènes, nous définissons plusieurs tests de non-conformité par échantillonnage postérieur prédictif pour étudier divers aspects de l’évolution moléculaire à partir de cartographies stochastiques. Ces tests montrent que le modèle à mélanges covarion utilisé avec une loi gamma est capable de refléter adéquatement les variations de substitutions tant à l’intérieur d’un site qu’entre les sites. Notre recherche permet de décrire de façon détaillée l’hétérotachie dans des données réelles et donne des pistes à suivre pour de futurs modèles hétérotaches. Les tests de non conformité par échantillonnage postérieur prédictif fournissent des outils de diagnostic pour évaluer les modèles en détails. De plus, nos deux études révèlent la non spécificité des modèles hétérogènes et, en conséquence, la présence d’interactions entre différents modèles hétérogènes. Nos études suggèrent fortement que les données contiennent différents caractères hétérogènes qui devraient être pris en compte simultanément dans les analyses phylogénétiques.

Estimation des corrélations phylogénétiques entre paramètres d'évolution moléculaire et Traits d'histoire de vie

Depuis quelques années, l'évolution moléculaire cherche à caractériser les variations et l'intensité de la sélection grâce au rapport entre taux de substitution synonyme et taux de substitution non-synonyme (dN/dS). Cette mesure, dN/dS, a permis d'étudier l'histoire de la variation de l'intensité de la sélection au cours du temps ou de détecter des épisodes de la sélection positive. Les liens entre sélection et variation de taille efficace interfèrent cependant dans ces mesures. Les méthodes comparatives, quant a elle, permettent de mesurer les corrélations entre caractères quantitatifs le long d'une phylogénie. Elles sont également utilisées pour tester des hypothèses sur l'évolution corrélée des traits d'histoire de vie, mais pour être employées pour étudier les corrélations entre traits d'histoire de vie, masse, taux de substitution ou dN/dS. Nous proposons ici une approche combinant une méthode comparative basée sur le principe des contrastes indépendants et un modèle d'évolution moléculaire, dans un cadre probabiliste Bayésien. Intégrant, le long d'une phylogénie, sur les reconstructions ancestrales des traits et et de dN/dS nous estimons les covariances entre traits ainsi qu'entre traits et paramètres du modèle d'évolution moléculaire. Un modèle hiérarchique, a été implémenté dans le cadre du logiciel coevol, publié au cours de cette maitrise. Ce modèle permet l'analyse simultané de plusieurs gènes sans perdre la puissance donnée par l'ensemble de séquences. Un travail deparallélisation des calculs donne la liberté d'augmenter la taille du modèle jusqu'à l'échelle du génome. Nous étudions ici les placentaires, pour lesquels beaucoup de génomes complets et de mesures phénotypiques sont disponibles. À la lumière des théories sur les traits d'histoire de vie, notre méthode devrait permettre de caractériser l'implication de groupes de gènes dans les processus biologique liés aux phénotypes étudiés.

Cytochrome-b variation in Apis mellifera samples and its association with COI-COII patterns

Five Mbo I (Mbo-A, Mbo-M, Mbo-C(1), Mbo-C(2) and Mbo-C(3)) and Hinf I (Hinf-1 to Hinf-5) patterns were observed in Apis mellifera samples after restriction of a 485 bp fragment of the mitochondrial cytochrome-b (cyt-b) gene. Associating the cyt-b Restriction fragment length polymorphism (RFLP) pattern of each sample to its respective previously established COI-COII (Dra I sites) pattern, five restriction patterns (Mbo-C(1), Mbo-C(2), Mbo-C(3), Hinf-1 and Hinf-4) were observed in samples of maternal origin associated to the evolutionary branch C. No deletions or insertions were observed and the nucleotide substitution rate was estimated at 5.4%. Higher nucleotide diversity was observed among the branch C-haplotypes when compared with A and M lineages. Further studies are needed to confirm if the cyt-b + COI-COII haplotypes help to assign certain phylogeographic patterns to the branch C and to clarify phylogenetic relationships among A. mellifera subspecies.

Genus Paracoccidioides: Species Recognition and Biogeographic Aspects

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

Spin reversal and magnetization jumps in ErMexMn1-xO3 perovskites (Me = Ni, Co)

The erbium-based manganite ErMnO3 has been partially substituted at the manganese site by the transition-metal elements Ni and Co. The perovskite orthorhombic structure is found from x(Ni) = 0.2-0.5 in the nickel-based solid solution ErNixMn1-xO3, while it can be extended up to x(Co) = 0.7 in the case of cobalt, provided that the synthesis is performed under oxygenation conditions to favor the presence of Co3+. Presence of different magnetic entities (i.e., Er3+, Ni2+, Co2+, Co3+, Mn3+, and Mn4+) leads to quite unusual magnetic properties, characterized by the coexistence of antiferromagnetic and ferromagnetic interactions. In ErNixMn1-xO3, a critical concentration x(crit)(Ni) = 1/3 separates two regimes: spin-canted AF interactions predominate at x < x(crit), while the ferromagnetic behavior is enhanced for x > x(crit). Spin reversal phenomena are present both in the nickel- and cobalt-based compounds. A phenomenological model based on two interacting sublattices, coupled by an antiferromagnetic exchange interaction, explains the inversion of the overall magnetic moment at low temperatures. In this model, the ferromagnetic transition-metal lattice, which orders at T-c, creates a strong local field at the erbium site, polarizing the Er moments in a direction opposite to the applied field. At low temperatures, when the contribution of the paramagnetic erbium sublattice, which varies as T-1, gets larger than the ferromagnetic contribution, the total magnetic moment changes its sign, leading to an overall ferrimagnetic state. The half-substituted compound ErCo0.50Mn0.50O3 was studied in detail, since the magnetization loops present two well-identified anomalies: an intersection of the magnetization branches at low fields, and magnetization jumps at high fields. The influence of the oxidizing conditions was studied in other compositions close to the 50/50 = Mn/Co substitution rate. These anomalies are clearly connected to the spin inversion phenomena and to the simultaneous presence of Co2+ and Co3+ magnetic moments. Dynamical aspects should be considered to well identify the high-field anomaly, since it depends on the magnetic field sweep rate. (C) 2006 Elsevier B.V. All rights reserved.

Utilização do coque verde de petróleo da Petrobras em substituição parcial ao carvão mineral das misturas injetadas em altos fornos

The Petrobras produced green petroleum coke (GPC) is a carbon rich fuel, virtually ash-free, with low sulfur content, and is a fuel suitable to replace metallurgical coke in blast furnaces. The GPC was tested in a pulverized coal injection simulator built in the Volta Redonda research center. It presented a low burning efficiency due to low volatile material content and high substitution rate by the carbon content. The tests were carried out in blast furnaces with ≤ 50% Petrobras GPC in the coal blends. The injected coal/CVP mixtures produced no negative side effects in the blast furnace grinding systems, pneumatic conveying, or operating process. The mixture burning process inside the blast furnace showed a decrease in fuel consumption, with a significant reduction in metallurgical coke consumption. The industrial-scale tests of the GPC mixtures did not reach the 70% maximum for lack of the GPC feedstock, which is necessary to continue with standard coal mixtures.

Molecular phylogenetic inferences on the position of the Mollusca within the Lophotrochozoa

Die phylogenetische Position der Mollusken innerhalb der Trochozoa sowie die interne Evolution der Klassen der Mollusca sind weitgehend unbekannt und wurden in meiner Arbeit anhand molekularer Merkmale untersucht. Phylogenomische Analysen zeigten in der Vergangenheit eine gute Auflösung für ursprüngliche Speziationsereignisse. Daher wurden hier drei neue EST Datensätze generiert: für Sipunculus nudus (Sipuncula), Barentsia elongata (Kamptozoa) und Lepidochitona cinerea, (Polyplacophora, Mollusca). Zusätzlich wurden gezielt Gene verschiedener Mollusken mittels RT-PCR amplifiziert. rnSowohl Kamptozoen als auch Sipunculiden wurden aufgrund morphologischer Kriterien bisher als mögliche Schwestergruppe der Mollusken gehandelt, aber die hier erzielten Ergebnisse zur Evolution der Hämerythrine, Gen-Anordnungen der mitochondrialen Genome und phylogenetische Analysen der ribosomalen und der mitochondriellen Proteine stützen diese Hypothese nicht. Die Position der Kamptozoa erwies sich hier generell als unbeständig; phylogenomische Analysen deuten eine Nähe zu den Bryozoen an, aber diese Position wird stark durch die Auswahl der Taxa beeinflusst. Dagegen weisen meine Analysen klar auf eine nähere Beziehung zwischen Annelida und Sipuncula hin. Die ribosomalen Proteine zeigen Sipuncula (und Echiura) sogar als Subtaxa der Anneliden. Wie den Mollusken fehlt den Sipunculiden jegliche Segmentierung und meine Ergebnisse legen hier die Möglichkeit des Verlusts dieses Merkmals innerhalb der Anneliden bei den Sipunculiden nahe. Innerhalb der Mollusken wurden die Solenogastren bereits als Schwestergruppe aller rezenten Mollusken vorgeschlagen. Im Rahmen meiner Arbeit wurden von drei verschiedenen Solenogastren-Arten die ersten zuverlässigen 18S rRNA-Sequenzen ermittelt, und es zeigte sich, dass alle bisher veröffentlichten 18S-Sequenzen dieser Molluskenklasse höchst unvollständig oder fehlerhaft sind. rnRibosomale Proteine sind gute phylogenetische Marker und hier wurden die Auswahl und Anzahl dieser Gene für phylogenetische Analysen optimiert. Über Sonden-basierte Detektion wurde eine sampling-Strategie getestet, die im Vergleich mit standard-phylogenomischen Ansätzen zukünftige molekulare Stammbaumrekonstruktionen mit größerem Taxonsampling ermöglicht.rn

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

Molecular evolution of primate immunodeficiency viruses and hepatitis delta virus

Primate immunodeficiency viruses, or lentiviruses (HIV-1, HIV-2, and SIV), and hepatitis delta virus (HDV) are RNA viruses characterized by rapid evolution. Infection by primate immunodeficiency viruses usually results in the development of acquired immunodeficiency syndrome (AIDS) in humans and AIDS-like illnesses in Asian macaques. Similarly, hepatitis delta virus infection causes hepatitis and liver cancer in humans. These viruses are heterogeneous within an infected patient and among individuals. Substitution rates in the virus genomes are high and vary in different lineages and among sites. Methods of phylogenetic analysis were applied to study the evolution of primate lentiviruses and the hepatitis delta virus. The following results have been obtained: (1) The substitution rate varies among sites of primate lentivirus genes according to the two parameter gamma distribution, with the shape parameter $\alpha$ being close to 1. (2) Primate immunodeficiency viruses fall into species-specific lineages. Therefore, viral transmissions across primate species are not as frequent as suggested by previous authors. (3) Primate lentiviruses have acquired or lost their pathogenicity several times in the course of evolution. (4) Evidence was provided for multiple infections of a North American patient by distinct HIV-1 strains of the B subtype. (5) Computer simulations indicate that the probability of committing an error in testing HIV transmission depends on the number of virus sequences and their length, the divergence times among sequences, and the model of nucleotide substitution. (6) For future investigations of HIV-1 transmissions, using longer virus sequences and avoiding the use of distant outgroups is recommended. (7) Hepatitis delta virus strains are usually related according to the geographic region of isolation. (8) Evolution of HDV is characterized by the rate of synonymous substitution being lower than the nonsynonymous substitution rate and the rate of evolution of the noncoding region. (9) There is a strong preference for G and C nucleotides at the third codon positions of the HDV coding region. ^

Long-term stability of contour augmentation in the esthetic zone: histologic and histomorphometric evaluation of 12 human biopsies 14 to 80 months after augmentation.

BACKGROUND Contour augmentation around early-placed implants (Type 2 placement) using autogenous bone chips combined with deproteinized bovine bone mineral (DBBM) and a collagen barrier membrane has been documented to predictably provide esthetically satisfactory clinical outcomes. In addition, recent data from cone beam computed tomography studies have shown the augmented volume to be stable long-term. However, no human histologic data are available to document the tissue reactions to this bone augmentation procedure. METHODS Over an 8-year period, 12 biopsies were harvested 14 to 80 months after implant placement with simultaneous contour augmentation in 10 patients. The biopsies were subjected to histologic and histomorphometric analysis. RESULTS The biopsies consisted of 32.0% ± 9.6% DBBM particles and 40.6% ± 14.6% mature bone. 70.3% ± 14.5% of the DBBM particle surfaces were covered with bone. On the remaining surface, multinucleated giant cells with varying intensity of tartrate-resistant acid phosphatase staining were regularly present. No signs of inflammation were visible, and no tendency toward a decreasing volume fraction of DBBM over time was observed. CONCLUSIONS The present study confirms previous findings that osseointegrated DBBM particles do not tend to undergo substitution over time. This low substitution rate may be the reason behind the clinically and radiographically documented long-term stability of contour augmentation using a combination of autogenous bone chips, DBBM particles, and a collagen membrane.