Search NASA⌕ Search

SEARCH · Search NASA

Results for “sequence divergence”

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

DIVA/DeviceEditor (DIVA) v6.0.0

The DIVA software interfaces a process in which researchers design their DNA with a web-based graphical user interface (DeviceEditor), submit their designs to a central queue, and a few weeks later receive their sequence-verified clonal constructs. Each researcher independently designs the DNA to be constructed with a web-based BioCAD tool, and presses a button to submit their designs to a central queue. Researchers have web-based access to their DNA design queues, and can track the progress of their submitted designs as they progress from "evaluation", to "waiting for reagents", to "in progress", to "complete". Researchers access their completed constructs through the central DNA repository. Along the way, all DNA construction success/failure rates are captured in a central database. Once a design has been submitted to the queue, a small number of dedicated staff evaluate the design for feasibility and provide feedback to the responsible researcher if the design is either unreasonable (e.g., encompasses a combinatorial library of a billion constructs) or small design changes could significantly facilitate the downstream implementation process. The dedicated staff then use DNA assembly design automation software to optimize the DNA construction process for the design, leveraging existing parts from the DNA repository where possible and ordering synthetic DNA where necessary. Once all requisite process inputs are available, the design progresses from "waiting for reagents" to "in progress" in the design queue. Human-readable and machine-parseable DNA construction protocols output by the DNA assembly design automation software are then executed by the dedicated staff exploiting lab automation devices wherever possible. Since the all employed DNA construction methods are sequence-agnostic, standardized (utilize the same enzymatic master mixes and reaction conditions), completely independent DNA construction tasks can be aggregated into the same multi-well plates and pursued in parallel. The resulting sets of cloned constructs can then be screened by high-throughput next-gen sequencing platforms for sequence correctness. A combination of long read-length (e.g., PacBio) and paired-end read platforms (e.g., Illumina) would be exploited depending the particular task at hand (e.g., PacBio might be sufficient to screen a set of pooled constructs with significant gene divergence). Post sequence verification, designs for which at least one correct clone was identified will progress to a "complete" status, while designs for which no correct clones were identified will progress to a "failure" status. Depending on the failure mode (e.g., no transformants), and how many prior attempts/variations of assembly protocol have been already made for a given design, subsequent attempts may be made or the design can progress to a "permanent failure" state. All success and failure rate information will be captured during the process, including at which stage a given clonal construction procedure failed (e.g., no PCR product) and what the exact failure was (e.g. assembly piece 2 missing). This success/failure rate data can be leveraged to refine the DNA assembly design process.

Plahar, Hector↗
DOE OSTI2016

DIVA V2.0

The DIVA software interfaces a process in which researchers design their DNA with a web-based graphical user interface, submit their designs to a central queue, and a few weeks later receive their sequence-verified clonal constructs. Each researcher independently designs the DNA to be constructed with a web-based BioCAD tool, and presses a button to submit their designs to a central queue. Researchers have web-based access to their DNA design queues, and can track the progress of their submitted designs as they progress from "evaluation", to "waiting for reagents", to "in progress", to "complete". Researchers access their completed constructs through the central DNA repository. Along the way, all DNA construction success/failure rates are captured in a central database. Once a design has been submitted to the queue, a small number of dedicated staff evaluate the design for feasibility and provide feedback to the responsible researcher if the design is either unreasonable (e.g., encompasses a combinatorial library of a billion constructs) or small design changes could significantly facilitate the downstream implementation process. The dedicated staff then use DNA assembly design automation software to optimize the DNA construction process for the design, leveraging existing parts from the DNA repository where possible and ordering synthetic DNA where necessary. SynTrack software manages the physical locations and availability of the various requisite reagents and process inputs (e.g., DNA templates). Once all requisite process inputs are available, the design progresses from "waiting for reagents" to "in progress" in the design queue. Human-readable and machine-parseable DNA construction protocols output by the DNA assembly design automation software are then executed by the dedicated staff exploiting lab automation devices wherever possible. Since the all employed DNA construction methods are sequence-agnostic, standardized (utilize the same enzymatic master mixes and reaction conditions), completely independent DNA construction tasks can be aggregated into the same multi-well plates and pursued in parallel. The resulting sets of cloned constructs can then be screened by high-throughput next-gen sequencing platforms for sequence correctness. A combination of long read-length (e.g., PacBio) and paired-end read platforms (e.g., Illumina) would be exploited depending the particular task at hand (e.g., PacBio might be sufficient to screen a set of pooled constructs with significant gene divergence). Post sequence verification, designs for which at least one correct clone was identified will progress to a "complete" status, while designs for which no correct clones wereidentified will progress to a "failure" status. Depending on the failure mode (e.g., no transformants), and how many prior attempts/variations of assembly protocol have been already made for a given design, subsequent attempts may be made or the design can progress to a "permanent failure" state. All success and failure rate information will be captured during the process, including at which stage a given clonal construction procedure failed (e.g., no PCR product) and what the exact failure was (e.g. assembly piece 2 missing). This success/failure rate data can be leveraged to refine the DNA assembly design process.

