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

Characterization of the mitochondrial Huso huso genome and new aspects of its organization in the presence of tandem repeats in 12S rRNA

BACKGROUND

The sturgeon group has been economically significant worldwide due to caviar production. Sturgeons consist of 27 species in the world. Mitogenome data could be used to infer genetic diversity and investigate the evolutionary history of sturgeons. A limited number of complete mitogenomes in this family were sequenced. Here, we annotated the mitochondrial Huso huso genome, which revealed new aspects of this species.

RESULTS

In this species, the mitochondrial genome consisted of 13 genes encoding proteins, 22tRNA and 2rRNA, and two non-coding regions that followed other vertebrates. In addition, H. huso had a pseudo-tRNA-Glu between ND6 and Cytb and a 52-nucleotide tandem repeat with two replications in 12S rRNA. This duplication event is probably related to the slipped strand during replication, which could remain in the strand due to mispairing during replication. Furthermore, an 82 bp repeat sequence with three replications was observed in the D-loop control region, which is usually visible in different species. Regulatory elements were also seen in the control region of the mitochondrial genome, which included termination sequences and conserved regulatory blocks. Genomic compounds showed the highest conservation in rRNA and tRNA, while protein-encoded genes and nonencoded regions had the highest divergence. The mitochondrial genome was phylogenetically assayed using 12 protein-encoding genes.

CONCLUSIONS

In H. huso sequencing, we identified a distinct genome organization relative to other species that have never been reported. In recent years, along with the advancement in sequencing identified more genome rearrangements. However, it is an essential aspect of researching the evolution of the mitochondrial genome that needs to be recognized.

Species Identification of Caviar Based on Multiple DNA Barcoding

This study aimed to explore the applicability of DNA barcoding for assessing the authenticity of caviar on the Chinese market. A set of universal COI primers and two sets of designed primers based on COI and D-loop genes were used to identify maternal species of samples from 21 batches of caviar. The results showed that the PCR products from three sets of primers had more than 98% similarity to the sequences in database. The COI gene could not distinguish sturgeons with closed genetic relationships, while D-loop gene could effectively improve the accuracy of DNA barcoding and was more suitable to the identification of interspecific sturgeon than the COI gene. The neighbor-joining dendrogram further confirmed the applicability and accuracy of COI and D-loop genes in identifying maternal relatives of caviar (Acipenser baerii/Acipenser gueldenstaedtii/Acipenser schrenckii/Huso dauricus/Huso huso). Despite the limitations of mitochondrial DNA in identifying hybrid sturgeon species, the presence of counterfeit caviar of non-sturgeon ingredients could be excluded. All the caviar samples were identified successfully as sturgeon species, but the mislabeling rate of species was 33.4%, indicating that there were illegal phenomena such as disorderly labeling, mislabeling, and adulteration on the market.

Mitogenome-wise codon usage pattern from comparative analysis of the first mitogenome of Blepharipa sp. (Muga uzifly) with other Oestroid flies

Uziflies (Family: Tachinidae) are dipteran endoparasites of sericigenous insects which cause major economic loss in the silk industry globally. Here, we are presenting the first full mitogenome of Blepharipa sp. (Acc: KY644698, 15,080 bp, A + T = 78.41%), a dipteran parasitoid of Muga silkworm (Antheraea assamensis) found in the Indian states of Assam and Meghalaya. This study has confirmed that Blepharipa sp. mitogenome gene content and arrangement is similar to other Tachinidae and Sarcophagidae flies of Oestroidea superfamily, typical of ancestral Diptera. Although, Calliphoridae and Oestridae flies have undergone tRNA translocation and insertion, forming unique intergenic spacers (IGS) and overlapping regions (OL) and a few of them (IGS, OL) have been conserved across Oestroidea flies. The Tachinidae mitogenomes exhibit more AT content and AT biased codons in their protein-coding genes (PCGs) than the Oestroidea counterpart. About 92.07% of all (3722) codons in PCGs of this new species have A/T in their 3rd codon position. The high proportion of AT and repeats in the control region (CR) affects sequence coverage, resulting in a short CR (Blepharipa sp.: 168 bp) and a smaller tachinid mitogenome. Our research unveils those genes with a high AT content had a reduced effective number of codons, leading to high codon usage bias. The neutrality test shows that natural selection has a stronger influence on codon usage bias than directed mutational pressure. This study also reveals that longer PCGs (e.g., nad5, cox1) have a higher codon usage bias than shorter PCGs (e.g., atp8, nad4l). The divergence rates increase nonlinearly as AT content at the 3rd codon position increases and higher rate of synonymous divergence than nonsynonymous divergence causes strong purifying selection. The phylogenetic analysis explains that Blepharipa sp. is well suited in the family of insectivorous tachinid maggots. It's possible that biased codon usage in the Tachinidae family reduces the effective number of codons, and purifying selection retains the core functions in their mitogenome, which could help with efficient metabolism in their endo-parasitic life style and survival strategy.

Chromosome-Level Genome Assembly of the Bioluminescent Cardinalfish Siphamia tubifer: An Emerging Model for Symbiosis Research

