Search NASA⌕ Search

SEARCH · Search NASA

Results for “Sequence Alignment”

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 217 records · Page 12

Vascular defense responses in rice: peroxidase accumulation in xylem parenchyma cells and xylem wall thickening

The rice bacterial blight pathogen Xanthomonas oryzae pv. oryzae is a vascular pathogen that elicits a defensive response through interaction with metabolically active rice cells. In leaves of 12-day-old rice seedlings, the exposed pit membrane separating the xylem lumen from the associated parenchyma cells allows contact with bacterial cells. During resistant responses, the xylem secondary walls thicken within 48 h and the pit diameter decreases, effectively reducing the area of pit membrane exposed for access by bacteria. In susceptible interactions and mock-inoculated controls, the xylem walls do not thicken within 48 h. Xylem secondary wall thickening is developmental and, in untreated 65-day-old rice plants, the size of the pit also is reduced. Activity and accumulation of a secreted cationic peroxidase, PO-C1, were previously shown to increase in xylem vessel walls and lumen. Peptide-specific antibodies and immunogold-labeling were used to demonstrate that PO-C1 is produced in the xylem parenchyma and secreted to the xylem lumen and walls. The timing of the accumulation is consistent with vessel secondary wall thickening. The PO-C1 gene is distinct but shares a high level of similarity with previously cloned pathogen-induced peroxidases in rice. PO-C1 gene expression was induced as early as 12 h during resistant interactions and peaked between 18 and 24 h after inoculation. Expression during susceptible interactions was lower than that observed in resistant interactions and was undetectable after infiltration with water, after mechanical wounding, or in mature leaves. These data are consistent with a role for vessel secondary wall thickening and peroxidase PO-C1 accumulation in the defense response in rice to X. oryzae pv. oryzae.

NASA Discipline Plant Biology↗

Multiplex detection and identification of viral, bacterial, and protozoan pathogens in human blood and plasma using an expanded high-density resequencing microarray platform

Introduction: Nucleic acid tests for blood donor screening have improved the safety of the blood supply; however, increasing numbers of emerging pathogen tests are burdensome. Multiplex testing platforms are a potential solution. Methods: The Blood Borne Pathogen Resequencing Microarray Expanded (BBP-RMAv.2) can perform multiplex detection and identification of 80 viruses, bacteria and parasites. This study evaluated pathogen detection in human blood or plasma. Samples spiked with selected pathogens, each with one of 6 viruses, 2 bacteria and 5 protozoans were tested on this platform. The nucleic acids were extracted, amplified using multiplexed sets of primers, and hybridized to a microarray. The reported sequences were aligned to a database to identify the pathogen. To directly compare the microarray to an emerging molecular approach, the amplified nucleic acids were also submitted to nanopore next generation sequencing (NGS). Results: The BBP-RMAv.2 detected viral pathogens at a concentration as low as 100 copies/ml and a range of concentrations from 1,000 to 100,000 copies/ml for all the spiked pathogens. Coded specimens were identified correctly demonstrating the effectiveness of the platform. The nanopore sequencing correctly identified most samples and the results of the two platforms were compared. Discussion: These results indicated that the BBP-RMAv.2 could be employed for multiplex detection with potential for use in blood safety or disease diagnosis. The NGS was nearly as effective at identifying pathogens in blood and performed better than BBP-RMAv.2 at identifying pathogen-negative samples.

59 BASIC BIOLOGICAL SCIENCES↗

The prokaryote-to-eukaryote transition reflected in the evolution of the V/F/A-ATPase catalytic and proteolipid subunits

Changes in the primary and quarternary structure of vacuolar and archaeal type ATPases that accompany the prokaryote-to-eukaryote transition are analyzed. The gene encoding the vacuolar-type proteolipid of the V-ATPase from Giardia lamblia is reported. Giardia has a typical vacuolar ATPase as observed from the common motifs shared between its proteolipid subunit and other eukaryotic vacuolar ATPases, suggesting that the former enzyme works as a hydrolase in this primitive eukaryote. The phylogenetic analyses of the V-ATPase catalytic subunit and the front and back halves of the proteolipid subunit placed Giardia as the deepest branch within the eukaryotes. Our phylogenetic analysis indicated that at least two independent duplication and fusion events gave rise to the larger proteolipid type found in eukaryotes and in Methanococcus. The spatial distribution of the conserved residues among the vacuolar-type proteolipids suggest a zipper-type interaction among the transmembrane helices and surrounding subunits of the V-ATPase complex. Important residues involved in the function of the F-ATP synthase proteolipid have been replaced during evolution in the V-proteolipid, but in some cases retained in the archaeal A-ATPase. Their possible implication in the evolution of V/F/A-ATPases is discussed.

