Gene Expression Complex Networks: Synthesis, Identification, and Analysis

Thanks to recent advances in molecular biology, allied to an ever increasing amount of experimental data, the functional state of thousands of genes can now be extracted simultaneously by using methods such as cDNA microarrays and RNA-Seq. Particularly important related investigations are the modeling and identification of gene regulatory networks from expression data sets. Such a knowledge is fundamental for many applications, such as disease treatment, therapeutic intervention strategies and drugs design, as well as for planning high-throughput new experiments. Methods have been developed for gene networks modeling and identification from expression profiles. However, an important open problem regards how to validate such approaches and its results. This work presents an objective approach for validation of gene network modeling and identification which comprises the following three main aspects: (1) Artificial Gene Networks (AGNs) model generation through theoretical models of complex networks, which is used to simulate temporal expression data; (2) a computational method for gene network identification from the simulated data, which is founded on a feature selection approach where a target gene is fixed and the expression profile is observed for all other genes in order to identify a relevant subset of predictors; and (3) validation of the identified AGN-based network through comparison with the original network. The proposed framework allows several types of AGNs to be generated and used in order to simulate temporal expression data. The results of the network identification method can then be compared to the original network in order to estimate its properties and accuracy. Some of the most important theoretical models of complex networks have been assessed: the uniformly-random Erdos-Renyi (ER), the small-world Watts-Strogatz (WS), the scale-free Barabasi-Albert (BA), and geographical networks (GG). The experimental results indicate that the inference method was sensitive to average degree k variation, decreasing its network recovery rate with the increase of k. The signal size was important for the inference method to get better accuracy in the network identification rate, presenting very good results with small expression profiles. However, the adopted inference method was not sensible to recognize distinct structures of interaction among genes, presenting a similar behavior when applied to different network topologies. In summary, the proposed framework, though simple, was adequate for the validation of the inferred networks by identifying some properties of the evaluated method, which can be extended to other inference methods.

Functional microarray analysis suggests repressed cell-cell signaling and cell survival-related modules inhibit progression of head and neck squamous cell carcinoma

Background: Cancer shows a great diversity in its clinical behavior which cannot be easily predicted using the currently available clinical or pathological markers. The identification of pathways associated with lymph node metastasis (N+) and recurrent head and neck squamous cell carcinoma (HNSCC) may increase our understanding of the complex biology of this disease. Methods: Tumor samples were obtained from untreated HNSCC patients undergoing surgery. Patients were classified according to pathologic lymph node status (positive or negative) or tumor recurrence (recurrent or non-recurrent tumor) after treatment (surgery with neck dissection followed by radiotherapy). Using microarray gene expression, we screened tumor samples according to modules comprised by genes in the same pathway or functional category. Results: The most frequent alterations were the repression of modules in negative lymph node (N0) and in non-recurrent tumors rather than induction of modules in N+ or in recurrent tumors. N0 tumors showed repression of modules that contain cell survival genes and in non-recurrent tumors cell-cell signaling and extracellular region modules were repressed. Conclusions: The repression of modules that contain cell survival genes in N0 tumors reinforces the important role that apoptosis plays in the regulation of metastasis. In addition, because tumor samples used here were not microdissected, tumor gene expression data are represented together with the stroma, which may reveal signaling between the microenvironment and tumor cells. For instance, in non-recurrent tumors, extracellular region module was repressed, indicating that the stroma and tumor cells may have fewer interactions, which disable metastasis development. Finally, the genes highlighted in our analysis can be implicated in more than one pathway or characteristic, suggesting that therapeutic approaches to prevent tumor progression should target more than one gene or pathway, specially apoptosis and interactions between tumor cells and the stroma.

The Penicillium chrysogenum aclA gene encodes a broad-substrate-specificity acyl-coenzyme A ligase involved in activation of adipic acid, a side-chain precursor for cephem antibiotics

Activation of the cephalosporin side-chain precursor to the corresponding CoA-thioester is an essential step for its incorporation into the P-lactam backbone. To identify an acyl-CoA ligase involved in activation of adipate, we searched in the genome database of Penicillium chrysogenum for putative structural genes encoding acyl-CoA ligases. Chemostat-based transcriptome analysis was used to identify the one presenting the highest expression level when cells were grown in the presence of adipate. Deletion of the gene renamed aclA, led to a 32% decreased specific rate of adipate consumption and a threefold reduction of adipoyl-6-aminopenicillanic acid levels, but did not affect penicillin V production. After overexpression in Escherichia coli, the purified protein was shown to have a broad substrate range including adipate. Finally, protein-fusion with cyan-fluorescent protein showed co-localization with microbody-borne acyl-transferase. Identification and functional characterization of aclA may aid in developing future metabolic engineering strategies for improving the production of different cephalosporins. (C) 2009 Elsevier Inc. All rights reserved.

