Robust identification of orthologues and paralogues for microbial pan-genomics using GET_HOMOLOGUES: a case study of pIncA/C plasmids

Methods for Pangenomic Core Detection

Computational pangenomics deals with the joint analysis of all genomic sequences of a species. It has already been successfully applied to various tasks in many research areas. Further advances in DNA sequencing technologies constantly let more and more genomic sequences become available for many species, leading to an increasing attractiveness of pangenomic studies. At the same time, larger datasets also pose new challenges for data structures and algorithms that are needed to handle the data. Efficient methods oftentimes make use of the concept of k-mers.Core detection is a common way of analyzing a pangenome. The pangenome's core is defined as the subset of genomic information shared among all individual members. Classically, it is not only determined on the abstract level of genes but can also be described on the sequence level.In this chapter, we provide an overview of k-mer-based methods in the context of pangenomics studies. We first revisit existing software solutions for k-mer counting and k-mer set representation. Afterward, we describe the usage of two k-mer-based approaches, Pangrowth and Corer, for pangenomic core detection.

Tulathromycin metaphylaxis increases nasopharyngeal isolation of multidrug resistant Mannheimia haemolytica in stocker heifers

Bovine respiratory disease (BRD) is a leading cause of disease in feedlot and stocker calves with Mannheimia haemolytica (MH) as one of the most common etiologies. One of the most effective means of controlling BRD is through metaphylaxis, which involves administering antimicrobials to all animals at high risk of developing BRD. However, increasing prevalence of multidrug resistant (MDR) MH may reduce efficacy of metaphylaxis due to decreased susceptibility to drugs used for metaphylaxis. Primarily, this study aimed to determine the effect of tulathromycin metaphylaxis and subsequent BRD treatment on antimicrobial resistance (AMR) in MH isolated from stocker calves. Secondary objectives included evaluating the effect of metaphylaxis and treatment for BRD on animal health and comparing the genetic relationship of MH isolated. Crossbred beef heifers (n = 331, mean weight = 232, SD = 17.8 kg) at high risk for BRD were randomly assigned to receive tulathromycin metaphylaxis (META, n = 167) or not (NO META, n = 164). Nasopharyngeal swabs were collected for MH isolation, antimicrobial susceptibility testing and whole genome sequencing at arrival and 3 (WK3) and 10 (WK10) weeks later. Mixed-effects logistic regression was used to identify risk factors for isolation of MH and MDR MH (resistant to ≥3 antimicrobial drug classes) at 3 and 10 weeks, BRD morbidity, and crude mortality. Animals in the META group had higher odds of isolation of MDR MH at 3 weeks [OR (95% CI) = 13.08 (5-30.9), p < 0.0001] and 10 weeks [OR (95% CI) = 5.92 (1.34-26.14), p = 0.019] after arrival. There was no difference in risk of isolation of any MH (resistant or susceptible) between META and NO META groups at all timepoints. Animals in the NO META group had 3 times higher odds of being treated for BRD [WK3: OR (95% CI) = 3.07 (1.70-5.52), p = 0.0002; WK10: OR (95% CI) = 2.76 (1.59-4.80), p = 0.0002]. Antimicrobial resistance genes found within isolates were associated with integrative conjugative element (ICE) genes. Tulathromycin metaphylaxis increased risk of isolation of MDR MH and in this population, the increase in MDR MH appeared to be associated with ICE containing antimicrobial resistance genes for multiple antimicrobial classes. This may have important implications for future efficacy of antimicrobials for control and treatment of BRD.

The Raphidiopsis (= Cylindrospermopsis) raciborskii pangenome updated: Two new metagenome-assembled genomes from the South American clade

Two Raphidiopsis (=Cylindrospermopsis) raciborskii metagenome-assembled genomes (MAGs) were recovered from two freshwater metagenomic datasets sampled in 2011 and 2012 in Pampulha Lake, a hypereutrophic, artificial, shallow reservoir, located in the city of Belo Horizonte (MG), Brazil. Since the late 1970s, the lake has undergone increasing eutrophication pressure, due to wastewater input, leading to the occurrence of frequent cyanobacterial blooms. The major difference observed between PAMP2011 and PAMP2012 MAGs was the lack of the saxitoxin gene cluster in PAMP2012, which also presented a smaller genome, while PAMP2011 presented the complete sxt cluster and all essential proteins and clusters. The pangenome analysis was performed with all Raphidiopsis/Cylindrospermopsis genomes available at NCBI to date, with the addition of PAMP2011 and PAMP2012 MAGs (All33 subset), but also without the South American strains (noSA subset), and only among the South American strains (SA10 and SA8 subsets). We observed a substantial increase in the core genome size for the 'noSA' subset, in comparison to 'All33' subset, and since the core genome reflects the closeness among the pangenome members, the results strongly suggest that the conservation level of the essential gene repertoire seems to be affected by the geographic origin of the strains being analyzed, supporting the existence of a distinct SA clade. The Raphidiopsis pangenome comprised a total of 7943 orthologous protein clusters, and the two new MAGs increased the pangenome size by 11%. The pangenome based phylogenetic relationships among the 33 analyzed genomes showed that the SA genomes clustered together with 99% bootstrap support, reinforcing the metabolic particularity of the Raphidiopsis South American clade, related to its saxitoxin producing unique ability, while also indicating a different evolutionary history due to its geographic isolation.

Functional annotation and comparative genomics analysis of Balamuthia mandrillaris reveals potential virulence-related genes

Balamuthia mandrillaris is a pathogenic protozoan that causes a rare but almost always fatal infection of the central nervous system and, in some cases, cutaneous lesions. Currently, the genomic data for this free-living amoeba include the description of several complete mitochondrial genomes. In contrast, two complete genomes with draft quality are available in GenBank, but none of these have a functional annotation. In the present study, the complete genome of B. mandrillaris isolated from a freshwater artificial lagoon was sequenced and assembled, obtaining an assembled genome with better assembly quality parameter values than the currently available genomes. Afterward, the genome mentioned earlier, along with strains V039 and 2046, were subjected to functional annotation. Finally, comparative genomics analysis was performed, and it was found that homologous genes in the core genome potentially involved in the virulence of Acanthamoeba spp. and Trypanosoma cruzi. Moreover, eleven of fifteen genes were identified in the three strains described as potential target genes to develop new treatment approaches for B. mandrillaris infections. These results describe proteins in this protozoan's complete genome and help prioritize which target genes could be used to develop new treatments.

Whither the genus Caldicellulosiruptor and the order Thermoanaerobacterales: phylogeny, taxonomy, ecology, and phenotype

