Search

In plants, cellular function is orchestrated by three distinct genomes located within the nucleus, mitochondrion, and plastid. These genomes are interdependent, requiring tightly coordinated maintenance and expression. Plastids host several multisubunit protein complexes encoded by both the plastid and nuclear genomes. To investigate plastid-nuclear coevolution, this study focused on plastid ribosomal protein genes that are encoded by both plastid and nuclear genomes from 50 taxa across 15 of the 22 early branching major clades of the legume subfamily Papilionoideae. Comparative analysis of substitution rates was conducted across five gene sets: nuclear-encoded plastid-targeted ribosomal protein genes (NuCpRP), nuclear-encoded cytosol-targeted ribosomal genes (NuCyRP), other nuclear-encoded plastid-targeted genes that are not involved in ribosomes (NuCpOT), plastid-encoded ribosomal protein genes (CpRP) and plastid-encoded photosynthesis genes (CpPS). 1 Elevated nonsynonymous substitution rates (d N ) and ratios of nonsynonymous to synonymous substitution rates (d N /d S ; ω) were observed in both CpRP and NuCpRP compared to the other gene sets. Significant differences in d N for CpRP and NuCpRP were found between the papilionoid 50-kb inversion clade and other legumes. Using coevolution statistics and evolutionary rate covariation, strong signals of cytonuclear coevolution were identified, where nonsynonymous substitutions in CpRP and NuCpRP genes co-occur along the same branches of the Papilionoideae phylogeny. Increased ω in a few CpRP genes was due to intensified positive selection whereas most of the CpRP and NuCpRP increased ω was caused by relaxed purifying selection. This pattern not only underscores the role of cytonuclear incompatibility in driving speciation but also highlights its constraints on the genetic enhancement of papilionoid crop species., (Copyright © 2024 Elsevier Inc. All rights reserved.)

Genisteae, a tribe in the subfamily Papilionoideae (Fabaceae), is characterized by the production of quinolizidine alkaloids that confer pest resistance in most of its members. Many relationships at the generic level remain unresolved due largely to a lack of modern attempts to reconstruct the phylogeny. Previous studies with limited taxon sampling and only a few molecular loci indicated the presence of three clades within the tribe: the Lupinus clade, the Cytisus-Genista complex and the Argyrolobium group. There are also two historical genera, Teline and Chamaecytisus, that have been reclassified over the years with some controversy. Species from Teline are currently classified in Genista, and Chamaecytisus species are placed inside of Cytisus. Sellocharis is another genus with vague placement inside of Genisteae near Anarthrophyllum, based mostly on morphology and cytology. Representative taxa from 24 of the 25 genera of Genisteae, along with species of historical genera Teline and Chamaecytisus, were sampled and utilized in a phylogenomic investigation using both plastid and nuclear data to resolve relationships at the generic level. Whole genomic DNA was sequenced and complete plastomes were assembled and annotated. Low-copy nuclear genes were retrieved from the genomic DNA sequences using a mapping-assembly-scaffold approach. Phylogenetic analyses using maximum likelihood, Bayesian and coalescence methods resulted in fully resolved and strongly supported trees for both nuclear and plastid data that show four major clades inside of Genisteae: Cytisus-Genista complex, Lupinus clade, Argyrolobium group and the novel Anarthrophyllum group. The resulting phylogenetic trees also supported the transfer of the Argyrolobium group from Crotalarieae to Genisteae, and the reclassification of Teline inside of Genista, both of which were previously suggested in literature. The phylogenetic trees also determined a placement for Sellocharis in the new Anarthrophyllum group. Although both nuclear and plastid trees were congruent with respect to the monophyly and relationships among the four major clades, incongruence was detected within some of the major clades and the potential causes are discussed., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Elsevier Inc. All rights reserved.)

The water-lily clade represents the second earliest-diverging branch of angiosperms. Most of its species belong to Nymphaeaceae, of which the "core Nymphaeaceae"-comprising the genera Euryale, Nymphaea and Victoria-is the most diverse clade. Despite previous molecular phylogenetic studies on the core Nymphaeaceae, various aspects of their evolutionary relationships have remained unresolved. The length-variable introns and intergenic spacers are known to contain most of the sequence variability within the water-lily plastomes. Despite the challenges with multiple sequence alignment, any new molecular phylogenetic investigation on the core Nymphaeaceae should focus on these noncoding plastome regions. For example, a new plastid phylogenomic study on the core Nymphaeaceae should generate DNA sequence alignments of all plastid introns and intergenic spacers based on the principle of conserved sequence motifs. In this investigation, we revisit the phylogenetic history of the core Nymphaeaceae by employing such an approach. Specifically, we use a plastid phylogenomic analysis strategy in which all coding and noncoding partitions are separated and then undergo software-driven DNA sequence alignment, followed by a motif-based alignment inspection and adjustment. This approach allows us to increase the reliability of the character base compared to the default practice of aligning complete plastomes through software algorithms alone. Our approach produces significantly different phylogenetic tree reconstructions for several of the plastome regions under study. The results of these reconstructions underscore that Nymphaea is paraphyletic in its current circumscription, that each of the five subgenera of Nymphaea is monophyletic, and that the subgenus Nymphaea is sister to all other subgenera of Nymphaea. Our results also clarify many evolutionary relationships within the Nymphaea subgenera Brachyceras, Hydrocallis and Nymphaea. In closing, we discuss whether the phylogenetic reconstructions obtained through our motif-based alignment adjustments are in line with morphological evidence on water-lily evolution., (© 2024 The Authors. Cladistics published by John Wiley & Sons Ltd on behalf of Willi Hennig Society.)

The North American Thermopsideae (Fabaceae: Papilionoideae), a monophyletic group comprising the North American endemic genus Baptisia, and the paraphyletic Eurasian-North American disjunct Thermopsis, is nested within the tribe Sophoreae. Previous phylogenetic studies have identified two East Asian taxa within the North American Thermopsideae, suggesting two independent dispersal events between North America-East Asia. More recent studies have also placed a third taxon, Vuralia turcica, an endemic species from Turkey, among the North American Thermopsideae. The presence of three geographically distant Eurasian taxa within a relatively young clade of North American origin is unprecedented among papilionoid legumes, and the biogeographic implications of this observation are not clear. To investigate this matter, 1540 low-copy nuclear genes and complete plastomes were obtained from 36 taxa across the core genistoids, including 26 newly sequenced taxa. Nuclear and plastome based maximum likelihood (ML) and ASTRAL analyses were conducted based on varying degrees of taxon coverage and read mapping consensus threshold values. Additional analyses were performed to estimate divergence times and to reconstruct biogeographic history. The results strongly support a previously undetected Old World clade, presently composed of V. turcica and T. chinensis, which diverged from the ancestor of the North American lineage during the mid to late Miocene. A single and recent North America-East Asia dispersal involving T. lupinoides is reported. Furthermore, the traditional inclusion of the genus Ammopiptanthus among Thermopsideae is not supported, and the monotypic generic status of Vuralia is called into question. A relatively high degree of cytonuclear discordance is reported within each sub-clade of the North American Thermopsideae. This finding is likely attributable to the high degree of interspecific hybridization reported within these groups and raises the need for more rigorous genome-scale testing to better delimit species within each of the reticulating subclades. Subjects: Biodiversity, Biogeography, Evolutionary Studies, Genetics, Plant Science., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024. Published by Elsevier Inc.)

