Liquid Glycan Array.

The M13 phage platform is a stable and monodisperse nanoscale carrier, which can be modified with different molecules by chemical conjugation strategies. Here, we describe M13 phage acylated on pVIII protein with a dibenzocyclooctyne reacting with azido glycan to yield 30-1500 copy numbers of glycan per phage and monitored by MALDI-TOF spectrometry to generate multivalent glycoconjugates that contain desired densities of glycans. We prepared the liquid glycan arrays (LiGA) such that both the structure and density of glycans were encoded in the DNA of the bacteriophage. The LiGA can be used to validate the binding properties of glycans to purified lectins and explore the effect of glycan density on such binding. From a mixture of multivalent glycan probes, LiGAs can also identify the glycoconjugates with optimal avidity necessary for binding to lectins on living cells in vitro and live animals in vivo.

Chemoenzymatic synthesis of genetically-encoded multivalent liquid N-glycan arrays.

Cellular glycosylation is characterized by chemical complexity and heterogeneity, which is challenging to reproduce synthetically. Here we show chemoenzymatic synthesis on phage to produce a genetically-encoded liquid glycan array (LiGA) of complex type N-glycans. Implementing the approach involved by ligating an azide-containing sialylglycosyl-asparagine to phage functionalized with 50-1000 copies of dibenzocyclooctyne. The resulting intermediate can be trimmed by glycosidases and extended by glycosyltransferases yielding a phage library with different N-glycans. Post-reaction analysis by MALDI-TOF MS allows rigorous characterization of N-glycan structure and mean density, which are both encoded in the phage DNA. Use of this LiGA with fifteen glycan-binding proteins, including CD22 or DC-SIGN on cells, reveals optimal structure/density combinations for recognition. Injection of the LiGA into mice identifies glycoconjugates with structures and avidity necessary for enrichment in specific organs. This work provides a quantitative evaluation of the interaction of complex N-glycans with GBPs in vitro and in vivo.

GlyNet: a multi-task neural network for predicting protein-glycan interactions.

Advances in diagnostics, therapeutics, vaccines, transfusion, and organ transplantation build on a fundamental understanding of glycan-protein interactions. To aid this, we developed GlyNet, a model that accurately predicts interactions (relative binding strengths) between mammalian glycans and 352 glycan-binding proteins, many at multiple concentrations. For each glycan input, our model produces 1257 outputs, each representing the relative interaction strength between the input glycan and a particular protein sample. GlyNet learns these continuous values using relative fluorescence units (RFUs) measured on 599 glycans in the Consortium for Functional Glycomics glycan arrays and extrapolates these to RFUs from additional, untested glycans. GlyNet's output of continuous values provides more detailed results than the standard binary classification models. After incorporating a simple threshold to transform such continuous outputs the resulting GlyNet classifier outperforms those standard classifiers. GlyNet is the first multi-output regression model for predicting protein-glycan interactions and serves as an important benchmark, facilitating development of quantitative computational glycobiology.

Isolation of a Human Betaretrovirus from Patients with Primary Biliary Cholangitis.

A human betaretrovirus (HBRV) has been linked with the autoimmune liver disease, primary biliary cholangitis (PBC), and various cancers, including breast cancer and lymphoma. HBRV is closely related to the mouse mammary tumor virus, and represents the only exogenous betaretrovirus characterized in humans to date. Evidence of infection in patients with PBC has been demonstrated through the identification of proviral integration sites in lymphoid tissue, the major reservoir of infection, as well as biliary epithelium, which is the site of the disease process. Accordingly, we tested the hypothesis that patients with PBC harbor a transmissible betaretrovirus by co-cultivation of PBC patients' lymph node homogenates with the HS578T breast cancer line. Because of the low level of HBRV replication, betaretrovirus producing cells were subcloned to optimize viral isolation and production. Evidence of infection was provided by electron microscopy, RT-PCR, in situ hybridization, cloning of the HBRV proviral genome and demonstration of more than 3400 integration sites. Further evidence of viral transmissibility was demonstrated by infection of biliary epithelial cells. While HBRV did not show a preference for integration proximal to specific genomic features, analyses of common insertion sites revealed evidence of integration proximal to cancer associated genes. These studies demonstrate the isolation of HBRV with features similar to mouse mammary tumor virus and confirm that patients with PBC display evidence of a transmissible viral infection.