The order Thermoanaerobacterales currently consists of fermentative anaerobic bacteria, including the genus Caldicellulosiruptor. Caldicellulosiruptor are represented by thirteen species; all, but one, have closed genome sequences. Interest in these extreme thermophiles has been motivated not only by their high optimal growth temperatures (≥70°C), but also by their ability to hydrolyze polysaccharides including, for some species, both xylan and microcrystalline cellulose. Caldicellulosiruptor species have been isolated from geographically diverse thermal terrestrial environments located in New Zealand, China, Russia, Iceland and North America. Evidence of their presence in other terrestrial locations is apparent from metagenomic signatures, including volcanic ash in permafrost. Here, phylogeny and taxonomy of the genus Caldicellulosiruptor was re-examined in light of new genome sequences. Based on genome analysis of 15 strains, a new order, Caldicellulosiruptorales, is proposed containing the family Caldicellulosiruptoraceae, consisting of two genera, Caldicellulosiruptor and Anaerocellum. Furthermore, the order Thermoanaerobacterales also was re-assessed, using 91 genome-sequenced strains, and should now include the family Thermoanaerobacteraceae containing the genera Thermoanaerobacter, Thermoanaerobacterium, Caldanaerobacter, the family Caldanaerobiaceae containing the genus Caldanaerobius, and the family Calorimonaceae containing the genus Calorimonas. A main outcome of ANI/AAI analysis indicates the need to reclassify several previously designated species in the Thermoanaerobacterales and Caldicellulosiruptorales by condensing them into strains of single species. Comparative genomics of carbohydrate-active enzyme inventories suggested differentiating phenotypic features, even among strains of the same species, reflecting available nutrients and ecological roles in their native biotopes.

Metabolic capabilities are highly conserved among human nasal-associated Corynebacterium species in pangenomic analyses

Corynebacterium species are globally ubiquitous in human nasal microbiota across the lifespan. Moreover, nasal microbiota profiles typified by higher relative abundances of Corynebacterium are often positively associated with health. Among the most common human nasal Corynebacterium species are C. propinquum, C. pseudodiphtheriticum, C. accolens, and C. tuberculostearicum. Based on the prevalence of these species, at least two likely coexist in the nasal microbiota of 82% of adults. To gain insight into the functions of these four species, we identified genomic, phylogenomic, and pangenomic properties and estimated the functional protein repertoire and metabolic capabilities of 87 distinct human nasal Corynebacterium strain genomes: 31 from Botswana and 56 from the U.S. C. pseudodiphtheriticum had geographically distinct clades consistent with localized strain circulation, whereas some strains from the other species had wide geographic distribution across Africa and North America. All four species had similar genomic and pangenomic structures. Gene clusters assigned to all COG metabolic categories were overrepresented in the persistent (core) compared to the accessory genome of each species indicating limited strain-level variability in metabolic capacity. Moreover, core metabolic capabilities were highly conserved among the four species indicating limited species-level metabolic variation. Strikingly, strains in the U.S. clade of C. pseudodiphtheriticum lacked genes for assimilatory sulfate reduction present in the Botswanan clade and in the other studied species, indicating a recent, geographically related loss of assimilatory sulfate reduction. Overall, the minimal species and strain variability in metabolic capacity implies coexisting strains might have limited ability to occupy distinct metabolic niches.

Genomic insight of sulfate reducing bacterial genus Desulfofaba reveals their metabolic versatility in biogeochemical cycling

BACKGROUND

Sulfate-reducing bacteria (SRB) drive the ocean sulfur and carbon cycling. They constitute a diverse phylogenetic and physiological group and are widely distributed in anoxic marine environments. From a physiological viewpoint, SRB's can be categorized as complete or incomplete oxidizers, meaning that they either oxidize their carbon substrate completely to CO₂ or to a stoichiometric mix of CO₂ and acetate. Members of Desulfofabaceae family are incomplete oxidizers, and within that family, Desulfofaba is the only genus with three isolates that are classified into three species. Previous physiological experiments revealed their capability of respiring oxygen.

RESULTS

Here, we sequenced the genomes of three isolates in Desulfofaba genus and reported on a genomic comparison of the three species to reveal their metabolic potentials. Based on their genomic contents, they all could oxidize propionate to acetate and CO₂. We confirmed their phylogenetic position as incomplete oxidizers based on dissimilatory sulfate reductase (DsrAB) phylogeny. We found the complete pathway for dissimilatory sulfate reduction, but also different key genes for nitrogen cycling, including nitrogen fixation, assimilatory nitrate/nitrite reduction, and hydroxylamine reduction to nitrous oxide. Their genomes also contain genes that allow them to cope with oxygen and oxidative stress. They have genes that encode for diverse central metabolisms for utilizing different substrates with the potential for more strains to be isolated in the future, yet their distribution is limited.

CONCLUSIONS

Results based on marker gene search and curated metagenome assembled genomes search suggest a limited environmental distribution of this genus. Our results reveal a large metabolic versatility within the Desulfofaba genus which establishes their importance in biogeochemical cycling of carbon in their respective habitats, as well as in the support of the entire microbial community through releasing easily degraded organic matters.

Complete Genome Sequences of Two Thermophilic Indigenous Bacteria Isolated from Wheat Straw, Thermoclostridium stercorarium subsp. Strain RKWS1 and Thermoanaerobacter sp. Strain RKWS2

Reported here are complete genome sequences for two anaerobic, thermophilic bacteria isolated from wheat straw, i.e., the (hemi)cellulolytic Thermoclostridium stercorarium subspecies strain RKWS1 (3,029,933 bp) and the hemicellulolytic Thermoanaerobacter species strain RKWS2 (2,827,640 bp). Discovery of indigenous thermophiles in plant biomass suggests that high-temperature microorganisms are more ubiquitous than previously thought.

CoreGenes5.0: An Updated User-Friendly Webserver for the Determination of Core Genes from Sets of Viral and Bacterial Genomes

The determination of core genes in viral and bacterial genomes is crucial for a better understanding of their relatedness and for their classification. CoreGenes5.0 is an updated user-friendly web-based software tool for the identification of core genes in and data mining of viral and bacterial genomes. This tool has been useful in the resolution of several issues arising in the taxonomic analysis of bacteriophages and has incorporated many suggestions from researchers in that community. The webserver displays result in a format that is easy to understand and allows for automated batch processing, without the need for any user-installed bioinformatics software. CoreGenes5.0 uses group protein clustering of genomes with one of three algorithm options to output a table of core genes from the input genomes. Previously annotated "unknown genes" may be identified with homologues in the output. The updated version of CoreGenes is able to handle more genomes, is faster, and is more robust, providing easier analysis of custom or proprietary datasets. CoreGenes5.0 is accessible at coregenes.org, migrating from a previous site.

