Complete mitochondrial genomes of the human follicle mites Demodex brevis and D. folliculorum: novel gene arrangement,truncated tRNA genes,and ancient divergence between species

Background

Follicle mites of the genus Demodex are found on a wide diversity of mammals, including humans; surprisingly little is known, however, about the evolution of this association. Additional sequence information promises to facilitate studies of Demodex variation within and between host species. Here we report the complete mitochondrial genome sequences of two species of Demodex known to live on humans—Demodex brevis and D. folliculorum—which are the first such genomes available for any member of the genus. We analyzed these sequences to gain insight into the evolution of mitochondrial genomes within the Acariformes. We also used relaxed molecular clock analyses, based on alignments of mitochondrial proteins, to estimate the time of divergence between these two species.

Results

Both Demodex genomes shared a novel gene order that differs substantially from the ancestral chelicerate pattern, with transfer RNA (tRNA) genes apparently having moved much more often than other genes. Mitochondrial tRNA genes of both species were unusually short, with most of them unable to encode tRNAs that could fold into the canonical cloverleaf structure; indeed, several examples lacked both D- and T-arms. Finally, the high level of sequence divergence observed between these species suggests that these two lineages last shared a common ancestor no more recently than about 87 mya.

Conclusions

Among Acariformes, rearrangements involving tRNA genes tend to occur much more often than those involving other genes. The truncated tRNA genes observed in both Demodex species would seem to require the evolution of extensive tRNA editing capabilities and/or coevolved interacting factors. The molecular machinery necessary for these unusual tRNAs to function might provide an avenue for developing treatments of skin disorders caused by Demodex. The deep divergence time estimated between these two species sets a lower bound on the time that Demodex have been coevolving with their mammalian hosts, and supports the hypothesis that there was an early split within the genus Demodex into species that dwell in different skin microhabitats.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1124) contains supplementary material, which is available to authorized users.     

There appear to be conserved constraints on the distribution of nucleotide sequences in cellular genomes

Summary The data from a genomic library can be sorted into the frequencies of every possible tetranucleotide in the sequence. This tabulation, a short sequence distribution, contains the frequency of occurrence of the 256 tetranucleotides and thus seems to serve as a vehicle for averaging sequence information. Two such distributions can be readily compared by correlation. Reported here are correlations (Spearmanr _s) of the distributions from all of the genomic libraries in GenBank 44.0 with sizes equal to or larger than that ofSalmonella typhimurium, except for the data for mouse and humans. All of the organisms examined showed highly significant correlations between the two DNA strands (not the complementarity expected from base pairing). Of 155 comparisons between libraries, 132 showed significant correlations at the 99% confidence level. Application of the correlation coefficients as a similarity matrix clustered most organisms in a phenogram in a pattern consistent with other hypotheses. This suggests a highly conserved pattern underlying all other genetic information in cellular DNA and affecting both DNA strands, perhaps caused by interaction with conserved factors necessary for DNA packaging.     

The chloroplast genomes of conifers lack one of the rRNA-encoding inverted repeats

Summary Chloroplast DNA from species of five different conifer genera was extracted and studied by Southern blot analysis. For all these species, hybridization with heterologous probes specific for 16 S and 23 S rDNA detected only one chloroplast DNA fragment per enzyme digest. This observation suggests that the 16 S and 23 S rRNA genes are not duplicated in these genomes. The unique 16 S rDNA-containing BamHI fragment from Pinus contorta Dougl. was cloned and restriction mapped. Apart from the 16 S rRNA gene, this fragment also contained the psbC and psbD genes. It is concluded that the chloroplast genomes of a wide taxonomic range of conifers lack one of the inverted repeat elements and that a dislocation of the psbDC gene cluster has occurred in P. contorta.     

The serine repeat antigen (SERA) gene family phylogeny in Plasmodium: the impact of GC content and reconciliation of gene and species trees

