Mitochondrial control region sequence variation within the argali wild sheep (Ovis ammon): evolution and conservation relevance

The phylogenetic relationship of several subspecies of Ovis ammon were analyzed by comparing DNA sequences within the entire mitochondrial D-loop region. Five putative subspecies of ammon (dalai-lamae, darwini, hodgsoni, sairensis and adamerzi) were sampled from four provinces in China [Xinjiang, Qinghai, Gansu and Xizang (Tibet)] and two (servertzovi and nigrimontana) from Uzbekistan. The argalis sampled represent most of the currently recognized putative Subspecies of argali. Analysis of mtDNA sequences revealed high variability within ammon (7.7%), ranging from 2.4 to 11.5%. MaxiMUM-Parsimony tree indicated that nigrimontana from Uzbekistan diverged First, followed by severtzovi from Uzbekistan. The dispersal of argalis into China gave rise to three clades, suggesting that the argali originated in Western Asia and then dispersed throughout the central Asian highlands on a southeastward course. Among the Chinese argalis, mtDNA analysis places dalailamae genetically closer to hodgsoni than to darwini. Severtzovi and.. nigrimontana are two separate subspecies and genetically distinct from the Chinese argali.

Reconstruction of late glacial and Holocene climate evolution in southern China from geolipids and pollen in the Dingnan peat sequence

IEECAS SKLLQG

Genomic organization, nucleotide sequence analysis of the core histone genes cluster in Chlamys farreri and molecular evolution assessment of the H2A and H2B

This work represents the nucleotide sequence of the core histone gene cluster from scallop Chlamys farreri. The tandemly repeated unit of 5671 bp containing a copy of the four core histone genes H4, H2B, H2A and H3 was amplified and identified by the techniques of homology cloning and genomic DNA walking. All the histone genes in the cluster had the structures in their 3' flanking region which related to the evolution of histone gene expression patterns throughout the cell cycle, including two different termination signals, the hairpin structure and at least one AATAAA polyadenylation signal. In their 5' region, the transcription initiation sites with a conserved sequence of 5'-PyATTCPu-3' known as the CAP site were present in all genes except to H2B, generally 37-45 bp upstream of the start code. Canonical TATA and CAAT boxes were identified only in certain histone genes. In the case of the promoters of H2B and H2A genes, there was a 5'-GATCC-3' element, which had been found to be essential to start transcription at the appropriate site. After this element, in the promoter of H2B, there was another sequence, 5'-GGATCGAAACGTTC-3', which was similar to the consensus sequence of 5'-GGAATAAACGTATTC-3' corresponding to the H2B-specific promoter element. The presence of enhancer sequences (5'-TGATATATG-3') was identified from the H4 and H3 genes, matching perfectly with the consensus sequence defined for histone genes. There were several slightly more complex repetitive DNA in the intergene regions. The presence of the series of conserved sequences and reiterated sequences was consistent with the view that mollusc histone gene cluster arose by duplicating of an ancestral precursor histone gene, the birth-and-death evolution model with strong purifying selection enabled the histone cluster less variation and more conserved function. Meanwhile, the H2A and the H2B were demonstrated to be potential good marks for phylogenetic analysis. All the results will be contributed to the characterization of repeating histone gene families in molluscs.

Analysis of Sequence Periodicity in E. coli Proteins: Empirical Investigation of the ?Duplication and Divergence? Theory of Protein Evolution

Gatherer, D., and McEwan, N.R. (2003). Analysis of sequence periodicity in E. coli proteins: empirical investigation of the 'duplication and divergence' theory of protein evolution. Journal of Molecular Evolution 57, 149-158. RAE2008

Sequence divergence of measles virus haemagglutinin during natural evolution and adaptation to cell culture

Phylogenetic analysis of the sequence of the H gene of 75 measles virus (MV) strains (32 published and 43 new sequences) was carried out. The lineage groups described from comparison of the nucleotide sequences encoding the C-terminal regions of the N protein of MV were the same as those derived from the H gene sequences in almost all cases. The databases document a number of distinct genotype switches that have occurred in Madrid (Spain). Well-documented is the complete replacement of lineage group C2, the common European genotype at that time, with that of group D3 around the autumn of 1993. No further isolations of group C2 took place in Madrid after this time. The rate of mutation of the H gene sequences of MV genotype D3 circulating in Madrid from 1993 to 1996 was very low (5 x 10(-4) per annum for a given nucleotide position). This is an order of magnitude lower than the rates of mutation observed in the HN genes of human influenza A viruses. The ratio of expressed over silent mutations indicated that the divergence was not driven by immune selection in this gene. Variations in amino acid 117 of the H protein (F or L) may be related to the ability of some strains to haemagglutinate only in the presence of salt. Adaptation of MV to different primate cell types was associated with very small numbers of mutations in the H gene. The changes could not be predicted when virus previously grown in human B cell lines was adapted to monkey Vero cells. In contrast, rodent brain-adapted viruses displayed a lot of amino acid sequence variation from normal MV strains. There was no convincing evidence for recombination between MV genotypes.

