Search NASA⌕ Search

SEARCH · Search NASA

Results for “sequencing experiments”

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 325 records · Page 18

Counting Kmers for Biological Sequences at Large Scale

Counting the abundance of all the distinct kmers in biological sequence data is a fundamental step in bioinformatics. These applications include de novo genome assembly, error correction, etc. With the development of sequencing technology, the sequence data in a single project can reach Petabyte-scale or Terabyte-scale nucleotides. Counting demand for the abundance of these sequencing data is beyond the memory and computing capacity of single computing node, and how to process it efficiently is a challenge on a high-performance computing cluster. As such, we propose SWAPCounter, a highly scalable distributed approach for kmer counting. This approach is embedded with an MPI streaming I/O module for loading huge data set at high speed, and a counting bloom filter module for both memory and communication efficiency. By overlapping all the counting steps, SWAPCounter achieves high scalability with high parallel efficiency. In this work, the experimental results indicate that SWAPCounter has competitive performance with two other tools on shared memory environment, KMC 2 , and MSPKmerCounter. Moreover, SWAPCounter also shows the highest scalability under strong scaling experiments. In our experiment on Cetus supercomputer, SWAPCounter scales to 32,768 cores with 79% parallel efficiency (using 2048 cores as baseline) when processing 4 TB sequence data of 1000 Genomes. The source code of SWAPCounter is publicly available at https://github.com/mengjintao/SWAPCounter.

59 BASIC BIOLOGICAL SCIENCES↗

Microbial Diversity in the Deep Sandstones of a Carbon Sequestration Well

A lack of data on microbes in the deep subsurface hinders our ability to predict how these microbial communities will respond to carbon sequestration. Previous studies predicted how microbial communities may shift due to an influx of carbon dioxide (CO2) by 1) examining waters with naturally-occurring elevated CO2 concentration; and/or 2) via experiments involving exposure of laboratory microbial strains to elevated concentrations of CO2 relative to their native subsurface environment. To identify microbes naturally found in carbon sequestration target reservoirs, water samples were collected from three potential target reservoirs in the Powder River Basin, Wyoming which include the Hulett, Lakota, and Minnelusa sandstones. Water samples were collected at depths of 8,060 ft (Lakota), 8,330 ft (Hulett), 9,380 ft (Minnelusa), 9,463 ft (Minnelusa), and 9,544 ft (Minnelusa) for microbial DNA analysis as well as water chemistry and stable isotope analyses to look for any microbial trends that could be related to water chemistry alone. Identification of microbial communities present in the sampled water was determined via 16S rRNA sequencing; these data were compared to 16S rRNA isolated from drilling mud to distinguish native from any microbes introduced via drilling processes. Enrichment experiments reveal that the collected water samples contain active (e.g., growing) microbes and how the microbial community may adapt upon the introduction of CO2.

Wawrousek, K↗

Preliminary Design of Critical Experiments Involving Commercially Available B 4 C Neutron Absorber Plates with Low-Enriched UO 2 Fuel

The International Criticality Safety Benchmark Evaluation Project (ICSBEP) is an initiative to provide high-quality benchmark data in a standardized format for criticality safety analysts to use to validate calculational tools and nuclear data. In the last released 2022 version, the ICSBEP handbook contains over 5,000 critical experiment descriptions, results, and associated models. This paper describes the critical configurations proposed and calculations performed as part of the conceptual design phase of Integral Experiment Request (IER) 554 (IER-554). The goal of IER-554 is to add to the benchmark experiments available to the criticality safety community by designing critical experiments that can analyze how adding commercially available Boralcan plates to an assembly of low-enriched UO2 fuel rods affects the effective neutron multiplication factor (k eff ). Boralcan, a neutron absorber product developed by Rio Tinto, is made of B4C in an Al 1100 matrix and is commonly used for criticality suppression in fuel storage pools. In the current version of the ICSBEP handbook, only a few dated critical experiments involve B4C materials, and none involve Boralcan in a thin plate shape. The experiments designed as part of IER-554 are intended to be performed at the Sandia Pulsed Reactor Facility/Critical Experiments (SPRF/CX) apparatus at Sandia National Laboratories. SPRF/CX is a well-characterized assembly considered trustworthy by the benchmarking community because of the numerous high-quality evaluations with very low experimental uncertainties included in the ICSBEP handbook. Figure 1 illustrates the assembly for a particular configuration from LCT-078, one of the published benchmarks in the ICSBEP handbook. The experiments will be moderated and reflected by light water at ambient atmospheric pressure. Before the plates are inserted in the critical assembly, they will be characterized with x-ray computed tomography (XCT) to identify the sizes and distribution of the B 4 C powder particles inside the plates. The preliminary design calculations were performed with SCALE 6.3.0 using the KENO V.a sequence for the criticality calculations, the CE-TSUNAMI-3D sequence for the sensitivity and uncertainty studies, and the ENDF/B-VII.1 Continuous Energy cross section library. All the calculation results presented have 10 pcm statistical uncertainties. A set of critical experiments with Boralcan plates could increase the confidence of the criticality safety community in its modeling methods when using this type of neutron absorber material. Experiments would also help validate k eff calculations, and the industry could use these validations to change the B loading credit limits from the US Nuclear Regulatory Commission standard review plan for dry cask storage of spent nuclear fuel. This paper summarizes only part of the analysis documented in the preliminary design report.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Constant-time 2D and 3D through-bond correlation NMR spectroscopy of solids under 60 kHz MAS