The bioluminescent symbiosis involving the sea urchin cardinalfish Siphamia tubifer and the luminous bacterium Photobacterium mandapamensis is an emerging vertebrate model for the study of microbial symbiosis. However, little genetic data are available for the host, limiting the scope of research that can be implemented with this association. We present a chromosome-level genome assembly for S. tubifer using a combination of PacBio HiFi sequencing and Hi-C technologies. The final assembly was 1.2 Gb distributed on 23 chromosomes and contained 32,365 protein coding genes with a BUSCO score of 99%. A comparison of the S. tubifer genome to that of another nonluminous species of cardinalfish revealed a high degree of synteny, whereas a comparison to a more distant relative in the sister order Gobiiformes revealed the fusion of two chromosomes in the cardinalfish genomes. The complete mitogenome of S. tubifer was also assembled, and an inversion in the vertebrate WANCY tRNA genes as well as heteroplasmy in the length of the control region were discovered. A phylogenetic analysis based on whole the mitochondrial genome indicated that S. tubifer is divergent from the rest of the cardinalfish family, highlighting the potential role of the bioluminescent symbiosis in the initial divergence of Siphamia. This high-quality reference genome will provide novel opportunities for the bioluminescent S. tubifer-P. mandapamensis association to be used as a model for symbiosis research.

Molluscan mitochondrial genomes break the rules

The first animal mitochondrial genomes to be sequenced were of several vertebrates and model organisms, and the consistency of genomic features found has led to a 'textbook description'. However, a more broad phylogenetic sampling of complete animal mitochondrial genomes has found many cases where these features do not exist, and the phylum Mollusca is especially replete with these exceptions. The characterization of full mollusc mitogenomes required considerable effort involving challenging molecular biology, but has created an enormous catalogue of surprising deviations from that textbook description, including wide variation in size, radical genome rearrangements, gene duplications and losses, the introduction of novel genes, and a complex system of inheritance dubbed 'doubly uniparental inheritance'. Here, we review the extraordinary variation in architecture, molecular functioning and intergenerational transmission of molluscan mitochondrial genomes. Such features represent a great potential for the discovery of biological history, processes and functions that are novel for animal mitochondrial genomes. This provides a model system for studying the evolution and the manifold roles that mitochondria play in organismal physiology, and many ways that the study of mitochondrial genomes are useful for phylogeny and population biology. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

Rapid divergence, molecular evolution, and morphological diversification of coastal host-parasite systems from southern Brazil

This study assessed the role of historical processes on the geographic isolation, molecular evolution, and morphological diversification of host-parasite populations from the southern Brazilian coast. Adult specimens of Scleromystax barbatus and Scleromystax macropterus were collected from the sub-basin of the Nhundiaquara River and the sub-basin of the Paranaguá Bay, state of Paraná, Brazil. Four species of Gyrodactylus were recovered from the body surface of both host species. Morphometric analysis of Gyrodactylus spp. and Scleromystax spp. indicated that subpopulations of parasites and hosts could be distinguished from different sub-basins and locations, but the degree of morphological differentiation seems to be little related to geographic distance between subpopulations. Phylogenetic relationships based on DNA sequences of Gyrodactylus spp. and Scleromystax spp. allowed distinguishing lineages of parasites and hosts from different sub-basins. However, the level of genetic structuring of parasites was higher in comparison to host species. Evidence of positive selection in mtDNA sequences is likely associated with local adaptation of lineages of parasites and hosts. A historical demographic analysis revealed that populations of Gyrodactylus and Scleromystax have expanded in the last 250 000 years. The genetic variation of parasites and hosts is consistent with population-specific selection, population expansions, and recent evolutionary co-divergence.

Highly Resolved Phylogenetic Relationships within Order Acipenseriformes According to Novel Nuclear Markers

Order Acipenseriformes contains 27 extant species distributed across the northern hemisphere, including so-called "living fossil" species of garfish and sturgeons. Previous studies have focused on their mitochondrial genetics and have rarely used nuclear genetic data, leaving questions as to their phylogenetic relationships. This study aimed to utilize a bioinformatics approach to screen for candidate single-copy nuclear genes, using transcriptomic data from sturgeon species and genomic data from the spotted gar, Lepisosteus oculatus. We utilized nested polymerase chain reaction (PCR) and degenerate primers to identify nuclear protein-coding (NPC) gene markers to determine phylogenetic relationships among the Acipenseriformes. We identified 193 nuclear single-copy genes, selected from 1850 candidate genes with at least one exon larger than 700 bp. Forty-three of these genes were used for primer design and development of 30 NPC markers, which were sequenced for at least 14 Acipenseriformes species. Twenty-seven NPC markers were found completely in 16 species. Gene trees according to Bayesian inference (BI) and maximum likelihood (ML) were calculated based on the 30 NPC markers (20,946 bp total). Both gene and species trees produced very similar topologies. A molecular clock model estimated the divergence time between sturgeon and paddlefish at 204.1 Mya, approximately 10% later than previous estimates based on cytochrome b data (184.4 Mya). The successful development and application of NPC markers provides a new perspective and insight for the phylogenetic relationships of Acipenseriformes. Furthermore, the newly developed nuclear markers may be useful in further studies on the conservation, evolution, and genomic biology of this group.

Mind the gap! The mitochondrial control region and its power as a phylogenetic marker in echinoids

BACKGROUND

In Metazoa, mitochondrial markers are the most commonly used targets for inferring species-level molecular phylogenies due to their extremely low rate of recombination, maternal inheritance, ease of use and fast substitution rate in comparison to nuclear DNA. The mitochondrial control region (CR) is the main non-coding area of the mitochondrial genome and contains the mitochondrial origin of replication and transcription. While sequences of the cytochrome oxidase subunit 1 (COI) and 16S rRNA genes are the prime mitochondrial markers in phylogenetic studies, the highly variable CR is typically ignored and not targeted in such analyses. However, the higher substitution rate of the CR can be harnessed to infer the phylogeny of closely related species, and the use of a non-coding region alleviates biases resulting from both directional and purifying selection. Additionally, complete mitochondrial genome assemblies utilizing next generation sequencing (NGS) data often show exceptionally low coverage at specific regions, including the CR. This can only be resolved by targeted sequencing of this region.

