Search

Macroalgae provide important habitats for microorganisms, with surfaces that are rich in nutrients and a microbial niche for isolation of different bacterial groups. In this study, we evaluated the production of enzymes (amylases, lipases, cellulases and agarases) and siderophores by strains isolated from the macroalga Ulva lactuca growing in the Santa Marta region of the Colombian Caribbean Sea. We taxonomically identified a subset of the bacteria and tested 207 bacterial isolates. Of these, 25 (12%) and 121 (58%) produced agarases and siderophores, respectively. Out of the isolated bacteria with amylolytic (97), cellulolytic (118) and lipolytic (26) activity, 31 (32%), 78 (66%) and 4 (15%) strains, respectively, produced a zone of clearance >1.5 cm, indicating substantial activity. Enzymatic activities and siderophore production were statistically different between years of sampling, and principal component analysis showed grouping for samples from the same year. These activities and production were recorded in bacterial strains belonging to the genera Vibrio, Pseudomonas and Bacillus. The results show that marine bacterial cultures isolated from the macroalga U. lactuca are a source of enzymes and siderophores that may have potential for biotechnological and industrial processes.

In this article, we introduce a mixture inverse Gaussian (MIG) model based on the Student-t distribution and apply it to bacterium-based protein production for food industry. This model is mainly useful to describe data that follow positively skewed distributions and accommodate atypical observations in a better way than its classical version. Specifically, we present a characterization of the MIG-t distribution. In addition, we carry out a hazard analysis of this distribution centered mainly on its hazard rate. Furthermore, we discuss the maximum likelihood method, which produces–in this case–robust parameter estimates. Moreover, to evaluate the potential influence of atypical observations, we produce a diagnostic analysis for the model. Finally, we apply the obtained results to novel bacterium-based protein production data and statistically compare two types of protein producers using the likelihood ratio test based on the MIG-t model as an alternative methodology to the procedures available until now. This fact is very important, since the evaluation of protein production using both constructions allows practitioners to choose the most productive one before the bacterial culture is scaled to an industrial level.En este artículo, introducimos un modelo Gaussiano inverso (MIG) mezclado basado en la distribución t-Student y lo aplicamos a la producción de proteínas basada en bacterias para la industria de alimentos. Este modelo es especialmente útil para describir datos que siguen una distribución con sesgo positivo ya que permite acomodar observaciones atípicas de mejor forma que su versión clásica. Específicamente, presentamos una caracterización de la distribución MIG-t y realizamos un análisis de confiabilidad de esta distribución centrado principalmente en la tasa de fallas. También, discutimos el método de verosimilitud máxima, el cual proporciona en este caso estimaciones robustas de los parámetros del modelo. Con el fin de evaluar la influencia potencial de observaciones atípicas, proponemos un análisis de diagnóstico para la distribución. Finalmente, aplicamos los resultados obtenidos al análisis de datos nuevos de producción de proteína basada en bacterias utilizada en la industria de alimentos y comparamos estadísticamente dos tipos de bacterias productoras usando la prueba de razón de verosimilitudes basada en el modelo MIG-t como una metodología alternativa a los procedimientos disponibles a la fecha. Este punto es muy importante, ya que la evaluación de producción de proteínas usando dos construcciones distintas permite a los investigadores escoger el tipo más productivo antes de proceder al cultivo industrial a gran escala.

The transcriptomic analysis of microarray and RNA-Seq datasets followed our own bioinformatic pipeline to identify a transcriptional regulatory network of lung cancer. Twenty-six transcription factors are dysregulated and co-expressed in most of the lung cancer and pulmonary arterial hypertension datasets, which makes them the most frequently dysregulated transcription factors. Co-expression, gene regulatory, coregulatory, and transcriptional regulatory networks, along with fibration symmetries, were constructed to identify common connection patterns, alignments, main regulators, and target genes in order to analyze transcription factor complex formation, as well as its synchronized co-expression patterns in every type of lung cancer. The regulatory function of the most frequently dysregulated transcription factors over lung cancer deregulated genes was validated with ChEA3 enrichment analysis. A Kaplan-Meier plotter analysis linked the dysregulation of the top transcription factors with lung cancer patients' survival. Our results indicate that lung cancer has unique and common deregulated genes and transcription factors with pulmonary arterial hypertension, co-expressed and regulated in a coordinated and cooperative manner by the transcriptional regulatory network that might be associated with critical biological processes and signaling pathways related to the acquisition of the hallmarks of cancer, making them potentially relevant tumor biomarkers for lung cancer early diagnosis and targets for the development of personalized therapies against lung cancer.

