Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sequence analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13

FTO regulates myoblast proliferation by controlling CCND1 expression in an m6A-YTHDF2-dependent manner

Highlights: • Silencing FTO delays cell cycle and inhibits cell proliferation. • FTO regulates cell proliferation depending on its m{sup 6}A demethylation activity. • Silencing FTO suppresses CCND1 expression via YTHDF2-mediated mRNA degradation. • Knockdown YTHDF2 eliminates the cell cycle arrest in FTO-deficient cells. N6-Methyladenosine (m{sup 6}A) modification is the most abundant chemical modification in mRNA, and it participates in various biological processes, such as cell differentiation and proliferation. However, little is known about the function of m{sup 6}A demethylase fat mass and obesity-associated (FTO) in myoblast proliferation. Here, we demonstrated that knockdown of FTO can significantly inhibit myoblast proliferation and promote apoptosis. RNA sequencing analysis revealed that a lot of downregulated genes in FTO knockdown cells are associated with cell cycle and apoptosis. Furthermore, silencing FTO drastically decreased cyclin D1 (CCND1) expression through YTHDF2-mediated mRNA degradation, thereby delaying the progression of G1 phase, and leading to impaired myoblast proliferation. These findings unraveled that FTO regulates myoblast proliferation by controlling CCND1 expression in an m{sup 6}A-YTHDF2-dependent manner, which highlights the critical roles of m{sup 6}A modification in myoblast proliferation.

60 APPLIED LIFE SCIENCES↗

Streptomyces wadayamensis MppP Is a Pyridoxal 5'-Phosphate-Dependent L -Arginine α-Deaminase, γ-Hydroxylase in the Enduracididine Biosynthetic Pathway

l-Enduracididine (l-End) is a nonproteinogenic amino acid found in a number of bioactive peptides, including the antibiotics teixobactin, enduracidin, and mannopeptimycin. The potent activity of these compounds against antibiotic-resistant pathogens like MRSA and their novel mode of action have garnered considerable interest for the development of these peptides into clinically relevant antibiotics. Here, this goal has been hampered, at least in part, by the fact that l-End is difficult to synthesize and not currently commercially available. We have begun to elucidate the biosynthetic pathway of this unusual building block. In mannopeptimycin-producing strains, like Streptomyces wadayamensis, l-End is produced from l-Arg by the action of three enzymes: MppP, MppQ, and MppR. Herein, we report the structural and functional characterization of MppP. This pyridoxal 5'-phosphate (PLP)-dependent enzyme was predicted to be a fold type I aminotransferase on the basis of sequence analysis. We show that MppP is actually the first example of a PLP-dependent hydroxylase that catalyzes a reaction of l-Arg with dioxygen to yield a mixture of 2-oxo-4-hydroxy-5-guanidinovaleric acid and 2-oxo-5-guanidinovaleric acid in a 1.7:1 ratio. The structure of MppP with PLP bound to the catalytic lysine residue (Lys221) shows that, while the tertiary structure is very similar to those of the well-studied aminotransferases, there are differences in the arrangement of active site residues around the cofactor that likely account for the unusual activity of this enzyme. The structure of MppP with the substrate analogue d-Arg bound shows how the enzyme binds its substrate and indicates why d-Arg is not a substrate. On the basis of this work and previous work with MppR, we propose a plausible biosynthetic scheme for l-End.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Genomic insights into the origin, domestication and diversification of Brassica juncea

Despite early domestication around 3000 BC, the evolutionary history of the ancient allotetraploid species Brassica juncea (L.) Czern & Coss remains uncertain. Here, we report a chromosome-scale de novo assembly of a yellow-seeded B. juncea genome by integrating long-read and short-read sequencing, optical mapping and Hi-C technologies. Nuclear and organelle phylogenies of 480 accessions worldwide supported that B. juncea is most likely a single origin in West Asia, 8,000–14,000 years ago, via natural interspecific hybridization. Subsequently, new crop types evolved through spontaneous gene mutations and introgressions along three independent routes of eastward expansion. Selective sweeps, genome-wide trait associations and tissue-specific RNA-sequencing analysis shed light on the domestication history of flowering time and seed weight, and on human selection for morphological diversification in this versatile species. Our data provide a comprehensive insight into the origin and domestication and a foundation for genomics-based breeding of B. juncea.

