InterPro protein classification

The Secretome of Adult Murine Hookworms Is Shaped by Host Expression of STAT6

Co-evolutionary adaptation of hookworms with their mammalian hosts has been selected for immunoregulatory excretory/secretory (E/S) products. However, it is not known whether, or if so, how host immunological status impacts the secreted profile of hematophagous adult worms. This study interrogated the impact of host Signal transducer and activator of transcription 6 (STAT6) expression during the experimental evolution of hookworms through the sequential passage of the life cycle in either STAT6 deficient or WT C57BL/6 mice. Proteomic analysis of E/S products by LC-MS showed increased abundance of 15 proteins, including myosin-3, related to muscle function, and aconitate hydratase, related to iron homeostasis. However, most E/S proteins (174 of 337 unique identities) were decreased, including those in the Ancylostoma-secreted protein (ASP) category, and metallopeptidases. Several identified proteins are established immune-modulators such as fatty acid-binding protein homologue, cystatin, and acetylcholinesterase. Enrichment analysis of InterPro functional categories showed down-regulation of Cysteine-rich secretory proteins, Antigen 5, and Pathogenesis-related 1 proteins (CAP), Astacin-like metallopeptidase, Glycoside hydrolase, and Transthyretin-like protein groups in STAT6 KO-adapted worms. Taken together, these data indicate that in an environment lacking Type 2 immunity, hookworms alter their secretome by reducing immune evasion proteins- and increasing locomotor- and feeding-associated proteins.

Long-Read Genome Assembly and Gene Model Annotations for the Rodent Malaria Parasite Plasmodium yoelii 17XNL

Malaria causes over 600 thousand fatalities each year, with most cases attributed to the human-infectious Plasmodium falciparum species. Many rodent-infectious Plasmodium species, like Plasmodium berghei and Plasmodium yoelii, have been used as model species that can expedite studies of this pathogen. P. yoelii is an especially good model for investigating the mosquito and liver stages of parasite development because key attributes closely resemble those of P. falciparum. Because of its importance, in 2002 the 17XNL strain of P. yoelii was the first rodent malaria parasite to be sequenced. While a breakthrough effort, the assembly consisted of >5000 contiguous sequences that adversely impacted the annotated gene models. While other rodent malaria parasite genomes have been sequenced and annotated since then, including the related P. yoelii 17X strain, the 17XNL strain has not. As a result, genomic data for 17X has become the de facto reference genome for the 17XNL strain while leaving open questions surrounding possible differences between the 17XNL and 17X genomes. In this work, we present a high-quality genome assembly for P. yoelii 17XNL using PacBio DNA sequencing. In addition, we use Nanopore and Illumina RNA sequencing of mixed blood stages to create complete gene models that include coding sequences, alternate isoforms, and UTR designations. A comparison of the 17X and this new 17XNL assembly revealed biologically meaningful differences between the strains due to the presence of coding sequence variants. Taken together, our work provides a new genomic framework for studies with this commonly used rodent malaria model species.

GUD-VE visualization tool for physicochemical properties of proteins

The physicochemical properties of primary sequences of proteins helps in determining both the structure and biological functions. The sequence analysis of the proteins and nucleic acids is most fundamental element of bioinformatics. Without these elements, it is impossible to gain insight deeper molecular and biochemical mechanisms. For this purpose, the computational methods like bioinformatics tools assist experts and novices alike in resolving issues relating to protein analysis. Similarly, this proposed work, for the graphical user interface (GUI) based prediction and visualization through the computations-based method done on Jupyter Notebook with tkinter package which allows the creation of a program on a local host platform and accessed by the programmer.•When it is queried with a protein sequence, it predicts physicochemical parameters of the peptides.•Users can choose to visualize the findings acquired either anonymously or on the user-specified email address and compare the biophysical properties of one protein with other using amino acids (AA) sequences. The aim of this paper is to meet the requirements of experimentalists, not just hardcore bioinformaticians related to biophysical properties prediction and comparison with other proteins. The code for it has been uploaded on GitHub (an online repository of codes) in private mode.

Genome-wide functional characterization of Canavalia rosea cysteine-rich trans-membrane module (CrCYSTM) genes to reveal their potential protective roles under extreme abiotic stress

Cysteine-rich transmembrane module (CYSTM) proteins constitute small molecular protein families and have been identified across eukaryotes, including yeast, humans, and several plant species. Plant CYSTMs play vital roles in growth regulation, development, phytohormone signal transduction, pathogen defense, environmental stress response, and even heavy metal binding and detoxification. Canavalia rosea (Sw.) DC is a perennial halophyte with great semi-arid and saline-alkali tolerance. In this study, the CrCYSTM family including 10 members were identified in the C. rosea genome, with the purpose of clarifying the possible roles of CrCYSTMs in C. rosea plants development and stress resistance. The phylogenetic relationships, exon-intron structure, domain structure, chromosomal localization, and putative cis-acting elements in promoter regions were predicted and analyzed. Transcriptome analysis combined with quantitative reverse transcription PCR showed that different CrCYSTM members exhibited varied expression patterns in different tissues and under different abiotic stress challenges. In addition, several CrCYSTMs were cloned and functionally characterized for their roles in abiotic stress tolerance with yeast expression system. Overall, these findings provide a foundation for functionally characterizing plant CYSTMs to unravel their possible roles in the adaptation of C. rosea to tropical coral reefs. Our results also lay the foundation for further research on the roles of plant CYSTM genes in abiotic stress signaling, especially for heavy metal detoxification.

Integrative omics analyses of the ligninolytic Rhodosporidium fluviale LM-2 disclose catabolic pathways for biobased chemical production

BACKGROUND

Lignin is an attractive alternative for producing biobased chemicals. It is the second major component of the plant cell wall and is an abundant natural source of aromatic compounds. Lignin degradation using microbial oxidative enzymes that depolymerize lignin and catabolize aromatic compounds into central metabolic intermediates is a promising strategy for lignin valorization. However, the intrinsic heterogeneity and recalcitrance of lignin severely hinder its biocatalytic conversion. In this context, examining microbial degradation systems can provide a fundamental understanding of the pathways and enzymes that are useful for lignin conversion into biotechnologically relevant compounds.

RESULTS

Lignin-degrading catabolism of a novel Rhodosporidium fluviale strain LM-2 was characterized using multi-omic strategies. This strain was previously isolated from a ligninolytic microbial consortium and presents a set of enzymes related to lignin depolymerization and aromatic compound catabolism. Furthermore, two catabolic routes for producing 4-vinyl guaiacol and vanillin were identified in R. fluviale LM-2.

CONCLUSIONS

The multi-omic analysis of R. fluviale LM-2, the first for this species, elucidated a repertoire of genes, transcripts, and secreted proteins involved in lignin degradation. This study expands the understanding of ligninolytic metabolism in a non-conventional yeast, which has the potential for future genetic manipulation. Moreover, this work unveiled critical pathways and enzymes that can be exported to other systems, including model organisms, for lignin valorization.

Long-Read Genome Assembly and Gene Model Annotations for the Rodent Malaria Parasite Plasmodium yoelii 17XNL

