Gene Synteny Research Articles

Analysis of bZIP transcription factors in Rhododendron simsii and functional study of RsbZIP6 in regulating anthocyanin biosynthesis

The basic leucine zipper (bZIP) transcription factors play a critical role in various plant biological processes, including anthocyanin biosynthesis. This study focuses on Rhododendron simsii, a notable ornamental species with insufficiently explored bZIP transcription factors. We identified 66 bZIP transcription factors in the R. simsii genome and conducted comprehensive bioinformatics analyses to determine their gene localization, phylogenetic relationships, grouping, gene/protein structure, duplication events, synteny, and expression profiles. Our analysis identified RsbZIP6, a homolog of HY5 known to influence anthocyanin biosynthesis in many plants, as a potential regulator of this pathway. We cloned the complete coding sequence of RsbZIP6, which encodes a 170-amino acid protein spanning 510 bp. Subcellular localization analysis verified the nuclear presence of the RsbZIP6 protein. RT-qPCR analysis revealed the highest expression of RsbZIP6 in petals, which correlated with anthocyanin accumulation. Transgenic experiments indicated that overexpressing RsbZIP6 in Arabidopsis enhanced anthocyanin accumulation by upregulating genes involved in anthocyanin biosynthesis (4CL, CHS, CHI, DFR, F3H, F3’H, ANS and UF3GT). Our findings enhance understanding of the bZIP transcription factor family in R. simsii and underscore the vital role of RsbZIP6 in anthocyanin biosynthesis, providing insights for future genetic enhancement strategies.

Identifying Calmodulin and Calmodulin-like Protein Members in Canavalia rosea and Exploring Their Potential Roles in Abiotic Stress Tolerance.

Calmodulins (CaMs) and calmodulin-like proteins (CMLs) belong to families of calcium-sensors that act as calcium ion (Ca2+) signal-decoding proteins and regulate downstream target proteins. As a tropical halophyte, Canavalia rosea shows great resistance to multiple abiotic stresses, including high salinity/alkalinity, extreme drought, heat, and intense sunlight. However, investigations of calcium ion signal transduction involved in the stress responses of C. rosea are limited. The CaM and CML gene families have been identified and characterized in many other plant species. Nevertheless, there is limited available information about these genes in C. rosea. In this study, a bioinformatic analysis, including the gene structures, conserved protein domains, phylogenetic relationships, chromosome distribution, and gene synteny, was comprehensively performed to identify and characterize CrCaMs and CrCMLs. A spatio-temporal expression assay in different organs and environmental conditions was then conducted using the RNA sequencing technique. Additionally, several CrCaM and CrCML members were then cloned and functionally characterized using the yeast heterogeneous expression system, and some of them were found to change the tolerance of yeast to heat, salt, alkalinity, and high osmotic stresses. The results of this study provide a foundation for understanding the possible roles of the CrCaM and CrCML genes, especially for halophyte C. rosea's natural ecological adaptability for its native habitats. This study also provides a theoretical basis for further study of the physiological and biochemical functions of plant CaMs and CMLs that are involved in tolerance to multiple abiotic stresses.

Insertion sequence elements and unique symmetrical genomic regions mediate chromosomal inversions in Streptococcus pyogenes.

Chromosomal inversions are a phenomenon in many bacterial species, often across the axis of replication. Inversions have been shown to alter gene expression, changing persistence of colonisation and infection following environmental stresses. In Streptococcus pyogenes, inversions have been reported. However, frequency and molecular markers of inversions have not been systematically examined. Here, 249 complete S.pyogenes genomes were analysed using a pangenomic core gene synteny framework to identify sequences associated with inversions. 47% of genomes (118/249) contained at least one inversion, from 23 unique inversion locations. Chromosomal locations enabling inversions were usually associated with mobile elements (insertion sequences n=9 and prophages n=7). Two insertion sequences, IS1548 and IS1239, accounted for>80% of insertion sequences and were the only insertion sequences associated with inversions. The most observed inversion location (n=104 genomes, 88% of genomes with an inversion) occurs between two conserved regions encoding rRNAs, tRNAsand sigma factor genes. The regions are symmetrically placed around the origin of replication forming a unique chromosomal structure inS. pyogenes, relative to other streptococci. Cataloging of thechromosomal location and frequency of inversions can direct dissection of phenotypic changes following chromosomal inversions. The framework used here can be transferred to other bacterial species to characterise chromosomal inversions.

