Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “genomic methods”

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 271 records · Page 15

Intrahost SARS-CoV-2 k-mer Identification Method (iSKIM) for Rapid Detection of Mutations of Concern Reveals Emergence of Global Mutation Patterns

Despite unprecedented global sequencing and surveillance of SARS-CoV-2, timely identification of the emergence and spread of novel variants of concern (VoCs) remains a challenge. Several million raw genome sequencing runs are now publicly available. We sought to survey these datasets for intrahost variation to study emerging mutations of concern. We developed iSKIM (“intrahost SARS-CoV-2 k-mer identification method”) to relatively quickly and efficiently screen the many SARS-CoV-2 datasets to identify intrahost mutations belonging to lineages of concern. Certain mutations surged in frequency as intrahost minor variants just prior to, or while lineages of concern arose. The Spike N501Y change common to several VoCs was found as a minor variant in 834 samples as early as October 2020. This coincides with the timing of the first detected samples with this mutation in the Alpha/B.1.1.7 and Beta/B.1.351 lineages. Using iSKIM, we also found that Spike L452R was detected as an intrahost minor variant as early as September 2020, prior to the observed rise of the Epsilon/B.1.429/B.1.427 lineages in late 2020. iSKIM rapidly screens for mutations of interest in raw data, prior to genome assembly, and can be used to detect increases in intrahost variants, potentially providing an early indication of novel variant spread.

60 APPLIED LIFE SCIENCES↗

LOGAN: High-Performance GPU-Based X-Drop Long-Read Alignment

Pairwise sequence alignment is one of the most computationally intensive kernels in genomic data analysis, accounting for more than 90% of the runtime for key bioinformatics applications. This method is particularly expensive for third-generation sequences due to the high computational cost of analyzing sequences of length between 1Kb and 1Mb. Given the quadratic overhead of exact pairwise algorithms for long alignments, the community primarily relies on approximate algorithms that search only for high-quality alignments and stop early when one is not found. In this work, we present the first GPU optimization of the popular X-drop alignment algorithm, that we named LOGAN. Results show that our high-performance multi-GPU implementation achieves up to 181.6 GCUPS and speed-ups up to 6.6× and 30.7× using 1 and 6 NVIDIA Tesla V100, respectively, over the state-of-the-art software running on two IBM Power9 processors using 168 CPU threads, with equivalent accuracy. We also demonstrate a 2.3× LOGAN speed-up versus ksw2, a state-of-art vectorized algorithm for sequence alignment implemented in minimap2, a long-read mapping software. Furthermore, to highlight the impact of our work on a real-world application, we couple LOGAN with a many-to-many long-read alignment software called BELLA, and demonstrate that our implementation improves the overall BELLA runtime by up to 10.6×. Finally, we adapt the Roofline model for LOGAN and demonstrate that our implementation is near optimal on the NVIDIA Tesla V100s.

97 MATHEMATICS AND COMPUTING↗

synDNA—a Synthetic DNA Spike-in Method for Absolute Quantification of Shotgun Metagenomic Sequencing

The synDNAs designed in this study enable accurate and reproducible measurements of absolute amount and fold changes of bacterial species in complex microbial communities. The method proposed here is versatile and promising as it can be applied to bacterial communities or genomic features like genes and operons, in addition to being easily adaptable by other research groups at a low cost.

59 BASIC BIOLOGICAL SCIENCES↗

The Exometabolome of Xylella fastidiosa in Contact with Paraburkholderia phytofirmans Supernatant Reveals Changes in Nicotinamide, Amino Acids, Biotin, and Plant Hormones