Plasmodium falciparum is the parasite responsible for the most acute form of malaria in humans. Recently, the serine repeat antigen (SERA) in P. falciparum has attracted attention as a potential vaccine and drug target, and it has been shown to be a member of a large gene family. To clarify the relationships among the numerous P. falciparum SERAs and to identify orthologs to SERA5 and SERA6 in Plasmodium species affecting rodents, gene trees were inferred from nucleotide and amino acid sequence data for 33 putative SERA homologs in seven different species. (A distance method for nucleotide sequences that is specifically designed to accommodate differing GC content yielded results that were largely compatible with the amino acid tree. Standard-distance and maximum-likelihood methods for nucleotide sequences, on the other hand, yielded gene trees that differed in important respects.) To infer the pattern of duplication, speciation, and gene loss events in the SERA gene family history, the resulting gene trees were then "reconciled" with two competing Plasmodium species tree topologies that have been identified by previous phylogenetic studies. Parsimony of reconciliation was used as a criterion for selecting a gene tree/species tree pair and provided (1) support for one of the two species trees and for the core topology of the amino acid-derived gene tree, (2) a basis for critiquing fine detail in a poorly resolved region of the gene tree, (3) a set of predicted "missing genes" in some species, (4) clarification of the relationship among the P. falciparum SERA, and (5) some information about SERA5 and SERA6 orthologs in the rodent malaria parasites. Parsimony of reconciliation and a second criterion--implied mutational pattern at two key active sites in the SERA proteins-were also seen to be useful supplements to standard "bootstrap" analysis for inferred topologies.     

Organization and evolution of the mitochondrial genome of yeast

Summary The mitochondrial genome of yeast (S. cerevisiae orS. carlsbergensis) appears to be formed by 60–70 genetic units, each one of which is formed by (1) a GC-rich sequence, possibly having a regulatory role; (2) a gene, and (3) an AT-rich spacer, which probably is not transcribed. Recombination in this genome appears to underlie a number of important phenomena. The organization of the mitochondrial genome of yeast and these recombinational events are discussed in relationship with the organization and evolution of the nuclear genome of eukaryotes.     

Protein domain analysis in the era of complete genomes

Domains present one of the most useful levels at which to understand protein function, and domain family-based analysis has had a profound impact on the study of individual proteins. Protein domain discovery has been progressing steadily over the past 30 years. What are the realistically achievable goals of sequence-based domain analysis, and how far off are they for the sequences encoded in eukaryotic genomes? Here we address some of the issues involved in better coverage of sequence-based domain annotation, and the integration of these results within the wider context of genomes, structures and function.     

Migration of Apicomplexa across biological barriers: the Toxoplasma and Plasmodium rides

The invasive stages of Apicomplexa parasites, called zoites, have been largely studied in in vitro systems, with a special emphasis on their unique gliding and host cell invasive capacities. In contrast, the means by which these parasites reach their destination in their hosts are still poorly understood. We summarize here our current understanding of the cellular basis of in vivo parasitism by two well-studied Apicomplexa zoites, the Toxoplasma tachyzoite and the Plasmodium sporozoite. Despite being close relatives, these two zoites use different strategies to reach their goal and establish infection.     

A highly conserved microsatellite in the dystrophin gene of diverse mammalian species

The presence of a CA repeat within the 3'-untranslated region (UTR) of the dystrophin gene has been reported previously in several species. Because microsatellites showing high cross-species homology can be conveniently used as markers in those species for which detailed linkage maps have not yet been developed, we evaluated whether the CA repeat could be amplified from a wide variety of mammalian species. Using a single pair of canine-specific oligonucleotide primers, we successfully amplified the 3'-UTR from 18 different carnivore and six additional species (human, chimpanzee, goat, cow, rabbit and mouse) and show conservation of the CA repeat in the dystrophin gene from a wide range of evolutionarily diverse mammalian species.     

A highly polymorphic region in the X chromosome of Drosophila

