Core Gene Families Research Articles

Chromosome-level genome assembly of two cultivated Jujubes

Jujube (Ziziphus jujuba Mill.) is a valuable tree species with economic, edible, medicinal, and ecological conservation benefits. In this study, we used PacBio HiFi sequencing and Hi-C technology to assemble chromosome-level genomes of two cultivated Jujubes, namely ‘Lingwuchangzao’ and ‘Shiguang’. The genome sizes of ‘Lingwuchangzao’ and ‘Shiguang’ were 385.66 Mb and 394.12 Mb, respectively, with contig N50 sizes of 30.62 Mb and 4.30 Mb. These genomes contained 31,082 and 31,015 protein-coding genes, with repeat element contents of 42.11% and 42.33%, respectively. Phylogenetic analysis revealed that ‘Lingwuchangzao’ was sister to ‘Shiguang’ and followed by ‘Dongzao’. Additionally, comparative analysis of gene families among ‘Lingwuchangzao’, ‘Shiguang’, ‘Dongzao’, ‘Junzao’, and ‘Suanzao’ identified 15,988 (57.98%) core gene families, 11,191 (40.59%) disposable gene families, and 394 (1.43%) private gene families. Overall, the assembly of the genomes of these two cultivated Jujube varieties provides valuable genetic information for Jujube evolution, functional genomics research, and molecular breeding.

Comparative Genomics of Limosilactobacillus reuteri YLR001 Reveals Genetic Diversity and Probiotic Properties.

To gain deeper insights into the genomic characteristics of Limosilactobacillus reuteri (L. reuteri) YLR001 and uncover its probiotic properties, in the current study, a comprehensive analysis of its whole genome was conducted, explicitly exploring the genetic variations associated with different host organisms. The genome of YLR001 consisted of a circular 2,242,943 bp chromosome with a GC content of 38.84%, along with three circular plasmids (24,864, 38, 926, and 132,625 bp). Among the 2183 protein-coding sequences (CDSs), the specific genes associated with genetic adaptation and stress resistance were identified. We predicted the function of COG protein genes and analyzed the KEGG pathways. Comparative genome analysis revealed that the pan-genome contained 5207 gene families, including 475 core gene families and 941 strain-specific genes. Phylogenetic analysis revealed distinct host specificity among 20 strains of L. reuteri, highlighting substantial genetic diversity across different hosts. This study enhanced our comprehension of the genetic diversity of L. reuteri YLR001, demonstrated its potential probiotic characteristics, and established more solid groundwork for future applications.

A comprehensive comparative genomic analysis revealed that plant growth promoting traits are ubiquitous in strains of Stenotrophomonas.

Stenotrophomonas strains, which are often described as plant growth promoting (PGP) bacteria, are ubiquitous in many environments. A total of 213 genomes of strains of Stenotrophomonas were analyzed using comparative genomics to better understand the ecological roles of these bacteria in the environment. The pan-genome of the 213 strains of Stenotrophomonas consists of 27,186 gene families, including 710 core gene families, 11,039 unique genes and 15,437 accessory genes. Nearly all strains of Stenotrophomonas harbor the genes for GH3-family cellulose degradation and GH2- and GH31-family hemicellulose hydrolase, as well as intact glycolysis and tricarboxylic acid cycle pathways. These abilities suggest that the strains of this genus can easily obtain carbon and energy from the environment. The Stenotrophomonas strains can respond to oxidative stress by synthesizing catalase, superoxide dismutase, methionine sulfoxide reductase, and disulfide isomerase, as well as managing their osmotic balance by accumulating potassium and synthesizing compatible solutes, such as betaine, trehalose, glutamate, and proline. Each Stenotrophomonas strain also contains many genes for resistance to antibiotics and heavy metals. These genes that mediate stress tolerance increase the ability of Stenotrophomonas strains to survive in extreme environments. In addition, many functional genes related to attachment and plant colonization, growth promotion and biocontrol were identified. In detail, the genes associated with flagellar assembly, motility, chemotaxis and biofilm formation enable the strains of Stenotrophomonas to effectively colonize host plants. The presence of genes for phosphate-solubilization and siderophore production and the polyamine, indole-3-acetic acid, and cytokinin biosynthetic pathways confer the ability to promote plant growth. These strains can produce antimicrobial compounds, chitinases, lipases and proteases. Each Stenotrophomonas genome contained 1-9 prophages and 17-60 genomic islands, and the genes related to antibiotic and heavy metal resistance and the biosynthesis of polyamines, indole-3-acetic acid, and cytokinin may be acquired by horizontal gene transfer. This study demonstrates that strains of Stenotrophomonas are highly adaptable for different environments and have strong potential for use as plant growth-promoting bacteria.

