The role SWISS-PROT and TrEMBL play in the genome research environment

Transcriptome sequencing of Cocos nucifera leaves in response to Rhynchophorus ferrugineus infestation

Red palm weevil (RPW, Rhynchophorus ferrugineus) is an invasive pest of palms. In China, coconut (Cocos nucifera) production is being significantly affected by the RPW attack. To develop a long-term RPW control strategy, host-plant resistance is the most sustainable option. In this regard, the availability of transcriptome sequencing data from RPW-infected coconut plants can be highly useful. Therefore, the present study assessed coconut leaf physiological responses and transcriptional changes after different days of RPW attack i.e., 5, 10, 15, 20, and 25 days after infestation (DAI). A comparison of physiological data indicated that populations with the higher number of RPW insects i.e., population C (15 males +21 females) and D (20 males +28 females) triggered higher antioxidant enzyme activities. We used this data to study the transcriptomic responses on 5 and 20 DAI. Of the 38,432 detected transcripts, 3,984, 1,981, 3,925, and 2,257 were differentially expressed in CK (control/no RPW)_vs._C (5 DAI), CK_vs._D (5 DAI), CK_vs._C (20 DAI), and CK_vs._D (20 DAI), respectively. These transcripts were enriched in plant-pathogen interaction, phenylpropanoid/flavonoid biosynthesis, amino sugar and nucleotide sugar metabolism, plant hormone signal transduction, mitogen-activated protein kinase, and reactive oxygen scavenging pathway. We discuss these results and present several candidate genes to be manipulated for developing a sustainable strategy to control RPW attack regarding host-plant resistance. Furthermore, these findings provide a basis for developing effective early and late RPW attack detection strategies.

Combined Transcriptome and Metabolome Profiling Provide Insights into Cold Responses in Rapeseed (Brassica napus L.) Genotypes with Contrasting Cold-Stress Sensitivity

Low temperature is a major environmental factor, which limits rapeseed (Brassica napus L.) growth, development, and productivity. So far, the physiological and molecular mechanisms of rapeseed responses to cold stress are not fully understood. Here, we explored the transcriptome and metabolome profiles of two rapeseed genotypes with contrasting cold responses, i.e., XY15 (cold-sensitive) and GX74 (cold-tolerant). The global metabolome profiling detected 545 metabolites in siliques of both genotypes before (CK) and after cold-stress treatment (LW). The contents of several sugar metabolites were affected by cold stress with the most accumulated saccharides being 3-dehydro-L-threonic acid, D-xylonic acid, inositol, D-mannose, D-fructose, D-glucose, and L-glucose. A total of 1943 and 5239 differentially expressed genes were identified from the transcriptome sequencing in XY15CK_vs_XY15LW and GX74CK_vs_GX74LW, respectively. We observed that genes enriched in sugar metabolism and biosynthesis-related pathways, photosynthesis, reactive oxygen species scavenging, phytohormone, and MAPK signaling were highly expressed in GX74LW. In addition, several genes associated with cold-tolerance-related pathways, e.g., the CBF-COR pathway and MAPK signaling, were specifically expressed in GX74LW. Contrarily, genes in the above-mentioned pathways were mostly downregulated in XY15LW. Thus, our results indicate the involvement of these pathways in the differential cold-stress responses in XY15 and GX74.

Genome-Wide DNA Methylation and Hydroxymethylation Changes Revealed Epigenetic Regulation of Neuromodulation and Myelination in Yak Hypothalamus

Both 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC) are important epigenetic modifications in neurodevelopment. However, there is little research examining the genome-wide patterns of 5mC and 5hmC in brain regions of animals under natural high-altitude conditions. We used oxidative reduced representation bisulfite sequencing (oxRRBS) to determine the 5mC and 5hmC sites in the brain, brainstem, cerebellum, and hypothalamus of yak and cattle. We reported the first map of genome-wide DNA methylation and hydroxymethylation in the brain, brainstem, cerebellum, and hypothalamus of yak (living at high altitudes) and cattle. Overall, we found striking differences in 5mC and 5hmC between the hypothalamus and other brain regions in both yak and cattle. Genome-wide profiling revealed that 5mC level decreased and 5hmC level increased in the hypothalamus than in other regions. Furthermore, we identified differentially methylated regions (DMRs) and differentially hydroxymethylated regions (DhMRs), most of which overlapped with each other. Interestingly, transcriptome results for these brain regions also showed distinctive gene levels in the hypothalamus. Finally, differentially expressed genes (DEGs) regulated by DMRs and DhMRs may play important roles in neuromodulation and myelination. Overall, our results suggested that mediation of 5mC and 5hmC on epigenetic regulation may broadly impact the development of hypothalamus and its biological functions.

Methylome and transcriptome profiles in three yak tissues revealed that DNA methylation and the transcription factor ZGPAT co-regulate milk production

Background

Domestic yaks play an indispensable role in sustaining the livelihood of Tibetans and other ethnic groups on the Qinghai-Tibetan Plateau (QTP), by providing milk and meat. They have evolved numerous physiological adaptations to high-altitude environment, including strong blood oxygen transportation capabilities and high metabolism. The roles of DNA methylation and gene expression in milk production and high-altitudes adaptation need further exploration.

