Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “protein function prediction”

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 487 records · Page 27

Hierarchical, rotation‐equivariant neural networks to select structural models of protein complexes

Abstract Predicting the structure of multi‐protein complexes is a grand challenge in biochemistry, with major implications for basic science and drug discovery. Computational structure prediction methods generally leverage predefined structural features to distinguish accurate structural models from less accurate ones. This raises the question of whether it is possible to learn characteristics of accurate models directly from atomic coordinates of protein complexes, with no prior assumptions. Here we introduce a machine learning method that learns directly from the 3D positions of all atoms to identify accurate models of protein complexes, without using any precomputed physics‐inspired or statistical terms. Our neural network architecture combines multiple ingredients that together enable end‐to‐end learning from molecular structures containing tens of thousands of atoms: a point‐based representation of atoms, equivariance with respect to rotation and translation, local convolutions, and hierarchical subsampling operations. When used in combination with previously developed scoring functions, our network substantially improves the identification of accurate structural models among a large set of possible models. Our network can also be used to predict the accuracy of a given structural model in absolute terms. The architecture we present is readily applicable to other tasks involving learning on 3D structures of large atomic systems.

Eismann, Stephan↗

A minimal complex of KHNYN and zinc-finger antiviral protein binds and degrades single-stranded RNA

Detecting viral infection is a key role of the innate immune system. The genomes of some RNA viruses have a high CpG dinucleotide content relative to most vertebrate cell RNAs, making CpGs a molecular marker of infection. The human zinc-finger antiviral protein (ZAP) recognizes CpG, mediates clearance of the foreign CpG-rich RNA, and causes attenuation of CpG-rich RNA viruses. While ZAP binds RNA, it lacks enzymatic activity that might be responsible for RNA degradation and thus requires interacting cofactors for its function. One of these cofactors, KHNYN, has a predicted nuclease domain. Using biochemical approaches, we found that the KHNYN NYN domain is a single-stranded RNA ribonuclease that does not have sequence specificity and digests RNA with or without CpG dinucleotides equivalently in vitro. We show that unlike most KH domains, the KHNYN KH domain does not bind RNA. Indeed, a crystal structure of the KH region revealed a double-KH domain with a negatively charged surface that accounts for the lack of RNA binding. Rather, the KHNYN C-terminal domain (CTD) interacts with the ZAP RNA-binding domain (RBD) to provide target RNA specificity. We define a minimal complex composed of the ZAP RBD and the KHNYN NYN-CTD and use a fluorescence polarization assay to propose a model for how this complex interacts with a CpG dinucleotide-containing RNA. In the context of the cell, this module would represent the minimum ZAP and KHNYN domains required for CpG-recognition and ribonuclease activity essential for attenuation of viruses with clusters of CpG dinucleotides.

Yeoh, Zoe C. (ORCID:0000000226949068)↗

An integrative analysis of tumor proteomic and phosphoproteomic profiles to examine the relationships between kinase activity and phosphorylation

Phosphorylation of proteins is a key way cells regulate function, both at the individual protein level and at the level of signaling pathways. Kinases are responsible for phosphorylation of substrates, generally on serine, threonine, or tyrosine residues. Though particular sequence patterns can be identified that dictate whether a residue will be phosphorylated by a specific kinase, these patterns are not highly predictive of phosphorylation. The availability of large scale proteomic and phosphoproteomic datasets generated using mass-spectrometry-based approaches provides an opportunity to study the important relationship between kinase activity, substrate specificity, and phosphorylation. In this study, we analyze relationships between protein abundance and phosphopeptide abundance across more than 150 tumor samples and show that phosphorylation at specific phosphosites is not well correlated with overall kinase abundance. However, individual kinases show a clear and statistically significant difference in correlation between known phosphosite targets for that kinase and randomly selected phosphosites. We further investigate relationships between phosphorylation of known activating or inhibitory sites on kinases and phosphorylation of their target phosphosites. Combined with motif-based analysis, this approach can predict novel kinase targets and show which subsets of a kinase’s target repertoire are specifically active in one condition versus another.

post-translational modification, kinase, phosphory↗

Proteomic Assessment of Fluid Shifts and Association with Visual Impairment and Intracranial Pressure in Twin Astronauts

