Mitochondrial DNA length variation and heteroplasmy in populations of white sturgeon (Acipenser transmontanus)

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.

Mitochondrial sequence-based evolutionary analysis of riverine-swamp hybrid buffaloes of India indicates novel maternal differentiation and domestication patterns

In this study, mitochondrial D-loop sequence data on riverine, swamp and hybrid buffaloes from India have been generated and compared with other reported Indian riverine, Chinese and Bangladeshi swamp buffalo populations. Sequence analysis revealed the presence of 132 haplotypes, with a haplotype diversity of 0.9611 ± 0.0045 and a nucleotide diversity of 0.04801 ± 0.00126. For the first time, the existence of riverine-swamp hybrids among the Indian Chilika buffalo population has been recorded, having 49 chromosomes, which was also confirmed by mitochondrial haplotype sharing between Chilika and Indian swamp as well as Chinese swamp buffalo populations in the network analysis. Phylogenetic analysis documents the sharing of reported pre-domestication haplogroups 'SA1', 'SA2', 'SA3' and 'SB1' between the Chilika and swamp buffalo populations of India, China and Bangladesh, an indication of the migration of swamp buffaloes towards Bangladesh and adjoining lower parts of India and north towards Chinese domestication sites. The results have also been supplemented by multidimension scaling, grouping Indian and Chinese swamp buffaloes more closely together with Bangladeshi buffaloes, but into a separate quadrant, whereas Chilika grouped away from other riverine as well as swamp buffaloes. These findings thus confirm the previous reports that the northeast region of India, close to the Indo-China border, is the point of evolution of swamp buffaloes with multiple sites of domestication.

The spatio-temporal features of chicken mitochondrial ND2 gene heteroplasmy and the effects of nutrition factors on this gene

Mitochondrial heterogeneity is the presence of two or more types of mitochondrial (mt)DNA in the same individual/tissue/cell. It is closely related to animal health and disease. ND2 is a protein-coding gene in mtDNA, which participates in mitochondrial respiratory chain and oxidative phosphorylation. In previous studies, we observed that the mt.A5703T and mt.T5727G sites in the ND2 gene were the heteroplasmic variation sites. We used pyrophosphate sequencing technology to examine chicken mt.A5703T and mt.T5727G heteroplasmic sites in the ND2 gene, in different tissues and at different development stages in chickens. We also investigated whether nutritional factors could affect the mt.A5703T and mt.T5727G heteroplasmy. Our results showed that chicken mt.A5703T and mt.T5727G heteroplasmy had clear spatio-temporal specificities, which varied between tissues/development stages. The mtDNA heterogeneity was relatively stable upon nutrition intervention, 30% dietary energy restriction (from 18 to 48 days old) and different types of dietary fats (at 5% concentration, from 1 to 42 days old) did not change the breast muscle heteroplasmy of broilers at the mt.A5703T and mt.T5727G sites. In addition, multiple potential heteroplasmic sites were detected by clone sequencing in the ND2 region, which potentially reflected abundant heteroplasmy in the chicken mitochondrial genome. These results provide an important reference for further research on heteroplasmy in chicken mitochondria.

Genetic diversity and relationship of Dulong chickens using mitochondrial DNA control region

The genetic structure and evolutionary relationship of Dulong chicken with other native Chinese species remained unclear. In this study, the mitochondrial control region was analyzed in total of 343 samples comprising 59 from Dulong chicken and 284 from 8 other Chinese local breeds revealed 51 mutation sites that defined 42 haplotypes. The maximum genetic variation was observed between the Shimian caoke and Pengxian yellow chickens. Phylogenetic analysis revealed that these local chickens mainly scatter in two southwestern clades. Dulong chickens have close relationship with other native chicken. Finding of this study suggests a single matrilineal lineage of indigenous Dulong chickens.

Islands and hybrid zones: combining the knowledge from "Natural Laboratories" to explain phylogeographic patterns of the European brown hare

Background

The aim of the study was to use hybrid populations as well as island populations of the European brown hare (Lepus europaeus) to explore the effect of evolutionary events, such as the post-deglaciation translocations, spontaneous and human-mediated, local adaptation and the genetic drift in the shaping of the phylogeographic patterns of the species. For this purpose, we used molecular markers, both nuclear and mitochondrial, that are indicative for local adaptation as well as neutral markers to elucidate the patterns of population differentiation based on geographic isolation and the clade of origin. To broaden our analysis, we included data from our previous studies concerning mainland populations, to explore the genetic differentiation in the base of the geographic origin (mainland/island) of the populations.

Results

Our results suggest that local adaptation shapes the differentiation in both genomes, favoring specific alleles in nuclear genes (e.g. DQA) or haplotypes in mtDNA (e.g. Control Region, CR). mtDNA variation was found to be in a higher level and was able to give a phylogeographic signal for the populations. Furthermore, the degree of variation was influenced not only by the geographic origin, but also by the clade of origin, since specific island populations of Anatolian origin showed a greater degree of variation compared to specific mainland populations of the European clade. Concerning the hybrid population, we confirmed the existence of both clades in the territory and we provided a possible explanation for the lack of introgression between the clades.