RESULTS

Here we provide novel sequence data for the echinoid mitochondrial control region in over 40 species across the echinoid phylogenetic tree. We demonstrate the advantages of directly targeting the CR and adjacent tRNAs to facilitate complementing low coverage NGS data from complete mitochondrial genome assemblies. Finally, we test the performance of this region as a phylogenetic marker both in the lab and in phylogenetic analyses, and demonstrate its superior performance over the other available mitochondrial markers in echinoids.

CONCLUSIONS

Our target region of the mitochondrial CR (1) facilitates the first thorough investigation of this region across a wide range of echinoid taxa, (2) provides a tool for complementing missing data in NGS experiments, and (3) identifies the CR as a powerful, novel marker for phylogenetic inference in echinoids due to its high variability, lack of selection, and high compatibility across the entire class, outperforming conventional mitochondrial markers.

The complete mitochondrial genome of the Giant Manta ray, Manta birostris

The complete mitochondrial genome of the giant manta ray (Manta birostris), consists of 18,075 bp with rich A + T and low G content. Gene organization and length is similar to other species of ray. It comprises of 13 protein-coding genes, 2 rRNAs genes, 23 tRNAs genes and 1 non-coding sequence, and the control region. We identified an AT tandem repeat region, similar to that reported in Mobula japanica.

Current Occurrence of the Atlantic Sturgeon Acipenser oxyrinchus in Northern Spain: A New Prospect for Sturgeon Conservation in Western Europe

Acipenser oxyrinchus is considered extirpated in Europe, but numerous breeding populations still exist on the Atlantic coast of North America. An adult female A. oxyrinchus, 2500 mm total length and 120 kg wet weight, was accidentally fished on 24 November 2010 near the coast of Gijón, Asturias, Spain. The fish was identified by its morphological pattern as well as by mitochondrial and nuclear DNA analyses. Because the sturgeon was found far away from any known breeding area, it was considered a stray or vagrant specimen. It certainly has a natural origin, but its eventual birthplace could not be determined. Because its current occurrence was unknown in southwestern Europe until now, the species is not cataloged or protected in this area. Therefore, the residual European stocks of A. oxyrinchus ought to be listed as Critically Endangered (CR) according to the IUCN categories. Likewise, it is imperative for southwestern European countries with an historic or recent occurrence of A. oxyrinchus to protect the species through domestic and international legislation. The present sympatric occurrence of A. sturio and A. oxyrinchus raises new challenges about key questions, such as the species selection for restoration program in European countries. Accurate monitoring is mandatory to obtain appropriate information for an assessment of the current occurrence of A. oxyrinchus in southwestern Europe.

Selections that isolate recombinant mitochondrial genomes in animals

Homologous recombination is widespread and catalyzes evolution. Nonetheless, its existence in animal mitochondrial DNA is questioned. We designed selections for recombination between co-resident mitochondrial genomes in various heteroplasmic Drosophila lines. In four experimental settings, recombinant genomes became the sole or dominant genome in the progeny. Thus, selection uncovers occurrence of homologous recombination in Drosophila mtDNA and documents its functional benefit. Double-strand breaks enhanced recombination in the germline and revealed somatic recombination. When the recombination partner was a diverged Drosophila melanogaster genome or a genome from a different species such as Drosophila yakuba, sequencing revealed long continuous stretches of exchange. In addition, the distribution of sequence polymorphisms in recombinants allowed us to map a selected trait to a particular region in the Drosophila mitochondrial genome. Thus, recombination can be harnessed to dissect function and evolution of mitochondrial genome.

Assembly and variation analyses of Clarias batrachus mitogenome retrieved from WGS data and its phylogenetic relationship with other catfishes

Whole genome sequencing (WGS) using next generation sequencing technologies paves the way to sequence the mitochondrial genomes with greater ease and lesser time. Here, we used the WGS data of Clarias batrachus, generated from Roche 454 and Ion Torrent sequencing platforms, to assemble the complete mitogenome using both de novo and reference based approaches. Both the methods yielded almost similar results and the best assembled mitogenome was of 16,510 bp size (GenBank Acc. No. KM259918). The mitogenome annotation resulted in 13 coding genes, 22 tRNA genes, 2 rRNA genes and one control region, and the gene order was found to be identical with other catfishes. Variation analyses between assembled and the reference (GenBank Acc. No. NC_023923) mitogenome revealed 51 variations. The phylogenetic analysis of coding DNA sequences and tRNA supports the monophyly of catfishes. Two SSRs were identified in C. batrachus mitogenome, out of which one was unique to this species. Based on the relative rate of gene evolution, protein coding mitochondrial genes were found to evolve at a much faster pace than the d-loop, which in turn are followed by the rRNAs; the tRNAs showed wide variability in the rate of sequence evolution, and on average evolve the slowest. Among the coding genes, ND2 evolves most rapidly. The variations present in the coding regions of the mitogenome and their comparative analyses with other catfish species may be useful in species conservation and management programs.

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.

Site heteroplasmy in the mitochondrial cytochrome b gene of the sterlet sturgeon Acipenser ruthenus

