Intraindividual sequence diversity and a hierarchical approach to the study of mitochondrial DNA mutations

Complete Mitochondrial Genome of Phytophthora nicotianae and Identification of Molecular Markers for the Oomycetes

Phytophthora nicotianae is one of the most destructive plant pathogens affecting a variety of plants, causing black shank of tobacco, among several other devastating diseases. Herein, we assembled the mitochondrial genome of P. nicotianae and analyzed its gene content and genome structure, performed comparative mitochondrial genomics analysis, and assessed phylogenetic relationships among oomycetes species. The circular mitogenome is 37,561 bp long, with 38 protein-coding genes, 25 transfer RNA (tRNA) genes, and 2 ribosomal RNA genes (rrnl and rrns). The mitochondrial genome showed a biased A/T usage versus G/C. The overall gene content and size of the P. nicotianae mitogenome are identical to those of other published Phytophthora mitogenomes. Interestingly, collinearity analysis using an existing ∼10 k inversion region (including 11 genes and 8 tRNAs) revealed that Phytophthora andina, Phytophthora infestans, Phytophthora mirabilis, Phytophthora ipomoeae, and Phytophthora phaseoli differed from Phytophthora nicotianae, Phytophthora sojae, Phytophthora ramorum, and Phytophthora polonica. Moreover, inverted repeat regions were found to be absent among species of the Peronosporales when compared with species from the Pythiales and Saprolegniales. A phylogenomic investigation based on 29 protein-coding genes demonstrated that Phytophthora is monophyletic, and placed P. nicotianae close to the clade including P. mirabilis, P. ipomoeae, P. andina, P. infestans, and P. phaseoli. Furthermore, we discovered six new candidate DNA molecular markers (rpl6, atp8, nad11, rps2, rps3, and rps4) based on these mitogenomes that would be suitable for species identification in the oomycetes, which have the same identification level as the whole mitogenome and ribosomal DNA sequences. These new molecular markers can not only provide a quick preview of the species without mitogenome information, but will also help to gain better understanding of the oomycetes pathogens and developing treatment or monitoring strategies.

Mitochondrial DNA assessment of Phytophthora infestans isolates from potato and tomato in Ethiopia reveals unexpected diversity

Mitochondrial DNA (mtDNA) haplotypes were determined using restriction fragment length polymorphism (RFLP) for P. infestans sampled from 513 foliar lesions of late blight found on potato and tomato in different regions of Ethiopia. Among the four reported mitochondrial haplotypes of Phytophthora infestans, Ia, Ib and IIb were detected in 93 % of the samples analyzed but the vast majority of these were Ia. The remaining 7 % represented a previously unreported haplotype. DNA sequencing of this new haplotype also confirmed a single base nucleotide substitution that resulted in loss of EcoRI restriction site and gain of two additional MspI sites in cox1 and atp1 genes, respectively. There were 28 polymorphic sites among all nucleotide sequences including five reference isolates. Sites with alignment gaps were observed in P4 with one nucleotide deletion in 11 Ethiopian isolates. None of the reference sequence produced frame-shifts, with the exception of the 3-nucleotide deletion in the P4 region by Phytophthora andina, a feature that can be used to distinguish the new Ethiopian isolates from P. andina. While a distinguishing molecular data presented here clearly separated them from P. infestans, 7 % of the isolates that share this feature formed an important component of the late blight pathogen causing disease on Solanum tuberosum in Ethiopia. Thus, these Ethiopian isolates could represent a novel Phytophthora species reported for the first time here.

Evolutionary history of the mitochondrial genome in Mycosphaerella populations infecting bread wheat, durum wheat and wild grasses

Plant pathogens emerge in agro-ecosystems following different evolutionary mechanisms over different time scales. Previous analyses based on sequence variation at six nuclear loci indicated that Mycosphaerella graminicola diverged from an ancestral population adapted to wild grasses during the process of wheat domestication approximately 10,500 years ago. We tested this hypothesis by conducting coalescence analyses based on four mitochondrial loci using 143 isolates that included four closely related pathogen species originating from four continents. Pathogen isolates from bread and durum wheat were included to evaluate the emergence of specificity towards these hosts in M. graminicola. Although mitochondrial and nuclear genomes differed greatly in degree of genetic variability, their coalescence was remarkably congruent, supporting the proposed origin of M. graminicola through host tracking. The coalescence analysis was unable to trace M. graminicola host specificity through recent evolutionary time, indicating that the specificity towards durum or bread wheat emerged following the domestication of the pathogen on wheat.

Mitochondrial genome sequences and molecular evolution of the Irish potato famine pathogen, Phytophthora infestans