Extension of the core map of common bean with EST-SSR, RGA, AFLP, and putative functional markers

Microsatellites and gene-derived markers are still underrepresented in the core molecular linkage map of common bean compared to other types of markers. In order to increase the density of the core map, a set of new markers were developed and mapped onto the RIL population derived from the `BAT93` x `Jalo EEP558` cross. The EST-SSR markers were first characterized using a set of 24 bean inbred lines. On average, the polymorphism information content was 0.40 and the mean number of alleles per locus was 2.7. In addition, AFLP and RGA markers based on the NBS-profiling method were developed and a subset of the mapped RGA was sequenced. With the integration of 282 new markers into the common bean core map, we were able to place markers with putative known function in some existing gaps including regions with QTL for resistance to anthracnose and rust. The distribution of the markers over 11 linkage groups is discussed and a newer version of the common bean core linkage map is proposed.

Cloning, characterization and functional analysis of a 1-FEH cDNA from Vernonia herbacea (Vell.) Rusby

Variations in the inulin contents have been detected in rhizophores of Vernonia herbacea during the phenological cycle. These variations indicate the occurrence of active inulin synthesis and depolymerization throughout the cycle and a role for this carbohydrate as a reserve compound. 1-Fructan exohydrolase (1-FEH) is the enzyme responsible for inulin depolymerization, and its activity has been detected in rhizophores of sprouting plants. Defoliation and low temperature are enhancer conditions of this 1-FEH activity. The aim of the present work was the cloning of this enzyme. Rhizophores were collected from plants induced to sprout, followed by storage at 5C. A full length 1-FEH cDNA sequence was obtained by PCR and inverse PCR techniques, and expressed in Pichia pastoris. Cold storage enhances FEH gene expression. Vh1-FEH was shown to be a functional 1-FEH, hydrolyzing predominantly -2,1 linkages, sharing high identity with chicory FEH sequences, and its activity was inhibited by 81 in the presence of 10 mM sucrose. In V. herbacea, low temperature and sucrose play a role in the control of fructan degradation. This is the first study concerning the cloning and functional analysis of a 1-FEH cDNA of a native species from the Brazilian Cerrado. Results will contribute to understanding the role of fructans in the establishment of a very successful fructan flora of the Brazilian Cerrado, subjected to water limitation and low temperature during winter.

Functional characterization of the Aspergillus fumigatus PHO80 homologue

Phosphate is an ion that is essential for fungal growth. The systems for inorganic phosphate (Pi) acquisition in eukaryotic cells (PHO) have been characterized as a low-affinity (that assures a supply of Pi at normal or high external Pi concentrations) and a high-affinity (activated in response to Pi starvation). Here, as an initial step to understand the PHO pathway in Aspergillus fumigatus, we characterized the PH080 homologue, PhoB(PHO80). We show that the Delta phoB(PHO80) mutant has a polar growth defect (i.e., a delayed germ tube emergence) and, by phenotypic and phosphate uptake analyses, establish a link between PhoB(PHO80), calcineurin and calcium metabolism. Microarray hybridizations carried out with RNA obtained from wild-type and Delta phoB(PHO80) mutant cells identify Afu4g03610 (phoD(PHO84)), Afu7g06350 (phoE(PHO89)), Afu4g06020 (phoC(PHO81)), and Afu2g09040 (vacuolar transporter Vtc4) as more expressed both in the Delta phoB(PHO80) mutant background and under phosphate-limiting conditions of 0.1 mM P-i. Epifluorescence microscopy revealed accumulation of poly-phosphate in Delta phoB(PHO80) vacuoles, which was independent of extracellular phosphate concentration. Surprisingly, a phoD(PHO84) deletion mutant is indistinguishable phenotypically from the corresponding wild-type strain. mRNA analyses suggest that protein kinase A absence supports the expression of PHO genes in A. fumigatus. Furthermore, Delta phoB(PHO80) and Delta phoD(PHO84) mutant are fully virulent in a murine low dose model for invasive aspergillosis. (C) 2008 Elsevier Inc. All rights reserved.

Functional characterization of the Aspergillus fumigatus CRZ1 homologue, CrzA

