Inverted Repeat Sequences Research Articles

The Complete Chloroplast Genome of Meconopsis simplicifolia and Its Genetic Comparison to Other Meconopsis Species.

Background: Chloroplasts, due to their high conservation and lack of recombination, serve as important genetic resources for the classification and evolutionary analysis of closely related species that are difficult to distinguish based on their morphological features. Meconopsis simplicifolia (M. simplicifolia), an endangered herb within the Meconopsis genus, has demonstrated therapeutic potential in treating various diseases. However, the highly polymorphic morphology of this species poses a challenge for accurate identification. Methods: In this study, the complete chloroplast genome of M. simplicifolia was sequenced and assembled using Illumina sequencing technology. Simple sequence repeats (SSRs) and repetitive sequences were characterized. In addition, a comparative analysis was conducted with the chloroplast genomes of six other Meconopsis species. Results: The chloroplast genome of M. simplicifolia has a quadripartite circular structure with a total length of 152,772 bp. It consists of a large single-copy region of 83,824 bp and a small single-copy region of 17,646 bp, separated by a pair of inverted repeat sequences (IRa and IRb, 25,651 bp). The genome contains 131 genes, 33 SSRs, and 27 long repetitive sequences. Comparative analysis with six other chloroplast genomes of Meconopsis revealed that M. simplicifolia is closely related to M. betonicifolia and that the rpl2 (ribosomal protein L2) gene in the IRb region has been deleted. This deletion is of significant importance for future taxonomic studies of M. simplicifolia. Conclusions: This study provides a valuable reference for the identification of M. simplicifolia and contributes to a deeper understanding of the phylogeny and evolution of the Meconopsis genus.

Flagellar switch inverted repeat impacts flagellar invertibility and varies Clostridioides difficile RT027/MLST1 virulence.

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.

The complete chloroplast genome of Lindernia crustacea (L.) F. Muell 1882 (Linderniaceae) and its phylogenetic analysis

Lindernia crustacea (L.) F. Muell 1882, a species in the Linderniaceae family, holds traditional medicinal value in China. To investigate its genetic diversity, we assembled, annotated, and characterized the first complete chloroplast genome of L. crustacea using Illumina sequencing data and various bioinformatics tools. The genome is 153,647 bp in length, with a GC content of 37.6%. It exhibits a typical quadripartite structure, consisting of a large single-copy region (LSC) of 85,411 bp, a small single-copy region (SSC) of 18,724 bp, and two inverted repeat sequences (IRa and IRb) of 25,816 bp each. The genome was predicted to contain 131 genes, including 87 protein-coding genes, 36 tRNA genes, and eight rRNA genes. Phylogenomic analysis indicated that L. crustacea is closely related to L. stricta. These findings provide a foundation for further research on the evolution and potential medicinal applications of the Linderniaceae family.

Identification of new banana endogenous virus sequences highlights the hallmark gene encoded by retroviruses integrated in banana genomes

Endogenous pararetrovirus sequences (EPRVs) originated from DNA viruses of the family Caulimoviridae are widely present in plant genomes. Banana streak viruses (BSVs) are a group of circular double-stranded DNA viruses in the genus Badnavirus of the family Caulimoviridae. Banana endogenous virus sequences (BEVs) derived from the ancestral genes of badnaviruses and fixed in the genomes of various bananas. However, the genomic characteristics of BEVs remain unknown. In this study, we identified 2 new variants of BEVs GZ5 and GZ13 by sequences analyses, Southern blot, and fluorescent in situ hybridization (FISH). BEV GZ5 had one copy of integration in the BB genome of bananas, while BEV GZ13 was only present in the genome of the variety Dajiao. Importantly, BEV GZ5 contained a complete gene of reverse transcriptase (RT) and ribonuclease H (RNase H) (RT/RNase H). In addition, a 340-bp inverted repeat sequence partially overlapping with RNase H was found upstream and downstream of BEV GZ5. However, the amino acid sequences of BEV GZ5 had deletions and mutations compared with BSVs. The bioinformatics analyses showed that BEV GZ5 protein composed of 412 amino acids with a molecular weight of 47.37 kDa and an isoelectric point of 9.40. Leucine, isoleucine, and lysine (Lys) were the main amino acids of BEV GZ5 protein. The analyses revealed that BEV GZ5 protein contained 35 potential phosphorylation sites. Additionally, it was a hydrophilic protein without a signal peptide and transmembrane region. The secondary structure of BEV GZ5 protein consisted of 37.26% α-helix, followed by 36.25% random coil. To our knowledge, this is the first report that novel BEVs with the complete gene of RT/RNase H has been characterized, which provide a basis for further exploration the function and integration mechanism of BEVs in bananas.

