Proteomic analysis of coronary sinus serum reveals leucine-rich α2-glycoprotein as a novel biomarker of ventricular dysfunction and heart failure

BACKGROUND: Heart failure (HF) prevention strategies require biomarkers that identify disease manifestation. Increases in B-type natriuretic peptide (BNP) correlate with increased risk of cardiovascular events and HF development. We hypothesize that coronary sinus serum from a high BNP hypertensive population reflects an active pathological process and can be used for biomarker exploration. Our aim was to discover differentially expressed disease-associated proteins that identify patients with ventricular dysfunction and HF.

METHODS AND RESULTS: Coronary sinus serum from 11 asymptomatic, hypertensive patients underwent quantitative differential protein expression analysis by 2-dimensional difference gel electrophoresis. Proteins were identified using mass spectrometry and then studied by enzyme-linked immunosorbent assay in sera from 40 asymptomatic, hypertensive patients and 105 patients across the spectrum of ventricular dysfunction (32 asymptomatic left ventricular diastolic dysfunction, 26 diastolic HF, and 47 systolic HF patients). Leucine-rich α2-glycoprotein (LRG) was consistently overexpressed in high BNP serum. LRG levels correlate significantly with BNP in hypertensive, asymptomatic left ventricular diastolic dysfunction, diastolic HF, and systolic HF patient groups (P≤0.05). LRG levels were able to identify HF independent of BNP. LRG correlates with coronary sinus serum levels of tumor necrosis factor-α (P=0.009) and interleukin-6 (P=0.021). LRG is expressed in myocardial tissue and correlates with transforming growth factor-βR1 (P<0.001) and α-smooth muscle actin (P=0.025) expression.

CONCLUSIONS: LRG was identified as a serum biomarker that accurately identifies patients with HF. Multivariable modeling confirmed that LRG is a stronger identifier of HF than BNP and this is independent of age, sex, creatinine, ischemia, β-blocker therapy, and BNP.

Characterisation of a secretory serine protease inhibitor (SjB6) from Schistosoma japonicum

BACKGROUND: Proteins belonging to the serine protease inhibitor (serpin) superfamily play essential physiological roles in many organisms. In pathogens, serpins are thought to have evolved specifically to limit host immune responses by interfering with the host immune-stimulatory signals. Serpins are less well characterised in parasitic helminths, although some are thought to be involved in mechanisms associated with host immune modulation. In this study, we cloned and partially characterised a secretory serpin from Schistosoma japonicum termed SjB6, these findings provide the basis for possible functional roles.

METHODS: SjB6 gene was identified through database mining of our previously published microarray data, cloned and detailed sequence and structural analysis and comparative modelling carried out using various bioinformatics and proteomics tools. Gene transcriptional profiling was determined by real-time PCR and the expression of native protein determined by immunoblotting. An immunological profile of the recombinant protein produced in insect cells was determined by ELISA.

RESULTS: SjB6 contains an open reading frame of 1160 base pairs that encodes a protein of 387 amino acid residues. Detailed sequence analysis, comparative modelling and structural-based alignment revealed that SjB6 contains the essential structural motifs and consensus secondary structures typical of inhibitory serpins. The presence of an N-terminal signal sequence indicated that SjB6 is a secretory protein. Real-time data indicated that SjB6 is expressed exclusively in the intra-mammalian stage of the parasite life cycle with its highest expression levels in the egg stage (p < 0.0001). The native protein is approximately 60 kDa in size and recombinant SjB6 (rSjB6) was recognised strongly by sera from rats experimentally infected with S. japonicum.

CONCLUSIONS: The significantly high expression of SjB6 in schistosome eggs, when compared to other life cycle stages, suggests a possible association with disease pathology, while the strong reactivity of sera from experimentally infected rats against rSjB6 suggests that native SjB6 is released into host tissue and induces an immune response. This study presents a comprehensive demonstration of sequence and structural-based analysis of a secretory serpin from a trematode and suggests SjB6 may be associated with important functional roles in S. japonicum, particularly in parasite modulation of the host microenvironment.

Primary pulmonary myxoid sarcoma with EWSR1-CREB1 fusion: a new tumor entity.