Malaria causes over 200 million infections and over 600 thousand fatalities each year, with most cases attributed to a human-infectious Plasmodium species, Plasmodium falciparum . Many rodent-infectious Plasmodium species, like Plasmodium berghei, Plasmodium chabaudi , and Plasmodium yoelii , have been used as genetically tractable model species that can expedite studies of this pathogen. In particular, P. yoelii is an especially good model for investigating the mosquito and liver stages of parasite development because key attributes closely resemble those of P. falciparum . Because of its importance to malaria research, in 2002 the 17XNL strain of P. yoelii was the first rodent malaria parasite to be sequenced. While sequencing and assembling this genome was a breakthrough effort, the final assembly consisted of >5000 contiguous sequences that impacted the creation of annotated gene models. While other important rodent malaria parasite genomes have been sequenced and annotated since then, including the related P. yoelii 17X strain, the 17XNL strain has not. As a result, genomic data for 17X has become the de facto reference genome for the 17XNL strain while leaving open questions surrounding possible differences between the 17XNL and 17X genomes. In this work, we present a high-quality genome assembly for P. yoelii 17XNL using HiFi PacBio long-read DNA sequencing. In addition, we use Nanopore long-read direct RNA-seq and Illumina short-read sequencing of mixed blood stages to create complete gene models that include not only coding sequences but also alternate transcript isoforms, and 5' and 3' UTR designations. A comparison of the 17X and this new 17XNL assembly revealed biologically meaningful differences between the strains due to the presence of coding sequence variants. Taken together, our work provides a new genomic and gene expression framework for studies with this commonly used rodent malaria model species.

Phylotranscriptomic Analyses of Mycoheterotrophic Monocots Show a Continuum of Convergent Evolutionary Changes in Expressed Nuclear Genes From Three Independent Nonphotosynthetic Lineages

Mycoheterotrophy is an alternative nutritional strategy whereby plants obtain sugars and other nutrients from soil fungi. Mycoheterotrophy and associated loss of photosynthesis have evolved repeatedly in plants, particularly in monocots. Although reductive evolution of plastomes in mycoheterotrophs is well documented, the dynamics of nuclear genome evolution remains largely unknown. Transcriptome datasets were generated from four mycoheterotrophs in three families (Orchidaceae, Burmanniaceae, Triuridaceae) and related green plants and used for phylogenomic analyses to resolve relationships among the mycoheterotrophs, their relatives, and representatives across the monocots. Phylogenetic trees based on 602 genes were mostly congruent with plastome phylogenies, except for an Asparagales + Liliales clade inferred in the nuclear trees. Reduction and loss of chlorophyll synthesis and photosynthetic gene expression and relaxation of purifying selection on retained genes were progressive, with greater loss in older nonphotosynthetic lineages. One hundred seventy-four of 1375 plant benchmark universally conserved orthologous genes were undetected in any mycoheterotroph transcriptome or the genome of the mycoheterotrophic orchid Gastrodia but were expressed in green relatives, providing evidence for massively convergent gene loss in nonphotosynthetic lineages. We designate this set of deleted or undetected genes Missing in Mycoheterotrophs (MIM). MIM genes encode not only mainly photosynthetic or plastid membrane proteins but also a diverse set of plastid processes, genes of unknown function, mitochondrial, and cellular processes. Transcription of a photosystem II gene (psb29) in all lineages implies a nonphotosynthetic function for this and other genes retained in mycoheterotrophs. Nonphotosynthetic plants enable novel insights into gene function as well as gene expression shifts, gene loss, and convergence in nuclear genomes.

A conserved Neurite Outgrowth and Guidance motif with biomimetic potential in neuronal Cell Adhesion Molecules

The discovery of conserved protein motifs can, in turn, unveil important regulatory signals, and when properly designed, synthetic peptides derived from such motifs can be used as biomimetics for biotechnological and therapeutic purposes. We report here that specific Ig-like repeats from the extracellular domains of neuronal Cell Adhesion Molecules share a highly conserved Neurite Outgrowth and Guidance (NOG) motif, which mediates homo- and heterophilic interactions crucial in neural development and repair. Synthetic peptides derived from the NOG motif of such proteins can boost neuritogenesis, and this potential is also retained by peptides with recombinant sequences, when fitting the NOG sequence pattern. The NOG motif discovery not only provides one more tile to the complex puzzle of neuritogenesis, but also opens the route to new neural regeneration strategies via a tunable biomimetic toolbox.

Screening Balanites aegyptiaca for inhibitors against putative drug targets in Microsporum gypseum - Subtractive proteome, docking and simulation approach

Microsporum gypseum is a keratinophilic fungi grouped under dermatophytes infecting skin, hair and nail portions in human and animals causing tinea corporis, tinea facei and tinea capitis. As both human and fungi are eukaryotes, the available drugs for treating dermatophytes produce some side effects due to drug interaction with human also. Apart from this, the gut microbiota has a very big role in the health of human which should not be affected by the drugs. Hence this study focused on finding a target which is unique and essential to M. gypseum and non-homologous to human and gut microbiota, non-homologous to human domain architecture, highly interacting with other proteins, sub-cellular localization of proteins and non-druggability analysis of the targets using subtractive proteomics approach which resulted with 3 novel drug targets from M. gypseum which were modeled using I-TASSER, refined by ModRefiner and validated by PROCHECK. Further these targets were docked with compounds identified through LC-MS of fractioned methanol extract of B. aegyptiaca fruit pulp using Glide module and the stability of the docked complex was analyzed by molecular dynamics simulation using Desmond module of Schrodinger. Cyanidin-3-O-rhamnoside had better interaction with all the targets and Taurocholic acid had better result with ECCP which suggests the multi-targeting potency of these two compounds against M. gypseum which has to be confirmed by in vitro and in vivo studies.

Genome analysis of a halophilic bacterium Halomonas malpeensis YU-PRIM-29T reveals its exopolysaccharide and pigment producing capabilities

Halomonas malpeensis strain YU-PRIM-29T is a yellow pigmented, exopolysaccharide (EPS) producing halophilic bacterium isolated from the coastal region. To understand the biosynthesis pathways involved in the EPS and pigment production, whole genome analysis was performed. The complete genome sequencing and the de novo assembly were carried out using Illumina sequencing and SPAdes genome assembler (ver 3.11.1) respectively followed by detailed genome annotation. The genome consists of 3,607,821 bp distributed in 18 contigs with 3337 protein coding genes and 53% of the annotated CDS are having putative functions. Gene annotation disclosed the presence of genes involved in ABC transporter-dependent pathway of EPS biosynthesis. As the ABC transporter-dependent pathway is also implicated in the capsular polysaccharide (CPS) biosynthesis, we employed extraction protocols for both EPS (from the culture supernatants) and CPS (from the cells) and found that the secreted polysaccharide i.e., EPS was predominant. The EPS showed good emulsifying activities against the petroleum hydrocarbons and its production was dependent on the carbon source supplied. The genome analysis also revealed genes involved in industrially important metabolites such as zeaxanthin pigment, ectoine and polyhydroxyalkanoate (PHA) biosynthesis. To confirm the genome data, we extracted these metabolites from the cultures and successfully identified them. The pigment extracted from the cells showed the distinct UV-Vis spectra having characteristic absorption peak of zeaxanthin (λmax 448 nm) with potent antioxidant activities. The ability of H. malpeensis strain YU-PRIM-29T to produce important biomolecules makes it an industrially important bacterium.

Proteome-Wide Analysis of Lysine 2-Hydroxyisobutyrylation in the Phytopathogenic Fungus Botrytis cinerea

Posttranslational modifications (PTMs) of the whole proteome have become a hot topic in the research field of epigenetics, and an increasing number of PTM types have been identified and shown to play significant roles in different cellular processes. Protein lysine 2-hydroxyisobutyrylation (K_hib) is a newly detected PTM, and the 2-hydroxyisobutyrylome has been identified in several species. Botrytis cinerea is recognized as one of the most destructive pathogens due to its broad host distribution and very large economic losses; thus the many aspects of its pathogenesis have been continuously studied. However, distribution and function of K_hib in this phytopathogenic fungus are not clear. In this study, a proteome-wide analysis of K_hib in B. cinerea was performed, and 5,398 K_hib sites on 1,181 proteins were identified. Bioinformatics analysis showed that the 2-hydroxyisobutyrylome in B. cinerea contains both conserved proteins and novel proteins when compared with K_hib proteins in other species. Functional classification, functional enrichment and protein interaction network analyses showed that K_hib proteins are widely distributed in cellular compartments and involved in diverse cellular processes. Significantly, 37 proteins involved in different aspects of regulating the pathogenicity of B. cinerea were detected as K_hib proteins. Our results provide a comprehensive view of the 2-hydroxyisobutyrylome and lay a foundation for further studying the regulatory mechanism of K_hib in both B. cinerea and other plant pathogens.

