Reference-guided de novo assembly approach improves genome reconstruction for related species

Production and characterization of bacterial cellulose synthesized by Enterobacter chuandaensis strain AEC using Phoenix dactylifera and Musa acuminata

Bacterial cellulose (BC) is a biopolymer synthesized extracellularly by certain bacteria through the polymerization of glucose monomers. This study aimed to produce BC using Enterobacter chuandaensis with fruit extracts from Phoenix dactylifera (D) and Musa acuminata (M) as carbon sources. Attenuated Total Reflectance-Fourier Transform Infrared Spectroscopy (ATR-FTIR) showed characteristic cellulose vibrations, while X-ray diffraction (XRD) identified distinct peaks at 15.34°, 19.98°, 22.58°, and 34.6°, confirming the cellulose structure. Whole-genome sequencing of E. chuandaensis identified key genes involved in BC production. The BC produced then exhibited a molecular weight of 1,857,804 g/mol, with yields of 2.8 g/L and 2.5 g/L for treatments D and M, respectively. The crystallinity index of the purified BC was 74.1, and 13C NMR analysis confirmed the dominant cellulose Iα crystalline form. The BC showed high biocompatibility in cytotoxicity assays, with cell viability between 92% and 100%, indicating its potential for use in biomedical applications. This investigation represents the first report of BC production by E. chuandaensis, which promises a new avenue for sustainable and efficient BC synthesis using fruit extracts as carbon sources.

Production and analysis of synthesized bacterial cellulose by Enterococcus faecalis strain AEF using Phoenix dactylifera and Musa acuminata fruit extracts

Bacterial cellulose (BC) is a highly versatile biopolymer renowned for its exceptional mechanical strength, water retention, and biocompatibility. These properties make it a valuable material for various industrial and biomedical applications. In this study, Enterococcus faecalis synthesized extracellular BC, utilizing Phoenix dactylifera and Musa acuminata fruit extracts as sustainable carbon sources. LC-MS analysis identified glucose as the primary carbohydrate in these extracts, providing a suitable substrate for BC production. Scanning Electron Microscopy (SEM) revealed a network of BC nanofibers on Congo red agar plates. ATR-FTIR spectroscopy confirmed the presence of characteristic cellulose functional groups, further supporting BC synthesis. X-ray diffraction (XRD) analysis indicated a high crystallinity index of 71%, consistent with the cellulose I structure, as evidenced by peaks at 16.22°, 21.46°, 22.52°, and 34.70°. Whole-genome sequencing of E. faecalis identified vital genes involved in BC biosynthesis, including bcsA, bcsB, diguanylate cyclase (DGC), and 6-phosphofructokinase (pfkA). Antibiotic susceptibility tests revealed resistance to cefotaxime, ceftazidime, and ceftriaxone, while susceptibility to imipenem was observed. Quantitative assessment demonstrated that higher concentrations of fruit extracts (5.0-20 mg/mL) significantly enhanced BC production. Cytotoxicity testing via the MTT assay confirmed excellent biocompatibility with NIH/3T3 fibroblast cells, showing high cell viability (97-105%). Unlike commonly studied Gram-negative bacteria like Acetobacter xylinum for BC production, this research focuses on Gram-positive Enterococcus faecalis and utilizes Phoenix dactylifera and Musa acuminata fruit extracts as carbon sources. This approach offers a sustainable and promising avenue for BC production.

A metagenomic approach to demystify the anaerobic digestion black box and achieve higher biogas yield: a review

The increasing reliance on fossil fuels and the growing accumulation of organic waste necessitates the exploration of sustainable energy alternatives. Anaerobic digestion (AD) presents one such solution by utilizing secondary biomass to produce biogas while reducing greenhouse gas emissions. Given the crucial role of microbial activity in anaerobic digestion, a deeper understanding of the microbial community is essential for optimizing biogas production. While metagenomics has emerged as a valuable tool for unravelling microbial composition and providing insights into the functional potential in biodigestion, it falls short of interpreting the functional and metabolic interactions, limiting a comprehensive understanding of individual roles in the community. This emphasizes the significance of expanding the scope of metagenomics through innovative tools that highlight the often-overlooked, yet crucial, role of microbiota in biomass digestion. These tools can more accurately elucidate microbial ecological fitness, shared metabolic pathways, and interspecies interactions. By addressing current limitations and integrating metagenomics with other omics approaches, more accurate predictive techniques can be developed, facilitating informed decision-making to optimize AD processes and enhance biogas yields, thereby contributing to a more sustainable future.

Soil microbiome characterization and its future directions with biosensing

Soil microbiome characterization is typically achieved with next-generation sequencing (NGS) techniques. Metabarcoding is very common, and meta-omics is growing in popularity. These techniques have been instrumental in microbiology, but they have limitations. They require extensive time, funding, expertise, and computing power to be effective. Moreover, these techniques are restricted to controlled laboratory conditions; they are not applicable in field settings, nor can they rapidly generate data. This hinders using NGS as an environmental monitoring tool or an in-situ checking device. Biosensing technology can be applied to soil microbiome characterization to overcome these limitations and to complement NGS techniques. Biosensing has been used in biomedical applications for decades, and many successful commercial products are on the market. Given its previous success, biosensing has much to offer soil microbiome characterization. There is a great variety of biosensors and biosensing techniques, and a few in particular are better suited for soil field studies. Aptamers are more stable than enzymes or antibodies and are more ready for field-use biosensors. Given that any microbiome is complex, a multiplex sensor will be needed, and with large, complicated datasets, machine learning might benefit these analyses. If the signals from the biosensors are optical, a smartphone can be used as a portable optical reader and potential data-analyzing device. Biosensing is a rich field that couples engineering and biology, and applying its toolset to help advance soil microbiome characterization would be a boon to microbiology more broadly.

Identification of distinct genotypes in circulating RSV A strains based on variants in the virus replication-associated genes

Respiratory syncytial virus (RSV) is a common cause of respiratory infection that often leads to hospitalization of infected younger children and older adults. RSV is classified into two strains, A and B, each with several subgroups or genotypes. One issue with the definition of these subgroups is the lack of a unified method of identification or genotyping. We propose that genotyping strategies based on the genes coding for replication-associated proteins could provide critical information on the replication capacity of the distinct subgroups, while clearly distinguishing genotypes. Here, we analyzed the virus replication-associated genes N, P, M2, and L from de novo assembled RSV A sequences obtained from 31 newly sequenced samples from hospitalized patients in Philadelphia and 78 additional publicly available sequences from different geographic locations within the United States. In-depth analysis and annotation of variants in the replication-associated proteins identified the polymerase protein L as a robust target for genotyping RSV subgroups. Importantly, our analysis revealed non-synonymous variations in L that were consistently accompanied by conserved changes in its co-factor P or the M2-2 protein, suggesting associations and interactions between specific domains of these proteins. Similar associations were seen among sequences of the related human metapneumovirus. These results highlight L as an alternative to other RSV genotyping targets and demonstrate the value of in-depth analyses and annotations of RSV sequences as it can serve as a foundation for subsequent in vitro and clinical studies on the efficiency of the polymerase and fitness of different virus isolates.IMPORTANCEGiven the historical heterogeneity of respiratory syncytial virus (RSV) and the disease it causes, there is a need to understand the properties of the circulating RSV strains each season. This information would benefit from an informative and consensus method of genotyping the virus. Here, we carried out a variant analysis that shows a pattern of specific variations among the replication-associated genes of RSV A across different seasons. Interestingly, these variation patterns, which were also seen in human metapneumovirus sequences, point to previously defined interactions of domains within these genes, suggesting co-variation in the replication-associated genes. Our results also suggest a genotyping strategy that can prove to be particularly important in understanding the genotype-phenotype correlation in the era of RSV vaccination, where selective pressure on the virus to evolve is anticipated. More importantly, the categorization of pneumoviruses based on these patterns may be of prognostic value.

Taking advantage of reference-guided assembly in a slowly-evolving lineage: Application to Testudo graeca

BACKGROUND

Obtaining de novo chromosome-level genome assemblies greatly enhances conservation and evolutionary biology studies. For many research teams, long-read sequencing technologies (that produce highly contiguous assemblies) remain unaffordable or unpractical. For the groups that display high synteny conservation, these limitations can be overcome by a reference-guided assembly using a close relative genome. Among chelonians, tortoises (Testudinidae) are considered one of the most endangered taxa, which calls for more genomic resources. Here we make the most of high synteny conservation in chelonians to produce the first chromosome-level genome assembly of the genus Testudo with one of the most iconic tortoise species in the Mediterranean basin: Testudo graeca.