A stacked Late Quaternary fluvio-periglacial sequence from the Axe valley, southern England with implications for landscape evolution and Palaeolithic archaeology

The current model of mid-latitude late Quaternary terrace sequences, is that they are uplift-driven but climatically controlled terrace staircases, relating to both regional-scale crustal and tectonic factors, and palaeohydrological variations forced by quasi-cyclic climatic conditions in the 100 K world (post Mid Pleistocene Transition). This model appears to hold for the majority of the river valleys draining into the English Channel which exhibit 8–15 terrace levels over approximately 60–100 m of altitudinal elevation. However, one valley, the Axe, has only one major morphological terrace and has long-been regarded as anomalous. This paper uses both conventional and novel stratigraphical methods (digital granulometry and terrestrial laser scanning) to show that this terrace is a stacked sedimentary sequence of 20–30 m thickness with a quasi-continuous (i.e. with hiatuses) pulsed, record of fluvial and periglacial sedimentation over at least the last 300–400 K yrs as determined principally by OSL dating of the upper two thirds of the sequence. Since uplift has been regional, there is no evidence of anomalous neotectonics, and climatic history must be comparable to the adjacent catchments (both of which have staircase sequences) a catchment-specific mechanism is required. The Axe is the only valley in North West Europe incised entirely into the near-horizontally bedded chert (crypto-crystalline quartz) and sand-rich Lower Cretaceous rocks creating a buried valley. Mapping of the valley slopes has identified many large landslide scars associated with past and present springs. It is proposed that these are thaw-slump scars and represent large hill-slope failures caused by Vauclausian water pressures and hydraulic fracturing of the chert during rapid permafrost melting. A simple 1D model of this thermokarstic process is used to explore this mechanism, and it is proposed that the resultant anomalously high input of chert and sand into the valley during terminations caused pulsed aggradation until the last termination. It is also proposed that interglacial and interstadial incision may have been prevented by the over-sized and interlocking nature of the sub-angular chert clasts until the Lateglacial when confinement of the river overcame this immobility threshold. One result of this hydrogeologically mediated valley evolution was to provide a sequence of proximal Palaeolithic archaeology over two MIS cycles. This study demonstrates that uplift tectonics and climate alone do not fully determine Quaternary valley evolution and that lithological and hydrogeological conditions are a fundamental cause of variation in terrestrial Quaternary records and landform evolution.

The Genome Sequence of Taurine Cattle: A Window to Ruminant Biology and Evolution

To understand the biology and evolution of ruminants, the cattle genome was sequenced to about sevenfold coverage. The cattle genome contains a minimum of 22,000 genes, with a core set of 14,345 orthologs shared among seven mammalian species of which 1217 are absent or undetected in noneutherian (marsupial or monotreme) genomes. Cattle-specific evolutionary breakpoint regions in chromosomes have a higher density of segmental duplications, enrichment of repetitive elements, and species-specific variations in genes associated with lactation and immune responsiveness. Genes involved in metabolism are generally highly conserved, although five metabolic genes are deleted or extensively diverged from their human orthologs. The cattle genome sequence thus provides a resource for understanding mammalian evolution and accelerating livestock genetic improvement for milk and meat production.

Comparative DNA sequence and methylation analyses of orthologous genes in humans and non-human primates: Genetic and epigenetic footprints of evolution

