Partial Genomic Data Research Articles

Detection of SARS-COV-2 variants and their proportions in wastewater samples using next-generation sequencing in Finland

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) variants may have different characteristics, e.g., in transmission, mortality, and the effectiveness of vaccines, indicating the importance of variant detection at the population level. Wastewater-based surveillance of SARS-CoV-2 RNA fragments has been shown to be an effective way to monitor the COVID-19 pandemic at the population level. Wastewater is a complex sample matrix affected by environmental factors and PCR inhibitors, causing insufficient coverage in sequencing, for example. Subsequently, results where part of the genome does not have sufficient coverage are not uncommon. To identify variants and their proportions in wastewater over time, we utilized next-generation sequencing with the ARTIC Network's primer set and bioinformatics pipeline to evaluate the presence of variants in partial genome data. Based on the wastewater data from November 2021 to February 2022, the Delta variant was dominant until mid-December in Helsinki, Finland’s capital, and thereafter in late December 2022 Omicron became the most common variant. At the same time, the Omicron variant of SARS-CoV-2 outcompeted the previous Delta variant in Finland in new COVID-19 cases. The SARS-CoV-2 variant findings from wastewater are in agreement with the variant information obtained from the patient samples when visually comparing trends in the sewerage network area. This indicates that the sequencing of wastewater is an effective way to monitor temporal and spatial trends of SARS-CoV-2 variants at the population level.

Partial genomic characterization of Chromobacterium piscinae from India reveals multi drug resistance.

Species of genus Chromobacterium have been isolated from diverse geographical settings, which exhibits significant metabolic flexibility as well as biotechnological and pathogenic properties. This study describes the isolation, characterization, draft assembly, and detailed sequence analysis of Chromobacterium piscinae strain W1B-CG-NIBSM isolated from water samples from multi use community pond. The organism was characterized by biochemical tests, Matrix Assisted Laser Desorption/Ionization Time of Flight Mass Spectrometry (MALDI TOF-MS) and partial genome sequencing. The partial genomic data ofChromobacterium pisciane isolate W1B NIBSM strainwas submitted to GenBank with Bio project number PRJNA803347 and accession no CP092474. An integrated genome analysis of Chromobacterium piscinae has been accomplished with PATRIC which indicates good quality genome. DNA sequencing using the illumina HiSeq 4000 system generated total length of 4,155,481bp with 63 contig with G + C content is 62.69%. This partial genome contains 4,126 protein-coding sequences (CDS), 27 repeats region and 78 transfer RNA (tRNA) genes as well as 3 ribosomal RNA (rRNA) genes. The genomic annotation of Chromobacterium W1B depicts 2,925 proteins with functional assignments and 1201 hypothetical proteins. A repertoire of specialty genes implicated in antibiotic resistance (45 genes), drug target (6 genes), Transporter (3 genes) and virulence factor (10 genes). The genomic analysis reveals the adaptability, displays metabolic varied pathways and shows specific structural complex and various virulence factors which makes this strain multi drug resistant. The isolate was found to be highly resistant to β-lactam antibiotics whereas it showed sensitivity towards aminoglycosides and fluoroquinolone antibiotics. Hence, the recovery of Chromobacterium piscinae from community pond evidenced for uncertain hidden source of public health hazard. To the best of authors knowledge this is first report of isolation and genomic description of C. piscinae from India.

Genomic features of Mycoplasma bovis subtypes currently circulating in France