Pan-Genome Analysis and Secondary Metabolic Pathway Mining of Biocontrol Bacterium Brevibacillus brevis

Brevibacillus brevis is one of the most common biocontrol strains with broad applications in the prevention and control of plant diseases and insect pests. In order to deepen our understanding of B. brevis genomes, describe their characteristics comprehensively, and mine secondary metabolites, we retrieved the genomic sequences of nine B. brevis strains that had been assembled into complete genomes from the NCBI database. These genomic sequences were analyzed using phylogenetic analysis software, pan-genome analysis software, and secondary metabolite mining software. Results revealed that the genome size of B. brevis strains ranged from 6.16 to 6.73 Mb, with GC content ranging from 47.0% to 54.0%. Phylogenetic analysis classified the nine B. brevis strains into three branches. The analyses of ANI and dDDH showed that B. brevis NEB573 had the potential to become a new species of Brevibacillus and needed further research in the future. The pan-genome analysis identified 10032 gene families, including 3257 core gene families, 3112 accessory gene families, and 3663 unique gene families. In addition, 123 secondary metabolite biosynthetic gene clusters of 20 classes were identified in the genomes of nine B. brevis strains. The major types of biosynthetic gene clusters were non-ribosomal peptide synthase (NRPS) and transAT polyketide synthase (transAT-PKS). Furthermore, a large number of untapped secondary metabolites were identified in B. brevis. In summary, this study elucidated the pan-genome characteristics of the biocontrol bacterium B. brevis and identified its secondary metabolites, providing valuable insights for its further development and utilization.

Genomic evolution and insights into agronomic trait innovations of Sesamum species

Sesame is an ancient oilseed crop with high oil content and quality. However, the evolutionary history and genetic mechanisms of its valuable agronomic traits remain unclear. Here, we report chromosome-scale genomes of cultivated sesame (Sesamum indicum L.) and six wild Sesamum species, representing all three karyotypes within this genus. Karyotyping and genome-based phylogenic analysis revealed the evolutionary route of Sesamum species from n = 13 to n = 16 and revealed that allotetraploidization occurred in the wild species Sesamum radiatum. Early divergence of the Sesamum genus (48.5–19.7 million years ago) during the Tertiary period and its ancient phylogenic position within eudicots were observed. Pan-genome analysis revealed 9164 core gene families in the 7 Sesamum species. These families are significantly enriched in various metabolic pathways, including fatty acid (FA) metabolism and FA biosynthesis. Structural variations in SiPT1 and SiDT1 within the phosphatidyl ethanolamine-binding protein gene family lead to the genomic evolution of plant-architecture and inflorescence-development phenotypes in Sesamum. A genome-wide association study (GWAS) of an interspecific population and genome comparisons revealed a long terminal repeat insertion and a sequence deletion in DIR genes of wild Sesamum angustifolium and cultivated sesame, respectively; both variations independently cause high susceptibility to Fusarium wilt disease. A GWAS of 560 sesame accessions combined with an overexpression study confirmed that the NAC1 and PPO genes play an important role in upregulating oil content of sesame. Our study provides high-quality genomic resources for cultivated and wild Sesamum species and insights that can improve molecular breeding strategies for sesame and other oilseed crops.

Pan-Chromosome and Comparative Analysis of Agrobacterium fabrum Reveal Important Traits Concerning the Genetic Diversity, Evolutionary Dynamics, and Niche Adaptation of the Species.