Welche genetische Unterschiede machen uns verschieden von unseren nächsten Verwandten, den Schimpansen, und andererseits so ähnlich zu den Schimpansen? Was wir untersuchen und auch verstehen wollen, ist die komplexe Beziehung zwischen den multiplen genetischen und epigenetischen Unterschieden, deren Interaktion mit diversen Umwelt- und Kulturfaktoren in den beobachteten phänotypischen Unterschieden resultieren. Um aufzuklären, ob chromosomale Rearrangements zur Divergenz zwischen Mensch und Schimpanse beigetragen haben und welche selektiven Kräfte ihre Evolution geprägt haben, habe ich die kodierenden Sequenzen von 2 Mb umfassenden, die perizentrischen Inversionsbruchpunkte flankierenden Regionen auf den Chromosomen 1, 4, 5, 9, 12, 17 und 18 untersucht. Als Kontrolle dienten dabei 4 Mb umfassende kollineare Regionen auf den rearrangierten Chromosomen, welche mindestens 10 Mb von den Bruchpunktregionen entfernt lagen. Dabei konnte ich in den Bruchpunkten flankierenden Regionen im Vergleich zu den Kontrollregionen keine höhere Proteinevolutionsrate feststellen. Meine Ergebnisse unterstützen nicht die chromosomale Speziationshypothese für Mensch und Schimpanse, da der Anteil der positiv selektierten Gene (5,1% in den Bruchpunkten flankierenden Regionen und 7% in den Kontrollregionen) in beiden Regionen ähnlich war. Durch den Vergleich der Anzahl der positiv und negativ selektierten Gene per Chromosom konnte ich feststellen, dass Chromosom 9 die meisten und Chromosom 5 die wenigsten positiv selektierten Gene in den Bruchpunkt flankierenden Regionen und Kontrollregionen enthalten. Die Anzahl der negativ selektierten Gene (68) war dabei viel höher als die Anzahl der positiv selektierten Gene (17). Eine bioinformatische Analyse von publizierten Microarray-Expressionsdaten (Affymetrix Chip U95 und U133v2) ergab 31 Gene, die zwischen Mensch und Schimpanse differentiell exprimiert sind. Durch Untersuchung des dN/dS-Verhältnisses dieser 31 Gene konnte ich 7 Gene als negativ selektiert und nur 1 Gen als positiv selektiert identifizieren. Dieser Befund steht im Einklang mit dem Konzept, dass Genexpressionslevel unter stabilisierender Selektion evolvieren. Die meisten positiv selektierten Gene spielen überdies eine Rolle bei der Fortpflanzung. Viele dieser Speziesunterschiede resultieren eher aus Änderungen in der Genregulation als aus strukturellen Änderungen der Genprodukte. Man nimmt an, dass die meisten Unterschiede in der Genregulation sich auf transkriptioneller Ebene manifestieren. Im Rahmen dieser Arbeit wurden die Unterschiede in der DNA-Methylierung zwischen Mensch und Schimpanse untersucht. Dazu wurden die Methylierungsmuster der Promotor-CpG-Inseln von 12 Genen im Cortex von Menschen und Schimpansen mittels klassischer Bisulfit-Sequenzierung und Bisulfit-Pyrosequenzierung analysiert. Die Kandidatengene wurden wegen ihrer differentiellen Expressionsmuster zwischen Mensch und Schimpanse sowie wegen Ihrer Assoziation mit menschlichen Krankheiten oder dem genomischen Imprinting ausgewählt. Mit Ausnahme einiger individueller Positionen zeigte die Mehrzahl der analysierten Gene keine hohe intra- oder interspezifische Variation der DNA-Methylierung zwischen den beiden Spezies. Nur bei einem Gen, CCRK, waren deutliche intraspezifische und interspezifische Unterschiede im Grad der DNA-Methylierung festzustellen. Die differentiell methylierten CpG-Positionen lagen innerhalb eines repetitiven Alu-Sg1-Elements. Die Untersuchung des CCRK-Gens liefert eine umfassende Analyse der intra- und interspezifischen Variabilität der DNA-Methylierung einer Alu-Insertion in eine regulatorische Region. Die beobachteten Speziesunterschiede deuten darauf hin, dass die Methylierungsmuster des CCRK-Gens wahrscheinlich in Adaption an spezifische Anforderungen zur Feinabstimmung der CCRK-Regulation unter positiver Selektion evolvieren. Der Promotor des CCRK-Gens ist anfällig für epigenetische Modifikationen durch DNA-Methylierung, welche zu komplexen Transkriptionsmustern führen können. Durch ihre genomische Mobilität, ihren hohen CpG-Anteil und ihren Einfluss auf die Genexpression sind Alu-Insertionen exzellente Kandidaten für die Förderung von Veränderungen während der Entwicklungsregulation von Primatengenen. Der Vergleich der intra- und interspezifischen Methylierung von spezifischen Alu-Insertionen in anderen Genen und Geweben stellt eine erfolgversprechende Strategie dar.

THEORETICAL STUDIES ON THE METHODS OF RECONSTRUCTING PHYLOGENETIC TREES FROM DNA SEQUENCE DATA (MOLECULAR EVOLUTION, HOMINOID EVOLUTION, COMPUTER SIMULATION)

