The mitochondrial genome of the chimpanzee louse,Pediculus schaeffi: insights into the process of mitochondrial genome fragmentation in the blood-sucking lice of great apes

Background

Blood-sucking lice in the genera Pediculus and Pthirus are obligate ectoparasites of great apes. Unlike most bilateral animals, which have 37 mitochondrial (mt) genes on a single circular chromosome, the sucking lice of humans have extensively fragmented mt genomes. The head louse, Pediculus capitis, and the body louse, Pe. humanus, have their 37 mt genes on 20 minichromosomes. The pubic louse, Pthirus pubis, has its 34 mt genes known on 14 minichromosomes. To understand the process of mt genome fragmentation in the sucking lice of great apes, we sequenced the mt genome of the chimpanzee louse, Pe. schaeffi, and compared it with the three human lice.

Results

We identified all of the 37 mt genes typical of bilateral animals in the chimpanzee louse; these genes are on 18 types of minichromosomes. Seventeen of the 18 minichromosomes of the chimpanzee louse have the same gene content and gene arrangement as their counterparts in the human head louse and the human body louse. However, five genes, cob, trnS₁, trnN, trnE and trnM, which are on three minichromosomes in the human head louse and the human body louse, are together on one minichromosome in the chimpanzee louse.

Conclusions

Using the human pubic louse, Pt. pubis, as an outgroup for comparison, we infer that a single minichromosome has fragmented into three in the lineage leading to the human head louse and the human body louse since this lineage diverged from the chimpanzee louse ~6 million years ago. Our results provide insights into the process of mt genome fragmentation in the sucking lice in a relatively fine evolutionary scale.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1843-3) contains supplementary material, which is available to authorized users.     

Heteroplasmy in the Mitochondrial Genomes of Human Lice and Ticks Revealed by High Throughput Sequencing

The typical mitochondrial (mt) genomes of bilateral animals consist of 37 genes on a single circular chromosome. The mt genomes of the human body louse, Pediculus humanus, and the human head louse, Pediculus capitis, however, are extensively fragmented and contain 20 minichromosomes, with one to three genes on each minichromosome. Heteroplasmy, i.e. nucleotide polymorphisms in the mt genome within individuals, has been shown to be significantly higher in the mt cox1 gene of human lice than in humans and other animals that have the typical mt genomes. To understand whether the extent of heteroplasmy in human lice is associated with mt genome fragmentation, we sequenced the entire coding regions of all of the mt minichromosomes of six human body lice and six human head lice from Ethiopia, China and France with an Illumina HiSeq platform. For comparison, we also sequenced the entire coding regions of the mt genomes of seven species of ticks, which have the typical mitochondrial genome organization of bilateral animals. We found that the level of heteroplasmy varies significantly both among the human lice and among the ticks. The human lice from Ethiopia have significantly higher level of heteroplasmy than those from China and France (P_t<0.05). The tick, Amblyomma cajennense, has significantly higher level of heteroplasmy than other ticks (P_t<0.05). Our results indicate that heteroplasmy level can be substantially variable within a species and among closely related species, and does not appear to be determined by single factors such as genome fragmentation.     

基于太平洋甲胁虱裂化线粒体基因组推测甲胁虱属祖先线粒体核型

[目的]动物典型的单一染色体线粒体基因组在甲胁虱属Hoplopleura已裂化成多个线粒体微环染色体.本研究旨在通过测定太平洋甲胁虱Hoplopleura pacifica的线粒体基因组来推测甲胁虱属祖先线粒体核型.[方法]利用Illumina HiSeq X Ten高通量测序技术对太平洋甲胁虱裂化线粒体基因组进行测定...     

Evolutionary history of mammalian sucking lice (Phthiraptera: Anoplura)

Background  

Sucking lice (Phthiraptera: Anoplura) are obligate, permanent ectoparasites of eutherian mammals, parasitizing members of 12 of the 29 recognized mammalian orders and approximately 20% of all mammalian species. These host specific, blood-sucking insects are morphologically adapted for life on mammals: they are wingless, dorso-ventrally flattened, possess tibio-tarsal claws for clinging to host hair, and have piercing mouthparts for feeding. Although there are more than 540 described species of Anoplura and despite the potential economical and medical implications of sucking louse infestations, this study represents the first attempt to examine higher-level anopluran relationships using molecular data. In this study, we use molecular data to reconstruct the evolutionary history of 65 sucking louse taxa with phylogenetic analyses and compare the results to findings based on morphological data. We also estimate divergence times among anopluran taxa and compare our results to host (mammal) relationships.     

吸虱亚目非典型线粒体基因组研究进展