RESULTS

We used high-quality, paired-end Illumina sequences to build a reference-guided assembly with the chromosome-level reference of Gopherus evgoodei. We reconstructed a 2.29 Gb haploid genome with a scaffold N50 of 107.598 Mb and 5.37% gaps. We sequenced 25,998 protein-coding genes, and identified 41.2% of the assembly as repeats. Demographic history reconstruction based on the genome revealed two events (population decline and recovery) that were consistent with previously suggested phylogeographic patterns for the species. This outlines the value of such reference-guided assemblies for phylogeographic studies.

CONCLUSIONS

Our results highlight the value of using close relatives to produce de novo draft assemblies in species where such resources are unavailable. Our annotated genome of T. graeca paves the way to delve deeper into the species' evolutionary history and provides a valuable resource to enhance direct conservation efforts on their threatened populations.

Reclassification of Botryococcus braunii chemical races into separate species based on a comparative genomics analysis

The colonial green microalga Botryococcus braunii is well known for producing liquid hydrocarbons that can be utilized as biofuel feedstocks. B. braunii is taxonomically classified as a single species made up of three chemical races, A, B, and L, that are mainly distinguished by the hydrocarbons produced. We previously reported a B race draft nuclear genome, and here we report the draft nuclear genomes for the A and L races. A comparative genomic study of the three B. braunii races and 14 other algal species within Chlorophyta revealed significant differences in the genomes of each race of B. braunii. Phylogenomically, there was a clear divergence of the three races with the A race diverging earlier than both the B and L races, and the B and L races diverging from a later common ancestor not shared by the A race. DNA repeat content analysis suggested the B race had more repeat content than the A or L races. Orthogroup analysis revealed the B. braunii races displayed more gene orthogroup diversity than three closely related Chlamydomonas species, with nearly 24-36% of all genes in each B. braunii race being specific to each race. This analysis suggests the three races are distinct species based on sufficient differences in their respective genomes. We propose reclassification of the three chemical races to the following species names: Botryococcus alkenealis (A race), Botryococcus braunii (B race), and Botryococcus lycopadienor (L race).

Genomes of the Orestias pupfish from the Andean Altiplano shed light on their evolutionary history and phylogenetic relationships within Cyprinodontiformes

BACKGROUND

To unravel the evolutionary history of a complex group, a comprehensive reconstruction of its phylogenetic relationships is crucial. This requires meticulous taxon sampling and careful consideration of multiple characters to ensure a complete and accurate reconstruction. The phylogenetic position of the Orestias genus has been estimated partly on unavailable or incomplete information. As a consequence, it was assigned to the family Cyprindontidae, relating this Andean fish to other geographically distant genera distributed in the Mediterranean, Middle East and North and Central America. In this study, using complete genome sequencing, we aim to clarify the phylogenetic position of Orestias within the Cyprinodontiformes order.

RESULTS

We sequenced the genome of three Orestias species from the Andean Altiplano. Our analysis revealed that the small genome size in this genus (~ 0.7 Gb) was caused by a contraction in transposable element (TE) content, particularly in DNA elements and short interspersed nuclear elements (SINEs). Using predicted gene sequences, we generated a phylogenetic tree of Cyprinodontiformes using 902 orthologs extracted from all 32 available genomes as well as three outgroup species. We complemented this analysis with a phylogenetic reconstruction and time calibration considering 12 molecular markers (eight nuclear and four mitochondrial genes) and a stratified taxon sampling to consider 198 species of nearly all families and genera of this order. Overall, our results show that phylogenetic closeness is directly related to geographical distance. Importantly, we found that Orestias is not part of the Cyprinodontidae family, and that it is more closely related to the South American fish fauna, being the Fluviphylacidae the closest sister group.

CONCLUSIONS

The evolutionary history of the Orestias genus is linked to the South American ichthyofauna and it should no longer be considered a member of the Cyprinodontidae family. Instead, we submit that Orestias belongs to the Orestiidae family, as suggested by Freyhof et al. (2017), and that it is the sister group of the Fluviphylacidae family, distributed in the Amazonian and Orinoco basins. These two groups likely diverged during the Late Eocene concomitant with hydrogeological changes in the South American landscape.

De novo transcriptome assembly of a lipoxygenase knock-down strain in the diatom Pseudo-nitzschia arenysensis

Diatoms are microalgae that live in marine and freshwater environments and are responsible for about 20% of the world's carbon fixation. Population dynamics of these cells is finely regulated by intricate signal transduction systems, in which oxylipins are thought to play a relevant role. These are oxygenated fatty acids whose biosynthesis is initiated by a lipoxygenase enzyme (LOX) and are widely distributed in all phyla, including diatoms. Here, we present a de novo transcriptome obtained from the RNA-seq performed in the diatom species Pseudo-nitzschia arenysensis, using both a wild-type and a LOX-silenced strain, which will represent a reliable reference for comparative analyses within the Pseudo-nitzschia genus and at a broader taxonomic scale. Moreover, the RNA-seq data can be interrogated to go deeper into the oxylipins metabolic pathways.

Identification of distinct genotypes in circulating RSV A strains based on variants on the virus replication-associated genes

Respiratory syncytial virus is a common cause of respiratory infection that often leads to hospitalization of infected younger children and older adults. RSV is classified into two strains, A and B, each with several subgroups or genotypes. One issue with the definition of these subgroups is the lack of a unified method of identification or genotyping. We propose that genotyping strategies based on the genes coding for replication-associated proteins could provide critical information on the replication capacity of the distinct subgroup, while clearly distinguishing genotypes. Here, we analyzed the virus replication-associated genes N, P, M2, and L from de novo assembled RSV A sequences obtained from 31 newly sequenced samples from hospitalized patients in Philadelphia and 78 additional publicly available sequences from different geographic locations within the US. In-depth analysis and annotation of the protein variants in L and the other replication-associated proteins N, P, M2-1, and M2-2 identified the polymerase protein L as a robust target for genotyping RSV subgroups. Importantly, our analysis revealed non-synonymous variations in L that were consistently accompanied by conserved changes in its co-factor P or the M2-2 protein, suggesting associations and interactions between specific domains of these proteins. These results highlight L as an alternative to other RSV genotyping targets and demonstrate the value of in-depth analyses and annotations of RSV sequences as it can serve as a foundation for subsequent in vitro and clinical studies on the efficiency of the polymerase and fitness of different virus isolates.

Haplotype-specific assembly of shattered chromosomes in esophageal adenocarcinomas

The epigenetic landscape of cancer is regulated by many factors, but primarily it derives from the underlying genome sequence. Chromothripsis is a catastrophic localized genome shattering event that drives, and often initiates, cancer evolution. We characterized five esophageal adenocarcinoma organoids with chromothripsis using long-read sequencing and transcriptome and epigenome profiling. Complex structural variation and subclonal variants meant that haplotype-aware de novo methods were required to generate contiguous cancer genome assemblies. Chromosomes were assembled separately and scaffolded using haplotype-resolved Hi-C reads, producing accurate assemblies even with up to 900 structural rearrangements. There were widespread differences between the chromothriptic and wild-type copies of chromosomes in topologically associated domains, chromatin accessibility, histone modifications, and gene expression. Differential epigenome peaks were most enriched within 10 kb of chromothriptic structural variants. Alterations in transcriptome and higher-order chromosome organization frequently occurred near differential epigenetic marks. Overall, chromothripsis reshapes gene regulation, causing coordinated changes in epigenetic landscape, transcription, and chromosome conformation.

MCPtaggR: R package for accurate genotype calling in reduced representation sequencing data by eliminating error-prone markers based on genome comparison

