Genomic and Transcriptomic Comparisons of the Twig Blight Pathogen, Passalora sequoiae, with Mycosphaerellaceae Foliar and Conifer Pathogens.

Passalora sequoiae is a foliar pathogen to conifer tree species. In this study, we conducted whole-genome and transcriptome analyses on isolates of P. sequoiae collected from symptomatic Leyland cypress leaves from a Christmas tree farm in Mississippi. The objectives for this research were to elucidate the pathogenicity mechanisms of P. sequoiae by characterizing the genome and transcriptome and possibly identify unique and shared predicted genes in comparison with non-conifer/canker and foliar pathogens in the family Mycosphaerellaceae. P. sequoiae was found to be similar to other foliar Mycosphaerellaceae pathogens and likely represents a hemibiotrophic lifestyle based on comparisons across pathogens. The genome and in planta transcriptome highlighted some unique features of P. sequoiae: the significant presence of chitin synthases and fructose-degrading carbohydrate-degrading enzymes, trans-AT PKS genes, and antibiotic gene clusters that were unique to P. sequoiae compared with the other Mycosphaerellaceae species genomes. Several transcripts that were highly expressed in planta were identified as effectors, yet the functions were not characterized. These targets provide ample resources to continue to characterize pathogen-conifer host interactions in conifer foliar pathogens. Furthermore, this research helps build genomic resources for an important plant pathogen on Leyland cypress that will further our ability to develop novel management practices that could begin with breeding for resistance.

ORTHOSCOPE*: A Phylogenetic Pipeline to Infer Gene Histories from Genome-Wide Data.

Comparative genome-scale analyses of protein-coding gene sequences are employed to examine evidence for whole-genome duplication and horizontal gene transfer. For this purpose, an orthogroup should be delineated to infer evolutionary history regarding each gene, and results of all orthogroup analyses need to be integrated to infer a genome-scale history. An orthogroup is a set of genes descended from a single gene in the last common ancestor of all species under consideration. However, such analyses confront several problems: 1) Analytical pipelines to infer all gene histories with methods comparing species and gene trees are not fully developed, and 2) without detailed analyses within orthogroups, evolutionary events of paralogous genes in the same orthogroup cannot be distinguished for genome-wide integration of results derived from multiple orthogroup analyses. Here I present an analytical pipeline, ORTHOSCOPE* (star), to infer evolutionary histories of animal/plant genes from genome-scale data. ORTHOSCOPE* estimates a tree for a specified gene, detects speciation/gene duplication events that occurred at nodes belonging to only one lineage leading to a species of interest, and then integrates results derived from gene trees estimated for all query genes in genome-wide data. Thus, ORTHOSCOPE* can be used to detect species nodes just after whole-genome duplications as a first step of comparative genomic analyses. Moreover, by examining the presence or absence of genes belonging to species lineages with dense taxon sampling available from the ORTHOSCOPE web version, ORTHOSCOPE* can detect genes lost in specific lineages and horizontal gene transfers. This pipeline is available at https://github.com/jun-inoue/ORTHOSCOPE_STAR.

Whole-Genome Relationships among Francisella Bacteria of Diverse Origins Define New Species and Provide Specific Regions for Detection.

Francisella tularensis is a highly virulent zoonotic pathogen that causes tularemia and, because of weaponization efforts in past world wars, is considered a tier 1 biothreat agent. Detection and surveillance of F. tularensis may be confounded by the presence of uncharacterized, closely related organisms. Through DNA-based diagnostics and environmental surveys, novel clinical and environmental Francisella isolates have been obtained in recent years. Here we present 7 new Francisella genomes and a comparison of their characteristics to each other and to 24 publicly available genomes as well as a comparative analysis of 16S rRNA and sdhA genes from over 90 Francisella strains. Delineation of new species in bacteria is challenging, especially when isolates having very close genomic characteristics exhibit different physiological features-for example, when some are virulent pathogens in humans and animals while others are nonpathogenic or are opportunistic pathogens. Species resolution within Francisella varies with analyses of single genes, multiple gene or protein sets, or whole-genome comparisons of nucleic acid and amino acid sequences. Analyses focusing on single genes (16S rRNA, sdhA), multiple gene sets (virulence genes, lipopolysaccharide [LPS] biosynthesis genes, pathogenicity island), and whole-genome comparisons (nucleotide and protein) gave congruent results, but with different levels of discrimination confidence. We designate four new species within the genus; Francisella opportunistica sp. nov. (MA06-7296), Francisella salina sp. nov. (TX07-7308), Francisella uliginis sp. nov. (TX07-7310), and Francisella frigiditurris sp. nov. (CA97-1460). This study provides a robust comparative framework to discern species and virulence features of newly detected Francisella bacteria. IMPORTANCE: DNA-based detection and sequencing methods have identified thousands of new bacteria in the human body and the environment. In most cases, there are no cultured isolates that correspond to these sequences. While DNA-based approaches are highly sensitive, accurately assigning species is difficult without known near relatives for comparison. This ambiguity poses challenges for clinical cases, disease epidemics, and environmental surveillance, for which response times must be short. Many new Francisella isolates have been identified globally. However, their species designations and potential for causing human disease remain ambiguous. Through detailed genome comparisons, we identified features that differentiate F. tularensis from clinical and environmental Francisella isolates and provide a knowledge base for future comparison of Francisella organisms identified in clinical samples or environmental surveys.

