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

Genetic diversity of two populations of the tufted puffin Fratercula cirrhata (Pallas, 1769)

The tufted puffin Fratercula cirrhata (Charadriiformes: Alcidae) is distributed throughout the boreal and low Arctic areas of the North Pacific, from California, USA to Hokkaido, Japan. Few studies have investigated the genetic diversity of this species. Therefore, we analyzed the genetic diversity of two captive populations using nucleotide sequences of two mitochondrial loci (COX1 and D-loop) and one nuclear locus (RHBG). We sequenced these loci for birds from Tokyo Sea Life Park (Kasai Rinkai Suizokuen), originally from Alaska, and birds from Aqua World Oarai, originally from far eastern Russia. We found five COX1 haplotypes and 17 D-loop haplotypes for the mitochondrial loci, and obtained 14 predicted haplotypes for the nuclear RHBG locus. The major haplotypes of all three loci occurred in individuals from both populations. Thus, there were no clear genetic differences between the populations with respect to these three loci. Although the breeding range of the tufted puffin covers the boreal and low Arctic from California to Hokkaido, our results suggest that the species has not genetically diverged within its breeding range.

New view on the organization and evolution of Palaeognathae mitogenomes poses the question on the ancestral gene rearrangement in Aves

BACKGROUND

Bird mitogenomes differ from other vertebrates in gene rearrangement. The most common avian gene order, identified first in Gallus gallus, is considered ancestral for all Aves. However, other rearrangements including a duplicated control region and neighboring genes have been reported in many representatives of avian orders. The repeated regions can be easily overlooked due to inappropriate DNA amplification or genome sequencing. This raises a question about the actual prevalence of mitogenomic duplications and the validity of the current view on the avian mitogenome evolution. In this context, Palaeognathae is especially interesting because is sister to all other living birds, i.e. Neognathae. So far, a unique duplicated region has been found in one palaeognath mitogenome, that of Eudromia elegans.

RESULTS

Therefore, we applied an appropriate PCR strategy to look for omitted duplications in other palaeognaths. The analyses revealed the duplicated control regions with adjacent genes in Crypturellus, Rhea and Struthio as well as ND6 pseudogene in three moas. The copies are very similar and were subjected to concerted evolution. Mapping the presence and absence of duplication onto the Palaeognathae phylogeny indicates that the duplication was an ancestral state for this avian group. This feature was inherited by early diverged lineages and lost two times in others. Comparison of incongruent phylogenetic trees based on mitochondrial and nuclear sequences showed that two variants of mitogenomes could exist in the evolution of palaeognaths. Data collected for other avian mitogenomes revealed that the last common ancestor of all birds and early diverging lineages of Neoaves could also possess the mitogenomic duplication.

CONCLUSIONS

The duplicated control regions with adjacent genes are more common in avian mitochondrial genomes than it was previously thought. These two regions could increase effectiveness of replication and transcription as well as the number of replicating mitogenomes per organelle. In consequence, energy production by mitochondria may be also more efficient. However, further physiological and molecular analyses are necessary to assess the potential selective advantages of the mitogenome duplications.

Mitochondrial DNA intra-individual variation in a bumblebee species: A challenge for evolutionary studies and molecular identification

Mitochondrial DNA (mtDNA) regions have been widely used as molecular markers in evolutionary studies and species identification. However, the presence of heteroplasmy and NUMTs may represent obstacles. Heteroplasmy is a state where an organism has different mitochondrial haplotypes. NUMTs are nuclear pseudogenes originating from mtDNA sequences transferred to nuclear DNA. Evidences of heteroplasmy were already verified in the bumblebee Bombus morio in an earlier study. The present work investigated in more detail the presence of intra-individual haplotypes variation in this species. Heteroplasmy was detected in individuals from all the ten sampled locations, with an average of six heteroplasmic haplotypes per individual. In addition, some of these heteroplasmic haplotypes were shared among individuals from different locations, suggesting the existence of stable heteroplasmy in B. morio. These results demonstrated that heteroplasmy is likely to affect inferences based on mtDNA analysis, especially in phylogenetic, phylogeographic and population genetics studies. In addition, NUMTs were also detected. These sequences showed divergence of 2.7% to 12% in relation to the mitochondrial haplotypes. These levels of divergence could mislead conclusions in evolutionary studies and affect species identification through DNA barcoding.

