Icarus: visualizer for de novo assembly evaluation

Bacterial nitrite production oxidizes Fe(II) bioremediating acidic abandoned coal mine drainage

Passive remediation systems (PRSs) treating either acidic or neutral abandoned coal mine drainage (AMD) are colonized by bacteria that can bioremediate iron (Fe) through chemical cycling. Due to the low pH in acidic AMD, iron oxidation from soluble Fe(II) to precipitated Fe(III) is mainly directed by microbial oxidation. Less well described are biotic reactions that lead to iron remediation through abiotic secondary reactions. We describe here iron oxidation in acidic AMD that is mediated by the bacterial reduction of nitrate to nitrite followed by the geochemical oxidation of Fe(II). Within an acidic PRS, 4,560 bacteria cultured from the microbial community were screened for their ability to oxidize iron and to perform nitrate-dependent iron oxidation (NDFO). Iron oxidation in the culturable community was observed in every pond of the system, ranging from 2.1% to 11.4%, and NDFO was observed in every pond, ranging from 1.4% to 6.0% of the culturable bacteria. Five NDFO isolates were purified and identified as Paraburkholderia spp. One of our isolates, Paraburkholderia sp. AV18 was shown to drive NDFO through the bacterial production of nitrite that in turn chemically oxidizes Fe(II) (nitrate reduction-iron oxidation; NRIO). AV18 expressed nitrate reductase, napA, concurrent to nitrite production. Burkholderiales are found by 16S rRNA gene sequencing in every pond of the PRS. The frequency of NDFO metabolism in the culturable microbial community and abundance of Burkholderiales in the PRS suggest nitrite producers contribute to the bioremediation of iron in acidic AMD and may be an unharnessed opportunity to increase iron bioremediation in acidic conditions.

IMPORTANCE

Our study sheds light on a poorly defined biogeochemical interaction, nitrate-dependent iron oxidation (NDFO), that has been described in several environments. We show that bacterial nitrate reduction produces nitrite, which can chemically oxidize ferrous iron, leading to insoluble ferric iron. We show that bacteria capable of the nitrate reduction-iron oxidation (NRIO) reactions are prevalent throughout multiple passive remediation systems that treat acidic coal mine drainage, indicating this may be a widespread mechanism for iron removal under acidic conditions. In acidic coal mine remediation, iron precipitation has been shown to be solely bacterially mediated, and NRIO provides a simple mechanism for aerobic oxidation of iron in these conditions.

Complete genome sequence of Planococcus koreensis isolated from soil in Fort Collins, Colorado

The complete genome of Planococcus koreensis was obtained using Nanopore MinION sequencing after isolation from soil in Colorado. The assembled genome contains one circular contig with 3,519,105 bp, 3,606 genes, 419 pseudogenes, and 47.62% guanine-cytosine content. This discovery provides a fully assembled P. koreensis genome available at the National Center for Biotechnology Information.

Draft genome sequences of 58 isolates of Campylobacter jejuni and Campylobacter coli from domestic ducks

We present the draft genomic sequences of 58 isolates of Campylobacter spp. (32 Campylobacter coli and 26 Campylobacter jejuni) isolated from domestic ducks (Anas platyrhynchos domesticus). This provides a substantial increase in the number of publicly available genomic sequences of Campylobacter isolated from farm-raised ducks in the United States.

Three new Pseudogymnoascus species (Pseudeurotiaceae, Thelebolales) described from Antarctic soils

The genus Pseudogymnoascus includes several species frequently isolated from extreme environments worldwide, including cold environments such as Antarctica. This study describes three new species of Pseudogymnoascus-P.russus sp. nov., P.irelandiae sp. nov., and P.ramosus sp. nov.-isolated from Antarctic soils. These species represent the first Pseudogymnoascus taxa to be formally described from Antarctic soil samples, expanding our understanding of fungal biodiversity in this extreme environment. Microscopic descriptions of asexual structures from living cultures, along with measurements of cultural characteristics and growth on various media types at different temperatures, identify three distinct new species. In addition, phylogenetic analyses based on five gene regions (ITS, LSU, MCM7, RPB2, TEF1) and whole-genome proteomes place these new species within three distinct previously described clades: P.irelandiae in clade K, P.ramosus in clade Q, and P.russus in clade B. These results provide further evidence of the extensive undescribed diversity of Pseudogymnoascus in high-latitude soils. This study contributes to the growing body of knowledge on Antarctic mycology and the broader ecology of psychrophilic and psychrotolerant fungi.

Chloroplast Genome of Arisaema takesimense: Comparative Genomics and Phylogenetic Insights into the Arisaema

Arisaema takesimense (Araceae) is a unique species found exclusively in Ulleung Island of Korea. This study presents the complete chloroplast (cp.) genome of A. takesimense, which comprises 174, 361 base pairs and exhibits a typical tetrad structure. The genome encodes 112 unique genes, including 78 protein-coding genes (CDS), 30 tRNA genes, and 4 rRNA genes. In this study, a total of 49 long repeats and 139 simple sequence repeats (SSRs) were identified, predominantly located in intergenic spacer regions (IGS). Additionally, several hotspot regions, including trnS-G, accD-psaI, ndhF and rps15-ycf1, were identified, which are commonly shared among Araceae species. The analysis of these repeats revealed species-specific SSR types and hotspot regions that can be utilized for population genetic studies and species identification. A comparative genomic analysis of eleven Arisaema taxa revealed that the large single copy region (LSC) exhibits the most variability, with non-coding genes displaying more variation than coding genes. The borders between the LSC-IR-SSC regions in Arisaema taxa were generally well-preserved, and there were notable exceptions in the positions of LSC/IRa, LSC/IRb and SSC/IRb junctions for A. takesimense, A. ringens, and A. nepenthoides. A phylogenetic analysis based on the cp. genome revealed a close relationship between A. takesimense and A. bockii. The outcomes of this study substantially increase the genomic resources available for Araceae, serving as a valuable resource for species identification and evaluating intraspecific diversity within the Arisaema genus.

Representative Metagenomes of Mesophilic Biogas Reactor Across South Korea

Biogas production through the anaerobic digestion (AD) of organic waste plays a crucial role in promoting sustainability and closing the carbon cycle. Over the past decade, this has driven global research on biogas-producing microbiomes, leading to significant advances in our understanding of microbial diversity and metabolic pathways within AD plants. However, substantial knowledge gaps persist, particularly in understanding the specific microbial communities involved in biogas production in countries such as South Korea. The present dataset addresses one of these gaps by providing comprehensive information on the metagenomes of five full-scale mesophilic biogas reactors in South Korea. From 110 GB of raw DNA sequences, 401 metagenome-assembled genomes (MAGs) were created, which include 42,301 annotated genes. Of these, 187 MAGs (46.7%) were classified as high-quality based on Minimum Information about Metagenome-Assembled Genome (MIMAG) standards. The data presented here contribute to a broader understanding of biogas-specific microbial communities and offers a significant resource for future studies and advancements in sustainable biogas production.

Developing pangenomes for large and complex plant genomes and their representation formats

BACKGROUND

The development of pangenomes has revolutionized genomic studies by capturing the complete genetic diversity within a species. Pangenome assembly integrates data from multiple individuals to construct a comprehensive genomic landscape, revealing both core and accessory genomic elements. This approach enables the identification of novel genes, structural variations, and gene presence-absence variations, providing insights into species evolution, adaptation, and trait variation. Representing pangenomes requires innovative visualization formats that effectively convey the complex genomic structures and variations.

AIM

This review delves into contemporary methodologies and recent advancements in constructing pangenomes, particularly in plant genomes. It examines the structure of pangenome representation, including format comparison, conversion, visualization techniques, and their implications for enhancing crop improvement strategies.

KEY SCIENTIFIC CONCEPTS OF REVIEW