Conclusion

Our results indicate that the Quaternary’s climatic oscillations played a major role in the shaping of the phylogeographic patterns of the species, by isolating populations in the distinct refugia, where they adapted and differentiate in allopatry, leading to genome incompatibilities observed nowadays.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1354-y) contains supplementary material, which is available to authorized users.

Genetic diversity of Jiangsu native chicken breeds assessed with the mitochondrial DNA D-loop region

1. The objective of this study was to determine the origin and evolution of chickens from 5 native breeds that are traditionally raised in Jiangsu Province. 2. To address this question, the complete mitochondrial DNA D-loop sequence of 149 chickens from 5 native breeds of Jiangsu Province was analysed. 3. Sequence read lengths of the native breeds were 1231 to 1232 bp, with a single-base deletion from the 859 bp site in the 1231 bp haplotype. A total of 33 variable sites that defined 19 haplotypes were identified. The average haplotype diversity and nucleotide diversity were 0.862 ± 0.017 and 0.00591 ± 0.00135. 4. Phylogenetic analysis showed that genetic structure of the mtDNA haplotypes of Jiangsu chickens are distributed across 5 clades (haplogroups): Clades A, B, C, D, and E. However, most of the individuals characterised in this study belonged to clades A and B. 5. The results of this study indicate that Jiangsu chicken populations have relatively low nucleotide and haplotype diversity and likely share 5 common maternal lineages.

Genetic Diversity and Phylogenetic Evolution of Tibetan Sheep Based on mtDNA D-Loop Sequences

The molecular and population genetic evidence of the phylogenetic status of the Tibetan sheep (Ovis aries) is not well understood, and little is known about this species' genetic diversity. This knowledge gap is partly due to the difficulty of sample collection. This is the first work to address this question. Here, the genetic diversity and phylogenetic relationship of 636 individual Tibetan sheep from fifteen populations were assessed using 642 complete sequences of the mitochondrial DNA D-loop. Samples were collected from the Qinghai-Tibetan Plateau area in China, and reference data were obtained from the six reference breed sequences available in GenBank. The length of the sequences varied considerably, between 1031 and 1259 bp. The haplotype diversity and nucleotide diversity were 0.992±0.010 and 0.019±0.001, respectively. The average number of nucleotide differences was 19.635. The mean nucleotide composition of the 350 haplotypes was 32.961% A, 29.708% T, 22.892% C, 14.439% G, 62.669% A+T, and 37.331% G+C. Phylogenetic analysis showed that all four previously defined haplogroups (A, B, C, and D) were found in the 636 individuals of the fifteen Tibetan sheep populations but that only the D haplogroup was found in Linzhou sheep. Further, the clustering analysis divided the fifteen Tibetan sheep populations into at least two clusters. The estimation of the demographic parameters from the mismatch analyses showed that haplogroups A, B, and C had at least one demographic expansion in Tibetan sheep. These results contribute to the knowledge of Tibetan sheep populations and will help inform future conservation programs about the Tibetan sheep native to the Qinghai-Tibetan Plateau.

Mitochondrial DNA variability to explore the relationship complexity of Schizothoracine (Teleostei: Cyprinidae)

Despite numerous studies on the taxonomy of a highly complex group of schizothoracine (snow trouts), with over five recognized species from Kashmir, India (Schizothorax niger, Schizothorax esocinus, Schizothorax plagiostomus, Schizothorax curvifrons and Schizothorax labiatus) based on traditional morphological data, the relationships between these species is poorly understood and the taxonomic validity is still under debate. To resolve the evolutionary relationships among these species, we sequenced mitochondrial fragments, including 16Sr RNA, Cytb and the D-loop. Separate analyses of 16S and Cytb showed intermixing of the species and 16S was found more conserved than Cytb. The D-loop was found highly variable and showed length variation between and within species. Length variation was observed in di-nucleotide (TA)n microsatellite repeats with a variable number of repeat units (n = 7-14) that did not show heteroplasmy. Central conserved sequence blocks (CSBs) in D-loop sequences were found comparable to other vertebrate species. All phylogenetic reconstructions recovered the focal taxa as a monophyletic clade within the schizothoracines. Analyses with combined mitochondrial data sets showed close genetic relationships of all the five species. In addition to a close relationship between S. niger and S. curvifrons, two distinct groupings of S. ecoscinus and S. plagiostomus were supported by all the analyses. This study gives an insight into molecular phylogeny of the species and improves our understanding of historical and taxonomic relationships derived from morphological and ecological studies.

Length variation and sequence divergence in mitochondrial control region of Schizothoracine (Teleostei: Cyperinidae) species

