Gonzalo Vera

XGAP: a uniform and extensible data model and software platform for genotype and phenotype experiments

We present an extensible software model for the genotype and phenotype community, XGAP. Readers can download a standard XGAP (http://www.xgap.org) or auto-generate a custom version using MOLGENIS with programming interfaces to R-software and web-services or user interfaces for biologists. XGAP has simple load formats for any type of genotype, epigenotype, transcript, protein, metabolite or other phenotype data. Current functionality includes tools ranging from eQTL analysis in mouse to genome-wide association studies in humans.

R/parallel – speeding up bioinformatics analysis with R

BackgroundR is the preferred tool for statistical analysis of many bioinformaticians due in part to the increasing number of freely available analytical methods. Such methods can be quickly reused and adapted to each particular experiment. However, in experiments where large amounts of data are generated, for example using high-throughput screening devices, the processing time required to analyze data is often quite long. A solution to reduce the processing time is the use of parallel computing technologies. Because R does not support parallel computations, several tools have been developed to enable such technologies. However, these tools require multiple modications to the way R programs are usually written or run. Although these tools can finally speed up the calculations, the time, skills and additional resources required to use them are an obstacle for most bioinformaticians.ResultsWe have designed and implemented an R add-on package, R/parallel, that extends R by adding user-friendly parallel computing capabilities. With R/parallel any bioinformatician can now easily automate the parallel execution of loops and benefit from the multicore processor power of todays desktop computers. Using a single and simple function, R/parallel can be integrated directly with other existing R packages. With no need to change the implemented algorithms, the processing time can be approximately reduced N-fold, N being the number of available processor cores.ConclusionR/parallel saves bioinformaticians time in their daily tasks of analyzing experimental data. It achieves this objective on two fronts: first, by reducing development time of parallel programs by avoiding reimplementation of existing methods and second, by reducing processing time by speeding up computations on current desktop computers. Future work is focused on extending the envelope of R/parallel by interconnecting and aggregating the power of several computers, both existing office computers and computing clusters.

miRNA Expression Profile Analysis in Kidney of Different Porcine Breeds

microRNAs (miRNAs) are important post-transcriptional regulators in eukaryotes that target mRNAs repressing their expression. The uncertain process of pig domestication, with different origin focuses, and the selection process that commercial breeds suffered, have generated a wide spectrum of breeds with clear genetic and phenotypic variability. The aim of this work was to define the miRNAs expression profile in kidney of several porcine breeds. Small RNA libraries from kidney were elaborated and high-throughput sequenced with the 454 Genome Sequencer FLX (Roche). Pigs used were classified into three groups: the European origin group (Iberian breed and European Wild Boar ancestor), European commercial breeds (Landrace, Large White and Piétrain breeds) and breeds with Asian origin (Meishan and Vietnamese breeds). A total of 229 miRNAs were described in the pig kidney miRNA profile, including 110 miRNAs out of the 257 previously described pig miRNAs and 119 orthologous miRNAs. The most expressed miRNAs in pig kidney microRNAome were Hsa-miR-200b-3p, Ssc-miR-125b and Ssc-miR-23b. Moreover, 5 novel porcine miRNAs and 3 orthologous miRNAs could be validated through RT-qPCR. miRNA sequence variation was determined in 116 miRNAs, evidencing the presence of isomiRs. 125 miRNAs were differentially expressed between breed groups. The identification of breed-specific miRNAs, which could be potentially associated to certain phenotypes, is becoming a new tool for the study of the genetic variability underlying complex traits and furthermore, it adds a new layer of complexity to the interesting process of pig evolution.

The role of viral and host microRNAs in the Aujeszky's disease virus during the infection process.

Porcine production is a primary market in the world economy. Controlling swine diseases in the farm is essential in order to achieve the sector necessities. Aujeszky’s disease is a viral condition affecting pigs and is endemic in many countries of the world, causing important economic losses in the swine industry. microRNAs (miRNAs) are non-coding RNAs which modulates gene expression in animals, plants and viruses. With the aim of understanding miRNA roles during the Aujeszky’s disease virus [ADV] (also known as suid herpesvirus type 1 [SuHV-1]) infection, the expression profiles of host and viral miRNAs were determined through deep sequencing in SuHV-1 infected porcine cell line (PK-15) and in an animal experimental SuHV-1 infection with virulent (NIA-3) and attenuated (Begonia) strains. In the in vivo approach miR-206, miR-133a, miR-133b and miR-378 presented differential expression between virus strains infection. In the in vitro approach, most miRNAs were down-regulated in infected groups. miR-92a and miR-92b-3p were up-regulated in Begonia infected samples. Functional analysis of all this over expressed miRNAs during the infection revealed their association in pathways related to viral infection processes and immune response. Furthermore, 8 viral miRNAs were detected by stem loop RT-qPCR in both in vitro and in vivo approaches, presenting a gene regulatory network affecting 59 viral genes. Most described viral miRNAs were related to Large Latency Transcript (LLT) and to viral transcription activators EP0 and IE180, and also to regulatory genes regarding their important roles in the host – pathogen interaction during viral infection.

affyGG: computational protocols for genetical genomics with Affymetrix arrays

MOTIVATION Affymetrix arrays use multiple probes per gene to measure mRNA abundances. Standard software takes averages over probes. Important information may be lost if polymorphisms in the mRNA affect the hybridization of individual probes. RESULTS We present custom software to analyze genetical genomics experiments in human, mouse and other organisms: (i) an R package providing functions for QTL analysis at the individual probe level and (ii) Perl scripts providing custom tracks in the UCSC Genome Browser to check for sequence polymorphisms in probe regions. AVAILABILITY http://gbic.biol.rug.nl/supplementary.

Identification of microRNAs in PCV2 subclinically infected pigs by high throughput sequencing

Porcine circovirus type 2 (PCV2) is the essential etiological infectious agent of PCV2-systemic disease and has been associated with other swine diseases, all of them collectively known as porcine circovirus diseases. MicroRNAs (miRNAs) are a new class of small non-coding RNAs that regulate gene expression post-transcriptionally. miRNAs play an increasing role in many biological processes. The study of miRNA-mediated host-pathogen interactions has emerged in the last decade due to the important role that miRNAs play in antiviral defense. The objective of this study was to identify the miRNA expression pattern in PCV2 subclinically infected and non-infected pigs. For this purpose an experimental PCV2 infection was carried out and small-RNA libraries were constructed from tonsil and mediastinal lymph node (MLN) of infected and non-infected pigs. High throughput sequencing determined differences in miRNA expression in MLN between infected and non-infected while, in tonsil, a very conserved pattern was observed. In MLN, miRNA 126-3p, miRNA 126-5p, let-7d-3p, mir-129a and mir-let-7b-3p were up-regulated whereas mir-193a-5p, mir-574-5p and mir-34a down-regulated. Prediction of functional analysis showed that these miRNAs can be involved in pathways related to immune system and in processes related to the pathogenesis of PCV2, although functional assays are needed to support these predictions. This is the first study on miRNA gene expression in pigs infected with PCV2 using a high throughput sequencing approach in which several host miRNAs were differentially expressed in response to PCV2 infection.

designGG: an R-package and web tool for the optimal design of genetical genomics experiments

BackgroundHigh-dimensional biomolecular profiling of genetically different individuals in one or more environmental conditions is an increasingly popular strategy for exploring the functioning of complex biological systems. The optimal design of such genetical genomics experiments in a cost-efficient and effective way is not trivial.ResultsThis paper presents designGG, an R package for designing optimal genetical genomics experiments. A web implementation for designGG is available at http://gbic.biol.rug.nl/designGG. All software, including source code and documentation, is freely available.ConclusionDesignGG allows users to intelligently select and allocate individuals to experimental units and conditions such as drug treatment. The user can maximize the power and resolution of detecting genetic, environmental and interaction effects in a genome-wide or local mode by giving more weight to genome regions of special interest, such as previously detected phenotypic quantitative trait loci. This will help to achieve high power and more accurate estimates of the effects of interesting factors, and thus yield a more reliable biological interpretation of data. DesignGG is applicable to linkage analysis of experimental crosses, e.g. recombinant inbred lines, as well as to association analysis of natural populations.

Evaluation of the capability of the PCV2 genome to encode miRNAs: lack of viral miRNA expression in an experimental infection

Porcine circovirus type 2 (PCV2) is a ssDNA virus causing PCV2-systemic disease (PCV2-SD), one of the most important diseases in swine. MicroRNAs (miRNAs) are a new class of small non-coding RNAs that regulate gene expression post-transcriptionally. Viral miRNAs have recently been described and the number of viral miRNAs has been increasing in the past few years. In this study, small RNA libraries were constructed from two tissues of subclinically PCV2 infected pigs to explore if PCV2 can encode viral miRNAs. The deep sequencing data revealed that PCV2 does not express miRNAs in an in vivo subclinical infection.

Differential expression of porcine microRNAs in African swine fever virus infected pigs: a proof-of-concept study

BackgroundAfrican swine fever (ASF) is a re-expanding devastating viral disease currently threatening the pig industry worldwide. MicroRNAs are a class of 17–25 nucleotide non- coding RNAs that have been shown to have critical functions in a wide variety of biological processes, such as cell differentiation, cell cycle regulation, carcinogenesis, apoptosis, regulation of immunity as well as in viral infections by cleavage or translational repression of mRNAs. Nevertheless, there is no information about miRNA expression in an ASFV infection.MethodsIn this proof-of-concept study, we have analyzed miRNAs expressed in spleen and submandibular lymph node of experimentally infected pigs with a virulent (E75) or its derived attenuated (E75CV1) ASFV strain, as well as, at different times post-infection with the virulent strain, by high throughput sequencing of small RNA libraries.ResultsSpleen presented a more differential expression pattern than lymph nodes in an ASFV infection. Of the most abundant miRNAs, 12 were differentially expressed in both tissues at two different times in infected animals with the virulent strain. Of these, miR-451, miR-145-5p, miR-181a and miR-122 presented up-regulation at late times post-infection while miR-92a, miR-23a, miR-92b-3p, miR-126-5p, miR-126-3p, miR-30d, miR-23b and miR-92c showed down-regulation. Of the 8 differentially expressed miRNAs identified at the same time post-infection in infected animals with the virulent strain compared with animals infected with its attenuated strain, miR-126-5p, miR-92c, miR-92a, miR-30e-5p and miR-500a-5p presented up-regulation whereas miR-125b, miR-451 and miR-125a were down-regulated. All these miRNAs have been shown to be associated with cellular genes involved in pathways related to the immune response, virus-host interactions as well as with several viral genes.ConclusionThe study of miRNA expression will contribute to a better understanding of African swine fever virus pathogenesis, essential in the development of any disease control strategy.

ATLS - A parallel loop scheduling scheme for dynamic environments

We here present ATLS, a self scheduling scheme designed for execution of parallel loops in dynamic environments of non-dedicated networked computers. Since grid and volunteer systems based on desktop computers are proving themselves as real and powerful alternatives for parallel computing, new scheduling schemes, better suited for these environments, are required. Our proposal, by tracking several performance change ratios at runtime, is able to properly adjust the load distribution using no prior information of the loop features nor the involved processors. The results obtained during the experiments performed to validate ATLS show that it is possible to improve former contributions of well-known parallel loop scheduling schemes in dynamic environments. The implementation of the scheduler has been done for the R language but, as it is exposed, it can easily be adapted to any other language and parallel loop based application.