Molecular characterization and evolution of a duck mitochondrial genome

Reduced Nitric Oxide Synthase Involvement in Aigamo Duck Basilar Arterial Relaxation

The basilar arterial endothelium mediates blood vessel relaxation partly through the release of nitric oxide (NO). Apoptosis of cerebrovascular endothelial cells is linked to a high mortality rate in chickens infected with the highly pathogenic avian influenza virus, but interestingly, ducks exhibit a greater resistance to this virus. In this study, we examined the responsiveness of duck basilar arteries (BAs) to various vasoactive substances, including 5-hydroxytryptamine (5-HT), histamine (His), angiotensin (Ang) II, noradrenaline (NA), acetylcholine (ACh), and avian bradykinin ornithokinin (OK), aiming to characterize the receptor subtypes involved and the role of endothelial NO in vitro. Our findings suggest that arterial contraction is mediated with 5-HT1 and H1 receptors, while relaxation is induced with β3-adrenergic and M3 receptors. Additionally, OK elicited a biphasic response in duck BAs, and Ang II had no effect. Endothelial NO appears to be crucial in relaxation mediated with M3 and OK receptors but not β3-adrenergic receptors in the duck BA. The reduced endothelial NO involvement in the receptor-mediated relaxation response in duck BAs represents a clear difference from the corresponding response reported in chicken BAs. This physiological difference may explain the differences in lethality between ducks and chickens when vascular endothelial cells are infected with the virus.

Holarctic phylogeographic structure of Eurasian wigeon (Mareca penelope)

The Eurasian wigeon (Mareca penelope) is one of the most numerous migrant species of waterfowl in the Palearctic. Annually, significant part of the world’s wigeon population makes seasonal flights over distances from tens to thousands or more kilometers. According to different estimates based on banding data, five geographic populations of the species were described in the Palearctic. However, distinct borders between the populations have not been identified. At the same time, no phylogeographic studies have been carried out for the complete native range of wigeon so far. In addition to the fundamental importance of such a study, knowledge of the genetic structure of populations is necessary for the development of measures to increase the number of and preserve this valuable game species. The aim of our work was a phylogeographic analysis of the wigeon across its vast native range in the Palearctic including ducks wintering in North America. We examined genetic diversity and differentiation of wigeon populations identified with banding data, phylogenetic relationships of mtDNA haplotypes and demographic history of populations and species as a whole by sequencing a 661 base-pair 5’-fragment of the mitochondrial control region from 195 individual ducks collected throughout the Palearctic and Nearctic. Genetic diversity was high in all studied populations. A reconstruction of haplotypes phylogeny revealed the absence of geographic structure in the data. Nonetheless, analysis of molecular variance (AMOVA) identified two groups of populations: EuropeanSiberian and East Asian. The former included wigeons from Europe, Siberia and the Atlantic coast of North America, and the latter comprised ducks from Russian Far East, Kamchatka Peninsula, Chukotka Autonomous District, the Aleutian Islands, Alaska, and the Pacific coast of North America.

A new insight into the classification of dusky thrush complex: bearings on the phylogenetic relationships within the Turdidae

Dusky thrush complex comprises of two sister species breeding in SC Siberia, which is the member of thrush Turdus from Turdidae. The phylogenetic resolution of Dusky thrush complex remains controversial, and a detailed research is still necessary. In this research, we determined the complete mtDNAs of both species, and estimated phylogenetic trees based on the mtDNA alignment of these and 21 other Turdidae species, to clarify the taxa status of the Dusky thrush complex. The squenced lengths of these three mitochondrial genomes were 16,737, 16,788 and 16,750 bp. The mtDNAs are circular molecules, containing the 37 typical genes, with an identical gene order and arrangement as those of other Turdidae. The ATG and TAA, respectively, are observed the most commonly start and stop codon. Most of the tRNA could be folded into the canonical cloverleaf secondary structure except for tRNA^Ser (AGY) and tRNA^Leu (CUN), which lose 'DHU' arm. The control region presented a higher A + T content than the average value for the whole mitogenome. The phylogenetic trees reconstructed by the concatenated nucleotide sequences of mtDNA genes (Cyt b, ND2 and COI) indicate the Dusky thrush complex cannot be divided into two species, but the relationships between Dusky thrush subspecies still need additional study. This study improves our understanding of mitogenomic structure and evolution of the Dusky thrush complex, which can provide further insights into our understanding of phylogeny and taxonomy in Turdidae.

