Predicting functional gene links from phylogenetic-statistical analyses of whole genomes

An important element of the developing field of proteomics is to understand protein-protein interactions and other functional links amongst genes. Across-species correlation methods for detecting functional links work on the premise that functionally linked proteins will tend to show a common pattern of presence and absence across a range of genomes. We describe a maximum likelihood statistical model for predicting functional gene linkages. The method detects independent instances of the correlated gain or loss of pairs of proteins on phylogenetic trees, reducing the high rates of false positives observed in conventional across-species methods that do not explicitly incorporate a phylogeny. We show, in a dataset of 10,551 protein pairs, that the phylogenetic method improves by up to 35% on across-species analyses at identifying known functionally linked proteins. The method shows that protein pairs with at least two to three correlated events of gain or loss are almost certainly functionally linked. Contingent evolution, in which one gene's presence or absence depends upon the presence of another, can also be detected phylogenetically, and may identify genes whose functional significance depends upon its interaction with other genes. Incorporating phylogenetic information improves the prediction of functional linkages. The improvement derives from having a lower rate of false positives and from detecting trends that across-species analyses miss. Phylogenetic methods can easily be incorporated into the screening of large-scale bioinformatics datasets to identify sets of protein links and to characterise gene networks.

A novel alpha-galactosidase from I'ifidobacterium bifidum with transgalactosylating properties: gene molecular cloning and heterologous expression

A genomic library of Bifidobacterium bifidum (NCIMB 41171) DNA was constructed in Escherichia coli RA11r (melA(-)B(+)) and one alpha-galactosidase encoding gene was isolated. Conceptual translation combined with insertional mutagenesis analysis indicated an open reading frame (ORF) of 759 amino acid (aa) residues encoding an alpha-galactosidase (named as MelA) of 82.8 kDa. Partial purification and characterisation showed that the enzyme had an apparent native molecular mass of a parts per thousand 243 kDa and a subunit size of a parts per thousand 85 kDa. The enzyme belongs to glycosyl hydrolases 36 family with high aa sequence similarities (a parts per thousand 73%) to other known alpha-galactosidases of bifidobacterial origin. Under optimum pH conditions for activity (pH 6.0) and high melibiose concentration (40% w/v), the enzyme was able to form oligosaccharides with degree of polymerisation (DP) a parts per thousand yen3 at higher concentration than DP = 2, with a total yield of 20.5% (w/w).

Variation in the FFAR1 gene modifies BMI, body composition and plasma lipids but not plasma insulin or Beta-cell function in overweight subjects

Background FFAR1 receptor is a long chain fatty acid G-protein coupled receptor which is expressed widely, but found in high density in the pancreas and central nervous system. It has been suggested that FFAR1 may play a role in insulin sensitivity, lipotoxicity and is associated with type 2 diabetes. Here we investigate the effect of three common SNPs of FFAR1 (rs2301151; rs16970264; rs1573611) on pancreatic function, BMI, body composition and plasma lipids. Methodology/Principal Findings For this enquiry we used the baseline RISCK data, which provides a cohort of overweight subjects at increased cardiometabolic risk with detailed phenotyping. The key findings were SNPs of the FFAR1 gene region were associated with differences in body composition and lipids, and the effects of the 3 SNPs combined were cumulative on BMI, body composition and total cholesterol. The effects on BMI and body fat were predominantly mediated by rs1573611 (1.06 kg/m2 higher (P = 0.009) BMI and 1.53% higher (P = 0.002) body fat per C allele). Differences in plasma lipids were also associated with the BMI-increasing allele of rs2301151 including higher total cholesterol (0.2 mmol/L per G allele, P = 0.01) and with the variant A allele of rs16970264 associated with lower total (0.3 mmol/L, P = 0.02) and LDL (0.2 mmol/L, P<0.05) cholesterol, but also with lower HDL-cholesterol (0.09 mmol/L, P<0.05) although the difference was not apparent when controlling for multiple testing. There were no statistically significant effects of the three SNPs on insulin sensitivity or beta cell function. However accumulated risk allele showed a lower beta cell function on increasing plasma fatty acids with a carbon chain greater than six. Conclusions/Significance Differences in body composition and lipids associated with common SNPs in the FFAR1 gene were apparently not mediated by changes in insulin sensitivity or beta-cell function.