Summary Many cloned regions of the Drosophila genome show minimal variation between strains in overall sequence arrangement. While restriction site polymorphisms occur and the location of transposable elements may vary from one strain to another, such changes appear to be relatively minor variations, superimposed on overall genome stability. In contrast to this general situation, we describe here a segment of the X chromosome that is highly polymorphic in four strains of D. melanogaster and in D. simulans. The strains differ in the presence and extent of a short duplication and the presence of repetitive DNA. These results suggest that different regions of the genome may be subject to different evolutionary constraints, with some regions being particularly prone to extensive changes, even within a single species.     

Achlya mitochondrial DNA: gene localization and analysis of inverted repeats

Summary Mitochondrial DNA from four strains of the oomycete Achlya has been compared and nine gene loci mapped, including that of the ribosomal protein gene, var1. Examination of the restriction enzyme site maps showed the presence of four insertions relative to a map common to all four strains. All the insertions were found in close proximity to genic regions. The four strains also cotained the inverted repeat first observed in A. ambisexualis (Hudspeth et al. 1983), allowing an examination by analysis of retained restriction sites of the evolutionary stability of repeated DNA sequences relative to single copy sequences. Although the inverted repeat is significantly more stable than single copy sequences, more detailed analysis indicated that this stability is limited to the portion encoding the ribosomal RNA genes. Thus, the apparent evolutionary stability of the repeat does not appear to derive from the inverted repeat structure per se.Abbreviations ATPase 6, 9 genes for ATPase subunits 6 and 9 - COI, II, III genes for cytochrome oxidase subunits 1, 2, and 3 - COB gene for apocytochrome b - L-, S-RNA genes for the mitochondrial large and small ribosomal RNAs - mtDNA mitochondrial DNA - var1 gene for the S. cerevisiae mitochondrially, encoded ribosomal protein - m.u. map units - bp base pairs - kb kilobase pairs     

The mitochondrial DNA of Dictyostelium discoideum: complete sequence, gene content and genome organization

We present an overview of the gene content and organization of the mitochondrial genome of Dictyostelium discoideum. The mitochondria genome consists of 55,564 bp with an A + T content of 72.6%. The identified genes include those for two ribosomal RNAs (rnl and rns), 18 tRNAs, ten subunits of the NADH dehydrogenase complex (nad1, 2, 3, 4, 4L, 5, 6, 7, 9 and 11), apocytochrome b (cytb), three subunits of the cytochrome oxidase (cox1/2 and 3), four subunits of the ATP synthase complex (atp1, 6, 8 and 9), 15 ribosomal proteins, and five other ORFs, excluding intronic ORFs. Notable features of D. discoideum mtDNA include the following. (1) All genes are encoded on the same strand of the DNA and a universal genetic code is used. (2) The cox1 gene has no termination codon and is fused to the downstream cox2 gene. The 13 genes for ribosomal proteins and four ORF genes form a cluster 15.4 kb long with several gene overlaps. (3) The number of tRNAs encoded in the genome is not sufficient to support the synthesis of mitochondrial protein. (4) In total, five group I introns reside in rnl and cox1/2, and three of those in cox1/2 contain four free-standing ORFs. We compare the genome to other sequenced mitochondrial genomes, particularly that of Acanthamoeba castellanii. Received: 5 July 1999 / Accepted: 17 January 2000     

A first glimpse into the pattern and scale of gene transfer in Apicomplexa

Reports of plant-like and bacterial-like genes for a number of parasitic organisms, most notably those within the Apicomplexa and Kinetoplastida, have appeared in the literature over the last few years. Among the apicomplexan organisms, following discovery of the apicomplexan plastid (apicoplast), the discovery of plant-like genes was less surprising although the extent of transfer and the relationship of transferred genes to the apicoplast remained unclear. We used new genome sequence data to begin a systematic examination of the extent and origin of transferred genes in the Apicomplexa combined with a phylogenomic approach to detect potential gene transfers in four apicomplexan genomes. We have detected genes of algal nuclear, chloroplast (cyanobacterial) and proteobacterial origin. Plant-like genes were detected in species not currently harbouring a plastid (e.g. Cryptosporidium parvum) and putatively transferred genes were detected that appear to be unrelated to the function of the apicoplast. While the mechanism of acquisition for many of the identified genes is not certain, it appears that some were most likely acquired via intracellular gene transfer from an algal endosymbiont while others may have been acquired via horizontal gene transfer.     