The mitochondrial genomes of haplotypes of the Irish potato famine pathogen, Phytophthora infestans, were sequenced. The genome sizes were 37,922, 39,870 and 39,840 bp for the type Ia, IIa and IIb mitochondrial DNA (mtDNA) haplotypes, respectively. The mitochondrial genome size for the type Ib haplotype, previously sequenced by others, was 37,957 bp. More than 90% of the genome contained coding regions. The GC content was 22.3%. A total of 18 genes involved in electron transport, 2 RNA-encoding genes, 16 ribosomal protein genes and 25 transfer RNA genes were coded on both strands with a conserved arrangement among the haplotypes. The type I haplotypes contained six unique open reading frames (ORFs) of unknown function while the type II haplotypes contained 13 ORFs of unknown function. Polymorphisms were observed in both coding and non-coding regions although the highest variation was in non-coding regions. The type I haplotypes (Ia and Ib) differed by only 14 polymorphic sites, whereas the type II haplotypes (IIa and IIb) differed by 50 polymorphic sites. The largest number (152) of polymorphic sites was found between the type IIb and Ia haplotypes. A large spacer flanked by the genes coding for tRNA-Tyr (trnY) and the small subunit RNA (rns) contained the largest number of polymorphic sites and corresponds to the region where a large indel that differentiates type II from type I haplotypes is located. The size of this region was 785, 2,666 and 2,670 bp in type Ia, IIa and IIb haplotypes, respectively. Among the four haplotypes, 81 mutations were identified. Phylogenetic and coalescent analysis revealed that although the type I and II haplotypes shared a common ancestor, they clearly formed two independent lineages that evolved independently. The type II haplotypes diverged earlier than the type I haplotypes. Thus our data do not support the previous hypothesis that the type II lineages evolved from the type I lineages. The type I haplotypes diverged more recently and the mutations associated with the evolution of the Ia and Ib types were identified.

Fluorescence-based directed termination PCR: direct mutation characterization without sequencing

We describe a fluorescence-based directed termination PCR (fluorescent DT-PCR) that allows accurate determination of actual sequence changes without dideoxy DNA sequencing. This is achieved using near infrared dye-labeled primers and performing two PCR reactions under low and unbalanced dNTP concentrations. Visualization of resulting termination fragments is accomplished with a dual dye Li-cor DNA sequencer. As each DT-PCR reaction generates two sets of terminating fragments, a pair of complementary reactions with limiting dATP and dCTP collectively provide information on the entire sequence of a target DNA, allowing an accurate determination of any base change. Blind analysis of 78 mutants of the supF reporter gene using fluorescent DT-PCR not only correctly determined the nature and position of all types of substitution mutations in the supF gene, but also allowed rapid scanning of the signature sequences among identical mutations. The method provides simplicity in the generation of terminating fragments and 100% accuracy in mutation characterization. Fluorescent DT-PCR was successfully used to generate a UV-induced spectrum of mutations in the supF gene following replication on a single plate of human DNA repair-deficient cells. We anticipate that the automated DT-PCR method will serve as a cost-effective alternative to dideoxy sequencing in studies involving large-scale analysis for nucleotide sequence changes.

Mutational spectrum analysis of RNase H(35) deficient Saccharomyces cerevisiae using fluorescence-based directed termination PCR

Mutational spectrum analysis has become an informative genetic tool to understand those protein functions involved in mutation avoidance pathways since specific types of mutations are often associated with particular protein defects involved in DNA replication and repair. In this study, we describe a novel, fluorescence-based procedure for direct determination of deletions and insertions with 100% accuracy. We performed two complementary directed termination PCR with near infrared dye-labeled primers, followed by visualization of termination fragments using an automated Li-cor DNA sequencer. This method is used for rapid analysis of mutational spectra generated in nuclease-defective strains of Saccharomyces cerevisiae to elucidate the role of RNase H(35) in RNA primer removal during DNA replication and in mutation avoidance. Strains deficient in RNH35 displayed a distinct spontaneous mutation spectrum of deletions characterized by a unique 4 bp deletion in a lys2-Bgl allele. This was in sharp contrast to strains deficient in rad27 that displayed duplication mutations. Further analysis of mutations in a rnh35/rad27 double mutant revealed a mixed spectrum. These results indicate that RNase H(35) may participate in a redundant pathway in Okazaki fragment processing and that mutational spectra caused by protein deficiencies may be more intermediate-specific than pathway-specific.

Terminal branch haplotype analysis: a novel approach to investigate newly arisen variants of mitochondrial DNA in natural populations

The discrimination of recent mutational derivatives from ancestral variation is a critical antecedent to any effort which aims to identify the factors modulating the rates of origin and persistence of new mutants. We propose that newly arisen mtDNA variants, which we designate as terminal branch haplotypes (TBHs), can be recognized by joint sequencing and phylogenetic analysis. This study examined mtDNA diversity in natural populations of the brown bullhead (Ameiurus nebulosus) from four heavily contaminated sites and three relatively pristine locations. While sequence analysis of the mtDNA D-loop region revealed that TBHs were prevalent in these populations, contaminant exposure appeared to play a minor role in their generation. Instead, most TBHs likely arose due to spontaneous mutations with variation in their incidence among sites reflecting the impact of demographic factors.