Circular Genome Research Articles

Snow viruses and their implications on red snow algal blooms.

Algal blooms can give snowmelt a red color, reducing snow albedo and creating a runaway effect that accelerates snow melting. The occurrence of red snow is predicted to grow in polar and subpolar regions with increasing global temperatures. We hypothesize that these algal blooms affect virus-bacteria interactions in snow, with potential effects on snowmelt dynamics. A genomic analysis of double-stranded DNA virus communities in red and white snow from the Whistler region of British Columbia, Canada, identified 792 putative viruses infecting bacteria. The most abundant putative snow viruses displayed low genomic similarity with known viruses. We recovered the complete circular genomes of nine putative viruses, two of which were classified as temperate. Putative snow viruses encoded genes involved in energy metabolisms, such as NAD+ synthesis and salvage pathways. In model phages, these genes facilitate increased viral particle production and lysis rates. The frequency of temperate phages was positively correlated with microbial abundance in the snow samples. These results suggest the increased frequency of temperate virus-bacteria interactions as microbial densities increase during snowmelt. We propose that this virus-bacteria dynamic may facilitate the red snow algae growth stimulated by bacteria.IMPORTANCEMicrobial communities in red snow algal blooms contribute to intensifying snowmelt rates. The role of viruses in snow during this environmental shift, however, has yet to be elucidated. Here, we characterize novel viruses extracted from snow viral metagenomes and define the functional capacities of snow viruses in both white and red snow. These results are contextualized using the composition and functions observed in the bacterial communities from the same snow samples. Together, these data demonstrate the energy metabolism performed by viruses and bacteria in a snow algal bloom, as well as expand the overall knowledge of viral genomes in extreme environments.

Codon usage bias of goose circovirus and its adaptation to host

Goose circovirus (GoCV), a potential immunosuppressive virus possessing a circular single-stranded DNA genome, is widely distributed in both domesticated and wild geese. This virus infection causes significant economic losses in the waterfowl industry. The codon usage patterns of viruses reflect the evolutionary history and genetic architecture, allowing them to adapt quickly to changes in the external environment, particularly to their hosts. In this study, we retrieved the coding sequences (Rep and Cap) and the genome of GoCV from GenBank, conducting comprehensive research to explore the codon usage patterns in 144 GoCV strains. The overall codon usage of the GoCV strains was relatively similar and exhibited a slight bias. The effective number of codons (ENC) indicated a low overall extent of codon usage bias (CUB) in GoCV. Combined with the base composition and relative synonymous codon usage (RSCU) analysis, the results revealed a bias toward A- and G-ending codons in the overall codon usage. Analysis of the ENC-GC3s plot and neutrality plot suggested that natural selection plays an important role in shaping the codon usage pattern of GoCV, with mutation pressure having a minor influence. Furthermore, the correlations between ENC and relative indices, as well as correspondence analysis (COA), showed that hydrophobicity and geographical distribution also contribute to codon usage variation in GoCV, suggesting the possible involvement of natural selection. In conclusion, GoCV exhibits comparatively slight CUB, with natural selection being the major factor shaping the codon usage pattern of GoCV. Our research contributes to a deeper understanding of GoCV evolution and its host adaptation, providing valuable insights for future basic studies and vaccine design related to GoCV.

The complete genome sequence of unculturable Mycoplasma faucium obtained through clinical metagenomic next-generation sequencing.

Diagnosing Mycoplasma faucium poses challenges, and it's unclear if its rare isolation is due to infrequent occurrence or its fastidious nutritional requirements. This study analyzes the complete genome sequence of M. faucium, obtained directly from the pus of a sternum infection in a lung transplant patient using metagenomic sequencing. Genome analysis revealed limited therapeutic options for the M. faucium infection, primarily susceptibility to tetracyclines. Three classes of mobile genetic elements were identified: two new insertion sequences, a new prophage (phiUMCG-1), and a species-specific variant of a mycoplasma integrative and conjugative element (MICE). Additionally, a Type I Restriction-Modification system was identified, featuring 5'-terminally truncated hsdS pseudogenes with overlapping repeats, indicating the potential for forming alternative hsdS variants through recombination. This study represents the first-ever acquisition of a complete circularized bacterial genome directly from a patient sample obtained from invasive infection of a primary sterile site using culture-independent, PCR-free clinical metagenomics.

