Modeling gene expression regulatory networks with the sparse vector autoregressive model

Abstract Background To understand the molecular mechanisms underlying important biological processes, a detailed description of the gene products networks involved is required. In order to define and understand such molecular networks, some statistical methods are proposed in the literature to estimate gene regulatory networks from time-series microarray data. However, several problems still need to be overcome. Firstly, information flow need to be inferred, in addition to the correlation between genes. Secondly, we usually try to identify large networks from a large number of genes (parameters) originating from a smaller number of microarray experiments (samples). Due to this situation, which is rather frequent in Bioinformatics, it is difficult to perform statistical tests using methods that model large gene-gene networks. In addition, most of the models are based on dimension reduction using clustering techniques, therefore, the resulting network is not a gene-gene network but a module-module network. Here, we present the Sparse Vector Autoregressive model as a solution to these problems. Results We have applied the Sparse Vector Autoregressive model to estimate gene regulatory networks based on gene expression profiles obtained from time-series microarray experiments. Through extensive simulations, by applying the SVAR method to artificial regulatory networks, we show that SVAR can infer true positive edges even under conditions in which the number of samples is smaller than the number of genes. Moreover, it is possible to control for false positives, a significant advantage when compared to other methods described in the literature, which are based on ranks or score functions. By applying SVAR to actual HeLa cell cycle gene expression data, we were able to identify well known transcription factor targets. Conclusion The proposed SVAR method is able to model gene regulatory networks in frequent situations in which the number of samples is lower than the number of genes, making it possible to naturally infer partial Granger causalities without any a priori information. In addition, we present a statistical test to control the false discovery rate, which was not previously possible using other gene regulatory network models.

Cloning and characterization of a novel alternatively spliced transcript of the human CHD7 putative helicase

Abstract Background The CHD7 (Chromodomain Helicase DNA binding protein 7) gene encodes a member of the chromodomain family of ATP-dependent chromatin remodeling enzymes. Mutations in the CHD7 gene are found in individuals with CHARGE, a syndrome characterized by multiple birth malformations in several tissues. CHD7 was identified as a binding partner of PBAF complex (Polybromo and BRG Associated Factor containing complex) playing a central role in the transcriptional reprogramming process associated to the formation of multipotent migratory neural crest, a transient cell population associated with the genesis of various tissues. CHD7 is a large gene containing 38 annotated exons and spanning 200 kb of genomic sequence. Although genes containing such number of exons are expected to have several alternative transcripts, there are very few evidences of alternative transcripts associated to CHD7 to date indicating that alternative splicing associated to this gene is poorly characterized. Findings Here, we report the cloning and characterization by experimental and computational studies of a novel alternative transcript of the human CHD7 (named CHD7 CRA_e), which lacks most of its coding exons. We confirmed by overexpression of CHD7 CRA_e alternative transcript that it is translated into a protein isoform lacking most of the domains displayed by the canonical isoform. Expression of the CHD7 CRA_e transcript was detected in normal liver, in addition to the DU145 human prostate carcinoma cell line from which it was originally isolated. Conclusions Our findings indicate that the splicing event associated to the CHD7 CRA_e alternative transcript is functional. The characterization of the CHD7 CRA_e novel isoform presented here not only sets the basis for more detailed functional studies of this isoform, but, also, contributes to the alternative splicing annotation of the CHD7 gene and the design of future functional studies aimed at the elucidation of the molecular functions of its gene products.

Cloning and characterization of bifunctional enzyme farnesyl diphosphate/geranylgeranyl diphosphate synthase from Plasmodium falciparum

