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.     

Description of Babesia duncani n.sp. (Apicomplexa: Babesiidae) from humans and its differentiation from other piroplasms

The morphologic, ultrastructural and genotypic characteristics of Babesia duncani n.sp. are described based on the characterization of two isolates (WA1, CA5) obtained from infected human patients in Washington and California. The intraerythrocytic stages of the parasite are morphologically indistinguishable from Babesia microti, which is the most commonly identified cause of human babesiosis in the USA. Intraerythrocytic trophozoites of B. duncani n.sp. are round to oval, with some piriform, ring and ameboid forms. Division occurs by intraerythrocytic schizogony, which results in the formation of merozoites in tetrads (syn. Maltese cross or quadruplet forms). The ultrastructural features of trophozoites and merozoites are similar to those described for B. microti and Theileria spp. However, intralymphocytic schizont stages characteristic of Theileria spp. have not been observed in infected humans. In phylogenetic analyses based on sequence data for the complete18S ribosomal RNA gene, B. duncani n.sp. lies in a distinct clade that includes isolates from humans, dogs and wildlife in the western United States but separate from Babesia sensu stricto, Theileria spp. and B. microti. ITS2 sequence analysis of the B. duncani n.sp. isolates (WA1, CA5) show that they are phylogenetically indistinguishable from each other and from two other human B. duncani-type parasites (CA6, WA2 clone1) but distinct from other Babesia and Theileria species sequenced. This analysis provides robust molecular support that the B. duncani n.sp. isolates are monophyletic and the same species. The morphologic characteristics together with the phylogenetic analysis of two genetic loci support the assertion that B. duncani n.sp. is a distinct species from other known Babesia spp. for which morphologic and sequence information are available.     

Insights into the demographic history of African Pygmies from complete mitochondrial genomes

Pygmy populations are among the few hunter-gatherers currently living in sub-Saharan Africa and are mainly represented by two groups, Eastern and Western, according to their current geographical distribution. They are scattered across the Central African belt and surrounded by Bantu-speaking farmers, with whom they have complex social and economic interactions. To investigate the demographic history of Pygmy groups, a population approach was applied to the analysis of 205 complete mitochondrial DNA (mtDNA) sequences from ten central African populations. No sharing of maternal lineages was observed between the two Pygmy groups, with haplogroup L1c being characteristic of the Western group but most of Eastern Pygmy lineages falling into subclades of L0a, L2a, and L5. Demographic inferences based on Bayesian coalescent simulations point to an early split among the maternal ancestors of Pygmies and those of Bantu-speaking farmers (～ 70,000 years ago [ya]). Evidence for population growth in the ancestors of Bantu-speaking farmers has been observed, starting ～ 65,000 ya, well before the diffusion of Bantu languages. Subsequently, the effective population size of the ancestors of Pygmies remained constant over time and ～ 27,000 ya, coincident with the Last Glacial Maximum, Eastern and Western Pygmies diverged, with evidence of subsequent migration only among the Western group and the Bantu-speaking farmers. Western Pygmies show signs of a recent bottleneck 4,000-650 ya, coincident with the diffusion of Bantu languages, whereas Eastern Pygmies seem to have experienced a more ancient decrease in population size (20,000-4,000 ya). In conclusion, the results of this first attempt at analyzing complete mtDNA sequences at the population level in sub-Saharan Africa not only support previous findings but also offer new insights into the demographic history of Pygmy populations, shedding new light on the ancient peopling of the African continent.     

Evolution of the apicoplast and its hosts: from heterotrophy to autotrophy and back again

The photosynthetic origin of apicomplexan parasites was proposed upon the discovery of a reduced non-photosynthetic plastid termed the apicoplast in their cells. Although it is clear that the apicoplast has evolved through a secondary endosymbiosis, its particular origin within the red or green plastid lineage remains controversial. The recent discovery of Chromera velia, the closest known photosynthetic relative to apicomplexan parasites, sheds new light on the evolutionary history of alveolate plastids. Here we review our knowledge on the evolutionary history of Apicomplexa and particularly their plastids, with a focus on the pathway by which they evolved from free-living heterotrophs through photoautotrophs to omnipresent obligatory intracellular parasites. New sequences from C. velia (histones H2A, H2B; GAPDH, TufA) and phylogenetic analyses are also presented and discussed here.     

