Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “software differences”

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 1,603 records · Page 89
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

SU-E-I-31: Differences Observed in Radiation Doses Across 2 Similar CT Scanners From Adult Brain-Neck CT Angiography

Purpose: The aim of this study is to evaluate the difference in radiation doses from adult Brain-Neck CT angiography (CTA) between two CT scanners. Methods: We collected CT dose index data (CTDIvol, DLP) from adult Brain-Neck CTA performed with two CT scanners (Sensation 64 (S64) and Definition AS (AS), Siemens Healthcare) performed at two of our facilities from Jan 1st to Dec 31th, 2014. X-ray dose management software (Radmetrics, Bayer Healthcare) was used to mine these data. All exams were performed with Tube Current Modulation (Care Dose 4D), tube voltage of 120 kVp, quality reference mAs of 300, beam collimation of 64*0.6 mm. The rotation time was set to 0.5 sec for S64 and 1.0 sec for AS. We also scanned an anthropomorphic skull and chest phantom under routine Brain-Neck CTA protocol with the two scanners and extracted the tube current values from the raw projection data. Results: The mean CTDIvol and DLP in Brain-Neck CTA was 72 mGy and 2554 mGy*cm for AS, which was substantially larger than the mean values of 46 mGy and 1699 mGy*cm for S64. The maximum tube current was 583 mA for most cases on the S64 while the maximum was 666 mA for AS even though the rotation time set for AS was 1.0 sec. Measurements obtained with the anthropomorphic phantom showed that the tube current reached 583 mA at the shoulder region for S64 while it reached to 666 mA for AS. Conclusion: The results of this study showed that substantially different CT doses can Result from Brain-Neck CTA protocols even when similar scanners and similar settings are used. Though both scanners have a similar maximum mA rating, differences in mA were observed through the shoulders, resulting in substantially different CTDIvol values.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-E-T-01: 2-D Characterization of DLG Among All MLC Leaf Pairs

Purpose: The aim of this study is to evaluate the variation of dosimetric leaf-gap (DLG) along the travel path of each MLC leaf pair. This study evaluates whether the spatial variations in DLG could cause dose differences between TPS-calculated and measured dose. Methods: The 6MV DLG values were measured for all leaf pairs in the direction of leaf motion using a 2-D diode array and 0.6cc ion chamber. These measurements were performed on two Varian Linacs, employing the Millennium 120-leaf MLC and a 2-D-DLG variation map was created via in-house software. Several test plans were created with sweeping MLC fields using constant gaps from 2mm to 10mm and corrected for 2-D variation utilizing in-house software. Measurements were performed utilizing the MapCHECK at 5.0cm depth for plans with and without the 2-D DLG correction and compared to the TPS calculated dose via gamma analysis (3%/3mm). Results: The measured DLGs for the middle 40 MLC leaf pairs (0.5cm width) were very similar along the central superior-inferior axis, with maximum variation of 0.2mm. The outer 20 MLC leaf pairs (1.0cm width) have DLG values from 0.32mm (mean) to 0.65mm (maximum) lower than the central leaf-pair, depending on off-axis distance. Gamma pass rates for the 2mm, 4mm, and 6mm sweep plans increased by 23.2%, 28.7%, and 26.0% respectively using the 2-D-DLG correction. The most improved dose points occur in areas modulated by the 1.0cm leaf-pairs. The gamma pass rate for the 10mm sweep plan increased by only 7.7%, indicating that the 2D variation becomes less significant for dynamic plans with larger MLC gaps. Conclusion: Fluences residing significantly off-axis with narrow sweeping gaps may exhibit significant variations from planned dose due to large differences between the true DLG exhibited by the 1.0cm leaf-pairs versus the constant DLG value utilized by the TPS for dose calculation.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

GEO2D - Two-Dimensional Computer Model of a Ground Source Heat Pump System

