Contrast-FEL—A Test for Differences in Selective Pressures at Individual Sites among Clades and Sets of Branches

A number of evolutionary hypotheses can be tested by comparing selective pressures among sets of branches in a phylogenetic tree. When the question of interest is to identify specific sites within genes that may be evolving differently, a common approach is to perform separate analyses on subsets of sequences and compare parameter estimates in a post hoc fashion. This approach is statistically suboptimal and not always applicable. Here, we develop a simple extension of a popular fixed effects likelihood method in the context of codon-based evolutionary phylogenetic maximum likelihood testing, Contrast-FEL. It is suitable for identifying individual alignment sites where any among the K≥2 sets of branches in a phylogenetic tree have detectably different ω ratios, indicative of different selective regimes. Using extensive simulations, we show that Contrast-FEL delivers good power, exceeding 90% for sufficiently large differences, while maintaining tight control over false positive rates, when the model is correctly specified. We conclude by applying Contrast-FEL to data from five previously published studies spanning a diverse range of organisms and focusing on different evolutionary questions.     

Metabolism of 20-hydroxyeicosatetraenoic acid by cyclooxygenase. Formation and identification of novel endothelium-dependent vasoconstrictor metabolites

We recently demonstrated that 20-hydroxyeicosatetraenoic acid (20-HETE) constricts rat aortic rings. The contractile response was partially dependent on the presence of endothelium and was abolished by pretreatment of the rings with either indomethacin or the endoperoxide/thromboxane receptor antagonist, SQ29548. Addition of GSH or SnCl2 to the organ bath diminished the contractile response of 20-HETE, whereas preincubation of the rings with a thromboxane synthase inhibitor did not affect the 20-HETE induced contractions. Short time incubation (2 min) of 20-HETE with ram seminal vesicle microsomes in the presence of p-hydroxymercurybenzoate yielded metabolites which migrated similarly on thin layer chromatography to the known arachidonate endoperoxides prostaglandin (PG) G2 and PGH2 and possess vasoconstrictory properties. The vasoconstriction was dose-dependent with a half-life of approximately 6.3 +/- 0.6 min. Addition of SQ29548 to the aortic ring bath 1 min after metabolite elicited vasoconstriction produced immediate relaxation. Furthermore, pretreatment of the rings with SQ29548 totally abolished the contraction. SnCl2 reduction of the metabolites produced in incubation of rat seminal vesicles with 20-HETE and p-hydroxymercurybenzoate resulted in a single radioactive peak which was further identified by gas chromatography/mass spectrometry as 20-hydroxy-PGF2 alpha. The inhibitory effect of SQ29548, the appearance of labile metabolites with a half-life of approximately 6 min and the production of 20-hydroxy-PGF2 alpha by SnCl2 reduction clearly indicate that the vasoconstrictor metabolites of 20-HETE are the labile endoperoxides of 20-HETE, 20-hydroxy-PGG2, and 20-hydroxy-PGH2.     

MORPHOLOGICAL EVOLUTION IN MUROID RODENTS I. CONSERVATIVE PATTERNS OF CRANIOMETRIC COVARIANCE AND THEIR ONTOGENETIC BASIS IN THE NEOTROPICAL GENUS ZYGODONTOMYS

