Two new avian mitochondrial genomes (penguin and goose) and a summary of bird and reptile mitogenomic features

Dynamic nucleotide mutation gradients and control region usage in squamate reptile mitochondrial genomes

Gradients of nucleotide bias and substitution rates occur in vertebrate mitochondrial genomes due to the asymmetric nature of the replication process. The evolution of these gradients has previously been studied in detail in primates, but not in other vertebrate groups. From the primate study, the strengths of these gradients are known to evolve in ways that can substantially alter the substitution process, but it is unclear how rapidly they evolve over evolutionary time or how different they may be in different lineages or groups of vertebrates. Given the importance of mitochondrial genomes in phylogenetics and molecular evolutionary research, a better understanding of how asymmetric mitochondrial substitution gradients evolve would contribute key insights into how this gradient evolution may mislead evolutionary inferences, and how it may also be incorporated into new evolutionary models. Most snake mitochondrial genomes have an additional interesting feature, 2 nearly identical control regions, which vary among different species in the extent that they are used as origins of replication. Given the expanded sampling of complete snake genomes currently available, together with 2 additional snakes sequenced in this study, we reexamined gradient strength and CR usage in alethinophidian snakes as well as several lizards that possess dual CRs. Our results suggest that nucleotide substitution gradients (and corresponding nucleotide bias) and CR usage is highly labile over the approximately 200 m.y. of squamate evolution, and demonstrates greater overall variability than previously shown in primates. The evidence for the existence of such gradients, and their ability to evolve rapidly and converge among unrelated species suggests that gradient dynamics could easily mislead phylogenetic and molecular evolutionary inferences, and argues strongly that these dynamics should be incorporated into phylogenetic models.

The complete mitochondrial genomes of the whistling duck (Dendrocygna javanica) and black swan (Cygnus atratus): dating evolutionary divergence in Galloanserae

Galloanserae is an ancient and diverse avian group, for which comprehensive molecular evidence relevant to phylogenetic analysis in the context of molecular chronology is lacking. In this study, we present two additional mitochondrial genome sequences of Galloanserae (the whistling duck, Dendrocygna javanica, and the black swan, Cygnus atratus) to broaden the scope of molecular phylogenetic reconstruction. The lengths of the whistling duck's and black swan's mitochondrial genomes are 16,753 and 16,748 bases, respectively. Phylogenetic analyses suggest that Dendrocygna is more likely to be in a basal position of the branch consisting of Anatinae and Anserinae, an affiliation that does not conform to its traditional classification. Bayesian approaches were employed to provide a rough timescale for Galloanserae evolution. In general, a narrow range of 95% confidence intervals gave younger estimates than those based on limited genes and estimated that at least two lineages originated before the Coniacian epoch around 90 MYA, well before the Cretaceous-Tertiary boundary. In addition, these results, which were compatible with estimates from fossil evidence, also imply that the origin of numerous genera in Anseriformes took place in the late Oligocene to early Miocene. Taken together, the results presented here provide a working framework for future research on Galloanserae evolution, and they underline the utility of whole mitochondrial genome sequences for the resolution of deep divergence.

Toward resolving deep neoaves phylogeny: data, signal enhancement, and priors

We report three developments toward resolving the challenge of the apparent basal polytomy of neoavian birds. First, we describe improved conditional down-weighting techniques to reduce noise relative to signal for deeper divergences and find increased agreement between data sets. Second, we present formulae for calculating the probabilities of finding predefined groupings in the optimal tree. Finally, we report a significant increase in data: nine new mitochondrial (mt) genomes (the dollarbird, New Zealand kingfisher, great potoo, Australian owlet-nightjar, white-tailed trogon, barn owl, a roadrunner [a ground cuckoo], New Zealand long-tailed cuckoo, and the peach-faced lovebird) and together they provide data for each of the six main groups of Neoaves proposed by Cracraft J (2001). We use his six main groups of modern birds as priors for evaluation of results. These include passerines, cuckoos, parrots, and three other groups termed "WoodKing" (woodpeckers/rollers/kingfishers), "SCA" (owls/potoos/owlet-nightjars/hummingbirds/swifts), and "Conglomerati." In general, the support is highly significant with just two exceptions, the owls move from the "SCA" group to the raptors, particularly accipitrids (buzzards/eagles) and the osprey, and the shorebirds may be an independent group from the rest of the "Conglomerati". Molecular dating mt genomes support a major diversification of at least 12 neoavian lineages in the Late Cretaceous. Our results form a basis for further testing with both nuclear-coding sequences and rare genomic changes.

