The InterPro protein families database: the classification resource after 15 years

Drought responsive gene expression regulatory divergence between upland and lowland ecotypes of a perennial C4 grass

Climatic adaptation is an example of a genotype-by-environment interaction (G×E) of fitness. Selection upon gene expression regulatory variation can contribute to adaptive phenotypic diversity; however, surprisingly few studies have examined how genome-wide patterns of gene expression G×E are manifested in response to environmental stress and other selective agents that cause climatic adaptation. Here, we characterize drought-responsive expression divergence between upland (drought-adapted) and lowland (mesic) ecotypes of the perennial C4 grass,Panicum hallii, in natural field conditions. Overall, we find that cis-regulatory elements contributed to gene expression divergence across 47% of genes, 7.2% of which exhibit drought-responsive G×E. While less well-represented, we observe 1294 genes (7.8%) with transeffects.Trans-by-environment interactions are weaker and much less common than cis G×E, occurring in only 0.7% oft rans-regulated genes. Finally, gene expression heterosis is highly enriched in expression phenotypes with significant G×E. As such, modes of inheritance that drive heterosis, such as dominance or overdominance, may be common among G×E genes. Interestingly, motifs specific to drought-responsive transcription factors are highly enriched in the promoters of genes exhibiting G×E and transregulation, indicating that expression G×E and heterosis may result from the evolution of transcription factors or their binding sites.P. hallii serves as the genomic model for its close relative and emerging biofuel crop, switchgrass (Panicum virgatum). Accordingly, the results here not only aid in the discovery of the genetic mechanisms that underlie local adaptation but also provide a foundation to improve switchgrass yield under water-limited conditions.

Biochemical and Spectroscopic Characterization of a Radical S-Adenosyl-L-methionine Enzyme Involved in the Formation of a Peptide Thioether Cross-Link

Peptide-derived natural products are a class of metabolites that afford the producing organism a selective advantage over other organisms in their biological niche. While the polypeptide antibiotics produced by the nonribosomal polypeptide synthetases (NRPS) are the most widely recognized, the ribosomally synthesized and post-translationally modified peptides (RiPPs) are an emerging group of natural products with diverse structures and biological functions. Both the NRPS derived peptides and the RiPPs undergo extensive post-translational modifications to produce structural diversity. Here we report the first characterization of the six cysteines in forty-five (SCIFF) [Haft, D. H. and Basu M. K. (2011) J. Bacteriol. 193, 2745-2755] peptide maturase Tte1186, which is a member of the radical S-adenosyl-l-methionine (SAM) superfamily. Tte1186 catalyzes the formation of a thioether cross-link in the peptide Tte1186a encoded by an orf located upstream of the maturase, under reducing conditions in the presence of SAM. Tte1186 contains three [4Fe-4S] clusters that are indispensable for thioether cross-link formation; however, only one cluster catalyzes the reductive cleavage of SAM. Mechanistic imperatives for the reaction catalyzed by the thioether forming radical SAM maturases will be discussed.

Chromosomal-Level Assembly of the Asian Seabass Genome Using Long Sequence Reads and Multi-layered Scaffolding

We report here the ~670 Mb genome assembly of the Asian seabass (Lates calcarifer), a tropical marine teleost. We used long-read sequencing augmented by transcriptomics, optical and genetic mapping along with shared synteny from closely related fish species to derive a chromosome-level assembly with a contig N50 size over 1 Mb and scaffold N50 size over 25 Mb that span ~90% of the genome. The population structure of L. calcarifer species complex was analyzed by re-sequencing 61 individuals representing various regions across the species' native range. SNP analyses identified high levels of genetic diversity and confirmed earlier indications of a population stratification comprising three clades with signs of admixture apparent in the South-East Asian population. The quality of the Asian seabass genome assembly far exceeds that of any other fish species, and will serve as a new standard for fish genomics.

Caenorhabditis elegans PAQR-2 and IGLR-2 Protect against Glucose Toxicity by Modulating Membrane Lipid Composition

In spite of the worldwide impact of diabetes on human health, the mechanisms behind glucose toxicity remain elusive. Here we show that C. elegans mutants lacking paqr-2, the worm homolog of the adiponectin receptors AdipoR1/2, or its newly identified functional partner iglr-2, are glucose intolerant and die in the presence of as little as 20 mM glucose. Using FRAP (Fluorescence Recovery After Photobleaching) on living worms, we found that cultivation in the presence of glucose causes a decrease in membrane fluidity in paqr-2 and iglr-2 mutants and that genetic suppressors of this sensitivity act to restore membrane fluidity by promoting fatty acid desaturation. The essential roles of paqr-2 and iglr-2 in the presence of glucose are completely independent from daf-2 and daf-16, the C. elegans homologs of the insulin receptor and its downstream target FoxO, respectively. Using bimolecular fluorescence complementation, we also show that PAQR-2 and IGLR-2 interact on plasma membranes and thus may act together as a fluidity sensor that controls membrane lipid composition.

Topology based identification and comprehensive classification of four-transmembrane helix containing proteins (4TMs) in the human genome

BACKGROUND

Membrane proteins are key components in a large spectrum of diverse functions and thus account for the major proportion of the drug-targeted portion of the genome. From a structural perspective, the α-helical transmembrane proteins can be categorized into major groups based on the number of transmembrane helices and these groups are often associated with specific functions. When compared to the well-characterized seven-transmembrane containing proteins (7TM), other TM groups are less explored and in particular the 4TM group. In this study, we identify the complete 4TM complement from the latest release of the human genome and assess the 4TM structure group as a whole. We functionally characterize this dataset and evaluate the resulting groups and ubiquitous functions, and furthermore describe disease and drug target involvement.

RESULTS

We classified 373 proteins, which represents ~7 % of the human membrane proteome, and includes 69 more proteins than our previous estimate. We have characterized the 4TM dataset based on functional, structural, and/or evolutionary similarities. Proteins that are involved in transport activity constitute 37 % of the dataset, 23 % are receptor-related, and 13 % have enzymatic functions. Intriguingly, proteins involved in transport are more than double the 15 % of transporters in the entire human membrane proteome, which might suggest that the 4TM topological architecture is more favored for transporting molecules over other functions. Moreover, we found an interesting exception to the ubiquitous intracellular N- and C-termini localization that is found throughout the entire membrane proteome and 4TM dataset in the neurotransmitter gated ion channel families. Overall, we estimate that 58 % of the dataset has a known association to disease conditions with 19 % of the genes possibly involved in different types of cancer.

CONCLUSIONS

We provide here the most robust and updated classification of the 4TM complement of the human genome as a platform to further understand the characteristics of 4TM functions and to explore pharmacological opportunities.

Ligand-binding specificity and promiscuity of the main lignocellulolytic enzyme families as revealed by active-site architecture analysis

Biomass can be converted into sugars by a series of lignocellulolytic enzymes, which belong to the glycoside hydrolase (GH) families summarized in CAZy databases. Here, using a structural bioinformatics method, we analyzed the active site architecture of the main lignocellulolytic enzyme families. The aromatic amino acids Trp/Tyr and polar amino acids Glu/Asp/Asn/Gln/Arg occurred at higher frequencies in the active site architecture than in the whole enzyme structure. And the number of potential subsites was significantly different among different families. In the cellulase and xylanase families, the conserved amino acids in the active site architecture were mostly found at the −2 to +1 subsites, while in β-glucosidase they were mainly concentrated at the −1 subsite. Families with more conserved binding amino acid residues displayed strong selectivity for their ligands, while those with fewer conserved binding amino acid residues often exhibited promiscuity when recognizing ligands. Enzymes with different activities also tended to bind different hydroxyl oxygen atoms on the ligand. These results may help us to better understand the common and unique structural bases of enzyme-ligand recognition from different families and provide a theoretical basis for the functional evolution and rational design of major lignocellulolytic enzymes.