NASA Discipline Exobiology↗

Structure, function, and evolution of bacterial ATP-binding cassette systems

The ATP-binding cassette (ABC) systems constitute one of the largest superfamilies of paralogous sequences. All ABC systems share a highly conserved ATP-hydrolyzing domain or protein (the ABC; also referred to as a nucleotide-binding domain [NBD]) that is unequivocally characterized by three short sequence motifs (Fig. 1): these are the Walker A and Walker B motifs, indicative of the presence of a nucleotide-binding site, and the signature motif, unique to ABC proteins, located upstream of the Walker B motif (426). Other motifs diagnostic of ABC proteins are also indicated in Fig. 1. The biological significance of these motifs is discussed in Structure, Function, and Dynamics of the ABC. ABC systems are widespread among living organisms and have been detected in all genera of the three kingdoms of life, with remarkable conservation in the primary sequence of the cassette and in the organization of the constitutive domains or subunits (203, 420). ABC systems couple the energy of ATP hydrolysis to an impressively large variety of essential biological phenomena, comprising not only transmembrane (TM) transport, for which they are best known, but also several non-transport-related processes, such as translation elongation (62) and DNA repair (174). Although ABC systems deserve much attention because they are involved in severe human inherited diseases (107), they were first discovered and characterized in detail in prokaryotes, as early as the 1970s (13, 148, 238, 468). The most extensively analyzed systems were the high-affinity histidine and maltose uptake systems of Salmonella enterica serovar Typhimurium and Escherichia coli. Over 2 decades ago, after the completion of the nucleotide sequences encoding these transporters in the respective laboratories of Giovanna Ames and Maurice Hofnung, Hiroshi Nikaido and colleagues noticed that the two systems displayed a global similarity in the nature of their components and, moreover, that the primary sequences of MalK and HisP, the proteins suspected to energize these transporters, shared as much as 32% identity in amino acid residues when their sequences were aligned (171). Later, it was found that several bacterial proteins involved in uptake of nutrients, export of toxins, cell division, bacterial nodulation of plants, and DNA repair displayed the same similarity in their sequences (127, 196). This led to the notion that the conserved protein, which had been shown to bind ATP (198, 201), would probably energize the systems mentioned above by coupling the energy of ATP hydrolysis to transport. The latter was demonstrated with the maltose and histidine transporters by use of isolated membrane vesicles (105, 379) and purified transporters reconstituted into proteoliposomes (30, 98). The determination of the sequence of the first eukaryotic protein strongly similar to these bacterial transporters (the P-glycoprotein, involved in resistance of cancer cells to multiple drugs) (169, 179) demonstrated that these proteins were not restricted to prokaryotes. Two names, 'traffic ATPases' (15) and the more accepted name 'ABC transporters' (193, 218), were proposed for members of this new superfamily. ABC systems can be divided into three main functional categories, as follows. Importers mediate the uptake of nutrients in prokaryotes. The nature of the substrates that are transported is very wide, including mono- and oligosaccharides, organic and inorganic ions, amino acids, peptides, ironsiderophores, metals, polyamine cations, opines, and vitamins. Exporters are involved in the secretion of various molecules, such as peptides, lipids, hydrophobic drugs, polysaccharides, and proteins, including toxins such as hemolysin. The third category of systems is apparently not involved in transport, with some members being involved in translation of mRNA and in DNA repair. Despite the large, diverse population of substrates handled and the difference in the polarity of transport, importers and exporters share a common organization made of two hydrophobic membrane-spanning or integral membrane (IM) domains and two hydrophilic domains carrying the ABC peripherally associated with the IM domains on the cytosolic side of the membrane (26). In importers, these four domains are almost always independent polypeptide chains that come together to form a multimeric complex. In most exporters, including the E. coli hemolysin exporter HlyB, the N-terminal IM and the C-terminal ABC domains are fused as a single polypeptide chain (IM-ABC). An inverted organization in which the IM domain is C-terminal with respect to the ABC domain (ABC-IM) exists, such as in the MacB protein, involved in macrolide resistance in E. coli. No IM domain partners have been identified for ABC proteins falling into the third category, and these proteins consist of two ABCs fused together (ABC2).