Pathogenic LRRK2 mutations do not alter gene expression in cell model systems or human brain tissue

Point mutations in LRRK2 cause autosomal dominant Parkinson's disease. Despite extensive efforts to determine the mechanism of cell death in patients with LRRK2 mutations, the aetiology of LRRK2 PD is not well understood. To examine possible alterations in gene expression linked to the presence of LRRK2 mutations, we carried out a case versus control analysis of global gene expression in three systems: fibroblasts isolated from LRRK2 mutation carriers and healthy, non-mutation carrying controls; brain tissue from G2019S mutation carriers and controls; and HEK293 inducible LRRK2 wild type and mutant cell lines. No significant alteration in gene expression was found in these systems following correction for multiple testing. These data suggest that any alterations in basal gene expression in fibroblasts or cell lines containing mutations in LRRK2 are likely to be quantitatively small. This work suggests that LRRK2 is unlikely to play a direct role in modulation of gene expression, although it remains possible that this protein can influence mRNA expression under pathogenic cicumstances.

Distribution, gene sequence and expression in vivo of the plasmid encoded fimbrial antigen of Salmonella serotype Enteritidis

The pefA gene which encoded the serotype associated plasmid (SAP) mediated fimbrial major subunit antigen of Salmonella enterica serotype Typhimurium shared genetic identity with 128 of 706 salmonella isolates as demonstrated by dot (colony) hybridization. Seventy-seven of 113 isolates of Typhimurium and individual isolates of serotypes Bovis-morbificans, Cholerae-suis and Enteritidis phage type 9b hybridized pefA strongly, whereas 48 isolates of Enteritidis hybridized pefA weakly and one Enteritidis isolate of phage type 14b failed to hybridize. Individual isolates of 294 serotypes and 247 individual isolates of serotype Dublin did not hybridize pefA. Southern hybridization of plasmids extracted from Enteritidis demonstrated that the pefA gene probe hybridized strongly an atypical SAP of 80 kb in size harboured by one Enteritidis isolate of phage-type 9b, whereas the typical SAP of 58 kb in size harboured by 48 Enteritidis isolates hybridized weakly. One Enteritidis isolate of phage type 14b which failed to hybridize pefA in dot (colony) hybridization experiments was demonstrated to be plasmid free. A cosmid library of Enteritidis phage type 4 expressed in Escherichia coli K12 was screened by hybridization for the presence of pef sequences. Recombinant clones which were deduced to harbour the entire pef operon elaborated a PEF-like fimbrial structure at the cell surface. The PEF-like fimbrial antigen was purified from one cosmid clone and used in western blot experiments with sera from chickens infected with Enteritidis phage-type 4. Seroconversion to the fimbrial antigen was observed which indicated that the Enteritidis PEF-like fimbrial structure was expressed at some stage during infection. Nucleotide sequence analysis demonstrated that the pefA alleles of Typhimurium and Enteritidis phage-type 4 shared 76% DNA nucleotide and 82% deduced amino acid sequence identity.

Brucella melitensis 16M: characterisation of the galE gene and mouse immunisation studies with a galE deficient mutant

The galE gene of Streptomyces lividans was used to probe a cosmid library harbouring Brucella melitensis 16M DNA and the nucleotide sequence of a 2.5 kb ClaI fragment which hybridised was determined. An open reading frame encoding a predicted polypeptide with significant homology to UDP-galactose-4-epimerases of Brucella arbortus strain 2308 and other bacterial species was identified. DNA sequences flanking the B. melitensis galE gene shared no identity with other gal genes and, as for B. abortus, were located adjacent to a mazG homologue. A plasmid which encoded the B. melitensis galE open reading frame complemented a galE mutation in Salmonella typhimurium LB5010, as shown by the restoration of smooth lipopolysaccharide (LPS) biosynthesis, sensitivity to phage P22 infection and restoration of UDP-galactose-4-epimerase activity. The galE gene on the B. melitensis 16M chromosome was disrupted by insertional inactivation and these mutants lacked UDP-galactose-4-epimerase activity but no discernible differences in LPS structure between parent and the mutants were observed. One B. melitensis 16M galE mutant, Bm92, was assessed for virulence in CD-1 and BALB/c mice and displayed similar kinetics of invasion and persistence in tissues compared with the parent bacterial strain. CD-1 mice immunised with B. melitensis 16M galE were protected against B. melitensis 16M challenge. Crown Copyright (C) 1999 Published by Elsevier Science B.V.