This file contains a zipped file that contains many files required to run GEO2D. GEO2D is a computer code for simulating ground source heat pump (GSHP) systems in two-dimensions. GEO2D performs a detailed finite difference simulation of the heat transfer occurring within the working fluid, the tube wall, the grout, and the ground. Both horizontal and vertical wells can be simulated with this program, but it should be noted that the vertical wall is modeled as a single tube. This program also models the heat pump in conjunction with the heat transfer occurring. GEO2D simulates the heat pump and ground loop as a system. Many results are produced by GEO2D as a function of time and position, such as heat transfer rates, temperatures and heat pump performance. On top of this information from an economic comparison between the geothermal system simulated and a comparable air heat pump systems or a comparable gas, oil or propane heating systems with a vapor compression air conditioner. The version of GEO2D in the attached file has been coupled to the DOE heating and cooling load software called ENERGYPLUS. This is a great convenience for the user because heating and cooling loads are an input to GEO2D. GEO2D is a user friendly program that uses a graphical user interface for inputs and outputs. These make entering data simple and they produce many plotted results that are easy to understand. In order to run GEO2D access to MATLAB is required. If this program is not available on your computer you can download the program MCRInstaller.exe, the 64 bit version, from the MATLAB website or from the Geothermal Data Repository. This is a free download which will enable you to run GEO2D.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Evaluation of the EarthSHAB Stratospheric Solar Hot Air Balloon Flight Prediction Model Using Balloon Trajectory Data

Abstract The heliotrope is a solar balloon design which is constructed out of painter’s plastic, and the exterior is coated in charcoal powder. Darkening the plastic gives the balloon a high solar absorptance, which allows it to ascend into the lower stratosphere and float for hours at a time. The balloons have previously been used to lift scientific instruments into the stratosphere to study chemical explosions, earthquakes, and stratospheric aerosols. They have also been proposed as a platform for planetary exploration. Flight predictions are crucial to preflight planning to reduce safety risks and meet flight objectives. However, there exists a wide range of possible flight paths due to varying environmental conditions and solar balloon configurations. EarthSHAB is one such software that was designed to support flight planning using the weather forecasts and balloon properties to predict the flight path of a solar balloon. We compare EarthSHAB-simulated flight paths to a set of observed flight paths for the 3.5-m diameter heliotrope design called the “Cloudskimmer.” Using the criteria that the modeled paths must fall within 5% of the observations to be considered successful, we found that EarthSHAB successfully predicted the Cloudskimmer ascent rate and average float altitude 10% and 90% of the time, respectively. We also found that the average difference in the observed and predicted landing locations was 97 km and landing times were 54 ± 38 min. Significant deviations between the observed and predicted ascent rates and excursions at float were found to be associated with heavy payloads and convective cloud development, respectively. Significance Statement Solar balloons are used to lift scientific instruments into the lower stratosphere for hours at a time to study chemical explosions, earthquakes, stratospheric aerosols, and more. The flight paths of solar balloons can be difficult to predict due to variability in their design and surrounding environment. We evaluate the accuracy of EarthSHAB, a software that predicts the altitude profile and horizontal trajectory of a solar balloon using inputs such as the weather, balloon size, and balloon mass. Our results suggest that for the balloon design used in this study, EarthSHAB is best suited for modeling the behavior of balloons with lightweight payloads that do not fly within or directly above clouds.