吸虱是真兽类哺乳动物体表的专性寄生虫,在全世界广泛分布,物种数量高达540种。近年来随着分子生物学的飞速发展,在NCBI中已收录15种吸虱线粒体基因组序列,其独特的非典型线粒体基因组发生剧烈的裂化现象,形成数目不等的微环染色体。本文综述了15种吸虱线粒体基因组的结构和组成、RNA基因、非编码区以及对吸虱祖先线粒体核型推测的方法和结果。探讨了不同种属间以及与其它昆虫的差异,提出今后对吸虱亚目线粒体基因组研究的展望。     

Evolution of multipartite mitochondrial genomes in the booklice of the genus Liposcelis (Psocoptera)

Background

The genus Liposcelis (Psocoptera: Troctomorpha) has more than 120 species with a worldwide distribution and they pose a risk for global food security. The organization of mitochondrial (mt) genomes varies between the two species of booklice investigated in the genus Liposcelis. Liposcelis decolor has its mt genes on a single chromosome, like most other insects; L. bostrychophila, however, has a multipartite mt genome with genes on two chromosomes.

Results

To understand how multipartite mt genome organization evolved in the genus Liposcelis, we sequenced the mt genomes of L. entomophila and L. paeta in this study. We found that these two species of booklice also have multipartite mt genomes, like L. bostrychophila, with the mt genes we identified on two chromosomes. Numerous pseudo mt genes and non-coding regions were found in the mt genomes of these two booklice, and account for 30% and 10% respectively of the entire length we sequenced. In L. bostrychophila, the mt genes are distributed approximately equally between the two chromosomes. In L. entomophila and L. paeta, however, one mt chromosome has most of the genes we identified whereas the other chromosome has largely pseudogenes and non-coding regions. L. entomophila and L. paeta differ substantially from each other and from L. bostrychophila in gene content and gene arrangement in their mt chromosomes.

Conclusions

Our results indicate unusually fast evolution in mt genome organization in the booklice of the genus Liposcelis, and reveal different patterns of mt genome fragmentation among L. bostrychophila, L. entomophila and L. paeta.

Electronic supplementary material

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

The Complete Mitochondrial Genome of Gossypium hirsutum and Evolutionary Analysis of Higher Plant Mitochondrial Genomes

Background

Mitochondria are the main manufacturers of cellular ATP in eukaryotes. The plant mitochondrial genome contains large number of foreign DNA and repeated sequences undergone frequently intramolecular recombination. Upland Cotton (Gossypium hirsutum L.) is one of the main natural fiber crops and also an important oil-producing plant in the world. Sequencing of the cotton mitochondrial (mt) genome could be helpful for the evolution research of plant mt genomes.

Methodology/Principal Findings

We utilized 454 technology for sequencing and combined with Fosmid library of the Gossypium hirsutum mt genome screening and positive clones sequencing and conducted a series of evolutionary analysis on Cycas taitungensis and 24 angiosperms mt genomes. After data assembling and contigs joining, the complete mitochondrial genome sequence of G. hirsutum was obtained. The completed G.hirsutum mt genome is 621,884 bp in length, and contained 68 genes, including 35 protein genes, four rRNA genes and 29 tRNA genes. Five gene clusters are found conserved in all plant mt genomes; one and four clusters are specifically conserved in monocots and dicots, respectively. Homologous sequences are distributed along the plant mt genomes and species closely related share the most homologous sequences. For species that have both mt and chloroplast genome sequences available, we checked the location of cp-like migration and found several fragments closely linked with mitochondrial genes.

Conclusion

The G. hirsutum mt genome possesses most of the common characters of higher plant mt genomes. The existence of syntenic gene clusters, as well as the conservation of some intergenic sequences and genic content among the plant mt genomes suggest that evolution of mt genomes is consistent with plant taxonomy but independent among different species.     

A novel mitochondrial genome architecture in thrips (Insecta: Thysanoptera): extreme size asymmetry among chromosomes and possible recent control region duplication

Background

Multipartite mitochondrial genomes are very rare in animals but have been found previously in two insect orders with highly rearranged genomes, the Phthiraptera (parasitic lice), and the Psocoptera (booklice/barklice).

Results

We provide the first report of a multipartite mitochondrial genome architecture in a third order with highly rearranged genomes: Thysanoptera (thrips). We sequenced the complete mitochondrial genomes of two divergent members of the Scirtothrips dorsalis cryptic species complex. The East Asia 1 species has the single circular chromosome common to animals while the South Asia 1 species has a genome consisting of two circular chromosomes. The fragmented South Asia 1 genome exhibits extreme chromosome size asymmetry with the majority of genes on the large, 14.28 kb, chromosome and only nad6 and trnC on the 0.92 kb mini-circle chromosome. This genome also features paralogous control regions with high similarity suggesting a very recent origin of the nad6 mini-circle chromosome in the South Asia 1 cryptic species.

