Organelle DNA Research Articles

Implications of mutation of organelle genomes for organelle function and evolution.

Organelle genomes undergo more variation, including that resulting from damage, than eukaryotic nuclear genomes, or bacterial genomes, under the same conditions. Recent advances in characterizing the changes to genomes of chloroplasts and mitochondria of Zea mays should, when applied more widely, help our understanding of how damage to organelle genomes relates to how organelle function is maintained through the life of individuals and in succeeding generations. Understanding of the degree of variation in the changes to organelle DNA and its repair among photosynthetic organisms might help to explain the variations in the rate of nucleotide substitution among organelle genomes. Further studies of organelle DNA variation, including that due to damage and its repair might also help us to understand why the extent of DNA turnover in the organelles is so much greater than that in their bacterial (cyanobacteria for chloroplasts, proteobacteria for mitochondria) relatives with similar rates of production of DNA-damaging reactive oxygen species. Finally, from the available data, even the longest-lived organelle-encoded proteins, and the RNAs needed for their synthesis, are unlikely to maintain organelle function for much more than a week after the complete loss of organelle DNA.

Genetic diversity in Japanese aromatic rice (Oryza sativa L.) as revealed by nuclear and organelle DNA markers

Aromatic rice (Oryza sativa L.) cultivated in Japan is regionally differentiated by geographical distribution and characteristics. We aimed to characterize the lineage of Japanese aromatic rice using DNA markers. Based on analyses with nuclear SSR markers, we found that Japanese aromatic rice cultivars belong, with one exception, to japonica but showed some differences from authentic japonica and were divided into two clades that were distributed in western and eastern Japan, respectively. Further analyses with organelle markers showed that most of the cultivars in eastern Japan had cytoplasm characterized by tropical japonica, whereas most of those in western Japan had cytoplasm characterized by temperate japonica. We postulate that the ancestor of the cultivars in eastern Japan differs from those of the cultivars in western Japan, and that the two groups may have been separately introduced from Taiwan into Japan. The cytoplasm of aromatic rice cultivars in western Japan may have originated from tropical japonica and been substituted into the cytoplasm of temperate japonica through hybridization between tropical japonica as a male parent and temperate japonica as a female parent.

Mitochondrial and plastid genome architecture: Reoccurring themes, but significant differences at the extremes

Mitochondrial and plastid genomes show a wide array of architectures, varying immensely in size, structure, and content. Some organelle DNAs have even developed elaborate eccentricities, such as scrambled coding regions, nonstandard genetic codes, and convoluted modes of posttranscriptional modification and editing. Here, we compare and contrast the breadth of genomic complexity between mitochondrial and plastid chromosomes. Both organelle genomes have independently evolved many of the same features and taken on similar genomic embellishments, often within the same species or lineage. This trend is most likely because the nuclear-encoded proteins mediating these processes eventually leak from one organelle into the other, leading to a high likelihood of processes appearing in both compartments in parallel. However, the complexity and intensity of genomic embellishments are consistently more pronounced for mitochondria than for plastids, even when they are found in both compartments. We explore the evolutionary forces responsible for these patterns and argue that organelle DNA repair processes, mutation rates, and population genetic landscapes are all important factors leading to the observed convergence and divergence in organelle genome architecture.

Organelle DNA rearrangement mapping reveals U-turn-like inversions as a major source of genomic instability in Arabidopsis and humans.

Failure to maintain organelle genome stability has been linked to numerous phenotypes, including variegation and cytosolic male sterility (CMS) in plants, as well as cancer and neurodegenerative diseases in mammals. Here we describe a next-generation sequencing approach that precisely maps and characterizes organelle DNA rearrangements in a single genome-wide experiment. In addition to displaying global portraits of genomic instability, it surprisingly unveiled an abundance of short-range rearrangements in Arabidopsis thaliana and human organelles. Among these, short-range U-turn-like inversions reach 25% of total rearrangements in wild-type Arabidopsis plastids and 60% in human mitochondria. Furthermore, we show that replication stress correlates with the accumulation of this type of rearrangement, suggesting that U-turn-like rearrangements could be the outcome of a replication-dependent mechanism. We also show that U-turn-like rearrangements are mostly generated using microhomologies and are repressed in plastids by Whirly proteins WHY1 and WHY3. A synergistic interaction is also observed between the genes for the plastid DNA recombinase RECA1 and those encoding plastid Whirly proteins, and the triple mutant why1why3reca1 accumulates almost 60 times the WT levels of U-turn-like rearrangements. We thus propose that the process leading to U-turn-like rearrangements may constitute a RecA-independent mechanism to restart stalled forks. Our results reveal that short-range rearrangements, and especially U-turn-like rearrangements, are a major factor of genomic instability in organelles, and this raises the question of whether they could have been underestimated in diseases associated with mitochondrial dysfunction.

