Search

Repeats can mediate rearrangements and recombination in plant mitochondrial genomes and plastid genomes. While repeat accumulations are linked to heightened evolutionary rates and complex structures in specific lineages, debates persist regarding the extent of their influence on sequence and structural evolution. In this study, 75 Plantago plastomes were analyzed to investigate the relationships between repeats, nucleotide substitution rates, and structural variations. Extensive repeat accumulations were associated with significant rearrangements and inversions in the large inverted repeats (IRs), suggesting that repeats contribute to rearrangement hotspots. Repeats caused infrequent recombination that potentially led to substoichiometric shifting, supported by long-read sequencing. Repeats were implicated in elevating evolutionary rates by facilitating localized hypermutation, likely through DNA damage and repair processes. This study also observed a decrease in nucleotide substitution rates for loci translocating into IRs, supporting the role of biased gene conversion in maintaining lower substitution rates. Combined with known parallel changes in mitogenomes, it is proposed that potential dysfunction in nuclear-encoded genes associated with DNA replication, recombination, and repair may drive the evolution of Plantago organellar genomes. These findings contribute to understanding how repeats impact organellar evolution and stability, particularly in rapidly evolving plant lineages. [ABSTRACT FROM AUTHOR]

Abstract Background Scaffolding is an intermediate stage of fragment assembly. It consists in orienting and ordering the contigs obtained by the assembly of the sequencing reads. In the general case, the problem has been largely studied with the use of distances data between the contigs. Here we focus on a dedicated scaffolding for the chloroplast genomes. As these genomes are small, circular and with few specific repeats, numerous approaches have been proposed to assemble them. However, their specificities have not been sufficiently exploited. Results We give a new formulation for the scaffolding in the case of chloroplast genomes as a discrete optimisation problem, that we prove the decision version to be $$\mathcal{NP}$$ NP -Complete. We take advantage of the knowledge of chloroplast genomes and succeed in expressing the relationships between a few specific genomic repeats in mathematical constraints. Our approach is independent of the distances and adopts a genomic regions view, with the priority on scaffolding the repeats first. In this way, we encode the structural haplotype issue in order to retrieve several genome forms that coexist in the same chloroplast cell. To solve exactly the optimisation problem, we develop an integer linear program that we implement in Python3 package khloraascaf. We test it on synthetic data to investigate its performance behaviour and its robustness against several chosen difficulties. Conclusions We succeed to model biological knowledge on genomic structures to scaffold chloroplast genomes. Our results suggest that modelling genomic regions is sufficient for scaffolding repeats and is suitable for finding several solutions corresponding to several genome forms.

Background: Scaffolding is an intermediate stage of fragment assembly. It consists in orienting and ordering the contigs obtained by the assembly of the sequencing reads. In the general case, the problem has been largely studied with the use of distances data between the contigs. Here we focus on a dedicated scaffolding for the chloroplast genomes. As these genomes are small, circular and with few specific repeats, numerous approaches have been proposed to assemble them. However, their specificities have not been sufficiently exploited. Results: We give a new formulation for the scaffolding in the case of chloroplast genomes as a discrete optimisation problem, that we prove the decision version to be NP -Complete. We take advantage of the knowledge of chloroplast genomes and succeed in expressing the relationships between a few specific genomic repeats in mathematical constraints. Our approach is independent of the distances and adopts a genomic regions view, with the priority on scaffolding the repeats first. In this way, we encode the structural haplotype issue in order to retrieve several genome forms that coexist in the same chloroplast cell. To solve exactly the optimisation problem, we develop an integer linear program that we implement in Python3 package khloraascaf. We test it on synthetic data to investigate its performance behaviour and its robustness against several chosen difficulties. Conclusions: We succeed to model biological knowledge on genomic structures to scaffold chloroplast genomes. Our results suggest that modelling genomic regions is sufficient for scaffolding repeats and is suitable for finding several solutions corresponding to several genome forms. [ABSTRACT FROM AUTHOR]

The insect mitogenome is typically a compact circular molecule with highly conserved gene contents. Nonetheless, mitogenome structural variations have been reported in specific taxa, and gene rearrangements, usually the tRNAs, occur in different lineages. Because synapomorphies of mitogenome organizations can provide information for phylogenetic inferences, comparative analyses of mitogenomes have been given increasing attention. However, most studies use a very few species to represent the whole genus, tribe, family, or even order, overlooking potential variations at lower taxonomic levels, which might lead to some incorrect inferences. To provide new insights into mitogenome organizations and their implications for phylogenetic inference, this study conducted comparative analyses for mitogenomes of three social bee tribes (Meliponini, Bombini, and Apini) based on the phylogenetic framework with denser taxonomic sampling at the species and population levels. Comparative analyses revealed that mitogenomes of Apini and Bombini are the typical type, while those of Meliponini show diverse variations in mitogenome sizes and organizations. Large inverted repeats (IRs) cause significant gene rearrangements of protein coding genes (PCGs) and rRNAs in Indo-Malay/Australian stingless bee species. Molecular evolution analyses showed that the lineage with IRs have lower dN/dS ratios for PCGs than lineages without IRs, indicating potential effects of IRs on the evolution of mitochondrial genes. The finding of IRs and different patterns of gene rearrangements suggested that Meliponini is a hotspot in mitogenome evolution. Unlike conserved PCGs and rRNAs whose rearrangements were found only in the mentioned lineages within Meliponini, tRNA rearrangements are common across all three tribes of social bees, and are significant even at the species level, indicating that comprehensive sampling is needed to fully understand the patterns of tRNA rearrangements, and their implications for phylogenetic inference. [ABSTRACT FROM AUTHOR]

Bacterial plasmids are important mobile genetic elements which often carry specific genes important for bacterial successful survival under various inhospitable environmental conditions. Most of the previous research has focused on large plasmids providing these beneficial traits to their host cells. In this study, small cryptic plasmid pALK1 (3 051 bp) was isolated from metallotolerant and alkalitollerant strain Acinetobacter sp. K1. The plasmid encodes Rep_3 family replication protein and MobM mobilization protein but none pALK1-like plasmids were detected in other Acinetobacter strains of environmental and animal origin. The secondary structure of the pALK1 plasmid is characterized by the complexity of multiple sets of direct and inverted repeats in its nucleotide sequence. Comparative genomics was used to hypothesize the biological functions of these repeats in Acinetobacter spp., whereas several similar plasmids with a related organization of direct repeats and palindromes are known in this genus. [ABSTRACT FROM AUTHOR]

This work aims to describe the observed enrichment of inverted repeats in the human genome; and to identify and describe, with detailed length profiles, the regions with significant and relevant enriched occurrence of inverted repeats. The enrichment is assessed and tested with a recently proposed measure (z-scores based measure). We simulate a genome using an order 7 Markov model trained with the data from the real genome. The simulated genome is used to establish the critical values which are used as decision thresholds to identify the regions with significant enriched concentrations. Several human genome regions are highly enriched in the occurrence of inverted repeats. This is observed in all the human chromosomes. The distribution of inverted repeat lengths varies along the genome. The majority of the regions with severely exaggerated enrichment contain mainly short length inverted repeats. There are also regions with regular peaks along the inverted repeats lengths distribution (periodic regularities) and other regions with exaggerated enrichment for long lengths (less frequent). However, adjacent regions tend to have similar distributions.