Characterization of Mycoviruses in Armillaria ostoyae and A. cepistipes in the Czech Republic.

Members of the genus Armillaria are widespread forest pathogens against which effective protection has not yet been developed. Due to their longevity and the creation of large-scale cloning of Armillaria individuals, the use of mycoviruses as biocontrol agents (BCAs) against these pathogens could be an effective alternative. This work describes the detection and characterization of viruses in Armillaria spp. collected in the Czech Republic through the application of stranded total RNA sequencing. A total of five single-stranded RNA viruses were detected in Armillaria ostoyae and A. cepistipes, including viruses of the family Tymoviridae and four viruses belonging to the recently described "ambivirus" group with a circular ambisense genome arrangement. Both hammerhead (HHRz) and hairpin (HpRz) ribozymes were detected in all the ambiviricot sequences. Armillaria viruses were compared through phylogenetic analysis and confirmed their specific host by direct RT-PCR. One virus appears to infect both Armillaria species, suggesting the occurrence of interspecies transmission in nature.

First telomere-to-telomere gapless assembly of the rice blast fungus Pyricularia oryzae

Rice blast caused by Pyricularia oryzae (syn., Magnaporthe oryzae) was one of the most destructive diseases of rice throughout the world. Genome assembly was fundamental to genetic variation identification and critically impacted the understanding of its ability to overcome host resistance. Here, we report a gapless genome assembly of rice blast fungus P. oryzae strain P131 using PacBio, Illumina and high throughput chromatin conformation capture (Hi-C) sequencing data. This assembly contained seven complete chromosomes (43,237,743 bp) and a circular mitochondrial genome (34,866 bp). Approximately 14.31% of this assembly carried repeat sequences, significantly greater than its previous assembled version. This assembly had a 99.9% complement in BUSCO evaluation. A total of 14,982 genes protein-coding genes were predicted. In summary, we assembled the first telomere-to-telomere gapless genome of P. oryzae, which would be a valuable genome resource for future research on the genome evolution and host adaptation.

Comparative Analysis of Viromes Identified in Multiple Macrofungi.

Macrofungi play important roles in the soil elemental cycle of terrestrial ecosystems. Fungal viruses are common in filamentous fungi, and some of them can affect the growth and development of hosts. However, the composition and evolution of macrofungal viruses are understudied. In this study, ninety strains of Trametes versicolor, Coprinellus micaceus, Amanita strobiliformis, and Trametes hirsuta were collected in China. Four mixed pools were generated by combining equal quantities of total RNA from each strain, according to the fungal species, and then subjected to RNA sequencing. The sequences were assembled, annotated, and then used for phylogenetic analysis. Twenty novel viruses or viral fragments were characterized from the four species of macrofungi. Based on the phylogenetic analysis, most of the viral contigs were classified into ten viral families or orders: Barnaviridae, Benyviridae, Botourmiaviridae, Deltaflexiviridae, Fusariviridae, Hypoviridae, Totiviridae, Mitoviridae, Mymonaviridae, and Bunyavirales. Of these, ambi-like viruses with circular genomes were widely distributed among the studied species. Furthermore, the number and overall abundance of viruses in these four species of macrofungi (Basidiomycota) were found to be much lower than those in broad-host phytopathogenic fungi (Ascomycota: Sclerotinia sclerotiorum, and Botrytis cinerea). By employing metatranscriptomic analysis in this study, for the first time, we demonstrated the presence of multiple mycoviruses in Amanita strobiliformis, Coprinellus micaceus, Trametes hirsute, and Trametes versicolor, significantly contributing to research on mycoviruses in macrofungi.