Nuclear and organellar genomes can evolve at vastly different rates despite occupying the same cell. In most bilaterian animals, mitochondrial DNA (mtDNA) evolves faster than nuclear DNA, whereas this trend is generally reversed in plants. However, in some exceptional angiosperm clades, mtDNA substitution rates have increased up to 5,000-fold compared with closely related lineages. The mechanisms responsible for this acceleration are generally unknown. Because plants rely on homologous recombination to repair mtDNA damage, we hypothesized that mtDNA copy numbers may predict evolutionary rates, as lower copy numbers may provide fewer templates for such repair mechanisms. In support of this hypothesis, we found that copy number explains 47% of the variation in synonymous substitution rates of mtDNA across 60 diverse seed plant species representing ~300 million years of evolution. Copy number was also negatively correlated with mitogenome size, which may be a cause or consequence of mutation rate variation. Both relationships were unique to mtDNA and not observed in plastid DNA. These results suggest that homologous recombinational repair plays a role in driving mtDNA substitution rates in plants and may explain variation in mtDNA evolution more broadly across eukaryotes. Our findings also contribute to broader questions about the relationships between mutation rates, genome size, selection efficiency, and the drift-barrier hypothesis., Competing Interests: Competing interests statement:The authors declare no competing interest.

Although plastid genome (plastome) structure is highly conserved across most seed plants, investigations during the past two decades have revealed several disparately related lineages that experienced substantial rearrangements. Most plastomes contain a large inverted repeat and two single-copy regions, and a few dispersed repeats; however, the plastomes of some taxa harbour long repeat sequences (>300 bp). These long repeats make it challenging to assemble complete plastomes using short-read data, leading to misassemblies and consensus sequences with spurious rearrangements. Single-molecule, long-read sequencing has the potential to overcome these challenges, yet there is no consensus on the most effective method for accurately assembling plastomes using long-read data. We generated a pipeline, plastid Genome Assembly Using Long-read data (ptGAUL), to address the problem of plastome assembly using long-read data from Oxford Nanopore Technologies (ONT) or Pacific Biosciences platforms. We demonstrated the efficacy of the ptGAUL pipeline using 16 published long-read data sets. We showed that ptGAUL quickly produces accurate and unbiased assemblies using only ~50× coverage of plastome data. Additionally, we deployed ptGAUL to assemble four new Juncus (Juncaceae) plastomes using ONT long reads. Our results revealed many long repeats and rearrangements in Juncus plastomes compared with basal lineages of Poales. The ptGAUL pipeline is available on GitHub: https://github.com/Bean061/ptgaul., (© 2023 John Wiley & Sons Ltd.)

Cyperaceae, the second largest family in the monocot order Poales, comprises >5500 species and includes the genus Eleocharis with ∼ 250 species. A previous study of complete plastomes of two Eleocharis species documented extensive structural heteroplasmy, gene order changes, high frequency of dispersed repeats along with gene losses and duplications. To better understand the phylogenetic distribution of gene and intron content as well as rates and patterns of sequence evolution within and between mitochondrial and plastid genomes of Eleocharis and Cyperaceae, an additional 29 Eleocharis organelle genomes were sequenced and analyzed. Eleocharis experienced extensive gene loss in both genomes while loss of introns was mitochondria-specific. Eleocharis has higher rates of synonymous (dS) and nonsynonymous (dN) substitutions in the plastid and mitochondrion than most sampled angiosperms, and the pattern was distinct from other eudicot lineages with accelerated rates. Several clades showed higher dS and dN in mitochondrial genes than in plastid genes. Furthermore, nucleotide substitution rates of mitochondrial genes were significantly accelerated on the branch leading to Cyperaceae compared to most angiosperms. Mitochondrial genes of Cyperaceae exhibited dramatic loss of RNA editing sites and a negative correlation between RNA editing and dS values was detected among angiosperms. Mutagenic retroprocessing and dysfunction of DNA replication, repair and recombination genes are the most likely cause of striking rate accelerations and loss of edit sites and introns in Eleocharis and Cyperaceae organelle genomes., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 Elsevier Inc. All rights reserved.)

The tomato ( Solanum lycopersicum L.) is considered one of the most important vegetable crops globally, both agronomically and economically; however, its fruit development regulation network is still unclear. The transcription factors serve as master regulators, activating many genes and/or metabolic pathways throughout the entire plant life cycle. In this study, we identified the transcription factors that are coordinated with TCP gene family regulation in early fruit development by making use of the high-throughput sequencing of RNA (RNAseq) technique. A total of 23 TCP -encoding genes were found to be regulated at various stages during the growth of the fruit. The expression patterns of five TCPs were consistent with those of other transcription factors and genes. There are two unique subgroups of this larger family: class I and class II TCPs. Others were directly associated with the growth and/or ripening of fruit, while others were involved in the production of the hormone auxin. Moreover, it was discovered that TCP18 had an expression pattern that was similar to that of the ethylene-responsive transcription factor 4 (ERF4). Tomato fruit set and overall development are under the direction of a gene called auxin response factor 5 (ARF5). TCP15 revealed an expression that was in sync with this gene. This study provides insight into the potential processes that help in acquiring superior fruit qualities by accelerating fruit growth and ripening.

Medicago truncatula is a model legume that has been extensively investigated in diverse subdisciplines of plant science. Medicago littoralis can interbreed with M. truncatula and M. italica; these three closely related species form a clade, i.e. TLI clade. Genetic studies have indicated that M. truncatula accessions are heterogeneous but their taxonomic identities have not been verified. To elucidate the phylogenetic position of diverse M. truncatula accessions within the genus, we assembled 54 plastid genomes (plastomes) using publicly available next-generation sequencing data and conducted phylogenetic analyses using maximum likelihood. Five accessions showed high levels of plastid DNA polymorphism. Three of these highly polymorphic accessions contained sequences from both M. truncatula and M. littoralis. Phylogenetic analyses of sequences placed some accessions closer to distantly related species suggesting misidentification of source material. Most accessions were placed within the TLI clade and maximally supported the interrelationships of three subclades. Two Medicago accessions were placed within a M. italica subclade of the TLI clade. Plastomes with a 45-kb (rpl20-ycf1) inversion were placed within the M. littoralis subclade. Our results suggest that the M. truncatula accession genome pool represents more than one species due to possible mistaken identities and gene flow among closely related species., (© 2022. The Author(s).)

Salt stress negatively impacts crop production worldwide. Genetic diversity among barley ( Hordeum vulgare ) landraces adapted to adverse conditions should provide a valuable reservoir of tolerance genes for breeding programs. To identify molecular and biochemical differences between barley genotypes, transcriptomic and antioxidant enzyme profiles along with several morpho-physiological features were compared between salt-tolerant (Boulifa) and salt-sensitive (Testour) genotypes subjected to salt stress. Decreases in biomass, photosynthetic parameters, and relative water content were low in Boulifa compared to Testour. Boulifa had better antioxidant protection against salt stress than Testour, with greater antioxidant enzymes activities including catalase, superoxide dismutase, and guaiacol peroxidase. Transcriptome assembly for both genotypes revealed greater accumulation of differentially expressed transcripts in Testour compared to Boulifa, emphasizing the elevated transcriptional response in Testour following salt exposure. Various salt-responsive genes, including the antioxidant catalase 3, the osmoprotectant betaine aldehyde dehydrogenase 2, and the transcription factors MYB20 and MYB41 , were induced only in Boulifa. By contrast, several genes associated with photosystems I and II, and light receptor chlorophylls A and B, were more repressed in Testour. Co-expression network analysis identified specific gene modules correlating with differences in genotypes and morpho-physiological traits. Overall, salinity-induced differential transcript accumulation underlies the differential morpho-physiological response in both genotypes and could be important for breeding salt tolerance in barley.

Plant nuclear genomes harbor sequence elements derived from the organelles (mitochondrion and plastid) through intracellular gene transfer (IGT). Nuclear genomes also show a dramatic range of repeat content, suggesting that any sequence can be readily amplified. These two aspects of plant nuclear genomes are well recognized but have rarely been linked. Through investigation of 31 Medicago taxa we detected exceptionally high post-IGT amplification of mitochondrial (mt) DNA sequences containing rps10 in the nuclear genome of Medicago polymorpha and closely related species. The amplified sequences were characterized as tandem arrays of five distinct repeat motifs (2157, 1064, 987, 971, and 587 bp) that have diverged from the mt genome (mitogenome) in the M. polymorpha nuclear genome. The mt rps10-like arrays were identified in seven loci (six intergenic and one telomeric) of the nuclear chromosome assemblies and were the most abundant tandem repeat family, representing 1.6-3.0% of total genomic DNA, a value approximately three-fold greater than the entire mitogenome in M. polymorpha. Compared to a typical mt gene, the mt rps10-like sequence coverage level was 691.5-7198-fold higher in M. polymorpha and closely related species. In addition to the post-IGT amplification, our analysis identified the canonical telomeric repeat and the species-specific satellite arrays that are likely attributable to an ancestral chromosomal fusion in M. polymorpha. A possible relationship between chromosomal instability and the mt rps10-like tandem repeat family in the M. polymorpha clade is discussed., (© 2022 Society for Experimental Biology and John Wiley & Sons Ltd.)

