The complete mitochondrial genome of the leafminer Liriomyza sativae (Diptera: Agromyzidae): great difference in the A+T-rich region compared to Liriomyza trifolii

The complete mitochondrial genome sequence of Liriomyza sativae Blanchard (15,551 bp) was determined and analyzed in this study. The circular genome contained 37 genes including 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes and an A + T-rich region. The initiation codons of COI and ND1 were ‘ATCA’ and ‘GTG’, respectively. ND2 gene used the truncated termination codon ‘T’. All the tRNA genes had the typical cloverleaf secondary structures except for tRNA^Ser(AGN) gene, which was found with the absence of a DHU arm. In addition, a tRNA-like secondary structure (tRNA^Met) was found in the A + T-rich region. The great difference was that the length of L. sativae A + T-rich region was 597 bp shorter than that of Liriomyza trifolii (Burgess). Meanwhile, some minor differences such as ‘TATA’ block were also observed in L. sativae in contrast to ‘TACA’ block in L. trifolii. There were also some essential structure elements such as ‘TATA’ block, ‘G(A)_nT’ block, poly-T stretch and stem-and-loop structure in the A + T-rich region of L. sativae mitochondrial genome.     

The complete mitochondrial genome sequence and gene organization of the mud crab (Scylla paramamosain) with phylogenetic consideration

The complete mitochondrial genome is of great importance for better understanding the genome-level characteristics and phylogenetic relationships among related species. In the present study, we determined the complete mitochondrial genome DNA sequence of the mud crab (Scylla paramamosain) by 454 deep sequencing and Sanger sequencing approaches. The complete genome DNA was 15,824 bp in length and contained a typical set of 13 protein-coding genes, 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes and a putative control region (CR). Of 37 genes, twenty-three were encoded by the heavy strand (H-strand), while the other ones were encoded by light strand (L-strand). The gene order in the mitochondrial genome was largely identical to those obtained in most arthropods, although the relative position of gene tRNA^His differed from other arthropods. Among 13 protein-coding genes, three (ATPase subunit 6 (ATP6), NADH dehydrogenase subunits 1 (ND1) and ND3) started with a rare start codon ATT, whereas, one gene cytochrome c oxidase subunit I (COI) ended with the incomplete stop codon TA. All 22 tRNAs could fold into a typical clover-leaf secondary structure, with the gene sizes ranging from 63 to 73 bp. The phylogenetic analysis based on 12 concatenated protein-coding genes showed that the molecular genetic relationship of 19 species of 11 genera was identical to the traditional taxonomy.     

The characterization of complete mitochondrial genome and phylogenetic relationship within Rattus genus (Rodentia: Muridae)

The genus Rattus is one of the main pest genus of rodent. Most species of the genus carry all kinds of pathogenic bacteria to human being. They are traditionally considered to be a least understood group. The complete mitochondrial genome of the White-Footed Indochinese Rat, Rattus nitidus was determined in this study. The characterization of mitochondrial genomes of Rattus genus was also analyzed based on comprehensive comparison. The result of evolutionary patterns of protein-coding genes (PCGs) suggested purifying selection was the predominant evolutionary forces in the mitochondrial genomes of Rattus genus. The NADH dehydrogenase 4 gene (ND4) showed a highly elevated Ka/Ks ratio compared to the other protein-coding genes, which indicated ND4 was most likely under relaxed selection pressure. Phylogenetic analysis provided a well-supported outline of Rattus genus, and revealed two groups in the genus. R. nitidus had a sister relationship with R. norvegicus.     

Complete mitochondrial genomes from four subspecies of common chaffinch (Fringilla coelebs): New inferences about mitochondrial rate heterogeneity,neutral theory,and phylogenetic relationships within the order Passeriformes