The mitochondrial genomes of three skippers: Insights into the evolution of the family Hesperiidae (Lepidoptera)

We sequenced the mitogenomes of Astictopterus jama, Isoteinon lamprospilus and Notocrypta curvifascia to obtain further insight into the mitogenomic architecture evolution and performed phylogenetic reconstruction using 29 Hesperiidae mitogenome sequences. The complete mitogenome sequences of A. jama, I. lamprospilus and N. curvifascia are 15,430, 15,430 and 15,546 bp in size, respectively. All contain 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes, and an A + T-rich region. Nucleotide composition is A + T biased, and the majority of the protein-coding genes exhibit a negative AT-skew, which is reflected in the nucleotide composition, codon, and amino acid usage. The A + T-rich region is comprised of nonrepetitive sequences, including the motif ATAGA followed by a poly-T stretch, a microsatellite-like element next to the ATTTA motif, and a poly-A adjacent to tRNAs. Although most genes evolve under a strong purifying selection, the entire nad gene family (especially nad6) exhibits somewhat relaxed purifying selection, and atp8, evolving under a highly relaxed selection, is an outlier in the family Hesperiidae. Several different approaches relatively consistently indicated that nad6, atp8 and nad4 are comparatively fast-evolving genes in this family, which may have implications for future phylogenetic, population genetics and species diagnostics studies. For phylogenetic analyses of Hesperiidae, we tested a few datasets, and found that the one comprising all 37 genes produced the highest node support, indicating that the inclusion of RNAs improves the phylogenetic signal. Results indicate that subfamilies Euschemoninae, Heteropterinae, and Coeliadinae are monophyletic with strong nodal support, but Pyrginae and Eudaminae are paraphyletic. Finally, we confirm that A. jama and I. lamprospilus are close relatives.     

A reevaluation of rice mitochondrial evolution based on the complete sequence of male-fertile and male-sterile mitochondrial genomes

Plant mitochondrial genomes have features that distinguish them radically from their animal counterparts: a high rate of rearrangement, of uptake and loss of DNA sequences, and an extremely low point mutation rate. Perhaps the most unique structural feature of plant mitochondrial DNAs is the presence of large repeated sequences involved in intramolecular and intermolecular recombination. In addition, rare recombination events can occur across shorter repeats, creating rearrangements that result in aberrant phenotypes, including pollen abortion, which is known as cytoplasmic male sterility (CMS). Using next-generation sequencing, we pyrosequenced two rice (Oryza sativa) mitochondrial genomes that belong to the indica subspecies. One genome is normal, while the other carries the wild abortive-CMS. We find that numerous rearrangements in the rice mitochondrial genome occur even between close cytotypes during rice evolution. Unlike maize (Zea mays), a closely related species also belonging to the grass family, integration of plastid sequences did not play a role in the sequence divergence between rice cytotypes. This study also uncovered an excellent candidate for the wild abortive-CMS-encoding gene; like most of the CMS-associated open reading frames that are known in other species, this candidate was created via a rearrangement, is chimeric in structure, possesses predicted transmembrane domains, and coopted the promoter of a genuine mitochondrial gene. Our data give new insights into rice mitochondrial evolution, correcting previous reports.     

Genome sequence and transcriptome of Sorbus pohuashanensis provide insights into population evolution and leaf sunburn response

