A High-Quality Blue Whale Genome, Segmental Duplications, and Historical Demography.

The blue whale, Balaenoptera musculus, is the largest animal known to have ever existed, making it an important case study in longevity and resistance to cancer. To further this and other blue whale-related research, we report a reference-quality, long-read-based genome assembly of this fascinating species. We assembled the genome from PacBio long reads and utilized Illumina/10×, optical maps, and Hi-C data for scaffolding, polishing, and manual curation. We also provided long read RNA-seq data to facilitate the annotation of the assembly by NCBI and Ensembl. Additionally, we annotated both haplotypes using TOGA and measured the genome size by flow cytometry. We then compared the blue whale genome with other cetaceans and artiodactyls, including vaquita (Phocoena sinus), the world's smallest cetacean, to investigate blue whale's unique biological traits. We found a dramatic amplification of several genes in the blue whale genome resulting from a recent burst in segmental duplications, though the possible connection between this amplification and giant body size requires further study. We also discovered sites in the insulin-like growth factor-1 gene correlated with body size in cetaceans. Finally, using our assembly to examine the heterozygosity and historical demography of Pacific and Atlantic blue whale populations, we found that the genomes of both populations are highly heterozygous and that their genetic isolation dates to the last interglacial period. Taken together, these results indicate how a high-quality, annotated blue whale genome will serve as an important resource for biology, evolution, and conservation research.

Chromosome level genome assembly of the Etruscan shrew Suncus etruscus.

Suncus etruscus is one of the world's smallest mammals, with an average body mass of about 2 grams. The Etruscan shrew's small body is accompanied by a very high energy demand and numerous metabolic adaptations. Here we report a chromosome-level genome assembly using PacBio long read sequencing, 10X Genomics linked short reads, optical mapping, and Hi-C linked reads. The assembly is partially phased, with the 2.472 Gbp primary pseudohaplotype and 1.515 Gbp alternate. We manually curated the primary assembly and identified 22 chromosomes, including X and Y sex chromosomes. The NCBI genome annotation pipeline identified 39,091 genes, 19,819 of them protein-coding. We also identified segmental duplications, inferred GO term annotations, and computed orthologs of human and mouse genes. This reference-quality genome will be an important resource for research on mammalian development, metabolism, and body size control.

Decreased plasma cartilage acidic protein 1 in COVID-19.

Cartilage acidic protein-1 (CRTAC1) is produced by several cell types, including Type 2 alveolar epithelial (T2AE) cells that are targeted by SARS-CoV2. Plasma CRTAC1 is known based on proteomic surveys to be low in patients with severe COVID-19. Using an ELISA, we found that patients treated for COVID-19 in an ICU almost uniformly had plasma concentrations of CRTAC1 below those of healthy controls. Magnitude of decrease in CRTAC1 distinguished COVID-19 from other causes of acute respiratory decompensation and correlated with established metrics of COVID-19 severity. CRTAC1 concentrations below those of controls were found in some patients a year after hospitalization with COVID-19, long COVID after less severe COVID-19, or chronic obstructive pulmonary disease. Decreases in CRTAC1 in severe COVID-19 correlated (r = 0.37, p = 0.0001) with decreases in CFP (properdin), which interacts with CRTAC1. Thus, decreases of CRTAC1 associated with severe COVID-19 may result from loss of production by T2AE cells or co-depletion with CFP. Determination of significance of and reasons behind decreased CRTAC1 concentration in a subset of patients with long COVID will require analysis of roles of preexisting lung disease, impact of prior acute COVID-19, age, and other confounding variables in a larger number of patients.

A haplotype-resolved genome assembly of the Nile rat facilitates exploration of the genetic basis of diabetes.