Multiple test corrections are a fundamental step in the analysis of differentially expressed genes, as the number of tests performed would otherwise inflate the false discovery rate (FDR). Recent methods for P-value correction involve a regression model in order to include covariates that are informative of the power of the test. Here, we present Progressive proportions plot (Prog-Plot), a visual tool to identify the functional relationship between the covariate and the proportion of P-values consistent with the null hypothesis. The relationship between the proportion of P-values and the covariate to be included is needed, but there are no available tools to verify it. The approach presented here aims at having an objective way to specify regression models instead of relying on prior knowledge., Competing Interests: Competing interests The authors declare no competing or financial interests., (© 2023. Published by The Company of Biologists Ltd.)

The use of a new bioinformatics pipeline allowed the identification of deregulated transcription factors (TFs) coexpressed in lung cancer that could become biomarkers of tumor establishment and progression. A gene regulatory network (GRN) of lung cancer was created with the normalized gene expression levels of differentially expressed genes (DEGs) from the microarray dataset GSE19804. Moreover, coregulatory and transcriptional regulatory network (TRN) analyses were performed for the main regulators identified in the GRN analysis. The gene targets and binding motifs of all potentially implicated regulators were identified in the TRN and with multiple alignments of the TFs' target gene sequences. Six transcription factors (E2F3, FHL2, ETS1, KAT6B, TWIST1, and RUNX2) were identified in the GRN as essential regulators of gene expression in non-small-cell lung cancer (NSCLC) and related to the lung tumoral process. Our findings indicate that RUNX2 could be an important regulator of the lung cancer GRN through the formation of coregulatory complexes with other TFs related to the establishment and progression of lung cancer. Therefore, RUNX2 could become an essential biomarker for developing diagnostic tools and specific treatments against tumoral diseases in the lung after the experimental validation of its regulatory function.

The bioinformatic pipeline previously developed in our research laboratory is used to identify potential general and specific deregulated tumor genes and transcription factors related to the establishment and progression of tumoral diseases, now comparing lung cancer with other two types of cancer. Twenty microarray datasets were selected and analyzed separately to identify hub differentiated expressed genes and compared to identify all the deregulated genes and transcription factors in common between the three types of cancer and those unique to lung cancer. The winning DEGs analysis allowed to identify an important number of TFs deregulated in the majority of microarray datasets, which can become key biomarkers of general tumors and specific to lung cancer. A coexpression network was constructed for every dataset with all deregulated genes associated with lung cancer, according to DAVID's tool enrichment analysis, and transcription factors capable of regulating them, according to oPOSSUM´s tool. Several genes and transcription factors are coexpressed in the networks, suggesting that they could be related to the establishment or progression of the tumoral pathology in any tissue and specifically in the lung. The comparison of the coexpression networks of lung cancer and other types of cancer allowed the identification of common connectivity patterns with deregulated genes and transcription factors correlated to important tumoral processes and signaling pathways that have not been studied yet to experimentally validate their role in lung cancer. The Kaplan-Meier estimator determined the association of thirteen deregulated top winning transcription factors with the survival of lung cancer patients. The coregulatory analysis identified two top winning transcription factors networks related to the regulatory control of gene expression in lung and breast cancer. Our transcriptomic analysis suggests that cancer has an important coregulatory network of transcription factors related to the acquisition of the hallmarks of cancer. Moreover, lung cancer has a group of genes and transcription factors unique to pulmonary tissue that are coexpressed during tumorigenesis and must be studied experimentally to fully understand their role in the pathogenesis within its very complex transcriptomic scenario. Therefore, the downstream bioinformatic analysis developed was able to identify a coregulatory metafirm of cancer in general and specific to lung cancer taking into account the great heterogeneity of the tumoral process at cellular and population levels.