Abstract Background Isoprenoids are the most diverse and abundant group of natural products. In Plasmodium falciparum, isoprenoid synthesis proceeds through the methyl erythritol diphosphate pathway and the products are further metabolized by farnesyl diphosphate synthase (FPPS), turning this enzyme into a key branch point of the isoprenoid synthesis. Changes in FPPS activity could alter the flux of isoprenoid compounds downstream of FPPS and, hence, play a central role in the regulation of a number of essential functions in Plasmodium parasites. Methods The isolation and cloning of gene PF3D7_18400 was done by amplification from cDNA from mixed stage parasites of P. falciparum. After sequencing, the fragment was subcloned in pGEX2T for recombinant protein expression. To verify if the PF3D7_1128400 gene encodes a functional rPfFPPS protein, its catalytic activity was assessed using the substrate [4-14C] isopentenyl diphosphate and three different allylic substrates: dimethylallyl diphosphate, geranyl diphosphate or farnesyl diphosphate. The reaction products were identified by thin layer chromatography and reverse phase high-performance liquid chromatography. To confirm the product spectrum formed of rPfFPPS, isoprenic compounds were also identified by mass spectrometry. Apparent kinetic constants KM and Vmax for each substrate were determined by Michaelis–Menten; also, inhibition assays were performed using risedronate. Results The expressed protein of P. falciparum FPPS (rPfFPPS) catalyzes the synthesis of farnesyl diphosphate, as well as geranylgeranyl diphosphate, being therefore a bifunctional FPPS/geranylgeranyl diphosphate synthase (GGPPS) enzyme. The apparent KM values for the substrates dimethylallyl diphosphate, geranyl diphosphate and farnesyl diphosphate were, respectively, 68 ± 5 μM, 7.8 ± 1.3 μM and 2.06 ± 0.4 μM. The protein is expressed constitutively in all intra-erythrocytic stages of P. falciparum, demonstrated by using transgenic parasites with a haemagglutinin-tagged version of FPPS. Also, the present data demonstrate that the recombinant protein is inhibited by risedronate. Conclusions The rPfFPPS is a bifunctional FPPS/GGPPS enzyme and the structure of products FOH and GGOH were confirmed mass spectrometry. Plasmodial FPPS represents a potential target for the rational design of chemotherapeutic agents to treat malaria.

Etablierung von Expressionsystemen für Gene der Indolalkaloid-Biosynthese unter besonderer Berücksichtigung von Cytochrom P450-Enzymen

Etablierung von Expressionsystemen für Gene der Indolalkaloid-Biosynthese unter besonderer Berücksichtigung von Cytochrom P450-Enzymen In der vorliegenden Arbeit wurden Enzyme aus der Arzneipflanze Rauvolfia serpentina bearbeitet. Es wurde versucht, das an der Biosynthese des Alkaloids Ajmalin beteiligte Cytochrom P450-Enzym Vinorin-Hydroxylase heterolog und funktionell zu exprimieren. Ein zunächst unvollständiger, unbekannter Cytochrom P450-Klon konnte komplettiert und eindeutig mittels heterologer Expression in sf9-Insektenzellen als Cinnamoyl-Hydroxylase identifiziert werden. Die Tauglichkeit des Insektenzellsystems für die Untersuchung der Vinorin-Hydroxylase ist auf Grund der deacetylierenden Wirkung der Insektenzellen auf das Substrat Vinorin nicht gegeben. Im Rahmen des Homology Cloning Projektes konnten mehrere Volllängenklone und diverse Teilsequenzen von neuen Cytochrom P450-Klonen ermittelt werden. Ausserdem wurde durch das unspezifische Binden eines degenerierten Primers ein zusätzlicher Klon gefunden, der der Gruppe der löslichen Reduktasen zugeordnet werden konnte. Diese putative Reduktase wurde auf die Aktivität von mehreren Schlüsselenzymen der Ajmalin-Biosynthese durch heterologe Expression in E.coli und anschliessende HPLC-gestützte Aktivitätstests ohne Erfolg geprüft. Bedingt durch die Untauglichkeit des Insektenzellsystems für die Identifizierung der Vinorin-Hydroxylase, wurde ein neuartiges Modul-gestütztes, pflanzliches Expressionsystem etabliert, um vorhandene P450-Volllängenklone auf Vinorin- Hydroxylaseaktivität testen zu können. Die Funktionalität des Systems konnte durch die heterologe Expression der Polyneuridinaldehyd Esterase bestätigt werden. Trotzdem war es bis jetzt nicht möglich, die Cinnamoyl-Hydroxylase als Kontrollenzym für das pflanzliche System oder aber die gesuchte Vinorin- Hydroxylase in aktiver Form zu exprimieren.

An unusual splice defect in the mitofusin 2 gene (MFN2) is associated with degenerative axonopathy in Tyrolean Grey cattle