Analyses of craniodental measurement data from 15 wild-collected population samples of the Neotropical muroid rodent genus Zygodontomys reveal consistent patterns of relative variability and correlation that suggest a common latent structure. Eigenanalysis of each sample covariance matrix of logarithms yields a first principal component that accounts for a large fraction of the total variance. Variances of subsequent sample principal components are much smaller, and the results of bootstrap resampling together with asymptotic statistics suggest that characteristic roots of the covariance matrix after the first are seldom distinct. The coefficients of normalized first principal components are strikingly similar from sample to sample: inner products of these vectors reveal an average between-sample correlation of 0.989, and the mean angle of divergence is only about eight degrees. Since first principal component coefficients identify the same contrasts among variables as comparisons of relative variability and correlation, we conclude that a single factor accounts for most of the common latent determination of these sample dispersions. Analyses of variance based on toothwear (a coarse index of age) and sex in the wild-collected samples, and on known age and sex in a captive-bred population, reveal that specimen scores on sample first principal components are age- and sex-dependent; residual sample dispersion, however, is essentially unaffected by age, sex, or age × sex interaction. The sample first principal component therefore reflects the covariance among measured dimensions induced by general growth, and its coefficients are interpretable as exponents of postnatal growth allometry. Path-analytic models that incorporate prior knowledge of the equivalent allometric effects of general growth within these samples can be used to decompose the between-sample variance by factors corresponding to other ontogenetic mechanisms of form change. The genetic or environmental determinants of differences in sample mean phenotypes induced by such mechanisms, however, can be demonstrated only by experiment.     

Evolutionary origin of human and primate malarias: evidence from the circumsporozoite protein gene

We have analyzed the conserved regions of the gene coding for the circumsporozoite protein (CSP) in 12 species of Plasmodium, the malaria parasite. The closest evolutionary relative of P. falciparum, the agent of malignant human malaria, is P. reichenowi, a chimpanzee parasite. This is consistent with the hypothesis that P. falciparum is an ancient human parasite, associated with humans since the divergence of the hominids from their closest hominoid relatives. Three other human Plasmodium species are each genetically indistinguishable from species parasitic to nonhuman primates; that is, for the DNA sequences included in our analysis, the differences between species are not greater than the differences between strains of the human species. The human P. malariae is indistinguishable from P. brasilianum, and P. vivax is indistinguishable from P. simium; P. brasilianum and P. simium are parasitic to New World monkeys. The human P. vivax-like is indistinguishable from P. simiovale, a parasite of Old World macaques. We conjecture that P. malariae, P. vivax, and P. vivax-like are evolutionarily recent human parasites, the first two at least acquired only within the last several thousand years, and perhaps within the last few hundred years, after the expansion of human populations in South America following the European colonizations. We estimate the rate of evolution of the conserved regions of the CSP gene as 2.46 x 10(-9) per site per year. The divergence between the P. falciparum and P. reichenowi lineages is accordingly dated 8.9 Myr ago. The divergence between the three lineages leading to the human parasites is very ancient, about 100 Myr old between P. malariae and P. vivax (and P. vivax-like) and about 165 Myr old between P. falciparum and the other two.      

The nitrogen cycle in a ‘Terra Firme’ rainforest on oxisol in the Amazon territory of Venezuela

Standing stocks and fluxes of nitrogen, including nitrogen fixation and denitrification, were measured in a tropical rainforest on Oxisol in the Amazon Territory of Venezuela. The standing stock of nitrogen was comparable to that of temperate forests, but was higher than that in an adjacent forest on Spodosol. Fluxes were higher than in forests in the temperate zone, but lower than in another tropical forest on more fertile soil. Even though nitrogen was abundant, this does not mean that nitrogen could not be limiting to agriculture if the forest is cleared and the land cultivated. The nitrogen fixing and nitrogen conserving mechanisms are dependent upon the structure of the undisturbed forest, and destruction of the forest would eventually decrease the input of nitrogen to the soil.     

Distribution patterns of Dikarya in arid and semiarid soils of Baja California,Mexico

Approximately four-fifths of the land area of Baja California (BC) in Mexico are occupied by arid and semiarid soils, the mycobiota of which is virtually uncharacterized. In the first culture-independent study of the mycobiota of BC, we collected soil from five different locations in the region and constructed a Dikarya-specific gene library for the ITS region of nuclear ribosomal DNA. Clones were analyzed by RFLP, were sequenced for phylogenetic analyses, and diversity and similarity indices were calculated. The ascomycete Penicillium dipodomyicola was the most frequent fungus found in soil at the most arid location studied, and the basidiomycete Coprinellus radians was the most frequent at the location receiving the highest rainfall. Other frequent members of the soil mycobiota were identified as Alternaria spp., Ceratobasidium sp., Coniozyma leucospermi, Nematoctonus robustus, Penicillium griseofulvum, Tulostoma kotlabae and uncultured members of the Dikarya. Several sequences were identified as those of uncultured fungi, one of which was previously reported from other hot deserts. Arid soils and the transitional zones between arid and semiarid soils had the most similar fungal diversity, with the former soils having a community from which basidiomycetes were absent, and the soil receiving the highest precipitation having a community dominated by basidiomycetes.     