BACKGROUND: The Nile rat (Avicanthis niloticus) is an important animal model because of its robust diurnal rhythm, a cone-rich retina, and a propensity to develop diet-induced diabetes without chemical or genetic modifications. A closer similarity to humans in these aspects, compared to the widely used Mus musculus and Rattus norvegicus models, holds the promise of better translation of research findings to the clinic. RESULTS: We report a 2.5 Gb, chromosome-level reference genome assembly with fully resolved parental haplotypes, generated with the Vertebrate Genomes Project (VGP). The assembly is highly contiguous, with contig N50 of 11.1 Mb, scaffold N50 of 83 Mb, and 95.2% of the sequence assigned to chromosomes. We used a novel workflow to identify 3613 segmental duplications and quantify duplicated genes. Comparative analyses revealed unique genomic features of the Nile rat, including some that affect genes associated with type 2 diabetes and metabolic dysfunctions. We discuss 14 genes that are heterozygous in the Nile rat or highly diverged from the house mouse. CONCLUSIONS: Our findings reflect the exceptional level of genomic resolution present in this assembly, which will greatly expand the potential of the Nile rat as a model organism.

Reference genome and demographic history of the most endangered marine mammal, the vaquita.

The vaquita is the most critically endangered marine mammal, with fewer than 19 remaining in the wild. First described in 1958, the vaquita has been in rapid decline for more than 20 years resulting from inadvertent deaths due to the increasing use of large-mesh gillnets. To understand the evolutionary and demographic history of the vaquita, we used combined long-read sequencing and long-range scaffolding methods with long- and short-read RNA sequencing to generate a near error-free annotated reference genome assembly from cell lines derived from a female individual. The genome assembly consists of 99.92% of the assembled sequence contained in 21 nearly gapless chromosome-length autosome scaffolds and the X-chromosome scaffold, with a scaffold N50 of 115 Mb. Genome-wide heterozygosity is the lowest (0.01%) of any mammalian species analysed to date, but heterozygosity is evenly distributed across the chromosomes, consistent with long-term small population size at genetic equilibrium, rather than low diversity resulting from a recent population bottleneck or inbreeding. Historical demography of the vaquita indicates long-term population stability at less than 5,000 (Ne) for over 200,000 years. Together, these analyses indicate that the vaquita genome has had ample opportunity to purge highly deleterious alleles and potentially maintain diversity necessary for population health.

Effects of PHENYLALANINE AMMONIA LYASE (PAL) knockdown on cell wall composition, biomass digestibility, and biotic and abiotic stress responses in Brachypodium.

The phenylpropanoid pathway in plants synthesizes a variety of structural and defence compounds, and is an important target in efforts to reduce cell wall lignin for improved biomass conversion to biofuels. Little is known concerning the trade-offs in grasses when perturbing the function of the first gene family in the pathway, PHENYLALANINE AMMONIA LYASE (PAL). Therefore, PAL isoforms in the model grass Brachypodium distachyon were targeted, by RNA interference (RNAi), and large reductions (up to 85%) in stem tissue transcript abundance for two of the eight putative BdPAL genes were identified. The cell walls of stems of BdPAL-knockdown plants had reductions of 43% in lignin and 57% in cell wall-bound ferulate, and a nearly 2-fold increase in the amounts of polysaccharide-derived carbohydrates released by thermochemical and hydrolytic enzymic partial digestion. PAL-knockdown plants exhibited delayed development and reduced root growth, along with increased susceptibilities to the fungal pathogens Fusarium culmorum and Magnaporthe oryzae. Surprisingly, these plants generally had wild-type (WT) resistances to caterpillar herbivory, drought, and ultraviolet light. RNA sequencing analyses revealed that the expression of genes associated with stress responses including ethylene biosynthesis and signalling were significantly altered in PAL knocked-down plants under non-challenging conditions. These data reveal that, although an attenuation of the phenylpropanoid pathway increases carbohydrate availability for biofuel, it can adversely affect plant growth and disease resistance to fungal pathogens. The data identify notable differences between the stress responses of these monocot pal mutants versus Arabidopsis (a dicot) pal mutants and provide insights into the challenges that may arise when deploying phenylpropanoid pathway-altered bioenergy crops.

Dicholine succinate, the neuronal insulin sensitizer, normalizes behavior, REM sleep, hippocampal pGSK3 beta and mRNAs of NMDA receptor subunits in mouse models of depression.