Sorbus pohuashanensis (Hance) Hedl. is a potential horticulture and medicinal plant, but its genomic and genetic backgrounds remain unknown. Here, we sequence and assemble the S. pohuashanensis reference genome using PacBio long reads. Based on the new reference genome, we resequence a core collection of 22 Sorbus spp. samples, which are divided into 2 groups (G1 and G2) based on phylogenetic and PCA analyses. These phylogenetic clusters are highly consistent with their classification based on leaf shape. Natural hybridization between the G1 and G2 groups is evidenced by a sample (R21) with a highly heterozygous genotype. Nucleotide diversity (π) analysis shows that G1 has a higher diversity than G2 and that G2 originated from G1. During the evolution process, the gene families involved in photosynthesis pathways expanded and the gene families involved in energy consumption contracted. RNA-seq data suggests that flavonoid biosynthesis and heat-shock protein (HSP)-heat-shock factor (HSF) pathways play important roles in protection against sunburn. This study provides new insights into the evolution of Sorbus spp. genomes. In addition, the genomic resources, and the identified genetic variations, especially those related to stress resistance, will help future efforts to produce and breed Sorbus spp.     

The sequence,organization, and evolution of the Locusta migratoria mitochondrial genome

The sequencing of the cloned Locusta migratoria mitochondrial genome has been completed. The sequence is 15,722 by in length and contains 75.3% A+T, the lowest value in any of the five insect mitochondrial sequences so far determined. The protein coding genes have a similar A+T content (74.1%) but are distinguished by a high cytosine content at the third codon position. The gene content and organization are the same as in Drosophila yakuba except for a rearrangement of the two tRNA genes tRNA^lys and tRNA^asp. The A+T-rich region has a lower A+T nucleotide content than in other insects, and this is largely due to the presence of two G+C-rich 155-bp repetitive sequences at the 5 end of this section and the beginning of the adjacent small rRNA gene. The sizes of the large and small rRNA genes are 1,314 and 827 bp, respectively, and both sequences can be folded to form secondary structures similar to those previously predicted for Drosophila. The tRNA genes have also been modeled and these show a strong resemblance to the dipteran tRNAs, all anticodons apparently being conserved between the two species. A comparison of the protein coding nucleotide sequences of the locust DNA with the homologous sequences of five other arthropods (Drosophila yakuba, Anopheles quadrimaculatus, Anopheles gambiae, Apis mellifera, and Artemia franciscana) was performed. The amino acid composition of the encoded proteins in Locusta is similar to that of Drosophila, with a Dayhoff distance twice that of the distance between the fruit fly and the mosquitoes. A phylogenetic analysis revealed the locust genes to be more similar to those of the Dipterans than to those of the honeybee at both the nucleotide and amino acid levels. A comparative analysis of tRNA orders, using crustacean mtDNAs as outgroups, supported this. This high level of divergence in the Apis genome has been noted elsewhere and is possibly an effect of directional mutation pressure having resulted in an accelerated pattern of sequence evolution. If the general assumption that the Holometabola are monophyletic holds, then these results emphasize the difficulties of reconstructing phylogenies that include lineages with variable substitution rates and base composition biases. The need to exercise caution in using information about tRNA gene orders in phylogenetic analysis is also illustrated. However, if the honeybee sequence is excluded, the correspondence between the other five arthropod sequences supports the findings of previous studies which have endorsed the use of mtDNA sequences for studies of phylogeny at deep levels of taxonomy when mutation rates are equivalent. Correspondence to: P.K. Flook     

Complete DNA sequences of the mitochondrial genomes of the pathogenic yeasts Candida orthopsilosis and Candida metapsilosis: insight into the evolution of linear DNA genomes from mitochondrial telomere mutants

We determined complete mitochondrial DNA sequences of the two yeast species, Candida orthopsilosis and Candida metapsilosis, and compared them with the linear mitochondrial genome of their close relative, C.parapsilosis. Mitochondria of all the three species harbor compact genomes encoding the same set of genes arranged in the identical order. Differences in the length of these genomes result mainly from the presence/absence of introns. Multiple alterations were identified also in the sequences of the ribosomal and transfer RNAs, and proteins. However, the most striking feature of C.orthopsilosis and C.metapsilosis is the existence of strains differing in the molecular form of the mitochondrial genome (circular-mapping versus linear). Their analysis opens a unique window for understanding the role of mitochondrial telomeres in the stability and evolution of molecular architecture of the genome. Our results indicate that the circular-mapping mitochondrial genome derived from the linear form by intramolecular end-to-end fusions. Moreover, we suggest that the linear mitochondrial genome evolved from a circular-mapping form present in a common ancestor of the three species and, at the same time, the emergence of mitochondrial telomeres enabled the formation of linear monomeric DNA forms. In addition, comparison of isogenic C.metapsilosis strains differing in the form of the organellar genome suggests a possibility that, under some circumstances, the linearity and/or the presence of telomeres provide a competitive advantage over a circular-mapping mitochondrial genome.     