59 BASIC BIOLOGICAL SCIENCES↗

Lyophilization of ASFV vaccine candidate ASFV-G-ΔI177L offers long term stability

Abstract For over a century African swine fever (ASF) has been causing outbreaks leading to devastating losses for the swine industry. The current pandemic of ASF has shown no signs of stopping and continues to spread causing outbreaks in additional countries. Currently control relies mostly on culling infected farms, and strict biosecurity procedures. Recently a vaccine, ASFV-G-ΔI177L was approved for use in Vietnam. In this study we evaluate the long-term stability of lyophilized ASFV-G-ΔI177L. Understanding the stability of different formulations of vaccines is information necessary for deployment of vaccines to ASF outbreak areas, particularly those that do not have a reliable well established cold chain to ensure conservation of vaccine quality. In this report, we determined that ASFV-G-ΔI177L, when lyophilized under specific conditions, is stable for up to one year at 4 °C, with similar vaccine titers after storage. Next-generation sequencing analysis also determined that lyophilization and long-term storage under these conditions had no effect on the genome of ASFV as the genome remained genetically identical to the original non-lyophilized form.

Science & Technology - Other Topics↗

Short-term early exposure to thirdhand cigarette smoke increases lung cancer incidence in mice

Exposure to thirdhand smoke (THS) is a recently described health concern that arises in many indoor environments. However, the carcinogenic potential of THS, a critical consideration in risk assessment, remains untested. Here we investigated the effects of short-term early exposure to THS on lung carcinogenesis in A/J mice. Forty weeks after THS exposure from 4 to 7 weeks of age, the mice had increased incidence of lung adenocarcinoma, tumor size and, multiplicity, compared with controls. In vitro studies using cultured human lung cancer cells showed that THS exposure induced DNA double-strand breaks and increased cell proliferation and colony formation. RNA sequencing analysis revealed that THS exposure induced endoplasmic reticulum stress and activated p53 signaling. Activation of the p53 pathway was confirmed by an increase in its targets p21 and BAX. These data indicate that early exposure to THS is associated with increased lung cancer risk.

59 BASIC BIOLOGICAL SCIENCES↗

Identification and grafting of a unique peptide-binding site in the Fab framework of monoclonal antibodies

Capitalizing on their extraordinary specificity, monoclonal antibodies (mAbs) have become one of the most reengineered classes of biological molecules. A major goal in many of these engineering efforts is to add new functionality to the parental mAb, including the addition of cytotoxins and imaging agents for medical applications. Herein, we present a unique peptide-binding site within the central cavity of the fragment antigen binding framework region of the chimeric, anti-epidermal growth factor receptor mAb cetuximab. We demonstrate through diffraction methods, biophysical studies, and sequence analysis that this peptide, a meditope, has moderate affinity for the Fab, is specific to cetuximab (i.e., does not bind to human IgGs), and has no significant effect on antigen binding. We further demonstrate by diffraction studies and biophysical methods that the meditope binding site can be grafted onto the anti-human epidermal growth factor receptor 2 mAb trastuzumab, and that the antigen binding affinity of the grafted trastuzumab is indistinguishable from the parental mAb. Lastly, we demonstrate a bivalent meditope variant binds specifically and stably to antigen-bearing cells only in the presence of the meditope-enabled mAbs. Collectively, this finding and the subsequent characterization and engineering efforts indicate that this unique interface could serve as a noncovalent “linker” for any meditope-enabled mAb with applications in multiple mAb-based technologies including diagnostics, imaging, and therapeutic delivery.

59 BASIC BIOLOGICAL SCIENCES↗

Structural basis of recognition of farnesylated and methylated KRAS4b by PDEδ