Pan-genome evolution and its association with divergence of metabolic functions in Bifidobacterium genus

Previous studies were mainly focused on genomic evolution and diversity of type species of Bifidobacterium genus due to their health-promoting effect on host. However, those studies were mainly based on species-level taxonomic resolution, adaptation, and characterization of carbohydrate metabolic features of the bifidobacterial species. Here, a comprehensive analysis of the type strain genome unveils the association of pan-genome evolution with the divergence of metabolic function of the Bifidobacterium genus. This study has also demonstrated that horizontal gene transfer, as well as genome expansion and reduction events, leads to the divergence of metabolic functions in Bifidobacterium genus. Furthermore, the genome-based search of probiotic traits among all the available bifidobacterial type strains gives hints on type species, that could confer health benefits to nutrient-deficient individuals. Altogether, the present study provides insight into the developments of genomic evolution, functional divergence, and potential probiotic type species of the Bifidobacterium genus.

The divergent genome of Scorpion Group 1 (SG1) intracellular bacterium from the venom glands of Vaejovis smithi (Scorpiones: Vaejovidae)

Scorpions were among the first animals on land around 430 million years ago. Like many arachnids, scorpions have evolved complex venoms used to paralyze their prey and for self-defense. Here we sequenced and analyzed the metagenomic DNA from venom glands from Vaejovis smithi scorpions. A metagenome-assembled genome (MAG) of 624,025 bp was obtained corresponding to the previously reported Scorpion Group 1 (SG1). The SG1 genome from venom glands had a low GC content (25.8%) characteristic of reduced genomes, many hypothetical genes and genes from the reported minimal set of bacterial genes. Phylogenomic reconstructions placed the uncultured SG1 distant from other reported bacteria constituting a taxonomic novelty. By PCR we detected SG1 in all tested venom glands from 30 independent individuals. Microscopically, we observed SG1 inside epithelial cells from the venom glands using FISH and its presence in scorpion embryos suggested that SG1 is transferred from mother to offspring.

A Novel Potent Crystalline Chitin Decomposer: Chitin Deacetylase from Acinetobacter schindleri MCDA01

Chitosan is a functional ingredient that is widely used in food chemistry as an emulsifier, flocculant, antioxidant, or preservative. Chitin deacetylases (CDAs) can catalyze the hydrolysis of acetyl groups, making them useful in the clean production of chitosan. However, the high inactivity of crystalline chitin catalyzed by CDAs has been regarded as the technical bottleneck of crystalline chitin deacetylation. Here, we mined the AsCDA gene from the genome of Acinetobacter schindleri MCDA01 and identified a member of the uraD_N-term-dom superfamily, which was a novel chitin deacetylase with the highest deacetylation activity. The AsCDA gene was expressed in Escherichia coli BL21 by IPTG induction, whose activity to colloidal chitin, α-chitin, and β-chitin reached 478.96 U/mg, 397.07 U/mg, and 133.27 U/mg, respectively. In 12 h, the enzymatic hydrolysis of AsCDA removed 63.05% of the acetyl groups from α-chitin to prepare industrial chitosan with a degree of deacetylation higher than 85%. AsCDA, as a potent chitin decomposer in the production of chitosan, plays a positive role in the upgrading of the chitosan industry and the value-added utilization of chitin biological resources.

Genome-based reclassification of Azospirillum brasilense Az39 as the type strain of Azospirillum argentinense sp. nov

Strain Az39T of Azospirillum is a diazotrophic plant growth-promoting bacterium isolated in 1982 from the roots of wheat plants growing in Marcos Juárez, Córdoba, Argentina. It produces indole-3-acetic acid in the presence of l-tryptophan as a precursor, grows at 20–38 °C (optimal 38 °C), and the cells are curved or spiral-shaped, with diameters ranging from 0.5–0.9 to 1.8–2.2 µm. They contain C16 : 0, C18 : 0 and C18 : 1  ω7c/ω6c as the main fatty acids. Phylogenetic analysis of its 16S rRNA gene sequence confirmed that this strain belongs to the genus Azospirillum , showing a close relationship with Azospirillum baldaniorum Sp245T, Azospirillum brasilense Sp7T and Azospirillum formosense CC-Nfb-7T. Housekeeping gene analysis revealed that Az39T, together with five strains of the genus (Az19, REC3, BR 11975, MTCC4035 and MTCC4036), form a cluster apart from A. baldaniorum Sp245T, A. brasilense Sp7T and A. formosense CC-Nfb-7T. Average nucleotide identity (ANI) and digital DNA–DNA hybridization (dDDH) between Az39T and the aforementioned type strains revealed values below 96 %, the circumscription limit for the species delineation (ANI: 95.3, 94.1 and 94.0 %; dDDH: 62.9, 56.3 and 55.6 %). Furthermore, a phylogeny evaluation of the core proteome, including 809 common shared proteins, showed an independent grouping of Az39T, Az19, REC3, BR 11975, MTCC4035 and MTCC4036. The G+C content in the genomic DNA of these six strains varied from 68.3 to 68.5 %. Based on the combined phylogenetic, genomic and phenotypic characterization presented here, we consider that strain Az39T, along with strains Az19, REC3, BR 11975, MTCC4035 and MTCC4036, are members of a new Azospirillum species, for which the name Azospirillum argentinense sp. nov. is proposed. The type strain is Az39T (=LBPCV39T=BR 148428T=CCCT 22.01T).

Comparative Genomics and Pathogenicity Analysis of Two Bacterial Symbionts of Entomopathogenic Nematodes: The Role of the GroEL Protein in Virulence

Bacteria of the genera Xenorhabdus and Photorhabdus are symbionts of entomopathogenic nematodes. Despite their close phylogenetic relationship, they show differences in their pathogenicity and virulence mechanisms in target insects. These differences were explored by the analysis of the pangenome, as it provides a framework for characterizing and defining the gene repertoire. We performed the first pangenome analysis of 91 strains of Xenorhabdus and Photorhabdus; the analysis showed that the Photorhabdus genus has a higher number of genes associated with pathogenicity. However, biological tests showed that whole cells of X. nematophila SC 0516 were more virulent than those of P. luminescens HIM3 when both were injected into G. mellonella larvae. In addition, we cloned and expressed the GroEL proteins of both bacteria, as this protein has been previously indicated to show insecticidal activity in the genus Xenorhabdus. Among these proteins, Cpn60-Xn was found to be the most toxic at all concentrations tested, with an LC50 value of 102.34 ng/larva. Sequence analysis suggested that the Cpn60-Xn toxin was homologous to Cpn60-Pl; however, Cpn60-Xn contained thirty-five differentially substituted amino acid residues that could be responsible for its insecticidal activity.