Metagenomics Versus Metatranscriptomics of the Murine Gut Microbiome for Assessing Microbial Metabolism During Inflammation.

Shotgun metagenomics studies have improved our understanding of microbial population dynamics and have revealed significant contributions of microbes to gut homeostasis. They also allow in silico inference of the metagenome. While they link the microbiome with metabolic abnormalities associated with disease phenotypes, they do not capture microbial gene expression patterns that occur in response to the multitude of stimuli that constantly ambush the gut environment. Metatranscriptomics closes that gap, but its implementation is more expensive and tedious. We assessed the metabolic perturbations associated with gut inflammation using shotgun metagenomics and metatranscriptomics. Shotgun metagenomics detected changes in abundance of bacterial taxa known to be SCFA producers, which favors gut homeostasis. Bacteria in the phylum Firmicutes were found at decreased abundance, while those in phyla Bacteroidetes and Proteobacteria were found at increased abundance. Surprisingly, inferring the coding capacity of the microbiome from shotgun metagenomics data did not result in any statistically significant difference, suggesting functional redundancy in the microbiome or poor resolution of shotgun metagenomics data to profile bacterial pathways, especially when sequencing is not very deep. Obviously, the ability of metatranscriptomics libraries to detect transcripts expressed at basal (or simply low) levels is also dependent on sequencing depth. Nevertheless, metatranscriptomics informed about contrasting roles of bacteria during inflammation. Functions involved in nutrient transport, immune suppression and regulation of tissue damage were dramatically upregulated, perhaps contributed by homeostasis-promoting bacteria. Functions ostensibly increasing bacteria pathogenesis were also found upregulated, perhaps as a consequence of increased abundance of Proteobacteria. Bacterial protein synthesis appeared downregulated. In summary, shotgun metagenomics was useful to profile bacterial population composition and taxa relative abundance, but did not inform about differential gene content associated with inflammation. Metatranscriptomics was more robust for capturing bacterial metabolism in real time. Although both approaches are complementary, it is often not possible to apply them in parallel. We hope our data will help researchers to decide which approach is more appropriate for the study of different aspects of the microbiome.

Genetically encoded multivalent liquid glycan array displayed on M13 bacteriophage.

The central dogma of biology does not allow for the study of glycans using DNA sequencing. We report a liquid glycan array (LiGA) platform comprising a library of DNA 'barcoded' M13 virions that display 30-1,500 copies of glycans per phage. A LiGA is synthesized by acylation of the phage pVIII protein with a dibenzocyclooctyne, followed by ligation of azido-modified glycans. Pulldown of the LiGA with lectins followed by deep sequencing of the barcodes in the bound phage decodes the optimal structure and density of the recognized glycans. The LiGA is target agnostic and can measure the glycan-binding profile of lectins, such as CD22, on cells in vitro and immune cells in a live mouse. From a mixture of multivalent glycan probes, LiGAs identify the glycoconjugates with optimal avidity necessary for binding to lectins on living cells in vitro and in vivo.

Access to RNA-sequencing data from 1,173 plant species: The 1000 Plant transcriptomes initiative (1KP).

BACKGROUND: The 1000 Plant transcriptomes initiative (1KP) explored genetic diversity by sequencing RNA from 1,342 samples representing 1,173 species of green plants (Viridiplantae). FINDINGS: This data release accompanies the initiative's final/capstone publication on a set of 3 analyses inferring species trees, whole genome duplications, and gene family expansions. These and previous analyses are based on de novo transcriptome assemblies and related gene predictions. Here, we assess their data and assembly qualities and explain how we detected potential contaminations. CONCLUSIONS: These data will be useful to plant and/or evolutionary scientists with interests in particular gene families, either across the green plant tree of life or in more focused lineages.

Impact of sequencing depth and technology on de novo RNA-Seq assembly.