Reduced representation sequencing (RRS) offers cost-effective, high-throughput genotyping platforms such as genotyping-by-sequencing (GBS). RRS reads are typically mapped onto a reference genome. However, mapping reads harbouring mismatches against the reference can potentially result in mismapping and biased mapping, leading to the detection of error-prone markers that provide incorrect genotype information. We established a genotype-calling pipeline named mappable collinear polymorphic tag genotyping (MCPtagg) to achieve accurate genotyping by eliminating error-prone markers. MCPtagg was designed for the RRS-based genotyping of a population derived from a biparental cross. The MCPtagg pipeline filters out error-prone markers prior to genotype calling based on marker collinearity information obtained by comparing the genome sequences of the parents of a population to be genotyped. A performance evaluation on real GBS data from a rice F2 population confirmed its effectiveness. Furthermore, our performance test using a genome assembly that was obtained by genome sequence polishing on an available genome assembly suggests that our pipeline performs well with converted genomes, rather than necessitating de novo assembly. This demonstrates its flexibility and scalability. The R package, MCPtaggR, was developed to provide functions for the pipeline and is available at https://github.com/tomoyukif/MCPtaggR.

A novel computational pipeline for var gene expression augments the discovery of changes in the Plasmodium falciparum transcriptome during transition from in vivo to short-term in vitro culture

The pathogenesis of severe Plasmodium falciparum malaria involves cytoadhesive microvascular sequestration of infected erythrocytes, mediated by P. falciparum erythrocyte membrane protein 1 (PfEMP1). PfEMP1 variants are encoded by the highly polymorphic family of var genes, the sequences of which are largely unknown in clinical samples. Previously, we published new approaches for var gene profiling and classification of predicted binding phenotypes in clinical P. falciparum isolates (Wichers et al., 2021), which represented a major technical advance. Building on this, we report here a novel method for var gene assembly and multidimensional quantification from RNA-sequencing that outperforms the earlier approach of Wichers et al., 2021, on both laboratory and clinical isolates across a combination of metrics. Importantly, the tool can interrogate the var transcriptome in context with the rest of the transcriptome and can be applied to enhance our understanding of the role of var genes in malaria pathogenesis. We applied this new method to investigate changes in var gene expression through early transition of parasite isolates to in vitro culture, using paired sets of ex vivo samples from our previous study, cultured for up to three generations. In parallel, changes in non-polymorphic core gene expression were investigated. Modest but unpredictable var gene switching and convergence towards var2csa were observed in culture, along with differential expression of 19% of the core transcriptome between paired ex vivo and generation 1 samples. Our results cast doubt on the validity of the common practice of using short-term cultured parasites to make inferences about in vivo phenotype and behaviour.

The Utilization of Reference-Guided Assembly and In Silico Libraries Improves the Draft Genome of Clarias batrachus and Culter alburnus

Long-read sequencing technologies can generate highly contiguous genome assemblies compared to short-read methods. However, their higher cost often poses a significant barrier. To address this, we explore the utilization of mapping-based genome assembly and reference-guided assembly as cost-effective alternative approaches. We assess the efficacy of these approaches in improving the contiguity of Clarias batrachus and Culter alburnus draft genomes. Our findings demonstrate that employing an iterative mapping strategy leads to a reduction in assembly errors. Specifically, after three iterations, the Mismatches per 100 kbp value for the C. batrachus genome decreased from 2447.20 to 2432.67, reaching a minimum of 2422.67 after two iterations. Additionally, the N50 value for the C. batrachus genome increased from 362,143 to 1,315,126 bp, with a maximum of 1,315,403 bp after two iterations. Furthermore, we achieved Mismatches per 100 kbp values of 3.70 for the reference-guided assembly of C. batrachus and 0.34 for C. alburnus. Correspondingly, the N50 value for the C. batrachus and C. alburnus genomes increased from 362,143 bp and 3,686,385 bp to 2,026,888 bp and 43,735,735 bp, respectively. Finally, we successfully utilized the improved C. batrachus and C. alburnus genomes to compare genome studies using the combined approach of Ragout and Ragtag. Through a comprehensive comparative analysis of mapping-based and reference-guided genome assembly methods, we shed light on the specific contributions of reference-guided assembly in reducing assembly errors and improving assembly continuity and integrity. These advancements establish reference-guided assembly and the utilization of in silico libraries as a promising and suitable approach for comparative genomics studies.

A chromosome-level genome assembly of the Korean crossbred pig Nanchukmacdon (Sus scrofa)

As plentiful high-quality genome assemblies have been accumulated, reference-guided genome assembly can be a good approach to reconstruct a high-quality assembly. Here, we present a chromosome-level genome assembly of the Korean crossbred pig called Nanchukmacdon (the NCMD assembly) using the reference-guided assembly approach with short and long reads. The NCMD assembly contains 20 chromosome-level scaffolds with a total size of 2.38 Gbp (N50: 138.77 Mbp). Its BUSCO score is 93.1%, which is comparable to the pig reference assembly, and a total of 20,588 protein-coding genes, 8,651 non-coding genes, and 996.14 Mbp of repetitive elements are annotated. The NCMD assembly was also used to close many gaps in the pig reference assembly. This NCMD assembly and annotation provide foundational resources for the genomic analyses of pig and related species.

The Draft Genome of the "Golden Tide" Seaweed, Sargassum horneri: Characterization and Comparative Analysis

Sargassum horneri, a prevalent species of brown algae found along the coast of the northwest Pacific Ocean, holds significant importance as a valuable source of bioactive compounds. However, its rapid growth can lead to the formation of a destructive "golden tide", causing severe damage to the local economy and coastal ecosystems. In this study, we carried out de novo whole-genome sequencing of S. horneri using next-generation sequencing to unravel the genetic information of this alga. By utilizing a reference-guided de novo assembly pipeline with a closely related species, we successfully established a final assembled genome with a total length of 385 Mb. Repetitive sequences made up approximately 30.6% of this genome. Among the identified putative genes, around 87.03% showed homology with entries in the NCBI non-redundant protein database, with Ectocarpus siliculosus being the most closely related species for approximately one-third of these genes. One gene encoding an alkaline phosphatase family protein was found to exhibit positive selection, which could give a clue for the formation of S. horneri golden tides. Additionally, we characterized putative genes involved in fucoidan biosynthesis metabolism, a significant pathway in S. horneri. This study represents the first genome-wide characterization of a S. horneri species, providing crucial insights for future investigations, such as ecological genomic analyses.

A systematic comparison of human mitochondrial genome assembly tools

BACKGROUND

Mitochondria are the cell organelles that produce most of the chemical energy required to power the cell's biochemical reactions. Despite being a part of a eukaryotic host cell, the mitochondria contain a separate genome whose origin is linked with the endosymbiosis of a prokaryotic cell by the host cell and encode independent genomic information throughout their genomes. Mitochondrial genomes accommodate essential genes and are regularly utilized in biotechnology and phylogenetics. Various assemblers capable of generating complete mitochondrial genomes are being continuously developed. These tools often use whole-genome sequencing data as an input containing reads from the mitochondrial genome. Till now, no published work has explored the systematic comparison of all the available tools for assembling human mitochondrial genomes using short-read sequencing data. This evaluation is required to identify the best tool that can be well-optimized for small-scale projects or even national-level research.

RESULTS

In this study, we have tested the mitochondrial genome assemblers for both simulated datasets and whole genome sequencing (WGS) datasets of humans. For the highest computational setting of 16 computational threads with the simulated dataset having 1000X read depth, MitoFlex took the least execution time of 69 s, and IOGA took the longest execution time of 1278 s. NOVOPlasty utilized the least computational memory of approximately 0.098 GB for the same setting, whereas IOGA utilized the highest computational memory of 11.858 GB. In the case of WGS datasets for humans, GetOrganelle and MitoFlex performed the best in capturing the SNPs information with a mean F1-score of 0.919 at the sequencing depth of 10X. MToolBox and NOVOPlasty performed consistently across all sequencing depths with a mean F1 score of 0.897 and 0.890, respectively.

CONCLUSIONS

Based on the overall performance metrics and consistency in assembly quality for all sequencing data, MToolBox performed the best. However, NOVOPlasty was the second fastest tool in execution time despite being single-threaded, and it utilized the least computational resources among all the assemblers when tested on simulated datasets. Therefore, NOVOPlasty may be more practical when there is a significant sample size and a lack of computational resources. Besides, as long-read sequencing gains popularity, mitochondrial genome assemblers must be developed to use long-read sequencing data.

Evolutionary dynamics of genome size and content during the adaptive radiation of Heliconiini butterflies