Conclusions

Thysanoptera, along with the other minor paraenopteran insect orders should be considered models for rapid mitochondrial genome evolution, including fragmentation. Continued use of these models will facilitate a greater understanding of recombination and other mitochondrial genome evolutionary processes across eukaryotes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1672-4) contains supplementary material, which is available to authorized users.     

The complete mitochondrial genomes of six heterodont bivalves (Tellinoidea and Solenoidea): variable gene arrangements and phylogenetic implications

Background

Taxonomy and phylogeny of subclass Heterodonta including Tellinoidea are long-debated issues and a complete agreement has not been reached yet. Mitochondrial (mt) genomes have been proved to be a powerful tool in resolving phylogenetic relationship. However, to date, only ten complete mitochondrial genomes of Heterodonta, which is by far the most diverse major group of Bivalvia, have been determined. In this paper, we newly sequenced the complete mt genomes of six species belonging to Heterodonta in order to resolve some problematical relationships among this subclass.

Principal Findings

The complete mt genomes of six species vary in size from 16,352 bp to 18,182. Hairpin-like secondary structures are found in the largest non-coding regions of six freshly sequenced mt genomes, five of which contain tandem repeats. It is noteworthy that two species belonging to the same genus show different gene arrangements with three translocations. The phylogenetic analysis of Heterodonta indicates that Sinonovacula constricta, distant from the Solecurtidae belonging to Tellinoidea, is as a sister group with Solen grandis of family Solenidae. Besides, all five species of Tellinoidea cluster together, while Sanguinolaria diphos has closer relationship with Solecurtus divaricatus, Moerella iridescens and Semele scaba rather than with Sanguinolaria olivacea.

Conclusions/Significance

By comparative study of gene order rearrangements and phylogenetic relationships of the five species belonging to Tellinoidea, our results support that comparisons of mt gene order rearrangements, to some extent, are a useful tool for phylogenetic studies. Based on phylogenetic analyses of multiple protein-coding genes, we prefer classifying the genus Sinonovacula within the superfamily Solenoidea and not the superfamily Tellinoidea. Besides, both gene order and sequence data agree that Sanguinolaria (Psammobiidae) is not monophyletic. Nevertheless, more studies based on more mt genomes via combination of gene order and phylogenetic analysis are needed to further understand the phylogenetic relationships in subclass Heterodonta.     

Chimeric mitochondrial minichromosomes of the human body louse, Pediculus humanus: evidence for homologous and non-homologous recombination

The mitochondrial (mt) genome of the human body louse, Pediculus humanus, consists of 18 minichromosomes. Each minichromosome is 3 to 4 kb long and has 1 to 3 genes. There is unequivocal evidence for recombination between different mt minichromosomes in P. humanus. It is not known, however, how these minichromosomes recombine. Here, we report the discovery of eight chimeric mt minichromosomes in P. humanus. We classify these chimeric mt minichromosomes into two groups: Group I and Group II. Group I chimeric minichromosomes contain parts of two different protein-coding genes that are from different minichromosomes. The two parts of protein-coding genes in each Group I chimeric minichromosome are joined at a microhomologous nucleotide sequence; microhomologous nucleotide sequences are hallmarks of non-homologous recombination. Group II chimeric minichromosomes contain all of the genes and the non-coding regions of two different minichromosomes. The conserved sequence blocks in the non-coding regions of Group II chimeric minichromosomes resemble the "recombination repeats" in the non-coding regions of the mt genomes of higher plants. These repeats are essential to homologous recombination in higher plants. Our analyses of the nucleotide sequences of chimeric mt minichromosomes indicate both homologous and non-homologous recombination between minichromosomes in the mitochondria of the human body louse.     

Diversity of 23S rRNA Genes within Individual Prokaryotic Genomes

Background

The concept of ribosomal constraints on rRNA genes is deduced primarily based on the comparison of consensus rRNA sequences between closely related species, but recent advances in whole-genome sequencing allow evaluation of this concept within organisms with multiple rRNA operons.

Methodology/Principal Findings