Analysis of the genomic sequences and metabolites of Serratia surfactantfaciens sp. nov. YD25T that simultaneously produces prodigiosin and serrawettin W2.

BACKGROUND: Gram-negative bacteria of the genus Serratia are potential producers of many useful secondary metabolites, such as prodigiosin and serrawettins, which have potential applications in environmental bioremediation or in the pharmaceutical industry. Several Serratia strains produce prodigiosin and serrawettin W1 as the main bioactive compounds, and the biosynthetic pathways are co-regulated by quorum sensing (QS). In contrast, the Serratia strain, which can simultaneously produce prodigiosin and serrawettin W2, has not been reported. This study focused on analyzing the genomic sequence of Serratia sp. strain YD25T isolated from rhizosphere soil under continuously planted burley tobacco collected from Yongding, Fujian province, China, which is unique in producing both prodigiosin and serrawettin W2. RESULTS: A hybrid polyketide synthases (PKS)-non-ribosomal peptide synthetases (NRPS) gene cluster putatively involved in biosynthesis of antimicrobial serrawettin W2 was identified in the genome of YD25T, and its biosynthesis pathway was proposed. We found potent antimicrobial activity of serrawettin W2 purified from YD25T against various pathogenic bacteria and fungi as well as antitumor activity against Hela cells. Subsequently, comparative genomic analyses were performed among a total of 133 Serratia species. The prodigiosin biosynthesis gene cluster in YD25T belongs to the type I pig cluster, which is the main form of pig-encoding genes existing in most of the pigmented Serratia species. In addition, a complete autoinducer-2 (AI-2) system (including luxS, lsrBACDEF, lsrGK, and lsrR) as a conserved bacterial operator is found in the genome of Serratia sp. strain YD25T. Phylogenetic analysis based on concatenated Lsr and LuxS proteins revealed that YD25T formed an independent branch and was clearly distant from the strains that solely produce either prodigiosin or serrawettin W2. The Fe (III) ion reduction assay confirmed that strain YD25T could produce an AI-2 signal molecule. Phylogenetic analysis using the genomic sequence of YD25T combined with phylogenetic and phenotypic analyses support this strain as a member of a novel and previously uncharacterized Serratia species. CONCLUSION: Genomic sequence and metabolite analysis of Serratia surfactantfaciens YD25T indicate that this strain can be further explored for the production of useful metabolites. Unveiling the genomic sequence of S. surfactantfaciens YD25T benefits the usage of this unique strain as a model system for studying the biosynthesis regulation of both prodigiosin and serrawettin W2 by the QS system.

Comparative study of the aldehyde dehydrogenase (ALDH) gene superfamily in the glycophyte Arabidopsis thaliana and Eutrema halophytes.