The Ty1 retrotransposon harbors a DNA region that performs dual functions as both a gene silencing and chromatin insulator

In various eukaryotic kingdoms, long terminal repeat (LTR) retrotransposons repress transcription by infiltrating heterochromatin generated within their elements. In contrast, the budding yeast LTR retrotransposon Ty1 does not itself undergo transcriptional repression, although it is capable of repressing the transcription of the inserted genes within it. In this study, we identified a DNA region within Ty1 that exerts its silencing effect via sequence orientation. We identified a DNA region within the Ty1 group-specific antigen (GAG) gene that causes gene silencing, termed GAG silencing (GAGsi), in which the silent chromatin in the GAGsi region is created by euchromatin-specific histone modifications. A characteristic inverted repeat (IR) sequence is present at the 5' end of this region, forming a chromatin boundary between promoter-specific chromatin upstream of the IR sequence and silent chromatin downstream of the IR sequence. In addition, Esc2 and Rad57, which are involved in DNA repair, were required for GAGsi silencing. Finally, the chromatin boundary was required for the transcription of Ty1 itself. Thus, the GAGsi sequence contributes to the creation of a chromatin environment that promotes Ty1 transcription.

A self-complementary AAV proviral plasmid that reduces cross-packaging and ITR promoter activity in AAV vector preparations

Adeno-associated viral vectors (AAV) are a leading delivery system for gene therapy in animal models and humans. With several FDA-approved AAV gene therapies on the market, issues related to vector manufacturing have become increasingly important. In this study, we focused on potentially toxic DNA contaminants that can arise from AAV proviral plasmids, the raw materials required for manufacturing recombinant AAV in eukaryotic cells. Typical AAV proviral plasmids are circular DNAs containing a therapeutic gene cassette flanked by natural AAV inverted terminal repeat (ITR) sequences, and a plasmid backbone carrying prokaryotic sequences required for plasmid replication and selection in bacteria. While the majority of AAV particles package the intended therapeutic payload, some capsids instead package the bacterial sequences located on the proviral plasmid backbone. Since ITR sequences also have promoter activity, potentially toxic bacterial open reading frames can be produced in vivo, thereby representing a safety risk. In this study, we describe a new AAV proviral plasmid for vector manufacturing that (1) significantly decreases cross-packaged bacterial sequences; (2) increases correctly packaged AAV payloads; and (3) blunts ITR-driven transcription of cross-packaged material to avoid expressing potentially toxic bacterial sequences. This system may help improve the safety of AAV vector products.

Dynamic stem–loop extension by Pol θ and templated insertion during DNA repair

Theta-mediated end-joining (TMEJ) is critical for survival of cancer cells when other DNA double-stranded break repair pathways are impaired. Human DNA polymerase theta (Pol θ) can extend single-stranded DNA oligonucleotides, but little is known about preferred substrates and mechanism. We show that Pol θ can extend both single-stranded DNA and RNA substrates by unimolecular stem loop synthesis initiated by only two 3ʹ terminal base-pairs. Given sufficient time, Pol θ uses alternative pairing configurations that greatly expand the repertoire of sequence outcomes. Further primer-template adjustments yield low-fidelity outcomes when the nucleotide pool is imbalanced. Unimolecular stem loop synthesis competes with bimolecular end-joining, even when a longer terminal microhomology for end-joining is available. Both reactions are partially suppressed by the ssDNA binding protein RPA. Protein-primer grasp residues that are specific to Pol θ are needed for rapid stem-loop synthesis. The ability to perform stem-loop synthesis from a minimally paired primer is rare amongst human DNA polymerases but we show that human DNA polymerases Pol η and Pol λ can catalyze related reactions. Using purified human Pol θ, we reconstituted in vitro TMEJ incorporating an insertion arising from a stem loop extension. These activities may help explain TMEJ repair events that include inverted repeat sequences.