Two homologous genes encoding interleukin (IL)-34 in the common carp (Cyprinus carpio L.): Roles in inflammatory modulation and anti-bacterial defense

In mammals, interleukin 34 (IL-34) is a ligand for macrophage colony-stimulating factor receptor (M-CSFR), promoting inflammatory responses and inducing the synthesis and secretion of various cytokines. However, studies on its function in lower vertebrates is limited, and its evolutionary relationship with homologous molecules in mammals remains unclear. In this study, two IL-34-encoding genes were cloned and identified in common carp (Cyprinus carpio L.), designated as CcIL-34A and CcIL-34B, with an amino acid sequence similarity of 77.7 %. Gene synteny analysis revealed that the IL-34 gene loci are relatively conserved, and both are located downstream of SF3B3. The expression patterns of CcIL-34s were analyzed using qRT-PCR, and this showed that they are expressed across all tested tissues, with higher levels in the liver, spleen, and head kidney and lower levels in the gills and intestines. Following infection with Aeromonas hydrophila, the mRNA expression levels of CcIL-34s in the gills, head kidney, intestines, and spleen were significantly upregulated. Immunofluorescence was also employed to assess changes in CcIL-34 protein expression, showing a significant increase in carp spleens 24 h after A. hydrophila infection, suggesting that CcIL-34s contribute to host defense against this bacterium. To investigate the immunological function of IL-34 in vivo, pc-CcIL-34A and pc-CcIL-34B eukaryotic expression plasmids were constructed and injected intramuscularly into fish. Five days after injection, the expression levels of inflammation-related cytokines in the head kidney and spleen were significantly altered. Furthermore, 24 h post-A. hydrophila infection, the bacterial loads in the liver, spleen, and kidneys were significantly reduced. Ten days post-infection, the survival rates in the groups with CcIL-34A and CcIL-34B overexpression were 40 % and 36.7 %, respectively, compared to 16.7 % in the control group. These findings suggest that CcIL-34s are involved in modulating inflammatory responses, enhancing the immune response, and improving survival rates in fish following bacterial infection, thus supporting the potential use of IL-34 molecules in aquaculture.

Chromosome-level subgenome-aware de novo assembly provides insight into Saccharomyces bayanus genome divergence after hybridization.

Interspecies hybridization is prevalent in various eukaryotic lineages and plays important roles in phenotypic diversification, adaptation, and speciation. To better understand the changes that occurred in the different subgenomes of a hybrid species and how they facilitate adaptation, we have completed chromosome-level de novo assemblies of all chromosomes for a recently formed hybrid yeast, Saccharomyces bayanus strain CBS380, using Oxford Nanopore Technologies' MinION long-read sequencing. We characterize the S. bayanus genome and compare it with its parent species, Saccharomyces uvarum and Saccharomyces eubayanus, and other S. bayanus genomes to better understand genome evolution after a relatively recent hybridization event. We observe multiple recombination events between the subgenomes in each chromosome, followed by loss of heterozygosity (LOH) in nine chromosome pairs. In addition to maintaining nearly all gene content and synteny from its parental genomes, S. bayanus has acquired many genes from other yeast species, primarily through the introgression of Saccharomyces cerevisiae, such as those involved in the maltose metabolism. Finally, the patterns of recombination and LOH suggest an allotetraploid origin of S. bayanus The gene acquisition and rapid LOH in the hybrid genome probably facilitated its adaptation to maltose brewing environments and mitigated the maladaptive effect of hybridization. This paper describes the first in-depth study using long-read sequencing technology of an S. bayanus hybrid genome, which may serve as an excellent reference for future studies of this important yeast and other yeast strains.

Characterization and comparative analysis of sericin protein 150 in Bombyx mori

