Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mutation testing”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Characterization of a nuclear localization signal in the foot-and-mouth disease virus polymerase

We have experimentally tested whether the MRKTKLAPT sequence in FMDV 3D protein (residues 16 to 24) can act as a nuclear localization signal (NLS). Mutants with substitutions in two basic residues within this sequence, K18E and K20E, were generated. A decreased nuclear localization was observed in transiently expressed 3D and its precursor 3CD, suggesting a role of K18 and K20 in nuclear targeting. Fusion of MRKTKLAPT to the green fluorescence protein (GFP) increased the nuclear localization of GFP, which was not observed when GFP was fused to the 3D mutated sequences. These results indicate that the sequence MRKTKLAPT can be functionally considered as a NLS. When introduced in a FMDV full length RNA replacements K18E and K20E led to production of revertant viruses that replaced the acidic residues introduced (E) by K, suggesting that the presence of lysins at positions 18 and 20 of 3D is essential for virus multiplication. - Highlights: • The FMDV 3D polymerase contains a nuclear localization signal. • Replacements K18E and K20E decrease nuclear localization of 3D and its precursor 3CD. • Fusion of the MRKTKLAPT 3D motif to GFP increases the nuclear localization of GFP. • Replacements K18E and K20E abolish the ability of MRKTKLAPT to relocate GFP. • RNAs harboring replacements K18E and K20E lead to recovery of revertant FMDVs.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Empirical evidence that glucan-interacting amino acid side chains within the transmembrane channel collectively facilitate cellulose synthase function

The fundamental mechanism of cellulose synthesis is widely conserved across Kingdoms and depends on cellulose synthases, which are processive, dual-function, family 2 glycosyltransferases (GT-2). These enzymes polymerize glucose on the cytoplasmic side of the plasma membrane and export the glucan chain to the cell surface through an integral transmembrane (TM) channel. Structural studies of active plant cellulose synthases (CESAs) have revealed interactions between the nascent glucan chain and the side chains of polar, charged, and aromatic amino acid residues that line the TM channel. However, the functional consequences of modifying these side chains have not been tested in vivo in CESAs or other processive GT-2s. To test this, we used an established in vivo assay based on genetic complementation of CESA5 in the moss, Physcomitrium patens. For accurate prediction of glucan-interacting amino acid residues, we generated a complete homotrimeric molecular model of PpCESA5 using a combination of homology and de novo modeling. All-atom molecular dynamics-based analyses of contact metrics and interaction energy identified 23 amino acid residues with high propensity to interact with the nascent glucan chain within the TM channel or on the apoplastic surface of PpCESA5. Mutating any one of 18 of these amino acid residues to alanine, thereby removing their side chains, abolished or impaired CESA function, with the strongest effects observed upon the loss of charged amino acid side chains. This provides direct evidence to support the hypothesis that multiple amino acid residues collectively maintain a smooth energy landscape within the TM channel to facilitate glucan translocation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

Associations of ATM Polymorphisms With Survival in Advanced Esophageal Squamous Cell Carcinoma Patients Receiving Radiation Therapy

Purpose: To investigate whether single nucleotide polymorphisms (SNPs) in the ataxia telangiectasia mutated (ATM) gene are associated with survival in patients with esophageal squamous cell carcinoma (ESCC) receiving radiation therapy or chemoradiation therapy or surgery only. Methods and Materials: Four tagSNPs of ATM were genotyped in 412 individuals with clinical stage III or IV ESCC receiving radiation therapy or chemoradiation therapy, and in 388 individuals with stage I, II, or III ESCC treated with surgery only. Overall survival time of ESCC among different genotypes was estimated by Kaplan-Meier plot, and the significance was examined by log-rank test. The hazard ratios (HRs) and 95% confidence intervals (CIs) for death from ESCC among different genotypes were computed by a Cox proportional regression model. Results: We found 2 SNPs, rs664143 and rs664677, associated with survival time of ESCC patients receiving radiation therapy. Individuals with the rs664143A allele had poorer median survival time compared with the rs664143G allele (14.0 vs 20.0 months), with the HR for death being 1.45 (95% CI 1.12-1.89). Individuals with the rs664677C allele also had worse median survival time than those with the rs664677T allele (14.0 vs 23.5 months), with the HR of 1.57 (95% CI 1.18-2.08). Stratified analysis showed that these associations were present in both stage III and IV cancer and different radiation therapy techniques. Significant associations were also found between the SNPs and locosregional progression or progression-free survival. No association between these SNPs and survival time was detected in ESCC patients treated with surgery only. Conclusion: These results suggest that the ATM polymorphisms might serve as independent biomarkers for predicting prognosis in ESCC patients receiving radiation therapy.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2015