BACKGROUND: RNA-Seq data is inherently nonuniform for different transcripts because of differences in gene expression. This makes it challenging to decide how much data should be generated from each sample. How much should one spend to recover the less expressed transcripts? The sequencing technology used is another consideration, as there are inevitably always biases against certain sequences. To investigate these effects, we first looked at high-depth libraries from a set of well-annotated organisms to ascertain the impact of sequencing depth on de novo assembly. We then looked at libraries sequenced from the Universal Human Reference RNA (UHRR) to compare the performance of Illumina HiSeq and MGI DNBseq™ technologies. RESULTS: On the issue of sequencing depth, the amount of exomic sequence assembled plateaued using data sets of approximately 2 to 8 Gbp. However, the amount of genomic sequence assembled did not plateau for many of the analyzed organisms. Most of the unannotated genomic sequences are single-exon transcripts whose biological significance will be questionable for some users. On the issue of sequencing technology, both of the analyzed platforms recovered a similar number of full-length transcripts. The missing "gap" regions in the HiSeq assemblies were often attributed to higher GC contents, but this may be an artefact of library preparation and not of sequencing technology. CONCLUSIONS: Increasing sequencing depth beyond modest data sets of less than 10 Gbp recovers a plethora of single-exon transcripts undocumented in genome annotations. DNBseq™ is a viable alternative to HiSeq for de novo RNA-Seq assembly.

Insights into the Evolution of Hydroxyproline-Rich Glycoproteins from 1000 Plant Transcriptomes.

The carbohydrate-rich cell walls of land plants and algae have been the focus of much interest given the value of cell wall-based products to our current and future economies. Hydroxyproline-rich glycoproteins (HRGPs), a major group of wall glycoproteins, play important roles in plant growth and development, yet little is known about how they have evolved in parallel with the polysaccharide components of walls. We investigate the origins and evolution of the HRGP superfamily, which is commonly divided into three major multigene families: the arabinogalactan proteins (AGPs), extensins (EXTs), and proline-rich proteins. Using motif and amino acid bias, a newly developed bioinformatics pipeline, we identified HRGPs in sequences from the 1000 Plants transcriptome project (www.onekp.com). Our analyses provide new insights into the evolution of HRGPs across major evolutionary milestones, including the transition to land and the early radiation of angiosperms. Significantly, data mining reveals the origin of glycosylphosphatidylinositol (GPI)-anchored AGPs in green algae and a 3- to 4-fold increase in GPI-AGPs in liverworts and mosses. The first detection of cross-linking (CL)-EXTs is observed in bryophytes, which suggests that CL-EXTs arose though the juxtaposition of preexisting SPn EXT glycomotifs with refined Y-based motifs. We also detected the loss of CL-EXT in a few lineages, including the grass family (Poaceae), that have a cell wall composition distinct from other monocots and eudicots. A key challenge in HRGP research is tracking individual HRGPs throughout evolution. Using the 1000 Plants output, we were able to find putative orthologs of Arabidopsis pollen-specific GPI-AGPs in basal eudicots.

Occurrence, structure, and evolution of nitric oxide synthase-like proteins in the plant kingdom.

Nitric oxide (NO) signaling regulates various physiological processes in both animals and plants. In animals, NO synthesis is mainly catalyzed by NO synthase (NOS) enzymes. Although NOS-like activities that are sensitive to mammalian NOS inhibitors have been detected in plant extracts, few bona fide plant NOS enzymes have been identified. We searched the data set produced by the 1000 Plants (1KP) international consortium for the presence of transcripts encoding NOS-like proteins in over 1000 species of land plants and algae. We also searched for genes encoding NOS-like enzymes in 24 publicly available algal genomes. We identified no typical NOS sequences in 1087 sequenced transcriptomes of land plants. In contrast, we identified NOS-like sequences in 15 of the 265 algal species analyzed. Even if the presence of NOS enzymes assembled from multipolypeptides in plants cannot be conclusively discarded, the emerging data suggest that, instead of generating NO with evolutionarily conserved NOS enzymes, land plants have evolved finely regulated nitrate assimilation and reduction processes to synthesize NO through a mechanism different than that in animals.

