Improved quality metrics for association and reproducibility in chromatin accessibility data using mutual information.

BACKGROUND: Correlation metrics are widely utilized in genomics analysis and often implemented with little regard to assumptions of normality, homoscedasticity, and independence of values. This is especially true when comparing values between replicated sequencing experiments that probe chromatin accessibility, such as assays for transposase-accessible chromatin via sequencing (ATAC-seq). Such data can possess several regions across the human genome with little to no sequencing depth and are thus non-normal with a large portion of zero values. Despite distributed use in the epigenomics field, few studies have evaluated and benchmarked how correlation and association statistics behave across ATAC-seq experiments with known differences or the effects of removing specific outliers from the data. Here, we developed a computational simulation of ATAC-seq data to elucidate the behavior of correlation statistics and to compare their accuracy under set conditions of reproducibility. RESULTS: Using these simulations, we monitored the behavior of several correlation statistics, including the Pearson's R and Spearman's [Formula: see text] coefficients as well as Kendall's [Formula: see text] and Top-Down correlation. We also test the behavior of association measures, including the coefficient of determination R[Formula: see text], Kendall's W, and normalized mutual information. Our experiments reveal an insensitivity of most statistics, including Spearman's [Formula: see text], Kendall's [Formula: see text], and Kendall's W, to increasing differences between simulated ATAC-seq replicates. The removal of co-zeros (regions lacking mapped sequenced reads) between simulated experiments greatly improves the estimates of correlation and association. After removing co-zeros, the R[Formula: see text] coefficient and normalized mutual information display the best performance, having a closer one-to-one relationship with the known portion of shared, enhanced loci between simulated replicates. When comparing values between experimental ATAC-seq data using a random forest model, mutual information best predicts ATAC-seq replicate relationships. CONCLUSIONS: Collectively, this study demonstrates how measures of correlation and association can behave in epigenomics experiments. We provide improved strategies for quantifying relationships in these increasingly prevalent and important chromatin accessibility assays.

Nonclassical autophagy activation pathways are essential for production of infectious Influenza A virus in vitro.

Autophagy is a critical mechanism deployed by eukaryotic cells in response to stress, including viral infection, to boost the innate antimicrobial responses. However, an increasing number of pathogens hijack the autophagic machinery to facilitate their own replication. Influenza A virus (IAV), responsible for several global pandemics, has an intricate dependence on autophagy for successful replication in mammalian cells. To elucidate key chokepoints in the host stress responses facilitating IAV replication, we constructed a meta-transcriptome of IAV and host gene expression dynamics during early (1-3 hpi), mid (4-6 hpi), and late (8-12 hpi) stages of the viral replication cycle at two multiplicities of infection (MOI): 1 and 5. We supplemented the global transcriptome study with phosphoproteomic analysis of stress-activated protein kinase (SAPK/JNK) signaling in lung carcinoma (predominantly used as an in vitro model of IAV replication) and normal human bronchial epithelial cells. We report significant differences in the activation profiles of autophagy regulating genes upon IAV infection at the two MOI as well as divergent dependence on ULK1 signaling within the normal and cancer cells. Regardless of the cell model, JNK-Thr187 signaling was crucial for the production of infectious viral particles.

RNA-Binding Protein RBP-P Is Required for Glutelin and Prolamine mRNA Localization in Rice Endosperm Cells.

In developing rice (Oryza sativa) endosperm, mRNAs of the major storage proteins, glutelin and prolamine, are transported and anchored to distinct subdomains of the cortical endoplasmic reticulum. RNA binding protein RBP-P binds to both glutelin and prolamine mRNAs, suggesting a role in some aspect of their RNA metabolism. Here, we show that rice lines expressing mutant RBP-P mislocalize both glutelin and prolamine mRNAs. Different mutant RBP-P proteins exhibited varying degrees of reduced RNA binding and/or protein-protein interaction properties, which may account for the mislocalization of storage protein RNAs. In addition, partial loss of RBP-P function conferred a broad phenotypic variation ranging from dwarfism, chlorophyll deficiency, and sterility to late flowering and low spikelet fertility. Transcriptome analysis highlighted the essential role of RBP-P in regulating storage protein genes and several essential biological processes during grain development. Overall, our data demonstrate the significant roles of RBP-P in glutelin and prolamine mRNA localization and in the regulation of genes important for plant growth and development through its RNA binding activity and cooperative regulation with interacting proteins.