Sturgeons are fish species with a complex biology. They are also characterized by complex aspects including polyploidization and easiness of hybridization. As with most of the Ponto-Caspian sturgeons, the populations of Acipenser ruthenus from the Danube have declined drastically during the last decades. This is the first report on mitochondrial point heteroplasmy in the cytochrome b gene of this species. The 1141 bp sequence of the cytb gene in wild sterlet sturgeon individuals from the Lower Danube was determined, and site heteroplasmy evidenced in three of the 30 specimens collected. Two nucleotide sequences were identified in these heteroplasmic individuals. The majority of the heteroplasmic sites are synonymous and do not modify the sequence of amino acids in cytochrome B protein. To date, several cases of point heteroplasmy have been reported in animals, mostly due to paternal leakage of mtDNA. The presence of specific point heteroplasmic sites might be interesting for a possible correlation with genetically distinct groups in the Danube River.

Complete mitochondrial DNA sequence of the ark shell Scapharca broughtonii: an ultra-large metazoan mitochondrial genome

The complete mitochondrial (mt) genome of the ark shell Scapharca broughtonii was determined using long PCR and a genome walking sequencing strategy with genus-specific primers. The S. broughtonii mt genome (GenBank accession number AB729113) contained 12 protein-coding genes (the atp8 gene is missing, as in most bivalves), 2 ribosomal RNA genes, and 42 transfer tRNA genes, in a length of 46,985 nucleotides for the size of mtDNA with only one copy of the heteroplasmic tandem repeat (HTR) unit. Moreover the S. broughtonii mt genome shows size variation; these genomes ranged in size from about 47 kb to about 50 kb because of variation in the number of repeat sequences in the non-coding region. The mt-genome of S. broughtonii is, to date, the longest reported metazoan mtDNA sequence. Sequence duplication in non-coding region and the formation of HTR arrays were two of the factors responsible for the ultra-large size of this mt genome. All the tRNA genes were found within the S. broughtonii mt genome, unlike the other bivalves usually lacking one or more tRNA genes. Twelve additional specimens were used to analyze the patterns of tandem repeat arrays by PCR amplification and agarose electrophoresis. Each of the 12 specimens displayed extensive heteroplasmy and had 8-10 length variants. The motifs of the HTR arrays are about 353-362 bp and the number of repeats ranges from 1 to 11.

Structure and organization of the mitochondrial DNA control region with tandemly repeated sequence in the Amazon ornamental fish

Tandemly repeated sequences are a common feature of vertebrate mitochondrial DNA control regions. However, questions still remain about their mode of evolution and function. To better understand patterns of variation in length and to explore the existence of previously described domain, we have characterized the control region structure of the Amazonian ornamental fish Nannostomus eques and Nannostomus unifasciatus. The control region ranged from 1121 to 1142 bp in length and could be separated into three domains: the domain associated with the extended terminal associated sequences, the central conserved domain, and the conserved sequence blocks domain. In the first domain, we encountered a sequence repeated 10 times in tandem (variable number tandem repeat (VNTR)) that could adopt an "inverted repetitions" type structural conformation. The results suggest that the VNTR pattern encountered in both N. eques and N. unifasciatus is consistent with the prerequisites of the illegitimate elongation model in which the unequal pairing of the chains near the 5'-end of the control region favors the formation of repetitions.

The complete mitochondrial genome of the taimen, Hucho taimen, and its unusual features in the control region

The whole mitochondrial genome of Hucho taimen was firstly sequenced and characterized. The genome is 16,833 bp in length and contains 13 protein-coding genes, 2 ribosomal RNAs, 22 transfer RNAs, and a noncoding control region. Twelve protein-coding genes on the heavy strand showed that the content of A+T was higher than that of G+C, whereas the nd6 protein-coding gene on the light strand displayed an opposite pattern. We described the secondary structure of the origin of light strand (oriL) replication and found that the conserved 5'-GCCGG-3' sequence motif is variable in H. taimen and some other salmonids. We conclude that the control region is variable in length and represents the high A+T content, compared with other mitochondrial control regions available in Salmonidae and other non-salmonids. Additionally, another interesting feature of H. taimen mitogenome is that a T-type mononucleotide microsatellite and an 82 bp tandem repeat were identified in the control region.

Mitochondrial DNA STR analysis as a tool for studying the green sea turtle (Chelonia mydas) populations: the Mediterranean Sea case study

The Mediterranean population of the green sea turtle Chelonia mydas is critically endangered. Genetic analysis of this population using the ordinary haplotyping system, based on sequence analysis of a segment of the mitochondrial DNA (mtDNA) D-loop (control region), revealed very little variation. The most common haplotype, CM-A13, was observed in all but three individuals in hundreds of samples in previous studies. In search for a more informative marker we sequenced the 3' of the mitochondrial control region which contains an AT-rich microsatellite. We found a unique pattern that consists of four AT short tandem repeats (STRs) with varying copy numbers. This allowed us to construct a new haplotyping system composed of four different STR sizes for each mtDNA sequence. Our new mitochondrial STR (mtSTR) haplotyping approach revealed 33 different haplotypes within the nesting and stranded sea turtles along the Mediterranean Israeli seashore. The Israeli coast nesting females had 10 different haplotypes that can be used for monitoring and conservation purposes. The mtSTR haplotyping system can clearly assist in fingerprinting of individual turtles. Moreover, it can be used for estimating phylogenetic distances within populations. This case study shows that the mtSTR haplotyping is applicable for the study of global green sea turtle populations and could also be considered as markers of genetic variability in other species.

Variation in repeat length and heteroplasmy of the mitochondrial DNA control region along a core-edge gradient in the eastern spadefoot toad (Pelobates syriacus)