ToxoDB: Functional Genomics Resource for Toxoplasma and Related Organisms

ToxoDB is a free online resource that provides access to genomic and functional genomic data. All data is made available through an intuitive queryable interface that enables scientists to build in silico experiments and develop testable hypothesis. The resource contains 32 fully sequenced and annotated genomes, with genomic sequence from multiple strains available for variant detection and copy number variation analysis. In addition to genomic sequence data, ToxoDB contains numerous functional genomic datasets including microarray, RNAseq, proteomics, ChIP-seq, and phenotypic data. In addition, results from a number of whole-genome analyses are incorporated including mapping to orthology clusters which allows users to leverage phylogenetic relationships in their analyses. Integration of primary data is made possible through a private galaxy interface and custom export tools that allow users to interrogate their own results in the context of all other data in the database.

The Class III Peroxidase (POD) Gene Family in Cassava: Identification, Phylogeny, Duplication, and Expression

The class III peroxidase (POD) enzymes participate in plant development, hormone signaling, and stress responses. However, little is known about the POD family in cassava. Here, we identified 91 cassava POD genes (MePODs) and classified them into six subgroups using phylogenetic analysis. Conserved motif analysis demonstrated that all MePOD proteins have typical peroxidase domains, and gene structure analysis showed that MePOD genes have between one and nine exons. Duplication pattern analysis suggests that tandem duplication has played a role in MePOD gene expansion. Comprehensive transcriptomic analysis revealed that MePOD genes in cassava are involved in the drought response and postharvest physiological deterioration. Several MePODs underwent transcriptional changes after various stresses and related signaling treatments were applied. In sum, we characterized the POD family in cassava and uncovered the transcriptional control of POD genes in response to various stresses and postharvest physiological deterioration conditions. These results can be used to identify potential target genes for improving the stress tolerance of cassava crops.

A highly expressed intestinal cysteine protease of Ancylostoma ceylanicum protects vaccinated hamsters from hookworm infection

BACKGROUND

Human hookworms (Necator americanus, Ancylostoma duodenale, and Ancylostoma ceylanicum) are intestinal blood-feeding parasites that infect ~500 million people worldwide and are among the leading causes of iron-deficiency anemia in the developing world. Drugs are useful against hookworm infections, but hookworms rapidly reinfect people, and the parasites can develop drug resistance. Therefore, having a hookworm vaccine would be of tremendous benefit.

METHODOLOGY/PRINCIPAL FINDINGS

We investigated the vaccine efficacy in outbred Syrian hamsters of three A. ceylanicum hookworm antigen candidates from two classes of proteins previously identified as promising vaccine candidates. These include two intestinally-enriched, putatively secreted cathepsin B cysteine proteases (AceyCP1, AceyCPL) and one small Kunitz-type protease inhibitor (AceySKPI3). Recombinant proteins were produced in Pichia pastoris, and adsorbed to Alhydrogel. Recombinant AceyCPL (rAceyCPL)/Alhydrogel and rAceySKPI3/Alhydrogel induced high serum immunoglobulin G (IgG) titers in 8/8 vaccinates, but were not protective. rAceyCP1/Alhydrogel induced intermediate serum IgG titers in ~60% of vaccinates in two different trials. rAceyCP1 serum IgG responders had highly significantly decreased hookworm burdens, fecal egg counts and clinical pathology compared to Alhydrogel controls and nonresponders. Protection was highly correlated with rAceyCP1 serum IgG titer. Antisera from rAceyCP1 serum IgG responders, but not nonresponders or rAceyCPL/Alhydrogel vaccinates, significantly reduced adult A. ceylanicum motility in vitro. Furthermore, rAceyCP1 serum IgG responders had canonical Th2-specific recall responses (IL4, IL5, IL13) in splenocytes stimulated ex vivo.

CONCLUSIONS/SIGNIFICANCE

These findings indicate that rAceyCP1 is a promising vaccine candidate and validates a genomic/transcriptomic approach to human hookworm vaccine discovery.

Expression of putative circadian clock components in the arbuscular mycorrhizal fungus Rhizoglomus irregulare

Arbuscular mycorrhizal fungi (AMF) are obligatory plant symbionts that live underground, so few studies have examined their response to light. Responses to blue light by other fungi can be mediated by White Collar-1 (WC-1) and WC-2 proteins. These wc genes, together with the frequency gene (frq), also form part of the endogenous circadian clock. The clock mechanism has never been studied in AMF, although circadian growth of their hyphae in the field has been reported. Using both genomic and transcriptomic data, we have found homologs of wc-1, wc-2, and frq and related circadian clock genes in the arbuscular mycorrhizal fungus Rhizoglomus irregulare (synonym Rhizophagus irregularis). Gene expression of wc-1, wc-2, and frq was analyzed using RT-qPCR on RNA extracted from germinating spores and from fungal material cultivated in vitro with transformed carrot roots. We found that all three core clock genes were expressed in both pre- and post-mycorrhizal stages of R. irregulare growth. Similar to the model fungus Neurospora crassa, the core circadian oscillator gene frq was induced by brief light stimulation. The presence of circadian clock and output genes in R. irregulare opens the door to the study of circadian clocks in the fungal partner of plant-AMF symbiosis. Our finding also provides new insight into the evolution of the circadian frq gene in fungi.

Distinct mutations with different inheritance mode caused similar retinal dystrophies in one family: a demonstration of the importance of genetic annotations in complicated pedigrees

Background

Retinitis pigmentosa (RP) is the most common form of inherited retinal dystrophy presenting remarkable genetic heterogeneity. Genetic annotations would help with better clinical assessments and benefit gene therapy, and therefore should be recommended for RP patients. This report reveals the disease causing mutations in two RP pedigrees with confusing inheritance patterns using whole exome sequencing (WES).

Methods

Twenty-five participants including eight patients from two families were recruited and received comprehensive ophthalmic evaluations. WES was applied for mutation identification. Bioinformatics annotations, intrafamilial co-segregation tests, and in silico analyses were subsequently conducted for mutation verification.

Results

All patients were clinically diagnosed with RP. The first family included two siblings born to parents with consanguineous marriage; however, no potential pathogenic variant was found shared by both patients. Further analysis revealed that the female patient carried a recurrent homozygous C8ORF37 p.W185*, while the male patient had hemizygous OFD1 p.T120A. The second family was found to segregate mutations in two genes, TULP1 and RP1. Two patients born to consanguineous marriage carried homozygous TULP1 p.R419W, while a recurrent heterozygous RP1 p.L762Yfs*17 was found in another four patients presenting an autosomal dominant inheritance pattern. Crystal structural analysis further indicated that the substitution from arginine to tryptophan at the highly conserved residue 419 of TULP1 could lead to the elimination of two hydrogen bonds between residue 419 and residues V488 and S534. All four genes, including C8ORF37, OFD1, TULP1 and RP1, have been previously implicated in RP etiology.

Conclusions

Our study demonstrates the coexistence of diverse inheritance modes and mutations affecting distinct disease causing genes in two RP families with consanguineous marriage. Our data provide novel insights into assessments of complicated pedigrees, reinforce the genetic complexity of RP, and highlight the need for extensive molecular evaluations in such challenging families with diverse inheritance modes and mutations.

Electronic supplementary material

The online version of this article (10.1186/s12967-018-1522-7) contains supplementary material, which is available to authorized users.

New reference genome sequences of hot pepper reveal the massive evolution of plant disease-resistance genes by retroduplication

Background

Transposable elements are major evolutionary forces which can cause new genome structure and species diversification. The role of transposable elements in the expansion of nucleotide-binding and leucine-rich-repeat proteins (NLRs), the major disease-resistance gene families, has been unexplored in plants.

