Characterization of the organellar genomes of Mazzaella laminarioides and Mazzaella membranacea (Gigartinaceae, Rhodophyta).

Mazzaella, a genus with no genomic resources available, has extensive distribution in the cold waters of the Pacific, where they represent ecologically and economically important species. In this study, we aimed to sequence, assemble, and annotate the complete mitochondrial and chloroplast genomes from two Mazzaella spp. and characterize the intraspecific variation among them. We report for the first time seven whole organellar genomes (mitochondria: OR915856, OR947465, OR947466, OR947467, OR947468, OR947469, OR947470; chloroplast: OR881974, OR909680, OR909681, OR909682, OR909683, OR909684, OR909685) obtained through high-throughput sequencing for six M. laminarioides sampled from three Chilean regions and one M. membranacea. Sequenced Mazzaella mitogenomes have identical gene number, gene order, and genome structure. The same results were observed for assembled plastomes. A total of 52 genes were identified in mitogenomes, and a total of 235 genes were identified in plastomes. Although the M. membranacea plastome included a full-length pbsA gene, in all M. laminarioides samples, the pbsA gene was split in three open reading frames (ORFs). Within M. laminarioides, we observed important plastome lineage-specific variations, such as the pseudogenization of the two hypothetical protein-coding genes, ycf23 and ycf45. Nonsense mutations in the ycf23 and ycf45 genes were only detected in the northern lineage. These results are consistent with phylogenetic reconstructions and divergence time estimation using concatenated coding sequences that not only support the monophyly of M. laminarioides but also underscore that the three M. laminarioides lineages are in an advanced stage of divergence. These new results open the question of the existence of still undisclosed species in M. laminarioides.

Organellar-genome analyses from the lycophyte genus Isoetes L. show one of the highest frequencies of RNA editing in land plants.

RNA editing is a post-transcriptional process that challenges the central dogma of molecular biology by modifying RNA sequences, introducing nucleotide changes at specific sites, and generating functional diversity beyond the genomic code, especially when it concerns organellar transcripts. In plants, this phenomenon is widespread, but its extent varies significantly among species and organellar genomes. Among land plants, the heterosporous lycophytes (i.e., Isoetes and Selaginella) stand out for their exceptionally high numbers of RNA-editing sites, despite their morphological stasis and ancient lineage. In this study, we explore the complete set of organellar protein-coding genes in the aquatic plant group Isoetes, providing a detailed analysis of RNA editing in both the mitochondrial and plastid genomes. Our findings reveal a remarkable abundance of RNA editing, particularly in the mitochondrial genome, with thousands of editing sites identified. Interestingly, the majority of these edits result in non-silent substitutions, suggesting a role in fine-tuning protein structure and function. Furthermore, we observe a consistent trend of increased hydrophobicity in membrane-bound proteins, supporting the notion that RNA editing may confer a selective advantage by preserving gene functionality in Isoetes. The conservation of highly edited RNA sequences over millions of years underscores the evolutionary significance of RNA editing. Additionally, the study sheds light on the dynamic nature of RNA editing, with shared editing sites reflecting common ancestry whereas exclusive edits matching more recent radiation events within the genus. This work advances our understanding of the intricate interplay between RNA editing, adaptation, and evolution in land plants and highlights the unique genomic features of Isoetes, providing a foundation for further investigations into the functional consequences of RNA editing in this enigmatic plant lineage.

Chloroplast genome assembly of Serjania erecta Raldk: comparative analysis reveals gene number variation and selection in protein-coding plastid genes of Sapindaceae.