Peripheral populations are those situated at the distribution margins of a species and are often subjected to more extreme abiotic and biotic conditions than those at the core. Here, we hypothesized that shorter repeat length and fewer heteroplasmic mitochondrial DNA (mtDNA) copies, which are associated with more efficient mitochondrial function, may be related to improved survival under extreme environmental conditions. We sampled eastern spadefoot toads (mostly as tadpoles) from 43 rain pools distributed along a 300-km gradient from core to edge of the species' distribution. We show that mean pool tandem repeat length and heteroplasmy increase from edge to core, even after controlling for body size. We evaluate several alternative hypotheses and propose the Fisher hypothesis as the most likely explanation. However, additional sequential sampling and experimental studies are required to determine whether selection under extreme conditions, or alternative mechanisms, could account for the gradient in heteroplasmy and repeat length in the mtDNA control region.

Managing polyploidy in ex situ conservation genetics: the case of the critically endangered Adriatic sturgeon (Acipenser naccarii)

While the current expansion of conservation genetics enables to address more efficiently the management of threatened species, alternative methods for genetic relatedness data analysis in polyploid species are necessary. Within this framework, we present a standardized and simple protocol specifically designed for polyploid species that can facilitate management of genetic diversity, as exemplified by the ex situ conservation program for the tetraploid Adriatic sturgeon Acipenser naccarii. A critically endangered endemic species of the Adriatic Sea tributaries, its persistence is strictly linked to the ex situ conservation of a single captive broodstock currently decimated to about 25 individuals, which represents the last remaining population of Adriatic sturgeon of certain wild origin. The genetic variability of three F1 broodstocks available as future breeders was estimated based on mitochondrial and microsatellite information and compared with the variability of the parental generation. Genetic data showed that the F1 stocks have only retained part of the genetic variation present in the original stock due to the few parent pairs used as founders. This prompts for the urgent improvement of the current F1 stocks by incorporating new founders that better represent the genetic diversity available. Following parental allocation based on band sharing values, we set up a user-friendly tool for selection of candidate breeders according to relatedness between all possible parent-pairs that secures the use of non-related individuals. The approach developed here could also be applied to other endangered tetraploid sturgeon species overexploited for caviar production, particularly in regions lacking proper infrastructure and/or expertise.

Repetitive sequences in the lamprey mitochondrial DNA control region and speciation of Lethenteron

The sequence of the mitochondrial DNA control region was examined in four species of lamprey in the genus Lethenteron. The 3' half of the control region contains highly variable repeat sequences, showing variation in both copy number and nucleotide sequence, even within local populations. Detailed analyses of the sequences of the repeats allowed us to deduce that slipped-strand mispairing during DNA replication, accompanied by a high rate of substitutions and indels, was primarily responsible for the variation in the repeats. We also found that some cases might be better explained by gene conversion, due to intermolecular recombination. Based on the observed variable nature of the mitochondrial control region, we searched for molecular markers in mitochondrial DNA, because there are few fixed genetic markers for distinguishing between Lethenteron japonicum and Lethenteron kessleri. However, we found no reliable markers in the control region. No fixed substitution was observed in intron sequences of the nuclear gene SoxD. Thus, these two species likely diverged quite recently and may possess only a limited number of fixed genetic loci.

Structure and conservation of tandem repeats in the mitochondrial DNA control region of the Least Brook lamprey (Lampetra aepyptera)

Tandemly repeated sequences are a common feature of vertebrate mitochondrial DNA control regions. However, questions still remain about their mode of evolution, function, and phylogenetic distribution. We report phylogenetic and geographic patterns of variation of control region repeat sequence and number in a nonparasitic lamprey, Lampetra aepyptera. A survey of populations from throughout the species' range revealed remarkably low repeat sequence polymorphism but some interpopulation variation in repeat number. The high sequence similarity extended to repeats observed in other species in the genus Lampetra and other lamprey genera. The very low levels of variation suggest a high copy turnover. Our data are consistent with the illegitimate elongation model of repeat gain and loss and further suggest that repeat change occurs at internal copies. However, the limited variation across some species of lamprey suggests that functional constraints may further limit variation.

Tracing the first steps of American sturgeon pioneers in Europe

BACKGROUND

A Baltic population of Atlantic sturgeon was founded approximately 1,200 years ago by migrants from North America, but after centuries of persistence, the population was extirpated in the 1960s, mainly as a result of over-harvest and habitat alterations. As there are four genetically distinct groups of Atlantic sturgeon inhabiting North American rivers today, we investigated the genetic provenance of the historic Baltic population by ancient DNA analyses using mitochondrial and nuclear markers.

RESULTS

The phylogeographic signal obtained from multilocus microsatellite DNA genotypes and mitochondrial DNA control region haplotypes, when compared to existing baseline datasets from extant populations, allowed for the identification of the region-of-origin of the North American Atlantic sturgeon founders. Moreover, statistical and simulation analyses of the multilocus genotypes allowed for the calculation of the effective number of individuals that originally founded the European population of Atlantic sturgeon. Our findings suggest that the Baltic population of A. oxyrinchus descended from a relatively small number of founders originating from the northern extent of the species' range in North America.

CONCLUSION

These results demonstrate that the most northerly distributed North American A. oxyrinchus colonized the Baltic Sea approximately 1,200 years ago, suggesting that Canadian specimens should be the primary source of broodstock used for restoration in Baltic rivers. This study illustrates the great potential of patterns obtained from ancient DNA to identify population-of-origin to investigate historic genotype structure of extinct populations.

Rapid evolution of the intergenic T-P spacer in the mtDNA of Arctic cod Arctogadus glacialis