Association between Chloroplast and Mitochondrial DNA sequences in Chinese Prunus genotypes (Prunus persica, Prunus domestica, and Prunus avium)

BackgroundThe nuclear DNA is conventionally used to assess the diversity and relatedness among different species, but variations at the DNA genome level has also been used to study the relationship among different organisms. In most species, mitochondrial and chloroplast genomes are inherited maternally; therefore it is anticipated that organelle DNA remains completely associated. Many research studies were conducted simultaneously on organelle genome. The objectives of this study was to analyze the genetic relationship between chloroplast and mitochondrial DNA in three Chinese Prunus genotypes viz., Prunus persica, Prunus domestica, and Prunus avium.ResultsWe investigated the genetic diversity of Prunus genotypes using simple sequence repeat (SSR) markers relevant to the chloroplast and mitochondria. Most of the genotypes were genetically similar as revealed by phylogenetic analysis. The Y2 Wu Xing (Cherry) and L2 Hong Xin Li (Plum) genotypes have a high similarity index (0.89), followed by Zi Ye Li (0.85), whereas; L1 Tai Yang Li (plum) has the lowest genetic similarity (0.35). In case of cpSSR, Hong Tao (Peach) and L1 Tai Yang Li (Plum) genotypes demonstrated similarity index of 0.85 and Huang Tao has the lowest similarity index of 0.50. The mtSSR nucleotide sequence analysis revealed that each genotype has similar amplicon length (509 bp) except M5Y1 i.e., 505 bp with CCB256 primer; while in case of NAD6 primer, all genotypes showed different sizes. The MEHO (Peach), MEY1 (Cherry), MEL2 (Plum) and MEL1 (Plum) have 586 bps; while MEY2 (Cherry), MEZI (Plum) and MEHU (Peach) have 585, 584 and 566 bp, respectively. The CCB256 primer showed highly conserved sequences and minute single polymorphic nucleotides with no deletion or mutation. The cpSSR (ARCP511) microsatellites showed the harmonious amplicon length. The CZI (Plum), CHO (Peach) and CL1 (Plum) showed 182 bp; whileCHU (Peach), CY2 (Cherry), CL2 (Plum) and CY1 (Cherry) showed 181 bp amplicon lengths.ConclusionsThese results demonstrated high conservation in chloroplast and mitochondrial genome among Prunus species during the evolutionary process. These findings are valuable to study the organelle DNA diversity in different species and genotypes of Prunus to provide in depth insight in to the mitochondrial and chloroplast genomes.

Identification of Doris verrucosa mollusc via mitochondrial 16S rDNA

The biological and molecular knowledge of the marine life is of enormous importance in discovering new classes of natural products and in improving management and sustainable utilization of these useful genetic resources. Nonetheless, in this frame, genomic and organelle DNA isolation has been reported to be a limiting step, since downstream applications of molecular biology, such as restriction enzyme digestion, polymerase chain reaction (PCR) and DNA sequencing are hampered by the presence of compounds that are concurrently extracted. In this work we compared different DNA extraction techniques in three marine organisms and tested the downstream applications. One of the species utilized in this work was the Nudibranch Doris verrucosa L., a species with a well-known pharmaceutical interest. The phylogeny of this mollusc has long been discussed on the basis of morphological characters. As a result, D. verrucosa has often been misidentified with other sister species. Mitochondrial 16S rDNA sequence was chosen to better solve this concern and to accurately settle D. verrucosa into the right position. A partial sequence of the mt 16S rDNA of the Nudibranch D. verrucosa was obtained for the first time.The results here presented will provide knowledge useful to a better management and use of marine genetic resources, as well as in resolving taxonomic and identification issues currently open in these marine invertebrates.

A comparative analysis of insertional effects in genetically engineered plants: considerations for pre-market assessments.