Results

We report two high-quality de novo genomes (Capsicum baccatum and C. chinense) and an improved reference genome (C. annuum) for peppers. Dynamic genome rearrangements involving translocations among chromosomes 3, 5, and 9 were detected in comparison between C. baccatum and the two other peppers. The amplification of athila LTR-retrotransposons, members of the gypsy superfamily, led to genome expansion in C. baccatum. In-depth genome-wide comparison of genes and repeats unveiled that the copy numbers of NLRs were greatly increased by LTR-retrotransposon-mediated retroduplication. Moreover, retroduplicated NLRs are abundant across the angiosperms and, in most cases, are lineage-specific.

Conclusions

Our study reveals that retroduplication has played key roles for the massive emergence of NLR genes including functional disease-resistance genes in pepper plants.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1341-9) contains supplementary material, which is available to authorized users.

New reference genome sequences of hot pepper reveal the massive evolution of plant disease-resistance genes by retroduplication

BACKGROUND

Transposable elements are major evolutionary forces which can cause new genome structure and species diversification. The role of transposable elements in the expansion of nucleotide-binding and leucine-rich-repeat proteins (NLRs), the major disease-resistance gene families, has been unexplored in plants.

RESULTS

We report two high-quality de novo genomes (Capsicum baccatum and C. chinense) and an improved reference genome (C. annuum) for peppers. Dynamic genome rearrangements involving translocations among chromosomes 3, 5, and 9 were detected in comparison between C. baccatum and the two other peppers. The amplification of athila LTR-retrotransposons, members of the gypsy superfamily, led to genome expansion in C. baccatum. In-depth genome-wide comparison of genes and repeats unveiled that the copy numbers of NLRs were greatly increased by LTR-retrotransposon-mediated retroduplication. Moreover, retroduplicated NLRs are abundant across the angiosperms and, in most cases, are lineage-specific.

CONCLUSIONS

Our study reveals that retroduplication has played key roles for the massive emergence of NLR genes including functional disease-resistance genes in pepper plants.

Mechanistic underpinnings of dehydration stress in the American dog tick revealed through RNA-Seq and metabolomics

Ticks are obligate blood feeders but spend the majority of their lifetime off-host where they must contend with a multitude of environmental stresses. Survival under desiccating conditions is a determinant for habitats where ticks can become established, and water-balance characteristics of ticks have been extensively studied. However, little is known about the underlying aspects associated with dehydration stress in ticks. In this study, we examined the response of male American dog ticks, Dermacentor variabilis, to dehydration using a combined transcriptomics and metabolomics approach. During dehydration, 497 genes were differentially expressed, including an up-regulation of stress-response and protein-catabolism genes and concurrent down-regulation of several energetically expensive biological processes. Accumulation of several metabolites, including specific amino acids, glycerol and gamma aminobutyric acid (GABA), and transcript shifts in the associated pathways for generating these metabolites indicated congruence between changes in the metabolome and gene expression. Ticks treated with exogenous glycerol and GABA demonstrated altered water-balance characteristics; specifically, increased water absorption at high relative humidity. Finally, we observed changes in locomotor activity in response to dehydration, but this change was not influenced by the accumulation of GABA. Overall, the responses to dehydration by these ticks were similar to those observed in other dehydration-tolerant arthropods, but several molecular and behavioral responses are distinct from those associated with other taxa.

Probing the influence of carboxyalkyl groups on the molecular flexibility and the charge density of apigenin derivatives

Apigenin is an important flavonoids due to its antidiabetic bioactivity. It was reported experimentally that the 7-substituent derivative of apigenin has higher biological activity than 4'- and 5-substituted derivatives while introducing sole carboxyalkyl group -(CH₂)₇COOH into the parent structure. Molecular docking studies indicated that the other two derivatives have lower binding affinities than the 7-substituent derivative (-7.52 kcal mol^-1), which is considered to be a better inhibitor than the parent molecule. Almost all of the carbon atoms and oxygen atoms are coplaner for all three molecules in solution phase, however, all carboxyalkyl groups bend inside into the parent molecules in the active site, and the jagged geometries of the carbon chains are destroyed correspondingly. In addition, most of the electron densities of the chemical bonds for all molecules are decreased, especially the 7-substituent derivative. In contrast, most of the Laplacian values for three molecules are increased in the active site, which suggests that the charge densities at the bond critical point (bcp) are much more depleted than the solution phase. Dipole moments of derivatives are all increased in the active site, suggesting strong intermolecular interactions. After interacting with the S. cerevisiae α-glucosidase, only the 7-substituent derivative has the lowest energy gap ΔE _HOMO-LUMO, which indicates the lowest stability and the highest inhibition activity. Graphical abstract Probing the influence of carboxyalkyl groups on the molecular flexibility and the charge density of apigenin derivatives.

Safety assessment of Staphylococcus phages of the family Myoviridae based on complete genome sequences

Staphylococcus phages of the Myoviridae family have a wide host range and potential applications in phage therapy. In this report, safety assessments of these phages were conducted based on their complete genome sequences. The complete genomes of Staphylococcus phages of the Myoviridae family were analyzed, and the Open Reading Frame (ORFs) were compared with a pool of virulence and antibiotic resistance genes using the BLAST algorithm. In addition, the lifestyle of the phages (virulent or temperate) was also confirmed using PHACTS. The results showed that all phages were lytic and did not contain resistance or virulence genes based on bioinformatic analyses, excluding the possibility that they could be vectors for the dissemination of these undesirable genes. These findings suggest that the phages are safe at the genome level. The SceD-like transglycosylase, which is a biomarker for vancomycin-intermediate strains, was widely distributed in the phage genomes. Approximately 70% of the ORFs encoded in the phage genomes have unknown functions; therefore, their roles in the antibiotic resistance and virulence of Staphylococcus aureus are still unknown and require consideration before use in phage therapy.

Maize pan-transcriptome provides novel insights into genome complexity and quantitative trait variation

Gene expression variation largely contributes to phenotypic diversity and constructing pan-transcriptome is considered necessary for species with complex genomes. However, the regulation mechanisms and functional consequences of pan-transcriptome is unexplored systematically. By analyzing RNA-seq data from 368 maize diverse inbred lines, we identified almost one-third nuclear genes under expression presence and absence variation, which tend to play regulatory roles and are likely regulated by distant eQTLs. The ePAV was directly used as “genotype” to perform GWAS for 15 agronomic phenotypes and 526 metabolic traits to efficiently explore the associations between transcriptomic and phenomic variations. Through a modified assembly strategy, 2,355 high-confidence novel sequences with total 1.9 Mb lengths were found absent within reference genome. Ten randomly selected novel sequences were fully validated with genomic PCR, including another two NBS_LRR candidates potentially affect flavonoids and disease-resistance. A simulation analysis suggested that the pan-transcriptome of the maize whole kernel is approaching a maximum value of 63,000 genes, and through developing two test-cross populations and surveying several most important yield traits, the dispensable genes were shown to contribute to heterosis. Novel perspectives and resources to discover maize quantitative trait variations were provided to better understand the kernel regulation networks and to enhance maize breeding.

DYNC2LI1 mutations broaden the clinical spectrum of dynein-2 defects

Skeletal ciliopathies are a heterogeneous group of autosomal recessive osteochondrodysplasias caused by defects in formation, maintenance and function of the primary cilium. Mutations in the underlying genes affect the molecular motors, intraflagellar transport complexes (IFT), or the basal body. The more severe phenotypes are caused by defects of genes of the dynein-2 complex, where mutations in DYNC2H1, WDR34 and WDR60 have been identified. In a patient with a Jeune-like phenotype we performed exome sequencing and identified compound heterozygous missense and nonsense mutations in DYNC2LI1 segregating with the phenotype. DYNC2LI1 is ubiquitously expressed and interacts with DYNC2H1 to form the dynein-2 complex important for retrograde IFT. Using DYNC2LI1 siRNA knockdown in fibroblasts we identified a significantly reduced cilia length proposed to affect cilia function. In addition, depletion of DYNC2LI1 induced altered cilia morphology with broadened ciliary tips and accumulation of IFT-B complex proteins in accordance with retrograde IFT defects. Our results expand the clinical spectrum of ciliopathies caused by defects of the dynein-2 complex.