Tyrolean Grey cattle represent a local breed with a population size of approximately 5000 registered cows. In 2003, a previously unknown neurological disorder was recognized in Tyrolean Grey cattle. The clinical signs of the disorder are similar to those of bovine progressive degenerative myeloencephalopathy (weaver syndrome) in Brown Swiss cattle but occur much earlier in life. The neuropathological investigation of an affected calf showed axonal degeneration in the central nervous system (CNS) and femoral nerve. The pedigrees of the affected calves suggested a monogenic autosomal recessive inheritance. We localized the responsible mutation to a 1.9 Mb interval on chromosome 16 by genome-wide association and haplotype mapping. The MFN2 gene located in this interval encodes mitofusin 2, a mitochondrial membrane protein. A heritable human axonal neuropathy, Charcot-Marie-Tooth disease-2A2 (CMT2A2), is caused by MFN2 mutations. Therefore, we considered MFN2 a positional and functional candidate gene and performed mutation analysis in affected and control Tyrolean Grey cattle. We did not find any non-synonymous variants. However, we identified a perfectly associated silent SNP in the coding region of exon 20 of the MFN2 gene. This SNP is located within a putative exonic splice enhancer (ESE) and the variant allele leads to partial retention of the entire intron 19 and a premature stop codon in the aberrant MFN2 transcript. Thus we have identified a highly unusual splicing defect, where an exonic single base exchange leads to the retention of the preceding intron. This splicing defect represents a potential explanation for the observed degenerative axonopathy. Marker assisted selection can now be used to eliminate degenerative axonopathy from Tyrolean Grey cattle.

Microphthalmia in Texel sheep is associated with a missense mutation in the paired-like homeodomain 3 (PITX3) gene

Microphthalmia in sheep is an autosomal recessive inherited congenital anomaly found within the Texel breed. It is characterized by extremely small or absent eyes and affected lambs are absolutely blind. For the first time, we use a genome-wide ovine SNP array for positional cloning of a Mendelian trait in sheep. Genotyping 23 cases and 23 controls using Illumina's OvineSNP50 BeadChip allowed us to localize the causative mutation for microphthalmia to a 2.4 Mb interval on sheep chromosome 22 by association and homozygosity mapping. The PITX3 gene is located within this interval and encodes a homeodomain-containing transcription factor involved in vertebrate lens formation. An abnormal development of the lens vesicle was shown to be the primary event in ovine microphthalmia. Therefore, we considered PITX3 a positional and functional candidate gene. An ovine BAC clone was sequenced, and after full-length cDNA cloning the PITX3 gene was annotated. Here we show that the ovine microphthalmia phenotype is perfectly associated with a missense mutation (c.338G>C, p.R113P) in the evolutionary conserved homeodomain of PITX3. Selection against this candidate causative mutation can now be used to eliminate microphthalmia from Texel sheep in production systems. Furthermore, the identification of a naturally occurring PITX3 mutation offers the opportunity to use the Texel as a genetically characterized large animal model for human microphthalmia.

Expression profiling of stem cell-related genes in neoadjuvant-treated gastric cancer: a NOTCH2, GSK3B and β-catenin gene signature predicts survival

Cancer stem cell (CSC) based gene expression signatures are associated with prognosis in various tumour types and CSCs are suggested to be particularly drug resistant. The aim of our study was first, to determine the prognostic significance of CSC-related gene expression in residual tumour cells of neoadjuvant-treated gastric cancer (GC) patients. Second, we wished to examine, whether expression alterations between pre- and post-therapeutic tumour samples exist, consistent with an enrichment of drug resistant tumour cells. The expression of 44 genes was analysed in 63 formalin-fixed, paraffin embedded tumour specimens with partial tumour regression (10-50% residual tumour) after neoadjuvant chemotherapy by quantitative real time PCR low-density arrays. A signature of combined GSK3B(high), β-catenin (CTNNB1)(high) and NOTCH2(low) expression was strongly correlated with better patient survival (p<0.001). A prognostic relevance of these genes was also found analysing publically available gene expression data. The expression of 9 genes was compared between pre-therapeutic biopsies and post-therapeutic resected specimens. A significant post-therapeutic increase in NOTCH2, LGR5 and POU5F1 expression was found in tumours with different tumour regression grades. No significant alterations were observed for GSK3B and CTNNB1. Immunohistochemical analysis demonstrated a chemotherapy-associated increase in the intensity of NOTCH2 staining, but not in the percentage of NOTCH2. Taken together, the GSK3B, CTNNB1 and NOTCH2 expression signature is a novel, promising prognostic parameter for GC. The results of the differential expression analysis indicate a prominent role for NOTCH2 and chemotherapy resistance in GC, which seems to be related to an effect of the drugs on NOTCH2 expression rather than to an enrichment of NOTCH2 expressing tumour cells.

