Unbiased estimation of symmetrical directional mutation pressure from protein-coding DNA

Mitochondrial genomes and Doubly Uniparental Inheritance: new insights from Musculista senhousia sex-linked mitochondrial DNAs (Bivalvia Mytilidae)

BACKGROUND

Doubly Uniparental Inheritance (DUI) is a fascinating exception to matrilinear inheritance of mitochondrial DNA (mtDNA). Species with DUI are characterized by two distinct mtDNAs that are inherited either through females (F-mtDNA) or through males (M-mtDNA). DUI sex-linked mitochondrial genomes share several unusual features, such as additional protein coding genes and unusual gene duplications/structures, which have been related to the functionality of DUI. Recently, new evidence for DUI was found in the mytilid bivalve Musculista senhousia. This paper describes the complete sex-linked mitochondrial genomes of this species.

RESULTS

Our analysis highlights that both M and F mtDNAs share roughly the same gene content and order, but with some remarkable differences. The Musculista sex-linked mtDNAs have differently organized putative control regions (CR), which include repeats and palindromic motifs, thought to provide sites for DNA-binding proteins involved in the transcriptional machinery. Moreover, in male mtDNA, two cox2 genes were found, one (M-cox2b) 123bp longer.

CONCLUSIONS

The complete mtDNA genome characterization of DUI bivalves is the first step to unravel the complex genetic signals allowing Doubly Uniparental Inheritance, and the evolutionary implications of such an unusual transmission route in mitochondrial genome evolution in Bivalvia. The observed redundancy of the palindromic motifs in Musculista M-mtDNA may have a role on the process by which sperm mtDNA becomes dominant or exclusive of the male germline of DUI species. Moreover, the duplicated M-COX2b gene may have a different, still unknown, function related to DUI, in accordance to what has been already proposed for other DUI species in which a similar cox2 extension has been hypothesized to be a tag for male mitochondria.

Nucleotide Content Gradients in Maternally and Paternally Inherited Mitochondrial Genomes of the Mussel Mytilus

Several studies have shown that in vertebrate mtDNAs the nucleotide content at fourfold degenerate sites is well correlated with the site's time of exposure to the single-strand state, as predicted from the asymmetrical model of mtDNA replication. Here we examine whether the same explanation may hold for the regional variation in nucleotide content in the maternal and paternal mtDNAs of the mussel Mytilus galloprovincialis. The origin of replication of the heavy strand (O(H)) of these genomes has been previously established. A systematic search of the two genomes for sequences that are likely to act as the origin of replication of the light strand (O(L)) suggested that the most probable site lies within the ND3 gene. By adopting this O(L) position we calculated times of exposure for 0(FD) (nondegenerate), 2(FD) (twofold degenerate), and 4(FD) (fourfold degenerate) sites of the protein-coding part of the genome and for the rRNA, tRNA and noncoding parts. The presence of thymine and absence of guanine at 4(FD) sites was highly correlated with the presumed time of exposure. Such an effect was not found for the 2(FD) sites, the rRNA, the tRNA, or the noncoding parts. There was a trend for a small increase in cytosine at 0(FD) sites with exposure time, which is explicable as the result of biased usage of 4(FD) codons. The same analysis was applied to a recently sequenced mitochondrial genome of Mytilus trossulus and produced similar results. These results are consistent with the asymmetrical model of replication and suggest that guanine oxidation due to single-strand exposure is the main cause of regional variation of nucleotide content in Mytilus mitochondrial genomes.

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

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

The Complete Maternal and Paternal Mitochondrial Genomes of the Mediterranean Mussel Mytilus galloprovincialis: Implications for the Doubly Uniparental Inheritance Mode of mtDNA

