Aedes aegypti on Madeira Island (Portugal): genetic variation of a recently introduced dengue vector

The increasing population of Aedes aegypti mosquitoes on Madeira Island (Portugal) resulted in the first autochthonous dengue outbreak, which occurred in October 2012. Our study establishes the first genetic evaluation based on the mitochondrial DNA (mtDNA) genes [cytochrome oxidase subunit I (COI) and NADH dehydrogenase subunit 4 (ND4)] and knockdown resistance ( kdr ) mutations exploring the colonisation history and the genetic diversity of this insular vector population. We included mosquito populations from Brazil and Venezuela in the analysis as putative geographic sources. The Ae. aegyptipopulation from Madeira showed extremely low mtDNA genetic variability, with a single haplotype for COI and ND4. We also detected the presence of two important kdr mutations and the quasi-fixation of one of these mutations (F1534C). These results are consistent with a unique recent founder event that occurred on the island of Ae. aegyptimosquitoes that carry kdr mutations associated with insecticide resistance. Finally, we also report the presence of the F1534C kdr mutation in the Brazil and Venezuela populations. To our knowledge, this is the first time this mutation has been found in South American Ae. aegypti mosquitoes. Given the present risk of Ae. aegypti re-invading continental Europe from Madeira and the recent dengue outbreaks on the island, this information is important to plan surveillance and control measures.

Molecular characterization of fire ants, Solenopsis spp., from Brazil based on analysis of mtDNA gene cytochrome oxidase I

Species from the Solenopsis saevissima (Smith) (Hymenoptera: Formicidae) species group are native to South America and have a cosmopolitan distribution because they have been accidentally introduced in many countries around the world. In Brazil, they have a wide distribution, including urban areas. The present study was conducted to investigate the characterization of Solenopsis genus populations associated with urban/human interference sites in Brazil by analyzing the mitochondrial gene cytochrome oxidase I and estimating the degree of relatedness of these populations to make inferences about their phylogeny and also observe the patterns of mitochondrial haplotype (mitotype) distribution across their range. The results revealed complete geographical coherence and polyphyly for the Solenopsis invicta Buren and Solenopsis saevissima species groups, which confirms the diversity of the genera. It also suggests the possibility that reproductively-isolated populations occur, resulting in the evolutionary process of speciation. No predominant haplotype was found in the populations analyzed, but some were more prevalent.

Genetic variability and demographic history of Heliothis virescens (Lepidoptera: Noctuidae) populations from Brazil inferred by mtDNA sequences

Intra-and inter-population genetic variability and the demographic history of Heliothis virescens (F.) populations were evaluated by using mtDNA markers (coxI, coxII and nad6) with samples from the major cotton-and soybean-producing regions in Brazil in the growing seasons 2007/08, 2008/09 and 2009/10. AMOVA indicated low and non-significant genetic structure, regardless of geographical scale, growing season or crop, with most of genetic variation occurring within populations. Clustering analyzes also indicated low genetic differentiation. The haplotype network obtained with combined datasets resulted in 35 haplotypes, with 28 exclusive occurrences, four of them sampled only from soybean fields. The minimum spanning network showed star-shaped structures typical of populations that underwent a recent demographic expansion. The recent expansion was supported by other demographic analyzes, such as the Bayesian skyline plot, the unimodal distribution of paired differences among mitochondrial sequences, and negative and significant values of neutrality tests for the Tajima's D and Fu's F-S parameters. In addition, high values of haplotype diversity ((H) over cap) and low values of nucleotide diversity (pi), combined with a high number of low frequency haplotypes and values of theta(pi)<theta(W), suggested a recent demographic expansion of H. virescens populations in Brazil. This demographic event could be responsible for the low genetic structure currently found; however, haplotypes present uniquely at the same geographic regions and from one specific host plant suggest an initial differentiation among H. virescens populations within Brazil.

LOW GENETIC DIVERSITY IN Wolbachia-INFECTED Culex quinquefasciatus (DIPTERA: CULICIDAE) FROM BRAZIL AND ARGENTINA

Culex quinquefasciatus is a vector of human pathogens, including filarial nematodes and several viruses. Although its epidemiological relevance is known to vary across geographical regions, an understanding of its population genetic structure is still incipient. In light of this, we evaluated the genetic diversity of Cx. quinquefasciatus and Cx. pipiens x Cx. quinquefasciatus hybrids collected from nine localities in Brazil and one site in Argentina. We used mitochondrial genes cox1 and nd4, along with the coxA and wsp genes of the maternally-inherited Wolbachia endosymbiont. The nd4 fragment was invariant between samples, whilst cox1 exhibited four haplotypes that separated two types of Cx. quinquefasciatus, one clustered in southern Brazil. Low sequence diversity was generally observed, being discussed. Both Brazilian and Argentinian mosquitoes were infected with a single Wolbachia strain. As reported in previous studies with these populations, cox1 and nd4 diversity is not congruent with the population structure revealed by nuclear markers or alar morphology. Future Cx. quinquefasciatus research should, if possible, evaluate mtDNA diversity in light of other markers.