Microbial competition within plant tissues affects invading pathogens’ fitness. Metabolomics is a great tool for studying their biochemical interactions by identifying accumulated metabolites. Xylella fastidiosa, a Gram-negative bacterium causing Pierce’s disease (PD) in grapevines, secretes various virulence factors including cell wall-degrading enzymes, adhesion proteins, and quorum-sensing molecules. These factors, along with outer membrane vesicles, contribute to its pathogenicity. Previous studies demonstrated that co-inoculating X. fastidiosa with the Paraburkholderia phytofirmans strain PsJN suppressed PD symptoms. Here, we further investigated the interaction between the phytopathogen and the endophyte by analyzing the exometabolome of wild-type X. fastidiosa and a diffusible signaling factor (DSF) mutant lacking quorum sensing, cultivated with 20% P. phytofirmans spent media. Liquid chromatography–mass spectrometry (LC-MS) and the Method for Metabolite Annotation and Gene Integration (MAGI) were used to detect and map metabolites to genomes, revealing a total of 121 metabolites, of which 25 were further investigated. These metabolites potentially relate to host adaptation, virulence, and pathogenicity. Notably, this study presents the first comprehensive profile of X. fastidiosa in the presence of a P. phytofirmans spent media. The results highlight that P. phytofirmans and the absence of functional quorum sensing affect the ratios of glutamine to glutamate (Gln:Glu) in X. fastidiosa. Additionally, two compounds with plant metabolism and growth properties, 2-aminoisobutyric acid and gibberellic acid, were downregulated when X. fastidiosa interacted with P. phytofirmans. These findings suggest that P. phytofirmans-mediated disease suppression involves modulation of the exometabolome of X. fastidiosa, impacting plant immunity.

60 APPLIED LIFE SCIENCES↗

Assessment of the Reversion to Virulence and Protective Efficacy in Pigs Receiving the Live Attenuated Classical Swine Fever Recombinant Vaccine Candidate FlagT4G

Background/Objectives: Control of classical swine fever virus (CSFV) in endemic countries relies on vaccination using live attenuated vaccines (LAVs). Most of these LAVs do not allow for the differentiation of vaccinated animals from infected animals (DIVA) based on their serological response. FlagT4G vaccine is a novel candidate that confers robust protective immunity early after vaccination and shows DIVA capabilities. Methods: This report presents the characterization of FlagT4G virus in terms of the stability of its genomic and attenuated phenotypes assessed by a reversion to virulence protocol, as well as its protective efficacy by determining the minimal protective dose. Results: Results presented here demonstrate that after five consecutive passages in groups of 5-week-old susceptible domestic pigs, FlagT4G virus remains genetically stable, and its attenuated phenotype remains unaltered. In terms of efficacy, FlagT4G virus induced solid protection against the intranasal challenge with 105 tissue culture infectious dose (TCID50) of virulent field isolate Brescia virus, even with a vaccine dose as low as 102 TCID50. Conclusions: Results presented here indicate that the FlagT4G vaccine may be a useful tool for CSFV control.

Immunology↗

Critical Assessment of Metagenome Interpretation: the second round of challenges

Abstract Evaluating metagenomic software is key for optimizing metagenome interpretation and focus of the Initiative for the Critical Assessment of Metagenome Interpretation (CAMI). The CAMI II challenge engaged the community to assess methods on realistic and complex datasets with long- and short-read sequences, created computationally from around 1,700 new and known genomes, as well as 600 new plasmids and viruses. Here we analyze 5,002 results by 76 program versions. Substantial improvements were seen in assembly, some due to long-read data. Related strains still were challenging for assembly and genome recovery through binning, as was assembly quality for the latter. Profilers markedly matured, with taxon profilers and binners excelling at higher bacterial ranks, but underperforming for viruses and Archaea. Clinical pathogen detection results revealed a need to improve reproducibility. Runtime and memory usage analyses identified efficient programs, including top performers with other metrics. The results identify challenges and guide researchers in selecting methods for analyses.

54 ENVIRONMENTAL SCIENCES↗

Genome-Scale Transcription-Translation Mapping Reveals Features of Zymomonas mobilis Transcription Units and Promoters

Efforts to rationally engineer synthetic pathways in Zymomonas mobilis are impeded by a lack of knowledge and tools for predictable and quantitative programming of gene regulation at the transcriptional, posttranscriptional, and posttranslational levels. With the detailed functional characterization of the Z. mobilis genome presented in this work, we provide crucial knowledge for the development of synthetic genetic parts tailored to Z. mobilis . This information is vital as researchers continue to develop Z. mobilis for synthetic biology applications. Our methods and statistical analyses also provide ways to rapidly advance the understanding of poorly characterized bacteria via empirical data that enable the experimental validation of sequence-based prediction for genome characterization and annotation.