Earlier comparative studies have illuminated novel gene sequences, copy number variations, and presence-absence variations across diverse crop species. The concept of a pan-genome, which captures multiple genetic variations from a broad spectrum of genotypes, offers a holistic perspective of a species' genetic makeup. However, constructing a pan-genome for plants with larger genomes poses challenges, including managing vast genome sequence data and comprehending the genetic variations within the germplasm. To address these challenges, researchers have explored cost-effective alternatives to encapsulate species diversity in a single assembly known as a pangenome. This involves reducing the volume of genome sequences while focusing on genetic variations. With the growing prominence of the pan-genome concept in plant genomics, several software tools have emerged to facilitate pangenome construction. This review sheds light on developing and utilizing software tools tailored for constructing pan-genomes in plants. It also discusses representation formats suitable for downstream analyses, offering valuable insights into the genetic landscape and evolutionary dynamics of plant species. In summary, this review underscores the significance of pan-genome construction and representation formats in resolving the genetic architecture of plants, particularly those with complex genomes. It provides a comprehensive overview of recent advancements, aiding in exploring and understanding plant genetic diversity.

Mitochondrial Genome and Phylogenetic Analysis of the Narrownose Smooth-Hound Shark Mustelus schmitti Springer, 1939

Southern Brazil is home to a large biodiversity of elasmobranchs from the Brazilian coast. Several genera and species of small sharks of the Triakidae family live in this marine environment. Studies on these shark species are scarce, with few genetic data and little information on animal population structures. The present study aimed to sequence the complete mitochondrial genome (mtDNA) of the endangered species Mustelus schmitti (narrownose smooth-hound shark) and to perform a phylogenetic analysis of the Triakidae family. The mtDNA sequenced here was 16,764 bp long and possessed the usual 13 mitochondrial protein coding genes (PCGs), 22 tRNAs, two rRNAs (12S and 16S) and a large D-loop DNA sequence, presenting an overall organization similar to other species from the genus Mustelus. Phylogenetic analyses were performed using a dataset containing this new mtDNA and 59 other mitochondrial genomes of the Carcharhiniformes species (including 14 from the Triakidae family), using the Maximum Likelihood (ML) method. All the species of the Triakidae family were clustered into a monophyletic topology group. In addition, polyphyly was observed in Galeorhinus galeus, Hemiatrakis japanica, Triakis megalopterus and Triakis semifasciata. In conclusion, this study contributes to a deeper understanding of the genetic diversity of sharks and represents an important step towards the conservation of these endangered animals.

Unlocking the growth-promoting and antagonistic power: A comprehensive whole genome study on Bacillus velezensis strains

Bacillus species are extensively documented as plant growth-promoting rhizobacteria, contributing significantly to the enhancement of soil fertility, nutrient recycling, and the control of phytopathogens. Utilizing them as biocontrol agents represents an environmentally friendly strategy, particularly within the rhizospheric community. This study presents the comprehensive genome sequences of three B. velezensis strains (LGMB12, LGMB319, and LGMB426) which were previously isolated from root samples of maize (Zea mays L.), along with a type strain FZB42. The research assesses the capability of the three strains for antagonizing fungi, specifically Fusarium graminearum, Fusarium verticillioides, Colletotrichum graminicola, and Stenocarpella sp. In paired cultures involving maize fungi, treatments containing bacteria B. velezensis exhibited statistically significant differences compared to both negative and positive treatments in terms of antagonism. Furthermore, genome mining techniques were employed to explore their inherent antagonistic potential. The assembly revealed that strains LGMB12, LGMB319, LGMB426, and FZB42 exhibit genome sizes of 4,187,541 bp, 4,244,954 bp, 3,976,537 bp, and 3,990,518 respectively. Their respective G + C content stands at 46.42 %, 46.50 %, 46.51 %, and 46.38 %. Moreover, the genomes present multiple gene clusters responsible for the synthesis of secondary metabolites and carbohydrate-active enzymes (CAZymes). These clusters highlight a diverse array of antibacterial and antifungal properties, complemented by numerous plant growth-promoting genes. These results highlight the potential of B. velezensis LGMB12, LGMB319, and LGMB426 strains as biocontrol and plant growth promotion agents, being promising candidates for further studies in agricultural production, including field trials.

Prediction of post-PCV13 pneumococcal evolution using invasive disease data enhanced by inverse-invasiveness weighting

After introducing pneumococcal conjugate vaccines (PCVs), serotype replacement occurred in Streptococcus pneumoniae. Predicting which pneumococcal strains will become common in carriage after vaccination can enhance vaccine design, public health interventions, and understanding of pneumococcal evolution. Invasive pneumococcal isolates were collected during 1998-2018 by the Active Bacterial Core surveillance (ABCs). Carriage data from Massachusetts (MA) and Southwest United States were used to calculate weights. Using pre-vaccine data, serotype-specific inverse-invasiveness weights were defined as the ratio of the proportion of the serotype in carriage to the proportion in invasive data. Genomic data were processed under bioinformatic pipelines to define genetically similar sequence clusters (i.e., strains), and accessory genes (COGs) present in 5-95% of isolates. Weights were applied to adjust observed strain proportions and COG frequencies. The negative frequency-dependent selection (NFDS) model predicted strain proportions by calculating the post-vaccine strain composition in the weighted invasive disease population that would best match pre-vaccine COG frequencies. Inverse-invasiveness weighting increased the correlation of COG frequencies between invasive and carriage data in linear or logit scale for pre-vaccine, post-PCV7, and post-PCV13; and between different epochs in the invasive data. Weighting the invasive data significantly improved the NFDS model's accuracy in predicting strain proportions in the carriage population in the post-PCV13 epoch, with the adjusted R2 increasing from 0.254 before weighting to 0.545 after weighting. The weighting system adjusted invasive disease data to better represent the pneumococcal carriage population, allowing the NFDS mechanism to predict strain proportions in carriage in the post-PCV13 epoch. Our methods enrich the value of genomic sequences from invasive disease surveillance.IMPORTANCEStreptococcus pneumoniae, a common colonizer in the human nasopharynx, can cause invasive diseases including pneumonia, bacteremia, and meningitis mostly in children under 5 years or older adults. The PCV7 was introduced in 2000 in the United States within the pediatric population to prevent disease and reduce deaths, followed by PCV13 in 2010, PCV15 in 2022, and PCV20 in 2023. After the removal of vaccine serotypes, the prevalence of carriage remained stable as the vacated pediatric ecological niche was filled with certain non-vaccine serotypes. Predicting which pneumococcal clones, and which serotypes, will be most successful in colonization after vaccination can enhance vaccine design and public health interventions, while also improving our understanding of pneumococcal evolution. While carriage data, which are collected from the pneumococcal population that is competing to colonize and transmit, are most directly relevant to evolutionary studies, invasive disease data are often more plentiful. Previously, evolutionary models based on negative frequency-dependent selection (NFDS) on the accessory genome were shown to predict which non-vaccine strains and serotypes were most successful in colonization following the introduction of PCV7. Here, we show that an inverse-invasiveness weighting system applied to invasive disease surveillance data allows the NFDS model to predict strain proportions in the projected carriage population in the post-PCV13/pre-PCV15 and pre-PCV20 epoch. The significance of our research lies in using a sample of invasive disease surveillance data to extend the use of NFDS as an evolutionary mechanism to predict post-PCV13 population dynamics. This has shown that we can correct for biased sampling that arises from differences in virulence and can enrich the value of genomic data from disease surveillance and advance our understanding of how NFDS impacts carriage population dynamics after both PCV7 and PCV13 vaccination.

Draft genome sequence of Paraclostridium benzoelyticum strain YGCA10, isolated from sediment of the Etoliko lagoon in Greece

Paraclostridium benzoelyticum is a Gram-positive, heterotrophic, and obligately anaerobic bacterial species. This report describes the draft genome sequence of P. benzoelyticum strain YGCA10, isolated from lagoon sediment under anaerobic culture conditions.

Long-read de novo genome assembly of Gulf toadfish (Opsanus beta)

BACKGROUND