Serjania erecta Raldk is an essential genetic resource due to its anti-inflammatory, gastric protection, and anti-Alzheimer properties. However, the genetic and evolutionary aspects of the species remain poorly known. Here, we sequenced and assembled the complete chloroplast genome of S. erecta and used it in a comparative analysis within the Sapindaceae family. S. erecta has a chloroplast genome (cpDNA) of 159,297 bp, divided into a Large Single Copy region (LSC) of 84,556 bp and a Small Single Copy region (SSC) of 18,057 bp that are surrounded by two Inverted Repeat regions (IRa and IRb) of 28,342 bp. Among the 12 species used in the comparative analysis, S. erecta has the fewest long and microsatellite repeats. The genome structure of Sapindaceae species is relatively conserved; the number of genes varies from 128 to 132 genes, and this variation is associated with three main factors: (1) Expansion and retraction events in the size of the IRs, resulting in variations in the number of rpl22, rps19, and rps3 genes; (2) Pseudogenization of the rps2 gene; and (3) Loss or duplication of genes encoding tRNAs, associated with the duplication of trnH-GUG in X. sorbifolium and the absence of trnT-CGU in the Dodonaeoideae subfamily. We identified 10 and 11 mutational hotspots for Sapindaceae and Sapindoideae, respectively, and identified six highly diverse regions (tRNA-Lys - rps16, ndhC - tRNA-Val, petA - psbJ, ndhF, rpl32 - ccsA, and ycf1) are found in both groups, which show potential for the development of DNA barcode markers for molecular taxonomic identification of Serjania. We identified that the psaI gene evolves under neutrality in Sapindaceae, while all other chloroplast genes are under strong negative selection. However, local positive selection exists in the ndhF, rpoC2, ycf1, and ycf2 genes. The genes ndhF and ycf1 also present high nucleotide diversity and local positive selection, demonstrating significant potential as markers. Our findings include providing the first chloroplast genome of a member of the Paullinieae tribe. Furthermore, we identified patterns in variations in the number of genes and selection in genes possibly associated with the family's evolutionary history.

A comprehensive genus-level phylogeny and biogeographical history of the Lythraceae based on whole plastome sequences.

BACKGROUND AND AIMS: The Lythraceae are a mainly subtropical to tropical family of the order Myrtales with 28 currently accepted genera and approximately 600 species. There is currently no well-supported phylogenetic and biogeographical hypothesis of the Lythraceae incorporating all currently accepted genera, which we sought to provide. METHODS: Plastomes of representative species of 18 distinct Lythraceae genera were sequenced and annotated. Together with existing sequences, plastomes of all 28 currently accepted genera in the Lythraceae were brought together for the first time. The plastomes were aligned and a Bayesian phylogenetic hypothesis was produced. We then conducted a time-calibrated Bayesian analysis and a biogeographical analysis. KEY RESULTS: Plastome-based Bayesian and maximum-likelihood phylogenetic trees are generally congruent with recent nuclear phylogenomic data and resolve two deeply branching major clades in the Lythraceae. One major clade concentrates shrubby and arboreal South American and African genera that inhabit seasonally dry environments, with larger, often winged seeds, adapted to dispersal by the wind. The second major clade concentrates North American, Asian, African and several near-cosmopolitan herbaceous, shrubby and arboreal genera, often inhabiting humid or aquatic environments, with smaller seeds possessing structures that facilitate dispersal by water. CONCLUSIONS: We hypothesize that the Lythraceae dispersed early in the Late Cretaceous from South American to North American continents, with subsequent expansion in the Late Cretaceous of a North American lineage through Laurasia to Africa via a boreotropical route. Two later expansions of South American clades to Africa in the Palaeocene and Eocene, respectively, are also hypothesized. Transoceanic dispersal in the family is possibly facilitated by adaptations to aquatic environments that are common to many extant genera of the Lythraceae, where long-distance dispersal and vicariance may be invoked to explain several remarkable disjunct distributions in Lythraceae clades.

Phylogenomics and plastome evolution of a Brazilian mycoheterotrophic orchid, Pogoniopsis schenckii.