Molecular cloning and characterization of equine thymic stromal lymphopoietin

Thymic stromal lymphopoietin (TSLP) is a novel cytokine that plays a central role in T helper 2 (Th2) cell differentiation and allergic inflammation. It is predominantly expressed by epithelial cells, and its expression is increased in patients with atopic dermatitis and asthma. Mice overexpressing TSLP in the skin develop allergic dermatitis and mice overexpressing TSLP in lungs develop asthma-like disease. However, it is not known whether TSLP plays an important role in equine allergies. Therefore, we cloned and sequenced the complete translated region of equine TSLP gene and measured its expression in various tissues. The equine TSLP gene is organized in 4 exons and encodes a protein of 143 amino acids, which has 62% amino acid identity with human TSLP.

The bovine aristaless-like homeobox 4 (ALX4) as a candidate gene for syndactyly

The ALX4 (aristaless-like homeobox 4) gene encodes a paired-type homeodomain transcriptional activator and plays a major role in anterior-posterior pattern formation during limb development. Here, the cloning, genomic structure and expression of the bovine ortholog of the ALX4 gene are reported. The bovine ALX4 gene consists of four exons and is located on BTA15q28-->q29 in a region syntenic to HSA11p11.2. The transcribed ALX4 mRNA encodes a 397-amino-acid protein showing a paired-type homeodomain and a C-terminal stretch of amino acids known as the OAR- or aristaless domain. The predicted protein shares 92.5% identity to human and mouse ALX4 proteins and all three species share almost complete identity in the conserved domains. ALX4 expression was detected by reverse transcriptase polymerase chain reaction in bovine fetal limb bones. The ALX4 gene was evaluated as a candidate gene for bovine syndactyly which has been mapped on the telomeric region of cattle chromosome 15. Sequencing of the four exons with flanking sequences of the bovine ALX4 gene from a panel of 14 affected animals belonging to German Holstein, German Fleckvieh and crossbreds, and 27 unaffected individuals from German Holstein revealed five silent SNPs within the coding region out of eleven SNPs in total. Four SNPs were polymorphic in the affected animals, but in comparison to the genotyped unaffected individuals the genotype distribution showed no evidence for an association to the phenotype. Therefore our data indicate that the ALX4 gene can probably be excluded as candidate gene for bovine syndactyly in the examined animals.

Detection of prion protein gene polymorphisms in Baltic breeds of sheep

A total of 167 sheep belonging to the Estonian whiteheaded mutton, Estonian blackheaded mutton, Lithuanian coarsewool native, Lithuanian blackface and Latvian darkheaded mutton breeds, and a population of sheep kept isolated on the Estonian island of Ruhnu, were sequence-analysed for polymorphisms in the prion protein (PrP) gene, to determine their genotype and the allele frequencies of polymorphisms in PrP known to confer resistance to scrapie. A 939 base pair fragment of exon 3 from the PrP gene was amplified by pcr and analysed by direct sequencing. For animals showing polymorphism at two nucleotide positions, both haplotypes of these double-heterozygous genotypes were further verified by pcr cloning and sequence analysis. Known polymorphisms were observed at codons 136, 154 and 171, and six different haplotypes (arr, ahq, arh, ahr, arq and vrq) were determined. On the basis of these polymorphisms, the six populations of sheep possessed the resistant arr haplotype at different frequencies. The high-risk arq haplotype occurred in high frequencies in all six populations, but vrq, the haplotype carrying the highest risk, occurred at low frequencies and in only three of the populations.

