Protein variety and functional diversity: Swiss-Prot annotation in its biological context

HPC-T-Annotator: an HPC tool for de novo transcriptome assembly annotation

BACKGROUND

The availability of transcriptomic data for species without a reference genome enables the construction of de novo transcriptome assemblies as alternative reference resources from RNA-Seq data. A transcriptome provides direct information about a species' protein-coding genes under specific experimental conditions. The de novo assembly process produces a unigenes file in FASTA format, subsequently targeted for the annotation. Homology-based annotation, a method to infer the function of sequences by estimating similarity with other sequences in a reference database, is a computationally demanding procedure.

RESULTS

To mitigate the computational burden, we introduce HPC-T-Annotator, a tool for de novo transcriptome homology annotation on high performance computing (HPC) infrastructures, designed for straightforward configuration via a Web interface. Once the configuration data are given, the entire parallel computing software for annotation is automatically generated and can be launched on a supercomputer using a simple command line. The output data can then be easily viewed using post-processing utilities in the form of Python notebooks integrated in the proposed software.

CONCLUSIONS

HPC-T-Annotator expedites homology-based annotation in de novo transcriptome assemblies. Its efficient parallelization strategy on HPC infrastructures significantly reduces computational load and execution times, enabling large-scale transcriptome analysis and comparison projects, while its intuitive graphical interface extends accessibility to users without IT skills.

ι-Carrageenan catabolism is initiated by key sulfatases in the marine bacterium Pseudoalteromonas haloplanktis LL1

Marine bacteria contribute substantially to cycle macroalgae polysaccharides in marine environments. Carrageenans are the primary cell wall polysaccharides of red macroalgae. The carrageenan catabolism mechanism and pathways are still largely unclear. Pseudoalteromonas is a representative bacterial genus that can utilize carrageenan. We previously isolated the strain Pseudoalteromonas haloplanktis LL1 that could grow on ι-carrageenan but produce no ι-carrageenase. Here, through a combination of bioinformatic, biochemical, and genetic analyses, we determined that P. haloplanktis LL1 processed a desulfurization-depolymerization sequential pathway for ι-carrageenan utilization, which was initiated by key sulfatases PhSulf1 and PhSulf2. PhSulf2 acted as an endo/exo-G4S (4-O-sulfation-β-D-galactopyranose) sulfatase, while PhSulf1 was identified as a novel endo-DA2S sulfatase that could function extracellularly. Because of the unique activity of PhSulf1 toward ι-carrageenan rather than oligosaccharides, P. haloplanktis LL1 was considered to have a distinct ι-carrageenan catabolic pathway compared to other known ι-carrageenan-degrading bacteria, which mainly employ multifunctional G4S sulfatases and exo-DA2S (2-O-sulfation-3,6-anhydro-α-D-galactopyranose) sulfatase for sulfate removal. Furthermore, we detected widespread occurrence of PhSulf1-encoding gene homologs in the global ocean, indicating the prevalence of such endo-acting DA2S sulfatases as well as the related ι-carrageenan catabolism pathway. This research provides valuable insights into the enzymatic processes involved in carrageenan catabolism within marine ecological systems.IMPORTANCECarrageenan is a type of linear sulfated polysaccharide that plays a significant role in forming cell walls of marine algae and is found extensively distributed throughout the world's oceans. To the best of our current knowledge, the ι-carrageenan catabolism in marine bacteria either follows the depolymerization-desulfurization sequential process initiated by ι-carrageenase or starts from the desulfurization step catalyzed by exo-acting sulfatases. In this study, we found that the marine bacterium Pseudoalteromonas haloplanktis LL1 processes a distinct pathway for ι-carrageenan catabolism employing a specific endo-acting DA2S-sulfatase PhSulf1 and a multifunctional G4S sulfatase PhSulf2. The unique PhSulf1 homologs appear to be widely present on a global scale, indicating the indispensable contribution of the marine bacteria containing the distinct ι-carrageenan catabolism pathway. Therefore, this study would significantly enrich our understanding of the molecular mechanisms underlying carrageenan utilization, providing valuable insights into the intricate roles of marine bacteria in polysaccharide cycling in marine environments.

Computational Protein Design - Where it goes?

Proteins have been playing a critical role in the regulation of diverse biological processes related to human life. With the increasing demand, functional proteins are sparse in this immense sequence space. Therefore, protein design has become an important task in various fields, including medicine, food, energy, materials, etc. Directed evolution has recently led to significant achievements. Molecular modification of proteins through directed evolution technology has significantly advanced the fields of enzyme engineering, metabolic engineering, medicine, and beyond. However, it is impossible to identify desirable sequences from a large number of synthetic sequences alone. As a result, computational methods, including data-driven machine learning and physics-based molecular modeling, have been introduced to protein engineering to produce more functional proteins. This review focuses on recent advances in computational protein design, highlighting the applicability of different approaches as well as their limitations.

Post-feeding transcriptomics reveals essential genes expressed in the midgut of the desert locust

The digestive tract constitutes an important interface between an animal's internal and external environment. In insects, available gut transcriptome studies are mostly exploratory or look at changes upon infection or upon exposure to xenobiotics, mainly performed in species belonging to holometabolan orders, such as Diptera, Lepidoptera or Coleoptera. By contrast, studies focusing on gene expression changes after food uptake and during digestion are underrepresented. We have therefore compared the gene expression profiles in the midgut of the desert locust, Schistocerca gregaria, between three different time points after feeding, i.e., 24 h (no active digestion), 10 min (the initial stage of feeding), and 2 h (active food digestion). The observed gene expression profiles were consistent with the polyphagous herbivorous lifestyle of this hemimetabolan (orthopteran) species. Our study reveals the upregulation of 576 genes 2 h post-feeding. These are mostly predicted to be associated with digestive physiology, such as genes encoding putative digestive enzymes or nutrient transporters, as well as genes putatively involved in immunity or in xenobiotic metabolism. The 10 min time point represented an intermediate condition, suggesting that the S. gregaria midgut can react rapidly at the transcriptional level to the presence of food. Additionally, our study demonstrated the critical importance of two transcripts that exhibited a significant upregulation 2 h post-feeding: the vacuolar-type H(+)-ATPase and the sterol transporter Niemann-Pick 1b protein, which upon RNAi-induced knockdown resulted in a marked increase in mortality. Their vital role and accessibility via the midgut lumen may make the encoded proteins promising insecticidal target candidates, considering that the desert locust is infamous for its huge migrating swarms that can devastate the agricultural production in large areas of Northern Africa, the Middle East, and South Asia. In conclusion, the transcriptome datasets presented here will provide a useful and promising resource for studying the midgut physiology of S. gregaria, a socio-economically important pest species.

Population transcriptogenomics highlights impaired metabolism and small population sizes in tree frogs living in the Chernobyl Exclusion Zone

BACKGROUND

Individual functional modifications shape the ability of wildlife populations to cope with anthropogenic environmental changes. But instead of adaptive response, human-altered environments can generate a succession of deleterious functional changes leading to the extinction of the population. To study how persistent anthropogenic changes impacted local species' population status, we characterised population structure, genetic diversity and individual response of gene expression in the tree frog Hyla orientalis along a gradient of radioactive contamination around the Chernobyl nuclear power plant.

RESULTS

We detected lower effective population size in populations most exposed to ionizing radiation in the Chernobyl Exclusion Zone that is not compensated by migrations from surrounding areas. We also highlighted a decreased body condition of frogs living in the most contaminated area, a distinctive transcriptomics signature and stop-gained mutations in genes involved in energy metabolism. While the association with dose will remain correlational until further experiments, a body of evidence suggests the direct or indirect involvement of radiation exposure in these changes.

CONCLUSIONS

Despite ongoing migration and lower total dose rates absorbed than at the time of the accident, our results demonstrate that Hyla orientalis specimens living in the Chernobyl Exclusion Zone are still undergoing deleterious changes, emphasizing the long-term impacts of the nuclear disaster.