Establishing connectivity and proximity of nuclei is an important step in elucidating the structure and dynamics of molecules in solids using magic angle spinning (MAS) NMR spectroscopy. Although recent studies have successfully demonstrated the feasibility of proton-detected multidimensional solid-state NMR experiments under ultrafast-MAS frequencies and obtaining high-resolution spectral lines of protons, assignment of proton resonances is a major challenge. In this study, we first re-visit and demonstrate the feasibility of 2D constant-time uniform-sign cross-peak correlation (CTUC-COSY) NMR experiment on rigid solids under ultrafast-MAS conditions, where the sensitivity of the experiment is enhanced by the reduced spin-spin relaxation rate and the use of low radio-frequency power for heteronuclear decoupling during the evolution intervals of the pulse sequence. In addition, we experimentally demonstrate the performance of a proton-detected pulse sequence to obtain a 3D {sup 1}H/{sup 13}C/{sup 1}H chemical shift correlation spectrum by incorporating an additional cross-polarization period in the CTUC-COSY pulse sequence to enable proton chemical shift evolution and proton detection in the incrementable t{sub 1} and t{sub 3} periods, respectively. In addition to through-space and through-bond {sup 13}C/{sup 1}H and {sup 13}C/{sup 13}C chemical shift correlations, the 3D {sup 1}H/{sup 13}C/{sup 1}H experiment also provides a COSY-type {sup 1}H/{sup 1}H chemical shift correlation spectrum, where only the chemical shifts of those protons, which are bonded to two neighboring carbons, are correlated. By extracting 2D F1/F3 slices ({sup 1}H/{sup 1}H chemical shift correlation spectrum) at different {sup 13}C chemical shift frequencies from the 3D {sup 1}H/{sup 13}C/{sup 1}H spectrum, resonances of proton atoms located close to a specific carbon atom can be identified. Overall, the through-bond and through-space homonuclear/heteronuclear proximities determined from the 3D {sup 1}H/{sup 13}C/{sup 1}H experiment would be useful to study the structure and dynamics of a variety of chemical and biological solids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A multi-stage stochastic programming model for adaptive biomass processing operation under uncertainty

Variations of physical and chemical characteristics of biomass reduce equipment utilization and increase operational costs of biomass processing. Biomass processing facilities use sensors to monitor the changes in biomass characteristics. Integrating sensory data into the operational decisions in biomass processing will increase its flexibility to the changing biomass conditions. In this paper, we propose a multi-stage stochastic programming model that minimizes the expected operational costs by identifying the initial inventory level and creating an operational decision policy for equipment speed settings. These policies take the sensory information data and the current biomass inventory level as inputs to dynamically adjust inventory levels and equipment settings according to the changes in the biomass' characteristics. We ensure that a prescribed target reactor utilization is consistently achieved by penalizing the violation of the target reactor feeding rate. A case study is developed using real-world data collected at Idaho National Laboratory's biomass processing facility. We show the value of multi-stage stochastic programming from an extensive computational experiment. Our sensitivity analysis indicates that updating the infeed rate of the system, the processing speed of equipment, and bale sequencing based on the moisture level of biomass improves the processing rate of the reactor and reduces operating costs.

09 BIOMASS FUELS↗

Multistage coupling of independent laser-plasma accelerators