PREMISE: Pogoniopsis likely represents an independent photosynthesis loss in orchids. We use phylogenomic data to better identify the phylogenetic placement of this fully mycoheterotrophic taxon, and investigate its molecular evolution. METHODS: We performed likelihood analysis of plastid and mitochondrial phylogenomic data to localize the position of Pogoniopsis schenckii in orchid phylogeny, and investigated the evolution of its plastid genome. RESULTS: All analyses place Pogoniopsis in subfamily Epidendroideae, with strongest support from mitochondrial data, which also place it near tribe Sobralieae with moderately strong support. Extreme rate elevation in Pogoniopsis plastid genes broadly depresses branch support; in contrast, mitochondrial genes are only mildly rate elevated and display very modest and localized reductions in bootstrap support. Despite considerable genome reduction, including loss of photosynthesis genes and multiple translation apparatus genes, gene order in Pogoniopsis plastomes is identical to related autotrophs, apart from moderately shifted inverted repeat (IR) boundaries. All cis-spliced introns have been lost in retained genes. Two plastid genes (accD, rpl2) show significant strengthening of purifying selection. A retained plastid tRNA gene (trnE-UUC) of Pogoniopsis lacks an anticodon; we predict that it no longer functions in translation but retains a secondary role in heme biosynthesis. CONCLUSIONS: Slowly evolving mitochondrial genes clarify the placement of Pogoniopsis in orchid phylogeny, a strong contrast with analysis of rate-elevated plastome data. We documented the effects of the novel loss of photosynthesis: for example, despite massive gene loss, its plastome is fully colinear with other orchids, and it displays only moderate shifts in selective pressure in retained genes.

Bactris gasipaes Kunth var. gasipaes complete plastome and phylogenetic analysis.

Bactris gasipaes var. gasipaes (Arecaceae, Palmae) is an economically and socially important plant species for populations across tropical South and Central America. It has been domesticated from its wild variety, B. gasipaes var. chichagui, since pre-Columbian times. In this study, we sequenced the plastome of the cultivated variety, B. gasipaes Kunth var. gasipaes and compared it with the published plastome of the wild variety. The chloroplast sequence obtained was 156,580 bp. The cultivated chloroplast sequence was conserved compared to the wild type sequence with 99.8% of nucleotide identity. We did, however, identify multiple Single Nucleotide Variants (SNVs), insertions, microsatellites and a resolved region of missing nucleotides. A SNV in one of the core barcode markers (matK) was detected between the wild and cultivated accessions. Phylogenetic analysis was carried out across the Arecaceae family and compared to previous reports, resulting in an identical topology. This study is a step forward in understanding the genome evolution of this species.

New plastome structural rearrangements discovered in core Tillandsioideae (Bromeliaceae) support recently adopted taxonomy.

Full plastome sequences for land plants have become readily accessible thanks to the development of Next Generation Sequencing (NGS) techniques and powerful bioinformatic tools. Despite this vast amount of genomic data, some lineages remain understudied. Full plastome sequences from the highly diverse (>1,500 spp.) subfamily Tillandsioideae (Bromeliaceae, Poales) have been published for only three (i.e., Guzmania, Tillandsia, and Vriesea) out of 22 currently recognized genera. Here, we focus on core Tillandsioideae, a clade within subfamily Tillandsioideae, and explore the contribution of individual plastid markers and data categories to inform deep divergences of a plastome phylogeny. We generated 37 high quality plastome assemblies and performed a comparative analysis in terms of plastome structure, size, gene content and order, GC content, as well as number and type of repeat motifs. Using the obtained phylogenetic context, we reconstructed the evolution of these plastome attributes and assessed if significant shifts on the evolutionary traits' rates have occurred in the evolution of the core Tillandsioideae. Our results agree with previously published phylogenetic hypotheses based on plastid data, providing stronger statistical support for some recalcitrant nodes. However, phylogenetic discordance with previously published nuclear marker-based hypotheses was found. Several plastid markers that have been consistently used to address phylogenetic relationships within Tillandsioideae were highly informative for the retrieved plastome phylogeny and further loci are here identified as promising additional markers for future studies. New lineage-specific plastome rearrangements were found to support recently adopted taxonomic groups, including large inversions, as well as expansions and contractions of the inverted repeats. Evolutionary trait rate shifts associated with changes in size and GC content of the plastome regions were found across the phylogeny of core Tillandsioideae.

Uncovering haplotype diversity in cultivated Mexican vanilla species.