Lepidopteran silk is a complex mixture of proteins, consisting mainly of fibroins and sericins. Sericins are a small family of highly divergent proteins that serve as adhesives and coatings for silk fibers. So far, five genes encoding sericin proteins have been identified in Bombyx mori. Having previously identified sericin protein 150 (SP150) as a major sericin-like protein in the cocoons of the pyralid moths Galleria mellonella and Ephestia kuehniella, we describe the identification of its homolog in B. mori. Our refined gene model shows that it consists of four exons and a long open reading frame with a conserved motif, CXCXCX, at the C-terminus, reminiscent of the structure observed in a class of mucin proteins. Notably, despite a similar expression pattern, both mRNA and protein levels of B. mori SP150 were significantly lower than those of its pyralid counterpart. We also discuss the synteny of homologous genes on corresponding chromosomes in different moth species and the possible phylogenetic relationships between SP150 and certain mucin-like proteins. Our results improve our understanding of silk structure and the evolutionary relationships between adhesion proteins in the silk of different lepidopteran species.

Characterisation of the complete chloroplast genome of Solanum tuberosum cv. White Lady

AbstractPotato (Solanum tuberosum) is considered worldwide as one of the most important non-cereal food crops. As a result of its adaptability and worldwide production area, potato displays a vast phenotypical variability as well as genomic diversity. Chloroplast genomes have long been a core issue in plant molecular evolution and phylogenetic studies, and have an important role in revealing photosynthetic mechanisms, metabolic regulations and the adaptive evolution of plants. We sequenced the complete chloroplast genome of the Hungarian cultivar White Lady, which is 155 549 base pairs (bp) in length and is characterised by the typical quadripartite structure composed of a large- and small single-copy region (85 991 bp and 18 374 bp, respectively) interspersed by two identical inverted repeats (25 592 bp). The genome consists of 127 genes of which 82 are protein-coding, eight are ribosomal RNAs and 37 are transfer RNAs. The overall gene content and distribution of the genes on the White Lady chloroplast was the same as found in other potato chloroplasts. The alignment of S. tuberosum chloroplast genome sequences resulted in a highly resolved tree, with 10 out of the 13 nodes recovered having bootstrap values over 90%. By comparing the White Lady chloroplast genome with available S. tuberosum sequences we found that gene content and synteny are highly conserved. The new chloroplast sequence can support further studies of genetic diversity, resource conservation, evolution and applied agricultural research. The new sequence can support further potato genetic diversity and evolutionary studies, resource conservation, and also applied agricultural research.

The first Chromosomal-level genome assembly of Sageretia thea using Nanopore long reads and Pore-C technology

Sageretia thea, a notable species within the mock buckthorn genus, is recognized for its intriguing biogeographical distribution and diverse medicinal properties. Despite this significance, genomic studies on S. thea are still in the nascent stages. We present the first chromosome-level genome assembly of S. thea that was generated using a combination of Oxford Nanopore long-read and Illumina short-read sequencing technologies complemented by Pore-C chromatin conformation capture. The genome assembly had a size of 197.8 Mb with 12 chromosomal scaffolds and a scaffold N50 length of 15.9 Mb. A total of 25,434 protein-coding genes were identified and functionally annotated, and the gene model indicated 96.5% complete eukaryotic BUSCOs. Additionally, orthologous gene profiling and synteny analysis were performed to elucidate the evolutionary relationships within the Rhamnaceae family and Rosales. This high-quality chromosomal genome is the first genomic view of S. thea, which will serve as the basis for future studies on its biological and medicinal properties, and evolutionary history.

Genome-Wide Identification and Characterization of the GASA Gene Family in Medicago truncatula, and Expression Patterns under Abiotic Stress and Hormone Treatments.