The complete chloroplast genome of Meconopsis torquata (Papaveraceae), a traditional Tibetan medicine

Meconopsis torquata Prain 1906, a national second-class rare and endangered plant, is reported here for the first time for its complete chloroplast genome. The genome is 153,290 bp in length, comprising a large single-copy region (LSC, 83,918 bp), a small single-copy region (SSC, 17,740 bp), and two inverted repeat sequences (IRa and IRb, each 25,816 bp). The overall GC content is 38.7%, with the IR region having the highest content (43.1%). The genome is annotated with 112 unique genes, including 4 rRNA genes, 29 tRNA genes, and 79 protein-coding genes. Analysis of codon usage bias reveals that codons ending in A/T account for 96.7% of those with a Relative Synonymous Codon Usage (RSCU) value above 1. This predominance of A/T-ending codons might be indicative of M. torquata adaptation to high-altitude environments. Phylogenetic analysis reveals a close kinship between M. torquata and M. pinnatifolia and M. paniculata, indicating that the ancestral groups of these species might have a complex evolutionary history. This study uncovers the genetic characteristics and adaptive evolution of M. torquata, offering a new perspective in understanding the phylogenetic relationships within the genus. The findings not only provide a solid theoretical foundation for the conservation and sustainable use of this rare and endangered species but also offer significant scientific support for the conservation of biodiversity.

The complete chloroplast genome sequence of leibnitzia anandria (linnaeus) turczaninow

Leibnitzia anandria is a perennial herbaceous plant with medicinal properties, and the entire plant can be used in traditional medicine. Leibnitzia anandria was once classified under the genus Gerbera Cass., but was reclassified under Leibnitzia Cass. recently. In this study, using the GeneLab M sequencing technology of the Genemind platform, we have sequenced, assembled, and analyzed the complete chloroplast genome of Leibnitzia anandria for the first time. The genome is 154168 bp in length, consisting of a large single-copy region(LSC, 80166 bp), a small single-copy region(SSC, 18202 bp), and a pair of inverted repeat sequences(IR, 27900 bp). We have predicted and annotated a total of 133 genes, including 88 protein-coding genes, 37 tRNA-coding genes, and 8 rRNA-coding genes. The results of the phylogenetic analysis indicate that Leibnitzia anandria and Leibnitzia nepalensis, as well as the closely related Gerbera plant, clustered into a separate clade, rather than grouping together with the other plants belonging to the tribe Mutisieae. This study provides new information for the phylogeny research of Leibnitzia anandria, contributing to a better understanding of its taxonomy and evolution.

Chloroplast genomes of Caragana tibetica and Caragana turkestanica: structures and comparative analysis

BackgroundThe genus Caragana encompasses multiple plant species that possess medicinal and ecological value. However, some species of Caragana are quite similar in morphology, so identifying species in this genus based on their morphological characteristics is considerably complex. In our research, illumina paired-end sequencing was employed to investigate the genetic organization and structure of Caragana tibetica and Caragana turkestanica, including the previously published chloroplast genome sequence of 7 Caragana plants.ResultsThe lengths of C. tibetica and C. turkestanica chloroplast genomes were 128,433 bp and 129,453 bp, respectively. The absence of inverted repeat sequences in these two species categorizes them under the inverted repeat loss clade (IRLC). They encode 110 and 111 genes (4 /4 rRNA genes, 30 /31tRNA genes, and 76 /76 protein-coding genes), respectively. Comparison of the chloroplast genomes of C. tibetica and C. turkestanica with 7 other Caragana species revealed a high overall sequence similarity. However, some divergence was observed between certain intergenic regions (matK-rbcL, psbD-psbM, atpA-psbI, and etc.). Nucleotide diversity (π) analysis revealed the detection of five highly likely variable regions, namely rps2-atpI, accD-psaI-ycf4, cemA-petA, psbN-psbH and rpoA-rps11. Phylogenetic analysis revealed that C. tibetica’s sister species is Caragana jubata, whereas C. turkestanica’s closest relative is Caragana arborescens.ConclusionsThe present study provides worthwhile information about the chloroplast genomes of C. tibetica and C. turkestanica, which aids in the identification and classification of Caragana species.