Selective sets of mRNAs localize to extracellular paramural bodies in a rice glup6 mutant.

The transport of rice glutelin storage proteins to the storage vacuoles requires the Rab5 GTPase and its related guanine nucleotide exchange factor (Rab5-GEF). Loss of function of these membrane vesicular trafficking factors results in the initial secretion of storage proteins and later their partial engulfment by the plasma membrane to form an extracellular paramural body (PMB), an aborted endosome complex. Here, we show that in the rice Rab5-GEF mutant glup6, glutelin RNAs are specifically mislocalized from their normal location on the cisternal endoplasmic reticulum (ER) to the protein body-ER, and are also apparently translocated to the PMBs. We substantiated the association of mRNAs with this aborted endosome complex by RNA-seq of PMBs purified by flow cytometry. Two PMB-associated groups of RNA were readily resolved: those that were specifically enriched in this aborted complex and those that were highly expressed in the cytoplasm. Examination of the PMB-enriched RNAs indicated that they were not a random sampling of the glup6 transcriptome but, instead, encompassed only a few functional mRNA classes. Although specific autophagy is also an alternative mechanism, our results support the view that RNA localization may co-opt membrane vesicular trafficking, and that many RNAs that share function or intracellular location are co-transported in developing rice seeds.

Genome Sequence and Transcriptome Analyses of Chrysochromulina tobin: Metabolic Tools for Enhanced Algal Fitness in the Prominent Order Prymnesiales (Haptophyceae).

Haptophytes are recognized as seminal players in aquatic ecosystem function. These algae are important in global carbon sequestration, form destructive harmful blooms, and given their rich fatty acid content, serve as a highly nutritive food source to a broad range of eco-cohorts. Haptophyte dominance in both fresh and marine waters is supported by the mixotrophic nature of many taxa. Despite their importance the nuclear genome sequence of only one haptophyte, Emiliania huxleyi (Isochrysidales), is available. Here we report the draft genome sequence of Chrysochromulina tobin (Prymnesiales), and transcriptome data collected at seven time points over a 24-hour light/dark cycle. The nuclear genome of C. tobin is small (59 Mb), compact (∼ 40% of the genome is protein coding) and encodes approximately 16,777 genes. Genes important to fatty acid synthesis, modification, and catabolism show distinct patterns of expression when monitored over the circadian photoperiod. The C. tobin genome harbors the first hybrid polyketide synthase/non-ribosomal peptide synthase gene complex reported for an algal species, and encodes potential anti-microbial peptides and proteins involved in multidrug and toxic compound extrusion. A new haptophyte xanthorhodopsin was also identified, together with two "red" RuBisCO activases that are shared across many algal lineages. The Chrysochromulina tobin genome sequence provides new information on the evolutionary history, ecology and economic importance of haptophytes.

A pangenomic analysis of the Nannochloropsis organellar genomes reveals novel genetic variations in key metabolic genes.

BACKGROUND: Microalgae in the genus Nannochloropsis are photosynthetic marine Eustigmatophytes of significant interest to the bioenergy and aquaculture sectors due to their ability to efficiently accumulate biomass and lipids for utilization in renewable transportation fuels, aquaculture feed, and other useful bioproducts. To better understand the genetic complement that drives the metabolic processes of these organisms, we present the assembly and comparative pangenomic analysis of the chloroplast and mitochondrial genomes from Nannochloropsis salina CCMP1776. RESULTS: The chloroplast and mitochondrial genomes of N. salina are 98.4% and 97% identical to their counterparts in Nannochloropsis gaditana. Comparison of the Nannochloropsis pangenome to other algae within and outside of the same phyla revealed regions of significant genetic divergence in key genes that encode proteins needed for regulation of branched chain amino synthesis (acetohydroxyacid synthase), carbon fixation (RuBisCO activase), energy conservation (ATP synthase), protein synthesis and homeostasis (Clp protease, ribosome). CONCLUSIONS: Many organellar gene modifications in Nannochloropsis are unique and deviate from conserved orthologs found across the tree of life. Implementation of secondary and tertiary structure prediction was crucial to functionally characterize many proteins and therefore should be implemented in automated annotation pipelines. The exceptional similarity of the N. salina and N. gaditana organellar genomes suggests that N. gaditana be reclassified as a strain of N. salina.