CHEN, JOANNA↗
DOE OSTI2022

Soils and sediments host Thermoplasmata archaea encoding novel copper membrane monooxygenases (CuMMOs)

Copper membrane monooxygenases (CuMMOs) play critical roles in the global carbon and nitrogen cycles. Organisms harboring these enzymes perform the first, and rate limiting, step in aerobic oxidation of ammonia, methane, or other simple hydrocarbons. Within archaea, only organisms in the order Nitrososphaerales (Thaumarchaeota) encode CuMMOs, which function exclusively as ammonia monooxygenases. From grassland and hillslope soils and aquifer sediments, we identified 20 genomes from distinct archaeal species encoding divergent CuMMO sequences. These archaea are phylogenetically clustered in a previously unnamed Thermoplasmatota order, herein named the Ca. Angelarchaeales. The CuMMO proteins in Ca. Angelarchaeales are more similar in structure to those in Nitrososphaerales than those of bacteria, and contain all functional residues required for general monooxygenase activity. Ca. Angelarchaeales genomes are significantly enriched in blue copper proteins (BCPs) relative to sibling lineages, including plastocyanin-like electron carriers and divergent nitrite reductase-like (nirk) 2-domain cupredoxin proteins co-located with electron transport machinery. Ca. Angelarchaeales also encode significant capacity for peptide/amino acid uptake and degradation and share numerous electron transport mechanisms with the Nitrososphaerales. Ca. Angelarchaeales are detected at high relative abundance in some of the environments where their genomes originated from. While the exact substrate specificities of the novel CuMMOs identified here have yet to be determined, activity on ammonia is possible given their metabolic and ecological context. The identification of an archaeal CuMMO outside of the Nitrososphaerales significantly expands the known diversity of CuMMO enzymes in archaea and suggests previously unaccounted organisms contribute to critical global nitrogen and/or carbon cycling functions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

A novel bacterial protein family that catalyses nitrous oxide reduction

Nitrous oxide (N 2 O), a driver of global warming and climate change, has reached unprecedented concentrations in Earth’s atmosphere. Current N 2 O sources outpace N 2 O sinks, emphasizing the need for comprehensive understanding of processes that consume N 2 O. Microbes that express the enzyme N 2 O reductase (N 2 OR) convert N 2 O to climate change-neutral dinitrogen (N 2 ). Known N 2 ORs belong to the canonical clade I and clade II NosZ reductases and are considered key enzymes for N 2 O reduction. Here we report a previously unrecognized protein family with a role in N 2 O reduction, clade III lactonase-type N 2 OR (L-N 2 OR), which diverges in sequence from canonical NosZ but conserves three-dimensional protein structural features. Integrated physiological, metagenomic, proteomic and structural modelling studies demonstrate that L-N 2 ORs catalyse N 2 O reduction. L-N 2 OR genes occur in several phyla, predominantly in uncultured taxa with broad geographic distribution. Our findings expand the known diversity of N 2 ORs and implicate previously unrecognized taxa (for example, Nitrospinota) in N 2 O consumption. In conclusion, the expansion of N 2 OR diversity and the identification of a novel type of catalyst for N 2 O reduction advances the understanding of N 2 O sinks, has implications for greenhouse gas emission and climate change modelling, and expands opportunities for innovative biotechnologies aimed at curbing N 2 O emissions.