BACKGROUND AND AIMS: Stresses such as drought or salinity induce the generation of reactive oxygen species, which subsequently cause excessive accumulation of aldehydes in plant cells. Aldehyde dehydrogenases (ALDHs) are considered as 'aldehyde scavengers' to eliminate toxic aldehydes caused by oxidative stress. The completion of the genome sequencing projects of the halophytes Eutrema parvulum and E. salsugineum has paved the way to explore the relationships and the roles of ALDH genes in the glycophyte Arabidopsis thaliana and halophyte model plants. METHODS: Protein sequences of all plant ALDH families were used as queries to search E. parvulum and E. salsugineum genome databases. Evolutionary analyses compared the phylogenetic relationships of ALDHs from A. thaliana and Eutrema. Expression patterns of several stress-associated ALDH genes were investigated under different salt conditions using reverse transcription-PCR. Putative cis-elements in the promoters of ALDH10A8 from A. thaliana and E. salsugineum were compared in silico. KEY RESULTS: Sixteen and 17 members of ten ALDH families were identified from E. parvulum and E. salsugineum genomes, respectively. Phylogenetic analysis of ALDH protein sequences indicated that Eutrema ALDHs are closely related to those of Arabidopsis, and members within these species possess nearly identical exon-intron structures. Gene expression analysis under different salt conditions showed that most of the ALDH genes have similar expression profiles in Arabidopsis and E. salsugineum, except for ALDH7B4 and ALDH10A8. In silico analysis of promoter regions of ALDH10A8 revealed different distributions of cis-elements in E. salsugineum and Arabidopsis. CONCLUSIONS: Genomic organization, copy number, sub-cellular localization and expression profiles of ALDH genes are conserved in Arabidopsis, E. parvulum and E. salsugineum. The different expression patterns of ALDH7B4 and ALDH10A8 in Arabidopsis and E. salsugineum suggest that E. salsugineum uses modified regulatory pathways, which may contribute to salinity tolerance.

Decreased susceptibility to tigecycline in Acinetobacter baumannii mediated by a mutation in trm encoding SAM-dependent methyltransferase.

OBJECTIVES: Acinetobacter baumannii is an important opportunistic pathogen and multidrug-resistant isolate. Although tigecycline is a potent antibiotic for treating infections with multidrug-resistant isolates, resistance is becoming a problem. This study aimed to explore the mechanism of tigecycline resistance in A. baumannii. METHODS: A serial passage experiment was performed to collect isolates selected by tigecycline. The expression of efflux pumps was quantified in the final selected isolate, 19606-T8. The whole genome of 19606-T8 was sequenced and the putative mutations were confirmed using PCR and Sanger sequencing. A complementation experiment was performed to evaluate the contribution of the mutations to decreased susceptibility to tigecycline. The significance of a deletion mutation was further investigated in terms of growth rate and antibiotic susceptibilities. RESULTS: We collected serial isolates by selective pressure of tigecycline, and designated them 19606-T1 to 19606-T8. The efflux pumps AdeABC, AdeFGH and AdeIJK were not overexpressed in 19606-T8, which had decreased susceptibility to tigecycline. Isolate 19606-T8 carried one deletion mutation in trm and three non-synonymous substitutions in msbA (A84V), lolA (P91L) and filC (N168K). The deletion mutation in trm (encoding S-adenosyl-L-methionine-dependent methyltransferase) resulted in decreased susceptibility to tigecycline as well as to minocycline and doxycycline. In complementation experiments, the MICs of tigecycline, minocycline and doxycycline in a tigecycline-resistant isolate were restored by complementation with wild-type trm. CONCLUSIONS: Given that the deletion mutation in trm was associated with decreased susceptibility to tigecycline and that a wild-type trm plasmid could restore the susceptibility, trm is considered to play an important role in decreased susceptibility to tigecycline in A. baumannii.

New Genera and Species of Caulobacter and Brevundimonas Bacteriophages Provide Insights into Phage Genome Evolution.

Previous studies have identified diverse bacteriophages that infect Caulobacter vibrioides strain CB15 ranging from small RNA phages to four genera of jumbo phages. In this study, we focus on 20 bacteriophages whose genomes range from 40 to 60 kb in length. Genome comparisons indicated that these diverse phages represent six Caulobacter phage genera and one additional genus that includes both Caulobacter and Brevundimonas phages. Within species, comparisons revealed that both single base changes and inserted or deleted genetic material cause the genomes of closely related phages to diverge. Among genera, the basic gene order and the orientation of key genes were retained with most of the observed variation occurring at ends of the genomes. We hypothesize that the nucleotide sequences of the ends of these phage genomes are less important than the need to maintain the size of the genome and the stability of the corresponding mRNAs.

Genome-wide analysis reveals the emergence of multidrug resistant Stenotrophomonas acidaminiphila strain SINDOREI isolated from a patient with sepsis.