Schizothoracine fish commonly called snow trouts inhibit the entire network of snow and spring fed cool waters of Kashmir, India. Over 10 species reported earlier, only five species have been found, these include Schizothorax niger, Schizothorax esocinus, Schizothorax plagiostomus, Schizothorax curvifrons and Schizothorax labiatus. The relationship between these species is contradicting. To understand the evolutionary relation of these species, we examined the sequence information of mitochondrial D-loop of 25 individuals representing five species. Sequence alignment showed D-loop region highly variable and length variation was observed in di-nucleotide (TA)n microsatellite between and within species. Interestingly, all these species have (TA)n microsatellite not associated with longer tandem repeats at the 3' end of the mitochondrial control region and do not show heteroplasmy. Our analysis also indicates the presence of four conserved sequence blocks (CSB), CSB-D, CSB-1, CSB-II and CSB-III, four (Termination Associated Sequence) TAS motifs and 15bp pyrimidine block within the mitochondrial control region, that are highly conserved within genus Schizothorax when compared with other species. The phylogenetic analysis carried by Maximum likelihood (ML), Neighbor Joining (NJ) and Bayesian inference (BI) generated almost identical results. The resultant BI tree showed a close genetic relationship of all the five species and supports two distinct grouping of S. esocinus species. Besides the species relation, the presence of length variation in tandem repeats is attributed to differences in predicting the stability of secondary structures. The role of CSBs and TASs, reported so far as main regulatory signals, would explain the conservation of these elements in evolution.

Geographic variation within a tandemly repeated mitochondrial DNA D-loop region of a North American freshwater fish, Pylodictis olivaris

The present study reports the distribution of a 35-bp mitochondrial DNA (mtDNA) D-loop tandemly repeated sequence in the populations of a North American freshwater catfish, Pylodictis olivaris, and the important role of a past geological event in the phylogeographic pattern of this species. A total of 330 individuals of flathead catfish, representing 34 drainages throughout the species' native range in the United States, were collected. While more than 70% of individuals sampled from the Southeastern Gulf Coast drainages were characterized by the presence of a 35-bp mtDNA D-loop tandem repeat proximal to the 5' end, more than 95% of samples from the Mississippi River and its tributaries, as well as from the drainages of the Southwest Gulf Coast region, lack this tandem repeat. Concomitantly, phylogenetic analyses revealed the existence of two distinct matrilineal lineages (lineage I and II) of P. olivaris, which were estimated to have diverged from a common ancestor sometime between 0.70 and 2.05myr ago. While one lineage is comprised of samples from the Mississippi River and its tributaries and rivers draining to the Southwest Gulf Coast, the other lineage is comprised of samples from the Southeastern Gulf Coast drainages. Each lineage also has two sub-lineages, which also showed geographic specificity.

An examination of the origin and evolution of additional tandem repeats in the mitochondrial DNA control region of Japanese sika deer (Cervus Nippon)

Tandem repeat units are only detected in the left domain of the mitochondrial DNA control region in sika deer. Previous studies showed that Japanese sika deer have more tandem repeat units than its cousins from the Asian continent and Taiwan, which often have only three repeat units. To determine the origin and evolution of these additional repeat units in Japanese sika deer, we obtained the sequence of repeat units from an expanded dataset of the control region from all sika deer lineages. The functional constraint is inferred to act on the first repeat unit because this repeat has the least sequence divergence in comparison to the other units. Based on slipped-strand mispairing mechanisms, the illegitimate elongation model could account for the addition or deletion of these additional repeat units in the Japanese sika deer population. We also report that these additional repeat units could be occurring in the internal positions of tandem repeat regions, possibly via coupling with a homogenization mechanism within and among these lineages. Moreover, the increased number of repeat units in the Japanese sika deer population could reflect a balance between mutation and selection, as well as genetic drift.

Genetic variation and phylogenetic relationship between three serow species of the genus Capricornis based on the complete mitochondrial DNA control region sequences

The molecular evidence of phylogenetic status regarding the Formosan serow (Capricornis swinhoei) is not robust and little is known about the genetic diversity of the Sumatran serow (Capricornis sumatraensis), which partly is due to the hardness in sample collection. Here we determined the sequences of the complete mitochondrial DNA control region (1,014 bp) of 19 Sumatran-serow individuals. Nine new haplotypes were defined based on 78 variable sites. Combined analysis with other 32 haplotypes downloaded from the public database, including 1 Sumatran-serow, 11 Formosan-serow and 20 Japanese-serow (Capricornis crispus) haplotypes, a relatively high level of nucleotide diversity was first observed in Sumatran serow (π = 0.0249). By comparative analysis with structural consensus sequences from other mammals, we have identified central, left and right domains and depicted the putative functional structure, including extend termination associated sequences and conserve sequence blocks, in mtDNA control region. The alignment of mtDNA control region revealed that both Sumatran and Japanese serow have two tandem repeats (TRs), but three TRs in Formosan serow. Phylogenetic analyses revealed that the Formosan serow is distinct species with the Japanese serow, but a sister group with the Sumatran serow. The divergence time estimated among three serow species revealed that Pleistocene climate changes and the uplift of Qinghai-Tibetan plateau might play an important role in the genetic differentiation of the serows. These results mainly provide the convinced evidence on the genetic relationship between three serow species.

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.

Genetic diversity and demographical history of Coilia ectenes (Clupeiformes: Engraulidae) inferred from the complete control region sequences of mitochondrial DNA

