Predominant contribution of cis‐regulatory divergence in the evolution of mouse alternative splicing

Divergence of alternative splicing represents one of the major driving forces to shape phenotypic diversity during evolution. However, the extent to which these divergences could be explained by the evolving cis‐regulatory versus trans‐acting factors remains unresolved. To globally investigate the relative contributions of the two factors for the first time in mammals, we measured splicing difference between C57BL/6J and SPRET/EiJ mouse strains and allele‐specific splicing pattern in their F1 hybrid. Out of 11,818 alternative splicing events expressed in the cultured fibroblast cells, we identified 796 with significant difference between the parental strains. After integrating allele‐specific data from F1 hybrid, we demonstrated that these events could be predominately attributed to cis‐regulatory variants, including those residing at and beyond canonical splicing sites. Contrary to previous observations in Drosophila, such predominant contribution was consistently observed across different types of alternative splicing. Further analysis of liver tissues from the same mouse strains and reanalysis of published datasets on other strains showed similar trends, implying in general the predominant contribution of cis‐regulatory changes in the evolution of mouse alternative splicing.     

Rapid evolution of prostatic protein PSP94 suggested by sequence divergence between rhesus monkey and human cDNAs

Comparison of the protein coding region of mRNA for the prostatic secretory protein PSP94 in human (hPSP94) with that in rhesus monkey (rmPSP94) indicates that, for the most part, its sequence has evolved with few constraints and at a relatively fast rate. Interestingly, half of the 22 residue differences between the two species involve charge changes, reflected by the acidic pI (5.4) of hPSP94 and the basic pI (10.6) of rmPSP94. However, the 10 cysteines and 5 of the 6 prolines of PSP94 were unaffected, suggesting that the three-dimensional conformations of the human and the monkey proteins may be similar. Rapid evolution of this gene might explain the apparent absence in nonprimates of homologous sequences detectable by hybridization.     

A novel comparative method for identifying shifts in the rate of character evolution on trees

Evolutionary biologists since Darwin have been fascinated by differences in the rate of trait-evolutionary change across lineages. Despite this continued interest, we still lack methods for identifying shifts in evolutionary rates on the growing tree of life while accommodating uncertainty in the evolutionary process. Here we introduce a Bayesian approach for identifying complex patterns in the evolution of continuous traits. The method (auteur) uses reversible-jump Markov chain Monte Carlo sampling to more fully characterize the complexity of trait evolution, considering models that range in complexity from those with a single global rate to potentially ones in which each branch in the tree has its own independent rate. This newly introduced approach performs well in recovering simulated rate shifts and simulated rates for datasets nearing the size typical for comparative phylogenetic study (i.e., ≥64 tips). Analysis of two large empirical datasets of vertebrate body size reveal overwhelming support for multiple-rate models of evolution, and we observe exceptionally high rates of body-size evolution in a group of emydid turtles relative to their evolutionary background. auteur will facilitate identification of exceptional evolutionary dynamics, essential to the study of both adaptive radiation and stasis.     

Intra- and interspecific variation of the CCR5 gene in higher primates

We have evaluated the molecular evolution of the chemokine receptor CCR5 in primates. The chemokine receptor CCR5 serves as a major co-receptor for human immunodeficiency virus/simian immunodeficiency virus (HIV/SIV) infection. Knowledge of evolution of the CCR5 molecule and selection on the CCR5 gene may shed light on its functional role. The comparison of differences between intraspecific polymorphisms and interspecific fixed substitutions provides useful information regarding modes of selection during the course of evolution. There is marked polymorphism in the CCR5 gene sequence within different primate species, whereas sequence divergence between different species is small. By using contingency tests, we compared synonymous (SS) and nonsynonymous (NS) CCR5 mutations occurring within and between a broad range of primates. Our results demonstrate that CCR5 evolution did not follow expectations of strict neutrality at the level of the whole gene. The proportion of NS to SS at the intraspecific level was significantly higher than that observed at the interspecific level. These results suggest that most CCR5 NS polymorphisms are slightly deleterious. However, at domains more closely correlated with its known biological functions, there was no obvious evidence to support deviation from neutrality.     