BackgroundMycoplasma (M.) bovis is a major etiological agent of bovine respiratory disease, which is the most economically costly disease of cattle worldwide. Cattle disease surveillance on M. bovis is increasingly using gene-based techniques, such as multilocus sequence typing (MLST), or genome-based techniques such as core genome MLST that both require only partial genomic data. However, accurate up-to-date surveillance also demands complete, circular genomes that can be used as reference to track the evolution of the different lineages. Yet, in France, two of the main subtypes currently circulating still have no representing genome in public databases. Here, to address this gap, we provide and compare three new complete M. bovis genomes obtained from recent clinical isolates that represent major subtypes circulating in France and Europe.ResultsGenomes were obtained using a hybrid assembly strategy (Illumina and Nanopore) with fine-tuning of settings and inputs used in the Unicycler assembly pipeline, such as size selection of reads and quality trimming of the FASTQ files. The main characteristics and synteny of the genomes were compared. The three genomes mainly differed by their content in terms of mobile genetic elements, i.e. integrative conjugative elements (ICE) and insertion sequences (IS), a feature that impacts their structure. For instance, strain L15527, representing subtype3 (st3), harbours an exceptionally high number of ICEs, which results in a bigger-sized genome than all those previously described and could be associated with the propensity of st3 to gain and fix mutations through chromosomal transfer mechanisms. In contrast, strain F9160, of st1, is very close to the PG45 type strain isolated in 1961 in the USA, and harbours a huge number of IS. These features may be associated with an evolution towards a host-restricted state or in a “closed” host or environment reservoir until a recent re-emergence.ConclusionsWhole-genome comparison of the three French M. bovis subtypes provides valuable resources for future studies combining epidemiology, phylogenetic data, and phylodynamic methods.

Redefining Chlorobotryaceae as one of the principal and most diverse lineages of eustigmatophyte algae

Eustigmatophyceae is one of the ∼17 classes of the vast algal phylum Ochrophyta. Over the last decade, the eustigmatophytes emerged as an expansive group that has grown from the initially recognized handful of species to well over 200 genetically distinct entities (potential species). Yet the majority of eustigs, remain represented by unidentified strains, or even only metabarcode sequences obtained from environmental samples. Moreover, the formal classification of the group has not yet been harmonized with the recently uncovered diversity and phylogenetic relationships within the class. Here we make a major step towards resolving this issue by addressing the diversity, phylogeny and classification of one of the most prominent eustigmatophyte clades previously informally called the “Eustigmataceae group”. We obtained 18S rDNA and rbcL gene sequences from four new strains from the “Eustigmataceae group”, and from several additional eustig strains, and performed the most comprehensive phylogenetic analyses of Eustigmatophyceae to date. Our results of these analyses confirm the monophyly of the “Eustigmataceae group” and define its major subclades. We also sequenced plastid genomes of five “Eustigmataceae group” strains to not only improve our understanding of the plastid gene content evolution in eustigs, but also to obtain a robustly resolved eustigmatophyte phylogeny. With this new genomic data, we have solidified the view of the “Eustigmataceae group” as a well-defined family level clade. Crucially, we also have firmly established the genus Chlorobotrys as a member of the “Eustigmataceae group”. This new molecular evidence, together with a critical analysis of the literature going back to the 19th century, provided the basis to radically redefine the historical concept of the family Chlorobotryaceae as the formal taxonomic rubric corresponding to the “Eustigmataceae group”. With this change, the family names Eustigmataceae and Characiopsidaceae are reduced to synonymy with the Chlorobotryaceae, with the latter having taxonomic priority. We additionally studied in detail the morphology and ultrastructure of two Chlorobotryaceae members, which we describe as Neustupella aerophytica gen. et sp. nov. and Lietzensia polymorpha gen. et sp. nov. Finally, our analyses of partial genomic data from several Chlorobotryaceae representatives identified genes for hallmark flagellar proteins in all of these strains. The presence of the flagellar proteins strongly suggests that zoosporogenesis is a common trait of the family and also occurs in the members never observed to produce flagellated stages. Altogether, our work paints a rich picture of one of the most diverse principal lineages of eustigmatophyte algae.

Evolutionary Rates are Correlated Between Buchnera Endosymbionts andthe Mitochondrial Genomes of Their Aphid Hosts.