To understand the evolution of coscinodiscophycean diatoms, plastome sequences of six coscinodiscophycean diatom species were constructed and analyzed in this study, doubling the number of constructed plastome sequences in Coscinodiscophyceae (radial centrics). The platome sizes varied substantially in Coscinodiscophyceae, ranging from 119.1 kb of Actinocyclus subtilis to 135.8 kb of Stephanopyxis turris. Plastomes in Paraliales and Stephanopyxales tended to be larger than those in Rhizosoleniales and Coscinodiacales, which were due to the expansion of the inverted repeats (IRs) and to the marked increase of the large single copy (LSC). Phylogenomic analysis indicated that Paralia and Stephanopyxis clustered tightly to form the Paraliales-Stephanopyxales complex, which was sister to the Rhizosoleniales-Coscinodiscales complex. The divergence time between Paraliales and Stephanopyxales was estimated at 85 MYA in the middle Upper Cretaceous, indicating that Paraliales and Stephanopyxales appeared later than Coscinodiacales and Rhizosoleniales according to their phylogenetic relationships. Frequent losses of housekeeping protein-coding genes (PCGs) were observed in these coscinodiscophycean plastomes, indicating that diatom plastomes showed an ongoing reduction in gene content during evolution. Two acpP genes (acpP1 and acpP2) detected in diatom plastomes were found to be originated from an early gene duplication event occurred in the common progenitor after diatom emergence, rather than multiple independent gene duplications occurring in different lineages of diatoms. The IRs in Stephanopyxis turris and Rhizosolenia fallax-imbricata exhibited a similar trend of large expansion to the small single copy (SSC) and slightly small contraction from the LSC, which eventually led to the conspicuous increase in IR size. Gene order was highly conserved in Coscinodiacales, while multiple rearrangements were observed in Rhizosoleniales and between Paraliales and Stephanopyxales. Our results greatly expanded the phylogenetic breadth in Coscinodiscophyceae and gained novel insights into the evolution of plastomes in diatoms. [ABSTRACT FROM AUTHOR]

Background: Aging in postmitotic tissues is associated with clonal expansion of somatic mitochondrial deletions, the origin of which is not well understood. Such deletions are often flanked by direct nucleotide repeats, but this alone does not fully explain their distribution. Here, we hypothesized that the close proximity of direct repeats on single-stranded mitochondrial DNA (mtDNA) might play a role in the formation of deletions. Results: By analyzing human mtDNA deletions in the major arc of mtDNA, which is single-stranded during replication and is characterized by a high number of deletions, we found a non-uniform distribution with a "hot spot" where one deletion breakpoint occurred within the region of 6–9 kb and another within 13–16 kb of the mtDNA. This distribution was not explained by the presence of direct repeats, suggesting that other factors, such as the spatial proximity of these two regions, can be the cause. In silico analyses revealed that the single-stranded major arc may be organized as a large-scale hairpin-like loop with a center close to 11 kb and contacting regions between 6–9 kb and 13–16 kb, which would explain the high deletion activity in this contact zone. The direct repeats located within the contact zone, such as the well-known common repeat with a first arm at 8470–8482 bp (base pair) and a second arm at 13,447–13,459 bp, are three times more likely to cause deletions compared to direct repeats located outside of the contact zone. A comparison of age- and disease-associated deletions demonstrated that the contact zone plays a crucial role in explaining the age-associated deletions, emphasizing its importance in the rate of healthy aging. Conclusions: Overall, we provide topological insights into the mechanism of age-associated deletion formation in human mtDNA, which could be used to predict somatic deletion burden and maximum lifespan in different human haplogroups and mammalian species. [ABSTRACT FROM AUTHOR]

Cissus quadrangularis L., a member of the Vitaceae family, is an important medicinal plant with widespread application in Indian traditional medicines. C. quadrangularis L. whole chloroplast genome of 160,404 bp was assembled using a genome skimming approach from the whole genome library. The assembled chloroplast genome contained a large single-copy region (88,987 bp), a small single-copy region (18,621 bp), and pairs of inverted repeat regions (26,398 bp). It also comprised 133 genes, including 37 tRNAs, eight rRNAs, and 88 protein-coding genes. Aside from that, we annotated three genes atpH, petB, and psbL, as well as one duplicated copy of the ycf1 gene in C. quadrangularis L. that had previously been missing from the annotation of compared Cissus chloroplast genomes. Five divergent hotspot regions such as petA_psbJ (0.1237), rps16_trnQ-UUG (0.0913), psbC_trnS-UGA (0.0847), rps15_ycf1 (0.0788), and rps2_rpoC2 (0.0788) were identified in the investigation that could aid in future species discrimination. Surprisingly, we found the overlapping genes ycf1 and ndhF on the IRb/SSC junction, rarely seen in angiosperms. The results of the phylogenetic study showed that the genomes of the Cissus species under study formed a single distinct clade. The detailed annotations given in this study could be useful in the future for genome annotations of Cissus species. The current findings of the study have the potential to serve as a useful resource for future research in the field of population genetics and the evolutionary relationships in the Cissus genus. [ABSTRACT FROM AUTHOR]

ABSTRACT Noncanonical secondary structures in nucleic acids have been studied intensively in recent years. Important biological roles of cruciform structures formed by inverted repeats (IRs) have been demonstrated in diverse organisms, including humans. Using Palindrome analyser, we analyzed IRs in all accessible bacterial genome sequences to determine their frequencies, lengths, and localizations. IR sequences were identified in all species, but their frequencies differed significantly across various evolutionary groups. We detected 242,373,717 IRs in all 1,565 bacterial genomes. The highest mean IR frequency was detected in the Tenericutes (61.89 IRs/kbp) and the lowest mean frequency was found in the Alphaproteobacteria (27.08 IRs/kbp). IRs were abundant near genes and around regulatory, tRNA, transfer-messenger RNA (tmRNA), and rRNA regions, pointing to the importance of IRs in such basic cellular processes as genome maintenance, DNA replication, and transcription. Moreover, we found that organisms with high IR frequencies were more likely to be endosymbiotic, antibiotic producing, or pathogenic. On the other hand, those with low IR frequencies were far more likely to be thermophilic. This first comprehensive analysis of IRs in all available bacterial genomes demonstrates their genomic ubiquity, nonrandom distribution, and enrichment in genomic regulatory regions. IMPORTANCE Our manuscript reports for the first time a complete analysis of inverted repeats in all fully sequenced bacterial genomes. Thanks to the availability of unique computational resources, we were able to statistically evaluate the presence and localization of these important regulatory sequences in bacterial genomes. This work revealed a strong abundance of these sequences in regulatory regions and provides researchers with a valuable tool for their manipulation.