59 BASIC BIOLOGICAL SCIENCES↗

Disentangling Sources of Gene Tree Discordance in Phylogenomic Data Sets: Testing Ancient Hybridizations in Amaranthaceae s.l

Gene tree discordance in large genomic data sets can be caused by evolutionary processes such as incomplete lineage sorting and hybridization, as well as model violation, and errors in data processing, orthology inference, and gene tree estimation. Species tree methods that identify and accommodate all sources of conflict are not available, but a combination of multiple approaches can help tease apart alternative sources of conflict. Here, using a phylotranscriptomic analysis in combination with reference genomes, we test a hypothesis of ancient hybridization events within the plant family Amaranthaceae s.l. that was previously supported by morphological, ecological, and Sanger-based molecular data. The data set included seven genomes and 88 transcriptomes, 17 generated for this study. We examined gene-tree discordance using coalescent-based species trees and network inference, gene tree discordance analyses, site pattern tests of introgression, topology tests, synteny analyses, and simulations. We found that a combination of processes might have generated the high levels of gene tree discordance in the backbone of Amaranthaceae s.l. Furthermore, we found evidence that three consecutive short internal branches produce anomalous trees contributing to the discordance. Overall, our results suggest that Amaranthaceae s.l. might be a product of an ancient and rapid lineage diversification, and remains, and probably will remain, unresolved. This work highlights the potential problems of identifiability associated with the sources of gene tree discordance including, in particular, phylogenetic network methods. Our results also demonstrate the importance of thoroughly testing for multiple sources of conflict in phylogenomic analyses, especially in the context of ancient, rapid radiations. We provide several recommendations for exploring conflicting signals in such situations.

59 BASIC BIOLOGICAL SCIENCES↗

Novel mutation leading to splice donor loss in a conserved site of DMD gene causes Duchenne muscular dystrophy with cryptorchidism

Background As one of the most common congenital abnormalities in male births, cryptorchidism has been found to have a polygenic aetiology according to previous studies of common variants. However, little is known about genetic predisposition of rare variants for cryptorchidism, since rare variants have larger effective size on diseases than common variants. Methods In this study, a cohort of 115 Chinese probands with cryptorchidism was analysed using whole-genome sequencing, alongside 19 parental controls and 2136 unaffected men. Additionally, CRISPR-Cas9 editing of a conserved variant was performed in a mouse model, with MRI screening used to observe the phenotype. Results In 30 of 115 patients (26.1%), we identified four novel genes ( ARSH , DMD , MAGEA4 and SHROOM2 ) affecting at least five unrelated patients and four known genes ( USP9Y , UBA1 , BCORL1 and KDM6A ) with the candidate rare pathogenic variants affecting at least two cases. Burden tests of rare variants revealed the genome-wide significances for newly identified genes (p<2.5×10 −6 ) under the Bonferroni correction. Surprisingly, novel and known genes were mainly found on X chromosome (seven on X and one on Y) and all rare X-chromosomal segregating variants exhibited a maternal inheritance rather than de novo origin. CRISPR-Cas9 mouse modelling of a splice donor loss variant in DMD (NC_000023.11:g.32454661C>G), which resides in a conserved site across vertebrates, replicated bilateral cryptorchidism phenotypes, confirmed by MRI at 4 and 10 weeks. The movement tests further revealed symptoms of Duchenne muscular dystrophy (DMD) in transgenic mice. Conclusion Our results revealed the role of the DMD gene mutation in causing cryptorchidism. The results also suggest that maternal-X inheritance of pathogenic defects could have a predominant role in the development of cryptorchidism.

Chen, Jianhai (ORCID:0000000300932003)↗

Host phylogeny matters: Examining sources of variation in infection risk by blood parasites across a tropical montane bird community in India