A noncoding intergenic spacer has previously been reported in mtDNA of Gadiformes. Here we present sequence information from two other cod species and variation within three species to clarify the evolution of this region. A general feature of the T-P spacer is high variation and folding into two or three hairpins. The variation among species both in structure of the region and sequence variation reflects the phylogenetic relationship of the species. A unique pattern is found within Arctic cod, Arctogadus glacialis, in which tandem repeat motifs result in new stable secondary structures. There is large variation in size of the region both within (heteroplasmy) and among individuals. A duplicated insertion is found in Greenland cod, Gadus ogac, at the same position as a corresponding duplication in Atlantic cod, Gadus morhua.

The complete mitochondrial genome of the Senegal sole, Solea senegalensis Kaup. Comparative analysis of tandem repeats in the control region among soles

The complete nucleotide sequence of the mitochondrial genome for the Senegal sole Solea senegalensis Kaup was determined. The mitochondrial DNA was 16,659 base pairs (bp) in length. Sequence features of the 13 protein-coding genes, two ribosomal RNAs and 22 transfer RNAs are described. The non-coding control region (1017 bp) was compared with those of the closely related soles Solea solea and Solea lascaris. The typical conservative blocks were identified. A cluster of 42 and 22 tandemly arrayed repeats was detected near the 3' end of control region in S. solea and S. lascaris, respectively. On the contrary, only two (93.8% of haplotypes) or three copies (6.2%) of an 8-bp repeated sequence motif was found in S. senegalensis. Phylogenetic analysis showed that 7 out of 9 of haplotypes bearing three copies grouped in a separate cluster. Possible mechanisms influencing the evolution of control region among soles are discussed.

Comparison of karyotypes of Acipenser oxyrinchus and A. sturio by chromosome banding and fluorescent in situ hybridization

A highly debated problem in Acipenseridae taxonomy is whether Acipenser oxyrinchus (North American Atlantic sturgeon) and A. sturio (European Atlantic sturgeon) are true species: a detailed comparison of their karyotypes could provide relevant information. Here we describe for the first time the karyotype of A. oxyrinchus (2n = 121 +/- 3), and its features, among which the constitutive heterochromatin, revealed by C-banding technique, the distribution of telomeric regions, and the 5S rRNA genes, detected by FISH. The results reveal that A. oxyrinchus and A. sturio karyotypes and features are quite similar. Moreover, comparing the results obtained through hybridization by FISH with HindIII and PstI satellite DNA in these and in other sturgeon species, no hybridization signals are detected in A. sturio and A. oxyrinchus, while A. stellatus and A. gueldenstaedtii show hybridization. Thus A. sturio and A. oxyrinchus appear very similar from a cytogenetic point of view: these and molecular data repeatedly point out that A. sturio and A. oxyrinchus represent a sister clade in comparison to all other sturgeon species up to now studied.

Anatomical, hormonal and histological descriptions of captive Russian sturgeon (Acipenser gueldenstaedtii) with intersex gonads

Sturgeons are known throughout the world as the source of black caviar. Their declining populations in their native habitats, mainly the Caspian Sea, due to over-fishing for meat and caviar production, destruction of their spawning grounds and water pollution, have led to their introduction into aquaculture in areas with suitable conditions, including Israel. Recently, we noticed an unusual phenomenon in these normally gonochoristic species. Several 5-year-old female sturgeons were found to have one or more testicular sections in each of their two gonads, forming an intersexual gender. Further examination of other fish from the same age group revealed 14% fish with intersex gonads among a population of 5000 fish that had been pre-selected as females. This phenomenon has not been found however in other age groups of Russian sturgeons, cultured at the same facility. Sturgeons are a generally gonochoristic species, and hermaphroditism is only very infrequently observed under natural or normal breeding conditions. Moreover, these rare cases have all been from polluted habitats. The present work is the first description of fish containing intersex gonads in Russian sturgeon (Acipenser gueldenstaedtii). We describe the phenomenon anatomically and histologically, and examine plasma steroid levels and pituitary gonadotropin gene expression by comparing fish with intersex gonads with normal females and males of the same age group. Intersex gonads were typical female ovaries with one or more white testicular components embedded in each. The testis components were not uniform in size or location among the two gonads of each fish or among different fish, and they showed marked differences in distribution. The ovarian component of the intersex gonad was at the pre-vitellogenic stage as in normal females, and the testis component contained spermatids and mature spermatozoa as in normal males of the same age. However, in terms of estradiol and 11-ketotestosterone plasma levels, as well as of pituitary gonadotropin (betaLH and betaFSH) gene-expression levels, the fish with intersex gonads were more similar to the normal males than to the normal females, even though the testis part of the intersex gonad was smaller than the ovarian part. To examine the possibility that the fish containing intersex gonads were hybrids, phylogenetic trees were constructed from the consensus sequences of Cytochrome b and control region (D-loop) genes. Results indicated no differences between the fish with intersex gonads and normal males or females of the same age group. However, statistically significant differences were found between different age groups of Russian sturgeon, as well as of white sturgeons (A. transmontanus), grown under the same culture conditions.

Unusual intraindividual variation of the nuclear 18S rRNA gene is widespread within the Acipenseridae

Significant intraindividual variation in the sequence of the 18S rRNA gene is unusual in animal genomes. In a previous study, multiple 18S rRNA gene sequences were observed within individuals of eight species of sturgeon from North America but not in the North American paddlefish, Polyodon spathula, in two species of Polypterus (Polypterus delhezi and Polypterus senegalus), in other primitive fishes (Erpetoichthys calabaricus, Lepisosteus osseus, Amia calva) or in a lungfish (Protopterus sp.). These observations led to the hypothesis that this unusual genetic characteristic arose within the Acipenseriformes after the presumed divergence of the sturgeon and paddlefish families. In the present study, a survey of nearly all Eurasian acipenseriform species was conducted to examine 18S rDNA variation. Intraindividual variation was not found in the polyodontid species, the Chinese paddlefish, Psephurus gladius, but variation was detected in all Eurasian acipenserid species. The comparison of sequences from two major segments of the 18S rRNA gene and identification of sites where insertion/deletion events have occurred are placed in the context of evolutionary relationships within the Acipenseriformes and the evolution of rDNA variation in this group.