Spreading Status of Hepatitis B Virus in Patients Attending Shone Primary Hospital from 2013-2015 E.C, Hadiya Zone, Central Ethiopia

Hepatitis virus is an inflammation of liver grouped among a member of hepadnavirus family with partially double stranded circular DNA genome, whose envelope contains protein called surface antigen. Due to existence of transmission methods within day to day activity, due to lack of knowledge, misuse of drugs of HBV and carelessness of health professionals we committed assessment of spreading status of HBV in Shone primary Hospital, to imprint knowledge among society, to recommend health professionals and concerning bodies too. This study was based on stratified and purposeful sampling. It proceeded with sex (male and female) and age group (<5, 6-15, 16-25, 26-35 and >36) as variables. During this study secondary data from 2013-2015E.C was used from the Hospital’s record office and laboratory registration books. As a result, age groups of <5, 6-15, 16-25, 26-35 and >36 were founded with 9.65%, 14.9%, 31.6 %, 24.44% and 18.41% among 114 positive patients in 2013 E.C, 9.4%, 12.75%, 41.6%, 27.5% and 8.75% among 149 positive patients in 2014 E.C and 8.5%, 16.33%, 34%, 23.53% and 17.64% among 153 positive patients in 2015 E.C respectively. And over all of 18.96% among patients examined for HBV were positive. 2.53% average increment for successive years and 5.52% and 2.2% annual increments was seen among females. Generally, lack of knowledge, unsafe blood transfusion, sexual intercourse, misuse of drugs and others were present among dwellers. So, knowledge imprinting, planning, wise use of drugs, punishing illegal, traditional peoples working with blood and contact without permission was needed to eradicate spread of disease.

Mitochondrial genome variation and intergenomic sequence transfers in Hevea species.

Among the Hevea species, rubber tree (Hevea brasiliensis) is the most important source of natural rubber. In previous studies, we sequenced the complete nuclear and chloroplast genomes of Hevea species, providing an invaluable resource for studying their phylogeny, disease resistance, and breeding. However, given that plant mitochondrial genomes are more complex and more difficult to assemble than that of the other organelles, little is known about their mitochondrial genome, which limits the comprehensive understanding of Hevea genomic evolution. In this study, we sequenced and assembled the mitochondrial genomes of four Hevea species. The four mitochondrial genomes had consistent GC contents, codon usages and AT skews. However, there were significant differences in the genome lengths and sequence repeats. Specifically, the circular mitochondrial genomes of the four Hevea species ranged from 935,732 to 1,402,206 bp, with 34-35 unique protein-coding genes, 35-38 tRNA genes, and 6-13 rRNA genes. In addition, there were 17,294-46,552 bp intergenomic transfer fragments between the chloroplast and mitochondrial genomes, consisting of eight intact genes (psaA, rrn16S, tRNA-Val, rrn5S, rrn4.5S, tRNA-Arg, tRNA-Asp, and tRNA-Asn), intergenic spacer regions and partial gene sequences. The evolutionary position of Hevea species, crucial for understanding its adaptive strategies and relation to other species, was verified by phylogenetic analysis based on the protein-coding genes in the mitochondrial genomes of 21 Malpighiales species. The findings from this study not only provide valuable insights into the structure and evolution of the Hevea mitochondrial genome but also lay the foundation for further molecular, evolutionary studies, and genomic breeding studies on rubber tree and other Hevea species, thereby potentially informing conservation and utilization strategies.

Determination of the chromosomal position effects for plug-and-play application in the Myxococcus xanthus chassis cells