59 BASIC BIOLOGICAL SCIENCES↗

Scalable Parallel Methods for Analyzing Metagenomics Data at Extreme Scale

The field of bioinformatics and computational biology is currently experiencing a data revolution. The exciting prospect of making fundamental biological discoveries is fueling the rapid development and deployment of numerous cost-effective, high-throughput next-generation sequencing technologies. The result is that the DNA and protein sequence repositories are being bombarded with new sequence information. Databases are continuing to report a Moore’s law-like growth trajectory in their database sizes, roughly doubling every 18 months. In what seems to be a paradigm-shift, individual projects are now capable of generating billions of raw sequence data that need to be analyzed in the presence of already annotated sequence information. While it is clear that data-driven methods, such as sequencing homology detection, are becoming the mainstay in the field of computational life sciences, the algorithmic advancements essential for implementing complex data analytics at scale have mostly lagged behind. Sequence homology detection is central to a number of bioinformatics applications including genome sequencing and protein family characterization. Given millions of sequences, the goal is to identify all pairs of sequences that are highly similar (or “homologous”) on the basis of alignment criteria. While there are optimal alignment algorithms to compute pairwise homology, their deployment for large-scale is currently not feasible; instead, heuristic methods are used at the expense of quality. In this dissertation, we present the design and evaluation of a parallel implementation for conducting optimal homology detection on distributed memory supercomputers. Our approach uses a combination of techniques from asynchronous load balancing (viz. work stealing, dynamic task counters), data replication, and exact-matching filters to achieve homology detection at scale. Results for a collection of 2.56M sequences show parallel efficiencies of ~75-100% on up to 8K cores, representing a time-to-solution of 33 seconds. We extend this work with a detailed analysis of single-node sequence alignment performance using the latest CPU vector instruction set extensions. Preliminary results reveal that current sequence alignment algorithms are unable to fully utilize widening vector registers.

59 BASIC BIOLOGICAL SCIENCES↗

The expression and function of the achaete-scute genes in Tribolium castaneum reveals conservation and variation in neural pattern formation and cell fate specification

The study of achaete-scute (ac/sc) genes has recently become a paradigm to understand the evolution and development of the arthropod nervous system. We describe the identification and characterization of the ac/sc genes in the coleopteran insect species Tribolium castaneum. We have identified two Tribolium ac/sc genes - achaete-scute homolog (Tc-ASH) a proneural gene and asense (Tc-ase) a neural precursor gene that reside in a gene complex. Focusing on the embryonic central nervous system we find that Tc-ASH is expressed in all neural precursors and the proneural clusters from which they segregate. Through RNAi and misexpression studies we show that Tc-ASH is necessary for neural precursor formation in Tribolium and sufficient for neural precursor formation in Drosophila. Comparison of the function of the Drosophila and Tribolium proneural ac/sc genes suggests that in the Drosophila lineage these genes have maintained their ancestral function in neural precursor formation and have acquired a new role in the fate specification of individual neural precursors. Furthermore, we find that Tc-ase is expressed in all neural precursors suggesting an important and conserved role for asense genes in insect nervous system development. Our analysis of the Tribolium ac/sc genes indicates significant plasticity in gene number, expression and function, and implicates these modifications in the evolution of arthropod neural development.

Non-NASA Center↗

Expression of the ctenophore Brain Factor 1 forkhead gene ortholog (ctenoBF-1) mRNA is restricted to the presumptive mouth and feeding apparatus: implications for axial organization in the Metazoa