Exosomes carry molecules of great biological and clinical interest, such as miRNAs. The contents of exosomes vary between healthy controls and cancer patients. Therefore, miRNAs and other molecules transported in exosomes are considered a potential source of diagnostic and prognostic biomarkers in cancer. Many miRNAs have been detected in recent years. Consequently, a substantial amount of miRNA-related data comparing patients and healthy individuals is available, which contributes to a better understanding of the initiation, development, malignancy, and metastasis of cancer using non-invasive sampling procedures. However, a re-analysis of available ncRNA data is rare. This study used available data about miRNAs in exosomes comparing healthy individuals and cancer patients to identify possible global changes related to the presence of cancer. A robust transcriptomic analysis identified two common miRNAs (miR-495-3p and miR-543) deregulated in five cancer datasets. They had already been implicated in different cancers but not reported in exosomes circulating in blood. The study also examined their target genes and the implications of these genes for functional processes.

Noncoding RNAs (ncRNAs) play prominent roles in the regulation of gene expression via their interactions with other biological molecules such as proteins and nucleic acids. Although much of our knowledge about how these ncRNAs operate in different biological processes has been obtained from experimental findings, computational biology can also clearly substantially boost this knowledge by suggesting possible novel interactions of these ncRNAs with other molecules. Computational predictions are thus used as an alternative source of new insights through a process of mutual enrichment because the information obtained through experiments continuously feeds through into computational methods. The results of these predictions in turn shed light on possible interactions that are subsequently validated experimentally. This review describes the latest advances in databases, bioinformatic tools, and new in silico strategies that allow the establishment or prediction of biological interactions of ncRNAs, particularly miRNAs and lncRNAs. The ncRNA species described in this work have a special emphasis on those found in humans, but information on ncRNA of other species is also included.

Aim: This study aims to investigate similarities and differences using lncRNA and mRNA coexpression network analysis in African ancestry (AA) and European ancestry (EA) among prostate cancer (PCa) patients., Methods: We performed weighted gene coexpression network analysis of the expression from 49 of AA and 49 of EA to identify lncRNAs-mRNAs., Results: 27 lncRNAs and 36 mRNAs were highly expressed in patients of AA. Two mRNAs and their antisense lncRNAs were expressed. Additionally, seven mRNAs were DE or coexpressed and had an impact on survival., Conclusion: We present a list of lncRNAs and mRNAs that were DE and coexpressed when comparing patients of AA and EA, and these data are a resource for future studies to understand the role of lncRNAs., Competing Interests: Financial & competing interests disclosure This study was supported by the School of Medicine and Health Sciences and the Faculty of Natural Sciences at the Universidad del Rosario. The authors have no other relevant affiliations or financial involvement with any organization or entity with a financial interest in or financial conflict with the subject matter or materials discussed in the manuscript apart from those disclosed. No writing assistance was utilized in the production of this manuscript., (© 2021 César Payán-Gómez.)

Microbial communities live on macroalgal surfaces. The identity and abundance of the bacteria making these epiphytic communities depend on the macroalgal host and the environmental conditions. Macroalgae rely on epiphytic bacteria for basic functions (spore settlement, morphogenesis, growth, and protection against pathogens). However, these marine bacterial-macroalgal associations are still poorly understood for macroalgae inhabiting the Colombian Caribbean. This study aimed at characterizing the epiphytic bacterial community from macroalgae of the species Ulva lactuca growing in La Punta de la Loma (Santa Marta, Colombia). We conducted a 16S rRNA gene sequencing-based study of these microbial communities sampled twice a year between 2014 and 2016. Within these communities, the Proteobacteria, Bacterioidetes, Cyanobacteria, Deinococcus-Thermus and Actinobacteria were the most abundant phyla. At low taxonomic levels, we found high variability among epiphytic bacteria from U. lactuca and bacterial communities associated with macroalgae from Germany and Australia. We observed differences in the bacterial community composition across years driven by abundance shifts of Rhodobacteraceae Hyphomonadaceae, and Flavobacteriaceae, probably caused by an increase of seawater temperature. Our results support the need for functional studies of the microbiota associated with U. lactuca, a common macroalga in the Colombian Caribbean Sea.