Historical Shifts in Brazilian P. falciparum Population Structure and Drug Resistance Alleles

Previous work suggests that Brazilian Plasmodium falciparum has limited genetic diversity and a history of bottlenecks, multiple reintroductions due to human migration, and clonal expansions. We hypothesized that Brazilian P. falciparum would exhibit clonal structure. We examined isolates collected across two decades from Amapá, Rondônia, and Pará state (n = 190). By examining more microsatellites markers on more chromosomes than previous studies, we hoped to define the extent of low diversity, linkage disequilibrium, bottlenecks, population structure, and parasite migration within Brazil. We used retrospective genotyping of samples from the 1980s and 1990s to explore the population genetics of SP resistant dhfr and dhps alleles. We tested an existing hypothesis that the triple mutant dhfr mutations 50R/51I/108N and 51I/108N/164L developed in southern Amazon from a single origin of common or similar parasites. We found that Brazilian P. falciparum had limited genetic diversity and isolation by distance was rejected, which suggests it underwent bottlenecks followed by migration between sites. Unlike Peru, there appeared to be gene flow across the Brazilian Amazon basin. We were unable to divide parasite populations by clonal lineages and pairwise FST were common. Most parasite diversity was found within sites in the Brazilian Amazon, according to AMOVA. Our results challenge the hypothesis that triple mutant alleles arose from a single lineage in the Southern Amazon. SP resistance, at both the double and triple mutant stages, developed twice and potentially in different regions of the Brazilian Amazon. We would have required samples from before the 1980s to describe how SP resistance spread across the basin or describe the complex internal migration of Brazilian parasites after the colonization efforts of past decades. The Brazilian Amazon basin may have sufficient internal migration for drug resistance reported in any particular region to rapidly spread to other parts of basin under similar drug pressure.     

Real-Time Loop-Mediated Isothermal Amplification (RealAmp) for the Species-Specific Identification of Plasmodium vivax

Plasmodium vivax infections remain a major source of malaria-related morbidity and mortality. Early and accurate diagnosis is an integral component of effective malaria control programs. Conventional molecular diagnostic methods provide accurate results but are often resource-intensive, expensive, have a long turnaround time and are beyond the capacity of most malaria-endemic countries. Our laboratory has recently developed a new platform called RealAmp, which combines loop-mediated isothermal amplification (LAMP) with a portable tube scanner real-time isothermal instrument for the rapid detection of malaria parasites. Here we describe new primers for the detection of P. vivax using the RealAmp method. Three pairs of amplification primers required for this method were derived from a conserved DNA sequence unique to the P. vivax genome. The amplification was carried out at 64°C using SYBR Green or SYTO-9 intercalating dyes for 90 minutes with the tube scanner set to collect fluorescence signals at 1-minute intervals. Clinical samples of P. vivax and other human-infecting malaria parasite species were used to determine the sensitivity and specificity of the primers by comparing with an 18S ribosomal RNA-based nested PCR as the gold standard. The new set of primers consistently detected laboratory-maintained isolates of P. vivax from different parts of the world. The primers detected P. vivax in the clinical samples with 94.59% sensitivity (95% CI: 87.48–98.26%) and 100% specificity (95% CI: 90.40–100%) compared to the gold standard nested-PCR method. The new primers also proved to be more sensitive than the published species-specific primers specifically developed for the LAMP method in detecting P. vivax.