The mitochondrial genome of Atrijuglans hetaohei Yang (Lepidoptera: Gelechioidea) and related phylogenetic analyses

Complete mitochondrial genome sequences are of great importance for better understanding the genome-level characteristics and phylogenetic relationships among related species. In this study, the complete mitochondrial genome of Atrijuglans hetaohei Yang is sequenced and analyzed, which is 15,379bp in length (GenBank: KT581634) and contains a typical set of 13 protein-coding genes, 22 tRNA genes, two rRNA genes and a non-coding region (control region). Except for cox1 gene that is initiated by CGA codon, all protein-coding genes start with ATN codons and end with the stop codon T, TA or TAA. All tRNAs have a typical clover-leaf secondary structure, except for trnS1, of which the DHU arm could not form a stable stem-loop structure. The secondary structure of rrnL and rrnS consists of 49 helices and 33 helices, respectively. Phylogenetic analyses of the complete mitochondrial genome sequences and of the amino acid sequences for 13 mitochondrial protein-coding genes among related species support the view that A. hetaohei is more closely related to the Gelechioidea than Yponomeutoidea. This result is consistent with a previous classification based on morphology.

The mitochondrial genome of the Saunders's gull Chroicocephalus saundersi (Charadriiformes: Laridae) and a higher phylogeny of shorebirds (Charadriiformes)

The complete mitogenome of Chroicocephalus saundersi was characterized and compared with the 6 published Charadriiformes mitogenomes. The mitogenome of C. saundersi is a closed circular molecule 16,739 bp in size, and contains 37 genes and a control region. The AT and GC skews are positive and negative, respectively, and in agreement with those of the other Charadriiformes mitogenomes. The mitogenome of C. saundersi contains 3 start codons (ATG, GTG, and ATT), 4 stop codons (TAA, TAG, AGG, and AGA), and an incomplete stop codon (T-) in 13 PCGs. A codon usage analysis of all available Charadriiformes mitogenomes showed that the ATG (78%) and TAA (50.5%) were the most common start codon and stop codon, respectively. An unusual start codon, ATT, is commonly found in the ND3s of Charadriiformes mitogenomes, whereas the more common start codons, ATC and ATA, are rarely found. In all the Laridae species, one extra cytosine was inserted at position 174 in ND3. The control region of C. saundersi is 1180-bp long, with a nucleotide composition of 30.2% A, 28.6% T, 27.3% C, and 14.0% G. Variable numbers of tandem repeats (VNTRs) with nine copies of the 10 bp repeat sequence (AACAACAAAC) are found within the CSB domain of the control region. The ML/BI analyses, based on the amino acids of the 13 mitochondrial PCGs, strongly support the monophyly of the order Charadriiformes, with the suborder Lari considered sister to the Scolopaci, which is in turn a sister group to the suborder Charadrii.

Comparative mitochondrial genomics and phylogenetic relationships of the Crossoptilon species (Phasianidae, Galliformes)

BACKGROUND

Phasianidae is a family of Galliformes containing 38 genera and approximately 138 species, which is grouped into two tribes based on their morphological features, the Pheasants and Partridges. Several studies have attempted to reconstruct the phylogenetic relationships of the Phasianidae, but many questions still remain unaddressed, such as the taxonomic status and phylogenetic relationships among Crossoptilon species. The mitochondrial genome (mitogenome) has been extensively used to infer avian genetic diversification with reasonable resolution. Here, we sequenced the entire mitogenomes of three Crossoptilon species (C. harmani, C. mantchuricum and C. crossoptilon) to investigate their evolutionary relationship among Crossoptilon species.

RESULTS

The complete mitogenomes of C. harmani, C. mantchuricum and C. crossoptilon are 16682 bp, 16690 bp and 16680 bp in length, respectively, encoding a standard set of 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes, and a putative control region. C. auritum and C. mantchuricum are more closely related genetically, whereas C. harmani is more closely related to C. crossoptilon. Crossoptilon has a closer relationship with Lophura, and the following phylogenetic relationship was reconstructed: ((Crossoptilon + Lophura) + (Phasianus + Chrysolophus)). The divergence time between the clades C. harmani-C. crossoptilon and C. mantchuricum-C. auritum is consistent with the uplift of the Tibetan Plateau during the Tertiary Pliocene. The Ka/Ks analysis showed that atp8 gene in the Crossoptilon likely experienced a strong selective pressure in adaptation to the plateau environment.