We describe whole mitochondrial genome sequences from four subspecies of the common chaffinch (Fringilla coelebs), and compare them to 31 publicly available mitochondrial genome sequences from other Passeriformes. Rates and patterns of mitochondrial gene evolution are analyzed at different taxonomic levels within this avian order, and evidence is adduced for and against the nearly neutral theory of molecular evolution and the role of positive selection in shaping genetic variation of this small but critical genome. We find evidence of mitochondrial rate heterogeneity in birds as in other vertebrates, likely due to differences in mutational pressure across the genome. Unlike in gadine fish and some of the human mitochondrial work we do not observe strong support for the nearly neutral theory of molecular evolution; instead evidence from molecular clocks, distribution of dN/dS ratios at different levels of the taxonomic hierarchy and in different lineages, McDonald–Kreitman tests within Fringillidae, and site-specific tests of selection within Passeriformes, all point to a role for positive selection, especially for the complex I NADH dehydrogenase genes. The protein-coding mitogenome phylogeny of the order Passeriformes is broadly consistent with previously-reported molecular findings, but provides support for a sister relationship between the superfamilies Muscicapoidea and Passeroidea on a short basal internode of the Passerida where relationships have been difficult to resolve. An unexpected placement of the Paridae (represented by Hume's groundpecker) within the Muscicapoidea was observed. Consistent with other molecular studies the mtDNA phylogeny reveals paraphyly within the Muscicapoidea and a sister relationship of Fringilla with Carduelis rather than Emberiza.     

Active transcription of the pseudogene for subunit 7 of the NADH dehydrogenase in Marchantia polymorpha mitochondria

A pseudogene, ψnad7, which has significant sequence similarity (66.7% amino acid identity) with the bovine nuclear gene for a 49 kDa subunit of the NADH dehydrogenase (NADH:ubiquinone oxidoreductase, EC 1.6.99.3), has been identified on the mitochondrial genome of the liverwort Marchantia polymorpha. The predicted coding region, which includes six termination codons, is actively transcribed into RNA molecules of 16 and 9.6 kb in length, but RNA splicing products were not detected in the liverwort mitochondria. Genomic DNA blot analysis and RNA blot analysis using poly(A)⁺ RNA suggest that a structurally related nuclear gene encodes the mitochondrial ND7 polypeptide. These results imply that this ψnad7 is a relic of a gene transfer event from the mitochondrial genome into the nuclear genome during mitochondrial evolution in M. polymorpha.     

Cloning and sequence analysis of a mitochondrial gene cluster encoding two NADH dehydrogenase subunits and seven tRNA molecules from Trichoderma reesei

A mitochondrial gene cluster, encoding proteins homologous to NADH dehydrogenase subunits II and III (ND2 and ND3) and seven tRNAs, from Trichoderma reesei QM9414 was cloned and sequenced. These genes are clustered tandemly on the mitochondrial genome of QM9414. Phylogenetic analysis showed that ND2 and ND3 were most closely related to the mitochondrial ND subunits II (71% identity) and III (70% identity) from Podospora anserine. Northern dot blot analysis showed that the nd2 and nd3 genes are actively transcribed in the T. reesei mitochondria.     

Conspicuous males suffer higher predation risk: visual modelling and experimental evidence from lizards

Colour pattern variation is a striking and widespread phenomenon. Differential predation risk between individuals is often invoked to explain colour variation, but empirical support for this hypothesis is equivocal. We investigated differential conspicuousness and predation risk in two species of Australian rock dragons, Ctenophorus decresii andC. vadnappa . To humans, the coloration of males of these species varies between ‘bright’ and ‘dull’. Visual modelling based on objective colour measurements and the spectral sensitivities of avian visual pigments showed that dragon colour variants are differentially conspicuous to the visual system of avian predators when viewed against the natural background. We conducted field experiments to test for differential predation risk, using plaster models of ‘bright’ and ‘dull’ males. ‘Bright’ models were attacked significantly more often than ‘dull’ models suggesting that differential conspicuousness translates to differential predation risk in the wild. We also examined the influence of natural geographical range on predation risk. Results from 22 localities suggest that predation rates vary according to whether predators are familiar with the prey species. This study is among the first to demonstrate both differential conspicuousness and differential predation risk in the wild using an experimental protocol. Copyright 2003 Published by Elsevier Ltd on behalf of The Association for the Study of Animal Behaviour.      

Characterization of the complete mitochondrial genome and a set of polymorphic microsatellite markers through next-generation sequencing for the brown brocket deer Mazama gouazoubira

The complete mitochondrial genome of the brown brocket deer Mazama gouazoubira and a set of polymorphic microsatellite markers were identified by 454-pyrosequencing. De novo genome assembly recovered 98% of the mitochondrial genome with a mean coverage of 9-fold. The mitogenome consisted of 16,356 base pairs that included 13 protein-coding genes, two ribosomal subunit genes, 22 transfer RNAs and the control region, as found in other deer. The genetic divergence between the mitogenome described here and a previously published report was ∼0.5%, with the control region and ND5 gene showing the highest intraspecific variation. Seven polymorphic loci were characterized using 15 unrelated individuals; there was moderate genetic variation across most loci (mean of 5.6 alleles/locus, mean expected heterozygosity = 0.70), with only one locus deviating significantly from Hardy-Weinberg equilibrium, probably because of null alleles. Marker independence was confirmed with tests for linkage disequilibrium. The genetic variation of the mitogenome and characterization of microsatellite markers will provide useful tools for assessing the phylogeography and population genetic patterns in M. gouazoubira, particularly in the context of habitat fragmentation in South America.     

Variation in Coding (NADH Dehydrogenase Subunits 2, 3, and 6) and Noncoding Intergenic Spacer Regions of the Mitochondrial Genome in Octocorallia (Cnidaria: Anthozoa)

Low rates of evolution in cnidarian mitochondrial genes such as COI and 16S rDNA have hindered molecular systematic studies in this important invertebrate group. We sequenced fragments of 3 mitochondrial protein-coding genes (NADH dehydrogenase subunits ND2, ND3 and ND6) as well as the COI-COII intergenic spacer, the longest noncoding region found in the octocoral mitochondrial genome, to determine if any of these regions contain levels of variation sufficient for reconstruction of phylogenetic relationships among genera of the anthozoan subclass Octocorallia. Within and between the soft coral families Alcyoniidae and Xeniidae, sequence divergence in the genes ND2 (539 bp), ND3 (102 bp), and ND6 (444 bp) ranged from 0.5% to 12%, with the greatest pairwise distances between the 2 families. The COI-COII intergenic spacer varied in length from 106 to 122 bp, and pairwise sequence divergence values ranged from 0% to 20.4%. Phylogenetic trees constructed using each region separately were poorly resolved. Better phylogenetic resolution was obtained in a combined analysis using all 3 protein-coding regions (1085 bp total). Although relationships among some pairs of species and genera were well supported in the combined analysis, the base of the alcyoniid family tree remained an unresolved polytomy. We conclude that variation in the NADH subunit coding regions is adequate to resolve phylogenetic relationships among families and some genera of Octocorallia, but insufficient for most species - or population-level studies. Although the COI-COII intergenic spacer exhibits greater variability than the protein-coding regions and may contain useful species-specific markers, its short length limits its phylogenetic utility.     

Positive and purifying selection in mitochondrial genomes of a bird with mitonuclear discordance

Diversifying selection on metabolic pathways can reduce intraspecific gene flow and promote population divergence. An opportunity to explore this arises from mitonuclear discordance observed in an Australian bird Eopsaltria australis. Across >1500 km, nuclear differentiation is low and latitudinally structured by isolation by distance, whereas two highly divergent, parapatric mitochondrial lineages (>6.6% in ND2) show a discordant longitudinal geographic pattern and experience different climates. Vicariance, incomplete lineage sorting and sex‐biased dispersal were shown earlier to be unlikely drivers of the mitonuclear discordance; instead, natural selection on a female‐linked trait was the preferred hypothesis. Accordingly, here we tested for signals of positive, divergent selection on mitochondrial genes in E. australis. We used codon models and physicochemical profiles of amino acid replacements to analyse complete mitochondrial genomes of the two mitochondrial lineages in E. australis, its sister species Eopsaltria griseogularis, and outgroups. We found evidence of positive selection on at least five amino acids, encoded by genes of two oxidative phosphorylation pathway complexes NADH dehydrogenase (ND4 and ND4L) and cytochrome bc₁ (cyt‐b) against a background of widespread purifying selection on all mitochondrial genes. Three of these amino acid replacements were fixed in ND4 of the geographically most widespread E. australis lineage. The other two replacements were fixed in ND4L and cyt‐b of the geographically more restricted E. australis lineage. We discuss whether this selection may reflect local environmental adaptation, a by‐product of other selective processes, or genetic incompatibilities, and propose how these hypotheses can be tested in future.     

Complete mitochondrial genome of <Emphasis Type="Italic">Ampittia dioscorides</Emphasis> (Lepidoptera: Hesperiidae) and its phylogenetic analysis

The complete mitochondrial genome of Ampittia dioscorides (Lepidoptera: Hesperiidae) was determined. The sequenced genome is a circular molecule of 15313 bp, containing 13 protein-coding genes (PCGs), 22 tRNA genes, two rRNA genes, and an A + T-rich region. The gene arrangements and transcribing directions are identical to those in most of the reported lepidopteran mitogenomes. The base composition of the whole genome and genes or regions are also similar to those in other lepidopteran species. All the PCGs are initiated by typical ATN codons; the exception being COI, which begins with a CGA codon. Eight genes (ND2, ATPase8, ATPase6, COIII, ND5, ND4L, ND6, and Cytb) end with a TAA stop codon, and two genes (ND1 and ND3) end with TAG. The remaining three genes (COI and COII, which end with TA-, and ND4, which ends with T-) have incomplete stop codons. All tRNAs have the typical clover-leaf structure of mitochondrial tRNAs, with the exception of tRNA^Ser(AGY). On the basis of the concatenated nucleotide and amino acid sequences of the 13 PCGs and wingless gene of 22 butterfly species, maximum parsimony (MP) and Bayesian inference (BI) trees were constructed, respectively. Both MP and BI trees had the same topological structure: ((((Nymphalidae + Danaidae) + Lycaenidae) + Pieridae) + Papilionidae) + Hesperiidae). The results provide support for Hesperiidae as a superfamily-level taxon.     