Genomic and Transcriptomic Features of Response to Anti-PD-1 Therapy in Metastatic Melanoma

PD-1 immune checkpoint blockade provides significant clinical benefits for melanoma patients. We analyzed the somatic mutanomes and transcriptomes of pretreatment melanoma biopsies to identify factors that may influence innate sensitivity or resistance to anti-PD-1 therapy. We find that overall high mutational loads associate with improved survival, and tumors from responding patients are enriched for mutations in the DNA repair gene BRCA2. Innately resistant tumors display a transcriptional signature (referred to as the IPRES, or innate anti-PD-1 resistance), indicating concurrent up-expression of genes involved in the regulation of mesenchymal transition, cell adhesion, extracellular matrix remodeling, angiogenesis, and wound healing. Notably, mitogen-activated protein kinase (MAPK)-targeted therapy (MAPK inhibitor) induces similar signatures in melanoma, suggesting that a non-genomic form of MAPK inhibitor resistance mediates cross-resistance to anti-PD-1 therapy. Validation of the IPRES in other independent tumor cohorts defines a transcriptomic subset across distinct types of advanced cancer. These findings suggest that attenuating the biological processes that underlie IPRES may improve anti-PD-1 response in melanoma and other cancer types.

Global phylogeography and evolutionary history of Shigella dysenteriae type 1

Together with plague, smallpox and typhus, epidemics of dysentery have been a major scourge of human populations for centuries(1). A previous genomic study concluded that Shigella dysenteriae type 1 (Sd1), the epidemic dysentery bacillus, emerged and spread worldwide after the First World War, with no clear pattern of transmission(2). This is not consistent with the massive cyclic dysentery epidemics reported in Europe during the eighteenth and nineteenth centuries(1,3,4) and the first isolation of Sd1 in Japan in 1897(5). Here, we report a whole-genome analysis of 331 Sd1 isolates from around the world, collected between 1915 and 2011, providing us with unprecedented insight into the historical spread of this pathogen. We show here that Sd1 has existed since at least the eighteenth century and that it swept the globe at the end of the nineteenth century, diversifying into distinct lineages associated with the First World War, Second World War and various conflicts or natural disasters across Africa, Asia and Central America. We also provide a unique historical perspective on the evolution of antibiotic resistance over a 100-year period, beginning decades before the antibiotic era, and identify a prevalent multiple antibiotic-resistant lineage in South Asia that was transmitted in several waves to Africa, where it caused severe outbreaks of disease.

GO annotation in InterPro: why stability does not indicate accuracy in a sea of changing annotations

The removal of annotation from biological databases is often perceived as an indicator of erroneous annotation. As a corollary, annotation stability is considered to be a measure of reliability. However, diverse data-driven events can affect the stability of annotations in both primary protein sequence databases and the protein family databases that are built upon the sequence databases and used to help annotate them. Here, we describe some of these events and their consequences for the InterPro database, and demonstrate that annotation removal or reassignment is not always linked to incorrect annotation by the curator.

Database URL:http://www.ebi.ac.uk/interpro

An integrative data analysis platform for gene set analysis and knowledge discovery in a data warehouse framework

Data analysis is one of the most critical and challenging steps in drug discovery and disease biology. A user-friendly resource to visualize and analyse high-throughput data provides a powerful medium for both experimental and computational biologists to understand vastly different biological data types and obtain a concise, simplified and meaningful output for better knowledge discovery. We have previously developed TargetMine, an integrated data warehouse optimized for target prioritization. Here we describe how upgraded and newly modelled data types in TargetMine can now survey the wider biological and chemical data space, relevant to drug discovery and development. To enhance the scope of TargetMine from target prioritization to broad-based knowledge discovery, we have also developed a new auxiliary toolkit to assist with data analysis and visualization in TargetMine. This toolkit features interactive data analysis tools to query and analyse the biological data compiled within the TargetMine data warehouse. The enhanced system enables users to discover new hypotheses interactively by performing complicated searches with no programming and obtaining the results in an easy to comprehend output format.

Database URL:http://targetmine.mizuguchilab.org

Evolution of Pentameric Ligand-Gated Ion Channels: Pro-Loop Receptors

Pentameric ligand-gated ion channels (pLGICs) are ubiquitous neurotransmitter receptors in Bilateria, with a small number of known prokaryotic homologues. Here we describe a new inventory and phylogenetic analysis of pLGIC genes across all kingdoms of life. Our main finding is a set of pLGIC genes in unicellular eukaryotes, some of which are metazoan-like Cys-loop receptors, and others devoid of Cys-loop cysteines, like their prokaryotic relatives. A number of such “Cys-less” receptors also appears in invertebrate metazoans. Together, those findings draw a new distribution of pLGICs in eukaryotes. A broader distribution of prokaryotic channels also emerges, including a major new archaeal taxon, Thaumarchaeota. More generally, pLGICs now appear nearly ubiquitous in major taxonomic groups except multicellular plants and fungi. However, pLGICs are sparsely present in unicellular taxa, suggesting a high rate of gene loss and a non-essential character, contrasting with their essential role as synaptic receptors of the bilaterian nervous system. Multiple alignments of these highly divergent sequences reveal a small number of conserved residues clustered at the interface between the extracellular and transmembrane domains. Only the “Cys-loop” proline is absolutely conserved, suggesting the more fitting name “Pro loop” for that motif, and “Pro-loop receptors” for the superfamily. The infered molecular phylogeny shows a Cys-loop and a Cys-less clade in eukaryotes, both containing metazoans and unicellular members. This suggests new hypotheses on the evolutionary history of the superfamily, such as a possible origin of the Cys-loop cysteines in an ancient unicellular eukaryote. Deeper phylogenetic relationships remain uncertain, particularly around the split between bacteria, archaea, and eukaryotes.

Genomic and Transcriptomic Features of Response to Anti-PD-1 Therapy in Metastatic Melanoma

PD-1 immune checkpoint blockade provides significant clinical benefits for melanoma patients. We analyzed the somatic mutanomes and transcriptomes of pretreatment melanoma biopsies to identify factors that may influence innate sensitivity or resistance to anti-PD-1 therapy. We find that overall high mutational loads associate with improved survival, and tumors from responding patients are enriched for mutations in the DNA repair gene BRCA2. Innately resistant tumors display a transcriptional signature (referred to as the IPRES, or innate anti-PD-1 resistance), indicating concurrent up-expression of genes involved in the regulation of mesenchymal transition, cell adhesion, extracellular matrix remodeling, angiogenesis, and wound healing. Notably, mitogen-activated protein kinase (MAPK)-targeted therapy (MAPK inhibitor) induces similar signatures in melanoma, suggesting that a non-genomic form of MAPK inhibitor resistance mediates cross-resistance to anti-PD-1 therapy. Validation of the IPRES in other independent tumor cohorts defines a transcriptomic subset across distinct types of advanced cancer. These findings suggest that attenuating the biological processes that underlie IPRES may improve anti-PD-1 response in melanoma and other cancer types.

PhyloPro2.0: a database for the dynamic exploration of phylogenetically conserved proteins and their domain architectures across the Eukarya

PhyloPro is a database and accompanying web-based application for the construction and exploration of phylogenetic profiles across the Eukarya. In this update article, we present six major new developments in PhyloPro: (i) integration of Pfam-A domain predictions for all proteins; (ii) new summary heatmaps and detailed level views of domain conservation; (iii) an interactive, network-based visualization tool for exploration of domain architectures and their conservation; (iv) ability to browse based on protein functional categories (GOSlim); (v) improvements to the web interface to enhance drill down capability from the heatmap view; and (vi) improved coverage including 164 eukaryotes and 12 reference species. In addition, we provide improved support for downloading data and images in a variety of formats. Among the existing tools available for phylogenetic profiles, PhyloPro provides several innovative domain-based features including a novel domain adjacency visualization tool. These are designed to allow the user to identify and compare proteins with similar domain architectures across species and thus develop hypotheses about the evolution of lineage-specific trajectories. Database URL: http://www.compsysbio.org/phylopro/.

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.