PREMISE: Although vanilla is one of the best-known spices, there is only a limited understanding of its biology and genetics within Mexico, where its cultivation originated and where phenotypic variability is high. This study aims to augment our understanding of vanilla's genetic resources by assessing species delimitation and genetic, geographic, and climatic variability within Mexican cultivated vanilla. METHODS: Using nuclear and plastid DNA sequence data from 58 Mexican samples collected from three regions and 133 ex situ accessions, we assessed species monophyly using phylogenetic analyses and genetic distances. Intraspecific genetic variation was summarized through the identification of haplotypes. Within the primarily cultivated species, Vanilla planifolia, haplotype relationships were further verified using plastome and rRNA gene sequences. Climatic niche and haplotype composition were assessed across the landscape. RESULTS: Three species (Vanilla planifolia, V. pompona, and V. insignis) and 13 haplotypes were identified among Mexican vanilla. Within V. planifolia haplotypes, hard phylogenetic incongruences between plastid and nuclear sequences suggest past hybridization events. Eight haplotypes consisted exclusively of Mexican samples. The dominant V. planifolia haplotype occurred throughout all three regions as well as outside of its country of origin. Haplotype richness was found to be highest in regions around Papantla and Chinantla. CONCLUSIONS: Long histories of regional cultivation support the consideration of endemic haplotypes as landraces shaped by adaptation to local conditions and/or hybridization. Results may aid further genomic investigations of vanilla's genetic resources and ultimately support the preservation of genetic diversity within the economically important crop.

Putting small and big pieces together: a genome assembly approach reveals the largest Lamiid plastome in a woody vine.

The plastid genome of flowering plants generally shows conserved structural organization, gene arrangement, and gene content. While structural reorganizations are uncommon, examples have been documented in the literature during the past years. Here we assembled the entire plastome of Bignonia magnifica and compared its structure and gene content with nine other Lamiid plastomes. The plastome of B. magnifica is composed of 183,052 bp and follows the canonical quadripartite structure, synteny, and gene composition of other angiosperms. Exceptionally large inverted repeat (IR) regions are responsible for the uncommon length of the genome. At least four events of IR expansion were observed among the seven Bignoniaceae species compared, suggesting multiple expansions of the IRs over the SC regions in the family. A comparison with 6,231 other complete plastomes of flowering plants available on GenBank revealed that the plastome of B. magnifica is the longest Lamiid plastome described to date. The newly generated plastid genome was used as a source of selected genes. These genes were combined with orthologous regions sampled from other species of Bignoniaceae and all gene alignments concatenated to infer a phylogeny of the family. The tree recovered is consistent with known relationships within the Bignoniaceae.

The Barley Chloroplast Mutator (cpm) Mutant: All Roads Lead to the Msh1 Gene.

The barley chloroplast mutator (cpm) is a nuclear gene mutant that induces a wide spectrum of cytoplasmically inherited chlorophyll deficiencies. Plastome instability of cpm seedlings was determined by identification of a particular landscape of polymorphisms that suggests failures in a plastome mismatch repair (MMR) protein. In Arabidopsis, MSH genes encode proteins that are in charge of mismatch repair and have anti-recombination activity. In this work, barley homologs of these genes were identified, and their sequences were analyzed in control and cpm mutant seedlings. A substitution, leading to a premature stop codon and a truncated MSH1 protein, was identified in the Msh1 gene of cpm plants. The relationship between this mutation and the presence of chlorophyll deficiencies was established in progenies from crosses and backcrosses. These results strongly suggest that the mutation identified in the Msh1 gene of the cpm mutant is responsible for the observed plastome instabilities. Interestingly, comparison of mutant phenotypes and molecular changes induced by the barley cpm mutant with those of Arabidopsis MSH1 mutants revealed marked differences.

The genome of a thorny species: comparative genomic analysis among South and North American Cactaceae.

MAIN CONCLUSION: The first South American cactus nuclear genome assembly associated with comparative genomic analyses provides insights into nuclear and plastid genomic features, such as size, transposable elements, and metabolic processes related to cactus development. Here, we assembled the partial genome, plastome, and transcriptome of Cereus fernambucensis (Cereeae, Cactaceae), a representative species of the South American core Cactoideae. We accessed other genomes and transcriptomes available for cactus species to compare the heterozygosity level, genome size, transposable elements, orthologous genes, and plastome structure. These estimates were obtained from the literature or using the same pipeline adopted for C. fermabucensis. In addition to the C. fernambucensis plastome, we also performed de novo plastome assembly of Pachycereus pringlei, Stenocereus thurberi, and Pereskia humboldtii based on the sequences available in public databases. We estimated a genome size of ~ 1.58 Gb for C. fernambucensis, the largest genome among the compared species. The genome heterozygosity was 0.88% in C. fernambucensis but ranged from 0.36 (Carnegiea gigantea) to 17.4% (Lophocereus schottii) in the other taxa. The genome lengths of the studied cacti are constituted by a high amount of transposable elements, ranging from ~ 57 to ~ 67%. Putative satellite DNAs are present in all species, excepting C. gigantea. The plastome of C. fernambucensis was ~ 104 kb, showing events of translocation, inversion, and gene loss. We observed a low number of shared unique orthologs, which may suggest gene duplication events and the simultaneous expression of paralogous genes. We recovered 37 genes that have undergone positive selection along the Cereus branch that are associated with different metabolic processes, such as improving photosynthesis during drought stress and nutrient absorption, which may be related to the adaptation to xeric areas of the Neotropics.