Enrichment or depletion? The impact of stool pretreatment on metaproteomic characterization of the human gut microbiota

To date, most metaproteomic studies of the gut microbiota employ stool sample pretreatment methods to enrich for microbial components. However, a specific investigation aimed at assessing if, how, and to what extent this may impact on the final taxonomic and functional results is still lacking. Here, stool replicates were either pretreated by differential centrifugation (DC) or not centrifuged. Protein extracts were then processed by filter-aided sample preparation, single-run LC, and high-resolution MS, and the metaproteomic data were compared by spectral counting. DC led to a higher number of identifications, a significantly richer microbial diversity, as well as to reduced information on the nonmicrobial components (host and food) when compared to not centrifuged. Nevertheless, dramatic differences in the relative abundance of several gut microbial taxa were also observed, including a significant change in the Firmicutes/Bacteroidetes ratio. Furthermore, some important microbial functional categories, including cell surface enzymes, membrane-associated proteins, extracellular proteins, and flagella, were significantly reduced after DC. In conclusion, this work underlines that a critical evaluation is needed when selecting the appropriate stool sample processing protocol in the context of a metaproteomic study, depending on the specific target to which the research is aimed. All MS data have been deposited in the ProteomeXchange with identifier PXD001573 (http://proteomecentral.proteomexchange.org/dataset/PXD001573).

The genome and transcriptome of the zoonotic hookworm Ancylostoma ceylanicum identify infection-specific gene families

Hookworms infect over 400 million people, stunting and impoverishing them. Sequencing hookworm genomes and finding which genes they express during infection should help in devising new drugs or vaccines against hookworms. Unlike other hookworms, Ancylostoma ceylanicum infects both humans and other mammals, providing a laboratory model for hookworm disease. We determined an A. ceylanicum genome sequence of 313 Mb, with transcriptomic data throughout infection showing expression of 30,738 genes. Approximately 900 genes were upregulated during early infection in vivo, including ASPRs, a cryptic subfamily of activation-associated secreted proteins (ASPs). Genes downregulated during early infection included ion channels and G protein-coupled receptors; this downregulation was observed in both parasitic and free-living nematodes. Later, at the onset of heavy blood feeding, C-lectin genes were upregulated along with genes for secreted clade V proteins (SCVPs), encoding a previously undescribed protein family. These findings provide new drug and vaccine targets and should help elucidate hookworm pathogenesis.

A bioinformatic survey of distribution, conservation, and probable functions of LuxR solo regulators in bacteria

LuxR solo transcriptional regulators contain both an autoinducer binding domain (ABD; N-terminal) and a DNA binding Helix-Turn-Helix domain (HTH; C-terminal), but are not associated with a cognate N-acyl homoserine lactone (AHL) synthase coding gene in the same genome. Although a few LuxR solos have been characterized, their distributions as well as their role in bacterial signal perception and other processes are poorly understood. In this study we have carried out a systematic survey of distribution of all ABD containing LuxR transcriptional regulators (QS domain LuxRs) available in the InterPro database (IPR005143), and identified those lacking a cognate AHL synthase. These LuxR solos were then analyzed regarding their taxonomical distribution, predicted functions of neighboring genes and the presence of complete AHL-QS systems in the genomes that carry them. Our analyses reveal the presence of one or multiple predicted LuxR solos in many proteobacterial genomes carrying QS domain LuxRs, some of them harboring genes for one or more AHL-QS circuits. The presence of LuxR solos in bacteria occupying diverse environments suggests potential ecological functions for these proteins beyond AHL and interkingdom signaling. Based on gene context and the conservation levels of invariant amino acids of ABD, we have classified LuxR solos into functionally meaningful groups or putative orthologs. Surprisingly, putative LuxR solos were also found in a few non-proteobacterial genomes which are not known to carry AHL-QS systems. Multiple predicted LuxR solos in the same genome appeared to have different levels of conservation of invariant amino acid residues of ABD questioning their binding to AHLs. In summary, this study provides a detailed overview of distribution of LuxR solos and their probable roles in bacteria with genome sequence information.

Comparative transcriptome analyses reveal core parasitism genes and suggest gene duplication and repurposing as sources of structural novelty

The origin of novel traits is recognized as an important process underlying many major evolutionary radiations. We studied the genetic basis for the evolution of haustoria, the novel feeding organs of parasitic flowering plants, using comparative transcriptome sequencing in three species of Orobanchaceae. Around 180 genes are upregulated during haustorial development following host attachment in at least two species, and these are enriched in proteases, cell wall modifying enzymes, and extracellular secretion proteins. Additionally, about 100 shared genes are upregulated in response to haustorium inducing factors prior to host attachment. Collectively, we refer to these newly identified genes as putative “parasitism genes.” Most of these parasitism genes are derived from gene duplications in a common ancestor of Orobanchaceae and Mimulus guttatus, a related nonparasitic plant. Additionally, the signature of relaxed purifying selection and/or adaptive evolution at specific sites was detected in many haustorial genes, and may play an important role in parasite evolution. Comparative analysis of gene expression patterns in parasitic and nonparasitic angiosperms suggests that parasitism genes are derived primarily from root and floral tissues, but with some genes co-opted from other tissues. Gene duplication, often taking place in a nonparasitic ancestor of Orobanchaceae, followed by regulatory neofunctionalization, was an important process in the origin of parasitic haustoria.

Comparative analyses of nonpathogenic, opportunistic, and totally pathogenic mycobacteria reveal genomic and biochemical variabilities and highlight the survival attributes of Mycobacterium tuberculosis

Mycobacterial evolution involves various processes, such as genome reduction, gene cooption, and critical gene acquisition. Our comparative genome size analysis of 44 mycobacterial genomes revealed that the nonpathogenic (NP) genomes were bigger than those of opportunistic (OP) or totally pathogenic (TP) mycobacteria, with the TP genomes being smaller yet variable in size—their genomic plasticity reflected their ability to evolve and survive under various environmental conditions. From the 44 mycobacterial species, 13 species, representing TP, OP, and NP, were selected for genomic-relatedness analyses. Analysis of homologous protein-coding genes shared between Mycobacterium indicus pranii (NP), Mycobacterium intracellulare ATCC 13950 (OP), and Mycobacterium tuberculosis H37Rv (TP) revealed that 4,995 (i.e., ~95%) M. indicaus pranii proteins have homology with M. intracellulare, whereas the homologies among M. indicus pranii, M. intracellulare ATCC 13950, and M. tuberculosis H37Rv were significantly lower. A total of 4,153 (~79%) M. indicus pranii proteins and 4,093 (~79%) M. intracellulare ATCC 13950 proteins exhibited homology with the M. tuberculosis H37Rv proteome, while 3,301 (~82%) and 3,295 (~82%) M. tuberculosis H37Rv proteins showed homology with M. indicus pranii and M. intracellulare ATCC 13950 proteomes, respectively. Comparative metabolic pathway analyses of TP/OP/NP mycobacteria showed enzymatic plasticity between M. indicus pranii (NP) and M. intracellulare ATCC 13950 (OP), Mycobacterium avium 104 (OP), and M. tuberculosis H37Rv (TP). Mycobacterium tuberculosis seems to have acquired novel alternate pathways with possible roles in metabolism, host-pathogen interactions, virulence, and intracellular survival, and by implication some of these could be potential drug targets.

The complete sequence analysis of Mycobacterium indicus pranii, a novel species of Mycobacterium shown earlier to have strong immunomodulatory properties and currently in use for the treatment of leprosy, places it evolutionarily at the point of transition to pathogenicity. With the purpose of establishing the importance of M. indicus pranii in providing insight into the virulence mechanism of tuberculous and nontuberculous mycobacteria, we carried out comparative genomic and proteomic analyses of 44 mycobacterial species representing nonpathogenic (NP), opportunistic (OP), and totally pathogenic (TP) mycobacteria. Our results clearly placed M. indicus pranii as an ancestor of the M. avium complex. Analyses of comparative metabolic pathways between M. indicus pranii (NP), M. tuberculosis (TP), and M. intracellulare (OP) pointed to the presence of novel alternative pathways in M. tuberculosis with implications for pathogenesis and survival in the human host and identification of new drug targets.

AnimalTFDB 2.0: a resource for expression, prediction and functional study of animal transcription factors

Transcription factors (TFs) are key regulators for gene expression. Here we updated the animal TF database AnimalTFDB to version 2.0 (http://bioinfo.life.hust.edu.cn/AnimalTFDB/). Using the improved prediction pipeline, we identified 72 336 TF genes, 21 053 transcription co-factor genes and 6502 chromatin remodeling factor genes from 65 species covering main animal lineages. Besides the abundant annotations (basic information, gene model, protein functional domain, gene ontology, pathway, protein interaction, ortholog and paralog, etc.) in the previous version, we made several new features and functions in the updated version. These new features are: (i) gene expression from RNA-Seq for nine model species, (ii) gene phenotype information, (iii) multiple sequence alignment of TF DNA-binding domains, and the weblogo and phylogenetic tree based on the alignment, (iv) a TF prediction server to identify new TFs from input sequences and (v) a BLAST server to search against TFs in AnimalTFDB. A new nice web interface was designed for AnimalTFDB 2.0 allowing users to browse and search all data in the database. We aim to maintain the AnimalTFDB as a solid resource for TF identification and studies of transcription regulation and comparative genomics.

Basis for substrate recognition and distinction by matrix metalloproteinases

Genomic sequencing and structural genomics produced a vast amount of sequence and structural data, creating an opportunity for structure-function analysis in silico [Radivojac P, et al. (2013) Nat Methods 10(3):221-227]. Unfortunately, only a few large experimental datasets exist to serve as benchmarks for function-related predictions. Furthermore, currently there are no reliable means to predict the extent of functional similarity among proteins. Here, we quantify structure-function relationships among three phylogenetic branches of the matrix metalloproteinase (MMP) family by comparing their cleavage efficiencies toward an extended set of phage peptide substrates that were selected from ∼ 64 million peptide sequences (i.e., a large unbiased representation of substrate space). The observed second-order rate constants [k(obs)] across the substrate space provide a distance measure of functional similarity among the MMPs. These functional distances directly correlate with MMP phylogenetic distance. There is also a remarkable and near-perfect correlation between the MMP substrate preference and sequence identity of 50-57 discontinuous residues surrounding the catalytic groove. We conclude that these residues represent the specificity-determining positions (SDPs) that allowed for the expansion of MMP proteolytic function during evolution. A transmutation of only a few selected SDPs proximal to the bound substrate peptide, and contributing the most to selectivity among the MMPs, is sufficient to enact a global change in the substrate preference of one MMP to that of another, indicating the potential for the rational and focused redesign of cleavage specificity in MMPs.

Xylem sap proteomics

Proteomic analysis of xylem sap has recently become a major field of interest to understand several biological questions related to plant development and responses to environmental clues. The xylem sap appears as a dynamic fluid undergoing changes in its proteome upon abiotic and biotic stresses. Unlike cell compartments which are amenable to purification in sufficient amount prior to proteomic analysis, the xylem sap has to be collected in particular conditions to avoid contamination by intracellular proteins and to obtain enough material. A model plant like Arabidopsis thaliana is not suitable for such an analysis because efficient harvesting of xylem sap is difficult. The analysis of the xylem sap proteome also requires specific procedures to concentrate proteins and to focus on proteins predicted to be secreted. Indeed, xylem sap proteins appear to be synthesized and secreted in the root stele or to originate from dying differentiated xylem cells. This chapter describes protocols to collect xylem sap from Brassica species and to prepare total and N-glycoprotein extracts for identification of proteins by mass spectrometry analyses and bioinformatics.

Functional genomic analyses of two morphologically distinct classes of Drosophila sensory neurons: post-mitotic roles of transcription factors in dendritic patterning

Background

Neurons are one of the most structurally and functionally diverse cell types found in nature, owing in large part to their unique class specific dendritic architectures. Dendrites, being highly specialized in receiving and processing neuronal signals, play a key role in the formation of functional neural circuits. Hence, in order to understand the emergence and assembly of a complex nervous system, it is critical to understand the molecular mechanisms that direct class specific dendritogenesis.

Methodology/Principal Findings

We have used the Drosophila dendritic arborization (da) neurons to gain systems-level insight into dendritogenesis by a comparative study of the morphologically distinct Class-I (C-I) and Class-IV (C-IV) da neurons. We have used a combination of cell-type specific transcriptional expression profiling coupled to a targeted and systematic in vivo RNAi functional validation screen. Our comparative transcriptomic analyses have revealed a large number of differentially enriched/depleted gene-sets between C-I and C-IV neurons, including a broad range of molecular factors and biological processes such as proteolytic and metabolic pathways. Further, using this data, we have identified and validated the role of 37 transcription factors in regulating class specific dendrite development using in vivo class-specific RNAi knockdowns followed by rigorous and quantitative neurometric analysis.

Conclusions/Significance

This study reports the first global gene-expression profiles from purified Drosophila C-I and C-IV da neurons. We also report the first large-scale semi-automated reconstruction of over 4,900 da neurons, which were used to quantitatively validate the RNAi screen phenotypes. Overall, these analyses shed global and unbiased novel insights into the molecular differences that underlie the morphological diversity of distinct neuronal cell-types. Furthermore, our class-specific gene expression datasets should prove a valuable community resource in guiding further investigations designed to explore the molecular mechanisms underlying class specific neuronal patterning.

Large-scale gene function analysis with the PANTHER classification system

The PANTHER (protein annotation through evolutionary relationship) classification system (http://www.pantherdb.org/) is a comprehensive system that combines gene function, ontology, pathways and statistical analysis tools that enable biologists to analyze large-scale, genome-wide data from sequencing, proteomics or gene expression experiments. The system is built with 82 complete genomes organized into gene families and subfamilies, and their evolutionary relationships are captured in phylogenetic trees, multiple sequence alignments and statistical models (hidden Markov models or HMMs). Genes are classified according to their function in several different ways: families and subfamilies are annotated with ontology terms (Gene Ontology (GO) and PANTHER protein class), and sequences are assigned to PANTHER pathways. The PANTHER website includes a suite of tools that enable users to browse and query gene functions, and to analyze large-scale experimental data with a number of statistical tests. It is widely used by bench scientists, bioinformaticians, computer scientists and systems biologists. In the 2013 release of PANTHER (v.8.0), in addition to an update of the data content, we redesigned the website interface to improve both user experience and the system's analytical capability. This protocol provides a detailed description of how to analyze genome-wide experimental data with the PANTHER classification system.

Comparative transcriptomics of infectious spores from the fungal pathogen Histoplasma capsulatum reveals a core set of transcripts that specify infectious and pathogenic states

Histoplasma capsulatum is a fungal pathogen that infects both healthy and immunocompromised hosts. In regions where it is endemic, H. capsulatum grows in the soil and causes respiratory and systemic disease when inhaled by humans. An interesting aspect of H. capsulatum biology is that it adopts specialized developmental programs in response to its environment. In the soil, it grows as filamentous chains of cells (mycelia) that produce asexual spores (conidia). When the soil is disrupted, conidia aerosolize and are inhaled by mammalian hosts. Inside a host, conidia germinate into yeast-form cells that colonize immune cells and cause disease. Despite the ability of conidia to initiate infection and disease, they have not been explored on a molecular level. We developed methods to purify H. capsulatum conidia, and we show here that these cells germinate into filaments at room temperature and into yeast-form cells at 37°C. Conidia internalized by macrophages germinate into the yeast form and proliferate within macrophages, ultimately lysing the host cells. Similarly, infection of mice with purified conidia is sufficient to establish infection and yield viable yeast-form cells in vivo. To characterize conidia on a molecular level, we performed whole-genome expression profiling of conidia, yeast, and mycelia from two highly divergent H. capsulatum strains. In parallel, we used homology and protein domain analysis to manually annotate the predicted genes of both strains. Analyses of the resultant data defined sets of transcripts that reflect the unique molecular states of H. capsulatum conidia, yeast, and mycelia.

Evolutionary, structural and functional interplay of the IκB family members

A primary level of control for nuclear factor kappa B (NF-κB) is effected through its interactions with the inhibitor protein, inhibitor of kappa B (IκB). Several lines of evidence confirm the existence of multiple forms of IκB that appear to regulate NF-κB by distinct mechanisms. Therefore, we performed a comprehensive bioinformatics analysis to understand the evolutionary history and intrinsic functional diversity of IκB family members. Phylogenetic relationships were constructed to trace the evolution of the IκB family genes. Our phylogenetic analysis revealed 10 IκB subfamily members that clustered into 5 major clades. Since the ankyrin (ANK) domain appears to be more ancient than the Rel homology domain (RHD), our phylogenetic analysis suggests that some undefined ancestral set of ANK repeats acquired an RHD before any duplication and was later duplicated and then diverged into the different IκB subfamilies. Functional analysis identified several functionally divergent sites in the ANK repeat domains (ARDs) and revealed that this region has undergone strong purifying selection, suggesting its functional importance in IκB genes. Structural analysis showed that the major variations in the number of ANK repeats and high conformational changes in the finger loop ARD region contribute to the differing binding partner specificities, thereby leading to distinct IκB functions. In summary, our study has provided useful information about the phylogeny and structural and functional divergence of the IκB family. Additionally, we identified a number of amino acid sites that contribute to the predicted functional divergence of these proteins.

Molecular evolution and structural features of IRAK family members

The interleukin-1 receptor-associated kinase (IRAK) family comprises critical signaling mediators of the TLR/IL-1R signaling pathways. IRAKs are Ser/Thr kinases. There are 4 members in the vertebrate genome (IRAK1, IRAK2, IRAKM, and IRAK4) and an IRAK homolog, Pelle, in insects. IRAK family members are highly conserved in vertebrates, but the evolutionary relationship between IRAKs in vertebrates and insects is not clear. To investigate the evolutionary history and functional divergence of IRAK members, we performed extensive bioinformatics analysis. The phylogenetic relationship between IRAK sequences suggests that gene duplication events occurred in the evolutionary lineage, leading to early vertebrates. A comparative phylogenetic analysis with insect homologs of IRAKs suggests that the Tube protein is a homolog of IRAK4, unlike the anticipated protein, Pelle. Furthermore, the analysis supports that an IRAK4-like kinase is an ancestral protein in the metazoan lineage of the IRAK family. Through functional analysis, several potentially diverged sites were identified in the common death domain and kinase domain. These sites have been constrained during evolution by strong purifying selection, suggesting their functional importance within IRAKs. In summary, our study highlighted the molecular evolution of the IRAK family, predicted the amino acids that contributed to functional divergence, and identified structural variations among the IRAK paralogs that may provide a starting point for further experimental investigations.

Genomics and localization of the Arabidopsis DHHC-cysteine-rich domain S-acyltransferase protein family

Protein lipid modification of cysteine residues, referred to as S-palmitoylation or S-acylation, is an important secondary and reversible modification that regulates membrane association, trafficking, and function of target proteins. This enzymatic reaction is mediated by protein S-acyl transferases (PATs). Here, the phylogeny, genomic organization, protein topology, expression, and localization pattern of the 24 PAT family members from Arabidopsis (Arabidopsis thaliana) is described. Most PATs are expressed at ubiquitous levels and tissues throughout the development, while few genes are expressed especially during flower development preferentially in pollen and stamen. The proteins display large sequence and structural variations but exhibit a common protein topology that is preserved in PATs from various organisms. Arabidopsis PAT proteins display a complex targeting pattern and were detected at the endoplasmic reticulum, Golgi, endosomal compartments, and the vacuolar membrane. However, most proteins were targeted to the plasma membrane. This large concentration of plant PAT activity to the plasma membrane suggests that the plant cellular S-acylation machinery is functionally different compared with that of yeast (Saccharomyces cerevisiae) and mammalians.

Genome-wide identification and analysis of late embryogenesis abundant (LEA) genes in Prunus mume

Late embryogenesis abundant (LEA) proteins play important roles in plant desiccation tolerance. In this study, 30 LEA genes were identified from Chinese plum (Prunus mume) through genome-wide analysis. The PmLEA genes are distributed on all Chinese plum chromosomes except chromosome 3. Twelve (40 %) and five PmLEA genes are arranged in tandem and segmental duplications, respectively. The PmLEA genes could be divided into eight groups (LEA_1, LEA_2, LEA_3, LEA_4, LEA_5, PvLEA18, dehydrin and seed maturation protein). Ten gene conversion events were observed and most of them (70 %) were identified in dehydrin group. Most PmLEA genes were highly expressed in flower (22/30) and up-regulated by ABA treatment (19/30).

Genome-wide identification and analysis of late embryogenesis abundant (LEA) genes in Prunus mume

Late embryogenesis abundant (LEA) proteins play important roles in plant desiccation tolerance. In this study, 30 LEA genes were identified from Chinese plum (Prunus mume) through genome-wide analysis. The PmLEA genes are distributed on all Chinese plum chromosomes except chromosome 3. Twelve (40 %) and five PmLEA genes are arranged in tandem and segmental duplications, respectively. The PmLEA genes could be divided into eight groups (LEA_1, LEA_2, LEA_3, LEA_4, LEA_5, PvLEA18, dehydrin and seed maturation protein). Ten gene conversion events were observed and most of them (70 %) were identified in dehydrin group. Most PmLEA genes were highly expressed in flower (22/30) and up-regulated by ABA treatment (19/30).

Sifting through genomes with iterative-sequence clustering produces a large, phylogenetically diverse protein-family resource

BACKGROUND

New computational resources are needed to manage the increasing volume of biological data from genome sequencing projects. One fundamental challenge is the ability to maintain a complete and current catalog of protein diversity. We developed a new approach for the identification of protein families that focuses on the rapid discovery of homologous protein sequences.

RESULTS

We implemented fully automated and high-throughput procedures to de novo cluster proteins into families based upon global alignment similarity. Our approach employs an iterative clustering strategy in which homologs of known families are sifted out of the search for new families. The resulting reduction in computational complexity enables us to rapidly identify novel protein families found in new genomes and to perform efficient, automated updates that keep pace with genome sequencing. We refer to protein families identified through this approach as "Sifting Families," or SFams. Our analysis of ~10.5 million protein sequences from 2,928 genomes identified 436,360 SFams, many of which are not represented in other protein family databases. We validated the quality of SFam clustering through statistical as well as network topology-based analyses.

CONCLUSIONS

We describe the rapid identification of SFams and demonstrate how they can be used to annotate genomes and metagenomes. The SFam database catalogs protein-family quality metrics, multiple sequence alignments, hidden Markov models, and phylogenetic trees. Our source code and database are publicly available and will be subject to frequent updates (http://edhar.genomecenter.ucdavis.edu/sifting_families/).

Identification and characterization of microRNAs and their target genes in Brassica oleracea

The microRNAs are a new class of small non-coding endogenous RNAs with lengths of approximately ~21 nt. MicroRNAs perform their biological function via the degradation of the target mRNAs or by inhibiting protein translation. Until recently, only limited numbers of miRNAs were identified in Brassica oleracea, a vegetable widely cultivated around the world. In present study, 193 potential miRNA candidates were identified from 17 expressed sequence tag (ESTs) and 152 genome survey sequences (GSSs) in B. oleracea. These miRNA candidates were classified into 70 families using a well-defined comparative genome-based computational analysis. Most miRNAs belong to the miRNA169, miR5021, miR156 and miR158 families. Of these, 36 miRNA families are firstly found in Brassica species. Around 1393 B. oleracea genes were predicted as candidate targets of 175 miRNAs. The mutual relationship between miRNAs and the candidate target genes was verified by checking differentially expression levels using quantitative real-time polymerase chain reaction (qRT-PCR) and 5' RLM-RACE analyses. These target genes participate in multiple biological and metabolic processes, including signal transduction, stress response, and plant development. Gene Ontology analysis shows that the 818, 514, and 265 target genes are involved in molecular functions, biological processes, and cellular component respectively. The Kyoto Encyclopedia of Genes and Genomes (KEGG) Pathway enrichment analysis suggests that these miRNAs might regulate 186 metabolic pathways, including those of lipid, energy, starch and sucrose, fatty acid and nitrogen.

Identification of conserved microRNAs and their targets in Chinese cabbage (Brassica rapa subsp. pekinensis)

The microRNAs (miRNAs) are a new class of small nonprotein-coding RNAs that have been identified to regulate gene expression at the post-transcriptional level by targeting mRNAs for degradation or by inhibiting protein translation. Until now, thousands of miRNAs have been identified in many plants species. However, only 23 miRNAs have been reported from the microRNA database in Chinese cabbage ( Brassica rapa subsp. pekinensis ), one of the most widely cultivated vegetables in China and East Asia. In the present study, 168 potential miRNAs, derived from 22 EST and 119 GSS sequences in Chinese cabbage were identified and classified into 38 miRNA families by well-defined computational analysis, in which most belonged to the miRNA1533, miRNA156, and miRNA2911 families. Totally, there are 129 identified miRNAs potentially targeting 1386 Chinese cabbage EST genes, which play roles in multiple biological and metabolic processes including metabolism, cell growth, signal transduction, stress response, and plant development. Gene ontology analysis, based on these target proteins, showed that 688, 532, and 287 genes were involved in molecular functions, biological processes, and cellular components, respectively. KEGG pathway analysis demonstrated that these miRNAs participated in 214 metabolism pathways, including, amongst others, plant–pathogen interaction, fatty acid metabolism, amino acid metabolism, nitrogen metabolism, plant hormone signal transduction.

Interactome-wide prediction of short, disordered protein interaction motifs in humans

Many of the specific functions of intrinsically disordered protein segments are mediated by Short Linear Motifs (SLiMs) interacting with other proteins. Well known examples include SLiMs that interact with 14-3-3, PDZ, SH2, SH3, and WW domains but the true extent and diversity of SLiM-mediated interactions is largely unknown. Here, we attempt to expand our knowledge of human SLiMs by applying in silico SLiM prediction to the human interactome. Combining data from seven different interaction databases, we analysed approximately 6000 protein-centred and 1600 domain-centred human interaction datasets of 3+ unrelated proteins that interact with a common partner. Results were placed in context through comparison to randomised datasets of similar size and composition. The search returned thousands of evolutionarily conserved, intrinsically disordered occurrences of hundreds of significantly enriched recurring motifs, including many that have never been previously identified (). In addition to True Positive results for at least 25 different known SLiMs, a striking number of "off-target" proteins/domains also returned significantly enriched known motifs. Often, this was due to the non-independence of the datasets, with many proteins sharing interaction partners or contributing interactions to multiple domain datasets. The majority of these motif classes, however, were also found to be significantly enriched in one or more randomised datasets. This highlights the need for care when interpreting motif predictions of this nature but also raises the possibility that SLiM occurrences may be successfully identified independently of interaction data. Although not as compositionally biased as previous studies, patterns matching known SLiMs tended to cluster into a few large groups of similar sequence, while novel predictions tended to be more distinctive and less abundant. Whether this is due to ascertainment bias or a true functional composition bias of SLiMs is not clear and warrants further investigation.

Phylogenetic-based propagation of functional annotations within the Gene Ontology consortium

The goal of the Gene Ontology (GO) project is to provide a uniform way to describe the functions of gene products from organisms across all kingdoms of life and thereby enable analysis of genomic data. Protein annotations are either based on experiments or predicted from protein sequences. Since most sequences have not been experimentally characterized, most available annotations need to be based on predictions. To make as accurate inferences as possible, the GO Consortium's Reference Genome Project is using an explicit evolutionary framework to infer annotations of proteins from a broad set of genomes from experimental annotations in a semi-automated manner. Most components in the pipeline, such as selection of sequences, building multiple sequence alignments and phylogenetic trees, retrieving experimental annotations and depositing inferred annotations, are fully automated. However, the most crucial step in our pipeline relies on software-assisted curation by an expert biologist. This curation tool, Phylogenetic Annotation and INference Tool (PAINT) helps curators to infer annotations among members of a protein family. PAINT allows curators to make precise assertions as to when functions were gained and lost during evolution and record the evidence (e.g. experimentally supported GO annotations and phylogenetic information including orthology) for those assertions. In this article, we describe how we use PAINT to infer protein function in a phylogenetic context with emphasis on its strengths, limitations and guidelines. We also discuss specific examples showing how PAINT annotations compare with those generated by other highly used homology-based methods.

Long-chain N-acyl amino acid synthases are linked to the putative PEP-CTERM/exosortase protein-sorting system in Gram-negative bacteria

Clones that encode the biosynthesis of long-chain N-acyl amino acids are frequently recovered from activity-based screens of soil metagenomic libraries. Members of a diverse set of enzymes referred to as N-acyl amino acid synthases are responsible for the production of all metagenome-derived N-acyl amino acids characterized to date. Based on the frequency at which N-acyl amino acid synthase genes have been identified from metagenomic samples, related genes are expected to be common throughout the global bacterial metagenome. Homologs of metagenome-derived N-acyl amino acid synthase genes are scarce, however, within the sequenced genomes of cultured bacterial species. Toward the goal of understanding the role(s) played by N-acyl amino acids in environmental bacteria, we looked for conserved genetic features that are positionally linked to metagenome-derived N-acyl amino acid synthase genes. This analysis revealed that N-acyl amino acid synthase genes are frequently found adjacent to genes predicted to encode PEP-CTERM motif-containing proteins and, in some cases, other conserved elements of the PEP-CTERM/exosortase system. Although relatively little is known about the PEP-CTERM/exosortase system, its core components are believed to represent the putative Gram-negative equivalent of the LPXTG/sortase protein-sorting system of Gram-positive bacteria. During the course of this investigation, we were able to provide evidence that an uncharacterized family of hypothetical acyltransferases, which had previously been linked to the PEP-CTERM/exosortase system by bioinformatics, is a new family of N-acyl amino acid synthases that is widely distributed among the PEP-CTERM/exosortase system-containing Proteobacteria.

Disease and phenotype data at Ensembl

Biological databases are an important resource for the life sciences community. Accessing the hundreds of databases supporting molecular biology and related fields is a daunting and time-consuming task. Integrating this information into one access point is a necessity for the life sciences community, which includes researchers focusing on human disease. Here we discuss the Ensembl genome browser, which acts as a single entry point with Graphical User Interface to data from multiple projects, including OMIM, dbSNP, and the NHGRI GWAS catalog. Ensembl provides a comprehensive source of annotation for the human genome, along with other species of biomedical interest. In this unit, we explore how to use the Ensembl genome browser in example queries related to human genetic diseases. Support protocols demonstrate quick sequence export using the BioMart tool.