Agrobacterium fabrum has been critical for the development of plant genetic engineering and agricultural biotechnology due to its ability to transform eukaryotic cells. However, the gene composition, evolutionary dynamics, and niche adaptation of this species is still unknown. Therefore, we established a comparative genomic analysis based on a pan-chromosome data set to evaluate the genetic diversity of A. fabrum. Here, 25 A. fabrum genomes were selected for analysis by core genome phylogeny combined with the average nucleotide identity (ANI), amino acid identity (AAI), and in silico DNA-DNA hybridization (DDH) values. An open pan-genome of A. fabrum exhibits genetic diversity with variable accessorial genes as evidenced by a consensus pan-genome of 12 representative genomes. The genomic plasticity of A. fabrum is apparent in its putative sequences for mobile genetic elements (MGEs), limited horizontal gene transfer barriers, and potentially horizontally transferred genes. The evolutionary constraints and functional enrichment in the pan-chromosome were measured by the Clusters of Orthologous Groups (COG) categories using eggNOG-mapper software, and the nonsynonymous/synonymous rate ratio (dN/dS) was determined using HYPHY software. Comparative analysis revealed significant differences in the functional enrichment and the degree of purifying selection between the core genome and non-core genome. We demonstrate that the core gene families undergo stronger purifying selection but have a significant bias to contain one or more positively selected sites. Furthermore, although they shared similar genetic diversity, we observed significant differences between chromosome 1 (Chr I) and the chromid in their functional features and evolutionary constraints. We demonstrate that putative genetic elements responsible for plant infection, ecological adaptation, and speciation represent the core genome, highlighting their importance in the adaptation of A. fabrum to plant-related niches. Our pan-chromosome analysis of A. fabrum provides comprehensive insights into the genetic properties, evolutionary patterns, and niche adaptation of the species. IMPORTANCE Agrobacterium spp. live in diverse plant-associated niches such as soil, the rhizosphere, and vegetation, which are challenged by multiple stressors such as diverse energy sources, plant defenses, and microbial competition. They have evolved the ability to utilize diverse resources, escape plant defenses, and defeat competitors. However, the underlying genetic diversity and evolutionary dynamics of Agrobacterium spp. remain unexplored. We examined the phylogeny and pan-genome of A. fabrum to define intraspecies evolutionary relationships. Our results indicate an open pan-genome and numerous MGEs and horizontally transferred genes among A. fabrum genomes, reflecting the flexibility of the chromosomes and the potential for genetic exchange. Furthermore, we observed significant differences in the functional features and evolutionary constraints between the core and accessory genomes and between Chr I and the chromid, respectively.

The Pan-Genomic Analysis of Corynebacterium striatum Revealed its Genetic Characteristics as an Emerging Multidrug-Resistant Pathogen.

Corynebacterium striatum is a Gram-positive bacterium that is straight or slightly curved and non-spore-forming. Although it was originally believed to be a part of the normal microbiome of human skin, a growing number of studies have identified it as a cause of various chronic diseases, bacteremia, and respiratory infections. However, despite its increasing importance as a pathogen, the genetic characteristics of the pathogen population, such as genomic characteristics and differences, the types of resistance genes and virulence factors carried by the pathogen and their distribution in the population are poorly understood. To address these knowledge gaps, we conducted a pan-genomic analysis of 314 strains of C. striatum isolated from various tissues and geographic locations. Our analysis revealed that C. striatum has an open pan-genome, comprising 5692 gene families, including 1845 core gene families, 2362 accessory gene families, and 1485 unique gene families. We also found that C. striatum exhibits a high degree of diversity across different sources, but strains isolated from skin tissue are more conserved. Furthermore, we identified 53 drug resistance genes and 42 virulence factors by comparing the strains to the drug resistance gene database (CARD) and the pathogen virulence factor database (VFDB), respectively. We found that these genes and factors are widely distributed among C. striatum, with 77.7% of strains carrying 2 or more resistance genes and displaying primary resistance to aminoglycosides, tetracyclines, lincomycin, macrolides, and streptomycin. The virulence factors are primarily associated with pathogen survival within the host, iron uptake, pili, and early biofilm formation. In summary, our study provides insights into the population diversity, resistance genes, and virulence factors ofC. striatum from different sources. Our findings could inform future research and clinical practices in the diagnosis, prevention, and treatment of C. striatum-associated diseases.

Pan-Genomes Provide Insights into the Genetic Basis of Auricularia heimuer Domestication.

