Genome Assembly Research Articles

Bioinformatics analysis of clinical isolates of Escherichia coli from Nigeria and other countries for quinolone and fluoroquinolone resistance

Quinolones were the drugs of choice for the treatment of bacterial infections, especially, infections caused by Gram negative bacteria. Unfortunately, these drugs have been resisted by the microbial agents including Escherichia coli, known to be the leading cause of Urinary Tract Infection (UTI) globally. This study therefore was aimed at detecting all the genes involved in Quinolone resistance by the E. coli pathogen isolated from Nigeria and from other geographical regions, using robust techniques. Twenty-three sequence data files of Escherichia coli from various countries of the world were retrieved from the National Centre for biotechnology Information (NCBI) database and sent for genome assembly for processing of the short reads into long reads. The outcome was config. fasta files which were comprehensively annotated and characterized for genetic functions and mechanisms. A total of 208 antibiotic resistance genes were detected, out of which 27(13.0%) were linked to quinolone resistance and 14(6.7%) to multi- drug resistance. The result of this study significantly implicated many genes in quinolone resistance; notably were the efflux pump genes and their high percentage abundance. We recommend in-depth study of the genes for their expression capabilities, also the structure and features of the efflux pump genes to enable proper redesigning of drugs by integrating anti efflux pump substances that will selectively prevent the expression of the genes for antibiotic resistance, without any harm to the host, or that can destabilize the positive regulation of the operon for antibiotic resistance.

Draft genome sequences of three Lactobacillus crispatus strains, isolated from the female urinary tract.

Lactobacillus crispatus is a frequent member of the female urogenital microbiota. Here, we present the draft genome assemblies of three L. crispatus strains: UMB4356, UMB5661, and UMB6244. All strains were isolated from voided urine samples from females with type 2 diabetes.

Complete genome sequences of Clostridioides difficile surveillance isolates representing the top 10 ribotypes from the Emerging Infections Program, United States, 2016.

Ten Clostridioides difficile isolates representing the top 10 ribotypes collected in 2016 through the Emerging Infections Program underwent long-read sequencing to obtain high-quality reference genome assemblies. These isolates are publicly available through the CDC & FDA Antibiotic Resistance Isolate Bank.

Sensitive Method To Analyze Cell Surface GPI-Anchored Proteins Using DNA Hybridization Chain Reaction-Mediated Signal Amplification.

GPI-anchored proteins (GPI-APs) are ubiquitous and essential but exist in low abundances on the cell surface, making their analysis and investigation especially challenging. To tackle the problem, a new method to detect and study GPI-APs based upon GPI metabolic engineering and DNA-facilitated fluorescence signal amplification was developed. In this context, cell surface GPI-APs were metabolically engineered using azido-inositol derivatives to introduce an azido group. This allowed GPI-AP coupling with alkyne-functionalized multifluorophore DNA assemblies generated by hybridization chain reaction (HCR). It was demonstrated that this approach could significantly improve the detection limit and sensitivity of GPI-APs, thereby enabling various biological studies, including the investigation of live cells. This new, enhanced GPI-AP detection method has been utilized to successfully explore GPI-AP engineering, analyze GPI-APs, and profile GPI-AP expression in different cells.

First draft reference genome and annotation of the alternative oil species Physaria fendleri.

In the wake of increasing demand for renewable energy sources, plant-based sources including alternative oilseeds have come to the forefront of interest. Hydroxy fatty acids (HFAs), produced in a few oilseed species, are important chemical feedstocks for industrial applications. An integrated approach was taken to assemble the first draft genome of the alternative HFA producer Physaria fendleri (n = 6), an outcrossing species with high heterozygosity. Both de novo transcriptome assemblies and genome assemblies were produced with public and generated sequencing reads. Resulting intermediate assemblies were then scaffolded and patched with multiple data sources, followed by super-scaffolding onto a masked genome of Camelina laxa (n = 6). Despite a current lack of available resources for the physical mapping of genomic scaffolds of Physaria fendleri, topography of the genome with respect to repeat and gene content was preserved at the scaffold level and not significantly lost via super-scaffolding. Read representation, gene and genome completion statistics, and annotation results illustrated the creation of a functional draft genome and a tool for future research on alternative oil species.

A chromosome-level genome assembly of the disco clam, Ctenoides ales.