Sequence analysis of the Rhop-3 gene of Plasmodium yoelii

The 110 kDa/Rhop-3 rhoptry protein of Plasmodium falciparum is non-covalently associated with two other proteins, the 140 kDa Rhop-1 and the 130 kDa Rhop-2. cDNAs encoding Rhop-3 from Plasmodium yoelii were isolated using rhoptry-specific antisera from Plasmodium falciparum, P. yoelii, and Plasmodium chabaudi. The cDNAs encoded peptides with partial homology to the C-terminal region (residues 541-861) of P. falciparum Rhop-3. Core regions of homology to the P. falciparum gene will be useful in determining the biological role of Rhop-3 and its potential as a vaccine candidate for malaria.     

Non-coding RNA gene families in the genomes of anopheline mosquitoes

Background

Only a small fraction of the mosquito species of the genus Anopheles are able to transmit malaria, one of the biggest killer diseases of poverty, which is mostly prevalent in the tropics. This diversity has genetic, yet unknown, causes. In a further attempt to contribute to the elucidation of these variances, the international “Anopheles Genomes Cluster Consortium” project (a.k.a. “16 Anopheles genomes project”) was established, aiming at a comprehensive genomic analysis of several anopheline species, most of which are malaria vectors. In the frame of the international consortium carrying out this project our team studied the genes encoding families of non-coding RNAs (ncRNAs), concentrating on four classes: microRNA (miRNA), ribosomal RNA (rRNA), small nuclear RNA (snRNA), and in particular small nucleolar RNA (snoRNA) and, finally, transfer RNA (tRNA).

Results

Our analysis was carried out using, exclusively, computational approaches, and evaluating both the primary NGS reads as well as the respective genome assemblies produced by the consortium and stored in VectorBase; moreover, the results of RNAseq surveys in cases in which these were available and meaningful were also accessed in order to obtain supplementary data, as were “pre-genomic era” sequence data stored in nucleic acid databases. The investigation included the identification and analysis, in most species studied, of ncRNA genes belonging to several families, as well as the analysis of the evolutionary relations of some of those genes in cross-comparisons to other members of the genus Anopheles.

Conclusions

Our study led to the identification of members of these gene families in the majority of twenty different anopheline taxa. A set of tools for the study of the evolution and molecular biology of important disease vectors has, thus, been obtained.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1038) contains supplementary material, which is available to authorized users.     

Evolution of the deep-sea gulper eel mitochondrial genomes: large-scale gene rearrangements originated within the eels

Recent studies have demonstrated that deviations from the typical vertebrate mitochondrial gene order are more frequent than initially thought. Such deviations, however, are minor, with inversions and/or translocations of a few genes being involved and tandem duplication of the gene regions followed by deletions of genes having been invoked as mechanisms originating in such novel gene order. During the course of molecular phylogenetic studies on the Elopomorpha (eels and their allies), we found that mitochondrial genomes (mitogenomes) from the two deep-sea gulper eels, Eurypharynx pelecanoides (Eurypharyngidae) and Saccopharynx lavenbergi (Saccopharyngidae), exhibit an identical gene order which greatly differs from that of any other vertebrates. Phylogenetic analysis using the mitogenomic data from 59 species of fish not only confirmed a single origin of such a gene order with confidence but also indicated that it had been derived from the typical vertebrate gene order. Detailed comparisons of the gulper eel gene order with that of typical vertebrates suggested that occurrence of a single step, large-scale duplication of gene region extending >12 kb, followed by deletions of genes in a common ancestor of the two species, most parsimoniously accounts for this unusual gene arrangement.     

Merozoite surface protein 2 of Plasmodium falciparum: expression, structure, dynamics, and fibril formation of the conserved N-terminal domain