During genetic engineering, DNA is inserted into a plant’s genome, and such insertions are often accompanied by the insertion of additional DNA, deletions and/or rearrangements. These genetic changes are collectively known as insertional effects, and they have the potential to give rise to unintended traits in plants. In addition, there are many other genetic changes that occur in plants both spontaneously and as a result of conventional breeding practices. Genetic changes similar to insertional effects occur in plants, namely as a result of the movement of transposable elements, the repair of double-strand breaks by non-homologous end-joining, and the intracellular transfer of organelle DNA. Based on this similarity, insertional effects should present a similar level of risk as these other genetic changes in plants, and it is within the context of these genetic changes that insertional effects must be considered. Increased familiarity with genetic engineering techniques and advances in molecular analysis techniques have provided us with a greater understanding of the nature and impact of genetic changes in plants, and this can be used to refine pre-market assessments of genetically engineered plants and food and feeds derived from genetically engineered plants.

Early and delayed long-term transcriptional changes and short-term transient responses during cold acclimation in olive leaves

Low temperature severely affects plant growth and development. To overcome this constraint, several plant species from regions having a cool season have evolved an adaptive response, called cold acclimation. We have studied this response in olive tree (Olea europaea L.) cv. Picual. Biochemical stress markers and cold-stress symptoms were detected after the first 24 h as sagging leaves. After 5 days, the plants were found to have completely recovered. Control and cold-stressed plants were sequenced by Illumina HiSeq 1000 paired-end technique. We also assembled a new olive transcriptome comprising 157,799 unigenes and found 6,309 unigenes differentially expressed in response to cold. Three types of response that led to cold acclimation were found: short-term transient response, early long-term response, and late long-term response. These subsets of unigenes were related to different biological processes. Early responses involved many cold-stress-responsive genes coding for, among many other things, C-repeat binding factor transcription factors, fatty acid desaturases, wax synthesis, and oligosaccharide metabolism. After long-term exposure to cold, a large proportion of gene down-regulation was found, including photosynthesis and plant growth genes. Up-regulated genes after long-term cold exposure were related to organelle fusion, nucleus organization, and DNA integration, including retrotransposons.

On the adaptive value of cytoplasmic genomes in plants.

Is DNA variation maintained in organelle genomes selectively neutral? The answer to this question has important implications for many aspects of ecology and evolution. While traditionally the answer has been 'yes', recent studies in animals have shown that, on the contrary, mitochondrial DNA polymorphism is frequently adaptive. In plants, however, the neutrality assumption has not been strongly challenged. Here, we begin with a critical evaluation of arguments in favour of this long-held view. We then discuss the latest empirical evidence for the opposing prediction that sequence variation in plant cytoplasmic genomes is frequently adaptive. While outstanding research progress is being made towards understanding this fundamental topic, we highlight the need for studies that combine information ranging from field experiments to physiology to molecular evolutionary biology. Such an interdisciplinary approach provides a means for determining the frequency, drivers and evolutionary significance of adaptive organelle DNA variation.

Remodeling of the AAA+ ATPase p97 by the UBX Adaptor Protein ASPL

The hexameric mammalian AAA+ ATPase p97, also known as VCP (valosin-containing protein; CDC48 in yeast), is a very abundant cytosolic protein and serves a wide variety of cellular functions. p97 is a central component in endoplasmic reticulum-associated degradation (ERAD) of proteins where it delivers ubiquitinated ERAD substrates to the proteasome. In addition, cellular roles of p97 in organelle membrane fusion, mitosis, DNA repair and suppression of apoptosis have been described. These different functions are linked to the binding of adaptor proteins to p97. Many of these adaptors contain ubiquitin regulatory X (UBX) domains. ASPL (alveolar soft part sarcoma locus, also known as TUG) was recenty identified as a p97 adaptor protein. As shown by crystal structure analysis, ASPL uses a substantially extended UBX domain for binding to the N domain of p97 where a lariat-like, mostly α-helical extension wraps around one subunit of p97. By this binding ASPL triggers the dissociation of functional p97 hexamers and the formation of p97:ASPL heterotetramers with 2:2 stoichiometry, leading to inactivation of the AAA+ ATPase. The p97-ASPL interaction in the heterotetramer is very tight, but p97 hexamer dissociation and heterotetramer formation may be suppressed by single-site mutations at p97-ASPL interfaces. p97 hexamer dissociation and p97-ASPL heterotetramer formation are linked to reduced ATPase activity of p97, cellular accumulation of ERAD substrates and apoptosis induction. To the best of our knowledge, this is the first time that the structural basis for adaptor protein-induced inactivation by hexamer dissociation of p97 and, indeed, any AAA+ ATPase has been demonstrated. This observation has far reaching implications for AAA+ ATPase-regulated processes.