The bivalve subclass Pteriomorphia, which includes the economically important scallops, oysters, mussels, and ark clams, exhibits extreme ecological, morphological, and behavioral diversity. Among this diversity are five morphologically distinct eye types, making Pteriomorphia an excellent setting to explore the molecular basis for the evolution of novel traits. Of pteriomorphian bivalves, Limida is the only order lacking genomic resources, greatly limiting the potential phylogenomic analyses related to eyes and phototransduction. Here, we present a limid genome assembly, the disco clam, Ctenoides ales, which is characterized by invaginated eyes, exceptionally long tentacles, and a flashing light display. This genome assembly was constructed with PacBio long reads and Dovetail Omni-CTM proximity-ligation sequencing. The final assembly is ∼2.3Gb and over 99% of the total length is contained in 18 pseudomolecule scaffolds. We annotated 41,064 protein coding genes and report a BUSCO completeness of 91.9% for metazoa_obd10. Additionally, we report a complete and annotated mitochondrial genome, which also had been lacking from Limida. The ∼20Kb mitogenome has 12 protein coding genes, 22 tRNAs, 2 rRNA genes, and a 1,589 bp duplicated sequence containing the origin of replication. The C. ales nuclear genome size is substantially larger than other pteriomorphian genomes, mainly accounted for by transposable element sequences. We inventoried the genome for opsins, the signaling proteins that initiate phototransduction, and found that, unlike its closest eyed-relatives, the scallops, C. ales lacks duplication of the rhabdomeric Gq-protein coupled opsin that is typically used for invertebrate vision. In fact, C. ales has uncharacteristically few opsins relative to the other pteriomorphian families, all of which have unique expansions of xenopsins, a recently discovered opsin subfamily. This chromosome-level assembly, along with the mitogenome, will be valuable resources for comparative genomics and phylogenetics in bivalves and particularly for the understudied but charismatic limids.

Complete Mitogenome of “Pumpo” (Bos taurus), a Top Bull from a Peruvian Genetic Nucleus, and Its Phylogenetic Analysis

The mitochondrial genome of Pumpo (Bos taurus), a prominent breed contributing to livestock farming, was sequenced using the Illumina HiSeq 2500 platform. Assembly and annotation of the mitochondrial genome were achieved through a multifaceted approach employing bioinformatics tools such as Trim Galore, SPAdes, and Geseq, followed by meticulous manual inspection. Additionally, analyses covering tRNA secondary structure and codon usage bias were conducted for comprehensive characterization. The 16,341 base pair mitochondrial genome comprises 13 protein-coding genes, 22 tRNA genes, and 2 rRNA genes. Phylogenetic analysis places Pumpo within a clade predominantly composed of European cattle, reflecting its prevalence in Europe. This comprehensive study underscores the importance of mitochondrial genome analysis in understanding cattle evolution and highlights the potential of genetic improvement programs in livestock farming, thus contributing to enhanced livestock practices.

Haplotype-resolved genome assembly and implementation of VitExpress, an open interactive transcriptomic platform for grapevine

Haplotype-resolved genome assemblies were produced for Chasselas and Ugni Blanc, two heterozygous Vitis vinifera cultivars by combining high-fidelity long-read sequencing and high-throughput chromosome conformation capture (Hi-C). The telomere-to-telomere full coverage of the chromosomes allowed us to assemble separately the two haplo-genomes of both cultivars and revealed structural variations between the two haplotypes of a given cultivar. The deletions/insertions, inversions, translocations, and duplications provide insight into the evolutionary history and parental relationship among grape varieties. Integration of de novo single long-read sequencing of full-length transcript isoforms (Iso-Seq) yielded a highly improved genome annotation. Given its higher contiguity, and the robustness of the IsoSeq-based annotation, the Chasselas assembly meets the standard to become the annotated reference genome for V. vinifera. Building on these resources, we developed VitExpress, an open interactive transcriptomic platform, that provides a genome browser and integrated web tools for expression profiling, and a set of statistical tools (StatTools) for the identification of highly correlated genes. Implementation of the correlation finder tool for MybA1, a major regulator of the anthocyanin pathway, identified candidate genes associated with anthocyanin metabolism, whose expression patterns were experimentally validated as discriminating between black and white grapes. These resources and innovative tools for mining genome-related data are anticipated to foster advances in several areas of grapevine research.

Klumpy: A tool to evaluate the integrity of long-read genome assemblies and illusive sequence motifs.

The improvement and decreasing costs of third-generation sequencing technologies has widened the scope of biological questions researchers can address with de novo genome assemblies. With the increasing number of reference genomes, validating their integrity with minimal overhead is vital for establishing confident results in their applications. Here, we present Klumpy, a tool for detecting and visualizing both misassembled regions in a genome assembly and genetic elements (e.g. genes) of interest in a set of sequences. By leveraging the initial raw reads in combination with their respective genome assembly, we illustrate Klumpy's utility by investigating antifreeze glycoprotein (afgp) loci across two icefishes, by searching for a reported absent gene in the northern snakehead fish, and by scanning the reference genomes of a mudskipper and bumblebee for misassembled regions. In the two former cases, we were able to provide support for the noncanonical placement of an afgp locus in the icefishes and locate the missing snakehead gene. Furthermore, our genome scans were able identify an unmappable locus in the mudskipper reference genome and identify a putative repetitive element shared among several species of bees.