Using the 23S rRNA gene as an example, we analyzed the diversity among individual rRNA genes within a genome. Of 184 prokaryotic species containing multiple 23S rRNA genes, diversity was observed in 113 (61.4%) genomes (mean 0.40%, range 0.01%–4.04%). Significant (1.17%–4.04%) intragenomic variation was found in 8 species. In 5 of the 8 species, the diversity in the primary structure had only minimal effect on the secondary structure (stem versus loop transition). In the remaining 3 species, the diversity significantly altered local secondary structure, but the alteration appears minimized through complex rearrangement. Intervening sequences (IVS), ranging between 9 and 1471 nt in size, were found in 7 species. IVS in Deinococcus radiodurans and Nostoc sp. encode transposases. T. tengcongensis was the only species in which intragenomic diversity >3% was observed among 4 paralogous 23S rRNA genes.

Conclusions/Significance

These findings indicate tight ribosomal constraints on individual 23S rRNA genes within a genome. Although classification using primary 23S rRNA sequences could be erroneous, significant diversity among paralogous 23S rRNA genes was observed only once in the 184 species analyzed, indicating little overall impact on the mainstream of 23S rRNA gene-based prokaryotic taxonomy.     

Rapid sequencing of the bamboo mitochondrial genome using Illumina technology and parallel episodic evolution of organelle genomes in grasses

Background

Compared to their counterparts in animals, the mitochondrial (mt) genomes of angiosperms exhibit a number of unique features. However, unravelling their evolution is hindered by the few completed genomes, of which are essentially Sanger sequenced. While next-generation sequencing technologies have revolutionized chloroplast genome sequencing, they are just beginning to be applied to angiosperm mt genomes. Chloroplast genomes of grasses (Poaceae) have undergone episodic evolution and the evolutionary rate was suggested to be correlated between chloroplast and mt genomes in Poaceae. It is interesting to investigate whether correlated rate change also occurred in grass mt genomes as expected under lineage effects. A time-calibrated phylogenetic tree is needed to examine rate change.

Methodology/Principal Findings

We determined a largely completed mt genome from a bamboo, Ferrocalamus rimosivaginus (Poaceae), through Illumina sequencing of total DNA. With combination of de novo and reference-guided assembly, 39.5-fold coverage Illumina reads were finally assembled into scaffolds totalling 432,839 bp. The assembled genome contains nearly the same genes as the completed mt genomes in Poaceae. For examining evolutionary rate in grass mt genomes, we reconstructed a phylogenetic tree including 22 taxa based on 31 mt genes. The topology of the well-resolved tree was almost identical to that inferred from chloroplast genome with only minor difference. The inconsistency possibly derived from long branch attraction in mtDNA tree. By calculating absolute substitution rates, we found significant rate change (∼4-fold) in mt genome before and after the diversification of Poaceae both in synonymous and nonsynonymous terms. Furthermore, the rate change was correlated with that of chloroplast genomes in grasses.

Conclusions/Significance

Our result demonstrates that it is a rapid and efficient approach to obtain angiosperm mt genome sequences using Illumina sequencing technology. The parallel episodic evolution of mt and chloroplast genomes in grasses is consistent with lineage effects.     

Comparative molecular cytogenetic analyses of a major tandemly repeated DNA family and retrotransposon sequences in cultivated jute Corchorus species (Malvaceae)

Background and Aims

The cultivated jute species Corchorus olitorius and Corchorus capsularis are important fibre crops. The analysis of repetitive DNA sequences, comprising a major part of plant genomes, has not been carried out in jute but is useful to investigate the long-range organization of chromosomes. The aim of this study was the identification of repetitive DNA sequences to facilitate comparative molecular and cytogenetic studies of two jute cultivars and to develop a fluorescent in situ hybridization (FISH) karyotype for chromosome identification.

Methods

A plasmid library was generated from C. olitorius and C. capsularis with genomic restriction fragments of 100–500 bp, which was complemented by targeted cloning of satellite DNA by PCR. The diversity of the repetitive DNA families was analysed comparatively. The genomic abundance and chromosomal localization of different repeat classes were investigated by Southern analysis and FISH, respectively. The cytosine methylation of satellite arrays was studied by immunolabelling.

Key Results

Major satellite repeats and retrotransposons have been identified from C. olitorius and C. capsularis. The satellite family CoSat I forms two undermethylated species-specific subfamilies, while the long terminal repeat (LTR) retrotransposons CoRetro I and CoRetro II show similarity to the Metaviridea of plant retroelements. FISH karyotypes were developed by multicolour FISH using these repetitive DNA sequences in combination with 5S and 18S–5·8S–25S rRNA genes which enable the unequivocal chromosome discrimination in both jute species.

Conclusions

The analysis of the structure and diversity of the repeated DNA is crucial for genome sequence annotation. The reference karyotypes will be useful for breeding of jute and provide the basis for karyotyping homeologous chromosomes of wild jute species to reveal the genetic and evolutionary relationship between cultivated and wild Corchorus species.     