Time dependency of molecular rate estimates and systematic overestimation of recent divergence times

Studies of molecular evolutionary rates have yielded a wide range of rate estimates for various genes and taxa. Recent studies based on population-level and pedigree data have produced remarkably high estimates of mutation rate, which strongly contrast with substitution rates inferred in phylogenetic (species-level) studies. Using Bayesian analysis with a relaxed-clock model, we estimated rates for three groups of mitochondrial data: avian protein-coding genes, primate protein-coding genes, and primate d-loop sequences. In all three cases, we found a measurable transition between the high, short-term (< 1-2 Myr) mutation rate and the low, long-term substitution rate. The relationship between the age of the calibration and the rate of change can be described by a vertically translated exponential decay curve, which may be used for correcting molecular date estimates. The phylogenetic substitution rates in mitochondria are approximately 0.5% per million years for avian protein-coding sequences and 1.5% per million years for primate protein-coding and d-loop sequences. Further analyses showed that purifying selection offers the most convincing explanation for the observed relationship between the estimated rate and the depth of the calibration. We rule out the possibility that it is a spurious result arising from sequence errors, and find it unlikely that the apparent decline in rates over time is caused by mutational saturation. Using a rate curve estimated from the d-loop data, several dates for last common ancestors were calculated: modern humans and Neandertals (354 ka; 222-705 ka), Neandertals (108 ka; 70-156 ka), and modern humans (76 ka; 47-110 ka). If the rate curve for a particular taxonomic group can be accurately estimated, it can be a useful tool for correcting divergence date estimates by taking the rate decay into account. Our results show that it is invalid to extrapolate molecular rates of change across different evolutionary timescales, which has important consequences for studies of populations, domestication, conservation genetics, and human evolution.     

Cloning of two splice variants of the rice PTS1 receptor, OsPex5pL and OsPex5pS, and their functional characterization using pex5-deficient yeast and Arabidopsis

Using the rice PEX14 cDNA as a bait in a yeast two-hybrid assay, two splice variants of the type I peroxisomal targeting signal (PTS1) receptor, OsPex5pL and OsPex5pS, were cloned from a pathogen-treated rice leaf cDNA library. The proteins were produced from a single gene by alternative splicing, which generated a full-length variant, OsPEX5L, and a variant that lacked exon 7, OsPEX5S. OsPex5pL contained 11 copies of the pentapeptide motif WXXXF/Y in its N-terminus, and seven tetratricopeptide repeats in its C-terminus. Expression of OsPEX5L and OsPEX5S predominantly occurred in leaf tissues, and was induced by various stresses, such as exposure to the pathogen Magnaporthe grisea, and treatment with fungal elicitor, methyl viologen, NaCl or hydrogen peroxide. The Arabidopsis T-DNA insertional pex5 mutant, Atpex5, which does not germinate in the absence of sucrose and was resistant to indole-3-butyric acid (IBA), was perfectly rescued by over-expression of OsPex5pL, but not by OsPex5pS. Using transient expression of OsPex5pL and OsPex5pS in the Atpex5 mutant, we show that OsPex5pL translocates both PTS1- and PTS2-containing proteins into the peroxisome by interacting with OsPex7p, whereas OsPex5pS is involved only in PTS1-dependent import in Arabidopsis.     

Male-biased mutation, sex linkage, and the rate of adaptive evolution

An interaction between sex-linked inheritance and sex-biased mutation rates may affect the rate of adaptive evolution. Males have much higher mutation rates than females in several vertebrate and plant taxa. When evolutionary rates are limited by the supply of favorable new mutations, then genes will evolve faster when located on sex chromosomes that spend more time in males. For mutations with additive effects, Y-linked genes evolve fastest, followed by Z-linked genes, autosomal genes, X-linked genes, and finally W-linked and cytoplasmic genes. This ordering can change when mutations show dominance. The predicted differences in substitution rates may be detectable at the molecular level. Male-biased mutation could cause adaptive changes to accumulate more readily on certain kinds of chromosomes and favor animals with Z-W sex determination to have rapidly evolving male sexual displays.     