Stenotrophomonas acidaminiphila, the most recent reported species in genus Stenotrophomonas, is a relatively rare bacteria and is an aerobic, glucose non-fermentative, Gram-negative bacterium. However, little information of S. acidaminiphila is known to cause human infections. In this research, we firstly reported a multidrug-resistant strain S. acidaminiphila SINDOREI isolated from the blood of a patient with sepsis, who was dead of infection eventually. The whole genome of strain SINDOREI was sequenced, and genome comparisons were performed among six closely related S. acidaminiphila strains. The core genes (2,506 genes) and strain-specific genes were identified, respectively, to know about the strain-level diversity in six S. acidaminiphila stains. The presence of a unique gene (narG) and essential genes involved in biofilm formation in strain SINDOREI are important for the pathogenesis of infections. Strain SINDOREI was resistant to trimethoprim/sulfamethoxazole, ciprofloxacin, ofloxacin, cefepime, ceftazidime, and aztreonam. Several common and specific antibiotic resistance genes were identified in strain SINDOREI. The presence of two sul genes and exclusive determinants GES-1, aadA3, qacL, and cmlA5 is responsible for the resistance to multidrug. The virulence factors and resistance determinants can show the relationship between the phenotype and genotype and afford potential therapeutic strategies for infections.

Genome Features and Secondary Metabolites Biosynthetic Potential of the Class Ktedonobacteria.

The prevalence of antibiotic resistance and the decrease in novel antibiotic discovery in recent years necessitates the identification of potentially novel microbial resources to produce natural products. Ktedonobacteria, a class of deeply branched bacterial lineage in the ancient phylum Chloroflexi, are ubiquitous in terrestrial environments and characterized by their large genome size and complex life cycle. These characteristics indicate Ktedonobacteria as a potential active producer of bioactive compounds. In this study, we observed the existence of a putative "megaplasmid," multiple copies of ribosomal RNA operons, and high ratio of hypothetical proteins with unknown functions in the class Ktedonobacteria. Furthermore, a total of 104 antiSMASH-predicted putative biosynthetic gene clusters (BGCs) for secondary metabolites with high novelty and diversity were identified in nine Ktedonobacteria genomes. Our investigation of domain composition and organization of the non-ribosomal peptide synthetase and polyketide synthase BGCs further supports the concept that class Ktedonobacteria may produce compounds structurally different from known natural products. Furthermore, screening of bioactive compounds from representative Ktedonobacteria strains resulted in the identification of broad antimicrobial activities against both Gram-positive and Gram-negative tested bacterial strains. Based on these findings, we propose the ancient, ubiquitous, and spore-forming Ktedonobacteria as a versatile and promising microbial resource for natural product discovery.

Comparative Analysis Highlights Variable Genome Content of Wheat Rusts and Divergence of the Mating Loci.

Three members of the Puccinia genus, Pucciniatriticina (Pt), Pstriiformis f.sp. tritici (Pst), and Pgraminis f.sp. tritici (Pgt), cause the most common and often most significant foliar diseases of wheat. While similar in biology and life cycle, each species is uniquely adapted and specialized. The genomes of Pt and Pst were sequenced and compared to that of Pgt to identify common and distinguishing gene content, to determine gene variation among wheat rust pathogens, other rust fungi, and basidiomycetes, and to identify genes of significance for infection. Pt had the largest genome of the three, estimated at 135 Mb with expansion due to mobile elements and repeats encompassing 50.9% of contig bases; in comparison, repeats occupy 31.5% for Pst and 36.5% for Pgt We find all three genomes are highly heterozygous, with Pst [5.97 single nucleotide polymorphisms (SNPs)/kb] nearly twice the level detected in Pt (2.57 SNPs/kb) and that previously reported for Pgt Of 1358 predicted effectors in Pt, 784 were found expressed across diverse life cycle stages including the sexual stage. Comparison to related fungi highlighted the expansion of gene families involved in transcriptional regulation and nucleotide binding, protein modification, and carbohydrate degradation enzymes. Two allelic homeodomain pairs, HD1 and HD2, were identified in each dikaryotic Puccinia species along with three pheromone receptor (STE3) mating-type genes, two of which are likely representing allelic specificities. The HD proteins were active in a heterologous Ustilago maydis mating assay and host-induced gene silencing (HIGS) of the HD and STE3 alleles reduced wheat host infection.