Abstract Background Identifying patterns and drivers of infection risk among host communities is crucial to elucidate disease dynamics and predict infectious disease risk in wildlife populations. Blood parasites of the genera Plasmodium and Haemoproteus are a diverse group of vector-borne protozoan parasites that affect bird populations globally. Despite their widespread distribution and exceptional diversity, factors underlying haemosporidian infection risk in wild bird communities remain poorly understood. While some studies have examined variation in avian haemosporidian risk, researchers have primarily focused on host ecological traits without considering host phylogenetic relationships. In this study, we employ a phylogenetically informed approach to examine the association between host ecological traits and haemosporidian infection risk in endemic bird communities in the Western Ghats Sky Islands. Methods We used parasite sequence data based on partial mitochondrial cytochrome b gene, that was amplified from genomic DNA extracted from 1177 birds (28 species) across the Western Ghats to assess infection of birds with haemosporidian parasites. We employed a Bayesian phylogenetic mixed effect modelling approach to test whether haemosporidian infection risk was affected by seven species-specific and four individual-level ecological predictors. We also examined the effect of host phylogenetic relationships on the observed patterns of variation in haemosporidian infection risk by estimating phylogenetic signal. Results Our study shows that host ecological traits and host phylogeny differentially influence infection risk by Plasmodium (generalist parasite) and Haemoproteus (specialist parasite). For Plasmodium , we found that sociality, sexual dimorphism and foraging strata were important ecological predictors. For Haemoproteus , patterns of infection risk among host species were associated with sociality, species elevation and individual body condition. Interestingly, variance in infection risk explained by host phylogeny was higher for Haemoproteus parasites compared to Plasmodium . Conclusions Our study highlights that while host ecological traits promoting parasite exposure and host susceptibility are important determinants of infection risk, host phylogeny also contributes substantially to predicting patterns of haemosporidian infection risk in multi-host communities. Importantly, infection risk is driven by joint contributions of host ecology and host phylogeny and studying these effects together could increase our ability to better understand the drivers of infection risk and predict future disease threats. Graphical abstract

59 BASIC BIOLOGICAL SCIENCES↗

Case Study: Can you find Delftia?

Delftia is a genus of bacteria with a bunch of cool features! The best studied species, Delftia acidovorans can produce gold nanoparticles from gold ions in solution. This bacteria has been found living in biofolms with Cupriavidis metallidurans on gold nuggets. It is also found in soil, in sinks and in rhizospheres of different plants where it promotes their growth. Delftia acidovorans forms gold nuggets by producing a short nonribosomal peptide called delftibactin. The 16 genes responsible for the production of delftibactin are called the del cluster (delA-delP). These genes were originally discovered in Delftia acidovorans SPH-1, but our research shows that the del cluster appears to be present across the genus. We'll start with a number of raw reads, clean them up and do a little taxonomy to figure out what we might be able to assemble. Next, we'll assemble the reads using 4 different methods and pick out the best assembly. Then we'll take those reads and sort them by genome and assess the quality of those genomes. We'll select the high quality genomes to annotate and insert into phylogenetic trees to find relatives.

59 BASIC BIOLOGICAL SCIENCES↗

Direct genome-scale screening of Gluconobacter oxydans B58 for rare earth element bioleaching

Abstract The transition to a sustainable energy economy will require an enormous increase in the supply of rare earth elements (REEs). Bioleaching offers a promising alternative to conventional hydrometallurgical methods for REE extraction from low-grade ores. However, exploiting this potential remains challenging due to large gaps in our understanding of the genetics involved, and inadequate biological tools to address them. We generated a highly non-redundant whole-genome knockout collection for the bioleaching microbe Gluconobacter oxydans B58, reducing redundancy by 85% compared to the previous best collection. This new collection was directly screened for bioleaching neodymium from a synthetic monazite powder, identifying 89 genes important for bioleaching, 68 of which have not previously been associated with this mechanism. We conducted bench-scale experiments to validate the extraction efficiency of promising strains: 8 demonstrated significant increases in extraction by up to 111% (δ GO_1598 , disruption of the gene encoding the orotate phosphoribosyltransferase enzyme PyrE), and one strain significantly reduced it by 97% (δ GO_1096 , disruption of the gene encoding the GTP-binding protein TypA). Notable changes in pH were only observed for 3 strains, suggesting an important role for non-acid mechanisms in bioleaching. These findings provide valuable insights into further enhancing REE-bioleaching by G. oxydans through genetic engineering.