The study of genomic structural evolution associated with accelerated evolutionary rates that result in avoidance of meltdown and increase biodiversity is becoming ever more possible as the number of available plastomes increases. To more comprehensively analyze rate heterogeneity among monocots and within Poaceae, we sequenced plastomes from four Poaceae species, combined them with publicly available data from ~200 plastomes, and conducted comparative analyses to quantify the pattern of rate heterogeneity between different lineages, functional groups, and periods of evolutionary time. We compared structural differences across the Poaceae to quantify how changes in plastome size correspond to different genomic subunits and the evolution of IR–SC junction boundaries. The substitution rates among ancestral Poaceae were inferred to be exceptionally rapid compared to other monocots but slowed after divergence into extant lineages, which could not be sufficiently explained by positive selection. As such, rapid rates in the ancestral lineage leading to Poaceae might be more closely linked to large‐scale structural changes like the loss of ycf1 and ycf2. The total increase in plastome size across Poaceae was positively correlated with the total length of intergenic spacers, tandem repeats, and dispersed repeats as well as large single copy, and inverted repeats (IRs). The continuous evolution of IR–SC junction boundaries was asynchronous with sizes of total genome and subunits across Poaceae. Future work is needed to better understand what factors in ancestral Poaceae evolved to harness such rapid rates of plastome evolution, avoid a mutational meltdown, and escape the stagnation of strong purifying selection as well as if these factors could be utilized to synthetically control rates. Highlights: Extensive sampling covering all monocotyledonous orders and subfamilies of Poaceae to compare different levels of plastome evolutionary rates heterogeneity.Comprehensively revealed the clade‐specific patterns of plastomes size variations among Poaceae subfamilies.Confirmed the distances of genes to IR–SC boundaries were effective for phylogeny, species identification, and maternal tracing of Poaceae species. [ABSTRACT FROM AUTHOR]

Non-B nucleic acids structures have arisen as key contributors to genetic variation in SARS-CoV-2. Herein, we investigated the presence of defining spike protein mutations falling within inverted repeats (IRs) for 18 SARS-CoV-2 variants, discussed the potential roles of G-quadruplexes (G4s) in SARS-CoV-2 biology, and identified potential pseudoknots within the SARS-CoV-2 genome. Surprisingly, there was a large variation in the number of defining spike protein mutations arising within IRs between variants and these were more likely to occur in the stem region of the predicted hairpin stem-loop secondary structure. Notably, mutations implicated in ACE2 binding and propagation (e.g., ΔH69/V70, N501Y, and D614G) were likely to occur within IRs, whilst mutations involved in antibody neutralization and reduced vaccine efficacy (e.g., T19R, ΔE156, ΔF157, R158G, and G446S) were rarely found within IRs. We also predicted that RNA pseudoknots could predominantly be found within, or next to, 29 mutations found in the SARS-CoV-2 spike protein. Finally, the Omicron variants BA.2, BA.4, BA.5, BA.2.12.1, and BA.2.75 appear to have lost two of the predicted G4-forming sequences found in other variants. These were found in nsp2 and the sequence complementary to the conserved stem-loop II-like motif (S2M) in the 3′ untranslated region (UTR). Taken together, non-B nucleic acids structures likely play an integral role in SARS-CoV-2 evolution and genetic diversity. [ABSTRACT FROM AUTHOR]

The current monkeypox virus (MPXV) strain differs from the strain arising in 2018 by 50+ single nucleotide polymorphisms (SNPs) and is mutating much faster than expected. The cytidine deaminase apolipoprotein B messenger RNA editing enzyme, catalytic subunit B (APOBEC3) was hypothesized to be driving this increased mutation. APOBEC has recently been identified to preferentially mutate cruciform DNA secondary structures formed by inverted repeats (IRs). IRs were recently identified as hot spots for mutation in severe acute respiratory syndrome coronavirus 2, and we aimed to identify whether IRs were also hot spots for mutation within MPXV genomes. We found that MPXV genomes were replete with IR sequences. Of the 50+ SNPs identified in the 2022 outbreak strain, 63.9% of these were found to have arisen within IR regions in the 2018 reference strain (MT903344.1). Notably, IR sequences found in the 2018 reference strain were significantly lost over time, with an average of 32.5% of these sequences being conserved in the 2022 MPXV genomes. This evidence was highly indicative that mutations were arising within IRs. This data provides further support to the hypothesis that APOBEC may be driving MPXV mutation and highlights the necessity for greater surveillance of IRs of MPXV genomes to detect new mutations. [ABSTRACT FROM AUTHOR]

Extrachromosomal circular DNA (eccDNA) is an important fraction of the genome. Recent studies proved that eccDNA plays important roles in genetic variation, aging and environmental adaptation, which have drawn wide attention. However, the characteristics of eccDNA in cattle remain unclear. Here, we studied eccDNAs from 676 cattle of 58 breeds using whole genome sequencing datasets. In total, 47 355 high‐confidence eccDNAs were identified and covered 4.6% of the cattle autosomes in length. Similarly to other species, the cattle eccDNA preferentially located in the genic and repeat sequences. Cattle eccDNAs contained complete sequences of 661 genes, which were significantly (p < 0.05) enriched in immunity‐related functions. The eccDNA was further proved to have inverted repeats on the boundaries, which contained a high proportion of A/T and ranged from 4 to 17 bp. Interestingly, we successfully separated animals according to their geographical distributions and their tissues where DNA was isolated. This implied possible roles for eccDNA in cattle selection and tissue development. Our study supplies basic knowledges on eccDNAs in cattle, which will promote understanding of extrachromosomal DNA. [ABSTRACT FROM AUTHOR]

Research pertaining to SARS-CoV-2 is in full swing to understand the origin and evolution of this deadly virus that can lead to its rapid detection. To achieve this, atypical genomic sequences which may be unique to SARS-CoV-2 or Coronaviridae family in general may be investigated. Such sequences in virus genomes may be responsible for target prediction, replication, defence mechanisms and viral packaging. This fact has motivated us to explore the different types of repeats such as palindromes, mirror repeats and inverted repeats in SARS-CoV-2, MERS-CoV and SARS-CoV-1. For this purpose, the respective reference sequence of SARS-CoV-2, MERS-CoV and SARS-CoV-1 is divided into descriptors of sequences of length ${k}$ using ${k}$ -mer technique. Thereafter, these descriptors are represented as a collection of tokens which are subsequently used for the identification of palindrome, mirror repeat and inverted repeat in the respective reference sequence. The highest number of palindromes, mirror repeats and inverted repeats are identified for descriptor length 10. As a result, for palindromes such values are 38, 42 and 33 and for mirror repeats they are 52, 38 and 33 for SARS-CoV-2, MERS-CoV and SARS-CoV-1 respectively. For inverted repeats, with a descriptor length 10 and intervening length 5, the values are 59, 56 and 70 respectively. Moreover, the identified repeats are then searched for in 108246, 291 and 340 SARS-CoV-2, MERS-CoV and SARS-CoV-1 virus sequences respectively to find the population coverage of such repeats. It surpasses 99% in most cases and even 100% for some. Furthermore, GC contents which mostly lie between 20%-50% are evaluated for these repeats as well in order to understand their binding efficacy.