Laser-plasma accelerators (LPAs) are capable of accelerating charged particles to very high energies in very compact structures. In theory, therefore, they offer advantages over conventional, large-scale particle accelerators. However, the energy gain in a single-stage LPA can be limited by laser diffraction, dephasing, electron-beam loading and laser-energy depletion. The problem of laser diffraction can be addressed by using laser-pulse guiding and preformed plasma waveguides to maintain the required laser intensity over distances of many Rayleigh lengths; dephasing can be mitigated by longitudinal tailoring of the plasma density; and beam loading can be controlled by proper shaping of the electron beam. To increase the beam energy further, it is necessary to tackle the problem of the depletion of laser energy, by sequencing the accelerator into stages, each powered by a separate laser pulse. In this work, we present results from an experiment that demonstrates such staging. Two LPA stages were coupled over a short distance (as is needed to preserve the average acceleration gradient) by a plasma mirror. Stable electron beams from a first LPA were focused to a twenty-micrometre radius-by a discharge capillary-based active plasma lens-into a second LPA, such that the beams interacted with the wakefield excited by a separate laser. Staged acceleration by the wakefield of the second stage is detected via an energy gain of 100 megaelectronvolts for a subset of the electron beam. Changing the arrival time of the electron beam with respect to the second-stage laser pulse allowed us to reconstruct the temporal wakefield structure and to determine the plasma density. Our results indicate that the fundamental limitation to energy gain presented by laser depletion can be overcome by using staged acceleration, suggesting a way of reaching the electron energies required for collider applications.

43 PARTICLE ACCELERATORS↗

Systematic and stochastic influences on the performance of the MinION nanopore sequencer across a range of nucleotide bias

Emerging sequencing technologies are allowing us to characterize environmental, clinical and laboratory samples with increasing speed and detail, including real-time analysis and interpretation of data. One example of this is being able to rapidly and accurately detect a wide range of pathogenic organisms, both in the clinic and the field. Genomes can have radically different GC content however, such that accurate sequence analysis can be challenging depending upon the technology used. Here, we have characterized the performance of the Oxford MinION nanopore sequencer for detection and evaluation of organisms with a range of genomic nucleotide bias. We have diagnosed the quality of base-calling across individual reads and discovered that the position within the read affects base-calling and quality scores. Finally, we have evaluated the performance of the current state-of-the-art neural network-based MinION basecaller, characterizing its behavior with respect to systemic errors as well as context- and sequence-specific errors. Overall, we present a detailed characterization the capabilities of the MinION in terms of generating high-accuracy sequence data from genomes with a wide range of nucleotide content. This study provides a framework for designing the appropriate experiments that are the likely to lead to accurate and rapid field-forward diagnostics.

59 BASIC BIOLOGICAL SCIENCES↗

Parallelized Kalman-Filter-Based Reconstruction of Particle Tracks on Many-Core Architectures

Faced with physical and energy density limitations on clock speed, contemporary microprocessor designers have increasingly turned to on-chip parallelism for performance gains. Algorithms should accordingly be designed with ample amounts of fine-grained parallelism if they are to realize the full performance of the hardware. This requirement can be challenging for algorithms that are naturally expressed as a sequence of small-matrix operations, such as the Kalman filter methods widely in use in high-energy physics experiments. In the High-Luminosity Large Hadron Collider (HL-LHC), for example, one of the dominant computational problems is expected to be finding and fitting charged-particle tracks during event reconstruction; today, the most common track-finding methods are those based on the Kalman filter. Experience at the LHC, both in the trigger and offline, has shown that these methods are robust and provide high physics performance. Previously we reported the significant parallel speedups that resulted from our efforts to adapt Kalman-filter-based tracking to many-core architectures such as Intel Xeon Phi. Here we report on how effectively those techniques can be applied to more realistic detector configurations and event complexity.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

The Progenitor Dependence of Core-collapse Supernovae from Three-dimensional Simulations with Progenitor Models of 12–40 M ⊙