Complete mitochondrial genome of the <Emphasis Type="Italic">Pleuronichthys lighti</Emphasis> (Pleuronectiformes,Pleuronectidae) with phylogenetic consideration

Mitochondrial genome has been used to shed light on many fields of both basic and applied research, including the study of molecular evolution. The complete mitochondrial genome sequence of 17368 bp nucleotides from the Pleuronichthys lighti was determined. It was a circular double-stranded DNA molecule with identical set of 22 transfer RNA genes, 2 ribosomal RNA genes, 13 protein-coding genes as well as a non-coding control region. Stand asymmetry in the nucleotide composition was reflected in the codon usage of genes oriented in opposite directions. In the control region, we identified the extended termination associated sequence domain, the central conserved sequence block domain and the conserved sequence block domain, and two complete repeat region. They were “TTACAATA” and “TGTTGTAA”, respectively. All known 12 mitochondrial genomes of Pleuronectinae fishes were downloaded and analyzed; there were 5570 variable sites in the consensus sequences of 15241 base pairs, calculation of total sites were 35.5%. The highest sequence divergence was 50% (ATP8) and the Kimura-2-parameter genetic distance was 0.235 (ND6), whereas the COIII had the lowest sequence divergence (28.8%) and genetic distance (0.128); the protein coding genes were mainly acted by purifying selection which was detected by selection tests. Analysis of confidence and the information content for per nucleotide revealed ND5, ATP6, COI and ND4 genes were suitable molecular markers for phylogenetic study of Pleuronectinae fishes. Phylogenetic analysis using Bayesian computational algorithms based on COI genes provided support for the taxonomic status of P. lighti, which was consistent with the traditional taxonomy.     

