Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “genomic mutations”

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 253 records · Page 14

Using intrahost single nucleotide variant data to predict SARS-CoV-2 detection cycle threshold values

Over the last four years, each successive wave of the COVID-19 pandemic has been caused by variants with mutations that improve the transmissibility of the virus. Despite this, we still lack tools for predicting clinically important features of the virus. In this study, we show that it is possible to predict the PCR cycle threshold (Ct) values from clinical detection assays using sequence data. Ct values often correspond with patient viral load and the epidemiological trajectory of the pandemic. Using a collection of 36,335 high quality genomes, we built models from SARS-CoV-2 intrahost single nucleotide variant (iSNV) data, computing XGBoost models from the frequencies of A, T, G, C, insertions, and deletions at each position relative to the Wuhan-Hu-1 reference genome. Our best model had an R 2 of 0.604 [0.593–0.616, 95% confidence interval] and a Root Mean Square Error (RMSE) of 5.247 [5.156–5.337], demonstrating modest predictive power. Overall, we show that the results are stable relative to an external holdout set of genomes selected from SRA and are robust to patient status and the detection instruments that were used. This study highlights the importance of developing modeling strategies that can be applied to publicly available genome sequence data for use in disease prevention and control.

COVID19↗

Quantifying the local adaptive landscape of a nascent bacterial community

Abstract The fitness effects of all possible mutations available to an organism largely shape the dynamics of evolutionary adaptation. Yet, whether and how this adaptive landscape changes over evolutionary times, especially upon ecological diversification and changes in community composition, remains poorly understood. We sought to fill this gap by analyzing a stable community of two closely related ecotypes (“L” and “S”) shortly after they emerged within the E. coli Long-Term Evolution Experiment (LTEE). We engineered genome-wide barcoded transposon libraries to measure the invasion fitness effects of all possible gene knockouts in the coexisting strains as well as their ancestor, for many different, ecologically relevant conditions. We find consistent statistical patterns of fitness effect variation across both genetic background and community composition, despite the idiosyncratic behavior of individual knockouts. Additionally, fitness effects are correlated with evolutionary outcomes for a number of conditions, possibly revealing shifting patterns of adaptation. Together, our results reveal how ecological and epistatic effects combine to shape the adaptive landscape in a nascent ecological community.

59 BASIC BIOLOGICAL SCIENCES↗

Detecting macroevolutionary genotype–phenotype associations using error-corrected rates of protein convergence

On macroevolutionary timescales, extensive mutations and phylogenetic uncertainty mask the signals of genotype–phenotype associations underlying convergent evolution. To overcome this problem, we extended the widely used framework of non-synonymous to synonymous substitution rate ratios and developed the novel metric ω C , which measures the error-corrected convergence rate of protein evolution. While ω C distinguishes natural selection from genetic noise and phylogenetic errors in simulation and real examples, its accuracy allows an exploratory genome-wide search of adaptive molecular convergence without phenotypic hypothesis or candidate genes. Using gene expression data, we explored over 20 million branch combinations in vertebrate genes and identified the joint convergence of expression patterns and protein sequences with amino acid substitutions in functionally important sites, providing hypotheses on undiscovered phenotypes. We further extended our method with a heuristic algorithm to detect highly repetitive convergence among computationally non-trivial higher-order phylogenetic combinations. Our approach allows bidirectional searches for genotype–phenotype associations, even in lineages that diverged for hundreds of millions of years.

59 BASIC BIOLOGICAL SCIENCES↗

The genome of the polyextremophilic yeast, Naganishia friedmannii, reveals adaptations involved in stress response pathways, carbohydrate metabolism expansion, and a limited DNA repair repertoire

