Discordant evolution of organellar genomes in peas (Pisum L.).

Plastids and mitochondria have their own small genomes, which do not undergo meiotic recombination and may have evolutionary fates different from each other and that of the nuclear genome. For the first time, we sequenced mitochondrial genomes of pea (Pisum L.) from 42 accessions mostly representing diverse wild germplasm from throughout the wild pea geographical range. Six structural types of the pea mitochondrial genome were revealed. From the same accessions, plastid genomes were sequenced. Phylogenetic trees based on the plastid and mitochondrial genomes were compared. The topologies of these trees were highly discordant, implying not less than six events of hybridisation between diverged wild peas in the past, with plastids and mitochondria differently inherited by the descendants. Such discordant inheritance of organelles could have been driven by plastid-nuclear incompatibility, which is known to be widespread in crosses involving wild peas and affects organellar inheritance. The topology of the phylogenetic tree based on nucleotide sequences of a nuclear gene, His5, encoding a histone H1 subtype, corresponded to the current taxonomy and resembled that based on the plastid genome. Wild peas (Pisum sativum subsp. elatius s.l.) inhabiting Southern Europe were shown to be of hybrid origin, resulting from crosses of peas related to those presently inhabiting the eastern Mediterranean in a broad sense. These results highlight the roles of hybridisation and cytonuclear conflict in shaping plant microevolution.

Cryptic divergences in the genus Pisum L. (peas), as revealed by phylogenetic analysis of plastid genomes.

Organellar genomes may shed light on complicated patterns of plant evolution at inter- and intraspecies level. Primary structure of plastid genomes sequenced in this study and taken from public databases was characterised and compared in 22 diverse, mostly wild representatives of the genus Pisum (peas). Phylogenetic trees reconstructed via Bayesian approach on the basis of entire plastid genomes resembled those reconstructed on the basis of a nuclear gene His5 coding for a minor histone H1 subtype. They reveal Pisum fulvum as an early divergence of the genus but do not support other taxonomical subdivisions. The positions of three accessions, classified as P. sativum subsp. elatius (the wild subspecies of the common pea), appeared quite unexpected. On the entire plastid genome tree, two accessions, from the Black Sea area of Turkey and Georgia, clustered with representatives of another species, P. fulvum, while the other, from Greece, was the first divergence of the P. sativum branch. We suppose these unusual plastid genomes to be ancient lineages ascending to a 'missing link' between P. fulvum and P. sativum, represented by accession Pe 013 from Turkey. Accessions with common pea appearance but deeply diverged plastids could occur through occasional crossing of diverged pea lines in the past and biparental plastid inheritance, both events being possible in peas.

Nuclear-cytoplasmic conflict in pea (Pisum sativum L.) is associated with nuclear and plastidic candidate genes encoding acetyl-CoA carboxylase subunits.

In crosses of wild and cultivated peas (Pisum sativum L.), nuclear-cytoplasmic incompatibility frequently occurs manifested as decreased pollen fertility, male gametophyte lethality, sporophyte lethality. High-throughput sequencing of plastid genomes of one cultivated and four wild pea accessions differing in cross-compatibility was performed. Candidate genes for involvement in the nuclear-plastid conflict were searched in the reconstructed plastid genomes. In the annotated Medicago truncatula genome, nuclear candidate genes were searched in the portion syntenic to the pea chromosome region known to harbor a locus involved in the conflict. In the plastid genomes, a substantial variability of the accD locus represented by nucleotide substitutions and indels was found to correspond to the pattern of cross-compatibility among the accessions analyzed. Amino acid substitutions in the polypeptides encoded by the alleles of a nuclear locus, designated as Bccp3, with a complementary function to accD, fitted the compatibility pattern. The accD locus in the plastid genome encoding beta subunit of the carboxyltransferase of acetyl-coA carboxylase and the nuclear locus Bccp3 encoding biotin carboxyl carrier protein of the same multi-subunit enzyme were nominated as candidate genes for main contribution to nuclear-cytoplasmic incompatibility in peas. Existence of another nuclear locus involved in the accD-mediated conflict is hypothesized.

Divergence and population traits in evolution of the genus Pisum L. as reconstructed using genes of two histone H1 subtypes showing different phylogenetic resolution.

Two histone H1 subtype genes, His7 and His5, were sequenced in a set of 56 pea accessions. Phylogenetic reconstruction based on concatenated His5 and His7 sequences had three main clades. First clade corresponded to Pisum fulvum, the next divergence separated a clade inside Pisum sativum in the broad sense that did not correspond strictly to any proposed taxonomical subdivisions. According to our estimations, the earliest divergence separating P. fulvum occurred 1.7±0.4MYA. The other divergence with high bootstrap support that separated two P. sativum groups took place approximately 1.3±0.3MYA. Thus, the main divergences in the genus took place either in late Pliocene or in early Pleistocene, the time of onset of the profound climate cooling in the northern hemisphere. The ω=K(a)/K(s) ratio was 2.5 times higher for His5 sequences than for His7. Thus, His7 gene, coding for a unique subtype specific for actively growing tissues, might have evolved under stricter evolutionary constraints than His5, that codes for a minor H1 subtype with less specific expression pattern. For this reason phylogenetic reconstructions separately obtained from His5 sequences resolved tree topology much better than those obtained from His7 sequences. Computational estimation of population dynamic parameters in the genus Pisum L. from His5-His7 sequences using IMa2 software revealed a decrease of effective population size on the early stage of Pisum evolution.