Heliconius butterflies, a speciose genus of Müllerian mimics, represent a classic example of an adaptive radiation that includes a range of derived dietary, life history, physiological and neural traits. However, key lineages within the genus, and across the broader Heliconiini tribe, lack genomic resources, limiting our understanding of how adaptive and neutral processes shaped genome evolution during their radiation. Here, we generate highly contiguous genome assemblies for nine Heliconiini, 29 additional reference-assembled genomes, and improve 10 existing assemblies. Altogether, we provide a dataset of annotated genomes for a total of 63 species, including 58 species within the Heliconiini tribe. We use this extensive dataset to generate a robust and dated heliconiine phylogeny, describe major patterns of introgression, explore the evolution of genome architecture, and the genomic basis of key innovations in this enigmatic group, including an assessment of the evolution of putative regulatory regions at the Heliconius stem. Our work illustrates how the increased resolution provided by such dense genomic sampling improves our power to generate and test gene-phenotype hypotheses, and precisely characterize how genomes evolve.

Khairuddin RF, Mahmood Y. Genome assembly composition of the String "ACGT" array: a review of data structure accuracy and performance challenges

Background

The development of sequencing technology increases the number of genomes being sequenced. However, obtaining a quality genome sequence remains a challenge in genome assembly by assembling a massive number of short strings (reads) with the presence of repetitive sequences (repeats). Computer algorithms for genome assembly construct the entire genome from reads in two approaches. The de novo approach concatenates the reads based on the exact match between their suffix-prefix (overlapping). Reference-guided approach orders the reads based on their offsets in a well-known reference genome (reads alignment). The presence of repeats extends the technical ambiguity, making the algorithm unable to distinguish the reads resulting in misassembly and affecting the assembly approach accuracy. On the other hand, the massive number of reads causes a big assembly performance challenge.

Method

The repeat identification method was introduced for misassembly by prior identification of repetitive sequences, creating a repeat knowledge base to reduce ambiguity during the assembly process, thus enhancing the accuracy of the assembled genome. Also, hybridization between assembly approaches resulted in a lower misassembly degree with the aid of the reference genome. The assembly performance is optimized through data structure indexing and parallelization. This article's primary aim and contribution are to support the researchers through an extensive review to ease other researchers' search for genome assembly studies. The study also, highlighted the most recent developments and limitations in genome assembly accuracy and performance optimization.

Results

Our findings show the limitations of the repeat identification methods available, which only allow to detect of specific lengths of the repeat, and may not perform well when various types of repeats are present in a genome. We also found that most of the hybrid assembly approaches, either starting with de novo or reference-guided, have some limitations in handling repetitive sequences as it is more computationally costly and time intensive. Although the hybrid approach was found to outperform individual assembly approaches, optimizing its performance remains a challenge. Also, the usage of parallelization in overlapping and reads alignment for genome assembly is yet to be fully implemented in the hybrid assembly approach.

Conclusion

We suggest combining multiple repeat identification methods to enhance the accuracy of identifying the repeats as an initial step to the hybrid assembly approach and combining genome indexing with parallelization for better optimization of its performance.

Occurrences of similar viral diversity in campus wastewater and reclaimed water of a university dormitory

Water reuse from wastewater sources still remain some critical safety concerns associated with treacherous contaminants like pathogenic viruses. In this study, viral diversities in campus wastewater (CWW) and its reclaimed water (RCW) recycled for toilet flushing and garden irrigation of a university dormitory were assessed using metagenomic sequencing for acquisition of more background data. Results suggested majority (>80%) of gene sequences within assembled contigs predicted by open reading frame (ORF) finder were no-hit yet believed to be novel/unrevealed viral genomic information whereas hits matched bacteriophages (i.e., mainly Myoviridae, Podoviridae, and Siphoviridae families) were predominant in both CWW and RCW samples. Moreover, few pathogenic viruses (<1%) related to infections of human skin (e.g., Molluscum contagiosum virus, MCV), digestion system (e.g., hepatitis C virus, HCV), and gastrointestinal tract (e.g., human norovirus, HuNoV) were also noticed raising safety concerns about application of reclaimed waters. Low-affinity interactions of particular viral exterior proteins (e.g., envelope glycoproteins or spike proteins) for disinfectant ligand (e.g., chlorite) elucidated treatment limitations of current sewage processing systems even with membrane bioreactor and disinfectant contactor. Revolutionary disinfection approaches together with routine monitoring and new regulations are prerequisite to secure pathogen-correlated water quality for safer reuse of reclaimed waters.

GALA: a computational framework for de novo chromosome-by-chromosome assembly with long reads

High-quality genome assembly has wide applications in genetics and medical studies. However, it is still very challenging to achieve gap-free chromosome-scale assemblies using current workflows for long-read platforms. Here we report on GALA (Gap-free long-read Assembly tool), a computational framework for chromosome-based sequencing data separation and de novo assembly implemented through a multi-layer graph that identifies discordances within preliminary assemblies and partitions the data into chromosome-scale scaffolding groups. The subsequent independent assembly of each scaffolding group generates a gap-free assembly likely free from the mis-assembly errors which usually hamper existing workflows. This flexible framework also allows us to integrate data from various technologies, such as Hi-C, genetic maps, and even motif analyses to generate gap-free chromosome-scale assemblies. As a proof of principle we de novo assemble the C. elegans genome using combined PacBio and Nanopore sequencing data and a rice cultivar genome using Nanopore sequencing data from publicly available datasets. We also demonstrate the proposed method's applicability with a gap-free assembly of the human genome using PacBio high-fidelity (HiFi) long reads. Thus, our method enables straightforward assembly of genomes with multiple data sources and overcomes barriers that at present restrict the application of de novo genome assembly technology.

Alkaloid production and response to natural adverse conditions in Peganum harmala: in silico transcriptome analyses

Peganum harmala is a valuable wild plant that grows and survives under adverse conditions and produces pharmaceutical alkaloid metabolites. Using different assemblers to develop a transcriptome improves the quality of assembled transcriptome. In this study, a concrete and accurate method for detecting stress-responsive transcripts by comparing stress-related gene ontology (GO) terms and public domains was designed. An integrated transcriptome for P. harmala including 42 656 coding sequences was created by merging de novo assembled transcriptomes. Around 35 000 transcripts were annotated with more than 90% resemblance to three closely related species of Citrus, which confirmed the robustness of the assembled transcriptome; 4853 stress-responsive transcripts were identified. CYP82 involved in alkaloid biosynthesis showed a higher number of transcripts in P. harmala than in other plants, indicating its diverse alkaloid biosynthesis attributes. Transcription factors (TFs) and regulatory elements with 3887 transcripts comprised 9% of the transcriptome. Among the TFs of the integrated transcriptome, cystein2/histidine2 (C2H2) and WD40 repeat families were the most abundant. The Kyoto Encyclopedia of Genes and Genomes (KEGG) MAPK (mitogen-activated protein kinase) signaling map and the plant hormone signal transduction map showed the highest assigned genes to these pathways, suggesting their potential stress resistance. The P. harmala whole-transcriptome survey provides important resources and paves the way for functional and comparative genomic studies on this plant to discover stress-tolerance-related markers and response mechanisms in stress physiology, phytochemistry, ecology, biodiversity, and evolution. P. harmala can be a potential model for studying adverse environmental cues and metabolite biosynthesis and a major source for the production of various alkaloids.

Reference-based read clustering improves the de novo genome assembly of microbial strains

Constructing accurate microbial genome assemblies is necessary to understand genetic diversity in microbial genomes and its functional consequences. However, it still remains as a challenging task especially when only short-read sequencing technologies are used. Here, we present a new read-clustering algorithm, called RBRC, for improving de novo microbial genome assembly, by accurately estimating read proximity using multiple reference genomes. The performance of RBRC was confirmed by simulation-based evaluation in terms of assembly contiguity and the number of misassemblies, and was successfully applied to existing fungal and bacterial genomes by improving the quality of the assemblies without using additional sequencing data. RBRC is a very useful read-clustering algorithm that can be used (i) for generating high-quality genome assemblies of microbial strains when genome assemblies of related strains are available, and (ii) for upgrading existing microbial genome assemblies when the generation of additional sequencing data, such as long reads, is difficult.

Next Generation and Other Sequencing Technologies in Diagnostic Microbiology and Infectious Diseases