Pangenome Analysis of Plant Transcripts and Coding Sequences

The pangenome of a species is the sum of the genomes of its individuals. As coding sequences often represent only a small fraction of each genome, analyzing the pangene set can be a cost-effective strategy for plants with large genomes or highly heterozygous species. Here, we describe a step-by-step protocol to analyze plant pangene sets with the software GET_HOMOLOGUES-EST . After a short introduction, where the main concepts are illustrated, the remaining sections cover the installation and typical operations required to analyze and annotate pantranscriptomes and gene sets of plants. The recipes include instructions on how to call core and accessory genes, how to compute a presence-absence pangenome matrix, and how to identify and analyze private genes, present only in some genotypes. Downstream phylogenetic analyses are also discussed.

The open pan-genome architecture and virulence landscape of Mycobacterium bovis

Animal tuberculosis (TB) is an emergent disease caused by Mycobacterium bovis, one of the animal-adapted ecotypes of the Mycobacterium tuberculosis complex (MTC). In this work, whole-genome comparative analyses of 70 M. bovis were performed to gain insights into the pan-genome architecture. The comparison across M. bovis predicted genome composition enabled clustering into the core- and accessory-genome components, with 2736 CDS for the former, while the accessory moiety included 3897 CDS, of which 2656 are restricted to one/two genomes only. These analyses predicted an open pan-genome architecture, with an average of 32 CDS added by each genome and show the diversification of discrete M. bovis subpopulations supported by both core- and accessory-genome components. The functional annotation of the pan-genome classified each CDS into one or several COG (Clusters of Orthologous Groups) categories, revealing ‘transcription’ (total average CDSs, n=258), ‘lipid metabolism and transport’ (n=242), ‘energy production and conversion’ (n=214) and ‘unknown function’ (n=876) as the most represented. The closer analysis of polymorphisms in virulence-related genes in a restrict group of M. bovis from a multi-host system enabled the identification of clade-monomorphic non-synonymous SNPs, illustrating clade-specific virulence landscapes and correlating with disease severity. This first comparative pan-genome study of a diverse collection of M. bovis encompassing all clonal complexes indicates a high percentage of accessory genes and denotes an open, dynamic non-conservative pan-genome structure, with high evolutionary potential, defying the canons of MTC biology. Furthermore, it shows that M. bovis can shape its virulence repertoire, either by acquisition and loss of genes or by SNP-based diversification, likely towards host immune evasion, adaptation and persistence.

Genomic Stability and Genetic Defense Systems in Dolosigranulum pigrum, a Candidate Beneficial Bacterium from the Human Microbiome

Dolosigranulum pigrum is positively associated with indicators of health in multiple epidemiological studies of human nasal microbiota. Knowledge of the basic biology of D. pigrum is a prerequisite for evaluating its potential for future therapeutic use; however, such data are very limited. To gain insight into D. pigrum's chromosomal structure, pangenome, and genomic stability, we compared the genomes of 28 D. pigrum strains that were collected across 20 years. Phylogenomic analysis showed closely related strains circulating over this period and closure of 19 genomes revealed highly conserved chromosomal synteny. Gene clusters involved in the mobilome and in defense against mobile genetic elements (MGEs) were enriched in the accessory genome versus the core genome. A systematic analysis for MGEs identified the first candidate D. pigrum prophage and insertion sequence. A systematic analysis for genetic elements that limit the spread of MGEs, including restriction modification (RM), CRISPR-Cas, and deity-named defense systems, revealed strain-level diversity in host defense systems that localized to specific genomic sites, including one RM system hot spot. Analysis of CRISPR spacers pointed to a wealth of MGEs against which D. pigrum defends itself. These results reveal a role for horizontal gene transfer and mobile genetic elements in strain diversification while highlighting that in D. pigrum this occurs within the context of a highly stable chromosomal organization protected by a variety of defense mechanisms. IMPORTANCE Dolosigranulum pigrum is a candidate beneficial bacterium with potential for future therapeutic use. This is based on its positive associations with characteristics of health in multiple studies of human nasal microbiota across the span of human life. For example, high levels of D. pigrum nasal colonization in adults predicts the absence of Staphylococcus aureus nasal colonization. Also, D. pigrum nasal colonization in young children is associated with healthy control groups in studies of middle ear infections. Our analysis of 28 genomes revealed a remarkable stability of D. pigrum strains colonizing people in the United States across a 20-year span. We subsequently identified factors that can influence this stability, including genomic stability, phage predators, the role of MGEs in strain-level variation, and defenses against MGEs. Finally, these D. pigrum strains also lacked predicted virulence factors. Overall, these findings add additional support to the potential for D. pigrum as a therapeutic bacterium.

Electrolocation? The evidence for redox-mediated taxis in Shewanella oneidensis

Shewanella oneidensis is a dissimilatory metal reducing bacterium and model for extracellular electron transfer (EET), a respiratory mechanism in which electrons are transferred out of the cell. In the last 10 years, migration to insoluble electron acceptors for EET has been shown to be nonrandom and tactic, seemingly in the absence of molecular or energy gradients that typically allow for taxis. As the ability to sense, locate, and respire electrodes has applications in bioelectrochemical technology, a better understanding of taxis in S. oneidensis is needed. While the EET conduits of S. oneidensis have been studied extensively, its taxis pathways and their interplay with EET are not yet understood, making investigation into taxis phenomena nontrivial. Since S. oneidensis is a member of an EET-encoding clade, the genetic circuitry of taxis to insoluble acceptors may be conserved. We performed a bioinformatic analysis of Shewanella genomes to identify S. oneidensis chemotaxis orthologs conserved in the genus. In addition to the previously reported core chemotaxis gene cluster, we identify several other conserved proteins in the taxis signaling pathway. We present the current evidence for the two proposed models of EET taxis, "electrokinesis" and flavin-mediated taxis, and highlight key areas in need of further investigation.

A Novel Family of Acinetobacter Mega-Plasmids Are Disseminating Multi-Drug Resistance Across the Globe While Acquiring Location-Specific Accessory Genes