Ctenophores are thoroughly modern animals whose ancestors are derived from a separate evolutionary branch than that of other eumetazoans. Their major longitudinal body axis is the oral-aboral axis. An apical sense organ, called the apical organ, is located at the aboral pole and contains a highly innervated statocyst and photodetecting cells. The apical organ integrates sensory information and controls the locomotory apparatus of ctenophores, the eight longitudinal rows of ctene/comb plates. In an effort to understand the developmental and evolutionary organization of axial properties of ctenophores we have isolated a forkhead gene from the Brain Factor 1 (BF-1) family. This gene, ctenoBF-1, is the first full-length nuclear gene reported from ctenophores. This makes ctenophores the most basal metazoan (to date) known to express definitive forkhead class transcription factors. Orthologs of BF-1 in vertebrates, Drosophila, and Caenorhabditis elegans are expressed in anterior neural structures. Surprisingly, in situ hybridizations with ctenoBF-1 antisense riboprobes show that this gene is not expressed in the apical organ of ctenophores. CtenoBF-1 is expressed prior to first cleavage. Transcripts become localized to the aboral pole by the 8-cell stage and are inherited by ectodermal micromeres generated from this region at the 16- and 32-cell stages. Expression in subsets of these cells persists and is seen around the edge of the blastopore (presumptive mouth) and in distinct ectodermal regions along the tentacular poles. Following gastrulation, stomodeal expression begins to fade and intense staining becomes restricted to two distinct domains in each tentacular feeding apparatus. We suggest that the apical organ is not homologous to the brain of bilaterians but that the oral pole of ctenophores corresponds to the anterior pole of bilaterian animals.

Non-NASA Center↗

Classification of bacterial plasmid and chromosome derived sequences using machine learning

Plasmids are important genetic elements that facilitate horizonal gene transfer between bacteria and contribute to the spread of virulence and antimicrobial resistance. Most bacterial genome sequences in the public archives exist in draft form with many contigs, making it difficult to determine if a contig is of chromosomal or plasmid origin. Using a training set of contigs comprising 10,584 chromosomes and 10,654 plasmids from the PATRIC database, we evaluated several machine learning models including random forest, logistic regression, XGBoost, and a neural network for their ability to classify chromosomal and plasmid sequences using nucleotide k-mers as features. Based on the methods tested, a neural network model that used nucleotide 6-mers as features that was trained on randomly selected chromosomal and plasmid subsequences 5kb in length achieved the best performance, outperforming existing out-of-the-box methods, with an average accuracy of 89.38% ± 2.16% over a 10-fold cross validation. The model accuracy can be improved to 92.08% by using a voting strategy when classifying holdout sequences. In both plasmids and chromosomes, subsequences encoding functions involved in horizontal gene transfer—including hypothetical proteins, transporters, phage, mobile elements, and CRISPR elements—were most likely to be misclassified by the model. This study provides a straightforward approach for identifying plasmid-encoding sequences in short read assemblies without the need for sequence alignment-based tools.

59 BASIC BIOLOGICAL SCIENCES↗

SCOPe: improvements to the structural classification of proteins – extended database to facilitate variant interpretation and machine learning