Outcrossing rates and organelle inheritance estimated from two natural populations of the Japanese endemic conifer Sciadopitys verticillata

The Japanese endemic conifer Sciadopitys verticillata is one of the most phylogenetically isolated species of all plants. Occurring in small and scattered populations, the species is currently classified as Near Threatened by the International Union for Conservation of Nature and Natural Resources (IUCN) and as Vulnerable in three prefectures of Japan. This study investigated two major factors that should impact the genetic structure of the species at both the nuclear and organelle DNA level, the mating system and the inheritance of both the chloroplast and mitochondrial genomes. The mating system is crucial to determining the degree of outcrossing of plant species and thus should have a key role in shaping the species' population level genetic diversity and gene flow between populations but as yet has not been studied in S. verticillata. Nine mother trees and their seedling progeny from two natural populations were genotyped using genetic markers from three plant genomes (eight nuclear microsatellites and DNA sequence for the chloroplast and mitochondria). Using a maximum likelihood method implemented in the software MLTR, the study found an outcrossing rate in the seedling stage of 0.49 and 0.79 for Aburazaka and Mount Shirotori populations, respectively, and an average of 0.66 at the species level. These outcrossing rates were low for conifers and therefore may have potential deleterious implications for the conservation of the species. The test of organelle inheritance supported paternal transmission of both the chloroplast and mitochondria consistent with previous microscopic evidence.

Nucleases in higher plants and their possible involvement in DNA degradation during leaf senescence

During leaf senescence, macromolecules such as proteins and lipids are known to be degraded for redistribution into upper tissues. Similarly, nucleic acids appear to undergo fragmentation or degradation during senescence, but the physiological role of nucleic acid degradation, particularly of genomic DNA degradation, remains unclear. To date, more than a dozen of plant deoxyribonucleases have been reported, whereas it remains to be verified whether any of them degrade DNA during leaf senescence. This review summarizes current knowledge related to the plant nucleases that are induced developmentally or in a tissue-specific manner and are known to degrade DNA biochemically. Of these, several endonucleases (BFN1, CAN1, and CAN2) and an exonuclease (DPD1) in Arabidopsis seem to act in leaf senescence because they were shown to be inducible at the transcript level. This review specifically examines DPD1, which is dual-targeted to chloroplasts and mitochondria. Results show that, among the exonuclease family to which DPD1 belongs, DPD1 expression is extraordinary when estimated using a microarray database. DPD1 is the only example among the nucleases in which DNA degradation has been confirmed in vivo in pollen by mutant analysis. These data imply a significant role of organelle DNA degradation during leaf senescence and implicate DPD1 as a potential target for deciphering nucleotide salvage in plants.

English

In this study, multiple displacement amplification (MDA) was used for transgene detection, and proved to effectively improve the sensitivity of most fluorescence polymerase chain (PCR). Five events of genetically modified plants were chosen for detection in different plant species by MDA and real-time PCRs. Individual primer/probe systems based on each transgene target were used. The extracted DNA of each sample was diluted and used as a template for MDA. All extracted DNA and diluted mdaDNA were detected by real-time PCR. The results were evaluated with respect to DNA quantity, DNA purity and difference of Ct values of the isogenetic DNA extract before and after MDA (D-values). The MDA amplifying effects were largely different among genes and plants. Organelle DNA has more copies than genomic DNA by MDA. Key words: Multiple displacement amplification, genetically modified organisms, detection, events.

Evolutionary history of Purple cone spruce (Picea purpurea) in the Qinghai-Tibet Plateau: homoploid hybrid origin and Pleistocene expansion.