Resolving Phylogenetic Relationships within Passeriformes Based on Mitochondrial Genes and Inferring the Evolution of Their Mitogenomes in Terms of Duplications

Mitochondrial genes are placed on one molecule, which implies that they should carry consistent phylogenetic information. Following this advantage, we present a well-supported phylogeny based on mitochondrial genomes from almost 300 representatives of Passeriformes, the most numerous and differentiated Aves order. The analyses resolved the phylogenetic position of paraphyletic Basal and Transitional Oscines. Passerida occurred divided into two groups, one containing Paroidea and Sylvioidea, whereas the other, Passeroidea and Muscicapoidea. Analyses of mitogenomes showed four types of rearrangements including a duplicated control region (CR) with adjacent genes. Mapping the presence and absence of duplications onto the phylogenetic tree revealed that the duplication was the ancestral state for passerines and was maintained in early diverged lineages. Next, the duplication could be lost and occurred independently at least four times according to the most parsimonious scenario. In some lineages, two CR copies have been inherited from an ancient duplication and highly diverged, whereas in others, the second copy became similar to the first one due to concerted evolution. The second CR copies accumulated over twice as many substitutions as the first ones. However, the second CRs were not completely eliminated and were retained for a long time, which suggests that both regions can fulfill an important role in mitogenomes. Phylogenetic analyses based on CR sequences subjected to the complex evolution can produce tree topologies inconsistent with real evolutionary relationships between species. Passerines with two CRs showed a higher metabolic rate in relation to their body mass.

Evidence for a duplicated mitochondrial region in Audubon's shearwater based on MinION sequencing