Abstract The Structural Classification of Proteins—extended (SCOPe, https://scop.berkeley.edu) knowledgebase aims to provide an accurate, detailed, and comprehensive description of the structural and evolutionary relationships amongst the majority of proteins of known structure, along with resources for analyzing the protein structures and their sequences. Structures from the PDB are divided into domains and classified using a combination of manual curation and highly precise automated methods. In the current release of SCOPe, 2.08, we have developed search and display tools for analysis of genetic variants we mapped to structures classified in SCOPe. In order to improve the utility of SCOPe to automated methods such as deep learning classifiers that rely on multiple alignment of sequences of homologous proteins, we have introduced new machine-parseable annotations that indicate aberrant structures as well as domains that are distinguished by a smaller repeat unit. We also classified structures from 74 of the largest Pfam families not previously classified in SCOPe, and we improved our algorithm to remove N- and C-terminal cloning, expression and purification sequences from SCOPe domains. SCOPe 2.08-stable classifies 106 976 PDB entries (about 60% of PDB entries).

59 BASIC BIOLOGICAL SCIENCES↗

Data from: Generation of transcript assemblies and identification of single nucleotide polymorphisms from seven lowland and upland cultivars of switchgrass

Switchgrass (Panicum virgatum) is a North American perennial prairie species that has been used as a rangeland and forage crop and has recently been targeted as a potential biofuel feedstock species. Switchgrass, which occurs as tetraploid and octoploid forms, is classified into lowland or upland ecotypes that differ in growth phenotypes and adaptation to distinct habitats. Using RNA-sequencing reads derived from crown, young shoot and leaf tissues, we generated sequence data from seven switchgrass cultivars, three lowland and four upland, to enable comparative analyses between switchgrass cultivars and to identify single nucleotide polymorphisms (SNPs) for use in breeding and genetic analysis. We also generated individual transcript assemblies for each of the cultivars. Transcript data indicate that subgenomes of octoploid switchgrass are not substantially different from subgenomes of tetraploids as expected for an autopolyploid origin of switchgrass octoploids. Using RNA-sequencing reads aligned to the switchgrass Release 0 AP13 reference genome, we identified 1,305,976 high-confidence SNPs. Of these SNPs, 438,464 were unique to lowland cultivars, but only 12,002 were found in all lowlands. Conversely, 723,678 SNPs were unique to upland cultivars with only 34,665 observed in all uplands. Comparison of our high-confidence transcriptome-derived SNPs with SNPs previously identified in a genotyping-by-sequencing study of an association panel revealed limited overlap between the two methods, highlighting the utility of transcriptome-based SNP discovery in augmenting genome diversity polymorphism datasets. The transcript and SNP data described here provide a useful resource for switchgrass gene annotation and marker-based analyses of the switchgrass genome.

09 BIOMASS FUELS↗

Functional characteristics of the calcium modulated proteins seen from an evolutionary perspective

We have constructed dendrograms relating 173 EF-hand proteins of known amino acid sequence. We aligned all of these proteins by their EF-hand domains, omitting interdomain regions. Initial dendrograms were computed by minimum mutation distance methods. Using these as starting points, we determined the best dendrogram by the method of maximum parsimony, scored by minimum mutation distance. We identified 14 distinct subfamilies as well as 6 unique proteins that are perhaps the sole representatives of other subfamilies. This information is given in tabular form. Within subfamilies one can easily align interdomain regions. The resulting dendrograms are very similar to those computed using domains only. Dendrograms constructed using pairs of domains show general congruence. However, there are enough exceptions to caution against an overly simple scheme in which one pair of gene duplications leads from one domain precurser to a four domain prototype from which all other forms evolved. The ability to bind calcium was lost and acquired several times during evolution. The distribution of introns does not conform to the dendrogram based on amino acid sequences. The rates of evolution appear to be much slower within subfamilies, especially within calmodulin, than those prior to the definition of subfamily.

Kretsinger, R. H.↗

Analysis of genomic signatures associated with Variovorax endosphere colonization

This repository contains the analysis code and supporting datasets associated with the study “Genomic signatures in Variovorax enabling colonization of the Populus endosphere.” Beals DG, Carper DL, Hochanadel LH, Jawdy SS, Klingeman DM, Piatkowski BT, Weston DJ, Doktycz MJ, Pelletier DA. 2026. Genomic signatures in Variovorax enabling colonization of the Populus endosphere. mSystems 11:e01605-25. https://doi.org/10.1128/msystems.01605-25 The scripts are organized sequentially (01–07) and document the workflows used for: Sequence-read alignment and feature counting Orthogroup and KEGG Ortholog annotation Count normalization Statistical analysis and aggregation Generation of manuscript figures and tables Repository contents The uncompressed files are the finalized, formatted datasets used to generate the figures and tables reported in the study, including the supplemental CSV files referenced in the manuscript. The accompanying ZIP archive contains the complete codebase and example data_input/ and data_output/ directories illustrating the organization and execution of the analytical workflow. Individual scripts identify the corresponding manuscript analyses and figure panels. Raw sequencing data Raw sequencing reads are available through the NCBI Sequence Read Archive under BioProject accession PRJNA1322484.

Beals, Delaney [ORNL] (ORCID:0000000306274574)↗

Structures of TOG1 and TOG2 from the human microtubule dynamics regulator CLASP1

Tubulin-binding TOG domains are found arrayed in a number of proteins that regulate microtubule dynamics. While much is known about the structure and function of TOG domains from the XMAP215 microtubule polymerase family, less in known about the TOG domain array found in animal CLASP family members. The animal CLASP TOG array promotes microtubule pause, potentiates rescue, and limits catastrophe. How structurally distinct the TOG domains of animal CLASP are from one another, from XMAP215 family TOG domains, and whether a specific order of structurally distinct TOG domains in the TOG array is conserved across animal CLASP family members is poorly understood. We present the x-ray crystal structures of Homo sapiens (H.s.) CLASP1 TOG1 and TOG2. The structures of H.s. CLASP1 TOG1 and TOG2 are distinct from each other and from the previously determined structure of Mus musculus (M.m.) CLASP2 TOG3. Comparative analyses of CLASP family TOG domain structures determined to date across species and paralogs supports a conserved CLASP TOG array paradigm in which structurally distinct TOG domains are arrayed in a specific order. H.s. CLASP1 TOG1 bears structural similarity to the free-tubulin binding TOG domains of the XMAP215 family but lacks many of the key tubulin-binding determinants found in XMAP215 family TOG domains. This aligns with studies that report that animal CLASP family TOG1 domains cannot bind free tubulin or microtubules. In contrast, animal CLASP family TOG2 and TOG3 domains have reported microtubule-binding activity but are structurally distinct from the free-tubulin binding TOG domains of the XMAP215 family. H.s. CLASP1 TOG2 has a convex architecture, predicted to engage a hyper-curved tubulin state that may underlie its ability to limit microtubule catastrophe and promote rescue. M.m. CLASP2 TOG3 has unique structural elements in the C-terminal half of its α-solenoid domain that our modeling studies implicate in binding to laterally-associated tubulin subunits in the microtubule lattice in a mode similar to, yet distinct from those predicted for the XMAP215 family TOG4 domain. The potential ability of the animal CLASP family TOG3 domain to engage lateral tubulin subunits may underlie the microtubule rescue activity ascribed to the domain. These findings highlight the structural diversity of TOG domains within the CLASP family TOG array and provide a molecular foundation for understanding CLASP-dependent effects on microtubule dynamics.

59 BASIC BIOLOGICAL SCIENCES↗

Standardized Residue Numbering and Secondary Structure Nomenclature in the Class D β-Lactamases

Over 1370 class D β-lactamases are currently known, and they pose a serious threat to the effective treatment of many infectious diseases, particularly in some pathogenic bacteria where evolving carbapenemase activity has been reported. Detailed understanding of their molecular biology, enzymology, and structural biology are critically important, but the lack of a standardized residue numbering scheme and inconsistent secondary structure annotation has made comparative analyses sometimes difficult and cumbersome. Compounding this, in the post-AlphaFold world where we currently find ourselves, an extraordinary wealth of detailed structural information on these enzymes is literally at our fingertips; therefore it is vitally important that a standard numbering system is in place to facilitate the accurate and straightforward analysis of their structures. In conclusion, here we present a residue numbering and secondary structure scheme for the class D enzymes based on the sequence and structure of OXA-48 and apply it to test targets to demonstrate the ease with which it can be used.

59 BASIC BIOLOGICAL SCIENCES↗

Proteome-wide analysis of HIV-specific naive and memory CD4 + T cells in unexposed blood donors

The preexisting HIV-1–specific T cell repertoire must influence both the immunodominance of T cells after infection and immunogenicity of vaccines. We directly compared two methods for measuring the preexisting CD4 + T cell repertoire in healthy HIV-1–negative volunteers, the HLA-peptide tetramer enrichment and T cell library technique, and show high concordance (r = 0.989). Using the library technique, we examined whether naive, central memory, and/or effector memory CD4 + T cells specific for overlapping peptides spanning the entire HIV-1 proteome were detectable in 10 HLA diverse, HIV-1–unexposed, seronegative donors. HIV-1–specific cells were detected in all donors at a mean of 55 cells/million naive cells and 38.9 and 34.1 cells/million in central and effector memory subsets. Remarkably, peptide mapping showed most epitopes recognized by naive (88%) and memory (56%) CD4 + T cells had been previously reported in natural HIV-1 infection. Furthermore, 83% of epitopes identified in preexisting memory subsets shared epitope length matches (8–12 amino acids) with human microbiome proteins, suggestive of a possible cross-reactive mechanism. These results underline the power of a proteome-wide analysis of peptide recognition by human T cells for the identification of dominant antigens and provide a baseline for optimizing HIV-1–specific helper cell responses by vaccination.

59 BASIC BIOLOGICAL SCIENCES↗

The repABC plasmids with quorum-regulated transfer systems in members of the Rhizobiales divide into two structurally and separately evolving groups

The large repABC plasmids of the order Rhizobiales with Class I quorum-regulated conjugative transfer systems often define the nature of the bacterium that harbors them. These otherwise diverse plasmids contain a core of highly conserved genes for replication and conjugation raising the question of their evolutionary relationships. In an analysis of 18 such plasmids these elements fall into two organizational classes, Group I and Group II, based on the sites at which cargo DNA is located. Cladograms constructed from proteins of the transfer and quorum-sensing components indicated that those of the Group I plasmids, while coevolving, have diverged from those coevolving proteins of the Group II plasmids. Moreover, within these groups the phylogenies of the proteins usually occupy similar, if not identical, tree topologies. Remarkably, such relationships were not seen among proteins of the replication system; although RepA and RepB coevolve, RepC does not. Nor do the replication proteins coevolve with the proteins of the transfer and quorum-sensing systems. Functional analysis was mostly consistent with phylogenies. TraR activated promoters from plasmids within its group, but not between groups and dimerized with TraR proteins from within but not between groups. However, oriT sequences, which are highly conserved, were processed by the transfer system of plasmids regardless of group. Here, we conclude that these plasmids diverged into two classes based on the locations at which cargo DNA is inserted, that the quorum-sensing and transfer functions are coevolving within but not between the two groups, and that this divergent evolution extends to function.

59 BASIC BIOLOGICAL SCIENCES↗

Tandem repeats in giant archaeal Borg elements undergo rapid evolution and create new intrinsically disordered regions in proteins

Borgs are huge, linear extrachromosomal elements associated with anaerobic methane-oxidizing archaea. Striking features of Borg genomes are pervasive tandem direct repeat (TR) regions. Here, we present six new Borg genomes and investigate the characteristics of TRs in all ten complete Borg genomes. We find that TR regions are rapidly evolving, recently formed, arise independently, and are virtually absent in host Methanoperedens genomes. Flanking partial repeats and A-enriched character constrain the TR formation mechanism. TRs can be in intergenic regions, where they might serve as regulatory RNAs, or in open reading frames (ORFs). TRs in ORFs are under very strong selective pressure, leading to perfect amino acid TRs (aaTRs) that are commonly intrinsically disordered regions. Proteins with aaTRs are often extracellular or membrane proteins, and functionally similar or homologous proteins often have aaTRs composed of the same amino acids. We propose that Borg aaTR-proteins functionally diversify Methanoperedens and all TRs are crucial for specific Borg–host associations and possibly cospeciation.

59 BASIC BIOLOGICAL SCIENCES↗

A taxonomy of bacterial microcompartment loci constructed by a novel scoring method

Bacterial microcompartments (BMCs) are proteinaceous organelles involved in both autotrophic and heterotrophic metabolism. All BMCs share homologous shell proteins but differ in their complement of enzymes; these are typically encoded adjacent to shell protein genes in genetic loci, or operons. To enable the identification and prediction of functional (sub)types of BMCs, we developed LoClass, an algorithm that finds putative BMC loci and inventories, weights, and compares their constituent pfam domains to construct a locus similarity network and predict locus (sub)types. In addition to using LoClass to analyze sequences in the Non-redundant Protein Database, we compared predicted BMC loci found in seven candidate bacterial phyla (six from single-cell genomic studies) to the LoClass taxonomy. Together, these analyses resulted in the identification of 23 different types of BMCs encoded in 30 distinct locus (sub)types found in 23 bacterial phyla. These include the two carboxysome types and a divergent set of metabolosomes, BMCs that share a common catalytic core and process distinct substrates via specific signature enzymes. Furthermore, many Candidate BMCs were found that lack one or more core metabolosome components, including one that is predicted to represent an entirely new paradigm for BMC-associated metabolism, joining the carboxysome and metabolosome. By placing these results in a phylogenetic context, we provide a framework for understanding the horizontal transfer of these loci, a starting point for studies aimed at understanding the evolution of BMCs. This comprehensive taxonomy of BMC loci, based on their constituent protein domains, foregrounds the functional diversity of BMCs and provides a reference for interpreting the role of BMC gene clusters encoded in isolate, single cell, and metagenomic data. Many loci encode ancillary functions such as transporters or genes for cofactor assembly; this expanded vocabulary of BMC-related functions should be useful for design of genetic modules for introducing BMCs in bioengineering applications.

59 BASIC BIOLOGICAL SCIENCES↗