We present clinicopathologic data on 10 pulmonary myxoid sarcomas, which are defined by distinctive histomorphologic features and characterized by a recurrent fusion gene, that appear to represent a distinct tumor entity at this site. The patients [7 female, 3 male; aged 27 to 67 y (mean, 45 y)] presented with local or systemic symptoms (n=5), symptoms from cerebral metastasis (1), or incidentally (2). Follow-up of 6 patients showed that 1 with brain metastasis died shortly after primary tumor resection, 1 developed a renal metastasis but is alive and well, and 4 are disease free after 1 to 15 years. All tumors involved pulmonary parenchyma, with a predominant endobronchial component in 8 and ranged from 1.5 to 4 cm. Microscopically, they were lobulated and composed of cords of polygonal, spindle, or stellate cells within myxoid stroma, morphologically reminiscent of extraskeletal myxoid chondrosarcoma. Four cases showed no or minimal atypia, 6 showed focal pleomorphism, and 5 had necrosis. Mitotic indices varied, with most tumors not exceeding 5/10 high-power fields. Tumors were immunoreactive for only vimentin and weakly focal for epithelial membrane antigen. Of 9 tumors, 7 were shown to harbor a specific EWSR1-CREB1 fusion by reverse transcription-polymerase chain reaction and direct sequencing, with 7 of 10 showing EWSR1 rearrangement by fluorescence in situ hybridization. This gene fusion has been described previously in 2 histologically and behaviorally different sarcomas: clear cell sarcoma-like tumors of the gastrointestinal tract and angiomatoid fibrous histiocytomas; however, this is a novel finding in tumors with the morphology we describe and that occur in the pulmonary region.

Gene mapping and expression analysis of 16q loss of heterozygosity identifies WWOX and CYLD as being important in determining clinical outcome in multiple myeloma.

We performed fluorescent in situ hybridization (FISH) for 16q23 abnormalities in 861 patients with newly diagnosed multiple myeloma and identified deletion of 16q [del(16q)] in 19.5%. In 467 cases in which demographic and survival data were available, del(16q) was associated with a worse overall survival (OS). It was an independent prognostic marker and conferred additional adverse survival impact in cases with the known poor-risk cytogenetic factors t(4;14) and del(17p). Gene expression profiling and gene mapping using 500K single-nucleotide polymorphism (SNP) mapping arrays revealed loss of heterozygosity (LOH) involving 3 regions: the whole of 16q, a region centered on 16q12 (the location of CYLD), and a region centered on 16q23 (the location of the WW domain-containing oxidoreductase gene WWOX). CYLD is a negative regulator of the NF-kappaB pathway, and cases with low expression of CYLD were used to define a "low-CYLD signature." Cases with 16q LOH or t(14;16) had significantly reduced WWOX expression. WWOX, the site of the translocation breakpoint in t(14;16) cases, is a known tumor suppressor gene involved in apoptosis, and we were able to generate a "low-WWOX signature" defined by WWOX expression. These 2 genes and their corresponding pathways provide an important insight into the potential mechanisms by which 16q LOH confers poor prognosis.

Insights into the multistep transformation of MGUS to myeloma using microarray expression analysis.

To define specific pathways important in the multistep transformation process of normal plasma cells (PCs) to monoclonal gammopathy of uncertain significance (MGUS) and multiple myeloma (MM), we have applied microarray analysis to PCs from 5 healthy donors (N), 7 patients with MGUS, and 24 patients with newly diagnosed MM. Unsupervised hierarchical clustering using 125 genes with a large variation across all samples defined 2 groups: N and MGUS/MM. Supervised analysis identified 263 genes differentially expressed between N and MGUS and 380 genes differentially expressed between N and MM, 197 of which were also differentially regulated between N and MGUS. Only 74 genes were differentially expressed between MGUS and MM samples, indicating that the differences between MGUS and MM are smaller than those between N and MM or N and MGUS. Differentially expressed genes included oncogenes/tumor-suppressor genes (LAF4, RB1, and disabled homolog 2), cell-signaling genes (RAS family members, B-cell signaling and NF-kappaB genes), DNA-binding and transcription-factor genes (XBP1, zinc finger proteins, forkhead box, and ring finger proteins), and developmental genes (WNT and SHH pathways). Understanding the molecular pathogenesis of MM by gene expression profiling has demonstrated sequential genetic changes from N to malignant PCs and highlighted important pathways involved in the transformation of MGUS to MM.