Abstract Background Ziziphus hajarensis is an endemic plant species well-distributed in the Western Hajar mountains of Oman. Despite its potential medicinal uses, little is known regarding its genomic architecture, phylogenetic position, or evolution. Here we sequenced and analyzed the entire chloroplast (cp) genome of Z. hajarensis to understand its genetic organization, structure, and phylogenomic disposition among Rhamnaceae species. Results The results revealed the genome of Z. hajarensis cp comprised 162,162 bp and exhibited a typical quadripartite structure, with a large single copy (LSC) region of 895,67 bp, a small single copy (SSC) region of 19,597 bp and an inverted repeat (IR) regions of 26,499 bp. In addition, the cp genome of Z. hajarensis comprises 126 genes, including 82 protein-coding genes, eight rRNA genes, and 36 tRNA genes. Furthermore, the analysis revealed 208 microsatellites, 96.6% of which were mononucleotides. Similarly, a total of 140 repeats were identified, including 11 palindromic, 24 forward, 14 reverse, and 104 tandem repeats. The whole cp genome comparison of Z. hajarensis and nine other species from family Rhamnaceae showed an overall high degree of sequence similarity, with divergence among some intergenic spacers. Comparative phylogenetic analysis based on the complete cp genome, 66 shared genes and matK gene revealed that Z. hajarensis shares a clade with Z. jujuba and that the family Rhamnaceae is the closest family to Barbeyaceae and Elaeagnaceae. Conclusion All the genome features such as genome size, GC content, genome organization and gene order were highly conserved compared to the other related genomes. The whole cp genome of Z. hajarensis gives fascinating insights and valuable data that may be used to identify related species and reconstruct the phylogeny of the species.

Egyptian henbane (Hyoscyamus muticus L. Mant. Pl. 1767) is an important medicinal plant of the family Solanaceae. Its complete chloroplast (cp) genome was assembled using Illumina high-throughput sequencing technology to contribute to its conservation genetics studies. Here, we report the complete sequence of the chloroplast genome of H. muticus. The cp genome was 156,271 bp in length with an asymmetric base composition (30.9% A, 18.9% C, 18.6% G and 31.6% T). It encodes 132 genes comprising 87 protein-coding genes, 29 tRNA genes, and 8 rRNA genes. The overall GC content of the H. muticus chloroplast genome was 37.5%, and the corresponding values in the large single-copy (LSC), the small single-copy (SSC), and the inverted repeat (IR) regions were 35.5%, 31.5%, and 43.0%, respectively. The complete chloroplast genome sequence was deposited to the GenBank (NCBI, Accession number: MZ450974). The maximum likelihood phylogenetic analysis showed that the H. muticus and H. niger were clustered into one clade with strong support values, indicating their closer relationship.

The development of next-generation sequencing technology and the increasing amount of sequencing data have brought the bioinformatic tools used in genome assembly into focus. The final step of the process is genome annotation, which works on assembled genome sequences to identify the location of genome features. In the case of organelle genomes, specialized annotation tools are used to identify organelle genes and structural features. Numerous annotation tools target chloroplast sequences. Most chloroplast DNA genomes have a quadripartite structure caused by two copies of a large inverted repeat. We investigated the strategies of six annotation tools (Chloë, Chloroplot, GeSeq, ORG.Annotate, PGA, Plann) for identifying inverted repeats and analyzed their success using publicly available complete chloroplast sequences of taxa belonging to the asterid and rosid clades. The annotation tools use two different approaches to identify inverted repeats, using existing general search tools or implementing stand-alone solutions. The chloroplast sequences studied show that there are different types of imperfections in the assembled data and that each tool performs better on some sequences than the others. [ABSTRACT FROM AUTHOR]

Abstract Background In most flowering plants, the plastid genome exhibits a quadripartite genome structure, comprising a large and a small single copy as well as two inverted repeat regions. Thousands of plastid genomes have been sequenced and submitted to public sequence repositories in recent years. The quality of sequence annotations in many of these submissions is known to be problematic, especially regarding annotations that specify the length and location of the inverted repeats: such annotations are either missing or portray the length or location of the repeats incorrectly. However, many biological investigations employ publicly available plastid genomes at face value and implicitly assume the correctness of their sequence annotations. Results We introduce airpg, a Python package that automatically assesses the frequency of incomplete or incorrect annotations of the inverted repeats among publicly available plastid genomes. Specifically, the tool automatically retrieves plastid genomes from NCBI Nucleotide under variable search parameters, surveys them for length and location specifications of inverted repeats, and confirms any inverted repeat annotations through self-comparisons of the genome sequences. The package also includes functionality for automatic identification and removal of duplicate genome records and accounts for taxa that genuinely lack inverted repeats. A survey of the presence of inverted repeat annotations among all plastid genomes of flowering plants submitted to NCBI Nucleotide until the end of 2020 using airpg, followed by a statistical analysis of potential associations with record metadata, highlights that release year and publication status of the genome records have a significant effect on the frequency of complete and equal-length inverted repeat annotations. Conclusion The number of plastid genomes on NCBI Nucleotide has increased dramatically in recent years, and many more genomes will likely be submitted over the next decade. airpg enables researchers to automatically access and evaluate the inverted repeats of these plastid genomes as well as their sequence annotations and, thus, contributes to increasing the reliability of publicly available plastid genomes. The software is freely available via the Python package index at http://pypi.python.org/pypi/airpg .

Abstract Background Angelica L. (family Apiaceae) is an economically important genus comprising ca. One hundred ten species. Angelica species are found on all continents of the Northern Hemisphere, and East Asia hosts the highest number of species. Morphological characters such as fruit anatomy, leaf morphology and subterranean structures of Angelica species show extreme diversity. Consequently, the taxonomic classification of Angelica species is complex and remains controversial, as the classifications proposed by previous studies based on morphological data and molecular data are highly discordant. In addition, the phylogenetic relationships of major clades in the Angelica group, particularly in the Angelica s. s. clade, remain unclear. Chloroplast (cp) genome sequences have been widely used in phylogenetic studies and for evaluating genetic diversity. Results In this study, we sequenced and assembled 28 complete cp genomes from 22 species, two varieties and two cultivars of Angelica. Combined with 36 available cp genomes in GenBank from representative clades of the subfamily Apioideae, the characteristics and evolutionary patterns of Angelica cp genomes were studied, and the phylogenetic relationships of Angelica species were resolved. The Angelica cp genomes had the typical quadripartite structure including a pair of inverted repeats (IRs: 5836–34,706 bp) separated by a large single-copy region (LSC: 76,657–103,161 bp) and a small single-copy region (SSC: 17,433–21,794 bp). Extensive expansion and contraction of the IR region were observed among cp genomes of Angelica species, and the pattern of the diversification of cp genomes showed high consistency with the phylogenetic placement of Angelica species. Species of Angelica were grouped into two major clades, with most species grouped in the Angelica group and A. omeiensis and A. sinensis grouped in the Sinodielsia with Ligusticum tenuissimum. Conclusions Our results further demonstrate the power of plastid phylogenomics in enhancing the phylogenetic reconstructions of complex genera and provide new insights into plastome evolution across Angelica L.