CONCLUSIONS

C. auritum with C. mantchuricum and C. harmani with C. crossoptilon form two pairs of sister groups. The genetic distance between C. harmani and C. crossoptilon is far less than the interspecific distance and is close to the intraspecific distance of Crossoptilon, indicating that C. harmani is much more closely related to C. crossoptilon. Our mito-phylogenomic analysis supports the monophyly of Crossoptilon and its closer relationship with Lophura. The uplift of Tibetan Plateau is suggested to impact the divergence between C. harmani-C. crossoptilon clade and C. mantchuricum-C. auritum clade during the Tertiary Pliocene. Atp8 gene in the Crossoptilon species might have experienced a strong selective pressure for adaptation to the plateau environment.

The complete mitochondrial genomes of sixteen ardeid birds revealing the evolutionary process of the gene rearrangements

BACKGROUND

The animal mitochondrial genome is generally considered to be under selection for both compactness and gene order conservation. As more mitochondrial genomes are sequenced, mitochondrial duplications and gene rearrangements have been frequently identified among diverse animal groups. Although several mechanisms of gene rearrangement have been proposed thus far, more observational evidence from major taxa is needed to validate specific mechanisms. In the current study, the complete mitochondrial DNA of sixteen bird species from the family Ardeidae was sequenced and the evolution of mitochondrial gene rearrangements was investigated. The mitochondrial genomes were then used to review the phylogenies of these ardeid birds.

RESULTS

The complete mitochondrial genome sequences of the sixteen ardeid birds exhibited four distinct mitochondrial gene orders in which two of them, named as "duplicate tRNA(Glu)-CR" and "duplicate tRNAThr-tRNA(Pro) and CR", were newly discovered. These gene rearrangements arose from an evolutionary process consistent with the tandem duplication--random loss model (TDRL). Additionally, duplications in these gene orders were near identical in nucleotide sequences within each individual, suggesting that they evolved in concert. Phylogenetic analyses of the sixteen ardeid species supported the idea that Ardea ibis, Ardea modesta and Ardea intermedia should be classified as genus Ardea, and Ixobrychus flavicollis as genus Ixobrychus, and indicated that within the subfamily Ardeinae, Nycticorax nycticorax is closely related to genus Egretta and that Ardeola bacchus and Butorides striatus are closely related to the genus Ardea.

CONCLUSIONS

The duplicate tRNAThr-CR gene order is found in most ardeid lineages, suggesting this gene order is the ancestral pattern within these birds and persisted in most lineages via concerted evolution. In two independent lineages, when the concerted evolution stopped in some subsections due to the accumulation of numerous substitutions and deletions, the duplicate tRNAThr-CR gene order was transformed into three other gene orders. The phylogenetic trees produced from concatenated rRNA and protein coding genes have high support values in most nodes, indicating that the mitochondrial genome sequences are promising markers for resolving the phylogenetic issues of ardeid birds when more taxa are added.

Construction and identification of a cDNA library for use in the yeast two-hybrid system from duck embryonic fibroblast cells post-infected with duck enteritis virus