The evolution of bacterial endosymbiont genomes is strongly influenced by host-driven selection. Factors affecting host genome evolution will potentially affect endosymbiont genomes in similar ways. One potential outcome is correlations in molecular rates between the genomes of the symbiotic partners. Recently, we presented the first evidence of such correlations between the mitochondrial genomes of cockroaches and the genomes of their endosymbiont (Blattabacterium cuenoti). Here we investigate whether similar patterns are found in additional host-symbiont partners. We use partial genome data from multiple strains of the bacterial endosymbionts Buchnera aphidicola and Sulcia muelleri, and the mitochondrial genomes of their sap-feeding insect hosts. Both endosymbionts show phylogenetic congruence with the mitochondria of their hosts, a result that is expected due to their identical mode of inheritance. We compared root-to-tip distances and branch lengths of phylogenetically independent species pairs. Both analyses showed a highly significant correlation of molecular rates between thegenomes of Buchnera and themitochondrial genomes of their hosts. A similar correlation was detected between Sulcia and their hosts, but was not statistically significant. Our results indicate that evolutionary rate correlations between hosts and long-term symbionts may be a widespread phenomenon.

Genome characterization of cetaceanpox virus from a managed Indo-Pacific bottlenose dolphin (Tursiops aduncus)

Cetaceanpox viruses (CePVs) are associated with a cutaneous disease in cetaceans often referred to as “tattoo” lesions. To date, only partial genomic data are available for CePVs, and thus, they remain unclassified members of the subfamily Chordopoxvirinae within the family Poxviridae. Herein, we describe the first complete CePV genome sequenced from the tattoo lesion of a managed Indo-Pacific bottlenose dolphin (Tursiops aduncus), using next-generation sequencing. The T. aduncus CePV genome (CePV-TA) was determined to encode 120 proteins, including eight genes unique to the CePV-TA and five genes predicted to function as immune-evasion genes. The results of CePV-TA genetic analyses supported the creation of a new chordopoxvirus genus for CePVs. The complete sequencing of a CePV represents an important first step in unraveling the evolutionary relationship and taxonomy of CePVs, and significantly increases our understanding of the genomic characteristics of these chordopoxviruses.

Evolutionary rates are correlated between cockroach symbionts and mitochondrial genomes.

Bacterial endosymbionts evolve under strong host-driven selection. Factors influencing host evolution might affect symbionts in similar ways, potentially leading to correlations between the molecular evolutionary rates of hosts and symbionts. Although there is evidence of rate correlations between mitochondrial and nuclear genes, similar investigations of hosts and symbionts are lacking. Here, we demonstrate a correlation in molecular rates between the genomes of an endosymbiont (Blattabacterium cuenoti) and the mitochondrial genomes of their hosts (cockroaches). We used partial genome data for multiple strains of B. cuenoti to compare phylogenetic relationships and evolutionary rates for 55 cockroach/symbiont pairs. The phylogenies inferred for B. cuenoti and the mitochondrial genomes of their hosts were largely congruent, as expected from their identical maternal and cytoplasmic mode of inheritance. We found a correlation between evolutionary rates of the two genomes, based on comparisons of root-to-tip distances and on comparisons of the branch lengths of phylogenetically independent species pairs. Our results underscore the profound effects that long-term symbiosis can have on the biology of each symbiotic partner.

Molecular Identification of a Novel Hantavirus in Malaysian Bronze Tube-Nosed Bats (Murina aenea).

In the past ten years, several novel hantaviruses were discovered in shrews, moles, and bats, suggesting the dispersal of hantaviruses in many animal taxa other than rodents during their evolution. Interestingly, the coevolutionary analyses of most recent studies have raised the possibility that nonrodents may have served as the primordial mammalian host and harboured the ancestors of rodent-borne hantaviruses as well. The aim of our study was to investigate the presence of hantaviruses in bat lung tissue homogenates originally collected for taxonomic purposes in Malaysia in 2015. Hantavirus-specific nested RT-PCR screening of 116 samples targeting the L segment of the virus has revealed the positivity of two lung tissue homogenates originating from two individuals, a female and a male of the Murina aenea bat species collected at the same site and sampling occasion. Nanopore sequencing of hantavirus positive samples resulted in partial genomic data from S, M, and L genome segments. The obtained results indicate molecular evidence for hantaviruses in the M. aenea bat species. Sequence analysis of the PCR amplicon and partial genome segments suggests that the identified virus may represent a novel species in the Mobatvirus genus within the Hantaviridae family. Our results provide additional genomic data to help extend our knowledge about the evolution of these viruses.