The maternal (F) and paternal (M) mitochondrial genomes of the mussel Mytilus galloprovincialis have diverged by about 20% in nucleotide sequence but retained identical gene content and gene arrangement and similar nucleotide composition and codon usage bias. Both lack the ATPase8 subunit gene, have two tRNAs for methionine and a longer open-reading frame for cox3 than seen in other mollusks. Between the F and M genomes, tRNAs are most conserved followed by rRNAs and protein-coding genes, even though the degree of divergence varies considerably among the latter. Divergence at nad3 is exceptionally low most likely because this gene includes the origin of transcription of the lagging strand (O(L)). Noncoding regions are the least conserved with the notable exception of the central domain of the main control region and a segment of another noncoding region immediately following nad3. The amino acid divergence (14%) of the two genomes is smaller than in two other pairs of conspecific genomes that are available in GenBank, that of the clam Venerupis philippinarum (34%) and of the fresh water mussel Inversidens japanensis (50%), suggesting that doubly uniparental inheritance of mtDNA emerged at different times in the three species or that there has been a relatively recent replacement of the male genome by the female in the Mytilus line. The latter hypothesis is supported from phylogenetic and population studies of Mytilidae. That the M genome contains a full complement of genes with no premature termination codons argues against it being a selfish element that rides with the sperm. It is shorter than the F by 118 bp, which apparently cannot account for the postulated replicative advantage of this genome over the F in male gonads. The high similarity of the two genomes explains why the F genome may assume the role of the M genome, but it does not exclude the possibility that for this to happen some M-specific sequences must be transferred on to the F genome by means of recombination. If such sequences exist they would most likely be located in noncoding regions.

Complete mitochondrial DNA genome sequences of extinct birds: ratite phylogenetics and the vicariance biogeography hypothesis

The ratites have stimulated much debate as to how such large flightless birds came to be distributed across the southern continents, and whether they are a monophyletic group or are composed of unrelated lineages that independently lost the power of flight. Hypotheses regarding the relationships among taxa differ for morphological and molecular data sets, thus hindering attempts to test whether plate tectonic events can explain ratite biogeography. Here, we present the complete mitochondrial DNA genomes of two extinct moas from New Zealand, along with those of five extant ratites (the lesser rhea, the ostrich, the great spotted kiwi, the emu and the southern cassowary and two tinamous from different genera. The non-stationary base composition in these sequences violates the assumptions of most tree-building methods. When this bias is corrected using neighbour-joining with log-determinant distances and non-homogeneous maximum likelihood, the ratites are found to be monophlyletic, with moas basal, as in morphological trees. The avian sequences also violate a molecular clock, so we applied a non-parametric rate smoothing algorithm, which minimizes ancestor-descendant local rate changes, to date nodes in the tree. Using this method, most of the major ratite lineages fit the vicariance biogeography hypothesis, the exceptions being the ostrich and the kiwi, which require dispersal to explain their present distribution.

A novel mitochondrial gene order in the crinoid echinoderm Florometra serratissima

The complete nucleotide sequence of the mitochondrial genome of the crinoid Florometra serratissima has been determined. It is a circular DNA molecule, 16,005 bp in length, containing the genes for 13 proteins, small and large ribosomal RNAs, and 22 transfer RNAs (tRNAs). Three regions of unassigned sequence (UAS) greater than 73 bp have been located. The largest, UAS I, is 432 bp long and exhibits sequence similarity to the putative mitochondrial control regions seen in other animals. UAS II (77 bp) and UAS III (73 bp) are located between the 5' ends of coding sequences and may play roles as bidirectional promoters. Analyses of nucleotide composition revealed that the major peptide-encoding strand is high in T and low in C. This bias is reflected in a specific pattern of codon usage. Molecular phylogenetic analyses based on cytochrome c oxidase (COI, COII, and COIII) amino acid and nucleotide sequences did not resolve all the relationships between echinoderm classes. The overall animal mitochondrial gene content has been maintained in the crinoid, but there is extensive rearrangement with respect to both the echinoid and the asteroid mtDNA gene maps. Florometra serratissima has a novel genome organization in a segment containing most of the tRNA genes, large and small rRNA genes, and the NADH dehydrogenase subunit 1 and 2 genes. Potential pathways and mechanisms for gene rearrangements between mitochondrial gene maps of echinoderm classes and vertebrates are discussed as indicators of early deuterostome phylogeny.