The mitochondrial and chloroplast genomes of the haptophyte Chrysochromulina tobin contain unique repeat structures and gene profiles.

BACKGROUND: Haptophytes are widely and abundantly distributed in both marine and freshwater ecosystems. Few genomic analyses of representatives within this taxon have been reported, despite their early evolutionary origins and their prominent role in global carbon fixation. RESULTS: The complete mitochondrial and chloroplast genome sequences of the haptophyte Chrysochromulina tobin (Prymnesiales) provide insight into the architecture and gene content of haptophyte organellar genomes. The mitochondrial genome (~34 kb) encodes 21 protein coding genes and contains a complex, 9 kb tandem repeat region. Similar to other haptophytes and rhodophytes, but not cryptophytes or stramenopiles, the mitochondrial genome has lost the nad7, nad9 and nad11 genes. The ~105 kb chloroplast genome encodes 112 protein coding genes, including ycf39 which has strong structural homology to NADP-binding nitrate transcriptional regulators; a divergent 'CheY-like' two-component response regulator (ycf55) and Tic/Toc (ycf60 and ycf80) membrane transporters. Notably, a zinc finger domain has been identified in the rpl36 ribosomal protein gene of all chloroplasts sequenced to date with the exception of haptophytes and cryptophytes--algae that have gained (via lateral gene transfer) an alternative rpl36 lacking the zinc finger motif. The two C. tobin chloroplast ribosomal RNA operon spacer regions differ in tRNA content. Additionally, each ribosomal operon contains multiple single nucleotide polymorphisms (SNPs)--a pattern observed in rhodophytes and cryptophytes, but few stramenopiles. Analysis of small (<200 bp) chloroplast encoded tandem and inverted repeats in C. tobin and 78 other algal chloroplast genomes show that repeat type, size and location are correlated with gene identity and taxonomic clade. CONCLUSION: The Chrysochromulina tobin organellar genomes provide new insight into organellar function and evolution. These are the first organellar genomes to be determined for the prymnesiales, a taxon that is present in both oceanic and freshwater systems and represents major primary photosynthetic producers and contributors to global ecosystem stability.

Multi-omics analysis reveals the dynamic interplay between Vero host chromatin structure and function during vaccinia virus infection.

The genome folds into complex configurations and structures thought to profoundly impact its function. The intricacies of this dynamic structure-function relationship are not well understood particularly in the context of viral infection. To unravel this interplay, here we provide a comprehensive investigation of simultaneous host chromatin structural (via Hi-C and ATAC-seq) and functional changes (via RNA-seq) in response to vaccinia virus infection. Over time, infection significantly impacts global and local chromatin structure by increasing long-range intra-chromosomal interactions and B compartmentalization and by decreasing chromatin accessibility and inter-chromosomal interactions. Local accessibility changes are independent of broad-scale chromatin compartment exchange (~12% of the genome), underscoring potential independent mechanisms for global and local chromatin reorganization. While infection structurally condenses the host genome, there is nearly equal bidirectional differential gene expression. Despite global weakening of intra-TAD interactions, functional changes including downregulated immunity genes are associated with alterations in local accessibility and loop domain restructuring. Therefore, chromatin accessibility and local structure profiling provide impactful predictions for host responses and may improve development of efficacious anti-viral counter measures including the optimization of vaccine design.

Human Coronavirus Infection Reorganizes Spatial Genomic Architecture in Permissive Lung Cells.