We present a first study of the progenitor star dependence of the three-dimensional (3D) neutrino mechanism of core-collapse supernovae. We employ full 3D general-relativistic multi-group neutrino radiation-hydrodynamics and simulate the postbounce evolutions of progenitors with zero-age main sequence masses of 12, 15, 20, 27, and 40 M ⊙ . All progenitors, with the exception of the 12 M ⊙ star, experience shock runaway by the end of their simulations. In most cases, a strongly asymmetric explosion will result. We demonstrate three qualitatively distinct evolutions that imply a complex dependence of explosion dynamics on progenitor density structure, neutrino heating, and 3D flow. (1) Progenitors with massive cores, shallow density profiles, and high post-core-bounce accretion rates experience very strong neutrino heating and neutrino-driven turbulent convection, leading to early shock runaway. Accretion continues at a high rate, likely leading to black hole formation. (2) Intermediate progenitors experience neutrino-driven, turbulence-aided explosions triggered by the arrival of density discontinuities at the shock. These occur typically at the silicon/silicon–oxygen shell boundary. (3) Progenitors with small cores and density profiles without strong discontinuities experience shock recession and develop the 3D standing-accretion shock instability (SASI). Shock runaway ensues late, once declining accretion rate, SASI, and neutrino-driven convection create favorable conditions. These differences in explosion times and dynamics conclude in a non-monotonic relationship between progenitor and compact remnant mass.

79 ASTRONOMY AND ASTROPHYSICS↗

PROPERTIES OF CARBON–OXYGEN WHITE DWARFS FROM MONTE CARLO STELLAR MODELS

We investigate properties of carbon–oxygen white dwarfs with respect to the composite uncertainties in the reaction rates using the stellar evolution toolkit, Modules for Experiments in Stellar Astrophysics (MESA) and the probability density functions in the reaction rate library STARLIB. These are the first Monte Carlo stellar evolution studies that use complete stellar models. Focusing on 3 M{sub ⊙} models evolved from the pre main-sequence to the first thermal pulse, we survey the remnant core mass, composition, and structure properties as a function of 26 STARLIB reaction rates covering hydrogen and helium burning using a Principal Component Analysis and Spearman Rank-Order Correlation. Relative to the arithmetic mean value, we find the width of the 95% confidence interval to be ΔM{sub 1TP} ≈ 0.019 M{sub ⊙} for the core mass at the first thermal pulse, Δt{sub 1TP} ≈ 12.50 Myr for the age, Δlog(T{sub c}/K) ≈ 0.013 for the central temperature, Δlog(ρ{sub c}/g cm{sup −3}) ≈ 0.060 for the central density, ΔY{sub e,c} ≈ 2.6 × 10{sup −5} for the central electron fraction, ΔX{sub c}({sup 22}Ne) ≈ 5.8 × 10{sup −4}, ΔX{sub c}({sup 12}C) ≈ 0.392, and ΔX{sub c}({sup 16}O) ≈ 0.392. Uncertainties in the experimental {sup 12}C(α,γ){sup 16}O, triple-α, and {sup 14}N(p,γ){sup 15}O reaction rates dominate these variations. We also consider a grid of 1–6 M{sub ⊙} models evolved from the pre main-sequence to the final white dwarf to probe the sensitivity of the initial–final mass relation to experimental uncertainties in the hydrogen and helium reaction rates.

79 ASTRONOMY AND ASTROPHYSICS↗

Multi-frequency spin manipulation using rapidly tunable superconducting coplanar waveguide microresonators

In this work, we demonstrate the use of frequency-tunable superconducting NbTiN coplanar waveguide microresonators for multi-frequency pulsed electron spin resonance (ESR) experiments. By applying a bias current to the center pin, the resonance frequency (~7.6 GHz) can be continuously tuned by as much as 95 MHz in 270 ns without a change in the quality factor of 3000 at 2 K. We demonstrate the ESR performance of our resonators by measuring donor spin ensembles in silicon and show that adiabatic pulses can be used to overcome magnetic field inhomogeneities and microwave power limitations due to the applied bias current. We take advantage of the rapid tunability of these resonators to manipulate both phosphorus and arsenic spins in a single pulse sequence, demonstrating pulsed double electron-electron resonance. Our NbTiN resonator design is useful for multi-frequency pulsed ESR and should also have applications in experiments where spin ensembles are used as quantum memories.

42 ENGINEERING↗

Proton-detected 3D {sup 1}H/{sup 13}C/{sup 1}H correlation experiment for structural analysis in rigid solids under ultrafast-MAS above 60 kHz