Medicago truncatula is a key model plant for studying legume plants, particularly alfalfa (Medicago sativa), due to its well-defined genetic background. Plant-specific GASA (Gibberellic Acid Stimulated Arabidopsis) genes play various roles in plant growth and development, abiotic stress, and hormone responses. However, limited information is available on GASA research in Medicago. In this study, 26 MtGASAs were identified and analyzed for its structure, evolution, and expressions. Sequence alignments and phylogeny revealed that 26 MtGASAs containing conserved GASA domains were classified into three clades. The chromosomal locations and gene synteny revealed segmental and tandem repetition evolution. Analysis of cis-regulatory elements indicates that family members likely influence various hormone signaling pathways and stress-related mechanisms. Moreover, the RNA-seq and qRT-PCR analyses revealed that 26 MtGASAs were extensively involved in abiotic stresses and hormone responses. Notably, seven MtGASA genes (MtGASA1, 10, 12, 17, 23, 25 and 26) were all dramatically activated by NaCl and Mannitol treatments, and four MtGASAs (MtGASA7, 10, 23 and 24) were significant activated by GA3, PBZ, ABA, and MeJA treatments. Collectively, this study is the first to identify and describe GASA genes in Medicago on a genome-wide scale. The results establish a basis for functional characterization, showing that these proteins are essential in responding to various abiotic stresses and hormonal signals.

SynGAP: a synteny-based toolkit for gene structure annotation polishing

Genome sequencing has become a routine task for biologists, but the challenge of gene structure annotation persists, impeding accurate genomic and genetic research. Here, we present a bioinformatics toolkit, SynGAP (Synteny-based Gene structure Annotation Polisher), which uses gene synteny information to accomplish precise and automated polishing of gene structure annotation of genomes. SynGAP offers exceptional capabilities in the improvement of gene structure annotation quality and the profiling of integrative gene synteny between species. Furthermore, an expression variation index is designed for comparative transcriptomics analysis to explore candidate genes responsible for the development of distinct traits observed in phylogenetically related species.

The hare syphilis agent is related to, but distinct from, the treponeme causing rabbit syphilis.

The treponemes infecting lagomorphs include Treponema paraluisleporidarum ecovar Cuniculus (TPeC) and ecovar Lepus (TPeL), infecting rabbits and hares, respectively. In this study, we described the first complete genome sequence of TPeL, isolate V3603-13, from an infected mountain hare (Lepus timidus) in Sweden. In addition, we determined 99.0% of the genome sequence of isolate V246-08 (also from an infected mountain hare, Sweden) and 31.7% of the genome sequence of isolate Z27 A77/78 (from a European hare, Lepus europeaus, The Netherlands). The TPeL V3603-13 genome had considerable gene synteny with the TPeC Cuniculi A genome and with the human pathogen T. pallidum, which causes syphilis (ssp. pallidum, TPA), yaws (ssp. pertenue, TPE) and endemic syphilis (ssp. endemicum, TEN). Compared to the TPeC Cuniculi A genome, TPeL V3603-13 contained four insertions and 11 deletions longer than three nucleotides (ranging between 6 and2,932 nts). In addition, there were 25 additional indels, from one to three nucleotides long, altogether spanning 36 nts. The number of single nucleotide variants (SNVs) between TPeC Cuniculi A and TPeL V3603-13 were represented by 309 nucleotide differences. Major proteome coding differences between TPeL and TPeC were found in the tpr gene family, and (predicted) genes coding for outer membrane proteins, suggesting that these components are essential for host adaptation in lagomorph syphilis. The phylogeny revealed that the TPeL sample from the European brown hare was more distantly related to TPeC Cuniculi A than V3603-13 and V246-08.

Genomic insights into cytokinin oxidase/dehydrogenase (CKX) gene family, identification, phylogeny and synteny analysis for its possible role in regulating seed number in Pigeonpea (Cajanus cajan (L.) Millsp.)

Cytokinin oxidase/dehydrogenase (CKX) regulates cytokinin levels in plants which are vital for plant growth and development. However, there is a paucity of evidence regarding their role in controlling embryo/seed development in pigeonpea. This comprehensive study provides information on the identification and characterization of CKX genes in pigeonpea. A genome-wide analysis identified 18 CKX genes, each with distinct structure, expression patterns, and possible diverse functions. Domain analysis revealed the presence of the sequences including FAD and CK-Binding domain, and subcellular localization analysis showed that almost 50 % of them reside within the nucleus. They were observed to be located unevenly on chromosome numbers 2, 4, 6, 7, and 11 with a majority of them present on the scaffolds. The 8 homologous pairs and various orthologous gene pairs provided further insights into their evolution pattern. Further, SNP/Indels variation in CKX genes and haplotype groups among contrasting genotypes for SNPP (seed number per pod) were analyzed. Spatiotemporal expression analysis revealed the significant expression pattern of CcCKX15, CcCKX17, and CcCKX2 in genotypes carrying low SNPP reiterating their possible role as negative regulators. These genes can be potential targets to undertake seed and biomass improvement in pigeonpea.