Hybridization and introgression can play an important role in speciation. Here, we examine their roles in the origin and evolution of Picea purpurea, a diploid spruce species occurring on the Qinghai-Tibet Plateau (QTP). Phylogenetic relationships and ecological differences between this species and its relatives, P.schrenkiana, P.likiangensis and P.wilsonii, are unclear. To clarify them, we surveyed sequence variation within and between them for 11 nuclear loci, three chloroplast (cp) and two mitochondrial (mt) DNA fragments, and examined their ecological requirements using ecological niche modelling. Initial analyses based on 11 nuclear loci rejected a close relationship between P.schrenkiana and P.purpurea. BP&P tests and ecological niche modelling indicated substantial divergence between the remaining three species and supported the species status of P.purpurea, which contained many private alleles as expected for a well-established species. Sequence variation for cpDNA and mtDNA suggested a close relationship between P.purpurea and P.wilsonii, while variation at the nuclear se1364 gene suggested P.purpurea was more closely related to P.likiangensis. Analyses of genetic divergence, Bayesian clustering and model comparison using approximate Bayesian computation (ABC) of nuclear (nr) DNA variation all supported the hypothesis that P.purpurea originated by homoploid hybrid speciation from P.wilsonii and P.likiangensis. The ABC analysis dated the origin of P.purpurea at the Pleistocene, and the estimated hybrid parameter indicated that 69% of its nuclear composition was contributed by P.likiangensis and 31% by P.wilsonii. Our results further suggested that during or immediately following its formation, P.purpurea was subject to organelle DNA introgression from P.wilsonii such that it came to possess both mtDNA and cpDNA of P.wilsonii. The estimated parameters indicated that following its origin, P.purpurea underwent an expansion during/after the largest Pleistocene glaciation recorded for the QTP.

Massive difference in synonymous substitution rates among mitochondrial, plastid, and nuclear genes of Phaeocystis algae

We are just beginning to understand how mutation rates differ among mitochondrial, plastid, and nuclear genomes. In most seed plants the mitochondrial mutation rate is estimated to be lower than those of the plastid and nucleus, whereas in the red alga Porphyra the opposite is true, and in certain green algae all three genomes appear to have similar rates of mutation. Relative rate statistics of organelle vs nuclear genes, however, are lacking for lineages that acquired their plastids through secondary endosymbiosis, but recent organelle DNA analyses suggest that they may differ drastically from what is observed in lineages with primary plastids, such as green plants and red algae. Here, by measuring synonymous nucleotide substitutions, we approximate the relative mutation rates within the haptophyte genus Phaeocystis, which has a red-algal-derived, secondary plastid. Synonymous-site divergence data indicate that for Phaeocystis antarctica and P. globosa the mitochondrial mutation rate is 10 and 3 times that of the plastid and nucleus, respectively. This differs drastically from relative rate estimates for primary-plastid-bearing lineages and presents a much more dynamic view of organelle vs nuclear mutation rates across the eukaryotic domain.

Proteomic Analysis of Eggs from Mytilus edulis Females Differing in Mitochondrial DNA Transmission Mode

Many bivalves have an unusual mechanism of mitochondrial DNA (mtDNA) inheritance called doubly uniparental inheritance (DUI) in which distinctly different genomes are inherited through the female (F genome) and male (M genome) lineages. In fertilized eggs that will develop into male embryos, the sperm mitochondria remain in an aggregation, which is believed to be delivered to the primordial germ cells and passed to the next generation through the sperm. In fertilized eggs that will develop into female embryos, the sperm mitochondria are dispersed throughout the developing embryo and make little if any contribution to the next generation. The frequency of embryos with the aggregated or dispersed mitochondrial type varies among females. Previous models of DUI have predicted that maternal nuclear factors cause molecular differences among unfertilized eggs from females producing embryos with predominantly dispersed or aggregated mitochondria. We test this hypothesis using females of each of the two types from a natural population. We have found small, yet detectable, differences of the predicted type at the proteome level. We also provide evidence that eggs of females giving the dispersed pattern have consistently lower expression for different proteasome subunits than eggs of females giving the aggregated pattern. These results, combined with those of an earlier study in which we used hatchery lines of Mytilus, and with a transcriptomic study in a clam that has the DUI system of mtDNA transmission, reinforce the hypothesis that the ubiquitin-proteasome system plays a key role in the mechanism of DUI and sex determination in bivalves. We also report that eggs of females giving the dispersed pattern have higher expression for arginine kinase and enolase, enzymes involved in energy production, whereas ferritin, which is involved in iron homeostasis, has lower expression. We discuss these results in the context of genetic models for DUI and suggest experimental methods for further understanding the role of these proteins in DUI.

Patterns of Genomic Integration of Nuclear Chloroplast DNA Fragments in Plant Species

