Mugsy: fast multiple alignment of closely related whole genomes

The physical biogeography of Fusobacterium nucleatum in health and disease

Fusobacterium nucleatum (Fn) is an oral commensal inhabiting the human gingival plaque that is rarely found in the gut. However, in colorectal cancer (CRC), Fn can be isolated from stool samples and detected in metagenomes. We hypothesized that ecological characteristics of the gut are altered by disease, enabling Fn to colonize. Multiple genomically distinct populations of Fn exist, but their ecological preferences are unstudied. We identified six well-separated populations in 133 Fn genomes and used simulated metagenomes to demonstrate sensitive detection of populations in human oral and gut metagenomes. In 9,560 samples from 11 studies, Fn population C2 animalis is elevated in gut metagenomes from CRC and Crohn's disease patients and is observed more frequently in CRC stool samples than in the gingiva. Polymorphum, the most prevalent gingival Fn population, is significantly increased in Crohn's stool samples; this effect was significantly stronger in male hosts than in female. We find polymorphum genomes are enriched for biosynthetic gene clusters and fluoride exporters, while C2 animalis are high in iron transporters. Fn populations thus associate with specific clinical and demographic phenotypes and harbor distinct functional features. Ecological differences in closely related groups of bacteria inform microbiome impacts on human health.

IMPORTANCE

Fusobacterium nucleatum is a bacterium normally found in the gingiva. F. nucleatum generally does not colonize the healthy gut, but is observed in approximately a third of colorectal cancer (CRC) patient guts. F. nucleatum's presence in the gut during CRC has been linked to worse prognosis and increased tumor proliferation. Here, we describe the population structure of F. nucleatum in oral and gut microbiomes. We report substantial diversity in gene carriage among six distinct populations of F. nucleatum and identify population disease and body-site preferences. We find the C2 animalis population is more common in the CRC gut than in the gingiva and is enriched for iron transporters, which support gut colonization in known pathogens. We find that C2 animalis is also enriched in Crohn's disease and type 2 diabetes, suggesting ecological commonalities between the three diseases. Our work shows that closely related bacteria can have different associations with human physiology.

Whole-genome automated assembly pipeline for Chlamydia trachomatis strains from reference, in vitro and clinical samples using the integrated CtGAP pipeline

Whole genome sequencing (WGS) is pivotal for the molecular characterization of Chlamydia trachomatis (Ct)-the leading bacterial cause of sexually transmitted infections and infectious blindness worldwide. Ct WGS can inform epidemiologic, public health and outbreak investigations of these human-restricted pathogens. However, challenges persist in generating high-quality genomes for downstream analyses given its obligate intracellular nature and difficulty with in vitro propagation. No single tool exists for the entirety of Ct genome assembly, necessitating the adaptation of multiple programs with varying success. Compounding this issue is the absence of reliable Ct reference strain genomes. We, therefore, developed CtGAP-Chlamydia trachomatisGenome Assembly Pipeline-as an integrated 'one-stop-shop' pipeline for assembly and characterization of Ct genome sequencing data from various sources including isolates, in vitro samples, clinical swabs and urine. CtGAP, written in Snakemake, enables read quality statistics output, adapter and quality trimming, host read removal, de novo and reference-guided assembly, contig scaffolding, selective ompA, multi-locus-sequence and plasmid typing, phylogenetic tree construction, and recombinant genome identification. Twenty Ct reference genomes were also generated. Successfully validated on a diverse collection of 363 samples containing Ct, CtGAP represents a novel pipeline requiring minimal bioinformatics expertise with easy adaptation for use with other bacterial species.

Re-evaluation of boundaries of Streptococcus mitis and Streptococcus oralis and demonstration of multiple later synonyms of Streptococcus mitis, Streptococcus oralis and Streptococcus thalassemiae: description of Streptococcus mitis subsp. carlssonii subsp. nov. and emended description of Streptococcus mitis

The commensal species Streptococcus mitis and Streptococcus oralis are genetically diverse to a degree that challenges traditional definitions of species. This causes automatic identification based on DNA sequences or cellular extract profiles problematic. Based on an initial analysis of 266 genomes, we subjected a subset of 100 representative genomes to detailed phylogenetic, pairwise distance and gene pattern analyses. S. mitis and S. oralis constitute a continuum of clones that are genetically unique. To recognize most isolates as separate species is biologically and practically meaningless. We recommend bending the proposed similarity borders to accommodate the biological reality. Accordingly, we conclude that Streptococcus toyakuensis, Streptococcus chosunensis, Streptococcus gwangjuensis, Streptococcus humanilactis and Streptococcus hohhotensis are later heterotypic synonyms of S. mitis. Type strains of effectively but not validly published 'Streptococcus shenyangsis', 'Streptococcus symci' and 'Streptococcus vulneris' belong in S. mitis. Streptococcus parapneumoniae and Streptococcus nakanonensis are later synonyms of Streptococcus thalassemiae. Streptococcus downii is a later synonym of Streptococcus oralis subsp. dentisani, and the type of 'Streptococcus halitosis' belongs in Streptococcus oralis subsp. tigurinus. The genome sequence of the type of the recently proposed 'Streptococcus bouchesdurhonensis' is based on a mixed culture. Phylogenetic results and the pattern of presence/absence of accessory genes support the distinction of two subspecies of S. mitis, i.e. Sreptococcus mitis subsp. mitis subsp. nov. (type strain is NCTC 12261T) and Sreptococcus mitis subsp. carlssonii subsp. nov. (type strain is SK608=CCUG 55085T=LMG 33510T). The special population structure of the Streptococcus mitis-pneumoniae-pseudopneumoniae-thalassemiae complex renders automated classification of isolates based on average nucleotide identity or digital DNA-DNA hybridization values problematic. As an alternative, for initial taxonomic assignment, we present a whole-genome phylogeny-based method that enables phylogenetic comparison of new isolates in the context of a set of 117 well-characterized reference strains assigned to the Mitis/Sanguinis group.

Comparative Genomic Analysis of an Apiotrichum cacaoliposimilis Strain Isolated from a Patient with Urinary Tract Infection

Opportunistic infections caused by fungi, particularly those occurring in immunocompromised patients, are considered challenging worldwide. Therefore, a comprehensive understanding of pathogenic fungi is necessary. The present study reports the isolation of a strain of Apiotrichum cacaoliposimilis, which is difficult to detect using conventional clinical assays, from the sterile urine samples of a patient with a urinary tract infection. Sanger sequencing of the internal transcribed spacer regions confirmed the genus of the microbe, while whole-genome sequencing yielded the initial genome assembly of A. cacaoliposimilis. A total of 7,161 predicted proteincoding genes were mapped using multiple databases, including Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, non-redundant protein database, Pathogen-Host Interactions Database, and Comprehensive Antibiotic Resistance Database. The phenotypic data, biochemical reactions, and antimicrobial susceptibility analyses were conducted to reveal the metabolic properties, virulence, and drug resistance profile of the isolated A. cacaoliposimilis. The rank-sum test revealed the differences in the intergeneric distribution of the highly virulent genes UgeB and Pem2. In addition, other genes exhibited significant overlap in terms of virulence factors with the clinical isolate Apiotrichum mycotoxinivorans GMU1709. Fortunately, similar to most fungi belonging to the Apiotrichum genus, the isolate investigated in the present study was also sensitive to the drug voriconazole (MIC = 0.06 μg/ml). In summary, the phylogenetic placement, potential pathogenic genes, drug sensitivity patterns, and morphological characteristics of the isolated A. cacaoliposimilis were determined precisely in the present study.

A stepwise guide for pangenome development in crop plants: an alfalfa (Medicago sativa) case study

BACKGROUND

The concept of pangenomics and the importance of structural variants is gaining recognition within the plant genomics community. Due to advancements in sequencing and computational technology, it has become feasible to sequence the entire genome of numerous individuals of a single species at a reasonable cost. Pangenomes have been constructed for many major diploid crops, including rice, maize, soybean, sorghum, pearl millet, peas, sunflower, grapes, and mustards. However, pangenomes for polyploid species are relatively scarce and are available in only few crops including wheat, cotton, rapeseed, and potatoes.

MAIN BODY

In this review, we explore the various methods used in crop pangenome development, discussing the challenges and implications of these techniques based on insights from published pangenome studies. We offer a systematic guide and discuss the tools available for constructing a pangenome and conducting downstream analyses. Alfalfa, a highly heterozygous, cross pollinated and autotetraploid forage crop species, is used as an example to discuss the concerns and challenges offered by polyploid crop species. We conducted a comparative analysis using linear and graph-based methods by constructing an alfalfa graph pangenome using three publicly available genome assemblies. To illustrate the intricacies captured by pangenome graphs for a complex crop genome, we used five different gene sequences and aligned them against the three graph-based pangenomes. The comparison of the three graph pangenome methods reveals notable variations in the genomic variation captured by each pipeline.

CONCLUSION