Marker Development for Phylogenomics: The Case of Orobanchaceae, a Plant Family with Contrasting Nutritional Modes.

Phylogenomic approaches, employing next-generation sequencing (NGS) techniques, have revolutionized systematic and evolutionary biology. Target enrichment is an efficient and cost-effective method in phylogenomics and is becoming increasingly popular. Depending on availability and quality of reference data as well as on biological features of the study system, (semi-)automated identification of suitable markers will require specific bioinformatic pipelines. Here, we established a highly flexible bioinformatic pipeline, BaitsFinder, to identify putative orthologous single copy genes (SCGs) and to construct bait sequences in a single workflow. Additionally, this pipeline has been constructed to be able to cope with challenging data sets, such as the nutritionally heterogeneous plant family Orobanchaceae. To this end, we used transcriptome data of differing quality available for four Orobanchaceae species and, as reference, SCG data from monkeyflower (Erythranthe guttata, syn. Mimulus g.; 1,915 genes) and tomato (Solanum lycopersicum; 391 genes). Depending on whether gaps were permitted in initial blast searches of the four Orobanchaceae species against the reference, our pipeline identified 1,307 and 981 SCGs with average length of 994 bp and 775 bp, respectively. Automated bait sequence construction (using 2× tiling) resulted in 38,170 and 21,856 bait sequences, respectively. In comparison to the recently published MarkerMiner 1.0 pipeline BaitsFinder identified about 1.6 times as many SCGs (of at least 900 bp length). Skipping steps specific to analyses of Orobanchaceae, BaitsFinder was successfully used in a group of non-parasitic plants (three Asteraceae species and, as reference, SCG data from Arabidopsis thaliana based on previously compiled SCGs). Thus, BaitsFinder is expected to be broadly applicable in groups, where only transcriptomes or partial genome data of differing quality are available.

Partitioning the effects of isolation by distance, environment, and physical barriers on genomic divergence between parapatric threespine stickleback.

Genetic divergence between populations is shaped by a combination of drift, migration, and selection, yielding patterns of isolation-by-distance (IBD) and isolation-by-environment (IBE). Unfortunately, IBD and IBE may be confounded when comparing divergence across habitat boundaries. For instance, parapatric lake and stream threespine stickleback (Gasterosteus aculeatus) may have diverged due to selection against migrants (IBE), or mere spatial separation (IBD). To quantitatively partition the strength of IBE and IBD, we used recently developed population genetic software (BEDASSLE) to analyze partial genomic data from three lake-stream clines on Vancouver Island. We find support for IBD within each of three outlet streams (unlike prior studies of lake-stream stickleback). In addition, we find evidence for IBE (controlling for geographic distance): the genetic effect of habitat is equivalent to geographic separation of ∼1.9 km of IBD. Remarkably, of our three lake-stream pairs, IBE is strongest where migration between habitats is easiest. Such microgeographic genetic divergence would require exceptionally strong divergent selection, which multiple experiments have failed to detect. Instead, we propose that nonrandom dispersal (e.g., habitat choice) contributes to IBE. Supporting this conclusion, we show that the few migrants between habitats are a nonrandom subset of the phenotype distribution of the source population.

Returning to more finished genomes.

Genomic data have become commonplace in most branches of the biological sciences and have fundamentally altered the way research is conducted. However, the predominance of short-read sequence data from second-generation sequencing technologies has commonly resulted in fragmented and partial genomic data characteristics. In this opinion, I will highlight how long, unbiased reads from single molecule, real-time (SMRT) sequencing now allow for a return to more contiguous and comprehensive views of genomes.

Reconstruction of Local Biochemical Reaction Network Based on Human Chromosome 9 Sequence Data

Metabolic networks are complex and highly interconnected, thus systems-level computational approaches are required to elucidate and to understand metabolic genotype-phenotype relationships. This paper has manually reconstructed the local human metabolic network based on DNA sequence data of human chromosome nine. Herein the paper describes the reconstruction process and discusses how the resulting chromosome-scale (or local) network differs from genome-scale ones. The underestimated results have revealed many gaps in the current understanding of human metabolism that require future experimental investigation. They also suggest possible problems arising from local reconstruction based on partial genome data. The study suggests further applications enabled by reconstruction of human metabolic network. The establishment of this network represents a step toward genome-scale human systems biology.