Genome-wide analysis identifies gain and loss/change of function within the small multigenic insecticidal Albumin 1 family of Medicago truncatula

BACKGROUND

Albumin 1b peptides (A1b) are small disulfide-knotted insecticidal peptides produced by Fabaceae (also called Leguminosae). To date, their diversity among this plant family has been essentially investigated through biochemical and PCR-based approaches. The availability of high-quality genomic resources for several fabaceae species, among which the model species Medicago truncatula (Mtr), allowed for a genomic analysis of this protein family aimed at i) deciphering the evolutionary history of A1b proteins and their links with A1b-nodulins that are short non-insecticidal disulfide-bonded peptides involved in root nodule signaling and ii) exploring the functional diversity of A1b for novel bioactive molecules.

RESULTS

Investigating the Mtr genome revealed a remarkable expansion, mainly through tandem duplications, of albumin1 (A1) genes, retaining nearly all of the same canonical structure at both gene and protein levels. Phylogenetic analysis revealed that the ancestral molecule was most probably insecticidal giving rise to, among others, A1b-nodulins. Expression meta-analysis revealed that many A1b coding genes are silent and a wide tissue distribution of the A1 transcripts/peptides within plant organs. Evolutionary rate analyses highlighted branches and sites with positive selection signatures, including two sites shown to be critical for insecticidal activity. Seven peptides were chemically synthesized and folded in vitro, then assayed for their biological activity. Among these, AG41 (aka MtrA1013 isoform, encoded by the orphan TA24778 contig.), showed an unexpectedly high insecticidal activity. The study highlights the unique burst of diversity of A1 peptides within the Medicago genus compared to the other taxa for which full-genomes are available: no A1 member in Lotus, or in red clover to date, while only a few are present in chick pea, soybean or pigeon pea genomes.

CONCLUSION

The expansion of the A1 family in the Medicago genus is reminiscent of the situation described for another disulfide-rich peptide family, the "Nodule-specific Cysteine-Rich" (NCR), discovered within the same species. The oldest insecticidal A1b toxin was described from the Sophorae, dating the birth of this seed-defense function to more than 58 million years, and making this model of plant/insect toxin/receptor (A1b/insect v-ATPase) one of the oldest known.

GenSeed-HMM: A Tool for Progressive Assembly Using Profile HMMs as Seeds and its Application in Alpavirinae Viral Discovery from Metagenomic Data

This work reports the development of GenSeed-HMM, a program that implements seed-driven progressive assembly, an approach to reconstruct specific sequences from unassembled data, starting from short nucleotide or protein seed sequences or profile Hidden Markov Models (HMM). The program can use any one of a number of sequence assemblers. Assembly is performed in multiple steps and relatively few reads are used in each cycle, consequently the program demands low computational resources. As a proof-of-concept and to demonstrate the power of HMM-driven progressive assemblies, GenSeed-HMM was applied to metagenomic datasets in the search for diverse ssDNA bacteriophages from the recently described Alpavirinae subfamily. Profile HMMs were built using Alpavirinae-specific regions from multiple sequence alignments (MSA) using either the viral protein 1 (VP1; major capsid protein) or VP4 (genome replication initiation protein). These profile HMMs were used by GenSeed-HMM (running Newbler assembler) as seeds to reconstruct viral genomes from sequencing datasets of human fecal samples. All contigs obtained were annotated and taxonomically classified using similarity searches and phylogenetic analyses. The most specific profile HMM seed enabled the reconstruction of 45 partial or complete Alpavirinae genomic sequences. A comparison with conventional (global) assembly of the same original dataset, using Newbler in a standalone execution, revealed that GenSeed-HMM outperformed global genomic assembly in several metrics employed. This approach is capable of detecting organisms that have not been used in the construction of the profile HMM, which opens up the possibility of diagnosing novel viruses, without previous specific information, constituting a de novo diagnosis. Additional applications include, but are not limited to, the specific assembly of extrachromosomal elements such as plastid and mitochondrial genomes from metagenomic data. Profile HMM seeds can also be used to reconstruct specific protein coding genes for gene diversity studies, and to determine all possible gene variants present in a metagenomic sample. Such surveys could be useful to detect the emergence of drug-resistance variants in sensitive environments such as hospitals and animal production facilities, where antibiotics are regularly used. Finally, GenSeed-HMM can be used as an adjunct for gap closure on assembly finishing projects, by using multiple contig ends as anchored seeds.

Regulation of amino-acid metabolism controls flux to lipid accumulation in Yarrowia lipolytica

Yarrowia lipolytica is a promising microbial cell factory for the production of lipids to be used as fuels and chemicals, but there are few studies on regulation of its metabolism. Here we performed the first integrated data analysis of Y. lipolytica grown in carbon and nitrogen limited chemostat cultures. We first reconstructed a genome-scale metabolic model and used this for integrative analysis of multilevel omics data. Metabolite profiling and lipidomics was used to quantify the cellular physiology, while regulatory changes were measured using RNAseq. Analysis of the data showed that lipid accumulation in Y. lipolytica does not involve transcriptional regulation of lipid metabolism but is associated with regulation of amino-acid biosynthesis, resulting in redirection of carbon flux during nitrogen limitation from amino acids to lipids. Lipid accumulation in Y. lipolytica at nitrogen limitation is similar to the overflow metabolism observed in many other microorganisms, e.g. ethanol production by Sacchromyces cerevisiae at nitrogen limitation.

An Ancient Fingerprint Indicates the Common Ancestry of Rossmann-Fold Enzymes Utilizing Different Ribose-Based Cofactors

Nucleoside-based cofactors are presumed to have preceded proteins. The Rossmann fold is one of the most ancient and functionally diverse protein folds, and most Rossmann enzymes utilize nucleoside-based cofactors. We analyzed an omnipresent Rossmann ribose-binding interaction: a carboxylate side chain at the tip of the second β-strand (β2-Asp/Glu). We identified a canonical motif, defined by the β2-topology and unique geometry. The latter relates to the interaction being bidentate (both ribose hydroxyls interacting with the carboxylate oxygens), to the angle between the carboxylate and the ribose, and to the ribose's ring configuration. We found that this canonical motif exhibits hallmarks of divergence rather than convergence. It is uniquely found in Rossmann enzymes that use different cofactors, primarily SAM (S-adenosyl methionine), NAD (nicotinamide adenine dinucleotide), and FAD (flavin adenine dinucleotide). Ribose-carboxylate bidentate interactions in other folds are not only rare but also have a different topology and geometry. We further show that the canonical geometry is not dictated by a physical constraint--geometries found in noncanonical interactions have similar calculated bond energies. Overall, these data indicate the divergence of several major Rossmann-fold enzyme classes, with different cofactors and catalytic chemistries, from a common pre-LUCA (last universal common ancestor) ancestor that possessed the β2-Asp/Glu motif.

Proteomic profiling of maize opaque endosperm mutants reveals selective accumulation of lysine-enriched proteins