A high-quality chromosome-level Eutrema salsugineum genome, an extremophile plant model

BACKGROUND

Eutrema salsugineum (2n = 14), a halophyte in the family Brassicaceae, is an attractive model to study abiotic stress tolerance in plants. Two versions of E. salsugineum genomes that previously reported were based on relatively short reads; thus, the repetitive regions were difficult to characterize.

RESULTS

We report the sequencing and assembly of the E. salsugineum (Shandong accession) genome using long-read sequencing and chromosome conformation capture data. We generated Oxford Nanopore long reads at high depth (> 60X) of genome coverage with additional short reads for error correction. The new assembly has a total size of 295.5 Mb with 52.8% repetitive sequences, and the karyotype of E. salsugineum is consistent with the ancestral translocation Proto-Calepineae Karyotype structure in both order and orientation. Compared with previous assemblies, this assembly has higher contiguity, especially in the centromere region. Based on this new assembly, we predicted 25,399 protein-coding genes and identified the positively selected genes associated with salt and drought stress responses.

CONCLUSION

The new genome assembly will provide a valuable resource for future genomic studies and facilitate comparative genomic analysis with other plants.

Evolutionary Conservation and Transcriptome Analyses Attribute Perenniality and Flowering to Day-Length Responsive Genes in Bulbous Barley (Hordeum bulbosum)

Rapid population growth and dramatic climatic turnovers are challenging global crop production. These challenges are spurring plant breeders to enhance adaptation and sustainability of major crops. One intriguing approach is to turn annual systems into perennial ones, yet long-term classical breeding efforts to induce perenniality have achieved limited success. Here, we report the results of our investigation of the genetic basis of bulb formation in the nonmodel organism Hordeum bulbosum, a perennial species closely related to barley. To identify candidate genes that regulate bulb formation in H. bulbosum, we applied two complementary approaches. First, we explored the evolutionary conservation of expressed genes among annual Poaceae species. Next, we assembled a reference transcriptome for H. bulbosum and conducted a differential expression (DE) analysis before and after stimulating bulb initiation. Low conservation was identified in genes related to perenniality in H. bulbosum compared with other species, including bulb development and sugar accumulation genes. We also inspected these genes using a DE analysis, which enabled identification of additional genes responsible for bulb initiation and flowering regulation. We propose a molecular model for the regulation of bulb formation involving storage organ development and starch biosynthesis genes. The high conservation observed along a major part of the pathway between H. bulbosum and barley suggests a potential for the application of biotechnological techniques to accelerate breeding toward perenniality in barley.

Influence of bisphenol A on growth and metabolism of Vicia faba ssp. minor seedlings depending on lighting conditions

The effect of one of anthropogenic pollutants, i.e., 4,4'-isopropylidenediphenol, called 2,2-bis (4-hydroxyphenyl) propane (BPA), at 30 and 120 mg L^-1 concentrations in the darkness (DK) or dark/light (DK/LT) on growth and selected elements of metabolism of seedlings and leaf discs of Vicia faba ssp. minor was studied. Treatment with 120 mg L^-1 BPA had greater effects which were reflected by increase in the number of necrotic changes in roots and stems as well as in leaf discs and reduction of the length of roots DK and DK/LT, and volume of roots in the DK group. However, minimal and no influence on the fresh and dry weight of roots and stems in plants growing under both types of lighting conditions were observed. In both DK and DK/LT groups these effects were correlated with reduced amounts of storage and cell wall-bound sugars as well as of proteins while in the DK/LT additionally with reduced soluble sugar levels in the roots and increased amounts of hydrogen peroxide and phenols in roots and stems as well as in treatment solutions, where these compounds were released. We suggest that endogenous phenols and BPA can be metabolised in roots and stems to quinones. It seems that TB-1,4-BQ, is the one of that of the five studied quinones. We expect that the results of this paper will help to answer the following question: does the phytomeliorative and phytosanitative V. faba ssp. minor plant is enough to be resistant on negative effects, and to be useful to reduce increasing amount of BPA in the environment?

Systematic analysis of HD-ZIP transcription factors in sesame genome and gene expression profiling of SiHD-ZIP class I entailing drought stress responses at early seedling stage

BACKGROUND

Sesame is an ancient oilseed crop, known for its high oil content and quality. Its sensitivity to drought at early seedling stage is one of the limiting factors affecting its world-wide growth and productivity. Among plant specific transcription factors, the association of HD-ZIPs with sesame drought responses at early seedling stage is not well-established yet and is very important to develop our molecular understanding on sesame drought tolerance.

METHODS AND RESULTS

In this study, total 61 sesame HD-ZIP proteins were identified, based on their protein sequence homology with Arabidopsis and protein domain(s) architecture prediction, followed by their phylogenetic, conserved domain(s) motifs and gene structure analyses to classify them into four classes (HD-ZIP Class I-IV). HD-ZIP Class I was also subdivided into four subgroups: α (SiHZ25, SiHZ43, SiHZ9 and SiHZ16), β1 (SiHZ10, SiHZ30, SiHZ32 and SiHZ26), β2 (SiHZ42 and SiHZ45) and γ (SiHZ17, SiHZ7 and SiHZ35) by a comparative phylogenetic analysis of sesame with Arabidopsis and maize. Afterwards, twenty-one days old sesame seedlings were exposed to drought stress by withholding water for 7 days (when soil moisture content reduced to ~16%) and gene expression of HD-ZIP Class I (13 members) was performed in well- watered (control) and drought stressed seedlings. The gene expression analysis showed that the expressions of SiHZ7 (6.8 fold) and SiHZ35 (2.6 fold) from γ subgroup were significantly high in drought seedlings.

CONCLUSIONS

This study is useful in demonstrating the role of SiHD-ZIP Class I in sesame drought responses at early seedling stage and to develop its novel drought tolerant varieties.

Degeneracy measures in biologically plausible random Boolean networks

Background

Degeneracy—the ability of structurally different elements to perform similar functions—is a property of many biological systems. Highly degenerate systems show resilience to perturbations and damage because the system can compensate for compromised function due to reconfiguration of the underlying network dynamics. Degeneracy thus suggests how biological systems can thrive despite changes to internal and external demands. Although degeneracy is a feature of network topologies and seems to be implicated in a wide variety of biological processes, research on degeneracy in biological networks is mostly limited to weighted networks. In this study, we test an information theoretic definition of degeneracy on random Boolean networks, frequently used to model gene regulatory networks. Random Boolean networks are discrete dynamical systems with binary connectivity and thus, these networks are well-suited for tracing information flow and the causal effects. By generating networks with random binary wiring diagrams, we test the effects of systematic lesioning of connections and perturbations of the network nodes on the degeneracy measure.

Results

Our analysis shows that degeneracy, on average, is the highest in networks in which ~ 20% of the connections are lesioned while 50% of the nodes are perturbed. Moreover, our results for the networks with no lesions and the fully-lesioned networks are comparable to the degeneracy measures from weighted networks, thus we show that the degeneracy measure is applicable to different networks.

Conclusions

Such a generalized applicability implies that degeneracy measures may be a useful tool for investigating a wide range of biological networks and, therefore, can be used to make predictions about the variety of systems’ ability to recover function.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04601-5.

Differential gene expression associated with a floral scent polymorphism in the evening primrose Oenothera harringtonii (Onagraceae)

Background

Plant volatiles play an important role in both plant-pollinator and plant-herbivore interactions. Intraspecific polymorphisms in volatile production are ubiquitous, but studies that explore underlying differential gene expression are rare. Oenothera harringtonii populations are polymorphic in floral emission of the monoterpene (R)-(−)-linalool; some plants emit (R)-(−)-linalool (linalool+ plants) while others do not (linalool- plants). However, the genes associated with differential production of this floral volatile in Oenothera are unknown. We used RNA-Seq to broadly characterize differential gene expression involved in (R)-(−)-linalool biosynthesis. To identify genes that may be associated with the polymorphism for this trait, we used RNA-Seq to compare gene expression in six different Oenothera harringtonii tissues from each of three linalool+ and linalool- plants.