Comprising 501 genera and around 14,000 species, Papilionoideae is not only the largest subfamily of Fabaceae (Leguminosae; legumes), but also one of the most extraordinarily diverse clades among angiosperms. Papilionoids are a major source of food and forage, are ecologically successful in all major biomes, and display dramatic variation in both floral architecture and plastid genome (plastome) structure. Plastid DNA-based phylogenetic analyses have greatly improved our understanding of relationships among the major groups of Papilionoideae, yet the backbone of the subfamily phylogeny remains unresolved. In this study, we sequenced and assembled 39 new plastomes that are covering key genera representing the morphological diversity in the subfamily. From 244 total taxa, we produced eight datasets for maximum likelihood (ML) analyses based on entire plastomes and/or concatenated sequences of 77 protein-coding sequences (CDS) and two datasets for multispecies coalescent (MSC) analyses based on individual gene trees. We additionally produced a combined nucleotide dataset comprising CDS plus matK gene sequences only, in which most papilionoid genera were sampled. A ML tree based on the entire plastome maximally supported all of the deep and most recent divergences of papilionoids (223 out of 236 nodes). The Swartzieae, ADA (Angylocalyceae, Dipterygeae, and Amburaneae), Cladrastis, Andira, and Exostyleae clades formed a grade to the remainder of the Papilionoideae, concordant with nine ML and two MSC trees. Phylogenetic relationships among the remaining five papilionoid lineages (Vataireoid, Dermatophyllum , Genistoid s.l., Dalbergioid s.l., and Baphieae + Non-Protein Amino Acid Accumulating or NPAAA clade) remained uncertain, because of insufficient support and/or conflicting relationships among trees. Our study fully resolved most of the deep nodes of Papilionoideae, however, some relationships require further exploration. More genome-scale data and rigorous analyses are needed to disentangle phylogenetic relationships among the five remaining lineages., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Choi, Cardoso, de Queiroz, de Lima, Lee, Ruhlman, Jansen and Wojciechowski.)

Plant specific mitoviruses in the 'genus' Mitovirus (Narnaviridae) and their integrated sequences (non-retroviral endogenous RNA viral elements or NERVEs) have been recently identified in various plant lineages. However, the sparse phylogenetic coverage of complete plant mitochondrial genome (mitogenome) sequences and the non-conserved nature of mitochondrial intergenic regions have hindered comparative studies on mitovirus NERVEs in plants. In this study, 10 new mitogenomes were sequenced from legumes (Fabaceae). Based on comparative genomic analysis of 27 total mitogenomes, we identified mitovirus NERVEs and transposable elements across the family. All legume mitogenomes included NERVEs and total NERVE length varied from ca. 2 kb in the papilionoid Trifolium to 35 kb in the mimosoid Acacia. Most of the NERVE integration sites were in highly variable intergenic regions, however, some were positioned in six cis-spliced mitochondrial introns. In the Acacia mitogenome, there were L1-like transposon sequences including an almost full-length copy with target site duplications (TSDs). The integration sites of NERVEs in four introns showed evidence of L1-like retrotransposition events. Phylogenetic analysis revealed that there were multiple instances of precise deletion of NERVEs between TSDs. This study provides clear evidence that a L1-like retrotransposition mechanism has a long history of contributing to the integration of viral RNA into plant mitogenomes while microhomology-mediated deletion can restore the integration site., (Copyright © 2021 Elsevier Inc. All rights reserved.)

The plastid genome (plastome), while surprisingly constant in gene order and content across most photosynthetic angiosperms, exhibits variability in several unrelated lineages. During the diversification history of the legume family Fabaceae, plastomes have undergone many rearrangements, including inversions, expansion, contraction and loss of the typical inverted repeat (IR), gene loss and repeat accumulation in both shared and independent events. While legume plastomes have been the subject of study for some time, most work has focused on agricultural species in the IR-lacking clade (IRLC) and the plant model Medicago truncatula. The subfamily Papilionoideae, which contains virtually all of the agricultural legume species, also comprises most of the plastome variation detected thus far in the family. In this study three non-papilioniods were included among 34 newly sequenced legume plastomes, along with 33 publicly available sequences, to assess plastome structural evolution in the subfamily. In an effort to examine plastome variation across the subfamily, approximately 20% of the sampling represents the IRLC with the remainder selected to represent the early-branching papilionoid clades. A number of IR-related and repeat-mediated changes were identified and examined in a phylogenetic context. Recombination between direct repeats associated with ycf2 resulted in intraindividual plastome heteroplasmy. Although loss of the IR has not been reported in legumes outside of the IRLC, one genistoid taxon was found to completely lack the typical plastome IR. The role of the IR and non-IR repeats in the progression of plastome change is discussed., (© 2021 Society for Experimental Biology and John Wiley & Sons Ltd.)

Barley is characterized by a rich genetic diversity, making it an important model for studies of salinity response with great potential for crop improvement. Moreover, salt stress severely affects barley growth and development, leading to substantial yield loss. Leaf and root transcriptomes of a salt-tolerant Tunisian landrace (Boulifa) exposed to 2, 8, and 24 h salt stress were compared with pre-exposure plants to identify candidate genes and pathways underlying barley's response. Expression of 3585 genes was upregulated and 5586 downregulated in leaves, while expression of 13,200 genes was upregulated and 10,575 downregulated in roots. Regulation of gene expression was severely impacted in roots, highlighting the complexity of salt stress response mechanisms in this tissue. Functional analyses in both tissues indicated that response to salt stress is mainly achieved through sensing and signaling pathways, strong transcriptional reprograming, hormone osmolyte and ion homeostasis stabilization, increased reactive oxygen scavenging, and activation of transport and photosynthesis systems. A number of candidate genes involved in hormone and kinase signaling pathways, as well as several transcription factor families and transporters, were identified. This study provides valuable information on early salt-stress-responsive genes in roots and leaves of barley and identifies several important players in salt tolerance.

Salt stress is considered one of the most devastating environmental stresses, affecting barley growth and leading to significant yield loss. Hence, there is considerable interest in investigating the most effective traits that determine barley growth under salt stress. The objective of this study was to elucidate the contribution of osmotic and oxidative stress components in leaves and roots growth under salt stress. Two distinct barley ( Hordeum vulgare ) salt-stress tolerant genotypes, Barrage Malleg (BM, tolerant) and Saouef (Sf, sensitive), were subjected to 200 mM NaCl at early vegetative stages. Stressed and control leaves and roots tissue were assessed for several growth traits, including fresh and dry weight and plant length, as well as the content of osmoprotectants proline and soluble sugars. In addition, malondialdehyde content and activities of superoxide dismutase (SOD), catalase (CAT) and ascorbate peroxidase (APX), as well as their corresponding gene expression patterns, were investigated. The results showed better performance of BM over Sf for leaf dry weight (LDW), root dry weight (RDW) and root length (RL). The salt-tolerant genotype (BM) had better osmoprotection against salt stress compared with the salt-sensitive genotype (Sf), with a higher accumulation of proline and soluble sugars in leaves and roots and a stronger antioxidant system as evidenced by higher activities of SOD, CAT and APX and more abundant Cu/Zn-SOD transcripts, especially in roots. Stepwise regression analysis indicated that under salt stress the most predominant trait of barley growth was Cu/Zn-SOD gene expression level, suggesting that alleviating oxidative stress and providing cell homeostasis is the first priority., (© The Author(s) 2021. Published by Oxford University Press on behalf of the Annals of Botany Company.)