Functional and topological diversity of LOV domain photoreceptors.

Light-oxygen-voltage sensitive (LOV) flavoproteins are ubiquitous photoreceptors that mediate responses to environmental cues. Photosensory inputs are transduced into signaling outputs via structural rearrangements in sensor domains that consequently modulate the activity of an effector domain or multidomain clusters. Establishing the diversity in effector function and sensor-effector topology will inform what signaling mechanisms govern light-responsive behaviors across multiple kingdoms of life and how these signals are transduced. Here, we report the bioinformatics identification of over 6,700 candidate LOV domains (including over 4,000 previously unidentified sequences from plants and protists), and insights from their annotations for ontological function and structural arrangements. Motif analysis identified the sensors from ∼42 million ORFs, with strong statistical separation from other flavoproteins and non-LOV members of the structurally related Per-aryl hydrocarbon receptor nuclear translocator (ARNT)-Sim family. Conserved-domain analysis determined putative light-regulated function and multidomain topologies. We found that for certain effectors, sensor-effector linker length is discretized based on both phylogeny and the preservation of α-helical heptad repeats within an extended coiled-coil linker structure. This finding suggests that preserving sensor-effector orientation is a key determinant of linker length, in addition to ancestry, in LOV signaling structure-function. We found a surprisingly high prevalence of effectors with functions previously thought to be rare among LOV proteins, such as regulators of G protein signaling, and discovered several previously unidentified effectors, such as lipases. This work highlights the value of applying genomic and transcriptomic technologies to diverse organisms to capture the structural and functional variation in photosensory proteins that are vastly important in adaptation, photobiology, and optogenetics.

Peptidomics of Circular Cysteine-Rich Plant Peptides: Analysis of the Diversity of Cyclotides from Viola tricolor by Transcriptome and Proteome Mining.

Cyclotides are plant-derived mini proteins. They are genetically encoded as precursor proteins that become post-translationally modified to yield circular cystine-knotted molecules. Because of this structural topology cyclotides resist enzymatic degradation in biological fluids, and hence they are considered as promising lead molecules for pharmaceutical applications. Despite ongoing efforts to discover novel cyclotides and analyze their biodiversity, it is not clear how many individual peptides a single plant specimen can express. Therefore, we investigated the transcriptome and cyclotide peptidome of Viola tricolor. Transcriptome mining enabled the characterization of cyclotide precursor architecture and processing sites important for biosynthesis of mature peptides. The cyclotide peptidome was explored by mass spectrometry and bottom-up proteomics using the extracted peptide sequences as queries for database searching. In total 164 cyclotides were discovered by nucleic acid and peptide analysis in V. tricolor. Therefore, violaceous plants at a global scale may be the source to as many as 150 000 individual cyclotides. Encompassing the diversity of V. tricolor as a combinatorial library of bioactive peptides, this commercially available medicinal herb may be a suitable starting point for future bioactivity-guided screening studies.

The Origin and Evolution of the Plant Cell Surface: Algal Integrin-Associated Proteins and a New Family of Integrin-Like Cytoskeleton-ECM Linker Proteins.

The extracellular matrix of scaly green flagellates consists of small organic scales consisting of polysaccharides and scale-associated proteins (SAPs). Molecular phylogenies have shown that these organisms represent the ancestral stock of flagellates from which all green plants (Viridiplantae) evolved. The molecular characterization of four different SAPs is presented. Three SAPs are type-2 membrane proteins with an arginine/alanine-rich short cytoplasmic tail and an extracellular domain that is most likely of bacterial origin. The fourth protein is a filamin-like protein. In addition, we report the presence of proteins similar to the integrin-associated proteins α-actinin (in transcriptomes of glaucophytes and some viridiplants), LIM-domain proteins, and integrin-associated kinase in transcriptomes of viridiplants, glaucophytes, and rhodophytes. We propose that the membrane proteins identified are the predicted linkers between scales and the cytoskeleton. These proteins are present in many green algae but are apparently absent from embryophytes. These proteins represent a new protein family we have termed gralins for green algal integrins. Gralins are absent from embryophytes. A model for the evolution of the cell surface proteins in Plantae is discussed.