Promoter regulatory mode evolution enhances the high multidrug resistance of tmexCD1-toprJ1.

Antibiotic resistance could rapidly emerge from acquiring the mobile antibiotic resistance genes, which are commonly evolved from an intrinsic gene. The emergence of the plasmid-borne mobilized efflux pump gene cluster tmexCD1-toprJ1 renders the last-resort antibiotic tigecycline ineffective, although its evolutionary mechanism remains unclear. In this study, we investigate the regulatory mechanisms of the progenitor NfxB-MexCD-OprJ, a chromosomally encoded operon that does not mediate antibiotic resistance in the wild-type version, and its homologs, TNfxB1-TMexCD1-TOprJ1 mediating high-level tigecycline resistance, and TNfxB3-TMexCD3-TOprJ1. Mechanistic studies demonstrated that in nfxB-mexCD-oprJ, MexCD expression was under a weaker promoter, PmexC and inhibited by a strong repressor NfxB. For tmexCD1-toprJ1, TMexCD1 was highly expressed owing to the presence of a strong promoter, PtmexC1, and an inactive suppressor, TNfxB1, with a T39R mutation that rendered it unable to bind to promoter DNA. In tnfxB3-tmexCD3-toprJ1b, TMexCD3 expression was intermediate because of the local regulator TNfxB3, which binds to two inverted repeat sequences of PtmexC. Additionally, TNfxB3 exhibited lower protein expression and weaker DNA binding affinity than its ancestor NfxB, together with their promoter activities difference explaining the different expression levels of tmexCD-toprJ homologs. Distinct fitness burdens on these homologs-carrying bacteria were observed due to the corresponding expression level, which might be associated with their global prevalence. In summary, our data depict the mechanisms underlying the evolution and dissemination of an important mobile antibiotic resistance gene from an intrinsic chromosomal gene.IMPORTANCEAs antibiotic resistance seriously challenges global health, tigecycline is one of the few effective drugs in the pipeline against infections caused by multidrug-resistant pathogens. Our previous work identified a novel tigecycline resistance efflux pump gene cluster tmexCD1-toprJ1 in animals and humans, together with its various variants, a rising clinical concern. Herein, this study focused on how the local regulation modes of tmexCD1-toprJ1 evolved to a highly expressed efflux pump. Through comparative analysis between three tnfxB-tmexCD-toprJ homologs and their progenitor nfxB-mexCD-oprJ, modes, we demonstrated the evolutionary dynamics from a chromosomal silent gene to an active state. We found the de-repression of the local regulator and an increase of the promoter activity work together to promote a high production of drug efflux machines and enhance multidrug resistance. Our findings revealed that TMexCD1-TOprJ1 adopts a distinct evolutionary path to achieve higher multidrug resistance, urgently needing tight surveillance.

Characterization of the transcriptionally active form of dephosphorylated DctD complexed with dephospho-IIAGlc.

