Search NASA⌕ Search

SEARCH · Search NASA

Results for “Sequence Alignment”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 451 records · Page 25

Comparative Genome Structure, Secondary Metabolite, and Effector Coding Capacity across Cochliobolus Pathogens

The genomes of five Cochliobolus heterostrophus strains, two Cochliobolus sativus strains, three additional Cochliobolus species (Cochliobolus victoriae, Cochliobolus carbonum, Cochliobolus miyabeanus), and closely related Setosphaeria turcica were sequenced at the Joint Genome Institute (JGI). The datasets were used to identify SNPs between strains and species, unique genomic regions, core secondary metabolism genes, and small secreted protein (SSP) candidate effector encoding genes with a view towards pinpointing structural elements and gene content associated with specificity of these closely related fungi to different cereal hosts. Whole-genome alignment shows that three to five of each genome differs between strains of the same species, while a quarter of each genome differs between species. On average, SNP counts among field isolates of the same C. heterostrophus species are more than 25 higher than those between inbred lines and 50 lower than SNPs between Cochliobolus species. The suites of nonribosomal peptide synthetase (NRPS), polyketide synthase (PKS), and SSP encoding genes are astoundingly diverse among species but remarkably conserved among isolates of the same species, whether inbred or field strains, except for defining examples that map to unique genomic regions. Functional analysis of several strain-unique PKSs and NRPSs reveal a strong correlation with a role in virulence.

59 BASIC BIOLOGICAL SCIENCES↗

Probing the Consequences of Cubic Particle Shape and Applied Field on Colloidal Crystal Engineering with DNA

In a magnetic field, cubic Fe 3 O 4 nanoparticles exhibit assembly behavior that is a consequence of a competition between magnetic dipole-dipole and ligand interactions. In most cases, the interactions between short hydrophobic ligands dominate and dictate assembly outcome. To better tune the face-to-face interactions, cubic Fe3O4 nanoparticles were functionalized with DNA. Their assembly behaviors were investigated both with and without an applied magnetic field. Upon application of a field, the tilted orientation of cubes, enabled by the flexible DNA ligand shell, led to an unexpected crystallographic alignment of the entire superlattice, as opposed to just the individual particles, along the field direction as revealed by small and wide-angle X-ray scattering. This observation is dependent upon DNA length and sequence and cube dimensions. Taken together, these studies show how combining physical and chemical control can expand the possibilities of crystal engineering with DNA.

36 MATERIALS SCIENCE↗

Peptide supramolecular materials for therapeutics

Supramolecular assembly of peptide-based monomers into nanostructures offers many promising applications in advanced therapies. In this Tutorial Review, we introduce molecular designs to control the structure and potential biological function of supramolecular assemblies. Here, an emphasis is placed on peptide-based supramolecular nanostructures that are intentionally designed to signal cells, either directly through the incorporation of amino acid sequences that activate receptors or indirectly by recruiting native signals such as growth factors. Additionally, we describe the use and future potential of hierarchical structures, such as single molecules that assemble into nanoscale fibers which then align to form macroscopic strings; the strings can then serve as scaffolds for cell growth, proliferation, and differentiation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Expanded Classification of Hepatitis C Virus Into 7 Genotypes and 67 Subtypes: Updated Criteria and Genotype Assignment Web Resource