Reduced prolamin (zein) accumulation and defective endoplasmic reticulum (ER) body formation occurs in maize opaque endosperm mutants opaque2 (o2), floury2 (fl2), defective endosperm*B30 (DeB30), and Mucronate (Mc), whereas other opaque mutants such as opaque1 (o1) and floury1 (fl1) are normal in these regards. This suggests that other factors contribute to kernel texture. A liquid chromatography approach coupled with tandem mass spectrometry (LC-MS/MS) proteomics was used to compare non-zein proteins of nearly isogenic opaque endosperm mutants. In total, 2762 proteins were identified that were enriched for biological processes such as protein transport and folding, amino acid biosynthesis, and proteolysis. Principal component analysis and pathway enrichment suggested that the mutants partitioned into three groups: (i) Mc, DeB30, fl2 and o2; (ii) o1; and (iii) fl1. Indicator species analysis revealed mutant-specific proteins, and highlighted ER secretory pathway components that were enriched in selected groups of mutants. The most significantly changed proteins were related to stress or defense and zein partitioning into the soluble fraction for Mc, DeB30, o1, and fl1 specifically. In silico dissection of the most significantly changed proteins revealed novel qualitative changes in lysine abundance contributing to the overall lysine increase and the nutritional rebalancing of the o2 and fl2 endosperm.

The Genomes of Three Uneven Siblings: Footprints of the Lifestyles of Three Trichoderma Species

The genus Trichoderma contains fungi with high relevance for humans, with applications in enzyme production for plant cell wall degradation and use in biocontrol. Here, we provide a broad, comprehensive overview of the genomic content of these species for "hot topic" research aspects, including CAZymes, transport, transcription factors, and development, along with a detailed analysis and annotation of less-studied topics, such as signal transduction, genome integrity, chromatin, photobiology, or lipid, sulfur, and nitrogen metabolism in T. reesei, T. atroviride, and T. virens, and we open up new perspectives to those topics discussed previously. In total, we covered more than 2,000 of the predicted 9,000 to 11,000 genes of each Trichoderma species discussed, which is >20% of the respective gene content. Additionally, we considered available transcriptome data for the annotated genes. Highlights of our analyses include overall carbohydrate cleavage preferences due to the different genomic contents and regulation of the respective genes. We found light regulation of many sulfur metabolic genes. Additionally, a new Golgi 1,2-mannosidase likely involved in N-linked glycosylation was detected, as were indications for the ability of Trichoderma spp. to generate hybrid galactose-containing N-linked glycans. The genomic inventory of effector proteins revealed numerous compounds unique to Trichoderma, and these warrant further investigation. We found interesting expansions in the Trichoderma genus in several signaling pathways, such as G-protein-coupled receptors, RAS GTPases, and casein kinases. A particularly interesting feature absolutely unique to T. atroviride is the duplication of the alternative sulfur amino acid synthesis pathway.

Acylsugar Acylhydrolases: Carboxylesterase-Catalyzed Hydrolysis of Acylsugars in Tomato Trichomes

Glandular trichomes of cultivated tomato (Solanum lycopersicum) and many other species throughout the Solanaceae produce and secrete mixtures of sugar esters (acylsugars) on the plant aerial surfaces. In wild and cultivated tomato, these metabolites consist of a sugar backbone, typically glucose or sucrose, and two to five acyl chains esterified to various positions on the sugar core. The aliphatic acyl chains vary in length and branching and are transferred to the sugar by a series of reactions catalyzed by acylsugar acyltransferases. A phenotypic screen of a set of S. lycopersicum M82 × Solanum pennellii LA0716 introgression lines identified a dominant genetic locus on chromosome 5 from the wild relative that affected total acylsugar levels. Genetic mapping revealed that the reduction in acylsugar levels was consistent with the presence and increased expression of two S. pennellii genes (Sopen05g030120 and Sopen05g030130) encoding putative carboxylesterase enzymes of the α/β-hydrolase superfamily. These two enzymes, named ACYLSUGAR ACYLHYDROLASE1 (ASH1) and ASH2, were shown to remove acyl chains from specific positions of certain types of acylsugars in vitro. A survey of related genes in M82 and LA0716 identified another trichome-expressed ASH gene on chromosome 9 (M82, Solyc09g075710; LA0716, Sopen09g030520) encoding a protein with similar activity. Characterization of the in vitro activities of the SpASH enzymes showed reduced activities with acylsugars produced by LA0716, presumably contributing to the high-level production of acylsugars in the presence of highly expressed SpASH genes.

The genomic basis of parasitism in the Strongyloides clade of nematodes

Soil-transmitted nematodes, including the Strongyloides genus, cause one of the most prevalent neglected tropical diseases. Here we compare the genomes of four Strongyloides species, including the human pathogen Strongyloides stercoralis, and their close relatives that are facultatively parasitic (Parastrongyloides trichosuri) and free-living (Rhabditophanes sp. KR3021). A significant paralogous expansion of key gene families--families encoding astacin-like and SCP/TAPS proteins--is associated with the evolution of parasitism in this clade. Exploiting the unique Strongyloides life cycle, we compare the transcriptomes of the parasitic and free-living stages and find that these same gene families are upregulated in the parasitic stages, underscoring their role in nematode parasitism.

Venom gland transcriptome analyses of two freshwater stingrays (Myliobatiformes: Potamotrygonidae) from Brazil

Stingrays commonly cause human envenoming related accidents in populations of the sea, near rivers and lakes. Transcriptomic profiles have been used to elucidate components of animal venom, since they are capable of providing molecular information on the biology of the animal and could have biomedical applications. In this study, we elucidated the transcriptomic profile of the venom glands from two different freshwater stingray species that are endemic to the Paraná-Paraguay basin in Brazil, Potamotrygon amandae and Potamotrygon falkneri. Using RNA-Seq, we identified species-specific transcripts and overlapping proteins in the venom gland of both species. Among the transcripts related with envenoming, high abundance of hyaluronidases was observed in both species. In addition, we built three-dimensional homology models based on several venom transcripts identified. Our study represents a significant improvement in the information about the venoms employed by these two species and their molecular characteristics. Moreover, the information generated by our group helps in a better understanding of the biology of freshwater cartilaginous fishes and offers clues for the development of clinical treatments for stingray envenoming in Brazil and around the world. Finally, our results might have biomedical implications in developing treatments for complex diseases.

Structure-based Markov random field model for representing evolutionary constraints on functional sites

Background

Elucidating the cooperative mechanism of interconnected residues is an important component toward understanding the biological function of a protein. Coevolution analysis has been developed to model the coevolutionary information reflecting structural and functional constraints. Recently, several methods have been developed based on a probabilistic graphical model called the Markov random field (MRF), which have led to significant improvements for coevolution analysis; however, thus far, the performance of these models has mainly been assessed by focusing on the aspect of protein structure.

Results

In this study, we built an MRF model whose graphical topology is determined by the residue proximity in the protein structure, and derived a novel positional coevolution estimate utilizing the node weight of the MRF model. This structure-based MRF method was evaluated for three data sets, each of which annotates catalytic site, allosteric site, and comprehensively determined functional site information. We demonstrate that the structure-based MRF architecture can encode the evolutionary information associated with biological function. Furthermore, we show that the node weight can more accurately represent positional coevolution information compared to the edge weight. Lastly, we demonstrate that the structure-based MRF model can be reliably built with only a few aligned sequences in linear time.

Conclusions

The results show that adoption of a structure-based architecture could be an acceptable approximation for coevolution modeling with efficient computation complexity.

Identification of repetitive units in protein structures with ReUPred

Over the last decade, numerous studies have demonstrated the fundamental importance of tandem repeat (TR) proteins in many biological processes. A plethora of new repeat structures have also been solved. The recently published RepeatsDB provides information on TR proteins. However, a detailed structural characterization of repetitive elements is largely missing, as repeat unit annotation is manually curated and currently covers only 3 % of the bona fide TR proteins. Repeat Protein Unit Predictor (ReUPred) is a novel method for the fast automatic prediction of repeat units and repeat classification using an extensive Structure Repeat Unit Library (SRUL) derived from RepeatsDB. ReUPred uses an iterative structural search against the SRUL to find repetitive units. On a test set of solenoid proteins, ReUPred is able to correctly detect 92 % of the proteins. Unlike previous methods, it is also able to correctly classify solenoid repeats in 89 % of cases. It also outperforms two recent state-of-the-art methods for the repeat unit identification problem. The accurate prediction of repeat units increases the number of annotated repeat units by an order of magnitude compared to the sequence-based Pfam classification. ReUPred is implemented in Python for Linux and freely available from the URL: http://protein.bio.unipd.it/reupred/ .