Results

We performed genome-wide DNA methylome and transcriptome analyses of breast, lung, and biceps brachii muscle tissues from yaks of different ages. We identified 432,350 differentially methylated regions (DMRs) across the age groups within each tissue. The post-mature breast tissue had considerably more differentially methylated regions (155,957) than that from the three younger age groups. Hypomethylated genes with high expression levels might regulate milk production by influencing protein processing in the endoplasmic reticulum. According to weighted gene correlation network analysis, the “hub” gene ZGPAT was highly expressed in the post-mature breast tissue, indicating that it potentially regulates the transcription of 280 genes that influence protein synthesis, processing, and secretion. The tissue network analysis indicated that high expression of HIF1A regulates energy metabolism in the lung.

Conclusions

This study provides a basis for understanding the epigenetic mechanisms underlying milk production in yaks, and the results offer insight to breeding programs aimed at improving milk production.

Supplementary information

Supplementary information accompanies this paper at 10.1186/s12864-020-07151-3.

Dissemination of Genetic Acquisition/Loss Provides a Variety of Quorum Sensing Regulatory Properties in Pseudoalteromonas

A bstract: Quorum sensing (QS) enables single-celled bacteria to communicate with chemical signals in order to synchronize group-level bacterial behavior. Pseudoalteromonas are marine bacteria found in versatile environments, of which QS regulation for their habitat adaptation is extremely fragmentary. To distinguish genes required for QS regulation in Pseudoalteromonas, comparative genomics was deployed to define the pan-genomics for twelve isolates and previously-sequenced genomes, of which acyl-homoserine lactone (AHL)-based QS traits were characterized. Additionally, transposon mutagenesis was used to identify the essential QS regulatory genes in the selected Pseudoalteromonas isolate. A remarkable feature showed that AHL-based colorization intensity of biosensors induced by Pseudoalteromonas most likely correlates with QS regulators genetic heterogeneity within the genus. This is supported by the relative expression levels of two of the main QS regulatory genes (luxO and rpoN) analyzed in representative Pseudoalteromonas isolates. Notably, comprehensive QS regulatory schema and the working model proposed in Pseudoalteromonas seem to phylogenetically include the network architectures derived from Escherichia coli, Pseudomonas, and Vibrio. Several associated genes were mapped by transposon mutagenesis. Among them, a right origin-binding protein-encoding gene (robp) was functionally identified as a positive QS regulatory gene. This gene lies on a genomic instable region and exists in the aforementioned bioinformatically recruited QS regulatory schema. The obtained data emphasize that the distinctly- and hierarchically-organized mechanisms probably target QS association in Pseudoalteromonas dynamic genomes, thus leading to bacterial ability to accommodate their adaption fitness and survival advantages.

Transcriptome analysis of orange-spotted grouper (Epinephelus coioides) spleen in response to Singapore grouper iridovirus

Background

Orange-spotted grouper (Epinephelus coioides) is an economically important marine fish cultured in China and Southeast Asian countries. The emergence of infectious viral diseases, including iridovirus and betanodavirus, have severely affected food products based on this species, causing heavy economic losses. Limited available information on the genomics of E. coioides has hampered the understanding of the molecular mechanisms that underlie host-virus interactions. In this study, we used a 454 pyrosequencing method to investigate differentially-expressed genes in the spleen of the E. coioides infected with Singapore grouper iridovirus (SGIV).

Results

Using 454 pyrosequencing, we obtained abundant high-quality ESTs from two spleen-complementary DNA libraries which were constructed from SGIV-infected (V) and PBS-injected fish (used as a control: C). A total of 407,027 and 421,141 ESTs were produced in control and SGIV infected libraries, respectively. Among the assembled ESTs, 9,616 (C) and 10,426 (V) ESTs were successfully matched against known genes in the NCBI non-redundant (nr) database with a cut-off E-value above 10^-5. Gene ontology (GO) analysis indicated that "cell part", "cellular process" and "binding" represented the largest category. Among the 25 clusters of orthologous group (COG) categories, the cluster for "translation, ribosomal structure and biogenesis" represented the largest group in the control (185 ESTs) and infected (172 ESTs) libraries. Further KEGG analysis revealed that pathways, including cellular metabolism and intracellular immune signaling, existed in the control and infected libraries. Comparative expression analysis indicated that certain genes associated with mitogen-activated protein kinase (MAPK), chemokine, toll-like receptor and RIG-I signaling pathway were alternated in response to SGIV infection. Moreover, changes in the pattern of gene expression were validated by qRT-PCR, including cytokines, cytokine receptors, and transcription factors, apoptosis-associated genes, and interferon related genes.

Conclusion

This study provided abundant ESTs that could contribute greatly to disclosing novel genes in marine fish. Furthermore, the alterations of predicted gene expression patterns reflected possible responses of these fish to the virus infection. Taken together, our data not only provided new information for identification of novel genes from marine vertebrates, but also shed new light on the understanding of defense mechanisms of marine fish to viral pathogens.