Baculovirus superinfection: a probable restriction factor on the surface display of proteins for library screening

In addition to the expression of recombinant proteins, baculoviruses have been developed as a platform for the display of complex eukaryotic proteins on the surface of virus particles or infected insect cells. Surface display has been used extensively for antigen presentation and targeted gene delivery but is also a candidate for the display of protein libraries for molecular screening. However, although baculovirus gene libraries can be efficiently expressed and displayed on the surface of insect cells, target gene selection is inefficient probably due to super-infection which gives rise to cells expressing more than one protein. In this report baculovirus superinfection of Sf9 cells has been investigated by the use of two recombinant multiple nucleopolyhedrovirus carrying green or red fluorescent proteins under the control of both early and late promoters (vAcBacGFP and vAcBacDsRed). The reporter gene expression was detected 8 hours after the infection of vAcBacGFP and cells in early and late phases of infection could be distinguished by the fluorescence intensity of the expressed protein. Simultaneous infection with vAcBacGFP and vAcBacDsRed viruses each at 0.5 MOI resulted in 80% of infected cells coexpressing the two fluorescent proteins at 48 hours post infection (hpi), and subsequent infection with the two viruses resulted in similar co-infection rate. Most Sf9 cells were re-infectable within the first several hours post infection, but the reinfection rate then decreased to a very low level by 16 hpi. Our data demonstrate that Sf9 cells were easily super-infectable during baculovirus infection, and super-infection could occur simultaneously at the time of the primary infection or subsequently during secondary infection by progeny viruses. The efficiency of super-infection may explain the difficulties of baculovirus display library screening but would benefit the production of complex proteins requiring co-expression of multiple polypeptides.

Genome dynamics explain the evolution of flowering time CCT domain gene families in the Poaceae

Numerous CCT domain genes are known to control flowering in plants. They belong to the CONSTANS-like (COL) and PREUDORESPONSE REGULATOR (PRR) gene families, which in addition to a CCT domain possess B-box or response-regulator domains, respectively. Ghd7 is the most recently identified COL gene to have a proven role in the control of flowering time in the Poaceae. However, as it lacks B-box domains, its inclusion within the COL gene family, technically, is incorrect. Here, we show Ghd7 belongs to a larger family of previously uncharacterized Poaceae genes which possess just a single CCT domain, termed here CCT MOTIF FAMILY (CMF) genes. We molecularly describe the CMF (and related COL and PRR) gene families in four sequenced Poaceae species, as well as in the draft genome assembly of barley (Hordeum vulgare). Genetic mapping of the ten barley CMF genes identified, as well as twelve previously unmapped HvCOL and HvPRR genes, finds the majority map to colinear positions relative to their Poaceae orthologues. Combined inter-/intra-species comparative and phylogenetic analysis of CMF, COL and PRR gene families indicates they evolved prior to the monocot/dicot divergence ~200 mya, with Poaceae CMF evolution described as the interplay between whole genome duplication in the ancestral cereal, and subsequent clade-specific mutation, deletion and duplication events. Given the proven role of CMF genes in the modulation of cereals flowering, the molecular, phylogenetic and comparative analysis of the Poaceae CMF, COL and PRR gene families presented here provides the foundation from which functional investigation can be undertaken.

Construction of a Festuca pratensis BAC library for map-based cloning in Festulolium substitution lines