Roles of Tetratricopeptide Repeat Proteins in Biogenesis of the Photosynthetic Apparatus

Biosynthesis of the photosynthetic apparatus is a complex operation, which includes the concerted synthesis and assembly of lipids, pigments and metal cofactors, and dozens of proteins. Research conducted in recent years has shown that these processes, as well as the stabilization and repair of this molecular machinery, are facilitated by transiently acting regulatory proteins, many of which belong to the superfamily of helical repeat proteins. Here, we focus on one of its families in photoautotrophic model organisms, the tetratricopeptide repeat (TPR) proteins, which participate in almost all of these steps and are crucial for biogenesis of the thylakoid membrane.

The handedness-associated PCSK6 locus spans an intronic promoter regulating novel transcripts

We recently reported the association of the PCSK6 gene with handedness through a quantitative genome-wide association study (GWAS; P < 0.5 × 10⁻⁸) for a relative hand skill measure in individuals with dyslexia. PCSK6 activates Nodal, a morphogen involved in regulating left–right body axis determination. Therefore, the GWAS data suggest that the biology underlying the patterning of structural asymmetries may also contribute to behavioural laterality, e.g. handedness. The association is further supported by an independent study reporting a variable number tandem repeat (VNTR) within the same PCSK6 locus to be associated with degree of handedness in a general population cohort. Here, we have conducted a functional analysis of the PCSK6 locus combining further genetic analysis, in silico predictions and molecular assays. We have shown that the previous GWAS signal was not tagging a VNTR effect, suggesting that the two markers have independent effects. We demonstrated experimentally that one of the top GWAS-associated markers, rs11855145, directly alters the binding site for a nuclear factor. Furthermore, we have shown that the predicted regulatory region adjacent to rs11855415 acts as a bidirectional promoter controlling the expression of novel RNA transcripts. These include both an antisense long non-coding RNA (lncRNA) and a short PCSK6 isoform predicted to be coding. This is the first molecular characterization of a handedness-associated locus that supports the role of common variants in non-coding sequences in influencing complex phenotypes through gene expression regulation.

The Genetic Basis of Haploid Induction in Maize Identified with a Novel Genome-Wide Association Method

In vivo haploid induction (HI) triggered by pollination with special intraspecific genotypes, called inducers, is unique to Zea maysL. within the plant kingdom and has revolutionized maize breeding during the last decade. However, the molecular mechanisms underlying HI in maize are still unclear. To investigate the genetic basis of HI, we developed a new approach for genome-wide association studies (GWAS), termed conditional haplotype extension (CHE) test that allows detection of selective sweeps even under almost perfect confounding of population structure and trait expression. Here, we applied this test to identify genomic regions required for HI expression and dissected the combined support interval (50.34 Mb) of the QTL qhir1, detected in a previous study, into two closely linked genomic segments relevant for HI expression. The first, termed qhir11(0.54 Mb), comprises an already fine-mapped region but was not diagnostic for differentiating inducers and noninducers. The second segment, termed qhir12(3.97 Mb), had a haplotype allele common to all 53 inducer lines but not found in any of the 1482 noninducers. By comparing resequencing data of one inducer with 14 noninducers, we detected in the qhir12 region three candidate genes involved in DNA or amino acid binding, however, none for qhir11 We propose that the CHE test can be utilized in introgression breeding and different fields of genetics to detect selective sweeps in heterogeneous genetic backgrounds.

High frequency of +1 programmed ribosomal frameshifting in Euplotes octocarinatus

Programmed -1 ribosomal frameshifting (-1 PRF) has been identified as a mechanism to regulate the expression of many viral genes and some cellular genes. The slippery site of -1 PRF has been well characterized, whereas the +1 PRF signal and the mechanism involved in +1 PRF remain poorly understood. Previous study confirmed that +1 PRF is required for the synthesis of protein products in several genes of ciliates from the genus Euplotes. To accurately assess the frequency of genes requiring frameshift in Euplotes, the macronuclear genome and transcriptome of Euplotes octocarinatus were analyzed in this study. A total of 3,700 +1 PRF candidate genes were identified from 32,353 transcripts, and the gene products of these putative +1 PRFs were mainly identified as protein kinases. Furthermore, we reported a putative suppressor tRNA of UAA which may provide new insights into the mechanism of +1 PRF in euplotids. For the first time, our transcriptome-wide survey of +1 PRF in E. octocarinatus provided a dataset which serves as a valuable resource for the future understanding of the mechanism underlying +1 PRF.

Understanding the Role of the Master Regulator XYR1 in Trichoderma reesei by Global Transcriptional Analysis

We defined the role of the transcriptional factor—XYR1—in the filamentous fungus Trichoderma reesei during cellulosic material degradation. In this regard, we performed a global transcriptome analysis using RNA-Seq of the Δxyr1 mutant strain of T. reesei compared with the parental strain QM9414 grown in the presence of cellulose, sophorose, and glucose as sole carbon sources. We found that 5885 genes were expressed differentially under the three tested carbon sources. Of these, 322 genes were upregulated in the presence of cellulose, while 367 and 188 were upregulated in sophorose and glucose, respectively. With respect to genes under the direct regulation of XYR1, 30 and 33 are exclusive to cellulose and sophorose, respectively. The most modulated genes in the Δxyr1 belong to Carbohydrate-Active Enzymes (CAZymes), transcription factors, and transporters families. Moreover, we highlight the downregulation of transporters belonging to the MFS and ABC transporter families. Of these, MFS members were mostly downregulated in the presence of cellulose. In sophorose and glucose, the expression of these transporters was mainly upregulated. Our results revealed that MFS and ABC transporters could be new players in cellulose degradation and their role was shown to be carbon source-dependent. Our findings contribute to a better understanding of the regulatory mechanisms of XYR1 to control cellulase gene expression in T. reesei in the presence of cellulosic material, thereby potentially enhancing its application in several biotechnology fields.

In Silico Analysis of the Structural and Biochemical Features of the NMD Factor UPF1 in Ustilago maydis

The molecular mechanisms regulating the accuracy of gene expression are still not fully understood. Among these mechanisms, Nonsense-mediated Decay (NMD) is a quality control process that detects post-transcriptionally abnormal transcripts and leads them to degradation. The UPF1 protein lays at the heart of NMD as shown by several structural and functional features reported for this factor mainly for Homo sapiens and Saccharomyces cerevisiae. This process is highly conserved in eukaryotes but functional diversity can be observed in various species. Ustilago maydis is a basidiomycete and the best-known smut, which has become a model to study molecular and cellular eukaryotic mechanisms. In this study, we performed in silico analysis to investigate the structural and biochemical properties of the putative UPF1 homolog in Ustilago maydis. The putative homolog for UPF1 was recognized in the annotated genome for the basidiomycete, exhibiting 66% identity with its human counterpart at the protein level. The known structural and functional domains characteristic of UPF1 homologs were also found. Based on the crystal structures available for UPF1, we constructed different three-dimensional models for umUPF1 in order to analyze the secondary and tertiary structural features of this factor. Using these models, we studied the spatial arrangement of umUPF1 and its capability to interact with UPF2. Moreover, we identified the critical amino acids that mediate the interaction of umUPF1 with UPF2, ATP, RNA and with UPF1 itself. Mutating these amino acids in silico showed an important effect over the native structure. Finally, we performed molecular dynamic simulations for UPF1 proteins from H. sapiens and U. maydis and the results obtained show a similar behavior and physicochemical properties for the protein in both organisms. Overall, our results indicate that the putative UPF1 identified in U. maydis shows a very similar sequence, structural organization, mechanical stability, physicochemical properties and spatial organization in comparison to the NMD factor depicted for Homo sapiens. These observations strongly support the notion that human and fungal UPF1 could perform equivalent biological activities.