Results

Three clusters of differentially expressed genes were enriched for terpene synthase activity: two were characterized by tissue-specific upregulation and one by upregulation only in plants with flowers that produce (R)-(−)-linalool. A molecular phylogeny of all terpene synthases identified two putative (R)-(−)-linalool synthase transcripts in Oenothera harringtonii, a single allele of which is found exclusively in linalool+ plants.

Conclusions

By using a naturally occurring polymorphism and comparing different tissues, we were able to identify candidate genes putatively involved in the biosynthesis of (R)-(−)-linalool. Expression of these genes in linalool- plants, while low, suggests a regulatory polymorphism, rather than a population-specific loss-of-function allele. Additional terpene biosynthesis-related genes that are up-regulated in plants that emit (R)-(−)-linalool may be associated with herbivore defense, suggesting a potential economy of scale between plant reproduction and defense.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08370-6.

Locating ligand binding sites in G-protein coupled receptors using combined information from docking and sequence conservation

GPCRs (G-protein coupled receptors) are the largest family of drug targets and share a conserved structure. Binding sites are unknown for many important GPCR ligands due to the difficulties of GPCR recombinant expression, biochemistry, and crystallography. We describe our approach, ConDockSite, for predicting ligand binding sites in class A GPCRs using combined information from surface conservation and docking, starting from crystal structures or homology models. We demonstrate the effectiveness of ConDockSite on crystallized class A GPCRs such as the beta2 adrenergic and A2A adenosine receptors. We also demonstrate that ConDockSite successfully predicts ligand binding sites from high-quality homology models. Finally, we apply ConDockSite to predict the ligand binding sites on a structurally uncharacterized GPCR, GPER, the G-protein coupled estrogen receptor. Most of the sites predicted by ConDockSite match those found in other independent modeling studies. ConDockSite predicts that four ligands bind to a common location on GPER at a site deep in the receptor cleft. Incorporating sequence conservation information in ConDockSite overcomes errors introduced from physics-based scoring functions and homology modeling.

Whole-Genome and Transcriptome Sequencing of Phlebopus portentosus Reveals Its Associated Ectomycorrhizal Niche and Conserved Pathways Involved in Fruiting Body Development

Phlebopus portentosus (Berk. and Broome) Boedijin, a widely consumed mushroom in China and Thailand, is the first species in the order Boletaceae to have been industrially cultivated on a large scale. However, to date, the lignocellulose degradation system and molecular basis of fruiting body development in P. portentosus have remained cryptic. In the present study, genome and transcriptome sequencing of P. portentosus was performed during the mycelium (S), primordium (P), and fruiting body (F) stages. A genome of 32.74 Mb with a 48.92% GC content across 62 scaffolds was obtained. A total of 9,464 putative genes were predicted from the genome, of which the number of genes related to plant cell wall-degrading enzymes was much lower than that of some saprophytic mushrooms with specific ectomycorrhizal niches. Principal component analysis of RNA-Seq data revealed that the gene expression profiles at all three stages were different. The low expression of plant cell wall-degrading genes also confirmed the limited ability to degrade lignocellulose. The expression profiles also revealed that some conserved and specific pathways were enriched in the different developmental stages of P. portentosus. Starch and sucrose metabolic pathways were enriched in the mycelium stage, while DNA replication, the proteasome and MAPK signaling pathways may be associated with maturation. These results provide a new perspective for understanding the key pathways and hub genes involved in P. portentosus development.

The Serological Cross-Detection of Bat-Borne Hantaviruses: A Valid Strategy or Taking Chances?

Bats are hosts of a range of viruses, and their great diversity and unique characteristics that distinguish them from all other mammals have been related to the maintenance, evolution, and dissemination of these pathogens. Recently, very divergent hantaviruses have been discovered in distinct species of bats worldwide, but their association with human disease remains unclear. Considering the low success rates of detecting hantavirus RNA in bat tissues and that to date no hantaviruses have been isolated from bat samples, immunodiagnostic tools could be very helpful to understand pathogenesis, epidemiology, and geographic range of bat-borne hantaviruses. In this sense, we aimed to identify in silico immunogenic B-cell epitopes present on bat-borne hantaviruses nucleoprotein (NP) and verify if they are conserved among them and other selected members of Mammantavirinae, using a combination of (the three most used) different prediction algorithms, ELLIPRO, Discotope 2.0, and PEPITO server. To support our data, we in silico modeled 3D structures of NPs from representative members of bat-borne hantaviruses, using comparative and ab initio methods due to the absence of crystallographic structures of studied proteins or similar models in the Protein Data Bank. Our analysis demonstrated the antigenic complexity of the bat-borne hantaviruses group, showing a low sequence conservation of epitopes among members of its own group and a minor conservation degree in comparison to Orthohantavirus, with a recognized importance to public health. Our data suggest that the use of recombinant rodent-borne hantavirus NPs to cross-detect antibodies against bat- or shrew-borne viruses could underestimate the real impact of this virus in nature.

Single-cell transcriptome profiling of buffelgrass (Cenchrus ciliaris) eggs unveils apomictic parthenogenesis signatures

Apomixis, a type of asexual reproduction in angiosperms, results in progenies that are genetically identical to the mother plant. It is a highly desirable trait in agriculture due to its potential to preserve heterosis of F₁ hybrids through subsequent generations. However, no major crops are apomictic. Deciphering mechanisms underlying apomixis becomes one of the alternatives to engineer self-reproducing capability into major crops. Parthenogenesis, a major component of apomixis, commonly described as the ability to initiate embryo formation from the egg cell without fertilization, also can be valuable in plant breeding for doubled haploid production. A deeper understanding of transcriptional differences between parthenogenetic and sexual or non-parthenogenetic eggs can assist with pathway engineering. By conducting laser capture microdissection-based RNA-seq on sexual and parthenogenetic egg cells on the day of anthesis, a de novo transcriptome for the Cenchrus ciliaris egg cells was created, transcriptional profiles that distinguish the parthenogenetic egg from its sexual counterpart were identified, and functional roles for a few transcription factors in promoting natural parthenogenesis were suggested. These transcriptome data expand upon previous gene expression studies and will be a resource for future research on the transcriptome of egg cells in parthenogenetic and sexual genotypes.

Regulatory Potential of Long Non-Coding RNAs (lncRNAs) in Boar Spermatozoa with Good and Poor Freezability

Long non-coding RNAs (lncRNAs) are suggested to play an important role in the sperm biological processes. We performed de novo transcriptome assembly to characterize lncRNAs in spermatozoa, and to investigate the role of the potential target genes of the differentially expressed lncRNAs (DElncRNAs) in sperm freezability. We detected approximately 4007 DElncRNAs, which were differentially expressed in spermatozoa from boars classified as having good and poor semen freezability (GSF and PSF, respectively). Most of the DElncRNAs were upregulated in boars of the PSF group and appeared to significantly affect the sperm's response to the cryopreservation conditions. Furthermore, we predicted that the potential target genes were regulated by DElncRNAs in cis or trans. It was found that DElncRNAs of both freezability groups had potential cis- and trans-regulatory effects on different protein-coding genes, such as COX7A2L, TXNDC8 and SOX-7. Gene Ontology (GO) enrichment revealed that the DElncRNA target genes are associated with numerous biological processes, including signal transduction, response to stress, cell death (apoptosis), motility and embryo development. Significant differences in the de novo assembled transcriptome expression profiles of the DElncRNAs between the freezability groups were confirmed by quantitative real-time PCR analysis. This study reveals the potential effects of protein-coding genes of DElncRNAs on sperm functions, which could contribute to further research on their relevance in semen freezability.

Transcriptome sequencing and analysis reveals the molecular response to selenium stimuli in Pueraria lobata (willd.) Ohwi

