Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “structural genomics”

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 667 records · Page 37

Specialization Restricts the Evolutionary Paths Available to Yeast Sugar Transporters

Functional innovation at the protein level is a key source of evolutionary novelties. The constraints on functional innovations are likely to be highly specific in different proteins, which are shaped by their unique histories and the extent of global epistasis that arises from their structures and biochemistries. These contextual nuances in the sequence–function relationship have implications both for a basic understanding of the evolutionary process and for engineering proteins with desirable properties. Here, we have investigated the molecular basis of novel function in a model member of an ancient, conserved, and biotechnologically relevant protein family. These Major Facilitator Superfamily sugar porters are a functionally diverse group of proteins that are thought to be highly plastic and evolvable. By dissecting a recent evolutionary innovation in an α-glucoside transporter from the yeast Saccharomyces eubayanus, we show that the ability to transport a novel substrate requires high-order interactions between many protein regions and numerous specific residues proximal to the transport channel. To reconcile the functional diversity of this family with the constrained evolution of this model protein, we generated new, state-of-the-art genome annotations for 332 Saccharomycotina yeast species spanning ~400 My of evolution. By integrating phylogenetic and phenotypic analyses across these species, we show that the model yeast α-glucoside transporters likely evolved from a multifunctional ancestor and became subfunctionalized. The accumulation of additive and epistatic substitutions likely entrenched this subfunction, which made the simultaneous acquisition of multiple interacting substitutions the only reasonably accessible path to novelty.

59 BASIC BIOLOGICAL SCIENCES↗

Genome-Resolved Metagenomics and Detailed Geochemical Speciation Analyses Yield New Insights into Microbial Mercury Cycling in Geothermal Springs

Geothermal systems emit substantial amounts of aqueous, gaseous, and methylated mercury, but little is known about microbial influences on mercury speciation. Here, we report results from genome-resolved metagenomics and mercury speciation analysis of acidic warm springs in the Ngawha Geothermal Field (< 55°C, pH < 4.5), Northland Region, Aotearoa New Zealand. Our aim was to identify the microorganisms genetically equipped for mercury methylation, demethylation, or Hg(II) reduction to volatile Hg(0) in these springs. Dissolved total and methylated mercury concentrations in two adjacent springs with different mercury speciation ranked among the highest reported from natural sources (250 to 16,000 ng liter -1 and 0.5 to 13.9 ng liter -1 , respectively). Total solid mercury concentrations in spring sediments ranged from 1,274 to 7,000 μg g -1 . In the context of such ultrahigh mercury levels, the geothermal microbiome was unexpectedly diverse and dominated by acidophilic and mesophilic sulfur- and iron-cycling bacteria, mercury- and arsenic-resistant bacteria, and thermophilic and acidophilic archaea. By integrating microbiome structure and metagenomic potential with geochemical constraints, we constructed a conceptual model for biogeochemical mercury cycling in geothermal springs. The model includes abiotic and biotic controls on mercury speciation and illustrates how geothermal mercury cycling may couple to microbial community dynamics and sulfur and iron biogeochemistry.

59 BASIC BIOLOGICAL SCIENCES↗

Nodavirus RNA replication crown architecture reveals proto-crown precursor and viral protein A conformational switching

Positive-strand RNA viruses replicate their genomes in virus-induced membrane vesicles, and the resulting RNA replication complexes are a major target for virus control. Nodavirus studies first revealed viral RNA replication proteins forming a 12-fold symmetric “crown” at the vesicle opening to the cytosol, an arrangement recently confirmed to extend to distantly related alphaviruses. Using cryoelectron microscopy (cryo-EM), we show that mature nodavirus crowns comprise two stacked 12-mer rings of multidomain viral RNA replication protein A. Each ring contains an ~19 nm circle of C-proximal polymerase domains, differentiated by strikingly diverged positions of N-proximal RNA capping/membrane binding domains. The lower ring is a “proto-crown” precursor that assembles prior to RNA template recruitment, RNA synthesis, and replication vesicle formation. In this proto-crown, the N-proximal segments interact to form a toroidal central floor, whose 3.1 Å resolution structure reveals many mechanistic details of the RNA capping/membrane binding domains. In the upper ring, cryo-EM fitting indicates that the N-proximal domains extend radially outside the polymerases, forming separated, membrane-binding “legs.” The polymerase and N-proximal domains are connected by a long linker accommodating the conformational switch between the two rings and possibly also polymerase movements associated with RNA synthesis and nonsymmetric electron density in the lower center of mature crowns. The results reveal remarkable viral protein multifunctionality, conformational flexibility, and evolutionary plasticity and insights into (+)RNA virus replication and control.

59 BASIC BIOLOGICAL SCIENCES↗

Structure and flexibility of the yeast NuA4 histone acetyltransferase complex