Background: Analysis of patients with chromosomal abnormalities, including Turner syndrome and Klinefelter syndrome, has highlighted the importance of X-linked gene dosage as a contributing factor for disease susceptibility. Escape from X-inactivation and X-linked imprinting can result in transcriptional differences between normal men and women as well as in patients with sex chromosome abnormalities., Objective: To identify differentially expressed genes among patients with Turner (45,X) and Klinefelter (46,XXY) syndrome using bioinformatics analysis., Methodology: Two gene expression data sets of Turner (45,X) and Klinefelter syndrome (47,XXY) were obtained from the Gene Omnibus Expression (GEO) database of the National Center for Biotechnology Information (NCBI). Statistical analysis was performed using R Bioconductor libraries. Differentially expressed genes (DEGs) were determined using significance analysis of microarray (SAM). The functional annotation of the DEGs was performed with DAVID v6.8 (The Database for Annotation, Visualization, and Integrated Discovery)., Results: There are no genes over-expressed simultaneously in both diseases. However, when crossing the list of under-expressed genes for 45,X cells and the list of over-expressed genes for 47,XXY cells, there are 16 common genes: SLC25A6, AKAP17A, ASMTL, KDM5C, KDM6A, ATRX, CSF2RA, DHRSX, CD99, ZBED1, EIF1AX, MVB12B, SMC1A, P2RY8, DOCK7, DDX3X, eight of which are involved in the regulation of gene expression by epigenetic mechanisms, regulation of splicing processes and protein synthesis., Conclusion: Of the 16 identified as under-expressed in 45,X cells and over-expressed in 47,XXY cells, 14 are located in X chromosome and 2 in autosomal chromosome; 8 of these genes are involved in the regulation of gene expression: 5 genes are related to epigenetic mechanisms, 2 in regulation of splicing processes, and 1 in the protein synthesis process. Our results are limited by it being the product of a bioinformatic analysis from mRNA isolated from whole blood, this makes necessary further exploration of the relationships between these genes and Turner syndrome and Klinefelter syndrome in the future., (© 2020 The Authors. Molecular Genetics & Genomic Medicine published by Wiley Periodicals LLC.)

The incidence of patients under 55 years old diagnosed with Prostate Cancer (EO-PCa) has increased during recent years. The molecular biology of PCa cancer in this group of patients remains unclear. Here, we applied weighted gene coexpression network analysis of the expression of miRNAs from 24 EO-PCa patients (38-45 years) and 25 late-onset PCa patients (LO-PCa, 71-74 years) to identify key miRNAs in EO-PCa patients. In total, 69 differentially expressed miRNAs were identified. Specifically, 26 and 14 miRNAs were exclusively deregulated in young and elderly patients, respectively, and 29 miRNAs were shared. We identified 20 hub miRNAs for the network built for EO-PCa. Six of these hub miRNAs exhibited prognostic significance in relapse-free or overall survival. Additionally, two of the hub miRNAs were coexpressed with mRNAs of genes previously identified as deregulated in EO-PCa and in the most aggressive forms of PCa in African-American patients compared with Caucasian patients. These genes are involved in activation of immune response pathways, increased rates of metastasis and poor prognosis in PCa patients. In conclusion, our analysis identified miRNAs that are potentially important in the molecular pathology of EO-PCa. These genes may serve as biomarkers in EO-PCa and as possible therapeutic targets.

Background: Epidemiological and clinical evidence points cancer comorbidity with pulmonary chronic disease. The acquisition of some hallmarks of cancer by cells affected with lung pathologies as a cell adaptive mechanism to a shear stress, suggests that could be associated with the establishment of tumoral processes. Objective: To propose a bioinformatic pipeline for the identification of all deregulated genes and the transcriptional regulators (TFs) that are coexpressed during lung cancer establishment, and therefore could be important for the acquisition of the hallmarks of cancer. Methods: Ten microarray datasets (six of lung cancer, four of lung diseases) comparing normal and diseases-related lung tissue were selected to identify hub differentiated expressed genes (DEGs) in common between lung pathologies and lung cancer, along with transcriptional regulators through the utilization of specialized libraries from R language. DAVID bioinformatics tool for gene enrichment analyses was used to identify genes with experimental evidence associated to tumoral processes and signaling pathways. Coexpression networks of DEGs and TFs in lung cancer establishment were created with Coexnet library, and a survival analysis of the main hub genes was made. Results: Two hundred ten DEGs were identified in common between lung cancer and other lung diseases related to the acquisition of tumoral characteristics, which are coexpressed in a lung cancer network with TFs, suggesting that could be related to the establishment of the tumoral pathology in lung. The comparison of the coexpression networks of lung cancer and other lung diseases allowed the identification of common connectivity patterns (CCPs) with DEGs and TFs correlated to important tumoral processes and signaling pathways, that haven´t been studied to experimentally validate their role in the early stages of lung cancer. Some of the TFs identified showed a correlation between its expression levels and the survival of lung cancer patients. Conclusion: Our findings indicate that lung diseases share genes with lung cancer which are coexpressed in lung cancer, and might be able to explain the epidemiological observations that point to direct and inverse comorbid associations between some chronic lung diseases and lung cancer and represent a complex transcriptomic scenario., (Copyright © 2019 Otálora-Otálora, Florez, López-Kleine, Canas Arboleda, Grajales Urrego and Rojas.)