(1) A mathematical theory for computing the probabilities of various nucleotide configurations is developed, and the probability of obtaining the correct phylogenetic tree (model tree) from sequence data is evaluated for six phylogenetic tree-making methods (UPGMA, distance Wagner method, transformed distance method, Fitch-Margoliash's method, maximum parsimony method, and compatibility method). The number of nucleotides (m*) necessary to obtain the correct tree with a probability of 95% is estimated with special reference to the human, chimpanzee, and gorilla divergence. m* is at least 4,200, but the availability of outgroup species greatly reduces m* for all methods except UPGMA. m* increases if transitions occur more frequently than transversions as in the case of mitochondrial DNA. (2) A new tree-making method called the neighbor-joining method is proposed. This method is applicable either for distance data or character state data. Computer simulation has shown that the neighbor-joining method is generally better than UPGMA, Farris' method, Li's method, and modified Farris method on recovering the true topology when distance data are used. A related method, the simultaneous partitioning method, is also discussed. (3) The maximum likelihood (ML) method for phylogeny reconstruction under the assumption of both constant and varying evolutionary rates is studied, and a new algorithm for obtaining the ML tree is presented. This method gives a tree similar to that obtained by UPGMA when constant evolutionary rate is assumed, whereas it gives a tree similar to that obtained by the maximum parsimony tree and the neighbor-joining method when varying evolutionary rate is assumed. ^

Evolution of cis-regulatory elements in a repetitive sequence adjacent to a sea urchin aboral ectoderm -specific gene

The creation, preservation, and degeneration of cis-regulatory elements controlling developmental gene expression are fundamental genome-level evolutionary processes about which little is known. In this study, critical differences in cis-regulatory elements controlling the expression of the sea urchin aboral ectoderm-specific spec genes were identified and explored. In genomes of species within the Strongylocentrotidae family, multiple copies of a repetitive sequence element termed RSR were present, but RSRs were not detected in genomes of species outside Strongylocentrotidae. RSRs are invariably associated with spec genes, and in Strongylocentrotus purpuratus, the spec2a RSR functioned as a transcriptional enhancer displaying greater activity than RSRs from the spec1 or spec2c paralogs. Single base-pair differences at two cis-regulatory elements within the spec2a RSR greatly increased the binding affinities of four transcription factors: SpCCAAT-binding factor at one element and SpOtx, SpGoosecoid, and SpGATA-E at another. The cis-regulatory elements to which SpCCAAT-binding factor, SpOtx, SpGoosecoid, and SpGATA-E bound were recent evolutionary acquisitions that could act either to activate or repress transcription, depending on the cell type. These elements were found in the spec2a RSR ortholog in Strongylocentrotus pallidus but not in the RSR orthologs of Strongylocentrotus droebachiensis or Hemicentrotus pulcherrimus. These results indicate that spec genes exhibit a dynamic pattern of cis-regulatory element evolution while stabilizing selection preserves their aboral ectoderm expression domain. ^

Striking sequence similarity in inter- and intra-specific comparisons of class I SLG alleles from Brassica oleracea and Brassica campestris: Implications for the evolution and recognition mechanism

Self-incompatibility in Brassica is controlled by a single multi-allelic locus (S locus), which contains at least two highly polymorphic genes expressed in the stigma: an S glycoprotein gene (SLG) and an S receptor kinase gene (SRK). The putative ligand-binding domain of SRK exhibits high homology to the secretory protein SLG, and it is believed that SLG and SRK form an active receptor kinase complex with a self-pollen ligand, which leads to the rejection of self-pollen. Here, we report 31 novel SLG sequences of Brassica oleracea and Brassica campestris. Sequence comparisons of a large number of SLG alleles and SLG-related genes revealed the following points. (i) The striking sequence similarity observed in an inter-specific comparison (95.6% identity between SLG14 of B. oleracea and SLG25 of B. campestris in deduced amino acid sequence) suggests that SLG diversification predates speciation. (ii) A perfect match of the sequences in hypervariable regions, which are thought to determine S specificity in an intra-specific comparison (SLG8 and SLG46 of B. campestris) and the observation that the hypervariable regions of SLG and SRK of the same S haplotype were not necessarily highly similar suggests that SLG and SRK bind different sites of the pollen ligand and that they together determine S specificity. (iii) Comparison of the hypervariable regions of SLG alleles suggests that intragenic recombination, together with point mutations, has contributed to the generation of the high level of sequence variation in SLG alleles. Models for the evolution of SLG/SRK are presented.