Next-generation sequencing (NGS) technologies have become increasingly available for use in the clinical microbiology diagnostic environment. There are three main applications of these technologies in the clinical microbiology laboratory: whole genome sequencing (WGS), targeted metagenomics sequencing and shotgun metagenomics sequencing. These applications are being utilized for initial identification of pathogenic organisms, the detection of antimicrobial resistance mechanisms and for epidemiologic tracking of organisms within and outside hospital systems. In this review, we analyze these three applications and provide a comprehensive summary of how these applications are currently being used in public health, basic research, and clinical microbiology laboratory environments. In the public health arena, WGS is being used to identify and epidemiologically track food borne outbreaks and disease surveillance. In clinical hospital systems, WGS is used to identify multi-drug-resistant nosocomial infections and track the transmission of these organisms. In addition, we examine how metagenomics sequencing approaches (targeted and shotgun) are being used to circumvent the traditional and biased microbiology culture methods to identify potential pathogens directly from specimens. We also expand on the important factors to consider when implementing these technologies, and what is possible for these technologies in infectious disease diagnosis in the next 5 years.

Whole-Genome Inter-Sex Variation in Russian Sturgeon (Acipenser gueldenstaedtii)

The Russian sturgeon (Acipenser gueldenstaedtii, AG) is an endangered fish species increasingly raised on fish farms for black caviar. Understanding the process of sex determination in AG is, therefore, of scientific and commercial importance. AG lacks sexual dimorphism until sexual maturation and has a predominantly octoploid genome without a definite sex chromosome. A conserved short female-specific genomic sequence was recently described, leading to the development of a genetic sex marker. However, no biological function has been reported for this sequence. Thus, the mechanism of sex determination and the overall inter-sex genomic variation in AG are still unknown. To comprehensively analyze the inter-sex genomic variation and assess the overall inter-species variation between AG and A. ruthenus (AR, sterlet), a related tetraploid sturgeon species, we performed whole-genome sequencing on DNA from 10 fish-farm-raised adult AG (5 males and 5 females). We produced a partially assembled, ~2390 MBp draft genome for AG. We validated in AG the female-specific region previously described in AR. We identified ~2.8 million loci (SNP/indels) varying between the species, but only ~7400 sex-associated loci in AG. We mapped the sex-associated AG loci to the AR genome and identified 15 peaks of sex-associated variation (10 kb segments with 30 or more sex-associated variants), 1 of which matched the previously reported sex-variable region. Finally, we identified 14 known and predicted genes in proximity to these peaks. Our analysis suggests that one or more of these genes may have functional roles in sex determination and/or sexual differentiation in sturgeons. Further functional studies are required to elucidate these roles.

Validation of reference-assisted assembly using existing and novel Heliothine genomes

Chloridea subflexa and Chloridea virescens are a pair of closely related noctuid species exhibiting pheromone-based sexual isolation and divergent host plant preferences. We produced a novel Illumina short read C. subflexa genome assembly and an improved C. virescens genome assembly, which offer opportunities to study the genomic basis for evolutionarily important traits in this lepidopteran family with few genomic resources. We then examined the feasibility of reference-assisted assembly, an approach that leverages existing high quality genomic resources for genome improvement in closely related taxa and applied it to our Heliothine genomes. Our work demonstrates that reference-assisted assembly has the potential to enhance contiguity and completeness of existing insect genomic resources with minimal additional laboratory costs. We conclude by discussing both the potential and pitfalls of reference-assisted assembly according to the intended downstream assembly application.

Reference-Guided De Novo Genome Assembly of the Flour Beetle Tribolium freemani

The flour beetle Tribolium freemani is a sibling species of the model organism and important pest Tribolium castaneum. The two species are so closely related that they can produce hybrid progeny, but the genetic basis of their differences has not been revealed. In this work, we sequenced the T. freemani genome by applying PacBio HiFi technology. Using the well-assembled T. castaneum genome as a reference, we assembled 262 Mb of the T. freemani genomic sequence and anchored it in 10 linkage groups corresponding to nine autosomes and sex chromosome X. The assembly showed 99.8% completeness of conserved insect genes, indicating a high-quality reference genome. Comparison with the T. castaneum assembly revealed that the main differences in genomic sequence between the two sibling species come from repetitive DNA, including interspersed and tandem repeats. In this work, we also provided the complete assembled mitochondrial genome of T. freemani. Although the genome assembly needs to be ameliorated in tandemly repeated regions, the first version of the T. freemani reference genome and the complete mitogenome presented here represent useful resources for comparative evolutionary studies of related species and for further basic and applied research on different biological aspects of economically important pests.

HGGA: hierarchical guided genome assembler

BACKGROUND

De novo genome assembly typically produces a set of contigs instead of the complete genome. Thus additional data such as genetic linkage maps, optical maps, or Hi-C data is needed to resolve the complete structure of the genome. Most of the previous work uses the additional data to order and orient contigs.

RESULTS

Here we introduce a framework to guide genome assembly with additional data. Our approach is based on clustering the reads, such that each read in each cluster originates from nearby positions in the genome according to the additional data. These sets are then assembled independently and the resulting contigs are further assembled in a hierarchical manner. We implemented our approach for genetic linkage maps in a tool called HGGA.

CONCLUSIONS

Our experiments on simulated and real Pacific Biosciences long reads and genetic linkage maps show that HGGA produces a more contiguous assembly with less contigs and from 1.2 to 9.8 times higher NGA50 or N50 than a plain assembly of the reads and 1.03 to 6.5 times higher NGA50 or N50 than a previous approach integrating genetic linkage maps with contig assembly. Furthermore, also the correctness of the assembly remains similar or improves as compared to an assembly using only the read data.

Fast, Ungapped Reads Mapping Using Squid

Advances in Next Generation Sequencing technologies allow us to inspect and unlock the genome to a level of detail that was unimaginable only a few decades ago. Omics-based studies are casting a light on the patterns and determinants of disease conditions in populations, as well as on the influence of microbial communities on human health, just to name a few. Through increasing volumes of sequencing information, for example, it is possible to compare genomic features and analyze the modulation of the transcriptome under different environmental stimuli. Although protocols for NGS preparation are intended to leave little to no space for contamination of any kind, a noticeable fraction of sequencing reads still may not uniquely represent what was intended to be sequenced in the first place. If a natural consequence of a sequencing sample is to assess the presence of features of interest by mapping the obtained reads to a genome of reference, sometimes it is useful to determine the fraction of those that do not map, or that map discordantly, and store this information to a new file for subsequent analyses. Here we propose a new mapper, which we called Squid, that among other accessory functionalities finds and returns sequencing reads that match or do not match to a reference sequence database in any orientation. We encourage the use of Squid prior to any quantification pipeline to assess, for instance, the presence of contaminants, especially in RNA-Seq experiments.

A Single Amino Acid Substitution in Elongation Factor G Can Confer Low-Level Gentamicin Resistance in Neisseria gonorrhoeae

The continued emergence of Neisseria gonorrhoeae isolates which are resistant to first-line antibiotics has reinvigorated interest in alternative therapies such as expanded use of gentamicin (Gen). We hypothesized that expanded use of Gen promotes emergence of gonococci with clinical resistance to this aminoglycoside. To understand how decreased susceptibility of gonococci to Gen might develop, we selected spontaneous low-level Gen-resistant (Gen^R) mutants (Gen MIC = 32 μg/mL) of the Gen-susceptible strain FA19. Consequently, we identified a novel missense mutation in fusA, which encodes elongation factor G (EF-G), causing an alanine (A) to valine (V) substitution at amino acid position 563 in domain IV of EF-G; the mutant allele was termed fusA2. Transformation analysis showed that fusA2 could increase the Gen MIC by 4-fold. While possession of fusA2 did not impair either in vitro gonococcal growth or protein synthesis, it did result in a fitness defect during experimental infection of the lower genital tract in female mice. Through bioinformatic analysis of whole-genome sequences of 10,634 international gonococcal clinical isolates, other fusA alleles were frequently detected, but genetic studies revealed that they could not decrease Gen susceptibility in a similar manner to fusA2. In contrast to these diverse international fusA alleles, the fusA2-encoded A563V substitution was detected in only a single gonococcal clinical isolate. We hypothesize that the rare occurrence of fusA2 in N. gonorrhoeae clinical isolates is likely due to a fitness cost during infection, but compensatory mutations which alleviate this fitness cost could emerge and promote Gen^R in global strains.

THE CARBON FOOTPRINT OF BIOINFORMATICS