Increasing Sequence Diversity with Flexible Backbone Protein Design: The Complete Redesign of a Protein Hydrophobic Core

Protein design tests our understanding of protein stability and structure. Successful design methods should allow the exploration of sequence space not found in nature. However, when redesigning naturally occurring protein structures, most fixed backbone design algorithms return amino acid sequences that share strong sequence identity with wild-type sequences, especially in the protein core. This behavior places a restriction on functional space that can be explored and is not consistent with observations from nature, where sequences of low identity have similar structures. Here, we allow backbone flexibility during design to mutate every position in the core (38 residues) of a four-helix bundle protein. Only small perturbations to the backbone, 12 {angstrom}, were needed to entirely mutate the core. The redesigned protein, DRNN, is exceptionally stable (melting point >140C). An NMR and X-ray crystal structure show that the side chains and backbone were accurately modeled (all-atom RMSD = 1.3 {angstrom}).

Murphy, Grant S.↗
DOE OSTI2020

New candidates for regulated gene integrity revealed through precise mapping of integrative genetic elements

Abstract Integrative genetic elements (IGEs) are mobile multigene DNA units that integrate into and excise from host bacterial genomes. Each IGE usually targets a specific site within a conserved host gene, integrating in a manner that preserves target gene function. However, a small number of bacterial genes are known to be inactivated upon IGE integration and reactivated upon excision, regulating phenotypes of virulence, mutation rate, and terminal differentiation in multicellular bacteria. The list of regulated gene integrity (RGI) cases has been slow-growing because IGEs have been challenging to precisely and comprehensively locate in genomes. We present software (TIGER) that maps IGEs with unprecedented precision and without attB site bias. TIGER uses a comparative genomic, ping-pong BLAST approach, based on the principle that the IGE integration module (i.e. its int-attP region) is cohesive. The resultant IGEs from 2168 genomes, along with integrase phylogenetic analysis and gene inactivation tests, revealed 19 new cases of genes whose integrity is regulated by IGEs (including dut, eccCa1, gntT, hrpB, merA, ompN, prkA, tqsA, traG, yifB, yfaT and ynfE), as well as recovering previously known cases (in sigK, spsM, comK, mlrA and hlb genes). It also recovered known clades of site-promiscuous integrases and identified possible new ones.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2014

Comprehensive Detection of Genes Causing a Phenotype Using Phenotype Sequencing and Pathway Analysis

Discovering all the genetic causes of a phenotype is an important goal in functional genomics. We combine an experimental design for detecting independent genetic causes of a phenotype with a high-throughput sequencing analysis that maximizes sensitivity for comprehensively identifying them. Testing this approach on a set of 24 mutant strains generated for a metabolic phenotype with many known genetic causes, we show that this pathway-based phenotype sequencing analysis greatly improves sensitivity of detection compared with previous methods, and reveals a wide range of pathways that can cause this phenotype. We demonstrate our approach on a metabolic re-engineering phenotype, the PEP/OAA metabolic node in E. coli, which is crucial to a substantial number of metabolic pathways and under renewed interest for biofuel research. Out of 2157 mutations in these strains, pathway-phenoseq discriminated just five gene groups (12 genes) as statistically significant causes of the phenotype. Experimentally, these five gene groups, and the next two high-scoring pathway-phenoseq groups, either have a clear connection to the PEP metabolite level or offer an alternative path of producing oxaloacetate (OAA), and thus clearly explain the phenotype. These high-scoring gene groups also show strong evidence of positive selection pressure, compared with strictly neutral selection in the rest of the genome.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