The protein phosphatase calcineurin is an important mediator connecting calcium-dependent signalling to various cellular responses in multiple organisms. In fungi calcineurin acts largely through regulating Crz1p-like transcription factors. Here we characterize an Aspergillus fumigatus CRZ1 homologue, CrzA and demonstrate its mediation of cellular tolerance to increased concentrations of calcium and manganese. In addition to acute sensitivitiy to these ions, and decreased conidiation, the crzA null mutant suffers altered expression of calcium transporter mRNAs under high concentrations of calcium, and loss of virulence when compared with the corresponding complemented and wild-type strains. We use multiple expression analyses to probe the transcriptional basis of A. fumigatus calcium tolerance identifying several genes having calA and/or crzA dependent mRNA accumulation patterns. We also demonstrate that contrary to previous findings, the gene encoding the Aspergillus nidulans calcineurin subunit homologue, cnaA, is not essential and that the cnaA deletion mutant shares the morphological phenotypes observed in the corresponding A. fumigatus mutant, Delta calA. Exploiting the A. nidulans model system, we have linked calcineurin activity with asexual developmental induction, finding that CrzA supports appropriate developmental induction in a calcineurin and brlA-dependent manner in both species.

Functional characterization of the putative Aspergillus nidulans DNA damage binding protein homologue DdbA