Marecos, Sabrina (ORCID:0000000158529933)↗

Optimizing resource allocation in Miscanthus breeding via sparse testing designs for genomic prediction

Phenotyping high-biomass perennial crops is laborious and the rate of genetic gain in conventional perennial crop breeding programs is typically low. So, it is especially important to identify methods that produce efficiency gains in the breeding process. Miscanthus is a C4 perennial grass with favorable characteristics for producing biomass as a feedstock for biofuels and diverse bio-based products. Increasing biomass yield will increase profitability and environmental benefits, so it is a key target for Miscanthus breeding. In addition, the identification of well-adapted genotypes across a wide range of environmental conditions requires the establishment of multi-environment trials (METs). Sparse testing is a genomic prediction-based strategy that reduces the phenotyping costs in METs by selecting a subset of genotypes to evaluate in a subset of environments and then predicts the performance of the unobserved genotype-environment combinations. A Miscanthus sacchariflorus (MSA) population comprising 336 genotypes observed across three environments was analyzed implementing sparse testing designs. Three prediction models considering main effects (environments, genotypes, genomic) and interaction effects (genotype-by-environment; G×E interaction) were implemented for forecasting dry biomass yield (YDY), total culm (TCM), average internode length (AIL), and culm node number (CNN). Multiple calibration sets based on different compositions and sizes were considered to evaluate performance in terms of the predictive ability (PA) and the mean square error (MSE) for a fixed testing set size. The training set size ranged from 52 to 112 to predict a fixed set of 224 unobserved genotypes across all three environments. The results showed that the model accounting for G×E interaction consistently presented the highest PA and the lowest MSE: for CNN (PA: ~0.77, MSE: ~0.5) and YDY (PA: ~0.70, MSE: ~1.3) while for TCM and AIL these ranged from ~0.28 to 0.41 and ~1.3 to 4.3, respectively. Overall, varying training sets and allocation strategies did not affect PA and MSE, with 52 non-overlapping and 0 overlapping genotypes per environment as the optimal cost-effective allocation framework. This suggests that implementing sparse testing designs could significantly reduce phenotyping costs by fivefold, without compromising PA in breeding programs for perennial crops such as Miscanthus.

Miscanthus sacchariflorus (MSA)↗

Improved bacterial recombineering by parallelized protein discovery

Exploiting bacteriophage-derived homologous recombination processes has offered precise, multiplex editing of microbial genomes and the construction of billions of customized genetic variants in a single day. The techniques that enable this, multiplex automated genome engineering (MAGE) and directed evolution with random genomic mutations (DIvERGE), are however, currently limited to a handful of microorganisms for which single-stranded DNA-annealing proteins (SSAPs) that promote efficient recombineering have been identified. Thus, to enable genome-scale engineering in new hosts, efficient SSAPs must first be found. Here we present a high-throughput method for SSAP discovery that we call “serial enrichment for efficient recombineering” (SEER). By performing SEER in Escherichia coli to screen hundreds of putative SSAPs, we identify highly active variants PapRecT and CspRecT. CspRecT increases the efficiency of single-locus editing to as high as 50% and improves multiplex editing by 5- to 10-fold in E. coli , while PapRecT enables efficient recombineering in Pseudomonas aeruginosa , a concerning human pathogen. CspRecT and PapRecT are also active in other, clinically and biotechnologically relevant enterobacteria. We envision that the deployment of SEER in new species will pave the way toward pooled interrogation of genotype-to-phenotype relationships in previously intractable bacteria.

60 APPLIED LIFE SCIENCES↗