A proton-detected 3D {sup 1}H/{sup 13}C/{sup 1}H chemical shift correlation experiment is proposed for the assignment of chemical shift resonances, identification of {sup 13}C-{sup 1}H connectivities, and proximities of {sup 13}C-{sup 1}H and {sup 1}H-{sup 1}H nuclei under ultrafast magic-angle-spinning (ultrafast-MAS) conditions. Ultrafast-MAS is used to suppress all anisotropic interactions including {sup 1}H-{sup 1}H dipolar couplings, while the finite-pulse radio frequency driven dipolar recoupling (fp-RFDR) pulse sequence is used to recouple dipolar couplings among protons and the insensitive nuclei enhanced by polarization transfer technique is used to transfer magnetization between heteronuclear spins. The 3D experiment eliminates signals from non-carbon-bonded protons and non-proton-bonded carbons to enhance spectral resolution. The 2D (F1/F3) {sup 1}H/{sup 1}H and 2D {sup 13}C/{sup 1}H (F2/F3) chemical shift correlation spectra extracted from the 3D spectrum enable the identification of {sup 1}H-{sup 1}H proximity and {sup 13}C-{sup 1}H connectivity. In addition, the 2D (F1/F2) {sup 1}H/{sup 13}C chemical shift correlation spectrum, incorporated with proton magnetization exchange via the fp-RFDR recoupling of {sup 1}H-{sup 1}H dipolar couplings, enables the measurement of proximities between {sup 13}C and even the remote non-carbon-bonded protons. The 3D experiment also gives three-spin proximities of {sup 1}H-{sup 1}H-{sup 13}C chains. Experimental results obtained from powder samples of L-alanine and L-histidine ⋅ H{sub 2}O ⋅ HCl demonstrate the efficiency of the 3D experiment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Microbial legacies alter decomposition in response to simulated global change

Terrestrial ecosystem models assume that microbial communities respond instantaneously, or are immediately resilient, to environmental change. Here we tested this assumption by quantifying the resilience of a leaf litter community to changes in precipitation or nitrogen availability. By manipulating composition within a global change experiment, we decoupled the legacies of abiotic parameters versus that of the microbial community itself. After one rainy season, more variation in fungal composition could be explained by the original microbial inoculum than the litterbag environment (18% versus 5.5% of total variation). This compositional legacy persisted for 3 years, when 6% of the variability in fungal composition was still explained by the microbial origin. In contrast, bacterial composition was generally more resilient than fungal composition. Microbial functioning (measured as decomposition rate) was not immediately resilient to the global change manipulations; decomposition depended on both the contemporary environment and rainfall the year prior. Finally, using metagenomic sequencing, we showed that changes in precipitation, but not nitrogen availability, altered the potential for bacterial carbohydrate degradation, suggesting why the functional consequences of the two experiments may have differed. Predictions of how terrestrial ecosystem processes respond to environmental change may thus be improved by considering the legacies of microbial communities.

54 ENVIRONMENTAL SCIENCES↗

Fungal and bacterial growth variation due to drought and nitrogen addition experimental treatments. Loma Ridge Experimental Project. 2010-2012

Terrestrial ecosystem models assume that microbial communities respond instantaneously, or are immediately resilient, to environmental change. Here we tested this assumption by quantifying the resilience of a leaf litter community to changes in precipitation or nitrogen availability. By manipulating composition within a global change experiment, we decoupled the legacies of abiotic parameters versus that of the microbial community itself. After one rainy season, more variation in fungal composition could be explained by the original microbial inoculum than the litterbag environment (18% versus 5.5% of total variation). This compositional legacy persisted for 3 years, when 6% of the variability in fungal composition was still explained by the microbial origin. In contrast, bacterial composition was generally more resilient than fungal composition. Microbial functioning (measured as decomposition rate) was not immediately resilient to the global change manipulations; decomposition depended on both the contemporary environment and rainfall the year prior. Finally, using metagenomic sequencing, we showed that changes in precipitation, but not nitrogen availability, altered the potential for bacterial carbohydrate degradation, suggesting why the functional consequences of the two experiments may have differed. Predictions of how terrestrial ecosystem processes respond to environmental change may thus be improved by considering the legacies of microbial communities. This data package includes ten csv files (five data files and their corresponding data dictionaries) and one file-level metadata excel file. Data files contains information about which plots were exposed to treatments related to drought and nitrogen, information about litter bags reciprocal transplants manipulation for water input and nitrogen, detail information about water addition, precipitation records, and litter variables collected. Data dictionary files include detail explanation for each column in the data files. The file-level metadata file describes each file mentioned above. All the analyses were done using the R software.

54 ENVIRONMENTAL SCIENCES↗

A Deep Learning-Aided Workflow for Decoding the Stress Regime of Southern Nevada