Salicaceae Endophytes: Growth Promotion Potential in Rice (Oryza sativa L.) and Maize (Zea mays L.) and Bio-Control of Plant Pathogen

Thesis (Ph.D.)--University of Washington, 2016-08

Whole Genome Analysis of Gene Expression Reveals Coordinated Activation of Signaling and Metabolic Pathways during Pollen-Pistil Interactions in Arabidopsis

Plant reproduction depends on the concerted activation of many genes to ensure correct communication between pollen and pistil. Here, we queried the whole transcriptome of Arabidopsis (Arabidopsis thaliana) in order to identify genes with specific reproductive functions. We used the Affymetrix ATH1 whole genome array to profile wild-type unpollinated pistils and unfertilized ovules. By comparing the expression profile of pistils at 0.5, 3.5, and 8.0 h after pollination and applying a number of statistical and bioinformatics criteria, we found 1,373 genes differentially regulated during pollen-pistil interactions. Robust clustering analysis grouped these genes in 16 time-course clusters representing distinct patterns of regulation. Coregulation within each cluster suggests the presence of distinct genetic pathways, which might be under the control of specific transcriptional regulators. A total of 78% of the regulated genes were expressed initially in unpollinated pistil and/or ovules, 15% were initially detected in the pollen data sets as enriched or preferentially expressed, and 7% were induced upon pollination. Among those, we found a particular enrichment for unknown transcripts predicted to encode secreted proteins or representing signaling and cell wall-related proteins, which may function by remodeling the extracellular matrix or as extracellular signaling molecules. A strict regulatory control in various metabolic pathways suggests that fine-tuning of the biochemical and physiological cellular environment is crucial for reproductive success. Our study provides a unique and detailed temporal and spatial gene expression profile of in vivo pollen-pistil interactions, providing a framework to better understand the basis of the molecular mechanisms operating during the reproductive process in higher plants.

Molecular biology of prolactin and prolactin receptor in a marine teleost: the sea bream(Sparus aurata)

Tese de dout. em Biologia, especialidade de Biologia Molecular, Unidade de Ciências e Tecnologias dos Recursos Aquáticos, Univ. do Algarve

Protein expression of Myt272-3 recombinant clone and in silico prediction of a possible vaccine candidate against Mycobacterium tuberculosis

Purpose: To investigate the expression of Myt272-3 recombinant protein and also to predict a possible protein vaccine candidate against Mycobacterium tuberculosis . Methods: Myt272-3 protein was expressed in pET30a+-Myt272-3 clone. The purity of the protein was determined using Dynabeads® His-Tag Isolation & Pulldown. Protein sequence was analysed in silico using bioinformatics software for the prediction of allergenicity, antigenicity, MHC-I and MHC-II binding, and B-cell epitope binding. Results: The candidate protein was a non-allergen with 15.19 % positive predictive value. It was also predicted to be antigenic, with binding affinity to MHC-I and MHC-II, as well as B-cell epitope binding. Conclusion: The predicted results obtained in this study provide a guide for practical design of a new tuberculosis vaccine.

Integrated Omics for the global mapping of biomolecules in LDL

Circulating low density lipoproteins (LDL) are thought to play a crucial role in the onset and development of atherosclerosis, though the detailed molecular mechanisms responsible for their biological effects remain controversial. The complexity of biomolecules (lipids, glycans and protein) and structural features (isoforms and chemical modifications) found in LDL particles hampers the complete understanding of the mechanism underlying its atherogenicity. For this reason the screening of LDL for features discriminative of a particular pathology in search of biomarkers is of high importance. Three major biomolecule classes (lipids, protein and glycans) in LDL particles were screened using mass spectrometry coupled to liquid chromatography. Dual-polarity screening resulted in good lipidome coverage, identifying over 300 lipid species from 12 lipid sub-classes. Multivariate analysis was used to investigate potential discriminators in the individual lipid sub-classes for different study groups (age, gender, pathology). Additionally, the high protein sequence coverage of ApoB-100 routinely achieved (≥70%) assisted in the search for protein modifications correlating to aging and pathology. The large size and complexity of the datasets required the use of chemometric methods (Partial Least Square-Discriminant Analysis, PLS-DA) for their analysis and for the identification of ions that discriminate between study groups. The peptide profile from enzymatically digested ApoB-100 can be correlated with the high structural complexity of lipids associated with ApoB-100 using exploratory data analysis. In addition, using targeted scanning modes, glycosylation sites within neutral and acidic sugar residues in ApoB-100 are also being explored. Together or individually, knowledge of the profiles and modifications of the major biomolecules in LDL particles will contribute towards an in-depth understanding, will help to map the structural features that contribute to the atherogenicity of LDL, and may allow identification of reliable, pathology-specific biomarkers. This research was supported by a Marie Curie Intra-European Fellowship within the 7th European Community Framework Program (IEF 255076). Work of A. Rudnitskaya was supported by Portuguese Science and Technology Foundation, through the European Social Fund (ESF) and "Programa Operacional Potencial Humano - POPH".