RhoA and vigilin are candidates for immunohistochemical markers for epithelioid malignant mesothelioma

Diagnostic markers of malignant mesothelioma (MM) have been extensively investigated. Immunohistochemistry (IHC) markers, such as calretinin, have been used for pathologic diagnosis. However, more diagnostic markers are required to improve the specificity and sensitivity of pathologic diagnosis. This study proposed two proteins as diagnostic markers for epithelioid MM. One is RhoA, an MM mutation-susceptible locus-derived protein, and another is vigilin, a lung small cell carcinoma marker. IHC was performed using 93 MM (epithelioid, 71 cases; sarcomatoid, 13 cases; and biphasic, 9 cases), 64 lung adenocarcinoma (LAC), 60 lung squamous cell carcinoma (LSC), and 14 normal mesothelial (NM) tissues. The majority of epithelioid MM cases were positive for both RhoA and vigilin, whereas both IHCs showed lower stainability in biphasic and sarcomatoid MM. Besides, both IHCs showed significantly higher stainability for RhoA and vigilin in epithelioid MM than in LAC and LSC (p < 0.05). Chi-square tests showed that both RhoA and vigilin IHC positive rate in epithelioid MM was not significantly different from that of calretinin (p > 0.05). In the differential diagnosis of MM from lung cancer, the accuracy and specificity of RhoA, vigilin, and calretinin staining were almost equivalent. Further, H-score test showed that there was no significant difference between RhoA versus calretinin and vigilin versus calretinin in IHC positivity in epithelioid MM (p > 0.05). In conclusion, RhoA and vigilin may be candidates for immunohistochemical markers for epithelioid MM.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

On the Rapid Calculation of Binding Affinities for Antigen and Antibody Design and Affinity Maturation Simulations

The accurate and efficient calculation of protein-protein binding affinities is an essential component in antibody and antigen design and optimization, and in computer modeling of antibody affinity maturation. Such calculations remain challenging despite advances in computer hardware and algorithms, primarily because proteins are flexible molecules, and thus, require explicit or implicit incorporation of multiple conformational states into the computational procedure. The astronomical size of the amino acid sequence space further compounds the challenge by requiring predictions to be computed within a short time so that many sequence variants can be tested. In this study, we compare three classes of methods for antibody/antigen (Ab/Ag) binding affinity calculations: (i) a method that relies on the physical separation of the Ab/Ag complex in equilibrium molecular dynamics (MD) simulations, (ii) a collection of 18 scoring functions that act on an ensemble of structures created using homology modeling software, and (iii) methods based on the molecular mechanics-generalized Born surface area (MM-GBSA) energy decomposition, in which the individual contributions of the energy terms are scaled to optimize agreement with the experiment. When applied to a set of 49 antibody mutations in two Ab/HIV gp120 complexes, all of the methods are found to have modest accuracy, with the highest Pearson correlations reaching about 0.6. In particular, the most computationally intensive method, i.e., MD simulation, did not outperform several scoring functions. The optimized energy decomposition methods provided marginally higher accuracy, but at the expense of requiring experimental data for parametrization. Within each method class, we examined the effect of the number of independent computational replicates, i.e., modeled structures or reinitialized MD simulations, on the prediction accuracy. We suggest using about ten modeled structures for scoring methods, and about five simulation replicates for MD simulations as a rule of thumb for obtaining reasonable convergence. We anticipate that our study will be a useful resource for practitioners working to incorporate binding affinity calculations within their protein design and optimization process.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2013