The family Batrachoididae are a group of ecologically important teleost fishes with unique life histories, behavior, and physiology that has made them popular model organisms. Batrachoididae remain understudied in the realm of genomics, with only four reference genome assemblies available for the family, with three being highly fragmented and not up to current assembly standards. Among these is the Gulf toadfish, Opsanus beta, a model organism for serotonin physiology which has recently been bred in captivity.

RESULTS

Here we present a new, de novo genome and transcriptome assemblies for the Gulf toadfish using PacBio long read technology. The genome size of the final assembly is 2.1 gigabases, which is among the largest teleost genomes. This new assembly improves significantly upon the currently available reference for Opsanus beta with a final scaffold count of 62, of which 23 are chromosome scale, an N50 of 98,402,768, and a BUSCO completeness score of 97.3%. Annotation with ab initio and transcriptome-based methods generated 41,076 gene models. The genome is highly repetitive, with ~ 70% of the genome composed of simple repeats and transposable elements. Satellite DNA analysis identified potential telomeric and centromeric regions.

CONCLUSIONS

This improved assembly represents a valuable resource for future research using this important model organism and to teleost genomics more broadly.

Improving Beneficial Traits in Bacillus cabrialesii subsp. cabrialesii TE3T through UV-Induced Genomic Changes

It is essential to hunt for new technologies that promote sustainable practices for agroecosystems; thus, the bioprospecting of beneficial microorganisms complementing with mutation induction techniques to improve their genomic, metabolic, and functional traits is a promising strategy for the development of sustainable microbial inoculants. Bacillus cabrialesii subsp. cabrialesii strain TE3T, a previously recognized plant growth-promoting and biological control agent, was subjected to UV mutation induction to improve these agro-biotechnological traits. Dilutions were made which were spread on Petri dishes and placed under a 20 W UV lamp at 10-min intervals for 60 min. After the UV-induced mutation of this strain, 27 bacterial colonies showed morphological differences compared to the wild-type strain; however, only a strain named TE3T-UV25 showed an improvement in 53.6% of the biocontrol against Bipolaris sorokiniana vs. the wild-type strain, by competition of nutrient and space (only detected in the mutant strain), as well as diffusible metabolites. Furthermore, the ability to promote wheat growth was evaluated by carrying out experiments under specific greenhouse conditions, considering un-inoculated, strain TE3T, and strain TE3T-UV25 treatments. Thus, after 120 days, biometric traits in seedlings were quantified and statistical analyses were performed, which showed that strain TE3T-UV25 maintained its ability to promote wheat growth in comparison with the wild-type strain. On the other hand, using bioinformatics tools such as ANI, GGDC, and TYGS, the Overall Genome Relatedness Index (OGRI) and phylogenomic relationship of mutant strain TE3T-UV25 were performed, confirming that it changed its taxonomic affiliation from B. cabrialesii subsp. cabrialesii to Bacillus subtilis. In addition, genome analysis showed that the mutant, wild-type, and B. subtilis strains shared 3654 orthologous genes; however, a higher number of shared genes (3954) was found between the TE3T-UV25 mutant strain and B. subtilis 168, while the mutant strain shared 3703 genes with the wild-type strain. Genome mining was carried out using the AntiSMASH v7.0 web server and showed that mutant and wild-type strains shared six biosynthetic gene clusters associated with biocontrol but additionally, pulcherriminic acid cluster only was detected in the genome of the mutant strain and Rhizocticin A was exclusively detected in the genome of the wild-type strain. Finally, using the PlaBase tool, differences in the number of genes (17) associated with beneficial functions in agroecosystems were detected in the genome of the mutant vs. wild-type strain, such as biofertilization, bioremediation, colonizing plant system, competitive exclusion, phytohormone, plant immune response stimulation, putative functions, stress control, and biocontrol. Thus, the UV-induced mutation was a successful strategy to improve the bioactivity of B. cabrialesii subsp. cabrialesii TE3T related to the agro-biotecnology applications. The obtained mutant strain, B. subtilis TE3T-UV25, is a promising strain to be further studied as an active ingredient for the bioformulation of bacterial inoculants to migrate sustainable agriculture.

A metagenomic survey of the fecal microbiome of the African savanna elephant (Loxodonta africana)

The African savanna elephant (Loxodonta africana) is the largest terrestrial animal on Earth and is found primarily in Southern and Eastern Africa. It is a hindgut, colonic fermenter and subsists on a diet of raw plant materials found in its grazing area. In this study the bacterial, archaeal and fungal populations of seven African savanna elephant fecal metagenomes were first characterized using amplicon sequencing. On the genus level it was observed that the p-1088-a5 gut group in the bacteriome, Methanocorpusulum and Methanobrevibacter in the archaeome and Alternaria, Aurobasidium, Didymella and Preussia in the mycome, predominated. Subsequently, metagenomic shotgun sequencing was employed to identify possible functional pathways and carbohydrate-active enzymes (CAZymes). Carbohydrate catabolic pathways represented the main degradation pathways, and the fecal metagenome was enriched in the glycohydroside (GH) class of CAZymes. Additionally, the top GH families identified - GH43, GH2, GH13 and GH3 - are known to be associated with cellulytic, hemicellulytic and pectolytic activities. Finally, the CAZymes families identified in the African savanna elephant were compared with those found in the Asian elephant and it was demonstrated that there is a unique repository of CAZymes that could be leveraged in the biotechnological context such as the degradation of lignocellulose for the production of second-generation biofuels and energy.

Unveiling Agricultural Biotechnological Prospects: The Draft Genome Sequence of Stenotrophomonas geniculata LGMB417

Stenotrophomonas species are recognized as rhizobacteria that play a pivotal role in promoting plant growth by making substantial contributions to enhanced soil fertility, nutrient recycling, and phytopathogen control. Employing them as bioinputs constitutes an environmentally sound strategy, particularly within the rhizospheric community. This study revealed the draft genome sequence of Stenotrophomonas geniculata LGMB417, which was originally isolated from root samples of maize (Zea mays L.). This research assessed the potential of a bacterial strain at the molecular level through genome mining, aiming to identify genes with biotechnological significance for promoting plant growth and protection. The assembly findings indicate that strain LGMB417 possesses a genome size of 4,654,011 bp, with a G + C content of 66.50%. The draft genome sequence revealed the presence of gene clusters responsible for the synthesis of secondary metabolites and carbohydrate active enzymes (CAZymes), glycoside hydrolases (23), glycosyltransferases (18), carbohydrate esterases (5), polysaccharide lyases (2), carbohydrate-binding modules (2), and auxiliary activities (1). Several genes related to growth promotion were found in the genome, including those associated with phosphate transport and solubilization, nitrogen metabolism, siderophore production and iron transport, hormonal modulation, stress responses (such as to drought, temperature fluctuations, osmotic challenges, and oxidative conditions), and volatile organic compounds (VOCs). Subsequent phases will encompass investigations utilizing gene expression methodologies, with future explorations concentrating on facets pertinent to agricultural production, including comprehensive field studies.

Metagenome-assembled genome of withering syndrome causative agent, "Candidatus Xenohaliotis californiensis," from endangered white abalone (Haliotis sorenseni)

The genome of "Candidatus Xenohaliotis californiensis" was assembled from shotgun metagenomic sequencing of experimentally infected white abalone. Ninety-one percent genome completeness was achieved with low contamination. Sequencing this genome provides the opportunity to track pathogen evolution over time, conduct gene expression experiments, and study dynamics between this pathogen and its phage.

Strain-specific Growth Parameters are Important to Accurately Model Bacterial Growth on Baby Spinach in Simulation Models