Pangenome resources are proving invaluable by offering insights into core and dispensable genes, novel gene discovery, and genome-wide patterns of variation. Developing user-friendly online portals for linear pangenome visualization has made these resources accessible to the broader scientific and breeding community. However, challenges remain with graph-based pangenomes including compatibility with other tools, extraction of sequence for regions of interest, and visualization of genetic variation captured in pangenome graphs. These issues necessitate further refinement of tools and pipelines to effectively address the complexities of polyploid, highly heterozygous, and cross-pollinated species.

The effects of Pseudomonas strains isolated from Achnatherum inebrians on plant growth: A genomic perspective

Achnatherum inebrians is a perennial grass widely distributed in northwest China. Nearly all wild A. inebrians plants are infected by Epichloë endophytes. In this study, bacteria from the phyllosphere were isolated from leaves of both endophyte-free and endophyte-infected A. inebrians and sequenced for identification. Pseudomonas, comprising 48.12% of the culturable bacterial communities, was the most dominant bacterial genus. Thirty-four strains from 12 Pseudomonas species were used to inoculate A. inebrians seeds and plants. Results indicated that Epichloë significantly increased the diversity and richness index of the phyllosphere. Pseudomonas Sp1, Sp3, Sp5 and Sp7 had a significantly positive effect on plant growth and photosynthesis, whereas Sp10, Sp11 and Sp12 had a significantly negative effect. Whole-genome and pan-genome analysis suggested that the variability in the effects of Pseudomonas on A. inebrians was related to differences in genome composition and genomic islands.

OrthoPhyl-streamlining large-scale, orthology-based phylogenomic studies of bacteria at broad evolutionary scales