Conserved structural chemistry for incision activity in structurally non-homologous apurinic/apyrimidinic endonuclease APE1 and endonuclease IV DNA repair enzymes.

Non-coding apurinic/apyrimidinic (AP) sites in DNA form spontaneously and as DNA base excision repair intermediates are the most common toxic and mutagenic in vivo DNA lesion. For repair, AP sites must be processed by 5' AP endonucleases in initial stages of base repair. Human APE1 and bacterial Nfo represent the two conserved 5' AP endonuclease families in the biosphere; they both recognize AP sites and incise the phosphodiester backbone 5' to the lesion, yet they lack similar structures and metal ion requirements. Here, we determined and analyzed crystal structures of a 2.4 ? resolution APE1-DNA product complex with Mg(2+) and a 0.92 Nfo with three metal ions. Structural and biochemical comparisons of these two evolutionarily distinct enzymes characterize key APE1 catalytic residues that are potentially functionally similar to Nfo active site components, as further tested and supported by computational analyses. We observe a magnesium-water cluster in the APE1 active site, with only Glu-96 forming the direct protein coordination to the Mg(2+). Despite differences in structure and metal requirements of APE1 and Nfo, comparison of their active site structures surprisingly reveals strong geometric conservation of the catalytic reaction, with APE1 catalytic side chains positioned analogously to Nfo metal positions, suggesting surprising functional equivalence between Nfo metal ions and APE1 residues. The finding that APE1 residues are positioned to substitute for Nfo metal ions is supported by the impact of mutations on activity. Collectively, the results illuminate the activities of residues, metal ions, and active site features for abasic site endonucleases.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2014

Targeted Analysis of Whole Genome Sequence Data to Diagnose Genetic Cardiomyopathy

Background—Cardiomyopathy is highly heritable but genetically diverse. At present, genetic testing for cardiomyopathy uses targeted sequencing to simultaneously assess the coding regions of more than 50 genes. New genes are routinely added to panels to improve the diagnostic yield. With the anticipated $1000 genome, it is expected that genetic testing will shift towards comprehensive genome sequencing accompanied by targeted gene analysis. Therefore, we assessed the reliability of whole genome sequencing and targeted analysis to identify cardiomyopathy variants in 11 subjects with cardiomyopathy. Methods and Results—Whole genome sequencing with an average of 37× coverage was combined with targeted analysis focused on 204 genes linked to cardiomyopathy. Genetic variants were scored using multiple prediction algorithms combined with frequency data from public databases. This pipeline yielded 1-14 potentially pathogenic variants per individual. Variants were further analyzed using clinical criteria and/or segregation analysis. Three of three previously identified primary mutations were detected by this analysis. In six subjects for whom the primary mutation was previously unknown, we identified mutations that segregated with disease, had clinical correlates, and/or had additional pathological correlation to provide evidence for causality. For two subjects with previously known primary mutations, we identified additional variants that may act as modifiers of disease severity. In total, we identified the likely pathological mutation in 9 of 11 (82%) subjects. We conclude that these pilot data demonstrate that ~30-40× coverage whole genome sequencing combined with targeted analysis is feasible and sensitive to identify rare variants in cardiomyopathy-associated genes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

classLog: Logistic regression for the classification of genetic sequences

Introduction Sequencing and phylogenetic classification have become a common task in human and animal diagnostic laboratories. It is routine to sequence pathogens to identify genetic variations of diagnostic significance and to use these data in realtime genomic contact tracing and surveillance. Under this paradigm, unprecedented volumes of data are generated that require rapid analysis to provide meaningful inference. Methods We present a machine learning logistic regression pipeline that can assign classifications to genetic sequence data. The pipeline implements an intuitive and customizable approach to developing a trained prediction model that runs in linear time complexity, generating accurate output rapidly, even with incomplete data. Our approach was benchmarked against porcine respiratory and reproductive syndrome virus (PRRSv) and swine H1 influenza A virus (IAV) datasets. Trained classifiers were tested against sequences and simulated datasets that artificially degraded sequence quality at 0, 10, 20, 30, and 40%. Results When applied to a poor-quality sequence data, the classifier achieved between >85% to 95% accuracy for the PRRSv and the swine H1 IAV HA dataset and this increased to near perfect accuracy when using the full dataset. The model also identifies amino acid positions used to determine genetic clade identity through a feature selection ranking within the model. These positions can be mapped onto a maximum-likelihood phylogenetic tree, allowing for the inference of clade defining mutations. Discussion Our approach is implemented as a python package with code available at https://github.com/flu-crew/classLog .