A novel genus of Pectobacterium bacteriophages display broad host range by targeting several species of Danish soft rot isolates

The bacterial diseases black leg and soft rot in potatoes cause heavy losses of potatoes worldwide. Bacteria within the genus Pectobacteriaceae are the causative agents of black leg and soft rot. The use of antibiotics in agriculture is heavily regulated and no other effective treatment currently exists, but bacteriophages (phages) have shown promise as potential biocontrol agents. In this study we isolated soft rot bacteria from potato tubers and plant tissue displaying soft rot or black leg symptoms collected in Danish fields. We then used the isolated bacterial strains as hosts for phage isolation. Using organic waste, we isolated phages targeting different species within Pectobacterium. Here we focus on seven of these phages representing a new genus primarily targeting P. brasiliense; phage Ymer, Amona, Sabo, Abuela, Koroua, Taid and Pappous. TEM image of phage Ymer showed siphovirus morphotype, and the proposed Ymer genus belongs to the class Caudoviricetes, with double-stranded DNA genomes varying from 39 kb to 43 kb. In silico host range prediction using a CRISPR-Cas spacer database suggested both P. brasiliense, P. polaris and P. versatile as natural hosts for phages within the proposed Ymer genus. A following host range experiment, using 47 bacterial isolates from Danish tubers and plants symptomatic with soft rot or black leg disease verified the in silico host range prediction, as the genus as a group were able to infect all three Pectobacterium species. Phages did, however, primarily target P. brasiliense isolates and displayed differences in host range even within the species level. Two of the phages were able to infect two or more Pectobacterium species. Despite no nucleotide similarity with any phages in the NCBI database, the proposed Ymer genus did share some similarity at the protein level, as well as gene synteny, with currently known phages. None of the phages encoded integrases or other genes typically associated with lysogeny. Similarly, no virulence factors nor antimicrobial resistance genes were found, and combined with their ability to infect several soft rot-causing Pectobacterium species from Danish fields, demonstrates their potential as biocontrol agents against soft rot and black leg diseases in potatoes.

Phylogenetic incongruence in an Asiatic species complex of the genus Caryodaphnopsis (Lauraceae)

BackgroundCaryodaphnopsis, a group of tropical trees (ca. 20 spp.) in the family Lauraceae, has an amphi-Pacific disjunct distribution: ten species are distributed in Southeast Asia, while eight species are restricted to tropical rainforests in South America. Previously, phylogenetic analyses using two nuclear markers resolved the relationships among the five species from Latin America. However, the phylogenetic relationships between the species in Asia remain poorly known.ResultsHere, we first determined the complete mitochondrial genome (mitogenome), plastome, and the nuclear ribosomal cistron (nrDNA) sequences of C. henryi with lengths of 1,168,029 bp, 154,938 bp, and 6495 bp, respectively. We found 2233 repeats and 368 potential SSRs in the mitogenome of C. henryi and 50 homologous DNA fragments between its mitogenome and plastome. Gene synteny analysis revealed a mass of rearrangements in the mitogenomes of Magnolia biondii, Hernandia nymphaeifolia, and C. henryi and only six conserved clustered genes among them. In order to reconstruct relationships for the ten Caryodaphnopsis species in Asia, we created three datasets: one for the mitogenome (coding genes and ten intergenic regions), another for the plastome (whole genome), and the other for the nuclear ribosomal cistron. All of the 22 Caryodaphnopsis individuals were divided into four, five, and six different clades in the phylogenies based on mitogenome, plastome, and nrDNA datasets, respectively.ConclusionsThe study showed phylogenetic conflicts within and between nuclear and organellar genome data of Caryodaphnopsis species. The sympatric Caryodaphnopsis species in Hekou and Malipo SW China may be related to the incomplete lineage sorting, chloroplast capture, and/or hybridization, which mixed the species as a complex in their evolutionary history.