A review of the new HGNC gene family resource

The HUGO Gene Nomenclature Committee (HGNC) approves unique gene symbols and names for human loci. As well as naming genomic loci, we manually curate genes into family sets based on shared characteristics such as function, homology or phenotype. Each HGNC gene family has its own dedicated gene family report on our website, www.genenames.org. We have recently redesigned these reports to support the visualisation and browsing of complex relationships between families and to provide extra curated information such as family descriptions, protein domain graphics and gene family aliases. Here, we review how our gene families are curated and explain how to view, search and download the gene family data.

Transcriptome-Wide Identification and Expression Profiling Analysis of Chrysanthemum Trihelix Transcription Factors

Trihelix transcription factors are thought to feature a typical DNA-binding trihelix (helix-loop-helix-loop-helix) domain that binds specifically to the GT motif, a light-responsive DNA element. Members of the trihelix family are known to function in a number of processes in plants. Here, we characterize 20 trihelix family genes in the important ornamental plant chrysanthemum (Chrysanthemum morifolium). Based on transcriptomic data, 20 distinct sequences distributed across four of five groups revealed by a phylogenetic tree were isolated and amplified. The phylogenetic analysis also identified four pairs of orthologous proteins shared by Arabidopsis and chrysanthemum and five pairs of paralogous proteins in chrysanthemum. Conserved motifs in the trihelix proteins shared by Arabidopsis and chrysanthemum were analyzed using MEME, and further bioinformatic analysis revealed that 16 CmTHs can be targeted by 20 miRNA families and that miR414 can target 9 CmTHs. qPCR results displayed that most chrysanthemum trihelix genes were highly expressed in inflorescences, while 20 CmTH genes were in response to phytohormone treatments and abiotic stresses. This work improves our understanding of the various functions of trihelix gene family members in response to hormonal stimuli and stress.

Modulation of FAK and Src adhesion signaling occurs independently of adhesion complex composition

Integrin adhesion complexes (IACs) form mechanochemical connections between the extracellular matrix and actin cytoskeleton and mediate phenotypic responses via posttranslational modifications. Here, we investigate the modularity and robustness of the IAC network to pharmacological perturbation of the key IAC signaling components focal adhesion kinase (FAK) and Src. FAK inhibition using AZ13256675 blocked FAK(Y397) phosphorylation but did not alter IAC composition, as reported by mass spectrometry. IAC composition was also insensitive to Src inhibition using AZD0530 alone or in combination with FAK inhibition. In contrast, kinase inhibition substantially reduced phosphorylation within IACs, cell migration and proliferation. Furthermore using fluorescence recovery after photobleaching, we found that FAK inhibition increased the exchange rate of a phosphotyrosine (pY) reporter (dSH2) at IACs. These data demonstrate that kinase-dependent signal propagation through IACs is independent of gross changes in IAC composition. Together, these findings demonstrate a general separation between the composition of IACs and their ability to relay pY-dependent signals.

The genome of the seagrass Zostera marina reveals angiosperm adaptation to the sea

Seagrasses colonized the sea on at least three independent occasions to form the basis of one of the most productive and widespread coastal ecosystems on the planet. Here we report the genome of Zostera marina (L.), the first, to our knowledge, marine angiosperm to be fully sequenced. This reveals unique insights into the genomic losses and gains involved in achieving the structural and physiological adaptations required for its marine lifestyle, arguably the most severe habitat shift ever accomplished by flowering plants. Key angiosperm innovations that were lost include the entire repertoire of stomatal genes, genes involved in the synthesis of terpenoids and ethylene signalling, and genes for ultraviolet protection and phytochromes for far-red sensing. Seagrasses have also regained functions enabling them to adjust to full salinity. Their cell walls contain all of the polysaccharides typical of land plants, but also contain polyanionic, low-methylated pectins and sulfated galactans, a feature shared with the cell walls of all macroalgae and that is important for ion homoeostasis, nutrient uptake and O2/CO2 exchange through leaf epidermal cells. The Z. marina genome resource will markedly advance a wide range of functional ecological studies from adaptation of marine ecosystems under climate warming, to unravelling the mechanisms of osmoregulation under high salinities that may further inform our understanding of the evolution of salt tolerance in crop plants.

Comprehensive annotation of Glossina pallidipes salivary gland hypertrophy virus from Ethiopian tsetse flies: a proteogenomics approach

Glossina pallidipes salivary gland hypertrophy virus (GpSGHV; family Hytrosaviridae) can establish asymptomatic and symptomatic infection in its tsetse fly host. Here, we present a comprehensive annotation of the genome of an Ethiopian GpSGHV isolate (GpSGHV-Eth) compared with the reference Ugandan GpSGHV isolate (GpSGHV-Uga; GenBank accession number EF568108). GpSGHV-Eth has higher salivary gland hypertrophy syndrome prevalence than GpSGHV-Uga. We show that the GpSGHV-Eth genome has 190 291 nt, a low G+C content (27.9 %) and encodes 174 putative ORFs. Using proteogenomic and transcriptome mapping, 141 and 86 ORFs were mapped by transcripts and peptides, respectively. Furthermore, of the 174 ORFs, 132 had putative transcriptional signals [TATA-like box and poly(A) signals]. Sixty ORFs had both TATA-like box promoter and poly(A) signals, and mapped by both transcripts and peptides, implying that these ORFs encode functional proteins. Of the 60 ORFs, 10 ORFs are homologues to baculovirus and nudivirus core genes, including three per os infectivity factors and four RNA polymerase subunits (LEF4, 5, 8 and 9). Whereas GpSGHV-Eth and GpSGHV-Uga are 98.1 % similar at the nucleotide level, 37 ORFs in the GpSGHV-Eth genome had nucleotide insertions (n = 17) and deletions (n = 20) compared with their homologues in GpSGHV-Uga. Furthermore, compared with the GpSGHV-Uga genome, 11 and 24 GpSGHV ORFs were deleted and novel, respectively. Further, 13 GpSGHV-Eth ORFs were non-canonical; they had either CTG or TTG start codons instead of ATG. Taken together, these data suggest that GpSGHV-Eth and GpSGHV-Uga represent two different lineages of the same virus. Genetic differences combined with host and environmental factors possibly explain the differential GpSGHV pathogenesis observed in different G. pallidipes colonies.

Gene expression and adaptive evolution of ZBED1 in the hibernating greater horseshoe bat (Rhinolophus ferrumequinum)

Mammalian hibernators experience physiological extremes, e.g. ischemia, muscle disuse and hypothermia, which are lethal to non-hibernators, implying the existence of underlying mechanisms that allow hibernators to withstand these physiological extremes. Increased cell proliferation is suggested to be such a strategy, but its molecular basis remains unknown. In this study, we characterized the expression pattern of ZBED1 (zinc finger, BED-type containing 1), a transcription factor that plays a crucial role in regulating cell proliferation, in five tissues of the greater horseshoe bat (Rhinolophus ferrumequinum) during pre-hibernation, deep hibernation and post-hibernation. Moreover, we investigated the ZBED1 genetic divergence from individuals with variable hibernation phenotypes that cover all three known mtDNA lineages of the species. Expression analyses showed that ZBED1 is overexpressed only in brain and skeletal muscle, not in the other three tissues, suggesting an increased cell proliferation in these two tissues during deep hibernation. Evolutionary analyses showed that ZBED1 sequences were clustered into two well-supported clades with each one dominated by hibernating and non-hibernating individuals, respectively. Positive selection analyses further showed some positively selected sites and a divergent selection pressure among hibernating and non-hibernating groups of R. ferrumequinum. Our results suggest that ZBED1 as a potential candidate gene that regulates cell proliferation for hibernators to face physiological extremes during hibernation.