Zeller, Michael A.↗
DOE OSTI2016

Tunable allosteric library of caspase-3 identifies coupling between conserved water molecules and conformational selection

The native ensemble of caspases is described globally by a complex energy landscape where the binding of substrate selects for the active conformation, whereas targeting an allosteric site in the dimer interface selects an inactive conformation that contains disordered active-site loops. Mutations and posttranslational modifications stabilize high-energy inactive conformations, with mostly formed, but distorted, active sites. To examine the interconversion of active and inactive states in the ensemble, we used detection of related solvent positions to analyze 4,995 waters in 15 high-resolution (<2.0 Å) structures of wild-type caspase-3, resulting in 450 clusters with the most highly conserved set containing 145 water molecules. The data show that regions of the protein that contact the conserved waters also correspond to sites of posttranslational modifications, suggesting that the conserved waters are an integral part of allosteric mechanisms. To test this hypothesis, we created a library of 19 caspase-3 variants through saturation mutagenesis in a single position of the allosteric site of the dimer interface, and we show that the enzyme activity varies by more than four orders of magnitude. Altogether, our database consists of 37 high-resolution structures of caspase-3 variants, and we demonstrate that the decrease in activity correlates with a loss of conserved water molecules. Furthermore, the data show that the activity of caspase-3 can be fine-tuned through globally desolvating the active conformation within the native ensemble, providing a mechanism for cells to repartition the ensemble and thus fine-tune activity through conformational selection.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2018

Sclerostin neutralization unleashes the osteoanabolic effects of Dkk1 inhibition [Genetic disruption of compensatory WNT inhibitor expression reveals a context-dependent, highly osteoanabolic role for DKK1 inhibition in the skeleton]

The WNT pathway has become an attractive target for skeletal therapies. High-bone-mass phenotypes in patients with loss-of-function mutations in the LRP5/6 inhibitor Sost (sclerosteosis), or in its downstream enhancer region (van Buchem disease), highlight the utility of targeting Sost/sclerostin to improve bone properties. Sclerostin-neutralizing antibody is highly osteoanabolic in animal models and in human clinical trials, but antibody-based inhibition of another potent LRP5/6 antagonist, Dkk1, is largely inefficacious for building bone in the unperturbed adult skeleton. In this work, we show that conditional deletion of Dkk1 from bone also has negligible effects on bone mass. Dkk1 inhibition increases Sost expression, suggesting a potential compensatory mechanism that might explain why Dkk1 suppression lacks anabolic action. To test this concept, we deleted Sost from osteocytes in, or administered sclerostin neutralizing antibody to, mice with a Dkk1-deficient skeleton. A robust anabolic response to Dkk1 deletion was manifest only when Sost/sclerostin was impaired. Whole-body DXA scans, μCT measurements of the femur and spine, histomorphometric measures of femoral bone formation rates, and biomechanical properties of whole bones confirmed the anabolic potential of Dkk1 inhibition in the absence of sclerostin. Moreover, combined administration of sclerostin and Dkk1 antibody in WT mice produced a synergistic effect on bone gain that greatly exceeded individual or additive effects of the therapies, confirming the therapeutic potential of inhibiting multiple WNT antagonists for skeletal health. In conclusion, the osteoanabolic effects of Dkk1 inhibition can be realized if sclerostin upregulation is prevented. Anabolic therapies for patients with low bone mass might benefit from a strategy that accounts for the compensatory milieu of WNT inhibitors in bone tissue.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Potential for Applying Continuous Directed Evolution to Plant Enzymes: An Exploratory Study