Coilia ectenes is a commercially important fishery species. In this study, genetic diversity and population structure of C. ectenes were examined by using mitochondrial DNA control region sequences in 246 individuals sampled from 10 localities in China. One hundred and ninety-five polymorphic sites defined 184 distinct haplotypes, revealing a moderately high haplotype diversity (Hd) and a relatively low nucleotide diversity (π) in the 10 localities. An excess of unique haplotypes at most sample locations were detected, which might influence the genetic structure of the C. ectenes populations. Hd ranged from 0.939 to 1.000 and π ranged from 0.26% to 1.15%. The Dongting fish population had the highest π level. The genetic distances ranged from 0.26% to 1.03% within populations and from 0.56% to 4.90% between populations. The distances between the Fuzhou (FZ) population and other populations were mostly >4.8%. Neighbor-joining tree indicated distinct patterns of phylogeographic structure among haplotypes from FZ population and those from other populations. Analyses of molecular variance and F(st) statistics suggested that the divergence existed among populations from 10 localities, indicating that gene flow might be restricted among those regions, despite the wide dispersal. In addition, neutral tests and analysis of mismatch distribution suggested that C. ectenes might have undergone a population expansion. Our study revealed the extant population genetic diversity and structure of the C. ectenes, and was in favor of the related fishery management issues including fishery stock identification and conservation.

Population genetics of the cytoplasm and the units of selection on mitochondrial DNA in Drosophila melanogaster

Biological variation exists across a nested set of hierarchical levels from nucleotides within genes to populations within species to lineages within the tree of life. How selection acts across this hierarchy is a long-standing question in evolutionary biology. Recent studies have suggested that genome size is influenced largely by the balance of selection, mutation and drift in lineages with different population sizes. Here we use population cage and maternal transmission experiments to identify the relative strength of selection at an individual and cytoplasmic level. No significant trends were observed in the frequency of large (L) and small (S) mtDNAs across 14 generations in population cages. In all replicate cages, new length variants were observed in heteroplasmic states indicating that spontaneous length mutations occurred in these experimental populations. Heteroplasmic flies carrying L genomes were more frequent than those carrying S genomes suggesting an asymmetric mutation dynamic from larger to smaller mtDNAs. Mother-offspring transmission of heteroplasmy showed that the L mtDNA increased in frequency within flies both between and within generations despite sampling drift of the same intensity as occurred in population cages. These results suggest that selection for mtDNA size is stronger at the cytoplasmic than at the organismal level. The fixation of novel mtDNAs within and between species requires a transient intracellular heteroplasmic stage. The balance of population genetic forces at the cytoplasmic and individual levels governs the units of selection on mtDNA, and has implications for evolutionary inference as well as for the effects of mtDNA mutations on fitness, disease and aging.

Inheritance of microsatellite loci in the white sturgeon (Acipenser transmontanus)

Nine tetramer motif (GATA)n microsatellite systems were developed for use in the white sturgeon, Acipenser transmontanus. We report inheritance patterns for these nine systems, which range from one possible disomic system to tetrasomy and octosomy, with some systems containing null alleles. Because of the complex modes of inheritance underlying these systems and the highly duplicated nature of the genome, we propose each allele be scored as its own dominant marker, similar to AFLPs or RAPDs. The utility of this method is validated by the observation that individual alleles within a microsatellite system generally fit the expectation for independent transmission and fit the expected transmission frequency for single copy nuclear markers.

Population genetics of shortnose sturgeon Acipenser brevirostrum based on mitochondrial DNA control region sequences

Shortnose sturgeon is an anadromous North American acipenserid that since 1973 has been designated as federally endangered in US waters. Historically, shortnose sturgeon occurred in as many as 19 rivers from the St. John River, NB, to the St. Johns River, FL, and these populations ranged in census size from 10(1) to 10(4), but little is known of their population structure or levels of gene flow. We used the polymerase chain reaction (PCR) and direct sequence analysis of a 440 bp portion of the mitochondrial DNA (mtDNA) control region to address these issues and to compare haplotype diversity with population size. Twenty-nine mtDNA nucleotide-substitution haplotypes were revealed among 275 specimens from 11 rivers and estuaries. Additionally, mtDNA length variation (6 haplotypes) and heteroplasmy (2-5 haplotypes for some individuals) were found. Significant genetic differentiation (P < 0.05) of mtDNA nucleotide-substitution haplotypes and length-variant haplotypes was observed among populations from all rivers and estuaries surveyed with the exception of the Delaware River and Chesapeake Bay collections. Significant haplotype differentiation was even observed between samples from two rivers (Kennebec and Androscoggin) within the Kennebec River drainage. The absence of haplotype frequency differences between samples from the Delaware River and Chesapeake Bay reflects a probable current absence of spawning within the Chesapeake Bay system and immigration of fish from the adjoining Delaware River. Haplotypic diversity indices ranged between 0.817 and 0.641; no relationship (P > 0.05) was found between haplotype diversity and census size. Gene flow estimates among populations were often low (< 2.0), but were generally higher at the latitudinal extremes of their distribution. A moderate level of haplotype diversity and a high percentage (37.9%) of haplotypes unique to the northern, once-glaciated region suggests that northern populations survived the Pleistocene in a northern refugium. Analysis of molecular variance best supported a five-region hierarchical grouping of populations, but our results indicate that in almost all cases populations of shortnose sturgeon should be managed as separate units.