Chromosome-scale Genome Assembly and Characterization of Top-quality Japanese Green Tea Cultivar 'Seimei'.

Japanese green tea, an essential beverage in Japanese culture, is characterized by the initial steaming of freshly harvested leaves during production. This process efficiently inactivates endogenous enzymes such as polyphenol oxidases, resulting in the production of sencha, gyokuro, and matcha that preserves the vibrant green color of young leaves. Although genome sequences of several tea cultivars and germplasms have been published, no reference genome sequences are available for Japanese green tea cultivars. Here, we constructed a reference genome sequence of the cultivar 'Seimei', which is used to produce high-quality Japanese green tea. Using the PacBio HiFi and Hi-C technologies for chromosome-scale genome assembly, we obtained 15 chromosome sequences with a total genome size of 3.1 Gb and an N50 of 214.9 Mb. By analyzing the genomic diversity of 23 Japanese tea cultivars and lines, including the leading green tea cultivars 'Yabukita' and 'Saemidori', revealed several candidate genes that could be related to the characteristics of Japanese green tea. The reference genome of 'Seimei' and information on genomic diversity of Japanese green tea cultivars should provide crucial information for effective breeding of such cultivars in the future.

Genome assembly of M. spongiola and comparative genomics of the genus Morchella provide initial insights into taxonomy and adaptive evolution.

Morchella spongiola is a highly prized mushroom for its delicious flavor and medical value and is one of the most flourishing, representative, and dominant macrofungi in the Qilian Mountains of the Qinghai-Tibet Plateau subkingdoms (QTPs). However, the understanding of M. spongiola remains largely unknown, and its taxonomy is ambiguous. In this study, we redescribed a unique species of M. spongiola, i.e., micromorphology, molecular data, genomics, and comparative genomics, and the historical biogeography of M. spongiola were estimated for 182 single-copy homologous genes. A high-quality chromosome-level reference genome of M. spongiola M12-10 was obtained by combining PacBio HiFi data and Illumina sequencing technologies; it was approximately 57.1Mb (contig N50 of 18.14Mb) and contained 9775 protein-coding genes. Comparative genome analysis revealed considerable conservation and unique characteristics between M. spongiola M12-10 and 32 other Morchella species. Molecular phylogenetic analysis indicated that M. spongiola M12-10 is similar to the M. prava/Mes-7 present in sandy soil near rivers, differentiating from black morels ~ 43.06 Mya (million years ago), and diverged from M. parva/Mes-7 at approximately 12.85 Mya (in the Miocene epoch), which is closely related to the geological activities in the QTPs (in the Neogene). Therefore, M. spongiola is a unique species rather than a synonym of M. vulgaris/Mes-5, which has a distinctive grey-brown sponge-like ascomata. This genome of M. spongiola M12-10 is the first published genome sequence of the species in the genus Morchella from the QTPs, which could aid future studies on functional gene identification, germplasm resource management, and molecular breeding efforts, as well as evolutionary studies on the Morchella taxon in the QTPs.

Compression algorithm for colored de Bruijn graphs

A colored de Bruijn graph (also called a set of k-mer sets), is a set of k-mers with every k-mer assigned a set of colors. Colored de Bruijn graphs are used in a variety of applications, including variant calling, genome assembly, and database search. However, their size has posed a scalability challenge to algorithm developers and users. There have been numerous indexing data structures proposed that allow to store the graph compactly while supporting fast query operations. However, disk compression algorithms, which do not need to support queries on the compressed data and can thus be more space-efficient, have received little attention. The dearth of specialized compression tools has been a detriment to tool developers, tool users, and reproducibility efforts. In this paper, we develop a new tool that compresses colored de Bruijn graphs to disk, building on previous ideas for compression of k-mer sets and indexing colored de Bruijn graphs. We test our tool, called ESS-color, on various datasets, including both sequencing data and whole genomes. ESS-color achieves better compression than all evaluated tools and all datasets, with no other tool able to consistently achieve less than 44% space overhead. The software is available at http://github.com/medvedevgroup/ESSColor.

DNA assembly by multi-fragment digestion/ligation and homologous recombination

