Reconstruction of the vertebrate ancestral genome reveals dynamic genome reorganization in early vertebrates

Although several vertebrate genomes have been sequenced, little is known about the genome evolution of early vertebrates and how large-scale genomic changes such as the two rounds of whole-genome duplications (2R WGD) affected evolutionary complexity and novelty in vertebrates. Reconstructing the ancestral vertebrate genome is highly nontrivial because of the difficulty in identifying traces originating from the 2R WGD. To resolve this problem, we developed a novel method capable of pinning down remains of the 2R WGD in the human and medaka fish genomes using invertebrate tunicate and sea urchin genes to define ohnologs, i.e., paralogs produced by the 2R WGD. We validated the reconstruction using the chicken genome, which was not considered in the reconstruction step, and observed that many ancestral proto-chromosomes were retained in the chicken genome and had one-to-one correspondence to chicken microchromosomes, thereby confirming the reconstructed ancestral genomes. Our reconstruction revealed a contrast between the slow karyotype evolution after the second WGD and the rapid, lineage-specific genome reorganizations that occurred in the ancestral lineages of major taxonomic groups such as teleost fishes, amphibians, reptiles, and marsupials.     

Multilocus sequence typing of hospital-associated Enterococcus faecium from Brazil reveals their unique evolutionary history

We studied the genetic relationships between vancomycin-susceptible (n = 11) and -resistant Enterococcus faecium (VRE, n = 20) recovered from Brazil using a multilocus sequence typing (MLST) scheme. Grouping of allelic profiles revealed six clusters of related sequence types (STs) that differ in no more than two of the seven alleles. Of these, one cluster harbored 16 of the 20 isolates recovered during the first VRE outbreak in Brazil. The ampicillin and gentamicin resistance profiles were stable in the isolates that clustered within the groups I-III. Comparison with the allelic profiles of 139 E. faecium from different geographical regions and origins found in the international database http://www.mlst.net revealed that the Brazilian outbreak clone did not cluster in the previously named complex-17. This genetic complex contains hospital epidemic and clinical isolates recovered from different countries and continents. Twenty two of the 31 Brazilian isolates, including the VRE outbreak clone, clustered apart from the E. faecium isolates from the database, suggesting that these Brazilian isolates have a distinct evolutionary history.     

The complete sequence of the zebrafish (Danio rerio) mitochondrial genome and evolutionary patterns in vertebrate mitochondrial DNA.

We describe the complete sequence of the 16,596-nucleotide mitochondrial genome of the zebrafish (Danio rerio); contained are 13 protein genes, 22 tRNAs, 2 rRNAs, and a noncoding control region. Codon usage in protein genes is generally biased toward the available tRNA species but also reflects strand-specific nucleotide frequencies. For 19 of the 20 amino acids, the most frequently used codon ends in either A or C, with A preferred over C for fourfold degenerate codons (the lone exception was AUG: methionine). We show that rates of sequence evolution vary nearly as much within vertebrate classes as between them, yet nucleotide and amino acid composition show directional evolutionary trends, including marked differences between mammals and all other taxa. Birds showed similar compositional characteristics to the other nonmammalian taxa, indicating that the evolutionary trend in mammals is not solely due to metabolic rate and thermoregulatory factors. Complete mitochondrial genomes provide a large character base for phylogenetic analysis and may provide for robust estimates of phylogeny. Phylogenetic analysis of zebrafish and 35 other taxa based on all protein-coding genes produced trees largely, but not completely, consistent with conventional views of vertebrate evolution. It appears that even with such a large number of nucleotide characters (11,592), limited taxon sampling can lead to problems associated with extensive evolution on long phyletic branches.     

Gorilla genome structural variation reveals evolutionary parallelisms with chimpanzee

Structural variation has played an important role in the evolutionary restructuring of human and great ape genomes. Recent analyses have suggested that the genomes of chimpanzee and human have been particularly enriched for this form of genetic variation. Here, we set out to assess the extent of structural variation in the gorilla lineage by generating 10-fold genomic sequence coverage from a western lowland gorilla and integrating these data into a physical and cytogenetic framework of structural variation. We discovered and validated over 7665 structural changes within the gorilla lineage, including sequence resolution of inversions, deletions, duplications, and mobile element insertions. A comparison with human and other ape genomes shows that the gorilla genome has been subjected to the highest rate of segmental duplication. We show that both the gorilla and chimpanzee genomes have experienced independent yet convergent patterns of structural mutation that have not occurred in humans, including the formation of subtelomeric heterochromatic caps, the hyperexpansion of segmental duplications, and bursts of retroviral integrations. Our analysis suggests that the chimpanzee and gorilla genomes are structurally more derived than either orangutan or human genomes.     