Bacterial enhancer-binding proteins (bEBPs) acquire a transcriptionally active state via phosphorylation. However, transcriptional activation by the dephosphorylated form of bEBP has been observed in DctD, which belongs to Group I bEBP. The formation of a complex between dephosphorylated DctD (d-DctD) and dephosphorylated IIAGlc (d-IIAGlc) is a prerequisite for the transcriptional activity of d-DctD. In the present study, characteristics of the transcriptionally active complex composed of d-IIAGlc and phosphorylation-deficient DctD (DctDD57Q) of Vibrio vulnificus were investigated in its multimeric conformation and DNA-binding ability. DctDD57Q formed a homodimer that could not bind to the DNA. In contrast, when DctDD57Q formed a complex with d-IIAGlc in a 1:1 molar ratio, it produced two conformations: dimer and dodecamer of the complex. Only the dodecameric complex exhibited ATP-hydrolyzing activity and DNA-binding affinity. For successful DNA-binding and transcriptional activation by the dodecameric d-IIAGlc/DctDD57Q complex, extended upstream activator sequences were required, which encompass the nucleotide sequences homologous to the known DctD-binding site and additional nucleotides downstream. This is the first report to demonstrate the molecular characteristics of a dephosphorylated bEBP complexed with another protein to form a transcriptionally active dodecameric complex, which has an affinity for a specific DNA-binding sequence.IMPORTANCEResponse regulators belonging to the bacterial two-component regulatory system activate the transcription initiation of their regulons when they are phosphorylated by cognate sensor kinases and oligomerized to the appropriate multimeric states. Recently, it has been shown that a dephosphorylated response regulator, DctD, could activate transcription in a phosphorylation-independent manner in Vibrio vulnificus. The dephosphorylated DctD activated transcription as efficiently as phosphorylated DctD when it formed a complex with dephosphorylated form of IIAGlc, a component of the glucose-phosphotransferase system. Functional mimicry of this complex with the typical form of transcriptionally active phosphorylated DctD led us to study the molecular characteristics of this heterodimeric complex. Through systematic analyses, it was surprisingly determined that a multimer constituted with 12 complexes gained the ability to hydrolyze ATP and recognize specific upstream activator sequences containing a typical inverted-repeat sequence flanked by distinct nucleotides.

A shorter inverted repeat length in a hairpin RNA results in a stronger silencing efficacy

RNA interference (RNAi), particularly hairpin RNA (hpRNA) constructs, is a potent genetic manipulation tool, influencing traits that are crucial for crop improvement. Typical hpRNA constructs consist of inverted repeat sequences derived from the mRNA of the target gene, separated by a non-complementary spacer. This configuration initiates the formation of a double-stranded RNA structure due to the complementary nature of sense and antisense sequences in the transcribed RNA. In this research, two hpRNA constructs, namely hp-TFL1-S and hp-TFL1-L, with distinct lengths of inverted repeats, were designed to target TFL1, a pivotal regulator of both flowering timing and inflorescence structure. hp-TFL1-L construct encompasses all the 150 nucleotide sequences of the hp-TFL1-S construct and an extra 135 nucleotides. These constructs included a chalcone synthase intron as a spacer, resulting in hairpin double-stranded RNA molecules of 150 and 285 base pair after transcription. The PCR and qRT-PCR analysis results confirmed the accurate design of both constructs and their successful suppression of TFL1 mRNA expression in transgenic Arabidopsis. Interestingly, Arabidopsis lines carrying hp-TFL1-S and hp-TFL1-L exhibited distinct phenotypic features. hp-TFL1-S induced more pronounced phenotypic changes, resulting in earlier flowering compared to the hp-TFL1-L and control group. Moreover, hp-TFL1-S had a significant impact on inflorescence architecture, distinguishing it from both hp-TFL1-L and control plants. This study provides an example where inverted repeat length alters the efficacy of silencing, and in this case, surprisingly the shorter construct appears to have a stronger silencing efficacy. This argues multiple sized hpRNAs should be tested in RNAi-based gene silencing experiments.

Comparative Analysis of the Chloroplast Genomes of Eight Species of the Genus Lirianthe Spach with Its Generic Delimitation Implications.