Microarray technology is widely recognized as one of the most important tools when it comes to understanding genetic expression in biological processes. In light of the thousands of gene expression level measurements (including measurements across a number of conditions), identifying differentially expressed genes necessarily implies data mining or large-scale multiple testing procedures. To date, advances with regard to this field have been multivariate-descriptive or inferential-univariate in nature and therefore have important limitations regarding the biological validity of detected genes. In the present article, we present a new multivariate inferential method designed to detect active differentially expressed genes in gene expression data. The proposed method estimates false discovery rates using artificial components. Our method excels when applied to the most common gene expression data structures, providing new insights into differentially expressed genes. The method described herein was programmed in an R-Bioconductor package called acde that has been available since 2015.

Runt‑related transcription factor (RUNX) proteins belong to a transcription factors family known as master regulators of important embryonic developmental programs. In the last decade, the whole family has been implicated in the regulation of different oncogenic processes and signaling pathways associated with cancer. Furthermore, a suppressor tumor function has been also reported, suggesting the RUNX family serves key role in all different types of cancer. In this review, the known biological characteristics, specific regulatory abilities and experimental evidence of RUNX proteins will be analyzed to demonstrate their oncogenic potential and tumor suppressor abilities during oncogenic processes, suggesting their importance as biomarkers of cancer. Additionally, the importance of continuing with the molecular studies of RUNX proteins' and its dual functions in cancer will be underlined in order to apply it in the future development of specific diagnostic methods and therapies against different types of cancer.

Background: Allergic conditions have shown an increase in the past few decades. Hygiene, changes in lifestyle, diet and vitamin D have been blamed for this increase., Objective: To determine the role of environment, diet and vitamin D in atopic diseases in a population of Colombian Afro-descendant children., Methods: Cross-sectional, descriptive, observational study in 200 Afro-descendant children of rural and urban areas from northern Colombia. Lifestyle and diet were established by questionnaires, atopy was assessed with skin tests and total and specific IgE and vitamin D levels were determined by ELISA., Results: Results: Atopy was more prevalent in the urban population (24 % versus 7 %, p < 0.001). Blomia tropicalis was the most common allergen (80.6 %). The PCA indicated two dietary patterns. In the rural area children, the consumption of dairy products and fruit/vegetables was higher, and related with a protective effect on atopy (OR, 0.21 [p < 0.027] and OR, 0.11 [ p < 0.04]). A similar pattern was observed with total IgE and vitamin D levels (2450.73 versus 777.56 kU/mL [p < 0.01] and 59.32 ng/mL versus 31.14 ng/mL [p < 0.001], respectively)., Conclusion: Allergic conditions were less prevalent in rural area residents, possibly owing to higher consumption of unpasteurized dairy products and fruit/vegetables, as well as higher levels of vitamin D.

Co-expression networks may provide insights into the patterns of molecular interactions that underlie cellular processes. To obtain a better understanding of miRNA expression patterns in gastric adenocarcinoma and to provide markers that can be associated with histopathological findings, we performed weighted gene correlation network analysis (WGCNA) and compare it with a supervised analysis. Integrative analysis of target predictions and miRNA expression profiles in gastric cancer samples was also performed. WGCNA identified a module of co-expressed miRNAs that were associated with histological traits and tumor condition. Hub genes were identified based on statistical analysis and network centrality. The miRNAs 100, let-7c, 125b and 99a stood out for their association with the diffuse histological subtype. The 181 miRNA family and miRNA 21 highlighted for their association with the tumoral phenotype. The integrated analysis of miRNA and gene expression profiles showed the let-7 miRNA family playing a central role in the regulatory relationships., (Copyright © 2016 Elsevier Inc. All rights reserved.)

Gene expression data (microarrays and RNA-sequencing data) as well as other kinds of genomic data can be extracted from publicly available genomic data. Here, we explain how to apply multivariate cluster and classification methods on gene expression data. These methods have become very popular and are implemented in freely available software in order to predict the participation of gene products in a specific functional category of interest. Taking into account the availability of data and of these methods, every biological study should apply them in order to obtain knowledge on the organism studied and functional category of interest. A special emphasis is made on the nonlinear kernel classification methods.