In order to reveal the genetic variation signals of Auricularia heimuer that have occurred during their domestication and to find potential functional gene families, we constructed a monokaryotic pan-genome of A. heimuer representing four cultivated strains and four wild strains. The pan-genome contained 14,089 gene families, of which 67.56% were core gene families and 31.88% were dispensable gene families. We screened substrate utilization-related genes such as the chitinase gene ahchi1 of the glycoside hydrolase (GH) 18 family and a carbohydrate-binding module (CBM)-related gene from the dispensable families of cultivated populations. The genomic difference in the ahchi1 gene between the wild and cultivated genomes was caused by a 33 kb presence/absence variation (PAV). The detection rate of the ahchi1 gene was 93.75% in the cultivated population, significantly higher than that in the wild population (17.39%), indicating that it has been selected in cultivated strains. Principal component analysis (PCA) of the polymorphic markers in fragments near the ahchi1 gene was enriched in cultivated strains, and this was caused by multiple independent instances of artificial selection. We revealed for the first time the genetic basis of the ahchi1 gene in domestication, thereby providing a foundation for elucidating the potential function of the ahchi1 gene in the breeding of A. heimuer.

Comparative Genomic Analysis of Vibrio cincinnatiensis Provides Insights into Genetic Diversity, Evolutionary Dynamics, and Pathogenic Traits of the Species.

Vibrio cincinnatiensis is a poorly understood pathogenic Vibrio species, and the underlying mechanisms of its genetic diversity, genomic plasticity, evolutionary dynamics, and pathogenicity have not yet been comprehensively investigated. Here, a comparative genomic analysis of V. cincinnatiensis was constructed. The open pan-genome with a flexible gene repertoire exhibited genetic diversity. The genomic plasticity and stability were characterized by the determinations of diverse mobile genetic elements (MGEs) and barriers to horizontal gene transfer (HGT), respectively. Evolutionary divergences were exhibited by the difference in functional enrichment and selective pressure between the different components of the pan-genome. The evolution on the Chr I and Chr II core genomes was mainly driven by purifying selection. Predicted essential genes in V. cincinnatiensis were mainly found in the core gene families on Chr I and were subject to stronger evolutionary constraints. We identified diverse virulence-related elements, including the gene clusters involved in encoding flagella, secretion systems, several pili, and scattered virulence genes. Our results indicated the pathogenic potential of V. cincinnatiensis and highlighted that HGT events from other Vibrio species promoted pathogenicity. This pan-genome study provides comprehensive insights into this poorly understood species from the genomic perspective.

Starch content changes and metabolism-related gene regulation of Chinese cabbage synergistically induced by Plasmodiophora brassicae infection.

Clubroot is one of the major diseases adversely affecting Chinese cabbage (Brassica rapa) yield and quality. To precisely characterize the Plasmodiophora brassicae infection of Chinese cabbage, we developed a dual fluorescent staining method for simultaneously examining the pathogen, cell structures, and starch grains. The number of starch (amylopectin) grains increased in B. rapa roots infected by P. brassicae, especially from 14 to 21 days after inoculation. Therefore, the expression levels of 38 core starch metabolism genes were investigated by quantitative real-time PCR. Most genes related to starch synthesis were up-regulated at 7 days after P. brassicae inoculation, whereas the expression levels of starch degradation-related genes were increased at 14 days after inoculation. Then, genes encoding the core enzymes involved in starch metabolism were investigated by assessing their chromosomal distributions, structures, duplication events, and synteny among Brassica species. Genome comparisons indicated that 38 non-redundant genes belonging to six core gene families related to starch metabolism are highly conserved among Arabidopsis thaliana, B. rapa, Brassica nigra, and Brassica oleracea. Previous genome sequencing projects have revealed that P. brassicae obtained host nutrients by manipulating plant metabolism. Starch may serve as a carbon source for P. brassicae colonization, as indicated by histological observations and transcriptomic analysis. Results of this study may elucidate the evolution and expression of core starch metabolism genes and provide researchers with novel insights into the pathogenesis of clubroot in B. rapa.

Roseomonas oleicola sp. nov., isolated from an oil production mixture in Yumen Oilfield, and emended description of Roseomonas frigidaquae.