The transfer of organelle DNA fragments to the nuclear genome is frequently observed in eukaryotes. These transfers are thought to play an important role in gene and genome evolution of eukaryotes. In plants, such transfers occur from plastid to nuclear [nuclear plastid DNAs (NUPTs)] and mitochondrial to nuclear (nuclear mitochondrial DNAs) genomes. The amount and genomic organization of organelle DNA fragments have been studied in model plant species, such as Arabidopsis thaliana and rice. At present, publicly available genomic data can be used to conduct such studies in non-model plants. In this study, we analysed the amount and genomic organization of NUPTs in 17 plant species for which genome sequences are available. The amount and distribution of NUPTs varied among the species. We also estimated the distribution of NUPTs according to the time of integration (relative age) by conducting sequence similarity analysis between NUPTs and the plastid genome. The age distributions suggested that the present genomic constitutions of NUPTs could be explained by the combination of the rapidly eliminated deleterious parts and few but constantly existing less deleterious parts.

Rapid and inexpensive NaOH based direct PCR for amplification of nuclear and organelle DNA from ramie (Boehmeria nivea), a bast fibre crop containing complex polysaccharides

Most of the PCR based approaches in plant science rely on lengthy and expensive DNA isolation protocols, which often involve use of hazardous chemicals. Direct PCR methods save time and cost of experiments and also increase efficiency of PCR. We have compared three rapid DNA extraction processes for direct PCR in ramie (Boehmeria nivea), a fibre crop species containing high amount of gummy complex polysaccharides and developed modified protocols for direct PCR using NaOH as extraction buffer and Tris/Tris–HCl/Tris–EDTA as dilution buffer. These protocols were successful in amplification of nuclear DNA from leaf and stem tissues using ISSR and SSR markers and also chloroplast DNA amplification using primers for rbcL gene. Our results also show that nature and quantity of dilution buffer are important for increasing efficiency of direct PCR. The NaOH based methods are simpler, cheaper and economical compared to other direct PCR methods and work very well for tissues containing high amount of complex polysaccharides. The protocols are suitable for batch processing and high throughput genotyping within a short time period, which will have many applications in plant genomic researches.

Palindromic Genes in the Linear Mitochondrial Genome of the Nonphotosynthetic Green Alga Polytomella magna

Organelle DNA is no stranger to palindromic repeats. But never has a mitochondrial or plastid genome been described in which every coding region is part of a distinct palindromic unit. While sequencing the mitochondrial DNA of the nonphotosynthetic green alga Polytomella magna, we uncovered precisely this type of genic arrangement. The P. magna mitochondrial genome is linear and made up entirely of palindromes, each containing 1–7 unique coding regions. Consequently, every gene in the genome is duplicated and in an inverted orientation relative to its partner. And when these palindromic genes are folded into putative stem-loops, their predicted translational start sites are often positioned in the apex of the loop. Gel electrophoresis results support the linear, 28-kb monomeric conformation of the P. magna mitochondrial genome. Analyses of other Polytomella taxa suggest that palindromic mitochondrial genes were present in the ancestor of the Polytomella lineage and lost or retained to various degrees in extant species. The possible origins and consequences of this bizarre genomic architecture are discussed.

DNA marker based identification of spontaneous hybrids between Pinus mugo and P. sylvestris at the Lithuanian sea‐side

The objective was to identify spontaneous hybrids between P. mugo and P. sylvestris using organelle DNA markers in sympatric zones at the sea‐side spit of Kursiu Nerija in western Lithuania. A field inventory was carried out over the entire Lithuanian part of the spit and 203 individuals morphologically intermediate between P. sylvestris and P. mugo were tested for their male parent with chloroplast DNA PCR‐RFLP markers and for their female parent with mitochondrial DNA PCR markers. Unfortunately, the mitochondrial DNA Nad7 marker failed to identify the female parent species in our study. However, the chloroplast DNA PCR‐RFLP marker revealed that out of 203 tested individuals only 23 had a paternity different than that indicated by morphology. Of these, 13 individuals were morphologically identified as P. sylvestris but possessed cpDNA of P. mugo (putative hybrids with P. sylvestris (female) ×P. mugo (male parent), and 10 individuals morphologically identified as P. mugo possessed cpDNA of P. sylvestris and may be hybrids with P. mugo (female) ×P. sylvestris (male parent). The remainder of the 177 individuals identified in the field inventory were were considered as pure species. In conclusion, our study indicates ongoing spontaneous hybridization between P. mugo and P. syvestris in Kursiu Nerija. Human impact via seed transfer on altered hybridization rates as well as evolutionary consequences are discussed.