The complete chloroplast genome of Artocarpus altilis (Moraceae) and phylogenetic relationships.

The chloroplast (cp) is an essential organelle in higher plants. The genes of the plastome are well suited to infer phylogenetic relationships among taxa. In this study, we report the assembly of the cp genome of Artocarpus altilis and its phylogeny among species from Moraceae family. The cp genome of A. altilis was 160,822 base pair (bp) in length, comprising one large single-copy region of 88,692 bp, one small single-copy region of 19,290 bp, and a pair of inverted repeat regions (IRs) of 26,420 bp. A total of 113 different genes were predicted, including 79 protein-coding genes, 30 tRNA genes, and four rRNA genes. The phylogenetic analysis of 19 species belonging to the Moraceae family confirmed the phylogenetic proximity of the genus Artocarpus and Morus and the genetic similarity of A. camansi and A. altilis.

Extremely low nucleotide diversity among thirty-six new chloroplast genome sequences from Aldama (Heliantheae, Asteraceae) and comparative chloroplast genomics analyses with closely related genera.

Aldama (Heliantheae, Asteraceae) is a diverse genus in the sunflower family. To date, nearly 200 Asteraceae chloroplast genomes have been sequenced, but the plastomes of Aldama remain undescribed. Plastomes in Asteraceae usually show little sequence divergence, consequently, our hypothesis is that species of Aldama will be overall conserved. In this study, we newly sequenced 36 plastomes of Aldama and of five species belonging to other Heliantheae genera selected as outgroups (i.e., Dimerostemma asperatum, Helianthus tuberosus, Iostephane heterophylla, Pappobolus lanatus var. lanatus, and Tithonia diversifolia). We analyzed the structure and gene content of the assembled plastomes and performed comparative analyses within Aldama and with other closely related genera. As expected, Aldama plastomes are very conserved, with the overall gene content and orientation being similar in all studied species. The length of the plastome is also consistent and the junction between regions usually contain the same genes and have similar lengths. A large ∼20 kb and a small ∼3 kb inversion were detected in the Large Single Copy (LSC) regions of all assembled plastomes, similarly to other Asteraceae species. The nucleotide diversity is very low, with only 1,509 variable sites in 127,466 bp (i.e., 1.18% of the sites in the alignment of 36 Aldama plastomes, with one of the IRs removed, is variable). Only one gene, rbcL, shows signatures of positive selection. The plastomes of the selected outgroups feature a similar gene content and structure compared to Aldama and also present the two inversions in the LSC region. Deletions of different lengths were observed in the gene ycf2. Multiple SSRs were identified for the sequenced Aldama and outgroups. The phylogenetic analysis shows that Aldama is not monophyletic due to the position of the Mexican species A. dentata. All Brazilian species form a strongly supported clade. Our results bring new understandings into the evolution and diversity of plastomes at the species level.

Plastome genomics in South American maize landraces: chloroplast lineages parallel the geographical structuring of nuclear gene pools.