To explore and isolate genes related to duck embryonic fibroblast cells (DEFs) post-infected with duck enteritis virus (DEV), a cDNA library was established using SMART (Switching Mechanism At 5' end of the RNA Transcript) technique coupling with a homologous recombination method. The cells were harvested and total RNA was extracted at 48 h post infection. Then the mRNAs were purified and reverse transcribed to first-strand cDNAs using oligo (dT) primers (CDS III). Subsequently, long distance-PCR was performed, the double-stranded cDNAs were purified, and a transformation assay was carried out in that order. Eventually, a high qualitative library was successfully established according to an evaluation on quality. The transformation efficiency was about 2.33 × 10(6) transformants/4.34 μg pGADT7-Rec (>1.0 × 10(6)). The cell density of the library was 1.75 × 10(9) cells/mL (>2 × 10(7) cells/mL). The titer of the primary cDNA library and amplified cDNA library was 6.75 × 10(5) and 2.33 × 10(7) CFU/mL respectively. The numbers for the primary cDNA library and amplified cDNA library were 1.01 × 10(7) and 1.14 × 10(9), respectively, and the recombinant rate was 97.14 %. The sequence results of 27 randomly picked independent clones revealed the insert ranged from 0.323 to 2.017 kb with an average insert size of 0.807 kb. Full-length transcripts of DEV-CHv LORF3, UL26 and UL35 genes were acquired through sequence similarity analysis from the non-redundant nucleic acid or protein database. Five polyA sites were identified in the DEV-CHv genome. Also, a new transcript of 668 bp was found between the IRS gene and US1 gene of the DEV-CHv genome. Thus, we concluded that the constructed cDNA library will be a useful tool in proteomic analysis of interactions between the DEV and host DEFs, and discovery of biomarkers studies on the mechanism of DEV and subsequently exploitation original vaccines and antiviral drugs to prevent or cure diseases.

Determination and analysis of complete mitochondrial genome sequence of mallard (Anas platyrhychos)

The complete mitochondrial genome (mtDNA) of mallard (Anas platyrhychos) was determined by long and accurate polymerase chain reaction and with primer walking sequence method. The entire genome was 16,606 bp in length. Similar to the typical mtDNA of vertebrates, it contained 37 genes (13 protein-coding genes, 2 rRNA genes, and 22 tRNA genes) and a non-coding region (D-loop). The characteristics of the mitochondrial genome were analyzed and discussed in detail.

Complete sequence and gene organization of the mitochondrial genome of scaly-sided merganser (Mergus squamatus) and phylogeny of some Anatidae species

The scaly-sided merganser (Mergus squamatus) is an endangered bird species on the IUCN Red List with the estimated global population of less than 2,500 individuals at present. In the present study, we studied the complete mitochondrial genome (mtDNA) and the phylogenetic of M. squamatus by PCR amplification and GenBank data. The genome was 16,595 bp in length and contained 37 genes (13 protein coding genes, two rRNAs, and 22 tRNAs) and a non-coding control region (D-loop). All protein-coding genes of M. squamatus mtDNA start with a typical ATG codon, except ND1, COI, and COII uses GTG as their initial codon. TAA, T- and TAG as the terminate codon occurred very commonly in the sequence. All tRNA genes can be folded into canonical cloverleaf secondary structure except for tRNA(Ser) (AGY) and tRNA(Leu) (CUN), which lose ''DHU'' arm. The genome sequences had been deposited in GenBank under accession number HQ833701. Based on the concatenated nucleotide sequences of mtDNA genes (Cyt b and D-loop), we reconstructed phylogenetic trees and discussed the phylogenetic relationships among ten Anatidae species. The results are different from the present classification, and we support Lophodytes cucullatus and Mergullus albellus to be members of the genus Mergus.

Phylogenetic relationships and historical biogeography of neotropical parrots (Psittaciformes: Psittacidae: Arini) inferred from mitochondrial and nuclear DNA sequences

Previous hypotheses of phylogenetic relationships among Neotropical parrots were based on limited taxon sampling and lacked support for most internal nodes. In this study we increased the number of taxa (29 species belonging to 25 of the 30 genera) and gene sequences (6388 base pairs of RAG-1, cyt b, NADH2, ATPase 6, ATPase 8, COIII, 12S rDNA, and 16S rDNA) to obtain a stronger molecular phylogenetic hypothesis for this group of birds. Analyses of the combined gene sequences using maximum likelihood and Bayesian methods resulted in a well-supported phylogeny and indicated that amazons and allies are a sister clade to macaws, conures, and relatives, and these two clades are in turn a sister group to parrotlets. Key morphological and behavioral characters used in previous classifications were mapped on the molecular tree and were phylogenetically uninformative. We estimated divergence times of taxa using the molecular tree and Bayesian and penalized likelihood methods that allow for rate variation in DNA substitutions among sites and taxa. Our estimates suggest that the Neotropical parrots shared a common ancestor with Australian parrots 59 Mya (million of years ago; 95% credibility interval (CrI) 66, 51 Mya), well before Australia separated from Antarctica and South America, implying that ancestral parrots were widespread in Gondwanaland. Thus, the divergence of Australian and Neotropical parrots could be attributed to vicariance. The three major clades of Neotropical parrots originated about 50 Mya (95% CrI 57, 41 Mya), coinciding with periods of higher sea level when both Antarctica and South America were fragmented with transcontinental seaways, and likely isolated the ancestors of modern Neotropical parrots in different regions in these continents. The correspondence between major paleoenvironmental changes in South America and the diversification of genera in the clade of amazons and allies between 46 and 16 Mya suggests they diversified exclusively in South America. Conversely, ancestors of parrotlets and of macaws, conures, and allies may have been isolated in Antarctica and/or the southern cone of South America, and only dispersed out of these southern regions when climate cooled and Antarctica became ice-encrusted about 35 Mya. The subsequent radiation of macaws and their allies in South America beginning about 28 Mya (95% CrI 22, 35 Mya) coincides with the uplift of the Andes and the subsequent formation of dry, open grassland habitats that would have facilitated ecological speciation via niche expansion from forested habitats.

Identification of DNA of human origin based on amplification of human-specific mitochondrial cytochrome b region

Species-specific differences in a non-polymorphic region of the mitochondrial cytochrome b gene appear to be large enough to allow human-specific amplification of forensic DNA samples. We therefore developed a PCR-based method using newly designed primers to amplify a 157-bp portion of the human mitochondrial cytochrome b gene. The forward and reverse primers were designed to hybridize to regions of the human mitochondrial cytochrome b gene with sequences differing from those of chimpanzee by 26% (7 bp/27 bp) and 26% (6 bp/23 bp), respectively. Using this primer pair, we successfully amplified DNA extracted from blood samples of 48 healthy adults. All these human samples produced a single band of the expected size on agarose gel electrophoresis, and the sequence of the single band was shown to be identical to that of the target region (157 bp) by sequence analysis. On the other hand, no visible bands were amplified from DNA extracted from blood samples of animals including non-human primates (chimpanzee, gorilla, Japanese monkey, crab-eating monkey) and other species (cow, pig, dog, goat, rat, chicken and tuna). Thus, DNA producing a single band following PCR amplification using this primer pair can be reasonably interpreted as being of human origin. In addition, aged biological specimens comprising bloodstains, hair shafts and bones were successfully identified as being of human origin, illustrating the applicability of the present method to forensic specimens.

Phylogenies of the Frigatebirds (Fregatidae) and Tropicbirds (Phaethonidae), two divergent groups of the traditional order Pelecaniformes, inferred from mitochondrial DNA sequences

The frigatebirds (Fregatidae) and Tropicbirds (Phaethonidae) represent the most morphologically and behaviorally distinct members of the traditional Order Pelecaniformes. Using 1756bp of mitochondrial DNA sequence consisting of the 12S, ATPase-6, ATPase-8, and COI genes obtained from all extant species, we derive a completely resolved phylogeny for both groups. The inferred relationships among these species are robust to the method of phylogenetic estimation used, and all branches are well supported, in spite of the relatively recent radiation within the frigatebirds. The two families are not closely related either to each other, or to any other putative relatives (e.g., pelicans; Pelecanidae).

Structure and evolution of the avian mitochondrial control region

The structural and evolutionary characteristics of the mitochondrial control region were studied by using control region sequences of 68 avian species. The distribution of the variable nucleotide positions within the control region was found to be genus specific and not dependant on the level of divergence, as suggested before. Saturation was shown to occur at the level of divergence of 10% in pairwise comparisons of the control region sequences, as has also been reported for the third codon positions in ND2 and cytochrome b genes of mtDNA. The ratio of control region vs cytochrome b divergence in pairwise comparisons of the sequences was shown to vary from 0.13 to 21.65, indicating that the control region is not always the most variable region of the mtDNA, but also that there are differences in the rate of divergence among the lineages. Only two of the conserved sequence blocks localized earlier for other species, D box and CSB-1, were found to show a considerable amount of sequence conservation across the avian and mammalian sequences. Additionally, a novel avian-specific sequence block was found.

A molecular phylogeny of anseriformes based on mitochondrial DNA analysis

To study the phylogenetic relationships among Anseriformes, sequences for the complete mitochondrial control region (CR) were determined from 45 waterfowl representing 24 genera, i.e., half of the existing genera. To confirm the results based on CR analysis we also analyzed representative species based on two mitochondrial protein-coding genes, cytochrome b (cytb) and NADH dehydrogenase subunit 2 (ND2). These data allowed us to construct a robust phylogeny of the Anseriformes and to compare it with existing phylogenies based on morphological or molecular data. Chauna and Dendrocygna were identified as early offshoots of the Anseriformes. All the remaining taxa fell into two clades that correspond to the two subfamilies Anatinae and Anserinae. Within Anserinae Branta and Anser cluster together, whereas Coscoroba, Cygnus, and Cereopsis form a relatively weak clade with Cygnus diverging first. Five clades are clearly recognizable among Anatinae: (i) the Anatini with Anas and Lophonetta; (ii) the Aythyini with Aythya and Netta; (iii) the Cairinini with Cairina and Aix; (iv) the Mergini with Mergus, Bucephala, Melanitta, Callonetta, Somateria, and Clangula, and (v) the Tadornini with Tadorna, Chloephaga, and Alopochen. The Tadornini diverged early on from the Anatinae; then the Mergini and a large group that comprises the Anatini, Aythyini, Cairinini, and two isolated genera, Chenonetta and Marmaronetta, diverged. The phylogeny obtained with the control region appears more robust than the one obtained with mitochondrial protein-coding genes such as ND2 and cytb. This suggests that the CR is a powerful tool for bird phylogeny, not only at a small scale (i.e., relationships between species) but also at the family level. Whereas morphological analysis effectively resolved the split between Anatinae and Anserinae and the existence of some of the clades, the precise composition of the clades are different when morphological and molecular data are compared.

Sex-biased gene flow in spectacled eiders (Anatidae): inferences from molecular markers with contrasting modes of inheritance

Genetic markers that differ in mode of inheritance and rate of evolution (a sex-linked Z-specific microsatellite locus, five biparentally inherited microsatellite loci, and maternally inherited mitochondrial [mtDNA] sequences) were used to evaluate the degree of spatial genetic structuring at macro- and microgeographic scales, among breeding regions and local nesting populations within each region, respectively, for a migratory sea duck species, the spectacled eider (Somateria fisheri). Disjunct and declining breeding populations coupled with sex-specific differences in seasonal migratory patterns and life history provide a series of hypotheses regarding rates and directionality of gene flow among breeding populations from the Indigirka River Delta, Russia, and the North Slope and Yukon-Kuskokwim Delta, Alaska. The degree of differentiation in mtDNA haplotype frequency among breeding regions and populations within regions was high (phiCT = 0.189, P < 0.01; phiSC = 0.059, P < 0.01, respectively). Eleven of 17 mtDNA haplotypes were restricted to a single breeding region. Genetic differences among regions were considerably lower for nuclear DNA loci (sex-linked: phiST = 0.001, P > 0.05; biparentally inherited microsatellites: mean theta = 0.001, P > 0.05) than was observed for mtDNA. Using models explicitly designed for uniparental and biparentally inherited genes, estimates of spatial divergence based on nuclear and mtDNA data together with elements of the species' breeding ecology were used to estimate effective population size and degree of male and female gene flow. Differences in the magnitude and spatial patterns of gene correlations for maternally inherited and nuclear genes revealed that females exhibit greater natal philopatry than do males. Estimates of generational female and male rates of gene flow among breeding regions differed markedly (3.67 x 10(-4) and 1.28 x 10(-2), respectively). Effective population size for mtDNA was estimated to be at least three times lower than that for biparental genes (30,671 and 101,528, respectively). Large disparities in population sizes among breeding areas greatly reduces the proportion of total genetic variance captured by dispersal, which may accelerate rates of inbreeding (i.e., promote higher coancestries) within populations due to nonrandom pairing of males with females from the same breeding population.

Duplication and concerted evolution of the mitochondrial control region in the parrot genus Amazona

We report a duplication and rearrangement of the mitochondrial genome involving the control region of parrots in the genus Amazona. This rearrangement results in a gene order of cytochrome b/tRNA(Thr)/pND6/pGlu/CR1/tRNA(Pro)/NADH dehydrogenase 6/tRNA(Glu)/CR2/tRNA(Phe)/12s rRNA, where CR1 and CR2 refer to duplicate control regions, and pND6 and pGlu indicate presumed pseudogenes. In contrast to previous reports of duplications involving the control regions of birds, neither copy of the parrot control region shows any indications of degeneration. Rather, both copies contain many of the conserved sequence features typically found in avian control regions, including the goose hairpin, TASs, the F, C, and D boxes, conserved sequence box 1 (CSB1), and an apparent homolog to the mammalian CSB3. We conducted a phylogenetic analysis of homologous portions of the duplicate control regions from 21 individuals representing four species of Amazona (A. ochrocephala, A. autumnalis, A. farinosa, and A. amazonica) and Pionus chalcopterus. This analysis revealed that an individual's two control region copies (i.e., the paralogous copies) were typically more closely related to one another than to corresponding segments of other individuals (i.e., the orthologous copies). The average sequence divergence of the paralogous control region copies within an individual was 1.4%, versus a mean value of 4.1% between control region orthologs representing nearest phylogenetic neighbors. No differences were found between the paralogous copies in either the rate or the pattern in which the two copies accumulated base pair changes. This pattern suggests concerted evolution of the two control regions, perhaps through occasional gene conversion events. We estimated that gene conversion events occurred on average every 34,670 +/- 18,400 years based on pairwise distances between the paralogous control region sequences of each individual. Our results add to the growing body of work indicating that under some circumstances duplicated mitochondrial control regions are retained through evolutionary time rather than degenerating and being lost, presumably due to selection for a small mitochondrial genome.

Mitochondrial genome organization and vertebrate phylogenetics

With the advent of DNA sequencing techniques the organization of the vertebrate mitochondrial genome shows variation between higher taxonomic levels. The most conserved gene order is found in placental mammals, turtles, fishes, some lizards and Xenopus. Birds, other species of lizards, crocodilians, marsupial mammals, snakes, tuatara, lamprey, and some other amphibians and one species of fish have gene orders that are less conserved. The most probable mechanism for new gene rearrangements seems to be tandem duplication and multiple deletion events, always associated with tRNA sequences. Some new rearrangements seem to be typical of monophyletic groups and the use of data from these groups may be useful for answering phylogenetic questions involving vertebrate higher taxonomic levels. Other features such as the secondary structure of tRNA, and the start and stop codons of protein-coding genes may also be useful in comparisons of vertebrate mitochondrial genomes.

A repeat complex in the mitochondrial control region of Adélie penguins from Antarctica

We have determined the nucleotide sequence of the entire mitochondrial control region (CR) of the Adélie penguin (Pygoscelis adeliae) from Antarctica. Like in most other birds, this CR region is flanked by the gene nad6 and transfer (t)RNA trnE(uuc) at the 5' end and the gene rns and trnF(gaa) at the 3' end. Sequence analysis shows that the Adelie penguin CR contains many elements in common with other CRs including the termination associated sequences (TAS), conserved F, E, D, and C boxes, the conserved sequence block (CSB)-1, as well as the putative light and heavy strand promoters sites (LSP-HSP). We report an extraordinarily long avian control region (1758 bp) which can be attributed to the presence, at the 3' peripheral domain, of five 81-bp repeat sequences, each containing a putative LSP-HSP, followed by 30 tetranucleotide microsatellite repeat sequences consisting of (dC-dA-dA-dA)30. The microsatellite and the 81-bp repeat reside in an area known to be transcribed in other species.

Animal mitochondrial genomes

Animal mitochondrial DNA is a small, extrachromosomal genome, typically approximately 16 kb in size. With few exceptions, all animal mitochondrial genomes contain the same 37 genes: two for rRNAs, 13 for proteins and 22 for tRNAs. The products of these genes, along with RNAs and proteins imported from the cytoplasm, endow mitochondria with their own systems for DNA replication, transcription, mRNA processing and translation of proteins. The study of these genomes as they function in mitochondrial systems-'mitochondrial genomics'-serves as a model for genome evolution. Furthermore, the comparison of animal mitochondrial gene arrangements has become a very powerful means for inferring ancient evolutionary relationships, since rearrangements appear to be unique, generally rare events that are unlikely to arise independently in separate evolutionary lineages. Complete mitochondrial gene arrangements have been published for 58 chordate species and 29 non-chordate species, and partial arrangements for hundreds of other taxa. This review compares and summarizes these gene arrangements and points out some of the questions that may be addressed by comparing mitochondrial systems.

Sequences and comparisons of duck mitochondrial DNA control regions

In this study, the nucleotide sequences of the displacement loop (D-loop) and the ribosomal RNA genes of mitochondrial DNA (mtDNA) were determined from a representative each of two genera of ducks, Cainina muschata and Anas platyrhynchos. The duck mtDNA shows a specific gene order at 5' upstream of D-loop (5'ND6-tRNA(Glu)-D-loop3') that is identical to chick mtDNA but is different from that of mammalian or amphibian (5' cytochrome b-tRNA(Thr)-tRNA(Pro)-D-loop3'). Nucleotide diversity is greatest in the D-loop while being most conserved in the 12S rRNA gene, as indicated from a sequence comparison between duck and chick mtDNA. A consensus sequence in the D-loop region, which may play influential roles in the regulation of transcription and replication of mtDNA, was found in both CSB-1 and repeated sequences of birds. Sequences of four tRNA genes in this region are also reported. Among them, tRNA(Glu) shows the greatest sequence divergence when different order of birds are compared.

Continental breakup and the ordinal diversification of birds and mammals

The classical hypothesis for the diversification of birds and mammals proposes that most of the orders diverged rapidly in adaptive radiations after the Cretaceous/Tertiary (K/T) extinction event 65 million years ago. Evidence is provided by the near-absence of fossils representing modern orders before the K/T boundary. However, fossil-based estimates of divergence time are known to be conservative because of sampling biases, and some molecular/time estimates point to earlier divergences among orders. In an attempt to resolve this controversy, we have estimated times of divergence among avian and mammalian orders with a comprehensive set of genes that exhibit a constant rate of substitution. Here we report molecular estimates of divergence times that average about 50-90% earlier than those predicted by the classical hypothesis, and show that the timing of these divergences coincides with the Mesozoic fragmentation of emergent land areas. This suggests that continental breakup may have been an important mechanism in the ordinal diversification of birds and mammals.

The sequence and structure of the meadow grasshopper (Chorthippus parallelus) mitochondrial srRNA, ND2, COI, COII ATPase8 and 9 tRNA genes

The nucleotide sequence of the mitochondrial ND2, COI, COII, ATPase8, srRNA and nine tRNA genes have been sequenced from two individual of the meadow grasshopper Chorthippus parallelus. Comparisons are made to other insects for which the same regions are completely sequenced. Percentage A + T is found to be relatively low in C. parallelus though consistent with that of the other Orthopteran, Locusta migratoria. The relative number of substitutions observed in the different protein-coding genes was analysed between pairs of insect species sharing different levels of relatedness. A clear change in this rate was observed between the within-genus and between-genera comparisons. This change is interpreted in terms of the functional constraints acting on these four different genes. The patterns seem to result from an early saturation of COI and COII genes with synonymous substitutions, and a tolerance of ND2 and ATPase8 function to high levels of amino acid replacements. This analysis highlights a need for further sequence studies and comparisons between taxa of different levels of divergence in order to understand the patterns of mtDNA evolution on which many evolutionary investigations are based.

Variable numbers of simple tandem repeats make birds of the order ciconiiformes heteroplasmic in their mitochondrial genomes

We have analyzed a variable domain of the mitochondrial DNA control region of 18 avian species. Intra-individual length variation was identified and characterized in 15 species. The occurrence of heteroplasmy among species is phylogenetically consistent with a current classification of birds. Polymerase chain reaction amplifications, direct sequencing, and Southern analysis of mitochondrial DNA showed that the heteroplasmy is due to variable numbers of direct repeats in a tandem organization, located in the control region close to the tRNAPhe gene. The tandem repeats consist of short sequence motifs that vary in size from 4 to 32 base pairs between species. Sequence complexity of the repeat motifs was low, with almost exclusively Ts and Gs in the heavy-strand. Extensive variation in the copy number of the repeats was seen both intra-specifically and within individuals. This is the first report of mitochondrial heteroplasmy characterized at the sequence level in birds.