Pueraria lobata (willd.) Ohwi is a consumable selenium-enriched plant used for medicinal purposes. The molecular response to selenium (Se) stimuli in P. lobata is currently unknown. We used RNA-Seq to identify potential genes involved in selenite metabolism and analyzed their expression profiles. We obtained a total of 150,567 unigenes, of which 90,961 were annotated, including 16 structural genes, 14 sulfate transporters, and 13 phosphate transporters that may be involved in Se metabolism, and 33 candidate structural genes involved in isoflavone biosynthesis. The genes with a —foldchange— >2 and q value <0.05 after sodium selenite treatment were identified as differentially expressed genes (DEGs). We obtained a total of 4,246 DEGs, which were enriched in GO terms that included “response to stimulus”, “response to stress”, “signal transduction”, “response to abiotic stimulus”, and “response to chemical”. Of the 4,246 DEGs, one sulfate transporter and five phosphate transporter genes involved Se metabolism, and nine structural genes involved in isoflavone biosynthesis were up-regulated. The expression patterns of 10 DEGs were selected randomly and validated using qRT-PCR. The Pearson Correlation Coefficient (r) was 0.86, indicating the reliability of RNA-Seq results. 22 Reactive Oxygen Species (ROS) scavenging DEGs were found, 11 of which were up-regulated. 436, 624 transcription factors (TFs) correlated with structural genes were identified that may be involved in Se and isoflavone biosynthesis, respectively, using r (r > 0.7 or r <  − 0.7). 556 TFs were related to at least one sulfate and phosphate transporter. Our results provided a comprehensive description of gene expression and regulation in response to Se stimuli in P. lobata.

Utilization of Tissue Ploidy Level Variation in de Novo Transcriptome Assembly of Pinus sylvestris

Compared to angiosperms, gymnosperms lag behind in the availability of assembled and annotated genomes. Most genomic analyses in gymnosperms, especially conifer tree species, rely on the use of de novo assembled transcriptomes. However, the level of allelic redundancy and transcript fragmentation in these assembled transcriptomes, and their effect on downstream applications have not been fully investigated. Here, we assessed three assembly strategies for short-reads data, including the utility of haploid megagametophyte tissue during de novo assembly as single-allele guides, for six individuals and five different tissues in Pinus sylvestris. We then contrasted haploid and diploid tissue genotype calls obtained from the assembled transcriptomes to evaluate the extent of paralog mapping. The use of the haploid tissue during assembly increased its completeness without reducing the number of assembled transcripts. Our results suggest that current strategies that rely on available genomic resources as guidance to minimize allelic redundancy are less effective than the application of strategies that cluster redundant assembled transcripts. The strategy yielding the lowest levels of allelic redundancy among the assembled transcriptomes assessed here was the generation of SuperTranscripts with Lace followed by CD-HIT clustering. However, we still observed some levels of heterozygosity (multiple gene fragments per transcript reflecting allelic redundancy) in this assembled transcriptome on the haploid tissue, indicating that further filtering is required before using these assemblies for downstream applications. We discuss the influence of allelic redundancy when these reference transcriptomes are used to select regions for probe design of exome capture baits and for estimation of population genetic diversity.

Plasma proteomic and autoantibody profiles reveal the proteomic characteristics involved in longevity families in Bama, China

BACKGROUND

Chinese Bama Yao Autonomous County is a well-known longevity region in the world. In the past 30 years, population and genome studies were undertaken to investigate the secret of longevity and showed that longevity is the result of a combination of multiple factors, such as genetic, environmental and other causes. In this study, characteristics of the blood plasma proteomic and autoantibody profiles of people from Bama longevity family were investigated.

METHODS

Sixty-six plasma donors from Chinese Bama longevity area were recruited in this study. Thirty-three offsprings of longevous families were selected as case studies (Longevous group) and 33 ABO (blood type), age, and gender-matched subjects from non-longevous families were selected as controls (Normal group). Each group contains 3 biological replicates. Tandem mass tag-based proteomic technique was used to investigate the differentially expressed plasma proteins between the two groups. The auto-reactive IgG antibody profiles of the 3 pooled samples in each group were revealed by human proteome microarrays with 17,000 recombinant human proteins.

RESULTS

Firstly, 525 plasma proteins were quantified and 12 proteins were discovered differentially expressed between the two groups. Secondly, more than 500 proteins were recognized by plasma antibodies, 14 proteins ware differentially reacted with the autoantibodies in the two groups. Bioinformatics analysis showed some of the differential proteins and targeted autoantigens were involved in cancer, cardiovascular disease and immunity.

CONCLUSIONS

Proteomic and autoantibody profiles varied between the offspring of longevous and normal families which are from the same area and shared the same environmental factors. The identified differences were reported to be involved in several physiological and pathological pathways. The identified proteins will contribute to a better understanding of the proteomic characteristics of people from Bama longevous area and a revelation of the molecular mechanisms of longevity.

Metabolic reconstruction and experimental verification of glucose utilization in Desulfurococcus amylolyticus DSM 16532

Desulfurococcus amylolyticus DSM 16532 is an anaerobic and hyperthermophilic crenarchaeon known to grow on a variety of different carbon sources, including monosaccharides and polysaccharides. Furthermore, D. amylolyticus is one of the few archaea that are known to be able to grow on cellulose. Here, we present the metabolic reconstruction of D. amylolyticus’ central carbon metabolism. Based on the published genome, the metabolic reconstruction was completed by integrating complementary information available from the KEGG, BRENDA, UniProt, NCBI, and PFAM databases, as well as from available literature. The genomic analysis of D. amylolyticus revealed genes for both the classical and the archaeal version of the Embden-Meyerhof pathway. The metabolic reconstruction highlighted gaps in carbon dioxide-fixation pathways. No complete carbon dioxide-fixation pathway such as the reductive citrate cycle or the dicarboxylate-4-hydroxybutyrate cycle could be identified. However, the metabolic reconstruction indicated that D. amylolyticus harbors all genes necessary for glucose metabolization. Closed batch experimental verification of glucose utilization by D. amylolyticus was performed in chemically defined medium. The findings from in silico analyses and from growth experiments are discussed with respect to physiological features of hyperthermophilic organisms.

Comparative whole genome re-sequencing analysis in upland New Rice for Africa: insights into the breeding history and respective genome compositions

BACKGROUND

Increasing rice demand is one of the consequences of the steadily improving socio-economic status of the African countries. New Rice for Africa (NERICA), which are interspecific hybrids between Asian and African rice varieties, are one of successful breeding products utilizing biodiversity across the two different rice crop species. Upland NERICA varieties (NU) exhibit agronomic traits of value for the harsh eco-geography, including shorter duration, higher yield and stress tolerance, compared to local African varieties. However, the molecular basis of the traits in NU varieties is largely unknown.

RESULTS

Whole genome re-sequencing was performed of four NU lines (3, 4, 5, and 7) and for the parental Oryza sativa WAB56-104 and Oryza glaberrima CG14. The k-mer analysis predicted large genomes for the four NU lines, most likely inherited from WAB56-104. Approximately 3.1, 0.10, and 0.40 million single nucleotide polymorphisms, multi nucleotide polymorphisms, and short insertions/deletions were mined between the parental lines, respectively. Integrated analysis with another four NU lines (1, 2, 8, and 9) showed that the ratios of the donor CG14 allelic sites in the NU lines ranged from 1.3 to 9.8%. High resolution graphical genotype indicated genome-level similarities and common genetic events during the breeding process: five xyloglucan fucosyltransferase from O. glaberrima were introgressed in common. Segregation of genic segments revealed potential causal genes for some agronomic traits including grain shattering, awnness, susceptibility to bacterial leaf bright, and salt tolerance. Analysis of unmapped sequences against the reference cultivar Nipponbare indicated existence of unique genes for pathogen and abiotic stress resistance in the NU varieties.

CONCLUSIONS

The results provide understanding of NU genomes for rice improvement for Africa reinforcing local capacity for food security and insights into molecular events in breeding of interspecific hybrid crops.

Comparative Metatranscriptomics of Wheat Rhizosphere Microbiomes in Disease Suppressive and Non-suppressive Soils for Rhizoctonia solani AG8