The chromosomal position effect can significantly affect the transgene expression, which may provide an efficient strategy for the inauguration of alien genes in new hosts, but has been less explored rationally. The bacterium Myxococcus xanthus harbors a large circular high-GC genome, and the position effect in this chassis may result in a thousand-fold expression variation of alien natural products. In this study, we conducted transposon insertion at TA sites on the M. xanthus genome, and used enrichment and dilution indexes to respectively appraise high and low expression potentials of alien genes at insertion sites. The enrichment sites are characteristically distributed along the genome, and the dilution sites are overlapped well with the horizontal transfer genes. We experimentally demonstrated the enrichment sites as high expression integration sites (HEISs), and the dilution sites unsuitable for gene integration expression. This work highlights that HEISs are the plug-and-play sites for efficient expression of integrated genes.

The virome of the kitome: small circular virus-like genomes in laboratory reagents.

In the course of studying the virome of protozoan parasites, we identified small circular genomes resembling viruses, which turned out to be contaminants from an RNA purification kit. We report their sequences here so others can detect possible contamination in their samples by aligning them to these targets.

Complete chloroplast genome of Angelica hirsutiflora Liu et al. 1961 (Apiaceae)

Angelica hirsutiflora Liu et al.1961, is a perennial herb in the Apiaceae family that is endemic to Taiwan. In this study, the complete circular chloroplast genome of A. hirsutiflora was reconstructed and annotated using Illumina sequencing. The size of the chloroplast genome is 154,266 bp, consisting of two inverted repeats (IRs, 25,075 bp) separated by a large single-copy region (LSC, 86,569 bp) and a small single-copy region (SSC, 17,547 bp). The GC content of the chloroplast genome is 37.6%. There are 114 different genes in the chloroplast genome of A. hirsutiflora, including 80 protein-coding genes, 30 tRNA genes and four rRNA genes. A maximum-likelihood phylogenetic analysis showed that A. hirsutiflora forms a distinct clade, and separated from other species within the genus Angelica. This study provided insights into the evolutionary relationships among different species of Angelica.

Characterization of the complete chloroplast genome sequence of Lycium qingshuiheense (Solanaceae)

Lycium qingshuiheense is a typical drought and salt-alkali-tolerant plant, which has been added to the new species of Lycium in recent years. Here, we first sequenced the complete chloroplast genome of L. qingshuiheense to investigate its evolutionary relationship within the family Solanaceae. Results suggested that the circular complete chloroplast genome of L. qingshuiheense was 154,945 bp in length, including a large single-copy (LSC) of 85,930 bp, a small single-copy (SSC) of 18,203 bp, and two inverted repeats (IRs) of 25,406 bp. The GC content accounts for 37.90% and annotated 131 genes, including 86 protein-coding genes, eight rRNA genes, and 37 tRNA genes. A neighbor-joining phylogenetic tree revealed that L. qingshuiheense was a sister species to L. ruthenicum. Our study provides a new insight into the systematic evolution of Lycium in the Solanaceae family.

Genomoviruses in Liver Samples of Molossus molossus Bats.

CRESS-DNA encompasses a broad spectrum of viruses documented across diverse organisms such as animals, plants, diatoms, fungi, and marine invertebrates. Despite this prevalence, the full extent of these viruses' impact on the environment and their respective hosts remains incompletely understood. Furthermore, an increasing number of viruses within this category lack detailed characterization. This investigation focuses on unveiling and characterizing viruses affiliated with the Genomoviridae family identified in liver samples from the bat Molossus molossus. Leveraging viral metagenomics, we identified seven sequences (MmGmV-PA) featuring a circular DNA genome housing two ORFs encoding replication-associated protein (Rep) and capsid protein (Cap). Predictions based on conserved domains typical of the Genomoviridae family were established. Phylogenetic analysis revealed the segregation of these sequences into two clades aligning with the genera Gemycirculavirus (MmGmV-06-PA and MmGmV-07-PA) and Gemykibivirus (MmGmV-01-PA, MmGmV-02-PA, MmGmV-03-PA, MmGmV-05-PA, and MmGmV-09-PA). At the species level, pairwise comparisons based on complete nucleotide sequences indicated the potential existence of three novel species. In summary, our study significantly contributes to an enhanced understanding of the diversity of Genomoviridae within bat samples, shedding light on previously undiscovered viral entities and their potential ecological implications.