Based on Sima and Lu's system of the family Magnoliaceae, the genus Lirianthe Spach s. l. includes approximately 25 species, each with exceptional landscaping and horticultural or medical worth. Many of these plants are considered rare and are protected due to their endangered status. The limited knowledge of species within this genus and the absence of research on its chloroplast genome have greatly impeded studies on the relationship between its evolution and systematics. In this study, the chloroplast genomes of eight species from the genus Lirianthe were sequenced and analyzed, and their phylogenetic relationships with other genera of the family Magnoliaceae were also elucidated. The results showed that the chloroplast genome sizes of the eight Lirianthe species ranged from 159,548 to 159,833 bp. The genomes consisted of a large single-copy region, a small single-copy region, and a pair of inverted repeat sequences. The GC content was very similar across species. Gene annotation revealed that the chloroplast genomes contained 85 protein-coding genes, 37 tRNA genes, and 8 rRNA genes, totaling 130 genes. Codon usage analysis indicated that codon usage was highly conserved among the eight Lirianthe species. Repeat sequence analysis identified 42-49 microsatellite sequences, 16-18 tandem repeats, and 50 dispersed repeats, with microsatellite sequences being predominantly single-nucleotide repeats. DNA polymorphism analysis revealed 10 highly variable regions located in the large single-copy and small single-copy regions, among which rpl32-trnL, petA-psbJ, and trnH-psbA were the recommended candidate DNA barcodes for the genus Lirianthe species. The inverted repeat boundary regions show little variation between species and are generally conserved. The result of phylogenetic analysis confirmed that the genus Lirianthe s. l. is a monophyletic taxon and the most affinal to the genera, Talauma and Dugandiodendron, in Sima and Lu's system and revealed that the genus Lirianthe s. s. is paraphyletic and the genus Talauma s. l. polyphyletic in Xia's system, while Magnolia subsection Gwillimia is paraphyletic and subsection Blumiana polyphyletic in Figlar and Nooteboom's system. Morphological studies found noticeable differences between Lirianthe species in aspects including leaf indumentum, stipule scars, floral orientation, tepal number, tepal texture, and fruit dehiscence. In summary, this study elucidated the chloroplast genome evolution within Lirianthe and laid a foundation for further systematic and taxonomic research on this genus.

Transposition mechanism of ISApl1-the determinant of colistin resistance dissemination.

Multidrug-resistant Enterobacteriaceae, a prominent family of gram-negative pathogenic bacteria, causes a wide range of severe diseases. Strains carrying the mobile colistin resistance (mcr-1) gene show resistance to polymyxin, the last line of defense against multidrug-resistant gram-negative bacteria. However, the transmission of mcr-1 is not well understood. In this study, genomes of mcr-1-positive strains were obtained from the NCBI database, revealing their widespread distribution in China. We also showed that ISApl1, a crucial factor in mcr-1 transmission, is capable of self-transposition. Moreover, the self-cyclization of ISApl1 is mediated by its own encoded transposase. The electrophoretic mobility shift assay experiment validated that the transposase can bind to the inverted repeats (IRs) on both ends, facilitating the cyclization of ISApl1. Through knockout or shortening of IRs at both ends of ISApl1, we demonstrated that the cyclization of ISApl1 is dependent on the sequences of the IRs at both ends. Simultaneously, altering the ATCG content of the bases at both ends of ISApl1 can impact the excision rate by modifying the binding ability between IRs and ISAPL1. Finally, we showed that heat-unstable nucleoid protein (HU) can inhibit ISApl1 transposition by binding to the IRs and preventing ISAPL1 binding and expression. In conclusion, the regulation of ISApl1-self-circling is predominantly controlled by the inverted repeat (IR) sequence and the HU protein. This molecular mechanism deepens our comprehension of mcr-1 dissemination.

Characterization and mitigation of artifacts derived from NGS library preparation due to structure-specific sequences in the human genome