BACKGROUND AND AIMS: The number of plastome sequences has increased exponentially during the last decade. However, there is still little knowledge of the levels and distribution of intraspecific variation. The aims of this study were to estimate plastome diversity within Zea mays and analyse the distribution of haplotypes in connection with the landrace groups previously delimited for South American maize based on nuclear markers. METHODS: We obtained the complete plastomes of 30 South American maize landraces and three teosintes by means of next-generation sequencing (NGS) and used them in combination with data from public repositories. After quality filtering, the curated data were employed to search for single-nucleotide polymorphisms, indels and chloroplast simple sequence repeats. Exact permutational contingency tests were performed to assess associations between plastome and nuclear variation. Network and Bayesian phylogenetic analyses were used to infer evolutionary relationships among haplotypes. KEY RESULTS: Our analyses identified a total of 124 polymorphic plastome loci, with the intergenic regions psbE-rps18, petN-rpoB, trnL_UAG-ndhF and rpoC2-atpI exhibiting the highest marker densities. Although restricted in number, these markers allowed the discrimination of 27 haplotypes in a total of 51 Zea mays individuals. Andean and lowland South American landraces differed significantly in haplotype distribution. However, overall differentiation patterns were not informative with respect to subspecies diversification, as evidenced by the scattered distribution of maize and teosinte plastomes in both the network and Bayesian phylogenetic reconstructions. CONCLUSIONS: Knowledge of intraspecific plastome variation provides the framework for a more comprehensive understanding of evolutionary processes at low taxonomic levels and may become increasingly important for future plant barcoding efforts. Whole-plastome sequencing provided useful variability to contribute to maize phylogeographic studies. The structuring of haplotype diversity in the maize landraces examined here clearly reflects the distinction between the Andean and South American lowland gene pools previously inferred based on nuclear markers.

Chloroplast genome assembly of Handroanthus impetiginosus: comparative analysis and molecular evolution in Bignoniaceae.

MAIN CONCLUSION: Bignoniaceae species have conserved chloroplast structure, with hotspots of nucleotide diversity. Several genes are under positive selection, and can be targets for evolutionary studies. Bignoniaceae is one of the most species-rich family of woody plants in Neotropical seasonally dry forests. Here we report the assembly of Handroanthus impetiginosus chloroplast genome and evolutionary comparative analyses of ten Bignoniaceae species representing the genera for which whole-genome chloroplast sequences were available. The chloroplast genome of H. impetiginosus is 159,462 bp in size and has a similar structure compared to the other nine species. The total number of genes was slightly variable amongst the Bignoniaceae, ranging from 124 in H. impetiginosus to 144 in Anemopaegma acutifolium. The inverted repeat (IR) size was variable, ranging from 24,657 bp (Tecomaria capensis) to 40,481 bp (A. acutifolium), due to the contraction and retraction at its boundaries. However, gene boundaries were very similar among the ten species. We found 98 forward and palindromic dispersed repeats, and 85 simple sequence repeats (SSRs). In general, chloroplast sequences were highly conserved, with few nucleotide diversity hotspots in the genes accD, clpP, rpoA, ycf1, ycf2. The phylogenetic analysis based on 77 coding genes was highly consistent with Angiosperm Phylogeny Group (APG) IV. Our results also indicate that most genes are under negative selection or neutral evolution. We found no evidence of branch-site selection, implying that H. impetiginosus is not evolving faster than the other species analyzed, notwithstanding we found site positive selection signal in several genes. These genes can provide targets for evolutionary studies in Bignoniaceae and Lamiales species.

Insights Into Chloroplast Genome Evolution Across Opuntioideae (Cactaceae) Reveals Robust Yet Sometimes Conflicting Phylogenetic Topologies.

Chloroplast genomes (plastomes) are frequently treated as highly conserved among land plants. However, many lineages of vascular plants have experienced extensive structural rearrangements, including inversions and modifications to the size and content of genes. Cacti are one of these lineages, containing the smallest plastome known for an obligately photosynthetic angiosperm, including the loss of one copy of the inverted repeat (∼25 kb) and the ndh gene suite, but only a few cacti from the subfamily Cactoideae have been sufficiently characterized. Here, we investigated the variation of plastome sequences across the second-major lineage of the Cactaceae, the subfamily Opuntioideae, to address (1) how variable is the content and arrangement of chloroplast genome sequences across the subfamily, and (2) how phylogenetically informative are the plastome sequences for resolving major relationships among the clades of Opuntioideae. Our de novo assembly of the Opuntia quimilo plastome recovered an organelle of 150,347 bp in length with both copies of the inverted repeat and the presence of all the ndh gene suite. An expansion of the large single copy unit and a reduction of the small single copy unit was observed, including translocations and inversion of genes, as well as the putative pseudogenization of some loci. Comparative analyses among all clades within Opuntioideae suggested that plastome structure and content vary across taxa of this subfamily, with putative independent losses of the ndh gene suite and pseudogenization of genes across disparate lineages, further demonstrating the dynamic nature of plastomes in Cactaceae. Our plastome dataset was robust in resolving three tribes with high support within Opuntioideae: Cylindropuntieae, Tephrocacteae and Opuntieae. However, conflicting topologies were recovered among major clades when exploring different assemblies of markers. A plastome-wide survey for highly informative phylogenetic markers revealed previously unused regions for future use in Sanger-based studies, presenting a valuable dataset with primers designed for continued evolutionary studies across Cactaceae. These results bring new insights into the evolution of plastomes in cacti, suggesting that further analyses should be carried out to address how ecological drivers, physiological constraints and morphological traits of cacti may be related with the common rearrangements in plastomes that have been reported across the family.