Dissecting Molecular Evolution in the Highly Diverse Plant Clade Caryophyllales Using Transcriptome Sequencing.

Many phylogenomic studies based on transcriptomes have been limited to "single-copy" genes due to methodological challenges in homology and orthology inferences. Only a relatively small number of studies have explored analyses beyond reconstructing species relationships. We sampled 69 transcriptomes in the hyperdiverse plant clade Caryophyllales and 27 outgroups from annotated genomes across eudicots. Using a combined similarity- and phylogenetic tree-based approach, we recovered 10,960 homolog groups, where each was represented by at least eight ingroup taxa. By decomposing these homolog trees, and taking gene duplications into account, we obtained 17,273 ortholog groups, where each was represented by at least ten ingroup taxa. We reconstructed the species phylogeny using a 1,122-gene data set with a gene occupancy of 92.1%. From the homolog trees, we found that both synonymous and nonsynonymous substitution rates in herbaceous lineages are up to three times as fast as in their woody relatives. This is the first time such a pattern has been shown across thousands of nuclear genes with dense taxon sampling. We also pinpointed regions of the Caryophyllales tree that were characterized by relatively high frequencies of gene duplication, including three previously unrecognized whole-genome duplications. By further combining information from homolog tree topology and synonymous distance between paralog pairs, phylogenetic locations for 13 putative genome duplication events were identified. Genes that experienced the greatest gene family expansion were concentrated among those involved in signal transduction and oxidoreduction, including a cytochrome P450 gene that encodes a key enzyme in the betalain synthesis pathway. Our approach demonstrates a new approach for functional phylogenomic analysis in nonmodel species that is based on homolog groups in addition to inferred ortholog groups.

Independent optical excitation of distinct neural populations.

Optogenetic tools enable examination of how specific cell types contribute to brain circuit functions. A long-standing question is whether it is possible to independently activate two distinct neural populations in mammalian brain tissue. Such a capability would enable the study of how different synapses or pathways interact to encode information in the brain. Here we describe two channelrhodopsins, Chronos and Chrimson, discovered through sequencing and physiological characterization of opsins from over 100 species of alga. Chrimson's excitation spectrum is red shifted by 45 nm relative to previous channelrhodopsins and can enable experiments in which red light is preferred. We show minimal visual system-mediated behavioral interference when using Chrimson in neurobehavioral studies in Drosophila melanogaster. Chronos has faster kinetics than previous channelrhodopsins yet is effectively more light sensitive. Together these two reagents enable two-color activation of neural spiking and downstream synaptic transmission in independent neural populations without detectable cross-talk in mouse brain slice.

Data access for the 1,000 Plants (1KP) project.

The 1,000 plants (1KP) project is an international multi-disciplinary consortium that has generated transcriptome data from over 1,000 plant species, with exemplars for all of the major lineages across the Viridiplantae (green plants) clade. Here, we describe how to access the data used in a phylogenomics analysis of the first 85 species, and how to visualize our gene and species trees. Users can develop computational pipelines to analyse these data, in conjunction with data of their own that they can upload. Computationally estimated protein-protein interactions and biochemical pathways can be visualized at another site. Finally, we comment on our future plans and how they fit within this scalable system for the dissemination, visualization, and analysis of large multi-species data sets.

Cyclotide discovery in Gentianales revisited--identification and characterization of cyclic cystine-knot peptides and their phylogenetic distribution in Rubiaceae plants.