There are a staggering number of publicly available bacterial genome sequences (at writing, 2.0 million assemblies in NCBI's GenBank alone), and the deposition rate continues to increase. This wealth of data begs for phylogenetic analyses to place these sequences within an evolutionary context. A phylogenetic placement not only aids in taxonomic classification but informs the evolution of novel phenotypes, targets of selection, and horizontal gene transfer. Building trees from multi-gene codon alignments is a laborious task that requires bioinformatic expertise, rigorous curation of orthologs, and heavy computation. Compounding the problem is the lack of tools that can streamline these processes for building trees from large-scale genomic data. Here we present OrthoPhyl, which takes bacterial genome assemblies and reconstructs trees from whole genome codon alignments. The analysis pipeline can analyze an arbitrarily large number of input genomes (>1200 tested here) by identifying a diversity-spanning subset of assemblies and using these genomes to build gene models to infer orthologs in the full dataset. To illustrate the versatility of OrthoPhyl, we show three use cases: E. coli/Shigella, Brucella/Ochrobactrum and the order Rickettsiales. We compare trees generated with OrthoPhyl to trees generated with kSNP3 and GToTree along with published trees using alternative methods. We show that OrthoPhyl trees are consistent with other methods while incorporating more data, allowing for greater numbers of input genomes, and more flexibility of analysis.

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

BACKGROUND

The short-read whole-genome sequencing (WGS) approach has been widely applied to investigate the genomic variation in the natural populations of many plant species. With the rapid advancements in long-read sequencing and genome assembly technologies, high-quality genome sequences are available for a group of varieties for many plant species. These genome sequences are expected to help researchers comprehensively investigate any type of genomic variants that are missed by the WGS technology. However, multiple genome alignment (MGA) tools designed by the human genome research community might be unsuitable for plant genomes.

RESULTS

To fill this gap, we developed the AnchorWave-Cactus Multiple Genome Alignment (ACMGA) pipeline, which improved the alignment of repeat elements and could identify long (> 50 bp) deletions or insertions (INDELs). We conducted MGA using ACMGA and Cactus for 8 Arabidopsis (Arabidopsis thaliana) and 26 Maize (Zea mays) de novo assembled genome sequences and compared them with the previously published short-read variant calling results. MGA identified more single nucleotide variants (SNVs) and long INDELs than did previously published WGS variant callings. Additionally, ACMGA detected significantly more SNVs and long INDELs in repetitive regions and the whole genome than did Cactus. Compared with the results of Cactus, the results of ACMGA were more similar to the previously published variants called using short-read. These two MGA pipelines identified numerous multi-allelic variants that were missed by the WGS variant calling pipeline.

CONCLUSIONS

Aligning de novo assembled genome sequences could identify more SNVs and INDELs than mapping short-read. ACMGA combines the advantages of AnchorWave and Cactus and offers a practical solution for plant MGA by integrating global alignment, a 2-piece-affine-gap cost strategy, and the progressive MGA algorithm.

Fine mapping of major QTL qshgd1 for spontaneous haploid genome doubling in maize (Zea mays L.)

KEY MESSAGE

A large-effect QTL was fine mapped, which revealed 79 gene models, with 10 promising candidate genes, along with a novel inversion. In commercial maize breeding, doubled haploid (DH) technology is arguably the most efficient resource for rapidly developing novel, completely homozygous lines. However, the DH strategy, using in vivo haploid induction, currently requires the use of mutagenic agents which can be not only hazardous, but laborious. This study focuses on an alternative approach to develop DH lines-spontaneous haploid genome duplication (SHGD) via naturally restored haploid male fertility (HMF). Inbred lines A427 and Wf9, the former with high HMF and the latter with low HMF, were selected to fine-map a large-effect QTL associated with SHGD-qshgd1. SHGD alleles were derived from A427, with novel haploid recombinant groups having varying levels of the A427 chromosomal region recovered. The chromosomal region of interest is composed of 45 megabases (Mb) of genetic information on chromosome 5. Significant differences between haploid recombinant groups for HMF were identified, signaling the possibility of mapping the QTL more closely. Due to suppression of recombination from the proximity of the centromere, and a newly discovered inversion region, the associated QTL was only confined to a 25 Mb region, within which only a single recombinant was observed among ca. 9,000 BC1 individuals. Nevertheless, 79 gene models were identified within this 25 Mb region. Additionally, 10 promising candidate genes, based on RNA-seq data, are described for future evaluation, while the narrowed down genome region is accessible for straightforward introgression into elite germplasm by BC methods.

A Novel Bacillus Velezensis for Efficient Degradation of Zearalenone

Zearalenone (ZEN) is considered one of the most serious mycotoxins contaminating grains and their by-products, causing significant economic losses in the feed and food industries. Biodegradation pathways are currently considered the most efficient solution to remove ZEN contamination from foods. However, low degradation rates and vulnerability to environmental impacts limit the application of biodegradation pathways. Therefore, the main research objective of this article was to screen strains that can efficiently degrade ZEN and survive under harsh conditions. This study successfully isolated a new strain L9 which can efficiently degrade ZEN from 108 food ingredients. The results of sequence alignment showed that L9 is Bacillus velezensis. Meanwhile, we found that the L9 degradation rate reached 91.14% at 24 h and confirmed that the primary degradation mechanism of this strain is biodegradation. The strain exhibits resistance to high temperature, acid, and 0.3% bile salts. The results of whole-genome sequencing analysis showed that, it is possible that the strain encodes the key enzyme, such as chitinase, carboxylesterases, and lactone hydrolase, that work together to degrade ZEN. In addition, 227 unique genes in this strain are primarily involved in its replication, recombination, repair, and protective mechanisms. In summary, we successfully excavated a ZEN-degrading, genetically distinct strain of Bacillus velezensis that provides a solid foundation for the detoxification of feed and food contamination in the natural environment.

Tailless and filamentous prophages are predominant in marine Vibrio

Although tailed bacteriophages (phages) of the class Caudoviricetes are thought to constitute the most abundant and ecologically relevant group of phages that can integrate their genome into the host chromosome, it is becoming increasingly clear that other prophages are widespread. Here, we show that prophages derived from filamentous and tailless phages with genome sizes below 16 kb make up the majority of prophages in marine bacteria of the genus Vibrio. To estimate prophage prevalence unaffected by database biases, we combined comparative genomics and chemical induction of 58 diverse Vibrio cyclitrophicus isolates, resulting in 107 well-curated prophages. Complemented with computationally predicted prophages, we obtained 1158 prophages from 931 naturally co-existing strains of the family Vibrionaceae. Prophages resembling tailless and filamentous phages predominated, accounting for 80% of all prophages in V. cyclitrophicus and 60% across the Vibrionaceae. In our experimental model, prophages of all three viral realms actively replicated upon induction indicating their ability to transfer to new hosts. Indeed, prophages were rapidly gained and lost, as suggested by variable prophage content between closely related V. cyclitrophicus. Prophages related to filamentous and tailless phages were integrated into only three genomic locations and restored the function of their integration site. Despite their small size, they contained highly diverse accessory genes that may contribute to host fitness, such as phage defense systems. We propose that, like their well-studied tailed equivalent, tailless and filamentous temperate phages are active and highly abundant drivers of host ecology and evolution in marine Vibrio, which have been largely overlooked.

Evolution of Mycobacterium abscessus in the human lung: Cumulative mutations and genomic rearrangement of porin genes in patient isolates

BACKGROUND

Mycobacterium abscessus subspecies massiliense (M. massiliense) is increasingly recognized as an emerging bacterial pathogen, particularly in cystic fibrosis (CF) patients and CF centres' respiratory outbreaks. We characterized genomic and phenotypic changes in 15 serial isolates from two CF patients (1S and 2B) with chronic pulmonary M. massiliense infection leading to death, as well as four isolates from a CF centre outbreak in which patient 2B was the index case.

RESULTS

Comparative genomic analysis revealed the mutations affecting growth rate, metabolism, transport, lipids (loss of glycopeptidolipids), antibiotic susceptibility (macrolides and aminoglycosides resistance), and virulence factors. Mutations in 23S rRNA, mmpL4, porin locus and tetR genes occurred in isolates from both CF patients. Interestingly, we identified two different spontaneous mutation events at the mycobacterial porin locus: a fusion of two tandem porin paralogs in patient 1S and a partial deletion of the first porin paralog in patient 2B. These genomic changes correlated with reduced porin protein expression, diminished 14C-glucose uptake, slower bacterial growth rates, and enhanced TNF-α induction in mycobacteria-infected THP-1 human cells. Porin gene complementation of porin mutants partly restored 14C-glucose uptake, growth rate and TNF-α levels to those of intact porin strains.

CONCLUSIONS

We hypothesize that specific mutations accumulated and maintained over time in M. massiliense, including mutations shared among transmissible strains, collectively lead to more virulent, host adapted lineages in CF patients and other susceptible hosts.

TRS: a method for determining transcript termini from RNAtag-seq sequencing data

In bacteria, determination of the 3' termini of transcripts plays an essential role in regulation of gene expression, affecting the functionality and stability of the transcript. Several experimental approaches were developed to identify the 3' termini of transcripts, however, these were applied only to a limited number of bacteria and growth conditions. Here we present a straightforward approach to identify 3' termini from widely available RNA-seq data without the need for additional experiments. Our approach relies on the observation that the RNAtag-seq sequencing protocol results in overabundance of reads mapped to transcript 3' termini. We present TRS (Termini by Read Starts), a computational pipeline exploiting this property to identify 3' termini in RNAtag-seq data, and show that the identified 3' termini are highly reliable. Since RNAtag-seq data are widely available for many bacteria and growth conditions, our approach paves the way for studying bacterial transcription termination in an unprecedented scope.

CoreDetector: a flexible and efficient program for core-genome alignment of evolutionary diverse genomes

MOTIVATION

Whole genome alignment of eukaryote species remains an important method for the determination of sequence and structural variations and can also be used to ascertain the representative non-redundant core-genome sequence of a population. Many whole genome alignment tools were first developed for the more mature analysis of prokaryote species with few current tools containing the functionality to process larger genomes of eukaryotes as well as genomes of more divergent species. In addition, the functionality of these tools becomes computationally prohibitive due to the significant compute resources needed to handle larger genomes.

RESULTS

In this research, we present CoreDetector, an easy-to-use general-purpose program that can align the core-genome sequences for a range of genome sizes and divergence levels. To illustrate the flexibility of CoreDetector, we conducted alignments of a large set of closely related fungal pathogen and hexaploid wheat cultivar genomes as well as more divergent fly and rodent species genomes. In all cases, compared to existing multiple genome alignment tools, CoreDetector exhibited improved flexibility, efficiency, and competitive accuracy in tested cases.

AVAILABILITY AND IMPLEMENTATION

CoreDetector was developed in the cross platform, and easily deployable, Java language. A packaged pipeline is readily executable in a bash terminal without any external need for Perl or Python environments. Installation, example data, and usage instructions for CoreDetector are freely available from https://github.com/mfruzan/CoreDetector.

Maast: genotyping thousands of microbial strains efficiently

Existing single nucleotide polymorphism (SNP) genotyping algorithms do not scale for species with thousands of sequenced strains, nor do they account for conspecific redundancy. Here we present a bioinformatics tool, Maast, which empowers population genetic meta-analysis of microbes at an unrivaled scale. Maast implements a novel algorithm to heuristically identify a minimal set of diverse conspecific genomes, then constructs a reliable SNP panel for each species, and enables rapid and accurate genotyping using a hybrid of whole-genome alignment and k-mer exact matching. We demonstrate Maast's utility by genotyping thousands of Helicobacter pylori strains and tracking SARS-CoV-2 diversification.

Comparative sequence analysis of pPATH pathogenicity plasmids in Pantoea agglomerans gall-forming bacteria

Acquisition of the pathogenicity plasmid pPATH that encodes a type III secretion system (T3SS) and effectors (T3Es) has likely led to the transition of a non-pathogenic bacterium into the tumorigenic pathogen Pantoea agglomerans. P. agglomerans pv. gypsophilae (Pag) forms galls on gypsophila (Gypsophila paniculata) and triggers immunity on sugar beet (Beta vulgaris), while P. agglomerans pv. betae (Pab) causes galls on both gypsophila and sugar beet. Draft sequences of the Pag and Pab genomes were previously generated using the MiSeq Illumina technology and used to determine partial T3E inventories of Pab and Pag. Here, we fully assembled the Pab and Pag genomes following sequencing with PacBio technology and carried out a comparative sequence analysis of the Pab and Pag pathogenicity plasmids pPATHpag and pPATHpab. Assembly of Pab and Pag genomes revealed a ~4 Mbp chromosome with a 55% GC content, and three and four plasmids in Pab and Pag, respectively. pPATHpag and pPATHpab share 97% identity within a 74% coverage, and a similar GC content (51%); they are ~156 kb and ~131 kb in size and consist of 198 and 155 coding sequences (CDSs), respectively. In both plasmids, we confirmed the presence of highly similar gene clusters encoding a T3SS, as well as auxin and cytokinins biosynthetic enzymes. Three putative novel T3Es were identified in Pab and one in Pag. Among T3SS-associated proteins encoded by Pag and Pab, we identified two novel chaperons of the ShcV and CesT families that are present in both pathovars with high similarity. We also identified insertion sequences (ISs) and transposons (Tns) that may have contributed to the evolution of the two pathovars. These include seven shared IS elements, and three ISs and two transposons unique to Pab. Finally, comparative sequence analysis revealed plasmid regions and CDSs that are present only in pPATHpab or in pPATHpag. The high similarity and common features of the pPATH plasmids support the hypothesis that the two strains recently evolved into host-specific pathogens.

Phylotranscriptomics supports numerous polyploidization events and phylogenetic relationships in Nicotiana

Introduction

Nicotiana L. (Solanaceae) is of great scientific and economic importance, and polyploidization has been pivotal in shaping this genus. Despite many previous studies on the Nicotiana phylogenetic relationship and hybridization, evidence from whole genome data is still lacking.

Methods

In this study, we obtained 995 low-copy genes and plastid transcript fragments from the transcriptome datasets of 26 Nicotiana species, including all sections. We reconstructed the phylogenetic relationship and phylogenetic network of diploid species.

Results

The incongruence among gene trees showed that the formation of N. sylvestris involved incomplete lineage sorting. The nuclear-plastid discordance and nuclear introgression absence indicated that organelle capture from section Trigonophyllae was involved in forming section Petunioides. Furthermore, we analyzed the evolutionary origin of polyploid species and dated the time of hybridization events based on the analysis of PhyloNet, sequence similarity search, and phylogeny of subgenome approaches. Our results highly evidenced the hybrid origins of five polyploid sections, including sections Nicotiana, Repandae, Rusticae, Polydicliae, and Suaveolentes. Notably, we provide novel insights into the hybridization event of section Polydicliae and Suaveolentes. The section Polydicliae formed from a single hybridization event between maternal progenitor N. attenuata and paternal progenitor N. undulata; the N. sylvestris (paternal progenitor) and the N. glauca (maternal progenitor) were involved in the formation of section Suaveolentes.

Discussion

This study represents the first exploration of Nicotiana polyploidization events and phylogenetic relationships using the high-throughput RNA-seq approach. It will provide guidance for further studies in molecular systematics, population genetics, and ecological adaption studies in Nicotiana and other related species.

PanGraph: scalable bacterial pan-genome graph construction

The genomic diversity of microbes is commonly parameterized as SNPs relative to a reference genome of a well-characterized, but arbitrary, isolate. However, any reference genome contains only a fraction of the microbial pangenome, the total set of genes observed in a given species. Reference-based approaches are thus blind to the dynamics of the accessory genome, as well as variation within gene order and copy number. With the widespread usage of long-read sequencing, the number of high-quality, complete genome assemblies has increased dramatically. In addition to pangenomic approaches that focus on the variation in the sets of genes present in different genomes, complete assemblies allow investigations of the evolution of genome structure and gene order. This latter problem, however, is computationally demanding with few tools available that shed light on these dynamics. Here, we present PanGraph, a Julia-based library and command line interface for aligning whole genomes into a graph. Each genome is represented as a path along vertices, which in turn encapsulate homologous multiple sequence alignments. The resultant data structure succinctly summarizes population-level nucleotide and structural polymorphisms and can be exported into several common formats for either downstream analysis or immediate visualization.

Host Association and Spatial Proximity Shape but Do Not Constrain Population Structure in the Mutualistic Symbiont Xenorhabdus bovienii

To what extent are generalist species cohesive evolutionary units rather than a compilation of recently diverged lineages? We examine this question in the context of host specificity and geographic structure in the insect pathogen and nematode mutualist Xenorhabdus bovienii. This bacterial species partners with multiple nematode species across two clades in the genus Steinernema. We sequenced the genomes of 42 X. bovienii strains isolated from four different nematode species and three field sites within a 240-km² region and compared them to globally available reference genomes. We hypothesized that X. bovienii would comprise several host-specific lineages, such that bacterial and nematode phylogenies would be largely congruent. Alternatively, we hypothesized that spatial proximity might be a dominant signal, as increasing geographic distance might lower shared selective pressures and opportunities for gene flow. We found partial support for both hypotheses. Isolates clustered largely by nematode host species but did not strictly match the nematode phylogeny, indicating that shifts in symbiont associations across nematode species and clades have occurred. Furthermore, both genetic similarity and gene flow decreased with geographic distance across nematode species, suggesting differentiation and constraints on gene flow across both factors, although no absolute barriers to gene flow were observed across the regional isolates. Several genes associated with biotic interactions were found to be undergoing selective sweeps within this regional population. The interactions included several insect toxins and genes implicated in microbial competition. Thus, gene flow maintains cohesiveness across host associations in this symbiont and may facilitate adaptive responses to a multipartite selective environment. IMPORTANCE Microbial populations and species are notoriously hard to delineate. We used a population genomics approach to examine the population structure and the spatial scale of gene flow in Xenorhabdus bovienii, an intriguing species that is both a specialized mutualistic symbiont of nematodes and a broadly virulent insect pathogen. We found a strong signature of nematode host association, as well as evidence for gene flow connecting isolates associated with different nematode host species and collected from distinct study sites. Furthermore, we saw signatures of selective sweeps for genes involved with nematode host associations, insect pathogenicity, and microbial competition. Thus, X. bovienii exemplifies the growing consensus that recombination not only maintains cohesion but can also allow the spread of niche-beneficial alleles.

Two major chromosome evolution events with unrivaled conserved gene content in pomegranate

Pomegranate has a unique evolutionary history given that different cultivars have eight or nine bivalent chromosomes with possible crossability between the two classes. Therefore, it is important to study chromosome evolution in pomegranate to understand the dynamics of its population. Here, we de novo assembled the Azerbaijani cultivar "Azerbaijan guloyshasi" (AG2017; 2n = 16) and re-sequenced six cultivars to track the evolution of pomegranate and to compare it with previously published de novo assembled and re-sequenced cultivars. High synteny was observed between AG2017, Bhagawa (2n = 16), Tunisia (2n = 16), and Dabenzi (2n = 18), but these four cultivars diverged from the cultivar Taishanhong (2n = 18) with several rearrangements indicating the presence of two major chromosome evolution events. Major presence/absence variations were not observed as >99% of the five genomes aligned across the cultivars, while >99% of the pan-genic content was represented by Tunisia and Taishanhong only. We also revisited the divergence between soft- and hard-seeded cultivars with less structured population genomic data, compared to previous studies, to refine the selected genomic regions and detect global migration routes for pomegranate. We reported a unique admixture between soft- and hard-seeded cultivars that can be exploited to improve the diversity, quality, and adaptability of local pomegranate varieties around the world. Our study adds body knowledge to understanding the evolution of the pomegranate genome and its implications for the population structure of global pomegranate diversity, as well as planning breeding programs aiming to develop improved cultivars.

A Global Survey of Hypervirulent Aeromonas hydrophila (vAh) Identified vAh Strains in the Lower Mekong River Basin and Diverse Opportunistic Pathogens from Farmed Fish and Other Environmental Sources

Hypervirulent Aeromonas hydrophila (vAh) has emerged as the etiologic agent of epidemic outbreaks of motile Aeromonas septicemia (MAS) in high-density aquaculture of farmed carp in China and catfish in the United States, which has caused millions of tons of lost fish. We conducted a global survey to better understand the evolution, geographical distribution, and phylogeny of vAh. Aeromonas isolates were isolated from fish that showed clinical symptoms of MAS, and pure cultures were screened for the ability to utilize myo-inositol as the sole carbon source. A total of 113 myo-inositol-utilizing bacterial strains were included in this study, including additional strains obtained from previously published culture collections. Based on a gyrB phylogeny, this collection included 66 A. hydrophila isolates, 48 of which were vAh. This collection also included five new vAh isolates from diseased Pangas catfish (Pangasius pangasius) and striped catfish (Pangasianodon hypophthalmus) obtained in Cambodia and Vietnam, respectively. Genome sequences were generated from representative vAh and non-vAh isolates to evaluate the potential for lateral genetic transfer of the myo-inositol catabolism pathway. Phylogenetic analyses of each of the nine genes required for myo-inositol utilization revealed the close affiliation of vAh strains regardless of geographic origin and suggested lateral genetic transfer of this catabolic pathway from an Enterobacter species. Prediction of virulence factors was conducted to determine differences between vAh and non-vAh strains in terms of virulence and secretion systems. Core genome phylogenetic analyses on vAh isolates and Aeromonas spp. disease isolates (55 in total) were conducted to evaluate the evolutionary relationships among vAh and other Aeromonas sp. isolates, which supported the clonal nature of vAh isolates. IMPORTANCE This global survey of vAh brought together scientists that study fish disease to evaluate the evolution, geographical distribution, phylogeny, and hosts of vAh and other Aeromonas sp. isolates. In addition to vAh isolates from China and the United States, four new vAh isolates were isolated from the lower Mekong River basin in Cambodia and Vietnam, indicating the significant threat of vAh to modern aquaculture and the need for improved biosecurity to prevent vAh spread.

Mitochondrial genome characteristics and phylogenetic analysis of the medicinal and edible plant Mesona chinensis Benth

Mesona chinensis Benth (MCB) (or Platostoma palustre or Platostoma chinense) is an important edible and medicinal plant in China. However, the mitochondrial genome (mitogenome, or mtDNA) of MCB has not been characterized or reported yet. In this study, we first sequenced and characterized the complete mitogenome of MCB. The MCB mitogenome was 494,599 bp in length and encoded 59 genes containing 37 protein-coding genes (PCGs), 19 tRNAs, and 3 rRNAs. Gene transfer analysis revealed that a total of 12 transfer segments with more than 93% identity (total length of 25,427 bp) were detected in the MCB mitogenome. Simple sequence repeats (SSR) analysis showed that 212 simple sequence repeats (SSR) were identified. Repeat sequence analysis revealed 305 repeat sequences (158 forward and 147 palindromic repeats) ranging from 30 bp to 48,383 bp and the 30-39 bp repeats were the majority type. Relative synonymous codon usage (RSCU) analysis uncovered that in total, 9,947 codons were encoding the protein-coding genes (PCGs). Serine (909, 9.1%) and leucine (879, 8.8%) were the two most abundant amino acids, while terminator (32, .3%) was the least abundant amino acid. Ka/Ks analysis indicated that almost all genes were subject to purification selection, except ccmB. Analysis of Lamiaceae mitogenomes constitution revealed that atpB and atpE were unique to the Rotheca serrata and Salvia miltiorrhiza mitogenomes. mttB gene loss was unique to the Boea hygrometrica mitogenome. The core fragments of the Lamiaceae mitogenomes harbored a higher GC content than the specific and variable fragments. In addition, phylogenetic analysis revealed that MCB was closely related to Salvia miltiorrhiza based on the mitogenomes. The current study provided valuable genomic resources for understanding and utilizing this important medicinal plant in the future.

PRAWNS: compact pan-genomic features for whole-genome population genomics

MOTIVATION

Scientists seeking to understand the genomic basis of bacterial phenotypes, such as antibiotic resistance, today have access to an unprecedented number of complete and nearly complete genomes. Making sense of these data requires computational tools able to perform multiple-genome comparisons efficiently, yet currently available tools cannot scale beyond several tens of genomes.

RESULTS

We describe PRAWNS, an efficient and scalable tool for multiple-genome analysis. PRAWNS defines a concise set of genomic features (metablocks), as well as pairwise relationships between them, which can be used as a basis for large-scale genotype-phenotype association studies. We demonstrate the effectiveness of PRAWNS by identifying genomic regions associated with antibiotic resistance in Acinetobacter baumannii.

AVAILABILITY AND IMPLEMENTATION

PRAWNS is implemented in C++ and Python3, licensed under the GPLv3 license, and freely downloadable from GitHub (https://github.com/KiranJavkar/PRAWNS.git).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A new haplotype-resolved turkey genome to enable turkey genetics and genomics research

BACKGROUND

The domesticated turkey (Meleagris gallopavo) is a species of significant agricultural importance and is the second largest contributor, behind broiler chickens, to world poultry meat production. The previous genome is of draft quality and partly based on the chicken (Gallus gallus) genome. A high-quality reference genome of M. gallopavo is essential for turkey genomics and genetics research and the breeding industry.

RESULTS

By adopting the trio-binning approach, we were able to assemble a high-quality chromosome-level F1 assembly and 2 parental haplotype assemblies, leveraging long-read technologies and genome-wide chromatin interaction data (Hi-C). From a total of 40 chromosomes (2n = 80), we captured 35 chromosomes in a single scaffold, showing much improved genome completeness and continuity compared to the old assembly build. The 3 assemblies are of higher quality than the previous draft quality assembly and comparable to the chicken assemblies (GRCg7) shown by the largest contig N50 (26.6 Mb) and comparable BUSCO gene set completeness scores (96-97%). Comparative analyses confirm a previously identified large inversion of around 19 Mbp on the Z chromosome not found in other Galliformes. Structural variation between the parent haplotypes was identified, which poses potential new target genes for breeding.

CONCLUSIONS

We contribute a new high-quality turkey genome at the chromosome level, benefiting turkey genetics and other avian genomics research as well as the turkey breeding industry.

Genomes of Vibrio metoecus co-isolated with Vibrio cholerae extend our understanding of differences between these closely related species

BACKGROUND

Vibrio cholerae, the causative agent of cholera, is a well-studied species, whereas Vibrio metoecus is a recently described close relative that is also associated with human infections. The availability of V. metoecus genomes provides further insight into its genetic differences from V. cholerae. Additionally, both species have been co-isolated from a cholera-free brackish coastal pond and have been suggested to interact with each other by horizontal gene transfer (HGT).

RESULTS

The genomes of 17 strains from each species were sequenced. All strains share a large core genome (2675 gene families) and very few genes are unique to each species (< 3% of the pan-genome of both species). This led to the identification of potential molecular markers-for nitrite reduction, as well as peptidase and rhodanese activities-to further distinguish V. metoecus from V. cholerae. Interspecies HGT events were inferred in 21% of the core genes and 45% of the accessory genes. A directional bias in gene transfer events was found in the core genome, where V. metoecus was a recipient of three times (75%) more genes from V. cholerae than it was a donor (25%).

CONCLUSION

V. metoecus was misclassified as an atypical variant of V. cholerae due to their resemblance in a majority of biochemical characteristics. More distinguishing phenotypic assays can be developed based on the discovery of potential gene markers to avoid any future misclassifications. Furthermore, differences in relative abundance or seasonality were observed between the species and could contribute to the bias in directionality of HGT.

Whole genome resequencing and comparative genome analysis of three Puccinia striiformis f. sp. tritici pathotypes prevalent in India

Stripe rust disease of wheat, caused by Puccinia striiformis f. sp. tritici, (Pst) is one of the most serious diseases of wheat worldwide. In India, virulent stripe rust races have been constantly evolving in the North-Western Plains Zone leading to the failure of some of the most widely grown resistant varieties in the region. With the goal of studying the recent evolution of virulent races in this region, we conducted whole-genome re-sequencing of three prevalent Indian Pst pathotypes Pst46S119, Pst78S84 and Pst110S119. We assembled 58.62, 58.33 and 55.78 Mb of Pst110S119, Pst46S119 and Pst78S84 genome, respectively and found that pathotypes were highly heterozygous. Comparative phylogenetic analysis indicated the recent evolution of pathotypes Pst110S119 and Pst78S84 from Pst46S119. Pathogenicity-related genes classes (CAZyme, proteases, effectors, and secretome proteins) were identified and found to be under positive selection. Higher rate of gene families expansion were also observed in the three pathotypes. A strong association between the effector genes and transposable elements may be the source of the rapid evolution of these strains. Phylogenetic analysis differentiated the Indian races in this study from other known United States, European, African, and Asian races. Diagnostic markers developed for the identification of three Pst pathotypes will help tracking of yellow rust at farmers field and strategizing resistance gene deployment.

Repeat infections with chlamydia in women may be more transcriptionally active with lower responses from some immune genes

Chlamydia trachomatis, the most common bacterial sexually transmitted infection worldwide, is responsible for considerable health burden due to its significant sequelae. There are growing concerns about chlamydial treatment and management due to widely documented increasing burden of repeat infections. In the current study, a cohort study design of 305 women with urogenital chlamydial infections demonstrated that 11.8% of women experienced repeat infections after treatment with azithromycin. The chlamydial DNA load measured by quantitative PCR was higher in women who experienced a repeat infection (p = 0.0097) and repeat infection was associated with sexual contact. There was no genomic or phenotypic evidence of azithromycin resistance within the chlamydial isolates. During repeat infection, or repeat positive tests during follow up, vaginal chlamydial gene expression (ompA, euo, omcB, htrA, trpAB) was markedly higher compared to baseline, and two of the selected immune genes analyzed had significantly lower expression at the time of repeat infection. Overall, there are two implications of these results. The results could be generalized to all recent infections, or repeat positive events, and indicate that chlamydial infections are have higher transcriptional activity of select genes early in the infection in women. Alternatively, after azithromycin treatment, repeat infections of Chlamydia may be more transcriptionally active at certain genes, and there may be post-treatment immunological alterations that interplay into repeat exposures establishing an active infection. The potential that recent infections may involve a higher level of activity from the organism may have implications for management by more regular testing of the most at risk women to reduce the risk of sequelae.

Phylodynamics of a regional SARS-CoV-2 rapid spreading event in Colorado in late 2020

Since the initial reported discovery of SARS-CoV-2 in late 2019, genomic surveillance has been an important tool to understand its transmission and evolution. Here, we sought to describe the underlying regional phylodynamics before and during a rapid spreading event that was documented by surveillance protocols of the United States Air Force Academy (USAFA) in late October-November of 2020. We used replicate long-read sequencing on Colorado SARS-CoV-2 genomes collected July through November 2020 at the University of Colorado Anschutz Medical campus in Aurora and the United States Air Force Academy in Colorado Springs. Replicate sequencing allowed rigorous validation of variation and placement in a phylogenetic relatedness network. We focus on describing the phylodynamics of a lineage that likely originated in the local Colorado Springs community and expanded rapidly over the course of two months in an outbreak within the well-controlled environment of the United States Air Force Academy. Divergence estimates from sampling dates indicate that the SARS-CoV-2 lineage associated with this rapid expansion event originated in late October 2020. These results are in agreement with transmission pathways inferred by the United States Air Force Academy, and provide a window into the evolutionary process and transmission dynamics of a potentially dangerous but ultimately contained variant.

Dynamic Subspecies Population Structure of Vibrio cholerae in Dhaka, Bangladesh

Cholera has been endemic to the Ganges Delta for centuries. Although the causative agent, Vibrio cholerae, is autochthonous to coastal and brackish water, cholera occurs continually in Dhaka, the inland capital city of Bangladesh which is surrounded by fresh water. Despite the persistence of this problem, little is known about the environmental abundance and distribution of lineages of V. cholerae, the most important being the pandemic generating (PG) lineage consisting mostly of serogroup O1 strains. To understand spatial and temporal dynamics of PG lineage and other lineages belonging to the V. cholerae species in surface water in and around Dhaka City, we used qPCR and high-throughput amplicon sequencing. Seven different freshwater sites across Dhaka were investigated for six consecutive months, and physiochemical parameters were measured in situ. Total abundance of V. cholerae was found to be relatively stable throughout the 6-month sampling period, with 2 × 10⁵ to 4 × 10⁵ genome copies/L at six sites and around 5 × 10⁵ genome copies/L at the site located in the most densely populated part of Dhaka City. PG O1 V. cholerae was present in high abundance during the entire sampling period and composed between 24 and 92% of the total V. cholerae population, only showing occasional but sudden reductions in abundance. In instances where PG O1 lost its dominance, other lineages underwent a rapid expansion while the size of the total V. cholerae population remained almost unchanged. Intraspecies richness of V. cholerae was positively correlated with salinity, conductivity, and total dissolved solids (TDS), while it was negatively correlated with dissolved oxygen (DO) concentration in water. Interestingly, negative correlation was observed specifically between PG O1 and salinity, even though the changes in this variable were minor (0-0.8 ppt). Observations in this study suggest that at the subspecies level, population composition of naturally occurring V. cholerae can be influenced by fluctuations in environmental factors, which can lead to altered competition dynamics among the lineages.

The effect of a temperature-sensitive prophage on the evolution of virulence in an opportunistic bacterial pathogen

Viruses are key actors of ecosystems and have major impacts on global biogeochemical cycles. Prophages deserve particular attention as they are ubiquitous in bacterial genomes and can enter a lytic cycle when triggered by environmental conditions. We explored how temperature affects the interactions between prophages and other biological levels using an opportunistic pathogen, the bacterium Serratia marcescens, which harbours several prophages and that had undergone an evolution experiment under several temperature regimes. We found that the release of one of the prophages was temperature-sensitive and malleable to evolutionary changes. We further discovered that the virulence of the bacterium in an insect model also evolved and was positively correlated with phage release rates. We determined through analysis of genetic and epigenetic data that changes in the bacterial outer cell wall structure possibly explain this phenomenon. We hypothezise that the temperature-dependent phage release rate acted as a selection pressure on S. marcescens and that it resulted in modified bacterial virulence in the insect host. Our study system illustrates how viruses can mediate the influence of abiotic environmental changes to other biological levels and thus be involved in ecosystem feedback loops.

Vibrio tarriae sp. nov., a novel member of the Cholerae clade

A number of bacteria with close resemblance to Vibrio cholerae have been isolated over the years by the Centres for Disease Control and Prevention (CDC), which could not be assigned a proper taxonomic designation on the basis of the results from preliminary identification methods. Nine such isolates have been found to share 16S rRNA gene identity exceeding 99 % with V. cholerae, yet DNA–DNA hybridization (60.4–62.1 %) and average nucleotide identity values (94.4–95.1 %) were below the species cut-off, indicating a potentially novel species. Phylogenetic analysis of core genomes places this group of isolates in a monophyletic clade, within the ‘Cholerae clade’, but distinct from any other species. Extensive phenotypic characterization reveals unique biochemical properties that distinguish this novel species from V. cholerae . Comparative genomic analysis reveals a unique set of siderophore genes, indicating that iron acquisition strategies could be vital for the divergence of the novel species from a common ancestor with V. cholerae . On the basis of the genetic, phylogenetic and phenotypic differences observed, we propose that these isolates represent a novel species of the genus Vibrio , for which the name Vibrio tarriae sp. nov. is proposed. Strain 2521-89 T (= DSM 112461=CCUG 75318), isolated from lake water, is the type strain.

Roles of Species-Specific Legumains in Pathogenicity of the Pinewood Nematode Bursaphelenchus xylophilus

Peptidases are very important to parasites, which have central roles in parasite biology and pathogenesis. In this study, by comparative genome analysis, genome-wide peptidase diversities among plant-parasitic nematodes are estimated. We find that genes encoding cysteine peptidases in family C13 (legumain) are significantly abundant in pine wood nematodes Bursaphelenchus genomes, compared to those in other plant-parasitic nematodes. By phylogenetic analysis, a clade of B. xylophilus-specific legumain is identified. RT-qPCR detection shows that these genes are highly expressed at early stage during the nematode infection process. Utilizing transgene technology, cDNAs of three species-specific legumain were introduced into the Arabidopsis γvpe mutant. Functional complementation assay shows that these B. xylophilus legumains can fully complement the activity of Arabidopsis γVPE to mediate plant cell death triggered by the fungal toxin FB1. Secretory activities of these legumains are experimentally validated. By comparative transcriptome analysis, genes involved in plant cell death mediated by legumains are identified, which enrich in GO terms related to ubiquitin protein transferase activity in category molecular function, and response to stimuli in category biological process. Our results suggest that B. xylophilu-specific legumains have potential as effectors to be involved in nematode-plant interaction and can be related to host cell death.

Multiple genome alignment in the telomere-to-telomere assembly era

With the arrival of telomere-to-telomere (T2T) assemblies of the human genome comes the computational challenge of efficiently and accurately constructing multiple genome alignments at an unprecedented scale. By identifying nucleotides across genomes which share a common ancestor, multiple genome alignments commonly serve as the bedrock for comparative genomics studies. In this review, we provide an overview of the algorithmic template that most multiple genome alignment methods follow. We also discuss prospective areas of improvement of multiple genome alignment for keeping up with continuously arriving high-quality T2T assembled genomes and for unlocking clinically-relevant insights.

Insertion sequences and other mobile elements associated with antibiotic resistance genes in Enterococcus isolates from an inpatient with prolonged bacteraemia

Insertion sequences (ISs) and other transposable elements are associated with the mobilization of antibiotic resistance determinants and the modulation of pathogenic characteristics. In this work, we aimed to investigate the association between ISs and antibiotic resistance genes, and their role in the dissemination and modification of the antibiotic-resistant phenotype. To that end, we leveraged fully resolved Enterococcus faecium and Enterococcus faecalis genomes of isolates collected over 5 days from an inpatient with prolonged bacteraemia. Isolates from both species harboured similar IS family content but showed significant species-dependent differences in copy number and arrangements of ISs throughout their replicons. Here, we describe two inter-specific IS-mediated recombination events and IS-mediated excision events in plasmids of E. faecium isolates. We also characterize a novel arrangement of the ISs in a Tn1546-like transposon in E. faecalis isolates likely implicated in a vancomycin genotype–phenotype discrepancy. Furthermore, an extended analysis revealed a novel association between daptomycin resistance mutations in liaSR genes and a putative composite transposon in E. faecium, offering a new paradigm for the study of daptomycin resistance and novel insights into its dissemination. In conclusion, our study highlights the role ISs and other transposable elements play in the rapid adaptation and response to clinically relevant stresses such as aggressive antibiotic treatment in enterococci.

Genome-wide analysis of citrus TCP transcription factors and their responses to abiotic stresses

BACKGROUND

Citrus is one of the most important fruit crops in the world, and it is worthy to conduct more research on artificially controlling citrus plant growth and development to adapt to different cultivation patterns and environmental conditions. The plant-specific TEOSINTE BRANCHED1, CYCOLOIDEA, and PROLIFERATING CELL FACTORS (TCP) transcription factors are crucial regulators controlling plant growth and development, as well as responding to abiotic stresses. However, the information about citrus TCP transcription factors remains unclear.

RESULTS

In this study, twenty putative TCP genes (CsTCPs) with the TCP domain were explored from Citrus sinensis genome, of which eleven (CsTCP3, - 4, - 5, - 6, - 10, - 11, - 15, - 16, - 18, - 19, - 20), five (CsTCP1, - 2, - 7, - 9, - 13), and four genes (CsTCP8, - 12, - 14, - 17) were unevenly distributed on chromosomes and divided into three subclades. Cis-acting element analysis indicated that most CsTCPs contained many phytohormone- and environment-responsive elements in promoter regions. All of CsTCPs were predominantly expressed in vegetative tissues or organs (stem, leaf, thorn, and bud) instead of reproductive tissues or organs (flower, fruit, and seed). Combined with collinearity analysis, CsTCP3, CsTCP9, and CsTCP13 may take part in leaf development; CsTCP12 and CsTCP14 may function in shoot branching, leaf development, or thorn development; CsTCP15 may participate in the development of stem, leaf, or thorn. In mature leaf, transcript levels of two CsTCPs (CsTCP19, - 20) were significantly increased while transcript levels of eight CsTCPs (CsTCP2, - 5, - 6, - 7, - 8, - 9, - 10, - 13) were significantly decreased by shading; except for two CsTCPs (CsTCP11, - 19), CsTCPs' transcript levels were significantly influenced by low temperature; moreover, transcript levels of two CsTCPs (CsTCP11, - 12) were significantly increased while five CsTCPs' (CsTCP14, - 16, - 18, - 19, - 20) transcript levels were significantly reduced by drought.

CONCLUSIONS

This study provides significant clues for research on roles of CsTCPs in regulating citrus plant growth and development, as well as responding to abiotic stresses.

Sequence-based pangenomic core detection

One of the most basic kinds of analysis to be performed on a pangenome is the detection of its core, i.e., the information shared among all members. Pangenomic core detection is classically done on the gene level and many tools focus exclusively on core detection in prokaryotes. Here, we present a new method for sequence-based pangenomic core detection. Our model generalizes from a strict core definition allowing us to flexibly determine suitable core properties depending on the research question and the dataset under consideration. We propose an algorithm based on a colored de Bruijn graph that runs in linear time with respect to the number of k-mers in the graph. An implementation of our method is called Corer. Because of the usage of a colored de Bruijn graph, it works alignment-free, is provided with a small memory footprint, and accepts as input assembled genomes as well as sequencing reads.

•

Pangenomic core detection for large collections of prokaryotes or higher eukaryotes

•

Whole-genome analysis with assemblies or even read data as input

•

Alignment-free, linear time algorithm with small memory footprint

•

Variation tolerance and quorum for flexible core detection

A large-scale population based organelle pan-genomes construction and phylogeny analysis reveal the genetic diversity and the evolutionary origins of chloroplast and mitochondrion in Brassica napus L

BACKGROUND

Allotetraploid oilseed rape (Brassica napus L.) is an important worldwide oil-producing crop. The origin of rapeseed is still undetermined due to the lack of wild resources. Despite certain genetic architecture and phylogenetic studies have been done focus on large group of Brassica nuclear genomes, the organelle genomes information under global pattern is largely unknown, which provide unique material for phylogenetic studies of B. napus. Here, based on de novo assemblies of 1,579 B. napus accessions collected globally, we constructed the chloroplast and mitochondrial pan-genomes of B. napus, and investigated the genetic diversity, phylogenetic relationships of B. napus, B. rapa and B. oleracea.

RESULTS

Based on mitotype-specific markers and mitotype-variant ORFs, four main cytoplasmic haplotypes were identified in our groups corresponding the nap, pol, ole, and cam mitotypes, among which the structure of chloroplast genomes was more conserved without any rearrangement than mitochondrial genomes. A total of 2,092 variants were detected in chloroplast genomes, whereas only 326 in mitochondrial genomes, indicating that chloroplast genomes exhibited a higher level of single-base polymorphism than mitochondrial genomes. Based on whole-genome variants diversity analysis, eleven genetic difference regions among different cytoplasmic haplotypes were identified on chloroplast genomes. The phylogenetic tree incorporating accessions of the B. rapa, B. oleracea, natural and synthetic populations of B. napus revealed multiple origins of B. napus cytoplasm. The cam-type and pol-type were both derived from B. rapa, while the ole-type was originated from B. oleracea. Notably, the nap-type cytoplasm was identified in both the B. rapa population and the synthetic B. napus, suggesting that B. rapa might be the maternal ancestor of nap-type B. napus.

CONCLUSIONS

The phylogenetic results provide novel insights into the organelle genomic evolution of Brassica species. The natural rapeseeds contained at least four cytoplastic haplotypes, of which the predominant nap-type might be originated from B. rapa. Besides, the organelle pan-genomes and the overall variation data offered useful resources for analysis of cytoplasmic inheritance related agronomical important traits of rapeseed, which can substantially facilitate the cultivation and improvement of rapeseed varieties.

Anopheles mosquitoes reveal new principles of 3D genome organization in insects

Chromosomes are hierarchically folded within cell nuclei into territories, domains and subdomains, but the functional importance and evolutionary dynamics of these hierarchies are poorly defined. Here, we comprehensively profile genome organizations of five Anopheles mosquito species and show how different levels of chromatin architecture influence each other. Patterns observed on Hi-C maps are associated with known cytological structures, epigenetic profiles, and gene expression levels. Evolutionary analysis reveals conservation of chromatin architecture within synteny blocks for tens of millions of years and enrichment of synteny breakpoints in regions with increased genomic insulation. However, in-depth analysis shows a confounding effect of gene density on both insulation and distribution of synteny breakpoints, suggesting limited causal relationship between breakpoints and regions with increased genomic insulation. At the level of individual loci, we identify specific, extremely long-ranged looping interactions, conserved for ~100 million years. We demonstrate that the mechanisms underlying these looping contacts differ from previously described Polycomb-dependent interactions and clustering of active chromatin.

Pathogenomes and variations in Shiga toxin production among geographically distinct clones of Escherichia coli O113:H21

Infections with globally disseminated Shiga toxin-producing Escherichia coli (STEC) of the O113:H21 serotype can progress to severe clinical complications, such as hemolytic uremic syndrome (HUS). Two phylogeographically distinct clonal complexes have been established by multi locus sequence typing (MLST). Infections with ST-820 isolates circulating exclusively in Australia have caused severe human disease, such as HUS. Conversely, ST-223 isolates prevalent in the US and outside Australia seem to rarely cause severe human disease but are frequent contaminants. Following a genomic epidemiology approach, we wanted to gain insights into the underlying cause for this disparity. We examined the plasticity in the genome make-up and Shiga toxin production in a collection of 20 ST-820 and ST-223 strains isolated from produce, the bovine reservoir, and clinical cases. STEC are notorious for assembly into fragmented draft sequences when using short-read sequencing technologies due to the extensive and partly homologous phage complement. The application of long-read technology (LRT) sequencing yielded closed reference chromosomes and plasmids for two representative ST-820 and ST-223 strains. The established high-resolution framework, based on whole genome alignments, single nucleotide polymorphism (SNP)-typing and MLST, includes the chromosomes and plasmids of other publicly available O113:H21 sequences and allowed us to refine the phylogeographical boundaries of ST-820 and ST-223 complex isolates and to further identify a historic non-shigatoxigenic strain from Mexico as a quasi-intermediate. Plasmid comparison revealed strong correlations between the strains' featured pO113 plasmid genotypes and chromosomally inferred ST, which suggests coevolution of the chromosome and virulence plasmids. Our pathogenicity assessment revealed statistically significant differences in the Stx_2a-production capabilities of ST-820 as compared to ST-223 strains under RecA-induced Stx phage mobilization, a condition that mimics Stx-phage induction. These observations suggest that ST-820 strains may confer an increased pathogenic potential in line with the strain-associated epidemiological metadata. Still, some of the tested ST-223 cultures sourced from contaminated produce or the bovine reservoir also produced Stx at levels comparable to those of ST-820 isolates, which calls for awareness and for continued surveillance of this lineage.

Complete Genome Sequencing and Comparative Analysis of the Clinically-Derived Apiotrichum mycotoxinivorans Strain GMU1709

Over the past decade, Apiotrichum mycotoxinivorans has been recognized globally as a source of opportunistic infections. It is a yeast-like fungus, and its association as an uncommon pulmonary pathogen with cystic fibrosis patients has been previously reported. Immunocompromised patients are at the highest risk of A. mycotoxinivorans infections. Therefore, to investigate the genetic basis for the pathogenicity of A. mycotoxinivorans, we performed whole-genome sequencing and comparative genomic analysis of A. mycotoxinivorans GMU1709 that was isolated from sputum specimens of a pneumonia patient receiving cardiac repair surgery. The assembly of Oxford Nanopore reads from the GMU1709 strain and its subsequent correction using Illumina paired-end reads yielded a high-quality complete genome with a genome size of 30.5 Mb in length, which comprised six chromosomes and one mitochondrion. Subsequently, 8,066 protein-coding genes were predicted based on multiple pieces of evidence, including transcriptomes. Phylogenomic analysis indicated that A. mycotoxinivorans exhibited the closest evolutionary affinity to A. veenhuisii, and both the A. mycotoxinivorans strains and the formerly Trichosporon cutaneum ACCC 20271 strain occupied the same phylogenetic position. Further comparative analysis supported that the ACCC 20271 strain belonged to A. mycotoxinivorans. Comparisons of three A. mycotoxinivorans strains indicated that the differences between clinical and non-clinical strains in pathogenicity and drug resistance may be little or none. Based on the comparisons with strains of other species in the Trichosporonaceae family, we identified potential key genetic factors associated with A. mycotoxinivorans infection or pathogenicity. In addition, we also deduced that A. mycotoxinivorans had great potential to inactivate some antibiotics (e.g., tetracycline), which may affect the efficacy of these drugs in co-infection. In general, our analyses provide a better understanding of the classification and phylogeny of the Trichosporonaceae family, uncover the underlying genetic basis of A. mycotoxinivorans infections and associated drug resistance, and provide clues into potential targets for further research and the therapeutic intervention of infections.

Genome-Wide Reconstruction of Rediploidization Following Autopolyploidization across One Hundred Million Years of Salmonid Evolution

The long-term evolutionary impacts of whole-genome duplication (WGD) are strongly influenced by the ensuing rediploidization process. Following autopolyploidization, rediploidization involves a transition from tetraploid to diploid meiotic pairing, allowing duplicated genes (ohnologs) to diverge genetically and functionally. Our understanding of autopolyploid rediploidization has been informed by a WGD event ancestral to salmonid fishes, where large genomic regions are characterized by temporally delayed rediploidization, allowing lineage-specific ohnolog sequence divergence in the major salmonid clades. Here, we investigate the long-term outcomes of autopolyploid rediploidization at genome-wide resolution, exploiting a recent "explosion" of salmonid genome assemblies, including a new genome sequence for the huchen (Hucho hucho). We developed a genome alignment approach to capture duplicated regions across multiple species, allowing us to create 121,864 phylogenetic trees describing genome-wide ohnolog divergence across salmonid evolution. Using molecular clock analysis, we show that 61% of the ancestral salmonid genome experienced an initial "wave" of rediploidization in the late Cretaceous (85-106 Ma). This was followed by a period of relative genomic stasis lasting 17-39 My, where much of the genome remained tetraploid. A second rediploidization wave began in the early Eocene and proceeded alongside species diversification, generating predictable patterns of lineage-specific ohnolog divergence, scaling in complexity with the number of speciation events. Using gene set enrichment, gene expression, and codon-based selection analyses, we provide insights into potential functional outcomes of delayed rediploidization. This study enhances our understanding of delayed autopolyploid rediploidization and has broad implications for future studies of WGD events.

Whole genome sequence of bacteremic Clostridium tertium in a World War I soldier, 1914

•

This original article is the first report of isolate and culture of a bacterium from ancient human samples and dental pulp in particular.

•

The dental pulp is a mirror of the individual's infectious state at the time of death.

•

Ancient dental pulp culture yielded to the identification and isolation of the bacterium Clostridium tertium responsible for septicaemia during World War I,.

Background

Dental pulp, encapsulating a blood drop, could be used to diagnose pathogen bacteraemia in archaeological materials using DNA-based techniques. We questioned the viability of such ancient pathogens preserved in ancient dental pulp.

Methods

After meticulous decontamination of 32 teeth collected from 31 World War I soldiers exhumed in Spincourt, France, dental pulps were extracted and cultured under strict anaerobiosis. Colonies were identified by mass spectrometry and whole genome sequencing. Fluorescent in situ hybridisation (FISH) was used for the direct microscopic detection of pathogens of interest in the dental pulp. All the experimental procedures included negative controls, notably sediments in contact with individual SQ517 to ensure that results did not arise from contamination.

Findings

Clostridium tertium was detected by FISH in two dental pulp specimens taken from a 1914 soldier. After a two-day incubation period, both dental pulp samples grew colonies identified by mass spectrometry and genome sequencing as C. tertium; whereas negative controls remained free of C. tertium in all the observations, and no C. tertium was founded in sediments. Skeletal remains of this soldier exhibited two notches in the left tibia evocative of a cold steel wound, and a probably fatal unhealed bullet impact in the hip bone.

Interpretation

Data indicated the presence of C. tertium in the dental pulp at the time of the death of one World War I soldier, in 1914. This observation diagnosed C. tertium bacteraemia, with war wounds as the probable portal of entry for C. tertium. Our C. tertium strains ante-dated by three years, the princeps description of this deadly opportunistic pathogen.

Comparative Population Genomics of Cryptic Speciation and Adaptive Divergence in Bicknell's and Gray-Cheeked Thrushes (Aves: Catharus bicknelli and Catharus minimus)

Cryptic speciation may occur when reproductive isolation is recent or the accumulation of morphological differences between sister lineages is slowed by stabilizing selection preventing phenotypic differentiation. In North America, Bicknell's Thrush (Catharus bicknelli) and its sister species, the Gray-cheeked Thrush (Catharus minimus), are parapatrically breeding migratory songbirds, distinguishable in nature only by subtle differences in song and coloration, and were recognized as distinct species only in the 1990s. Previous molecular studies have estimated that the species diverged approximately 120,000-420,000 YBP and found very low levels of introgression despite their similarity and sympatry in the spring (prebreeding) migration. To further clarify the history, genetic divergence, genomic structure, and adaptive processes in C. bicknelli and C. minimus, we sequenced and assembled high-coverage reference genomes of both species and resequenced genomes from population samples of C. bicknelli, C. minimus, and two individuals of the Swainson's Thrush (Catharus ustulatus). The genome of C. bicknelli exhibits markedly higher abundances of transposable elements compared with other Catharus and chicken. Demographic and admixture analyses confirm moderate genome-wide differentiation (Fst ≈ 0.10) and limited gene flow between C. bicknelli and C. minimus, but suggest a more recent divergence than estimates based on mtDNA. We find evidence of rapid evolution of the Z-chromosome and elevated divergence consistent with natural selection on genomic regions near genes involved with neuronal processes in C. bicknelli. These genomes are a useful resource for future investigations of speciation, migration, and adaptation in Catharus thrushes.

Hybridization drives mitochondrial DNA degeneration and metabolic shift in a species with biparental mitochondrial inheritance

Mitochondrial DNA (mtDNA) is a cytoplasmic genome that is essential for respiratory metabolism. Although uniparental mtDNA inheritance is most common in animals and plants, distinct mtDNA haplotypes can coexist in a state of heteroplasmy, either because of paternal leakage or de novo mutations. mtDNA integrity and the resolution of heteroplasmy have important implications, notably for mitochondrial genetic disorders, speciation, and genome evolution in hybrids. However, the impact of genetic variation on the transition to homoplasmy from initially heteroplasmic backgrounds remains largely unknown. Here, we use Saccharomyces yeasts, fungi with constitutive biparental mtDNA inheritance, to investigate the resolution of mtDNA heteroplasmy in a variety of hybrid genotypes. We previously designed 11 crosses along a gradient of parental evolutionary divergence using undomesticated isolates of Saccharomyces paradoxus and Saccharomyces cerevisiae Each cross was independently replicated 48 to 96 times, and the resulting 864 hybrids were evolved under relaxed selection for mitochondrial function. Genome sequencing of 446 MA lines revealed extensive mtDNA recombination, but the recombination rate was not predicted by parental divergence level. We found a strong positive relationship between parental divergence and the rate of large-scale mtDNA deletions, which led to the loss of respiratory metabolism. We also uncovered associations between mtDNA recombination, mtDNA deletion, and genome instability that were genotype specific. Our results show that hybridization in yeast induces mtDNA degeneration through large-scale deletion and loss of function, with deep consequences for mtDNA evolution, metabolism, and the emergence of reproductive isolation.

Genome Sequence Analysis and Characterization of Shiga Toxin 2 Production by Escherichia coli O157:H7 Strains Associated With a Laboratory Infection

A laboratory-acquired E. coli O157:H7 infection with associated severe sequelae including hemolytic uremic syndrome occurred in an individual working in the laboratory with a mixture of nalidixic acid-resistant (Nal^R) O157:H7 mutant strains in a soil-biochar blend. The patient was hospitalized and treated with an intravenous combination of metronidazole and levofloxacin. The present study investigated the source of this severe laboratory acquired infection and further examined the influence of the antibiotics used during treatment on the expression and production of Shiga toxin. Genomes of two Stx_2a-and eae-positive O157:H7 strains isolated from the patient's stool were sequenced along with two pairs of the wt strains and their derived Nal^R mutants used in the laboratory experiments. High-resolution SNP typing determined the strains' individual genetic relatedness and unambiguously identified the two laboratory-derived Nal^R mutant strains as the source of the researcher's life-threatening disease, rather than a conceivable ingestion of unrelated O157:H7 isolates circulating at the same time. It was further confirmed that in sublethal doses, the antibiotics increased toxin expression and production. Our results support a simultaneous co-infection with clinical strains in the laboratory, which were the causative agents of previous O157:H7 outbreaks, and further that the administration of antibiotics may have impacted the outcome of the infection.

Vertical Inheritance Facilitates Interspecies Diversification in Biosynthetic Gene Clusters and Specialized Metabolites

While specialized metabolites are thought to mediate ecological interactions, the evolutionary processes driving chemical diversification, particularly among closely related lineages, remain poorly understood. Here, we examine the evolutionary dynamics governing the distribution of natural product biosynthetic gene clusters (BGCs) among 118 strains representing all nine currently named species of the marine actinobacterial genus Salinispora. While much attention has been given to the role of horizontal gene transfer (HGT) in structuring BGC distributions, we find that vertical descent facilitates interspecies BGC diversification over evolutionary timescales. Moreover, we identified a distinct phylogenetic signal among Salinispora species at both the BGC and metabolite level, indicating that specialized metabolism represents a conserved phylogenetic trait. Using a combination of genomic analyses and liquid chromatography–high-resolution tandem mass spectrometry (LC-MS/MS) targeting nine experimentally characterized BGCs and their small molecule products, we identified gene gain/loss events, constrained interspecies recombination, and other evolutionary processes associated with vertical inheritance as major contributors to BGC diversification. These evolutionary dynamics had direct consequences for the compounds produced, as exemplified by species-level differences in salinosporamide production. Together, our results support the concept that specialized metabolites, and their cognate BGCs, can represent phylogenetically conserved functional traits with chemical diversification proceeding in species-specific patterns over evolutionary time frames.

Pan-genome Analyses of 3 Strains of Inonotus obliquus and Prediction of Polysaccharide and Terpenoid Genes

Inonotus obliquus is a rare, edible and medicinal fungus that is widely used as a remedy for various diseases. Its main bioactive substances are polysaccharides and terpenoids. In this study, we characterized and investigated the pan-genome of three strains of I. obliquus. The genome sizes of JL01, HE, and NBRC8681 were 32.04, 29.04, and 31.78 Mb, respectively. There were 6 543 core gene families and 6 197 accessory gene families among the three strains, with 14 polysaccharide-related core gene families and seven accessory gene families. For terpenoids, there were 13 core gene families and 17 accessory gene families. Pan-genome sequencing of I. obliquus has improved our understanding of biological characteristics related to the biosynthesis of polysaccharides and terpenoids at the molecular level, which in turn will enable us to increase the production of polysaccharides and terpenoids by this mushroom.

Enhancer architecture-dependent multilayered transcriptional regulation orchestrates RA signaling-induced early lineage differentiation of ESCs

Signaling pathway-driven target gene transcription is critical for fate determination of embryonic stem cells (ESCs), but enhancer-dependent transcriptional regulation in these processes remains poorly understood. Here, we report enhancer architecture-dependent multilayered transcriptional regulation at the Halr1–Hoxa1 locus that orchestrates retinoic acid (RA) signaling-induced early lineage differentiation of ESCs. We show that both homeobox A1 (Hoxa1) and Hoxa adjacent long non-coding RNA 1 (Halr1) are identified as direct downstream targets of RA signaling and regulated by RARA/RXRA via RA response elements (RAREs). Chromosome conformation capture-based screens indicate that RA signaling promotes enhancer interactions essential for Hoxa1 and Halr1 expression and mesendoderm differentiation of ESCs. Furthermore, the results also show that HOXA1 promotes expression of Halr1 through binding to enhancer; conversely, loss of Halr1 enhances interaction between Hoxa1 chromatin and four distal enhancers but weakens interaction with chromatin inside the HoxA cluster, leading to RA signaling-induced Hoxa1 overactivation and enhanced endoderm differentiation. These findings reveal complex transcriptional regulation involving synergistic regulation by enhancers, transcription factors and lncRNA. This work provides new insight into intrinsic molecular mechanisms underlying ESC fate determination during RA signaling-induced early differentiation.

Rapid evolutionary turnover of mobile genetic elements drives bacterial resistance to phages

[Figure: see text].