Bioinformatic research relies on large-scale computational infrastructures which have a nonzero carbon footprint but so far, no study has quantified the environmental costs of bioinformatic tools and commonly run analyses. In this work, we estimate the carbon footprint of bioinformatics (in kilograms of CO₂ equivalent units, kgCO₂e) using the freely available Green Algorithms calculator (www.green-algorithms.org, last accessed 2022). We assessed 1) bioinformatic approaches in genome-wide association studies (GWAS), RNA sequencing, genome assembly, metagenomics, phylogenetics, and molecular simulations, as well as 2) computation strategies, such as parallelization, CPU (central processing unit) versus GPU (graphics processing unit), cloud versus local computing infrastructure, and geography. In particular, we found that biobank-scale GWAS emitted substantial kgCO₂e and simple software upgrades could make it greener, for example, upgrading from BOLT-LMM v1 to v2.3 reduced carbon footprint by 73%. Moreover, switching from the average data center to a more efficient one can reduce carbon footprint by approximately 34%. Memory over-allocation can also be a substantial contributor to an algorithm’s greenhouse gas emissions. The use of faster processors or greater parallelization reduces running time but can lead to greater carbon footprint. Finally, we provide guidance on how researchers can reduce power consumption and minimize kgCO₂e. Overall, this work elucidates the carbon footprint of common analyses in bioinformatics and provides solutions which empower a move toward greener research.

Comparative Genomics of Three Colletotrichum scovillei Strains and Genetic Analysis Revealed Genes Involved in Fungal Growth and Virulence on Chili Pepper

Colletotrichum scovillei causes anthracnose of chili pepper in many countries. Three strains of this pathogen, Coll-524, Coll-153, and Coll-365, show varied virulence on chili pepper. Among the three strains, Coll-365 showed significant defects in growth and virulence. To decipher the genetic variations among these strains and identify genes contributing to growth and virulence, comparative genomic analysis and gene transformation to show gene function were applied in this study. Compared to Coll-524, Coll-153, and Coll-365 had numerous gene losses including 32 candidate effector genes that are mainly exist in acutatum species complex. A cluster of 14 genes in a 34-kb genomic fragment was lost in Coll-365. Through gene transformation, three genes in the 34-kb fragment were identified to have functions in growth and/or virulence of C. scovillei. CsPLAA encoding a phospholipase A2-activating protein enhanced the growth of Coll-365. A combination of CsPLAA with one transcription factor CsBZTF and one C6 zinc finger domain-containing protein CsCZCP was found to enhance the pathogenicity of Coll-365. Introduction of CsGIP, which encodes a hypothetical protein, into Coll-365 caused a reduction in the germination rate of Coll-365. In conclusion, the highest virulent strain Coll-524 had more genes and encoded more pathogenicity related proteins and transposable elements than the other two strains, which may contribute to the high virulence of Coll-524. In addition, the absence of the 34-kb fragment plays a critical role in the defects of growth and virulence of strain Coll-365.

Sequencing and Reconstructing Helminth Mitochondrial Genomes Directly from Genomic Next-Generation Sequencing Data

We present a detailed method for extraction of high-molecular weight genomic DNA suitable for numerous DNA sequencing applications, and a straightforward in silico approach for reconstructing novel mitochondrial (mt) genomes directly from total genomic DNA extracts derived from next-generation sequencing (NGS) data sets. The in silico post-sequencing pipeline described is fast, accurate, and highly efficient, with modest memory requirements that can be performed using a standard desktop computer. The approach is particularly effective for obtaining mitochondrial genomes for species with little or no mitochondrial sequence information currently available and overcomes many of the limitations of traditional strategies. The described methodologies are also applicable for metagenomics sequencing from mixed or pooled samples containing multiple species and subsequent specific assembly of specific mitochondrial genomes.

Rapid screening and identification of viral pathogens in metagenomic data

BACKGROUND

Virus screening and viral genome reconstruction are urgent and crucial for the rapid identification of viral pathogens, i.e., tracing the source and understanding the pathogenesis when a viral outbreak occurs. Next-generation sequencing (NGS) provides an efficient and unbiased way to identify viral pathogens in host-associated and environmental samples without prior knowledge. Despite the availability of software, data analysis still requires human operations. A mature pipeline is urgently needed when thousands of viral pathogen and viral genome reconstruction samples need to be rapidly identified.

RESULTS

In this paper, we present a rapid and accurate workflow to screen metagenomics sequencing data for viral pathogens and other compositions, as well as enable a reference-based assembler to reconstruct viral genomes. Moreover, we tested our workflow on several metagenomics datasets, including a SARS-CoV-2 patient sample with NGS data, pangolins tissues with NGS data, Middle East Respiratory Syndrome (MERS)-infected cells with NGS data, etc. Our workflow demonstrated high accuracy and efficiency when identifying target viruses from large scale NGS metagenomics data. Our workflow was flexible when working with a broad range of NGS datasets from small (kb) to large (100 Gb). This took from a few minutes to a few hours to complete each task. At the same time, our workflow automatically generates reports that incorporate visualized feedback (e.g., metagenomics data quality statistics, host and viral sequence compositions, details about each of the identified viral pathogens and their coverages, and reassembled viral pathogen sequences based on their closest references).

CONCLUSIONS

Overall, our system enabled the rapid screening and identification of viral pathogens from metagenomics data, providing an important piece to support viral pathogen research during a pandemic. The visualized report contains information from raw sequence quality to a reconstructed viral sequence, which allows non-professional people to screen their samples for viruses by themselves (Additional file 1).

Comparative Genomics of Clinical and Environmental Isolates of Vibrio spp. of Colombia: Implications of Traits Associated with Virulence and Resistance

There is widespread concern about the increase in cases of human and animal infections caused by pathogenic Vibrio species due to the emergence of epidemic lineages. In Colombia, active surveillance by the National Institute of Health (INS) has confirmed the presence of Vibrio; however, in routine surveillance, these isolates are not genomically characterized. This study focused on the pangenome analysis of six Vibrio species: V. parahaemolyticus, V. vulnificus, V. alginolyticus, V. fluvialis, V. diabolicus and V. furnissii to determine the genetic architectures of potentially virulent and antimicrobial resistance traits. Isolates from environmental and clinical samples were genome sequenced, assembled and annotated. The most important species in public health were further characterized by multilocus sequence typing and phylogenomics. For V. parahaemolyticus, we found the virulent ST3 and ST120 genotypes. For V. vulnificus, we identified isolates belonging to lineages 1 and 2. Virulence gene homologues between species were found even in non-pathogenic species such as V. diabolicus. Annotations related to the mobilome, integrative mobile and conjugative elements and resistance genes were obtained from environmental and clinical isolates. This study contributes genomic information to the intensified surveillance program implemented by the INS to establish potential sources of vibriosis in Colombia.

Biosynthesis of Ditropolonyl Sulfide, an Antibacterial Compound Produced by Burkholderia cepacia Complex Strain R-12632

Burkholderia cepacia complex strain R-12632 produces ditropolonyl sulfide, an unusual sulfur-containing tropone, via a yet-unknown biosynthetic pathway. Ditropolonyl sulfide purified from a culture of strain R-12632 inhibits the growth of various Gram-positive and Gram-negative resistant bacteria, with MIC values as low as 16 μg/ml. In the present study, we used a transposon mutagenesis approach combined with metabolite analyses to identify the genetic basis for antibacterial activity of strain R-12632 against Gram-negative bacterial pathogens. Fifteen of the 8304 transposon mutants investigated completely lost antibacterial activity against Klebsiella pneumoniae LMG 2095. In these loss-of-activity mutants, nine genes were interrupted. Four of those genes were involved in assimilatory sulfate reduction, two were involved in phenylacetic acid (PAA) catabolism, and one was involved in glutathione metabolism. Via semipreparative fractionation and metabolite identification, it was confirmed that inactivation of the PAA degradation pathway or glutathione metabolism led to loss of ditropolonyl sulfide production. Based on earlier studies on the biosynthesis of tropolone compounds, the requirement for a functional PAA catabolic pathway for antibacterial activity in strain R-12632 indicated that this pathway likely provides the tropolone backbone for ditropolonyl sulfide. Loss of activity observed in mutants defective in assimilatory sulfate reduction and glutathione biosynthesis suggested that cysteine and glutathione are potential sources of the sulfur atom linking the two tropolone moieties. The demonstrated antibacterial activity of the unusual antibacterial compound ditropolonyl sulfide warrants further studies into its biosynthesis and biological role. IMPORTANCE Burkholderia bacteria are historically known for their biocontrol properties and have been proposed as a promising and underexplored source of bioactive specialized metabolites. Burkholderia cepacia complex strain R-12632 inhibits various Gram-positive and Gram-negative resistant pathogens and produces numerous specialized metabolites, among which is ditropolonyl sulfide. This unusual antimicrobial has been poorly studied and its biosynthetic pathway remains unknown. In the present study, we performed transposon mutagenesis of strain R-12632 and performed genome and metabolite analyses of loss-of-activity mutants to study the genetic basis for antibacterial activity. Our results indicate that phenylacetic acid catabolism, assimilatory sulfate reduction, and glutathione metabolism are necessary for ditropolonyl sulfide production. These findings contribute to understanding of the biosynthesis and biological role of this unusual antimicrobial.