Acinetobacter species are emerging as major nosocomial pathogens, aided by their ability to acquire resistance to all classes of antibiotics. A key factor leading to their multi-drug resistance phenotypes is the acquisition of a wide variety of mobile genetic elements, particularly large conjugative plasmids. Here, we characterize a family of 21 multi-drug resistance mega-plasmids in 11 different Acinetobacter species isolated from various locations across the globe. The plasmid family exhibits a highly dynamic and diverse accessory genome, including 221 antibiotic resistance genes (ARGs) that confer resistance to 13 classes of antibiotics. We show that plasmids isolated within the same geographic region are often evolutionarily divergent members of this family based on their core-genome, yet they exhibit a more similar accessory genome. Individual plasmids, therefore, can disseminate to different locations around the globe, where they then appear to acquire diverse sets of accessory genes from their local surroundings. Further, we show that plasmids from several geographic regions were enriched with location-specific functional traits. Together, our findings show that these mega-plasmids can transmit across species boundaries, have the capacity for global dissemination, can accumulate a diverse suite of location-specific accessory genes, and can confer multi-drug resistance phenotypes of significant concern for human health. We therefore highlight this previously undescribed plasmid family as a serious threat to healthcare systems worldwide. These findings also add to the growing concern that mega-plasmids are key disseminators of antibiotic resistance and require global surveillance.

Phylogenomic Analysis of Secondary Metabolism in the Toxic Cyanobacterial Genera Anabaena, Dolichospermum and Aphanizomenon

Cyanobacteria produce an array of toxins that pose serious health risks to humans and animals. The closely related diazotrophic genera, Anabaena, Dolichospermum, and Aphanizomenon, frequently form poisonous blooms in lakes and brackish waters around the world. These genera form a complex now termed the Anabaena, Dolichospermum, and Aphanizomenon (ADA) clade and produce a greater array of toxins than any other cyanobacteria group. However, taxonomic confusion masks the distribution of toxin biosynthetic pathways in cyanobacteria. Here we obtained 11 new draft genomes to improve the understanding of toxin production in these genera. Comparison of secondary metabolite pathways in all available 31 genomes for these three genera suggests that the ability to produce microcystin, anatoxin-a, and saxitoxin is associated with specific subgroups. Each toxin gene cluster was concentrated or even limited to a certain subgroup within the ADA clade. Our results indicate that members of the ADA clade encode a variety of secondary metabolites following the phylogenetic clustering of constituent species. The newly sequenced members of the ADA clade show that phylogenetic separation of planktonic Dolichospermum and benthic Anabaena is not complete. This underscores the importance of taxonomic revision of Anabaena, Dolichospermum, and Aphanizomenon genera to reflect current phylogenomic understanding.

Genomic and Metabolomic Analyses of Natural Products in Nodularia spumigena Isolated from a Shrimp Culture Pond

The bloom-forming cyanobacterium Nodularia spumigena CENA596 encodes the biosynthetic gene clusters (BGCs) of the known natural products nodularins, spumigins, anabaenopeptins/namalides, aeruginosins, mycosporin-like amino acids, and scytonemin, along with the terpenoid geosmin. Targeted metabolomics confirmed the production of these metabolic compounds, except for the alkaloid scytonemin. Genome mining of N. spumigena CENA596 and its three closely related Nodularia strains—two planktonic strains from the Baltic Sea and one benthic strain from Japanese marine sediment—revealed that the number of BGCs in planktonic strains was higher than in benthic one. Geosmin—a volatile compound with unpleasant taste and odor—was unique to the Brazilian strain CENA596. Automatic annotation of the genomes using subsystems technology revealed a related number of coding sequences and functional roles. Orthologs from the Nodularia genomes are involved in the primary and secondary metabolisms. Phylogenomic analysis of N. spumigena CENA596 based on 120 conserved protein sequences positioned this strain close to the Baltic Nodularia. Phylogeny of the 16S rRNA genes separated the Brazilian CENA596 strain from those of the Baltic Sea, despite their high sequence identities (99% identity, 100% coverage). The comparative analysis among planktic Nodularia strains showed that their genomes were considerably similar despite their geographically distant origin.

Fifty Aureobasidium pullulans genomes reveal a recombining polyextremotolerant generalist

The black yeast Aureobasidium pullulans is a textbook example of a generalistic and ubiquitous fungus thriving in a wide variety of environments. To investigate whether A. pullulans is a true generalist, or alternatively, whether part of its versatility can be attributed to intraspecific specialization masked by cryptic diversification undetectable by traditional phylogenetic analyses, we sequenced and analysed the genomes of 50 strains of A. pullulans from different habitats and geographic locations. No population structure was observed in the sequenced strains. Decay of linkage disequilibrium over shorter physical distances (<100 bp) than in many sexually reproducing fungi indicates a high level of recombination in the species. A homothallic mating locus was found in all of the sequenced genomes. Aureobasidium pullulans appears to have a homogeneous population genetics structure, which is best explained by good dispersal and high levels of recombination. This means that A. pullulans is a true generalist that can inhabit different habitats without substantial specialization to any of these habitats at the genomic level. Furthermore, in the future, the high level of A. pullulans recombination can be exploited for the identification of genomic loci that are involved in the many biotechnologically useful traits of this black yeast.

Population Genomics of an Obligately Halophilic Basidiomycete Wallemia ichthyophaga

Background

Wallemia ichthyophaga is a highly specialized basidiomycetous fungus. It is one of the most halophilic fungi ever described, only able to grow at low water activity. This specialization is thought to explain why it is only rarely isolated from nature.

Results

Genomes of 21 W. ichthyophaga strains were sequenced with PE150 reads on BGISEQ500 platform. The genomes shared high similarity with the reference genome of the species, they were all smaller than 10 Mbp and had a low number of predicted genes. Groups of strains isolated in the same location encompassed clones as well as very divergent strains. There was little concordance between phylogenies of predicted genes. Linkage disequilibrium of pairs of polymorphic loci decayed relatively quickly as a function of distance between the loci (LD decay distance 1270 bp). For the first time a putative mating-type locus was identified in the genomes of W. ichthyophaga.

Conclusion

Based on the comparison of W. ichthyophaga genomes it appears that some phylogenetic lineages of the species can persist in the same location over at least several years. Apart from this, the differences between the strains do not reflect the isolation habitat or geographic location. Together with results supporting the existence of (sexual) recombination in W. ichthyophaga, the presented results indicate that strains of W. ichthyophaga can form a single recombining population even between different habitats and over large geographical distances.

