Search NASASearch

SEARCH · Search NASA

Results for “Bioinformatics”

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 73 records · Page 4

rcsb-api : Python Toolkit for Streamlining Access to RCSB Protein Data Bank APIs

The Protein Data Bank (PDB) was founded in 1971 as the first open-access digital data resource in biology to serve as the single global archive for three-dimensional (3D) macromolecular structure data. Current PDB holdings exceed 230,000 experimentally determined structures of proteins, nucleic acids, viruses, and macromolecular machines. The RCSB Protein Data Bank RCSB.org research-focused web portal facilitates search, analyses, and visualization of every PDB structure along with more than one million Computed Structure Models from AlphaFold DB and the ModelArchive. It is powered by a set of publicly available Application Programming Interfaces (APIs) that both support RCSB.org users and provide programmatic access to PDB data. Given the breadth and levels of granularity encompassed in this rich data collection, efficiently accessing the information programmatically may be challenging for new users. RCSB PDB has developed a Python software package, rcsb-api , that facilitates easy and efficient use of RCSB PDB APIs within a Python environment. This software tool is designed to streamline access to the extensive corpus of data housed within the PDB, enabling researchers to search, retrieve, and analyze 3D biostructure data seamlessly. Its use will accelerate research in structural biology, molecular biology and biochemistry, drug discovery, and bioinformatics by providing more efficient tools for data integration and analysis. The new toolkit is available on GitHub (github.com/rcsb/py-rcsb-api) and published to the public Python package repository (PyPI) to foster wider usage and support basic and applied research in fundamental biology, biomedicine, and the energy sciences.

FAIR principles

Emerging protein sequencing technologies: proteomics without mass spectrometry?

Liquid chromatography-tandem mass spectrometry (LC-MS/MS) has been a leading method for proteomics for 30 years. Advantages provided by LC-MS/MS are offset by significant disadvantages, including cost. Recently, several non-mass spectrometric methods have emerged, but little information is available about their capacity to analyze the complex mixtures routine for mass spectrometry. Areas Covered: We review recent non-mass-spectrometric methods for sequencing proteins and peptides, including those using nanopores, sequencing by degradation, reverse translation, and short-epitope mapping, with comments on bioinformatics challenges, fundamental limitations, and areas where new technologies will be more or less competitive with LC-MS/MS. In addition to conventional literature searches, instrument vendor websites, patents, webinars, and preprints were also consulted to give a more up-to-date picture. Expert Opinion: Many new technologies are promising. However, demonstrations that they outperform mass spectrometry in terms of peptides and proteins identified have not yet been published, and astute observers note important disadvantages, especially relating to the dynamic range of single-molecule measurements of complex mixtures. Still, even if the performance of emerging methods proves inferior to LC-MS/MS, their low cost could create a different kind of revolution: a dramatic increase in the number of biology laboratories engaging in new forms of proteomics research.

59 BASIC BIOLOGICAL SCIENCES

Energy metric prediction for double insertion mutants via the RoseNet deep learning framework

Studying the structural and functional implications of protein mutations is an important task in computational biology and bioinformatics. We leverage our previously proposed RoseNet neural network architecture to predict energy metrics of proteins with double amino acid insertions or deletions (InDels). We train models on previously generated benchmark datasets containing the exhaustive double InDel mutations for three proteins, as well as an additional three proteins for which ∼145k random mutants, each with two InDels, have been generated. We expand on our previous work by evaluating three additional proteins and analyzing domain features that impact the prediction capabilities of RoseNet. These features include InDels into secondary structures and the solvent accessible surface area (SASA) scores of the residues. We uncover further evidence to support that RoseNet has a higher proficiency of generalizing to unseen residue combinations than unseen insertion positions. We also observe that RoseNet produces higher-quality predictions when inserting into a β-sheet over an α-helix. Additionally, when the insertions fall in an area of high SASA, RoseNet often displays better performance than inserting into areas of low SASA.

59 BASIC BIOLOGICAL SCIENCES

Mining Thermophile Photosynthesis Genes: A Synthetic Operon Expressing Chloroflexota Species Reaction Center Genes in Rhodobacter sphaeroides