Whole-genome assembly and annotation of northern wild rice, Zizania palustris L., supports a whole-genome duplication in the Zizania genus

Zizania palustris L. (northern wild rice, NWR) is an aquatic grass native to North America that is notable for its nutritious grain. This is an important species with ecological, cultural and agricultural significance, specifically in the Great Lakes region of the USA. Using flow cytometry, we first estimated the NWR genome size to be 1.8 Gb. Using long- and short-range sequencing, Hi-C scaffolding and RNA-seq data from eight tissues, we generated an annotated whole-genome de novo assembly of NWR. The assembly was 1.29 Gb in length, highly repetitive (approx. 76.0%) and contained 46 421 putative protein-coding genes. The expansion of retrotransposons within the genome and a whole-genome duplication (WGD) after the Zizania-Oryza speciation event have both led to an increase in the genome size of NWR in comparison with Oryza sativa L. and Zizania latifolia. Both events depict a genome rapidly undergoing change over a short evolutionary time. Comparative analyses revealed the conservation of large syntenic blocks between NWR and O. sativa, which were used to identify putative seed-shattering genes. Estimates of divergence times revealed that the Zizania genus diverged from Oryza approximately 26-30 million years ago (26-30 MYA), whereas NWR and Z. latifolia diverged from one another approximately 6-8 MYA. Comparative genomics confirmed evidence of a WGD in the Zizania genus and provided support that the event occurred prior to the NWR-Z. latifolia speciation event. This genome assembly and annotation provides a valuable resource for comparative genomics in the Oryzeae tribe and provides an important resource for future conservation and breeding efforts of NWR.

PromethION sequencing and assembly of the genome of Micropoecilia picta, a fish with a highly Degenerated Y chromosome

We here describe sequencing and assembly of both the autosomes and the sex chromosome in M. picta, the closest related species to the guppy, Poecilia reticulata. Poecilia ()Micropoecilia) picta is a close outgroup for studying the guppy, an important organism for studies in evolutionary ecology and in sex chromosome evolution. The guppy XY pair (LG12) has long been studied as a test case for the importance of sexually antagonistic variants in selection for suppressed recombination between Y and X chromosomes. The guppy Y chromosome is not degenerated, but appears to carry functional copies of all genes that are present on its X counterpart. The X chromosomes of M. picta (and its relative M. parae) are homologous to the guppy XY pair, but their Y chromosomes are highly degenerated, and no genes can be identified in the fully Y-linked region. A complete genome sequence of a M. picta male may therefore contribute to understanding how the guppy Y evolved. These fish species' genomes are estimated to be about 750 Mb, with high densities of repetitive sequences, suggesting that long-read sequencing is needed. We evaluated several assembly approaches, and used our results to investigate the extent of Y chromosome degeneration in this species.

Exposure to dexamethasone modifies transcriptomic responses of free-living stages of Strongyloides stercoralis

Hyperinfection and disseminated infection by the parasitic nematode Strongyloides stercoralis can be induced by iatrogenic administration of steroids and immunosuppression and lead to an elevated risk of mortality. Responses of free-living stages of S. stercoralis to the therapeutic corticosteroid dexamethasone (DXM) were investigated using RNA-seq transcriptomes of DXM-treated female and male worms. A total of 17,950 genes representing the transcriptome of these free-living adult stages were obtained, among which 199 and 263 were differentially expressed between DXM-treated females and DXM-treated males, respectively, compared with controls. According to Gene Ontology analysis, differentially expressed genes from DXM-treated females participate in developmental process, multicellular organismal process, cell differentiation, carbohydrate metabolic process and embryonic morphogenesis. Others are involved in signaling and signal transduction, including cAMP, cGMP-dependent protein kinase pathway, endocrine system, and thyroid hormone pathway, as based on Kyoto Encyclopedia of Genes and Genomes analysis. The novel findings warrant deeper investigation of the influence of DXM on growth and other pathways in this neglected tropical disease pathogen, particularly in a setting of autoimmune and/or allergic disease, which may require the clinical use of steroid-like hormones during latent or covert strongyloidiasis.

Underpinning the molecular programming attributing heat stress associated thermotolerance in tea (Camellia sinensis (L.) O. Kuntze)

The most daunting issue of global climate change is the deleterious impact of extreme temperatures on tea productivity and quality, which has resulted in a quest among researchers and growers. The current study aims to unravel molecular programming underpinning thermotolerance by characterizing heat tolerance and sensitivity response in 20 tea cultivars. The significantly higher negative influence of heat stress was recorded in a sensitive cultivar with reduced water retention (47%), chlorophyll content (33.79%), oxidation potential (32.48%), and increase in membrane damage (76.4%). Transcriptional profiling of most tolerant and sensitive cultivars identified 78 differentially expressed unigenes with chaperon domains, including low and high molecular weight heat shock protein (HSP) and heat shock transcription factors (HSFs) involved in heat shock response (HSR). Further, predicted transcriptional interactome network revealed their key role in thermotolerance via well-co-ordinated transcriptional regulation of aquaporins, starch metabolism, chlorophyll biosynthesis, calcium, and ethylene mediated plant signaling system. The study identified the key role of HSPs (CsHSP90) in regulating HSR in tea, wherein, structure-based molecular docking revealed the inhibitory role of geldanamycin (GDA) on CsHSP90 by blocking ATP binding site at N-terminal domain of predicted structure. Subsequently, GDA mediated leaf disc inhibitor assay further affirmed enhanced HSR with higher expression of CsHSP17.6, CsHSP70, HSP101, and CsHSFA2 genes in tea. Through the current study, efforts were made to extrapolate a deeper understanding of chaperons mediated regulation of HSR attributing thermotolerance in tea.

LoReTTA, a user-friendly tool for assembling viral genomes from PacBio sequence data

Long-read, single-molecule DNA sequencing technologies have triggered a revolution in genomics by enabling the determination of large, reference-quality genomes in ways that overcome some of the limitations of short-read sequencing. However, the greater length and higher error rate of the reads generated on long-read platforms make the tools used for assembling short reads unsuitable for use in data assembly and motivate the development of new approaches. We present LoReTTA (Long Read Template-Targeted Assembler), a tool designed for performing de novo assembly of long reads generated from viral genomes on the PacBio platform. LoReTTA exploits a reference genome to guide the assembly process, an approach that has been successful with short reads. The tool was designed to deal with reads originating from viral genomes, which feature high genetic variability, possible multiple isoforms, and the dominant presence of additional organisms in clinical or environmental samples. LoReTTA was tested on a range of simulated and experimental datasets and outperformed established long-read assemblers in terms of assembly contiguity and accuracy. The software runs under the Linux operating system, is designed for easy adaptation to alternative systems, and features an automatic installation pipeline that takes care of the required dependencies. A command-line version and a user-friendly graphical interface version are available under a GPLv3 license at https://bioinformatics.cvr.ac.uk/software/ with the manual and a test dataset.

Comparative Genomics Supports That Brazilian Bioethanol Saccharomyces cerevisiae Comprise a Unified Group of Domesticated Strains Related to Cachaça Spirit Yeasts

