Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Constructive Tree”

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 55 records · Page 3

Classification of Complete Proteomes of Different Organisms and Protein Sets Based on Their Protein Distributions in Terms of Some Key Attributes of Proteins

The existence of complete genome sequences makes it important to develop different approaches for classification of large-scale data sets and to make extraction of biological insights easier. Here, we propose an approach for classification of complete proteomes/protein sets based on protein distributions on some basic attributes. We demonstrate the usefulness of this approach by determining protein distributions in terms of two attributes: protein lengths and protein intrinsic disorder contents (ID). The protein distributions based on L and ID are surveyed for representative proteome organisms and protein sets from the three domains of life. The two-dimensional maps (designated as fingerprints here) from the protein distribution densities in the LD space defined by ln( L ) and ID are then constructed. The fingerprints for different organisms and protein sets are found to be distinct with each other, and they can therefore be used for comparative studies. As a test case, phylogenetic trees have been constructed based on the protein distribution densities in the fingerprints of proteomes of organisms without performing any protein sequence comparison and alignments. The phylogenetic trees generated are biologically meaningful, demonstrating that the protein distributions in the LD space may serve as unique phylogenetic signals of the organisms at the proteome level.

59 BASIC BIOLOGICAL SCIENCES↗

Disturbance legacies and climate jointly drive tree growth and mortality in an intensively studied boreal forest

How will regional growth and mortality change with even relatively small climate shifts, even independent of catastrophic disturbances? This question is particularly acute for the North American boreal forest, which is carbon-dense and subject The goals of this study were to combine dendrochronological sampling, inventory records, and machine-learning algorithms to understand how tree growth and death have changed at one highly studied site (Northern Old Black Spruce, NOBS) in the central Canadian boreal forest. Over the 1999-2012 inventory period, mean DBH increased even as stand density and basal area declined significantly from 41.3 to 37.5 m2 ha-1. Tree mortality averaged 1.4±0.6% yr-1, with most mortality occurring in medium-sized trees. A combined tree ring chronology constructed from 2001, 2004, and 2012 sampling showed several periods of extreme growth depression, with increased mortality lagging depressed growth by ~5 years. Minimum and maximum air temperatures exerted a negative influence on tree growth, while precipitation and climate moisture index had a positive effect; both current- and previous-year data exerted significant effects. Models based on these variables explained 23-44% of the ring-width variability. There have been at least one, and probably two, significant recruitment episodes since stand initiation, and we infer that past climate extremes led to significant NOBS mortality still visible in the current forest structure. These results imply that a combination of successional and demographic processes, along with mortality driven by abiotic factors, continue to affect the stand, with significant implications for our understanding of previous work at NOBS and the sustainable management of regional forests.

black spruce↗

Poplar: a phylogenomics pipeline

Motivation Generating phylogenomic trees from the genomic data is essential in understanding biological systems. Each step of this complex process has received extensive attention and has been significantly streamlined over the years. Given the public availability of data, obtaining genomes for a wide selection of species is straightforward. However, analyzing that data to generate a phylogenomic tree is a multistep process with legitimate scientific and technical challenges, often requiring a significant input from a domain-area scientist. Results We present Poplar, a new, streamlined computational pipeline, to address the computational logistical issues that arise when constructing the phylogenomic trees. It provides a framework that runs state-of-the-art software for essential steps in the phylogenomic pipeline, beginning from a genome with or without an annotation, and resulting in a species tree. Running Poplar requires no external databases. In the execution, it enables parallelism for execution for clusters and cloud computing. The trees generated by Poplar match closely with state-of-the-art published trees. The usage and performance of Poplar is far simpler and quicker than manually running a phylogenomic pipeline. Availability and implementation Freely available on GitHub at https://github.com/sandialabs/poplar. Implemented using Python and supported on Linux.