An intraspecific comparative analysis of character divergence between sympatric species

Although sympatric character divergence between closely related species has been described in a wide variety of taxa, the evolutionary processes responsible for generating these patterns are difficult to identify. One hypothesis that can explain sympatric differences is ecological character displacement: the sympatric origin of morphologically divergent phenotypes in response to selection caused by interspecific competition. Alternatively, populations may adapt to different conditions in allopatry, with sympatric distributions evolving through selective colonization and proliferation of ecologically compatible phenotypes. In this study, I characterize geographic variation within two sibling species of rocky-shore gastropods that have partially overlapping distributions in central California. In sympatry, both Nucella emarginata and N. ostrina show significant differences in shell shape and shell ornamentation that together suggest that where the two species co-exist, divergent phenotypes arose as an evolutionary consequence of competition. To examine the evolutionary origins of divergent characters in sympatry, I used a comparative method based on spatial autocorrelation to remove the portion of the phenotypic variance among populations that is explained by genetic distance (using mitochondrial DNA sequences and allozyme frequency data). Because the remaining portion of the phenotypic variance represents the independent divergence of individual populations, a significant sympatric difference in the corrected dataset provides evidence of true character displacement: significant sympatric character evolution that is independent of population history. After removal of genetic distance effects in Nucella, shell shape differences remain statistically significant in N. emarginata, providing evidence of significant sympatric character divergence. However, for external shell ornamentation in both species and shell shape in N. ostrina, the significance of sympatric differences is lost in the corrected dataset, indicating that colonization events and gene flow have played important roles in the evolutionary history of character divergence in sympatry. Although the absence of a widely dispersing planktonic larva in the life cycle of Nucella will promote local adaptation, the results here indicate that once advantageous traits arise, demographic processes, such as recurrent gene flow between established populations and extinction and recolonization, are important factors contributing to the geographic pattern of sympatric character divergence.     

Divergence population genetic analysis of hybridization between rhesus and cynomolgus macaques

The geographic ranges of rhesus ( Macaca mulatta ) and cynomolgus ( M. fascicularis ) macaques adjoin in Indochina where they appear to hybridize. We used published and newly generated DNA sequences from 19 loci spanning ~20 kb to test whether introgression has occurred between these macaque species. We studied introgression at the level of nuclear DNA and distinguished between incomplete lineage sorting of ancestral polymorphisms or interspecific gene flow. We implemented a divergence population genetics approach by fitting our data to an isolation model implemented in the software IMa. The model that posits no gene flow from the rhesus into the cynomolgus macaque was rejected ( P  = 1.99 × 10⁻⁸). Gene flow in this direction was estimated as 2 Nm ~1.2, while gene flow in the reverse direction was nonsignificantly different from zero ( P  = 0.16). The divergence time between species was estimated as ~1.3 million years. Balancing selection, a special case of incomplete sorting, was taken into consideration, as well as potential crossbreeding in captivity. Parameter estimates varied between analyses of subsets of data, although we still rejected isolation models. Geographic sampling of the data, where samples of cynomolgus macaques derived from Indochina were excluded, revealed a lost signature of gene flow, indicating that interspecific gene flow is restricted to mainland Indochina. Our results, in conjunction with those by others, justify future detailed analyses into the genetics of reproductive barriers and reticulate evolution in these two genome-enabled primates. Future studies of the natural hybridization between rhesus and cynomolgus macaques would expand the repertoire of systems available for speciation studies in primates.     

Rates of genome evolution and branching order from whole genome analysis