Low genetic diversity in Wolbachia-Infected Culex quinquefasciatus (Diptera: Culicidae) from Brazil and Argentina

Culex quinquefasciatus is a vector of human pathogens, including filarial nematodes and several viruses. Although its epidemiological relevance is known to vary across geographical regions, an understanding of its population genetic structure is still incipient. In light of this, we evaluated the genetic diversity of Cx. quinquefasciatus and Cx. pipiens x Cx. quinquefasciatus hybrids collected from nine localities in Brazil and one site in Argentina. We used mitochondrial genes cox1 and nd4, along with the coxA and wsp genes of the maternally-inherited Wolbachia endosymbiont. The nd4 fragment was invariant between samples, whilst cox1 exhibited four haplotypes that separated two types of Cx. quinquefasciatus, one clustered in southern Brazil. Low sequence diversity was generally observed, being discussed. Both Brazilian and Argentinian mosquitoes were infected with a single Wolbachia strain. As reported in previous studies with these populations, cox1 and nd4 diversity is not congruent with the population structure revealed by nuclear markers or alar morphology. Future Cx. quinquefasciatus research should, if possible, evaluate mtDNA diversity in light of other markers.

Molecular bases, pathogenic mechanisms and possible therapeutic approach in Leber's Hereditary Optic Neuropathy

Leber’s hereditary optic neuropathy (LHON) is a mitochondrial disease characterized by a rapid loss of central vision and optic atrophy, due to the selective degeneration of retinal ganglion cells. The age of onset is around 20, and the degenerative process is fast and usually the second eye becomes affected in weeks or months. Even if this pathology is well known and has been well characterized, there are still open questions on its pathophysiology, such as the male prevalence, the incomplete penetrance and the tissue selectivity. This maternally inherited disease is caused by mutations in mitochondrial encoded genes of NADH ubiquinone oxidoreductase (complex I) of the respiratory chain. The 90% of LHON cases are caused by one of the three common mitochondrial DNA mutations (11778/ND4, 14484/ND6 and 3460/ND1) and the remaining 10% is caused by rare pathogenic mutations, reported in literature in one or few families. Moreover, there is also a small subset of patients reported with new putative pathogenic nucleotide changes, which awaits to be confirmed. We here clarify some molecular aspects of LHON, mainly the incomplete penetrance and the role of rare mtDNA mutations or variants on LHON expression, and attempt a possible therapeutic approach using the cybrids cell model. We generated novel structural models for mitochondrial encoded complex I subunits and a conservation analysis and pathogenicity prediction have been carried out for LHON reported mutations. This in-silico approach allowed us to locate LHON pathogenic mutations in defined and conserved protein domains and can be a useful tool in the analysis of novel mtDNA variants with unclear pathogenic/functional role. Four rare LHON pathogenic mutations have been identified, confirming that the ND1 and ND6 genes are mutational hot spots for LHON. All mutations were previously described at least once and we validated their pathogenic role, suggesting the need for their screening in LHON diagnostic protocols. Two novel mtDNA variants with a possible pathogenic role have been also identified in two independent branches of a large pedigree. Functional studies are necessary to define their contribution to LHON in this family. It also been demonstrated that the combination of mtDNA rare polymorphic variants is relevant in determining the maternal recurrence of myoclonus in unrelated LHON pedigrees. Thus, we suggest that particular mtDNA backgrounds and /or the presence of specific rare mutations may increase the pathogenic potential of the primary LHON mutations, thereby giving rise to the extraocular clinical features characteristic of the LHON “plus” phenotype. We identified the first molecular parameter that clearly discriminates LHON affected individuals from asymptomatic carriers, the mtDNA copy number. This provides a valuable mechanism for future investigations on variable penetrance in LHON. However, the increased mtDNA content in LHON individuals was not correlated to the functional polymorphism G1444A of PGC-1 alpha, the master regulator of mitochondrial biogenesis, but may be due to gene expression of genes involved in this signaling pathway, such as PGC-1 alpha/beta and Tfam. Future studies will be necessary to identify the biochemical effects of rare pathogenic mutations and to validate the novel candidate mutations here described, in terms of cellular bioenergetic characterization of these variants. Moreover, we were not able to induce mitochondrial biogenesis in cybrids cell lines using bezafibrate. However, other cell line models are available, such as fibroblasts harboring LHON mutations, or other approaches can be used to trigger the mitochondrial biogenesis.