Plant evolution has produced enzymes that may not be optimal for maximizing yield and quality in today’s agricultural environments and plant biotechnology applications. By improving enzyme performance, it should be possible to alleviate constraints on yield and quality currently imposed by kinetic properties or enzyme instability. Enzymes can be optimized more quickly than naturally possible by applying directed evolution, which entails mutating a target gene in vitro and screening or selecting the mutated gene products for the desired characteristics. Continuous directed evolution is a more efficient and scalable version that accomplishes the mutagenesis and selection steps simultaneously in vivo via error-prone replication of the target gene and coupling of the host cell’s growth rate to the target gene’s function. However, published continuous systems require custom plasmid assembly, and convenient multipurpose platforms are not available. We discuss two systems suitable for continuous directed evolution of enzymes, OrthoRep in Saccharomyces cerevisiae and EvolvR in Escherichia coli, and our pilot efforts to adapt each system for high-throughput plant enzyme engineering. To test our modified systems, we used the thiamin synthesis enzyme THI4, previously identified as a prime candidate for improvement. Our adapted OrthoRep system shows promise for efficient plant enzyme engineering.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Binding of adenovirus species C hexon to prothrombin and the influence of hexon on vector properties in vitro and in vivo

The majority of adenovirus (Ad) vectors are based on human Ad type 5, which is a member of Ad species C. Species C also includes the closely-related types 1, 2, 6, 57 and 89. It is known that coagulation factors bind to Ad5 hexon and play a key role in the liver tropism of Ad5 vectors, but it is unclear how coagulation factors affect vectors derived from other species C Ads. We evaluated species C Ad vectors both in vitro and following intravenous injection in mice. To assess the impact of hexon differences, we constructed chimeric Ad5 vectors that contain the hexon hypervariable regions from other species C types, including vectors with hexon mutations that decreased coagulation factor binding. After intravenous injection into mice, vectors with Ad5 or Ad6 hexon had strong liver tropism, while vectors with chimeric hexon from other Ad types had weaker liver tropism due to inhibition by natural antibodies and complement. In addition, we discovered a novel ability of hexon to bind prothrombin, which is the most abundant coagulation factor in blood, and we found striking differences in the affinity of Ads for human, mouse and bovine coagulation factors. When compared to Ad5, vectors with non-Ad5 species C hexons had considerably higher affinity for both human and mouse prothrombin. Most of the vectors tested were strongly dependent on coagulation factors for liver transduction, but vectors with chimeric Ad6 hexon showed much less dependence on coagulation factors than other vectors. We found that in vitro neutralization experiments with mouse serum predicted in vivo behavior of Ad5 vectors, but in vitro experiments did not predict the in vivo behavior of vectors based on other Ad types. In sum, hexons from different human Ad species C viruses confer diverse properties on vectors, including differing abilities to target the liver.

Tian, Jie↗
DOE OSTI2014

Hybrid fusions show that inter-monomer electron transfer robustly supports cytochrome bc{sub 1} function in vivo

Highlights: • We used hybrid fusion bc{sub 1} complex to test inter-monomer electron transfer in vivo. • Cross-inactivated complexes were able to sustain photoheterotrophic growth. • Inter-monomer electron transfer supports catalytic cycle in vivo. • bc{sub 1} dimer is functional even when cytochrome b subunits come from different species. - Abstract: Electronic connection between Q{sub o} and Q{sub i} quinone catalytic sites of dimeric cytochrome bc{sub 1} is a central feature of the energy-conserving Q cycle. While both the intra- and inter-monomer electron transfers were shown to connect the sites in the enzyme, mechanistic and physiological significance of the latter remains unclear. Here, using a series of mutated hybrid cytochrome bc{sub 1}-like complexes, we show that inter-monomer electron transfer robustly sustains the function of the enzyme in vivo, even when the two subunits in a dimer come from different species. This indicates that minimal requirement for bioenergetic efficiency is to provide a chain of cofactors for uncompromised electron flux between the catalytic sites, while the details of protein scaffold are secondary.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Estimating the Timing of Early Simian-Human Immunodeficiency Virus Infections: a Comparison between Poisson Fitter and BEAST