Adaptive evolution and functional redesign of core metabolic proteins in snakes

BACKGROUND

Adaptive evolutionary episodes in core metabolic proteins are uncommon, and are even more rarely linked to major macroevolutionary shifts.

METHODOLOGY/PRINCIPAL FINDINGS

We conducted extensive molecular evolutionary analyses on snake mitochondrial proteins and discovered multiple lines of evidence suggesting that the proteins at the core of aerobic metabolism in snakes have undergone remarkably large episodic bursts of adaptive change. We show that snake mitochondrial proteins experienced unprecedented levels of positive selection, coevolution, convergence, and reversion at functionally critical residues. We examined Cytochrome C oxidase subunit I (COI) in detail, and show that it experienced extensive modification of normally conserved residues involved in proton transport and delivery of electrons and oxygen. Thus, adaptive changes likely altered the flow of protons and other aspects of function in CO, thereby influencing fundamental characteristics of aerobic metabolism. We refer to these processes as "evolutionary redesign" because of the magnitude of the episodic bursts and the degree to which they affected core functional residues.

CONCLUSIONS/SIGNIFICANCE

The evolutionary redesign of snake COI coincided with adaptive bursts in other mitochondrial proteins and substantial changes in mitochondrial genome structure. It also generally coincided with or preceded major shifts in ecological niche and the evolution of extensive physiological adaptations related to lung reduction, large prey consumption, and venom evolution. The parallel timing of these major evolutionary events suggests that evolutionary redesign of metabolic and mitochondrial function may be related to, or underlie, the extreme changes in physiological and metabolic efficiency, flexibility, and innovation observed in snake evolution.

Bird evolution: testing the Metaves clade with six new mitochondrial genomes

BACKGROUND

Evolutionary biologists are often misled by convergence of morphology and this has been common in the study of bird evolution. However, the use of molecular data sets have their own problems and phylogenies based on short DNA sequences have the potential to mislead us too. The relationships among clades and timing of the evolution of modern birds (Neoaves) has not yet been well resolved. Evidence of convergence of morphology remain controversial. With six new bird mitochondrial genomes (hummingbird, swift, kagu, rail, flamingo and grebe) we test the proposed Metaves/Coronaves division within Neoaves and the parallel radiations in this primary avian clade.

RESULTS

Our mitochondrial trees did not return the Metaves clade that had been proposed based on one nuclear intron sequence. We suggest that the high number of indels within the seventh intron of the beta-fibrinogen gene at this phylogenetic level, which left a dataset with not a single site across the alignment shared by all taxa, resulted in artifacts during analysis. With respect to the overall avian tree, we find the flamingo and grebe are sister taxa and basal to the shorebirds (Charadriiformes). Using a novel site-stripping technique for noise-reduction we found this relationship to be stable. The hummingbird/swift clade is outside the large and very diverse group of raptors, shore and sea birds. Unexpectedly the kagu is not closely related to the rail in our analysis, but because neither the kagu nor the rail have close affinity to any taxa within this dataset of 41 birds, their placement is not yet resolved.

CONCLUSION

Our phylogenetic hypothesis based on 41 avian mitochondrial genomes (13,229 bp) rejects monophyly of seven Metaves species and we therefore conclude that the members of Metaves do not share a common evolutionary history within the Neoaves.

Cost minimization of ribosomal frameshifts