A Repertory of Rearrangements and the Loss of an Inverted Repeat Region in Passiflora Chloroplast Genomes.

Chloroplast genomes (cpDNA) in angiosperms are usually highly conserved. Although rearrangements have been observed in some lineages, such as Passiflora, the mechanisms that lead to rearrangements are still poorly elucidated. In the present study, we obtained 20 new chloroplast genomes (18 species from the genus Passiflora, and Dilkea retusa and Mitostemma brevifilis from the family Passifloraceae) in order to investigate cpDNA evolutionary history in this group. Passiflora cpDNAs vary in size considerably, with ∼50 kb between shortest and longest. Large inverted repeat (IR) expansions were identified, and at the extreme opposite, the loss of an IR was detected for the first time in Passiflora, a rare event in angiosperms. The loss of an IR region was detected in Passiflora capsularis and Passiflora costaricensis, a species in which occasional biparental chloroplast inheritance has previously been reported. A repertory of rearrangements such as inversions and gene losses were detected, making Passiflora one of the few groups with complex chloroplast genome evolution. We also performed a phylogenomic study based on all the available cp genomes and our analysis implies that there is a need to reconsider the taxonomic classifications of some species in the group.

ECuADOR-Easy Curation of Angiosperm Duplicated Organellar Regions, a tool for cleaning and curating plastomes assembled from next generation sequencing pipelines.

BACKGROUND: With the rapid increase in availability of genomic resources offered by Next-Generation Sequencing (NGS) and the availability of free online genomic databases, efficient and standardized metadata curation approaches have become increasingly critical for the post-processing stages of biological data. Especially in organelle-based studies using circular chloroplast genome datasets, the assembly of the main structural regions in random order and orientation represents a major limitation in our ability to easily generate "ready-to-align" datasets for phylogenetic reconstruction, at both small and large taxonomic scales. In addition, current practices discard the most variable regions of the genomes to facilitate the alignment of the remaining coding regions. Nevertheless, no software is currently available to perform curation to such a degree, through simple detection, organization and positioning of the main plastome regions, making it a time-consuming and error-prone process. Here we introduce a fast and user friendly software ECuADOR, a Perl script specifically designed to automate the detection and reorganization of newly assembled plastomes obtained from any source available (NGS, sanger sequencing or assembler output). METHODS: ECuADOR uses a sliding-window approach to detect long repeated sequences in draft sequences, which then identifies the inverted repeat regions (IRs), even in case of artifactual breaks or sequencing errors and automates the rearrangement of the sequence to the widely used LSC-Irb-SSC-IRa order. This facilitates rapid post-editing steps such as creation of genome alignments, detection of variable regions, SNP detection and phylogenomic analyses. RESULTS: ECuADOR was successfully tested on plant families throughout the angiosperm phylogeny by curating 161 chloroplast datasets. ECuADOR first identified and reordered the central regions (LSC-Irb-SSC-IRa) for each dataset and then produced a new annotation for the chloroplast sequences. The process took less than 20 min with a maximum memory requirement of 150 MB and an accuracy of over 99%. CONCLUSIONS: ECuADOR is the sole de novo one-step recognition and re-ordination tool that provides facilitation in the post-processing analysis of the extra nuclear genomes from NGS data. The program is available at https://github.com/BiodivGenomic/ECuADOR/.