Farnesylation and carboxymethylation of KRAS4b (Kirsten rat sarcoma isoform 4b) are essential for its interaction with the plasma membrane where KRAS-mediated signaling events occur. Phosphodiesterase-δ (PDEδ) binds to KRAS4b and plays an important role in targeting it to cellular membranes. We solved structures of human farnesylated–methylated KRAS4b in complex with PDEδ in two different crystal forms. In these structures, the interaction is driven by the C-terminal amino acids together with the farnesylated and methylated C185 of KRAS4b that binds tightly in the central hydrophobic pocket present in PDEδ. In crystal form II, we see the full-length structure of farnesylated–methylated KRAS4b, including the hypervariable region. Crystal form I reveals structural details of farnesylated–methylated KRAS4b binding to PDEδ, and crystal form II suggests the potential binding mode of geranylgeranylated–methylated KRAS4b to PDEδ. We identified a 5-aa-long sequence motif (Lys-Ser-Lys-Thr-Lys) in KRAS4b that may enable PDEδ to bind both forms of prenylated KRAS4b. Structure and sequence analysis of various prenylated proteins that have been previously tested for binding to PDEδ provides a rationale for why some prenylated proteins, such as KRAS4a, RalA, RalB, and Rac1, do not bind to PDEδ. Comparison of all four available structures of PDEδ complexed with various prenylated proteins/peptides shows the presence of additional interactions due to a larger protein–protein interaction interface in KRAS4b–PDEδ complex. Here, this interface might be exploited for designing an inhibitor with minimal off-target effects.

59 BASIC BIOLOGICAL SCIENCES↗

Protein cost minimization promotes the emergence of coenzyme redundancy

Metabolism relies on a small class of molecules (coenzymes) that serve as universal donors and acceptors of key chemical groups and electrons. Although metabolic networks crucially depend on structurally redundant coenzymes [e.g., NAD(H) and NADP(H)] associated with different enzymes, the criteria that led to the emergence of this redundancy remain poorly understood. Our combination of modeling and structural and sequence analysis indicates that coenzyme redundancy may not be essential for metabolism but could rather constitute an evolved strategy promoting efficient usage of enzymes when biochemical reactions are near equilibrium. Our work suggests that early metabolism may have operated with fewer coenzymes and that adaptation for metabolic efficiency may have driven the rise of coenzyme diversity in living systems.

59 BASIC BIOLOGICAL SCIENCES↗

Structural insights into the polyphyletic origins of glycyl tRNA synthetases

Glycyl tRNA synthetase (GlyRS) provides a unique case among class II aminoacyl tRNA synthetases, with two clearly widespread types of enzymes: a dimeric (α 2 ) species present in some bacteria, archaea, and eukaryotes; and a heterotetrameric form (α 2 β 2 ) present in most bacteria. Although the differences between both types of GlyRS at the anticodon binding domain level are evident, the extent and implications of the variations in the catalytic domain have not been described, and it is unclear whether the mechanism of amino acid recognition is also dissimilar. Here, we show that the α-subunit of the α 2 β 2 GlyRS from the bacterium Aquifex aeolicus is able to perform the first step of the aminoacylation reaction, which involves the activation of the amino acid with ATP. The crystal structure of the α-subunit in the complex with an analog of glycyl adenylate at 2.8 Å resolution presents a conformational arrangement that properly positions the cognate amino acid. This work shows that glycine is recognized by a subset of different residues in the two types of GlyRS. Furthermore, a structural and sequence analysis of class II catalytic domains shows that bacterial GlyRS is closely related to alanyl tRNA synthetase, which led us to define a new subclassification of these ancient enzymes and to propose an evolutionary path of α 2 β 2 GlyRS, convergent with α 2 GlyRS and divergent from AlaRS, thus providing a possible explanation for the puzzling existence of two proteins sharing the same fold and function but not a common ancestor.

59 BASIC BIOLOGICAL SCIENCES↗

Validation of Continuous-Energy ENDF/B-VIII.0 16 O, 56 Fe, and 63,65 Cu Cross Sections for Nuclear Criticality Safety Applications