The soilborne fungus Rhizoctonia solani anastomosis group (AG) 8 is a major pathogen of grain crops resulting in substantial production losses. In the absence of resistant cultivars of wheat or barley, a sustainable and enduring method for disease control may lie in the enhancement of biological disease suppression. Evidence of effective biological control of R. solani AG8 through disease suppression has been well documented at our study site in Avon, South Australia. A comparative metatranscriptomic approach was applied to assess the taxonomic and functional characteristics of the rhizosphere microbiome of wheat plants grown in adjacent fields which are suppressive and non-suppressive to the plant pathogen R. solani AG8. Analysis of 12 rhizosphere metatranscriptomes (six per field) was undertaken using two bioinformatic approaches involving unassembled and assembled reads. Differential expression analysis showed the dominant taxa in the rhizosphere based on mRNA annotation were Arthrobacter spp. and Pseudomonas spp. for non-suppressive samples and Stenotrophomonas spp. and Buttiauxella spp. for the suppressive samples. The assembled metatranscriptome analysis identified more differentially expressed genes than the unassembled analysis in the comparison of suppressive and non-suppressive samples. Suppressive samples showed greater expression of a polyketide cyclase, a terpenoid biosynthesis backbone gene (dxs) and many cold shock proteins (csp). Non-suppressive samples were characterised by greater expression of antibiotic genes such as non-heme chloroperoxidase (cpo) which is involved in pyrrolnitrin synthesis, and phenazine biosynthesis family protein F (phzF) and its transcriptional activator protein (phzR). A large number of genes involved in detoxifying reactive oxygen species (ROS) and superoxide radicals (sod, cat, ahp, bcp, gpx1, trx) were also expressed in the non-suppressive rhizosphere samples most likely in response to the infection of wheat roots by R. solani AG8. Together these results provide new insight into microbial gene expression in the rhizosphere of wheat in soils suppressive and non-suppressive to R. solani AG8. The approach taken and the genes involved in these functions provide direction for future studies to determine more precisely the molecular interplay of plant-microbe-pathogen interactions with the ultimate goal of the development of management options that promote beneficial rhizosphere microflora to reduce R. solani AG8 infection of crops.

ILDR2 Is a Novel B7-like Protein That Negatively Regulates T Cell Responses

The B7-like protein family members play critical immunomodulatory roles and constitute attractive targets for the development of novel therapies for human diseases. We identified Ig-like domain-containing receptor (ILDR)2 as a novel B7-like protein with robust T cell inhibitory activity, expressed in immune cells and in immune-privileged and inflamed tissues. A fusion protein, consisting of ILDR2 extracellular domain with an Fc fragment, that binds to a putative counterpart on activated T cells showed a beneficial effect in the collagen-induced arthritis model and abrogated the production of proinflammatory cytokines and chemokines in autologous synovial-like cocultures of macrophages and cytokine-stimulated T cells. Collectively, these findings point to ILDR2 as a novel negative regulator for T cells, with potential roles in the development of immune-related diseases, including autoimmunity and cancer.

A Tissue-Mapped Axolotl De Novo Transcriptome Enables Identification of Limb Regeneration Factors

Mammals have extremely limited regenerative capabilities; however, axolotls are profoundly regenerative and can replace entire limbs. The mechanisms underlying limb regeneration remain poorly understood, partly because the enormous and incompletely sequenced genomes of axolotls have hindered the study of genes facilitating regeneration. We assembled and annotated a de novo transcriptome using RNA-sequencing profiles for a broad spectrum of tissues that is estimated to have near-complete sequence information for 88% of axolotl genes. We devised expression analyses that identified the axolotl orthologs of cirbp and kazald1 as highly expressed and enriched in blastemas. Using morpholino anti-sense oligonucleotides, we find evidence that cirbp plays a cytoprotective role during limb regeneration whereas manipulation of kazald1 expression disrupts regeneration. Our transcriptome and annotation resources greatly complement previous transcriptomic studies and will be a valuable resource for future research in regenerative biology.

Identification of regenerative roadblocks via repeat deployment of limb regeneration in axolotls

Axolotl salamanders are powerful models for understanding how regeneration of complex body parts can be achieved, whereas mammals are severely limited in this ability. Factors that promote normal axolotl regeneration can be examined in mammals to determine if they exhibit altered activity in this context. Furthermore, factors prohibiting axolotl regeneration can offer key insight into the mechanisms present in regeneration-incompetent species. We sought to determine if we could experimentally compromise the axolotl's ability to regenerate limbs and, if so, discover the molecular changes that might underlie their inability to regenerate. We found that repeated limb amputation severely compromised axolotls' ability to initiate limb regeneration. Using RNA-seq, we observed that a majority of differentially expressed transcripts were hyperactivated in limbs compromised by repeated amputation, suggesting that mis-regulation of these genes antagonizes regeneration. To confirm our findings, we additionally assayed the role of amphiregulin, an EGF-like ligand, which is aberrantly upregulated in compromised animals. During normal limb regeneration, amphiregulin is expressed by the early wound epidermis, and mis-expressing this factor lead to thickened wound epithelium, delayed initiation of regeneration, and severe regenerative defects. Collectively, our results suggest that repeatedly amputated limbs may undergo a persistent wound healing response, which interferes with their ability to initiate the regenerative program. These findings have important implications for human regenerative medicine.

Selected Bioinformatic Tools and MS (MALDI-TOF, PMF) Techniques Used in the Strategy for the Identification of Oat Proteins After 2-DE

Computer analysis of protein maps obtained from the separation of proteins with two-dimensional polyacrylamide gel electrophoresis (2-DE), in combination with mass spectrometry (MS) analysis and selected bioinformatic tools is used in the strategy for the identification of oat proteins. In proteomic research the most often used MS technique is the combination of ion sources: matrix-assisted laser desorption/ionization (MALDI) and the analyzer of the time of flight (TOF), i.e., MALDI-TOF MS.This chapter describes the possibilities of the use of selected bioinformatic tools (UniProtKB database, ProtParam, Compute pI/MW programs) for initial identification of separated oat proteins (especially prolamin fractions) with the 2-DE technique. Also the procedure of preparation of samples obtained from cut out protein spots for analysis with the MALDI-TOF MS and peptide mass fingerprinting (PMF) technique is presented.Among oat prolamins separated with the 2-DE technique (see Chapter 17 ), 13 protein spots are considered to be the most characteristic (range of MW 27.0-34.6 kDa, pI 5.7-7.6) for this fraction of proteins. Among them there are four protein spots (MW 27.0-28.0 kDa) and two spots (MW 31.4-32.1 kDa) which can correspond to avenins (Accession numbers (AC) in UniProtKB: L0L5I0, I4EP88, I4EP64, L0L4I8 and F2Q9W5, L0L6J0, respectively).

Exploring lateral genetic transfer among microbial genomes using TF-IDF

Many microbes can acquire genetic material from their environment and incorporate it into their genome, a process known as lateral genetic transfer (LGT). Computational approaches have been developed to detect genomic regions of lateral origin, but typically lack sensitivity, ability to distinguish donor from recipient, and scalability to very large datasets. To address these issues we have introduced an alignment-free method based on ideas from document analysis, term frequency-inverse document frequency (TF-IDF). Here we examine the performance of TF-IDF on three empirical datasets: 27 genomes of Escherichia coli and Shigella, 110 genomes of enteric bacteria, and 143 genomes across 12 bacterial and three archaeal phyla. We investigate the effect of k-mer size, gap size and delineation of groups on the inference of genomic regions of lateral origin, finding an interplay among these parameters and sequence divergence. Because TF-IDF identifies donor groups and delineates regions of lateral origin within recipient genomes, aggregating these regions by gene enables us to explore, for the first time, the mosaic nature of lateral genes including the multiplicity of biological sources, ancestry of transfer and over-writing by subsequent transfers. We carry out Gene Ontology enrichment tests to investigate which biological processes are potentially affected by LGT.