The plastid organelle is essential for many vital cellular processes and the growth and development of plants. The availability of a large number of complete plastid genomes could be effectively utilized to understand the evolution of the plastid genomes and phylogenetic relationships among plants. We comprehensively analyzed the plastid genomes of Viridiplantae comprising 3,654 taxa from 298 families and 111 orders and compared the genomic organizations in their plastid genomic DNA among major clades, which include gene gain/loss, gene copy number, GC content, and gene blocks. We discovered that some important genes that exhibit similar functions likely formed gene blocks, such as the psb family presumably showing co-occurrence and forming gene blocks in Viridiplantae. The inverted repeats (IRs) in plastid genomes have doubled in size across land plants, and their GC content is substantially higher than non-IR genes. By employing three different data sets [all nucleotide positions (nt123), only the first and second codon positions (nt12), and amino acids (AA)], our phylogenomic analyses revealed Chlorokybales + Mesostigmatales as the earliest-branching lineage of streptophytes. Hornworts, mosses, and liverworts forming a monophylum were identified as the sister lineage of tracheophytes. Based on nt12 and AA data sets, monocots, Chloranthales and magnoliids are successive sister lineages to the eudicots + Ceratophyllales clade. The comprehensive taxon sampling and analysis of different data sets from plastid genomes recovered well-supported relationships of green plants, thereby contributing to resolving some long-standing uncertainties in the plant phylogeny. [ABSTRACT FROM AUTHOR]

Members of the Peronosporaceae (Oomycota, Chromista), which currently consists of 25 genera and approximately 1,000 recognized species, are responsible for disease on a wide range of plant hosts. Molecular phylogenetic analyses over the last two decades have improved our understanding of evolutionary relationships within Peronosporaceae. To date, 16 numbered and three named clades have been recognized; it is clear from these studies that the current taxonomy does not reflect evolutionary relationships. Whole organelle genome sequences are an increasingly important source of phylogenetic information, and in this study, we present comparative and phylogenetic analyses of mitogenome sequences from 15 of the 19 currently recognized clades of Peronosporaceae, including 44 newly assembled sequences. Our analyses suggest strong conservation of mitogenome size and gene content across Peronosporaceae but, as previous studies have suggested, limited conservation of synteny. Specifically, we identified 28 distinct syntenies amongst the 71 examined isolates. Moreover, 19 of the isolates contained inverted or direct repeats, suggesting repeated sequences may be more common than previously thought. In terms of phylogenetic relationships, our analyses of 34 concatenated mitochondrial gene sequences resulted in a topology that was broadly consistent with previous studies. However, unlike previous studies concatenated mitochondrial sequences provided strong support for higher-level relationships within the family. [ABSTRACT FROM AUTHOR]

Abstract Background Co-expression of two distinct guide RNAs (gRNAs) has been used to facilitate the application of CRISPR/Cas9 system in fields such as large genomic deletion. The paired gRNAs are often placed adjacently in the same direction and expressed individually by two identical promoters, constituting direct repeats (DRs) which are susceptible to self-homologous recombination. As a result, the paired-gRNA plasmids cannot remain stable, which greatly prevents extensible applications of CRISPR/Cas9 system. Results To address this limitation, different DRs-involved paired-gRNA plasmids were designed and the events of recombination were characterized. Deletion between DRs occurred with high frequencies during plasmid construction and subsequent plasmid propagation. This recombination event was RecA-independent, which agreed with the replication slippage model. To increase plasmid stability, a reversed paired-gRNA plasmids (RPGPs) cloning strategy was developed by converting DRs to the more stable invert repeats (IRs), which completely eliminated DRs-induced recombination. Using RPGPs, rapid deletion of chromosome fragments up to 100 kb with an efficiency of 83.33% was achieved in Escherichia coli. Conclusions The RPGPs cloning strategy serves as a general solution to avoid plasmid RecA-independent recombination. It can be adapted to applications that rely on paired gRNAs or repeated genetic parts.

The plastid organelle is essential for many vital cellular processes and the growth and development of plants. The availability of a large number of complete plastid genomes could be effectively utilized to understand the evolution of the plastid genomes and phylogenetic relationships among plants. We comprehensively analyzed the plastid genomes of Viridiplantae comprising 3,654 taxa from 298 families and 111 orders and compared the genomic organizations in their plastid genomic DNA among major clades, which include gene gain/loss, gene copy number, GC content, and gene blocks. We discovered that some important genes that exhibit similar functions likely formed gene blocks, such as the psb family presumably showing co-occurrence and forming gene blocks in Viridiplantae. The inverted repeats (IRs) in plastid genomes have doubled in size across land plants, and their GC content is substantially higher than non-IR genes. By employing three different data sets [all nucleotide positions (nt123), only the first and second codon positions (nt12), and amino acids (AA)], our phylogenomic analyses revealed Chlorokybales + Mesostigmatales as the earliest-branching lineage of streptophytes. Hornworts, mosses, and liverworts forming a monophylum were identified as the sister lineage of tracheophytes. Based on nt12 and AA data sets, monocots, Chloranthales and magnoliids are successive sister lineages to the eudicots + Ceratophyllales clade. The comprehensive taxon sampling and analysis of different data sets from plastid genomes recovered well-supported relationships of green plants, thereby contributing to resolving some long-standing uncertainties in the plant phylogeny.

Clostridioides difficile RT027 strains cause infections that vary in severity from asymptomatic to lethal, but the molecular basis for this variability is poorly understood. Through comparative analyses of RT027 clinical isolates, we determined that isolates that exhibit greater variability in their flagellar gene expression exhibit greater virulence in vivo . C. difficile flagellar genes are phase-variably expressed due to the site-specific inversion of the flgB 5'UTR region, which reversibly generates ON vs. OFF orientations for the flagellar switch. We found that longer inverted repeat (IR) sequences in this switch region correlate with greater disease severity, with RT027 strains carrying 6A/6T IR sequences exhibiting greater phenotypic heterogeneity in flagellar gene expression (60%-75% ON) and causing more severe disease than those with shorter IRs (> 99% ON or OFF). Taken together, our results reveal that phenotypic heterogeneity in flagellar gene expression may contribute to the variable disease severity observed in C. difficile patients., Competing Interests: Declaration of interests None.

Phase variation is defined as the rapid and reversible switching of gene expression, and typically occurs in genes encoding surface features in small genome bacterial pathogens. Phase variation has evolved to provide an extra survival mechanism in bacteria that lack multiple 'sense-and-respond' gene regulation systems. Many bacterial pathogens also encode DNA methyltransferases that are phase-variable, controlling systems called 'phasevarions' (phase-variable regulons). This primer will summarize the current understanding of phase variation, describing the role of major phase-variable factors, and phasevarions, in bacterial pathobiology.

Background: Ziziphus hajarensis is an endemic plant species well-distributed in the Western Hajar mountains of Oman. Despite its potential medicinal uses, little is known regarding its genomic architecture, phylogenetic position, or evolution. Here we sequenced and analyzed the entire chloroplast (cp) genome of Z. hajarensis to understand its genetic organization, structure, and phylogenomic disposition among Rhamnaceae species. Results: The results revealed the genome of Z. hajarensis cp comprised 162,162 bp and exhibited a typical quadripartite structure, with a large single copy (LSC) region of 895,67 bp, a small single copy (SSC) region of 19,597 bp and an inverted repeat (IR) regions of 26,499 bp. In addition, the cp genome of Z. hajarensis comprises 126 genes, including 82 protein-coding genes, eight rRNA genes, and 36 tRNA genes. Furthermore, the analysis revealed 208 microsatellites, 96.6% of which were mononucleotides. Similarly, a total of 140 repeats were identified, including 11 palindromic, 24 forward, 14 reverse, and 104 tandem repeats. The whole cp genome comparison of Z. hajarensis and nine other species from family Rhamnaceae showed an overall high degree of sequence similarity, with divergence among some intergenic spacers. Comparative phylogenetic analysis based on the complete cp genome, 66 shared genes and matK gene revealed that Z. hajarensis shares a clade with Z. jujuba and that the family Rhamnaceae is the closest family to Barbeyaceae and Elaeagnaceae. Conclusion: All the genome features such as genome size, GC content, genome organization and gene order were highly conserved compared to the other related genomes. The whole cp genome of Z. hajarensis gives fascinating insights and valuable data that may be used to identify related species and reconstruct the phylogeny of the species. [ABSTRACT FROM AUTHOR]