Patterns of prehistoric human mobility in Polynesia indicated by mtDNA from the Pacific rat

Human settlement of Polynesia was a major event in world prehistory. Despite the vastness of the distances covered, research suggests that prehistoric Polynesian populations maintained spheres of continuing interaction for at least some period of time in some regions. A low level of genetic variation in ancestral Polynesian populations, genetic admixture (both prehistoric and post-European contact), and severe population crashes resulting from introduction of European diseases make it difficult to trace prehistoric human mobility in the region by using only human genetic and morphological markers. We focus instead on an animal that accompanied the ancestral Polynesians on their voyages. DNA phylogenies derived from mitochondrial control-region sequences of Pacific rats (Rattus exulans) from east Polynesia are presented. A range of specific hypotheses regarding the degree of interaction within Polynesia are tested. These include the issues of multiple contacts between central east Polynesia and the geographically distinct archipelagos of New Zealand and Hawaii. Results are inconsistent with models of Pacific settlement involving substantial isolation after colonization and confirm the value of genetic studies on commensal species for elucidating the history of human settlement.

Apparent mtDNA heteroplasmy in Alzheimer’s disease patients and in normals due to PCR amplification of nucleus-embedded mtDNA pseudogenes

In an unprecedented finding, Davis et al. [Davis, R. E., Miller, S., Herrnstadt, C., Ghosh, S. S., Fahy, E., Shinobu, L. A., Galasko, D., Thal, L. J., Beal, M. F., Howell, N. & Parker, W. D., Jr. (1997) Proc. Natl. Acad. Sci. USA 94, 4526–4531] used an unusual DNA isolation method to show that healthy adults harbor a specific population of mutated mitochondrial cytochrome c oxidase (COX) genes that coexist with normal mtDNAs. They reported that this heteroplasmic population was present at a level of 10–15% in the blood of normal individuals and at a significantly higher level (20–30%) in patients with sporadic Alzheimer’s disease. We provide compelling evidence that the DNA isolation method employed resulted in the coamplification of authentic mtDNA-encoded COX genes together with highly similar COX-like sequences embedded in nuclear DNA (“mtDNA pseudogenes”). We conclude that the observed heteroplasmy is an artifact.

Ancient mtDNA sequences in the human nuclear genome: A potential source of errors in identifying pathogenic mutations

Nuclear-localized mtDNA pseudogenes might explain a recent report describing a heteroplasmic mtDNA molecule containing five linked missense mutations dispersed over the contiguous mtDNA CO1 and CO2 genes in Alzheimer’s disease (AD) patients. To test this hypothesis, we have used the PCR primers utilized in the original report to amplify CO1 and CO2 sequences from two independent ρ° (mtDNA-less) cell lines. CO1 and CO2 sequences amplified from both of the ρ° cells, demonstrating that these sequences are also present in the human nuclear DNA. The nuclear pseudogene CO1 and CO2 sequences were then tested for each of the five “AD” missense mutations by restriction endonuclease site variant assays. All five mutations were found in the nuclear CO1 and CO2 PCR products from ρ° cells, but none were found in the PCR products obtained from cells with normal mtDNA. Moreover, when the overlapping nuclear CO1 and CO2 PCR products were cloned and sequenced, all five missense mutations were found, as well as a linked synonymous mutation. Unlike the findings in the original report, an additional 32 base substitutions were found, including two in adjacent tRNAs and a two base pair deletion in the CO2 gene. Phylogenetic analysis of the nuclear CO1 and CO2 sequences revealed that they diverged from modern human mtDNAs early in hominid evolution about 770,000 years before present. These data would be consistent with the interpretation that the missense mutations proposed to cause AD may be the product of ancient mtDNA variants preserved as nuclear pseudogenes.

Muscle-specific mutations accumulate with aging in critical human mtDNA control sites for replication

The recently discovered aging-dependent large accumulation of point mutations in the human fibroblast mtDNA control region raised the question of their occurrence in postmitotic tissues. In the present work, analysis of biopsied or autopsied human skeletal muscle revealed the absence or only minimal presence of those mutations. By contrast, surprisingly, most of 26 individuals 53 to 92 years old, without a known history of neuromuscular disease, exhibited at mtDNA replication control sites in muscle an accumulation of two new point mutations, i.e., A189G and T408A, which were absent or marginally present in 19 individuals younger than 34 years. These two mutations were not found in fibroblasts from 22 subjects 64 to 101 years of age (T408A), or were present only in three subjects in very low amounts (A189G). Furthermore, in several older individuals exhibiting an accumulation in muscle of one or both of these mutations, they were nearly absent in other tissues, whereas the most frequent fibroblast-specific mutation (T414G) was present in skin, but not in muscle. Among eight additional individuals exhibiting partial denervation of their biopsied muscle, four subjects >80 years old had accumulated the two muscle-specific point mutations, which were, conversely, present at only very low levels in four subjects ≤40 years old. The striking tissue specificity of the muscle mtDNA mutations detected here and their mapping at critical sites for mtDNA replication strongly point to the involvement of a specific mutagenic machinery and to the functional relevance of these mutations.