Chromatin conformation capture followed by next-generation sequencing in combination with large-scale polymer simulations (4DHiC) produces detailed information on genomic loci interactions, allowing for the interrogation of 3D spatial genomic structures. Here, Hi-C data was acquired from the infection of fetal lung fibroblast (MRC5) cells with α-coronavirus 229E (CoV229E). Experimental Hi-C contact maps were used to determine viral-induced changes in genomic architecture over a 48-hour time period following viral infection, revealing substantial alterations in contacts within chromosomes and in contacts between different chromosomes. To gain further structural insight and quantify the underlying changes, we applied the 4DHiC polymer simulation method to reconstruct the 3D genomic structures and dynamics corresponding to the Hi-C maps. The models successfully reproduced experimental Hi-C data, including the changes in contacts induced by viral infection. Our 3D spatial simulations uncovered widespread chromatin restructuring, including increased chromosome compactness and A-B compartment mixing arising from infection. Our model also suggests increased spatial accessibility to regions containing interferon-stimulated genes upon infection with CoV229E, followed by chromatin restructuring at later time points, potentially inducing the migration of chromatin into more compact regions. This is consistent with previously observed suppression of gene expression. Our spatial genomics study provides a mechanistic structural basis for changes in chromosome architecture induced by coronavirus infection in lung cells.

Trait drift in microalgae and applications for strain improvement.

Microalgae are increasingly used to generate a wide range of commercial products, and there is growing evidence that microalgae-based products can be produced sustainably. However, industrial production of microalgal biomass is not as developed as other biomanufacturing platform technologies. In addition, results of bench-scale research often fail to translate to large-scale or mass production systems. This disconnect may result from trait drift and evolution occurring, through time, in response to unique drivers in each environment, such as cultivation regimes, weather, and pests. Moreover, outdoor and indoor cultivation of microalgae has the potential to impose negative selection pressures, which makes the maintenance of desired traits a challenge. In this context, this review sheds the light on our current understanding of trait drift and evolution in microalgae. We delineate the basics of phenotype plasticity and evolution, with a focus on how microalgae respond under various conditions. In addition, we review techniques that exploit phenotypic plasticity and evolution for strain improvement in view of industrial commercial applications, highlighting associated advantages and shortcomings. Finally, we suggest future research directions and recommendations to overcome unwanted trait drift and evolution in microalgae cultivation.

Metabolism of Scenedesmus obliquus cultivated with raw plant substrates.

The potential benefits of adding raw, non-food, lignocellulosic plant material as a carbon source for mixotrophic growth of microalgae have previously been demonstrated. This approach has advantages over using traditional carbon sources like glucose or acetate due to wide-spread plant biomass availability and substrate recalcitrance to bacterial contamination. Here, we report the overall growth characteristics and explore the metabolic patterns of Scenedesmus obliquus cultured in the presence raw plant substrate. An initial screen of plant substrate candidates showed an increase in specific growth rate and biomass accumulation when S. obliquus was cultured in the presence of switchgrass or yard waste compared to media alone. We observed a near doubling of microalgal dry weight when S. obliquus was grown with 0.2% (w/v) switchgrass under ambient CO2. Scanning electron microscopy (SEM) of corn stem after S. obliquus cultivation exhibited substantial phloem degradation. Transcriptomic analyses of S. obliquus during mid- and late-log phase growth revealed a dynamic metabolic landscape within many KEGG pathways. Notably, differential expression was observed for several potential glycosyl hydrolases. We also investigated the influence of switchgrass on the growth of S. obliquus at 50 L volume in mini raceway ponds to determine the scalability of this approach.

The Genome Sequence of an Algal Strain of Nannochloropsis QH25.

Species of Nannochloropsis are single-celled Stramenopiles commonly used in microalgae-based technologies for the manufacturing of bioproducts. Nannochloropsis oceanica QH25 was isolated from an algal cultivation pond located in Imperial, Texas (USA). We used PacBio continuous long read (CLR) sequencing to produce a highly contiguous 29.34 Mb genome.

Genome of the cyanobacterium Microcoleus vaginatus FGP-2, a photosynthetic ecosystem engineer of arid land soil biocrusts worldwide.