BackgroundHybridization capture-based targeted next generation sequencing (NGS) is gaining importance in routine cancer clinical practice. DNA library preparation is a fundamental step to produce high-quality sequencing data. Numerous unexpected, low variant allele frequency calls were observed in libraries using sonication fragmentation and enzymatic fragmentation. In this study, we investigated the characteristics of the artifact reads induced by sonication and enzymatic fragmentation. We also developed a bioinformatic algorithm to filter these sequencing errors.ResultsWe used pairwise comparisons of somatic single nucleotide variants (SNVs) and insertions and deletions (indels) of the same tumor DNA samples prepared using both ultrasonic and enzymatic fragmentation protocols. Our analysis revealed that the number of artifact variants was significantly greater in the samples generated using enzymatic fragmentation than using sonication. Most of the artifacts derived from the sonication-treated libraries were chimeric artifact reads containing both cis- and trans-inverted repeat sequences of the genomic DNA. In contrast, chimeric artifact reads of endonuclease-treated libraries contained palindromic sequences with mismatched bases. Based on these distinctive features, we proposed a mechanistic hypothesis model, PDSM (pairing of partial single strands derived from a similar molecule), by which these sequencing errors derive from ultrasonication and enzymatic fragmentation library preparation. We developed a bioinformatic algorithm to generate a custom mutation “blacklist” in the BED region to reduce errors in downstream analyses.ConclusionsWe first proposed a mechanistic hypothesis model (PDSM) of sequencing errors caused by specific structures of inverted repeat sequences and palindromic sequences in the natural genome. This new hypothesis predicts the existence of chimeric reads that could not be explained by previous models, and provides a new direction for further improving NGS analysis accuracy. A bioinformatic algorithm, ArtifactsFinder, was developed and used to reduce the sequencing errors in libraries produced using sonication and enzymatic fragmentation.

Chloroplast genome analyses of Caragana arborescens and Caragana opulens

BackgroundNumerous species within the genus Caragana have high ecological and medicinal value. However, species identification based on morphological characteristics is quite complicated in the genus. To address this issue, we analyzed complete plastid genome data for the genus.ResultsWe obtained chloroplast genomes of two species, Caragana arborescens and Caragana opulens, using Illumina sequencing technology, with lengths of 129,473 bp and 132,815 bp, respectively. The absence of inverted repeat sequences in the two species indicated that they could be assigned to the inverted repeat-lacking clade (IRLC). The genomes included 111 distinct genes (4 rRNA genes, 31 tRNA genes, and 76 protein-coding genes). In addition, 16 genes containing introns were identified in the two genomes, the majority of which contained a single intron. Repeat analyses revealed 129 and 229 repeats in C. arborescens and C. opulens, respectively. C. arborescens and C. opulens genomes contained 277 and 265 simple sequence repeats, respectively. The two Caragana species exhibited similar codon usage patterns. rpl20-clpP, rps19-rpl2, and rpl23-ycf2 showed the highest nucleotide diversity (pi). In an analysis of sequence divergence, certain intergenic regions (matK-rbcL, psbM-petN, atpA-psbI, petA-psbL, psbE-petL, and rps7-rps12) were highly variable. A phylogenetic analysis showed that C. arborescens and C. opulens were related and clustered together with four other Caragana species. The genera Astragalus and Caragana were relatively closely related.ConclusionsThe present study provides valuable information about the chloroplast genomes of C. arborescens and C. opulens and lays a foundation for future phylogenetic research and molecular marker development.

Restoration of the Functional nif Gene Cluster by Complex Recombination Events during Heterocyst Development in the Nitrogen-Fixing Cyanobacterium Calothrix sp. NIES-4101.