Whole-genome sequence analysis of multidrug-resistant uropathogenic strains of Escherichia coli from Mexico

Background: Escherichia coli is the main bacterium associated with urinary tract infections (UTIs), including cystitis and pyelonephritis. Uropathogenic E. coli (UPEC) harbors numerous genes that encode diverse virulence factors contributing to its pathogenicity. The treatment of UTIs has become complicated due to the natural selection of E. coli strains that are multiresistant to several groups of antibiotics regularly used in clinical settings such as hospitals. Genomic reports of the global composition and distribution of the antibiotic resistance and virulence genes of these pathogenic strains are lacking in the Mexican population.

Purpose and methods: The aim of this study was to globally characterize the genomes of a group of UPEC strains by massive parallel sequencing to determine the prevalence and distribution of virulence and antibiotic resistance genes associated with different serotypes and phylogenetic groups.

Results: The strains exhibited 138-197 virulence genes and 29 antibiotic resistance genes related to antibiotics that are commonly used in clinical practice. 

Conclusions: These findings are relevant to the definition of new strategies for treating urinary tract infections in public hospitals and private practice. To further define the epidemiological distribution and composition of these virulence and antibiotic resistance genes, larger studies are needed.

Genomic Characterization of Jumbo Salmonella Phages That Effectively Target United Kingdom Pig-Associated Salmonella Serotypes

A common cause of human food poisoning is through ingestion of pork products contaminated with Salmonella spp. Worryingly multi-drug resistant (MDR) Salmonella strains have been isolated from pigs, which motivates the need for alternative antimicrobials. In this study isolation and characterization of 21 lytic Salmonella phages is described. All 21 phages, labeled as SPFM phages were shown to efficiently infect MDR Salmonella strains isolated from United Kingdom pigs and phages SPFM1, SPFM3, SPFM10, SPFM14, SPFM15, SPFM17, and SPFM19 could lyse 100% of strains tested. The phage genome sizes range from 233 to 242 Kb, which qualifies them as jumbo phages. All SPFM phage genomes are approximately 95% similar to each other by average nucleotide identity, they encode between 258-307 coding sequences and share 188 core genes. Phylogenetic analysis shows these phages are most similar to phages of the genus Seoulvirus and to further characterize phages within the genus, genes under positive selection were identified. Several of the genes under evolutionary selection pressure were predicted to encode for proteins that interact with bacteria. We describe the phenotypic and genetic characterization of this novel Salmonella phage set. As the phages efficiently kill MDR Salmonella strains, they may offer a promising alternative to antibiotics.

Genomic Evidence of Recombination in the Basidiomycete Wallemia mellicola

One of the most commonly encountered species in the small basidiomycetous sub-phylum Wallemiomycotina is Wallemia mellicola, a xerotolerant fungus with a widespread distribution. To investigate the population characteristics of the species, whole genomes of twenty-five strains were sequenced. Apart from identification of four strains of clonal origin, the distances between the genomes failed to reflect either the isolation habitat of the strains or their geographical origin. Strains from different parts of the world appeared to represent a relatively homogenous and widespread population. The lack of concordance between individual gene phylogenies and the decay of linkage disequilibrium indicated that W. mellicola is at least occasionally recombining. Two versions of a putative mating-type locus have been found in all sequenced genomes, each present in approximately half of the strains. W. mellicola thus appears to be capable of (sexual) recombination and shows no signs of allopatric speciation or specialization to specific habitats.

Insight into the genome and brackish water adaptation strategies of toxic and bloom-forming Baltic Sea Dolichospermum sp. UHCC 0315

The Baltic Sea is a shallow basin of brackish water in which the spatial salinity gradient is one of the most important factors contributing to species distribution. The Baltic Sea is infamous for its annual cyanobacterial blooms comprised of Nodularia spumigena, Aphanizomenon spp., and Dolichospermum spp. that cause harm, especially for recreational users. To broaden our knowledge of the cyanobacterial adaptation strategies for brackish water environments, we sequenced the entire genome of Dolichospermum sp. UHCC 0315, a species occurring not only in freshwater environments but also in brackish water. Comparative genomics analyses revealed a close association with Dolichospermum sp. UHCC 0090 isolated from a lake in Finland. The genome closure of Dolichospermum sp. UHCC 0315 unraveled a mixture of two subtypes in the original culture, and subtypes exhibited distinct buoyancy phenotypes. Salinity less than 3 g L-1 NaCl enabled proper growth of Dolichospermum sp. UHCC 0315, whereas growth was arrested at moderate salinity (6 g L-1 NaCl). The concentrations of toxins, microcystins, increased at moderate salinity, whereas RNA sequencing data implied that Dolichospermum remodeled its primary metabolism in unfavorable high salinity. Based on our results, the predicted salinity decrease in the Baltic Sea may favor toxic blooms of Dolichospermum spp.

Genome analysis of Salmonella enterica subsp. diarizonae isolates from invasive human infections reveals enrichment of virulence-related functions in lineage ST1256

BACKGROUND

Salmonella enterica subsp. diarizonae (IIIb) is frequently isolated from the environment, cold-blooded reptiles, sheep and humans; however only a few studies describe the isolation of this subspecies from invasive human infections. The factors contributing to this unusual behavior are currently unknown.

RESULTS

We report here the genome features of two diarizonae strains, SBO13 and SBO27, isolated from endocervical tissue collected post-abortion and from cerebrospinal fluid of a newborn child, respectively, in the city of Santa Cruz, Bolivia. Although isolated six years apart, SBO27 in 2008 and SBO13 in 2014, both strains belong to the same sequence type 1256 (ST1256) and show a high degree of genome conservation sharing more than 99% of their genes, including the conservation of a ~ 10 kb plasmid. A prominent feature of the two genomes is the presence of 24 genomic islands (GIs), in addition to 10 complete Salmonella pathogenicity islands (SPI) and fragments of SPI-7, a pathogenicity island first reported in the human-adapted serovar Typhi. Some of the GIs identified in SBO13 and SBO27 harbor genes putatively encoding auto-transporters involved in adhesion, lipopolysaccharide modifying enzymes, putative toxins, pili-related proteins, efflux pumps, and several putative membrane cation transport related-genes, among others. These two Bolivian isolates also share genes encoding the type-III secretion system effector proteins SseK2, SseK3 and SlrP with other diarizonae sequence types (ST) mainly-associated with infections in humans. The sseK2, sseK3 and slrP genes were either absent or showing frameshift mutations in a significant proportion of genomes from environmental diarizonae isolates.

CONCLUSIONS