The genomic evolution of visual opsin genes in amphibians

Among tetrapod (terrestrial) vertebrates, amphibians remain more closely tied to an amphibious lifestyle than amniotes, and their visual opsin genes may be adapted to this lifestyle. Previous studies have discussed physiological, morphological, and molecular changes in the evolution of amphibian vision. We predicted the locations of the visual opsin genes, their neighboring genes, and the tuning sites of the visual opsins, in 39 amphibian genomes. We found that all of the examined genomes lacked the Rh2 gene. The caecilian genomes have further lost the SWS1 and SWS2 genes; only the Rh1 and LWS genes were retained. The loss of the SWS1 and SWS2 genes in caecilians may be correlated with their cryptic lifestyles. The opsin gene syntenies were predicted to be highly similar to those of other bony vertebrates. Moreover, dual syntenies were identified in allotetraploid Xenopus laevis and X. borealis. Tuning site analysis showed that only some Caudata species might have UV vision. In addition, the S164A that occurred several times in LWS evolution might either functionally compensate for the Rh2 gene loss or fine-tuning visual adaptation. Our study provides the first genomic evidence for a caecilian LWS gene and a genomic viewpoint of visual opsin genes by reviewing the gains and losses of visual opsin genes, the rearrangement of syntenies, and the alteration of spectral tuning in the course of amphibians’ evolution.

The Description and Analysis of the Complete Genome of Dermacoccus barathri FBCC-B549 Strain.

Dermacoccus barathri is the first reported pathogen within the Dermacoccus genus to cause a catheter-related bloodstream infection, which occurred in 2015. In this study, the complete genome assembly of Dermacoccus barathri was constructed, and the complete genome of Dermacoccus barathri FBCC-B549 consists of a single chromosome (3,137,745 bp) without plasmids. The constructed genome of D. barathri was compared with those of two closely related species within the Dermacoccus genus. D. barathri exhibited a pattern similar to Dermacoccus abyssi in terms of gene clusters and synteny analysis. Contrary to previous studies, biosynthetic gene cluster (BGC) analysis for predicting secondary metabolites revealed the presence of the LAP biosynthesis pathway in the complete genome of D. barathri, predicting the potential synthesis of the secondary metabolite plantazolicin. Furthermore, an analysis to investigate the potential pathogenicity of D. barathri did not reveal any antibiotic resistance genes; however, nine virulence factors were identified in the Virulence Factor Database (VFDB). According to these matching results in the VFDB, despite identifying a few factors involved in biofilm formation, further research is required to determine the actual impact of D. barathri on pathogenicity. The complete genome of D. barathri is expected to serve as a valuable resource for future studies on D. barathri, which currently lack sufficient genomic sequence information.

The chromosome level whole genome sequence and the seconary matabolism gene cluster prediction of Fusarium meridionale, the pathogen causing maize ear rot.

Fusarium meridionale is one of the pathogens causing maize ear rot, it produce bioactive secondary metabolites may threaten humans food safty, however, the production mechanism of the secondary metabolites and their interaction with maize ear remains poorly understood. To facilitate related studies, we sequenced and assembled the genome of F. meridionale strain JX18-4. The size of F. meridionale JX18-4 genome is 37.11 Mbp, include four nuclear chromosome contigs that consists of 11920 coding genes and one mitochondrial contig. 95.64% gene synteny collinearity was found between the assembly and the reference genomes F. graminearum strain PH-1. Compared to the sequences of seconary matabolism gene clusters sequences reported previously, the stain JX18-4 was predicted potential producing 8 clusters, including nivalenol, zearalenone, aurofusarin, fusarielin, fusaristatin A, fusarin, fusarubin and butenolide. This study aims to reveal the molecular mechanism of secondary metabolites producing, and the genomic information of JX18-4 will provide resources for the study of biological control mechanisms and plant-microbe interactions.

Comparative genomics reveals high genetic similarity among strains of Salmonella enterica serovar Infantis isolated from multiple sources in Brazil.