Beyond annotation transfer by homology: novel protein-function prediction methods to assist drug discovery

Every entirely sequenced genome reveals 100 s to 1000 s of protein sequences for which the only annotation available is 'hypothetical protein'. Thus, in the human genome and in the genomes of pathogenic agents there could be 1000 s of potential, unexplored drug targets. Computational prediction of protein function can play a role in studying these targets. We shall review the challenges, research approaches and recently developed tools in the field of computational function-prediction and we will discuss the ways these issues can change the process of drug discovery.

CAMPO, SCR_FIND and CHC_FIND: a suite of web tools for computational structural biology

The identification of evolutionarily conserved features of protein structures can provide insights into their functional and structural properties. Three methods have been developed and implemented as WWW tools, CAMPO, SCR_FIND and CHC_FIND, to analyze evolutionarily conserved residues (ECRs), structurally conserved regions (SCRs) and conserved hydrophobic contacts (CHCs) in protein families and superfamilies, on the basis of their 3D structures and the homologous sequences available. The programs identify protein segments that conserve a similar main-chain conformation, compute residue-to-residue hydrophobic contacts involving only apolar atoms common to all the 3D structures analyzed and allow the identification of conserved amino-acid sites among protein structures and their homologous sequences. The programs also allow the visualization of SCRs, CHCs and ECRs directly on the superposed structures and their multiple structural and sequence alignments. Tools and tutorials explaining their usage are available at , and .

Role of MHC II affinity and molecular mimicry in defining anti-HER-2/neu MAb-3 linear peptide epitope

With the aid of computational biology, we have studied the possibility of predicting the peptides able to evoke humoral immune response by using as experimental model the human HER-2/neu breast cancer-associated antigen. We already demonstrated that HER-2/neu peptides, that are the target of humoral human and mouse immune responses, correspond to those sequences having a low degree of sequence similarity to host's proteome. Here we report that the linear peptide determinant of the anti-HER-2/neu MAb-3 is characterized by a low degree of sequence similarity to mouse proteome in combination with high binding potential to specific MHC II molecule.

Structural plasticity of thermophilic serine hydroxymethyltransferases

Serine hydroxymethyltransferase (SHMT) catalyzes the reversible cleavage of serine to form glycine and monocarbonic groups, essential in several biosynthetic pathways. The availability of crystallographic structures of SHMT from mesophilic organisms and information produced by the genomic projects prompted the analysis of the adaptation of SHMT to "extreme" environments, such as high temperatures, by exploitation of structural data from thermophilic organisms. The sequences of 10 thermophilic/hyperthermophilic SHMTs were multiply aligned to 53 mesophilic homologs and analyzed by a comparative approach, examining the amino acid compositions and preferred residue exchanges between mesophiles and extremophiles. The structural basis of the observed exchanges was further investigated through the application of homology modeling to the 10 extremophilic SHMTs. The results of this study indicate that, in SHMT, thermal stability can be achieved mainly through three strategies: (i) increased number of charged residues at the protein surface; (ii) increased hydrophobicity of the protein core; and (iii) substitution of thermolabile residues exposed to the solvent. Additional features of the archaeal SHMTs, for which no structural data are available yet, were also investigated to explain their quaternary assemblage and the interaction with modified folates.

Monoclonal and polyclonal humoral immune response to EC HER-2/NEU peptides with low similarity to the host's proteome

We are studying peptide immunogenicity as a function of the similarity level to the host's proteome. By using as a model the breast/prostate cancer-associated HER-2/neu antigen, we analyzed the monoclonal and polyclonal humoral immune responses against HER-2/neu peptide motifs not shared with the host proteome. We show here that (i) a mouse monoclonal antibody (MAb) raised against the extracellular domain (EC) of human HER-2/neu oncoprotein recognized a linear peptide motif endowed with low similarity level to the mouse proteome; (ii) likewise, human sera from breast/prostate cancer patients preferentially recognized peptide fragments from the EC of the HER-2/neu oncoprotein having sequences that are not present in the human proteome. Together with previous results obtained in other disease models (cervical cancer-associated HPV16 E7 oncoprotein and Pemphigus vulgaris auto-antigen desmoglein-3), the present data suggest that a low level of sequence similarity to the host's proteome might be an important factor in shaping the pool of B cell epitopes.

Individuation of monoclonal anti-HPV16 E7 antibody linear peptide epitope by computational biology

We applied computational biology to identify the linear amino acid sequence recognized by a mouse monoclonal antibody raised against the full length HPV16 E7 oncoprotein. Computer-assisted search for the epitopic peptide used two parameters: the capability of E7 peptides to bind to MHC class II molecules, and the similarity level of the oncoprotein sequence to the mouse proteome. We report that anti-E7 mAb recognized the peptide having both high binding potential to MHC II molecules and low level of molecular mimicry to mouse proteome. Peptide ability to bind to MHC II molecules appears a necessary but not sufficient condition to determine peptide immunodominance, by needing to be supported by a low degree of peptide similarity to the host's proteome.