Digital tools to predict produce shelf life have the potential to reduce food waste and improve consumer satisfaction. To address this need, we (i) performed an observational study on the microbial quality of baby spinach, (ii) completed growth experiments of bacteria that are representative of the baby spinach microbiota, and (iii) developed an initial simulation model of bacterial growth on baby spinach. Our observational data showed that the predominant genera found on baby spinach were Pseudomonas, Pantoea and Exiguobacterium. Rifampicin-resistant mutants (rifR mutants) of representative bacterial subtypes were subsequently generated to obtain strain-specific growth parameters on baby spinach. These experiments showed that: (i) it is difficult to select rifR mutants that do not have fitness costs affecting growth (9 of 15 rifR mutants showed substantial differences in growth, compared to their corresponding wild-type strain) and (ii) based on estimates from primary growth models, the mean (geometric) maximum population of rifR mutants on baby spinach (7.6 log10 CFU/g, at 6°C) appears lower than that of the spinach microbiota (9.6 log10 CFU/g, at 6°C), even if rifR mutants did not have substantial growth-related fitness costs. Thus, a simulation model, parameterized with the data obtained here as well as literature data on home refrigeration temperatures, underestimated bacterial growth on baby spinach. The root mean square error of the simulation's output, compared against data from the observational study, was 1.11 log10 CFU/g. Sensitivity analysis was used to identify key parameters (e.g., strain maximum population) that impact the simulation model's output, allowing for prioritization of future data collection to improve the simulation model. Overall, this study provides a roadmap for the development of models to predict bacterial growth on leafy vegetables with strain-specific parameters and suggests that additional data are required to improve these models.

Metagenomic sequencing of a Patescibacteria-containing enrichment from Zodletone spring in Oklahoma, USA

An enrichment of sulfidic sediments from Zodletone spring was sequenced as a metagenome. Draft genomes representing Cloacimonadota, Deltabacterota, Firmicutes, and Patescibacteria were binned and annotated and will aid functional genomics and cultivation efforts.

Draft genome sequences of two Micromonospora strains isolated from the root nodules of Alnus glutinosa

In this paper, the draft genomes of Micromonospora RTGN7 and RTP1Z1, derived from Alnus glutinosa root nodules, are reported. The assembly of RTGN7 is 6.6 Mbp, composed of 59 contigs, with an N50 of 321,872. RTP1Z1's assembly is 6.3 Mbp, composed of 151 contigs, with an N50 of 76,442 bp.

The gut metagenome harbors metabolic and antibiotic resistance signatures of moderate-to-severe asthma

Asthma is a common allergic airway disease that has been associated with the development of the human microbiome early in life. Both the composition and function of the infant gut microbiota have been linked to asthma risk, but functional alterations in the gut microbiota of older patients with established asthma remain an important knowledge gap. Here, we performed whole metagenomic shotgun sequencing of 95 stool samples from a cross-sectional cohort of 59 healthy and 36 subjects with moderate-to-severe asthma to characterize the metagenomes of gut microbiota in adults and children 6 years and older. Mapping of functional orthologs revealed that asthma contributes to 2.9% of the variation in metagenomic content even when accounting for other important clinical demographics. Differential abundance analysis showed an enrichment of long-chain fatty acid (LCFA) metabolism pathways, which have been previously implicated in airway smooth muscle and immune responses in asthma. We also observed increased richness of antibiotic resistance genes (ARGs) in people with asthma. Several differentially abundant ARGs in the asthma cohort encode resistance to macrolide antibiotics, which are often prescribed to patients with asthma. Lastly, we found that ARG and virulence factor (VF) richness in the microbiome were correlated in both cohorts. ARG and VF pairs co-occurred in both cohorts suggesting that virulence and antibiotic resistance traits are coselected and maintained in the fecal microbiota of people with asthma. Overall, our results show functional alterations via LCFA biosynthetic genes and increases in antibiotic resistance genes in the gut microbiota of subjects with moderate-to-severe asthma and could have implications for asthma management and treatment.

High-Quality Genome Assemblies of 4 Members of the Podospora anserina Species Complex

The filamentous fungus Podospora anserina is a model organism used extensively in the study of molecular biology, senescence, prion biology, meiotic drive, mating-type chromosome evolution, and plant biomass degradation. It has recently been established that P. anserina is a member of a complex of 7 closely related species. In addition to P. anserina, high-quality genomic resources are available for 2 of these taxa. Here, we provide chromosome-level annotated assemblies of the 4 remaining species of the complex, as well as a comprehensive data set of annotated assemblies from a total of 28 Podospora genomes. We find that all 7 species have genomes of around 35 Mb arranged in 7 chromosomes that are mostly collinear and less than 2% divergent from each other at genic regions. We further attempt to resolve their phylogenetic relationships, finding significant levels of phylogenetic conflict as expected from a rapid and recent diversification.

Annotated genome and transcriptome of the endangered Caribbean mountainous star coral (Orbicella faveolata) using PacBio long-read sequencing

Long-read sequencing is revolutionizing de-novo genome assemblies, with continued advancements making it more readily available for previously understudied, non-model organisms. Stony corals are one such example, with long-read de-novo genome assemblies now starting to be publicly available, opening the door for a wide array of 'omics-based research. Here we present a new de-novo genome assembly for the endangered Caribbean star coral, Orbicella faveolata, using PacBio circular consensus reads. Our genome assembly improved the contiguity (51 versus 1,933 contigs) and complete and single copy BUSCO orthologs (93.6% versus 85.3%, database metazoa_odb10), compared to the currently available reference genome generated using short-read methodologies. Our new de-novo assembled genome also showed comparable quality metrics to other coral long-read genomes. Telomeric repeat analysis identified putative chromosomes in our scaffolded assembly, with these repeats at either one, or both ends, of scaffolded contigs. We identified 32,172 protein coding genes in our assembly through use of long-read RNA sequencing (ISO-seq) of additional O. faveolata fragments exposed to a range of abiotic and biotic treatments, and publicly available short-read RNA-seq data. With anthropogenic influences heavily affecting O. faveolata, as well as its increasing incorporation into reef restoration activities, this updated genome resource can be used for population genomics and other 'omics analyses to aid in the conservation of this species.

Antimicrobial peptide glatiramer acetate targets Pseudomonas aeruginosa lipopolysaccharides to breach membranes without altering lipopolysaccharide modification

Antimicrobial peptides (AMPs) are key components of innate immunity across all domains of life. Natural and synthetic AMPs are receiving renewed attention in efforts to combat the antimicrobial resistance (AMR) crisis and the loss of antibiotic efficacy. The gram-negative pathogen Pseudomonas aeruginosa is one of the most concerning infecting bacteria in AMR, particularly in people with cystic fibrosis (CF) where respiratory infections are difficult to eradicate and associated with increased morbidity and mortality. Cationic AMPs exploit the negatively charged lipopolysaccharides (LPS) on P. aeruginosa to bind and disrupt bacterial membrane(s), causing lethal damage. P. aeruginosa modifies its LPS to evade AMP killing. Free-LPS is also a component of CF sputum and feeds pro-inflammatory cycles. Glatiramer acetate (GA) is a random peptide co-polymer-of glycine, lysine, alanine, tyrosine-used as a drug in treatment of multiple sclerosis (MS); we have previously shown GA to be an AMP which synergises with tobramycin against CF P. aeruginosa, functioning via bacterial membrane disruption. Here, we demonstrate GA's direct binding and sequestration/neutralisation of P. aeruginosa LPS, in keeping with GA's ability to disrupt the outer membrane. At CF-relevant LPS concentrations, however, membrane disruption by GA was not strongly inhibited. Furthermore, exposure to GA did not result in increased Lipid A modification of LPS or in increased gene expression of systems involved in AMP sensing and LPS modification. Therefore, despite the electrostatic targeting of LPS by GA as part of its activity, P. aeruginosa does not demonstrate LPS modification in its defence.

Draft genome sequences of five Mycobacterium strains, isolated from Alnus glutinosa root nodules

Mycobacterium is a clinically relevant genus of bacteria, with this paper reporting draft genomes of five Mycobacterium strains derived from Alnus glutinosa root nodules. The genome sizes of the isolates ranged from 6.1 to 6.9 Mbp, composed of 22-59 contigs. The N50 values ranged from 303,875 to 865,751 bp, presenting a GC% of 66.07%-66.96%.