Cyclotides are a unique class of ribosomally synthesized cysteine-rich miniproteins characterized by a head-to-tail cyclized backbone and three conserved disulfide-bonds in a knotted arrangement. Originally they were discovered in the coffee-family plant Oldenlandia affinis (Rubiaceae) and have since been identified in several species of the violet, cucurbit, pea, potato, and grass families. However, the identification of novel cyclotide-containing plant species still is a major challenge due to the lack of a rapid and accurate analytical workflow in particular for large sampling numbers. As a consequence, their phylogeny in the plant kingdom remains unclear. To gain further insight into the distribution and evolution of plant cyclotides, we analyzed ∼300 species of >40 different families, with special emphasis on plants from the order Gentianales. For this purpose, we have developed a refined screening methodology combining chemical analysis of plant extracts and bioinformatic analysis of transcript databases. Using mass spectrometry and transcriptome-mining, we identified nine novel cyclotide-containing species and their related cyclotide precursor genes in the tribe Palicoureeae. The characterization of novel peptide sequences underlines the high variability and plasticity of the cyclotide framework, and a comparison of novel precursor proteins from Carapichea ipecacuanha illustrated their typical cyclotide gene architectures. Phylogenetic analysis of their distribution within the Psychotria alliance revealed cyclotides to be restricted to Palicourea, Margaritopsis, Notopleura, Carapichea, Chassalia, and Geophila. In line with previous reports, our findings confirm cyclotides to be one of the largest peptide families within the plant kingdom and suggest that their total number may exceed tens of thousands.

Evaluating methods for isolating total RNA and predicting the success of sequencing phylogenetically diverse plant transcriptomes.

Next-generation sequencing plays a central role in the characterization and quantification of transcriptomes. Although numerous metrics are purported to quantify the quality of RNA, there have been no large-scale empirical evaluations of the major determinants of sequencing success. We used a combination of existing and newly developed methods to isolate total RNA from 1115 samples from 695 plant species in 324 families, which represents >900 million years of phylogenetic diversity from green algae through flowering plants, including many plants of economic importance. We then sequenced 629 of these samples on Illumina GAIIx and HiSeq platforms and performed a large comparative analysis to identify predictors of RNA quality and the diversity of putative genes (scaffolds) expressed within samples. Tissue types (e.g., leaf vs. flower) varied in RNA quality, sequencing depth and the number of scaffolds. Tissue age also influenced RNA quality but not the number of scaffolds ≥ 1000 bp. Overall, 36% of the variation in the number of scaffolds was explained by metrics of RNA integrity (RIN score), RNA purity (OD 260/230), sequencing platform (GAIIx vs HiSeq) and the amount of total RNA used for sequencing. However, our results show that the most commonly used measures of RNA quality (e.g., RIN) are weak predictors of the number of scaffolds because Illumina sequencing is robust to variation in RNA quality. These results provide novel insight into the methods that are most important in isolating high quality RNA for sequencing and assembling plant transcriptomes. The methods and recommendations provided here could increase the efficiency and decrease the cost of RNA sequencing for individual labs and genome centers.

A genome triplication associated with early diversification of the core eudicots.

BACKGROUND: Although it is agreed that a major polyploidy event, gamma, occurred within the eudicots, the phylogenetic placement of the event remains unclear. RESULTS: To determine when this polyploidization occurred relative to speciation events in angiosperm history, we employed a phylogenomic approach to investigate the timing of gene set duplications located on syntenic gamma blocks. We populated 769 putative gene families with large sets of homologs obtained from public transcriptomes of basal angiosperms, magnoliids, asterids, and more than 91.8 gigabases of new next-generation transcriptome sequences of non-grass monocots and basal eudicots. The overwhelming majority (95%) of well-resolved gamma duplications was placed before the separation of rosids and asterids and after the split of monocots and eudicots, providing strong evidence that the gamma polyploidy event occurred early in eudicot evolution. Further, the majority of gene duplications was placed after the divergence of the Ranunculales and core eudicots, indicating that the gamma appears to be restricted to core eudicots. Molecular dating estimates indicate that the duplication events were intensely concentrated around 117 million years ago. CONCLUSIONS: The rapid radiation of core eudicot lineages that gave rise to nearly 75% of angiosperm species appears to have occurred coincidentally or shortly following the gamma triplication event. Reconciliation of gene trees with a species phylogeny can elucidate the timing of major events in genome evolution, even when genome sequences are only available for a subset of species represented in the gene trees. Comprehensive transcriptome datasets are valuable complements to genome sequences for high-resolution phylogenomic analysis.