The filamentous cyanobacterium Microcoleus vaginatusis found in arid land soils worldwide. The genome of M. vaginatus strain FGP-2 allows exploration of genes involved in photosynthesis, desiccation tolerance, alkane production, and other features contributing to this organism's ability to function as a major component of biological soil crusts in arid lands.

A multi-omic characterization of temperature stress in a halotolerant Scenedesmus strain for algal biotechnology.

Microalgae efficiently convert sunlight into lipids and carbohydrates, offering bio-based alternatives for energy and chemical production. Improving algal productivity and robustness against abiotic stress requires a systems level characterization enabled by functional genomics. Here, we characterize a halotolerant microalga Scenedesmus sp. NREL 46B-D3 demonstrating peak growth near 25 °C that reaches 30 g/m2/day and the highest biomass accumulation capacity post cell division reported to date for a halotolerant strain. Functional genomics analysis revealed that genes involved in lipid production, ion channels and antiporters are expanded and expressed. Exposure to temperature stress shifts fatty acid metabolism and increases amino acids synthesis. Co-expression analysis shows that many fatty acid biosynthesis genes are overexpressed with specific transcription factors under cold stress. These and other genes involved in the metabolic and regulatory response to temperature stress can be further explored for strain improvement.

High-throughput phenotyping of cell-to-cell interactions in gel microdroplet pico-cultures.

Microbiomes exert significant influence on our planet's ecology. Elucidating the identities of individual microbes within these communities and how they interact is a vital research imperative. Using traditional plating and culturing methods, it is impractical to assess even a small fraction of the interactions that exist within microbial communities. To address this technology gap, we integrated gel microdroplet technology with microfluidics to generate millions of microdroplet cultures (MDs) that sequester individual cells for phenotyping MDs, facilitating rapid analysis and viable recovery using flow cytometry. Herein, we describe a validated high-throughput phenotyping pipeline that elucidates cell-to-cell interactions for millions of combinations of microorganisms. Through iterative co-culturing of an algae and a pool of environmentally sourced microbes, we successfully isolated bacteria that improved algal growth.

High-Quality Draft Genome Sequence of the Green Alga Tetraselmis striata (Chlorophyta) Generated from PacBio Sequencing.

A high-quality draft genome sequence of the microalgal species Tetraselmis striata was generated using PacBio sequencing. The assembled genome is 228 Mb, derived from 3,613 polished contigs at 84× coverage depth. This genome contains an average GC content of 57.9% and 48,906 predicted genes.

Development of a high-productivity, halophilic, thermotolerant microalga Picochlorum renovo.

Microalgae are promising biocatalysts for applications in sustainable fuel, food, and chemical production. Here, we describe culture collection screening, down-selection, and development of a high-productivity, halophilic, thermotolerant microalga, Picochlorum renovo. This microalga displays a rapid growth rate and high diel biomass productivity (34 g m-2 day-1), with a composition well-suited for downstream processing. P. renovo exhibits broad salinity tolerance (growth at 107.5 g L-1 salinity) and thermotolerance (growth up to 40 °C), beneficial traits for outdoor cultivation. We report complete genome sequencing and analysis, and genetic tool development suitable for expression of transgenes inserted into the nuclear or chloroplast genomes. We further evaluate mechanisms of halotolerance via comparative transcriptomics, identifying novel genes differentially regulated in response to high salinity cultivation. These findings will enable basic science inquiries into control mechanisms governing Picochlorum biology and lay the foundation for development of a microalga with industrially relevant traits as a model photobiology platform.

Detection of Abrin-Like and Prepropulchellin-Like Toxin Genes and Transcripts Using Whole Genome Sequencing and Full-Length Transcript Sequencing of Abrus precatorius.