The NuA4 protein complex acetylates histones H4 and H2A to activate both transcription and DNA repair. We report the 3.1-Å resolution cryo-electron microscopy structure of the central hub of NuA4, which flexibly tethers the histone acetyltransferase (HAT) and Trimer Independent of NuA4 involved in Transcription Interactions with Nucleosomes (TINTIN) modules. The hub contains the large Tra1 subunit and a core that includes Swc4, Arp4, Act1, Eaf1, and the C-terminal region of Epl1. Eaf1 stands out as the primary scaffolding factor that interacts with the Tra1, Swc4, and Epl1 subunits and contributes the conserved HSA helix to the Arp module. Using nucleosome-binding assays, we find that the HAT module, which is anchored to the core through Epl1, recognizes H3K4me3 nucleosomes with hyperacetylated H3 tails, while the TINTIN module, anchored to the core via Eaf1, recognizes nucleosomes that have hyperacetylated H2A and H4 tails. Together with the known interaction of Tra1 with site-specific transcription factors, our data suggest a model in which Tra1 recruits NuA4 to specific genomic sites then allowing the flexible HAT and TINTIN modules to select nearby nucleosomes for acetylation.

59 BASIC BIOLOGICAL SCIENCES↗

Label-free structural imaging of plant roots and microbes using third-harmonic generation microscopy

Root biology is pivotal in addressing global challenges including sustainable agriculture and climate change. However, roots have been relatively understudied among plant organs, partly due to the difficulties in imaging root structures in their natural environment. Here we used microfabricated ecosystems (EcoFABs) to establish growing environments with optical access and employed nonlinear multimodal microscopy of third-harmonic generation (THG) and three-photon fluorescence (3PF) to achieve label-free, in situ imaging of live roots and microbes at high spatiotemporal resolution. THG enabled us to observe key plant root structures including the vasculature, Casparian strips, dividing meristematic cells, and root cap cells, as well as subcellular features including nuclear envelopes, nucleoli, starch granules, and putative stress granules. THG from the cell walls of bacteria and fungi also provides label-free contrast for visualizing these microbes in the root rhizosphere. With simultaneously recorded 3PF signal, we demonstrated our ability to investigate root-microbe interactions by achieving single-bacterium tracking and subcellular imaging of fungal spores and hyphae in the rhizosphere.