Sturgeon genetics and cytogenetics: recent advancements and perspectives

The aim of this review is to introduce current knowledge in the field of sturgeon genetics. The first section deals with sturgeon cytogenetics, reviewing karyotype organization and polyploidization events during evolution of Acipenseriformes. The second section concerns the results of applications of molecular biology to studies of phylogenetic relationships between extant species, intraspecific analysis of wild populations and stocks for conservation purposes, together with characterization of molecular markers for species identification, relevant to forensic and conservation issues.

Genome duplication events and functional reduction of ploidy levels in sturgeon (Acipenser, Huso and Scaphirhynchus)

Sturgeon (order Acipenserformes) provide an ideal taxonomic context for examination of genome duplication events. Multiple levels of ploidy exist among these fish. In a novel microsatellite approach, data from 962 fish from 20 sturgeon species were used for analysis of ploidy in sturgeon. Allele numbers in a sample of individuals were assessed at six microsatellite loci. Species with approximately 120 chromosomes are classified as functional diploid species, species with approximately 250 chromosomes as functional tetraploid species, and with approximately 500 chromosomes as functional octaploids. A molecular phylogeny of the sturgeon was determined on the basis of sequences of the entire mitochondrial cytochrome b gene. By mapping the estimated levels of ploidy on this proposed phylogeny we demonstrate that (I) polyploidization events independently occurred in the acipenseriform radiation; (II) the process of functional genome reduction is nearly finished in species with approximately 120 chromosomes and more active in species with approximately 250 chromosomes and approximately 500 chromosomes; and (III) species with approximately 250 and approximately 500 chromosomes arose more recently than those with approximately 120 chromosomes. These results suggest that gene silencing, chromosomal rearrangements, and transposition events played an important role in the acipenseriform genome formation. Furthermore, this phylogeny is broadly consistent with previous hypotheses but reveals a highly supported oceanic (Atlantic-Pacific) subdivision within the Acipenser/Huso complex.

Restriction site mapping of the mitochondrial DNA of the four-wing flyingfish, Hirundichthys affinis

Mitochondrial DNA was isolated from samples of the four-wing flyingfish, Hirundichthys affinis, collected in Barbados in January 1996 and subjected to restriction enzyme analysis, using 13 restriction endonucleases which recognise hexanucleotide sequences, in single and double digests. The resulting restricted DNA fragments were used to map the 14 enzyme recognition sites of 6 endonucleases (7 had no sites) on the flyingfish mtDNA molecule for the first time. In addition, the mtDNA D-loop region was positioned on the restriction site map, for the first time, by selective restriction digestion of the mtDNA molecule followed by polymerase chain reaction (PCR) amplification of the resulting fragments using specific mtDNA D-loop primers. The size of the flyingfish mtDNA molecule (18 kb) was also determined.

Heteroplasmy in the mtDNA control region of sturgeon (Acipenser, Huso and Scaphirhynchus)

Data from 1238 fishes from 19 sturgeon species and 1 paddlefish were used to analyze heteroplasmy in sturgeon. Lengths of central repeat units ranged from 74 to 83 bp among sturgeon species. No repeat sequence was found in the paddlefish, Polyodon spathula. A general feature of the repeat units was the presence of termination associated sequence (TAS) motifs. About 50% of 138 interspecific mutations observed among the D-loop sequences are located 10 bp down- and upstream from these TAS motifs. Interestingly, most homoplasmic species showed deletions upstream to the TAS motifs, whereas deletions downstream to the TAS motifs observed in two species do not seem to preclude heteroplasmy. Calculations of secondary structures and thermal stabilities of repeat units showed DeltaG values for all heteroplasmic species to be <-8 and for most homoplasmic species DeltaG value to be >-8. Most heteroplasmic fishes had two and/or three repeat units. No homoplasmic sturgeon with >2 repeat units were observed. Molecular phylogeny based on the entire cytochrome b showed that heteroplasmy probably resulted from a single evolutionary event. Our data demonstrate that heteroplasmy is present in most sturgeon species and suggest that the thermal stability of the secondary structure of the repeat unit in combination with mutations downstream of the TAS sequences influences heteroplasmy.

Tandemly repeated sequences in mtDNA control region of whitefish, Coregonus lavaretus

Length variation of the mitochondrial DNA control region was observed with PCR amplification of a sample of 138 whitefish (Coregonus lavaretus). Nucleotide sequences of representative PCR products showed that the variation was due to the presence of an approximately 100-bp motif tandemly repeated two, three, or five times in the region between the conserved sequence block-3 (CSB-3) and the gene for phenylalanine tRNA. This is the first report on the tandem array composed of long repeat units in mitochondrial DNA of salmonids.

