Summary: The program MODELTEST uses log likelihood scores to establish the model of DNA evolution that best fits the data. Availability: The MODELTEST package, including the source code and some documentation is available at http://bioag.byu.edu/zoology/crandall―lab/modeltest.html. Contact: dp47@email.byu.edu.

/pdf/modeltest-testing-the-model-of-dna-substitution-2dmippvvyz.pdf

MODELTEST: testing the model of DNA substitution.

PhyML is a phylogeny software based on the maximum-likelihood principle. Early PhyML versions used a fast algorithm performing nearest neighbor interchanges to improve a reasonable starting tree topology. Since the original publication (Guindon S., Gascuel O. 2003. A simple, fast and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst. Biol. 52:696-704), PhyML has been widely used (>2500 citations in ISI Web of Science) because of its simplicity and a fair compromise between accuracy and speed. In the meantime, research around PhyML has continued, and this article describes the new algorithms and methods implemented in the program. First, we introduce a new algorithm to search the tree space with user-defined intensity using subtree pruning and regrafting topological moves. The parsimony criterion is used here to filter out the least promising topology modifications with respect to the likelihood function. The analysis of a large collection of real nucleotide and amino acid data sets of various sizes demonstrates the good performance of this method. Second, we describe a new test to assess the support of the data for internal branches of a phylogeny. This approach extends the recently proposed approximate likelihood-ratio test and relies on a nonparametric, Shimodaira-Hasegawa-like procedure. A detailed analysis of real alignments sheds light on the links between this new approach and the more classical nonparametric bootstrap method. Overall, our tests show that the last version (3.0) of PhyML is fast, accurate, stable, and ready to use. A Web server and binary files are available from http://www.atgc-montpellier.fr/phyml/.

https://www.zora.uzh.ch/id/eprint/155669/1/ZORA_NL_155669.pdf

New Algorithms and Methods to Estimate Maximum-Likelihood Phylogenies: Assessing the Performance of PhyML 3.0

A set of oligonucleotide primers capable of initiating enzymatic amplification (polymerase chain reaction) on a phylogenetically and taxonomically wide range of bacteria is described along with methods for their use and examples. One pair of primers is capable of amplifying nearly full-length 16S ribosomal DNA (rDNA) from many bacterial genera; the additional primers are useful for various exceptional sequences. Methods for purification of amplified material, direct sequencing, cloning, sequencing, and transcription are outlined. An obligate intracellular parasite of bovine erythrocytes, Anaplasma marginale, is used as an example; its 16S rDNA was amplified, cloned, sequenced, and phylogenetically placed. Anaplasmas are related to the genera Rickettsia and Ehrlichia. In addition, 16S rDNAs from several species were readily amplified from material found in lyophilized ampoules from the American Type Culture Collection. By use of this method, the phylogenetic study of extremely fastidious or highly pathogenic bacterial species can be carried out without the need to culture them. In theory, any gene segment for which polymerase chain reaction primer design is possible can be derived from a readily obtainable lyophilized bacterial culture.

16S ribosomal DNA amplification for phylogenetic study.

jModelTest is a new program for the statistical selection of models of nucleotide substitution based on "Phyml" (Guindon and Gascuel 2003. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 52:696-704.). It implements 5 different selection strategies, including "hierarchical and dynamical likelihood ratio tests," the "Akaike information criterion," the "Bayesian information criterion," and a "decision-theoretic performance-based" approach. This program also calculates the relative importance and model-averaged estimates of substitution parameters, including a model-averaged estimate of the phylogeny. jModelTest is written in Java and runs under Mac OSX, Windows, and Unix systems with a Java Runtime Environment installed. The program, including documentation, can be freely downloaded from the software section at http://darwin.uvigo.es.

/pdf/jmodeltest-phylogenetic-model-averaging-4hpp3p5d0f.pdf

jModelTest: Phylogenetic Model Averaging

Motivation: In 2001 and 2002, we published two papers (Bioinformatics, 17, 282--283, Bioinformatics, 18, 77--82) describing an ultrafast protein sequence clustering program called cd-hit. This program can efficiently cluster a huge protein database with millions of sequences. However, the applications of the underlying algorithm are not limited to only protein sequences clustering, here we present several new programs using the same algorithm including cd-hit-2d, cd-hit-est and cd-hit-est-2d. Cd-hit-2d compares two protein datasets and reports similar matches between them; cd-hit-est clusters a DNA/RNA sequence database and cd-hit-est-2d compares two nucleotide datasets. All these programs can handle huge datasets with millions of sequences and can be hundreds of times faster than methods based on the popular sequence comparison and database search tools, such as BLAST.