Merozoite surface protein 2 (MSP2) is a GPI-anchored protein on the surface of the merozoite stage of the malaria parasite Plasmodium falciparum. It is largely disordered in solution, but has a propensity to form amyloid-like fibrils under physiological conditions. The N-terminal conserved region (MSP2(1-25)) is part of the protease-resistant core of these fibrils. To investigate the structure and dynamics of this region, its ability to form fibrils, and the role of individual residues in these properties, we have developed a bacterial expression system that yields > or =10 mg of unlabeled or (15)N-labeled peptide per litre of culture. Two recombinant versions of MSP2(1-25), wild-type and a Y7A/Y16A mutant, have been produced. Detailed conformational analysis of the wild-type peptide and backbone (15)N relaxation data indicated that it contains beta-turn and nascent helical structures in the central and C-terminal regions. Residues 6-21 represent the most ordered region of the structure, although there is some flexibility around residues 8 and 9. The 10-residue sequence (MSP2(7-16)) (with two Tyr residues) was predicted to have a higher propensity for beta-aggregation than the 8-mer sequence (MSP2(8-15)), but there was no significant difference in conformation between MSP2(1-25) and [Y7A,Y16A]MSP2(1-25) and the rate of fibril formation was only slightly slower in the mutant. The peptide expression system described here will facilitate further mutational analyses to define the roles of individual residues in transient structural elements and fibril formation, and thus contribute to the further development of MSP2 as a malaria vaccine candidate.     

基于线粒体基因cyt b的鹟亚科部分鸟类系统发育

采用分子系统学方法对鹟亚科(Muscicapinae)6属31种鸟类的cytb基因序列992bp进行系统发生分析。以荒漠伯劳(Lanius isabellinus)和发冠卷尾(Dicrurus hottentottus)为外群,采用贝叶斯法(Bayesian,BI)、最大似然法(Maximum-likelihood,ML)和最大简约法(Maximumparsimony,MP)分别构建鹟亚科的系统发育树。结果支持:寿带属(Terpsiphone)、扇尾鹟属(Rhipidura)与方尾鹟属(Culicicapa)可从鹟亚科中移出,其中寿带属归入王鹟科(Monarchidae),扇尾鹟属与方尾属归入扇尾鹟科(Rhipiduridae);鹟属(Muscicapa)、仙鹟属(Niltava)为单系发生,并聚为姐妹群,亲缘关系较近;姬鹟属(Ficedula)并非单系发生,白眉姬鹟(Ficedulazanthopygia)在3种系统发生树中的位置差别较大,研究结果未能确定其分类地位;铜蓝(Muscicapa thalassina)与白腹蓝(Cyanoptila cyanomelana)亲缘关系较近,前者应从属中移出,后者应从姬属移出,共同归入仙属或列为仙属的姐妹属。上述结论解决了亚科部分有争议属、种间的进化关系,为亚科分类系统提供了DNA水平证据。     

Antioxidant properties of UCP1 are evolutionarily conserved in mammals and buffer mitochondrial reactive oxygen species

Mitochondrial uncoupling reduces reactive oxygen species (ROS) production and appears to be important for cellular signaling/protection, making it a focus for the treatment of metabolic and age-related diseases. Whereas the physiological role of uncoupling protein 1 (UCP1) of brown adipose tissue is established for thermogenesis, the function of UCP1 in the reduction of ROS in cold-exposed animals is currently under debate. Here, we investigated the role of UCP1 in mitochondrial ROS handling in the Lesser hedgehog tenrec (Echinops telfairi), a unique protoendothermic Malagasy mammal with recently identified brown adipose tissue (BAT). We show that the reduction of ROS by UCP1 activity also occurs in BAT mitochondria of the tenrec, suggesting that the antioxidative role of UCP1 is an ancient mammalian trait. Our analysis shows that the quantity of UCP1 displays strong control over mitochondrial hydrogen peroxide release, whereas other factors, such as mild cold, nonshivering thermogenesis, oxidative capacity, and mitochondrial respiration, do not correlate. Furthermore, hydrogen peroxide release from recoupled BAT mitochondria was positively associated with mitochondrial membrane potential. These findings led to a model of UCP1 controlling mitochondrial ROS release and, presumably, being controlled by high membrane potential, as proposed in the canonical model of “mild uncoupling”. Our study further promotes a conserved role for UCP1 in the prevention of oxidative stress, which was presumably established during evolution before UCP1 was physiologically integrated into nonshivering thermogenesis.     