Background: The clinical course of chronic lymphocytic leukemia (CLL) is highly variable; some patients follow an indolent course, but others progress to a more advanced stage. The mutational status of rearranged immunoglobulin heavy chain variable (IGVH) genes in CLL is a feature that is widely recognized for dividing patients into groups that are related to their prognoses. However, the regulatory programs associated with the IGVH statuses are poorly understood, and markers that can precisely predict survival outcomes have yet to be identified., Methods: In this study, (i) we reconstructed gene regulatory networks in CLL by applying an information-theoretic approach to the expression profiles of 5 cohorts. (ii) We applied master regulator analysis (MRA) to these networks to identify transcription factors (TFs) that regulate an IGVH mutational status signature. The IGVH mutational status signature was developed by searching for differentially expressed genes between the IGVH mutational statuses in numerous CLL cohorts. (iii) To evaluate the biological implication of the inferred regulators, prognostic values were determined using time to treatment (TTT) and overall survival (OS) in two different cohorts., Results: A robust IGVH expression signature was obtained, and various TFs emerged as regulators of the signature in most of the reconstructed networks. The TF targets expression profiles exhibited significant differences with respect to survival, which allowed the definition of a reduced profile with a high value for OS. TCF7 and its targets stood out for their roles in progression., Conclusion: TFs and their targets, which were obtained merely from inferred regulatory associations, have prognostic implications and reflect a regulatory context for prognosis.

Genomes from four Clostridium sp. strains considered to be mesophilic anaerobic bacteria, isolated from crop soil in Colombia, with a strong potential to produce alcohols like 1,3-propanediol, were analyzed. We present the draft genome of these strains, which will be useful for developing genetic engineering strategies., (Copyright © 2015 Rosas-Morales et al.)

Gene co-expression networks (GCNs) are graphic representations that depict the coordinated transcription of genes in response to certain stimuli. GCNs provide functional annotations of genes whose function is unknown and are further used in studies of translational functional genomics among species. In this work, a methodology for the reconstruction and comparison of GCNs is presented. This approach was applied using gene expression data that were obtained from immunity experiments in Arabidopsis thaliana, rice, soybean, tomato and cassava. After the evaluation of diverse similarity metrics for the GCN reconstruction, we recommended the mutual information coefficient measurement and a clustering coefficient-based method for similarity threshold selection. To compare GCNs, we proposed a multivariate approach based on the Principal Component Analysis (PCA). Branches of plant immunity that were exemplified by each experiment were analyzed in conjunction with the PCA results, suggesting both the robustness and the dynamic nature of the cellular responses. The dynamic of molecular plant responses produced networks with different characteristics that are differentiable using our methodology. The comparison of GCNs from plant pathosystems, showed that in response to similar pathogens plants could activate conserved signaling pathways. The results confirmed that the closeness of GCNs projected on the principal component space is an indicative of similarity among GCNs. This also can be used to understand global patterns of events triggered during plant immune responses.

Relationships between genes are best represented using networks constructed from information of different types, with metabolic information being the most valuable and widely used for genetic network reconstruction. Other types of information are usually also available, and it would be desirable to systematically include them in algorithms for network reconstruction. Here, we present an algorithm to construct a global metabolic network that uses all available enzymatic and metabolic information about the organism. We construct a global enzymatic network (GEN) with a total of 4226 nodes (EC numbers) and 42723 edges representing all known metabolic reactions. As an example we use microarray data for Arabidopsis thaliana and combine it with the metabolic network constructing a final gene interaction network for this organism with 8212 nodes (genes) and 4606,901 edges. All scripts are available to be used for any organism for which genomic data is available.

Microarray experiments are often conducted in order to compare gene expression between two conditions. Tests to detected mean differential expression of genes between conditions are conducted applying correction for multiple testing. Seldom, relationships between gene expression and microarray conditions are investigated in a multivariate approach. Here we propose determining the relationship between genes and conditions using a Principal Component Analysis (PCA) space and classifying genes to one of two biological conditions based on their position relative to a direction on the PC space representing each condition.

