Search NASA⌕ Search

SEARCH · Search NASA

Results for “genetic diversity”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Retracted: Monitoring Hatchery Broodstock Composition and Genetic Variation of Spring/Summer Chinook Salmon in the Columbia River Basin with Multigeneration Pedigrees

Abstract Hatchery production of Chinook Salmon Oncorhynchus tshawytscha in the Columbia River basin comprises most of the anadromous salmonid production in this region. Hatchery facilities and programs serve to mitigate for impacts to salmonids due to the construction and operation of hydropower dams and habitat impacts from development in addition to the conservation and restoration of natural populations. A genetic method referred to as parentage-based tagging (PBT) enables highly reliable detection of hatchery-origin fish and inference of multigeneration pedigrees. This study compiles 11 years of PBT data from nearly 125,000 interior stream-type Chinook Salmon from 24 spawning hatcheries located in tributaries of the mid- and upper Columbia River as well as the Salmon, Clearwater, and Grande Ronde River subbasins. Multigenerational pedigrees allowed investigation of the proportion of natural- and hatchery-origin broodstock (pNOB and pHOB, respectively) for each hatchery and enumeration of the scale of production between segregated and integrated programs. We then compared how pHOB and the scale of production influenced the number of stray fish observed, genetic diversity, relatedness, and age-class compositions within broodstocks. Over 91.0% of hatchery broodstock could be assigned back to their parents, and there was an overall rate of less than 1.0% of broodstock that were unintentionally incorporated into nonnatal hatchery programs. We evaluated 14 segregated programs with 0.0–10.0% pNOB, 7 intermediate programs (10.1–50.0% pNOB), and 3 integrated programs (>50.0% pNOB). There was no correlation between the scale of production or pNOB with the level of genetic diversity, but as production size increased, so did the effective number of breeders. This study demonstrates the utility of PBT as a monitoring tool for hatchery broodstocks, and results suggest that segregated and integrated programs have tradeoffs that generally align with their intended broodstock management purpose of providing fish for harvest and/or fish for supplementation or reintroduction.

Horn, Rebekah L.↗
DOE OSTI2020

Genomics of sorghum local adaptation to a parasitic plant

Host–parasite coevolution can maintain high levels of genetic diversity in traits involved in species interactions. In many systems, host traits exploited by parasites are constrained by use in other functions, leading to complex selective pressures across space and time. Here, we study genome-wide variation in the staple crop Sorghum bicolor (L.) Moench and its association with the parasitic weed Striga hermonthica (Delile) Benth., a major constraint to food security in Africa. We hypothesize that geographic selection mosaics across gradients of parasite occurrence maintain genetic diversity in sorghum landrace resistance. Suggesting a role in local adaptation to parasite pressure, multiple independent loss-of-function alleles at sorghum LOW GERMINATION STIMULANT 1 (LGS1) are broadly distributed among African landraces and geographically associated with S. hermonthica occurrence. However, low frequency of these alleles within S. hermonthica-prone regions and their absence elsewhere implicate potential trade-offs restricting their fixation. LGS1 is thought to cause resistance by changing stereochemistry of strigolactones, hormones that control plant architecture and below-ground signaling to mycorrhizae and are required to stimulate parasite germination. Consistent with trade-offs, we find signatures of balancing selection surrounding LGS1 and other candidates from analysis of genome-wide associations with parasite distribution. Experiments with CRISPR–Cas9-edited sorghum further indicate that the benefit of LGS1-mediated resistance strongly depends on parasite genotype and abiotic environment and comes at the cost of reduced photosystem gene expression. Our study demonstrates long-term maintenance of diversity in host resistance genes across smallholder agroecosystems, providing a valuable comparison to both industrial farming systems and natural communities.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

A haplotype‐resolved reference genome of Quercus alba sheds light on the evolutionary history of oaks