Comparative Genomics of Field Isolates of Mycobacterium bovis and M. caprae Provides Evidence for Possible Correlates with Bacterial Viability and Virulence

Mycobacteria of the Mycobacterium tuberculosis complex (MTBC) greatly affect humans and animals worldwide. The life cycle of mycobacteria is complex and the mechanisms resulting in pathogen infection and survival in host cells are not fully understood. Recently, comparative genomics analyses have provided new insights into the evolution and adaptation of the MTBC to survive inside the host. However, most of this information has been obtained using M. tuberculosis but not other members of the MTBC such as M. bovis and M. caprae. In this study, the genome of three M. bovis (MB1, MB3, MB4) and one M. caprae (MB2) field isolates with different lesion score, prevalence and host distribution phenotypes were sequenced. Genome sequence information was used for whole-genome and protein-targeted comparative genomics analysis with the aim of finding correlates with phenotypic variation with potential implications for tuberculosis (TB) disease risk assessment and control. At the whole-genome level the results of the first comparative genomics study of field isolates of M. bovis including M. caprae showed that as previously reported for M. tuberculosis, sequential chromosomal nucleotide substitutions were the main driver of the M. bovis genome evolution. The phylogenetic analysis provided a strong support for the M. bovis/M. caprae clade, but supported M. caprae as a separate species. The comparison of the MB1 and MB4 isolates revealed differences in genome sequence, including gene families that are important for bacterial infection and transmission, thus highlighting differences with functional implications between isolates otherwise classified with the same spoligotype. Strategic protein-targeted analysis using the ESX or type VII secretion system, proteins linking stress response with lipid metabolism, host T cell epitopes of mycobacteria, antigens and peptidoglycan assembly protein identified new genetic markers and candidate vaccine antigens that warrant further study to develop tools to evaluate risks for TB disease caused by M. bovis/M.caprae and for TB control in humans and animals.

Characterization of plasmids in a human clinical strain of Lactococcus garvieae

The present work describes the molecular characterization of five circular plasmids found in the human clinical strain Lactococcus garvieae 21881. The plasmids were designated pGL1-pGL5, with molecular sizes of 4,536 bp, 4,572 bp, 12,948 bp, 14,006 bp and 68,798 bp, respectively. Based on detailed sequence analysis, some of these plasmids appear to be mosaics composed of DNA obtained by modular exchange between different species of lactic acid bacteria. Based on sequence data and the derived presence of certain genes and proteins, the plasmid pGL2 appears to replicate via a rolling-circle mechanism, while the other four plasmids appear to belong to the group of lactococcal theta-type replicons. The plasmids pGL1, pGL2 and pGL5 encode putative proteins related with bacteriocin synthesis and bacteriocin secretion and immunity. The plasmid pGL5 harbors genes (txn, orf5 and orf25) encoding proteins that could be considered putative virulence factors. The gene txn encodes a protein with an enzymatic domain corresponding to the family actin-ADP-ribosyltransferases toxins, which are known to play a key role in pathogenesis of a variety of bacterial pathogens. The genes orf5 and orf25 encode two putative surface proteins containing the cell wall-sorting motif LPXTG, with mucin-binding and collagen-binding protein domains, respectively. These proteins could be involved in the adherence of L. garvieae to mucus from the intestine, facilitating further interaction with intestinal epithelial cells and to collagenous tissues such as the collagen-rich heart valves. To our knowledge, this is the first report on the characterization of plasmids in a human clinical strain of this pathogen.

The impact and control of malignant catarrhal fever in Tanzania