Comparative Genomics of Two Sequential Candida glabrata Clinical Isolates.

Candida glabrata is an important fungal pathogen which develops rapid antifungal resistance in treated patients. It is known that azole treatments lead to antifungal resistance in this fungal species and that multidrug efflux transporters are involved in this process. Specific mutations in the transcriptional regulator PDR1 result in upregulation of the transporters. In addition, we showed that the PDR1 mutations can contribute to enhance virulence in animal models. In this study, we were interested to compare genomes of two specific C. glabrata-related isolates, one of which was azole susceptible (DSY562) while the other was azole resistant (DSY565). DSY565 contained a PDR1 mutation (L280F) and was isolated after a time-lapse of 50 d of azole therapy. We expected that genome comparisons between both isolates could reveal additional mutations reflecting host adaptation or even additional resistance mechanisms. The PacBio technology used here yielded 14 major contigs (sizes 0.18-1.6 Mb) and mitochondrial genomes from both DSY562 and DSY565 isolates that were highly similar to each other. Comparisons of the clinical genomes with the published CBS138 genome indicated important genome rearrangements, but not between the clinical strains. Among the unique features, several retrotransposons were identified in the genomes of the investigated clinical isolates. DSY562 and DSY565 each contained a large set of adhesin-like genes (101 and 107, respectively), which exceed by far the number of reported adhesins (63) in the CBS138 genome. Comparison between DSY562 and DSY565 yielded 17 nonsynonymous SNPs (among which the was the expected PDR1 mutation) as well as small size indels in coding regions (11) but mainly in adhesin-like genes. The genomes contained a DNA mismatch repair allele of MSH2 known to be involved in the so-called hyper-mutator phenotype of this yeast species and the number of accumulated mutations between both clinical isolates is consistent with the presence of a MSH2 defect. In conclusion, this study is the first to compare genomes of C. glabrata sequential clinical isolates using the PacBio technology as an approach. The genomes of these isolates taken in the same patient at two different time points exhibited limited variations, even if submitted to the host pressure.

Complete genome sequence of Vibrio anguillarum strain NB10, a virulent isolate from the Gulf of Bothnia.

Vibrio anguillarum causes a fatal hemorrhagic septicemia in marine fish that leads to great economical losses in aquaculture world-wide. Vibrio anguillarum strain NB10 serotype O1 is a Gram-negative, motile, curved rod-shaped bacterium, isolated from a diseased fish on the Swedish coast of the Gulf of Bothnia, and is slightly halophilic. Strain NB10 is a virulent isolate that readily colonizes fish skin and intestinal tissues. Here, the features of this bacterium are described and the annotation and analysis of its complete genome sequence is presented. The genome is 4,373,835 bp in size, consists of two circular chromosomes and one plasmid, and contains 3,783 protein-coding genes and 129 RNA genes.

Description of Bacillus toyonensis sp. nov., a novel species of the Bacillus cereus group, and pairwise genome comparisons of the species of the group by means of ANI calculations.

Strain BCT-7112(T) was isolated in 1966 in Japan from a survey designed to obtain naturally occurring microorganisms as pure cultures in the laboratory for use as probiotics in animal nutrition. This strain, which was primarily identified as Bacillus cereus var toyoi, has been in use for more than 30 years as the active ingredient of the preparation TOYOCERIN(®), an additive for use in animal nutrition (e.g. swine, poultry, cattle, rabbits and aquaculture). Despite the fact that the strain was initially classified as B. cereus, it showed significant genomic differences from the type strains of the B. cereus group that were large enough (ANI values below 92%) to allow it to be considered as a different species within the group. The polyphasic taxonomic study presented here provides sufficient discriminative parameters to classify BCT-7112(T) as a new species for which the name Bacillus toyonensis sp. nov. is proposed, with BCT-7112(T) (=CECT 876(T); =NCIMB 14858(T)) being designated as the type strain. In addition, a pairwise comparison between the available genomes of the whole B. cereus group by means of average nucleotide identity (ANI) calculations indicated that besides the eight classified species (including B. toyonensis), additional genomospecies could be detected, and most of them also had ANI values below 94%. ANI values were on the borderline of a species definition only in the cases of representatives of B. cereus versus B. thuringiensis, and B. mycoides and B. weihenstephanensis.