The 2005 consensus proposal for the classification of hepatitis C virus (HCV) presented an agreed and uniform nomenclature for HCV variants and the criteria for their assignment into genotypes and subtypes. Since its publication, the available dataset of HCV sequences has vastly expanded through advancement in nucleotide sequencing technologies and an increasing focus on the role of HCV genetic variation in disease and treatment outcomes. The current study represents a major update to the previous consensus HCV classification, incorporating additional sequence information derived from over 1,300 (near-)complete genome sequences of HCV available on public databases in May 2013. Analysis resolved several nomenclature conflicts between genotype designations and using consensus criteria created a classification of HCV into seven confirmed genotypes and 67 subtypes. There are 21 additional complete coding region sequences of unassigned subtype. The study additionally describes the development of a Web resource hosted by the International Committee for Taxonomy of Viruses (ICTV) that maintains and regularly updates tables of reference isolates, accession numbers, and annotated alignments (http://talk.ictvonline.org/links/hcv/hcv-classification.htm). The Flaviviridae Study Group urges those who need to check or propose new genotypes or subtypes of HCV to contact the Study Group in advance of publication to avoid nomenclature conflicts appearing in the literature. While the criteria for assigning genotypes and subtypes remain unchanged from previous consensus proposals, changes are proposed in the assignment of provisional subtypes, subtype numbering beyond “w,” and the nomenclature of intergenotypic recombinant. Conclusion: This study represents an important reference point for the consensus classification of HCV variants that will be of value to researchers working in clinical and basic science fields.

59 BASIC BIOLOGICAL SCIENCES↗

Catalytic Upgrading of Carbohydrates in Waste Streams to Hydrocarbons: Paper Sludge to Fuel Project (PStF)

This research, funded by the Department of Energy Bioenergy Technologies Office (DOE BETO), aims to understand the barriers and assess opportunities for transforming carbohydrates in paper sludge, a solid lignocellulosic residue from the pulp and paper industry. The goal is to turn the paper sludge into a liquid hydrocarbon product that can be blended into jet fuel, both economically and sustainably. Research groups at North Carolina State University, National Renewable Energy Laboratory, and Yale University collaborated synergistically, leveraging their expertise in pulp and paper operations, biomass deconstruction, and catalytic upgrading to propose a pathway that efficiently captures the energy in paper sludge. Findings from this study could potentially contribute to advancing biomass conversion technologies, aligning with the efforts of the U.S. DOE BETO. This report documents the experimental and simulation results of a biochemical and catalytic pathway designed to transform paper sludge into a liquid hydrocarbon product. The process involves a sequence of seven steps, including ash removal, carbohydrate enzymatic hydrolysis, sugar dehydration, solvent recovery, aldol-condensation between furans and ketone, hydrogenation, and hydrodeoxygenation to obtain a hydrocarbon blend in the ~C10 range. The experimental efforts from the initiation of the project were guided by techno-economic analysis (TEA) and sensitivity analysis results including around seventy-eight operational parameters. This methodology facilitated the efficient use of resources over time. This study relies on detailed process simulations and TEA to determine the minimum fuel selling price (MFSP) for the hydrocarbon fuel product. Preliminary TEA results led to the evaluation of eight case studies considering alternative dehydration co-solvents, operational settings, and biorefinery layouts. Finally, the life cycle assessment of twenty-eight scenarios, comparing various dehydration co-solvents, fuel sources, chemical feedstock sources, side product utilization, and other process settings, was conducted. Landfilling scenarios with and without landfill gas recovery were also estimated, analyzed, and compared.

09 BIOMASS FUELS↗

Bioaerosols as indicators of central Arctic ice nucleating particle sources

The Arctic is warming at a rapid rate, with implications for microbial communities as the ecosystems change. Some microbes and biogenic materials can affect the persistence of long-lived mixed-phase clouds by serving as ice nucleating particles (INPs). The presence of INPs modulates the cloud phase, and long-term measurements are important to elucidate their seasonal sources and to predict future change. The Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition in 2019–2020 provided the first year-long measurements of bioaerosols and INPs in the central Arctic, with 3 d filters for amplicon sequencing and cumulative INP concentrations from −5 to −30 °C. Here, we investigated the INP seasonal cycle and its relation to the seasonal cycle of bacteria and eukaryotes. INPs were greatly elevated and compositionally similar in summer, aligning with a greater prevalence of local bioaerosol sources, but, despite this, a diverse mixture of sources (marine and terrestrial) was present all times. A common broader Arctic INP population is hypothesized for much of the year by comparable coincident data collected in Svalbard and a sensitivity of both the INPs and bioaerosols to large-scale events.

Barry, Kevin R. [Colorado State Univ., Fort Collin↗

Consolidated Nuclear Security, LLC, Excess Facility Disposition Program - 17289

Consolidated Nuclear Security, LLC (CNS), the management and operating contractor for the Y-12 National Security Complex (Y-12) in Oak Ridge, Tennessee, and the Pantex Plant (Pantex) in Amarillo, Texas, has developed an integrated strategy to address badly deteriorating excess facilities presenting risks to people, the environment, and the National Nuclear Security Administration (NNSA) mission. The comprehensive, integrated Excess Facilities Disposition Program (EFDP) incorporates NNSA scope for excess facility risk reduction, preparation of facilities to be transferred to the U.S. Department of Energy (DOE) Office of Environmental Management (EM) for demolition, and demolition of non-process contaminated facilities. The Program has been developed using the NNSA Office of Safety, Infrastructure, and Operations (NA-50) Program Management Plan as a guide. The Program has undertaken actions to address the most immediate risks associated with CNS excess facilities. At Y-12, new foam roofs have been installed on Buildings 9201-5 (Alpha 5), 9204-4 (Beta 4), and 9206 to minimize further facility degradation. Tanks and dikes outside of Alpha 5 and Beta 4 are also being disconnected, drained, and filled with concrete. In addition, Beta 4 will be getting a new electrical system as its Manhattan Project-era system will be replaced by a temporary 'construction power' setup that will provide as-needed electricity to the building that is safer, more reliable, and will ease demolition activities in the future. The EFDP is also focusing on the basement of Alpha 5 and is planning to remove the water from the aging facility's flooded basement, which in some areas contains up to 54 inches of water. The EFDP is also planning several efforts to remove excess materials from Alpha 5 and Beta 4, de-inventory the Alpha 5 Mercury House System and numerous transformers, and de-inventory Beta 4 hydraulic pumps and systems. Efforts are also underway to deactivate and decommission non-process-contaminated facilities at Y-12 and at Pantex. The EFDP has prioritized the performance of work to reduce risk while at the same time supporting the preparation of process-contaminated Y-12 facilities for transfer to EM. The Program aligns with the EM goal to have the NNSA scope completed in Alpha 5 and Beta 4 by the end of FY 2020, and Building 9206 completed by the end of FY 2024. These dates and project sequencing will support the Y-12-related DOE EM cleanup initiatives. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Deep RNA Sequencing Reveals Hidden Features and Dynamics of Early Gene Transcription in Paramecium bursaria Chlorella Virus 1

Paramecium bursaria chlorella virus 1 (PBCV-1) is the prototype of the genus Chlorovirus (family Phycodnaviridae) that infects the unicellular, eukaryotic green alga Chlorella variabilis NC64A. The 331-kb PBCV-1 genome contains 416 major open reading frames. A mRNA-seq approach was used to analyze PBCV-1 transcriptomes at 6 progressive times during the first hour of infection. The alignment of 17 million reads to the PBCV-1 genome allowed the construction of single-base transcriptome maps. Significant transcription was detected for a subset of 50 viral genes as soon as 7 min after infection. By 20 min post infection (p.i.), transcripts were detected for most PBCV-1 genes and transcript levels continued to increase globally up to 60 min p.i., at which time 41% or the poly (A+)-containing RNAs in the infected cells mapped to the PBCV-1 genome. For some viral genes, the number of transcripts in the latter time points (20 to 60 min p.i.) was much higher than that of the most highly expressed host genes. RNA-seq data revealed putative polyadenylation signal sequences in PBCV-1 genes that were identical to the polyadenylation signal AAUAAA of green algae. Several transcripts have an RNA fragment excised. However, the frequency of excision and the resulting putative shortened protein products suggest that most of these excision events have no functional role but are probably the result of the activity of misled splicesomes.

59 BASIC BIOLOGICAL SCIENCES↗

Genome mapping of quantitative trait loci (QTL) controlling domestication traits of intermediate wheatgrass ( Thinopyrum intermedium )

Allohexaploid (2n = 6x = 42) intermediate wheatgrass (Thinopyrum intermedium), abbreviated IWG, is an outcrossing perennial grass belonging to the tertiary gene pool of wheat. Perenniality would be valuable option for grain production, but attempts to introgress this complex trait from wheat-Thinopyrum hybrids have not been commercially successful. Efforts to breed IWG itself as a dual-purpose forage and grain crop have demonstrated useful progress and applications, but grain yields are significantly less than wheat. Therefore, genetic and physical maps have been developed to accelerate domestication of IWG. Herein, these maps were used to identify quantitative trait loci (QTLs) and candidate genes associated with IWG grain production traits in a family of 266 full-sib progenies derived from two heterozygous parents, M26 and M35. Transgressive segregation was observed for 17 traits related to seed size, shattering, threshing, inflorescence capacity, fertility, stem size, and flowering time. A total of 111 QTLs were detected in 36 different regions using 3826 genotype-by-sequence markers in 21 linkage groups. The most prominent QTL had a LOD score of 15 with synergistic effects of 29% and 22% over the family means for seed retention and percentage of naked seeds, respectively. Many QTLs aligned with one or more IWG gene models corresponding to 42 possible domestication orthogenes including the wheat Q and RHT genes. A cluster of seed-size and fertility QTLs showed possible alignment to a putative Z self-incompatibility gene, which could have detrimental grain-yield effects when genetic variability is low. Finally, these findings elucidate pathways and possible hurdles in the domestication of IWG.

59 BASIC BIOLOGICAL SCIENCES↗

Data Analysis of the 2020 Central Idaho Mainshock-Aftershock Sequence

In an effort to inform the Senior Seismic Hazard Analysis Committee for the Idaho National Laboratory, we provide an improved aftershock catalog related to the March 31, 2020, Mw6.5 Stanley, Idaho earthquake from picks related to a temporary network of two real-time and 15 non-telemetered seismometers within the epicentral area. From the permanent and temporary (XP) real-time network, the USGS cataloged 1,946 aftershocks between April 1, 2020 and October 31, 2020. To improve aftershock location and magnitudes, we manually picked arrival times of P and S waves from off-line stations in the XP temporary network, generated a new crustal velocity model, and independently relocated each event using the HypoDD double-difference earthquake algorithm. We created our new velocity model from existing broadband and active source seismic campaign data that were acquired near the epicentral region prior to the 2020 earthquake. We compare arrival time differences, epicentral locations and depths between aftershocks recorded with the two catalogs. We find the addition of local stations provides tighter aftershock clustering that suggests an improved aftershock locations. To detect lower magnitude events, we employed deep learning. Our method solves common problems associated with detecting many events that have a low signal-to-noise ratio. From the machine learning database, we detected more than 74,000 aftershocks. Based on the number of identified earthquakes and Gutenberg-Richter relationships derived from the USGS catalog, we estimate that we have reduced the completion magnitude for the Stanley earthquake sequence to below M1 using this machine learning approach. We located each aftershock with our new velocity model. Our new velocity model and picks suggests aftershocks occurred mostly at shallower depths than assessed in the USGS catalog. These aftershocks align along two linear trends that suggest the activation of two unnamed primary faults.

58 GEOSCIENCES↗

ORBIT AND STELLAR PROPERTIES OF THE YOUNG TRIPLE V807 Tau

We present new orbital measurements of the pre-main-sequence triple system, V807 Tau, using adaptive optics imaging at the Keck Observatory. We computed an orbit for the close pair, V807 Tau Ba-Bb, with a period of 12.312 {+-} 0.058 years and a semi-major axis of 38.59 {+-} 0.16 mas. By modeling the center of mass motion of the components in the close pair relative to the wide component, V807 Tau A, we measured a mass ratio of 0.843 {+-} 0.050 for Bb/Ba. Combined with the total mass from the relative orbit, we derived individual masses of M{sub Ba} = 0.564 {+-} 0.018 (d/140 pc){sup 3} M{sub Sun} and M{sub Bb} = 0.476 {+-} 0.017 (d/140 pc){sup 3} M{sub Sun} at an average distance of 140 pc to the Taurus star-forming region. We computed spectral energy distributions to determine the luminosities of the three components. We also measured their spectral types, effective temperatures, and rotational velocities based on spatially resolved spectra obtained at the Keck Observatory. If the rotational axes are aligned, then the projected rotational velocities imply that V807 Tau Ba and Bb are rotating much faster than V807 Tau A. The uncertainty in the stellar effective temperatures and distance to the system currently limit the comparison of our dynamical mass measurements with predictions based on evolutionary tracks for pre-main-sequence stars. We also report preliminary results from a program to map the 3.6 cm radio emission from V807 Tau using the Very Long Baseline Array. With continued monitoring, these observations will provide a precise parallax for placing the dynamical masses on an absolute scale.

79 ASTRONOMY AND ASTROPHYSICS↗

Laser wakefield acceleration driven by a discrete flying focus

Laser wakefield acceleration (LWFA) may enable the next generation of TeV-scale lepton colliders. Reaching such energies will likely require multiple LWFA stages to overcome limitations on the energy gain achievable in a single stage. The use of stages, however, introduces challenges such as alignment, adiabatic matching between stages, and a lower average accelerating gradient. Here, we propose a discrete flying focus that can deliver higher energy gain in a single stage, thereby reducing the number of stages required for a target energy. A sequence of laser pulses with staggered focal points and delays drives a plasma wave in which an electron beam experiences a near-constant accelerating gradient over distances beyond those attainable with a conventional pulse. Simulations demonstrate that a discrete flying focus with a total energy of 150 J can transfer 40 GeV per electron to a 50-pC beam in a single 30-cm stage, corresponding to 50 dephasing lengths.

laser wakefield acceleration↗

Further Evidence of Modified Spin-down in Sun-like Stars: Pileups in the Temperature–Period Distribution

We combine stellar surface rotation periods determined from NASA's Kepler mission with spectroscopic temperatures to demonstrate the existence of pileups at the long-period and short-period edges of the temperature–period distribution for main-sequence stars with temperatures exceeding ~5500 K. The long-period pileup is well described by a curve of constant Rossby number, with a critical value of Ro crit ≲ Ro ⊙ . The long-period pileup was predicted by van Saders et al. as a consequence of weakened magnetic braking, in which wind-driven angular momentum losses cease once stars reach a critical Rossby number. Stars in the long-period pileup are found to have a wide range of ages (~2–6 Gyr), meaning that, along the pileup, rotation period is strongly predictive of a star's surface temperature but weakly predictive of its age. The short-period pileup, which is also well described by a curve of constant Rossby number, is not a prediction of the weakened magnetic braking hypothesis but may instead be related to a phase of slowed surface spin-down due to core-envelope coupling. The same mechanism was proposed by Curtis et al. to explain the overlapping rotation sequences of low-mass members of differently aged open clusters. The relative dearth of stars with intermediate rotation periods between the short- and long-period pileups is also well described by a curve of constant Rossby number, which aligns with the period gap initially discovered by McQuillan et al. in M-type stars. These observations provide further support for the hypothesis that the period gap is due to stellar astrophysics, rather than a nonuniform star formation history in the Kepler field.

79 ASTRONOMY AND ASTROPHYSICS↗

Electron beam-induced nanopores in Bernal-stacked hexagonal boron nitride

Controlling the size and shape of nanopores in two-dimensional materials is a key challenge in applications such as DNA sequencing, sieving, and quantum emission in artificial atoms. Here, we investigate experimentally and theoretically triangular vacancies in (unconventional) Bernal-stacked AB-h-BN formed using a high-energy electron beam. Due to the geometric configuration of AB-h-BN, triangular pores in different layers are aligned, and their sizes are controlled by the duration of the electron irradiation. Interlayer covalent bonding at the vacancy edge is not favored, as opposed to what occurs in the more common AA'-stacked BN. A variety of monolayer, concentric and bilayer pores in bilayer AB-h-BN are observed in high-resolution transmission electron microscopy and characterized using ab initio simulations. Bilayer pores in AB-h-BN are commonly formed, and grow without breaking the bilayer character. Nanopores in AB-h-BN exhibit a wide range of electronic properties, ranging from half-metallic to non-magnetic and magnetic semiconducting. Therefore, because of the controllability of the pore size, the electronic structure is also highly controllable in these systems, and can potentially be tuned for particular applications.

2D materials↗

Comparative analysis of the transcriptome across distant species

The transcriptome is the readout of the genome. Identifying common features in it across distant species can reveal fundamental principles. To this end, the ENCODE and modENCODE consortia have generated large amounts of matched RNA-sequencing data for human, worm and fly. Uniform processing and comprehensive annotation of these data allow comparison across metazoan phyla, extending beyond earlier within-phylum transcriptome comparisons and revealing ancient, conserved features[1,2,3,4,5,6]. Specifically, we discovered co-expression modules shared across animals, many of which are enriched in developmental genes. Moreover, we used expression patterns to align the stages in worm and fly development, finding a novel pairing between worm embryo and fly pupae, in addition to the expected embryo-to-embryo and larvae-to-larvae pairings. Furthermore, we found that the extent of non-canonical, non-coding transcription is similar in each organism, per base-pair. Finally, we found in all three organisms the gene-expression levels, both coding and non-coding, can be quantitatively predicted from chromatin features at the promoter using a “universal model,” based on a single set of organism-independent parameters. Our comparison used the ENCODE-modENCODE RNA resource (Fig. ED1). This resource comprises: (1) deeply sequenced RNA-Seq data from many distinct samples from all three organisms; (2) comprehensive annotation of transcribed elements and (3) uniformly processed, standardized analysis files, focusing on non-coding transcription and expression patterns. Where practical, these datasets match comparable samples across organisms and to other types of functional genomics data. In total, the resource contains 575 different experiments containing >67B sequence reads. It encompasses many different RNA types, including poly(A)+, poly(A)- and ribosomal-RNA-depleted RNA and short and long RNA.

59 BASIC BIOLOGICAL SCIENCES↗

CPGminer: An Interactive Dashboard to Explore the Genomic Features and Taxonomy of Complete Prokaryotic Genomes

Prokaryotes, the earliest forms of life on Earth, play crucial roles in global biogeochemical processes in virtually all ecosystems. The ever-increasing amount of prokaryotic genome sequencing data provides a wealth of information to examine fundamental and applied questions through systematic genome comparison. Genomic features, such as genome size and GC content, and taxonomy-centric genomic features of complete prokaryotic genomes (CPGs) are crucial for various fields of microbial research and education, yet they are often overlooked. Additionally, creating systematically curated datasets that align with research concerns is an essential yet challenging task for wet-lab researchers. In this study, we introduce CPGminer, a user-friendly tool that allows researchers to quickly and easily examine the genomic features and taxonomy of CPGs and curate genome datasets. We also provide several examples to demonstrate its practical utility in addressing descriptive questions.

59 BASIC BIOLOGICAL SCIENCES↗

SPIN-ORBIT ALIGNMENT FOR THE CIRCUMBINARY PLANET HOST KEPLER-16 A

Kepler-16 is an eccentric low-mass eclipsing binary with a circumbinary transiting planet. Here, we investigate the angular momentum of the primary star, based on Kepler photometry and Keck spectroscopy. The primary star's rotation period is 35.1 {+-} 1.0 days, and its projected obliquity with respect to the stellar binary orbit is 1.{sup 0}6 {+-} 2.{sup 0}4. Therefore, the three largest sources of angular momentum-the stellar orbit, the planetary orbit, and the primary's rotation-are all closely aligned. This finding supports a formation scenario involving accretion from a single disk. Alternatively, tides may have realigned the stars despite their relatively wide separation (0.2 AU), a hypothesis that is supported by the agreement between the measured rotation period and the 'pseudosynchronous' period of tidal evolution theory. The rotation period, chromospheric activity level, and fractional light variations suggest a main-sequence age of 2-4 Gyr. Evolutionary models of low-mass stars can match the observed masses and radii of the primary and secondary stars to within about 3%.

79 ASTRONOMY AND ASTROPHYSICS↗

Model of metabolism and gene expression predicts proteome allocation in Pseudomonas putida

Abstract The genome-scale model of metabolism and gene expression (ME-model) forPseudomonas putidaKT2440,iPpu1676-ME, provides a comprehensive representation of biosynthetic costs and proteome allocation. Compared to a metabolic-only model,iPpu1676-ME significantly expands on gene expression, macromolecular assembly, and cofactor utilization, enabling accurate growth predictions without additional constraints. Multi-omics analysis using RNA sequencing and ribosomal profiling data revealed translational prioritization inP. putida, with core pathways, such as nicotinamide biosynthesis and queuosine metabolism, exhibiting higher translational efficiency, while secondary pathways displayed lower priority. Notably, the ME-model significantly outperformed the M-model in alignment with multi-omics data, thereby validating its predictive capacity. Thus,iPpu1676-ME offers valuable insights intoP. putida’s proteome allocation and presents a powerful tool for understanding resource allocation in this industrially relevant microorganism.

Mathematical & Computational Biology↗