Patterns of Genome-Wide Variation in Glossina fuscipes fuscipes Tsetse Flies from Uganda

The tsetse fly Glossina fuscipes fuscipes (Gff) is the insect vector of the two forms of Human African Trypanosomiasis (HAT) that exist in Uganda. Understanding Gff population dynamics, and the underlying genetics of epidemiologically relevant phenotypes is key to reducing disease transmission. Using ddRAD sequence technology, complemented with whole-genome sequencing, we developed a panel of ∼73,000 single-nucleotide polymorphisms (SNPs) distributed across the Gff genome that can be used for population genomics and to perform genome-wide-association studies. We used these markers to estimate genomic patterns of linkage disequilibrium (LD) in Gff, and used the information, in combination with outlier-locus detection tests, to identify candidate regions of the genome under selection. LD in individual populations decays to half of its maximum value (r(2) max/2) between 1359 and 2429 bp. The overall LD estimated for the species reaches r(2) max/2 at 708 bp, an order of magnitude slower than in Drosophila Using 53 infected (Trypanosoma spp.) and uninfected flies from four genetically distinct Ugandan populations adapted to different environmental conditions, we were able to identify SNPs associated with the infection status of the fly and local environmental adaptation. The extent of LD in Gff likely facilitated the detection of loci under selection, despite the small sample size. Furthermore, it is probable that LD in the regions identified is much higher than the average genomic LD due to strong selection. Our results show that even modest sample sizes can reveal significant genetic associations in this species, which has implications for future studies given the difficulties of collecting field specimens with contrasting phenotypes for association analysis.

Micropollutant degradation via extracted native enzymes from activated sludge

A procedure was developed to assess the biodegradation of micropollutants in cell-free lysates produced from activated sludge of a municipal wastewater treatment plant (WWTP). This proof-of-principle provides the basis for further investigations of micropollutant biodegradation via native enzymes in a solution of reduced complexity, facilitating downstream protein analysis. Differently produced lysates, containing a variety of native enzymes, showed significant enzymatic activities of acid phosphatase, β-galactosidase and β-glucuronidase in conventional colorimetric enzyme assays, whereas heat-deactivated controls did not. To determine the enzymatic activity towards micropollutants, 20 compounds were spiked to the cell-free lysates under aerobic conditions and were monitored via LC-ESI-MS/MS. The micropollutants were selected to span a wide range of different biodegradabilities in conventional activated sludge treatment via distinct primary degradation reactions. Of the 20 spiked micropollutants, 18 could be degraded by intact sludge under assay conditions, while six showed reproducible degradation in the lysates compared to the heat-deactivated negative controls: acetaminophen, N-acetyl-sulfamethoxazole (acetyl-SMX), atenolol, bezafibrate, erythromycin and 10,11-dihydro-10-hydroxycarbamazepine (10-OH-CBZ). The primary biotransformation of the first four compounds can be attributed to amide hydrolysis. However, the observed biotransformations in the lysates were differently influenced by experimental parameters such as sludge pre-treatment and the addition of ammonium sulfate or peptidase inhibitors, suggesting that different hydrolase enzymes were involved in the primary degradation, among them possibly peptidases. Furthermore, the transformation of 10-OH-CBZ to 9-CA-ADIN was caused by a biologically-mediated oxidation, which indicates that in addition to hydrolases further enzyme classes (probably oxidoreductases) are present in the native lysates. Although the full variety of indigenous enzymatic activity of the activated sludge source material could not be restored, experimental modifications, e.g. different lysate filtration, significantly enhanced specific enzyme activities (e.g. >96% removal of the antibiotic erythromycin). Therefore, the approach presented in this study provides the experimental basis for a further elucidation of the enzymatic processes underlying wastewater treatment on the level of native proteins.

Exploration of Structural and Functional Variations Owing to Point Mutations in α-NAGA

Schindler disease is a lysosomal storage disorder caused due to deficiency or defective activity of alpha-N-acetylgalactosaminidase (α-NAGA). Mutations in gene encoding α-NAGA cause wide range of diseases, characterized with mild to severe clinical features. Molecular effects of these mutations are yet to be explored in detail. Therefore, this study was focused on four missense mutations of α-NAGA namely, S160C, E325K, R329Q and R329W. Native and mutant structures of α-NAGA were analysed to determine geometrical deviations such as the contours of root mean square deviation, root mean square fluctuation, percentage of residues in allowed regions of Ramachandran plot and solvent accessible surface area, using conformational sampling technique. Additionally, global energy-minimized structures of native and mutants were further analysed to compute their intra-molecular interactions, hydrogen bond dilution and distribution of secondary structure. In addition, docking studies were also performed to determine variations in binding energies between native and mutants. The deleterious effects of mutants were evident due to variations in their active site residues pertaining to spatial conformation and flexibility, comparatively. Hence, variations exhibited by mutants, namely S160C, E325K, R329Q and R329W to that of native, consequently, lead to the detrimental effects causing Schindler disease. This study computationally explains the underlying reasons for the pathogenesis of the disease, thereby aiding future researchers in drug development and disease management.

BacHbpred: Support Vector Machine Methods for the Prediction of Bacterial Hemoglobin-Like Proteins

The recent upsurge in microbial genome data has revealed that hemoglobin-like (HbL) proteins may be widely distributed among bacteria and that some organisms may carry more than one HbL encoding gene. However, the discovery of HbL proteins has been limited to a small number of bacteria only. This study describes the prediction of HbL proteins and their domain classification using a machine learning approach. Support vector machine (SVM) models were developed for predicting HbL proteins based upon amino acid composition (AC), dipeptide composition (DC), hybrid method (AC + DC), and position specific scoring matrix (PSSM). In addition, we introduce for the first time a new prediction method based on max to min amino acid residue (MM) profiles. The average accuracy, standard deviation (SD), false positive rate (FPR), confusion matrix, and receiver operating characteristic (ROC) were analyzed. We also compared the performance of our proposed models in homology detection databases. The performance of the different approaches was estimated using fivefold cross-validation techniques. Prediction accuracy was further investigated through confusion matrix and ROC curve analysis. All experimental results indicate that the proposed BacHbpred can be a perspective predictor for determination of HbL related proteins. BacHbpred, a web tool, has been developed for HbL prediction.

Molecular dynamics study of naturally existing cavity couplings in proteins

Couplings between protein sub-structures are a common property of protein dynamics. Some of these couplings are especially interesting since they relate to function and its regulation. In this article we have studied the case of cavity couplings because cavities can host functional sites, allosteric sites, and are the locus of interactions with the cell milieu. We have divided this problem into two parts. In the first part, we have explored the presence of cavity couplings in the natural dynamics of 75 proteins, using 20 ns molecular dynamics simulations. For each of these proteins, we have obtained two trajectories around their native state. After applying a stringent filtering procedure, we found significant cavity correlations in 60% of the proteins. We analyze and discuss the structure origins of these correlations, including neighbourhood, cavity distance, etc. In the second part of our study, we have used longer simulations (≥100 ns) from the MoDEL project, to obtain a broader view of cavity couplings, particularly about their dependence on time. Using moving window computations we explored the fluctuations of cavity couplings along time, finding that these couplings could fluctuate substantially during the trajectory, reaching in several cases correlations above 0.25/0.5. In summary, we describe the structural origin and the variations with time of cavity couplings. We complete our work with a brief discussion of the biological implications of these results.

Coordinated gene expression of neuroinflammatory and cell signaling markers in dorsolateral prefrontal cortex during human brain development and aging

BACKGROUND

Age changes in expression of inflammatory, synaptic, and neurotrophic genes are not well characterized during human brain development and senescence. Knowing these changes may elucidate structural, metabolic, and functional brain processes over the lifespan, as well vulnerability to neurodevelopmental or neurodegenerative diseases.

HYPOTHESIS

Expression levels of inflammatory, synaptic, and neurotrophic genes in the human brain are coordinated over the lifespan and underlie changes in phenotypic networks or cascades.

METHODS