A pink, ovoid-shaped, Gram-stain-negative, strictly aerobic and motile bacterial strain, designated ROY-5-3T, was isolated from an oil production mixture from Yumen Oilfield in PR China. The strain grew at 4-42 °C (optimum, 30 °C), at pH 5-10 (optimum, 7) and with 0-5 % (w/v) NaCl (optimum, 0%). The results of phylogenetic analysis based on 16S rRNA gene sequences indicated that ROY-5-3T belongs to the genus Roseomonas and shared the highest pairwise similarities with Roseomonas frigidaquae CW67T (98.1%), Roseomonas selenitidurans BU-1T (97.8%), Roseomonas tokyonensis K-20T (97.7%) and Roseomonas stagni HS-69T (97.3%). The average nucleotide identity and digital DNA-DNA hybridization values between ROY-5-3T and other related type strains of Roseomonas species were less than 84.08 and 28.60 %, respectively, both below the species delineation threshold. Pan-genomic analysis showed that the novel isolate ROY-5-3T shared 3265 core gene families with the four closely related type strains in Roseomonas, and the number of strain-specific gene families was 513. The major fatty acids were identified as summed feature 8 (C18 : 1 ω6c/C18 : 1 ω7c), summed feature 3 (C16 : 1 ω6c/C16 : 1 ω7c) and C16 : 0. Strain ROY-5-3T contained Q-10 as the main ubiquinone and the genomic DNA G+C content was 69.8 mol%. The major polar lipids were diphosphatidylglycerol, phosphatidylcholine, phosphatidylethanolamine and phosphatidylglycerol. Based on the phylogenetic, morphological, physiological, chemotaxonomic and genome analyses, strain ROY-5-3T represents a novel species of the genus Roseomonas for which the name Roseomonas oleicola sp. nov. is proposed. The type strain is ROY-5-3T (=CGMCC 1.13459T =KCTC 82484T).

Draft genome of the blister beetle, Epicauta chinensis

Existence of cantharidin (CTD) in blister beetles is a significant ecological adaptive mechanism that defends against predators and regulates courtship and mating behaviors. To better understand CTD biosynthetic information as well as its biology and pharmacology, we assembled a genome of 151.88 Mb for Epicauta chinensis using PacBio sequencing technology. Gene annotation yielded 249,238 repeats, 527 non-coding RNAs and 12,520 protein-coding genes. Compared to other 11 insects, expansions of gene families in E. chinensis for most core gene families likely associated with environmental adaptation, such as chemoreception, immunity, and detoxification.We further annotated P450s and immune-related genes, a total of 117 putative P450s comprising 7 CYP2, 67 CYP3, 36 CYP4, and 7 mitochondrial P450s and 281 immune-related genes were identified. Comparative analysis of the insect immune repertoires indicated presence of immune genes detected only from Coleopteran insects such as MD2-like. This suggested a lineage-specific gene evolution for Coleopteran insects. Based on the gene family evolution analysis, we identified two probable candidate genes including CYP4TT1 and phytanoyl-CoA dioxygenase for CTD biosynthesis. The high-quality reference genome of E. chinensis provides the genetic basis for further investigation of CTD biosynthesis and in-depth studies of the development and evolution of blister beetles.

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

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

Capnocytophaga periodontitidis sp. nov., isolated from subgingival plaque of periodontitis patient.

Two carbon dioxide-requiring, gliding, Gram-stain-negative strains, designated p1a2T and 051621, were isolated from subgingival plaque in association with severe periodontitis. The 16S rRNA gene sequence analysis revealed that they represented members of the genus Capnocytophaga and had less than 96.4 % pairwise similarity with species with validly published names in this genus. The whole-genome sequences of those strains had less than 91.9 % average nucleotide identity and 48.4 % digital DNA-DNA hybridization values with the other type strains of species of the genus Capnocytophaga, both below the species delineation threshold. The results of pan-genomic analysis indicated that p1a2T and 051621 shared 765 core gene families with the other ten species in this genus, and the numbers of strain-specific gene families were 493 and 455, respectively. The major fatty acids were iso-C15 : 0 and C16 : 0. A combination of phenotypic, chemotaxonomic, phylogenetic and genotypic data clearly indicate that p1a2T and 051621 should be considered to represent a novel species of the genus Capnocytophaga, for which the name Capnocytophaga periodontitidis sp. nov. is proposed. The type strain is p1a2T (=CGMCC 1.17337T=JCM 34126T).

Pseudomonas nanhaiensis sp. nov., a lipase-producing bacterium isolated from deep-sea sediment of the South China Sea.