Circuits of cancer drivers revealed by convergent misregulation of transcription factor targets across tumor types

BACKGROUND

Large tumor genome sequencing projects have now uncovered a few hundred genes involved in the onset of tumorigenesis, or drivers, in some two dozen malignancies. One of the main challenges emerging from this catalog of drivers is how to make sense of their heterogeneity in most cancer types. This is key not only to understand how carcinogenesis appears and develops in these malignancies to be able to early diagnose them, but also to open up the possibility to employ therapeutic strategies targeting a driver protein to counteract the alteration of another connected driver.

METHODS

Here, I focus on driver transcription factors and their connection to tumorigensis in several tumor types through the alteration of the expression of their targets. First, I explore their involvement in tumorigenesis as mutational drivers in 28 different tumor types. Then, I collect a list of downstream targets of the all driver transcription factors (TFs), and identify which of them exhibit a differential expression upon alterations of driver transcription factors.

RESULTS

I identify the subset of targets of each TF most likely mediating the tumorigenic effect of their driver alterations in each tumor type, and explore their overlap. Furthermore, I am able to identify other driver genes that cause tumorigenesis through the alteration of very similar sets of targets.

CONCLUSIONS

I thus uncover these circuits of connected drivers which cause tumorigenesis through the perturbation of overlapping cellular pathways in a pan-cancer manner across 15 malignancies. The systematic detection of these circuits may be key to propose novel therapeutic strategies indirectly targeting driver alterations in tumors.

Dynamin-2 is a novel NOS1β interacting protein and negative regulator in the collecting duct

Nitric oxide synthase 1 (NOS1)-derived nitric oxide (NO) production in collecting ducts is critical for maintaining fluid-electrolyte balance. Rat collecting ducts express both the full-length NOS1α and its truncated variant NOS1β, while NOS1β predominates in mouse collecting ducts. We reported that dynamin-2 (DNM2), a protein involved in excising vesicles from the plasma membrane, and NOS1α form a protein-protein interaction that promotes NO production in rat collecting ducts. NOS1β was found to be highly expressed in human renal cortical/medullary samples; hence, we tested the hypothesis that DNM2 is a positive regulator of NOS1β-derived NO production. COS7 and mouse inner medullary collecting duct-3 (mIMCD3) cells were transfected with NOS1β and/or DNM2. Coimmunoprecipitation experiments show that NOS1β and DNM2 formed a protein-protein interaction. DNM2 overexpression decreased nitrite production (index of NO) in both COS7 and mIMCD-3 cells by 50-75%. mIMCD-3 cells treated with a panel of dynamin inhibitors or DNM2 siRNA displayed increased nitrite production. To elucidate the physiological significance of IMCD DNM2/NOS1β regulation in vivo, flox control and CDNOS1 knockout mice were placed on a high-salt diet, and freshly isolated IMCDs were treated acutely with a dynamin inhibitor. Dynamin inhibition increased nitrite production by IMCDs from flox mice. This response was blunted (but not abolished) in collecting duct-specific NOS1 knockout mice, suggesting that DNM2 also negatively regulates NOS3 in the mouse IMCD. We conclude that DNM2 is a novel negative regulator of NO production in mouse collecting ducts. We propose that DNM2 acts as a "break" to prevent excess or potentially toxic NO levels under high-salt conditions.

Effects of Halide Ions on the Carbamidocyclophane Biosynthesis in Nostoc sp. CAVN2

In this study, the influence of halide ions on [7.7]paracyclophane biosynthesis in the cyanobacterium Nostoc sp. CAVN2 was investigated. In contrast to KI and KF, supplementation of the culture medium with KCl or KBr resulted not only in an increase of growth but also in an up-regulation of carbamidocyclophane production. LC-MS analysis indicated the presence of chlorinated, brominated, but also non-halogenated derivatives. In addition to 22 known cylindrocyclophanes and carbamidocyclophanes, 27 putative congeners have been detected. Nine compounds, carbamidocyclophanes M-U, were isolated, and their structural elucidation by 1D and 2D NMR experiments in combination with HRMS and ECD analysis revealed that they are brominated analogues of chlorinated carbamidocyclophanes. Quantification of the carbamidocyclophanes showed that chloride is the preferably utilized halide, but incorporation is reduced in the presence of bromide. Evaluation of the antibacterial activity of 30 [7.7]paracyclophanes and related derivatives against selected pathogenic Gram-positive and Gram-negative bacteria exhibited remarkable effects especially against methicillin- and vancomycin-resistant staphylococci and Mycobacterium tuberculosis. For deeper insights into the mechanisms of biosynthesis, the carbamidocyclophane biosynthetic gene cluster in Nostoc sp. CAVN2 was studied. The gene putatively coding for the carbamoyltransferase has been identified. Based on bioinformatic analyses, a possible biosynthetic assembly is discussed.

Bacillus subtilis Spore Inner Membrane Proteome

The endospore is the dormant form of Bacillus subtilis and many other Firmicutes. By sporulation, these spore formers can survive very harsh physical and chemical conditions. Yet, they need to go through germination to return to their growing form. The spore inner membrane (IM) has been shown to play an essential role in triggering the initiation of germination. In this study, we isolated the IM of bacterial spores, in parallel with the isolation of the membrane of vegetative cells. With the use of GeLC-MS/MS, over 900 proteins were identified from the B. subtilis spore IM preparations. By bioinformatics-based membrane protein predictions, ca. one-third could be predicted to be membrane-localized. A large number of unique proteins as well as proteins common to the two membrane proteomes were identified. In addition to previously known IM proteins, a number of IM proteins were newly identified, at least some of which are likely to provide new insights into IM physiology, unveiling proteins putatively involved in spore germination machinery and hence putative germination inhibition targets.

Use of a Probabilistic Motif Search to Identify Histidine Phosphotransfer Domain-Containing Proteins

The wealth of newly obtained proteomic information affords researchers the possibility of searching for proteins of a given structure or function. Here we describe a general method for the detection of a protein domain of interest in any species for which a complete proteome exists. In particular, we apply this approach to identify histidine phosphotransfer (HPt) domain-containing proteins across a range of eukaryotic species. From the sequences of known HPt domains, we created an amino acid occurrence matrix which we then used to define a conserved, probabilistic motif. Examination of various organisms either known to contain (plant and fungal species) or believed to lack (mammals) HPt domains established criteria by which new HPt candidates were identified and ranked. Search results using a probabilistic motif matrix compare favorably with data to be found in several commonly used protein structure/function databases: our method identified all known HPt proteins in the Arabidopsis thaliana proteome, confirmed the absence of such motifs in mice and humans, and suggests new candidate HPts in several organisms. Moreover, probabilistic motif searching can be applied more generally, in a manner both readily customized and computationally compact, to other protein domains; this utility is demonstrated by our identification of histones in a range of eukaryotic organisms.

Assembly, Assessment, and Availability of De novo Generated Eukaryotic Transcriptomes

De novo assembly of a complete transcriptome without the need for a guiding reference genome is attractive, particularly where the cost and complexity of generating a eukaryote genome is prohibitive. The transcriptome should not however be seen as just a quick and cheap alternative to building a complete genome. Transcriptomics allows the understanding and comparison of spatial and temporal samples within an organism, and allows surveying of multiple individuals or closely related species. De novo assembly in theory allows the building of a complete transcriptome without any prior knowledge of the genome. It also allows the discovery of alternate splice forms of coding RNAs and also non-coding RNAs, which are often missed by proteomic approaches, or are incompletely annotated in genome studies. The limitations of the method are that the generation of a truly complete assembly is unlikely, and so we require some methods for the assessment of the quality and appropriateness of a generated transcriptome. Whilst no single consensus pipeline or tool is agreed as optimal, various algorithms, and easy to use software do exist making transcriptome generation a more common approach. With this expansion of data, questions still exist relating to how do we make these datasets fully discoverable, comparable and most useful to understand complex biological systems?