In the genome of the heterocystous cyanobacterium Calothrix sp. NIES-4101 (NIES-4101), the four genes essential for nitrogen fixation (nifB, nifH, nifD and nifK) are highly fragmented into 13 parts in a 350-kb chromosomal region, and four of these parts are encoded in the reverse strand. Such a complex fragmentation feature makes it difficult to restore the intact nifBHDK genes by the excision mechanism found in the nifD gene of the Anabaena sp. PCC 7120 heterocyst. To examine the nitrogen-fixing ability of NIES-4101, we confirmed that NIES-4101 grew well on a combined nitrogen-free medium and showed high nitrogenase activity, which strongly suggested that the complete nifBHDK genes are restored by a complex recombination process in heterocysts. Next, we resequenced the genome prepared from cells grown under nitrogen-fixing conditions. Two contigs covering the complete nifHDK and nifB genes were found by de novo assembly of the sequencing reads. In addition, the DNA fragments covering the nifBHDK operon were successfully amplified by PCR. We propose that the process of nifBHDK restoration occurs as follows. First, the nifD-nifK genes are restored by four excision events. Then, the complete nifH and nifB genes are restored by two excision events followed by two successive inversion events between the inverted repeat sequences and one excision event, forming the functional nif gene cluster, nifB-fdxN-nifS-nifU-nifH-nifD-nifK. All genes coding recombinases responsible for these nine recombination events are located close to the terminal repeat sequences. The restoration of the nifBHDK genes in NIES-4101 is the most complex genome reorganization reported in heterocystous cyanobacteria.

The complete chloroplast genome of Nekemias hypoglauca (Hance) J. Wen & Z. L. Nie 2014 (Family: Vitaceae) and its phylogenetic analysis

Nekemias is a perennial woody vine with nine species that had been originally placed in Ampelopsis. These species of Nekemias are economically and medically important. Limited information is available on the genomic characteristics of the chloroplasts of this genus. Nekemias hypoglauca (Hance) J. Wen & Z. L. Nie 2014 contains 131 unique genes (86 protein-coding genes, 8 rRNAs, and 37 tRNAs). The complete chloroplast sequence contains 162,976 bp. The large single-copy region contains 89,291 bp; the small single-copy region contains 19,063 bp, and a pair of inverted repeat sequences is composed of 27,311 bp. There are 84 simple sequence repeat (SSR) loci in the complete chloroplast genome of N. hypoglauca, with mononucleotide, dinucleotide, trinucleotide, tetranucleotide and hexanucleotide SSRs of 58, 9, 6, 10 and 1, respectively. A total of 337 repeats were identified, including 172 forward repeats, three reverse repeats and 163 palindromic repeats. A phylogenetic analysis based on the complete genome data of the chloroplasts of 10 plant species indicated the monophyly of Nekemias and determined the phylogenetic relationships of N. hypoglauca in Nekemias. This study provides a reference for further studies on the taxonomy, identification, origin and evolution of N. hypoglauca and Nekemias.

Genomic resistant determinants of multidrug-resistant Campylobacter spp. isolates in Peru

Antimicrobial resistant (AMR) Campylobacter is a global health threat; however, there is limited information on genomic determinants of resistance in low- and middle-income countries. We evaluated genomic determinants of AMR using a collection of whole genome sequenced Campylobacter jejuni and C. coli isolates from Iquitos, Peru. Campylobacter isolates from two paediatric cohort studies enriched with isolates that demonstrated resistance to ciprofloxacin and azithromycin were sequenced and mined for AMR determinants. The gyrA mutation leading to the Thr86Ile amino acid change was the only gyrA mutation associated with fluoroquinolone resistance identified. The A2075G mutation in 23S rRNA was present, but three other 23S rRNA mutations previously associated with macrolide resistance were not identified. A resistant-enhancing variant of the cmeABC efflux pump genotype (RE-cmeABC) was identified in 36.1% (35/97) of C. jejuni genomes and 17.9% (12/67) of C. coli genomes. Mutations identified in the CmeR-binding site, an inverted repeat sequence in the cmeABC promoter region that increases expression of the operon, were identified in 24/97 C. jejuni and 14/67 C. coli genomes. The presence of these variants, in addition to RE-cmeABC, was noted in 18 of the 24 C. jejuni and 9 of the 14 C. coli genomes. Both RE-cmeABC and mutations in the CmeR-binding site were strongly associated with the MDR phenotype in C. jejuni and C. coli. This is the first report of RE-cmeABC in Peru and suggests it is a major driver of resistance to the principal therapies used to treat human campylobacteriosis in this setting.