Whole-genome analysis of Alu repeat elements reveals complex evolutionary history

Alu repeats are the most abundant family of repeats in the human genome, with over 1 million copies comprising 10% of the genome. They have been implicated in human genetic disease and in the enrichment of gene-rich segmental duplications in the human genome, and they form a rich fossil record of primate and human history. Alu repeat elements are believed to have arisen from the replication of a small number of source elements, whose evolution over time gives rise to the 31 Alu subfamilies currently reported in Repbase Update. We apply a novel method to identify and statistically validate 213 Alu subfamilies. We build an evolutionary tree of these subfamilies and conclude that the history of Alu evolution is more complex than previous studies had indicated.     

Systematic classification of vertebrate chemokines based on conserved synteny and evolutionary history

The genes involved in host defences are known to undergo rapid evolution. Therefore, it is often difficult to assign orthologs in multigene families among various vertebrate species. Chemokines are a large family of small cytokines that orchestrate cell migration in health and disease. Herein, we have surveyed the genomes of 18 representative vertebrate species for chemokine genes and identified a total of 553 genes. We have determined their orthologous relationships and classified them in accordance with the current systematic chemokine nomenclature system. Our study reveals an interesting evolutionary history that gave origin and diversification to the vertebrate chemokine superfamily.     

Protein structure and evolutionary history determine sequence space topology

Understanding the observed variability in the number of homologs of a gene is a very important unsolved problem that has broad implications for research into coevolution of structure and function, gene duplication, pseudogene formation, and possibly for emerging diseases. Here, we attempt to define and elucidate some possible causes behind the observed irregularity in sequence space. We present evidence that sequence variability and functional diversity of a gene or fold family is influenced by quantifiable characteristics of the protein structure. These characteristics reflect the structural potential for sequence plasticity, i.e., the ability to accept mutation without losing thermodynamic stability. We identify a structural feature of a protein domain-contact density-that serves as a determinant of entropy in sequence space, i.e., the ability of a protein to accept mutations without destroying the fold (also known as fold designability). We show that (log) of average gene family size exhibits statistical correlation (R(2) > 0.9.) with contact density of its three-dimensional structure. We present evidence that the size of individual gene families are influenced not only by the designability of the structure, but also by evolutionary history, e.g., the amount of time the gene family was in existence. We further show that our observed statistical correlation between gene family size and contact density of the structure is valid on many levels of evolutionary divergence, i.e., not only for closely related sequence, but also for less-related fold and superfamily levels of homology.     

Mitochondrial genome variation and evolutionary history of Australian and New Guinean aborigines

To study the evolutionary history of the Australian and New Guinean indigenous peoples, we analyzed 101 complete mitochondrial genomes including populations from Australia and New Guinea as well as from Africa, India, Europe, Asia, Melanesia, and Polynesia. The genetic diversity of the Australian mitochondrial sequences is remarkably high and is similar to that found across Asia. This is in contrast to the pattern seen in previously described Y-chromosome data where an Australia-specific haplotype was found at high frequency. The mitochondrial genome data indicate that Australia was colonized between 40 and 70 thousand years ago, either by a single migration from a heterogeneous source population or by multiple movements of smaller groups occurring over a period of time. Some Australian and New Guinea sequences form clades, suggesting the possibility of a joint colonization and/or admixture between the two regions.     

FBR murine osteosarcoma virus. II. Nucleotide sequence of the provirus reveals that the genome contains sequences acquired from two cellular genes

The complete nucleotide sequence of the FBR proviral DNA has been determined. The provirus of 3791 nucleotides (specifying a genome of 3284 bases) encodes a single gag- fos fusion product of 554 amino acids. The fos portion of the gene lacks the sequences which code for the first 24 and the last 98 amino acids of the 380-amino acid mouse c- fos gene product. In addition, the coding region has sustained three in-frame deletions, one in the p30gag portion, and two in the fos region, as compared to the sequences of AKR-MLV and the c- fos gene, respectively. The gene product terminates in sequences, termed v-fox, that are present in uninfected mouse DNA at loci unrelated to the c- fos gene. The c-fox gene(s) is expressed as an abundant class of polyadenylated RNA in normal mouse tissues.     

Molecular characterization of a large group of Mucopolysaccharidosis type IIIC patients reveals the evolutionary history of the disease