Modelling the evolution of genomes with integrated external and internal functions

The genomes that organisms transmit between generations contain information about different kinds of functions. The genome with the "best" mix and number of genes for these functions is the one that natural selection favours. Here I introduce a new way to model simple organisms with genes for external and internal functions, and use it to study the evolution of genome size. The external functions are exemplified by resource use and the internal functions by mutation control (repair). It is shown that even with a suitable proportion of genes for mutation control, the genomes in the organisms do not forever incorporate genes that increase resource use. Instead they evolve towards an optimal genome of limited size. The optimal proportion of genes for mutation control is shown to have an upper limit given by the ease with which transmission accuracy is improved by adding extra genes for this purpose to the genome. The model illustrates how natural selection on genomes integrates systems for the transmission of genetic information with systems relating to the external adaptation of the organism. It also opens up for other, more detailed theoretical investigations of genome functions.     

Localisation of Plasmodium falciparum uroporphyrinogen III decarboxylase of the heme-biosynthetic pathway in the apicoplast and characterisation of its catalytic properties

Uroporphyrinogen decarboxylase (UROD) is a key enzyme in the heme-biosynthetic pathway and in Plasmodium falciparum it occupies a strategic position in the proposed hybrid pathway for heme biosynthesis involving shuttling of intermediates between different subcellular compartments in the parasite. In the present study, we demonstrate that an N-terminally truncated recombinant P. falciparum UROD (r(Δ)PfUROD) over-expressed and purified from Escherichia coli cells, as well as the native enzyme from the parasite were catalytically less efficient compared with the host enzyme, although they were similar in other enzyme parameters. Molecular modeling of PfUROD based on the known crystal structure of the human enzyme indicated that the protein manifests a distorted triose phosphate isomerase (TIM) barrel fold which is conserved in all the known structures of UROD. The parasite enzyme shares all the conserved or invariant amino acid residues at the active and substrate binding sites, but is rich in lysine residues compared with the host enzyme. Mutation of specific lysine residues corresponding to residues at the dimer interface in human UROD enhanced the catalytic efficiency of the enzyme and dimer stability indicating that the lysine rich nature and weak dimer interface of the wild-type PfUROD could be responsible for its low catalytic efficiency. PfUROD was localised to the apicoplast, indicating the requirement of additional mechanisms for transport of the product coproporphyrinogen to other subcellular sites for its further conversion and ultimate heme formation.     

三株新城疫广西分离株全基因组序列的测定与分析

根据GenBank上所公布的新城疫病毒(NDV)的全基因组序列,设计了8对引物,运用RT-PCR 方法获取了3株广西地方强毒株GX7/02、GX9/03和GX11/03的全基因组序列,并对其进行了比较分析.此3株病毒的全基因组序列均由15192个碱基组成,与GenBank公布的ZJ1、U.S/Largo/71、Italy/2736/00等7个毒株的全基因组序列长度相同,比LaSota、Clone-30和B1的全基因组序列多出6个核苷酸,此6个核苷酸位于np基因的非编码区内,相对与NDV毒株LaSota、Clone-30和B1序列的1647～1648nt 位.通过序列的比较分析,发现GX7/02、GX9/03和GX11/03与ZJ1毒株的同源性较高,而与LaSota、Clone-30和B1等毒株的同源性较低.     

Modeling within-host HIV-1 dynamics and the evolution of drug resistance: trade-offs between viral enzyme function and drug susceptibility