Malignant Catarrhal Fever (MCF), an often-lethal infectious disease, presents as a variable complex of lesions in susceptible ungulate species. The disease is caused by a -herpesvirus following transmission from an inapparent carrier host. Two major epidemiological forms exist: wildebeest-associated MCF (WA-MCF), in which the virus is transmitted to susceptible species by wildebeest calves less than approximately four months of age, and sheepassociated MCF (SA-MCF) in which the virus is spread by sheep (primarily adolescents). Due to the lack of an in-vitro propagation system for the causative agent of the more economically significant SA-MCF, and with the expectation that cross-protective immunity may be provided, vaccine development has focused on the more easily propagated alcelaphine herpesvirus-1 (AlHV-1) that causes WA-MCF. In 2008 a direct viral challenge trial showed that a novel vaccine, employing an attenuated AlHV-1 (atAlHV-1) `C5000 virus strain, protected British Friesian-Holstein (FH) cattle against an intranasal challenge with virulent AlHV-1 `C5000 virus. For cattle keeping people living near wildebeest calving areas in sub-Saharan Africa an effective vaccine would have value as it would release them from the costly annual disease avoidance strategy of having to move their herds away from the oncoming wildebeest. On the other hand, an effective vaccine will release herd owners from the need to avoid MCF, allowing them to graze their cattle alongside wildebeest on the highly nutritious pastures of the calving areas. As such conservationists have raised concerns that the development of a vaccine might lead to detrimental grazing competition. The principle objective of this study was to test the novel vaccine on Tanzanian shorthorn zebu cross cattle (SZC).We did this firstly using a natural challenge field trial (Chapter Two) which demonstrated that immunisation with the atAlHV-1 vaccine was well tolerated and induced an oro-nasopharyngeal AlHV-1-specific and -neutralising antibody response. This resulted in an immunity in SZC cattle that was partially protective and reduced naturally transmitted infection by 56%. We also demonstrated that non-fatal infections occurred with a much higher frequency than previously thought. Because the calculated efficacy of the vaccine was less than that seen in British FH cattle we wanted to determine whether host factors, particular to SZC cattle, had impacted the outcomes of the field trial. To do this we repeated the 2008 direct viral challenge trial using SZC cattle (Chapter Four). During this trial we also investigated whether the recombinant bacterial flagellin monomer (FliC), when used as an adjuvant, might improve the vaccine’s efficacy. The findings from this trial indicated that direct challenge with pathogenic AlHV-1 is effective at inducing MCF in SZC cattle and that FliC is not an appropriate adjuvant for this vaccine. Furthermore, with less control group cattle dying of MCF than expected we speculate that SZC cattle may have a degree of resistance to MCF that affords them protection from infection and developing fatal disease. In Chapter Three we investigated aspects of the epidemiology of MCF, specifically whether wildebeest placenta, long implicated by Maasai cattle owners as a source of MCF, might play a role in viral transmission. Additionally, through comparative sequence analysis, at two specific genes (A9.5 and ORF50) of wild-type and atAlHV-1, we investigated whether the `C5000 strain, the source of which was taken from Africa more than 40 years ago, was appropriate for vaccine development. The detection of AlHV-1 virus in approximately 50% of placentae indicated that infection can occur in-utero and that this tissue might play a role in disease transmission. And, despite describing three new alleles of the A9.5 gene (supporting previous evidence that this gene is polymorphic and encodes a secretory protein with interleukin-4 as the major homologue), the observation that the most frequently detected haplotypes, in both wild-type and attenuated AlHV-1, were identical suggests that AlHV-1 has a slow molecular clock and that the attenuated strain was appropriate for vaccine development. In Chapter Five we present the first quantitative assessment of the annual MCF avoidance costs that Maasai pastoralists incur. In particular we estimated that as a result of MCF avoidance 64% of the total daily milk yield during the MCF season was not available to be used by the 81% of the family unit remaining at the permanent boma. This represents an upper-bound loss of approximately 8% of a household0s annual income. Despite these considerable losses we concluded that, given an incidence of fatal MCF in cattle living in wildebeest calving areas of 5% to 10%, if herd owners were to stop trying to avoid MCF by allowing their cattle to graze alongside wildebeest, any gains made through increased availability of milk, improved body condition and reduced energy demands would be offset by an increase in MCF-incidence. With the development of an effective vaccine, however, this alternative strategy might become optimal. The overall conclusion we draw therefore is that, despite the substantial costs incurred each year avoiding MCF, the partial protection afforded by the novel vaccine strategy is not sufficient to warrant a wholesale change in disease avoidance strategy. Nonetheless, even the partial protection provided by this vaccine could be of value to protect animals that cannot be moved, for example where some of the herd remain at the boma to provide milk or where land-use changes make traditional disease avoidance difficult. Furthermore, the vaccine may offer a feasible solution to some of the current land-use challenges and conflicts, providing a degree of protection to valuable livestock where avoidance strategies are not possible, but with less risk of precipitating the potentially damaging environmental consequences, such as overgrazing of highly nutritious seasonal pastures, that might result if herd owners decide they no longer need to avoid wildebeest.

A genome editing approach to the study of Parvovirus B19

Parvovirus B19 (B19V) is a ssDNA virus, with a 5596 nt long genome encapsidated within an icosahedral capsid with a diameter of 22 nm. Viral proteins are subdivided into structural and non-structural: the main non-structural one is the NS1, while the 2 structural proteins VP1 and VP2 assemble originating the capsid shell. B19V tropism is mainly limited to erythroid progenitor cells (EPCs), however, virus can be detected in several districts persisting in tissues possibly lifelong. The virus can induce anemia and erythroid aplasia. Therapeutic strategies are only symptomatic, so the search for antivirals is strongly active, with screenings showing the activity in vitro of different compounds like hydroxyurea, cidofovir and brincidofovir. In the first project, a functional minigenome of B19V was developed, able to express only the NS1 protein. This minigenome proved able to replicate and express the NS1 at levels comparable to unmodified clones. Furthermore, the ability of this minigenome to complement the function of NS1-deficient genomes was demonstrated, thus providing a proof-of-concept of B19V genome editing possibility and, at the same time, a useful tool to study the NS1 protein also as an antiviral target. In the second project I addressed the interplay between B19V and the cellular restriction factor APOBEC3B (A3B), a cytidine deaminase acting on ssDNA, whose footprint on B19V genome was proved by a bioinformatic sequence analysis performed by the hosting lab. To understand whether A3B still exerts activity and a potential antiviral effect on B19V, the UT7/EpoS1 cells were transduced with lentiviral vectors to silence A3B expression, then used as a model to study viral behavior. No significant role of A3B on B19V was demonstrated, in agreement with the hypothesis of viral adaptation to this cellular restriction factor; anyway, virus ability to alter A3B expression would deserve further investigations.

Deep learning and embeddings for problems of computational biology

The development of Next Generation Sequencing promotes Biology in the Big Data era. The ever-increasing gap between proteins with known sequences and those with a complete functional annotation requires computational methods for automatic structure and functional annotation. My research has been focusing on proteins and led so far to the development of three novel tools, DeepREx, E-SNPs&GO and ISPRED-SEQ, based on Machine and Deep Learning approaches. DeepREx computes the solvent exposure of residues in a protein chain. This problem is relevant for the definition of structural constraints regarding the possible folding of the protein. DeepREx exploits Long Short-Term Memory layers to capture residue-level interactions between positions distant in the sequence, achieving state-of-the-art performances. With DeepRex, I conducted a large-scale analysis investigating the relationship between solvent exposure of a residue and its probability to be pathogenic upon mutation. E-SNPs&GO predicts the pathogenicity of a Single Residue Variation. Variations occurring on a protein sequence can have different effects, possibly leading to the onset of diseases. E-SNPs&GO exploits protein embeddings generated by two novel Protein Language Models (PLMs), as well as a new way of representing functional information coming from the Gene Ontology. The method achieves state-of-the-art performances and is extremely time-efficient when compared to traditional approaches. ISPRED-SEQ predicts the presence of Protein-Protein Interaction sites in a protein sequence. Knowing how a protein interacts with other molecules is crucial for accurate functional characterization. ISPRED-SEQ exploits a convolutional layer to parse local context after embedding the protein sequence with two novel PLMs, greatly surpassing the current state-of-the-art. All methods are published in international journals and are available as user-friendly web servers. They have been developed keeping in mind standard guidelines for FAIRness (FAIR: Findable, Accessible, Interoperable, Reusable) and are integrated into the public collection of tools provided by ELIXIR, the European infrastructure for Bioinformatics.