The complete mitochondrial genome of Leucoptera coffeella (Lepidoptera: Lyonetiidae) and phylogenetic relationships within the Yponomeutoidea superfamily

The coffee leaf miner (Leucoptera coffeella) is one of the major pests of coffee crops in the neotropical regions, and causes major economic losses. Few molecular data are available to identify this pest and advances in the knowledge of the genome of L. coffeella will contribute to improving pest identification and also clarify taxonomy of this microlepidoptera. L. coffeella DNA was extracted and sequenced using PacBio HiFi technology. Here we report the complete L. coffeella circular mitochondrial genome (16,407 bp) assembled using Aladin software. We found a total of 37 genes, including 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs) and an A + T rich-region and a D-loop. The L. coffeella mitochondrial gene organization is highly conserved with similarities to lepidopteran mitochondrial gene rearrangements (trnM-trnI-trnQ). We concatenated the 13 PCG to construct a phylogenetic tree and inferred the relationship between L. coffeella and other lepidopteran species. L. coffeella is found in the Lyonetiidae clade together with L. malifoliella and Lyonetia clerkella, both leaf miners. Interestingly, this clade is assigned in the Yponomeutoidea superfamily together with Gracillariidae, and both superfamilies displayed species with leaf-mining feeding habits.

Circular metagenome-assembled genome of Candidatus Patescibacteria recovered from anaerobic digestion sludge.

A single-contig, circular metagenome-assembled genome (cMAG) of Candidatus (Ca.) Patescibacteria was reconstructed from a mesophilic full-scale food waste treatment plant in Japan. The genome is of small size and lacks fundamental biosynthetic pathways. Taxonomic analysis using the Genome Taxonomy Database revealed that this cMAG belonged to the genus JAEZRQ01 (Ca. Parcubacteria).

Complete Mitogenome sequencing of the fish louse Argulus japonicus (Crustacea: Branchiura): Comparative analyses and phylogenetic implications.

The fish louse Argulus japonicus, a branchiuran crustacean of the Argulidae family, is attracting increasing attention because of its parasitic tendencies and significant health threats to global fish farming. The mitogenomes can yield a foundation for studying epidemiology, genetic diversity, and molecular ecology and therefore may be used to assist in the surveillance and control of A. japonicus. In this study, we sequenced and assembled the complete mitogenome of A. japonicus to shed light on its genetic and evolutionary blueprint. Our investigation indicated that the 15,045-bp circular genome of A. japonicus encodes 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), and 2 ribosomal RNAs (rRNAs) with significant AT and GC skews. Comparative genomics provided an evolutionary scenario for the genetic diversity of 13 PCGs: all were under purifying selection, with cox1 and nad6 having the lowest and highest evolutionary rates, respectively. Genome-wide phylogenetic trees established a close relationship between species of the families Argulidae (Arguloida) and Armilliferidae (Porocephalida) within Crustacea, and further, A. japonicus and Argulus americanus were determined to be more closely related to each other than to others within the family Argulidae. Single PCG-based phylogenies supported nad1 and nad6 as the best genetic markers for evolutionary and phylogenetic studies for branchiuran crustaceans due to their similar phylogenetic topologies with those of genome-based phylogenetic analyses. To sum up, these comprehensive mitogenomic data of A. japonicus and related species refine valuable marker resources and should contribute to molecular diagnostic methods, epidemiological investigations, and ecological studies of the fish ectoparasites in Crustacea.

Comparative Genomic Analysis of the Mutant Rhodotorula mucilaginosa JH-R23 Provides Insight into the High-Yield Carotenoid Mechanism