Gene transfers from mitochondria and plastids to the nucleus are an important process in the evolution of the eukaryotic cell. Plastid (pt) gene losses have been documented in multiple angiosperm lineages and are often associated with functional transfers to the nucleus or substitutions by duplicated nuclear genes targeted to both the plastid and mitochondrion. The plastid genome sequence of Euphorbia schimperi was assembled and three major genomic changes were detected, the complete loss of rpl32 and pseudogenization of rps16 and infA. The nuclear transcriptome of E. schimperi was sequenced to investigate the transfer/substitution of the rpl32 and rps16 genes to the nucleus. Transfer of plastid-encoded rpl32 to the nucleus was identified previously in three families of Malpighiales, Rhizophoraceae, Salicaceae and Passifloraceae. An E. schimperi transcript of pt SOD-1-RPL32 confirmed that the transfer in Euphorbiaceae is similar to other Malpighiales indicating that it occurred early in the divergence of the order. Ribosomal protein S16 (rps16) is encoded in the plastome in most angiosperms but not in Salicaceae and Passifloraceae. Substitution of the E. schimperi pt rps16 was likely due to a duplication of nuclear-encoded mitochondrial-targeted rps16 resulting in copies dually targeted to the mitochondrion and plastid. Sequences of RPS16-1 and RPS16-2 in the three families of Malpighiales (Salicaceae, Passifloraceae and Euphorbiaceae) have high sequence identity suggesting that the substitution event dates to the early divergence within Malpighiales.

Plastid inheritance in angiosperms is presumed to be largely maternal, with the potential to inherit plastids biparentally estimated for about 20% of species. In Passiflora , maternal, paternal and biparental inheritance has been reported; however, these studies were limited in the number of crosses and progeny examined. To improve the understanding of plastid transmission in Passiflora, the progeny of 45 interspecific crosses were analyzed in the three subgenera: Passiflora, Decaloba and Astrophea. Plastid types were assessed following restriction digestion of PCR amplified plastid DNA in hybrid embryos, cotyledons and leaves at different developmental stages. Clade-specific patterns of inheritance were detected such that hybrid progeny from subgenera Passiflora and Astrophea predominantly inherited paternal plastids with occasional incidences of maternal inheritance, whereas subgenus Decaloba showed predominantly maternal and biparental inheritance. Biparental plastid inheritance was also detected in some hybrids from subgenus Passiflora. Heteroplasmy due to biparental inheritance was restricted to hybrid cotyledons and first leaves with a single parental plastid type detectable in mature plants. This indicates that in Passiflora , plastid retention at later stages of plant development may not reflect the plastid inheritance patterns in embryos. Passiflora exhibits diverse patterns of plastid inheritance, providing an excellent system to investigate underlying mechanisms in angiosperms.

Geraniaceae organelle genomes have been shown to exhibit several highly unusual features compared to most other photosynthetic angiosperms. This includes massively rearranged plastomes with considerable size variation, extensive gene and intron loss, accelerated rates of nucleotide substitutions in both mitogenomes and plastomes, and biparental inheritance and cytonuclear incompatibility of the plastome. Most previous studies have focused on plastome evolution with mitogenome comparisons limited to only a few taxa or genes. In this study, mitogenomes and transcriptomes were examined for 27 species of Geraniales, including 13 species of Pelargonium. Extensive gene and intron losses were detected across the Geraniales with Pelargonium representing the most gene depauperate lineage in the family. Plotting these events on the Geraniaceae phylogenetic tree showed that gene losses occurred multiple times, whereas intron losses more closely reflected the relationships among taxa. In addition, P. australe acquired an intron by horizontal transfer. Comparisons of nucleotide substitution rates in Pelargonium showed that synonymous changes in nuclear genes were much lower than in mitochondrial genes. This is in contrast to the previously published studies that indicated that nuclear genes have 16 fold higher rates than mitochondrial genes across angiosperms. Elevated synonymous substitutions occurred for each mitochondrial gene in Pelargonium with the highest values 783 and 324 times higher than outgroups and other Geraniaceae, respectively. Pelargonium is one of four unrelated genera of angiosperms (Ajuga, Plantago and Silene) that have experienced highly accelerated nucleotide substitutions in mitogenomes. It is distinct from most angiosperms in also having elevated substitution rates in plastid genes but the cause of rate accelerations in Pelargonium plastomes and mitogenomes may be different., (Copyright © 2020 Elsevier Inc. All rights reserved.)

The plastid genome (plastome ) has proved a valuable source of data for evaluating evolutionary relationships among angiosperms. Through basic and applied approaches, plastid transformation technology offers the potential to understand and improve plant productivity, providing food, fiber, energy, and medicines to meet the needs of a burgeoning global population. The growing genomic resources available to both phylogenetic and biotechnological investigations is allowing novel insights and expanding the scope of plastome research to encompass new species. In this chapter, we present an overview of some of the seminal and contemporary research that has contributed to our current understanding of plastome evolution and attempt to highlight the relationship between evolutionary mechanisms and the tools of plastid genetic engineering.

Diatoms are the largest group of heterokont algae with more than 100,000 species. As one of the single-celled photosynthetic organisms that inhabit marine, aquatic and terrestrial ecosystems, diatoms contribute ~ 45% of global primary production. Despite their ubiquity and environmental significance, very few diatom plastid genomes (plastomes) have been sequenced and studied. This study explored patterns of nucleotide substitution rates of diatom plastids across the entire suite of plastome protein-coding genes for 40 taxa representing the major clades. The highest substitution rate was lineage-specific within the araphid 2 taxon Astrosyne radiata and radial 2 taxon Proboscia sp. Rate heterogeneity was also evident in different functional classes and individual genes. Similar to land plants, proteins genes involved in photosynthetic metabolism have lower synonymous and nonsynonymous substitutions rates than those involved in transcription and translation. Significant positive correlations were identified between substitution rates and measures of genomic rearrangements, including indels and inversions, which is a similar result to what was found in legume plants. This work advances the understanding of the molecular evolution of diatom plastomes and provides a foundation for future studies.