We used a large-scale microarray dataset from human prefrontal cortex, BrainCloud, to quantify age changes over the lifespan, divided into Development (0 to 21 years, 87 brains) and Aging (22 to 78 years, 144 brains) intervals, in transcription levels of 39 genes.

RESULTS

Gene expression levels followed different trajectories over the lifespan. Many changes were intercorrelated within three similar groups or clusters of genes during both Development and Aging, despite different roles of the gene products in the two intervals. During Development, changes were related to reported neuronal loss, dendritic growth and pruning, and microglial events; TLR4, IL1R1, NFKB1, MOBP, PLA2G4A, and PTGS2 expression increased in the first years of life, while expression of synaptic genes GAP43 and DBN1 decreased, before reaching plateaus. During Aging, expression was upregulated for potentially pro-inflammatory genes such as NFKB1, TRAF6, TLR4, IL1R1, TSPO, and GFAP, but downregulated for neurotrophic and synaptic integrity genes such as BDNF, NGF, PDGFA, SYN, and DBN1.

CONCLUSIONS

Coordinated changes in gene transcription cascades underlie changes in synaptic, neurotrophic, and inflammatory phenotypic networks during brain Development and Aging. Early postnatal expression changes relate to neuronal, glial, and myelin growth and synaptic pruning events, while late Aging is associated with pro-inflammatory and synaptic loss changes. Thus, comparable transcriptional regulatory networks that operate throughout the lifespan underlie different phenotypic processes during Aging compared to Development.

Metaproteomics to unravel major microbial players in leaf litter and soil environments: challenges and perspectives

Soil- and litter-borne microorganisms vitally contribute to biogeochemical cycles. However, changes in environmental parameters but also human interferences may alter species composition and elicit alterations in microbial activities. Soil and litter metaproteomics, implying the assignment of soil and litter proteins to specific phylogenetic and functional groups, has a great potential to provide essential new insights into the impact of microbial diversity on soil ecosystem functioning. This article will illuminate challenges and perspectives of current soil and litter metaproteomics research, starting with an introduction to an appropriate experimental design and state-of-the-art proteomics methodologies. This will be followed by a summary of important studies aimed at (i) the discovery of the major biotic drivers of leaf litter decomposition, (ii) metaproteomics analyses of rhizosphere-inhabiting microbes, and (iii) global approaches to study bioremediation processes. The review will be closed by a brief outlook on future developments and some concluding remarks, which should assist the reader to develop successful concepts for soil and litter metaproteomics studies.

Combining bioinformatics and MS-based proteomics: clinical implications

Clinical proteomics research aims at i) discovery of protein biomarkers for screening, diagnosis and prognosis of disease, ii) discovery of protein therapeutic targets for improvement of disease prevention, treatment and follow-up, and iii) development of mass spectrometry (MS)-based assays that could be implemented in clinical chemistry, microbiology or hematology laboratories. MS has been increasingly applied in clinical proteomics studies for the identification and quantification of proteins. Bioinformatics plays a key role in the exploitation of MS data in several aspects such as the generation and curation of protein sequence databases, the development of appropriate software for MS data treatment and integration with other omics data and the establishment of adequate standard files for data sharing. In this article, we discuss the main MS approaches and bioinformatics solutions that are currently applied to accomplish the objectives of clinical proteomic research.

HeatMapViewer: interactive display of 2D data in biology

Summary: The HeatMapViewer is a BioJS component that lays-out and renders two-dimensional (2D) plots or heat maps that are ideally suited to visualize matrix formatted data in biology such as for the display of microarray experiments or the outcome of mutational studies and the study of SNP-like sequence variants. It can be easily integrated into documents and provides a powerful, interactive way to visualize heat maps in web applications. The software uses a scalable graphics technology that adapts the visualization component to any required resolution, a useful feature for a presentation with many different data-points. The component can be applied to present various biological data types. Here, we present two such cases – showing gene expression data and visualizing mutability landscape analysis.

Availability:https://github.com/biojs/biojs; http://dx.doi.org/10.5281/zenodo.7706.

An attempt to construct a (general) mathematical framework to model biological "context-dependence"

Context-dependent nature of biological phenomena is well documented in every branch of biology. While there have been few previous attempts to (implicitly) model various (particular) facets of biological context-dependence, a formal and general mathematical construct to model the wide spectrum of context-dependence, eludes the students of biology. Such an objective model, from both 'bottom-up' as well as 'top-down' perspective, is proposed here to serve as the template to describe the various kinds of context-dependence that we encounter in different branches of biology. Interactions between biological contexts was found to be transitive but non-commutative. It is found that a hierarchical nature of dependence among the biological contexts models the emergent biological properties efficiently. Reasons for these findings are provided in a general model to describe biological reality. Scheme to algorithmically implement the hierarchic structure of organization of biological contexts was proposed with a construct named 'Context tree'. A 'Context tree' based analysis of context interactions among biophysical factors influencing protein structure was performed.

Antibody and antigen contact residues define epitope and paratope size and structure

A total of 111 Ag-Ab x-ray crystal structures of large protein Ag epitopes and paratopes were analyzed to inform the process of eliciting or selecting functional and therapeutic Abs. These analyses illustrate that Ab contact residues (CR) are distributed in three prominent CR regions (CRR) on L and H chains that overlap but do not coincide with Ab CDR. The number of Ag and Ab CRs per structure are overlapping and centered around 18 and 19, respectively. The CR span (CRS), a novel measure introduced in this article, is defined as the minimum contiguous amino acid sequence containing all CRs of an Ag or Ab and represents the size of a complete structural epitope or paratope, inclusive of CR and the minimum set of supporting residues required for proper conformation. The most frequent size of epitope CRS is 50-79 aa, which is similar in size to L (60-69) and H chain (70-79) CRS. The size distribution of epitope CRS analyzed in this study ranges from ~20 to 400 aa, similar to the distribution of independent protein domain sizes reported in the literature. Together, the number of CRs and the size of the CRS demonstrate that, on average, complete structural epitopes and paratopes are equal in size to each other and similar in size to intact protein domains. Thus, independent protein domains inclusive of biologically relevant sites represent the fundamental structural unit bound by, and useful for eliciting or selecting, functional and therapeutic Abs.

EuLoc: a web-server for accurately predict protein subcellular localization in eukaryotes by incorporating various features of sequence segments into the general form of Chou's PseAAC

The function of a protein is generally related to its subcellular localization. Therefore, knowing its subcellular localization is helpful in understanding its potential functions and roles in biological processes. This work develops a hybrid method for computationally predicting the subcellular localization of eukaryotic protein. The method is called EuLoc and incorporates the Hidden Markov Model (HMM) method, homology search approach and the support vector machines (SVM) method by fusing several new features into Chou's pseudo-amino acid composition. The proposed SVM module overcomes the shortcoming of the homology search approach in predicting the subcellular localization of a protein which only finds low-homologous or non-homologous sequences in a protein subcellular localization annotated database. The proposed HMM modules overcome the shortcoming of SVM in predicting subcellular localizations using few data on protein sequences. Several features of a protein sequence are considered, including the sequence-based features, the biological features derived from PROSITE, NLSdb and Pfam, the post-transcriptional modification features and others. The overall accuracy and location accuracy of EuLoc are 90.5 and 91.2 %, respectively, revealing a better predictive performance than obtained elsewhere. Although the amounts of data of the various subcellular location groups in benchmark dataset differ markedly, the accuracies of 12 subcellular localizations of EuLoc range from 82.5 to 100 %, indicating that this tool is much more balanced than other tools. EuLoc offers a high, balanced predictive power for each subcellular localization. EuLoc is now available on the web at http://euloc.mbc.nctu.edu.tw/.

Modeling of three dimensional structure of human alpha-fetoprotein complexed with diethylstilbestrol: docking and molecular dynamics simulation study