To develop an accurate and efficient protocol for multi-fragment assembly and multi-site mutagenesis, we integrated and optimized the common multi-fragment assembly methods and validated the established method by using fructose-1,6-diphosphatase 1 (FBP1) with 4 mutant sites. The fragments containing mutations were assembled by introducing mutant sites and Bsa I recognition sequences. After digestion/ligation, the ligated fragment was amplified with the primers containing overlap region to the linearized vector. The amplified fragment was ligated to the linearized vector and the ligation product was transformed into Escherichia coli. After screening and sequencing, the recombinant plasmid with 4 mutant sites was obtained. This protocol overcame the major defects of Gibson assembly and Golden Gate assembly, serving as an efficient solution for multi-fragment assembly and multi-site mutagenesis.

ACMGA: a reference-free multiple-genome alignment pipeline for plant species

BackgroundThe short-read whole-genome sequencing (WGS) approach has been widely applied to investigate the genomic variation in the natural populations of many plant species. With the rapid advancements in long-read sequencing and genome assembly technologies, high-quality genome sequences are available for a group of varieties for many plant species. These genome sequences are expected to help researchers comprehensively investigate any type of genomic variants that are missed by the WGS technology. However, multiple genome alignment (MGA) tools designed by the human genome research community might be unsuitable for plant genomes.ResultsTo fill this gap, we developed the AnchorWave-Cactus Multiple Genome Alignment (ACMGA) pipeline, which improved the alignment of repeat elements and could identify long (> 50 bp) deletions or insertions (INDELs). We conducted MGA using ACMGA and Cactus for 8 Arabidopsis (Arabidopsis thaliana) and 26 Maize (Zea mays) de novo assembled genome sequences and compared them with the previously published short-read variant calling results. MGA identified more single nucleotide variants (SNVs) and long INDELs than did previously published WGS variant callings. Additionally, ACMGA detected significantly more SNVs and long INDELs in repetitive regions and the whole genome than did Cactus. Compared with the results of Cactus, the results of ACMGA were more similar to the previously published variants called using short-read. These two MGA pipelines identified numerous multi-allelic variants that were missed by the WGS variant calling pipeline.ConclusionsAligning denovo assembled genome sequences could identify more SNVs and INDELs than mapping short-read. ACMGA combines the advantages of AnchorWave and Cactus and offers a practical solution for plant MGA by integrating global alignment, a 2-piece-affine-gap cost strategy, and the progressive MGA algorithm.

The reference genome of Macropodus opercularis (the paradise fish)

Amongst fishes, zebrafish (Danio rerio) has gained popularity as a model system over most other species and while their value as a model is well documented, their usefulness is limited in certain fields of research such as behavior. By embracing other, less conventional experimental organisms, opportunities arise to gain broader insights into evolution and development, as well as studying behavioral aspects not available in current popular model systems. The anabantoid paradise fish (Macropodus opercularis), an “air-breather” species has a highly complex behavioral repertoire and has been the subject of many ethological investigations but lacks genomic resources. Here we report the reference genome assembly of M. opercularis using long-read sequences at 150-fold coverage. The final assembly consisted of 483,077,705 base pairs (~483 Mb) on 152 contigs. Within the assembled genome we identified and annotated 20,157 protein coding genes and assigned ~90% of them to orthogroups.

Conifers concentrate large numbers of NLR immune receptor genes on one chromosome.

Nucleotide-binding domain and Leucine-rich Repeat (NLR) immune receptor genes form a major line of defence in plants, acting in both pathogen recognition and resistance machinery activation. NLRs are reported to form large gene clusters in limber pine (Pinus flexilis) but it is unknown how widespread this genomic architecture may be among the extant species of conifers (Pinophyta). We used comparative genomic analyses to assess patterns in the abundance, diversity and genomic distribution of NLR genes. Chromosome-level whole genome assemblies and high-density linkage maps in the Pinaceae, Cupressaceae, Taxaceae and other gymnosperms were scanned for NLR genes using existing and customised pipelines. Discovered genes were mapped across chromosomes and linkage groups and analysed phylogenetically for evolutionary history. Conifer genomes are characterised by dense clusters of NLR genes, highly localised on one chromosome. These clusters are rich in TNL-encoding genes, which seem to have formed through multiple tandem duplication events. In contrast to angiosperms and non-coniferous gymnosperms, genomic clustering of NLR genes is ubiquitous in conifers. NLR-dense genomic regions are likely to influence a large part of the plant's resistance, informing our understanding of adaptation to biotic stress and the development of genetic resources through breeding.

Chromosome-Scale, De Novo, Phased Genome Assemblies of Three Australian Limes: Citrus australasica, C. inodora, and C. glauca