Tandemly repeated sequences in the mitochondrial DNA control region and phylogeography of the Pike-Perches Stizostedion

DNA sequences from the mitochondrial DNA control region are used to test the phylogeographic relationships among the pike-perches, Stizostedion (Teleostei: Percidae) and to examine patterns of variation. Sequences reveal two types of variability: single nucleotide polymorphisms and 6 to 14 copies of 10- to 11-base-pair tandemly repeated sequences. Numbers of copies of the tandem repeats are found to evolve too rapidly to detect phylogenetic signal at any taxonomic level, even among populations. Sequence similarities of the tandem repeats among Stizostedion and other percids suggest concerted evolutionary processes. Predicted folding of the tandem repeats and their proximity to termination-associated sequences indicate that secondary structure mediates slipped-strand mispairing among the d-loop, heavy, and light strands. Neighbor-joining and maximum parsimony analyses of sequences indicate that the genus is divided into clades on the continents of North America and Eurasia. Calibrating genetic distances with divergence times supports the hypothesis that Stizostedion dispersed from Eurasia to North America across a North Pacific Beringial land bridge approximately 4 million years before present, near the beginning of the Pliocene Epoch. The North American S. vitreum and S. canadense appear separated by about 2.75 million years, and the Eurasian S. lucioperca and S. volgensis are diverged by about 1.8 million years, suggesting that speciation occurred during the late Pliocene Epoch.

Heteroplasmy, length and sequence variation in the mtDNA control regions of three percid fish species (Perca fluviatilis, Acerina cernua, Stizostedion lucioperca)

The nucleotide sequence of the control region and flanking tRNA genes of perch (Perca fluviatilis) mtDNA was determined. The organization of this region is similar to that of other vertebrates. A tandem array of 10-bp repeats, associated with length variation and heteroplasmy was observed in the 5' end. While the location of the array corresponds to that reported in other species, the length of the repeated unit is shorter than previously observed for tandem repeats in this region. The repeated sequence was highly similar to the Mt5 element which has been shown to specifically bind a putative D-loop DNA termination protein. Of 149 perch analyzed, 74% showed length variation heteroplasmy. Single-cell PCR on oocytes suggested that the high level of heteroplasmy is passively maintained by maternal transmission. The array was also observed in the two other percid species, ruffe (Acerina cernua) and zander (Stizostedion lucioperca). The array and the associated length variation heteroplasmy are therefore likely to be general features of percid mtDNAs. Among the perch repeats, the mutation pattern is consistent with unidirectional slippage, and statistical analyses supported the notion that the various haplotypes are associated with different levels of heteroplasmy. The variation in array length among and within species is ascribed to differences in predicted stability of secondary structures made between repeat units.

Evolution of repeated sequence arrays in the D-loop region of bat mitochondrial DNA

Analysis of mitochondrial DNA control region sequences from 41 species of bats representing 11 families revealed that repeated sequence arrays near the tRNA-Pro gene are present in all vespertilionine bats. Across 18 species tandem repeats varied in size from 78 to 85 bp and contained two to nine repeats. Heteroplasmy ranged from 15% to 63%. Fewer repeats among heteroplasmic than homoplasmic individuals in a species with up to nine repeats indicates selection may act against long arrays. A lower limit of two repeats and more repeats among heteroplasmic than homoplasmic individuals in two species with few repeats suggests length mutations are biased. Significant regressions of heteroplasmy, theta and pi, on repeat number further suggest that repeat duplication rate increases with repeat number. Comparison of vespertilionine bat consensus repeats to mammal control region sequences revealed that tandem repeats of similar size, sequence and number also occur in shrews, cats and bighorn sheep. The presence of two conserved protein-binding sequences in all repeat units indicates that convergent evolution has occurred by duplication of functional units. We speculate that D-loop region tandem repeats may provide signal redundancy and a primitive repair mechanism in the event of somatic mutations to these binding sites.

Evolution of repeated sequence arrays in the D-loop region of bat mitochondrial DNA

Analysis of mitochondrial DNA control region sequences from 41 species of bats representing 11 families revealed that repeated sequence arrays near the tRNA-Pro gene are present in all vespertilionine bats. Across 18 species tandem repeats varied in size from 78 to 85 bp and contained two to nine repeats. Heteroplasmy ranged from 15% to 63%. Fewer repeats among heteroplasmic than homoplasmic individuals in a species with up to nine repeats indicates selection may act against long arrays. A lower limit of two repeats and more repeats among heteroplasmic than homoplasmic individuals in two species with few repeats suggests length mutations are biased. Significant regressions of heteroplasmy, theta and pi, on repeat number further suggest that repeat duplication rate increases with repeat number. Comparison of vespertilionine bat consensus repeats to mammal control region sequences revealed that tandem repeats of similar size, sequence and number also occur in shrews, cats and bighorn sheep. The presence of two conserved protein-binding sequences in all repeat units indicates that convergent evolution has occurred by duplication of functional units. We speculate that D-loop region tandem repeats may provide signal redundancy and a primitive repair mechanism in the event of somatic mutations to these binding sites.