A bacterial lipase producing bacterium, designated SCS 2-3, was isolated from deep-sea sediment of the South China Sea. Phylogenetic analysis based on the 16S rRNA sequence revealed that strain SCS2-3 belonged to the genus Pseudomonas and had 98.56% similarity to P. xinjiangensis NRRL B-51270T as the closest relative strain. MLSA using four protein-coding genes (dnaK, gyrA, recA, and rpoB) showed strain SCS 2-3 to form a separate branch. ANI and in silico DDH values between strain SCS 2-3 and related type strains of Pseudomonas were less than 81.51% and 23.80%, respectively. Genome comparison showed that strain SCS 2-3 shared 1875 core gene families with other eight closely related type strains in Pseudomonas, and the number of strain-unique genes was 263. Through gene annotations, genes related to lipase were found in the genome. Furthermore, a combination of phenotypic, chemotaxonomic, phylogenetic and genotypic data clearly indicated that strain SCS 2-3 represents a novel species of the genus Pseudomonas, for which the name Pseudomonas nanhaiensis sp. nov. is proposed. The type strain is SCS 2-3T (= GDMCC 1.2219T = JCM 34440T).

Discovery of core gene families associated with liver metastasis in colorectal cancer and regulatory roles in tumor cell immune infiltration

In this study, we aimed to uncover genes that drive the pathogenesis of liver metastasis in colorectal cancer (CRC), and identify effective genes that could serve as potential therapeutic targets for treating with colorectal liver metastasis patients based on two GEO datasets. Several bioinformatics approaches were implemented. First, differential expression analysis screened out key differentially expressed genes (DEGs) across the two GEO datasets. Based on gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses, we identified the enrichment functions and pathways of the DEGs that were associated with liver metastasis in CRC. Second, immune infiltration analysis identified key immune signature gene sets associated with CRC liver metastasis, among which two key immune gene families (CD and CCL) identified as key DEGs were filtered by protein-protein interaction (PPI) network. Some of the members in these gene families were associated with disease free survival (DFS) or overall survival (OS) in two subtypes of CRC, namely COAD and READ. Finally, functional enrichment analysis of the two gene families and their neighboring genes revealed that they were closely associated with cytokine, leukocyte proliferation and chemotaxis. These results are valuable in comprehending the pathogenesis of liver metastasis in CRC, and are of seminal importance in understanding the role of immune tumor infiltration in CRC. Our study also identified potentially effective therapeutic targets for liver metastasis in CRC including CCL20, CCL24 and CD70.

Azospirillum oleiclasticum sp. nov, a nitrogen-fixing and heavy oil degrading bacterium isolated from an oil production mixture of Yumen Oilfield

Two nitrogen-fixing and heavy oil degrading strains, designated RWY-5-1-1T and ROY-1-1-2, were isolated from an oil production mixture from Yumen Oilfield in China. The 16S rRNA gene sequence showed they belong to Azospirillum and have less than 96.1 % pairwise similarity with each species in this genus. The average nucleotide identity and digital DNA-DNA hybridization values between them and other type strains of Azospirillum species were less than 75.69 % and 22.0 %, respectively, both below the species delineation threshold. Pan-genomic analysis showed that the novel isolate RWY-5-1-1T shared 2145 core gene families with other type strains in Azospirillum, and the number of strain-specific gene families was 1623, almost two times more than the number known from other species. Furthermore, genes related to nitrogenase, hydrocarbon degradation and biosurfactant production were found in the isolates’ genomes. Also, this strain was capable of reducing acetylene to ethylene at a rate of 22nmol ethylene h-1 (108 cells) and degrading heavy oil at a rate of 36.2 %. The major fatty acids and polar lipids were summed feature 8 (C18:1ω7c/C18:1ω6c), and phosphatidylethanolamine, diphosphatidylglycerol, phosphatidylglycerol, and phosphatidylcholine. Furthermore, a combination of phenotypic, chemotaxonomic, phylogenetic and genotypic data clearly indicated that strains RWY-5-1-1T and ROY-1-1-2 represent a novel species, for which the name Azospirillum oleiclasticum sp. nov. is proposed. The type strain is RWY-5-1-1T (=CGMCC 1.13426T =KCTC 72259 T). Azospirillum novel strains with the ability of heavy oil degradation associated with the promotion of plant growth has never been reported to date.

Genomic signatures of globally enhanced gene duplicate accumulation in the megadiverse higher Diptera fueling intralocus sexual conflict resolution.