Central insulin receptor-mediated signaling is attracting the growing attention of researchers because of rapidly accumulating evidence implicating it in the mechanisms of plasticity, stress response, and neuropsychiatric disorders including depression. Dicholine succinate (DS), a mitochondrial complex II substrate, was shown to enhance insulin-receptor mediated signaling in neurons and is regarded as a sensitizer of the neuronal insulin receptor. Compounds enhancing neuronal insulin receptor-mediated transmission exert an antidepressant-like effect in several pre-clinical paradigms of depression; similarly, such properties for DS were found with a stress-induced anhedonia model. Here, we additionally studied the effects of DS on several variables which were ameliorated by other insulin receptor sensitizers in mice. Pre-treatment with DS of chronically stressed C57BL6 mice rescued normal contextual fear conditioning, hippocampal gene expression of NMDA receptor subunit NR2A, the NR2A/NR2B ratio and increased REM sleep rebound after acute predation. In 18-month-old C57BL6 mice, a model of elderly depression, DS restored normal sucrose preference and activated the expression of neural plasticity factors in the hippocampus as shown by Illumina microarray. Finally, young naïve DS-treated C57BL6 mice had reduced depressive- and anxiety-like behaviors and, similarly to imipramine-treated mice, preserved hippocampal levels of the phosphorylated (inactive) form of GSK3 beta that was lowered by forced swimming in pharmacologically naïve animals. Thus, DS can ameliorate behavioral and molecular outcomes under a variety of stress- and depression-related conditions. This further highlights neuronal insulin signaling as a new factor of pathogenesis and a potential pharmacotherapy of affective pathologies.

Plant-derived antifungal agent poacic acid targets ß-1,3-glucan.

A rise in resistance to current antifungals necessitates strategies to identify alternative sources of effective fungicides. We report the discovery of poacic acid, a potent antifungal compound found in lignocellulosic hydrolysates of grasses. Chemical genomics using Saccharomyces cerevisiae showed that loss of cell wall synthesis and maintenance genes conferred increased sensitivity to poacic acid. Morphological analysis revealed that cells treated with poacic acid behaved similarly to cells treated with other cell wall-targeting drugs and mutants with deletions in genes involved in processes related to cell wall biogenesis. Poacic acid causes rapid cell lysis and is synergistic with caspofungin and fluconazole. The cellular target was identified; poacic acid localized to the cell wall and inhibited ß-1,3-glucan synthesis in vivo and in vitro, apparently by directly binding ß-1,3-glucan. Through its activity on the glucan layer, poacic acid inhibits growth of the fungi Sclerotinia sclerotiorum and Alternaria solani as well as the oomycete Phytophthora sojae. A single application of poacic acid to leaves infected with the broad-range fungal pathogen S. sclerotiorum substantially reduced lesion development. The discovery of poacic acid as a natural antifungal agent targeting ß-1,3-glucan highlights the potential side use of products generated in the processing of renewable biomass toward biofuels as a source of valuable bioactive compounds and further clarifies the nature and mechanism of fermentation inhibitors found in lignocellulosic hydrolysates.

Systems biology defines the biological significance of redox-active proteins during cellulose degradation in an aerobic bacterium.

Microbial depolymerization of plant cell walls contributes to global carbon balance and is a critical component of renewable energy. The genomes of lignocellulose degrading microorganisms encode diverse classes of carbohydrate modifying enzymes, although currently there is a paucity of knowledge on the role of these proteins in vivo. We report the comprehensive analysis of the cellulose degradation system in the saprophytic bacterium Cellvibrio japonicus. Gene expression profiling of C. japonicus demonstrated that three of the 12 predicted ß-1,4 endoglucanases (cel5A, cel5B, and cel45A) and the sole predicted cellobiohydrolase (cel6A) showed elevated expression during growth on cellulose. Targeted gene disruptions of all 13 predicted cellulase genes showed that only cel5B and cel6A were required for optimal growth on cellulose. Our analysis also identified three additional genes required for cellulose degradation: lpmo10B encodes a lytic polysaccharide monooxygenase (LPMO), while cbp2D and cbp2E encode proteins containing carbohydrate binding modules and predicted cytochrome domains for electron transfer. CjLPMO10B oxidized cellulose and Cbp2D demonstrated spectral properties consistent with redox function. Collectively, this report provides insight into the biological role of LPMOs and redox proteins in cellulose utilization and suggests that C. japonicus utilizes a combination of hydrolytic and oxidative cleavage mechanisms to degrade cellulose.