It has been long experimentally demonstrated that human alpha-fetoprotein (HAFP) has an ability to bind immobilized estrogens with the most efficiency for synthetic estrogen analog - diethylstilbestrol (DES). However, the question remains why the human AFP (HAFP), unlike rodent AFP, cannot bind free estrogens. Moreover, despite the fact that AFP was first discovered more than 50 years ago and is presently recognized as a "golden standard" among onco-biomarkers, its three-dimensional (3D) structure has not been experimentally solved yet. In this work using MODELLER program, we generated 3D model of HAFP on the basis of homology with human serum albumin (HSA) and Vitamin D-binding protein (VTDB) with subsequent molecular docking of DES to the model structure and molecular dynamics (MD) simulation study of the complex obtained. The model constructed has U-shaped structure in which a cavity may be distinguished. In this cavity the putative estrogen-binding site is localized. Validation by RMSD calculation and with the use of PROCHECK program showed good quality of the model and stability of extended region of four alpha-helical structures that contains putative hormone-binding residues. Data extracted from MD simulation trajectory allow proposing two types of interactions between amino acid residues of HAFP and DES molecule: (1) hydrogen bonding with involvement of residues S445, R452, and E551; (2) hydrophobic interactions with participation of L138, M448, and M548 residues. A suggestion is made that immobilization of the hormone using a long spacer provides delivery of the estrogen molecule to the binding site and, thereby, facilitates interaction between HAFP and the hormone.

The three M's: melanoma, microphthalmia-associated transcription factor and microRNA

Studies examining intratumor heterogeneity have indicated that several cancer types, including melanoma, can display phenotypic plasticity, corresponding to their capacity to undergo transient reversible cellular changes. Conceptual models constructed to explain the process of cancer propagation differ in their treatment of intratumor heterogeneity. Recent observations of reversible phenotypic heterogeneity in melanoma have led to the proposal of a novel 'phenotypic plasticity' model of cancer propagation. Microphthalmia-associated transcription factor (MITF), the melanocyte 'lineage-specific' transcription factor, has emerged as one of the central players in melanoma phenotypic plasticity. Here we discuss the conceptual models suggested to explain the relations between MITF and melanoma plasticity, in addition to the complex regulatory roles that MITF plays in melanocytes and melanoma development. Finally, we provide an in-depth literature survey of microRNAs (miRNAs) involved in MITF activity, melanoma propagation and metastasis, in addition to their potential use as agents of personalized therapy.

Influence of intramolecular interactions on conformational and dynamic properties of analogs of heptapeptide AFP(14-20)

Conformational and dynamic properties of proteins and peptides play an important role in their functioning. However, mechanisms that underlie this influence have not been fully elucidated. In the present work we computationally constructed analogs of heptapeptide AFP(14-20) (LDSYQCT) - one of the biologically active sites of human α-fetoprotein (AFP) - to study their conformational and dynamic properties using molecular dynamics simulation. Analogs were obtained by point substitutions of amino acid residues taking into account differences in their physicochemical properties and also on the basis of analysis of amino acid substitutions in the AFP(14-20)-like motifs revealed in different physiologically active proteins. It is shown that changes in conformational mobility of amino acid residues of analogs are due to disruption or arising of intramolecular interactions that, in turn, determine existence of steric restrictions during rotation around covalent bonds of the peptide backbone. Substitution of an amino acid by another one with significant difference in physicochemical properties may not lead to remarkable changes in conformational and dynamic properties of the peptide if intramolecular interactions remain unchanged.

A bioinformatics workflow for variant peptide detection in shotgun proteomics

Shotgun proteomics data analysis usually relies on database search. However, commonly used protein sequence databases do not contain information on protein variants and thus prevent variant peptides and proteins from been identified. Including known coding variations into protein sequence databases could help alleviate this problem. Based on our recently published human Cancer Proteome Variation Database, we have created a protein sequence database that comprehensively annotates thousands of cancer-related coding variants collected in the Cancer Proteome Variation Database as well as noncancer-specific ones from the Single Nucleotide Polymorphism Database (dbSNP). Using this database, we then developed a data analysis workflow for variant peptide identification in shotgun proteomics. The high risk of false positive variant identifications was addressed by a modified false discovery rate estimation method. Analysis of colorectal cancer cell lines SW480, RKO, and HCT-116 revealed a total of 81 peptides that contain either noncancer-specific or cancer-related variations. Twenty-three out of 26 variants randomly selected from the 81 were confirmed by genomic sequencing. We further applied the workflow on data sets from three individual colorectal tumor specimens. A total of 204 distinct variant peptides were detected, and five carried known cancer-related mutations. Each individual showed a specific pattern of cancer-related mutations, suggesting potential use of this type of information for personalized medicine. Compatibility of the workflow has been tested with four popular database search engines including Sequest, Mascot, X!Tandem, and MyriMatch. In summary, we have developed a workflow that effectively uses existing genomic data to enable variant peptide detection in proteomics.

Solution structure of 4'-phosphopantetheine - GmACP3 from Geobacter metallireducens: a specialized acyl carrier protein with atypical structural features and a putative role in lipopolysaccharide biosynthesis

GmACP3 from Geobacter metallireducens is a specialized acyl carrier protein (ACP) whose gene, gmet_2339, is located near genes encoding many proteins involved in lipopolysaccharide (LPS) biosynthesis, indicating a likely function for GmACP3 in LPS production. By overexpression in Escherichia coli, about 50% holo-GmACP3 and 50% apo-GmACP3 were obtained. Apo-GmACP3 exhibited slow precipitation and non-monomeric behavior by (15)N NMR relaxation measurements. Addition of 4'-phosphopantetheine (4'-PP) via enzymatic conversion by E. coli holo-ACP synthase resulted in stable >95% holo-GmACP3 that was characterized as monomeric by (15)N relaxation measurements and had no indication of conformational exchange. We have determined a high-resolution solution structure of holo-GmACP3 by standard NMR methods, including refinement with two sets of NH residual dipolar couplings, allowing for a detailed structural analysis of the interactions between 4'-PP and GmACP3. Whereas the overall four helix bundle topology is similar to previously solved ACP structures, this structure has unique characteristics, including an ordered 4'-PP conformation that places the thiol at the entrance to a central hydrophobic cavity near a conserved hydrogen-bonded Trp-His pair. These residues are part of a conserved WDSLxH/N motif found in GmACP3 and its orthologs. The helix locations and the large hydrophobic cavity are more similar to medium- and long-chain acyl-ACPs than to other apo- and holo-ACP structures. Taken together, structural characterization along with bioinformatic analysis of nearby genes suggests that GmACP3 is involved in lipid A acylation, possibly by atypical long-chain hydroxy fatty acids, and potentially is involved in synthesis of secondary metabolites.

Structural investigation of myo-inositol dehydrogenase from Bacillus subtilis: implications for catalytic mechanism and inositol dehydrogenase subfamily classification

Inositol dehydrogenase from Bacillus subtilis (BsIDH) is a NAD+-dependent enzyme that catalyses the oxidation of the axial hydroxy group of myo-inositol to form scyllo-inosose. We have determined the crystal structures of wild-type BsIDH and of the inactive K97V mutant in apo-, holo- and ternary complexes with inositol and inosose. BsIDH is a tetramer, with a novel arrangement consisting of two long continuous β-sheets, formed from all four monomers, in which the two central strands are crossed over to form the core of the tetramer. Each subunit in the tetramer consists of two domains: an N-terminal Rossmann fold domain containing the cofactor-binding site, and a C-terminal domain containing the inositol-binding site. Structural analysis allowed us to determine residues important in cofactor and substrate binding. Lys97, Asp172 and His176 are the catalytic triad involved in the catalytic mechanism of BsIDH, similar to what has been proposed for related enzymes and short-chain dehydrogenases. Furthermore, a conformational change in the nicotinamide ring was observed in some ternary complexes, suggesting hydride transfer to the si-face of NAD+. Finally, comparison of the structure and sequence of BsIDH with other putative inositol dehydrogenases allowed us to differentiate these enzymes into four subfamilies based on six consensus sequence motifs defining the cofactor- and substrate-binding sites.