The Rock Valley fault zone in southern Nevada has a notable history of seismic activity and is the site of a future direct comparison experiment of explosion and earthquake sources. This study aims to gain insight into regional tectonic processes by leveraging recent advances in seismic monitoring capabilities to elucidate the local stress regime. A crucial step in this investigation is the accurate determination of P-wave first-motion polarities, which play a vital role in resolving earthquake focal mechanisms of small earthquakes. Here, we deploy a deep learning-based method for automatic determination of first-motion polarities to vastly expand the polarity dataset beyond what has been reviewed by human analysts. By the integrating P-wave polarities with new measurements of S/P amplitude ratios, we obtain robust focal mechanism estimates for 1306 earthquakes with a local magnitude of 1 and above occurring between 2010 and 2023 in southern Nevada. We then use the focal mechanism catalog to examine the regional stress orientation, confirming an overall trans-tensional stress regime with smaller scale complexities illuminated by individual earthquake sequences. These findings demonstrate how detailed analyses of small earthquakes can provide fundamental information for understanding earthquake processes in the region and inform future experiments at the Nevada National Security Site.

58 GEOSCIENCES↗

Gate Set Tomography

Gate set tomography (GST) is a protocol for detailed, predictive characterization of logic operations (gates) on quantum computing processors. Early versions of GST emerged around 2012-13, and since then it has been refined, demonstrated, and used in a large number of experiments. This paper presents the foundations of GST in comprehensive detail. The most important feature of GST, compared to older state and process tomography protocols, is that it is calibration-free . GST does not rely on pre-calibrated state preparations and measurements. Instead, it characterizes all the operations in a gate set simultaneously and self-consistently, relative to each other. Long sequence GST can estimate gates with very high precision and efficiency, achieving Heisenberg scaling in regimes of practical interest. In this paper, we cover GST's intellectual history, the techniques and experiments used to achieve its intended purpose, data analysis, gauge freedom and fixing, error bars, and the interpretation of gauge-fixed estimates of gate sets. Our focus is fundamental mathematical aspects of GST, rather than implementation details, but we touch on some of the foundational algorithmic tricks used in the pyGSTi implementation.

97 MATHEMATICS AND COMPUTING↗

Hardware In the Loop for Demand Flexibility (HIL4DF) v1.0

The software package in question is a collection of simulation models in the Modelica language, representing a variety of mechanical system designs and envelope conditions related to LBL's FLEXLAB facility. The collection of models also features multiple controls sequences that can be simulated with the FLEXLAB model to simulate different demand flexibility scenarios. Additionally, this package will feature datasets from 3 experimental tests, used for calibration, validation and comparison against the Modelica models, this includes weather data that can be used to replicate different scenarios in simulation across the same weather conditions experienced in real experiments. Given FLEXLAB high level of instrumentation and available data, the models are calibrated across multiple measurement points, and thus results from the extension of this model to other climate zones or control sequences, would provide high level of confidence.

Huang, Weiping↗

Warmer incubation temperatures and later lay–orders lead to shorter telomere lengths in wood duck ( Aix sponsa ) ducklings

The environment that animals experience during development shapes phenotypic expression. In birds, two important aspects of the early-developmental environment are lay-order sequence and incubation. Later-laid eggs tend to produce weaker offspring, sometimes with compensatory mechanisms to accelerate their growth rate to catch-up to their siblings. Further, small decreases in incubation temperature slow down embryonic growth rates and lead to wide-ranging negative effects on many post-hatch traits. Recently, telomeres, non-coding DNA sequences at the end of chromosomes, have been recognized as a potential proxy for fitness because longer telomeres are positively related to lifespan and individual quality in many animals, including birds. Although telomeres appear to be mechanistically linked to growth rate, little is known about how incubation temperature and lay-order may influence telomere length. We incubated wood duck (Aix sponsa) eggs at two ecologically-relevant temperatures (34.9 and 36.2ºC) and measured telomere length at hatch and one week after. We found that ducklings incubated at the lower temperature had longer telomeres than those incubated at the higher temperature both at hatch and one week later. Further, we found that later-laid eggs produced ducklings with shorter telomeres than those laid early in the lay-sequence, although lay-order was not related to embryonic developmental rate. Furthermore, this study contributes to our broader understanding of how parental effects can affect telomere length early in life. More work is needed to determine if these effects on telomere length persist until adulthood, and if they are associated with effects on fitness in this precocial species.

59 BASIC BIOLOGICAL SCIENCES↗