Draft genome sequence of Streptomyces poriferorum RTGN2, a bacterial endophyte isolated from Alnus glutinosa root nodules

Herein is reported the draft genome sequence of Streptomyces poriferum RTGN2, a bacterial isolate of Alnus glutinosa root nodules, collected from Saltwell Park, Gateshead, United Kingdom. The assembly is 9.5 Mbp in size, composed of 187 contigs, with a N50 of 189,630 bp, presenting a GC content of 71.2%.

Genome sequences of four colistin-resistant ESKAPE bacterial strains isolated from patients within the same hospital

The genomes of four clinical Gram-negative ESKAPE bacterial strains highly resistant to the last-resort antibiotic colistin were sequenced and analyzed. The strains were found to carry multidrug-resistant genes besides colistin-resistant genes.

Draft genome sequence of Amycolatopsis camponoti RTGN1, a bacterial endophyte isolated from Alnus glutinosa root nodules

Here, we report the draft genome sequence of Amycolatopsis camponoti RTGN1, a bacterial endophyte of Alnus glutinosa root nodules, collected from Saltwell Park, United Kingdom. The genome is 11.9 Mbp in size, composed of 147 contigs, with an N50 of 179,211 bp and presenting a GC content of 70.9%.

Draft genome sequence of Thermococcus waiotapuensis WT1T, a thermophilic sulfur-dependent archaeon from the order Thermococcales

Thermococcus waiotapuensis WT1T is a thermophilic, peptide, and amino acid-fermenting archaeon from the order Thermococcales. It was isolated from Waiotapu, Aotearoa-New Zealand, and has a genome size of 1.80 Mbp. The genome contains 2,000 total genes, of which 1,913 encode proteins and 46 encode tRNA.

A widespread single amino acid mutation in AcrA reduces tigecycline susceptibility in Klebsiella pneumoniae

IMPORTANCE

Tigecycline, a glycecycline antibiotic with broad-spectrum activity against almost all Gram-positive and Gram-negative bacteria, is a highly concerned "last-resort" antibiotic. In addition to plasmid-hosted mobile tet(X) conferring high-level resistance to tigecycline, there are many reports suggesting increased expression of AcrAB-TolC efflux pump leads to tigecycline non-susceptibility. However, the role of mutations in AcrAB-TolC on tigecycline resistance has not been identified. This study reports a novel T188A mutation of the AcrA subunit of AcrAB-TolC complex in a clinical tigecycline-resistant Klebsiella pneumoniae strain and reveals the role of AcrA mutation on tigecycline resistance in K. pneumoniae. High prevalence of A188 type AcrA in hypervirulent multidrug-resistant K. pneumoniae indicates that mutations of the AcrAB-TolC complex may play a larger role in determining bacterial pathogenesis and antibiotic susceptibility than previously expected.

Prediction of post-PCV13 pneumococcal evolution using invasive disease data enhanced by inverse-invasiveness weighting

Background

After introduction of pneumococcal conjugate vaccines (PCVs), serotype replacement occurred in the population of Streptococcus pneumoniae. Predicting which pneumococcal clones and serotypes will become more common in carriage after vaccination can enhance vaccine design and public health interventions, while also improving our understanding of pneumococcal evolution. We sought to use invasive disease data to assess how well negative frequency-dependent selection (NFDS) models could explain pneumococcal carriage population evolution in the post-PCV13 epoch by weighting invasive data to approximate strain proportions in the carriage population.

Methods

Invasive pneumococcal isolates were collected and sequenced during 1998-2018 by the Active Bacterial Core surveillance (ABCs) from the Centers for Disease Control and Prevention (CDC). To predict the post-PCV13 population dynamics in the carriage population using a NFDS model, all genomic data were processed under a bioinformatic pipeline of assembly, annotation, and pangenome analysis to define genetically similar sequence clusters (i.e., strains) and a set of accessory genes present in 5% to 95% of the isolates. The NFDS model predicted the strain proportion by calculating the post-vaccine strain composition in the weighted invasive disease population that would best match pre-vaccine accessory gene frequencies. To overcome the biases of invasive disease data, serotype-specific inverse-invasiveness weights were defined as the ratio of the proportion of the serotype in the carriage data to the proportion in the invasive data, using data from 1998-2001 in the United States, before conjugate vaccine introduction. The weights were applied to adjust both the observed strain proportion and the accessory gene frequencies.

Results

Inverse-invasiveness weighting increased the correlation of accessory gene frequencies between invasive and carriage data with reduced residuals in linear or logit scale for pre-vaccine, post-PCV7, and post-PCV13. Similarly, weighting increased the correlation of accessory gene frequencies between different time periods in the invasive data. By weighting the invasive data, we were able to use the NFDS model to predict strain proportions in the carriage population in the post-PCV13 epoch, with the adjusted R-squared between predicted and observed strain proportions increasing from 0.176 to 0.544 after weighting.

Conclusions

The weighting system adjusted the invasive disease surveillance data to better represent the carriage population of S. pneumoniae. The NFDS mechanism predicted the strain proportions in the projected carriage population as estimated from the weighted invasive disease frequencies in the post-PCV13 epoch. Our methods enrich the value of genomic sequences from invasive disease surveillance, which is readily available, easy to collect, and of direct interest to public health.

Comparative genomic and metabolomic study of three Streptomyces sp. differing in biological activity

The Streptomyces genus is known to produce many specialized metabolites of value for medicine, but the potential of these metabolites in agronomy remains largely unexplored. In this study, we investigated three phylogenetically closely related Streptomyces strains (B5, B91, and B135) isolated from three distinct soil samples in Sudan. Despite belonging to the same species, these strains exhibited different ranges of Phytophthora infestans inhibition. The objective of this work was to identify the active compound(s) responsible for the inhibition of P. infestans and of other plant pathogens by comparing the genomes and metabolomes of the three strains which showed distinct activity patterns: B5 was the strongest inhibitor of oomycetes, B5 and B91 both inhibited most fungi and B135 was the only strain showing antibacterial activity. Our comparative genomic and metabolomic analysis identified borrelidin as the bioactive compound underlying B5's strong anti-oomycete activity and highlighted a few other metabolites as putative candidates underlying the strains' antifungal and antibacterial activities. This study illustrates the power of comparative genomics and metabolomics on phylogenetically closely related strains of differing activities to highlight bioactive compounds that could contribute to new sustainable crop protection strategies.

Metagenome-assembled genomes from sugarcane mill mud

The genomes of 11 bacteria and 3 archaea were assembled from metagenomic DNA extracted from sugarcane mill mud. These metagenome-assembled genomes ranged from 1.79 to 6.45 Mb, with 2,263 to 5,551 predicted proteins, 80.65% to 100% genome completeness, and 43.19% to 68.02% G+C content.

Draft genome sequence of Paenibacillus sp. strain RC67, an isolate from a long-term forest soil warming experiment in Petersham, Massachusetts

Paenibacillus sp. strain RC67 was isolated from the Harvard Forest long-term soil warming experiment. The assembled genome is a single contig with 7,963,753 bp and 99.4% completion. Genome annotation suggests that the isolate is of a novel bacterial species.

High-quality genomes of Paenibacillus spp. RC334 and RC343, isolated from a long-term forest soil warming experiment

Paenibacillus spp. RC334 and RC343 were isolated from heated soil in a long-term soil warming experiment. Both genomes were 5.98 Mb and assembled as a single contig. We describe the assembly and annotation of the two high-quality draft genomes for these isolates here.

A bacterial symbiont in the gill of the marine scallop Argopecten irradians irradians metabolizes dimethylsulfoniopropionate