Mucopolysaccharidosis type IIIC (MPSIIIC) is a severe, rare autosomal recessive disorder caused by variants in the heparan‐α‐glucosaminide N‐acetyltransferase (HGSNAT) gene which result in lysosomal accumulation of heparan sulfate. We analyzed clinical presentation, molecular defects and their haplotype context in 78 (27 novel) MPSIIIC cases from 22 countries, the largest group studied so far. We describe for the first time disease‐causing variants in the patients from Brazil, Algeria, Azerbaijan, and Iran, and extend their spectrum within Canada, Colombia, Turkey, and the USA. Six variants are novel: two missense, c.773A>T/p.N258I and c.1267G>T/p.G423W, a nonsense c.164T>A/p.L55*, a splice‐site mutation c.494?1G>A/p.[P165_L187delinsQSCYVTQAGVRWHHLGSLQALPPGFTPFSYLSLLSSWNC,P165fs], a deletion c.1348delG/p.(D450fs) and an insertion c.1479dupA/p.(Leu494fs). The missense HGSNAT variants lacked lysosomal targeting, enzymatic activity, and likely the correct folding. The haplotype analysis identified founder mutations, p.N258I, c.525dupT, and p.L55* in the Brazilian state of Paraiba, c.493+1G>A in Eastern Canada/Quebec, p.A489E in the USA, p.R384* in Poland, p.R344C and p.S518F in the Netherlands and suggested that variants c.525dupT, c.372?2G>A, and c.234+1G>A present in cis with c.564‐98T>C and c.710C>A rare single‐nucleotide polymorphisms, have been introduced by Portuguese settlers in Brazil. Altogether, our results provide insights into the origin, migration roots and founder effects of HGSNAT disease‐causing variants, and reveal the evolutionary history of MPSIIIC.     

Complete genome sequence of avian paramyxovirus type 3 reveals an unusually long trailer region

The complete genome sequence was determined for prototype parakeet/Netherlands/449/75 strain of avian paramyxovirus (APMV) serotype 3. The genome is 16,272 nucleotides (nt) in length, consisting of six non-overlapping genes in the order of 3'-N-P/V/W-M-F-HN-L-5', with intergenic regions of 31-63nt. APMV-3 genome follows the "rule of six" and is the largest among the avian paramyxoviruses reported to date, with a trailer region of 707nt, the longest in the family Paramyxoviridae. The cleavage site of F protein, A-R-P-R-G-R downward arrowL, does not conform to the preferred cleavage site of the ubiquitous cellular protease furin. Therefore, exogenous protease was needed for replication in vitro. Alignment and phylogenetic analysis of the predicted amino acid sequences of strain Netherlands proteins with the cognate proteins of viruses of all of the five genera of family Paramyxoviridae showed that APMV-3 strain Netherlands is more closely related to APMV-1 than APMV-6.     

Comparative analysis of genome sequences of three isolates of Orf virus reveals unexpected sequence variation

Orf virus (ORFV) is the type species of the Parapoxvirus genus. Here, we present the genomic sequence of the most well studied ORFV isolate, strain NZ2. The NZ2 genome is 138 kbp and contains 132 putative genes, 88 of which are present in all analyzed chordopoxviruses. Comparison of the NZ2 genome with the genomes of 2 other fully sequenced isolates of ORFV revealed that all 3 genomes carry each of the 132 genes, but there are substantial sequence variations between isolates in a significant number of genes, including 9 with inter-isolate amino acid sequence identity of only 38-79%. Each genome has an average of 64% G+C but each has a distinctive pattern of substantial deviation from the average within particular regions of the genome. The same pattern of variation was also seen in the genome of another parapoxvirus species and was clearly unlike the uniform patterns of G+C content seen in all other genera of chordopoxviruses. The availability of genomic sequences of three orf virus isolates allowed us to more accurately assess likely coding regions and thereby revise published data for 24 genes and to predict two previously unrecognized genes.     

The complete genome sequence of an enterovirus 76 isolate in China reveals a recombination event

Enterovirus 76 (EV76) is a new member of species Human Enterovirus A (HEV-A). So far, the only complete genome sequence of the prototype strain from France has been available. In this study, we determined the complete nucleotide sequence of strain 04360, isolated from an acute flaccid paralysis patient in Shandong province, China in 2004. Sequence analysis revealed 80.7-94.7% VP1 nucleotide identity to other EV76 strains and provided evidence of recombination with other types of HEV-A in the P2 and P3 coding regions.     

Decoding the human genome sequence