Properties of mRNA leading regions that modulate protein synthesis are little known (besides effects of their secondary structure). Here I explore how coding properties of leading regions may account for their disparate efficiencies. Trinucleotides that form off frame stop codons decrease costs of ribosomal slippages during protein synthesis: protein activity (as a proxy of gene expression, and as measured in experiments using artificial variants of 5' leading sequences of beta galactosidase in Escherichia coli) increases proportionally to the number of stop motifs in any frame in the 5' leading region. This suggests that stop codons in the 5' leading region, upstream of the recognized coding sequence, terminate eventual translations that sometimes start before ribosomes reach the mRNA's recognized start codon, increasing efficiency. This hypothesis is confirmed by further analyses: mRNAs with 5' leading regions containing in the same frame a start preceding a stop codon (in any frame) produce less enzymatic activity than those with the stop preceding the start. Hence coding properties, in addition to other properties, such as the secondary structure of the 5' leading region, regulate translation. This experimentally (a) confirms that within coding regions, off frame stops increase protein synthesis efficiency by early stopping frameshifted translation; (b) suggests that this occurs for all frames also in 5' leading regions and that (c) several alternative start codons that function at different probabilities should routinely be considered for all genes in the region of the recognized initiation codon. An unknown number of short peptides might be translated from coding and non-coding regions of RNAs.

The problem of rooting rapid radiations

There are many examples of groups (such as birds, bees, mammals, multicellular animals, and flowering plants) that have undergone a rapid radiation. In such cases, where there is a combination of short internal and long external branches, correctly estimating and rooting phylogenetic trees is known to be a difficult problem. In this simulation study, we tested the performances of different phylogenetic methods at estimating a tree that models a rapid radiation. We found that maximum likelihood, corrected and uncorrected neighbor-joining, and corrected and uncorrected parsimony, all suffer from biases toward specific tree topologies. In addition, we found that using a single-taxon outgroup to root a tree frequently disrupts an otherwise correct ingroup phylogeny. Moreover, for uncorrected parsimony, we found cases where several individual trees (in which the outgroup was placed incorrectly) were selected more frequently than the correct tree. Even for parameter settings where the correct tree was selected most frequently when using extremely long sequences, for sequences of up to 60,000 nucleotides the incorrectly rooted trees were each selected more frequently than the correct tree. For all the cases tested here, tree estimation using a two taxon outgroup was more accurate than when using a single-taxon outgroup. However, the ingroup was most accurately recovered when no outgroup was used.

Comparative mitochondrial genomics of snakes: extraordinary substitution rate dynamics and functionality of the duplicate control region

BACKGROUND

The mitochondrial genomes of snakes are characterized by an overall evolutionary rate that appears to be one of the most accelerated among vertebrates. They also possess other unusual features, including short tRNAs and other genes, and a duplicated control region that has been stably maintained since it originated more than 70 million years ago. Here, we provide a detailed analysis of evolutionary dynamics in snake mitochondrial genomes to better understand the basis of these extreme characteristics, and to explore the relationship between mitochondrial genome molecular evolution, genome architecture, and molecular function. We sequenced complete mitochondrial genomes from Slowinski's corn snake (Pantherophis slowinskii) and two cottonmouths (Agkistrodon piscivorus) to complement previously existing mitochondrial genomes, and to provide an improved comparative view of how genome architecture affects molecular evolution at contrasting levels of divergence.

RESULTS

We present a Bayesian genetic approach that suggests that the duplicated control region can function as an additional origin of heavy strand replication. The two control regions also appear to have different intra-specific versus inter-specific evolutionary dynamics that may be associated with complex modes of concerted evolution. We find that different genomic regions have experienced substantial accelerated evolution along early branches in snakes, with different genes having experienced dramatic accelerations along specific branches. Some of these accelerations appear to coincide with, or subsequent to, the shortening of various mitochondrial genes and the duplication of the control region and flanking tRNAs.

CONCLUSION

Fluctuations in the strength and pattern of selection during snake evolution have had widely varying gene-specific effects on substitution rates, and these rate accelerations may have been functionally related to unusual changes in genomic architecture. The among-lineage and among-gene variation in rate dynamics observed in snakes is the most extreme thus far observed in animal genomes, and provides an important study system for further evaluating the biochemical and physiological basis of evolutionary pressures in vertebrate mitochondria.

Higher-order phylogeny of modern birds (Theropoda, Aves: Neornithes) based on comparative anatomy. II. Analysis and discussion

In recent years, avian systematics has been characterized by a diminished reliance on morphological cladistics of modern taxa, intensive palaeornithogical research stimulated by new discoveries and an inundation by analyses based on DNA sequences. Unfortunately, in contrast to significant insights into basal origins, the broad picture of neornithine phylogeny remains largely unresolved. Morphological studies have emphasized characters of use in palaeontological contexts. Molecular studies, following disillusionment with the pioneering, but non-cladistic, work of Sibley and Ahlquist, have differed markedly from each other and from morphological works in both methods and findings. Consequently, at the turn of the millennium, points of robust agreement among schools concerning higher-order neornithine phylogeny have been limited to the two basalmost and several mid-level, primary groups. This paper describes a phylogenetic (cladistic) analysis of 150 taxa of Neornithes, including exemplars from all non-passeriform families, and subordinal representatives of Passeriformes. Thirty-five outgroup taxa encompassing Crocodylia, predominately theropod Dinosauria, and selected Mesozoic birds were used to root the trees. Based on study of specimens and the literature, 2954 morphological characters were defined; these characters have been described in a companion work, approximately one-third of which were multistate (i.e. comprised at least three states), and states within more than one-half of these multistate characters were ordered for analysis. Complete heuristic searches using 10 000 random-addition replicates recovered a total solution set of 97 well-resolved, most-parsimonious trees (MPTs). The set of MPTs was confirmed by an expanded heuristic search based on 10 000 random-addition replicates and a full ratchet-augmented exploration to ascertain global optima. A strict consensus tree of MPTs included only six trichotomies, i.e. nodes differing topologically among MPTs. Bootstrapping (based on 10 000 replicates) percentages and ratchet-minimized support (Bremer) indices indicated most nodes to be robust. Several fossil Neornithes (e.g. Dinornithiformes, Aepyornithiformes) were placed within the ingroup a posteriori either through unconstrained, heursitic searches based on the complete matrix augmented by these taxa separately or using backbone-constraints. Analysis confirmed the topology among outgroup Theropoda and achieved robust resolution at virtually all levels of the Neornithes. Findings included monophyly of the palaeognathous birds, comprising the sister taxa Tinamiformes and ratites, respectively, and the Anseriformes and Galliformes as monophyletic sister-groups, together forming the sister-group to other Neornithes exclusive of the Palaeognathae (Neoaves). Noteworthy inferences include: (i) the sister-group to remaining Neoaves comprises a diversity of marine and wading birds; (ii) Podicipedidae are the sister-group of Gaviidae, and not closely related to the Phoenicopteridae, as recently suggested; (iii) the traditional Pelecaniformes, including the shoebill (Balaeniceps rex) as sister-taxon to other members, are monophyletic; (iv) traditional Ciconiiformes are monophyletic; (v) Strigiformes and Falconiformes are sister-groups; (vi) Cathartidae is the sister-group of the remaining Falconiformes; (vii) Ralliformes (Rallidae and Heliornithidae) are the sister-group to the monophyletic Charadriiformes, with the traditionally composed Gruiformes and Turniciformes (Turnicidae and Mesitornithidae) sequentially paraphyletic to the entire foregoing clade; (viii) Opisthocomus hoazin is the sister-taxon to the Cuculiformes (including the Musophagidae); (ix) traditional Caprimulgiformes are monophyletic and the sister-group of the Apodiformes; (x) Trogoniformes are the sister-group of Coliiformes; (xi) Coraciiformes, Piciformes and Passeriformes are mutually monophyletic and closely related; and (xii) the Galbulae are retained within the Piciformes. Unresolved portions of the Neornithes (nodes having more than one most-parsimonious solution) comprised three parts of the tree: (a) several interfamilial nodes within the Charadriiformes; (b) a trichotomy comprising the (i) Psittaciformes, (ii) Columbiformes and (iii) Trogonomorphae (Trogoniformes, Coliiformes) + Passerimorphae (Coraciiformes, Piciformes, Passeriformes); and (c) a trichotomy comprising the Coraciiformes, Piciformes and Passeriformes. The remaining polytomies were among outgroups, although several of the highest-order nodes were only marginally supported; however, the majority of nodes were resolved and met or surpassed conventional standards of support. Quantitative comparisons with alternative hypotheses, examination of highly supportive and diagnostic characters for higher taxa, correspondences with prior studies, complementarity and philosophical differences with palaeontological phylogenetics, promises and challenges of palaeogeography and calibration of evolutionary rates of birds, and classes of promising evidence and future directions of study are reviewed. Homology, as applied to avian examples of apparent homologues, is considered in terms of recent theory, and a revised annotated classification of higher-order taxa of Neornithes and other closely related Theropoda is proposed. (c) 2007 The Linnean Society of London, Zoological Journal of the Linnean Society, 2007, 149, 1-95.

The mitochondrial genome of Cygnus columbianus, the Whistling Swan

The nucleotide sequence of the mitochondrial genome of the Whistling Swan, Cygnus columbianus, is reported. Many of the features common to avian mitochondrial genomes are present in C. columbianus and are described here. The gene order is the same as in Gallus gallus. The sequence of this mitochondrial genome allows relationships within the family Anatidae (swans, geese and ducks) to be reconsidered in the light of a large suite of mitochondrial characters. Protein-coding gene sequences of C. columbianus were concatenated to form a supergene, which was analyzed phylogenetically with similar constructs from previously published avian genomes. Relationships within Anatidae and between the Anatidae and the galliform birds were addressed. Three independent phylogenetic methods confirmed traditional classifications and the existence of the Galloanseres clade.

New candidate species most closely related to penguins

The phylogenetic position of the order Spenisciformes in Aves remains unclear despite several independent analyses based on morphological and molecular data. To address this issue, we determined the complete mtDNA sequence of rockhopper penguins. The mitochondrial genome, excluding the region from the D-loop to 12SrRNA, was also sequenced for petrel, albatross, frigatebird, loon and grebe, which previous studies suggest are related to penguins. A maximum likelihood analysis of the phylogenetic placement of penguins with 23 birds, including 17 species whose mtDNA sequences were previously reported, suggested that storks are the closest extant relatives of penguins, with 78% and 56% bootstrap supports, depending on the choice of outgroup species. Thus, ciconiiform birds constitute new candidates as the closest extant relatives of penguins (previously proposed candidates were either gaviiform, podicipediform, or procellariiform birds). In addition to this new evidence, our analysis gave evidence to some of ambiguous relationships in the avian tree: our analysis supported a basal split between passerines and other neoavians within Neoaves, and rejected the monophyly of Falconiformes as well as that of loons and grebes.

Early Penguin Fossils, Plus Mitochondrial Genomes, Calibrate Avian Evolution

Testing models of macroevolution, and especially the sufficiency of microevolutionary processes, requires good collaboration between molecular biologists and paleontologists. We report such a test for events around the Late Cretaceous by describing the earliest penguin fossils, analyzing complete mitochondrial genomes from an albatross, a petrel, and a loon, and describe the gradual decline of pterosaurs at the same time modern birds radiate. The penguin fossils comprise four naturally associated skeletons from the New Zealand Waipara Greensand, a Paleocene (early Tertiary) formation just above a well-known Cretaceous/Tertiary boundary site. The fossils, in a new genus (Waimanu), provide a lower estimate of 61-62 Ma for the divergence between penguins and other birds and thus establish a reliable calibration point for avian evolution. Combining fossil calibration points, DNA sequences, maximum likelihood, and Bayesian analysis, the penguin calibrations imply a radiation of modern (crown group) birds in the Late Cretaceous. This includes a conservative estimate that modern sea and shorebird lineages diverged at least by the Late Cretaceous about 74 +/- 3 Ma (Campanian). It is clear that modern birds from at least the latest Cretaceous lived at the same time as archaic birds including Hesperornis, Ichthyornis, and the diverse Enantiornithiformes. Pterosaurs, which also coexisted with early crown birds, show notable changes through the Late Cretaceous. There was a decrease in taxonomic diversity, and small- to medium-sized species disappeared well before the end of the Cretaceous. A simple reading of the fossil record might suggest competitive interactions with birds, but much more needs to be understood about pterosaur life histories. Additional fossils and molecular data are still required to help understand the role of biotic interactions in the evolution of Late Cretaceous birds and thus to test that the mechanisms of microevolution are sufficient to explain macroevolution.

Sequences from 14 mitochondrial genes provide a well-supported phylogeny of the Charadriiform birds congruent with the nuclear RAG-1 tree

Because of the difficulties of constructing a robust phylogeny for Charadriiform birds using morphological characters, recent studies have turned to DNA sequences to resolve the systematic uncertainties of family-level relationships in this group. However, trees constructed using nuclear genes or the mitochondrial Cytochrome b gene suggest deep-level relationships of shorebirds that differ from previous studies based on morphology or DNA-DNA hybridization distances. To test phylogenetic hypotheses based on nuclear genes (RAG-1, myoglobin intron-2) and single mitochondrial genes (Cytochrome b), approximately 13,000 bp of mitochondrial sequence was collected for one exemplar species of 17 families of Charadriiformes plus potential outgroups. Maximum likelihood and Bayesian analyses show that trees constructed from long mitochondrial sequences are congruent with the nuclear gene topologies [Chardrii (Lari, Scolopaci)]. Unlike short mitochondrial sequences (such as Cytochrome b alone), longer sequences yield a well-supported phylogeny for shorebirds across various taxonomic levels. Examination of substitution patterns among mitochondrial genes reveals specific genes (especially ND5, ND4, ND2, and COI) that are better suited for phylogenetic analyses among shorebird families because of their relatively homogeneous nucleotide composition among lineages, slower accumulation of substitutions at third codon positions, and phylogenetic utility in both closely and distantly related lineages. For systematic studies of birds in which family and generic levels are examined simultaneously, we recommend the use of both nuclear and mitochondrial sequences as the best strategy to recover relationships that most likely reflect the phylogenetic history of these lineages.

An unusual source of apparent mitochondrial heteroplasmy: duplicate mitochondrial control regions in Thalassarche albatrosses

Molecular ecologists, in search of suitable molecular markers, frequently PCR-amplify regions of mitochondrial DNA from total DNA extracts. This approach, although common, is prone to the co-amplification of nuclear copies of transposed DNA sequences (numts), which can then generate apparent mitochondrial sequence heteroplasmy. In this study we describe the discovery of apparent mitochondrial sequence heteroplasmy in Thalassarche albatrosses but eliminate the possibility of true sequence heteroplasmy and numts and instead reveal the source of the apparent heteroplasmy to be a duplicated control region. The two control regions align easily but are not identical in sequence or in length. Comparisons of functionally significant conserved sequence blocks do not provide evidence of degeneration in either duplicate. Phylogenetic analyses of domain I of both control region copies in five Thalassarche species indicate that they are largely evolving in concert; however, a short section within them is clearly evolving independently. To our knowledge this is the first time contrasting evolutionary patterns have been reported for duplicate control regions. Available evidence suggests that this duplication may be taxonomically widespread, so the results presented here should be considered in future evolutionary studies targeting the control region of all Procellariiformes and potentially other closely related avian groups.

Mitochondrial replication origin stability and propensity of adjacent tRNA genes to form putative replication origins increase developmental stability in Lizards

Secondary structure stability of mitochondrial origins of light-strand replication (OL) presumably reduces delayed formation of light-strand initiating replication forks on the heavy strand. Delayed replication initiation prolongs single strandedness of the heavy strand. More mutations accumulate during the prolonged time spent single stranded. Presumably, delayed replication initiation and excess mutations affect mitochondrial biochemical processes and ultimately morphological outcomes of development at the whole-organism level. This predicts that developmental stability increases with OL secondary structure stability and with formation of OL-like structures by the five tRNA genes flanking recognized OLs. Stable OLs and high percentages of OL-resembling secondary structures of adjacent tRNA genes (predicted by Mfold) correlate positively with developmental stability in three lizard families (Anguidae, Amphisbaenidae, and Polychrotidae). Accounting for effects of the regular OL, Sfold-predicted OL-like propensity of the entire tRNA gene cluster (not of individual genes) correlates with increased developmental stability in Anguidae, also across the entire free-energy range of Boltzmann's distribution of secondary structures. In the fossorial Amphisbaenidae, the OL-like structure-forming propensity of tRNA genes correlates positively with developmental stability for the distribution's sub-optimally stable regions, and negatively for its optimally stable regions, suggesting the thermoregulated functioning of OL vs. flanking tRNA genes as replication origins. Results for polychrotid tRNA genes are intermediate. Anguid tRNA genes possibly function in addition to the regular OL. Mitochondrial tRNA genes may thus frequently acquire and lose the alternative OL function, without sequence (gene) duplication and loss of their primary function.

The rise and fall of the CR1 subfamily in the lineage leading to penguins

The evolution of penguins has been investigated extensively, although inconclusively, by morphologists, biogeographers and molecular phylogeneticists. We investigated this issue using retroposon analysis of insertions of CR1, which is a member of the LINE (long interspersed element) family, in the genomes of penguins and penguin relatives. The retroposon method is a powerful tool for identifying monophyletic groups. Because retroposons often show different relative frequencies of retroposition during evolution, it is first necessary to identify a certain subgroup that was specifically active during the period when the species in question diverged. Hence, we systematically analyzed many CR1 members isolated from penguin and penguin-related genomes. These CR1s are divided into at least three distinct subgroups that share diagnostic nucleotide insertions and/or deletions, namely, penguin CR1 Sph I, Sph II type A and Sph II type B. The analysis of the inserted retroposons by PCR revealed that different CR1 subfamilies or types had amplified at different rates among different periods during penguin evolution. Namely, the penguin CR1 Sph I subfamily had higher rates of retroposition in a common ancestor of all orders examined in this study or at least in a common ancestor of all extant penguins, and the subfamily Sph II type A also had the same tendency. Therefore, these CR1 members can be used to elucidate the phylogenetic relationships of Sphenisciformes (penguins) among different avian orders. In contrast, the penguin CR1 Sph II type B subfamily had higher rates of retroposition just before and after the emergence of the extant genera in Spheniscidae, suggesting that they are useful for elucidating the intra-relationships among extant penguins. This is the first report for the characterization among the members of CR1 family in avian genomes excluding those of chickens. Hence, this work will be a cornerstone for elucidating the phylogenetic relationships in penguin evolution using the retroposon method.

Molecular architecture and rates of DNA substitutions of the mitochondrial control region of cracid birds

The control region (CR) plays an important role in replication and transcription in the mitochondrial genome. Its supposedly high rate of DNA substitution makes it a suitable marker for studies of population and closely related species. Three domains can be identified in CR, each having its own characteristics regarding base composition, pattern of variability and rate of DNA substitution. We sequenced the complete CR for 27 cracids, including all 11 genera to characterize its molecular organization, describe patterns of DNA substitution along the gene, and estimate absolute rates of DNA substitution. Our results show that cracid CR organization and patterns of DNA substitution are typical of other birds. Conserved blocks C and B, fused CSB2/3, and a putative segment for bidirectional replication not usually present in birds were identified in cracids. We also suggest a new delimitation for domains to allow better comparisons among CRs at different taxonomic ranking. Our estimates of absolute rates of DNA substitution show that, in general, CR evolves at a rate slower than that of two protein-coding genes (CR, 0.14%-0.3%; ND2, 0.37%-0.47%; and cytochrome b, 0.29%-0.35% per site per million years within genera). Within CR domains, rates within genera ranged from 0.05% to 0.8% per site per million years.

Convergent evolution of strigiform and caprimulgiform dark-activity is supported by phylogenetic analysis using the arylalkylamine N-acetyltransferase (Aanat) gene

Alternative hypotheses propose the sister order of owls (Strigiformes) to be either day-active raptors (Falconiformes) or dark-active nightjars and allies (Caprimulgiformes). In an effort to identify molecular characters distinguishing between these hypotheses we examined a gene, arylalkylamine N-acetyltransferase (Aanat), potentially associated with the evolution of avian dark-activity. Partial Aanat coding sequences, and two introns, were obtained from the genomic DNA of 16 species: Strigiformes (four species), Falconiformes (four species), Caprimulgiformes (five species), with outgroups: Ciconiiformes (one species), Passeriformes (one species), and Apterygiformes (one species). Phylogenetic trees derived from aligned, evolutionarily conserved Aanat regions did not consistently recover clades corresponding to orders Strigiformes and Falconiformes but did place a caprimulgiform clade more distant from the strigiform and falconiform species than the latter two groups are to each other. This finding was supported by spectral analysis. The taxonomic distribution of seven intronic indels is consistent with the Aanat derived phylogenetic trees and supports conventional family-level groupings within both Strigiformes and Caprimulgiformes. The phylogenetic analyses also indicate that Caprimulgiformes is a polyphyletic grouping. In conclusion the data support, but do not conclusively prove, the proposal that Falconiformes is the sister order to Strigiformes and therefore, that the dark-activity characteristic of Strigiformes and Caprimulgiformes arose by convergent evolution.

Phylogenetic expression of Galalpha1-4Gal on avian glycoproteins: glycan differentiation inscribed in the early history of modern birds

Glycoproteins containing Galalpha1-4Gal (galabiose) had been rarely found in vertebrates, except in a few species of birds and amphibians. We had previously reported that pigeon (Columba livia) egg white and serum glycoproteins are rich in N-glycans with Galalpha1-4Gal at nonreducing termini. To investigate the origin of Galalpha1-4Gal expression in avian evolution, we examined the presence of Galalpha1-4Gal glycoproteins in egg whites from 20 orders, 88 families, 163 genera, and 181 species of birds, as probed by Western blot with Griffonia simplicifolia-I lectin (terminal alpha-Gal/GalNAc-specific) and anti-P(1) mAb (Galalpha1-4Galbeta1-4GlcNAcbeta1-specific). One of the significant observations is the total absence of Galalpha1-4Gal glycoproteins in Struthioniformes (four species), Tinamiformes (three species), Craciformes (two species), Galliformes (14 species), and Anseriformes (10 species), which are phylogenetically separated from other orders at earlier stage of modern bird diversification (100-65 million years ago). The presence or absence of Galalpha1-4Gal glycoproteins in other avian orders varied by the species (104 species positive, and 44 species negative), even though some of them belong to the same order or family. Our results revealed that the expression of Galalpha1-4Gal glycoproteins is not rare among avians, and is correlated with the phylogeny. The expression was most likely differentiated at earlier stage of diversification in modern birds, but some birds might have lost the facility for the expression relatively recently.

Evolution of preproinsulin gene in birds

The coding region of the preproinsulin gene has been cloned and partly sequenced in a variety of marine and terrestrial birds (28 species). All genes showed the "ancestral" structure with a large intron-2. The size of intron-2 changed considerably during the evolution of birds (2.4-4.2kb). The hydrophobicity of signal peptides was conserved. Bird C-peptides were predicted to be 28 amino acids long, but circulating C-peptides would be only 26 amino acids long, with Passer as a possible exception. Bird C-peptides were found to lack the sequences identified in mammals as responsible for peptide bioactivity and the structure of the central part. In contrast, predicted insulin sequences were highly conserved. Only two types of analog were identified: the hypoactive form (GluA8), present only in Anseriformes and the hyperactive form (His A8), present in all other species. Based on 3'-nucleotide sequence analysis (extending into intron-2), birds appeared to be monophyletic. Five groups were clearly identified: Paleognathae, Galliformes, Anseriformes, Passeriformes, and Charadriiformes. Paleognathae were suggested as the basal group, supporting the traditional view of avian evolution. Subsequent branching identified a gallo-anserae group and a group containing all other Neognathae. Surprisingly, Columba livia (Columbiforme order) clustered with Galliformes. With represented species, Procellariiformes and possibly Ciconiiformes, and Pelicaniformes were suggested as paraphyletic, in agreement with conclusions from some studies based on mitochondrial DNA sequences.