The sequenced genome and the leaf transcriptome of a near relative of Abrus pulchellus and Abrus precatorius was analyzed to characterize the genetic basis of toxin gene expression. From the high-quality genome assembly, a total of 26 potential coding regions were identified that contain genes with abrin-like, pulchellin-like, and agglutinin-like homology, with full-length transcripts detected in leaf tissue for 9 of the 26 coding regions. All of the toxin-like genes were identified within only five isolated regions of the genome, with each region containing 1 to 16 gene variants within each genomic region (<1 Mbp). The Abrusprecatorius cultivar sequenced here contains genes which encode for proteins that are homologous to certain abrin and prepropulchellin genes previously identified, and we observed substantial diversity of genes and predicted gene products in Abrus precatorius and previously characterized toxins. This suggests diverse toxin repertoires within Abrus, potentially the results of rapid toxin evolution.

Gut microbiota and their putative metabolic functions in fragmented Bengal tiger population of Nepal.

Bengal tigers (Panthera tigris tigris) serve a pivotal role as an apex predator in forest ecosystems. To increase our knowledge on factors impacting the viability and health of this endangered species, we studied the gut microbiota in 32 individual Bengal tigers from three geographically separated areas (Chitwan National Park (CNP), Bardia National Park (BNP) and Suklaphanta Wildlife Reserve (SWR)) in Nepal, using noninvasive genetic sampling methods. Gut microbiota influence the immune system, impact various physiological functions, and modulates metabolic reactions, that ultimately impact the host health, behavior and development. Across the tiger populations in Nepal, we found significant differences in the composition of microbial communities based on their geographic locations. Specifically, we detected significant differences between CNP and the other two protected areas (CNP vs BNP: pseudo t = 1.944, P = 0.006; CNP vs SWR: pseudo t = 1.9942, P = 0.0071), but no differences between BNP and SWR. This mirrors what has been found for tiger gene flow in the same populations, suggesting gut microbiota composition and host gene flow may be linked. Furthermore, predictive metagenome functional content analysis (PICRUSt) revealed a higher functional enrichment and diversity for significant gut microbiota in the Chitwan tiger population and the lowest enrichment and diversity in Suklaphanta. The CNP tiger population contained higher proportions of microbiota that are associated with predicted functions relevant for metabolism of amino acid, lipid, xenobiotics biodegradation, terpenoides and polyketides than the SWR population. We conclude the tiger population structure, gut microbiota profile and associated functional metabolic categories are correlated, with geographically most separated CNP and SWR tiger population having the most distinct and different host genotype and microbiota profiles. Our work dramatically expands the understanding of tiger microbiota in wild populations and provides a valuable case study on how to investigate genetic diversity at different hierarchical levels, including hosts as well as their microbial communities.

Expression of a putative nitrite reductase and the reversible inhibition of nitrite-dependent respiration by nitric oxide in Nitrobacter winogradskyi Nb-255.

The nitrite oxidizing Alphaproteobacterium, Nitrobacter winogradskyi, primarily conserves energy from the oxidation of nitrite (NO(2)(-))to nitrate (NO(3)(-)) through aerobic respiration. Almost 20 years ago, NO-dependent NADH formation was reported to occur in both aerobic and anaerobic cell suspensions of N. winogradskyi strain 'agilis', suggesting that NO oxidation might contribute to energy conservation by Nitrobacter. Recently, the N. winogradskyi Nb-255 genome was found to contain a gene (Nwin_2648) that encodes a putative copper-containing nitrite reductase (NirK), which may reduce NO(2)(-) to NO. In this study, the putative nirK was found to be maximally transcribed under low O(2) (between zero and 4% O(2)) in the presence of NO(2)(-). Transcription of nirK was not detected under anaerobic conditions in the absence of NO(2)(-) or in the presence of NO(3)(-) and pyruvate. Although net production of NO could not be detected from either aerobically grown or anaerobically incubated cells, exogenous NO was consumed by viable cells and concomitantly inhibited NO(2)(-)-dependent O(2) uptake in a reversible, concentration dependent manner. Both NO(2(-)-dependent O(2) uptake and NO consumption were inhibited by 1 mM cyanide suggesting involvement of cytochrome oxidase with NO consumption. Abiotic consumption of NO was measured, yet, both the rates and kinetics of NO transformation in buffer alone, or by heat killed, or cyanide-treated cells differed from those of viable cells. In light of this new information, a modified model is proposed to explain how NirK and NO manage electron flux in Nitrobacter.