Gene duplication is an important source of evolutionary innovation. To explore the relative impact of gene duplication during the diversification of major insect model system lineages, we performed a comparative analysis of lineage-specific gene duplications in the fruit fly Drosophila melanogaster (Diptera: Brachycera), the mosquito Anopheles gambiae (Diptera: Culicomorpha), the red flour beetle Tribolium castaneum (Coleoptera), and the honeybee Apis mellifera (Hymenoptera). Focusing on close to 6,000 insect core gene families containing maximally six paralogs, we detected a conspicuously higher number of lineage-specific duplications in Drosophila (689) compared to Anopheles (315), Tribolium (386), and Apis (223). Based on analyses of sequence divergence, phylogenetic distribution, and gene ontology information, we present evidence that an increased background rate of gene duplicate accumulation played an exceptional role during the diversification of the higher Diptera (Brachycera), in part by providing enriched opportunities for intralocus sexual conflict resolution, which may have boosted speciation rates during the early radiation of the megadiverse brachyceran subclade Schizophora.

Pan-genomic analysis of the species Salmonella enterica: Identification of core essential and putative essential genes

BackgroundEssential genes are defined as the minimal gene set required to support bacterial life. In order to develop new antimicrobials to treat multidrug-resistant pathogens, such as serovars of Salmonella enterica, the identification of essential genes is crucial. MethodologyIn the present work, we hypothesize that essential genes within a group of evolutionary closely related organisms may be highly conserved. We, therefore, conducted an extensive comparative genomic analysis of 44 genome sequences representing 17 serovars of S. enterica to gain an improved understanding of conserved essential genes for its survival. ResultsPan-genome estimates indicate that the genus Salmonella displays an open pan-genome structure comprising a reservoir of 10,775 gene families. Of these, 2847, 4657, and 3271 constitute the core gene families (CGFs), dispensable gene families (DGFs), and strain-specific gene families (SSGFs), respectively. The pan-genome family tree based on the presence/absence of gene families is highly concordant with the 16S rRNA tree, though the former provides a more robust phylogenetic resolution. The Clusters of Orthologous Groups of proteins (COGs) database categorized the vast majority of the CGFs (40.9%) to metabolism, whereas a large proportion of the DGFs (70.6%) was uncharacterized. Homology analysis of the CGFs against the Database of essential genes (DEG) identified 1695 essential CGFs (E-CGFs). Of these, 687 are experimentally verified as essential in Salmonella, 1157 are identified in ≥2 species, 159 are conserved in ≥7 species, and 538 were present in at least one species. Thus, for the species, S. enterica 69%, 52%, and 31% of the genome are dedicated to the core, essential, and dispensable functions, respectively. ConclusionThe E-CGFs identified may serve as important targets for the development of novel antimicrobials, and their detailed analysis may shed new light on a better understanding of Salmonella's survival.

Reanalysis of Lactobacillus paracasei Lbs2 Strain and Large-Scale Comparative Genomics Places Many Strains into Their Correct Taxonomic Position.

Lactobacillus paracasei are diverse Gram-positive bacteria that are very closely related to Lactobacillus casei, belonging to the Lactobacillus casei group. Due to extreme genome similarities between L. casei and L. paracasei, many strains have been cross placed in the other group. We had earlier sequenced and analyzed the genome of Lactobacillus paracasei Lbs2, but mistakenly identified it as L. casei. We re-analyzed Lbs2 reads into a 2.5 MB genome that is 91.28% complete with 0.8% contamination, which is now suitably placed under L. paracasei based on Average Nucleotide Identity and Average Amino Acid Identity. We took 74 sequenced genomes of L. paracasei from GenBank with assembly sizes ranging from 2.3 to 3.3 MB and genome completeness between 88% and 100% for comparison. The pan-genome of 75 L. paracasei strains hold 15,945 gene families (21,5232 genes), while the core genome contained about 8.4% of the total genes (243 gene families with 18,225 genes) of pan-genome. Phylogenomic analysis based on core gene families revealed that the Lbs2 strain has a closer relationship with L. paracasei subsp. tolerans DSM20258. Finally, the in-silico analysis of the L. paracasei Lbs2 genome revealed an important pathway that could underpin the production of thiamin, which may contribute to the host energy metabolism.