Accurate estimation of any phylogeny is important as a framework for evolutionary analysis of form and function at all levels of organization from sequence to whole organism. Using alignments of nonrepetitive components of opossum, human, mouse, rat, and dog genomes we evaluated two alternative tree topologies for eutherian evolution. We show with very high confidence that there is a basal split between rodents (as represented by the mouse and rat) and a branch joining primates (as represented by humans) and carnivores (as represented by dogs), consistent with some but not the most widely accepted mammalian phylogenies. The result was robust to substitution model choice with equivalent inference returned from a spectrum of models ranging from a general time reversible model, a model that treated nucleotides as either purines and pyrimidines, and variants of these that incorporated rate heterogeneity among sites. By determining this particular branching order we are able to show that the rate of molecular evolution is almost identical in rodent and carnivore lineages and that sequences evolve approximately 11%-14% faster in these lineages than in the primate lineage. In addition by applying the chicken as outgroup the analyses suggested that the rate of evolution in all eutherian lineages is approximately 30% slower than in the opossum lineage. This pattern of relative rates is inconsistent with the hypothesis that generation time is an important determinant of substitution rates and, by implication, mutation rates. Possible factors causing rate differences between the lineages include differences in DNA repair and replication enzymology, and shifts in nucleotide pools. Our analysis demonstrates the importance of using multiple sequences from across the genome to estimate phylogeny and relative evolutionary rate in order to reduce the influence of distorting local effects evident even in relatively long sequences.     

The chloroplast genome of Phalaenopsis aphrodite (Orchidaceae): comparative analysis of evolutionary rate with that of grasses and its phylogenetic implications

Whether the Amborella/Amborella-Nymphaeales or the grass lineage diverged first within the angiosperms has recently been debated. Central to this issue has been focused on the artifacts that might result from sampling only grasses within the monocots. We therefore sequenced the entire chloroplast genome (cpDNA) of Phalaenopsis aphrodite, Taiwan moth orchid. The cpDNA is a circular molecule of 148,964 bp with a comparatively short single-copy region (11,543 bp) due to the unusual loss and truncation/scattered deletion of certain ndh subunits. An open reading frame, orf91, located in the complementary strand of the rrn23 was reported for the first time. A comparison of nucleotide substitutions between P. aphrodite and the grasses indicates that only the plastid expression genes have a strong positive correlation between nonsynonymous (Ka) and synonymous (Ks) substitutions per site, providing evidence for a generation time effect, mainly across these genes. Among the intron-containing protein-coding genes of the sampled monocots, the Ks of the genes are significantly correlated to transitional substitutions of their introns. We compiled a concatenated 61 protein-coding gene alignment for the available 20 cpDNAs of vascular plants and analyzed the data set using Bayesian inference, maximum parsimony, and neighbor-joining (NJ) methods. The analyses yielded robust support for the Amborella/Amborella-Nymphaeales-basal hypothesis and for the orchid and grasses together being a monophyletic group nested within the remaining angiosperms. However, the NJ analysis using Ka, the first two codon positions, or amino acid sequences, respectively, supports the monocots-basal hypothesis. We demonstrated that these conflicting angiosperm phylogenies are most probably linked to the transitional sites at all codon positions, especially at the third one where the strong base-composition bias and saturation effect take place.     

Theoretical Calculations on Zeolites: The Aluminium Substitution in Mordenite,Ferrierite and ZSM-5

Abstract

The experimental determination of Al siting in zeolites involves the use of multiple techniques. Great effort has been made on both, experimental and theoretical approaches. The present study presents a novel methodology for calculating Al/Si replacement energies. Simple semiempirical calculations were applied on Modenite, Ferrierite and ZSM-5 zeolites, resulting in good agreement with the experimetal evidences. We have found that the favored Al substitution sites are T3 and T4 in Mordenite, while T2 and T4 are in Ferrierite, and only the T9 site is favored in ZSM-5. The method presented is based on an average of partial Al/Si replacement energies, evaluated for all rings belonging to each T site, rather than in the calculation of a total replacement energy evaluated for only one representative aggregate.     

Cloning and sequence analysis of two monkey metallothionein cDNAs

We have isolated two metallothionein (MT) cDNA clones copied from the RNA of cadmium-resistant monkey kidney cells. The complete DNA sequences of these clones show that they encode two distinct MTs. One clone appears to represent monkey MT-II, as shown by its close homology to the human MT-II sequence, whereas the second may correspond to monkey MT-I or a related variant metallothionein. Conserved sequences were identified in both the 5′ and 3′ untranslated regions of these clones.