Nucleotide excision repair (NER) eliminates helix-distorting DNA base lesions. Seven XP-deficient genetic complementation groups (XPA to XPG) have already been identified in mammals, and their corresponding genes have been cloned. Hereditary defects in NER are associated with several diseases, including xeroderma pigmentosum (XP). UV-DDB (XPE) is formed by two associated subunits, DDB1 and DDB2. UV-DDB was identified biochemically as a protein factor that exhibits very strong and specific binding to ultraviolet (UV)-treated DNA. As a preliminary step to characterize the components of the NER in the filamentous fungus Aspergillus nidulans, here we identified a putative DDB1 homologue, DdbA. Deletion and expression analysis indicated that A. nidulans ddbA gene is involved in the DNA damage response, more specifically in the UV light response and 4-nitroquinoline oxide (4-NQO) sensitivity. Furthermore, the Delta ddbA strain cannot self-cross and expression analysis showed that ddbA can be induced by oxidative stress and is developmentally regulated in both asexual and sexual processes. The Delta ddbA mutation can genetically interact with uvsB(ATR), atmA(ATM), nkuA(KU70), H2AX-S129A (a replacement of the conserved serine in the C-terminal of H2AX with alanine), and cshB (a mutation in CSB Cockayne`s syndrome protein involved in the transcription-coupled repair subpathway of NER) mutations. Finally, to determine the DdbA cellular localization, we constructed a GFP:DdbA strain. In the presence and absence of DNA damage, DdbA was mostly detected in the nuclei, indicating that DdbA localizes to nuclei and its cellular localization is not affected by the cellular response to DNA damage induced by 4-NQO and UV light.

A functional matrix metalloproteinase (MMP)-9 polymorphism modifies plasma MMP-9 levels in subjects environmentally exposed to mercury

Mercury (Hg) exposure causes health problems including cardiovascular diseases. Although precise mechanisms have not been precisely defined yet, matrix metalloproteinases (MMPs) may be involved. The gene encoding MMP-9 presents genetic polymorphisms which affect the expression and activity level of this enzyme. Two polymorphisms in the promoter region [C(-1562)T and (CA)(n)] are functionally relevant, and are implicated in several diseases. This study aimed at examining how these polymorphisms affect the circulating MMP-9 levels and its endogenous inhibitor, the tissue inhibitor of metalloproteinase-1 (TIMP-1) in 266 subjects environmentally exposed to Hg. Blood and plasma Hg concentrations were determined by inductively coupled plasma-mass spectrometry (ICP-MS). MMP-9 and TIMP-1 concentrations were measured in plasma samples by gelatin zymography and ELISA, respectively. Genotypes for the C(-1562)T and the microsatellite (CA)(n) polymorphisms were determined. We found a positive association (P<0.05) between plasma Hg concentrations and MMP-9/TIMP-1 ratio (an index of net MMP-9 activity). When the subjects were divided into tertiles with basis on their plasma Hg concentrations, we found that the (CA)(n) polymorphism modified MMP-9 concentrations and MMP-9/TIMP-1 ratio in subjects with the lowest Hg concentrations (first tertile), with the highest MMP-9 levels being found in subjects with genotypes including alleles with 21 or more CA repeats (H alleles) (P<0.05). Conversely, this polymorphism had no effects on subjects with intermediate or high plasma Hg levels (second and third tertiles, respectively). The C(-1562)T polymorphism had no effects on MMP-9 levels. These findings suggest a significant interaction between the (CA)(n) polymorphism and low levels of Hg exposure, possibly increasing the risk of developing diseases in subjects with H alleles. (c) 2010 Elsevier B.V. All rights reserved.

X-linked myoclonic epilepsy with spasticity and intellectual disability - Mutation in the homeobox gene ARX

Objective: To describe a new syndrome of X-linked myoclonic epilepsy with generalized spasticity and intellectual disability (XMESID) and identify the gene defect underlying this disorder. Methods: The authors studied a family in which six boys over two generations had intractable seizures using a validated seizure questionnaire, clinical examination, and EEG studies. Previous records and investigations were obtained. Information on seizure disorders was obtained on 271 members of the extended family. Molecular genetic analysis included linkage studies and mutational analysis using a positional candidate gene approach. Results: All six affected boys had myoclonic seizures and TCS; two had infantile spasms, but only one had hypsarrhythmia. EEG studies show diffuse background slowing with slow generalized spike wave activity. All affected boys had moderate to profound intellectual disability. Hyperreflexia was observed in obligate carrier women. A late-onset progressive spastic ataxia in the matriarch raises the possibility of late clinical manifestations in obligate carriers. The disorder was mapped to Xp11.2-22.2 with a maximum lod score of 1.8. As recently reported, a missense mutation (1058C>T/P353L) was identified within the homeodomain of the novel human Aristaless related homeobox gene (ARX). Conclusions: XMESID is a rare X-linked recessive myoclonic epilepsy with spasticity and intellectual disability in boys. Hyperreflexia is found in carrier women. XMESID is associated with a missense mutation in ARX. This disorder is allelic with X-linked infantile spasms (ISSX; MIM 308350) where polyalanine tract expansions are the commonly observed molecular defect. Mutations of ARX are associated with a wide range of phenotypes; functional studies in the future may lend insights to the neurobiology of myoclonic seizures and infantile spasms.

Regulation of the gene encoding glutathione S-transferase M1 (GSTM1) by the Myb oncoprotein

The identification of Myb 'target' genes will not only aid in the understanding of how overexpression of Myb, or expression of activated forms of Myb, leads to cellular transformation but will also shed light on its role in normal cells. Using a combination of an estrogen-regulated Myb-transformed cell line (ERMYB) and PCR-based subtractive hybridization, we have identified the gene (GSTM1) encoding the detoxification enzyme glutathione S-transferase M1 as being transcriptionally upregulated by Myb. Functional analysis of the GSTM1 promoter using reporter assays indicated that both the DNA binding and transactivation domains of Myb were required for transcriptional activation. Mutational analysis of consensus Myb-binding sites (MBS) in the promoter and electrophoretic mobility gel shift analysis indicated that one of the three potential MBS can bind Myb protein, and is the primary site involved in the regulation of this promoter by Myb.

Identification of a minimal promoter sequence for the human N-acetyltransferase Type I gene that binds AP-1 (activator protein 1) and YY-1 (Yin and Yang 1)

Human N-acetyltransferase Type I (NAT1) catalyses the acetylation of many aromatic amine and hydrazine compounds and it has been implicated in the catabolism of folic acid. The enzyme is widely expressed in the body, although there are considerable differences in the level of activity between tissues. A search of the mRNA databases revealed the presence of several NAT1 transcripts in human tissue that appear to be derived from different promoters. Because little is known about NAT1 gene regulation, the present study was undertaken to characterize one of the putative promoter sequences of the NAT1 gene located just upstream of the coding region. We show with reverse-transcriptase PCR that mRNA transcribed from this promoter (Promoter 1) is present in a variety of human cell-lines, but not in quiescent peripheral blood mononuclear cells. Using deletion mutant constructs, we identified a 20 bp sequence located 245 bases upstream of the translation start site which was sufficient for basal NAT1 expression. It comprised an AP-1 (activator protein 1)-binding site, flanked on either side by a TCATT motif. Mutational analysis showed that the AP-1 site and the 3' TCATT sequence were necessary for gene expression, whereas the 5' TCATT appeared to attenuate promoter activity. Electromobility shift assays revealed two specific bands made up by complexes of c-Fos/Fra, c-Jun, YY-1 (Yin and Yang 1) and possibly Oct-1. PMA treatment enhanced expression from the NAT1 promoter via the AP-1-binding site. Furthermore, in peripheral blood mononuclear cells, PMA increased endogenous NAT1 activity and induced mRNA expression from Promoter I, suggesting that it is functional in vivo.

Functional and Evolutionary Insights from the Genomes of Three Parasitoid Nasonia Species

We report here genome sequences and comparative analyses of three closely related parasitoid wasps: Nasonia vitripennis, N. giraulti, and N. longicornis. Parasitoids are important regulators of arthropod populations, including major agricultural pests and disease vectors, and Nasonia is an emerging genetic model, particularly for evolutionary and developmental genetics. Key findings include the identification of a functional DNA methylation tool kit; hymenopteran-specific genes including diverse venoms; lateral gene transfers among Pox viruses, Wolbachia, and Nasonia; and the rapid evolution of genes involved in nuclear-mitochondrial interactions that are implicated in speciation. Newly developed genome resources advance Nasonia for genetic research, accelerate mapping and cloning of quantitative trait loci, and will ultimately provide tools and knowledge for further increasing the utility of parasitoids as pest insect-control agents.

A non-invasive method for silencing gene transcription in honeybees maintained under natural conditions

in the Apis mellifera post-genomic era, RNAi protocols have been used in functional approaches. However, sample manipulation and invasive methods such as injection of double-stranded RNA (dsRNA) can compromise physiology and survival. To circumvent these problems, we developed a non-invasive method for honeybee gene knockdown, using a well-established vitellogenin RNAi system as a model. Second instar larvae received dsRNA for vitellogenin (dsVg-RNA) in their natural diet. For exogenous control, larvae received dsRNA for GFP (dsGFP-RNA). Untreated larvae formed another control group. Around 60% of the treated larvae naturally developed until adult emergence when 0.5 mu g of dsVg-RNA or dsGFP-RNA was offered while no larvae that received 3.0 mu g of dsRNA reached pupal stages. Diet dilution did not affect the removal rates. Viability depends not only on the delivered doses but also on the internal conditions of colonies. The weight of treated and untreated groups showed no statistical differences. This showed that RNAi ingestion did not elicit drastic collateral effects. Approximately 90% of vitellogenin transcripts from 7-day-old workers were silenced compared to controls. A large number of samples are handled in a relatively short time and smaller quantities of RNAi molecules are used compared to invasive methods. These advantages culminate in a versatile and a cost-effective approach. (c) 2008 Elsevier Ltd. All rights reserved.

The p.A2215D Thyroglobulin Gene Mutation Leads to Deficient Synthesis and Secretion of the Mutated Protein and Congenital Hypothyroidism with Wide Phenotype Variation

Context: Thyroglobulin (TG) is a large glycoprotein and functions as a matrix for thyroid hormone synthesis. TG gene mutations give rise to goitrous congenital hypothyroidism (CH) with considerable phenotype variation. Objectives: The aim of the study was to report the genetic screening of 15 patients with CH due to TG gene mutations and to perform functional analysis of the p. A2215D mutation. Design: Clinical evaluation and DNA sequencing of the TG gene were performed in all patients. TG expression was analyzed in the goitrous tissue of one patient. Human cells were transfected with expression vectors containing mutated and wild-type human TG cDNA. Results: All patients had an absent rise of serum TG after stimulation with recombinant human TSH. Sequence analysis revealed three previously described mutations (p. A2215D, p. R277X, and g. IVS30 + 1G > T), and two novel mutations (p. Q2142X and g. IVS46-1G > A). Two known (g. IVS30 + 1G/p. A2215D and p. A2215D/p. R277X) and one novel (p. R277X/g. IVS46-1G > A) compound heterozygous constellations were also identified. Functional analysis indicated deficiency in TG synthesis, reduction of TG secretion, and retention of the mutant TG within the cell, leading to an endoplasmic reticulum storage disease, whereas small amounts of mutant TG were still secreted within the cell system. Conclusion: All studied patients were either homozygous or heterozygous for TG gene mutations. Two novel mutations have been detected, and we show that TG mutation p. A2215D promotes the retention of TG within the endoplasmic reticulum and reduces TG synthesis and secretion, causing mild hypothyroidism. In the presence of sufficient iodine supply, some patients with TG mutations are able to compensate the impaired hormonogenesis and generate thyroid hormone. (J Clin Endocrinol Metab 94: 2938-2944, 2009)