Recently completed cross-section evaluations sponsored in part by the Nuclear Criticality Safety Program were incorporated into the 2018 release of the ENDF/B-VIII.0 cross-section library. Evaluated isotopes of interest to the nuclear data and criticality safety community include 16 O, 56 Fe, and 63,65 Cu. For performance validation, benchmark models defined in the International Criticality Safety Benchmark Evaluation Project Handbook were selected based on energy-integrated k eff sensitivities to total cross sections of interest and compared with experimental values. Of the 102 benchmark configurations that were utilized, 63 are sensitive to 16 O, 32 sensitive to 63,65 Cu, and 25 sensitive to 56 Fe. Selected benchmarks were modeled in SCALE 6.2.3 Criticality Safety Analysis Sequence (CSAS) continuous-energy Monte Carlo k eff calculations with ENDF/B-VII.1, with a hybrid ENDF/B-VII.1 with ENDF/B-VIII.0 data substituted for individual isotopes of interest, and with ENDF/B-VIII.0. ENDF/B-VIII.0 showed improved agreement with experimental k eff for 56 Fe, 63 Cu, elemental copper, and full library substitution while producing lessened agreement for 16 O and 65 Cu. Furthermore, with full library and isotope-specific ENDF/B-VIII.0 performance, a best-case ENDF library was formed by excluding underperforming isotopes’ ENDF/B-VIII.0 data, reverting 16 O and 65 Cu cross sections to ENDF/B-VII.1. This resulted in the average relative deviation between calculated and experimental data improving from 1.45σ for the ENDF/B-VIII.0 library to 1.32σ for the best-case library, relative to benchmark uncertainty.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Continuous-Energy ENDF/B-VIII.0 Cross Section and SCALE 6.2.4 Performance for Nuclear Criticality Safety Applications: 1 H, C, 58,60 Ni, 182,183,184,186 W, 235,238 U, 239 Pu

As part of the nuclear data evaluation and validation cycle, the ENDF/B-VIII.0 cross-section library released in 2018 requires testing to determine areas of improvement and deterioration. Previous work by the authors investigated the performance of 16 O, 56 Fe, and 63,65 Cu cross sections, with this study acting as an extension of the prior work. In addition to the isotopes and nuclear criticality safety benchmarks of interest to the prior work, benchmarks from the International Criticality Safety Benchmark Evaluation Project Handbook were selected for their k eff sensitivity to 1 H, C, 58,60 Ni, 182,183,184,186 W, 235,238U, or 239 Pu cross sections and were modeled in the SCALE code system maintained by Oak Ridge National Laboratory. In total, 253 benchmark configurations were selected for their sensitivities and modeled using SCALE 6.2.4 Criticality Safety Analysis Sequences (CSAS) continuous-energy Monte Carlo k eff calculations. This collection includes and expands upon the 99 benchmarks in the prior work. The AMPX-processed ENDF/B-VIII.0 library was decomposed into individual ENDF/B-VIII.0 datum libraries for each isotope of interest. Doing so allowed for the individual substitution of an ENDF/B-VIII.0 cross section in the place of ENDF/B-VII.1, determining isotope-specific effects of ENDF/B-VIII.0 relative to ENDF/B-VII.1. Full library calculations with entirely ENDF/B-VII.1 data or entirely ENDF/B-VIII.0 data were also executed. As a measure of performance, the average relative deviation was determined as the ratio of the deviation between calculated and experimental keff to the propagated calculational and experimental uncertainty. With calculated full library and isotope-specific ENDF/B-VIII.0 k eff ’s, an optimized combination of data libraries was estimated and confirmed with SCALE calculations. This showed that reverting 239 Pu, 58 Ni, 16 O, and 65 Cu cross sections to ENDF/B-VII.1 resulted in improved performance relative to the full ENDF/B-VIII.0 library. Across all 253 benchmarks, the average relative deviation was 1.29σ for the full ENDF/B-VII.1 library, 1.17σ for the full ENDF/B-VIII.0 library, and 0.97σ for the optimized combination. The reversion of 239 Pu, 58 Ni, 16 O, and 65 Cu cross sections to ENDF/B-VII.1 in the 99 benchmarks of the prior work resulted in further improved experimental agreement compared to the previously reported improvement from 16 O and 65 Cu alone. Therefore, it is suggested that applications with significant sensitivities to 239 Pu, 58 Ni, 16 O, and 65 Cu consider their choice of nuclear data library.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Hypermut 3: identifying specific mutational patterns in a defined nucleotide context that allows multistate characters

Abstract Motivation The detection of APOBEC3F- and APOBEC3G-induced mutations in virus sequences is useful for identifying hypermutated sequences. These sequences are not representative of viral evolution and can therefore alter the results of downstream sequence analyses if included. We previously published the software Hypermut, which detects hypermutation events in sequences relative to a reference. Two versions of this method are available as a webtool. Neither of these methods consider multistate characters or gaps in the sequence alignment. Results Here, we present an updated, user-friendly web and command-line version of Hypermut with functionality to handle multistate characters and gaps in the sequence alignment. This tool allows for straightforward integration of hypermutation detection into sequence analysis pipelines. As with the previous tool, while the main purpose is to identify G to A hypermutation events, any mutational pattern and context can be specified. Availability and implementation Hypermut 3 is written in Python 3. It is available as a command-line tool at https://github.com/MolEvolEpid/hypermut3 and as a webtool at https://www.hiv.lanl.gov/content/sequence/HYPERMUT/hypermutv3.html.