Koning, Elizabeth [Sandia National Laboratories (S↗

Tensor networks, p -adic fields, and algebraic curves: arithmetic and the AdS 3 /CFT 2 correspondence

One of the many remarkable properties of conformal field theory in two dimensions is its connection to algebraic geometry. Since every compact Riemann surface is a projective algebraic curve, many constructions of interest in physics (which a priori depend on the analytic structure of the spacetime) can be formulated in purely algebraic language. This opens the door to interesting generalizations, obtained by taking another choice of field: for instance, the p-adics. Herein, we generalize the $\mathrm{AdS} / \mathrm{CFT}$ correspondence according to this principle; the result is a formulation of holography in which the bulk geometry is discrete—the Bruhat–Tits tree for $\mathrm{PGL} (2, \mathbb{Q}_p)$—but the group of bulk isometries nonetheless agrees with that of boundary conformal transformations and is not broken by discretization. We anticipate that this forms the natural geometric setting for tensor networks that have been proposed as models of bulk reconstruction via quantum error correcting codes; in certain cases, geodesics in the Bruhat–Tits tree reproduce those constructed using quantum error correction. Other aspects of holography also hold: Standard holographic results for massive free scalar fields in a fixed background carry over to the tree, whose vertical direction can be interpreted as a renormalization-group scale for modes in the boundary $\mathrm{CFT}$. Higher-genus bulk geometries (the $\mathrm{BTZ}$ black hole and its generalizations) can be understood straightforwardly in our setting, and the Ryu–Takayanagi formula for the entanglement entropy appears naturally

97 MATHEMATICS AND COMPUTING↗

Host-region parameters for an adjustable model for crustal earthquakes to facilitate the implementation of the backbone approach to building ground-motion logic trees in probabilistic seismic hazard analysis

The backbone approach to constructing a ground-motion logic tree for probabilistic seismic hazard analysis (PSHA) can address shortcomings in the traditional approach of populating the branches with multiple existing, or potentially modified, ground-motion models (GMMs) by rendering more transparent the relationship between branch weights and the resulting distribution of predicted accelerations. To capture epistemic uncertainty in a tractable manner, there are benefits in building the logic tree through the application of successive adjustments for differences in source, path, and site characteristics between the host region of the selected backbone GMM and the target region for which the PSHA is being conducted. The implementation of this approach is facilitated by selecting a backbone GMM that is amenable to such host-to-target adjustments for individual source, path, and site characteristics. The NGA-West2 GMM of Chiou and Youngs (CY14) has been identified as a highly adaptable model for crustal seismicity that is well suited to such adjustments. Rather than using generic source, path, and site characteristics assumed appropriate for the host region, the final suite of adjusted GMMs for the target region will be better constrained if the host-region parameters are defined specifically on the basis of their compatibility with the CY14 backbone GMM. To this end, making use of a recently developed crustal shear-wave velocity profile consistent with CY14, we present an inversion of the model to estimate the key source and path parameters, namely the stress parameter and the anelastic attenuation. With these outputs, the effort in constructing a ground-motion logic tree for any PSHA dealing with crustal seismicity can be focused primarily on the estimation of the target-region characteristics and their associated uncertainties. The inversion procedure can also be adapted for any application in which different constraints might be relevant.

Stafford, Peter J.↗

Optimal decision trees for categorical data via integer programming

Decision trees have been a very popular class of predictive models for decades due to their interpretability and good performance on categorical features. However, they are not always robust and tend to overfit the data. Additionally, if allowed to grow large, they lose interpretability. In this paper, we present a mixed integer programming formulation to construct optimal decision trees of a prespecified size. We take the special structure of categorical features into account and allow combinatorial decisions (based on subsets of values of features) at each node. Our approach can also handle numerical features via thresholding. Here we show that very good accuracy can be achieved with small trees using moderately-sized training sets. The optimization problems we solve are tractable with modern solvers.

97 MATHEMATICS AND COMPUTING↗

Regional surname affinity: A spatial network approach

Abstract Objective We investigate surname affinities among areas of modern‐day China, by constructing a spatial network, and making community detection. It reports a geographical genealogy of the Chinese population that is result of population origins, historical migrations, and societal evolutions. Materials and methods We acquire data from the census records supplied by China's National Citizen Identity Information System, including the surname and regional information of 1.28 billion registered Chinese citizens. We propose a multilayer minimum spanning tree (MMST) to construct a spatial network based on the matrix of isonymic distances, which is often used to characterize the dissimilarity of surname structure among areas. We use the fast unfolding algorithm to detect network communities. Results We obtain a 10‐layer MMST network of 362 prefecture nodes and 3,610 edges derived from the matrix of the Euclidean distances among these areas. These prefectures are divided into eight groups in the spatial network via community detection. We measure the partition by comparing the inter‐distances and intra‐distances of the communities and obtain meaningful regional ethnicity classification. Discussion The visualization of the resulting communities on the map indicates that the prefectures in the same community are usually geographically adjacent. The formation of this partition is influenced by geographical factors, historic migrations, trade and economic factors, as well as isolation of culture and language. The MMST algorithm proves to be effective in geo‐genealogy and ethnicity classification for it retains essential information about surname affinity and highlights the geographical consanguinity of the population.

59 BASIC BIOLOGICAL SCIENCES↗

Simulation Model of Microsphere Distribution for Selective Internal Radiation Therapy Agrees With Observations

Purpose: To perform a detailed analysis of microsphere distribution in biopsy material from a patient treated with {sup 90}Y-labeled resin spheres and characterize microsphere distribution in the hepatic artery tree, and to construct a novel dichotomous bifurcation model for microsphere deposits and evaluate its accuracy in simulating the observed microsphere deposits. Methods and Materials: Our virtual model consisted of arteries that successively branched into 2 new generations of arteries at 20 nodes. The artery diameter exponentially decreased from the lowest generation to the highest generation. Three variable parameters were optimized to obtain concordance between simulations and measure microsphere distributions: an artery coefficient of variation (ACV) for the diameter of all artery generations and the microsphere flow distribution at the nodes; a hepatic tree distribution volume (HDV) for the artery tree; and an artery diameter reduction (ADR) parameter. The model was tested against previously measured activity concentrations in 84 biopsies from the liver of 1 patient. In 16 of 84 biopsies, the microsphere distribution regarding cluster size and localization in the artery tree was determined via light microscopy of 30-μm sections (mean concentration, 14 microspheres/mg; distributions divided into 3 groups with mean microsphere concentrations of 4.6, 14, and 28 microspheres/mg). Results: Single spheres and small clusters were observed in terminal arterioles, whereas large clusters, up to 450 microspheres, were observed in larger arterioles. For 14 microspheres/mg, the optimized parameter values were ACV=0.35, HDV = 50 cm{sup 3}, and ADR=6 μm. For 4.6 microspheres/mg, ACV and ADR decreased to 0.26 and 0 μm, respectively, whereas HDV increased to 130 cm{sup 3}. The opposite trend was observed for 28 microspheres/mg: ACV = 0.49, HDV = 20 cm{sup 3}, and ADR = 8 μm. Conclusion: Simulations and measurements reveal that microsphere clusters are larger and more common in volumes with high microsphere concentrations and indicate that the spatial distribution of the artery tree must be considered in estimates of microsphere distributions.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

B-Trees and Cache-Oblivious B-Trees with Different-Sized Atomic Keys

Most B-tree articles assume that all $\textit{N}$ keys have the same size $\textit{K}$, that $\textit{f = B/K}$ keys fit in a disk block, and therefore that the search cost is $O(\text{log}_{f+1}N) $ block transfers. When keys have variable size, B-tree operations have no nontrivial performance guarantees, however. This article provides B-tree-like performance guarantees on dictionaries that contain keys of different sizes in a model in which keys must be stored and compared as opaque objects. The resulting atomic-key dictionaries exhibit performance bounds in terms of the average key size and match the bounds when all keys are the same size. Atomic-key dictionaries can be built with minimal modification to the B-tree structure, simply by choosing the pivot keys properly. Further, this article describes both static and dynamic atomic-key dictionaries. In the static case, if there are $\textit{N}$ keys with average size $\textit{K}$, the search cost is $O(\lceil K/B\rceil \text{log}_{1 + \lceil K/B\rceil} N)$ expected transfers. It is not possible to transform these expected bounds into worst-case bounds. The cost to build the tree is $\textit{O(NK)}$ operations and $\textit{O(NK/B)}$ transfers if all keys are presented in sorted order. If not, the cost is the sorting cost. For the dynamic dictionaries, the amortized cost to insert a key $\textit{κ}$ of arbitrary length at an arbitrary rank is dominated by the cost to search for $\textit{κ}$. Specifically, the amortized cost to insert a key $\textit{κ}$ of arbitrary length and random rank is $O(\lceil K/B\rceil log_{1+ \lceil B/K\rceil } N + |κ|/B)$ transfers. A dynamic-programming algorithm is shown for constructing a search tree with minimal expected cost. This article also gives a cache-oblivious static atomic-key B-tree, which achieves the same asymptotic performance as the static B-tree dictionary, mentioned previously. A cache-oblivious data structure or algorithm is not parameterized by the block size $\textit{B}$ or memory size $\textit{M}$ in the memory hierarchy; rather, it is universal, working simultaneously for all possible values of $\textit{B}$ or $\textit{M}$. On a machine with block size $\textit{B}$, if there are $\textit{N}$ keys with average size $\textit{K}$, search operations costs $O(\lceil K/B \rceil log_{1+ \lceil B/K \rceil} N)$ block transfers in expectation. This cache-oblivious layout can be built in $\textit{O(N}\text{log}\textit{(NK))}$ processor operations.

97 MATHEMATICS AND COMPUTING↗

Phylogenomics and diversification of Sordariomycetes

The Sordariomycetes is a specious, morphologically diverse, and widely distributed class of the phylum Ascomycota that forms a well-supported clade diverged from Leotiomycetes. Aside from their ecological significance as plant and human pathogens, saprobes, endophytes, and fungicolous taxa, species of Sordariomycetes produces a wide range of chemically novel and diverse metabolites used in important fields. Recent phylogenetic analyses derived from a small number of genes have considerably increased our understanding of the family, order, and subclass relationships within Sordariomycetes, but several important groups have not been resolved well. In addition, there are various paraphyletic or polyphyletic groups. Moreover, the criteria used to establish higher ranks remain highly variable across different studies. Therefore, the taxonomy of Sordariomycetes is in constant flux, remains poorly understood, and is subject to much controversy. Here, for the first time, we have assembled a phylogenetic dataset containing 638 genomes representing the 156 genera, 50 families, and 17 orders and 5 subclasses of Sordariomycetes. This data set is based on 1124 genes and results in a well-resolved phylogenomic tree. We further constructed an evolutionary timeline of Sordariomycetes diversification based on the genomic data sets. Our divergence time estimate results are inconsistent with previous studies, suggesting estimates of node ages are less precise and varied. Based on these results, we discuss the higher ranks of Sordariomycetes and empirically propose an unprecedented taxonomic framework for the class.

59 BASIC BIOLOGICAL SCIENCES↗

When can we detect lianas from space? Toward a mechanistic understanding of liana‐infested forest optics

Abstract Lianas, woody vines acting as structural parasites of trees, have profound effects on the composition and structure of tropical forests, impacting tree growth, mortality, and forest succession. Remote sensing could offer a powerful tool for quantifying the scale of liana infestation, provided the availability of robust detection methods. We analyze the consistency and global geographic specificity of spectral signals—reflectance across wavelengths—from liana‐infested tree crowns and forest stands, examining the underlying mechanisms of these signals. We compiled a uniquely comprehensive database, including leaf reflectance spectra from 5424 leaves, fine‐scale airborne reflectance data from 999 liana‐infested canopies, and coarse‐scale satellite reflectance data covering 775 ha of liana‐infested forest stands. To unravel the mechanisms of the liana spectral signal, we applied mechanistic radiative transfer models across scales, establishing a synthesis of the relative importance of different mechanisms, which we corroborate with field data on liana leaf chemistry and canopy structure. We find a consistent liana spectral signal at canopy and stand scales across globally distributed sites. This signature mainly arises at the canopy level due to direct effects of more horizontal leaf angles, resulting in a larger projected leaf area, and indirect effects from increased light scattering in the near and short‐wave infrared regions, linked to lianas' less costly leaf construction compared with trees on average. The existence of a consistent global spectral signal for lianas suggests that large‐scale quantification of liana infestation is feasible. However, because the traits responsible for the liana canopy‐reflectance signal are not exclusive to lianas, accurate large‐scale detection requires rigorously validated remote sensing methods. Our models highlight challenges in automated detection, such as potential misidentification due to leaf phenology, tree life history, topography, and climate, especially where the scale of liana infestation is less than a single remote sensing pixel. The observed cross‐site patterns also prompt ecological questions about lianas' adaptive similarities in optical traits across environments, indicating possible convergent evolution due to shared constraints on leaf biochemical and structural traits.

Environmental Sciences & Ecology↗

Challenges in leveraging existing human performance data for quantifying the IDHEAS HRA method

Our article documents an exploratory study for collecting and using human performance data to inform human error probability (HEP) estimates for a new human reliability analysis (HRA) method, the IntegrateD Human Event Analysis System (IDHEAS). The method was based on cognitive models and mechanisms underlying human behaviour and employs a framework of 14 crew failure modes (CFMs) to represent human failures typical for human performance in nuclear power plant (NPP) internal, at-power events [1]. A decision tree (DT) was constructed for each CFM to assess the probability of the CFM occurring in different contexts. Data needs for IDHEAS quantification are discussed. Then, the data collection framework and process is described and how the collected data were used to inform HEP estimation is illustrated with two examples. Next, five major technical challenges are identified for leveraging human performance data for IDHEAS quantification. Furthermore, these challenges reflect the data needs specific to IDHEAS. More importantly, they also represent the general issues with current human performance data and can provide insight for a path forward to support HRA data collection, use, and exchange for HRA method development, implementation, and validation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Phylogenetic relationships and phylogeography of relevant lineages within the complex Campanulaceae family in Macaronesia

Abstract Macaronesia has long been recognized as a natural model for studying evolutionary processes in plant diversification. Several studies have attempted to focus on single lineages, and few have covered the diversification of a family across all the archipelagos. We used a comprehensive sample to clarify the phylogenetic relationships and the biogeographic history of the Macaronesian Campanulaceae. Hypotheses related to the colonization of these archipelagos will be used to examine the diversification patterns of different lineages. We sequenced the ITS region and six cp DNA markers ( atpB , matK , petD , rbcL , trnL‐F , and psbA‐trnH ) from 10 Campanulaceae species, including seven endemic species in Macaronesia. The phylogeny of these taxa was reconstructed using maximum parsimony, maximum likelihood, and Bayesian inference. To study the relationships within each lineage, haplotype networks were calculated using NeighborNet and TCS algorithms. Moreover, data were combined with fossil information to construct time‐calibrated trees for the Macaronesian Campanulaceae species. The phylogenetic analyses are largely congruent with current taxon circumscriptions, and all the endemic genera formed monophyletic clades, namely Azorina in Azores; Musschia in Madeira; and Campanula in Cape Verde. The Azorina clade and the Cape Verde endemic Campanula may share a common ancestor in North Africa, and the divergence was dated ca. 12.3 million years ago (Mya). The divergence of the Musschia clade began in the Pliocene ca. 3.4 Mya. Moreover, several examples of intraspecific variation were revealed among the native species with a clear geographic structured patterns, suggesting that cryptic diversity might exist within the native Macaronesian Campanulaceae when compared to the close mainland taxa (e.g., Campanula erinus , Trachelium caeruleum ), but additional studies are needed to support the molecular data. This study highlights the power of combining data (e.g., phylogeny and divergence times, with species distribution data) for testing diversification hypotheses within the unique Macaronesian flora, providing useful information for future conservation efforts.

Menezes, Tiago↗

Gauged soft recursion: on-shell construction of Goldstone-gauge amplitudes

We present a new on-shell recursion relation for scattering amplitudes involving Nambu-Goldstone bosons with a gauged unbroken symmetry. A central challenge is that gauge interactions break Adler’s zero condition for charged scalars, invalidating the standard soft recursion. To overcome this, we introduce a “gauged soft recursion” that leverages the soft theorems of the gauge bosons themselves, combined with a novel decomposition of amplitudes into gauge-invariant components where Adler’s zero is partially restored. The formalism, which also incorporates internal gauge bosons via angular momentum constraints, enables the systematic construction of tree-level amplitudes with arbitrary numbers of Goldstone bosons and gauge bosons in both Abelian and non-Abelian theories, as we demonstrate with explicit examples.

Chiral Lagrangian↗

Epidural Analgesia Versus Patient-Controlled Analgesia for Pain Relief in Uterine Artery Embolization for Uterine Fibroids: A Decision Analysis

Purpose: This study was designed to compare the costs and effects of epidural analgesia (EDA) to those of patient-controlled intravenous analgesia (PCA) for postintervention pain relief in women having uterine artery embolization (UAE) for systematic uterine fibroids. Methods: Cost-effectiveness analysis (CEA) based on data from the literature by constructing a decision tree to model the clinical pathways for estimating the effects and costs of treatment with EDA and PCA. Literature on EDA for pain-relief after UAE was missing, and therefore, data on EDA for abdominal surgery were used. Outcome measures were compared costs to reduce one point in visual analogue score (VAS) or numeric rating scale (NRS) for pain 6 and 24 h after UAE and risk for complications. Results: Six hours after the intervention, the VAS was 3.56 when using PCA and 2.0 when using EDA. The costs for pain relief in women undergoing UAE with PCA and EDA were Euro-Sign 191 and Euro-Sign 355, respectively. The costs for EDA to reduce the VAS score 6 h after the intervention with one point compared with PCA were Euro-Sign 105 and Euro-Sign 179 after 24 h. The risk of having a complication was 2.45 times higher when using EDA. Conclusions: The results of this indirect comparison of EDA for abdominal surgery with PCA for UAE show that EDA would provide superior analgesia for post UAE pain at 6 and 24 h but with higher costs and an increased risk of complications.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Comparative genomics and transcriptomics analysis reveals evolution patterns of selection in the Salix phylogeny

Background: Willows are widely distributed in the northern hemisphere and have good adaptability to different living environment. The increasing of genome and transcriptome data provides a chance for comparative analysis to study the evolution patterns with the different origin and geographical distributions in the Salix phylogeny. Results: Transcript sequences of 10 Salicaceae species were downloaded from public databases. All pairwise of orthologues were identified by comparative analysis in these species, from which we constructed a phylogenetic tree and estimated the rate of diverse. Divergence times were estimated in the 10 Salicaceae using comparative transcriptomic analysis. All of the fast-evolving positive selection sequences were identified, and some cold-, drought-, light-, universal-, and heat- resistance genes were discovered. Conclusions: The divergence time of subgenus Vetrix and Salix was about 17.6–16.0 Mya during the period of Middle Miocene Climate Transition (21–14 Mya). Subgenus Vetrix diverged to migratory and resident groups when the climate changed to the cool and dry trend by 14 Mya. Cold- and light- stress genes were involved in positive selection among the resident Vetrix, and which would help them to adapt the cooling stage. Universal- stress genes exhibited positive selection among the migratory group and subgenus Salix. These data are useful for comprehending the adaptive evolution and speciation in the Salix lineage.

59 BASIC BIOLOGICAL SCIENCES↗

Systematic Identification of Gene Families for Use as “Markers” for Phylogenetic and Phylogeny-Driven Ecological Studies of Bacteria and Archaea and Their Major Subgroups

With the astonishing rate that genomic and metagenomic sequence data sets are accumulating, there are many reasons to constrain the data analyses. One approach to such constrained analyses is to focus on select subsets of gene families that are particularly well suited for the tasks at hand. Such gene families have generally been referred to as “marker” genes. We are particularly interested in identifying and using such marker genes for phylogenetic and phylogeny-driven ecological studies of microbes and their communities (e.g., construction of species trees, phylogenetic based assignment of metagenomic sequence reads to taxonomic groups, phylogeny-based assessment of alpha- and beta-diversity of microbial communities from metagenomic data). We therefore refer to these as PhyEco (for phylogenetic and phylogenetic ecology) markers. The dual use of these PhyEco markers means that we needed to develop and apply a set of somewhat novel criteria for identification of the best candidates for such markers. The criteria we focused on included universality across the taxa of interest, ability to be used to produce robust phylogenetic trees that reflect as much as possible the evolution of the species from which the genes come, and low variation in copy number across taxa. We describe here an automated protocol for identifying potential PhyEco markers from a set of complete genome sequences. The protocol combines rapid searching, clustering and phylogenetic tree building algorithms to generate protein families that meet the criteria listed above. We report here the identification of PhyEco markers for different taxonomic levels including 40 for “all bacteria and archaea”, 114 for “all bacteria (greatly expanding on the ~30 commonly used), and 100 s to 1000 s for some of the individual phyla of bacteria. This new list of PhyEco markers should allow much more detailed automated phylogenetic and phylogenetic ecology analyses of these groups than possible previously.

59 BASIC BIOLOGICAL SCIENCES↗

Genome collection processing for “Conserved upper thermal limits and small safety margins in soil copiotrophic bacteria”

We extracted the genomic DNA of 400 randomly selected isolates using a Quick-DNA Microprep Kit (Zymo Research D3020) according to the manufacturer’s protocol. We then submitted the extracted gDNA samples for short-read Illumina sequencing (200 Mbp) at SeqCoast Genomics (Portsmouth, NH, USA). After preprocessing the sequences using Trimmommatic (Bolger et al. 2014), we assembled the genomes using SPADES (Bankevich et al. 2012) and checked the quality of each assembly using QUAST (Gurevich et al. 2013). We processed the genome assemblies using a KBase (v1.4.0) pipeline (Allen et al. 2017; Arkin et al. 2018). Briefly, we used DRAM (v0.1.2) with default settings to annotate the genome assemblies. We then evaluated genome quality and possible contamination levels using CheckM (v1.0.18) (Parks et al. 2015) and retained genomes with completeness above 98% and contamination below 5% (n = 354), following the authors' guidelines. We then obtained taxonomic assignments for all remaining isolates using the Genome Taxonomy Database tool GTDB-Tk (v2.3.2, database version r214) (Chaumeil et al. 2019). We constructed a phylogenetic tree using the tool SpeciesTree (v2.2.0). We then trimmed the tree (using Trim SpeciesTree to GenomeSet- v1.4.0), retaining only tips within our collection with measured thermal performance.

59 BASIC BIOLOGICAL SCIENCES↗