Microbial lysis of dimethylsulfoniopropionate (DMSP) is a key step in marine organic sulfur cycling and has been recently demonstrated to play an important role in mediating interactions between bacteria, algae, and zooplankton. To date, microbes that have been found to lyse DMSP are largely confined to free-living and surface-attached bacteria. In this study, we report for the first time that a symbiont (termed "Rhodobiaceae bacterium HWgs001") in the gill of the marine scallop Argopecten irradians irradians can lyse and metabolize DMSP. Analysis of 16S rRNA gene sequences suggested that HWgs001 accounted for up to 93% of the gill microbiota. Microscopic observations suggested that HWgs001 lived within the gill tissue. Unlike symbionts of other bivalves, HWgs001 belongs to Alphaproteobacteria rather than Gammaproteobacteria, and no genes for carbon fixation were identified in its small genome. Moreover, HWgs001 was found to possess a dddP gene, responsible for the lysis of DMSP to acrylate. The enzymatic activity of dddP was confirmed using the heterologous expression, and in situ transcription of the gene in scallop gill tissues was demonstrated using reverse-transcription PCR. Together, these results revealed a taxonomically and functionally unique symbiont, which represents the first-documented DMSP-metabolizing symbiont likely to play significant roles in coastal marine ecosystems.

Genomic analysis of 61 Chlamydia psittaci strains reveals extensive divergence associated with host preference

BACKGROUND

Chlamydia (C.) psittaci, the causative agent of avian chlamydiosis and human psittacosis, is a genetically heterogeneous species. Its broad host range includes parrots and many other birds, but occasionally also humans (via zoonotic transmission), ruminants, horses, swine and rodents. To assess whether there are genetic markers associated with host tropism we comparatively analyzed whole-genome sequences of 61 C. psittaci strains, 47 of which carrying a 7.6-kbp plasmid.

RESULTS

Following clean-up, reassembly and polishing of poorly assembled genomes from public databases, phylogenetic analyses using C. psittaci whole-genome sequence alignment revealed four major clades within this species. Clade 1 represents the most recent lineage comprising 40/61 strains and contains 9/10 of the psittacine strains, including type strain 6BC, and 10/13 of human isolates. Strains from different non-psittacine hosts clustered in Clades 2- 4. We found that clade membership correlates with typing schemes based on SNP types, ompA genotypes, multilocus sequence types as well as plasticity zone (PZ) structure and host preference. Genome analysis also revealed that i) sequence variation in the major outer membrane porin MOMP can result in 3D structural changes of immunogenic domains, ii) past host change of Clade 3 and 4 strains could be associated with loss of MAC/perforin in the PZ, rather than the large cytotoxin, iii) the distinct phylogeny of atypical strains (Clades 3 and 4) is also reflected in their repertoire of inclusion proteins (Inc family) and polymorphic membrane proteins (Pmps).

CONCLUSIONS

Our study identified a number of genomic features that can be correlated with the phylogeny and host preference of C. psittaci strains. Our data show that intra-species genomic divergence is associated with past host change and includes deletions in the plasticity zone, structural variations in immunogenic domains and distinct repertoires of virulence factors.

Complete genome sequencing of Bacillus cabrialesii TE3T: A plant growth-promoting and biological control agent isolated from wheat (Triticum turgidum subsp. durum) in the Yaqui Valley

Bacillus cabrialesii TE3^T is a strictly aerobic and Gram-stain-positive plant growth-promoting bacterium, motile and catalase-positive. In addition, strain TE3^T was also recently described as a biological control agent. Here, we present the complete circularized genome of this type strain, as well as a whole genome analysis identifying genes of agricultural interest. Thus, a hybrid assembly method was performed using short-read sequencing through the Illumina MiSeq platform, and long-read sequencing through the MinION sequencing technology by Oxford Nanopore Technology (ONT). This assembly method showed a closed circular chromosome of 4,125,766 bp and 44.2% G + C content. The strain TE3^T genome annotation, based on the RAST platform, presented 4,282 Coding DNA sequences (CDS) distributed in 335 subsystems, from which 4 CDS are related to the promotion of plant growth and 28 CDS to biological control. Also, Prokka (Rapid Prokaryotic Genome Annotation) predicted a total of 119 RNAs composed of 87 tRNAs, 31 rRNA, and 1 tmRNA; and the PGAP (Prokaryotic Genome Annotation Pipeline) predicted a total of 4,212 genes (3,991 CDS). Additionally, seven putative biosynthetic gene clusters were identified by antiSMASH, such as Fengycin, Bacilysin, Subtilosin A, Bacillibactin, Bacillaene, Surfactin, and Rizocticin A, which are related to antimicrobial and antifungal properties, whose gene presence was further supported by the Prokaryotic Genome Annotation Pipeline (PGAP) annotation. Thus, the complete genome of Bacillus cabrialesii TE3^T showed promising bioactivities for the use of this type strain to bioformulate bacterial inoculants for sustainable agriculture.

Draft Genome Sequences of 14 Fungal Species from Alternaria Section Infectoriae

Isolates representing 14 Alternaria section Infectoriae species collected from diverse hosts and locations were shotgun sequenced. These genome assemblies and annotations will improve the study of taxonomy and biology in this group of ubiquitous fungi.

The gut metagenome harbors metabolic and antibiotic resistance signatures of moderate-to-severe asthma

Asthma is a common allergic airway disease that develops in association with the human microbiome early in life. Both the composition and function of the infant gut microbiota have been linked to asthma risk, but functional alterations in the gut microbiota of older patients with established asthma remain an important knowledge gap. Here, we performed whole metagenomic shotgun sequencing of 95 stool samples from 59 healthy and 36 subjects with moderate-to-severe asthma to characterize the metagenomes of gut microbiota in children and adults 6 years and older. Mapping of functional orthologs revealed that asthma contributes to 2.9% of the variation in metagenomic content even when accounting for other important clinical demographics. Differential abundance analysis showed an enrichment of long-chain fatty acid (LCFA) metabolism pathways which have been previously implicated in airway smooth muscle and immune responses in asthma. We also observed increased richness of antibiotic resistance genes (ARGs) in people with asthma. One differentially abundant ARG was a macrolide resistance marker, ermF, which significantly co-occurred with the Bacteroides fragilis toxin, suggesting a possible relationship between enterotoxigenic B. fragilis, antibiotic resistance, and asthma. Lastly, we found multiple virulence factor (VF) and ARG pairs that co-occurred in both cohorts suggesting that virulence and antibiotic resistance traits are co-selected and maintained in the fecal microbiota of people with asthma. Overall, our results show functional alterations via LCFA biosynthetic genes and increases in antibiotic resistance genes in the gut microbiota of subjects with moderate-to-severe asthma and could have implications for asthma management and treatment.

Comparative analysis of Penicillium genomes reveals the absence of a specific genetic basis for biocontrol in Penicillium rubens strain 212

Penicillium rubens strain 212 (PO212) is a filamentous fungus belonging to the division Ascomycete. PO212 acts as an effective biocontrol agent against several pathogens in a variety of horticultural crops including Fusarium oxysporum f.sp. lycopersici, causing vascular wilt disease in tomato plants. We assembled draft genomes of two P. rubens strains, the biocontrol agent PO212 and the soil isolate S27, which lacks biocontrol activity. We also performed comparative analyses of the genomic sequence of PO212 with that of the other P. rubens and P. chrysogenum strains. This is the first Penicillium strain with biocontrol activity whose genome has been sequenced and compared. PO212 genome size is 2,982 Mb, which is currently organized into 65 scaffolds and a total of 10,164 predicted Open Reading Frames (ORFs). Sequencing confirmed that PO212 belongs to P. rubens clade. The comparative analysis of the PO212 genome with the genomes of other P. rubens and Penicillium chrysogenum strains available in databases showed strong conservation among genomes, but a correlation was not found between these genomic data and the biocontrol phenotype displayed by PO212. Finally, the comparative analysis between PO212 and S27 genomes showed high sequence conservation and a low number of variations mainly located in ORF regions. These differences found in coding regions between PO212 and S27 genomes can explain neither the biocontrol activity of PO212 nor the absence of such activity in S27, opening a possible avenue toward transcriptomic and epigenetic studies that may shed light on this mechanism for fighting plant diseases caused by fungal pathogens. The genome sequences described in this study provide a useful novel resource for future research into the biology, ecology, and evolution of biological control agents.