In this study, the wild-type Rhodotorula mucilaginosa GDMCC 2.30 and its high carotenoid-producing mutant JH-R23, which was screened from the space mutation breeding treated wild type, were used as materials. Through whole-genome sequencing and resequencing analysis, the carotenoid metabolic pathway and mechanism of high carotenoid production in the mutant were explored. The R. mucilaginosa GDMCC 2.30 genome comprised 18 scaffolds and one circular mitochondrial genome with a total size of 20.31 Mb, a GC content of 60.52%, and encoding 7128 genes. The mitochondrial genome comprised 40,152 bp with a GC content of 40.59%. Based on functional annotations in the GO, KEGG, and other protein databases, nine candidate genes associated with carotenoid metabolic pathways, and candidate genes of the CrtS and CrtR homologous gene families were identified. The carotenoid metabolic pathway was inferred to start from sugar metabolism to the mevalonate pathway, as is common to most fungi, and the final product of the mevalonate pathway, geranylgeranyl diphosphate, is a precursor for various carotenoids, including β-carotene, lycopene, astaxanthin, and torularhodin, formed through the activity of crucial enzymes encoded by genes such as CrtI, CrtYB, CrtS, and CrtR. Resequencing analysis of the mutant JH-R23 detected mutations in the exons of four genes, including those encoding Gal83, 3-oxoacyl-reductase, p24 proteins, and GTPase. These mutations are interpreted to have an important impact on carotenoid synthesis by JH-R23.

Isolation and identification of a novel phage targeting clinical multidrug-resistant Corynebacterium striatum isolates.

Over the past decade, Corynebacterium striatum (C. striatum), an emerging multidrug-resistant (MDR) pathogen, has significantly challenged healthcare settings, especially those involving individuals with weakened immune systems. The rise of these superbugs necessitates innovative solutions. This study aimed to isolate and characterize bacteriophages targeting MDR-C. striatum. Utilizing 54 MDR-C. striatum isolates from a local hospital as target strains, samples were collected from restroom puddles for phage screening. Dot Plaque and Double-layer plate Assays were employed for screening. A novel temperate bacteriophage, named CSP1, was identified through a series of procedures, including purification, genome extraction, sequencing, and one-step growth curves. CSP1 possesses a 39,752 base pair circular double-stranded DNA genome with HK97-like structural proteins and potential for site-specific recombination. It represents a new species within the unclassified Caudoviricetes class, as supported by transmission electron microscopy, genomic evolutionary analysis, and collinearity studies. Notably, CSP1 infected and lysed 21 clinical MDR-C. striatum isolates, demonstrating a wide host range. The phage remained stable in conditions ranging from -40 to 55°C, pH 4 to 12, and in 0.9% NaCl buffer, showing no cytotoxicity. The identification of CSP1 as the first phage targeting clinical C. striatum strains opens new possibilities in bacteriophage therapy research, and the development of diagnostic and therapeutic tools against pathogenic bacteria.

Complete genome sequences of Lactobacillus acidophilus strain P42 and Limosilactobacillus reuteri strain P43 isolated from chicken cecum.

The gut microflora contains a diverse microbial population that is influenced by the host and the environment. We report the complete circular genome sequences of Lactobacillus acidophilus strain P42 and Limosilactobacillus reuteri strain P43 isolated from chicken cecal samples. P42 and P43 could potentially serve as poultry probiotic strains.

Genome sequence and annotation of Pyramidobacter sp. strain YE332, isolated from a cattle rumen fermentation of Leucaena leaf.

We report the 2.78-Mb circular genome sequence of Pyramidobacter sp. strain YE332, isolated from a fermentation of bovine rumen fluid, supplied with leaf material from Leucaena leucocephala cv. Cunningham. This genome sequence consists of 2,795,328 bp with 60% G + C content, 2,573 predicted coding DNA sequences, and 70 RNAs.