Cephaleuros is often known as an algal pathogen with 19 taxonomically valid species, some of which are responsible for red rust and algal spot diseases in vascular plants. No chloroplast genomes have yet been reported in this genus, and the limited genetic information is an obstacle to understanding the evolution of this genus. In this study, we sequenced six new Trentepohliales chloroplast genomes, including four Cephaleuros and two Trentepohlia. The chloroplast genomes of Trentepohliales are large compared to most green algae, ranging from 216 to 408 kbp. They encode between 93 and 98 genes and have a GC content of 26–36%. All new chloroplast genomes were circular-mapping and lacked a quadripartite structure, in contrast to the previously sequenced Trentepohlia odorata , which does have an inverted repeat. The duplicated trnD -GTC , petD , and atpA genes in C. karstenii may be remnants of the IR region and shed light on its reduction. Chloroplast genes of Trentepohliales show elevated rates of evolution, strong rearrangement dynamics and several genes display an alternative genetic code with reassignment of the UGA/UAG codon presumably coding for arginine. Our results present the first whole chloroplast genome of the genus Cephaleuros and enrich the chloroplast genome resources of Trentepohliales. [ABSTRACT FROM AUTHOR]

The Oedogoniales comprises the genera Oedogonium, Oedocladium and Bulbochaete which include more than 600 described species. However, only two chloroplast genome sequences have been reported to date, and little is known about the evolutionary relationships between these genera or their infrageneric classification. This group is economically important in astaxanthin production and the energy sector. We examined four Oedocladium strains from different locations in China. One was identified as O. prescottii, and the morphology and molecular data for the other three strains suggested that they were O. carolinianum. Our chloroplast (cp) genome sequences of O. prescottii (FACHB-2452) and O. carolinianum (FACHB-2453) were compared with previously reported Oedogoniales cp genomes. The cp genome of O. carolinianum FACHB-2453 was of similar size, structure and composition to that of O. carolinianum UTEX LB 1686, apart from two missing group I introns in the genes psbA and psbD. O. prescottii showed more differences with regard to genome size, structure, number of genes and introns, and repeats than the other Oedogoniales cp genomes. The inverted repeats (IRs) region of the O. prescottii cp genome had a 10-kb contraction. Synteny analysis of the four cp genomes showed a high degree of syntenic conservation in general; however, some structural variation at genus level was observed in Oedocladium. In phylogenetic analyses of Oedocladium using sequences of 18S rDNA, ITS and 59 cp protein-coding genes, species were placed on two separate branches within the Oedogonium clade, indicating that Oedocladium is a paraphyletic group. [ABSTRACT FROM AUTHOR]

Cephaleuros is often known as an algal pathogen with 19 taxonomically valid species, some of which are responsible for red rust and algal spot diseases in vascular plants. No chloroplast genomes have yet been reported in this genus, and the limited genetic information is an obstacle to understanding the evolution of this genus. In this study, we sequenced six new Trentepohliales chloroplast genomes, including four Cephaleuros and two Trentepohlia. The chloroplast genomes of Trentepohliales are large compared to most green algae, ranging from 216 to 408 kbp. They encode between 93 and 98 genes and have a GC content of 26–36%. All new chloroplast genomes were circular-mapping and lacked a quadripartite structure, in contrast to the previously sequenced Trentepohlia odorata, which does have an inverted repeat. The duplicated trnD-GTC, petD, and atpA genes in C. karstenii may be remnants of the IR region and shed light on its reduction. Chloroplast genes of Trentepohliales show elevated rates of evolution, strong rearrangement dynamics and several genes display an alternative genetic code with reassignment of the UGA/UAG codon presumably coding for arginine. Our results present the first whole chloroplast genome of the genus Cephaleuros and enrich the chloroplast genome resources of Trentepohliales.

Background: In most flowering plants, the plastid genome exhibits a quadripartite genome structure, comprising a large and a small single copy as well as two inverted repeat regions. Thousands of plastid genomes have been sequenced and submitted to public sequence repositories in recent years. The quality of sequence annotations in many of these submissions is known to be problematic, especially regarding annotations that specify the length and location of the inverted repeats: such annotations are either missing or portray the length or location of the repeats incorrectly. However, many biological investigations employ publicly available plastid genomes at face value and implicitly assume the correctness of their sequence annotations. Results: We introduce airpg, a Python package that automatically assesses the frequency of incomplete or incorrect annotations of the inverted repeats among publicly available plastid genomes. Specifically, the tool automatically retrieves plastid genomes from NCBI Nucleotide under variable search parameters, surveys them for length and location specifications of inverted repeats, and confirms any inverted repeat annotations through self-comparisons of the genome sequences. The package also includes functionality for automatic identification and removal of duplicate genome records and accounts for taxa that genuinely lack inverted repeats. A survey of the presence of inverted repeat annotations among all plastid genomes of flowering plants submitted to NCBI Nucleotide until the end of 2020 using airpg, followed by a statistical analysis of potential associations with record metadata, highlights that release year and publication status of the genome records have a significant effect on the frequency of complete and equal-length inverted repeat annotations. Conclusion: The number of plastid genomes on NCBI Nucleotide has increased dramatically in recent years, and many more genomes will likely be submitted over the next decade. airpg enables researchers to automatically access and evaluate the inverted repeats of these plastid genomes as well as their sequence annotations and, thus, contributes to increasing the reliability of publicly available plastid genomes. The software is freely available via the Python package index at http://pypi.python.org/pypi/airpg. [ABSTRACT FROM AUTHOR]

Egyptian henbane (Hyoscyamus muticus L. Mant. Pl. 1767) is an important medicinal plant of the family Solanaceae. Its complete chloroplast (cp) genome was assembled using Illumina high-throughput sequencing technology to contribute to its conservation genetics studies. Here, we report the complete sequence of the chloroplast genome of H. muticus. The cp genome was 156,271 bp in length with an asymmetric base composition (30.9% A, 18.9% C, 18.6% G and 31.6% T). It encodes 132 genes comprising 87 protein-coding genes, 29 tRNA genes, and 8 rRNA genes. The overall GC content of the H. muticus chloroplast genome was 37.5%, and the corresponding values in the large single-copy (LSC), the small single-copy (SSC), and the inverted repeat (IR) regions were 35.5%, 31.5%, and 43.0%, respectively. The complete chloroplast genome sequence was deposited to the GenBank (NCBI, Accession number: MZ450974). The maximum likelihood phylogenetic analysis showed that the H. muticus and H. niger were clustered into one clade with strong support values, indicating their closer relationship. [ABSTRACT FROM AUTHOR]