Introgression in Festulolium is a potentially powerful tool to isolate genes for a large number of traits which differ between Festuca pratensis Huds. and Lolium perenne L. Not only are hybrids between the two species fertile, but the two genomes can be distinguished by genomic in situ hybridisation and a high frequency of recombination occurs between homoeologous chromosomes and chromosome segments. By a programme of introgression and a series of backcrosses, L. perenne lines have been produced which contain small F. pratensis substitutions. This material is a rich source of polymorphic markers targeted towards any trait carried on the F. pratensis substitution not observed in the L. perenne background. We describe here the construction of an F. pratensis BAC library, which establishes the basis of a map-based cloning strategy in L. perenne. The library contains 49,152 clones, with an average insert size of 112 kbp, providing coverage of 2.5 haploid genome equivalents. We have screened the library for eight amplified fragment length polymorphism (AFLP) derived markers known to be linked to an F. pratensis gene introgressed into L. perenne and conferring a staygreen phenotype as a consequence of a mutation in primary chlorophyll catabolism. While for four of the markers it was possible to identify bacterial artificial chromosome (BAC) clones, the other four AFLPs were too repetitive to enable reliable identification of locus-specific BACs. Moreover, when the four BACs were partially sequenced, no obvious coding regions could be identified. This contrasted to BACs identified using cDNA sequences, when multiple genes were identified on the same BAC.

A maize bacterial artificial chromosome (BAC) library from the European flint inbred line F2

We report here the construction and characterisation of a BAC library from the maize flint inbred line F2, widely used in European maize breeding programs. The library contains 86,858 clones with an average insert size of approximately 90 kb, giving approximately 3.2-times genome coverage. High-efficiency BAC cloning was achieved through the use of a single size selection for the high-molecular-weight genomic DNA, and co-transformation of the ligation with yeast tRNA to optimise transformation efficiency. Characterisation of the library showed that less than 0.5% of the clones contained no inserts, while 5.52% of clones consisted of chloroplast DNA. The library was gridded onto 29 nylon filters in a double-spotted 8 × 8 array, and screened by hybridisation with a number of single-copy and gene-family probes. A 3-dimensional DNA pooling scheme was used to allow rapid PCR screening of the library based on primer pairs from simple sequence repeat (SSR) and expressed sequence tag (EST) markers. Positive clones were obtained in all hybridisation and PCR screens carried out so far. Six BAC clones, which hybridised to a portion of the cloned Rp1-D rust resistance gene, were further characterised and found to form contigs covering most of this complex resistance locus.

Modulation of stretch-induced myocyte remodeling and gene expression by nitric oxide: a novel role for lipoma preferred partner in myofibrillogenesis

Prolonged hemodynamic load as a result of hypertension eventually leads to maladaptive cardiac adaptation and heart failure. The signalling pathways that underlie these changes are still poorly understood. The adaptive response to mechanical load is mediated by mechanosensors which convert the mechanical stimuli into a biological response. We examined the effect of cyclic mechanical stretch on myocyte adaptation using neonatal rat ventricular myocytes with 10% (adaptive) or 20% (maladaptive) maximum strain, 1Hz for 48 hours to mimic in vivo mechanical stress. Cells were also treated with and without L-NAME, a general nitric oxide synthase (NOS) inhibitor to suppress NO production. Maladaptive 20% mechanical stretch led to a significant loss of intact sarcomeres which was rescued by LNAME (P<0.05, n≥5 cultures). We hypothesized that the mechanism was through NOinduced alteration of myocyte gene expression. L-NAME up-regulated the mechanosensing proteins Muscle LIM protein (MLP (by 100%, p<0.05, n=4 cultures)) and lipoma preferred partner, a novel cardiac protein (LPP (by 80%, p<0.05, n=4 cultures)). L-NAME also significantly altered the subcellular localisation of LPP and MLP in a manner that favoured growth and adaptation. These findings suggest that NO participates in stretch-mediated adaptation. The use of isoform selective NOS inhibitors indicated a complex interaction between iNOS and nNOS isoforms regulate gene expression. LPP knockdown by siRNA led to formation of α-actinin aggregates and Z-bodies showing that myofibrillogenesis was impaired. There was an up-regulation of E3 ubiquitin ligase (MUL1) by 75% (P<0.05, n=5 cultures). This indicates that NO contributes to stretch-mediated adaptation via the upregulation of proteins associated mechansensing and myofibrillogenesis, thereby presenting potential therapeutic targets during the progression of heart failure. Keywords: Mechanotransduction, heart failure, stretch, heart, hypertrophy