He, Guang 何广 [Univ. of Tennessee, Knoxville, TN (U↗
DOE OSTI2019

Structure of interferon-stimulated gene product 15 (ISG15) from the bat species Myotis davidii and the impact of interdomain ISG15 interactions on viral protein engagement

Bats have long been observed to be the hosts and the origin of numerous human diseases. Bats, like all mammals, rely on a number of innate immune mechanisms to combat invading pathogens, including the interferon type I, II and III responses. Ubiquitin-like interferon-stimulated gene product 15 (ISG15) is a key modulator of these interferon responses. Within these pathways, ISG15 can serve to stabilize host proteins modulating innate immune responses and act as a cytokine. Post-translational modifications of viral proteins introduced by ISG15 have also been observed to directly affect the function of numerous viral proteins. Unlike ubiquitin, which is virtually identical across all animals, comparison of ISG15s across species reveals that they are relatively divergent, with sequence identity dropping to as low as ~58% among mammals. In addition to serving as an obstacle to the zoonotic transmission of influenza, these ISG15 species–species differences have also long been shown to have an impact on the function of viral deISGylases. Recently, the structure of the first nonhuman ISG15, originating from mouse, suggested that the structures of human ISG15 may not be reflective of other species. Here, the structure of ISG15 from the bat species Myotis davidii solved to 1.37 Å resolution is reported. Comparison of this ISG15 structure with those from human and mouse not only underscores the structural impact of ISG15 species–species differences, but also highlights a conserved hydrophobic motif formed between the two domains of ISG15. Using the papain-like deISGylase from Severe acute respiratory syndrome coronavirus as a probe, the biochemical importance of this motif in ISG15–protein engagements was illuminated.

Langley, Caroline↗
DOE OSTI2021

Widespread misregulation of inter-species hybrid transcriptomes due to sex-specific and sex-chromosome regulatory evolution

When gene regulatory networks diverge between species, their dysfunctional expression in inter-species hybrid individuals can create genetic incompatibilities that generate the developmental defects responsible for intrinsic post-zygotic reproductive isolation. Both cis- and trans-acting regulatory divergence can be hastened by directional selection through adaptation, sexual selection, and inter-sexual conflict, in addition to cryptic evolution under stabilizing selection. Dysfunctional sex-biased gene expression, in particular, may provide an important source of sexually-dimorphic genetic incompatibilities. Here, we characterize and compare male and female/hermaphrodite transcriptome profiles for sibling nematode species Caenorhabditis briggsae and C. nigoni, along with allele-specific expression in their F 1 hybrids, to deconvolve features of expression divergence and regulatory dysfunction. Despite evidence of widespread stabilizing selection on gene expression, misexpression of sex-biased genes pervades F 1 hybrids of both sexes. This finding implicates greater fragility of male genetic networks to produce dysfunctional organismal phenotypes. Spermatogenesis genes are especially prone to high divergence in both expression and coding sequences, consistent with a “faster male” model for Haldane’s rule and elevated sterility of hybrid males. Moreover, underdominant expression pervades male-biased genes compared to female-biased and sex-neutral genes and an excess of cis-trans compensatory regulatory divergence for X-linked genes underscores a “large-X effect” for hybrid male expression dysfunction. Extensive regulatory divergence in sex determination pathway genes likely contributes to demasculinization of XX hybrids. The evolution of genetic incompatibilities due to regulatory versus coding sequence divergence, however, are expected to arise in an uncorrelated fashion. This study identifies important differences between the sexes in how regulatory networks diverge to contribute to sex-biases in how genetic incompatibilities manifest during the speciation process.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

Benchmarking viromics: an in silico evaluation of metagenome-enabled estimates of viral community composition and diversity

Background: Viral metagenomics (viromics) is increasingly used to obtain uncultivated viral genomes, evaluate community diversity, and assess ecological hypotheses. While viromic experimental methods are relatively mature and widely accepted by the research community, robust bioinformatics standards remain to be established. Here we used in silico mock viral communities to evaluate the viromic sequence-to-ecological-inference pipeline, including (i) read pre-processing and metagenome assembly, (ii) thresholds applied to estimate viral relative abundances based on read mapping to assembled contigs, and (iii) normalization methods applied to the matrix of viral relative abundances for alpha and beta diversity estimates. Results: Tools specifically designed for metagenomes, specifically metaSPAdes, MEGAHIT, and IDBA-UD, were the most effective at assembling viromes. Read pre-processing, such as partitioning, had virtually no impact on assembly output, but may be useful when hardware is limited. Viral populations with 2–5 × coverage typically assembled well, whereas lesser coverage led to fragmented assembly. Strain heterogeneity within populations hampered assembly, especially when strains were closely related (average nucleotide identity, or ANI ≥97%) and when the most abundant strain represented <50% of the population. Viral community composition assessments based on read recruitment were generally accurate when the following thresholds for detection were applied: (i) ≥10 kb contig lengths to define populations, (ii) coverage defined from reads mapping at ≥90% identity, and (iii) ≥75% of contig length with ≥1 × coverage. Finally, although data are limited to the most abundant viruses in a community, alpha and beta diversity patterns were robustly estimated (±10%) when comparing samples of similar sequencing depth, but more divergent (up to 80%) when sequencing depth was uneven across the dataset. In the latter cases, the use of normalization methods specifically developed for metagenomes provided the best estimates. Conclusions: These simulations provide benchmarks for selecting analysis cut-offs and establish that an optimized sample-to-ecological-inference viromics pipeline is robust for making ecological inferences from natural viral communities. Continued development to better accessing RNA, rare, and/or diverse viral populations and improved reference viral genome availability will alleviate many of viromics remaining limitations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

Divergence in Gene Regulation Contributes to Sympatric Speciation of Shewanella baltica Strains

Niche partitioning and sequence evolution drive genomic and phenotypic divergence, which ultimately leads to bacterial diversification. This study investigated the genomic composition of two Shewanella baltica clades previously identified through multilocus sequencing typing and recovered from the redox transition zone in the central Baltic Sea. Comparative genomic analysis revealed significantly higher interclade than intraclade genomic dissimilarity and that a subset of genes present in clade A were associated with potential adaptation to respiration of sulfur compounds present in the redox transition zone. The transcriptomic divergence between two representative strains of clades A and D, OS185 and OS195, was also characterized and revealed marked regulatory differences. We found that both the transcriptional divergence of shared genes and expression of strain-specific genes led to differences in regulatory patterns between strains that correlate with environmental redox niches. For example, under anoxic conditions of respiratory nitrate ammonification, OS185—the strain isolated from a nitrate-rich environment—upregulated nearly twice the number of shared genes upregulated by OS195—the strain isolated from an H 2 S-containing anoxic environment. On the other hand, OS195 showed stronger induction of strain-specific genes, especially those associated with sulfur compound respiration, under thiosulfate-reducing conditions. A positive association between the level of transcriptional divergence and the level of sequence divergence for shared genes was also noted. Our results provide further support for the hypothesis that genomic changes impacting transcriptional regulation play an important role in the diversification of ecologically distinct populations. IMPORTANCE: This study examined potential mechanisms through which co-occurring Shewanella baltica strains diversified to form ecologically distinct populations. At the time of isolation, the strains studied composed the major fraction of culturable nitrate-reducing communities in the Baltica Sea. Analysis of genomic content of 13 S. baltica strains from two clades representing different ecotypes demonstrated that one clade specifically possesses a number of genes that could favor successful adaptation to respire sulfur compounds in the portion of the water column from which these strains were isolated. In addition, transcriptional profiling of fully sequenced strains representative of these two clades, OS185 and OS195, under oxygen-, nitrate-, and thiosulfate-respiring conditions demonstrated that the strains exhibit relatively similar transcriptional responses during aerobic growth but more-distinct transcriptional responses under nitrate- and thiosulfate-respiring conditions. Results from this study provide insights into how genomic and gene regulatory diversification together impacted the redox specialization of the S. baltica strains.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

PacBio high‐throughput multi‐locus sequencing reveals high genetic diversity in mushroom‐forming fungi

Abstract Multi‐locus sequence data are widely used in fungal systematic and taxonomic studies to delimit species and infer evolutionary relationships. We developed and assessed the efficacy of a multi‐locus pooled sequencing method using PacBio long‐read high‐throughput sequencing. Samples included fresh and dried voucher specimens, cultures and archival DNA extracts of Agaricomycetes with an emphasis on the order Cantharellales. Of the 283 specimens sequenced, 93.6% successfully amplified at one or more loci with a mean of 3.3 loci amplified. Our method recovered multiple sequence variants representing alleles of rDNA loci and single copy protein‐coding genesrpb1,rpb2 andtef1. Within‐sample genetic variation differed by locus and taxonomic group, with the greatest genetic divergence observed among sequence variants ofrpb2 andtef1 from corticioid Cantharellales. Our method is a cost‐effective approach for generating accurate multi‐locus sequence data coupled with recovery of alleles from polymorphic samples and multi‐organism specimens. These results have important implications for understanding intra‐individual genomic variation among genetic loci commonly used in species delimitation of fungi.

Biochemistry & Molecular Biology↗
DOE OSTI2024

Functional diversification within the heme-binding split-barrel family

Due to neofunctionalization, a single fold can be identified in multiple proteins that have distinct molecular functions. Depending on the time that has passed since gene duplication and the number of mutations, the sequence similarity between functionally divergent proteins can be relatively high, eroding the value of sequence similarity as the sole tool for accurately annotating the function of uncharacterized homologs. Here, we combine bioinformatic approaches with targeted experimentation to reveal a large multifunctional family of putative enzymatic and nonenzymatic proteins involved in heme metabolism. This family (homolog of HugZ (HOZ)) is embedded in the “FMN-binding split barrel” superfamily and contains separate groups of proteins from prokaryotes, plants, and algae, which bind heme and either catalyze its degradation or function as nonenzymatic heme sensors. In prokaryotes these proteins are often involved in iron assimilation, whereas several plant and algal homologs are predicted to degrade heme in the plastid or regulate heme biosynthesis. In the plant Arabidopsis thaliana, which contains two HOZ subfamilies that can degrade heme in vitro (HOZ1 and HOZ2), disruption of AtHOZ1 (AT3G03890) or AtHOZ2A (AT1G51560) causes developmental delays, pointing to important biological roles in the plastid. In the tree Populus trichocarpa, a recent duplication event of a HOZ1 ancestor has resulted in localization of a paralog to the cytosol. Structural characterization of this cytosolic paralog and comparison to published homologous structures suggests conservation of heme-binding sites. This study unifies our understanding of the sequence-structure-function relationships within this multilineage family of heme-binding proteins and presents new molecular players in plant and bacterial heme metabolism.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2018

Cytosolic expression, solution structures, and molecular dynamics simulation of genetically encodable disulfide‐rich de novo designed peptides

Abstract Disulfide‐rich peptides represent an important protein family with broad pharmacological potential. Recent advances in computational methods have made it possible to design new peptides which adopt a stable conformation de novo. Here, we describe a system to produce disulfide‐rich de novo peptides using Escherichia coli as the expression host. The advantage of this system is that it enables production of uniformly 13 C‐ and 15 N‐labeled peptides for solution nuclear magnetic resonance (NMR) studies. This expression system was used to isotopically label two previously reported de novo designed peptides, and to determine their solution structures using NMR. The ensemble of NMR structures calculated for both peptides agreed well with the design models, further confirming the accuracy of the design protocol. Collection of NMR data on the peptides under reducing conditions revealed a dependency on disulfide bonds to maintain stability. Furthermore, we performed long‐time molecular dynamics (MD) simulations with tempering to assess the stability of two families of de novo designed peptides. Initial designs which exhibited a stable structure during simulations were more likely to adopt a stable structure in vitro , but attempts to utilize this method to redesign unstable peptides to fold into a stable state were unsuccessful. Further work is therefore needed to assess the utility of MD simulation techniques for de novo protein design.

Buchko, Garry W.↗
DOE OSTI2018

Correlation of structure, function and protein dynamics in GH7 cellobiohydrolases from Trichoderma atroviride, T. reesei and T. harzianum

The ascomycete fungus Trichoderma reesei is the predominant source of enzymes for industrial conversion of lignocellulose. Its glycoside hydrolase family 7 cellobiohydrolase (GH7 CBH) TreCel7A constitutes nearly half of the enzyme cocktail by weight and is the major workhorse in the cellulose hydrolysis process. The orthologs from Trichoderma atroviride (TatCel7A) and Trichoderma harzianum (ThaCel7A) show high sequence identity with TreCel7A, ~ 80%, and represent naturally evolved combinations of cellulose-binding tunnel-enclosing loop motifs, which have been suggested to influence intrinsic cellobiohydrolase properties, such as endo-initiation, processivity, and off-rate. The TatCel7A, ThaCel7A, and TreCel7A enzymes were characterized for comparison of function. The catalytic domain of TatCel7A was crystallized, and two structures were determined: without ligand and with thio-cellotriose in the active site. Initial hydrolysis of bacterial cellulose was faster with TatCel7A than either ThaCel7A or TreCel7A. In synergistic saccharification of pretreated corn stover, both TatCel7A and ThaCel7A were more efficient than TreCel7A, although TatCel7A was more sensitive to thermal inactivation. Structural analyses and molecular dynamics (MD) simulations were performed to elucidate important structure/function correlations. Moreover, reverse conservation analysis (RCA) of sequence diversity revealed divergent regions of interest located outside the cellulose-binding tunnel of Trichoderma spp. GH7 CBHs. We hypothesize that the combination of loop motifs is the main determinant for the observed differences in Cel7A activity on cellulosic substrates. Fine-tuning of the loop flexibility appears to be an important evolutionary target in Trichoderma spp., a conclusion supported by the RCA data. Our results indicate that, for industrial use, it would be beneficial to combine loop motifs from TatCel7A with the thermostability features of TreCel7A. Furthermore, one region implicated in thermal unfolding is suggested as a primary target for protein engineering.

09 BIOMASS FUELS↗
DOE OSTI2022

Borgs are giant genetic elements with potential to expand metabolic capacity

Anaerobic methane oxidation exerts a key control on greenhouse gas emissions1, yet factors that modulate the activity of microorganisms performing this function remain poorly understood. Here we discovered extraordinarily large, diverse DNA sequences that primarily encode hypothetical proteins through studying groundwater, sediments and wetland soil where methane production and oxidation occur. Four curated, complete genomes are linear, up to approximately 1 Mb in length and share genome organization, including replichore structure, long inverted terminal repeats and genome-wide unique perfect tandem direct repeats that are intergenic or generate amino acid repeats. We infer that these are highly divergent archaeal extrachromosomal elements with a distinct evolutionary origin. Gene sequence similarity, phylogeny and local divergence of sequence composition indicate that many of their genes were assimilated from methane-oxidizing Methanoperedens archaea. We refer to these elements as ‘Borgs’. We identified at least 19 different Borg types coexisting with Methanoperedens spp. in four distinct ecosystems. Borgs provide methane-oxidizing Methanoperedens archaea access to genes encoding proteins involved in redox reactions and energy conservation (for example, clusters of multihaem cytochromes and methyl coenzyme M reductase). These data suggest that Borgs might have previously unrecognized roles in the metabolism of this group of archaea, which are known to modulate greenhouse gas emissions, but further studies are now needed to establish their functional relevance.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

Methylamine-specific methyltransferase paralogs in Methanosarcina are functionally distinct despite frequent gene conversion

Abstract Sequenced archaeal genomes are mostly smaller and more streamlined than typical bacterial genomes; however, members of the Methanosarcina genus within the Euryarchaeaota are a significant exception, with M. acetivorans being the largest archaeal genome (5.8 Mbp) sequenced thus far. This finding is partially explained by extensive gene duplication within Methanosarcina spp. Significantly, the evolutionary pressures leading to gene duplication and subsequent genome expansion have not been well investigated, especially with respect to biological methane production (methanogenesis), which is the key biological trait of these environmentally important organisms. In this study, we address this question by specifically probing the functional evolution of two methylamine-specific methyltransferase paralogs in members of the Methanosarcina genus. Using the genetically tractable strain, M. acetivorans, we first show that the two paralogs have distinct cellular functions: one being required for methanogenesis from methylamine, the other for use of methylamine as a nitrogen source. Subsequently, through comparative sequence analyses, we show that functional divergence of paralogs is primarily mediated by divergent evolution of the 5′ regulatory region, despite frequent gene conversion within the coding sequence. This unique evolutionary paradigm for functional divergence of genes post-duplication underscores a divergent role for an enzyme singularly associated with methanogenic metabolism in other aspects of cell physiology.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2003

Conservation of Endo16 expression in sea urchins despite evolutionary divergence in both cis and trans-acting components of transcriptional regulation

Evolutionary changes in transcriptional regulation undoubtedly play an important role in creating morphological diversity. However, there is little information about the evolutionary dynamics of cis-regulatory sequences. This study examines the functional consequence of evolutionary changes in the Endo16 promoter of sea urchins. The Endo16 gene encodes a large extracellular protein that is expressed in the endoderm and may play a role in cell adhesion. Its promoter has been characterized in exceptional detail in the purple sea urchin, Strongylocentrotus purpuratus. We have characterized the structure and function of the Endo16 promoter from a second sea urchin species, Lytechinus variegatus. The Endo16 promoter sequences have evolved in a strongly mosaic manner since these species diverged approximately 35 million years ago: the most proximal region (module A) is conserved, but the remaining modules (B-G) are unalignable. Despite extensive divergence in promoter sequences, the pattern of Endo16 transcription is largely conserved during embryonic and larval development. Transient expression assays demonstrate that 2.2 kb of upstream sequence in either species is sufficient to drive GFP reporter expression that correctly mimics this pattern of Endo16 transcription. Reciprocal cross-species transient expression assays imply that changes have also evolved in the set of transcription factors that interact with the Endo16 promoter. Taken together, these results suggest that stabilizing selection on the transcriptional output may have operated to maintain a similar pattern of Endo16 expression in S. purpuratus and L. variegatus, despite dramatic divergence in promoter sequence and mechanisms of transcriptional regulation.

NASA Discipline Evolutionary Biology↗
DOE OSTI2018

A traditional evolutionary history of foot-and-mouth disease viruses in Southeast Asia challenged by analyses of non-structural protein coding sequences

Recombination of rapidly evolving RNA-viruses provides an important mechanism for diversifcation, spread, and emergence of new variants with enhanced ftness. Foot-and-mouth disease virus (FMDV) causes an important transboundary disease of livestock that is endemic to most countries in Asia and Africa. Maintenance and spread of FMDV are driven by periods of dominance of specifc viral lineages. Current understanding of the molecular epidemiology of FMDV lineages is generally based on the phylogenetic relationship of the capsid-encoding genes, with less attention to the process of recombination and evolution of non-structural proteins. In this study, the putative recombination breakpoints of FMDVs endemic to Southeast Asia were determined using full-open reading frame sequences. Subsequently, the lineages’ divergence times of recombination-free genome regions were estimated. These analyses revealed a close relationship between two of the earliest endemic viral lineages that appear unrelated when only considering the phylogeny of their capsid proteins. Contrastingly, one lineage, named O/CATHAY, known for having a particular host predilection (pigs) has evolved independently. Additionally, intra-lineage recombination occurred at diferent breakpoints compared to the inter-lineage process. These results provide new insights about FMDV recombination patterns and the evolutionary interdependence of FMDV serotypes and lineages.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Permanent draft genome of strain ESFC-1: ecological genomics of a newly discovered lineage of filamentous diazotrophic cyanobacteria

The nonheterocystous filamentous cyanobacterium, strain ESFC-1, is a recently described member of the order Oscillatoriales within the Cyanobacteria . ESFC-1 has been shown to be a major diazotroph in the intertidal microbial mat system at Elkhorn Slough, CA, USA. Based on phylogenetic analyses of the 16S RNA gene, ESFC-1 appears to belong to a unique, genus-level divergence; the draft genome sequence of this strain has now been determined. Here we report features of this genome as they relate to the ecological functions and capabilities of strain ESFC-1. The 5,632,035 bp genome sequence encodes 4914 protein-coding genes and 92 RNA genes. One striking feature of this cyanobacterium is the apparent lack of either uptake or bi-directional hydrogenases typically expected within a diazotroph. In addition, a large genomic island is found that contains numerous low GC-content genes and genes related to extracellular polysaccharide production and cell wall synthesis and maintenance.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2016

Permanent Draft Genome of Strain ESFC-1: Ecological Genomics of a Newly Discovered Lineage of Filamentous Diazotrophic Cyanobacteria

The nonheterocystous filamentous cyanobacterium, strain ESFC-1, is a recently described member of the order Oscillatoriales within the Cyanobacteria. ESFC-1 has been shown to be a major diazotroph in the intertidal microbial mat system at Elkhorn Slough, CA, USA. Based on phylogenetic analyses of the 16S RNA gene, ESFC-1 appears to belong to a unique, genus-level divergence; the draft genome sequence of this strain has now been determined. Here we report features of this genome as they relate to the ecological functions and capabilities of strain ESFC-1. The 5,632,035 bp genome sequence encodes 4914 protein-coding genes and 92 RNA genes. One striking feature of this cyanobacterium is the apparent lack of either uptake or bi-directional hydrogenases typically expected within a diazotroph. Additionally, a large genomic island is found that contains numerous low GC-content genes and genes related to extracellular polysaccharide production and cell wall synthesis and maintenance.

R Craig Everroad↗