There are many biological steps between viral infection of CD4(+) T cells and the production of HIV-1 virions. Here we incorporate an eclipse phase, representing the stage in which infected T cells have not started to produce new virus, into a simple HIV-1 model. Model calculations suggest that the quicker infected T cells progress from the eclipse stage to the productively infected stage, the more likely that a viral strain will persist. Long-term treatment effectiveness of antiretroviral drugs is often hindered by the frequent emergence of drug resistant virus during therapy. We link drug resistance to both the rate of progression of the eclipse phase and the rate of viral production of the resistant strain, and explore how the resistant strain could evolve to maximize its within-host viral fitness. We obtained the optimal progression rate and the optimal viral production rate, which maximize the fitness of a drug resistant strain in the presence of drugs. We show that the window of opportunity for invasion of drug resistant strains is widened for a higher level of drug efficacy provided that the treatment is not potent enough to eradicate both the sensitive and resistant virus.     

三株新城疫广西分离株全基因组序列的测定与分析

根据GenBank上所公布的新城疫病毒(NDV)的全基因组序列,设计了8对引物,运用RT-PCR方法获取了3株广西地方强毒株GX7/02、GX9/03和GX11/03的全基因组序列,并对其进行了比较分析。此3株病毒的全基因组序列均由15192个碱基组成,与GenBank公布的ZJ1、U.S/Largo/71、Italy/2736/00等7个毒株的全基因组序列长度相同,比LaSota、Clone-30和B1的全基因组序列多出6个核苷酸,此6个核苷酸位于np基因的非编码区内,相对与NDV毒株LaSota、Clone-30和B1序列的1647~1648nt位。通过序列的比较分析,发现GX7/02、GX9/03和GX11/03与ZJ1毒株的同源性较高,而与LaSota、Clone-30和B1等毒株的同源性较低。     

2009-2011年胆石症患者胆汁肠球菌菌群分布及菌种耐药性变迁

目的 了解2009年1月年2011年12月胆石症患者胆汁中肠球菌的分布及菌种耐药性变迁,指导临床合理应用抗菌药物.方法 采用法国生物梅里埃公司生产的VITEK-2细菌分析仪对胆汁标本中分离的肠球菌做鉴定及药敏试验.结果 2009年1月至2010年6月共分离到74株肠球菌;2010年7月至2011年12月共分离到115株肠球菌,排在前两位的均为粪肠球菌和屎肠球菌.药敏结果显示粪肠球菌对除红霉素和喹努普汀/达福普汀外的其他抗菌药物具有较高的抑菌活性;屎肠球菌对除喹努普汀/达福普汀、利奈唑胺、万古霉素外的其他抗菌药物具有较高的耐药率;首次检测到耐利奈唑胺的肠球菌.结论 胆石症患者胆汁中肠球菌感染以粪肠球菌和屎肠球菌为主,且多种细菌对抗菌药物均有不同程度的耐药,因此加强细菌耐药监测对临床合理使用抗菌药物有重要的参考意义.     

The complete chloroplast and mitochondrial genome sequences of Boea hygrometrica: insights into the evolution of plant organellar genomes

The complete nucleotide sequences of the chloroplast (cp) and mitochondrial (mt) genomes of resurrection plant Boea hygrometrica (Bh, Gesneriaceae) have been determined with the lengths of 153,493 bp and 510,519 bp, respectively. The smaller chloroplast genome contains more genes (147) with a 72% coding sequence, and the larger mitochondrial genome have less genes (65) with a coding faction of 12%. Similar to other seed plants, the Bh cp genome has a typical quadripartite organization with a conserved gene in each region. The Bh mt genome has three recombinant sequence repeats of 222 bp, 843 bp, and 1474 bp in length, which divide the genome into a single master circle (MC) and four isomeric molecules. Compared to other angiosperms, one remarkable feature of the Bh mt genome is the frequent transfer of genetic material from the cp genome during recent Bh evolution. We also analyzed organellar genome evolution in general regarding genome features as well as compositional dynamics of sequence and gene structure/organization, providing clues for the understanding of the evolution of organellar genomes in plants. The cp-derived sequences including tRNAs found in angiosperm mt genomes support the conclusion that frequent gene transfer events may have begun early in the land plant lineage.