Mitochondrial genetic markers have been extensively used to study the phylogenetics and phylogeography of many birds, including seabirds of the order Procellariiformes. Evidence suggests that part of the mitochondrial genome of Procellariiformes, especially albatrosses, is duplicated, but no DNA fragment covering the entire duplication has been sequenced. We sequenced the complete mitochondrial genome of a non-albatross species of Procellariiformes, Puffinus lherminieri (Audubon's shearwater) using the long-read MinION (ONT) technology. Two mitogenomes were assembled from the same individual, differing by 52 SNPs and in length. The shorter was 19 kb long while the longer was 21 kb, due to the presence of two identical copies of nad6, three tRNA, and two dissimilar copies of the control region (CR). Contrary to albatrosses, cob was not duplicated. We further detected a complex repeated region of undetermined length between the CR and 12S. Long-read sequencing suggests heteroplasmy and a novel arrangement within the duplicated region, indicating a complex evolution of the mitogenome in Procellariiformes.

Complete mitochondrial genome sequences of the northern spotted owl (Strix occidentalis caurina) and the barred owl (Strix varia; Aves: Strigiformes: Strigidae) confirm the presence of a duplicated control region

We report here the successful assembly of the complete mitochondrial genomes of the northern spotted owl (Strix occidentalis caurina) and the barred owl (S. varia). We utilized sequence data from two sequencing methodologies, Illumina paired-end sequence data with insert lengths ranging from approximately 250 nucleotides (nt) to 9,600 nt and read lengths from 100–375 nt and Sanger-derived sequences. We employed multiple assemblers and alignment methods to generate the final assemblies. The circular genomes of S. o. caurina and S. varia are comprised of 19,948 nt and 18,975 nt, respectively. Both code for two rRNAs, twenty-two tRNAs, and thirteen polypeptides. They both have duplicated control region sequences with complex repeat structures. We were not able to assemble the control regions solely using Illumina paired-end sequence data. By fully spanning the control regions, Sanger-derived sequences enabled accurate and complete assembly of these mitochondrial genomes. These are the first complete mitochondrial genome sequences of owls (Aves: Strigiformes) possessing duplicated control regions. We searched the nuclear genome of S. o. caurina for copies of mitochondrial genes and found at least nine separate stretches of nuclear copies of gene sequences originating in the mitochondrial genome (Numts). The Numts ranged from 226–19,522 nt in length and included copies of all mitochondrial genes except tRNA^Pro, ND6, and tRNA^Glu. Strix occidentalis caurina and S. varia exhibited an average of 10.74% (8.68% uncorrected p-distance) divergence across the non-tRNA mitochondrial genes.

Extreme variation in patterns of tandem repeats in mitochondrial control region of yellow-browed tits (Sylviparus modestus, Paridae)

To investigate the evolutionary pattern and origins of tandem repeats in the mitochondrial control region of the yellow-browed tit (Sylviparus modestus), the control region and another four mitochondrial loci from fifteen individuals were analyzed. A 117-bp tandem repeat unit that repeated once, twice or three times in different individuals was found, and a rarely reported arrangement for this tandem repeats region that a 5' imperfect copy at its downstream and a 3' imperfect copy at its upstream was observed. The haplotype network, phylogenetic trees, and ancestral state reconstruction of the combined dataset of five loci suggested multiple origins of the same repeat number. The turnover model via slipped-strand mispairing was introduced to interpret the results, because mispairing occurred so frequently that multiple origins of certain repeat number were observed. Insertion via recombination should be a better explanation for the origin of this tandem repeat unit, considering characteristics of the combined sequence of the 3' and 5' imperfect copy, including identification of its homolog in other passerines and its predicted secondary structure.

Lack of Structural Variation but Extensive Length Polymorphisms and Heteroplasmic Length Variations in the Mitochondrial DNA Control Region of Highly Inbred Crested Ibis, Nipponia nippon

The animal mitochondrial DNA (mtDNA) length polymorphism and heteroplasmy are accepted to be universal. Here we report the lack of structural variation but the presence of length polymorphism as well as heteroplasmy in mtDNA control region of an endangered avian species - the Crested Ibis (Nipponia nippon). The complete control region was directly sequenced while the distribution pattern and inheritance of the length variations were examined using both direct sequencing and genotyping of the PCR fragments from captive birds with pedigrees, wild birds and a historical specimen. Our results demonstrated that there was no structural variation in the control region, however, different numbers of short tandem repeats with an identical motif of CA3CA2CA3 at the 3'-end of the control region determined the length polymorphisms among and heteroplasmy within individual birds. There were one to three predominant fragments in every bird; nevertheless multiple minor fragments coexist in all birds. These extremely high polymorphisms were suggested to have derived from the 'replication slippage' of a perfect microsatellite evolution following the step-wise mutational model. The patterns of heteroplasmy were found to be shifted between generations and among siblings but rather stable between blood and feather samples. This study provides the first evidence of a very extensive mtDNA length polymorphism and heteroplasmy in the highly inbred Crested Ibis which carries an mtDNA genome lack of structural genetic diversity. The analysis of pedigreed samples also sheds light on the transmission of mtDNA length heteroplasmy in birds following the genetic bottleneck theory. Further research focusing on the generation and transmission of particular mtDNA heteroplasmy patterns in single germ line of Crested Ibis is encouraged by this study.

Comparative analysis of mitochondrial control region in polyploid hybrids of red crucian carp (Carassius auratus) x blunt snout bream (Megalobrama amblycephala)

The entire sequences of the mitochondrial (mt)DNA control region (CR) and portions of its flanking genes in the red crucian carp (RC) and blunt snout bream (BSB) as well as their polyploid hybrids (3nRB, 4nRB and 5nRB) were determined and subjected to a comparative analysis. The mtDNA-CRs of these five fish species ranged from 923 to 937 bp in length, they had the same flanking gene arrangement as other vertebrates and the pattern of nucleotide substitution bias was also similar to that in other vertebrates. Our data are consistent with the viewpoint of three domains [extended terminal associated sequence (ETAS domain), central conserved sequence block domain and conserved sequence block (CSB) domain] within the mtDNA-CR of mammals. On the basis our comparative analysis of the mtDNA-CRs of these five fish species, we were able to identify the consensus sequences of functional conserved units, including the ETAS, CSB-F, CSB-D, CSB-E, CSB1, CSB2 and CSB3 and putative promoter. The percentage of variable nucleotide positions (41.98%) in the central domain was lower than those in the ETAS and conserved domain (71.70 and 47.12%, respectively), suggesting that the central domain was the most conserved part of the mtDNA-CR. These results provide useful and important information for the further study of mtDNA-CR structure in fish. The sequence similarities of mtDNA-CR among the 3nRB, 4nRB, 5nRB hybrids and their respective female parents were higher than those among the 3nRB, 4nRB, 5nRB hybrids and their respective male parents, providing the direct evidence of stringent maternal inheritance of mtDNA-CR in the 3nRB, 4nRB and 5nRB hybrids.

Tandem duplication of mitochondrial DNA in the black-faced spoonbill, Platalea minor

Mitochondrial (mt) heteroplasmy in the control region (CR) of the black-faced spoonbill was investigated using LA-PCR. To avoid amplification of transpositioned nuclear genome fragment from mtDNA (numt), PCR product of the almost-complete mitochondrial genome was amplified using primers designed to anneal on the COIII gene. Then nested LA-PCR product was amplified between the cyt b and 12S rRNA genes using the almost-complete mitochondrial genome PCR product as a template. Nucleotide sequencing revealed tandem duplication composed of two units. The first contains cyt b-1, tRNA(Thr)-1, tRNA(Pro)-1, ND6-1, tRNA(Glu)-1 and CR1, and the second consists of cyt b-2, tRNA(Thr)-2, tRNA(Pro)-2, ND6-2, tRNA(Glu)-2 and CR2, followed by tRNA(Phe) and 12S rRNA. The duplicated cyt b-2 sequence coincided with 499 bp at the 3' end of cyt b-1. With the exception of the CR, the other genes in the duplicated sequence were identical to the original corresponding gene. Even though both CR1 and CR2 contain functional blocks, such as a poly-C site, a goose hairpin and a TAS structure in Domain I, the 3' end of CR1 was followed by a 112 bp sequence (non-coding region) that was not found in CR2 or in sequence homology analysis of similar genes. Meanwhile, CR2 ended in a complicated repeat sequence. The 5' franking region in the Domain I (Region A) and the 3' franking region in the Domain I (Region B) of the two CRs evolve in quite different manners: Region A was highly variable between CR1 and CR2 in the same individuals, while Region B was almost identical between them, which indicates concerted evolution.

Halibut mitochondrial genomes contain extensive heteroplasmic tandem repeat arrays involved in DNA recombination

Background

Halibuts are commercially important flatfish species confined to the North Pacific and North Atlantic Oceans. We have determined the complete mitochondrial genome sequences of four specimens each of Atlantic halibut (Hippoglossus hippoglossus), Pacific halibut (Hippoglossus stenolepis) and Greenland halibut (Reinhardtius hippoglossoides), and assessed the nucleotide variability within and between species.

Results

About 100 variable positions were identified within the four specimens in each halibut species, with the control regions as the most variable parts of the genomes (10 times that of the mitochondrial ribosomal DNA). Due to tandem repeat arrays, the control regions have unusually large sizes compared to most vertebrate mtDNAs. The arrays are highly heteroplasmic in size and consist mainly of different variants of a 61-bp motif. Halibut mitochondrial genomes lacking arrays were also detected.

Conclusion

The complexity, distribution, and biological role of the heteroplasmic tandem repeat arrays in halibut mitochondrial control regions are discussed. We conclude that the most plausible explanation for array maintenance includes both the slipped-strand mispairing and DNA recombination mechanisms.

Mechanisms of population differentiation in seabirds

Despite recent advances in population genetic theory and empirical research, the extent of genetic differentiation among natural populations of animals remains difficult to predict. We reviewed studies of geographic variation in mitochondrial DNA in seabirds to test the importance of various factors in generating population genetic and phylogeographic structure. The extent of population genetic and phylogeographic structure varies extensively among species. Species fragmented by land or ice invariably exhibit population genetic structure and most also have phylogeographic structure. However, many populations (26 of 37) display genetic structure in the absence of land, suggesting that other barriers to gene flow exist. In these populations, the extent of genetic structure is best explained by nonbreeding distribution: almost all species with two or more population-specific nonbreeding areas (or seasons) have phylogeographic structure, and all species that are resident at or near breeding colonies year-round have population genetic structure. Geographic distance between colonies and foraging range appeared to have a weak influence on the extent of population genetic structure, but little evidence was found for an effect of colony dispersion or population bottlenecks. In two species (Galapagos petrel, Pterodroma phaeopygia, and Xantus's murrelet, Synthliboramphus hypoleucus), population genetic structure, and even phylogeographic structure, exist in the absence of any recognizable physical or nonphysical barrier, suggesting that other selective or behavioural processes such as philopatry may limit gene flow. Retained ancestral variation may be masking barriers to dispersal in some species, especially at high latitudes. Allopatric speciation undoubtedly occurs in this group, but reproductive isolation also appears to have evolved through founder-induced speciation, and there is strong evidence that parapatric and sympatric speciation occur. While many questions remain unanswered, results of the present review should aid conservation efforts by enabling managers to predict the extent of population differentiation in species that have not yet been studied using molecular markers, and, thus, enable the identification of management units and evolutionary significant units for conservation.

Inheritance and recombination of mitochondrial genomes in plants, fungi and animals

It is generally assumed that mitochondrial genomes are uniparentally transmitted, homoplasmic and nonrecombining. However, these assumptions draw largely from early studies on animal mitochondrial DNA (mtDNA). In this review, we show that plants, animals and fungi are all characterized by episodes of biparental inheritance, recombination among genetically distinct partners, and selfish elements within the mitochondrial genome, but that the extent of these phenomena may vary substantially across taxa. We argue that occasional biparental mitochondrial transmission may allow organisms to achieve the best of both worlds by facilitating mutational clearance but continuing to restrict the spread of selfish genetic elements. We also show that methodological biases and disproportionately allocated study effort are likely to have influenced current estimates of the extent of biparental inheritance, heteroplasmy and recombination in mitochondrial genomes from different taxa. Despite these complications, there do seem to be discernible similarities and differences in transmission dynamics and likelihood of recombination of mtDNA in plant, animal and fungal taxa that should provide an excellent opportunity for comparative investigation of the evolution of mitochondrial genome dynamics.

Origin and evolution of tandem repeats in the mitochondrial DNA control region of shrikes (Lanius spp.)

The origin and evolution of a 128-bp tandem repeat in the mtDNA control region of shrikes (Lanius: Aves) were investigated. The tandem repeat is present in only two species, L. excubitor and L. ludovicianus. In contrast to the variation in repeat number in L. ludovicianus, all individuals of three subspecies of L. excubitor had three repeats. Comparative analysis suggests that a short direct repeat, and a secondary structure including the tandem repeat and a downstream inverted repeat, may be important in the origin of the tandem repeat by slipped-strand mispairing and its subsequent turnover. Homogenization of repeat sequences is most simply explained by expansion and contraction of the repeat array. Surprisingly, mtDNA sequences from L. excubitor were found to be paraphyletic with respect to L. ludovicianus. These results show the utility of a comparative analysis for insights into the evolutionary dynamics of mtDNA tandem repeats.

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.

Phylogeography of shy and white-capped albatrosses inferred from mitochondrial DNA sequences: implications for population history and taxonomy

The evolutionary relationship between shy (Thalassarche cauta) and white-capped (T. steadi) albatrosses was examined using mitochondrial control region sequences. Results were interpreted in the context of a recent and contentious taxonomic revision that recommended full species status for shy and white-capped albatrosses. Low sequence divergence between shy and white-capped albatrosses (1.8%) and between their close relatives, Salvin's and Chatham albatrosses (2.9%), was observed. Much higher sequence divergence was found between the shy/white-capped pair and the Salvin's/Chatham pair (7.0%). Phylogenetic analyses confirmed the separation of the shy/white-capped pair from the Salvin's/Chatham pair but did not provide species-level resolution. Phylogeographic analyses, including a nested clade analysis, FST estimates and an analysis of molecular variance, indicated unambiguous genetic structuring between shy and white-capped albatrosses, thus confirming the demographic isolation of the species, but showed little to no structuring within each species. The geographical distribution of mtDNA haplotypes and other evidence suggest that shy albatrosses arose through range expansion by white-capped albatrosses.

Characterization of the red knot (Calidris canutus) mitochondrial control region

We sequenced the complete mitochondrial control regions of 11 red knots (Calidris canutus). The control region is 1168 bp in length and is flanked by tRNA glutamate (glu) and the gene ND6 at its 5' end and tRNA phenylalanine (phe) and the gene 12S on its 3' end. The sequence possesses conserved sequence blocks F, E, D, C, CSB-1, and the bird similarity box (BSB), as expected for a mitochondrial copy. Flanking tRNA regions show correct secondary structure, and a relative rate test indicated no significant difference between substitution rates in the sequence we obtained versus the known mitochondrial sequence of turnstones (Charadriiformes: Scolopacidae). These characteristics indicate that the sequence is mitochondrial in origin. To confirm this, we sequenced the control region of a single individual using both purified mitochondrial DNA and genomic DNA. The sequences were identical using both methods. The sequence and methods presented in this paper may now serve as a reference for future studies using knot and other avian control regions. Furthermore, the discovery of five variable sites in 11 knots towards the 3' end of the control region, and the variability of this region in contrast to the more conserved central domain in the alignment between knots and other Charadriiformes, highlights the importance of this area as a source of variation for future studies in knots and other birds.

Dynamics in the evolution of sexual traits: losses and gains, radiation and convergence in yellow wagtails (Motacilla flava)

We analyse patterns of genetic diversity and song complexity in the Palaearctic yellow wagtail (Motacilla flava), a highly polytypic species complex. Mitochondrial and nuclear DNA show that the complex is polyphyletic, despite parallel plumage variation in western and eastern clades. In the western clade there is genetic structure among southern subspecies, haplotype diversity decreases with latitude, and northern subspecies show evidence of bottlenecking and rapid expansions, as expected from isolation in glacial refugia followed by postglacial colonization. However, northern subspecies, which have more divergent male plumages, lack genetic structure and sing simpler songs. Loss of song complexity and evolution of plumage in founder populations are consistent with the Kaneshiro model, which posits that variation among species is a consequence of founder-induced shifts in female preference leading to loss of ancestral male sexual traits. Our results suggest possible postglacial founder-effect mechanisms for the morhological diversification of the yellow wagtail complex.

Mitochondrial DNA sequence evolution and phylogeny of the Atlantic Alcidae, including the extinct great auk (Pinguinus impennis)

The Atlantic auk assemblage includes four extant species, razorbill (Alca torda), dovekie (Alle alle), common murre (Uria aalge), and thick-billed murre (U. lomvia), and one recently extinct species, the flightless great auk (Pinguinus impennis). To determine the phylogenetic relationships among the species, a contiguous 4.2-kb region of the mitochondrial genome from the extant species was amplified using PCR. This region included one ribosomal RNA gene, four transfer RNA genes, two protein-coding genes, the control region, and intergenic spacers. Sets of PCR primers for amplifying the same region from great auk were designed from sequences of the extant species. The authenticity of the great auk sequence was ascertained by alternative amplifications, cloning, and separate analyses in an independent laboratory. Phylogenetic analyses of the entire assemblage, made possible by the great auk sequence, fully resolved the phylogenetic relationships and split it into two primary lineages, Uria versus Alle, Alca, and Pinguinus. A sister group relationship was identified between Alca and Pinguinus to the exclusion of ALLE: Phylogenetically, the flightless great auk originated late relative to other divergences within the assemblage. This suggests that three highly divergent species in terms of adaptive specializations, Alca, Alle, and Pinguinus, evolved from a single lineage in the Atlantic Ocean, in a process similar to the initial adaptive radiation of alcids in the Pacific Ocean.

Systematics of Large White-Headed Gulls: Patterns of Mitochondrial DNA Variation in Western European Taxa

Although the large white-headed gull group (genus Larus) has long been a model in speciation studies, the systematic status and evolutionary relationships of many of its taxa remain unresolved. We used mitochondrial DNA control region and cytochrome-b gene sequences in an attempt to resolve some of those uncertainties. In contrast to previously published results based on nuclear markers, mitochondrial DNA was found to be strongly structured among species, indicating that mitochondrial gene flow is very low. Phylogenetic relationships remain largely unresolved, mainly because of the low amount of variation between species. Horizontal transfer of mitochondrial lineages is demonstrated or suspected between most taxa and obscured the reconstruction of the history of the group. The Mediterranean form michahellis was as differentiated from the other western European species as these are from each others, confirming it is neither conspecific with L. fuscus nor with L. argentatus. The forms fuscus and graellsii do not show any significant differences in haplotypes frequencies, arguing for their subspecific status.

Global relationships amongst black-browed and grey-headed albatrosses: analysis of population structure using mitochondrial DNA and microsatellites

The population structure of black-browed (Thalassarche melanophris and T. impavida) and grey-headed (T. chrysostoma) albatrosses was examined using both mitochondrial DNA (mtDNA) and microsatellite analyses. mtDNA sequences from 73 black-browed and 50 grey-headed albatrosses were obtained from five island groups in the Southern Ocean. High levels of sequence divergence were found in both taxa (0.55-7.20% in black-browed albatrosses and 2.10-3.90% in grey-headed albatrosses). Black-browed albatrosses form three distinct groups: Falklands, Diego Ramirez/South Georgia/Kerguelen, and Campbell Island (T. impavida). T. melanophris from Campbell Island contain birds from each of the three groups, indicating high levels of mixture and hybridization. In contrast, grey-headed albatrosses form one globally panmictic population. Microsatellite analyses on a larger number of samples using seven highly variable markers found similar population structure to the mtDNA analyses in both black-browed and grey-headed albatrosses. Differences in population structure between these two very similar and closely related species could be the result of differences in foraging and dispersal patterns. Breeding black-browed albatrosses forage mainly over continental shelves and migrate to similar areas when not breeding. Grey-headed albatrosses forage mainly at frontal systems, travelling widely across oceanic habitats outside the breeding season. Genetic analyses support the current classification of T. impavida as being distinct from T. melanophris, but would also suggest splitting T. melanophris into two groups: Falkland Islands, and Diego Ramirez/South Georgia/Kerguelen.

Genetic diversity and population history of two related seabird species based on mitochondrial DNA control region sequences

Geographical variation in two related seabird species, the razorbill (Alca torda) and common guillemot (Uria aalge), was investigated using sequence analysis of mitochondrial DNA (mtDNA) control regions. We determined the nucleotide sequence of the variable 5' segment of the control region in razorbills and common guillemots from breeding colonies across the Atlantic Ocean. The ecology and life history characteristics of razorbill and common guillemot are in many respects similar. They are both considered highly philopatric and have largely overlapping distributions in temperate and subarctic regions of the North Atlantic, yet the species were found to differ widely in the extent and spatial distribution of mtDNA variation. Moreover, the differences in genetic differentiation and diversity were in the opposite direction to that expected from a consideration of traditional classifications and current population sizes. Indices of genetic diversity were highest in razorbill and varied among colonies, as did genotype frequencies, suggestive of restrictions to gene flow. The distribution of genetic variation suggests that razorbills originated from a refugial population in the south-western Atlantic Ocean through sequential founder events and subsequent expansion in the east and north. In common guillemots, genetic diversity was low and there was a lack of geographical structure, consistent with a recent population bottleneck, expansion and gene flow. We suggest that the reduced level of genetic diversity and differentiation in the common guillemot is caused by an inherent propensity for repeated population bottlenecks and concomitantly unstable population structure related to their specialized feeding ecology.

Dynamic microsatellites in transcribed regions of gastropod mitochondrial 16S rDNA

Length variations of repetitive sequences in different AT-rich loop-coding regions of mitochondrial 16S rDNA in two gastropod species were discovered during intraspecific haplotype surveys. Examination of the discrete length variation of the basic repeat unit in a phylogenetic framework led to the conclusion of a microsatellite-like mutational dynamic. The observations suggest that the presence of a repetitive sequence structure alone is sufficient to trigger this dynamic.

Slow rate of evolution in the mitochondrial control region of gulls (Aves: Laridae)

We sequenced part of the mitochondrial control region and the cytochrome b gene in 72 specimens from 32 gull species (Laridae, Larini) and 2 outgroup representatives (terns: Laridae, Sternini). Our control region segment spanned the conserved central domain II and the usually hypervariable 3' domain III. Apart from some heteroplasmy at the 3' end of the control region, domain III was not more variable than domain II or the cytochrome b gene. Furthermore, variation in the tempo of evolution of domain III was apparent between phyletic species groups. The lack of variation of the gull control region could not be explained by an increase in the proportion of conserved sequences in these birds, and the gull control region showed an organization similar to those of other avian control regions studied to date. A novel invariant direct repeat was identified in domain II of gulls, and in domain III, two to three inverted, sometimes imperfect, repeats are able to form a significantly stable stem-and-loop structure. These putative secondary structures have not been reported before, and a comparison between species groups showed that they are more stable in the group with the more conserved control region. The unusually slow rate of evolution of control region part III of the gulls could thus be partly explained by the existence of secondary structures in domain III of these species.

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.

The colonization history and present-day population structure of the european great tit (Parus major major)

The colonization history and present-day population structure of the European subspecies of the great tit Parus major major were studied using mitochondrial control region sequences. One major haplotype was found in all but one of the eight sampled populations from Spain to northern Finland. The other haplotypes differed from the common one by just a few substitutions; the overall nucleotide diversity was 0.00187 and haplotype diversity 0.8633. No population structuring was detected. The mismatch distribution followed the expected distribution of an expanding population. The estimated time to the most recent common ancestor coincides with the last glacial period. The results suggest that P. m. major survived the last glacial period in a single isolated refuge probably by the Mediterranean Sea. This was followed by rapid colonization of the European continent and population growth. The most recent range expansion northwards is still occurring. Gene flow between the sampled populations is extensive. It is aided by juvenile dispersal, long-distance movements of juvenile flocks and partial migration in the northern parts of the great tit's range.

Precise sequence assignment of replication origin in the control region of chick mitochondrial DNA relative to 5' and 3' D-loop ends, secondary structure, DNA synthesis, and protein binding

The data reported identify for the first time the sequence of an avian mitochondrial heavy-strand replication origin, OH, located only about 12 nucleotides (nt) downstream from the conserved sequence block CSB-1, as well as the sequence of premature synthesis arrest of the 781 (+/-1) nt D-loop strand, only 6-7 nt downstream from a TAS-like (termination-associated) element. Both sites are associated with putative cruciform secondary structures. A major sequence-specific DNA-binding/cleavage site of a potential regulatory protein, the approximately 36-kDa aMDP1 (shown previously to stimulate mtDNA synthesis), is located about 90 nt upstream of OH. Correlated in vivo analysis of avian genome-length mtDNA replication provides missing evidence on the functional equivalence of D-loop origin with nascent initiation, and on the direction, asymmetry and temporal aspects of a full round of replication. The importance of the results to understanding the regulation of linked replication/transcription and the unusual sequence evolution of avian mtDNA is

Polymorphic simple sequence repeat regions in chloroplast genomes: applications to the population genetics of pines

Simple sequence repeats (SSRs), consisting of tandemly repeated multiple copies of mono-, di-, tri-, or tetranucleotide motifs, are ubiquitous in eukaryotic genomes and are frequently used as genetic markers, taking advantage of their length polymorphism. We have examined the polymorphism of such sequences in the chloroplast genomes of plants, by using a PCR-based assay. GenBank searches identified the presence of several (dA)n.(dT)n mononucleotide stretches in chloroplast genomes. A chloroplast (cp) SSR was identified in three pine species (Pinus contorta, Pinus sylvestris, and Pinus thunbergii) 312 bp upstream of the psbA gene. DNA amplification of this repeated region from 11 pine species identified nine length variants. The polymorphic amplified fragments were isolated and the DNA sequence was determined, confirming that the length polymorphism was caused by variation in the length of the repeated region. In the pines, the chloroplast genome is transmitted through pollen and this PCR assay may be used to monitor gene flow in this genus. Analysis of 305 individuals from seven populations of Pinus leucodermis Ant. revealed the presence of four variants with intrapopulational diversities ranging from 0.000 to 0.629 and an average of 0.320. Restriction fragment length polymorphism analysis of cpDNA on the same populations previously failed to detect any variation. Population subdivision based on cpSSR was higher (Gst = 0.22, where Gst is coefficient of gene differentiation) than that revealed in a previous isozyme study (Gst = 0.05). We anticipate that SSR loci within the chloroplast genome should provide a highly informative assay for the analysis of the genetic structure of plant populations.