Photosynthesis is the foundation of the vast majority of life systems, and is therefore the most important bioenergetic process on earth. The greatest diversity of photosynthetic systems is found in microorganisms. However, our understanding of the biophysical and biochemical processes that transduce light into chemical energy is derived from a relatively small subset of proteins from microbes that are amenable to cultivation, in contrast to the huge number of predicted proteins that catalyze the initial photochemical reactions deposited in databases, such as from metagenomics. We describe the use of a Rhodobacter sphaeroides laboratory strain for the expression of heterologous photosynthesis genes to demonstrate the feasibility of mining this resource, focusing on hot spring Chloroflexota gene sequences. Using a synthetic operon of genes, we produced a photochemically active complex of reaction center proteins in our biological system. We also present bioinformatic analyses of anoxygenic type II reaction center sequences from metagenomic samples collected from hot (42–90 °C) springs available through the JGI IMG database, to generate a resource of diverse sequences that are potentially adapted to photosynthesis at such temperatures. These data provide a view into the natural diversity of anoxygenic photosynthesis, through a lens focused on high-temperature environments. The approach we took to express such genes can be applied for potential biotechnology purposes as well as for studies of fundamental catalytic properties of these heretofore inaccessible protein complexes.

Chloroflexota

Helmet-Mounted Display Technology for Eva Training in NASA's Neutral Buoyancy Lab

Currently, during extravehicular activity (EVA) on the International Space Station (ISS), astronauts are in constant communication with mission control; however, this paradigm will significantly change during future NASA exploration missions due to limited and time-delayed communications. As exploration missions will thereby require an increase in crew autonomy, it is hypothesized that heads-up display (HUD) technology can supplement communications with mission control. Before introducing novel display technology to EVA operations, it must first be demonstrated and evaluated in an EVA training environment. The UC Davis Center for Spaceflight Research, in collaboration with NASA Johnson Space Center’s Human Physiology, Performance, Protection, and Operations (H-3PO) Lab has developed a low-cost and rapid development helmet-mounted display (HMD) for EVA training at NASA’s Neutral Buoyancy Lab (NBL). The first two phases of HMD testing (Phase I and II) at the NBL focused on technology and capability development and demonstration. Crew member feedback was obtained that evaluated aspects of HMD including its user satisfaction, visibility, readability, and usefulness. Overall, feedback indicated HMD was readable, un-obtrusive, and the available display modes had potential to be useful and enhance EVA training. HMD display modes currently available include: real-time biofeedback (metabolic rate), phase elapsed time (PET), ability to set timers, task procedural aids (e.g. Pistol Grip Tool settings) and custom messages. The goal for Phase III, which is currently underway, is to evaluate the effects that access to bioinformatics through HMD has on astronaut training and performance at the NBL. For this phase of testing, crewmembers will be asked to perform a go/no-go task while wearing HMD throughout the NBL run to characterize crew member state of situational awareness with respect to HMD. This information could be used to inform how best to display relatively important or supplemental, but not critical, data in the future.

Heads up display

Challenges and Opportunities in State‐of‐the‐Art Proteomics Analysis for Biomarker Development From Plasma Extracellular Vesicles

Extracellular vesicles (EVs) are membrane-bound particles secreted by cells, playing crucial roles in intercellular communication. The composition of EVs can undergo changes in response to stress and disease conditions, making them excellent biomarker candidates. However, extracting protein information from EVs can be challenging due to their low abundance in complex biofluids and copurification with contaminant proteins and particles. Techniques to enrich EVs have their strengths and limitations, without one being able to purify EVs to complete homogeneity. This can lead to compromised recovery rates and increased complexity, making data interpretation difficult. In this viewpoint article, we explore the concept that better characterization of EV composition, followed by quantification of EV proteins in complex samples, might be a more viable route for biomarker development. Mass spectrometers can provide reproducible deep coverage of the EV proteome, despite sample impurities. This paradigm shift presents opportunities to integrate advanced bioinformatics tools to refine the EV proteome landscape, identify novel biomarkers, and streamline validation processes in biomarker development. By focusing on leveraging technology rather than achieving absolute purity, this approach can transform current practices and open opportunities for robust biomarker discovery. Herein, we highlight not only such opportunities but also challenges to implement this concept.