The year 2000 is marked by the production of the sequence of the human genome. A 'working draft' of high quality sequence covering 90% of the genome has been determined and a quarter is in finished form, including the first two completed chromosomes. All sequence data from the project is made freely available to the community via the Internet, for further analysis and exploitation. The challenge which lies ahead is to decipher the information. Knowledge of the human genome sequence will enable us to understand how the genetic information determines the development, structure and function of the human body. We will be able to explore how variations within our DNA sequence cause disease, how they affect our interaction with our environment and ultimately to develop new and effective ways to improve human health.     

The Chlamydophila abortus genome sequence reveals an array of variable proteins that contribute to interspecies variation

The obligate intracellular bacterial pathogen Chlamydophila abortus strain S26/3 (formerly the abortion subtype of Chlamydia psittaci) is an important cause of late gestation abortions in ruminants and pigs. Furthermore, although relatively rare, zoonotic infection can result in acute illness and miscarriage in pregnant women. The complete genome sequence was determined and shows a high level of conservation in both sequence and overall gene content in comparison to other Chlamydiaceae. The 1,144,377-bp genome contains 961 predicted coding sequences, 842 of which are conserved with those of Chlamydophila caviae and Chlamydophila pneumoniae. Within this conserved Cp. abortus core genome we have identified the major regions of variation and have focused our analysis on these loci, several of which were found to encode highly variable protein families, such as TMH/Inc and Pmp families, which are strong candidates for the source of diversity in host tropism and disease causation in this group of organisms. Significantly, Cp. abortus lacks any toxin genes, and also lacks genes involved in tryptophan metabolism and nucleotide salvaging (guaB is present as a pseudogene), suggesting that the genetic basis of niche adaptation of this species is distinct from those previously proposed for other chlamydial species.     

The complete mitochondrial DNA sequence of the green alga Oltmannsiellopsis viridis: evolutionary trends of the mitochondrial genome in the Ulvophyceae

The mitochondrial genome displays a highly plastic architecture in the green algal division comprising the classes Prasinophyceae, Trebouxiophyceae, Ulvophyceae, and Chlorophyceae (Chlorophyta). The compact mitochondrial DNAs (mtDNAs) of Nephroselmis (Prasinophyceae) and Prototheca (Trebouxiophyceae) encode about 60 genes and have been ascribed an ‘ancestral’ pattern of evolution, whereas those of chlorophycean green algae are much more reduced in gene content and size. Although the mtDNA of the early-diverging ulvophyte Pseudendoclonium contains 57 conserved genes, it differs from ‘ancestral’ chlorophyte mtDNAs by its unusually large size (96 kb) and long intergenic spacers. To gain insights into the evolutionary trends of mtDNA in the Ulvophyceae, we have determined the complete mtDNA sequence of Oltmannsiellopsis viridis, an ulvophyte belonging to a distinct, early-diverging lineage. This 56,761 bp genome harbours 54 conserved genes, numerous repeated sequences, and only three introns. From our comparative analyses with Pseudendoclonium mtDNA, we infer that the mitochondrial genome of the last common ancestor of the two ulvophytes closely resembled that of the trebouxiophyte Prototheca in terms of gene content and gene density. Our results also provide strong evidence for the intracellular, interorganellar transfer of a group I intron and for two distinct events of intercellular, horizontal DNA transfer.Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

Comparative genome analyses of Arabidopsis spp.: inferring chromosomal rearrangement events in the evolutionary history of A. thaliana

Comparative genome analysis is a powerful tool that can facilitate the reconstruction of the evolutionary history of the genomes of modern-day species. The model plant Arabidopsis thaliana with its n = 5 genome is thought to be derived from an ancestral n = 8 genome. Pairwise comparative genome analyses of A. thaliana with polyploid and diploid Brassicaceae species have suggested that rapid genome evolution, manifested by chromosomal rearrangements and duplications, characterizes the polyploid, but not the diploid, lineages of this family. In this study, we constructed a low-density genetic linkage map of Arabidopsis lyrata ssp. lyrata (A. l. lyrata; n = 8, diploid), the closest known relative of A. thaliana (MRCA approximately 5 Mya), using A. thaliana-specific markers that resolve into the expected eight linkage groups. We then performed comparative Bayesian analyses using raw mapping data from this study and from a Capsella study to infer the number and nature of rearrangements that distinguish the n = 8 genomes of A. l. lyrata and Capsella from the n = 5 genome of A. thaliana. We conclude that there is strong statistical support in favor of the parsimony scenarios of 10 major chromosomal rearrangements separating these n = 8 genomes from A. thaliana. These chromosomal rearrangement events contribute to a rate of chromosomal evolution higher than previously reported in this lineage. We infer that at least seven of these events, common to both sets of data, are responsible for the change in karyotype and underlie genome reduction in A. thaliana.