Structure and evolution of the mitochondrial control region and flanking sequences in the European cave salamander Proteus anguinus

The European cave salamander Proteus anguinus Laurenti 1768 is one of the best-known subterranean animals, yet its evolutionary history and systematic relationships remain enigmatic. This is the first comprehensive study on molecular evolution within the taxon, using an mtDNA segment containing the control region (CR) and adjacent sequences. Two to seven tandem repeats of 24-32 bp were found in the intergenic spacer region (VNTR1), and three, four or six repeats, 59-77 bp each, in the 3' end of the CR (VNTR2). Different molecular mechanisms account for VNTR2 formation in different lineages of Proteus. The overall CR variation was lower than that of the spacer region, the 3' end of the cytb gene, or the tRNA genes. Individual genes and the concatenated non-repetitive sequences produced similar, well resolved maximum likelihood, Bayesian inference and parsimony trees. The numbers of repeat elements as well as the genealogy of the VNTR2 repeat units were mostly inconsistent with the groupings of the non-repetitive sequences. Different degrees of repeat array homogenization were detected in all major groups. Orthology was established for the first and the second VNTR2 elements of some populations. These two copies may therefore be used for analyses at the population level. The pattern of CR sequence variation points to strong genetic isolation of hydrographically separated populations. Genetic separation of the major groups of populations is incongruent with the current division into subspecies.

Minding the gap: frequency of indels in mtDNA control region sequence data and influence on population genetic analyses

Insertions and deletions (indels) result in sequences of various lengths when homologous gene regions are compared among individuals or species. Although indels are typically phylogenetically informative, occurrence and incorporation of these characters as gaps in intraspecific population genetic data sets are rarely discussed. Moreover, the impact of gaps on estimates of fixation indices, such as F(ST), has not been reviewed. Here, I summarize the occurrence and population genetic signal of indels among 60 published studies that involved alignments of multiple sequences from the mitochondrial DNA (mtDNA) control region of vertebrate taxa. Among 30 studies observing indels, an average of 12% of both variable and parsimony-informative sites were composed of these sites. There was no consistent trend between levels of population differentiation and the number of gap characters in a data block. Across all studies, the average influence on estimates of PhiST was small, explaining only an additional 1.8% of among population variance (range 0.0-8.0%). Studies most likely to observe an increase in PhiST with the inclusion of gap characters were those with < 20 variable sites, but a near equal number of studies with few variable sites did not show an increase. In contrast to studies at interspecific levels, the influence of indels for intraspecific population genetic analyses of control region DNA appears small, dependent upon total number of variable sites in the data block, and related to species-specific characteristics and the spatial distribution of mtDNA lineages that contain indels.

Evidence for recombination of mitochondrial DNA in triploid crucian carp

In this study, we report the complete mitochondrial DNA (mtDNA) sequences of the allotetraploid and triploid crucian carp and compare the complete mtDNA sequences between the triploid crucian carp and its female parent Japanese crucian carp and between the triploid crucian carp and its male parent allotetraploid. Our results indicate that the complete mtDNA nucleotide identity (98%) between the triploid crucian carp and its male parent allotetraploid was higher than that (93%) between the triploid crucian carp and its female parent Japanese crucian carp. Moreover, the presence of a pattern of identity and difference at synonymous sites of mitochondrial genomes between the triploid crucian carp and its parents provides direct evidence that triploid crucian carp possessed the recombination mtDNA fragment (12,759 bp) derived from the paternal fish. These results suggest that mtDNA recombination was derived from the fusion of the maternal and paternal mtDNAs. Compared with the haploid egg with one set of genome from the Japanese crucian carp, the diploid sperm with two sets of genomes from the allotetraploid could more easily make its mtDNA fuse with the mtDNA of the haploid egg. In addition, the triple hybrid nature of the triploid crucian carp probably allowed its better mtDNA recombination. In summary, our results provide the first evidence of mtDNA combination in polyploid fish.

Genomic organization and evolution of the 5S ribosomal DNA in the ancient fish sturgeon

Ribosomal DNA in sturgeon is informative when analyzed at the molecular level because it bears unique characteristics that are, to a certain extent, ancestral within vertebrates. In this paper, we examine the structure and the molecular evolution of the 5S ribosomal DNA (rDNA) region in 13 sturgeon species, comparing both the 5S ribosomal RNA (rRNA) genes and the non-transcribed spacer (NTS) sequences between the coding regions. We have found that different NTS and 5S gene variants are intermixed in the 5S rDNA arrays of the different sturgeon species and that all variants are ancestral, having been maintained over many millions of years. Using predictive models, we have found similar levels of sequence diversity in the coding regions, as well as in the non-coding region, but fixed interspecific differences are underrepresented for 5S genes. However, contrary to the expectations, we have not found fixed differences between NTS sequences when comparing many pairs of species. Specifically, when they belong to the same phylogeographic clade of the four into which the sturgeon is divided, but fixation of mutations and divergence is found between species belonging to different phylogeographic clades. Our results suggest that the evolution of the two parts of the 5S rDNA region cannot be explained exclusively as the outcome of a balance between mutational, homogenizing (i.e., gene conversion as a predominant force in sturgeon), and selective forces. Rather, they suggest that other factors (i.e., hybridization) might be superimposed over those forces and thus could to some extent be masking their effects.