Recent advances in genomic and post-genomic technologies have provided the opportunity to generate a previously unimaginable amount of information. However, biological knowledge is still needed to improve the understanding of complex mechanisms such as plant immune responses. Better knowledge of this process could improve crop production and management. Here, we used holistic analysis to combine our own microarray and RNA-seq data with public genomic data from Arabidopsis and cassava in order to acquire biological knowledge about the relationships between proteins encoded by immunity-related genes (IRGs) and other genes. This approach was based on a kernel method adapted for the construction of gene networks. The obtained results allowed us to propose a list of new IRGs. A putative function in the immunity pathway was predicted for the new IRGs. The analysis of networks revealed that our predicted IRGs are either well documented or recognized in previous co-expression studies. In addition to robust relationships between IRGs, there is evidence suggesting that other cellular processes may be also strongly related to immunity., (Copyright © 2013. Production and hosting by Elsevier Ltd.)

Techniques in molecular biology have permitted the gathering of an extremely large amount of information relating organisms and their genes. The current challenge is assigning a putative function to thousands of genes that have been detected in different organisms. One of the most informative types of genomic data to achieve a better knowledge of protein function is gene expression data. Based on gene expression data and assuming that genes involved in the same function should have a similar or correlated expression pattern, a function can be attributed to those genes with unknown functions when they appear to be linked in a gene co-expression network (GCN). Several tools for the construction of GCNs have been proposed and applied to plant gene expression data. Here, we review recent methodologies used for plant gene expression data and compare the results, advantages and disadvantages in order to help researchers in their choice of a method for the construction of GCNs.

The main objective of the present study was to reanalyse tomato expression data that was previously submitted to the Tomato Expression Database to dissect the resistance/defence genomic and metabolic responses of tomato to Phytophthora infestans under field conditions. Overrepresented gene sets belonging to chromosome 10 were identified using the Gene Set Enrichment Analysis, and we found that these genes tend to be located towards the end of the chromosome 10. An analysis of syntenic regions between Arabidopsis thaliana chromosomes and the tomato chromosome 10 allowed us to identify conserved regions in the two genomes. In addition to allowing for the identification of tomato candidate genes participating in resistance/defence in the field, this approach allowed us to investigate the relationships of the candidate genes with chromosomal position and participation in metabolic functions, thus offering more insight into the phenomena occurring during the infection process., (Copyright © 2013 Elsevier Inc. All rights reserved.)

Publicly available genomic data are a great source of biological knowledge that can be extracted when appropriate data analysis is used. Predicting the biological function of genes is of interest to understand molecular mechanisms of virulence and resistance in pathogens and hosts and is important for drug discovery and disease control. This is commonly done by searching for similar gene expression behavior. Here, we used publicly available Streptococcus pyogenes microarray data obtained during primate infection to identify genes that have a potential influence on virulence and Phytophtora infestance inoculated tomato microarray data to identify genes potentially implicated in resistance processes. This approach goes beyond co-expression analysis. We employed a quasi-likelihood model separated by primate gender/inoculation condition to model median gene expression of known virulence/resistance factors. Based on this model, an influence analysis considering time course measurement was performed to detect genes with atypical expression. This procedure allowed for the detection of genes potentially implicated in the infection process. Finally, we discuss the biological meaning of these results, showing that influence analysis is an efficient and useful alternative for functional gene prediction.

Since the introduction of genome sequencing techniques several methods for genomic data preprocessing and analysis have been published and applied to answer different biological questions. Rarely, multivariate methods have been used to extract knowledge about protein roles. Two of the most informative types of data are gene expression data (microarrays) and phylogenetic profiles indicating presence of genes in other organisms and therefore providing information about their co-evolution. Here we show that these two types of data, analyzed by means of principal component analysis and non parametric discriminant analysis, provide useful information about protein function and their participation in virulence processes.

Interesting biological information as, for example, gene expression data (microarrays), can be extracted from publicly available genomic data. As a starting point in order to narrow down the great possibilities of wet lab experiments, global high throughput data and available knowledge should be used to infer biological knowledge and emit biological hypothesis. Here, based on microarray data, we propose the use of cluster and classification methods that have become very popular and are implemented in freely available software in order to predict the participation in virulence mechanisms of different proteins coded by genes of the pathogen Streptococcus pyogenes. Confidence of predictions is based on classification errors of known genes and repetitive prediction by more than three methods. A special emphasis is done on the nonlinear kernel classification methods used. We propose a list of interesting candidates that could be virulence factors or that participate in the virulence process of S. pyogenes. Biological validations should start using this list of candidates as they show similar behavior to known virulence factors.