The Position of Gnetales among Seed Plants: Overcoming Pitfalls of Chloroplast Phylogenomics

The phylogenetic position of Gnetales is one of the most contentious issues in the seed plant systematics. To elucidate the Gnetales position, an improved amino acid substitution matrix was estimated based on 64 chloroplast (cp) genomes and was applied to cp genome data including all three lineages of Gnetales in maximum likelihood analyses of proteins. Although the initial analysis strongly supported the sister relation of Gnetales with Cryptomeria (Cupressophyta or non-Pinaceae conifers) (the "Gnecup" hypothesis), the support seems to be caused by a long-branch attraction (LBA) artifact. Indeed, by removing fastest evolving proteins that are most likely associated with the LBA, the support drastically declined. Furthermore, another analysis of partial genome data with dense taxon sampling of conifers showed that, in psbC, rpl2, and rps7 proteins, there are many parallel amino acid substitutions between the lineages leading to Gnetales and to Cryptomeria, and by further excluding these three genes, the sister relation of Gnetales with Pinaceae (the "Gnepine" hypothesis) became supported. Overall, our analyses indicate that the LBA and parallel substitutions cause a seriously biased inference of phylogenetic position of Gnetales with the cp genome data.

Using homolog groups to create a whole-genomic tree of free-living organisms: an update.

Genomic trees have been constructed based on the presence and absence of families of protein-encoding genes observed in 27 complete genomes, including genomes of 15 free-living organisms. This method does not rely on the identification of suspected orthologs in each genome, nor the specific alignment used to compare gene sequences because the protein-encoding gene families are formed by grouping any protein with a pairwise similarity score greater than a preset value. Because of this all inclusive grouping, this method is resilient to some effects of lateral gene transfer because transfers of genes are masked when the recipient genome already has a homolog (not necessarily an ortholog) of the incoming gene. Of 71 genes suspected to have been laterally transferred to the genome of Aeropyrum pernix, only approximately 7 to 15 represent genes where a lateral gene transfer appears to have generated homoplasy in our character dataset. The genomic tree of the 15 free-living taxa includes six different bacterial orders, six different archaeal orders, and two different eukaryotic kingdoms. The results are remarkably similar to results obtained by analysis of rRNA. Inclusion of the other 12 genomes resulted in a tree only broadly similar to that suggested by rRNA with at least some of the differences due to artifacts caused by the small genome size of many of these species. Very small genomes, such as those of the two Mycoplasma genomes included, fall to the base of the Bacterial domain, a result expected due to the substantial gene loss inherent to these lineages. Finally, artificial "partial genomes" were generated by randomly selecting ORFs from the complete genomes in order to test our ability to recover the tree generated by the whole genome sequences when only partial data are available. The results indicated that partial genomic data, when sampled randomly, could robustly recover the tree generated by the whole genome sequences.

Comparative genomics for identification of clone-specific sequence blocks in Streptococcus pneumoniae

The partial genome sequences of a serotype 3 and a serotype 2 pneumococcal strain were compared to the complete type 4 pneumococcal genome. Over 500 000 and 150 000 base pairs of the partial genome data, obtained from published patents, were analysed respectively. Global alignment showed that nearly the whole genome is highly conserved in accordance with data of multilocus sequence typing of housekeeping genes. The search for clone-specific genes revealed 17 new open reading frames in the type 3 strain, while no new open reading frame was detected in the type 2 strain. Allelic variation of genes was restricted by the use of crude sequence data, but still permitted identification of some new alleles and the observation that all surface proteins present in the partial genome data were highly conserved. In both strains we observed also a variety of chromosomal rearrangements and variations due to mobile genetic elements. All together, this comparative genomic approach gives a genome-based overview of strain relatedness and a prospective on what could be expected when sequencing other pneumococcal strains.