Lien, Jessica M. [Sandia National Laboratories, Al↗
DOE OSTI2019

Comparison of two software systems for quantification of myocardial blood flow in patients with hypertrophic cardiomyopathy

Backgorund: Quantification of myocardial blood flow (MBF) by positron emission tomography (PET) is important for investigation of angina in hypertrophic cardiomyopathy (HCM). Several software programs exist for MBF quantification, but they have been mostly evaluated in patients (with normal cardiac geometry), referred for evaluation of coronary artery disease (CAD). Software performance has not been evaluated in HCM patients who frequently have hyperdynamic LV function, LV outflow tract (LVOT) obstruction, small LV cavity size, and variation in the degree/location of LV hypertrophy. Aim: We compared results of MBF obtained using PMod, which permits manual segmentation, to those obtained by FDA-approved QPET software which has an automated segmentation algorithm. Methods: {sup 13}N-ammonia PET perfusion data were acquired in list mode at rest and during pharmacologic vasodilation, in 76 HCM patients and 10 non-HCM patients referred for evaluation of CAD (CAD group.) Data were resampled to create static, ECG-gated and 36-frame-dynamic images. Myocardial flow reserve (MFR) and MBF (in ml/min/g) were calculated using QPET and PMod softwares. Results: All HCM patients had asymmetric septal hypertrophy, and 50% had evidence of LVOT obstruction, whereas non-HCM patients (CAD group) had normal wall thickness and ejection fraction. PMod yielded significantly higher values for global and regional stress-MBF and MFR than for QPET in HCM. Reasonably fair correlation was observed for global rest-MBF, stress-MBF, and MFR using these two softwares (rest-MBF: r = 0.78; stress-MBF: r = 0.66.; MFR: r = 0.7) in HCM patients. Agreement between global MBF and MFR values improved when HCM patients with high spillover fractions (> 0.65) were excluded from the analysis (rest-MBF: r = 0.84; stress-MBF: r = 0.72; MFR: r = 0.8.) Regionally, the highest agreement between PMod and QPET was observed in the LAD territory (rest-MBF: r = 0.82, Stress-MBF: r = 0.68) where spillover fraction was the lowest. Unlike HCM patients, the non-HCM patients (CAD group) demonstrated excellent agreement in MBF/MFR values, obtained by the two softwares, when patients with high spillover fractions were excluded (rest-MBF: r = 0.95; stress-MBF: r = 0.92; MFR: r = 0.95). Conclusions: Anatomic characteristics specific to HCM hearts contribute to lower correlations between MBF/MFR values obtained by PMod and QPET, compared with non-HCM patients. These differences indicate that PMod and QPET cannot be used interchangeably for MBF/MFR analyses in HCM patients.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2017

Heterologous Gene Expression of N -Terminally Truncated Variants of LipPks1 Suggests a Functionally Critical Structural Motif in the N -terminus of Modular Polyketide Synthase

Streptomyces genomes have a high G + C content and typically use an ATG or GTG codon to initiate protein synthesis. Although gene-finding tools perform well in low GC genomes, it is known that the accuracy in predicting a translational start site (TSS) is much less for high GC genomes. LipPks1 is a Streptomyces-derived, well-characterized modular polyketide synthase (PKS). Here, using this enzyme as a model, we experimentally investigated the effects of alternative TSSs using a heterologous host, Streptomyces venezuelae. One of the TSSs employed boosted the protein level by 59-fold and the product yield by 23-fold compared to the originally annotated start codon. Interestingly, a structural model of the PKS indicated the presence of a structural motif in the N-terminus, which may explain the observed different protein levels together with a proline and arginine-rich sequence that may inhibit translational initiation. Lastly, this structure was also found in six other modular PKSs that utilize noncarboxylated starter substrates, which may guide the selection of optimal TSSs in conjunction with start-codon prediction software.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

Analysis of rupture area of aftershocks caused by twin earthquakes (Case study: 11 April 2012 earthquakes of Aceh-North Sumatra)

The 11 April 2012 earthquakes off-shore Aceh-North Sumatra are unique events for the history of Indonesian earthquake. It is unique because that they have similar magnitude, 8.5 Mw and 8.1 Mw; close to epicenter distance, similar strike-slip focal mechanism, and occuring in outer rise area. The purposes of this research are: (1) comparing area of earthquakes base on models and that of calculation, (2) fitting the shape and the area of earthquake rupture zones, (3) analyzing the relationship between rupture area and magnitude of the earthquakes. Rupture area of the earthquake fault are determined by using 4 different formulas, i.e. Utsu and Seki (1954), Wells and Coppersmith (1994), Ellsworth (2003), and Christophersen and Smith (2000). The earthquakes aftershock parameters are taken from PGN (PusatGempabumiNasional or National Earthquake Information Center) of BMKG (Indonesia Agency Meteorology Climatology and Geophysics). The aftershock epicenters are plotted by GMT’s software. After that, ellipse and rectangular models of aftershock spreading are made. The results show that: (1) rupture areas were calculated using magnitude relationship which are larger than the the aftershock distributions model, (2) the best fitting model for that earthquake aftershock distribution is rectangular associated with Utsu and Seki (1954) formula, (3) the larger the magnitude of the earthquake, the larger area of the fault.

58 GEOSCIENCES↗
DOE OSTI2023

Color Transformations of Photometric Measurements of Galaxies in Optical and Near-infrared Wide-field Imaging Surveys*

Abstract Over the past 2 decades, wide-field photometric surveys in optical and infrared domains reached a nearly all-sky coverage thanks to numerous observational facilities operating in both hemispheres. However, subtle differences among exact realizations of Johnson and SDSS photometric systems require one to convert photometric measurements into the same system prior to analysis of composite data sets originating from multiple surveys. It turns out that the published photometric transformations lead to substantial biases when applied to integrated photometry of galaxies from the corresponding catalogs. Here we present photometric transformations based on piece-wise linear approximations of integrated photometry of galaxies in the optical surveys SDSS, DECaLS, BASS, MzLS, DES, DELVE, KiDS, VST ATLAS, and the near-infrared surveys UKIDSS, UHS, VHS, and VIKING. We validate our transformations by constructing k -corrected color–magnitude diagrams of non-active galaxies and measuring the position and tightness of the “red sequence.” We also provide transformations for aperture magnitudes and show how they are affected by the image quality difference among the surveys. We present the implementation of the derived transformations in python and idl and also a web-based color transformation calculator for galaxies. By comparing DECaLS and DES, we identified systematic issues in DECaLS photometry for extended galaxies, which we attribute to the photometric software package used by DECaLS. As an application of our method, we compiled two multi-wavelength photometric catalogs for over 200,000 low- and intermediate-redshift galaxies originating from CfA FAST and Hectospec spectral archives.

Astronomy & Astrophysics↗
DOE OSTI2017

Development of a modular test system for the silicon sensor R&D of the ATLAS Upgrade

High Voltage CMOS sensors are a promising technology for tracking detectors in collider experiments. Extensive R&D studies are being carried out by the ATLAS Collaboration for a possible use of HV-CMOS in the High Luminosity LHC upgrade of the Inner Tracker detector. CaRIBOu (Control and Readout Itk BOard) is a modular test system developed to test Silicon based detectors. It currently includes five custom designed boards, a Xilinx ZC706 development board, FELIX (Front-End LInk eXchange) PCIe card and a host computer. A software program has been developed in Python to control the CaRIBOu hardware. CaRIBOu has been used in the testbeam of the HV-CMOS sensor AMS180v4 at CERN. Preliminary results have shown that the test system is very versatile. In conclusion, further development is ongoing to adapt to different sensors, and to make it available to various lab test stands.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Multi-angle Precession Electron Diffraction (MAPED): A Versatile Approach to 4D-STEM Precession

Precession of a converged beam during acquisition of a 4D-STEM dataset improves strain, orientation, and phase mapping accuracy by averaging over continuous angles of illumination. Precession experiments usually rely on integrated systems, where automatic alignments lead to fast, high-quality results. The dependence of these experiments on specific hardware and software is evident even when switching to nonintegrated detectors on a precession tool, as experimental set-up becomes challenging and time-consuming. Here, we introduce multi-angle precession electron diffraction (MAPED): a method to perform electron diffraction by collecting sequential 4D-STEM scans at different incident beam tilts. The multiple diffraction datasets are averaged together postacquisition, resulting in a single dataset that minimizes the impact of the curvature and orientation of the Ewald sphere relative to the crystal under study. Our results demonstrate that even four additional tilts improved measurement of material properties, namely strain and orientation, as compared to single-tilt 4D-STEM experiments. We show the versatility and flexibility of our MAPED approach with data collected on a number of microscopes with different hardware configurations and a variety of detectors.

4D-STEM↗
DOE OSTI2025

Three-dimensional modeling of hyphal fusion, branching, and nutrient transport in filamentous fungi

Fungi exhibit behaviors distinct from other microbes. Filamentous fungi grow by extending complex networks of branched filaments collectively referred to as the mycelium. These networks can expand over large distances and traverse low-nutrient areas by translocating nutrients through the filament network. This spatial characteristic makes filamentous fungi crucial for soil ecosystems, supporting stable microbial communities and promoting plant growth. However, simulating these behaviors is complex. The elongated nature of fungal compartments results in different mechanical interactions compared to the commonly modeled spherical bacteria. These detailed hyphal mechanics require specialized consideration and are often excluded from conventional fungal simulation packages. Additionally, the extensive fungal networks in nature demand computationally intensive simulations, necessitating high-performance algorithms. Therefore, realistic fungi simulations require specialized software. Here, we introduce a fungal modeling expansion to the high-performance biological modelling and interface exchange (bmx) software suite. bmx leverages adaptive mesh refinement in AMReX for chemical diffusion and incorporates a full mechanical model for bacterial cells, accelerated by GPUs. By extending bmx to model filamentous particles, we demonstrate the formation of complex filament networks through interactions like hyphal branching and fusion (anastomosis). We show that the networks produced match real-world fungal structures through various metrics. This work supports computational studies of fungal growth dynamics and can be adapted to investigate the growth of other filamentous structures in biology or materials science. The expanded-BMX package is open-sourced and is available online.

Cell mechanics↗
DOE OSTI2024

Radiocarbon of deep fine roots across four in lowland seasonal tropical forests in Panama

Objectives:Fine roots are key to ecosystem-scale nutrient, carbon (C), and water cycling, but our understanding of fine root traits variation within and among tropical forests, one of Earth’s most C-rich ecosystems, is limited. In 2022 and 2023, we aimed to explore differences in deep root characteristics among four lowland tropical forests in Panama, which vary in fertility and mean annual precipitation. We measured radiocarbon content (fraction modern [FM] and Δ14C) and δ13C of live fine roots at depths greater than 80 cm, up to 120 cm. The goal was to understand how deep root characteristics differ across these sites.Datasets included:The datasets provided include .csv and .xlsx files for radiocarbon content (fraction modern [FM] and Δ14C) and δ13C of live fine roots at depths greater than 80 cm, up to 120 cm collected from late 2017 to 2018 across four different forests. Additionally, a .kml file is provided with coordinates for all 32 plots included in the study across four forests (n = 8 plots per site). This dataset serves as baseline data before a throughfall exclusion experiment, Panama Rainforest Changes with Experimental Drying (PARCHED), was implemented. No special software is needed to open these files.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2013

FAST: A Fast Full System Cycle Accurate Computer Simulator (Early Career Principal Investigator Award Final Report)

This project is researching FAST, a methodology to build very fast, cycle-accurate full system computer simulators and building the first set of such simulators and the first set of tools to help construct those simulators. The methodology relies on a functional model that is a fast, full-system but not cycle-accurate simulator coupled with a timing model that is models the micro-architectural structure and arbitration of a computer system, but not its functionality. The way FAST simulators differ from other simulators partitioned in the same way is that the interface between the functional and timing model is optimized to minimize the need for round-trip communication. The optimized communication enables FAST timing models to be implemented in an FPGA and the functional model in software while still achieving extremely high performance. Our general strategy is to design a methodology and then prototype an example using that methodology to ensure the methodology is sound before committing to it. This strategy ensures that the selected methodology has been tested and provides an early example of the output of that methodology.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

REScheck Technical Support

PNNL has developed REScheck compliance materials for three different editions of the MEC (CABO 1992, 1993, and 1995) and all editions of the IECC (ICC 1998, 1999, 2003, 2006, 2007, 2009, 2012, 2015, and 2018). However, per DOE policy (dated March 18, 2014), only IECC editions 2009, 2012, 2015, and 2018 are presently supported in the software with the long-term intent to support only the three most recent code cycle editions. This report explains the methodology used to develop version 4.6.5 of the REScheck software in order to support compliance determination for IECC editions 2009, 2012, 2015, and 2018.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Advanced Monitoring and Control in the ANL METL Facility Using an Engineering Digital Twin

The potential benefits of using an engineering digital twin to achieve greater autonomy for monitoring and control functions in advanced reactors was investigated for the Mechanisms Engineering Test Loop (METL) liquid sodium facility at Argonne National Laboratory. The METL sodium purification system served as a representative system as it requires significant human in-the-loop interaction to accomplish its design function. The objective was to demonstrate how real-time operation could be automated while preserving oversight of the operator for ensuring that the system design functions are met. A digital twin model of the purification system was developed for both the cold trap purification loop and plugging meter diagnostic loop using information from the METL piping and instrumentation diagram (P&ID). Automated monitoring and diagnosis of component degradation in the METL facility was demonstrated in tests using the PRO-AID health monitoring software with the digital twin model incorporated in the library of components. Component failures were introduced and were successfully diagnosed in real time. These tests serve to demonstrate an advanced monitoring capability able to differentiate sensor degradation from component degradation, to generate a rank ordering of probabilities of different failure mechanisms that serves to circumvent the false alarm problem with current anomaly detection methods, and how facility monitoring can be transformed from anomaly detection to identification of a specific fault. Automated control of the purification system was demonstrated through simulations that exercised a model predictive controller designed using the digital twin model. Results of these simulations compared favorably with experimental data showing very good reference tracking response with negligible overshoot. In conclusion, these pilot tests and simulations successfully demonstrated the use of a digital twin for improved automation of monitoring and control. It was shown how the digital twin enables switching between control modes from cold trap operation where impurities are removed to plugging meter operation where impurity concentrations are measured.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

EcoSpec: Highly Equipped Tower-Based Hyperspectral and Thermal Infrared Automatic Remote Sensing System for Investigating Plant Responses to Environmental Changes

Despite an advanced ability to forecast ecosystem functions and climate at regional and global scales, little is known about relationships between local variations in water and carbon fluxes and large-scale phenomena. To enable data collection of local-scale ecosystem functions to support such investigations, we developed the EcoSpec system, a highly equipped remote sensing system that houses a hyperspectral radiometer (350–2500 nm) and five optical and infrared sensors in a compact tower. Its custom software controls the sequence and timing of movement of the sensors and system components and collects measurements at 12 locations around the tower. The data collected using the system was processed to remove sun-angle effects, and spectral vegetation indices computed from the data (i.e., the Normalized Difference Vegetation Index (NDVI), Normalized Difference Water Index (NDWI), Photochemical Reflectance Index (PRI), and Moisture Stress Index (MSI)) were compared with the fraction of photochemically active radiation (fPAR) and canopy temperature. The results showed that the NDVI, NDWI, and PRI were strongly correlated with fPAR; the MSI was correlated with canopy temperature at the diurnal scale. These correlations suggest that this type of near-surface remote sensing system would complement existing observatories to validate satellite remote sensing observations and link local and large-scale phenomena to improve our ability to forecast ecosystem functions and climate. The system is also relevant for precision agriculture to study crop growth, detect disease and pests, and compare traits of cultivars.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Results of Application Demonstration of Autonomous Offshore Power System Supporting Multiple Mobile and Static Assets

Over the last decade, several technology trends have dominated the innovation cycle within the terrestrial economy. Electrification, cloud computing, decarbonization, digitization, and autonomy are representative examples of macro trends that have had a rapid and fundamental impact; an impact felt not only in user-facing technologies, such as electric vehicles and software-as-a-service, but also in the infrastructure required to deliver those technologies. However, the scale of the impact and rollout has not been felt equally between the terrestrial and ocean economies. Here, the innovation cycle in the ocean economy is significantly slower. What accounts for the difference and are there solutions to this issue?

16 TIDAL AND WAVE POWER↗
DOE OSTI2022

Reekeekee- and roodoodooviruses, two different Microviridae clades constituted by the smallest DNA phages

Small circular single-stranded DNA viruses of the Microviridae family are both prevalent and diverse in all ecosystems. They usually harbor a genome between 4.3 and 6.3 kb, with a microvirus recently isolated from a marine Alphaproteobacteria being the smallest known genome of a DNA phage (4.248 kb). A subfamily, Amoyvirinae, has been proposed to classify this virus and other related small Alphaproteobacteria-infecting phages. Here, we report the discovery, in meta-omics data sets from various aquatic ecosystems, of sixteen complete microvirus genomes significantly smaller (2.991–3.692 kb) than known ones. Phylogenetic analysis reveals that these sixteen genomes represent two related, yet distinct and diverse, novel groups of microviruses—amoyviruses being their closest known relatives. We propose that these small microviruses are members of two tentatively named subfamilies Reekeekeevirinae and Roodoodoovirinae. As known microvirus genomes encode many overlapping and overprinted genes that are not identified by gene prediction software, we developed a new methodology to identify all genes based on protein conservation, amino acid composition, and selection pressure estimations. Surprisingly, only four to five genes could be identified per genome, with the number of overprinted genes lower than that in phiX174. These small genomes thus tend to have both a lower number of genes and a shorter length for each gene, leaving no place for variable gene regions that could harbor overprinted genes. Even more surprisingly, these two Microviridae groups had specific and different gene content, and major differences in their conserved protein sequences, highlighting that these two related groups of small genome microviruses use very different strategies to fulfill their lifecycle with such a small number of genes. The discovery of these genomes and the detailed prediction and annotation of their genome content expand our understanding of ssDNA phages in nature and are further evidence that these viruses have explored a wide range of possibilities during their long evolution.

59 BASIC BIOLOGICAL SCIENCES↗