Phenotypic and genomic characterization provide new insights into adaptation to environmental stressors and biotechnological relevance of mangrove Alcaligenes faecalis D334

Alcaligenes faecalis D334 was determined in this study as a salt-tolerant bacterium isolated from mangrove sediment. In response to 6% (w/v) NaCl, strain D334 produced the highest ectoines of 14.14 wt%. To understand adaptive features to mangrove environment, strain D334 was sequenced using Pacific BioScience platform, resulting in a circular chromosome of 4.23 Mb. Of note, D334 genome harbored 81 salt-responsive genes, among which two membrane-associated genes ompc and eric were absent in 3 selected A. faecalis genomes. Apart from that, a complete pathway for ectoine and 5-hydroxyectoine synthesis was predicted. To resist 40 mM H₂O₂, 46 genetic determinants contributing to oxidative stress response were employed. Moreover, two operons involved in polyhydroxyalkanoate (PHA) production were identified in the D334 genome, resulting in maximum PHA content of 5.03 ± 0.04 wt% and PHA concentration of 0.13 ± 0.001 g/L. A large flagellar biosynthesis operon contributing to swimming motility was found to be conserved in D334 and 8 other A. faecalis genomes. These findings shed light for the first time on the high versatility of A. faecalis D334 genome to adapt to mangrove lifestyle and the possibility to develop D334 as an industrial platform for PHA and 5-hydroxyectoine production.

NPGREAT: assembly of human subtelomere regions with the use of ultralong nanopore reads and linked-reads

BACKGROUND

Human subtelomeric DNA regulates the length and stability of adjacent telomeres that are critical for cellular function, and contains many gene/pseudogene families. Large evolutionarily recent segmental duplications and associated structural variation in human subtelomeres has made complete sequencing and assembly of these regions difficult to impossible for many loci, complicating or precluding a wide range of genetic analyses to investigate their function.

RESULTS

We present a hybrid assembly method, NanoPore Guided REgional Assembly Tool (NPGREAT), which combines Linked-Read data with mapped ultralong nanopore reads spanning subtelomeric segmental duplications to potentially overcome these difficulties. Linked-Read sets of DNA sequences identified by matches with 1-copy subtelomere sequence adjacent to segmental duplications are assembled and extended into the segmental duplication regions using Regional Extension of Assemblies using Linked-Reads (REXTAL). Mapped telomere-containing ultralong nanopore reads are then used to provide contiguity and correct orientation for matching REXTAL sequence contigs as well as identification/correction of any misassemblies. Our method was tested for a subset of representative subtelomeres with ultralong nanopore read coverage in the haploid human cell line CHM13. A 10X Linked-Read dataset from CHM13 was combined with ultralong nanopore reads from the same genome to provide improved subtelomere assemblies. Comparison of Nanopore-only assemblies using SHASTA with our NPGREAT assemblies in the distal-most subtelomere regions showed that NPGREAT produced higher-quality and more complete assemblies than SHASTA alone when these regions had low ultralong nanopore coverage (such as cases where large segmental duplications were immediately adjacent to (TTAGGG) tracts).

CONCLUSION

In genomic regions with large segmental duplications adjacent to telomeres, NPGREAT offers an alternative economical approach to improving assembly accuracy and coverage using linked-read datasets when more expensive HiFi datasets of 10-20 kb reads are unavailable.

Genome sequence and silkomics of the spindle ermine moth, Yponomeuta cagnagella, representing the early diverging lineage of the ditrysian Lepidoptera

Many lepidopteran species produce silk, cocoons, feeding tubes, or nests for protection from predators and parasites for caterpillars and pupae. Yet, the number of lepidopteran species whose silk composition has been studied in detail is very small, because the genes encoding the major structural silk proteins tend to be large and repetitive, making their assembly and sequence analysis difficult. Here we have analyzed the silk of Yponomeuta cagnagella, which represents one of the early diverging lineages of the ditrysian Lepidoptera thus improving the coverage of the order. To obtain a comprehensive list of the Y. cagnagella silk genes, we sequenced and assembled a draft genome using Oxford Nanopore and Illumina technologies. We used a silk-gland transcriptome and a silk proteome to identify major silk components and verified the tissue specificity of expression of individual genes. A detailed annotation of the major genes and their putative products, including their complete sequences and exon-intron structures is provided. The morphology of silk glands and fibers are also shown. This study fills an important gap in our growing understanding of the structure, evolution, and function of silk genes and provides genomic resources for future studies of the chemical ecology of Yponomeuta species.

Lactobacillus reuteri tryptophan metabolism promotes host susceptibility to CNS autoimmunity

BACKGROUND

Dysregulation of gut microbiota-associated tryptophan metabolism has been observed in patients with multiple sclerosis. However, defining direct mechanistic links between this apparent metabolic rewiring and individual constituents of the gut microbiota remains challenging. We and others have previously shown that colonization with the gut commensal and putative probiotic species, Lactobacillus reuteri, unexpectedly enhances host susceptibility to experimental autoimmune encephalomyelitis (EAE), a murine model of multiple sclerosis. To identify underlying mechanisms, we characterized the genome of commensal L. reuteri isolates, coupled with in vitro and in vivo metabolomic profiling, modulation of dietary substrates, and gut microbiota manipulation.

RESULTS

The enzymes necessary to metabolize dietary tryptophan into immunomodulatory indole derivatives were enriched in the L. reuteri genomes, including araT, fldH, and amiE. Moreover, metabolite profiling of L. reuteri monocultures and serum of L. reuteri-colonized mice revealed a depletion of kynurenines and production of a wide array of known and novel tryptophan-derived aryl hydrocarbon receptor (AhR) agonists and antagonists, including indole acetate, indole-3-glyoxylic acid, tryptamine, p-cresol, and diverse imidazole derivatives. Functionally, dietary tryptophan was required for L. reuteri-dependent EAE exacerbation, while depletion of dietary tryptophan suppressed disease activity and inflammatory T cell responses in the CNS. Mechanistically, L. reuteri tryptophan-derived metabolites activated the AhR and enhanced T cell production of IL-17.

CONCLUSIONS

Our data suggests that tryptophan metabolism by gut commensals, such as the putative probiotic species L. reuteri, can unexpectedly enhance autoimmunity, inducing broad shifts in the metabolome and immunological repertoire. Video Abstract.

Genomic trajectories of a near-extinction event in the Chatham Island black robin

BACKGROUND

Understanding the micro--evolutionary response of populations to demographic declines is a major goal in evolutionary and conservation biology. In small populations, genetic drift can lead to an accumulation of deleterious mutations, which will increase the risk of extinction. However, demographic recovery can still occur after extreme declines, suggesting that natural selection may purge deleterious mutations, even in extremely small populations. The Chatham Island black robin (Petroica traversi) is arguably the most inbred bird species in the world. It avoided imminent extinction in the early 1980s and after a remarkable recovery from a single pair, a second population was established and the two extant populations have evolved in complete isolation since then. Here, we analysed 52 modern and historical genomes to examine the genomic consequences of this extreme bottleneck and the subsequent translocation.

RESULTS

We found evidence for two-fold decline in heterozygosity and three- to four-fold increase in inbreeding in modern genomes. Moreover, there was partial support for temporal reduction in total load for detrimental variation. In contrast, compared to historical genomes, modern genomes showed a significantly higher realised load, reflecting the temporal increase in inbreeding. Furthermore, the translocation induced only small changes in the frequency of deleterious alleles, with the majority of detrimental variation being shared between the two populations.

CONCLUSION

Our results highlight the dynamics of mutational load in a species that recovered from the brink of extinction, and show rather limited temporal changes in mutational load. We hypothesise that ancestral purging may have been facilitated by population fragmentation and isolation on several islands for thousands of generations and may have already reduced much of the highly deleterious load well before human arrival and introduction of pests to the archipelago. The majority of fixed deleterious variation was shared between the modern populations, but translocation of individuals with low mutational load could possibly mitigate further fixation of high-frequency deleterious variation.