Salmonella enterica serovar Infantis (Salmonella Infantis) is a zoonotic, ubiquitous and foodborne pathogen of worldwide distribution. Despite Brazil's relevance as a major meat exporter, few studies were conducted to characterize strains of this serovar by genomic analyses in this country. Therefore, this study aimed to assess the diversity of 80 Salmonella Infantis strains isolated from veterinary, food and human sources in Brazil between 2013 and 2018 by comparative genomic analyses. Additional genomes of non-Brazilian countries (n=18) were included for comparison purposes in some analyses. Analyses of whole-genome multi-locus sequence typing (wgMLST), using PGAdb-builder, and of fragmented genomes, using Gegenees, were conducted to compare the 80 Brazilian strains to the 18 non-Brazilian genomes. Pangenome analyses and calculations were performed for all Salmonella Infantis genomes analyzed. The presence of prophages was determined using PHASTER for the 80 Brazilian strains. The genome plasticity using BLAST Ring Image Generator (BRIG) and gene synteny using Mauve were evaluated for 20 selected Salmonella Infantis genomes from Brazil and ten from non-Brazilian countries. Unique orthologous protein clusters were searched in ten selected Salmonella Infantis genomes from Brazil and ten from non-Brazilian countries. wgMLST and Gegenees showed a high genomic similarity among some Brazilian Salmonella Infantis genomes, and also the correlation of some clusters with non-Brazilian genomes. Gegenees also showed an overall similarity >91% among all Salmonella Infantis genomes. Pangenome calculations revealed an open pangenome for all Salmonella Infantis subsets analyzed and a high gene content in the core genomes. Fifteen types of prophages were detected among 97.5% of the Brazilian strains. BRIG and Mauve demonstrated a high structural similarity among the Brazilian and non-Brazilian isolates. Unique orthologous protein clusters related to biological processes, molecular functions, and cellular components were detected among Brazilian and non-Brazilian genomes. The results presented using different genomic approaches emphasized the significant genomic similarity among Brazilian Salmonella Infantis genomes analyzed, suggesting wide distribution of closely related genotypes among diverse sources in Brazil. The data generated contributed to novel information regarding the genomic diversity of Brazilian and non-Brazilian Salmonella Infantis in comparison. The different genetically related subtypes of Salmonella Infantis from Brazil can either occur exclusively within the country, or also in other countries, suggesting that some exportation of the Brazilian genotypes may have already occurred.

First Pangenome of Corynebacterium rouxii, a Potentially Toxigenic Species of Corynebacterium diphtheriae Complex

Corynebacterium rouxii is one of the recently described species of the Corynebacterium diphtheriae complex. As this species can potentially infect different hosts and harbor the tox gene, producing diphtheria toxin, we present its first pangenomic analysis in this work. A total of fifteen genomes deposited in online databases were included. After confirming the taxonomic position of the isolates by genomic taxonomy, the genomes were submitted to genomic plasticity, gene synteny, and pangenome prediction analyses. In addition, virulence and antimicrobial resistance genes were investigated. Finally, epidemiological data were obtained through molecular typing, clustering, and phylogenetic analysis. Our data demonstrated genetic diversity within the species with low synteny. However, the gene content is extensively conserved, and the pangenome is composed of 2606 gene families, of which 1916 are in the core genome and 80 are related to unique genes. Prophages, insertion sequences, and genomic islands were found. A type I-E CRISPR-Cas system was also detected. Besides the tox gene, determinants involved in adhesion and iron acquisition and two putative antimicrobial resistance genes were predicted. These findings provide valuable insight about this species’ pathogenicity, evolution, and diversity. In the future, our data can contribute to different areas, including vaccinology and epidemiology.

LiftoffTools: a toolkit for comparing gene annotations mapped between genome assemblies.

In 2020 we published Liftoff, which was the first standalone tool specifically designed for transferring gene annotations between genome assemblies of the same or closely related species. While the gene content is expected to be very similar in closely related genomes, the differences may be biologically consequential, and a computational method to extract all gene-related differences should prove useful in the analysis of such genomes. Here we present LiftoffTools, a toolkit to automate the detection and analysis of gene sequence variants, synteny, and gene copy number changes. We provide a description of the toolkit and an example of its use comparing genes mapped between two human genome assemblies.