Aromatic inhibitors derived from ammonia-pretreated lignocellulose hinder bacterial ethanologenesis by activating regulatory circuits controlling inhibitor efflux and detoxification.

Efficient microbial conversion of lignocellulosic hydrolysates to biofuels is a key barrier to the economically viable deployment of lignocellulosic biofuels. A chief contributor to this barrier is the impact on microbial processes and energy metabolism of lignocellulose-derived inhibitors, including phenolic carboxylates, phenolic amides (for ammonia-pretreated biomass), phenolic aldehydes, and furfurals. To understand the bacterial pathways induced by inhibitors present in ammonia-pretreated biomass hydrolysates, which are less well studied than acid-pretreated biomass hydrolysates, we developed and exploited synthetic mimics of ammonia-pretreated corn stover hydrolysate (ACSH). To determine regulatory responses to the inhibitors normally present in ACSH, we measured transcript and protein levels in an Escherichia coli ethanologen using RNA-seq and quantitative proteomics during fermentation to ethanol of synthetic hydrolysates containing or lacking the inhibitors. Our study identified four major regulators mediating these responses, the MarA/SoxS/Rob network, AaeR, FrmR, and YqhC. Induction of these regulons was correlated with a reduced rate of ethanol production, buildup of pyruvate, depletion of ATP and NAD(P)H, and an inhibition of xylose conversion. The aromatic aldehyde inhibitor 5-hydroxymethylfurfural appeared to be reduced to its alcohol form by the ethanologen during fermentation, whereas phenolic acid and amide inhibitors were not metabolized. Together, our findings establish that the major regulatory responses to lignocellulose-derived inhibitors are mediated by transcriptional rather than translational regulators, suggest that energy consumed for inhibitor efflux and detoxification may limit biofuel production, and identify a network of regulators for future synthetic biology efforts.

Altered lipid composition and enhanced nutritional value of Arabidopsis leaves following introduction of an algal diacylglycerol acyltransferase 2.

Enhancement of acyl-CoA-dependent triacylglycerol (TAG) synthesis in vegetative tissues is widely discussed as a potential avenue to increase the energy density of crops. Here, we report the identification and characterization of Chlamydomonas reinhardtii diacylglycerol acyltransferase type two (DGTT) enzymes and use DGTT2 to alter acyl carbon partitioning in plant vegetative tissues. This enzyme can accept a broad range of acyl-CoA substrates, allowing us to interrogate different acyl pools in transgenic plants. Expression of DGTT2 in Arabidopsis thaliana increased leaf TAG content, with some molecular species containing very-long-chain fatty acids. The acyl compositions of sphingolipids and surface waxes were altered, and cutin was decreased. The increased carbon partitioning into TAGs in the leaves of DGTT2-expressing lines had little effect on transcripts of the sphingolipid/wax/cutin pathway, suggesting that the supply of acyl groups for the assembly of these lipids is not transcriptionally adjusted. Caterpillars of the generalist herbivore Spodoptera exigua reared on transgenic plants gained more weight. Thus, the nutritional value and/or energy density of the transgenic lines was increased by ectopic expression of DGTT2 and acyl groups were diverted from different pools into TAGs, demonstrating the interconnectivity of acyl metabolism in leaves.

Acidic nuclear phosphoprotein 32kDa (ANP32)B-deficient mouse reveals a hierarchy of ANP32 importance in mammalian development.

The highly conserved ANP32 proteins are proposed to function in a broad array of physiological activities through molecular mechanisms as diverse as phosphatase inhibition, chromatin regulation, caspase activation, and intracellular transport. On the basis of previous analyses of mice bearing targeted mutations of Anp32a or Anp32e, there has been speculation that all ANP32 proteins play redundant roles and are dispensable for normal development. However, more recent work has suggested that ANP32B may in fact have functions that are not shared by other ANP32 family members. Here we report that ANP32B expression is associated with a poor prognosis in human breast cancer, consistent with the increased levels of Anp32b mRNA present in proliferating wild-type (WT) murine embryonic fibroblasts and stimulated WT B and T lymphocytes. Moreover, we show that, contrary to previous assumptions, Anp32b is very important for murine embryogenesis. In a mixed genetic background, ANP32B-deficient mice displayed a partially penetrant perinatal lethality that became fully penetrant in a pure C57BL/6 background. Surviving ANP32B-deficient mice showed reduced viability due to variable defects in various organ systems. Study of compound mutants lacking ANP32A, ANP32B, and/or ANP32E revealed previously hidden roles for ANP32A in mouse development that became apparent only in the complete absence of ANP32B. Our data demonstrate a hierarchy of importance for the mammalian Anp32 genes, with Anp32b being the most critical for normal development.