Stock structure and homing fidelity in Gulf of Mexico sturgeon (Acipenser oxyrinchus desotoi) based on restriction fragment length polymorphism and sequence analyses of mitochondrial DNA

Efforts have been proposed worldwide to restore sturgeon populations through the use of hatcheries to supplement natural reproduction and to reintroduce sturgeon where they have become extinct. We examined the population structure and inferred the extent of homing in the anadromous Gulf of Mexico (Gulf) sturgeon (Acipenser oxyrinchus desotoi). Restriction fragment length polymorphism and control region sequence analyses of mitochondrial DNA (mtDNA) were used to identify haplotypes of Gulf sturgeon specimens obtained from eight drainages spanning the subspecies' entire distribution from Louisiana to Florida. Significant differences in haplotype frequencies indicated substantial geographic structuring of populations. A minimum of four regional or river-specific populations were identified (from west to east): (1) Pearl River, LA and Pascagoula River, MS, (2) Escambia and Yellow rivers, FI, (3) Choctawbatchee River, FL and (4) Apalachicola Ochlockonee, and Suwannee rivers, FL. Estimates of maternally mediated gene flow between any pair of the four regional or river-specific stocks ranged between 0.15 to 1.2. Tandem repeats in the mtDNA control region of Gulf sturgeon were not perfectly conserved. This result, together with an absence of heteroplasmy and length variation in Gulf sturgeon mtDNA, indicates that the molecular mechanisms of mtDNA control region sequence evolution differ among acipenserids.

Length variation, heteroplasmy and sequence divergence in the mitochondrial DNA of four species of sturgeon (Acipenser)

The extent of mtDNA length variation and heteroplasmy as well as DNA sequences of the control region and two tRNA genes were determined for four North American sturgeon species: Acipenser transmontanus, A. medirostris, A. fulvescens and A. oxyrhnychus. Across the Continental Divide, a division in the occurrence of length variation and heteroplasmy was observed that was concordant with species biogeography as well as with phylogenies inferred from restriction fragment length polymorphisms (RFLP) of whole mtDNA and pairwise comparisons of unique sequences of the control region. In all species, mtDNA length variation was due to repeated arrays of 78-82-bp sequences each containing a D-loop strand synthesis termination associated sequence (TAS). Individual repeats showed greater sequence conservation within individuals and species rather than between species, which is suggestive of concerted evolution. Differences in the frequencies of multiple copy genomes and heteroplasmy among the four species may be ascribed to differences in the rates of recurrent mutation. A mechanism that may offset the high rate of mutation for increased copy number is suggested on the basis that an increase in the number of functional TAS motifs might reduce the frequency of successfully initiated H-strand replications.

Differences in crossover frequency and distribution among three sibling species of Drosophila

Comparisons of the genetic and cytogenetic maps of three sibling species of Drosophila reveal marked differences in the frequency and cumulative distribution of crossovers during meiosis. The maps for two of these species, Drosophila melanogaster and D. simulans, have previously been described, while this report presents new map data for D. mauritiana, obtained using a set of P element markers. A genetic map covering nearly the entire genome was constructed by estimating the recombination fraction for each pair of adjacent inserts. The P-based genetic map of mauritiana is approximately 1.8 times longer than the standard melanogaster map. It appears that mauritiana has higher recombination along the entire length of each chromosome, but the difference is greates in centromere-proximal regions of the autosomes. The mauritiana autosomes show little or no centromeric recombinational suppression, a characteristic that is prominent in melanogaster. D. simulans appears to be intermediate both in terms of total map length and intensity of the autosomal centromeric effect. These interspecific differences in recombination have important evolutionary implications for DNA sequence organization and variability. In particular, mauritiana is expected to differ from melanogaster in patterns and amounts of sequence variation and transposon insertions.

Unique features in the mitochondrial D-loop region of the European seabass Dicentrarchus labrax

We have cloned and sequenced the displacement-loop (D-loop) region of the mitochondrial DNA (mtDNA) from the European seabass Dicentrarchus labrax (Dl). This sequencing revealed the presence of four tandemly repeated elements (R1, R2, R3 and R4); the individual variation in mtDNA total length is entirely accounted for by their variable number. The individuals examined also possessed an imperfect copy of one of the tandem repeats (psi R2). At least one termination-associated sequence (TAS) is present in each of the repeats and in two copies 5' upstream from the tandem array as well. The alignment of the Dl D-loop region with D-loop sequences from four other Teleosts and one Chondrosteus showed the Dl sequence to be larger than that of other fish. The extraordinary length of the Dl D-loop sequence is also due to the 5' and 3' regions that are flanking the tandem array, the largest ones to date analyzed in fish. In this study, we also report the unique organization and localization of putative TAS and conserved-sequence block (CSB) elements, and the presence of a conserved 218-bp sequence in the Dl D-loop region.

Sex-limited mitochondrial DNA transmission in the marine mussel Mytilus edulis