Assessing the Application of a Genomic Network Analysis in Population Ecology: Inferring Patterns of Dispersal and Geographic Structure in the Emerging Pathogen, Coccidioides

A challenge in population ecology studies is identifying how to best group individuals into populations, especially when individual origin is unknown. Machine learning has improved upon traditional methods of identifying population structure and is more efficient at handling large, complex datasets. We demonstrate the applicability of a machine learning method to identify hierarchical population structure in an emerging pathogen, Coccidioides spp., the causative agent of Valley fever. We compared the network clusters to structure identified by traditional tools as a validation of the network performance. We used publicly available whole-genome data for 48 C. immitis and 102 C. posadasii, resulting in 168,211 genome-wide SNPs among the two species. The network analysis grouped samples into populations comparable to the literature for these species but also identified fine-scale geographic structure and travel-associated cases not reported thus far. Exploring different resolutions in the network made it easy to identify unique genotypes specific to California and possibly Nevada, as well as Phoenix- and Tucson-acquired infections in non-endemic areas, regardless of reported travel history. The present study provides a promising example of how a ML-based network analysis can improve our ability to understand pathogen ecology, group cases into populations and infer travel-associated infections.

59 BASIC BIOLOGICAL SCIENCES↗

Improving I/O-aware Workflow Scheduling via Data Flow Characterization and trade-off Analysis

The scientific computing paradigm has transitioned from compute-intensive to I/O-intensive and memory-intensive in the past decade, especially when data-driven science has become common practice. Numerous empirical I/O-aware scheduling optimizations have been developed by incorporating I/O capacity and bandwidth as constraints into scheduling. Unfortunately, there is a lack of data flow (I/O) characterization tool and an understanding of trade-offs between concurrency, locality, and I/O bandwidth. To bridge the gap, this work 1) presents a set of descriptors to characterize, organize, and visualize I/O profiles, including flow size, I/O bandwidth, and operation count, which group data flows by I/O types, tasks, and files; 2) proposes an I/O Roofline model-based trade-off analysis to find the optimal trade-off between flow operational intensity, concurrency, and flow performance. The I/O descriptors generate useful insights into complicated I/O behaviors, suggesting distinct concurrency, storage, and scheduling to be used by types, tasks, and files. The proposed trade-off analysis guides scheduling decisions that generate resource assignment with the best flow parallelism. We evaluate our I/O-aware scheduling methodology on a highly I/O-intensive workflow–1000 Genomes. The experimental results demonstrate speedups of up to 2.4× compared to the state-of-the- art methods.

Guo, Luanzheng [BATTELLE (PACIFIC NW LAB)]↗

A Systems Biology Approach to Energy Flow in H 2 Producing Microbial Communities: Microbial Mats