Many HIV prevention strategies are currently under consideration where it is highly informative to know the study participants’ times of infection. These can be estimated using viral sequence data sampled early in infection. However, there are several scenarios that, if not addressed, can skew timing estimates. These include multiple transmitted/founder (TF) viruses, APOBEC (apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like)-mediated mutational enrichment, and recombination. Here, we suggest a pipeline to identify these problems and resolve the biases that they introduce. We then compare two modeling strategies to obtain timing estimates from sequence data. The first, Poisson Fitter (PF), is based on a Poisson model of random accumulation of mutations relative to the TF virus (or viruses) that established the infection. The second uses a coalescence-based phylogenetic strategy as implemented in BEAST. The comparison is based on timing predictions using plasma viral RNA (cDNA) sequence data from 28 simian-human immunodeficiency virus (SHIV)-infected animals for which the exact day of infection is known. In this particular setting, based on nucleotide sequences from samples obtained in early infection, the Poisson method yielded more accurate, more precise, and unbiased estimates for the time of infection than did the explored implementations of BEAST. The inference of the time of infection is a critical parameter in testing the efficacy of clinical interventions in protecting against HIV-1 infection. For example, in clinical trials evaluating the efficacy of passively delivered antibodies (Abs) for preventing infections, accurate time of infection data are essential for discerning levels of the Abs required to confer protection, given the natural Ab decay rate in the human body. In such trials, genetic sequences from early in the infection are regularly sampled from study participants, generally prior to immune selection, when the viral population is still expanding and genetic diversity is low. In this particular setting of early viral growth, the Poisson method is superior to the alternative approach based on coalescent methods. This approach can also be applied in human vaccine trials, where accurate estimates of infection times help ascertain if vaccine-elicited immune protection wanes over time.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2013

Methods for Performing Crosses in Setaria viridis , a New Model System for the Grasses

Setaria viridis is an emerging model system for C 4 grasses. It is closely related to the bioenergy feed stock switchgrass and the grain crop foxtail millet. Recently, the 510 Mb genome of foxtail millet, S. italica, has been sequenced 1,2 and a 25x coverage genome sequence of the weedy relative S. viridis is in progress. S. viridis has a number of characteristics that make it a potentially excellent model genetic system including a rapid generation time, small stature, simple growth requirements, prolific seed production 3 and developed systems for both transient and stable transformation 4 . However, the genetics of S. viridis is largely unexplored, in part, due to the lack of detailed methods for performing crosses. To date, no standard protocol has been adopted that will permit rapid production of seeds from controlled crosses. The protocol presented here is optimized for performing genetic crosses in S. viridis, accession A10.1. We have employed a simple heat treatment with warm water for emasculation after pruning the panicle to retain 20-30 florets and labeling of flowers to eliminate seeds resulting from newly developed flowers after emasculation. After testing a series of heat treatments at permissive temperatures and varying the duration of dipping, we have established an optimum temperature and time range of 48 °C for 3-6 min. By using this method, a minimum of 15 crosses can be performed by a single worker per day and an average of 3-5 outcross progeny per panicle can be recovered. Therefore, an average of 45-75 outcross progeny can be produced by one person in a single day. Broad implementation of this technique will facilitate the development of recombinant inbred line populations of S. viridis X S. viridis or S. viridis X S. italica, mapping mutations through bulk segregant analysis and creating higher order mutants for genetic analysis.

54 ENVIRONMENTAL SCIENCES↗