Gene expression changes in phosphorus deficient potato (Solanum tuberosum L.) leaves and the potential for diagnostic gene expression markers

Background: There are compelling economic and environmental reasons to reduce our reliance on inorganic phosphate (Pi) fertilisers. Better management of Pi fertiliser applications is one option to improve the efficiency of Pi fertiliser use, whilst maintaining crop yields. Application rates of Pi fertilisers are traditionally determined from analyses of soil or plant tissues. Alternatively, diagnostic genes with altered expression under Pi limiting conditions that suggest a physiological requirement for Pi fertilisation, could be used to manage Pifertiliser applications, and might be more precise than indirect measurements of soil or tissue samples. Results: We grew potato (Solanum tuberosum L.) plants hydroponically, under glasshouse conditions, to control their nutrient status accurately. Samples of total leaf RNA taken periodically after Pi was removed from the nutrient solution were labelled and hybridised to potato oligonucleotide arrays. A total of 1,659 genes were significantly differentially expressed following Pi withdrawal. These included genes that encode proteins involved in lipid, protein, and carbohydrate metabolism, characteristic of Pi deficient leaves and included potential novel roles for genes encoding patatin like proteins in potatoes. The array data were analysed using a support vector machine algorithm to identify groups of genes that could predict the Pi status of the crop. These groups of diagnostic genes were tested using field grown potatoes that had either been fertilised or unfertilised. A group of 200 genes could correctly predict the Pi status of field grown potatoes. Conclusions: This paper provides a proof-of-concept demonstration for using microarrays and class prediction tools to predict the Pi status of a field grown potato crop. There is potential to develop this technology for other biotic and abiotic stresses in field grown crops. Ultimately, a better understanding of crop stresses may improve our management of the crop, improving the sustainability of agriculture.

Modelling negative feedback networks for activating transcription factor 3 predicts a dominant role for miRNAs in immediate early gene regulation

Activating transcription factor 3 (Atf3) is rapidly and transiently upregulated in numerous systems, and is associated with various disease states. Atf3 is required for negative feedback regulation of other genes, but is itself subject to negative feedback regulation possibly by autorepression. In cardiomyocytes, Atf3 and Egr1 mRNAs are upregulated via ERK1/2 signalling and Atf3 suppresses Egr1 expression. We previously developed a mathematical model for the Atf3-Egr1 system. Here, we adjusted and extended the model to explore mechanisms of Atf3 feedback regulation. Introduction of an autorepressive loop for Atf3 tuned down its expression and inhibition of Egr1 was lost, demonstrating that negative feedback regulation of Atf3 by Atf3 itself is implausible in this context. Experimentally, signals downstream from ERK1/2 suppress Atf3 expression. Mathematical modelling indicated that this cannot occur by phosphorylation of pre-existing inhibitory transcriptional regulators because the time delay is too short. De novo synthesis of an inhibitory transcription factor (ITF) with a high affinity for the Atf3 promoter could suppress Atf3 expression, but (as with the Atf3 autorepression loop) inhibition of Egr1 was lost. Developing the model to include newly-synthesised miRNAs very efficiently terminated Atf3 protein expression and, with a 4-fold increase in the rate of degradation of mRNA from the mRNA/miRNA complex, profiles for Atf3 mRNA, Atf3 protein and Egr1 mRNA approximated to the experimental data. Combining the ITF model with that of the miRNA did not improve the profiles suggesting that miRNAs are likely to play a dominant role in switching off Atf3 expression post-induction.

Identification of novel cholesteatoma-related gene expression signatures using full-genome microarrays

This study demonstrates that the expression profile of cholesteatoma is similar to a metastatic tumour and chronically inflamed tissue. Based on the investigated profiles we present novel protein-protein interaction and signal transduction networks, which include cholesteatoma-regulated transcripts and may be of great value for drug targeting and therapy development.