Global gene expression patterns in Clostridium thermocellum as determined by microarray analysis of chemostat cultures on cellulose or cellobiose.

A microarray study of chemostat growth on insoluble cellulose or soluble cellobiose has provided substantial new information on Clostridium thermocellum gene expression. This is the first comprehensive examination of gene expression in C. thermocellum under defined growth conditions. Expression was detected from 2,846 of 3,189 genes, and regression analysis revealed 348 genes whose changes in expression patterns were growth rate and/or substrate dependent. Successfully modeled genes included those for scaffoldin and cellulosomal enzymes, intracellular metabolic enzymes, transcriptional regulators, sigma factors, signal transducers, transporters, and hypothetical proteins. Unique genes encoding glycolytic pathway and ethanol fermentation enzymes expressed at high levels simultaneously with previously established maximal ethanol production were also identified. Ranking of normalized expression intensities revealed significant changes in transcriptional levels of these genes. The pattern of expression of transcriptional regulators, sigma factors, and signal transducers indicates that response to growth rate is the dominant global mechanism used for control of gene expression in C. thermocellum.

Design and analysis of quantitative differential proteomics investigations using LC-MS technology.

Liquid chromatography-mass spectrometry (LC-MS)-based proteomics is becoming an increasingly important tool in characterizing the abundance of proteins in biological samples of various types and across conditions. Effects of disease or drug treatments on protein abundance are of particular interest for the characterization of biological processes and the identification of biomarkers. Although state-of-the-art instrumentation is available to make high-quality measurements and commercially available software is available to process the data, the complexity of the technology and data presents challenges for bioinformaticians and statisticians. Here, we describe a pipeline for the analysis of quantitative LC-MS data. Key components of this pipeline include experimental design (sample pooling, blocking, and randomization) as well as deconvolution and alignment of mass chromatograms to generate a matrix of molecular abundance profiles. An important challenge in LC-MS-based quantitation is to be able to accurately identify and assign abundance measurements to members of protein families. To address this issue, we implement a novel statistical method for inferring the relative abundance of related members of protein families from tryptic peptide intensities. This pipeline has been used to analyze quantitative LC-MS data from multiple biomarker discovery projects. We illustrate our pipeline here with examples from two of these studies, and show that the pipeline constitutes a complete workable framework for LC-MS-based differential quantitation. Supplementary material is available at http://iec01.mie.utoronto.ca/~thodoros/Bukhman/.

Large-scale mapping of human protein-protein interactions by mass spectrometry.

Mapping protein-protein interactions is an invaluable tool for understanding protein function. Here, we report the first large-scale study of protein-protein interactions in human cells using a mass spectrometry-based approach. The study maps protein interactions for 338 bait proteins that were selected based on known or suspected disease and functional associations. Large-scale immunoprecipitation of Flag-tagged versions of these proteins followed by LC-ESI-MS/MS analysis resulted in the identification of 24,540 potential protein interactions. False positives and redundant hits were filtered out using empirical criteria and a calculated interaction confidence score, producing a data set of 6463 interactions between 2235 distinct proteins. This data set was further cross-validated using previously published and predicted human protein interactions. In-depth mining of the data set shows that it represents a valuable source of novel protein-protein interactions with relevance to human diseases. In addition, via our preliminary analysis, we report many novel protein interactions and pathway associations.

Differential analysis of membrane proteins in mouse fore- and hindbrain using a label-free approach.