59 BASIC BIOLOGICAL SCIENCES↗

Draft genome sequences of strains CBS6241 and CBS6242 of the basidiomycetous yeast Filobasidium floriforme

The Tremellomycetes are a species-rich group within the basidiomycete fungi; however, most analyses of this group to date have focused on pathogenic Cryptococcus species within the order Tremellales. Recent genome-assisted studies of other Tremellomycetes have identified interesting features with respect to biotechnological applications as well as the evolution of genes involved in mating and sexual development. Here, we report genome sequences of two strains of Filobasidium floriforme, a species from the order Filobasidiales, which branches basally to the Tremellales, Trichosporonales, and Holtermanniales. The assembled genomes of strains CBS6241 and CBS6242 are 27.4 Mb and 26.4 Mb in size, respectively, with 8314 and 7695 predicted protein-coding genes. Overall sequence identity at nucleic acid level between the strains is 97%. Among the predicted genes are pheromone precursor and pheromone receptor genes as well as two genes encoding homedomain (HD) transcription factors, which are predicted to be part of the mating type (MAT) locus. Sequence analysis indicates that CBS6241 and CBS6242 carry different alleles for both the pheromone/receptor genes as well as the HD transcription factors. Orthology inference identified 1482 orthogroups exclusively found in F. floriforme, some of which were involved in carbohydrate transport and metabolism. Subsequent CAZyme repertoire characterization identified 267 and 247 enzymes for CBS6241 and CBS6242, respectively, the second highest number of CAZymes among the analyzed Tremellomycete species. In addition, F. floriforme contains five CAZymes absent in other species and several plant-cell-wall degrading CAZymes with the highest copy number in Tremellomycota, indicating the biotechnological potential of this species.

59 BASIC BIOLOGICAL SCIENCES↗

Comparative Genomics of Saccharomyces cerevisiae Natural Isolates for Bioenergy Production

Lignocellulosic plant material is a viable source of biomass to produce alternative energy including ethanol and other biofuels. However, several factors—including toxic by products from biomass pretreatment and poor fermentation of xylose and other pentose sugars—currently limit the efficiency of microbial biofuel production. To begin to understand the genetic basis of desirable traits, we characterized three strains of Saccharomyces cerevisiae with robust growth in a pretreated lignocellulosic hydrolysate or tolerance to stress conditions relevant to industrial biofuel production, through genome and transcriptome sequencing analysis. All stress resistant strains were highly mosaic, suggesting that genetic admixture may contribute to novel allele combinations underlying these phenotypes. Strain-specific gene sets not found in the lab strain were functionally linked to the tolerances of particular strains. Furthermore, genes with signatures of evolutionary selection were enriched for functional categories important for stress resistance and included stress-responsive signaling factors. Comparison of the strains’ transcriptomic responses to heat and ethanol treatment—two stresses relevant to industrial bioethanol production—pointed to physiological processes that were related to particular stress resistance profiles. Many of the genotype-by-environment expression responses occurred at targets of transcription factors with signatures of positive selection, suggesting that these strains have undergone positive selection for stress tolerance. Our results generate new insights into potential mechanisms of tolerance to stresses relevant to biofuel production, including ethanol and heat, present a backdrop for further engineering, and provide glimpses into the natural variation of stress tolerance in wild yeast strains.

59 BASIC BIOLOGICAL SCIENCES↗

Occurrence and molecular characterization of Enterocytozoon bieneusi in wild and domestic animal species in Portugal