The comparative genomic study of two diarizonae strains isolated in Bolivia from human patients uncovered the presence of many genes putatively related to virulence. The statistically-significant acquisition of a unique combination of these functions by diarizonae strains isolated from humans may have impacted the ability of these isolates to successfully infect the human host.

Evidence of the Red-Queen Hypothesis from Accelerated Rates of Evolution of Genes Involved in Biotic Interactions in Pneumocystis

Pneumocystis species are ascomycete fungi adapted to live inside the lungs of mammals. These ascomycetes show extensive stenoxenism, meaning that each species of Pneumocystis infects a single species of host. Here, we study the effect exerted by natural selection on gene evolution in the genomes of three Pneumocystis species. We show that genes involved in host interaction evolve under positive selection. In the first place, we found strong evidence of episodic diversifying selection in Major surface glycoproteins (Msg). These proteins are located on the surface of Pneumocystis and are used for host attachment and probably for immune system evasion. Consistent with their function as antigens, most sites under diversifying selection in Msg code for residues with large relative surface accessibility areas. We also found evidence of positive selection in part of the cell machinery used to export Msg to the cell surface. Specifically, we found that genes participating in glycosylphosphatidylinositol (GPI) biosynthesis show an increased rate of nonsynonymous substitutions (dN) versus synonymous substitutions (dS). GPI is a molecule synthesized in the endoplasmic reticulum that is used to anchor proteins to membranes. We interpret the aforementioned findings as evidence of selective pressure exerted by the host immune system on Pneumocystis species, shaping the evolution of Msg and several proteins involved in GPI biosynthesis. We suggest that genome evolution in Pneumocystis is well described by the Red-Queen hypothesis whereby genes relevant for biotic interactions show accelerated rates of evolution.

The Genomic Basis of Intrinsic and Acquired Antibiotic Resistance in the Genus Serratia

Serratia marcescens, a member of the Enterobacteriaceae family, was long thought to be a non-pathogenic bacterium prevalent in environmental habitats. Together with other members of this genus, it has emerged in recent years as an opportunistic nosocomial pathogen causing various types of infections. One important feature of pathogens belonging to this genus is their intrinsic and acquired resistance to a variety of antibiotic families, including β-lactam, aminoglycosides, quinolones and polypeptide antibiotics. The aim of this study was to elucidate which genes participate in the intrinsic and acquired antibiotic resistance of this genus in order to determine the Serratia genus resistome. We performed phylogenomic and comparative genomic analyses using 32 Serratia spp. genomes deposited in the NCBI GenBank from strains isolated from different ecological niches and different lifestyles. S. marcescens strain SmUNAM836, which was previously isolated from a Mexican adult with obstructive pulmonary disease, was included in this study. The results show that most of the antibiotic resistance genes (ARGs) were found on the chromosome, and to a lesser degree, on plasmids and transposons acquired through horizontal gene transfer. Four strains contained the gyrA point mutation in codon Ser83 that confers quinolone resistance. Pathogenic and environmental isolates presented a high number of ARGs, especially genes associated with efflux systems. Pathogenic strains, specifically nosocomial strains, presented more acquired resistance genes than environmental isolates. We may conclude that the environment provides a natural reservoir for antibiotic resistance, which has been underestimated in the medical field.

GET_PHYLOMARKERS, a Software Package to Select Optimal Orthologous Clusters for Phylogenomics and Inferring Pan-Genome Phylogenies, Used for a Critical Geno-Taxonomic Revision of the Genus Stenotrophomonas

The massive accumulation of genome-sequences in public databases promoted the proliferation of genome-level phylogenetic analyses in many areas of biological research. However, due to diverse evolutionary and genetic processes, many loci have undesirable properties for phylogenetic reconstruction. These, if undetected, can result in erroneous or biased estimates, particularly when estimating species trees from concatenated datasets. To deal with these problems, we developed GET_PHYLOMARKERS, a pipeline designed to identify high-quality markers to estimate robust genome phylogenies from the orthologous clusters, or the pan-genome matrix (PGM), computed by GET_HOMOLOGUES. In the first context, a set of sequential filters are applied to exclude recombinant alignments and those producing anomalous or poorly resolved trees. Multiple sequence alignments and maximum likelihood (ML) phylogenies are computed in parallel on multi-core computers. A ML species tree is estimated from the concatenated set of top-ranking alignments at the DNA or protein levels, using either FastTree or IQ-TREE (IQT). The latter is used by default due to its superior performance revealed in an extensive benchmark analysis. In addition, parsimony and ML phylogenies can be estimated from the PGM. We demonstrate the practical utility of the software by analyzing 170 Stenotrophomonas genome sequences available in RefSeq and 10 new complete genomes of Mexican environmental S. maltophilia complex (Smc) isolates reported herein. A combination of core-genome and PGM analyses was used to revise the molecular systematics of the genus. An unsupervised learning approach that uses a goodness of clustering statistic identified 20 groups within the Smc at a core-genome average nucleotide identity (cgANIb) of 95.9% that are perfectly consistent with strongly supported clades on the core- and pan-genome trees. In addition, we identified 16 misclassified RefSeq genome sequences, 14 of them labeled as S. maltophilia, demonstrating the broad utility of the software for phylogenomics and geno-taxonomic studies. The code, a detailed manual and tutorials are freely available for Linux/UNIX servers under the GNU GPLv3 license at https://github.com/vinuesa/get_phylomarkers. A docker image bundling GET_PHYLOMARKERS with GET_HOMOLOGUES is available at https://hub.docker.com/r/csicunam/get_homologues/, which can be easily run on any platform.

Genome sequences of Chlorella sorokiniana UTEX 1602 and Micractinium conductrix SAG 241.80: implications to maltose excretion by a green alga

Green algae represent a key segment of the global species capable of photoautotrophic-driven biological carbon fixation. Algae partition fixed-carbon into chemical compounds required for biomass, while diverting excess carbon into internal storage compounds such as starch and lipids or, in certain cases, into targeted extracellular compounds. Two green algae were selected to probe for critical components associated with sugar production and release in a model alga. Chlorella sorokiniana UTEX 1602 - which does not release significant quantities of sugars to the extracellular space - was selected as a control to compare with the maltose-releasing Micractinium conductrix SAG 241.80 - which was originally isolated from an endosymbiotic association with the ciliate Paramecium bursaria. Both strains were subjected to three sequencing approaches to assemble their genomes and annotate their genes. This analysis was further complemented with transcriptional studies during maltose release by M. conductrix SAG 241.80 versus conditions where sugar release is minimal. The annotation revealed that both strains contain homologs for the key components of a putative pathway leading to cytosolic maltose accumulation, while transcriptional studies found few changes in mRNA levels for the genes associated with these established intracellular sugar pathways. A further analysis of potential sugar transporters found multiple homologs for SWEETs and tonoplast sugar transporters. The analysis of transcriptional differences revealed a lesser and more measured global response for M. conductrix SAG 241.80 versus C. sorokiniana UTEX 1602 during conditions resulting in sugar release, providing a catalog of genes that might play a role in extracellular sugar transport.