Dakup, Panshak P. [Pacific Northwest National Labo

Visualizing and analyzing 3D biomolecular structures using Mol* at RCSB.org: Influenza A H5N1 virus proteome case study

The easiest and often most useful way to work with experimentally determined or computationally predicted structures of biomolecules is by viewing their three-dimensional (3D) shapes using a molecular visualization tool. Mol* was collaboratively developed by RCSB Protein Data Bank (RCSB PDB, RCSB.org) and Protein Data Bank in Europe (PDBe, PDBe.org) as an open-source, web-based, 3D visualization software suite for examination and analyses of biostructures. It is capable of displaying atomic coordinates and related experimental data of biomolecular structures together with a variety of annotations, facilitating basic and applied research, training, education, and information dissemination. Across RCSB.org, the RCSB PDB research-focused web portal, Mol* has been implemented to support single-mouse-click atomic-level visualization of biomolecules (e.g., proteins, nucleic acids, carbohydrates) with bound cofactors, small-molecule ligands, ions, water molecules, or other macromolecules. RCSB.org Mol* can seamlessly display 3D structures from various sources, allowing structure interrogation, superimposition, and comparison. Using influenza A H5N1 virus as a topical case study of an important pathogen, we exemplify how Mol* has been embedded within various RCSB.org tools—allowing users to view polymer sequence and structure-based annotations integrated from trusted bioinformatics data resources, assess patterns and trends in groups of structures, and view structures of any size and compositional complexity. In addition to being linked to every experimentally determined biostructure and Computed Structure Model made available at RCSB.org, Standalone Mol* is freely available for visualizing any atomic-level or multi-scale biostructure at rcsb.org/3d-view.

3D biostructure

Cohort-based pan-cancer analysis and experimental studies reveal ISG15 gene as a novel biomarker for prognosis and immunotherapy efficacy prediction

Abstract ISG15, an interferon-stimulated ubiquitin-like protein, plays a multifaceted role in tumorigenesis and immune regulation. This study comprehensively evaluates ISG15 as a prognostic biomarker and predictor of immunotherapy response through pan-cancer bioinformatics analysis and experimental validation. By integrating multiomics data from TCGA, GEO, and clinical cohorts, we found that ISG15 is significantly overexpressed in multiple cancers and generally correlates with poor prognosis. Elevated ISG15 expression is associated with increased immune checkpoint gene expression, particularly PD-L1, and immune infiltration, notably M2-like tumor-associated macrophages. Immunohistochemistry and multiplexed immunofluorescence confirmed a strong positive correlation between ISG15, PD-L1, and M2-TAM infiltration in lung and gastric cancer samples. Functional analysis at the single-cell level revealed significant associations between ISG15 and tumor proliferation, angiogenesis, and immune suppression. Immunotherapy cohort analysis demonstrated that tumors with high ISG15 expression responded favorably to PD-L1 inhibitors but exhibited resistance to CTLA-4 blockade, findings further validated in lung cancer patients receiving anti-PD-1 therapy. These results suggest that ISG15 is a promising biomarker for prognosis and immunotherapy response prediction across cancers. Its integration into clinical decision-making may enhance personalized treatment strategies, improve immunotherapy outcomes, and provide new insights into the tumor immune microenvironment, cancer progression, and potential therapeutic targets for future drug development.

Immunology

The unique architecture of umbrella toxins permits a two-tiered molecular bet-hedging strategy for interbacterial antagonism

Bacteria exist in competitive and rapidly changing environments in which the nature of future threats cannot be easily predicted. Streptomyces coelicolor produces three antibacterial umbrella particles that harbor distinct polymorphic toxin domains and an overlapping set of six diversified lectins. Here, we show that the exquisite specificity of umbrella particles derives from lectin-mediated species-specific binding to previously undescribed hypervariable surface glycoconjugates. A cryo-electron microscopy (cryo-EM) structure of one such lectin in complex with its oligosaccharide substrate defines the molecular basis for targeting through the coordinated recognition of multiple glycan features. Biochemical and genetic studies of several target species, in conjunction with lectin-swapping experiments, support a model whereby S. coelicolor umbrella toxin diversification at the levels of lectin composition and toxin polymorphism represents a unique, two-tiered bet-hedging strategy. Bioinformatic analyses support this as a means by which the unusual architecture of umbrella toxins offers Streptomyces a generalizable strategy to antagonize an unpredictable array of competitors.

59 BASIC BIOLOGICAL SCIENCES

Expression, purification, and characterization of diacylated Lipo-YcjN from Escherichia coli

YcjN is a putative substrate binding protein expressed from a cluster of genes involved in carbohydrate import and metabolism in Escherichia coli. Here, we determine the crystal structure of YcjN to a resolution of 1.95 Å, revealing that its three-dimensional structure is similar to substrate binding proteins in subcluster D-I, which includes the well-characterized maltose binding protein. Furthermore, we found that recombinant overexpression of YcjN results in the formation of a lipidated form of YcjN that is posttranslationally diacylated at cysteine 21. Comparisons of size-exclusion chromatography profiles and dynamic light scattering measurements of lipidated and nonlipidated YcjN proteins suggest that lipidated YcjN aggregates in solution via its lipid moiety. Additionally, bioinformatic analysis indicates that YcjN-like proteins may exist in both Bacteria and Archaea, potentially in both lipidated and nonlipidated forms. Together, our results provide a better understanding of the aggregation properties of recombinantly expressed bacterial lipoproteins in solution and establish a foundation for future studies that aim to elucidate the role of these proteins in bacterial physiology.

Escherichia coli

Functional diversification within the heme-binding split-barrel family

Due to neofunctionalization, a single fold can be identified in multiple proteins that have distinct molecular functions. Depending on the time that has passed since gene duplication and the number of mutations, the sequence similarity between functionally divergent proteins can be relatively high, eroding the value of sequence similarity as the sole tool for accurately annotating the function of uncharacterized homologs. Here, we combine bioinformatic approaches with targeted experimentation to reveal a large multifunctional family of putative enzymatic and nonenzymatic proteins involved in heme metabolism. This family (homolog of HugZ (HOZ)) is embedded in the “FMN-binding split barrel” superfamily and contains separate groups of proteins from prokaryotes, plants, and algae, which bind heme and either catalyze its degradation or function as nonenzymatic heme sensors. In prokaryotes these proteins are often involved in iron assimilation, whereas several plant and algal homologs are predicted to degrade heme in the plastid or regulate heme biosynthesis. In the plant Arabidopsis thaliana, which contains two HOZ subfamilies that can degrade heme in vitro (HOZ1 and HOZ2), disruption of AtHOZ1 (AT3G03890) or AtHOZ2A (AT1G51560) causes developmental delays, pointing to important biological roles in the plastid. In the tree Populus trichocarpa, a recent duplication event of a HOZ1 ancestor has resulted in localization of a paralog to the cytosol. Structural characterization of this cytosolic paralog and comparison to published homologous structures suggests conservation of heme-binding sites. This study unifies our understanding of the sequence-structure-function relationships within this multilineage family of heme-binding proteins and presents new molecular players in plant and bacterial heme metabolism.

59 BASIC BIOLOGICAL SCIENCES

Antibacterial ADP-ribosyl cyclase toxins inhibit bacterial growth by rapidly depleting NAD(P) +

In metazoans, enzymes belonging to the bifunctional ADP-ribosyl cyclase/cyclic ADP-ribose (cADPr) hydrolase family regulate diverse cellular processes by synthesizing and hydrolyzing the intracellular second messenger cADPr, derived from the electron carrier NAD+. However, bacterial enzymes belonging to this family have not been characterized. Here, we identify a bacterial ADP-ribosyl cyclase that is associated with the type VII secretion system and functions as an antibacterial toxin. This enzyme, which we name Tac1, inhibits bacterial growth by rapidly hydrolyzing NAD+ and NADP+. We determine the X-ray crystal structure of Tac1 to a resolution of 1.4 Å, which reveals that this protein adopts the core catalytic fold of metazoan ADP-ribosyl cyclase enzymes such as CD38. Using a combination of biochemical and mutagenesis approaches, we identify catalytic residues within the active site of Tac1, which are responsible for the formation of a cADPr catalytic intermediate and subsequent hydrolysis of this intermediate into linear ADP-ribose. A bioinformatic analysis reveals that Tac1 is the founding member of a widespread family of bacterial ADP-ribosyl cyclase enzymes, many of which are associated with interbacterial conflict systems. We also identify enzymes in this family that are not associated with biological conflict systems and demonstrate that they produce cADPr as their major product rather than linear ADP-ribose, suggesting that these enzymes serve a biological function distinct from interbacterial antagonism. Together, these findings demonstrate that ADP-ribosyl cyclase/cADPr hydrolase enzymes function as toxins in diverse bacterial conflict systems and suggest that cADPr may play a previously overlooked role in bacterial physiology.

Colautti, Jake

Environmental Contributions to Proton Sharing in Protein Low-Barrier Hydrogen Bonds

Hydrogen bonds (H-bonds) are central to biomolecular structure and dynamics. Although H-bonds are typically characterized by well-defined proton positions, proton delocalization has been proposed to play a role in facilitating enzyme catalysis and allostery in some systems. Experimentally locating protons is difficult, hampering the study of proton mobility in H-bonds. We used neutron crystallography, atomic resolution X-ray bond length analysis, and large quantum mechanics/molecular mechanics-Born–Oppenheimer molecular dynamics (QM/MM-BOMD) simulations to comprehensively characterize the shared proton/deuteron in a Glu–Asp low-barrier hydrogen bond (LBHB) in the bacterial protein YajL that is a conventional H-bond in the homologous disease-associated human protein DJ-1. X-ray bond length analysis of protiated and perdeuterated DJ-1 and YajL shows no significant effect of deuteron substitution on these carboxylic acid-carboxylate H-bonds but does reveal an effect at the active site glutamic acid near a cysteine thiolate. Residues in an H-bonded network that might favor LBHB formation in YajL were interrogated by the mutation of homologous residues in DJ-1. A distal DJ-1 substitution increases proton delocalization in the Glu–Asp H-bond, demonstrating that mutations within extended H-bond networks can modulate proton transfer barriers in carboxylic acid-carboxylate H-bonds. In addition, proton mobility in the H-bond is correlated with dimer-spanning motions in the QM/MM-BOMD simulations of YajL and DJ-1. Our results show that proton delocalization can be tuned using combined bioinformatic, structural, and computational information, opening the possibility of using engineered proton delocalization as a probe of H-bonding environments and as a tool to test hypotheses about LBHB function.

Lin, Jiusheng [University of Nebraska, Lincoln, NE

Bottom-Up Simulation, Reconstruction, and Quantification of Macromolecule Sequences from Experimental Polymerizations

Motivated by the canonical sequence–structure–function paradigm, tools to characterize chemical patterning in natural biomacromolecules, from proteins to nucleic acids, have grown exponentially in recent years. However, analogous strategies for synthetic macromolecules remain in nascent stages, complicated by sequence polydispersity and analytical limitations. To address this, we have developed a comprehensive and open-source Python package, PRISM (polymer rate insights and sequence modeling), an end-to-end workflow that provides a path from experimental kinetics measurements to quantitative and qualitative metrics for describing chemical patterning in stochastic polymers. First, a numerical integration strategy was constructed to simulate and fit experimental data from reversible addition–fragmentation chain transfer (RAFT) polymerization kinetics, enabling the facile estimation of relevant reactivity ratios. These ratios were then used in a mechanism-specific stochastic kinetic simulation strategy to simulate sequence ensembles corresponding to model systems spanning experimental copolymers, classes of statistical polymers (e.g., alternating, block, and gradient), and multiblock copolymers. Lastly, inspired by sequence homology metrics from bioinformatics, we introduce visualization strategies and quantitative metrics to facilitate comparisons of different sequence ensembles. As the sequence–structure–function paradigm becomes increasingly central in de novo design of synthetic macromolecules, this toolkit provides a first step toward accurate and representative sequence description and featurization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Quantitative Dissection of Agrobacterium Virulence to Generate a Synthetic Ti Plasmid

Agrobacterium is not only a costly plant pathogen but is also an essential tool for plant transformation. Though Agrobacterium-mediated transformation (AMT) has been heavily studied, its polygenic nature and complex transcriptional regulation make identification of the genetic basis of transformational efficiency difficult through traditional genetic and bioinformatic approaches. Here, we use a bottom-up synthetic approach to systematically engineer the tumor-inducing plasmid (pTi), wherein the majority of virulence machinery is encoded. Using a validated toolkit to control Agrobacterium gene expression in planta, we perform a quantitative dissection of AMT to investigate the contributions of critical vir-genes at different expression levels. We construct a synthetic pTi capable of transient plant and stable fungal transformation and characterize bottlenecks and solutions for complex polygenic synthetic pTi designs. Our reductionist approach demonstrates how bottom-up engineering can be used to dissect and elucidate the genetic underpinnings of complex biological traits, laying the foundation for future engineering to establish full synthetic control over the critical process of AMT.

Agrobacterium-mediated transformation

Structural Insights into the Mechanism of a Polyketide Synthase Thiocysteine Lyase Domain

Polyketide synthases (PKSs) are renowned for the structural diversity of the polyketide natural products they produce, but sulfur-containing functionalities are rarely installed by PKSs. We previously characterized thiocysteine lyase (SH) domains involved in the biosynthesis of the leinamycin (LNM) family of natural products, exemplified by LnmJ-SH and guangnanmycin (GnmT-SH). Here we report a detailed investigation into the PLP-dependent reaction catalyzed by the SH domains, guided by a 1.8 Å resolution crystal structure of GnmT-SH. A series of elaborate substrate mimics were synthesized to answer specific questions garnered from the crystal structure and from the biosynthetic logic of the LNM family of natural products. Here, through a combination of bioinformatics, molecular modeling, in vitro assays, and mutagenesis, we have developed a detailed model of acyl carrier protein (ACP)-tethered substrate-SH, and interdomain interactions, that contribute to the observed substrate specificity. Comparison of the GnmT-SH structure with archetypical PLP-dependent enzyme structures revealed how Nature, via evolution, has modified a common protein structural motif to accommodate an ACP-tethered substrate, which is significantly larger than any of those previously characterized. Overall, this study demonstrates how PLP-dependent chemistry can be incorporated into the context of PKS assembly lines and sets the stage for engineering PKSs to produce sulfur-containing polyketides.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Systematic benchmarking demonstrates large language models have not reached the diagnostic accuracy of traditional rare-disease decision support tools

Large language models (LLMs) show promise in supporting differential diagnosis, but their performance is challenging to evaluate due to the unstructured nature of their responses, and their accuracy compared to existing diagnostic tools is not well characterized. To assess the current capabilities of LLMs to diagnose genetic diseases, we benchmarked these models on 5213 previously published case reports using the Phenopacket Schema, the Human Phenotype Ontology and Mondo disease ontology. Prompts generated from each phenopacket were sent to seven LLMs, including four generalist models and three LLMs specialized for medical applications. The same phenopackets were used as input to a widely used diagnostic tool, Exomiser, in phenotype-only mode. The best LLM ranked the correct diagnosis first in 23.6% of cases, whereas Exomiser did so in 35.5% of cases. While the performance of LLMs for supporting differential diagnosis has been improving, it has not reached the level of commonly used traditional bioinformatics tools. Future research is needed to determine the best approach to incorporate LLMs into diagnostic pipelines.

Reese, Justin T. [Lawrence Berkeley National Labor

Microbial species and intraspecies units exist and are maintained by ecological cohesiveness coupled to high homologous recombination

Abstract Recent genomic analyses have revealed that microbial communities are predominantly composed of persistent, sequence-discrete species and intraspecies units (genomovars), but the mechanisms that create and maintain these units remain unclear. By analyzing closely-related isolate genomes from the same or related samples and identifying recent recombination events using a novel bioinformatics methodology, we show that high ecological cohesiveness coupled to frequent-enough and unbiased (i.e., not selection-driven) horizontal gene flow, mediated by homologous recombination, often underlie these diversity patterns. Ecological cohesiveness was inferred based on greater similarity in temporal abundance patterns of genomes of the same vs. different units, and recombination was shown to affect all sizable segments of the genome (i.e., be genome-wide) and have two times or greater impact on sequence evolution than point mutations. These results were observed in bothSalinibacter ruber, an environmental halophilic organism, andEscherichia coli, the model gut-associated organism and an opportunistic pathogen, indicating that they may be more broadly applicable to the microbial world. Therefore, our results represent a departure compared to previous models of microbial speciation that invoke either ecology or recombination, but not necessarily their synergistic effect, and answer an important question for microbiology: what a species and a subspecies are.

Science & Technology - Other Topics