Sequence Divergence and Conservation in Genomes of Helicobacter cetorum Strains from a Dolphin and a Whale

Background and Objectives

Strains of Helicobacter cetorum have been cultured from several marine mammals and have been found to be closely related in 16 S rDNA sequence to the human gastric pathogen H. pylori, but their genomes were not characterized further.

Methods

The genomes of H. cetorum strains from a dolphin and a whale were sequenced completely using 454 technology and PCR and capillary sequencing.

Results

These genomes are 1.8 and 1.95 mb in size, some 7–26% larger than H. pylori genomes, and differ markedly from one another in gene content, and sequences and arrangements of shared genes. However, each strain is more related overall to H. pylori and its descendant H. acinonychis than to other known species. These H. cetorum strains lack cag pathogenicity islands, but contain novel alleles of the virulence-associated vacuolating cytotoxin (vacA) gene. Of particular note are (i) an extra triplet of vacA genes with ≤50% protein-level identity to each other in the 5′ two-thirds of the gene needed for host factor interaction; (ii) divergent sets of outer membrane protein genes; (iii) several metabolic genes distinct from those of H. pylori; (iv) genes for an iron-cofactored urease related to those of Helicobacter species from terrestrial carnivores, in addition to genes for a nickel co-factored urease; and (v) members of the slr multigene family, some of which modulate host responses to infection and improve Helicobacter growth with mammalian cells.

Conclusions

Our genome sequence data provide a glimpse into the novelty and great genetic diversity of marine helicobacters. These data should aid further analyses of microbial genome diversity and evolution and infection and disease mechanisms in vast and often fragile ocean ecosystems.     

虱目裂化线粒体基因组研究进展

虱目是哺乳类和鸟类体表的专性寄生虫。在虱科、阴虱科、长角鸟虱科和兽羽虱科的某些寄生虱种中发现了线粒体基因组裂化现象, 其线粒体基因组裂化成了多个环状的线粒体染色体, 如体虱(Pediculus humanus)、头虱(pediculus capitis)和阴虱(Pthirus pubis)的线粒体基因组分别裂化形成20个、20个和14个微环染色体。微环染色体可能是基因删除和同源重组的结果, 关于线粒体基因组裂化的具体原因和机制, 目前并不清楚, 推测可能是进化选择或随机遗传漂变的结果或与线粒体单链DNA结合蛋白的缺失有关。鉴于线粒体基因组裂化研究对于深入理解线粒体的起源和进化方面具有重要意义, 文章以虱目裂化线粒体基因组为主线, 列举了动物裂化线粒体基因组和裂化特征, 阐述了虱目裂化线粒体基因组的研究现状, 分析了虱目线粒体基因组裂化的类型、原因和机制, 并对该领域未来的研究方向进行了展望。     

Different dynamics of genome content shuffling among host-specificity groups of the symbiotic actinobacterium Frankia

Background

Frankia is a genus of soil actinobacteria forming nitrogen-fixing root-nodule symbiotic relationships with non-leguminous woody plant species, collectively called actinorhizals, from eight dicotyledonous families. Frankia strains are classified into four host-specificity groups (HSGs), each of which exhibits a distinct host range. Genome sizes of representative strains of Alnus, Casuarina, and Elaeagnus HSGs are highly diverged and are positively correlated with the size of their host ranges.

Results

The content and size of 12 Frankia genomes were investigated by in silico comparative genome hybridization and pulsed-field gel electrophoresis, respectively. Data were collected from four query strains of each HSG and compared with those of reference strains possessing completely sequenced genomes. The degree of difference in genome content between query and reference strains varied depending on HSG. Elaeagnus query strains were missing the greatest number (22–32%) of genes compared with the corresponding reference genome; Casuarina query strains lacked the fewest (0–4%), with Alnus query strains intermediate (14–18%). In spite of the remarkable gene loss, genome sizes of Alnus and Elaeagnus query strains were larger than would be expected based on total length of the absent genes. In contrast, Casuarina query strains had smaller genomes than expected.

Conclusions

The positive correlation between genome size and host range held true across all investigated strains, supporting the hypothesis that size and genome content differences are responsible for observed diversity in host plants and host plant biogeography among Frankia strains. In addition, our results suggest that different dynamics of shuffling of genome content have contributed to these symbiotic and biogeographic adaptations. Elaeagnus strains, and to a lesser extent Alnus strains, have gained and lost many genes to adapt to a wide range of environments and host plants. Conversely, rather than acquiring new genes, Casuarina strains have discarded genes to reduce genome size, suggesting an evolutionary orientation towards existence as specialist symbionts.

Electronic supplementary material

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