Complex mtDNA constitutes an approximate 620-kb insertion on Arabidopsis thaliana chromosome 2: Implication of potential sequencing errors caused by large-unit repeats

Previously conducted sequence analysis of Arabidopsis thaliana (ecotype Columbia-0) reported an insertion of 270-kb mtDNA into the pericentric region on the short arm of chromosome 2. DNA fiber-based fluorescence in situ hybridization analyses reveal that the mtDNA insert is 618 ± 42 kb, ≈2.3 times greater than that determined by contig assembly and sequencing analysis. Portions of the mitochondrial genome previously believed to be absent were identified within the insert. Sections of the mtDNA are repeated throughout the insert. The cytological data illustrate that DNA contig assembly by using bacterial artificial chromosomes tends to produce a minimal clone path by skipping over duplicated regions, thereby resulting in sequencing errors. We demonstrate that fiber-fluorescence in situ hybridization is a powerful technique to analyze large repetitive regions in the higher eukaryotic genomes and is a valuable complement to ongoing large genome sequencing projects.

Fialuridine and its metabolites inhibit DNA polymerase gamma at sites of multiple adjacent analog incorporation, decrease mtDNA abundance, and cause mitochondrial structural defects in cultured hepatoblasts.

The thymidine analog fialuridine deoxy-2-fluoro-beta-D-arabinofuranosyl)-5-iodouracil (FIAU) was toxic in trials for chronic hepatitis B infection. One mechanism postulated that defective mtDNA replication was mediated through inhibition of DNA polymerase-gamma (DNA pol-gamma), by FIAU triphosphate (FIALTP) or by triphosphates of FIAU metabolites. Inhibition kinetics and primer-extension analyses determined biochemical mechanisms of FIAU, 1-(2-deoxy-2-fluoro-beta-D-arabinofuranosyl) -5-methyluracil (FAU), 1-(2-deoxy-2-fluoro-beta-D-arabinofuranosyl)uracil triphosphate (TP) inhibition of DNA pol-gamma. dTMP incorporation by DNA pol-gamma was inhibited competitively by FIAUTP, FMAUTP, and FAUTP (K1=0.015, 0.03, and 1.0 microM, respectively). By using oliginucleotide template-primers. DNA pol-gamma incorporated each analog into DNA opposite a single adenosine efficiently without effects on DNA chain elongation. Incorporation of multiple adjacent analogs at positions of consecutive adenosines dramatically impaired chain elongation by DNA pol-gamma. Effects of FIAU, FMAU, and FAU on HepG2 cell mmtDNA abundance and ultrastructure were determined. After 14 days, mtDNA decreased by 30% with 20 microM FIAU or 20 microM FMAU and decreased less than 10% with 100 microM FAU. FIAU and FMAU disrupted mitochondria and caused accumulation of intracytoplasmic lipid droplets. Biochemical and cell biological findings suggest that FIAU and its metabolites inhibit mtDNA replication, most likely at positions of adenosine tracts, leading to decreased mtDNA and mitochondrial ultrastructural defects.

Geographical structuring in the mtDNA of Italians.

Geographical patterns of mtDNA variation were studied in 12 Italian samples (1072 individuals) by two different spatial autocorrelation methods. Separate analyses of the frequencies of 12 restriction morphs show North-South clines, differences between Sardinia and the mainland populations, and the effects of isolation by distance. A recently developed autocorrelation statistic summarizing molecular similarity at all sites (AIDA; autocorrelation index for DNA analysis) confirms the presence of a clinical pattern; differences between random pairs of haplotypes tend to increase with their geographical distance. The partition of gene diversity, however, reveals that most variability occurs within populations, whereas differences between populations are minor (GST = 0.057). When the data from the 12 samples are pooled, two descriptors of genetic variability (number of polymorphic sites and average sequence difference between pairs of individuals) do not behave as expected under neutrality. The presence of clinal patterns, Tajima's tests, and a simulation experiment agree in suggesting that population sizes increased rapidly in Italy and Sicily but not necessarily so in Sardinia. The distribution of pairwise sequence differences in the Italian peninsula (excluding Sardinia) permits a tentative location of the demographic increase between 8000 and 20,500 years ago. These dates are consistent with archaeological estimates of two distinct expansion processes, occurring, respectively, in the Neolithic and after the last glacial maximum in the Paleolithic. Conversely, there is no genetic evidence that such processes have had a major impact on the Sardinian population.