A Novel Mitochondrial Heteroplasmic C13806A Point Mutation Associated with Iranian Friedreich’s Ataxia

Friedreich’s ataxia (FRDA) is an autosomal recessive neurodegenerative disorder caused by decreased expression of the protein Frataxin. Frataxin deficiency leads to excessive free radical production and dysfunction of chain complexes. Mitochondrial DNA (mtDNA) could be considered a candidate modifier factor for FRDA disease, since mitochondrial oxidative stress is thought to be involved in the pathogenesis of this disease. It prompted us to focus on the mtDNA and monitor the nucleotide changes of genome which are probably the cause of respiratory chain defects and reduced ATP generation. We searched about 46% of the entire mitochondrial genome by temporal temperature gradient gel electrophoresis (TTGE) and DNA fragments showing abnormal banding patterns were sequenced for the identification of exact mutations. In 18 patients, for the first time, we detected 26 mtDNA mutations; of which 5 (19.2%) was novel and 21 (80.8%) have been reported in other diseases. Heteroplasmic C13806A polymorphisms were associated with Iranian FRDA patients (55.5%). Our results showed that NADH dehydrogenase (ND) genes mutations in FRDA samples were higher than normal controls (P < 0.001) and we found statistically significant inverse correlation (r = −0.8) between number of mutation in ND genes and age of onset in FRDA patients. It is possible that mutations in ND genes could constitute a predisposing factor which in combination with environmental risk factors affects age of onset and disease progression.     