Here we report the draft genome sequence of Naganishia friedmannii (formerly Cryptococcus friedmannii) isolate, a Basidiomycota yeast commonly found in some of the most extreme environments of the Earth's cryosphere. We isolated N. friedmannii strain Llullensis from soils at 6000 m above sea level on Volcán Llullaillaco, Argentina. The genome was 22.2 Mb with 6251 identified protein coding genes. Proteins known to be associated with thermal, osmotic, and radiation stress were identified in the genome. Comparative analysis with seven other Naganishia genomes revealed unique features underlying its polyextremophilic lifestyle. Naganishia friedmannii showed an expansion of genes involved in breaking down plant-derived carbohydrates, supporting the hypothesis that it survives at high elevations by metabolizing wind-deposited organic matter. Surprisingly, many genes involved in cell-cycle checkpoints and DNA repair were missing, as in several other Naganishia species. This extensive loss may be adaptive in extreme environments prone to abiotic stress, where a high mutation rate could generate advantageous traits, and reduced cell-cycle control may allow for faster reproduction that would be advantageous for rapid growth during brief periods of soil wetting following rare snow events.

Vimercati, Lara↗

Distant-acting Enhancers in Development and Disease

Multiple converging lines of evidence from experimental and human genetic studies indicate that both common and rare sequence variants involving enhancers play major roles in Mendelian and complex human disease phenotypes. We use a combination of sequence-based molecular approaches including chromatin mapping and synthesis-enabled modification of enhancer sequences, large-scale transgenic mouse studies (http://enhancer.lbl.gov), and CRISPR genome editing in the mouse model to study the in vivofunction of enhancers in developmental, evolutionary, and disease-related processes. Using examples from our ongoing work, I will illustrate how these methods provide insight into the function and evolution of distant-acting regulatory sequences and offer a starting point for understanding their role in human disease. I will discuss how improved in vivo mouse assays can be used to study candidate pathogenic human mutations in enhancers at scale, revealing the extent to which human genetic studies alone can be insufficient to identify truly causal mutations in enhancers. I will also show how such assays can be used to explore the sensitivity of enhancers that are extremely conserved ("ultraconserved") in evolution to mutations, revealing a remarkable functional robustness to sequence changes that contrasts with apparent evolutionary constraints.

Visel, A↗

Identification of ancestral gnathostome Gli3 enhancers with activity in mammals

Abnormal expression of the transcriptional regulator and hedgehog (Hh) signaling pathway effector Gli3 is known to trigger congenital disease, most frequently affecting the central nervous system (CNS) and the limbs. Accurate delineation of the genomic cis-regulatory landscape controlling Gli3 transcription during embryonic development is critical for the interpretation of noncoding variants associated with congenital defects. Here, we employed a comparative genomic analysis on fish species with a slow rate of molecular evolution to identify seven previously unknown conserved noncoding elements (CNEs) in Gli3 intronic intervals (CNE15–21). Transgenic assays in zebrafish revealed that most of these elements drive activities in Gli3 expressing tissues, predominantly the fins, CNS, and the heart. Intersection of these CNEs with human disease associated SNPs identified CNE15 as a putative mammalian craniofacial enhancer, with conserved activity in vertebrates and potentially affected by mutation associated with human craniofacial morphology. Finally, comparative functional dissection of an appendage-specific CNE conserved in slowly evolving fish (elephant shark), but not in teleost (CNE14/hs1586) indicates co-option of limb specificity from other tissues prior to the divergence of amniotes and lobe-finned fish. In conclusion, these results uncover a novel subset of intronic Gli3 enhancers that arose in the common ancestor of gnathostomes and whose sequence components were likely gradually modified in other species during the process of evolutionary diversification.

59 BASIC BIOLOGICAL SCIENCES↗

Machine learning identifies novel signatures of antifungal drug resistance in Saccharomycotina yeasts

Antifungal drug resistance is a major challenge in fungal infection management. Numerous genomic changes are known to contribute to acquired drug resistance in clinical isolates of specific pathogens, but whether they broadly explain natural resistance across entire lineages is unknown. We leveraged genomic, ecological, and phenotypic trait data from naturally sampled strains from nearly all known species in subphylum Saccharomycotina to examine the evolution of resistance to eight antifungal drugs. The phylogenetic distribution of drug resistance varied by drug; fluconazole resistance was widespread, while 5-fluorocytosine resistance was rare, except in Lipomycetales. A random forest algorithm trained on genomic data predicted drug-resistant yeasts with 54–75% accuracy. Fluconazole resistance was consistently predicted with the highest accuracy (75.2%). Furthermore, fluconazole resistance prediction accuracy was similar between models trained on genome-wide variation in the presence and number of InterPro protein annotations across Saccharomycotina (75.2%) and those trained on amino acid sequence alignment data of Erg11, a protein known to be involved in fluconazole resistance (74.3-74.9%). Interestingly, the top Erg11 residues for predicting fluconazole resistance across Saccharomycotina do not overlap with, are not spatially close to, and are less conserved than those previously linked to resistance in clinical isolates of Candida albicans. In silico deep mutational scanning of the C. albicans Erg11 protein reveals that amino acid variants implicated in clinical cases of resistance are almost universally destabilizing while variants in our most informative residues are energetically more neutral, explaining why the latter are much more common than the former in natural populations. Importantly, previous experimental analyses of C. albicans Erg11 have shown that amino acid variation in our most informative residues, despite having never been directly implicated in clinical cases, can directly contribute to resistance. Our results suggest that studies of natural resistance in yeast species never encountered in the clinic will yield a fuller understanding of antifungal drug resistance.

Harrison, Marie-Claire [Vanderbilt Univ., Nashvill↗

Genomic Features of Antimicrobial Resistance in Staphylococcus pseudintermedius Isolated from Dogs with Pyoderma in Argentina and the United States: A Comparative Study

Staphylococcus pseudintermedius is the most common opportunistic pathogen in dogs and methicillin resistance (MRSP) has been identified as an emerging problem in canine pyoderma. Here, we evaluated the antimicrobial resistance (AMR) features and phylogeny of S. pseudintermedius isolated from canine pyoderma cases in Argentina (n = 29) and the United States (n = 29). 62% of isolates showed multi-drug resistance. The AMR genes found: mecA, blaZ, ermB, dfrG, catA, tetM, aac(6')-aph(2"), in addition to tetK and lnuA (only found in U.S. isolates). Two point mutations were detected: grlA(S80I)-gyrA(S84L), and grlA(D84N)-gyrA(S84L) in one U.S. isolate. A mutation in rpoB (H481N) was found in two isolates from Argentina. SCCmec type III, SCCmec type V, ΨSCCmec 57395 were identified in the Argentinian isolates; and SCCmec type III, SCCmec type IVg, SCCmec type V, and SCCmec type VII variant in the U.S. cohort. Sequence type (ST) ST71 belonging to a dominant clone was found in isolates from both countries, and ST45 only in Argentinian isolates. This is the first study to comparatively analyze the population structure of canine pyoderma-associated S. pseudintermedius isolates in Argentina and in the U.S. It is important to maintain surveillance on S. pseudintermedius populations to monitor AMR and gain further understanding of its evolution and dissemination.

60 APPLIED LIFE SCIENCES↗

Comparative Genomics Applied to Systematically Assess Pathogenicity Potential in Shiga Toxin-Producing Escherichia coli O145:H28

Shiga toxin-producing Escherichia coli (STEC) O145:H28 can cause severe disease in humans and is a predominant serotype in STEC O145 environmental isolates. Here, comparative genomics was applied to a set of clinical and environmental strains to systematically evaluate the pathogenicity potential in environmental strains. While the core genes-based tree separated all O145:H28 strains from the non O145:H28 reference strains, it failed to segregate environmental strains from the clinical. In contrast, the accessory genes-based tree placed all clinical strains in the same clade regardless of their genotypes or serotypes, apart from the environmental strains. Loss-of-function mutations were common in the virulence genes examined, with a high frequency in genes related to adherence, autotransporters, and the type three secretion system. Distinct differences in pathogenicity islands LEE, OI-122, and OI-57, the acid fitness island, and the tellurite resistance island were detected between the O145:H28 and reference strains. A great amount of genetic variation was detected in O145:H28, which was mainly attributed to deletions, insertions, and gene acquisition at several chromosomal “hot spots”. Finally, our study demonstrated a distinct virulence gene repertoire among the STEC O145:H28 strains originating from the same geographical region and revealed unforeseen contributions of loss-of-function mutations to virulence evolution and genetic diversification in STEC.

59 BASIC BIOLOGICAL SCIENCES↗

Whole genome resequencing data from a collection of Clostridium Thermocellum strains

Clostridium thermocellum is an anaerobic thermophilic bacterium that natively ferments cellulose to ethanol and organic acids. This data set is a collection of whole genome resequencing data for several hundred strains of Clostridium thermocellum. It includes strains that have been engineered to increase ethanol production, strains that have been engineered to understand microbial physiology, and strains that have been adapted for desired phenotypes including increased ethanol tolerance. Resequencing data consists of paired Illumina reads, 100-150 bp on each end, with a ~500 bp insert size. One data file containing raw Illumina data (interleaved) is available for each strain. We also provide data describing the mutations identified in each strain, and distinguish between inherited and newly observed mutations. In addition to resequencing data, we also provide metadata describing the lineage of each strain, and any targeted genetic modifications.

resequencing bio energy fermentation↗

Rapid Construction of an Infectious Clone of the Zika Virus, Strain ZKC2

Zika virus (ZIKV) has had detrimental effects on global public health in recent years. This is because the management of the disease has been limited, in part because its pathogenic mechanisms are not yet completely understood. Infectious clones are an important tool that utilize reverse genetics; these can be used to modify the ZIKV genomic RNA at the DNA level. A homologous recombination clone was used to construct pWSK29, a low copy plasmid that contained sequences for a T7 promoter, the whole genome of ZIKV ZKC2 strain, and a hepatitis delta virus ribozyme. High fidelity PCR was then used to amplify the T7 transcription template. The transcript was then transfected into susceptible cells via lipofection to recover the ZIKV ZKC2 strain. Finally, the virulence of rZKC2 was evaluated both in vitro and in vivo . The rZKC2 was successfully obtained and it showed the same virulence as its parent, the ZIKV ZKC2 strain (pZKC2), both in vitro and in vivo . The 3730 (NS2A-D62G) mutation site was identified as being important, since it had significant impacts on rZKC2 recovery. The 4015 (NS2A, A157V) mutation may reduce virus production by increasing the interferon type I response. In this study, one of the earliest strains of ZIKV that was imported into China was used for infectious clone construction and one possible site for antiviral medication development was discovered. The use of homologous recombination clones, of PCR products as templates for T7 transcription, and of lipofection for large RNA transfection could increase the efficiency of infectious clone construction. Our infectious clone provides an effective tool which can be used to explore the life cycle and medical treatment of ZIKV.

Qin, Zhiran↗

Molecular principles of Piwi-mediated cotranscriptional silencing through the dimeric SFiNX complex

Nuclear Argonaute proteins, guided by their bound small RNAs to nascent target transcripts, mediate cotranscriptional silencing of transposons and repetitive genomic loci through heterochromatin formation. The molecular mechanisms involved in this process are incompletely understood. Here, we show that the SFiNX complex, a silencing mediator downstream from nuclear Piwi-piRNA complexes in Drosophila, facilitates cotranscriptional silencing as a homodimer. The dynein light chain protein Cut up/LC8 mediates SFiNX dimerization, and its function can be bypassed by a heterologous dimerization domain, arguing for a constitutive SFiNX dimer. Dimeric, but not monomeric SFiNX, is capable of forming molecular condensates in a nucleic acid-stimulated manner. Mutations that prevent SFiNX dimerization result in loss of condensate formation in vitro and the inability of Piwi to initiate heterochromatin formation and silence transposons in vivo. We propose that multivalent SFiNX-nucleic acid interactions are critical for heterochromatin establishment at piRNA target loci in a cotranscriptional manner.

59 BASIC BIOLOGICAL SCIENCES↗

Widespread transposon co-option in the Caenorhabditis germline regulatory network

The movement of selfish DNA elements can lead to widespread genomic alterations with potential to create novel functions. We show that transposon expansions in Caenorhabditis nematodes led to extensive rewiring of germline transcriptional regulation. We find that about one-third of Caenorhabditis elegans germline-specific promoters have been co-opted from two related miniature inverted repeat transposable elements (TEs), CERP2 and CELE2. These promoters are regulated by HIM-17, a THAP domain–containing transcription factor related to a transposase. Expansion of CERP2 occurred before radiation of the Caenorhabditis genus, as did fixation of mutations in HIM-17 through positive selection, whereas CELE2 expanded only in C. elegans. Through comparative analyses in Caenorhabditis briggsae, we find not only evolutionary conservation of most CERP2 co-opted promoters but also a substantial fraction that are species-specific. Our work reveals the emergence and evolutionary conservation of a novel transcriptional network driven by TE co-option with a major impact on regulatory evolution.

59 BASIC BIOLOGICAL SCIENCES↗

Data for Comparison of Genotyping Assays for Detection of Targeted CRISPR/Cas Mutagenesis in Highly Polyploid Sugarcane

Sugarcane ( Saccharum spp.) is an important biofuel feedstock and a leading source of global table sugar. Saccharum hybrid cultivars are highly polyploid (2n = 100–130), containing large numbers of functionally redundant hom(e)ologs in their genomes. Genome editing with sequence-specific nucleases holds tremendous promise for sugarcane breeding. However, identification of plants with the desired level of co-editing within a pool of primary transformants can be difficult. While DNA sequencing provides direct evidence of targeted mutagenesis, it is cost-prohibitive as a primary screening method in sugarcane and most other methods of identifying mutant lines have not been optimized for use in highly polyploid species. In this study, non-sequencing methods of mutant screening, including capillary electrophoresis (CE), Cas9 RNP assay, and high-resolution melt analysis (HRMA), were compared to assess their potential for CRISPR/Cas9-mediated mutant screening in sugarcane. These assays were used to analyze sugarcane lines containing mutations at one or more of six sgRNA target sites. All three methods distinguished edited lines from wild type, with co-mutation frequencies ranging from 2% to 100%. Cas9 RNP assays were able to identify mutant sugarcane lines with as low as 3.2% co-mutation frequency, and samples could be scored based on undigested band intensity. CE was highlighted as the most comprehensive assay, delivering precise information on both mutagenesis frequency and indel size to a 1 bp resolution across all six targets. This represents an economical and comprehensive alternative to sequencing-based genotyping methods which could be applied in other polyploid species.

Genomics↗

Multiplex knockout of trichome-regulating MYB duplicates in hybrid poplar using a single gRNA

As the focus for CRISPR/Cas-edited plants moves from proof-of-concept to real-world applications, precise gene manipulation will increasingly require concurrent multiplex editing for polygenic traits. A common approach for editing across multiple sites is to design one guide RNA (gRNA) per target; however, this complicates construct assembly and increases the possibility of off-target mutations. In this study, we utilized one gRNA to target MYB186, a known positive trichome regulator, as well as its paralogs MYB138 and MYB38 at a consensus site for mutagenesis in hybrid poplar (Populus tremula × P. alba INRA 717-1B4). Unexpected duplications of MYB186 and MYB138 resulted in eight alleles for the three targeted genes in the hybrid poplar. Deep sequencing and polymerase chain reaction analyses confirmed editing across all eight targets in nearly all of the resultant glabrous mutants, ranging from small indels to large genomic dropouts, with no off-target activity detected at four potential sites. This highlights the effectiveness of a single gRNA targeting conserved exonic regions for multiplex editing. Additionally, cuticular wax and whole-leaf analyses showed a complete absence of triterpenes in the trichomeless mutants, hinting at a previously undescribed role for the nonglandular trichomes of poplar.

59 BASIC BIOLOGICAL SCIENCES↗

Lineage frequency time series reveal elevated levels of genetic drift in SARS-CoV-2 transmission in England

Genetic drift in infectious disease transmission results from randomness of transmission and host recovery or death. The strength of genetic drift for SARS-CoV-2 transmission is expected to be high due to high levels of superspreading, and this is expected to substantially impact disease epidemiology and evolution. However, we don’t yet have an understanding of how genetic drift changes over time or across locations. Furthermore, noise that results from data collection can potentially confound estimates of genetic drift. To address this challenge, we develop and validate a method to jointly infer genetic drift and measurement noise from time-series lineage frequency data. Our method is highly scalable to increasingly large genomic datasets, which overcomes a limitation in commonly used phylogenetic methods. We apply this method to over 490,000 SARS-CoV-2 genomic sequences from England collected between March 2020 and December 2021 by the COVID-19 Genomics UK (COG-UK) consortium and separately infer the strength of genetic drift for pre-B.1.177, B.1.177, Alpha, and Delta. We find that even after correcting for measurement noise, the strength of genetic drift is consistently, throughout time, higher than that expected from the observed number of COVID-19 positive individuals in England by 1 to 3 orders of magnitude, which cannot be explained by literature values of superspreading. Our estimates of genetic drift suggest low and time-varying establishment probabilities for new mutations, inform the parametrization of SARS-CoV-2 evolutionary models, and motivate future studies of the potential mechanisms for increased stochasticity in this system.

60 APPLIED LIFE SCIENCES↗

Human XPG nuclease structure, assembly, and activities with insights for neurodegeneration and cancer from pathogenic mutations

Xeroderma pigmentosum group G (XPG) protein is both a functional partner in multiple DNA damage responses (DDR) and a pathway coordinator and structure-specific endonuclease in nucleotide excision repair (NER). Different mutations in the XPG gene ERCC5 lead to either of two distinct human diseases: Cancer-prone xeroderma pigmentosum (XP-G) or the fatal neurodevelopmental disorder Cockayne syndrome (XP-G/CS). To address the enigmatic structural mechanism for these differing disease phenotypes and for XPG’s role in multiple DDRs, here we determined the crystal structure of human XPG catalytic domain (XPGcat), revealing XPG-specific features for its activities and regulation. Furthermore, XPG DNA binding elements conserved with FEN1 superfamily members enable insights on DNA interactions. Notably, all but one of the known pathogenic point mutations map to XPGcat, and both XP-G and XP-G/CS mutations destabilize XPG and reduce its cellular protein levels. Mapping the distinct mutation classes provides structure-based predictions for disease phenotypes: Residues mutated in XP-G are positioned to reduce local stability and NER activity, whereas residues mutated in XP-G/CS have implied long-range structural defects that would likely disrupt stability of the whole protein, and thus interfere with its functional interactions. Combined data from crystallography, biochemistry, small angle X-ray scattering, and electron microscopy unveil an XPG homodimer that binds, unstacks, and sculpts duplex DNA at internal unpaired regions (bubbles) into strongly bent structures, and suggest how XPG complexes may bind both NER bubble junctions and replication forks. Collective results support XPG scaffolding and DNA sculpting functions in multiple DDR processes to maintain genome stability.

59 BASIC BIOLOGICAL SCIENCES↗

Identification of genetic interactions with priB links the PriA/PriB DNA replication restart pathway to double-strand DNA break repair in Escherichia coli

Collisions between DNA replication complexes (replisomes) and impediments such as damaged DNA or proteins tightly bound to the chromosome lead to premature dissociation of replisomes at least once per cell cycle in Escherichia coli. Left unrepaired, these events produce incompletely replicated chromosomes that cannot be properly partitioned into daughter cells. DNA replication restart, the process that reloads replisomes at prematurely terminated sites, is therefore essential in E. coli and other bacteria. Three replication restart pathways have been identified in E. coli: PriA/PriB, PriA/PriC, and PriC/Rep. A limited number of genetic interactions between replication restart and other genome maintenance pathways have been defined, but a systematic study placing replication restart reactions in a broader cellular context has not been performed. We have utilized transposon-insertion sequencing to identify new genetic interactions between DNA replication restart pathways and other cellular systems. Known genetic interactors with the priB replication restart gene (uniquely involved in the PriA/PriB pathway) were confirmed and several novel priB interactions were discovered. Targeted genetic and imaging-based experiments with priB and its genetic partners revealed significant double-strand DNA break accumulation in strains with mutations in dam, rep, rdgC, lexA, or polA. Modulating the activity of the RecA recombinase partially suppressed the detrimental effects of rdgC or lexA mutations in ΔpriB cells. Taken together, our results highlight roles for several genes in double-strand DNA break homeostasis and define a genetic network that facilitates DNA repair/processing upstream of PriA/PriB-mediated DNA replication restart in E. coli.

59 BASIC BIOLOGICAL SCIENCES↗