Mitochondrial DNA (mtDNA) was thought to be inherited maternally in animals, although paternal leakage has been reported in mice and Drosophila. Recently, direct evidence of extensive paternal inheritance of mtDNA has been found in the marine mussel Mytilus. We give evidence that whereas female mussels are homoplasmic for a genome that is transmitted to eggs, male mussels are heteroplasmic for this genome and for a second genome that is transmitted preferentially to sperm. The results provide support for the existence of separate male and female routes of mtDNA inheritance in mussels. The two genomes show a base sequence divergence exceeding 20% at three protein coding genes, consistent with long term maintenance of the heteroplasmic state. We propose that the two genomes differ in fitness in males and females, possibly as a result of interaction with nuclear genes.

Length variation in mitochondrial DNA of the minnow Cyprinella spiloptera

Length differences in animal mitochondrial DNA (mtDNA) are common, frequently due to variation in copy number of direct tandem duplications. While such duplications appear to form without great difficulty in some taxonomic groups, they appear to be relatively short-lived, as typical duplication products are geographically restricted within species and infrequently shared among species. To better understand such length variation, we have studied a tandem and direct duplication of approximately 260 bp in the control region of the cyprinid fish, Cyprinella spiloptera. Restriction site analysis of 38 individuals was used to characterize population structure and the distribution of variation in repeat copy number. This revealed two length variants, including individuals with two or three copies of the repeat, and little geographic structure among populations. No standard length (single copy) genomes were found and heteroplasmy, a common feature of length variation in other taxa, was absent. Nucleotide sequence of tandem duplications and flanking regions localized duplication junctions in the phenylalanine tRNA and near the origin of replication. The locations of these junctions and the stability of folded repeat copies support the hypothesized importance of secondary structures in models of duplication formation.

Polymorphic sequence in the D-loop region of equine mitochondrial DNA

The D-loop regions in equine mitochondrial DNA were cloned from three thoroughbred horses by polymerase chain reaction (PCR). The total number of bases in the D-loop region were 1114 bp, 1115 bp and 1146 bp. The equine D-loop region is A/T rich like many other mammalian D-loops. The large central conserved sequence block and small conserved sequence blocks 1, 2 and 3, that are common to other mammals, were observed. Between conserved sequence blocks 1 and 2 there were tandem repeats of an 8 bp equine-specific sequence TGTGCACC, and the number of tandem repeats differed among individual horses. The base composition in the unit of these repeats is G/C rich as are the short repeats in the D-loops of rabbit and pig. Comparing DNA sequences between horse and other mammals, the difference in the D-loop region length is mostly due to the difference in the number of DNA sequences at both extremities. The similarities of the DNA sequences are in the middle part of the D-loop. In comparison of the sequences among three thoroughbred horses, it was determined that the region between tRNA(Pro) and the large central conserved sequence block was the richest in variation. PCR primers in the D-loop region were designed and the expected maternal inheritance was confirmed by PCR-RFLP (restriction fragment length polymorphism).

Transfer of paternal mitochondrial DNA during fertilization of honeybee (Apis mellifera L.) eggs

Strict maternal inheritance of mitochondrial (mt) DNA is believed to be the rule in most eukaryotic organisms because of exclusion of paternal mitochondria from the egg cytoplasm during fertilization. In honeybees, polyspermic fertilization occurs, and many spermatozoa, including their mitochondria-rich flagellum, can completely penetrate the egg, thus allowing for a possibly high paternal leakage. In order to identify paternal mtDNA in honeybee eggs, restriction fragment length polymorphisms (RFLP) of different subspecies were used. Total DNA extracts of different developmental stages of an Apis mellifera carnica x Apis mellifera capensis hybrid brood were tested with a radioactively-labelled diagnostic mtDNA probe. Densitograms of autoradiographs indicated that the male contribution represents up to 27% of the total mitochondrial DNA in the fertilized eggs 12 h after oviposition. In subsequent developmental stages the portion of paternal mtDNA slowly decreased until hatching of the larvae when only traces were found. Although rapid disintegration of paternal mtDNA does not occur, the initially high paternal mitochondrial contribution is not maintained in the adult animal.

Evolution by DNA turnover in the control region of vertebrate mitochondrial DNA

The control region of animal mitochondrial DNA is heterogeneous, including both highly conserved and highly variable sequences. Within the variable regions, variable number tandem repeat sequences have been described for numerous species. Repeats at one location, just upstream of the origin of replication, show an unprecedented level of length variation in somatic tissue. Recent comparison of these sequences in different species indicates a pattern of DNA turnover acting at different rates and over motifs of various sizes.

The putative origin of heavy strand replication (oriH) in mitochondrial DNA is highly conserved among the teleost fishes

We have determined the nucleotide sequence of a 645 bp EcoRI-HindIII fragment from mitochondrial DNA of Arctic charr, Salvelinus alpinus. The sequence includes the tRNA(Phe) gene, part of the small subunit ribosomal RNA gene as well as the putative origin of heavy strand replication located at the major non-coding region, the D-loop containing region. Upon comparison of the Arctic charr sequence to mitochondrial DNA sequences from several distantly related teleost fishes, the origin of heavy strand replication was found to be highly conserved among the teleost fishes.