Gene losses in plastid genomes (plastomes) are often accompanied by functional transfer to the nucleus or substitution of an alternative nuclear-encoded gene. Despite the highly conserved gene content in plastomes of photosynthetic land plants, recent gene loss events have been documented in several disparate angiosperm clades. Among these lineages, Passiflora lacks several essential ribosomal genes, rps7, rps16, rpl20, rpl22, and rpl32, the two largest plastid genes, ycf1 and ycf2, and has a highly divergent rpoA. Comparative transcriptome analyses were performed to determine the fate of the missing genes in Passiflora. Putative functional transfers of rps7, rpl22, and rpl32 to nucleus were detected, with the nuclear transfer of rps7, representing a novel event in angiosperms. Plastid-encoded rps7 was transferred into the intron of a nuclear-encoded plastid-targeted thioredoxin m-type gene, acquiring its plastid transit peptide (TP). Plastid rpl20 likely experienced a novel substitution by a duplicated, nuclear-encoded mitochondrial-targeted rpl20 that has a similar gene structure. Additionally, among rosids, evidence for a third independent transfer of rpl22 in Passiflora was detected that gained a TP from a nuclear gene containing an organelle RNA recognition motif. Nuclear transcripts representing rpoA, ycf1, and ycf2 were not detected. Further analyses suggest that the divergent rpoA remains functional and that the gene is under positive or purifying selection in different clades. Comparative analyses indicate that alternative translocon and motor protein complexes may have substituted for the loss of ycf1 and ycf2 in Passiflora., (© The Author(s) 2020. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Premise: Despite the fast pace of exploration of the patterns and processes influencing Neotropical plant hyperdiversity, the taxa explored are mostly from large groups that are widely distributed, morphologically diverse, or economically important. Vochysiaceae is an example of an undersampled taxon, providing an excellent system for investigating Neotropical biogeography. We present a phylogenomics-based hypothesis of species relationships in Vochysiaceae to investigate its evolutionary history through space and time., Methods: We inferred a phylogeny for 122 species from Vochysiaceae and seven other families of Myrtales. Fossils from four myrtalean families were used to estimate the divergence times within Vochysiaceae. Historical biogeography was estimated using ancestral range probabilities and stochastic mapping., Results: Monophyly of all genera was supported except for Qualea, which was split by Ruizterania into two clades. Vochysiaceae originated ~100 mya, splitting into an Afrotropical and a Neotropical lineage ~50 mya, and its ancestral range is in the area currently occupied by the Cerrado., Conclusions: The most recent common ancestor of Vochysiaceae + Myrtaceae had a West Gondwanan distribution, supporting a South American + African ancestral range of Vochysiaceae. On a global scale, geographic range reduction was the principal biogeographic event. At a finer scale, initial range reduction was also important and the Cerrado region was the most ancestral area with multiple colonization events to the Amazon, Central America, and the Atlantic Forest. Colonization events occurred from open areas to forest vegetation, an unusual finding regarding the evolution of plants in the Neotropics., (© 2020 Botanical Society of America.)

Plastid genomes (plastomes) of land plants have a conserved quadripartite structure in a gene-dense unit genome consisting of a large inverted repeat that separates two single copy regions. Recently, alternative plastome structures were suggested in Geraniaceae and in some conifers and Medicago the coexistence of inversion isomers has been noted. In this study, plastome sequences of two Cyperaceae, Eleocharis dulcis (water chestnut) and Eleocharis cellulosa (gulf coast spikerush), were completed. Unlike the conserved plastomes in basal groups of Poales, these Eleocharis plastomes have remarkably divergent features, including large plastome sizes, high rates of sequence rearrangements, low GC content and gene density, gene duplications and losses, and increased repetitive DNA sequences. A novel finding among these features was the unprecedented level of heteroplasmy with the presence of multiple plastome structural types within a single individual. Illumina paired-end assemblies combined with PacBio single-molecule real-time sequencing, long-range polymerase chain reaction, and Sanger sequencing data identified at least four different plastome structural types in both Eleocharis species. PacBio long read data suggested that one of the four E. dulcis plastome types predominates., (© The Author(s) 2020. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

The genus Trifolium is the largest of the tribe Trifolieae in the subfamily Papilionoideae (Fabaceae). The paucity of mitochondrial genome (mitogenome) sequences has hindered comparative analyses among the three genomic compartments of the plant cell (nucleus, mitochondrion and plastid). We assembled four mitogenomes from the two subgenera ( Chronosemium and Trifolium ) of the genus. The four Trifolium mitogenomes were compact (294,911-348,724 bp in length) and contained limited repetitive (6.6-8.6%) DNA. Comparison of organelle repeat content highlighted the distinct evolutionary trajectory of plastid genomes in a subset of Trifolium species. Intracellular gene transfer (IGT) was analyzed among the three genomic compartments revealing functional transfer of mitochondrial rps1 to nuclear genome along with other IGT events. Phylogenetic analysis based on mitochondrial and nuclear rps1 sequences revealed that the functional transfer in Trifolieae was independent from the event that occurred in robinioid clade that includes genus Lotus . A novel, independent fission event of ccmFn in Trifolium was identified, caused by a 59 bp deletion. Fissions of this gene reported previously in land plants were reassessed and compared with Trifolium .

Obesity is a multifactorial disease caused by complex interactions between genes and dietary factors. Salt-rich diet is related to the development and progression of several chronic diseases including obesity. However, the molecular basis of how salt sensitivity genes (SSG) contribute to adiposity in obesity patients remains unexplored. In this study, we used the microarray expression data of visceral adipose tissue samples and constructed a complex protein-interaction network of salt sensitivity genes and their co-expressed genes to trace the molecular pathways connected to obesity. The Salt Sensitivity Protein Interaction Network (SSPIN) of 2691 differentially expressed genes and their 15474 interactions has shown that adipose tissues are enriched with the expression of 23 SSGs, 16 hubs and 84 bottlenecks (p = 2.52 x 10-16) involved in diverse molecular pathways connected to adiposity. Fifteen of these 23 SSGs along with 8 other SSGs showed a co-expression with enriched obesity-related genes (r ≥ 0.8). These SSGs and their co-expression partners are involved in diverse metabolic pathways including adipogenesis, adipocytokine signaling pathway, renin-angiotensin system, etc. This study concludes that SSGs could act as molecular signatures for tracing the basis of adipogenesis among obese patients. Integrated network centered methods may accelerate the identification of new molecular targets from the complex obesity genomics data., Competing Interests: The authors have declared that no competing interests exist.

Obesity is a multifactorial metabolic disorder characterized by low grade chronic inflammation. Rare and novel mutations in genes which are vital in several key pathways have been reported to alter the energy expenditure which regulates body weight. The TP53 or p53 gene plays a prominent role in regulating various metabolic activities such as glycolysis, lipolysis, and glycogen synthesis. Recent genome-wide association studies reported that tumor suppressor gene p53 variants play a critical role in the predisposition of type 2 diabetes and obesity. Till date, no reports are available from the Arabian population; hence the present study was intended to assess the association between p53 variants with risk of obesity development in the Saudi population. We have selected three p53 polymorphisms, rs1642785 (C > G), and rs9894946 (A > G), and rs1042522 (Pro72Arg; C > G) and assessed their association with obesity risk in the Saudi population. Phenotypic and biochemical parameters were also evaluated to check their association with p53 genotypes and obesity. Genotyping was carried out on 136 obese and 122 normal samples. We observed that there is significantly increased prevalence p52 Pro72Arg (rs1042522) polymorphism in obese persons when compared to controls at GG genotype in overall comparison (OR: 2.169, 95% CI: 1.086-4.334, p = 0.02716). Male obese subjects showed three-fold higher risk at GG genotype (OR: 3.275, 95% CI: 1.230-8.716, p = 0.01560) and two-fold risk at G allele (OR: 1.827, 95% CI: 1.128-2.958, p = 0.01388) of p53 variant Pro72Arg respectively. This variant has also shown significant influence on cholesterol, LDL level, and random insulin levels in obese subjects (p ≤ 0.05). In conclusion, p53 Pro72Arg variant is highly prevalent among obese individuals and may act as a genetic modifier for obesity development among Saudis., (© 2019 Production and hosting by Elsevier B.V. on behalf of King Saud University.)

Background: Organelle genome studies of Fabaceae, an economically and ecologically important plant family, have been biased towards the plastid genome (plastome). Thus far, less than 15 mitochondrial genome (mitogenome) sequences of Fabaceae have been published, all but four of which belong to the subfamily Papilionoideae, limiting the understanding of size variation and content across the family. To address this, four mitogenomes were sequenced and assembled from three different subfamilies (Cercidoideae, Detarioideae and Caesalpinioideae)., Results: Phylogenetic analysis based on shared mitochondrial protein coding regions produced a fully resolved and well-supported phylogeny that was completely congruent with the plastome tree. Comparative analyses suggest that two kinds of mitogenome expansions have occurred in Fabaceae. Size expansion of four genera (Tamarindus, Libidibia, Haematoxylum, and Leucaena) in two subfamilies (Detarioideae and Caesalpinioideae) occurred in relatively deep nodes, and was mainly caused by intercellular gene transfer and/or interspecific horizontal gene transfer (HGT). The second, more recent expansion occurred in the Papilionoideae as a result of duplication of native mitochondrial sequences. Family-wide gene content analysis revealed 11 gene losses, four (rps2, 7, 11 and 13) of which occurred in the ancestor of Fabaceae. Losses of the remaining seven genes (cox2, rpl2, rpl10, rps1, rps19, sdh3, sdh4) were restricted to specific lineages or occurred independently in different clades. Introns of three genes (cox2, ccmFc and rps10) showed extensive lineage-specific length variation due to large sequence insertions and deletions. Shared DNA analysis among Fabaceae mitogenomes demonstrated a substantial decay of intergenic spacers and provided further insight into HGT between the mimosoid clade of Caesalpinioideae and the holoparasitic Lophophytum (Balanophoraceae)., Conclusion: This study represents the most exhaustive analysis of Fabaceae mitogenomes so far, and extends the understanding the dynamic variation in size and gene/intron content. The four newly sequenced mitogenomes reported here expands the phylogenetic coverage to four subfamilies. The family has experienced multiple mitogenome size fluctuations in both ancient and recent times. The causes of these size variations are distinct in different lineages. Fabaceae mitogenomes experienced extensive size fluctuation by recruitment of exogenous DNA and duplication of native mitochondrial DNA.

Plastid genomes (plastomes) of photosynthetic angiosperms are for the most part highly conserved in their organization, mode of inheritance and rates of nucleotide substitution. A small number of distantly related lineages share a syndrome of features that deviate from this general pattern, including extensive genomic rearrangements, accelerated rates of nucleotide substitution, biparental inheritance and plastome-genome incompatibility. Previous studies of plastomes in Passiflora with limited taxon sampling suggested that the genus exhibits this syndrome. To examine this phenomenon further, 15 new plastomes from Passiflora were sequenced and combined with previously published data to examine the phylogenetic relationships, genome organization and evolutionary rates across all five subgenera and the sister genus Adenia. Phylogenomic analyses using 68 protein-coding genes shared by Passiflora generated a fully resolved and strongly supported tree that is congruent with previous phylogenies based on a few plastid and nuclear loci. This phylogeny was used to examine the distribution of plastome rearrangements across Passiflora. Multiple gene and intron losses and inversions were identified in Passiflora with some occurring in parallel and others that extended across the Passifloraceae. Furthermore, extensive expansions and contractions of the inverted repeat (IR) were uncovered and in some cases this resulted in exclusion of all ribosomal RNA genes from the IR. The most highly rearranged lineage was subgenus Decaloba, which experienced extensive IR expansion that incorporated up to 25 protein-coding genes usually located in large single copy region. Nucleotide substitution rate analyses of 68 protein-coding genes across the genus showed lineage- and locus-specific acceleration. Significant increase in dS, dN and dN/dS was detected for clpP across the genus and for ycf4 in certain lineages. Significant increases in dN and dN/dS for ribosomal subunits and plastid-encoded RNA polymerase genes were detected in the branch leading to the expanded IR-clade in subgenus Decaloba. This subgenus displays the syndrome of unusual features, making it an ideal system to investigate the dynamic evolution of angiosperm plastomes., (Copyright © 2019 Elsevier Inc. All rights reserved.)

The current classification of angiosperms is based primarily on concatenated plastid markers and maximum likelihood (ML) inference. This approach has been justified by the assumption that plastid DNA (ptDNA) is inherited as a single locus and that its individual genes produce congruent trees. However, structural and functional characteristics of ptDNA suggest that plastid genes may not evolve as a single locus and are experiencing different evolutionary forces. To examine this idea, we produced new complete plastid genome (plastome) sequences of 27 species and combined these data with publicly available sequences to produce a final dataset that includes 78 plastid genes for 89 species of rosids and five outgroups. We used four data matrices (i.e., gene, exon, codon-aligned, and amino acid) to infer species and gene trees using ML and multispecies coalescent (MSC) methods. Rosids include about one third of all angiosperms and their two major clades, fabids and malvids, were recovered in almost all analyses. However, we detected incongruence between species trees inferred with different matrices and methods and previously published plastid and nuclear phylogenies. We visualized and tested the significance of incongruence between gene trees and species trees. We then measured the distribution of phylogenetic signal across sites and genes supporting alternative placements of five controversial nodes at different taxonomic levels. Gene trees inferred with plastid data often disagree with species trees inferred using both ML (with unpartitioned or partitioned data) and MSC. Species trees inferred with both methods produced alternative topologies for a few taxa. Our results show that, in a phylogenetic context, plastid protein-coding genes may not be fully linked and behaving as a single locus. Furthermore, concatenated matrices may produce highly supported phylogenies that are discordant with individual gene trees. We also show that phylogenies inferred with MSC are accurate. We therefore emphasize the importance of considering variation in phylogenetic signal across plastid genes and the exploration of plastome data to increase accuracy of estimating relationships. We also support the use of MSC with plastome matrices in future phylogenomic investigations., (Copyright © 2019 Elsevier Inc. All rights reserved.)

The predominantly South-African plant genus Pelargonium L'Hér. (Geraniaceae) displays remarkable morphological diversity, several basic chromosome numbers as well as high levels of organelle genomic rearrangements, and represents the 7th largest Cape Floristic Region clade. In this study, we reconstructed a phylogenetic tree based on 74 plastome exons and nuclear rDNA ITS regions for 120 species, which represents 43% taxon coverage for Pelargonium. We also performed a dating analysis to examine the timing of the major radiations in the genus. Phylogenetic analyses of nucleotide, amino acid, and ITS alignments confirmed the previously-documented subgeneric split into five main clades ((C1,C2),(B(A1,A2))) although clade only A1 received low bootstrap support. Using calibration evidence from a range of sources the Pelargonium crown age was estimated to be 9.7 My old, much younger than previous estimates for the genus but similar to recent studies of other Cape Floristic lineages that are part of both Fynbos and Succulent Karoo biomes., (Copyright © 2019 Elsevier Inc. All rights reserved.)

Phylogenomics has become increasingly popular in recent years mostly due to the increased affordability of next generation sequencing techniques. Phylogenomics has sparked interest in multiple fields of research, including systematics, ecology, epidemiology, and even personalized medicine, agriculture and pharmacy. Despite this trend, it is usually difficult to learn and understand how the analyses were done, how the results were obtained, and most importantly, how to replicate the study. Here we present the data and all of the code utilized to perform phylogenomic inferences using plastome data: from raw data to extensive phylogenetic inference and accuracy assessment. The data presented here utilizes plastome sequences available on GenBank (accession numbers of 94 species are available below) and the code is also available at https://github.com/deisejpg/rosids. Gonçalves et al. is the research article associated with the data analyses presented here.

The diatom Phaeodactylum tricornutum has been used as a model for cell biologists and ecologists for over a century. We have incorporated several new raphid pennates into a three gene phylogenetic dataset (SSU, rbcL, psbC), and recover Gomphonemopsis sp. as sister to P. tricornutum with 100% BS support. This is the first time a close relative has been identified for P. tricornutum with robust statistical support. We test and reject a succession of hypotheses for other relatives. Our molecular data are statistically significantly incongruent with placement of either or both species among the Cymbellales, an order of diatoms with which both have been associated. We believe that further resolution of the phylogenetic position of P. tricornutum will rely more on increased taxon sampling than increased genetic sampling. Gomphonemopsis is a benthic diatom, and its phylogenetic relationship with P. tricornutum is congruent with the hypothesis that P. tricornutum is a benthic diatom with specific adaptations that lead to active recruitment into the plankton. We hypothesize that other benthic diatoms are likely to have similar adaptations and are not merely passively recruited into the plankton., Competing Interests: The authors have declared that no competing interests exist.

Background: Calotropis procera is a wild plant species in the family Apocynaceae that is able to grow in harsh, arid and heat stressed conditions. Understanding how this highly adapted plant persists in harsh environments should inform future efforts to improve the hardiness of crop and forage plant species. To study the plant response to droμght and osmotic stress, we treated plants with polyethylene glycol and NaCl and carried out transcriptomic and metabolomics measurements across a time-course of five days., Results: We identified a highly dynamic transcriptional response across the time-course including dramatic changes in inositol signaling, stress response genes and cytokinins. The resulting metabolome changes also involved sharp increases of myo-inositol, a key signaling molecule and elevated amino acid metabolites at later times., Conclusions: The data generated here provide a first glimpse at the expressed genome of C. procera, a plant that is exceptionally well adapted to arid environments. We demonstrate, through transcriptome and metabolome analysis that myo-inositol signaling is strongly induced in response to drought and salt stress and that there is elevation of amino acid concentrations after prolonged osmotic stress. This work should lay the foundations of future studies in adaptation to arid environments.

In plant evolution, intracellular gene transfer (IGT) is a prevalent, ongoing process. While nuclear and mitochondrial genomes are known to integrate foreign DNA via IGT and horizontal gene transfer (HGT), plastid genomes (plastomes) have resisted foreign DNA incorporation and only recently has IGT been uncovered in the plastomes of a few land plants. In this study, we completed plastome sequences for l0 crop species and describe a number of structural features including variation in gene and intron content, inversions, and expansion and contraction of the inverted repeat (IR). We identified a putative in cinnamon ( J. Presl) and other sequenced Lauraceae and an apparent functional transfer of to the nucleus of quinoa ( Willd.). In the orchard tree cashew ( L.), we report the insertion of an ∼6.7-kb fragment of mitochondrial DNA into the plastome IR. BLASTn analyses returned high identity hits to mitogenome sequences including an intact open reading frame. Using three plastome markers for five species of , we generated a phylogeny to investigate the distribution and timing of the insertion. Four species share the insertion, suggesting that this event occurred <20 million yr ago in a single clade in the genus. Our study extends the observation of mitochondrial to plastome IGT to include long-lived tree species. While previous studies have suggested possible mechanisms facilitating IGT to the plastome, more examples of this phenomenon, along with more complete mitogenome sequences, will be required before a common, or variable, mechanism can be elucidated., (Copyright © 2017 Crop Science Society of America.)

Geraniaceae have emerged as a model system for investigating the causes and consequences of variation in plastid and mitochondrial genomes. Incredible structural variation in plastid genomes (plastomes) and highly accelerated evolutionary rates have been reported in selected lineages and functional groups of genes in both plastomes and mitochondrial genomes (mitogenomes), and these phenomena have been implicated in cytonuclear incompatibility. Previous organelle genome studies have included limited sampling of Geranium, the largest genus in the family with over 400 species. This study reports on rates and patterns of nucleotide substitutions in plastomes and mitogenomes of 17 species of Geranium and representatives of other Geraniaceae. As detected across other angiosperms, substitution rates in the plastome are 3.5 times higher than the mitogenome in most Geranium. However, in the branch leading to Geranium brycei/Geranium incanum mitochondrial genes experienced significantly higher dN and dS than plastid genes, a pattern that has only been detected in one other angiosperm. Furthermore, rate accelerations differ in the two organelle genomes with plastomes having increased dN and mitogenomes with increased dS. In the Geranium phaeum/Geranium reflexum clade, duplicate copies of clpP and rpoA genes that experienced asymmetric rate divergence were detected in the single copy region of the plastome. In the case of rpoA, the branch leading to G. phaeum/G. reflexum experienced positive selection or relaxation of purifying selection. Finally, the evolution of acetyl-CoA carboxylase is unusual in Geraniaceae because it is only the second angiosperm family where both prokaryotic and eukaryotic ACCases functionally coexist in the plastid., (© The Author 2017. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Rhazya stricta is an evergreen shrub that is widely distributed across Western and South Asia, and like many other members of the Apocynaceae produces monoterpene indole alkaloids that have anti-cancer properties. This species is adapted to very harsh desert conditions making it an excellent system for studying tolerance to high temperatures and salinity. RNA-Seq analysis was performed on R. stricta exposed to severe salt stress (500 mM NaCl) across four time intervals (0, 2, 12 and 24 h) to examine mechanisms of salt tolerance. A large number of transcripts including genes encoding tetrapyrroles and pentatricopeptide repeat (PPR) proteins were regulated only after 12 h of stress of seedlings grown in controlled greenhouse conditions. Mechanisms of salt tolerance in R. stricta may involve the upregulation of genes encoding chaperone protein Dnaj6, UDP-glucosyl transferase 85a2, protein transparent testa 12 and respiratory burst oxidase homolog protein b. Many of the highly-expressed genes act on protecting protein folding during salt stress and the production of flavonoids, key secondary metabolites in stress tolerance. Other regulated genes encode enzymes in the porphyrin and chlorophyll metabolic pathway with important roles during plant growth, photosynthesis, hormone signaling and abiotic responses. Heme biosynthesis in R. stricta leaves might add to the level of salt stress tolerance by maintaining appropriate levels of photosynthesis and normal plant growth as well as by the participation in reactive oxygen species (ROS) production under stress. We speculate that the high expression levels of PPR genes may be dependent on expression levels of their targeted editing genes. Although the results of PPR gene family indicated regulation of a large number of transcripts under salt stress, PPR actions were independent of the salt stress because their RNA editing patterns were unchanged.

This study represents the most comprehensive plastome-wide comparison of nucleotide substitution rates across the three subfamilies of Fabaceae: Caesalpinioideae, Mimosoideae, and Papilionoideae. Caesalpinioid and mimosoid legumes have large, unrearranged plastomes compared with papilionoids, which exhibit varying levels of rearrangement including the loss of the inverted repeat (IR) in the IR-lacking clade (IRLC). Using 71 genes common to 39 legume taxa representing all the three subfamilies, we show that papilionoids consistently have higher nucleotide substitution rates than caesalpinioids and mimosoids, and rates in the IRLC papilionoids are generally higher than those in the IR-containing papilionoids. Unsurprisingly, this pattern was significantly correlated with growth habit as most papilionoids are herbaceous, whereas caesalpinioids and mimosoids are largely woody. Both nonsynonymous (dN) and synonymous (dS) substitution rates were also correlated with several biological features including plastome size and plastomic rearrangements such as the number of inversions and indels. In agreement with previous reports, we found that genes in the IR exhibit between three and fourfold reductions in the substitution rates relative to genes within the large single-copy or small single-copy regions. Furthermore, former IR genes in IR-lacking taxa exhibit accelerated rates compared with genes contained in the IR.

For species with minor inverted repeat (IR) boundary changes in the plastid genome (plastome), nucleotide substitution rates were previously shown to be lower in the IR than the single copy regions (SC). However, the impact of large-scale IR expansion/contraction on plastid nucleotide substitution rates among closely related species remains unclear. We included plastomes from 22 Pelargonium species, including eight newly sequenced genomes, and used both pairwise and model-based comparisons to investigate the impact of the IR on sequence evolution in plastids. Ten types of plastome organization with different inversions or IR boundary changes were identified in Pelargonium. Inclusion in the IR was not sufficient to explain the variation of nucleotide substitution rates. Instead, the rate heterogeneity in Pelargonium plastomes was a mixture of locus-specific, lineage-specific and IR-dependent effects. Our study of Pelargonium plastomes that vary in IR length and gene content demonstrates that the evolutionary consequences of retaining these repeats are more complicated than previously suggested., (© 2016 The Authors. New Phytologist © 2016 New Phytologist Trust.)

Premise of the Study: There is a misinterpretation in the literature regarding the variable orientation of the small single copy region of plastid genomes (plastomes). The common phenomenon of small and large single copy inversion, hypothesized to occur through intramolecular recombination between inverted repeats (IR) in a circular, single unit-genome, in fact, more likely occurs through recombination-dependent replication (RDR) of linear plastome templates. If RDR can be primed through both intra- and intermolecular recombination, then this mechanism could not only create inversion isomers of so-called single copy regions, but also an array of alternative sequence arrangements., Methods: We used Illumina paired-end and PacBio single-molecule real-time (SMRT) sequences to characterize repeat structure in the plastome of Monsonia emarginata (Geraniaceae). We used OrgConv and inspected nucleotide alignments to infer ancestral nucleotides and identify gene conversion among repeats and mapped long (>1 kb) SMRT reads against the unit-genome assembly to identify alternative sequence arrangements., Results: Although M. emarginata lacks the canonical IR, we found that large repeats (>1 kilobase; kb) represent ∼22% of the plastome nucleotide content. Among the largest repeats (>2 kb), we identified GC-biased gene conversion and mapping filtered, long SMRT reads to the M. emarginata unit-genome assembly revealed alternative, substoichiometric sequence arrangements., Conclusion: We offer a model based on RDR and gene conversion between long repeated sequences in the M. emarginata plastome and provide support that both intra-and intermolecular recombination between large repeats, particularly in repeat-rich plastomes, varies unit-genome structure while homogenizing the nucleotide sequence of repeats., (© 2017 Botanical Society of America.)

Several lineages of raphe-bearing diatoms possess a "stauros," which is a transverse, usually thickened area free of pores across the center of the valve. It has been suggested that this structure has evolved several times across the raphid diatoms, but we have noticed similarities beyond the stauros between two marine genera-Craspedostauros and Staurotropis-in the structure of their pore occlusions. We have isolated, cultured and extracted DNA from several strains of both genera to infer the phylogenetic relationship between these taxa, as well as test the suggested relationship of Craspedostauros to Achnanthes and Mastogloia based on plastid morphology. DNA sequence data (nuclear-encoded rRNA SSU, plastid-encoded rbcL and psbC) suggest that, except for Mastogloia, these genera are closely-related, though not sister taxa. The DNA phylogeny also suggests that the Mastogloiales are not monophyletic, with clades containing Achnanthes and Craspedostauros sister to clades containing taxa in the Bacillariales. Using evidence from molecular and morphological data, we describe the following new taxa: Craspedostauros alyoubii and C. paradoxa from the Red Sea and Guam, respectively; Staurotropis khiyamii and S. americana from the Red Sea and the Gulf of Mexico, respectively; and Dreuhlago cuneata n. gen., n. sp. from Guam., (Copyright © 2016 Elsevier GmbH. All rights reserved.)

Background: Climate change is predicted to be a serious threat to agriculture due to the need for crops to be able to tolerate increased heat stress. Desert plants have already adapted to high levels of heat stress so they make excellent systems for identifying genes involved in thermotolerance. Rhazya stricta is an evergreen shrub that is native to extremely hot regions across Western and South Asia, making it an excellent system for examining plant responses to heat stress. Transcriptomes of apical and mature leaves of R. stricta were analyzed at different temperatures during several time points of the day to detect heat response mechanisms that might confer thermotolerance and protection of the plant photosynthetic apparatus., Results: Biological pathways that were crosstalking during the day involved the biosynthesis of several heat stress-related compounds, including soluble sugars, polyols, secondary metabolites, phenolics and methionine. Highly downregulated leaf transcripts at the hottest time of the day (40-42.4 °C) included genes encoding cyclin, cytochrome p450/secologanin synthase and U-box containing proteins, while upregulated, abundant transcripts included genes encoding heat shock proteins (HSPs), chaperones, UDP-glycosyltransferase, aquaporins and protein transparent testa 12. The upregulation of transcripts encoding HSPs, chaperones and UDP-glucosyltransferase and downregulation of transcripts encoding U-box containing proteins likely contributed to thermotolerance in R. stricta leaf by correcting protein folding and preventing protein degradation. Transcription factors that may regulate expression of genes encoding HSPs and chaperones under heat stress included HSFA2 to 4, AP2-EREBP and WRKY27., Conclusion: This study contributed new insights into the regulatory mechanisms of thermotolerance in the wild plant species R. stricta, an arid land, perennial evergreen shrub common in the Arabian Peninsula and Indian subcontinent. Enzymes from several pathways are interacting in the biosynthesis of soluble sugars, polyols, secondary metabolites, phenolics and methionine and are the primary contributors to thermotolerance in this species.

Alkaloid accumulation in plants is activated in response to stress, is limited in distribution and specific alkaloid repertoires are variable across taxa. Rauvolfioideae (Apocynaceae, Gentianales) represents a major center of structural expansion in the monoterpenoid indole alkaloids (MIAs) yielding thousands of unique molecules including highly valuable chemotherapeutics. The paucity of genome-level data for Apocynaceae precludes a deeper understanding of MIA pathway evolution hindering the elucidation of remaining pathway enzymes and the improvement of MIA availability in planta or in vitro. We sequenced the nuclear genome of Rhazya stricta (Apocynaceae, Rauvolfioideae) and present this high quality assembly in comparison with that of coffee (Rubiaceae, Coffea canephora, Gentianales) and others to investigate the evolution of genome-scale features. The annotated Rhazya genome was used to develop the community resource, RhaCyc, a metabolic pathway database. Gene family trees were constructed to identify homologs of MIA pathway genes and to examine their evolutionary history. We found that, unlike Coffea, the Rhazya lineage has experienced many structural rearrangements. Gene tree analyses suggest recent, lineage-specific expansion and diversification among homologs encoding MIA pathway genes in Gentianales and provide candidate sequences with the potential to close gaps in characterized pathways and support prospecting for new MIA production avenues.

Plastids and mitochondria have many protein complexes that include subunits encoded by organelle and nuclear genomes. In animal cells, compensatory evolution between mitochondrial and nuclear-encoded subunits was identified and the high mitochondrial mutation rates were hypothesized to drive compensatory evolution in nuclear genomes. In plant cells, compensatory evolution between plastid and nucleus has rarely been investigated in a phylogenetic framework. To investigate plastid-nuclear coevolution, we focused on plastid ribosomal protein genes that are encoded by plastid and nuclear genomes from 27 Geraniales species. Substitution rates were compared for five sets of genes representing plastid- and nuclear-encoded ribosomal subunit proteins targeted to the cytosol or the plastid as well as nonribosomal protein controls. We found that nonsynonymous substitution rates (dN) and the ratios of nonsynonymous to synonymous substitution rates (ω) were accelerated in both plastid- (CpRP) and nuclear-encoded subunits (NuCpRP) of the plastid ribosome relative to control sequences. Our analyses revealed strong signals of cytonuclear coevolution between plastid- and nuclear-encoded subunits, in which nonsynonymous substitutions in CpRP and NuCpRP tend to occur along the same branches in the Geraniaceae phylogeny. This coevolution pattern cannot be explained by physical interaction between amino acid residues. The forces driving accelerated coevolution varied with cellular compartment of the sequence. Increased ω in CpRP was mainly due to intensified positive selection whereas increased ω in NuCpRP was caused by relaxed purifying selection. In addition, the many indels identified in plastid rRNA genes in Geraniaceae may have contributed to changes in plastid subunits., (© The Author 2016. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Background and Aims: Several unrelated lineages such as plastids, viruses and plasmids, have converged on quadripartite genomes of similar size with large and small single copy regions and a large inverted repeat (IR). Except for Erodium (Geraniaceae), saguaro cactus and some legumes, the plastomes of all photosynthetic angiosperms display this structure. The functional significance of the IR is not understood and Erodium provides a system to examine the role of the IR in the long-term stability of these genomes. We compared the degree of genomic rearrangement in plastomes of Erodium that differ in the presence and absence of the IR., Methods: We sequenced 17 new Erodium plastomes. Using 454, Illumina, PacBio and Sanger sequences, 16 genomes were assembled and categorized along with one incomplete and two previously published Erodium plastomes. We conducted phylogenetic analyses among these species using a dataset of 19 protein-coding genes and determined if significantly higher evolutionary rates had caused the long branch seen previously in phylogenetic reconstructions within the genus. Bioinformatic comparisons were also performed to evaluate plastome evolution across the genus., Key Results: Erodium plastomes fell into four types (Type 1-4) that differ in their substitution rates, short dispersed repeat content and degree of genomic rearrangement, gene and intron content and GC content. Type 4 plastomes had significantly higher rates of synonymous substitutions (dS) for all genes and for 14 of the 19 genes non-synonymous substitutions (dN) were significantly accelerated. We evaluated the evidence for a single IR loss in Erodium and in doing so discovered that Type 4 plastomes contain a novel IR., Conclusions: The presence or absence of the IR does not affect plastome stability in Erodium. Rather, the overall repeat content shows a negative correlation with genome stability, a pattern in agreement with other angiosperm groups and recent findings on genome stability in bacterial endosymbionts., (© The Author 2016. Published by Oxford University Press on behalf of the Annals of Botany Company. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)