SARS-CoV-2 is an intensively investigated virus from the order Nidovirales (Coronaviridae family) that causes COVID-19 disease in humans. Through enormous scientific effort, thousands of viral strains have been sequenced to date, thereby creating a strong background for deep bioinformatics studies of the SARS-CoV-2 genome. In this study, we inspected high-frequency mutations of SARS-CoV-2 and carried out systematic analyses of their overlay with inverted repeat (IR) loci and CpG islands. The main conclusion of our study is that SARS-CoV-2 hot-spot mutations are significantly enriched within both IRs and CpG island loci. This points to their role in genomic instability and may predict further mutational drive of the SARS-CoV-2 genome. Moreover, CpG islands are strongly enriched upstream from viral ORFs and thus could play important roles in transcription and the viral life cycle. We hypothesize that hypermethylation of these loci will decrease the transcription of viral ORFs and could therefore limit the progression of the disease. [ABSTRACT FROM AUTHOR]

Non-canonical nucleic acid structures play important roles in the regulation of molecular processes. Considering the importance of the ongoing coronavirus crisis, we decided to evaluate genomes of all coronaviruses sequenced to date (stated more broadly, the order Nidovirales) to determine if they contain non-canonical nucleic acid structures. We discovered much evidence of putative G-quadruplex sites and even much more of inverted repeats (IRs) loci, which in fact are ubiquitous along the whole genomic sequence and indicate a possible mechanism for genomic RNA packaging. The most notable enrichment of IRs was found inside 5′UTR for IRs of size 12+ nucleotides, and the most notable enrichment of putative quadruplex sites (PQSs) was located before 3′UTR, inside 5′UTR, and before mRNA. This indicates crucial regulatory roles for both IRs and PQSs. Moreover, we found multiple G-quadruplex binding motifs in human proteins having potential for binding of SARS-CoV-2 RNA. Non-canonical nucleic acids structures in Nidovirales and in novel SARS-CoV-2 are therefore promising druggable structures that can be targeted and utilized in the future.

BetP is an Na(+)-coupled betaine-specific transporter of the betaine-choline-carnitine (BCC) transporter family involved in the response to hyperosmotic stress. The crystal structure of BetP revealed an overall fold of two inverted structurally related repeats (LeuT-fold) that BetP shares with other sequence-unrelated Na(+)-coupled symporters. Numerous structures of LeuT-fold transporters in distinct conformational states have contributed substantially to our understanding of the alternating access mechanism of transport. Nevertheless, coupling of substrate and co-transported ion fluxes has not been structurally corroborated to the same extent. We converted BetP by a single-point mutation--glycine to aspartate--into an H(+)-coupled choline-specific transporter and solved the crystal structure of this mutant in complex with choline. The structure of BetP-G153D demonstrates a new inward-facing open conformation for BetP. Choline binding to a location close to the second, low-affinity sodium-binding site (Na2) of LeuT-fold transporters is facilitated by the introduced aspartate. Our data confirm the importance of a cation-binding site in BetP, playing a key role in a proposed molecular mechanism of Na(+) and H(+) coupling in BCC transporters.

Site-specific recombination is a DNA breaking and reconstructing process that plays important roles in various cellular pathways for both prokaryotes and eukaryotes. This process requires a site-specific recombinase and direct or inverted repeats. Some tyrosine site-specific recombinases catalyze DNA inversions and regulate subpopulation diversity and phase variation in many bacterial species. In Streptococcus pneumoniae, the PsrA tyrosine recombinase was shown to control DNA inversions in the three DNA methyltransferase hsdS genes of the type I restriction-modification cod locus. Such DNA inversions are mediated by three inverted repeats (IR1, IR2, and IR3). In this work, we purified an untagged form of the PsrA protein and studied its DNA-binding and catalytic features. Gel retardation assays showed that PsrA binds to linear and supercoiled DNAs, containing or not inverted repeats. Nevertheless, DNase I footprinting assays showed that, on linear DNAs, PsrA has a preference for sites that include an IR1 sequence (IR1.1 or IR1.2) and its boundary sequences. Furthermore, on supercoiled DNAs, PsrA was able to generate DNA inversions between specific inverted repeats (IR1, IR2, and IR3), which supports its ability to locate specific target sites. Unlike other site-specific recombinases, PsrA showed reliance on magnesium ions for efficient catalysis of IR1-mediated DNA inversions. We discuss that PsrA might find its specific binding sites on the bacterial genome by a mechanism that involves transitory non-specific interactions between protein and DNA.

Based on the nrDNA ITS sequence data, the Tordylieae tribe is recognized as monophyletic with three major lineages: the subtribe Tordyliinae, the Cymbocarpum clade, and the Lefebvrea clade. Recent phylogenomic investigations showed incongruence between the nuclear and plastid genome evolution in the tribe. To assess phylogenetic relations and structure evolution of plastomes in Tordylieae, we generated eleven complete plastome sequences using the genome skimming approach and compared them with the available data from this tribe and close relatives. Newly assembled plastomes had lengths ranging from 141,148 to 150,103 base pairs and contained 122–127 genes, including 79–82 protein-coding genes, 35–37 tRNAs, and 8 rRNAs. We observed substantial differences in the inverted repeat length and gene content, accompanied by a complex picture of multiple JLA and JLB shifts. In concatenated phylogenetic analyses, Tordylieae plastomes formed at least three not closely related lineages with plastomes of the Lefebvrea clade as a sister group to plastomes from the Selineae tribe. The newly obtained data have increased our knowledge on the range of plastome variability in Apiaceae.

Background: Angelica L. (family Apiaceae) is an economically important genus comprising ca. One hundred ten species. Angelica species are found on all continents of the Northern Hemisphere, and East Asia hosts the highest number of species. Morphological characters such as fruit anatomy, leaf morphology and subterranean structures of Angelica species show extreme diversity. Consequently, the taxonomic classification of Angelica species is complex and remains controversial, as the classifications proposed by previous studies based on morphological data and molecular data are highly discordant. In addition, the phylogenetic relationships of major clades in the Angelica group, particularly in the Angelica s. s. clade, remain unclear. Chloroplast (cp) genome sequences have been widely used in phylogenetic studies and for evaluating genetic diversity. Results: In this study, we sequenced and assembled 28 complete cp genomes from 22 species, two varieties and two cultivars of Angelica. Combined with 36 available cp genomes in GenBank from representative clades of the subfamily Apioideae, the characteristics and evolutionary patterns of Angelica cp genomes were studied, and the phylogenetic relationships of Angelica species were resolved. The Angelica cp genomes had the typical quadripartite structure including a pair of inverted repeats (IRs: 5836–34,706 bp) separated by a large single-copy region (LSC: 76,657–103,161 bp) and a small single-copy region (SSC: 17,433–21,794 bp). Extensive expansion and contraction of the IR region were observed among cp genomes of Angelica species, and the pattern of the diversification of cp genomes showed high consistency with the phylogenetic placement of Angelica species. Species of Angelica were grouped into two major clades, with most species grouped in the Angelica group and A. omeiensis and A. sinensis grouped in the Sinodielsia with Ligusticum tenuissimum. Conclusions: Our results further demonstrate the power of plastid phylogenomics in enhancing the phylogenetic reconstructions of complex genera and provide new insights into plastome evolution across Angelica L. [ABSTRACT FROM AUTHOR]