The influence of phylogeny and life history on telomere lengths and telomere rate of change among bird species: A meta‐analysis

Longevity is highly variable among animal species and has coevolved with other life‐history traits, such as body size and rates of reproduction. Telomeres, through their erosion over time, are one of the cell mechanisms that produce senescence at the cell level and might even have an influence on the rate of aging in whole organisms. However, uneroded telomeres are also risk factors of cell immortalization. The associations of telomere lengths, their rate of change, and life‐history traits independent of body size are largely underexplored for birds. To test associations of life‐history traits and telomere dynamics, we conducted a phylogenetic meta‐analysis using studies of 53 species of birds. We restricted analyses to studies that applied the telomere restriction fragment length (TRF) method, and examined relationships between mean telomere length at the chick (Chick TL) and adult (Adult TL) stages, the mean rate of change in telomere length during life (TROC), and life‐history traits. We examined 3 principal components of 12 life‐history variables that represented: body size (PC1), the slow–fast continuum of pace of life (PC2), and postfledging parental care (PC3). Phylogeny had at best a small‐to‐medium influence on Adult and Chick TL (r ² = .190 and .138, respectively), but a substantial influence on TROC (r ² = .688). Phylogeny strongly influenced life histories: PC1 (r ² = .828), PC2 (.838), and PC3 (.613). Adult TL and Chick TL were poorly associated with the life‐history variables. TROC, however, was negatively and moderate‐to‐strongly associated with PC2 (unadjusted r = −.340; with phylogenetic correction, r = −.490). Independent of body size, long‐lived species with smaller clutches, and slower embryonic rate of growth may exhibit less change in telomere length over their lifetimes. We suggest that telomere lengths may have diverged, even among closely avian‐related species, yet telomere dynamics are strongly linked to the pace of life.     

Molecular characterization and expression analysis of ubiquitin-activating enzyme E1 gene in Citrus reticulata

Ubiquitin-activating enzyme E1 (UBE1) catalyzes the first step in the ubiquitination reaction, which targets a protein for degradation via a proteasome pathway. UBE1 plays an important role in metabolic processes. In this study, full-length cDNA and DNA sequences of UBE1 gene, designated CrUBE1, were obtained from ‘Wuzishatangju’ (self-incompatible, SI) and ‘Shatangju’ (self-compatible, SC) mandarins. 5 amino acids and 8 bases were different in cDNA and DNA sequences of CrUBE1 between ‘Wuzishatangju’ and ‘Shatangju’, respectively. Southern blot analysis showed that there existed only one copy of the CrUBE1 gene in genome of ‘Wuzishatangju’ and ‘Shatangju’. The temporal and spatial expression characteristics of the CrUBE1 gene were investigated using semi-quantitative RT-PCR (SqPCR) and quantitative real-time PCR (qPCR). The expression level of the CrUBE1 gene in anthers of ‘Shatangju’ was approximately 10-fold higher than in anthers of ‘Wuzishatangju’. The highest expression level of CrUBE1 was detected in pistils at 7 days after self-pollination of ‘Wuzishatangju’, which was approximately 5-fold higher than at 0 h. To obtain CrUBE1 protein, the full-length cDNA of CrUBE1 genes from ‘Wuzishatangju’ and ‘Shatangju’ were successfully expressed in Pichia pastoris. Pollen germination frequency of ‘Wuzishatangju’ was significantly inhibited with increasing of CrUBE1 protein concentrations from ‘Wuzishatangju’.     