Cell wall proteome of sugarcane stems: comparison of a destructive and a non-destructive extraction method showed differences in glycoside hydrolases and peroxidases

BACKGROUND

Sugarcane has been used as the main crop for ethanol production for more than 40 years in Brazil. Recently, the production of bioethanol from bagasse and straw, also called second generation (2G) ethanol, became a reality with the first commercial plants started in the USA and Brazil. However, the industrial processes still need to be improved to generate a low cost fuel. One possibility is the remodeling of cell walls, by means of genetic improvement or transgenesis, in order to make the bagasse more accessible to hydrolytic enzymes. We aimed at characterizing the cell wall proteome of young sugarcane culms, to identify proteins involved in cell wall biogenesis. Proteins were extracted from the cell walls of 2-month-old culms using two protocols, non-destructive by vacuum infiltration vs destructive. The proteins were identified by mass spectrometry and bioinformatics.

RESULTS

A predicted signal peptide was found in 84 different proteins, called cell wall proteins (CWPs). As expected, the non-destructive method showed a lower percentage of proteins predicted to be intracellular than the destructive one (33% vs 44%). About 19% of CWPs were identified with both methods, whilst the infiltration protocol could lead to the identification of 75% more CWPs. In both cases, the most populated protein functional classes were those of proteins related to lipid metabolism and oxido-reductases. Curiously, a single glycoside hydrolase (GH) was identified using the non-destructive method whereas 10 GHs were found with the destructive one. Quantitative data analysis allowed the identification of the most abundant proteins.

CONCLUSIONS

The results highlighted the importance of using different protocols to extract proteins from cell walls to expand the coverage of the cell wall proteome. Ten GHs were indicated as possible targets for further studies in order to obtain cell walls less recalcitrant to deconstruction. Therefore, this work contributed to two goals: enlarge the coverage of the sugarcane cell wall proteome, and provide target proteins that could be used in future research to facilitate 2G ethanol production.

Integrating 400 million variants from 80,000 human samples with extensive annotations: towards a knowledge base to analyze disease cohorts

BACKGROUND

Data from a plethora of high-throughput sequencing studies is readily available to researchers, providing genetic variants detected in a variety of healthy and disease populations. While each individual cohort helps gain insights into polymorphic and disease-associated variants, a joint perspective can be more powerful in identifying polymorphisms, rare variants, disease-associations, genetic burden, somatic variants, and disease mechanisms.

DESCRIPTION

We have set up a Reference Variant Store (RVS) containing variants observed in a number of large-scale sequencing efforts, such as 1000 Genomes, ExAC, Scripps Wellderly, UK10K; various genotyping studies; and disease association databases. RVS holds extensive annotations pertaining to affected genes, functional impacts, disease associations, and population frequencies. RVS currently stores 400 million distinct variants observed in more than 80,000 human samples.

CONCLUSIONS

RVS facilitates cross-study analysis to discover novel genetic risk factors, gene-disease associations, potential disease mechanisms, and actionable variants. Due to its large reference populations, RVS can also be employed for variant filtration and gene prioritization.

AVAILABILITY

A web interface to public datasets and annotations in RVS is available at https://rvs.u.hpc.mssm.edu/.

Genome-wide analysis of tomato NF-Y factors and their role in fruit ripening

Background

Fruit ripening is a complex developmental process that depends on a coordinated regulation of numerous genes, including ripening-related transcription factors (TFs), fruit-related microRNAs, DNA methylation and chromatin remodeling. It is known that various TFs, such as MADS-domain, MYB, AP2/ERF and SBP/SPL family proteins play key roles in modulating ripening. However, little attention has been given to members of the large NF-Y TF family in this regard, although genes in this family are known to have important functions in regulating plant growth, development, and abiotic or biotic stress responses.

Results

In this study, the evolutionary relationship between Arabidopsis thaliana and tomato (Solanum lycopersicum) NF-Y genes was examined to predict similarities in function. Furthermore, through gene expression analysis, 13 tomato NF-Y genes were identified as candidate regulators of fruit ripening. Functional studies involving suppression of NF-Y gene expression using virus induced gene silencing (VIGS) indicated that five NF-Y genes, including two members of the NF-YB subgroup (Solyc06g069310, Solyc07g065500) and three members of the NF-YA subgroup (Solyc01g087240, Solyc08g062210, Solyc11g065700), influence ripening. In addition, subcellular localization analyses using NF-Y proteins fused to a green fluorescent protein (GFP) reporter showed that the three NF-YA proteins accumulated in the nucleus, while the two NF-YB proteins were observed in both the nucleus and cytoplasm.

Conclusions

In this study, we identified tomato NF-Y genes by analyzing the tomato genome sequence using bioinformatics approaches, and characterized their chromosomal distribution, gene structures, phylogenetic relationship and expression patterns. We also examined their biological functions in regulating tomato fruit via VIGS and subcellular localization analyses. The results indicated that five NF-Y transcription factors play roles in tomato fruit ripening. This information provides a platform for further investigation of their biological functions.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2334-2) contains supplementary material, which is available to authorized users.

3DIANA: 3D Domain Interaction Analysis: A Toolbox for Quaternary Structure Modeling

Electron microscopy (EM) is experiencing a revolution with the advent of a new generation of Direct Electron Detectors, enabling a broad range of large and flexible structures to be resolved well below 1 nm resolution. Although EM techniques are evolving to the point of directly obtaining structural data at near-atomic resolution, for many molecules the attainable resolution might not be enough to propose high-resolution structural models. However, accessing information on atomic coordinates is a necessary step toward a deeper understanding of the molecular mechanisms that allow proteins to perform specific tasks. For that reason, methods for the integration of EM three-dimensional maps with x-ray and NMR structural data are being developed, a modeling task that is normally referred to as fitting, resulting in the so called hybrid models. In this work, we present a novel application—3DIANA—specially targeted to those cases in which the EM map resolution is medium or low and additional experimental structural information is scarce or even lacking. In this way, 3DIANA statistically evaluates proposed/potential contacts between protein domains, presents a complete catalog of both structurally resolved and predicted interacting regions involving these domains and, finally, suggests structural templates to model the interaction between them. The evaluation of the proposed interactions is computed with DIMERO, a new method that scores physical binding sites based on the topology of protein interaction networks, which has recently shown the capability to increase by 200% the number of domain-domain interactions predicted in interactomes as compared to previous approaches. The new application displays the information at a sequence and structural level and is accessible through a web browser or as a Chimera plugin at http://3diana.cnb.csic.es.

ANISEED 2015: a digital framework for the comparative developmental biology of ascidians

Ascidians belong to the tunicates, the sister group of vertebrates and are recognized model organisms in the field of embryonic development, regeneration and stem cells. ANISEED is the main information system in the field of ascidian developmental biology. This article reports the development of the system since its initial publication in 2010. Over the past five years, we refactored the system from an initial custom schema to an extended version of the Chado schema and redesigned all user and back end interfaces. This new architecture was used to improve and enrich the description of Ciona intestinalis embryonic development, based on an improved genome assembly and gene model set, refined functional gene annotation, and anatomical ontologies, and a new collection of full ORF cDNAs. The genomes of nine ascidian species have been sequenced since the release of the C. intestinalis genome. In ANISEED 2015, all nine new ascidian species can be explored via dedicated genome browsers, and searched by Blast. In addition, ANISEED provides full functional gene annotation, anatomical ontologies and some gene expression data for the six species with highest quality genomes. ANISEED is publicly available at: http://www.aniseed.cnrs.fr.