Pan, Daisong [University of California, Berkeley, ↗

Classical Swine Fever Virus Structural Glycoprotein E2 Interacts with Host Protein ACADM during the Virus Infectious Cycle

The E2 glycoprotein is one of the four structural proteins of the classical swine fever virus (CSFV) particle. E2 has been shown to be involved in many virus functions, including adsorption to host cells, virus virulence and interaction with several host proteins. Using a yeast two-hybrid screen, we have previously shown that the CSFV E2 specifically interacts with swine host protein medium-chain-specific acyl-Coenzyme A dehydrogenase (ACADM), an enzyme that catalyzes the initial step of the mitochondrial fatty acid beta-oxidation pathway. Here, we show that interaction between ACADM and E2 also happens in swine cells infected with CSFV using two different procedures: coimmunoprecipitation and a proximity ligation assay (PLA). In addition, the amino acid residues in E2 critically mediating the interaction with ACADM, M49 and P130 were identified via a reverse yeast two-hybrid screen using an expression library composed of randomly mutated versions of E2. A recombinant CSFV, E2ΔACADMv, harboring substitutions at residues M49I and P130Q in E2, was developed via reverse genomics from the highly virulent Brescia isolate. E2ΔACADMv was shown to have the same kinetics growth in swine primary macrophages and SK6 cell cultures as the parental Brescia strain. Similarly, E2ΔACADMv demonstrated a similar level of virulence when inoculated to domestic pigs as the parental Brescia. Animals intranasally inoculated with 10 5 TCID 50 developed a lethal form of clinical disease with virological and hematological kinetics changes undistinguishable from those produced by the parental strain. Therefore, interaction between CSFV E2 and host ACADM is not critically involved in the processes of virus replication and disease production.

60 APPLIED LIFE SCIENCES↗

HydraGNN_Predictive_GFM_2026 - Ensemble of predictive graph foundation models for atomistic materials modeling

This release contains data and parameters of HydraGNN-based graph foundation models trained as a result of the work published in the pre-print "Exascale Multi-Task Graph Foundation Models for Imbalanced, Multi-Fidelity Atomistic Data" by M. Lupo Pasini et al. (https://arxiv.org/abs/2604.15380). We jointly train on 16 open first-principles datasets (544+ million structures covering 85+ elements) using a multi-task architecture with per-dataset heads and a scalable ADIOS2/DDStore data pipeline. On Frontier, we execute six large-scale DeepHyper hyperparameter optimization campaigns in FP64 and promote the top-performing message-passing models to sustained 2,048-node training, yielding a PaiNN-based lead model. The version of HydraGNN used to generate the outputs provided in this release is HydraGNN v5.0 (https://github.com/ORNL/HydraGNN/releases/tag/v5.0) The list of datasets used for the training of the graph foundation model is the following: 1) Alexandria [1] 2) ANI1x [2] 3) MPTrj [3] 4) Open Catalyst 2020 (OC20) [4] 5) Open Catalyst 2022 (OC22) [5] 6) Open Catalyst 2025 (OC25) [6] 7) Open Direct ir Capture 2023 (ODAC23) [7] 8) Open Materials 2024 (OMat24) [8] 9) Open Molecules 2025 (OMol25) [9] 10) OMol25-neutral (subset of OMol25 that contains only molecules with zero total charge) 11) OMol25-non-neutral (subset of OMol25 that contains only molecules with non-zero total charge) 12) Open Polymers 2026 (OPoly2026) [10] 13) Nabla2DFT [11] 14) QCML [12] 15) QM7X [reference 13] 16) transition1x [14] Dataset references: [1] J. Schmidt et al., “A dataset of 175k stable and metastable materials calculated with the PBEsol and SCAN functionals,” Scientific Data, vol. 9, p. 64, 2022. [2] J. S. Smith et al., “The ANI-1ccx and ANI-1x data sets, coupled-cluster and density functional theory properties for molecules,” Scientific Data, vol. 7, p. 134, 2020. [Online]. Available: https: //www.nature.com/articles/s41597-020-0473-z [3] A. Jain et al., “Commentary: The Materials Project: A materials genome approach to accelerating materials innovation,” APL Materials, vol. 1, no. 1, p. 011002, 07 2013. [Online]. Available: https://doi.org/10.1063/1.4812323 [4] L. Chanussot et al., “Open catalyst 2020 (oc20) dataset and community challenges,” ACS Catalysis, vol. 11, no. 10, pp. 6059–6072, 2021. [Online]. Available: https://doi.org/10.1021/acscatal.0c04525 [5] K. Tran et al., “Open catalyst 2022 (oc22) dataset and challenges for oxidation electrocatalysts,” ACS Catalysis, vol. 13, no. 5, pp. 3066–3084, 2023. [Online]. Available: https://doi.org/10.1021/acscatal.2c05426 [6] S. J. Sahoo et al., “The open catalyst 2025 (oc25) dataset and models for solid-liquid interfaces,” arXiv preprint arXiv:2509.17862, 2025. [Online]. Available: https://arxiv.org/abs/2509.17862 [7] A. Sriram et al., “The open DAC 2023 dataset and challenges for sorbent discovery in direct air capture,” ACS Central Science, vol. 10, no. 5, pp. 923–941, 2024. [8] L. Barroso-Luque et al., “Open materials 2024 (omat24) inorganic materials dataset and models,” 2024. [Online]. Available: https://arxiv.org/abs/2410.12771 [9] D. S. Levine et al., “The open molecules 2025 (OMol25) dataset, evaluations, and models,” 2025. [Online]. Available: https://arxiv.org/abs/2505.08762 [10] D. S. Levine et al., The open polymers 2026 (OPoly26) dataset and evaluations,” arXiv preprint arXiv:2512.23117, 2025. [Online]. Available: https://arxiv.org/abs/2512.23117 [11] K. Khrabrov et al., “Nabla2dft: A universal quantum chemistry dataset of drug-like molecules and a benchmark for neural network potentials,” in NeurIPS 2024 Datasets and Benchmarks Track, 2024. [Online]. Available: https://openreview.net/forum?id=ElUrNM9U8c [12] S. Ganscha et al., “The QCML dataset, quantum chemistry reference data from 33.5M DFT and 14.7B semi-empirical calculations,” Scientific Data, vol. 12, p. 406, 2025. [13] J. Hoja et al., “QM7-X, a comprehensive dataset of quantum-mechanical properties spanning the chemical space of small organic molecules,” Scientific Data, vol. 8, p. 43, 2021. [Online]. Available: https://www.nature.com/articles/s41597-021-00812-2 [14] M. Schreiner et al., “Transition1x - a dataset for building generalizable reactive machine learning potentials,” Scientific Data, vol. 9, p. 779, 2022. The folder "datasets_ADIOS2_format" contains the set of pre-processed datasets in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used for the development and training of GFMs in this work. The "datasets_ADIOS2_format" directory contains 2 sub-directories, one for the version "v1" of the datasets and one for the version "v2" of the datasets. The version "v1" of the datasets provides values of the total energy as they are extracted from the original data as it was released by the respective institutions. The version "v2" of the datasets provides values of the energy that have been realigned. The realignment was performed by training a linear regression model that predicts the total energy as a function of the chemical composition of the atomistic structure, and then subtract such prediction from the original value of the total energy. Both folders "v1" and "v2" contain 16 sub-directories, each corresponding to an ADIOS2-formatted dataset The folder "DeepHyper-results" contains the configurational files and model's parameters for all the 186 HPO trials that were successfully completed by the scalable hyperparameter optimization (HPO) runs on Frontier. The content of the folder "DeepHyper-results" I structured as follows: 1) task-list.txt: list of mpnn name, jobid, and deephyper task id 2) gfm_${MPNN}_${JOBID}_0.${TASKID}: run directory with checkpoint files 3) gfm_${MPNN}: deephyper summary directory (*.csv) for each specific MPNN type 4) deephyper-experiment-${JOBID}: output and error logs for each job The file "deephyper-sorted.csv" contains the details of each HydraGNN model built and tested by HPO, obtained by merging the (*.csv) filed from each HPO run executed. Out of all the HPO trials, we selected 10 to continue the training of the respective HydraGNN models. Due to limited computational budget available in the LRN070 allocation we could not complete the training till convergence for all these 10 selected models. The folder "models" contains multiple sub-folders, one per each HydraGNN model trained. Each model sub-folder contains the parameters of each HydraGNN model, with multiple checkpoint-restarts. The list of sub-folders are as follows: 1) multidataset_hpo-BEST1-fp64 2) multidataset_hpo-BEST2-fp64 3) multidataset_hpo-BEST3-fp64 4) multidataset_hpo-BEST4-fp64 5) multidataset_hpo-BEST5-fp64 6) multidataset_hpo-BEST6-fp64 7) multidataset_hpo-BEST7-fp64 8) multidataset_hpo-BEST8-fp64 9) multidataset_hpo-BEST9-fp64 10) multidataset_hpo-BEST10-fp64 Within each one of these folders, additional auxiliary log files are provided with descriptions about how the training proceeded. The lead PaiNN-model is contained inside "multidataset_hpo-BEST6-fp64". The file "mlp_branch_weights" contains the parameters of the multi-layer perceptron (MLP) used to reconcile the predictions of the 16 output decoding heads of the HydragNN architectures. The MLP takes in input the chemical composition of the atomistic structure and predicts averaging weights to linearly mix the predictions of each output decoding head toward consolidating them into a single one. The folder "1.1billion-structure-inference" contains 1.1 billion atomistic structures randomly generated. Each structures is associated with energy and forces predicted with the lead-PaiNN model combined with the MLP model for reconciliation of the multi-branch predictions generated by the 16 output decoding heads. The folder "1.1billion-structure-inference" contains 9,300 (*.tar.gz) subdirectories, one per Frontier compute node used to execute the inference at exascale. Once uncompressed, each (*.tar.gz) subdirectory contains an ADIOS2 (*.bp) file container, where each atomistic structure is stored as a PyTorch-Geometric Data object. The file "export_dataset_environment_variables.sh" contains the environment variables that need to be set before running the HydraGNN code to reproduce the results provided in this dataset release. The code that can be used to load the ADIOS2 files, load HydraGNN models, and run inference is available at: https://github.com/ORNL/HydraGNN/releases/tag/v5.0

36 MATERIALS SCIENCE↗

DOE new players Carbon cycle (2016-2021)

The overarching scientific goals of our multidisciplinary grant was to further expand understanding about the key microorganisms (players), metabolic strategies (processes), and interspecies relationships (interactions) involved in the formation and oxidation of methane in the environment. This research applied novel environmental metagenomics, transcriptomics, and proteomic techniques, state-of-the-art analytical imaging, stable isotope geochemistry, and reaction-transport modeling to address these goals and develop an ‘ecosystems level’ understanding of the factors which regulate microbial methane cycling in anoxic sedimentary ecosystems. For decades, it was believed that the obligate step in methanogenesis catalyzed by methyl coenzyme M reductase (Mcr) was limited to a specific branch of the archaeal Domain, formerly known as the Euryarchaeota. Less than a decade ago co-I Tyson’s team published a surprising metagenomic-based discovery of divergent Mcr genes in a novel uncultured phylum (Bathyarchaeota), catalyzing a major shift in thinking about the diversity of microorganisms that encode the potential for methane (or higher alkane) metabolism in anoxic environments (Evans et al., 2015). In our work here, we further expand on the groups of archaea harboring the genomic potential for methane or hydrocarbon metabolism using environmental metagenomics and new gene targeted bioinformatics techniques. We additionally advanced understanding about the terminal electron acceptors and metabolic potential supporting the anaerobic oxidation of methane (AOM) in terrestrial ecosystems, specifically focused on new lineages of ANME archaea capable of respiring manganese oxides with methane presumably using large extracellular multi-heme cytochrome complexes. New details about specific syntrophic mechanisms underlying the exchange of electrons during sulfate-coupled methane oxidation between ANME-2 archaea and their sulfate-reducing bacterial partners were also elucidated as part of this funded project. Through a series of experimental ‘omics and single cell stable isotope probing studies with incubated environmental sediment samples and a cultured model electrogenic microorganism combined with model-based predictions. Combined, this work provides strong support for the hypothesis of direct interspecies electron transfer (DIET) is the dominant syntrophic mechanism controlling the anaerobic oxidation of methane with sulfate over other proposed mechanisms and additionally illustrates important spatial constraints and the underlying physico-chemical factors influencing AOM syntrophic consortia structure for DIET and extracellular metal respiration. This collaborative multi-institutional project successfully advanced several of our milestone goals including the identification of new microbial players containing methyl coenzyme M reductases hypothesized to be central to methane or hydrocarbon cycling in anoxic environments and enhancing fundamental knowledge about the role extracellular electron transfer plays in the ecophysiology of methanotrophic archaea respiring metal oxides and in the physical and metabolic structuring of syntrophic interactions in methane-rich sedimentary ecosystems.

03 NATURAL GAS↗

Data for Spatial Analysis of Cell Patterning to Aid Genetic and Phenotypic Understanding of Grass Stomatal Density: A Case Study in Maize

Biological processes involve complex hierarchies where composite traits result from multiple component traits. However, holistically understanding of how sets of component traits interact to underpin genotype-to-phenotype relationships is generally lacking. Stomatal density (SD) is a tractable model system for exploring how high-throughput phenotyping (HTP) data could be exploited by a new spatial analysis approach to better understand a developmentally and functionally important trait. SD is a composite trait, resulting from various components related to cell identity and size, which are themselves governed by a series of spatio-developmental processes. Data from 192 recombinant inbred lines of maize [Zea mays (L.)] were analyzed by a new stomatal patterning phenotype (SPP) to (1) describe the average spatial probability distribution of the nearest neighboring stomata; (2) derive a core set of component traits related to cell size, cell packing, and positional probabilities; (3) build a structural equation model of component traits underlying SD; and (4) identify stomatal patterning quantitative trait loci (QTL). The core set of SPP-derived traits explained 74% of the variation in SD. Analyzing SPP component traits allowed some loci previously identified as generic SD QTL to be recognized as specific to lateral versus longitudinal elements of stomatal patterning. Therefore, this study highlights how novel insights can be gained by decomposing a composite trait (e.g., SD) into a set of component traits that were present in HTP data but not previously exploited.

AI/ML↗

Genetic control of photoprotection and photosystem II operating efficiency in plants

Summary Photoprotection against excess light via nonphotochemical quenching (NPQ) is indispensable for plant survival. However, slow NPQ relaxation under low light conditions can decrease yield of field‐grown crops up to 40%. Using semi‐high‐throughput assay, we quantified the kinetics of NPQ and photosystem II operating efficiency (ΦPSII) in a replicated field trial of more than 700 maize ( Zea mays ) genotypes across 2 yr. Parametrized kinetics data were used to conduct genome‐wide association studies. For six candidate genes involved in NPQ and ΦPSII kinetics in maize the loss of function alleles of orthologous genes in Arabidopsis ( Arabidopsis thaliana ) were characterized: two thioredoxin genes, and genes encoding a transporter in the chloroplast envelope, an initiator of chloroplast movement, a putative regulator of cell elongation and stomatal patterning, and a protein involved in plant energy homeostasis. Since maize and Arabidopsis are distantly related, we propose that genes involved in photoprotection and PSII function are conserved across vascular plants. The genes and naturally occurring functional alleles identified here considerably expand the toolbox to achieving a sustainable increase in crop productivity.

59 BASIC BIOLOGICAL SCIENCES↗

Genome-wide identification of novel flagellar motility genes in Pseudomonas syringae pv. tomato DC3000

Pseudomonas syringaepv.tomatoDC3000 (PstDC3000) is a plant pathogenic bacterium that possesses complicated motility regulation pathways including a typical chemotaxis system. A significant portion of our understanding about the genes functioning inPstDC3000 motility is based on comparison to other bacteria. This leaves uncertainty about whether gene functions are conserved, especially since specific regulatory modules can have opposite functions in sets ofPseudomonas. In this study, we used a competitive selection to enrich for mutants with altered swimming motility and used random barcode transposon-site sequencing (RB-TnSeq) to identify genes with significant roles in swimming motility. Besides many of the known or predicted chemotaxis and motility genes, our method identified PSPTO_0406 (dipA), PSPTO_1042 (chrR) and PSPTO_4229 (hypothetical protein) as novel motility regulators. PSPTO_0406 is a homolog ofdipA, a known cyclic di-GMP degrading enzyme inP. aeruginosa. PSPTO_1042 is part of an extracytoplasmic sensing system that controls gene expression in response to reactive oxygen species, suggesting that PSPTO_1042 may function as part of a mechanism that enablesPstDC3000 to alter motility when encountering oxidative stressors. PSPTO_4229 encodes a protein containing an HD-related output domain (HDOD), but with no previously identified functions. We found that deletion and overexpression of PSPTO_4229 both reduce swimming motility, suggesting that its function is sensitive to expression level. We used the overexpression phenotype to screen for nonsense and missense mutants of PSPTO_4229 that no longer reduce swimming motility and found a pair of conserved arginine residues that are necessary for motility suppression. Together these results provide a global perspective on regulatory and structural genes controlling flagellar motility inPstDC3000.

Microbiology↗

Describing patterns of familial cancer risk in subfertile men using population pedigree data

STUDY QUESTION Can we simultaneously assess risk for multiple cancers to identify familial multicancer patterns in families of azoospermic and severely oligozoospermic men? SUMMARY ANSWER Here, distinct familial cancer patterns were observed in the azoospermia and severe oligozoospermia cohorts, suggesting heterogeneity in familial cancer risk by both type of subfertility and within subfertility type. WHAT IS KNOWN ALREADY Subfertile men and their relatives show increased risk for certain cancers including testicular, thyroid, and pediatric. STUDY DESIGN, SIZE, DURATION A retrospective cohort of subfertile men (N = 786) was identified and matched to fertile population controls (N = 5674). Family members out to third-degree relatives were identified for both subfertile men and fertile population controls (N = 337 754). The study period was 1966–2017. Individuals were censored at death or loss to follow-up, loss to follow-up occurred if they left Utah during the study period. PARTICIPANTS/MATERIALS, SETTING, METHODS Azoospermic (0 × 10 6 /mL) and severely oligozoospermic (<1.5 × 10 6 /mL) men were identified in the Subfertility Health and Assisted Reproduction and the Environment cohort (SHARE). Subfertile men were age- and sex-matched 5:1 to fertile population controls and family members out to third-degree relatives were identified using the Utah Population Database (UPDB). Cancer diagnoses were identified through the Utah Cancer Registry. Families containing ≥10 members with ≥1 year of follow-up 1966–2017 were included (azoospermic: N = 426 families, 21 361 individuals; oligozoospermic: N = 360 families, 18 818 individuals). Unsupervised clustering based on standardized incidence ratios for 34 cancer phenotypes in the families was used to identify familial multicancer patterns; azoospermia and severe oligospermia families were assessed separately. MAIN RESULTS AND THE ROLE OF CHANCE Compared to control families, significant increases in cancer risks were observed in the azoospermia cohort for five cancer types: bone and joint cancers hazard ratio (HR) = 2.56 (95% CI = 1.48–4.42), soft tissue cancers HR = 1.56 (95% CI = 1.01–2.39), uterine cancers HR = 1.27 (95% CI = 1.03–1.56), Hodgkin lymphomas HR = 1.60 (95% CI = 1.07–2.39), and thyroid cancer HR = 1.54 (95% CI = 1.21–1.97). Among severe oligozoospermia families, increased risk was seen for three cancer types: colon cancer HR = 1.16 (95% CI = 1.01–1.32), bone and joint cancers HR = 2.43 (95% CI = 1.30–4.54), and testis cancer HR = 2.34 (95% CI = 1.60–3.42) along with a significant decrease in esophageal cancer risk HR = 0.39 (95% CI = 0.16–0.97). Thirteen clusters of familial multicancer patterns were identified in families of azoospermic men, 66% of families in the azoospermia cohort showed population-level cancer risks, however, the remaining 12 clusters showed elevated risk for 2-7 cancer types. Several of the clusters with elevated cancer risks also showed increased odds of cancer diagnoses at young ages with six clusters showing increased odds of adolescent and young adult (AYA) diagnosis [odds ratio (OR) = 1.96–2.88] and two clusters showing increased odds of pediatric cancer diagnosis (OR = 3.64–12.63). Within the severe oligozoospermia cohort, 12 distinct familial multicancer clusters were identified. All 12 clusters showed elevated risk for 1–3 cancer types. An increase in odds of cancer diagnoses at young ages was also seen in five of the severe oligozoospermia familial multicancer clusters, three clusters showed increased odds of AYA diagnosis (OR = 2.19–2.78) with an additional two clusters showing increased odds of a pediatric diagnosis (OR = 3.84–9.32). LIMITATIONS, REASONS FOR CAUTION Although this study has many strengths, including population data for family structure, cancer diagnoses and subfertility, there are limitations. First, semen measures are not available for the sample of fertile men. Second, there is no information on medical comorbidities or lifestyle risk factors such as smoking status, BMI, or environmental exposures. Third, all of the subfertile men included in this study were seen at a fertility clinic for evaluation. These men were therefore a subset of the overall population experiencing fertility problems and likely represent those with the socioeconomic means for evaluation by a physician. WIDER IMPLICATIONS OF THE FINDINGS This analysis leveraged unique population-level data resources, SHARE and the UPDB, to describe novel multicancer clusters among the families of azoospermic and severely oligozoospermic men. Distinct overall multicancer risk and familial multicancer patterns were observed in the azoospermia and severe oligozoospermia cohorts, suggesting heterogeneity in cancer risk by type of subfertility and within subfertility type. Describing families with similar cancer risk patterns provides a new avenue to increase homogeneity for focused gene discovery and environmental risk factor studies. Such discoveries will lead to more accurate risk predictions and improved counseling for patients and their families. STUDY FUNDING/COMPETING INTEREST(S) This work was funded by GEMS: Genomic approach to connecting Elevated germline Mutation rates with male infertility and Somatic health (Eunice Kennedy Shriver National Institute of Child Health and Human Development (NICHD): R01 HD106112). The authors have no conflicts of interest relevant to this work.

60 APPLIED LIFE SCIENCES↗

Base excision repair of the N -(2-deoxy- d - erythro -pentofuranosyl)-urea lesion by the hNEIL1 glycosylase

Abstract The N-(2-deoxy-d-erythro-pentofuranosyl)-urea DNA lesion forms following hydrolytic fragmentation of cis-5R,6S- and trans-5R,6R-dihydroxy-5,6-dihydrothymidine (thymine glycol, Tg) or from oxidation of 7,8-dihydro-8-oxo-deoxyguanosine (8-oxodG) and subsequent hydrolysis. It interconverts between α and β deoxyribose anomers. Synthetic oligodeoxynucleotides containing this adduct are efficiently incised by unedited (K242) and edited (R242) forms of the hNEIL1 glycosylase. The structure of a complex between the active site unedited mutant CΔ100 P2G hNEIL1 (K242) glycosylase and double-stranded (ds) DNA containing a urea lesion reveals a pre-cleavage intermediate, in which the Gly2 N-terminal amine forms a conjugate with the deoxyribose C1′ of the lesion, with the urea moiety remaining intact. This structure supports a proposed catalytic mechanism in which Glu3-mediated protonation of O4′ facilitates attack at deoxyribose C1′. The deoxyribose is in the ring-opened configuration with the O4′ oxygen protonated. The electron density of Lys242 suggests the ‘residue 242-in conformation’ associated with catalysis. This complex likely arises because the proton transfer steps involving Glu6 and Lys242 are hindered due to Glu6-mediated H-bonding with the Gly2 and the urea lesion. Consistent with crystallographic data, biochemical analyses show that the CΔ100 P2G hNEIL1 (K242) glycosylase exhibits a residual activity against urea-containing dsDNA.

59 BASIC BIOLOGICAL SCIENCES↗

ENVnet provides a global molecular resource of dissolved organic matter

Dissolved organic matter (DOM) is an important component of Earth's carbon cycle and one of the planet's most chemically diverse pools, yet the molecular structures of its constituents remain largely unresolved. This limitation has hindered our ability to link DOM composition to microbial processes and ecosystem function. Here we present ENVnet, a global molecular repository built from tandem mass spectrometry data collected across 13 terrestrial and aquatic environment types, including 419 newly generated samples that expand publicly available DOM metabolomics data and cover previously underrepresented environments. By computationally deconvolving chimeric mass spectra, a longstanding challenge in environmental metabolomics, we recover high-quality fragmentation data for >22,000 distinct molecular features (defined by a specific precursor mass and fragmentation pattern). Using ENVnet, we uncover conserved and environment-specific molecular patterns in DOM composition and underlying biogeochemical processes. We also use molecular features encoded in ENVnet to train predictive models of DOM persistence, allowing molecular-level assessment of microbial turnover in independent systems.

54 ENVIRONMENTAL SCIENCES↗

Machine learning in nuclear materials research

Nuclear materials are often demanded to function for extended time in extreme environments, including high radiation fluxes with associated transmutations, high temperature and temperature gradients, mechanical stresses, and corrosive coolants. They also have a wide range of microstructural and chemical makeups, resulting in multifaceted and often out-of-equilibrium interactions. Machine learning (ML) is increasingly being used to tackle these complex time-dependent interactions and aid researchers in developing models and making predictions, sometimes with better accuracy than traditional modeling that focuses on one or two parameters at a time. Conventional practices of acquiring new experimental data in nuclear materials research are often slow and expensive, limiting the opportunity for data-centric ML, but new methods are changing that paradigm. Here we review high-throughput computational and experimental data approaches, especially robotic experimentation and active learning that is based on Gaussian process and Bayesian optimization. We show ML examples in structural materials (e.g., reactor pressure vessel (RPV) alloys and radiation detecting scintillating materials) and highlight new techniques of high-throughput sample preparation and characterizations, and automated radiation/environmental exposures and real-time online diagnostics. Herein, this review suggests that ML models of material constitutive relations in plasticity, damage, and even electronic and optical responses to radiation are likely to become powerful tools as they develop. Finally, we speculate on how the recent trends of using natural language processing (NLP) to aid the collection and analysis of literature data, interpretable artificial intelligence (AI), and the use of streamlined scripting, database, workflow management, and cloud computing platforms that will soon make the utilization of ML techniques as commonplace as the spreadsheet curve-fitting practices of today.

36 MATERIALS SCIENCE↗

Leveraging structure-informed machine learning for fast steric zipper propensity prediction across whole proteomes

Predicting the amyloid fold and the propensity of peptide segments to adopt amyloid-like structures remain a challenge. However, recent progress has facilitated structure-based prediction of steric zipper propensity and the use of machine learning to accelerate the calculation of predictive models across many scientific areas. Leveraging these advances, we have developed a new approach for rapid proteome-wide assessment of zipper profiles that is informed by four million steric zipper predictions collected over ten years. This collection is used to build a machine learning model capable of rapidly predicting steric zipper propensity, and allowing for the assessment of zippers at both the protein and proteome level. Our predictions show enrichment for zipper forming segments in proteins involved in cell wall reorganization in yeast, highlighting a potential category of interest for experimental characterization. Overall, our predictive model allows for the exploration of amyloid formation across the tree of life and provides a tool for assessment of both novel and designed sequences for zipper density.

Biochemistry & Molecular Biology↗

Ultrastructure of the green alga Dunaliella salina strain CCAP19/18 (Chlorophyta) as investigated by quick-freeze deep-etch electron microscopy

The single-celled green alga Dunaliella salina is a model system for studies on stress biology, in particular regarding secondary carotenoid accumulation. Under non-stress conditions the cells are green, but under abiotic stress the cells turn orange, because they switch their metabolism and accumulate β-Carotene in globules in the chloroplast. For the first time, Quick-freeze deep-etch electron microscopy was used to visualize cellular structures in green and orange cells of D. salina strain CCAP19/18. This allowed us to present an in-depth analysis of the cellular ultrastructure describing and comparing the features of the two cell types. Our images illustrate the presence of a pericellular matrix for this strain of D. salina. The pericellular matrix was spongy and strands of unknown material anchored it into the plasma membrane. The cytoplasm contained a variety of vesicles, vacuoles, and acidocalcisomes. We could show for strain CCAP19/18 that cytoplasmic lipid bodies were often in close proximity to and sometimes in contact with the outer chloroplast envelope membrane and with the endoplasmic reticulum. Major visible differences between green and orange cells were in the chloroplast: the orange cells have greatly reduced amounts of thylakoid membranes and greatly increased numbers of β-Carotene globules. We showed that the β-Carotene globules often made point contacts with thylakoid membranes, and frequently laid side-by-side along the thylakoid membrane surface, providing support to studies that indicated exchange of molecules between β-Carotene globules and thylakoid membranes. A unique finding was the β-Carotene globule duplets, suggesting intermediate stages in β-Carotene globule morphogenesis. Overall, the β-Carotene globules appear to be similar to plant plastoglobuli. We provide this description of the cellular ultrastructure features as a resource in the context of the recent publication of the genome of D. salina strain CCAP19/18 to expand on the knowledge regarding this novel reference strain.

59 BASIC BIOLOGICAL SCIENCES↗

Quality-Control Mechanism for Telomerase RNA Folding in the Cell

Long non-coding RNAs can often fold into different conformations. Telomerase RNA, an essential component of the telomerase ribonucleoprotein (RNP) enzyme, must fold into a defined structure to fulfill its function with the protein catalytic subunit (TERT) and other accessory factors. However, the mechanism by which the correct folding of telomerase RNA is warranted in a cell is still unknown. Here we show that La-related protein Pof8 specifically recognizes the conserved pseudoknot region of telomerase RNA and instructs the binding of the Lsm2-8 complex to its mature 3' end, thus selectively protecting the correctly folded RNA from exonucleolytic degradation. In the absence of Pof8, TERT assembles with misfolded RNA and produces little telomerase activity. Therefore, Pof8 plays a key role in telomerase RNA folding quality control, ensuring that TERT only assembles with functional telomerase RNA to form active telomerase. Our finding reveals a mechanism for non-coding RNA folding quality control.

59 BASIC BIOLOGICAL SCIENCES↗