Does a pleiotropic gene explain deafness and blue irises in white cats?

The prevalence of deafness is high in cat populations in which the dominant white gene is segregating. The objective of this study was to investigate whether there is a gene that is responsible for deafness as well as for blue eyes and to establish a plausible mode of inheritance. For this purpose, data from an experimental colony with deaf cats were analyzed. The hearing status was determined by acoustically evoked brain stem responses (BAER). Complex segregation analyses were conducted to find out the most probable mode of inheritance using maximum likelihood procedures. The prevalence of deafness and partial hearing in the experimental colony was 67% and 29%, respectively. The results of the bivariate segregation analysis support the hypothesis of a pleiotropic major gene segregating for deafness and blue iris colour. The high heritability coefficients for both traits, 0.55 and 0.75 respectively, indicate that beside the major gene there is an important influence of polygenic effects.

Classification Using Generalized Partial Least Squares

The advances in computational biology have made simultaneous monitoring of thousands of features possible. The high throughput technologies not only bring about a much richer information context in which to study various aspects of gene functions but they also present challenge of analyzing data with large number of covariates and few samples. As an integral part of machine learning, classification of samples into two or more categories is almost always of interest to scientists. In this paper, we address the question of classification in this setting by extending partial least squares (PLS), a popular dimension reduction tool in chemometrics, in the context of generalized linear regression based on a previous approach, Iteratively ReWeighted Partial Least Squares, i.e. IRWPLS (Marx, 1996). We compare our results with two-stage PLS (Nguyen and Rocke, 2002A; Nguyen and Rocke, 2002B) and other classifiers. We show that by phrasing the problem in a generalized linear model setting and by applying bias correction to the likelihood to avoid (quasi)separation, we often get lower classification error rates.

Congenital hypotrichosis and partial anodontia in a crossbred beef calf

Clinical examination, skin biopsies, skull radiographs, and DNA analysis of a 2-day-old Red Angus-Charolais-Simmental cross bull calf confirmed the diagnosis of congenital hypotrichosis and anodontia defect (HAD), also called anhidrotic ectodermal dysplasia, which is a rare anomaly caused by a deletion in the bovine EDA gene on the X chromosome.

Molecular characterization of the porcine DNAL4 gene

The DNAL4 (dynein, axonemal, light polypeptide 4) gene encodes a light chain of dynein. Dyneins are motor proteins that contribute to axonal transport. Cloning and characterization of the porcine DNAL4 revealed a conserved organization with respect to the human ortholog. The porcine DNAL4 gene consists of 4 exons and codes for a peptide of 105 amino acids. The porcine DNAL4 gene is located on SSC5p15. Analysis of the naturally occurring variation of the DNAL4 gene in pigs from the Piétrain und Duroc breeds revealed five SNPs in non-coding regions of the gene.

Molecular characterisation of BSR4, a novel bradyzoite-specific gene from Neospora caninum

Here we present the identification and cloning of the NcBSR4 gene, the putative Neospora caninum orthologue to the Toxoplasma gondii TgBSR4 gene. To isolate NcBSR4, genome walking PCR was performed on N. caninum genomic DNA using the expressed sequence tag NcEST3c28h02.y1 sequence, which shares a 44% identity with the TgBSR4 gene, as a framework. Nucleotide sequencing of amplified DNA fragments revealed a single uninterrupted 1227 bp open reading frame that encodes a protein of 408 amino acids with 66% similarity to the TgBSR4 antigen. A putative 39-residue signal peptide was found at the NH2-terminus, followed by a hydrophilic region. At the COOH-terminus, a potential site for a glycosylphosphatidylinositol anchor was identified at amino acid 379. A polyclonal serum against recombinant NcBSR4 protein was raised in rabbits, and immunolabelling demonstrated stage-specific expression of the NcBSR4 antigen in N. caninum bradyzoites produced in vitro and in vivo. Furthermore, RT-PCR analysis showed a slight increase of NcBSR4 transcripts in bradyzoites generated during in vitro tachyzoite-to-bradyzoite stage-conversion, suggesting that this gene is specifically expressed at the bradyzoite stage and that its transcription relies on the switch to this stage.