Huanglongbing (HLB) is a severe citrus disease worldwide. Wild Australian limes like Citrus australasica, C. inodora, and C. glauca possess beneficial HLB resistance traits. Individual trees of the three taxa were extensively used in a breeding program for over a decade to introgress resistance traits into commercial-quality citrus germplasm. We generated high-quality, phased, de novo genome assemblies of the three Australian limes using PacBio long-read sequencing. The genome assembly sizes of the primary and alternate haplotypes were determined for C. australasica (337 Mb/335 Mb), C. inodora (304 Mb/299 Mb), and C. glauca (376 Mb/379 Mb). The nine chromosome-scale scaffolds included 86–91% of the genome sequences generated. The integrity and completeness of the assembled genomes were estimated to be at 97.2–98.8%. Gene annotation studies identified 25,461 genes in C. australasica, 27,665 in C. inodora, and 30,067 in C. glauca. Genes belonging to 118 orthogroups were specific to Australian lime genomes compared to other citrus genomes analyzed. Significantly fewer canonical resistance (R) genes were found in C. inodora and C. glauca (319 and 449, respectively) compared to C. australasica (576), C. clementina (579), and C. sinensis (651). Similar patterns were observed for other gene families associated with potential HLB resistance, including Phloem protein 2 (PP2) and Callose synthase (CalS) genes predicted in the Australian lime genomes. The genomic information on Australian limes developed in the present study will help understand the genetic basis of HLB resistance.

Bacterial dynamics of the plastisphere microbiome exposed to sub-lethal antibiotic pollution

BackgroundAntibiotics and microplastics are two major aquatic pollutants that have been associated to antibiotic resistance selection in the environment and are considered a risk to human health. However, little is known about the interaction of these pollutants at environmental concentrations and the response of the microbial communities in the plastisphere to sub-lethal antibiotic pollution. Here, we describe the bacterial dynamics underlying this response in surface water bacteria at the community, resistome and mobilome level using a combination of methods (next-generation sequencing and qPCR), sequencing targets (16S rRNA gene, pre-clinical and clinical class 1 integron cassettes and metagenomes), technologies (short and long read sequencing), and assembly approaches (non-assembled reads, genome assembly, bacteriophage and plasmid assembly).ResultsOur results show a shift in the microbial community response to antibiotics in the plastisphere microbiome compared to surface water communities and describe the bacterial subpopulations that respond differently to antibiotic and microplastic pollution. The plastisphere showed an increased tolerance to antibiotics and selected different antibiotic resistance bacteria (ARB) and antibiotic resistance genes (ARGs). Several metagenome assembled genomes (MAGs) derived from the antibiotic-exposed plastisphere contained ARGs, virulence factors, and genes involved in plasmid conjugation. These include Comamonas, Chryseobacterium, the opportunistic pathogen Stenotrophomonas maltophilia, and other MAGs belonging to genera that have been associated to human infections, such as Achromobacter. The abundance of the integron-associated ciprofloxacin resistance gene aac(6’)-Ib-cr increased under ciprofloxacin exposure in both freshwater microbial communities and in the plastisphere. Regarding the antibiotic mobilome, although no significant changes in ARG load in class 1 integrons and plasmids were observed in polluted samples, we identified three ARG-containing viral contigs that were integrated into MAGs as prophages.ConclusionsThis study illustrates how the selective nature of the plastisphere influences bacterial response to antibiotics at sub-lethal selective pressure. The microbial changes identified here help define the selective role of the plastisphere and its impact on the maintenance of environmental antibiotic resistance in combination with other anthropogenic pollutants. This research highlights the need to evaluate the impact of aquatic pollutants in environmental microbial communities using complex scenarios with combined stresses.3sW_WiGiki6wh61j1XWikeVideo

The genome sequence of a hoverfly Eristalinus aeneus (Scopoli, 1763)

We present a genome assembly from an individual female Eristalinus aeneus (a hoverfly; Arthropoda; Insecta; Diptera; Syrphidae). The genome sequence is 495.4 megabases in span. Most of the assembly is scaffolded into 6 chromosomal pseudomolecules. The mitochondrial genome has also been assembled and is 15.97 kilobases in length.

The genome sequence of the Common Carpet moth, Epirrhoe alternata (Müller, 1764)

We present a genome assembly from an individual male Epirrhoe alternata (the Common Carpet; Arthropoda; Insecta; Lepidoptera; Geometridae). The genome sequence is 358.5 megabases in span. Most of the assembly is scaffolded into 31 chromosomal pseudomolecules, including the Z sex chromosome. The mitochondrial genome has also been assembled and is 16.99 kilobases in length.