Ethanol production from sugarcane is a key renewable fuel industry in Brazil. Major drivers of this alcoholic fermentation are Saccharomyces cerevisiae strains that originally were contaminants to the system and yet prevail in the industrial process. Here we present newly sequenced genomes (using Illumina short-read and PacBio long-read data) of two monosporic isolates (H3 and H4) of the S. cerevisiae PE-2, a predominant bioethanol strain in Brazil. The assembled genomes of H3 and H4, together with 42 draft genomes of sugarcane-fermenting (fuel ethanol plus cachaça) strains, were compared against those of the reference S288C and diverse S. cerevisiae. All genomes of bioethanol yeasts have amplified SNO2(3)/SNZ2(3) gene clusters for vitamin B1/B6 biosynthesis, and display ubiquitous presence of a particular family of SAM-dependent methyl transferases, rare in S. cerevisiae. Widespread amplifications of quinone oxidoreductases YCR102C/YLR460C/YNL134C, and the structural or punctual variations among aquaporins and components of the iron homeostasis system, likely represent adaptations to industrial fermentation. Interesting is the pervasive presence among the bioethanol/cachaça strains of a five-gene cluster (Region B) that is a known phylogenetic signature of European wine yeasts. Combining genomes of H3, H4, and 195 yeast strains, we comprehensively assessed whole-genome phylogeny of these taxa using an alignment-free approach. The 197-genome phylogeny substantiates that bioethanol yeasts are monophyletic and closely related to the cachaça and wine strains. Our results support the hypothesis that biofuel-producing yeasts in Brazil may have been co-opted from a pool of yeasts that were pre-adapted to alcoholic fermentation of sugarcane for the distillation of cachaça spirit, which historically is a much older industry than the large-scale fuel ethanol production.

Fine-scale empirical data on niche divergence and homeolog expression patterns in an allopolyploid and its diploid progenitor species

Polyploidization is pervasive in plants, but little is known about the niche divergence of wild allopolyploids (species that harbor polyploid genomes originating from different diploid species) relative to their diploid progenitor species and the gene expression patterns that may underlie such ecological divergence. We conducted a fine-scale empirical study on habitat and gene expression of an allopolyploid and its diploid progenitors. We quantified soil properties and light availability of habitats of an allotetraploid Cardamine flexuosa and its diploid progenitors Cardamine amara and Cardamine hirsuta in two seasons. We analyzed expression patterns of genes and homeologs (homeologous gene copies in allopolyploids) using RNA sequencing. We detected niche divergence between the allopolyploid and its diploid progenitors along water availability gradient at a fine scale: the diploids in opposite extremes and the allopolyploid in a broader range between diploids, with limited overlap with diploids at both ends. Most of the genes whose homeolog expression ratio changed among habitats in C. flexuosa varied spatially and temporally. These findings provide empirical evidence for niche divergence between an allopolyploid and its diploid progenitor species at a fine scale and suggest that divergent expression patterns of homeologs in an allopolyploid may underlie its persistence in diverse habitats.

High-Level Antibiotic Tolerance of a Clinically Isolated Enterococcus faecalis Strain

Bacteria can survive antibiotic treatment both by acquiring antibiotic resistance genes and through mechanisms of tolerance that are based on phenotypic changes and the formation of metabolically inactive cells. Here, we report an Enterococcus faecalis strain (E. faecalis UM001B) that was isolated from a cystic fibrosis patient and had no increase in resistance but extremely high-level tolerance to ampicillin, vancomycin, and tetracycline. Specifically, the percentages of cells that survived 3.5-h antibiotic treatment (at 100 μg · ml^-1) were 25.4% ± 4.3% and 51.9% ± 4.0% for ampicillin and tetracycline, respectively; vancomycin did not exhibit any significant killing. Consistent with the changes in antibiotic susceptibility, UM001B was found to have reduced penetration of ampicillin and vancomycin and accumulation of tetracycline compared to the reference strain ATCC 29212. Based on whole-genome sequencing, four amino acid substitutions were identified in one of the tetracycline efflux pump repressors (TetRs), compared to ATCC 29212. Results of molecular simulations and experimental assays revealed that these mutations could lead to higher levels of tetracycline efflux activity. Consistently, replicating these mutations in Escherichia coli MG1655 increased its tolerance to tetracycline. Overall, these findings provide new insights into the development of multidrug tolerance in E. faecalis, which can facilitate future studies to better control enterococcal infections.IMPORTANCE Enterococcus faecalis represents a major group of pathogens causing nosocomial infections that are resistant to multiple classes of antibiotics. An important challenge associated with E. faecalis infection is the emergence of multidrug-tolerant strains, which have normal MICs but do not respond to antibiotic treatment. Here, we report a strain of E. faecalis that was isolated from a cystic fibrosis patient and demonstrated high-level tolerance to ampicillin, vancomycin, and tetracycline. Whole-genome sequencing revealed critical substitutions in one of the tetracycline efflux pump repressors that are consistent with the increased tolerance of E. faecalis UM001B to tetracycline. These findings provide new information about bacterial antibiotic tolerance and may help develop more effective therapeutics.

The Allelic Diversity of the Gibberellin Signaling Pathway Genes in Aegilops tauschii Coss

Gibberellin-insensitive reduced height genes are widely spread in modern wheat varieties, making them resistant to lodging under conditions of intensive farming. However, the limited diversity of these genes present in wheat germplasm can limit the adaptability of newly created cultivars to the changing climate. The diversity of the gibberellin signaling pathway genes involved in plant height control- Reduced height 1 (Rht-D1), Gibberellin-insensitive dwarf 1 (Gid1‑D) and Gibberellin-insensitive dwarf 2 (Gid2-D)-was studied in the diploid wild goatgrass Aegilops tauschii Coss., one of the ancestral species of the bread wheat (Triticum aestivum L.) and the donor of its D subgenome, using high-throughput sequencing. The examination of 24 Ae. tauschii accessions of different geographical origins revealed a large number of new alleles (haplotypes) not found in bread wheat varieties. Some of the detected polymorphisms lead to changes in the amino acid sequence of proteins. Four isoforms (amino acid sequence variants) were found for the RHT-D1 protein, and two isoforms-for the GID1 and GID2 proteins, each. An analysis of the co-occurrence frequencies of various isoforms of the three proteins showed that their combinations were not random in Ae. tauschii, which may indicate the functional significance of their differences. New alleles of the Rht-D1, Gid1-D, and Gid2-D genes are promising for introgression into bread wheat and studying their effect on plant height and adaptability.

Reference genome and demographic history of the most endangered marine mammal, the vaquita

The vaquita is the most critically endangered marine mammal, with fewer than 19 remaining in the wild. First described in 1958, the vaquita has been in rapid decline for more than 20 years resulting from inadvertent deaths due to the increasing use of large-mesh gillnets. To understand the evolutionary and demographic history of the vaquita, we used combined long-read sequencing and long-range scaffolding methods with long- and short-read RNA sequencing to generate a near error-free annotated reference genome assembly from cell lines derived from a female individual. The genome assembly consists of 99.92% of the assembled sequence contained in 21 nearly gapless chromosome-length autosome scaffolds and the X-chromosome scaffold, with a scaffold N50 of 115 Mb. Genome-wide heterozygosity is the lowest (0.01%) of any mammalian species analysed to date, but heterozygosity is evenly distributed across the chromosomes, consistent with long-term small population size at genetic equilibrium, rather than low diversity resulting from a recent population bottleneck or inbreeding. Historical demography of the vaquita indicates long-term population stability at less than 5,000 (Ne) for over 200,000 years. Together, these analyses indicate that the vaquita genome has had ample opportunity to purge highly deleterious alleles and potentially maintain diversity necessary for population health.

Twelve quick steps for genome assembly and annotation in the classroom

Eukaryotic genome sequencing and de novo assembly, once the exclusive domain of well-funded international consortia, have become increasingly affordable, thus fitting the budgets of individual research groups. Third-generation long-read DNA sequencing technologies are increasingly used, providing extensive genomic toolkits that were once reserved for a few select model organisms. Generating high-quality genome assemblies and annotations for many aquatic species still presents significant challenges due to their large genome sizes, complexity, and high chromosome numbers. Indeed, selecting the most appropriate sequencing and software platforms and annotation pipelines for a new genome project can be daunting because tools often only work in limited contexts. In genomics, generating a high-quality genome assembly/annotation has become an indispensable tool for better understanding the biology of any species. Herein, we state 12 steps to help researchers get started in genome projects by presenting guidelines that are broadly applicable (to any species), sustainable over time, and cover all aspects of genome assembly and annotation projects from start to finish. We review some commonly used approaches, including practical methods to extract high-quality DNA and choices for the best sequencing platforms and library preparations. In addition, we discuss the range of potential bioinformatics pipelines, including structural and functional annotations (e.g., transposable elements and repetitive sequences). This paper also includes information on how to build a wide community for a genome project, the importance of data management, and how to make the data and results Findable, Accessible, Interoperable, and Reusable (FAIR) by submitting them to a public repository and sharing them with the research community.