Genome sequences of lower Great Lakes Microcystis sp. reveal strain-specific genes that are present and expressed in western Lake Erie blooms

Blooms of the potentially toxic cyanobacterium Microcystis are increasing worldwide. In the Laurentian Great Lakes they pose major socioeconomic, ecological, and human health threats, particularly in western Lake Erie. However, the interpretation of "omics" data is constrained by the highly variable genome of Microcystis and the small number of reference genome sequences from strains isolated from the Great Lakes. To address this, we sequenced two Microcystis isolates from Lake Erie (Microcystis aeruginosa LE3 and M. wesenbergii LE013-01) and one from upstream Lake St. Clair (M. cf aeruginosa LSC13-02), and compared these data to the genomes of seventeen Microcystis spp. from across the globe as well as one metagenome and seven metatranscriptomes from a 2014 Lake Erie Microcystis bloom. For the publically available strains analyzed, the core genome is ~1900 genes, representing ~11% of total genes in the pan-genome and ~45% of each strain's genome. The flexible genome content was related to Microcystis subclades defined by phylogenetic analysis of both housekeeping genes and total core genes. To our knowledge this is the first evidence that the flexible genome is linked to the core genome of the Microcystis species complex. The majority of strain-specific genes were present and expressed in bloom communities in Lake Erie. Roughly 8% of these genes from the lower Great Lakes are involved in genome plasticity (rapid gain, loss, or rearrangement of genes) and resistance to foreign genetic elements (such as CRISPR-Cas systems). Intriguingly, strain-specific genes from Microcystis cultured from around the world were also present and expressed in the Lake Erie blooms, suggesting that the Microcystis pangenome is truly global. The presence and expression of flexible genes, including strain-specific genes, suggests that strain-level genomic diversity may be important in maintaining Microcystis abundance during bloom events.

KinFin: Software for Taxon-Aware Analysis of Clustered Protein Sequences

The field of comparative genomics is concerned with the study of similarities and differences between the information encoded in the genomes of organisms. A common approach is to define gene families by clustering protein sequences based on sequence similarity, and analyze protein cluster presence and absence in different species groups as a guide to biology. Due to the high dimensionality of these data, downstream analysis of protein clusters inferred from large numbers of species, or species with many genes, is nontrivial, and few solutions exist for transparent, reproducible, and customizable analyses. We present KinFin, a streamlined software solution capable of integrating data from common file formats and delivering aggregative annotation of protein clusters. KinFin delivers analyses based on systematic taxonomy of the species analyzed, or on user-defined, groupings of taxa, for example, sets based on attributes such as life history traits, organismal phenotypes, or competing phylogenetic hypotheses. Results are reported through graphical and detailed text output files. We illustrate the utility of the KinFin pipeline by addressing questions regarding the biology of filarial nematodes, which include parasites of veterinary and medical importance. We resolve the phylogenetic relationships between the species and explore functional annotation of proteins in clusters in key lineages and between custom taxon sets, identifying gene families of interest. KinFin can easily be integrated into existing comparative genomic workflows, and promotes transparent and reproducible analysis of clustered protein data.

Pan genome and CRISPR analyses of the bacterial fish pathogen Moritella viscosa

BACKGROUND

Winter-ulcer Moritella viscosa infections continue to be a significant burden in Atlantic salmon (Salmo salar L.) farming. M. viscosa comprises two main clusters that differ in genetic variation and phenotypes including virulence. Horizontal gene transfer through acquisition and loss of mobile genetic elements (MGEs) is a major driving force of bacterial diversification. To gain insight into genomic traits that could affect sublineage evolution within this bacterium we examined the genome sequences of twelve M. viscosa strains. Matches between M. viscosa clustered, regularly interspaced, short palindromic, repeats and associated cas genes (CRISPR-Cas) were analysed to correlate CRISPR-Cas with adaptive immunity against MGEs.

RESULTS

The comparative genomic analysis of M. viscosa isolates from across the North Atlantic region and from different fish species support delineation of M. viscosa into four phylogenetic lineages. The results showed that M. viscosa carries two distinct variants of the CRISPR-Cas subtype I-F systems and that CRISPR features follow the phylogenetic lineages. A subset of the spacer content match prophage and plasmid genes dispersed among the M. viscosa strains. Further analysis revealed that prophage and plasmid-like element distribution were reflected in the content of the CRISPR-spacer profiles.

CONCLUSIONS

Our data suggests that CRISPR-Cas mediated interactions with MGEs impact genome properties among M. viscosa, and that patterns in spacer and MGE distributions are linked to strain relationships.

Analysis of Plant Pan-Genomes and Transcriptomes with GET_HOMOLOGUES-EST, a Clustering Solution for Sequences of the Same Species

The pan-genome of a species is defined as the union of all the genes and non-coding sequences found in all its individuals. However, constructing a pan-genome for plants with large genomes is daunting both in sequencing cost and the scale of the required computational analysis. A more affordable alternative is to focus on the genic repertoire by using transcriptomic data. Here, the software GET_HOMOLOGUES-EST was benchmarked with genomic and RNA-seq data of 19 Arabidopsis thaliana ecotypes and then applied to the analysis of transcripts from 16 Hordeum vulgare genotypes. The goal was to sample their pan-genomes and classify sequences as core, if detected in all accessions, or accessory, when absent in some of them. The resulting sequence clusters were used to simulate pan-genome growth, and to compile Average Nucleotide Identity matrices that summarize intra-species variation. Although transcripts were found to under-estimate pan-genome size by at least 10%, we concluded that clusters of expressed sequences can recapitulate phylogeny and reproduce two properties observed in A. thaliana gene models: accessory loci show lower expression and higher non-synonymous substitution rates than core genes. Finally, accessory sequences were observed to preferentially encode transposon components in both species, plus disease resistance genes in cultivated barleys, and a variety of protein domains from other families that appear frequently associated with presence/absence variation in the literature. These results demonstrate that pan-genome analyses are useful to explore germplasm diversity.