BACKGROUND: Astronauts participating in long duration space missions are at an increased risk of physiological disruptions. The development of visual impairment and intracranial pressure (VIIP) syndrome is one of the leading health concerns for crew members on long-duration space missions; microgravity-induced fluid shifts and chronic elevated cabin CO2 may be contributing factors. By studying physiological and molecular changes in one identical twin during his 1-year ISS mission and his ground-based co-twin, this work extends a current NASA-funded investigation to assess space flight induced "Fluid Shifts" in association with the development of VIIP. This twin study uniquely integrates physiological and -omic signatures to further our understanding of the molecular mechanisms underlying space flight-induced VIIP. We are: (i) conducting longitudinal proteomic assessments of plasma to identify fluid regulation-related molecular pathways altered by long-term space flight; and (ii) integrating physiological and proteomic data with genomic data to understand the genomic mechanism by which these proteomic signatures are regulated. PURPOSE: We are exploring proteomic signatures and genomic mechanisms underlying space flight-induced VIIP symptoms with the future goal of developing early biomarkers to detect and monitor the progression of VIIP. This study is first to employ a male monozygous twin pair to systematically determine the impact of fluid distribution in microgravity, integrating a comprehensive set of structural and functional measures with proteomic, metabolomic and genomic data. This project has a broader impact on Earth-based clinical areas, such as traumatic brain injury-induced elevations of intracranial pressure, hydrocephalus, and glaucoma. HYPOTHESIS: We predict that the space-flown twin will experience a space flight-induced alteration in proteins and peptides related to fluid balance, fluid control and brain injury as compared to his pre-flight protein/peptide signatures. Conversely, the trajectory of these protein signatures will remain relatively constant in his ground based co-twin. METHODS: We are using proteomic and standard immunoelectrophoresis techniques to delineate the change in protein signatures throughout the course of a long duration space flight in relation to the development of VIIP. We are also applying a novel cell-based metaboloic organ system assay ("Organs on a Plate") to address how these circulating biomarkers affect physiological processes at the cellular and organ level which could result in VIIP symptoms. These molecular data will be correlated with physiological measures (eg. extra and intracellular fluid volume, vascular filling/flow patterns, MRI, and Optic Coherence Tomography. DISCUSSION: Pre- and in-flight data collection is in progress for the space-flown twin, and similar data have been obtained from the ground-based twin. Biosamples will be batch processed when received from ISS after the conclusion of the 1-year mission. Omic and Physiological measures from the twin astronauts will be compared to similar data being collected on twin subjects who participated in simulated microgravity study. bed rest study.

Rana, Brinda K.↗

Speciation Underpinned by Unexpected Molecular Diversity in the Mycorrhizal Fungal Genus Pisolithus

The mutualistic ectomycorrhizal (ECM) fungal genus Pisolithus comprises 19 species defined to date which colonize the roots of >50 hosts worldwide suggesting that substantial genomic and functional evolution occurred during speciation. To better understand this intra-genus variation, we undertook a comparative multi-omic study of nine Pisolithus species sampled from North America, South America, Asia, and Australasia. We found that there was a small core set of genes common to all species (13%), and that these genes were more likely to be significantly regulated during symbiosis with a host than accessory or species-specific genes. Thus, the genetic “toolbox” foundational to the symbiotic lifestyle in this genus is small. Transposable elements were located significantly closer to gene classes including effector-like small secreted proteins (SSPs). Poorly conserved SSPs were more likely to be induced by symbiosis, suggesting that they may be a class of protein that tune host specificity. The Pisolithus gene repertoire is characterized by divergent CAZyme profiles when compared with other fungi, both symbiotic and saprotrophic. This was driven by differences in enzymes associated with symbiotic sugar processing, although metabolomic analysis suggest that neither copy number nor expression of these genes is sufficient to predict sugar capture from a host plant or its metabolism in fungal hyphae. Our results demonstrate that intra-genus genomic and functional diversity within ECM fungi is greater than previously thought, underlining the importance of continued comparative studies within the fungal tree of life to refine our focus on pathways and evolutionary processes foundational to this symbiotic lifestyle.

59 BASIC BIOLOGICAL SCIENCES↗

Sugarcane Elongin C is involved in infection by sugarcane mosaic disease pathogens

Sugarcane (Saccharum sp. hybrid) provides the main source of sugar for humans. Sugarcane mosaic disease (SMD) is a major threat to sugarcane production. Currently, control of SMD is mainly dependent on breeding resistant cultivars through hybridization, which is time-consuming. Understanding the mechanism of viral infection may facilitate novel strategies to breed cultivars resistant to SMD and to control the disease. In this study, a wide interaction was detected between the viral VPg protein and host proteins. Several genes were screened from sugarcane cDNA library that could interact with Sugarcane streak mosaic virus VPg, including SceIF4E1 and ScELC. ScELC was predicted to be a cytoplasmic protein, but subcellular localization analysis showed it was distributed both in cytoplasmic and nuclear, and interactions were also detected between ScELC and VPg of SCMV or SrMV that reveal ScELC was widely used in the SMD pathogen infection process. ScELC and VPgs interacted in the nucleus, and may function to enhance the viral transcription rate. ScELC also interacted with SceIF4E2 both in the cytoplasm and nucleus, but not with SceIF4E1 and SceIF4E3. These results suggest that ScELC may be essential for the function of SceIF4E2, an isomer of eIF4E. - Highlights: • We cloned ScELC, SceIF4E1, SceIF4E2 and SceIF4E3 from sugarcane accession Badila. • We examined interactions among VPg, ScELC, SceIF4E1, SceIF4E2 and SceIF4E3. • We proofed that ScELC interacted with VPgs of SCMV, SrMV and SCSMV. • We proofed that ScELC interacted with SceIF4E2 but not SceIF4E1 or SceIF4E3.

60 APPLIED LIFE SCIENCES↗

A functional proteomics platform to reveal the sequence determinants of lysine methyltransferase substrate selectivity

Lysine methylation is a key regulator of histone protein function. Beyond histones, few connections have been made to the enzymes responsible for the deposition of these posttranslational modifications. Here, we debut a high-throughput functional proteomics platform that maps the sequence determinants of lysine methyltransferase (KMT) substrate selectivity without a priori knowledge of a substrate or target proteome. We demonstrate the predictive power of this approach for identifying KMT substrates, generating scaffolds for inhibitor design, and predicting the impact of missense mutations on lysine methylation signaling. By comparing KMT selectivity profiles to available lysine methylome datasets, we reveal a disconnect between preferred KMT substrates and the ability to detect these motifs using standard mass spectrometry pipelines. Collectively, our studies validate the use of this platform for guiding the study of lysine methylation signaling and suggest that substantial gaps exist in proteome-wide curation of lysine methylomes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The lethal cargo of Myxococcus xanthus outer membrane vesicles

Myxococcus xanthus is a bacterial micro-predator known for hunting other microbes in a wolf pack-like manner. Outer membrane vesicles (OMVs) are produced in large quantities by M. xanthus and have a highly organized structure in the extracellular milieu, sometimes occurring in chains that link neighboring cells within a biofilm. OMVs may be a vehicle for mediating wolf pack activity by delivering hydrolytic enzymes and antibiotics aimed at killing prey microbes. Here, both the protein and small molecule cargo of the OMV and membrane fractions of M. xanthus were characterized and compared. Our analysis indicates a number of proteins that are OMV-specific or OMV-enriched, including several with putative hydrolytic function. Secondary metabolite profiling of OMVs identifies 16 molecules, many associated with antibiotic activities. Several hydrolytic enzyme homologs were identified, including the protein encoded by MXAN_3564 (mepA), an M36 protease homolog. Genetic disruption of mepA leads to a significant reduction in extracellular protease activity suggesting MepA is part of the long-predicted (yet to date undetermined) extracellular protease suite of M. xanthus.

59 BASIC BIOLOGICAL SCIENCES↗

Protein Molecular Dynamics Simulations with Approximate QM: What Can We Learn?

Classical force fields are essential for computer simulations of proteins and are typically parameterized to reproduce secondary and tertiary structure of isolated proteins. However, while protein-protein interactions are ubiquitous in nature, they are not considered in parameterization efforts and are far less understood than isolated proteins. A better characterization of intermolecular interactions is widely recognized as a key to revolutionizing drug and therapeutic developments with high-throughput computational screening. Urgently needed is a critical assessment of the performance of modern protein force fields against first-principles electronic structure methods and experiments. In a daring step toward this goal, we here describe a comparison of peptide folding dynamics as predicted by a molecular mechanics force field on the one hand and by an approximate electronic structure quantum mechanical (QM) method based on density-functional tight-binding (DFTB) on the other. We further compare the dynamics from straightforward DFTB simulations with a near-linear scaling version of DFTB for massively parallel computation based on the fragment molecular orbital (FMO-DFTB) method. We illustrate differences between the phenomenology of the folding dynamics from these three methods for a small model peptide, as well as charge polarization and dynamic fluctuations, point out possible correlations and implications for force field developers, and discuss the lessons learned that might become applicable to future predictive high-throughput computer screening for personalized neoantigen cancer therapy.

Irle, Stephan↗

A nucleotide-dependent oligomerization of the Escherichia coli replication initiator DnaA requires residue His136 for remodeling of the chromosomal origin

Escherichia coli replication initiator protein DnaA binds ATP with high affinity but the amount of ATP required to initiate replication greatly exceeds the amount required for binding. Previously, we showed that ATP-DnaA, not ADP-DnaA, undergoes a conformational change at the higher nucleotide concentration, which allows DnaA oligomerization at the replication origin but the association state remains unclear. Here, we used Small Angle X-ray Scattering (SAXS) to investigate oligomerization of DnaA in solution. Whereas ADP-DnaA was predominantly monomeric, AMP–PNP–DnaA (a non-hydrolysable ATP-analog bound-DnaA) was oligomeric, primarily dimeric. Functional studies using DnaA mutants revealed that DnaA(H136Q) is defective in initiating replication in vivo. The mutant retains high-affinity ATP binding, but was defective in producing replication-competent initiation complexes. Docking of ATP on a structure of E. coli DnaA, modeled upon the crystallographic structure of Aquifex aeolicus DnaA, predicts a hydrogen bond between ATP and imidazole ring of His136, which is disrupted when Gln is present at position 136. SAXS performed on AMP–PNP–DnaA (H136Q) indicates that the protein has lost its ability to form oligomers. These results show the importance of high ATP in DnaA oligomerization and its dependence on the His136 residue.

59 BASIC BIOLOGICAL SCIENCES↗

Evidence for an unusual transmembrane configuration of AGG3, a class C Gγ subunit of Arabidopsis

Heterotrimeric G proteins are crucial for the perception of external signals and subsequent signal transduction in animal and plant cells. In both model systems, the complex is comprised of one Gα, one Gβ and one Gγ subunit. However, in addition to the canonical Gγ subunits (Class A), plants also possess two unusual, plant-specific classes of Gγ subunits (Classes B and C) not yet found in animals. These include Gγ subunits lacking the C-terminal CaaX motif (Class B) which is important for membrane anchoring of the protein, and thus give rise to a flexible subpopulation of Gβ/γ heterodimers that is not necessarily restricted to the plasma membrane. Even more interesting, plants also contain Class C Gγ subunits which are twice the size of canonical Gγs, with a predicted transmembrane domain, and a large cysteine-rich, extracellular C-terminus. However, neither the presence of the transmembrane domain nor the membrane topology has been unequivocally demonstrated. Finally, we provide compelling evidence that AGG3, a Class C Ggamma subunit of Arabidopsis, contains a functional transmembrane domain, which is sufficient but not essential for plasma membrane localization, and that the cysteine-rich C-terminus is extracellular.

59 BASIC BIOLOGICAL SCIENCES↗

Different dynamics and pathway of disulfide bonds reduction of two human defensins, a molecular dynamics simulation study

Human defensins are a class of antimicrobial peptides that are crucial components of the innate immune system. Both human α defensin type 5 (HD5) and human β defensin type 3 (hBD-3) have 6 cysteine residues which form 3 pairs of disulfide bonds in oxidizing condition. Disulfide bond linking is important to the protein structure stabilization, and the disulfide bond linking and breaking order have been shown to influence protein function. In this project, microsecond long molecular dynamics simulations were performed to study the structure and dynamics of HD5 and hBD-3 wildtype and analogs which have all 3 disulfide bonds released in reducing condition. The structure of hBD-3 was found to be more dynamic and flexible than HD5, based on RMSD, RMSF, and radius of gyration calculations. The disulfide bridge breaking order of HD5 and hBD-3 in reducing condition was predicted by two kinds of methods, which gave consistent results. It was found that the disulfide bonds breaking pathways for HD5 and hBD-3 are very different. The breaking of disulfide bonds can influence the dimer interface by making the dimer structure less stable for both kinds of defensin. In order to understand the difference in dynamics and disulfide bond breaking pathway, hydrophilic and hydrophobic accessible surface areas (ASA), buried surface area between cysteine pairs, entropy of cysteine pairs, and internal energy were calculated. Comparing to the wildtype, hBD-3 analog is more hydrophobic, while HD5 is more hydrophilic. For hBD-3, the disulfide breaking is mainly entropy driven, while other factors such as the solvation effects may take the major role in controlling HD5 disulfide breaking pathway. Proteins 2017; 85:665–681. © 2016 Wiley Periodicals, Inc.

Biochemistry & Molecular Biology↗

Transcriptome sequencing and differential gene expression analysis in Viola yedoensis Makino (Fam. Violaceae) responsive to cadmium (Cd) pollution

Viola yedoensis Makino is an important Chinese traditional medicine plant adapted to cadmium (Cd) pollution regions. Illumina sequencing technology was used to sequence the transcriptome of V. yedoensis Makino. We sequenced Cd-treated (VIYCd) and untreated (VIYCK) samples of V. yedoensis, and obtained 100,410,834 and 83,587,676 high quality reads, respectively. After de novo assembly and quantitative assessment, 109,800 unigenes were finally generated with an average length of 661 bp. We then obtained functional annotations by aligning unigenes with public protein databases including NR, NT, SwissProt, KEGG and COG. In addition, 892 differentially expressed genes (DEGs) were investigated between the two libraries of untreated (VIYCK) and Cd-treated (VIYCd) plants. Moreover, 15 randomly selected DEGs were further validated with qRT-PCR and the results were highly accordant with the Solexa analysis. This study firstly generated a successful global analysis of the V. yedoensis transcriptome and it will provide for further studies on gene expression, genomics, and functional genomics in Violaceae. - Highlights: • A de novo assembly generated 109,800 unigenes and 5,4479 of them were annotated. • 31,285 could be classified into 26 COG categories. • 263 biosynthesis pathways were predicted and classified into five categories. • 892 DEGs were detected and 15 of them were validated by qRT-PCR.

60 APPLIED LIFE SCIENCES↗

Structural and biochemical basis for regiospecificity of the flavonoid glycosyltransferase UGT95A1

Glycosylation is a predominant strategy plants use to fine-tune the properties of small molecule metabolites to affect their bioactivity, transport, and storage. It is also important in biotechnology and medicine as many glycosides are utilized in human health. Small molecule glycosylation is largely carried out by family 1 glycosyltransferases. Here, we report a structural and biochemical investigation of UGT95A1, a family 1 GT enzyme from Pilosella officinarum that exhibits a strong, unusual regiospecificity for the 3'-O position of flavonoid acceptor substrate luteolin. We obtained an apo crystal structure to help drive the analyses of a series of binding site mutants, revealing that while most residues are tolerant to mutations, key residues M145 and D464 are important for overall glycosylation activity. Interestingly, E347 is crucial for maintaining the strong preference for 3'-O glycosylation, while R462 can be mutated to increase regioselectivity. The structural determinants of regioselectivity were further confirmed in homologous enzymes. Our study also suggests that the enzyme contains large, highly dynamic, disordered regions. We showed that while most disordered regions of the protein have little to no implication in catalysis, the disordered regions conserved among investigated homologs are important to both the overall efficiency and regiospecificity of the enzyme. This report represents a comprehensive in-depth analysis of a family 1 GT enzyme with a unique substrate regiospecificity and may provide a basis for enzyme functional prediction and engineering.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Combining Geometric Pocket Detection and Desolvation Properties to Detect Putative Ligand Binding Sites on Proteins

The accurate identification of cavities that can bind ligands on the surface of proteins is of major importance for the characterization of the function of proteins based on its structure. In addition it can be helpful for rational structure-based drug design on target proteins of medical relevance and for evaluating the tendency of proteins to aggregate or oligomerize. A new approach termed dPredGB to detect and evaluate putativebinding cavities on protein surfaces has been developed. In contrast to existing prediction methods that are based on purely geometric features of binding sites or on possible direct interactions with a putativebinding partner the dPredGB approach combines rapid geometricdetection with an evaluation of the desolvationproperties of the putativebindingpocket. It has been tested on a variety of proteins known to bind ligands in bound and unbound conformations. The approach outperforms most available methods and offers also the spatial characterization of the desolvationproperties of a binding region. On a test set of proteins the method identifies in 69% of the unbound cases and 85% of the bound cases the known ligandbinding cavity as the top ranking prediction. Possibilities to improve the prediction performance even further are also discussed.

Environmental Molecular Sciences Laboratory↗

Aberrant FBXW7-mediated ubiquitination and degradation of ZMYND8 enhances tumor progression and stemness in bladder cancer

ZMYND8, an epigenetic regulator, was identified as a common oncogene across various tumors. However, little was reported about the association between ZMYND8 and bladder cancer. Besides, aberrant mechanisms that contribute to abnormal ZMYND8 expressions still remain unclear. In the current study, we first found that ZMYND8 protein levels were significantly elevated in Bca samples versus normal tissues, but not the mRNA levels. We then utilized the Cell Counting Kit-8 (CCK-8) assay, clone formation assay and transwell analysis to confirm that ZMYND8 could remarkably promote the tumor progression in vitro, including growth capacity and migration. Bioinformatic predictive analysis revealed that E3 ubiquitin ligase FBXW7 interacts directly with ZMYND8 and degrades ZMYND8 in a polyubiquitination manner. Low FBXW7 was a hazard factor for promoting and depending on accumulated ZMYND8 proteins to promote Bca progression. Gene set enrichment analysis (GSEA) further indicated that ZMYND8 was notably associated with stemness process, which was well functionally validated. Lastly, ZMYND8 deficiency was observed to inhibit tumor growth of Bca in vivo, revealing a promising translational significance in Bca treatment. In conclusion, our study for the first time provided evidence for a novel mechanism of FBXW7/ZMYND8 axis in Bca, providing therapeutic vulnerability for individualized cancer treatment.

60 APPLIED LIFE SCIENCES↗

Ecosystems and Networks Integrated with Genes and Molecular Assemblies (ENIGMA): Component 5: Imaging Protein Conformations, Shapes & Assemblies in Solution & Administration project (Final Scientific/Technical Report, Subcontract No. 6974584)

We set ambitious goals to examine microorganism communities and measure both their chemical input and output as a read out of specific biochemical activity. These scientific goals are driving the development of sophisticated algorithms to analyze large amounts of experimental measurements made using high throughput technologies to explain and predict how the environment influences biological function at multiple scales and how the microbial systems in tum modify the environment. By examining how bacteria communities rely on symbiotic metabolic relationships for survival and reproductive success, and how these relationships consequently affect their biochemical capabilities. This was accomplished using state-of-the-art mass spectrometry-based methods, and metabolic fingerprinting approaches with a high degree of chemical specificity and sensitivity. During this period the original effort transitioned and was consolidated into what is now ENIGMA, the efforts on technology development expanded to include more untargeted metabolomics and its application to organisms on a systems level.

59 BASIC BIOLOGICAL SCIENCES↗

Metapangenomics reveals depth-dependent shifts in metabolic potential for the ubiquitous marine bacterial SAR324 lineage

Abstract Background Oceanic microbiomes play a pivotal role in the global carbon cycle and are central to the transformation and recycling of carbon and energy in the ocean’s interior. SAR324 is a ubiquitous but poorly understood uncultivated clade of Deltaproteobacteria that inhabits the entire water column, from ocean surface waters to its deep interior. Although some progress has been made in elucidating potential metabolic traits of SAR324 in the dark ocean, very little is known about the ecology and the metabolic capabilities of this group in the euphotic and twilight zones. To investigate the comparative genomics, ecology, and physiological potential of the SAR324 clade, we examined the distribution and variability of key genomic features and metabolic pathways in this group from surface waters to the abyss in the North Pacific Subtropical Gyre, one of the largest biomes on Earth. Results We leveraged a pangenomic ecological approach, combining spatio-temporally resolved single-amplified genome, metagenomic, and metatranscriptomic datasets. The data revealed substantial genomic diversity throughout the SAR324 clade, with distinct depth and temporal distributions that clearly differentiated ecotypes. Phylogenomic subclade delineation, environmental distributions, genomic feature similarities, and metabolic capacities revealed strong congruence. The four SAR324 ecotypes delineated in this study revealed striking divergence from one another with respect to their habitat-specific metabolic potentials. The ecotypes living in the dark or twilight oceans shared genomic features and metabolic capabilities consistent with a sulfur-based chemolithoautotrophic lifestyle. In contrast, those inhabiting the sunlit ocean displayed higher plasticity energy-related metabolic pathways, supporting a presumptive photoheterotrophic lifestyle. In epipelagic SAR324 ecotypes, we observed the presence of two types of proton-pumping rhodopsins, as well as genomic, transcriptomic, and ecological evidence for active photoheterotrophy, based on xanthorhodopsin-like light-harvesting proteins. Conclusions Combining pangenomic and both metagenomic and metatranscriptomic profiling revealed a striking divergence in the vertical distribution, genomic composition, metabolic potential, and predicted lifestyle strategies of geographically co-located members of the SAR324 bacterial clade. The results highlight the utility of metapangenomic approaches employed across environmental gradients, to decipher the properties and variation in function and ecological traits of specific phylogenetic clades within complex microbiomes.

59 BASIC BIOLOGICAL SCIENCES↗