The ability to quantitatively compare protein levels across different regions of the brain to identify disease mechanisms remains a fundamental research challenge. It requires both a robust method to efficiently isolate proteins from small amounts of tissue and a differential technique that provides a sensitive and comprehensive analysis of these proteins. Here, we describe a proteomic approach for the quantitative mapping of membrane proteins between mouse fore- and hindbrain regions. The approach focuses primarily on a recently developed method for the fractionation of membranes and on-membrane protein digestion, but incorporates off-line SCX-fractionation of the peptide mixture and nano-LC-MS/MS analysis using an LTQ-FT-ICR instrument as part of the analytical method. Comparison of mass spectral peak intensities between samples, mapping of peaks to peptides and protein sequences, and statistical analysis were performed using in-house differential analysis software (DAS). In total, 1213 proteins were identified and 967 were quantified; 81% of the identified proteins were known membrane proteins and 38% of the protein sequences were predicted to contain transmembrane helices. Although this paper focuses primarily on characterizing the efficiency of this purification method from a typical sample set, for many of the quantified proteins such as glutamate receptors, GABA receptors, calcium channel subunits, and ATPases, the observed ratios of protein abundance were in good agreement with the known mRNA expression levels and/or intensities of immunostaining in rostral and caudal regions of murine brain. This suggests that the approach would be well-suited for incorporation in more rigorous, larger scale quantitative analysis designed to achieve biological significance.

Comparison of microarray-based mRNA profiling technologies for identification of psychiatric disease and drug signatures.

The gene expression profiles of human postmortem parietal and prefrontal cortex samples of normal controls and patients with bipolar disease, or human neuroblastoma flat (NBFL) cells treated with the mood-stabilizing drug, valproate, were used to compare the performance of Affymetrix oligonucleotide U133A GeneChips and Agilent Human 1 cDNA microarrays. Among those genes represented on both platforms, the oligo array identified 26-53% more differentially expressed genes compared to the cDNA array in the three experiments, when identical fold change and t-test criteria were applied. The increased sensitivity was primarily the result of more robust fold changes measured by the oligonucleotide system. Essentially all gene changes overlapping between the two platforms were co-directional, and ranged from 4 to 19% depending upon the amount of biological variability within and between the comparison groups. Q-PCR validation rates were virtually identical for the two platforms, with 23-24% validation in the prefrontal cortex experiment, and 56% for both platforms in the cell culture experiment. Validated genes included dopa decarboxylase, dopamine beta-hydroxylase, and dihydropyrimidinase-related protein 3, which were decreased in NBFL cells exposed to valproate, and spinocerebellar ataxia 7, which was increased in bipolar disease. The modest overlap but similar validation rates show that each microarray system identifies a unique set of differentially expressed genes, and thus the greatest information is obtained from the use of both platforms.

Electroconvulsive seizures regulate gene expression of distinct neurotrophic signaling pathways.

Electroconvulsive therapy (ECT) remains the treatment of choice for drug-resistant patients with depressive disorders, yet the mechanism for its efficacy remains unknown. Gene transcription changes were measured in the frontal cortex and hippocampus of rats subjected to sham seizures or to 1 or 10 electroconvulsive seizures (ECS), a model of ECT. Among the 3500-4400 RNA sequences detected in each sample, ECS increased by 1.5- to 11-fold or decreased by at least 34% the expression of 120 unique genes. The hippocampus produced more than three times the number of gene changes seen in the cortex, and many hippocampal gene changes persisted with chronic ECS, unlike in the cortex. Among the 120 genes, 77 have not been reported in previous studies of ECS or seizure responses, and 39 were confirmed among 59 studied by quantitative real time PCR. Another 19 genes, 10 previously unreported, changed by <1.5-fold but with very high significance. Multiple genes were identified within distinct pathways, including the BDNF-MAP kinase-cAMP-cAMP response element-binding protein pathway (15 genes), the arachidonic acid pathway (5 genes), and more than 10 genes in each of the immediate-early gene, neurogenesis, and exercise response gene groups. Neurogenesis, neurite outgrowth, and neuronal plasticity associated with BDNF, glutamate, and cAMP-protein kinase A signaling pathways may mediate the antidepressant effects of ECT in humans. These genes, and others that increase only with chronic ECS such as neuropeptide Y and thyrotropin-releasing hormone, may provide novel ways to select drugs for the treatment of depression and mimic the rapid effectiveness of ECT.