Chromosome-level reference genome for European flat oyster (Ostrea edulis L.)

The European flat oyster (Ostrea edulis L.) is a bivalve naturally distributed across Europe, which was an integral part of human diets for centuries, until anthropogenic activities and disease outbreaks severely reduced wild populations. Despite a growing interest in genetic applications to support population management and aquaculture, a reference genome for this species is lacking to date. Here, we report a chromosome-level assembly and annotation for the European Flat oyster genome, generated using Oxford Nanopore, Illumina, Dovetail OmniC™ proximity ligation and RNA sequencing. A contig assembly (N50: 2.38 Mb) was scaffolded into the expected karyotype of 10 pseudochromosomes. The final assembly is 935.13 Mb, with a scaffold-N50 of 95.56 Mb, with a predicted repeat landscape dominated by unclassified elements specific to O. edulis. The assembly was verified for accuracy and completeness using multiple approaches, including a novel linkage map built with ddRAD-Seq technology, comprising 4016 SNPs from four full-sib families (eight parents and 163 F1 offspring). Annotation of the genome integrating multitissue transcriptome data, comparative protein evidence and ab-initio gene prediction identified 35,699 protein-coding genes. Chromosome-level synteny was demonstrated against multiple high-quality bivalve genome assemblies, including an O. edulis genome generated independently for a French O. edulis individual. Comparative genomics was used to characterize gene family expansions during Ostrea evolution that potentially facilitated adaptation. This new reference genome for European flat oyster will enable high-resolution genomics in support of conservation and aquaculture initiatives, and improves our understanding of bivalve genome evolution.

Culture-Free Phylogenetic Analysis of Legionella pneumophila Using Targeted CRISPR/Cas9 Next-Generation Sequencing

Currently available methods for the laboratory investigation of Legionella pneumophila outbreaks require organism culture. The ability to sequence L. pneumophila directly from clinical samples would significantly reduce delays. Here, we develop a method for targeted next-generation sequencing (NGS) of selected L. pneumophila genes utilizing a CRISPR/Cas9-based target enrichment system. We determine the method's utility by typing cultured L. pneumophila isolates and subsequently apply the method directly to patient samples. We sequenced 10 L. pneumophila isolates by 2 methods, (i) whole-genome sequencing (WGS) and (ii) targeted (CRISPR/Cas9-based) finding low-abundance sequences by hybridization (FLASH)-NGS, sequencing 57 selected genes. The targeted NGS of 57 genes was more efficient than WGS, and phylogenetic analysis of the 57 genes yielded the same classification of the L. pneumophila isolates as that based on analysis of whole-genome data. Furthermore, targeted NGS of L. pneumophila performed directly on patient respiratory samples correctly classified the patients according to their corresponding cultured isolates. This provides proof of concept that targeted NGS can be used to sequence L. pneumophila directly from patient samples. Studies on a larger number of patient samples will further validate this method. Nonetheless, CRISPR/Cas9 targeted NGS methods have the potential to be widely applicable to microbial-outbreak investigations in the future, particularly in the context of difficult and slow-growing organisms. IMPORTANCE The bacterium Legionella pneumophila is responsible for outbreaks of serious and life-threatening pneumonia called Legionnaires' disease. There is a need for new molecular methods that allow investigation of Legionella outbreaks directly from patient samples, without the need for prior microbiological culture, which causes delays. Our study aims to address this problem. We have utilized a CRISPR/Cas9-based targeted next-generation sequencing (NGS) method that can be applied directly on human specimens. Furthermore, we show that analysis of the sequences of a small number of targeted genes offers the same classification of L. pneumophila as that based on data derived from the whole genome. Given the rising interest globally in sequencing pathogens directly from human samples, CRISPR/Cas9 targeted NGS methods have the potential to be widely applicable to microbial-outbreak investigations in the future, particularly in the context of difficult and slow-growing organisms.

Interspecies Transmission of CMY-2-Producing Escherichia coli Sequence Type 963 Isolates between Humans and Gulls in Australia

We have provided the first comprehensive genomic study of E. coli ST963 by analyzing various genomic and phenotypic data sets of isolates from Australian silver gulls and comparison with genomes from geographically dispersed regions of human and animal origin. Our study suggests the emergence of a specific bla CMY-2 -carrying E. coli ST963 clone in Australia that is widely spread across the continent by humans and birds.

Benchmarking and Assessment of Eight De Novo Genome Assemblers on Viral Next-Generation Sequencing Data, Including the SARS-CoV-2

Viral genomics has become crucial in clinical diagnostics and ecology, not to mention to stem the COVID-19 pandemic. Whole-genome sequencing (WGS) is pivotal in gaining an improved understanding of viral evolution, genomic epidemiology, infectious outbreaks, pathobiology, clinical management, and vaccine development. Genome assembly is one of the crucial steps in WGS data analyses. A series of different assemblers has been developed with the advent of high-throughput next-generation sequencing (NGS). Various studies have reported the evaluation of these assembly tools on distinct datasets; however, these lack data from viral origin. In this study, we performed a comparative evaluation and benchmarking of eight de novo assemblers: SOAPdenovo, Velvet, assembly by short sequences (ABySS), iterative De Bruijn graph assembler (IDBA), SPAdes, Edena, iterative virus assembler, and VICUNA on the viral NGS data from distinct Illumina (GAIIx, Hiseq, Miseq, and Nextseq) platforms. WGS data of diverse viruses, that is, severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), dengue virus 3, human immunodeficiency virus 1, hepatitis B virus, human herpesvirus 8, human papillomavirus 16, rhinovirus A, and West Nile virus, were utilized to assess these assemblers. Performance metrics such as genome fraction recovery, assembly lengths, NG50, N50, contig length, contig numbers, mismatches, and misassemblies were analyzed. Overall, three assemblers, that is, SPAdes, IDBA, and ABySS, performed consistently well, including for genome assembly of SARS-CoV-2. These assembly methods should be considered and recommended for future studies of viruses. The study also suggests that implementing two or more assembly approaches should be considered in viral NGS studies, especially in clinical settings. Taken together, the benchmarking of eight de novo genome assemblers reported in this study can inform future public health and ecology research concerning the viruses, the COVID-19 pandemic, and viral outbreaks.

Transposon activity, local duplications and propagation of structural variants across haplotypes drive the evolution of the Drosophila S2 cell line

BACKGROUND

Immortalized cell lines are widely used model systems whose genomes are often highly rearranged and polyploid. However, their genome structure is seldom deciphered and is thus not accounted for during analyses. We therefore used linked short- and long-read sequencing to perform haplotype-level reconstruction of the genome of a Drosophila melanogaster cell line (S2-DRSC) with a complex genome structure.

RESULTS

Using a custom implementation (that is designed to use ultra-long reads in complex genomes with nested rearrangements) to call structural variants (SVs), we found that the most common SV was repetitive sequence insertion or deletion (> 80% of SVs), with Gypsy retrotransposon insertions dominating. The second most common SV was local sequence duplication. SNPs and other SVs were rarer, but several large chromosomal translocations and mitochondrial genome insertions were observed. Haplotypes were highly similar at the nucleotide level but structurally very different. Insertion SVs existed at various haplotype frequencies and were unlinked on chromosomes, demonstrating that haplotypes have different structures and suggesting the existence of a mechanism that allows SVs to propagate across haplotypes. Finally, using public short-read data, we found that transposable element insertions and local duplications are common in other D. melanogaster cell lines.

CONCLUSIONS

The S2-DRSC cell line evolved through retrotransposon activity and vast local sequence duplications, that we hypothesize were the products of DNA re-replication events. Additionally, mutations can propagate across haplotypes (possibly explained by mitotic recombination), which enables fine-tuning of mutational impact and prevents accumulation of deleterious events, an inherent problem of clonal reproduction. We conclude that traditional linear homozygous genome representation conceals the complexity when dealing with rearranged and heterozygous clonal cells.