Identification, distribution and molecular evolution of the pacifastin gene family in Metazoa

Background

This study was motivated by the observation of unusual mitochondrial haplotype distributions and associated physiological differences between populations of the killifish Fundulus heteroclitus distributed along the Atlantic coast of North America. A distinct "northern" haplotype is fixed in all populations north of New Jersey, and does not appear south of New Jersey except in extreme upper-estuary fresh water habitats, and northern individuals are known to be more tolerant of hyposmotic conditions than southern individuals. Complete mitochondrial genomes were sequenced from individuals from northern coastal, southern coastal, and fresh water populations (and from out-groups). Comparative genomics approaches were used to test multiple evolutionary hypotheses proposed to explain among-population genome variation including directional selection and hybridization.

Results

Structure and organization of the Fundulus mitochondrial genome is typical of animals, yet subtle differences in substitution patterns exist among populations. No signals of directional selection or hybridization were detected. Mitochondrial genes evolve at variable rates, but all genes exhibit very low dN/dS ratios across all lineages, and the southern population harbors more synonymous polymorphism than other populations.

Conclusion

Evolution of mitochondrial genomes within Fundulus is primarily governed by interaction between strong purifying selection and demographic influences, including larger historical population size in the south. Though directional selection and hybridization hypotheses were not supported, adaptive processes may indirectly contribute to partitioning of variation between populations.     

Phylogenetic analyses of the harvest mouse,Micromys minutus (Rodentia: Muridae) based on the complete mitogenome sequences

Harvest Mouse (Micromys minutus) has a very wide range of distribution in Asia and Europe. However, the phylogenetic relationship of M. minutus is still uncertain. In this study, we determined the complete mitochondrial (mt) genome sequences of M. minutus, and used the complete mitochondrial genome sequences constructed the phylogenetic tree of Muroidea. The size of the genome is 16,232 bp in length and has a base composition of 33.6% A, 29.1% T, 24.8% C, and 12.5% G. The mitogenome structure was similar to that of typical vertebrate and other rodents' mitochondrial genomes, includes 13 protein-coding genes, 2 rRNA genes (12S rRNA and 16S rRNA), 22 tRNA genes, and 1 control region. We suggested a new initiation codon for ND5 (NADH dehydrogenase subunit), which has been never reported in the mitochondrial genome of vertebrate. The ML and BI phylogenetic trees, which based on the combination of the 12 protein-coding genes, supported strongly that the genus Micromys was represent an early offshoot within the Muridae with high support values (BI = 1.00, ML = 100).     

Genome‐wide single‐generation signatures of local selection in the panmictic European eel

Next‐generation sequencing and the collection of genome‐wide data allow identifying adaptive variation and footprints of directional selection. Using a large SNP data set from 259 RAD‐sequenced European eel individuals (glass eels) from eight locations between 34 and 64^oN, we examined the patterns of genome‐wide genetic diversity across locations. We tested for local selection by searching for increased population differentiation using F_ST‐based outlier tests and by testing for significant associations between allele frequencies and environmental variables. The overall low genetic differentiation found (F_ST = 0.0007) indicates that most of the genome is homogenized by gene flow, providing further evidence for genomic panmixia in the European eel. The lack of genetic substructuring was consistent at both nuclear and mitochondrial SNPs. Using an extensive number of diagnostic SNPs, results showed a low occurrence of hybrids between European and American eel, mainly limited to Iceland (5.9%), although individuals with signatures of introgression several generations back in time were found in mainland Europe. Despite panmixia, a small set of SNPs showed high genetic differentiation consistent with single‐generation signatures of spatially varying selection acting on glass eels. After screening 50 354 SNPs, a total of 754 potentially locally selected SNPs were identified. Candidate genes for local selection constituted a wide array of functions, including calcium signalling, neuroactive ligand–receptor interaction and circadian rhythm. Remarkably, one of the candidate genes identified is PERIOD, possibly related to differences in local photoperiod associated with the >30° difference in latitude between locations. Genes under selection were spread across the genome, and there were no large regions of increased differentiation as expected when selection occurs within just a single generation due to panmixia. This supports the conclusion that most of the genome is homogenized by gene flow that removes any effects of diversifying selection from each new generation.