Summary White oak ( Quercus alba ) is an abundant forest tree species across eastern North America that is ecologically, culturally, and economically important. We report the first haplotype‐resolved chromosome‐scale genome assembly of Q. alba and conduct comparative analyses of genome structure and gene content against other published Fagaceae genomes. We investigate the genetic diversity of this widespread species and the phylogenetic relationships among oaks using whole genome data. Despite strongly conserved chromosome synteny and genome size across Quercus , certain gene families have undergone rapid changes in size, including defense genes. Unbiased annotation of resistance (R) genes across oaks revealed that the overall number of R genes is similar across species – as are the chromosomal locations of R gene clusters – but, gene number within clusters is more labile. We found that Q. alba has high genetic diversity, much of which predates its divergence from other oaks and likely impacts divergence time estimations. Our phylogenetic results highlight widespread phylogenetic discordance across the genus. The white oak genome represents a major new resource for studying genome diversity and evolution in Quercus . Additionally, we show that unbiased gene annotation is key to accurately assessing R gene evolution in Quercus .

Larson, Drew A. [Department of Biology Indiana Uni↗
DOE OSTI2024

Climate adaptation and sustainability in switchgrass: exploring plant-microbe-soil interactions across continental scale environmental gradients

Less carbon-intensive energy sources are needed to reduce greenhouse gas emissions and their predicted role in climate change. There is growing interest in the potential of biofuels for meeting this need. A critical question is whether large-scale biofuel production can be sustainable over the time scales needed to mitigate our carbon debt from fossil fuel consumption. The carbon balance and ultimately the sustainability of biofuel feedstock production is the result of complex climate-coupled interactions between carbon fixation, sequestration, and release through combustion. Similarly, the long-term productivity of biofuels depends on the environmental factors limiting plant growth. These factors are often related to soil resources which involve complex interactions at the plant-microbe-soil interface impacting their availability and cycling. Our collaborative project addressed sustainable switchgrass (Panicum virgatum) production by exploring Plant Systems, Plant-Microbiome Interactions, and Ecosystem Processes through the integrating lens of Multi-Scale Modeling. Our research was based on detailed characterization of genetically diverse switchgrass genotypes planted in common gardens across a continental latitudinal gradient. The underlying theme of our Plant Systems research was the use of locally adapted plant material to explore plant function, to understand the mechanistic basis of environmental interactions, and to discover the plant genes important for adaptation and sustainability in the face of climate change. Our Plant-Microbiome Interaction project characterized the microbial communities associated with switchgrass using genomic tools. Our Ecosystem Processes research focused on carbon cycle responses at the ecosystem level using stand level plantings. Finally, our Multi-Scale Modeling helped to define conditions of a sustainable biofuel system and identify key tradeoffs between genetic diversity, productivity, and ecosystem services. Genome-wide association analyses were used to identify alleles that contribute to successful establishment and biomass production across North America. Together, our work provided a baseline analyses of the potential of switchgrass as a biofuel feedstock. Our project resulted in a number of successful outcomes. First, we were successful in collecting switchgrass germplasm across the species range, propagating the material, and establishing common garden experiments across the species range. In collaboration with DOE JGI, we successfully assembled the first tetraploid switchgrass genome and published this resource with an analyses of the genetic basis local adaptation from our gardens (Lowry et al. 2019, Lovell et al. 2021). The gardens were used to characterize the genetic architecture for a number of important plant phenotypes. Our project also conducted extensive sampling and sequencing to characterize the bacterial and fungal associates of switchgrass roots and leaves. We showed that host genotype, location, and harvesting practices can play a role in microbiome assembly (Singer et al. 2019 & 2022, Van Wallendael et al. 2020 & 2022, Edwards et al. 2023). Our ecosystem processes work created baseline dataset of carbon and nutrient cycling in realistic stand plantings of switchgrass. Data from this experiment provided new insight into the role of plant traits, phenology, and local environments in ecosystem processes like soil respiration, net-ecosystem exchange, and dynamics of soil and plant nutrients (Ricketts et al. 2023). Finally, our crop modelling experiments help to characterize the sensitivity of common modeling frameworks to parameters, identify key limiters of productivity across large geographic scales, and leverage patterns of local adaptation in prediction. Ultimately, these studies help to identify critical plant-microbe-soil traits that may be manipulated, through breeding or agronomic management, to improve the sustainability of biofuel feedstocks.

09 BIOMASS FUELS↗
DOE OSTI2021

Mitochondrial Genomes of the United States Distribution of Gray Fox (Urocyon cinereoargenteus) Reveal a Major Phylogeographic Break at the Great Plains Suture Zone

We examined phylogeographic structure in gray fox (Urocyon cinereoargenteus) across the United States to identify the location of secondary contact zone(s) between eastern and western lineages and investigate the possibility of additional cryptic intraspecific divergences. We generated and analyzed complete mitochondrial genome sequence data from 75 samples and partial control region mitochondrial DNA sequences from 378 samples to investigate levels of genetic diversity and structure through population- and individual-based analyses including estimates of divergence (FST and SAMOVA), median joining networks, and phylogenies. We used complete mitochondrial genomes to infer phylogenetic relationships and date divergence times of major lineages of Urocyon in the United States. Despite broad-scale sampling, we did not recover additional major lineages of Urocyon within the United States, but identified a deep east-west split (~0.8 million years) with secondary contact at the Great Plains Suture Zone and confirmed the Channel Island fox (Urocyon littoralis) is nested within U. cinereoargenteus. Genetic diversity declined at northern latitudes in the eastern United States, a pattern concordant with post-glacial recolonization and range expansion. Beyond the east-west divergence, morphologically-based subspecies did not form monophyletic groups, though unique haplotypes were often geographically limited. Gray foxes in the United States displayed a deep, cryptic divergence suggesting taxonomic revision is needed. Secondary contact at a common phylogeographic break, the Great Plains Suture Zone, where environmental variables show a sharp cline, suggests ongoing evolutionary processes may reinforce this divergence. Follow-up study with nuclear markers should investigate whether hybridization is occurring along the suture zone and characterize contemporary population structure to help identify conservation units. Comparative work on other wide-ranging carnivores in the region should test whether similar evolutionary patterns and processes are occurring.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Structural and genetic convergence of HIV-1 neutralizing antibodies in vaccinated non-human primates

A primary goal of HIV-1 vaccine development is the consistent elicitation of protective, neutralizing antibodies. While highly similar neutralizing antibodies (nAbs) have been isolated from multiple HIV-infected individuals, it is unclear whether vaccination can consistently elicit highly similar nAbs in genetically diverse primates. Here, we show in three outbred rhesus macaques that immunization with Env elicits a genotypically and phenotypically conserved nAb response. From these vaccinated macaques, we isolated four antibody lineages that had commonalities in immunoglobulin variable, diversity, and joining gene segment usage. Atomic-level structures of the antigen binding fragments of the two most similar antibodies showed nearly identical paratopes. The Env binding modes of each of the four vaccine-induced nAbs were distinct from previously known monoclonal HIV-1 neutralizing antibodies, but were nearly identical to each other. The similarities of these antibodies show that the immune system in outbred primates can respond to HIV-1 Env vaccination with a similar structural and genotypic solution for recognizing a particular neutralizing epitope. These results support rational vaccine design for HIV-1 that aims to reproducibly elicit, in genetically diverse primates, nAbs with specific paratope structures capable of binding conserved epitopes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Data for Genetic Variation in Zea mays Influences Microbial Nitrification and DeNitrification in Conventional Agroecosystems

Nitrogenous fertilizers provide a short-lived benefit to crops in agroecosystems, but stimulate nitrification and denitrification, processes that result in nitrate pollution, N2O production, and reduced soil fertility. Recent advances in plant microbiome science suggest that genetic variation in plants can modulate the composition and activity of rhizosphere N-cycling microorganisms. Here we attempted to determine whether genetic variation exists in Zea mays for the ability to influence the rhizosphere nitrifier and denitrifier microbiome under “real-world” conventional agricultural conditions. To capture an extensive amount of genetic diversity within maize we grew and sampled the rhizosphere microbiome of a diversity panel of germplasm that included ex-PVP inbreds ( Z. mays ssp. mays ), ex-PVP hybrids ( Z. mays ssp. may s), and teosinte ( Z. mays ssp. mexicana and Z. mays ssp. parviglumis ). From these samples, we characterized the microbiome, a suite of microbial genes involved in nitrification and denitrification and carried out N-cycling potential assays. Here we are showing that populations/genotypes of a single species can vary in their ecological interaction with denitrifers and nitrifers. Some hybrid and teosinte genotypes supported microbial communities with lower potential nitrification and potential denitrification activity in the rhizosphere, while inbred genotypes stimulated/did not inhibit these N-cycling activities. These potential differences translated to functional differences in N2O fluxes, with teosinte plots producing less GHG than maize plots. Taken together, these results suggest that Zea genetic variation can lead to changes in N-cycling processes that result in N leaching and N2O production, and thereby are selectable targets for crop improvement. Understanding the underlying genetic variation contributing to belowground microbiome N-cycling into our conventional agricultural system could be useful for sustainability.

Nitrogen↗
DOE OSTI2022

Molecular Epidemiology of HIV-1 in Ghana: Subtype Distribution, Drug Resistance and Coreceptor Usage

The greatest HIV-1 genetic diversity is found in West/Central Africa due to the pandemic’s origins in this region, but this diversity remains understudied. We characterized HIV-1 subtype diversity (from both sub-genomic and full-genome viral sequences), drug resistance and coreceptor usage in 103 predominantly (90%) antiretroviral-naive individuals living with HIV-1 in Ghana. Full-genome HIV-1 subtyping confirmed the circulating recombinant form CRF02_AG as the dominant (53.9%) subtype in the region, with the complex recombinant 06_cpx (4%) present as well. Unique recombinants, most of which were mosaics containing CRF02_AG and/or 06_cpx, made up 37% of sequences, while “pure” subtypes were rare (<6%). Pretreatment resistance to at least one drug class was observed in 17% of the cohort, with NNRTI resistance being the most common (12%) and INSTI resistance being relatively rare (2%). CXCR4-using HIV-1 sequences were identified in 23% of participants. Overall, our findings advance our understanding of HIV-1 molecular epidemiology in Ghana. Extensive HIV-1 genetic diversity in the region appears to be fueling the ongoing creation of novel recombinants, the majority CRF02_AG-containing, in the region. The relatively high prevalence of pretreatment NNRTI resistance but low prevalence of INSTI resistance supports the use of INSTI-based first-line regimens in Ghana.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Virulence factors and antimicrobial resistance profiles of Campylobacter isolates recovered from consecutively reused broiler litter

ABSTRACT Campylobacterinfections are a leading cause of bacterial diarrhea in humans globally. Infections are due to consumption of contaminated food products and are highly associated with chicken meat, with chickens being an important reservoir forCampylobacter. Here, we characterized the genetic diversity ofCampylobacter jejuni(C. jejuni) andCampylobacter coli(C. coli) detected in broiler chicken litter over three consecutive flocks and determined their antimicrobial resistance (ARM) and virulence factor (VF) profiles.Campylobacterwas detected in 9.38% (27/288) of litter samples collected. Antimicrobial susceptibility testing and whole genome sequencing were performed onC. jejuni(n= 39) andC. coli(n= 5) isolates.Campylobactervirulence factors differed within and across broiler houses but were explained by the broiler flock cohort raised on litter,Campylobacterspecies andCampylobactermultilocus sequence type (MLST). Virulence factors involved in the ability to invade and colonize host tissues and evade host defenses were present inC. jejuniisolates (ST-464) from flock cohorts 1 and 2 but absent inC. jejuniisolates (ST-48) from flock cohort 3.C. jejuniisolates from house three harbored a significantly higher proportion of virulence genes with functions related to glycosylation and immune evasion thanC. jejuniisolates from houses 1 and 2 (P< 0.01). AllC. jejuniisolates were susceptible to all antibiotics tested whileC. coli(n= 4) were resistant to tetracycline and harbored the tetracycline resistant ribosomal protection protein (TetO). Our results suggest that house environment and broiler management practices imposed selective pressures on virulence factors and antimicrobial resistance genes ofCampylobacter. IMPORTANCE Campylobacteris a leading cause of foodborne illness in the United States due to consumption of contaminated or mishandled food products, often associated with chicken meat.Campylobacteris common in the microbiota of avian and mammalian gut; however, acquisition of antimicrobial resistance genes (ARGs) and virulence factors (VFs) may result in strains that pose significant threat to public health. Although there are studies investigating the genetic diversity ofCampylobacterstrains isolated from post-harvest chicken samples, there are limited data on the genome characteristics of isolates recovered from preharvest broiler production. Here, we show thatCampylobacter jejuniandCampylobacter colidiffer in their carriage of antimicrobial resistance and virulence factors may also differ in their ability to persist in litter during consecutive grow-out of broiler flocks. We found that presence/absence of virulence factors needed for evasion of host defense mechanisms and gut colonization played an integral role in differentiatingCampylobacterstrains.

Microbiology↗
DOE OSTI2022

COMPILE: a GWAS computational pipeline for gene discovery in complex genomes

Abstract Background Genome-Wide Association Studies (GWAS) are used to identify genes and alleles that contribute to quantitative traits in large and genetically diverse populations. However, traits with complex genetic architectures create an enormous computational load for discovery of candidate genes with acceptable statistical certainty. We developed a streamlined computational pipeline for GWAS (COMPILE) to accelerate identification and annotation of candidate maize genes associated with a quantitative trait, and then matches maize genes to their closest rice and Arabidopsis homologs by sequence similarity. Results COMPILE executed GWAS using a Mixed Linear Model that incorporated, without compression, recent advancements in population structure control, then linked significant Quantitative Trait Loci (QTL) to candidate genes and RNA regulatory elements contained in any genome. COMPILE was validated using published data to identify QTL associated with the traits of α-tocopherol biosynthesis and flowering time, and identified published candidate genes as well as additional genes and non-coding RNAs. We then applied COMPILE to 274 genotypes of the maize Goodman Association Panel to identify candidate loci contributing to resistance of maize stems to penetration by larvae of the European Corn Borer ( Ostrinia nubilalis ). Candidate genes included those that encode a gene of unknown function, WRKY and MYB-like transcriptional factors, receptor-kinase signaling, riboflavin synthesis, nucleotide-sugar interconversion, and prolyl hydroxylation. Expression of the gene of unknown function has been associated with pathogen stress in maize and in rice homologs closest in sequence identity. Conclusions The relative speed of data analysis using COMPILE allowed comparison of population size and compression. Limitations in population size and diversity are major constraints for a trait and are not overcome by increasing marker density. COMPILE is customizable and is readily adaptable for application to species with robust genomic and proteome databases.

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

Low Mutation Rate and Atypical Mutation Spectrum in Prasinoderma coloniale : Insights From an Early Diverging Green Lineage

Mutations are the ultimate source of genetic diversity on which natural selection and genetic drift act, playing a crucial role in evolution and long-term adaptation. At the molecular level, the spontaneous mutation rate (µ), defined as the number of mutations per base per generation, thus determines the adaptive potential of a species. Through a mutation accumulation experiment, we estimate the mutation rate and spectrum in Prasinoderma coloniale, a phytoplankton species from an early-branching lineage within the Archaeplastida, characterized by an unusually high genomic guanine-cytosine (GC) content (69.8%). We find that P. coloniale has a very low total mutation rate of µ = 2.00 × 10 -10 . The insertion–deletion mutation rate is almost 5 times lesser than the single nucleotide mutation rate with µ ID = 3.40 × 10 -11 and µ SNM = 1.62 × 10 -10 . Prasinoderma coloniale also exhibits an atypical mutational spectrum: While essentially all other eukaryotes show a bias toward GC to AT mutations, no evidence of this AT-bias is observed in P. coloniale. Since cytosine methylation is known to be mutagenic, we hypothesized that this may result from an absence of C-methylation. Surprisingly, we found high levels of C-methylation (14% in 5mC, 25% in 5mCG contexts). Methylated cytosines did not show increased mutation rates compared with unmethylated ones, not supporting the prevailing notion that C-methylation universally leads to higher mutation rates. Overall, P. coloniale combines a GC-rich genome with a low mutation rate and original mutation spectrum, suggesting the almost universal AT-bias may not have been present in the ancestor of the green lineage.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Characterization of Senecavirus A Isolates Collected From the Environment of U.S. Sow Slaughter Plants

Vesicular disease caused by Senecavirus A (SVA) is clinically indistinguishable from foot-and-mouth disease (FMD) and other vesicular diseases of swine. When a vesicle is observed in FMD-free countries, a costly and time-consuming foreign animal disease investigation (FADI) is performed to rule out FMD. Recently, there has been an increase in the number of FADIs and SVA positive samples at slaughter plants in the U.S. The objectives of this investigation were to: (1) describe the environmental burden of SVA in sow slaughter plants; (2) determine whether there was a correlation between PCR diagnostics, virus isolation (VI), and swine bioassay results; and (3) phylogenetically characterize the genetic diversity of contemporary SVA isolates. Environmental swabs were collected from three sow slaughter plants (Plants 1-3) and one market-weight slaughter plant (Plant 4) between June to December 2020. Of the 426 samples taken from Plants 1-3, 304 samples were PCR positive and 107 were VI positive. There was no detection of SVA by PCR or VI at Plant 4. SVA positive samples were most frequently found in the summer (78.3% June-September, vs. 59.4% October-December), with a peak at 85% in August. Eighteen PCR positive environmental samples with a range of C t values were selected for a swine bioassay: a single sample infected piglets ( n = 2). A random subset of the PCR positive samples was sequenced; and phylogenetic analysis demonstrated co-circulation and divergence of two genetically distinct groups of SVA. These data demonstrate that SVA was frequently found in the environment of sow slaughter plants, but environmental persistence and diagnostic detection was not indicative of whether a sampled was infectious to swine. Consequently, a more detailed understanding of the epidemiology of SVA and its environmental persistence in the marketing chain is necessary to reduce the number of FADIs and aide in the development of control measures to reduce the spread of SVA.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Identification of pleiotropic loci mediating structural and non-structural carbohydrate accumulation within the sorghum bioenergy association panel using high-throughput markers

Molecular characterization of diverse germplasm can contribute to breeding programs by increasing genetic gain for sorghum [ Sorghum bicolor (L.) Moench] improvement. Identifying novel marker-trait associations and candidate genes enriches the existing genomic resources and can improve bioenergy-related traits using genomic-assisted breeding. In the current scenario, identifying the genetic loci underlying biomass and carbon partitioning is vital for ongoing efforts to maximize each carbon sink’s yield for bioenergy production. Here, we have processed a high-density genomic marker (22 466 550) data based on whole-genome sequencing (WGS) using a set of 365 accessions from the bioenergy association panel (BAP), which includes ~19.7 million (19 744 726) single nucleotide polymorphism (SNPs) and 2.7 million (~2 721 824) insertion deletions (indels). A set of high-quality filtered SNP (~5.48 million) derived markers facilitated the assessment of population structure, genetic diversity, and genome-wide association studies (GWAS) for various traits related to biomass and its composition using the BAP. The phenotypic traits for GWAS included seed color (SC), plant height (PH), days to harvest (DTH), fresh weight (FW), dry weight (DW), brix content % (BRX), neutral detergent fiber (NDF), acid detergent fiber (ADF), non-fibrous carbohydrate (NFC), and lignin content. Several novel loci and candidate genes were identified for bioenergy-related traits, and some well-characterized genes for plant height ( Dw1 and Dw2 ) and the YELLOW SEED1 locus ( Y1 ) were validated. We further performed a multi-variate adaptive shrinkage analysis to identify pleiotropic QTL, which resulted in several shared marker-trait associations among bioenergy and compositional traits. Significant marker-trait associations with pleiotropic effects can be used to develop molecular markers for trait improvement using a marker-assisted breeding approach. Significant nucleotide diversity and heterozygosity were observed between photoperiod-sensitive and insensitive individuals of the panel. This diverse bioenergy panel with genomic resources will provide an excellent opportunity for further genetic studies, including selecting parental lines for superior hybrid development to improve biomass-related traits in sorghum.

bioenergy association panel↗
DOE OSTI2023

Characterization of novel recombinant mycobacteriophages derived from homologous recombination between two temperate phages

Comparative analyses of mycobacteriophage genomes reveals extensive genetic diversity in genome organization and gene content, contributing to widespread mosaicism. We previously reported that the prophage of mycobacteriophage Butters (cluster N) provides defense against infection by Island3 (subcluster I1). To explore the anti-Island3 defense mechanism, we attempted to isolate Island3 defense escape mutants on a Butters lysogen, but only uncovered phages with recombinant genomes comprised of regions of Butters and Island3 arranged from left arm to right arm as Butters-Island3-Butters (BIBs). Recombination occurs within two distinct homologous regions that encompass lysin A, lysin B, and holin genes in one segment, and RecE and RecT genes in the other. Structural genes of mosaic BIB genomes are contributed by Butters while the immunity cassette is derived from Island3. Consequently, BIBs are morphologically identical to Butters (as shown by transmission electron microscopy) but are homoimmune with Island3. Recombinant phages overcome antiphage defense and silencing of the lytic cycle. We leverage this observation to propose a stratagem to generate novel phages for potential therapeutic use.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Can fractal dimensions objectivize gastropod shell morphometrics? A case study from Lake Lugu (SW China)

1. Morphometrics are fundamental for the analysis of size and shape in fossils, particularly because soft parts or DNA are rarely preserved and hard parts such as shells are commonly the only source of information. Geometric morphometrics, that is, landmark analysis, is well established for the description of shape but it exhibits a couple of shortcomings resulting from subjective choices during landmarking (number and position of landmarks) and from difficulties in resolving shape at the level of micro-sculpture. 2. With the aid of high-resolution 3D scanning technology and analyses of fractal dimensions, we test whether such shortcomings of linear and landmark morphometrics can be overcome. As a model group, we selected a clade of modern viviparid gastropods from Lake Lugu, with shells that show a high degree of sculptural variation. Linear and landmark analyses were applied to the same shells in order to establish the fractal dimensions. 3. The genetic diversity of the gastropod clade was assessed. The genetic results suggest that the gastropod clade represents a single species. The results of all morphometric methods applied are in line with the genetic results, which is that no specific morphotype could be delimited. Apart from this overall agreement, landmark and fractal dimension analyses do not correspond to each other but represent data sets with different information. Generally, the fractal dimension values quantify the roughness of the shell surface, the resolution of the 3D scans determining the level. In our approach, we captured the micro-sculpture but not the first-order sculptural elements, which explains that fractal dimension and landmark data are not in phase. 4. We can show that analyzing fractal dimensions of gastropod shells opens a window to more detailed information that can be considered in evolutionary and ecological contexts. We propose that using low-resolution 3D scans may successfully substitute landmark analyses because it overcomes the subjective landmarking. Analyses of 3D scans with higher resolution than used in this study will provide surface roughness information at the mineralogical level. We suggest that fractal dimension analyses of a combination of differently resolved 3D models will significantly improve the quality of shell morphometrics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A variant selection framework for genome graphs

Abstract Motivation Variation graph representations are projected to either replace or supplement conventional single genome references due to their ability to capture population genetic diversity and reduce reference bias. Vast catalogues of genetic variants for many species now exist, and it is natural to ask which among these are crucial to circumvent reference bias during read mapping. Results In this work, we propose a novel mathematical framework for variant selection, by casting it in terms of minimizing variation graph size subject to preserving paths of length α with at most δ differences. This framework leads to a rich set of problems based on the types of variants [e.g. single nucleotide polymorphisms (SNPs), indels or structural variants (SVs)], and whether the goal is to minimize the number of positions at which variants are listed or to minimize the total number of variants listed. We classify the computational complexity of these problems and provide efficient algorithms along with their software implementation when feasible. We empirically evaluate the magnitude of graph reduction achieved in human chromosome variation graphs using multiple α and δ parameter values corresponding to short and long-read resequencing characteristics. When our algorithm is run with parameter settings amenable to long-read mapping (α = 10 kbp, δ = 1000), 99.99% SNPs and 73% SVs can be safely excluded from human chromosome 1 variation graph. The graph size reduction can benefit downstream pan-genome analysis. Availability and implementation https://github.com/AT-CG/VF. Supplementary information Supplementary data are available at Bioinformatics online.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Genomic dissection of anthracnose resistance response in sorghum [Sorghum bicolor (L.) Moench]

Sorghum [Sorghum bicolor (L.) Moench] is the fifth most important grain crop behind maize, wheat, rice, and barley. Today, it is of interest as a source of fermentable sugars for the production of renewable fuels and chemicals, and as a source of biomass for co-firing. The productivity and profitability of sorghum are limited by several biotic constraints, most notably anthracnose caused by the fungal pathogen Colletotrichum sublineolum. The most cost-effective and environmentally benign strategy to control anthracnose is through the incorporation of resistance genes. Over the last three years, our research efforts have been directed to identify new sources of resistance in temperate adapted and tropical germplasm, and to delimited genomic regions associated with the observe anthracnose resistant response. Three biparental mapping populations derived from the resistant lines SC112-14, QL3 and IS18760 were evaluated for anthracnose resistance response in Texas, Georgia, Florida and Puerto Rico. In parallel, three high density recombination maps were constructed and used to identify resistant loci. Anthracnose resistant response in line SC112-14 is controlled by a major locus on chromosome 5. Segregation analysis of 1,500 progenies delimited the resistance locus on chromosome 5 to a 23-kb region harboring three candidate genes, including Sobic.005G17230 identified by GWAS of the sorghum association panel (SAP). The latter gene belongs to a family of genes encoding F-box proteins indicating that this resistance response involved in signaling cascades and transcriptional reprograming, rather than recognition of pathotype-associated molecular patterns. In contrast, anthracnose resistant response in lines QL3 and IS18760 is controlled by multiple small-effect genes. Greenhouse evaluation of a representative subset of the three mapping populations against nine pathotypes found that lines susceptible in the field could be resistant to a single pathotype in the greenhouse. Thus, the activation of a resistance response system by a single pathotype could not provide a broader resistance response against multiple pathotypes. The screening of 1,801 sweet sorghum accessions from the National Plant Germplasm System identified 654 accessions with Brix value larger than 10, which in turn was used to select a subset of 233 accessions for evaluation of anthracnose resistant response. Even though most of the accessions were not completely infected by anthracnose, 28 accessions were completely resistant against pathotypes from Texas, Georgia, Florida and Puerto Rico. Genotyping-by-sequencing analysis of this subset identified 157,843 single nucleotide polymorphisms. Population structure analysis of the subset based on a subset of 2,345 unlinked SNPs found that the genetic diversity could be divided into four populations. The genetic relatedness among accessions within populations suggests most of the resistant germplasm may contain few different resistance sources. These resistance sources present in sweet sorghum germplasm could expedite the development of new resistant sweet sorghum cultivars and hybrids by avoiding time-consuming introgression breeding approaches with non-sweet sorghums serving as donor of the resistance alleles.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

CMPLE: Correlation Modeling to Decode Photosynthesis Using the Minorize–Maximize Algorithm

In plant genomic experiments, correlations among various biological traits (phenotypes) give new insights into how genetic diversity may have tuned biological processes to enhance fitness under diverse conditions. Consequently, knowing how the correlations are affected by genetic (G) and environmental (E) factors helps develop climate-resilient plants. However, the current literature lacks any method for assessing the effect of predictors on pairwise correlations among multiple phenotypes together with easily interpretable model parameters. To address this need, we propose to model pairwise correlations directly in terms of G and E and develop a computationally efficient inference procedure. Two major novelties in our methodology are (1) the use of a composite pairwise likelihood method to avoid the positive definiteness restriction on the correlation matrix and (2) the use of a novel Minorize–Maximize (MM) algorithm for the efficient estimation of a large number of parameters. The proposed method shows excellent numerical performance on synthetic datasets. Here, the analysis of the motivating data on cowpea reveals that the rates of solar energy storage by photosynthesis (the aggregate trait) are differentially affected by different genetic loci through two distinct processes: “photoinhibition” which results from photodamage caused by excess light, and “photoprotection” which protects plants from photodamage but also results in energy loss.

Correlation modeling↗