Interpretation of variability of mitochondrial genomes in the species Aspergillus carbonarius

For interpretation of intraspecific polymorphism and the considerable differences in the size of mtDNAs among three groups of A. carbonarius, restriction maps were constructed from several enzymes. Functional maps were also developed to compare genome organisations and gene content. The appearance of various mtDNAs of A. carbonarius strains are different in size, but their gene content is almost identical. The 1.1 kb size difference between two closely related subgroups (1a, 1b) can be attributed to the presence or absence of an intron in cox2 gene. This phenomenon demonstrates that the migration of introns is possibly responsible for the development of variable mitochondrial genomes in nature. The striking differences in size and restriction patterns between two main mtDNA groups might derive from both the intronal variations and the altered intergenic organisation.     

The complete mitochondrial genomes of two ghost moths, Thitarodes renzhiensis and Thitarodes yunnanensis: the ancestral gene arrangement in Lepidoptera

ABSTRACT: BACKGROUND: Lepidoptera encompasses more than 160,000 described species that have been classified into 45-48 superfamilies. The previously determined Lepidoptera mitochondrial genomes (mitogenomes) are limited to six superfamilies of the most derived lepidopteran lineage Ditrysia. Compared with the ancestral insect gene order, these mitogenomes all contain a tRNA rearrangement. To gain new insights into Lepidoptera mitogenome evolution, we sequenced the mitogenomes of two ghost moths that belong to primitive lepidopteran lineages and conducted a comparative mitogenomic analysis across Lepidoptera. RESULTS: The mitogenomes of Thitarodes renzhiensis and T. yunnanensis are 16,173 bp and 15,814 bp long with an A+T content of 81.28% and 82.33%, respectively. Different tandem repeats in the A+T-rich region mainly account for the size difference between the two mitogenomes. Both mitogenomes include 13 protein-coding genes, 22 transfer RNA genes, and 2 ribosomal RNA genes. The 1,584-bp sequence from rrnS to nad2 was also determined for Thitarodes sp.QL, which has no repetitive sequence in the A+T-rich region. All three Thitarodes species possess the ancestral gene order with trnI-trnQ-trnM located between the A+T-rich region and nad2, which is different from the gene order trnM-trnI-trnQ in all previously sequenced Lepidoptera species. The formerly identified conserved elements of Lepidoptera mitogenomes (i.e. the motif 'ATAGA' and poly-T stretch in the A+T-rich region and the long intergenic spacer upstream of nad2) are absent in the Thitarodes mitogenomes. The phylogenetic analysis supports that Hepialoidea, represented by T. renzhiensis and T. yunnanensis, occupies a basal position in the currently sampled seven superfamilies. The relationships of the other six superfamilies are (((((Bombycoidea + Geometroidea) + Noctuoidea) + Pyraloidea) + Papilionoidea) + Tortricoidea). CONCLUSION: The mitogenomes of T. renzhiensis and T. yunnanensis exhibit unusual features compared with the previously determined Lepidoptera mitogenomes. Their ancestral gene order indicates that the tRNA rearrangement event occurred after Lepidoptera diverged from other holometabolous insect orders. Phylogenetic analysis based on mitogenome sequences is a power tool for addressing phylogenetic relationships among major Lepidoptera superfamilies. Characterization of the two ghost moth mitogenomes has enriched our knowledge of Lepidoptera mitogenomes and contributed to our understanding of the mechanisms underlying mitogenome evolution, especially gene rearrangements.