This report describes the results of work conducted at NASA Ames Research Center, in collaboration with investigators at Lawrence Livermore National Laboratory, during the time period 12/15/2014 to 12/15/2017. The work was an extension of a collaboration with LLNL that began in 2008. Work this reporting period focused primarily on tasks related to microbial mats from both the Elkhorn Slough and Baja California field sites, including the analysis and publication of data produced prior to 12/15/2014. During this reporting period, in addition to the microbial mat work, we also conducted a series of experiments with a particular algal-bacterial association (C. sorokiniana and A. brasilense). A number of aspects of this well-characterized algal-bacterial association enabled us not only to further our understanding of interactions between these particular organisms, but to apply that knowledge and techniques to our collection of microorganisms isolated from microbial mats. Microbial mat experiments focused both on intact mat communities and on cultures of organisms isolated from the mats (the Microbial Mat Living Library). In a series of co-culture experiments, we documented both positive and negative interactions amongst heterotrophic and autotrophic isolates of microbial mat microorganisms. One particularly important and novel mat cyanobacterium (ESFC-1) was studied in detail. We developed techniques to pursue an investigation of functional redundancy (the idea that closely related organisms perform identical roles in a community) in cyanobacterial autotrophs in the Elkhorn Slough microbial mat. In those experiments, we documented numerous instances of increases in growth rate between co-cultures of two organisms (relative to their growth rate in isolation) performing what would seem to be identical functional roles. Simplified microbial mat communities consisting of organisms from the Microbial Mat Living Library were constructed and studied in laboratory experiments. The “constructed mats” were grown in standard laboratory glassware and plasticware, but cultures were also combined with artificial materials in an effort to re-create the three-dimensional structure of naturally occurring microbial mats. These constructed mats re-created a number of functions attributable to natural microbial mats (N-fixation, H2-production) and simplified their study using molecular ecological techniques due to the smaller number of microbial “species” present. Cyanobacteria which migrated in response to vertical variations in light intensity (irradiance) obtained characteristic positions in constructed mats that were similar to their positions in natural mats. A detailed investigation of nitrogen cycling in both the Baja California and Elkhorn Slough microbial mats was initiated during this reporting period. The nitrogen cycle in the microbial mats has not been studied in depth, despite the fact that nitrogen usually limits growth in marine environments. The organisms conducting nitrogen cycling in both microbial mats were studied using the genomes available from earlier work on these mats, as well as with a targeted approach using molecular ecological methods developed here. Quantitative measurements of rates of nitrogen transformations were made using stable isotopically labeled nitrogen. The processes quantified and examined using molecular ecological methods included: nitrogen fixation, ammonification, ammonium and nitrate assimilation, denitrification, anaerobic oxidation of ammonium, and dissimilatory nitrate reduction to ammonium. In a series of experiments designed to examine a potential application for microbial mats, we demonstrated that both Elkhorn Slough and Baja California microbial mats were very effective in removing large quantities of both nitrate and ammonium from waste streams containing these compounds at levels where they would be considered to be pollutants. Considerable effort was expended during this reporting period to analyze, and prepare for publication, metagenomes and metatranscriptomes collected from two experimental manipulations of microbial mats. The first, a diel study, was conducted in November 2011, while the second, an investigation of the effects of sulfate on community composition, was conducted during this reporting period. For the earlier study, reference-based and reference-free methods were used to assess microbial ecology and genetic partitioning in these complex microbial systems, and that work is available as a pre-print. Four metagenomes from the same experiment were published. In the second, longer (nearly 200 day) manipulation, rRNA iTag libraries for cDNA and DNA were analyzed to determine the effects of sulfate concentration on microbial mat community composition over the course of the experiment. Work on the algal-bacterial association (Chlorella sorokiniana and Azospirillum brasilense) focused initially on demonstrating the transfer of carbon and nitrogen between the partners, using stable isotopic labeling and the LLNL NanoSIMS. Subsequently, the potential role of diazotrophy (nitrogen fixation) by the bacterial partner in supporting the growth of the association was assessed in a series of experiments. It appears to be the case that, although the bacterial partner demonstrably increases the fitness of the algae under nitrogen limitation, and cells of C. sorokiniana receive nitrogen originally fixed from atmospheric nitrogen by A. brasilense, the amount of nitrogen obtained by diazotrophy is too small to explain the increase in fitness.

08 HYDROGEN↗

GENESPACE tracks regions of interest and gene copy number variation across multiple genomes

The development of multiple chromosome-scale reference genome sequences in many taxonomic groups has yielded a high-resolution view of the patterns and processes of molecular evolution. Nonetheless, leveraging information across multiple genomes remains a significant challenge in nearly all eukaryotic systems. These challenges range from studying the evolution of chromosome structure, to finding candidate genes for quantitative trait loci, to testing hypotheses about speciation and adaptation. Here, we present GENESPACE, which addresses these challenges by integrating conserved gene order and orthology to define the expected physical position of all genes across multiple genomes. We demonstrate this utility by dissecting presence–absence, copy-number, and structural variation at three levels of biological organization: spanning 300 million years of vertebrate sex chromosome evolution, across the diversity of the Poaceae (grass) plant family, and among 26 maize cultivars. The methods to build and visualize syntenic orthology in the GENESPACE R package offer a significant addition to existing gene family and synteny programs, especially in polyploid, outbred, and other complex genomes.

59 BASIC BIOLOGICAL SCIENCES↗