Availability: http://cd-hit.org

Contact: [email protected]

Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences

Mammalian protein metabolism

The precision of ultraviolet absorption measurements in the Schmidt-Thannhauser procedure for nucleic acid estimation

The mRNAs for the heavy and light subunits of the iron-storage protein ferritin occur in cells largely as inactive ribonucleoprotein particles, which are recruited for translation when iron enters the cell. Cytoplasmic extracts from rat tissues and hepatoma cells were shown by an electrophoretic separation procedure to form RNA-protein complexes involving a highly conserved sequence in the 5' untranslated region of both ferritin heavy- and light-subunit mRNAs. The pattern of complex formation was affected by pretreatment of rats or cells with iron. Crosslinking by UV irradiation showed that the complexes contained an 87-kDa protein interacting with the conserved sequence of the ferritin mRNA. We propose that intracellular iron levels regulate ferritin synthesis by causing changes in specific protein binding to the conserved sequence in the ferritin heavy- and light-subunit mRNAs.

Cytoplasmic protein binds in vitro to a highly conserved sequence in the 5' untranslated region of ferritin heavy- and light-subunit mRNAs

The mRNAs for the heavy and light subunits of the iron-storage protein ferritin occur in cells largely as inactive ribonucleoprotein particles, which are recruited for translation when iron enters the cell. Cytoplasmic extracts from rat tissues and hepatoma cells were shown by an electro- phoretic separation procedure to form RNA-protein com- plexes involving a highly conserved sequence in the 5' untrans- lated region of both ferritin heavy- and light-subunit mRNAs. The pattern of complex formation was affected by pretreat- ment of rats or cells with iron. Crosslinking by UV irradiation showed that the complexes contained an 87-kDa protein inter- acting with the conserved sequence of the ferritin mRNA. We propose that intracellular iron levels regulate ferritin synthesis by causing changes in specific protein binding to the conserved sequence in the ferritin heavy- and light-subunit mRNAs. Femtin, an iron-storage protein found in animal, plant, fungal, and bacterial cells (1), has a protein shell (Mr = 500,000) consisting in vertebrates of 24 subunits of two kinds, heavy (H; Mr 20,000) and light (L; Mr 19,000).

Cytoplasmic protein binds in vitro to a highly conserved sequence in the 5' untranslated region of ferritin heavy- and light-subunit mRNAs (RNA-protein complexes/iron/translational control)

In previous studies, we showed that acute administration of iron to intact rats or to rat hepatoma cells in culture induces synthesis of the iron-storage protein ferritin by activating translation of inactive cytoplasmic ferritin mRNAs for both the heavy (H) and the light (L) subunits. In the course of activation, these ferritin mRNAs are recruited onto polysomes. To elucidate the structural features of these mRNAs involved in the translational response to iron, a chimera was constructed from the 5' and 3' untranslated regions (UTRs) of ferritin L subunit mRNA fused to the reading frame of the mRNA of bacterial chloramphenicol acetyltransferase (CAT). This chimera and deletion constructs derived from it were introduced into a rat hepatoma cell line by retrovirus-mediated gene transfer. The complete chimera showed increased CAT activity in response to iron enrichment of the medium, whereas deletion of the first 67 nucleotides of the 5' UTR, which contain a highly conserved sequence, caused loss of regulation by iron. Whereas cis-acting sequences located in the 5' flanking regions of many genes have been repeatedly implicated in modulating their transcriptional expression, we report here a specific regulatory translational sequence found within the 5' UTR of a eukaryotic mRNA.

https://www.pnas.org/content/pnas/84/23/8478.full.pdf

Hamish N. Munro

Papers

Mammalian protein metabolism

The precision of ultraviolet absorption measurements in the Schmidt-Thannhauser procedure for nucleic acid estimation

Cytoplasmic protein binds in vitro to a highly conserved sequence in the 5' untranslated region of ferritin heavy- and light-subunit mRNAs

Cytoplasmic protein binds in vitro to a highly conserved sequence in the 5' untranslated region of ferritin heavy- and light-subunit mRNAs (RNA-protein complexes/iron/translational control)

Iron regulates ferritin mRNA translation through a segment of its 5' untranslated region.