In genome engineering, sgRNAs define the genomic target to be modified in CRISPR/Cas9 system. Either for single gene editing or genome-wide screens, sgRNAs are cloned into plasmid vectors. During our performance of CRISPR/Cas9 gene knock out, we found that the central part of a sgRNA was inverted after transformation into Escherichia coli. Interestingly, the inverted portion was found to be flanked by inverted repeats, and sealing of nicks inside the plasmid could correct the inversion. This type of sgRNA recombination completely changed its original sequence and should be noted during sgRNA design and performance of CRISPR/Cas9. [ABSTRACT FROM AUTHOR]

Non-canonical nucleic acid structures play important roles in the regulation of molecular processes. Considering the importance of the ongoing coronavirus crisis, we decided to evaluate genomes of all coronaviruses sequenced to date (stated more broadly, the order Nidovirales) to determine if they contain non-canonical nucleic acid structures. We discovered much evidence of putative G-quadruplex sites and even much more of inverted repeats (IRs) loci, which in fact are ubiquitous along the whole genomic sequence and indicate a possible mechanism for genomic RNA packaging. The most notable enrichment of IRs was found inside 5′UTR for IRs of size 12+ nucleotides, and the most notable enrichment of putative quadruplex sites (PQSs) was located before 3′UTR, inside 5′UTR, and before mRNA. This indicates crucial regulatory roles for both IRs and PQSs. Moreover, we found multiple G-quadruplex binding motifs in human proteins having potential for binding of SARS-CoV-2 RNA. Non-canonical nucleic acids structures in Nidovirales and in novel SARS-CoV-2 are therefore promising druggable structures that can be targeted and utilized in the future. [ABSTRACT FROM AUTHOR]

The advent of novel high-throughput sequencing techniques has revealed that eukaryotic genomes are massively transcribed although only a small fraction of RNAs exhibits protein-coding capacity. In the last years, long noncoding RNAs (lncRNAs) have emerged as regulators of eukaryotic gene expression in a wide range of molecular mechanisms. Plant lncRNAs can be transcribed by alternative RNA polymerases, acting directly as long transcripts or can be processed into active small RNAs. Several lncRNAs have been recently shown to interact with chromatin, DNA or nuclear proteins to condition the epigenetic environment of target genes or modulate the activity of transcriptional complexes. In this review, we will summarize the recent discoveries about the actions of plant lncRNAs in the regulation of gene expression at the transcriptional level. [ABSTRACT FROM AUTHOR]

Background: Co-expression of two distinct guide RNAs (gRNAs) has been used to facilitate the application of CRISPR/Cas9 system in fields such as large genomic deletion. The paired gRNAs are often placed adjacently in the same direction and expressed individually by two identical promoters, constituting direct repeats (DRs) which are susceptible to self-homologous recombination. As a result, the paired-gRNA plasmids cannot remain stable, which greatly prevents extensible applications of CRISPR/Cas9 system. Results: To address this limitation, different DRs-involved paired-gRNA plasmids were designed and the events of recombination were characterized. Deletion between DRs occurred with high frequencies during plasmid construction and subsequent plasmid propagation. This recombination event was RecA-independent, which agreed with the replication slippage model. To increase plasmid stability, a reversed paired-gRNA plasmids (RPGPs) cloning strategy was developed by converting DRs to the more stable invert repeats (IRs), which completely eliminated DRs-induced recombination. Using RPGPs, rapid deletion of chromosome fragments up to 100 kb with an efficiency of 83.33% was achieved in Escherichia coli. Conclusions: The RPGPs cloning strategy serves as a general solution to avoid plasmid RecA-independent recombination. It can be adapted to applications that rely on paired gRNAs or repeated genetic parts. [ABSTRACT FROM AUTHOR]

Okazaki fragments that are formed during lagging strand DNA synthesis include an initiating primer consisting of both RNA and DNA. The RNA fragment must be removed before the fragments are joined. In Saccharomyces cerevisiae, a key player in this process is the structure-specific flap endonuclease, Rad27p (human homolog FEN1). To obtain a genomic view of the mutational consequence of loss of RAD27, a S. cerevisiae rad27Δ strain was subcultured for 25 generations and sequenced using Illumina paired-end sequencing. Out of the 455 changes observed in 10 colonies isolated the two most common types of events were insertions or deletions (INDELs) in simple sequence repeats (SSRs) and INDELs mediated by short direct repeats. Surprisingly, we also detected a previously neglected class of 21 template-switching events. These events were presumably generated by quasi-palindrome to palindrome correction, as well as palindrome elongation. The formation of these events is best explained by folding back of the stalled nascent strand and resumption of DNA synthesis using the same nascent strand as a template. Evidence of quasi-palindrome to palindrome correction that could be generated by template switching appears also in yeast genome evolution. Out of the 455 events, 55 events appeared in multiple isolates; further analysis indicates that these loci are mutational hotspots. Since Rad27 acts on the lagging strand when the leading strand should not contain any gaps, we propose a mechanism favoring intramolecular strand switching over an intermolecular mechanism. We note that our results open new ways of understanding template switching that occurs during genome instability and evolution.

Long inverted repeats (LIRs) are evolutionarily and functionally important structures in genomes because of their involvement in RNA interference, DNA recombination, and gene duplication. Identification of LIRs is highly complicated when mismatches and indels between the repeats are permitted. Long inverted repeat explorer (Lirex) was developed and introduced in this report. Written in Java, Lirex provides a user-friendly interface and allows users to specify LIR searching criteria, such as length of the region, as well as pattern and size of the repeats. Recombinogenic LIRs can be selected on the basis of mismatch rate and internal spacer size from identified LIRs. Lirex, as a cross-platform tool to identify LIRs in a genome, may assist in designing following experiments to explore the function of LIRs. Our tool can identify more LIRs than other LIR searching tools. Lirex is publicly available at http://124.16.219.129/Lirex.

Rheum likiangense Samuelsson (Polygonaceae) is an endangered alpine rhubarb in the Qinghai-Tibet Plateau. In this study we report the complete chloroplast genome sequence (plastome) of Rh. likiangense. The assembled plastome is 162,291 bp in length with 31,741 bp inverted repeat (IR) regions and 128 annotated genes, including 34 tRNA genes, 8 rRNA genes, and 86 protein-coding genes. Phylogenetic analyses based on the full plastome sequences suggest the close relationship of Rh. likiangense with Rh. acuminatum and Rh. nobile. The plastome reported here is highly useful for designing plastome SSR markers to conduct a further conservation genetic study of this endangered rhubarb.

We sequenced the complete mitogenome of a specimen of Cymbiola nobilis (Neogastropoda, Volutidae) from Phú Quốc Island in Vietnam. The mitogenome of this predatory sea snail is 16,314-bp long and comprises the typical set of genes encoding 13 proteins, 2 rRNAs, and 22 tRNAs. It features a complex control region that contains two copies of an inverted repeat sequence separated by an AT-rich region, a characteristic previously observed in a few other gastropods. Phylogenetic analysis based on all protein-coding genes revealed that the Volutidae form a strongly supported clade.