Abstract The phylum Microsporidia encompasses a diverse group of obligate, intracellular, and spore-forming organisms able to infect a wide range of animal hosts. Among them, Enterocytozoon bieneusi is the most frequently reported species in humans and animals. Little is known about the presence and epidemiology of E. bieneusi in wildlife. We investigated E. bieneusi occurrence and genetic diversity in wild and domestic mammals, through molecular-detection methods, from different regions across Portugal. A total of 756 samples were collected from 288, 242, and 226 wild carnivores, wild ungulates, and domestic animals, respectively. Overall, eight specimens were E. bieneusi-positive (1.1%, 8/756) obtained from five wild (Iberian lynx, Iberian wolf, red fox, stone marten, and wild boar) and one domestic (sheep) host. Nucleotide sequence analysis identified four genotypes of E. bieneusi, Type IV, Wildboar3, BEB6, and PtEbIX. Three of those genotypes belong to Groups 1 (Type IV and Wildboar3) and 2 (BEB6), which are known to contain genotypes capable of infecting a variety of hosts, including humans, highlighting their public health importance. PtEbIX belongs to the dog-specific Group 11. This study represents the first, largest, and most comprehensive molecular-based epidemiology survey carried out in Portugal in wild and domestic animals to date and the first worldwide identification of E. bieneusi in wolf species. Our study showed that wild carnivores and ungulates may act as reservoirs of zoonotic genotypes of E. bieneusi, establishing their role in maintaining the sylvatic cycle of this parasite while representing a potential source of infection for humans and domestic animals.

Figueiredo, Ana M. (ORCID:0000000226236340)↗

Microhomology-mediated end joining is activated in irradiated human cells due to phosphorylation-dependent formation of the XRCC1 repair complex

Microhomology-mediated end joining (MMEJ), an error-prone pathway for DNA double-strand break (DSB) repair, is implicated in genomic rearrangement and oncogenic transformation; however, its contribution to repair of radiation-induced DSBs has not been characterized. We used recircularization of a linearized plasmid with 3'-P-blocked termini, mimicking those at X-ray-induced strand breaks, to recapitulate DSB repair via MMEJ or nonhomologous end-joining (NHEJ). Sequence analysis of the circularized plasmids allowed measurement of relative activity of MMEJ versus NHEJ. While we predictably observed NHEJ to be the predominant pathway for DSB repair in our assay, MMEJ was significantly enhanced in preirradiated cells, independent of their radiation-induced arrest in the G2/M phase. MMEJ activation was dependent on XRCC1 phosphorylation by casein kinase 2 (CK2), enhancing XRCC1’s interaction with the end resection enzymes MRE11 and CtIP. Both endonuclease and exonuclease activities of MRE11 were required for MMEJ, as has been observed for homology-directed DSB repair (HDR). Furthermore, the XRCC1 co-immunoprecipitate complex (IP) displayed MMEJ activity in vitro, which was significantly elevated after irradiation. Our studies thus suggest that radiation-mediated enhancement of MMEJ in cells surviving radiation therapy may contribute to their radioresistance and could be therapeutically targeted.

59 BASIC BIOLOGICAL SCIENCES↗

Unique and conserved genome regions in Vibrio harveyi and related species in comparison with the shrimp pathogen Vibrio harveyi CAIM 1792

Vibrio harveyi CAIM 1792 is a marine bacterial strain that causes mortality in farmed shrimp in north-west Mexico, and the identification of virulence genes in this strain is important for understanding its pathogenicity. The goal of this study was to compare the V. harveyi CAIM 1792 genome with related genome sequences to determine their phylogenic relationship and explore unique regions in silico that differentiate this strain from other V. harveyi strains. Twenty-one newly sequenced genomes were compared in silico against the CAIM 1792 genome at nucleotidic and predicted proteome levels. The proteome of CAIM 1792 had higher similarity to those of other V. harveyi strains (78 %) than to those of the other closely related species Vibrio owensii (67 %), Vibrio rotiferianus (63 %) and Vibrio campbellii (59 %). Pan-genome ORFans trees showed the best fit with the accepted phylogeny based on DNA–DNA hybridization and multi-locus sequence analysis of 11 concatenated housekeeping genes. SNP analysis clustered 34/38 genomes within their accepted species. Laslty, the pangenomic and SNP trees showed that V. harveyi is the most conserved of the four species studied and V. campbellii may be divided into at least three subspecies, supported by intergenomic distance analysis. BLASTP atlases were created to identify unique regions among the genomes most related to V. harveyi CAIM 1792; these regions included genes encoding glycosyltransferases, specific type restriction modification systems and a transcriptional regulator, LysR, reported to be involved in virulence, metabolism, quorum sensing and motility.

59 BASIC BIOLOGICAL SCIENCES↗