Phylogeny of elasmobranchs based on LSU and SSU ribosomal RNA genes

The dominant view of the phylogeny of living elasmobranchs, based on morphological characters, is that batoids (skates and rays) are derived sharks, joined with saw sharks, and angel sharks in the clade Hypnosqualea [S. Shirai, Squalean Phylogeny: A New Framework of 'Squaloid' Sharks and Related Taxa, Hokkaido University Press, Sapporo, 1992]. By contrast, a recent molecular-phylogenetic study based on mitochondrial genes for 12S and 16S rRNA and tRNA valine [C.J. Douady et al., Mol. Phylogenet. Evol., 26 (2003) 215-221] supported the older view that batoids and sharks are separate lineages. Here, we tested these two different views using combined, nuclear large-subunit and small-subunit rRNA gene sequences ( approximately 5.3kb) from 22 elasmobranchs, two chimeras, and two bony fishes. We used maximum likelihood, maximum parsimony, minimum evolution, and Bayesian inference for tree reconstruction, and found the large-subunit rRNA gene to contain far more signal than the small-subunit gene for resolving this mostly Mesozoic radiation. Our findings matched those of in separating batoids from sharks and in statistically rejecting Hypnosqualea. The angel shark (Squatina) was the sister group to squaliforms (dogfish sharks), and our findings are consistent with the idea that "orbitostylic" sharks form a monophyletic group (squaliforms+the hexanchiform Chlamydoselachus+Squatina+Pristiophorus). In the galeomorph sharks, however, lamniforms grouped with orectolobiforms, opposing the widely accepted 'lamniform+carcharhiniform' grouping. A tree based on the mitochondrial gene for cytochrome b also supported a separation of sharks and batoids, in contrast to Hypnosqualea. Among elasmobranchs, variation in the evolutionary rates of the nuclear rRNA genes was higher than that of cytochrome b genes, mainly due to the relatively rapid evolution of rRNA in some carcharhiniforms. In conclusion, several different molecular studies now refute the Hypnosqualea hypothesis of elasmobranch interrelationships.

Restriction site heteroplasmy in the mitochondrial DNA of Brycon opalinus (Cuvier, 1819) (Characiformes, Characidae, Bryconiae)

Homoplasmy is a feature usually found in the mtDNA of higher animal taxa. On the other hand, the presence of two classes of mtDNA in the same cell or organism is rare and may appear in length or site variation. Data from mtDNA RFLP analysis of Brycon opalinus populations (Cuvier, 1819; Characiformes, Characidae, Bryconinae) revealed site heteroplasmy from endonuclease NheI digestion. Southern blotting hybridization was used to survey a total of 257 specimens with 24 restriction enzymes. Three different restriction fragment patterns of mtDNA were obtained from NheI digestion. Two individuals from hatchery broodstock were found to have two of them. NheI digests of heteroplasmic individuals yielded two fragments of approximately 1180 and 1260 bp. Despite the low frequency of this type of heteroplasmy in the whole B. opalinus population, the presence of site heteroplasmy in this species supports the evidence of this phenomenon in lower vertebrate groups.

Nonconcordant evolutionary history of maternal and paternal lineages in Adriatic sturgeon

Although analyses of intraspecific variability are an important prerequisite for species identification assays, only a few studies have focused on population genetics and historical biogeography of sturgeon species. Here we present the first study on genetic variability of the last remaining Adriatic sturgeon, Acipenser naccarii, derived from mitochondrial and nuclear DNA. Our mitochondrial DNA analyses arranged individuals into three distinguished mitochondrial DNA haplogroups (Po1, Po2 and Buna). Two haplogroups (Po1 and Buna) were correlated to geographical distribution, whereas the third (Po2) was not. It was, however, very closely related to one lineage of its Ponto-Caspian sister species, A. gueldenstaedtii. The distribution of nuclear markers (microsatellites and amplified fragment length polymorphism) was strongly correlated to geographical distribution. An assignment test based on nuclear data placed no specimen of A. naccarii to A. gueldenstaedtii and vice versa. Therefore, the presence of gueldenstaedtii-like haplotypes within the Po population is either the result of a postglacial introgression or an ancestral polymorphism and does not indicate a hybrid population. The most valuable tool for forensic species identification purposes is one diagnostic deletion separating all A. naccarii from A. gueldenstaedtii. As both A. naccarii populations are genetically differentiated, stocking of sturgeon from the Po River in Italy into waters of the Buna River would jeopardize the genetic differences between both populations and should thus be avoided.

The crocodilian mitochondrial control region: general structure, conserved sequences, and evolutionary implications

We present the first comprehensive analysis of the crocodilian control region. We have analyzed sequences from all three families of Crocodylia (Crocodylidae, Gavialidae, Alligatoridae), incorporating all genera except Paleosuchus and Melanosuchus. Within the control region of other vertebrates, several sequence motifs and their order appear to be conserved. Herein, we compare aligned crocodilian D-loop sequences to homologous sequences from other vertebrates ranging from fish to birds. Among other findings, we have discovered that while domain I tends to be shorter than the same region in mammals and birds, it contains sequences similar in structure to both the goose-hairpin and termination associated sequences (TAS). Domain II is highly conservative with regard to size among the taxa examined and contains several of the conserved sequence boxes characterized in other vertebrates. Domain III contains several interesting sequence motifs including tandemly repeated sequences, a long poly-A region in the Crocodylidae, and possible bidirection promoter sequences.