Network signatures based on gene pair expression ratios improve classification and the analysis of muscle-invasive urothelial cancer

Urothelial cancer (UC) is highly recurrent and can progress from non-invasive (NMIUC) to a more aggressive muscle-invasive (MIUC) subtype that invades the muscle tissue layer of the bladder. We present a proof of principle study that network-based features of gene pairs can be used to improve classifier performance and the functional analysis of urothelial cancer gene expression data. In the first step of our procedure each individual sample of a UC gene expression dataset is inflated by gene pair expression ratios that are defined based on a given network structure. In the second step an elastic net feature selection procedure for network-based signatures is applied to discriminate between NMIUC and MIUC samples. We performed a repeated random subsampling cross validation in three independent datasets. The network signatures were characterized by a functional enrichment analysis and studied for the enrichment of known cancer genes. We observed that the network-based gene signatures from meta collections of proteinprotein interaction (PPI) databases such as CPDB and the PPI databases HPRD and BioGrid improved the classification performance compared to single gene based signatures. The network based signatures that were derived from PPI databases showed a prominent enrichment of cancer genes (e.g., TP53, TRIM27 and HNRNPA2Bl). We provide a novel integrative approach for large-scale gene expression analysis for the identification and development of novel diagnostical targets in bladder cancer. Further, our method allowed to link cancer gene associations to network-based expression signatures that are not observed in gene-based expression signatures.

Disease-Gene Association Using a Genetic Algorithm

Understanding the relationship between genetic diseases and the genes associated with them is an important problem regarding human health. The vast amount of data created from a large number of high-throughput experiments performed in the last few years has resulted in an unprecedented growth in computational methods to tackle the disease gene association problem. Nowadays, it is clear that a genetic disease is not a consequence of a defect in a single gene. Instead, the disease phenotype is a reflection of various genetic components interacting in a complex network. In fact, genetic diseases, like any other phenotype, occur as a result of various genes working in sync with each other in a single or several biological module(s). Using a genetic algorithm, our method tries to evolve communities containing the set of potential disease genes likely to be involved in a given genetic disease. Having a set of known disease genes, we first obtain a protein-protein interaction (PPI) network containing all the known disease genes. All the other genes inside the procured PPI network are then considered as candidate disease genes as they lie in the vicinity of the known disease genes in the network. Our method attempts to find communities of potential disease genes strongly working with one another and with the set of known disease genes. As a proof of concept, we tested our approach on 16 breast cancer genes and 15 Parkinson's Disease genes. We obtained comparable or better results than CIPHER, ENDEAVOUR and GPEC, three of the most reliable and frequently used disease-gene ranking frameworks.

Gene and species trees of a Neotropical group of treefrogs: Genetic diversification in the Brazilian Atlantic Forest and the origin of a polyploid species

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Evolução e organização genômica do gene FOXL2 em vertebrados

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Reciprocal gene flow between conventional and genetically modified soybean cultivars

The objective of this work was to determine the reciprocal gene flow between two soybean cultivars, one tolerant and the other sensitive to glyphosate, as well as to use estimators to determine the outcrossing rate in the population and the number of hybrid seeds in the progeny. The experiment was composed of four blocks of 40 soybean rows, of which 20 rows of each cultivar (CD217 and CD219RR). At the R8 stage, five rows, distant 0.5, 1, 2, 4 and 8 m from the adjacent cultivar, were harvested, threshed and analyzed as for the occurrence of gene flow. As phenotypical markers, the trait color of flowers, hypocotil and pubescence, as well as the tolerance to glyphosate were used. The cultivars contrast for all the analyzed traits, each one conditioned by a single gene with two alleles, in a complete dominance interaction. In the tolerant cultivar progeny, the largest outcross rate was 0.27%, and in the sensitive cultivar progeny, 0.83% was identified; by the dilution effect hypothesis, the outcross rates in natural populations would be 0.104 and 0.388%, respectively. The reciprocal gene flow between CD217 and CD219RR cultivars is not the same in both directions. The proposed estimators are useful for determining the hybrid rates in seed samples.

Genetic epidemiology of gallbladder disease in Mexican Americans and cholesterol 7$\alpha$-hydroxylase gene sequence variation

Among Mexican Americans, the second largest minority group in the United States, the prevalence of gallbladder disease is markedly elevated. Previous data from both genetic admixture and family studies indicate that there is a genetic component to the occurrence of gallbladder disease in Mexican Americans. However, prior to this thesis no formal genetic analysis of gallbladder disease had been carried out nor had any contributing genes been identified.^ The results of complex segregation analysis in a sample of 232 Mexican American pedigrees documented the existence of a major gene having two alleles with age- and gender-specific effects influencing the occurrence of gallbladder disease. The estimated frequency of the allele increasing susceptibility was 0.39. The lifetime probabilities that an individual will be affected by gallbladder disease were 1.0, 0.54, and 0.00 for females of genotypes "AA", "Aa", and "aa", respectively, and 0.68, 0.30, and 0.00 for males, respectively. This analysis provided the first conclusive evidence for the existence of a common single gene having a large effect on the occurrence of gallbladder disease.^ Human cholesterol 7$\alpha$-hydroxylase is the rate-limiting enzyme in bile acid synthesis. The results of an association study in both a random sample and a matched case/control sample showed that there is a significant association between cholesterol 7$\alpha$-hydroxylase gene variation and the occurrence of gallbladder disease in Mexican Americans males but not in females. These data have implicated a specific gene, 7$\alpha$-hydroxylase, in the etiology of gallbladder disease in this population.^ Finally, I asked whether the inferred major gene from complex segregation analysis is genetically linked to the cholesterol 7$\alpha$-hydroxylase gene. Three pedigrees predicted to be informative for linkage analysis by virtue of supporting the major gene hypothesis and having parents with informative genotypes and multiple offspring were selected for this linkage analysis. In each of these pedigrees, the recombination fractions maximized at 0 with a positive, albeit low, LOD score. The results of this linkage analysis provide preliminary and suggestive evidence that the cholesterol 7$\alpha$-hydroxylase gene and the inferred gallbladder disease susceptibility gene are genetically linked. ^

The Crithidia fasciculata RNH1 gene encodes both nuclear and mitochondrial isoforms of RNase H

The Crithidia fasciculata RNH1 gene encodes an RNase H, an enzyme that specifically degrades the RNA strand of RNA–DNA hybrids. The RNH1 gene is contained within an open reading frame (ORF) predicted to encode a protein of 53.7 kDa. Previous work has shown that RNH1 expresses two proteins: a 38 kDa protein and a 45 kDa protein which is enriched in kinetoplast extracts. Epitope tagging of the C-terminus of the RNH1 gene results in localization of the protein to both the kinetoplast and the nucleus. Translation of the ORF beginning at the second in-frame methionine codon predicts a protein of 38 kDa. Insertion of two tandem stop codons between the first ATG codon and the second in-frame ATG codon of the ORF results in expression of only the 38 kDa protein and the protein localizes specifically to the nucleus. Mutation of the second methionine codon to a valine codon prevents expression of the 38 kDa protein and results in exclusive production of the 45 kDa protein and localization of the protein only in the kinetoplast. These results suggest that the kinetoplast enzyme results from processing of the full-length 53.7 kDa protein. The nuclear enzyme appears to result from translation initiation at the second in-frame ATG codon. This is the first example in trypanosomatids of the production of nuclear and mitochondrial isoforms of a protein from a single gene and is the only eukaryotic gene in the RNase HI gene family shown to encode a mitochondrial RNase H.

VIPP1, a nuclear gene of Arabidopsis thaliana essential for thylakoid membrane formation

The conversion of light to chemical energy by the process of photosynthesis is localized to the thylakoid membrane network in plant chloroplasts. Although several pathways have been described that target proteins into and across the thylakoids, little is known about the origin of this membrane system or how the lipid backbone of the thylakoids is transported and fused with the target membrane. Thylakoid biogenesis and maintenance seem to involve the flow of membrane elements via vesicular transport. Here we show by mutational analysis that deletion of a single gene called VIPP1 (vesicle-inducing protein in plastids 1) is deleterious to thylakoid membrane formation. Although VIPP1 is a hydrophilic protein it is found in both the inner envelope and the thylakoid membranes. In VIPP1 deletion mutants vesicle formation is abolished. We propose that VIPP1 is essential for the maintenance of thylakoids by a transport pathway not previously recognized.

Generation of broad-spectrum disease resistance by overexpression of an essential regulatory gene in systemic acquired resistance

The recently cloned NPR1 gene of Arabidopsis thaliana is a key regulator of acquired resistance responses. Upon induction, NPR1 expression is elevated and the NPR1 protein is activated, in turn inducing expression of a battery of downstream pathogenesis-related genes. In this study, we found that NPR1 confers resistance to the pathogens Pseudomonas syringae and Peronospora parasitica in a dosage-dependent fashion. Overexpression of NPR1 leads to enhanced resistance with no obvious detrimental effect on the plants. Thus, for the first time, a single gene is shown to be a workable target for genetic engineering of nonspecific resistance in plants.

Analysis of Promoter Activity for the Gene Encoding Pyruvate Orthophosphate Dikinase in Stably Transformed C4 Flaveria Species1

The C4 enzyme pyruvate orthophosphate dikinase is encoded by a single gene, Pdk, in the C4 plant Flaveria trinervia. This gene also encodes enzyme isoforms located in the chloroplast and in the cytosol that do not have a function in C4 photosynthesis. Our goal is to identify cis-acting DNA sequences that regulate the expression of the gene that is active in the C4 cycle. We fused 1.5 kb of a 5′ flanking region from the Pdk gene, including the entire 5′ untranslated region, to the uidA reporter gene and stably transformed the closely related C4 species Flaveria bidentis. β-Glucuronidase (GUS) activity was detected at high levels in leaf mesophyll cells. GUS activity was detected at lower levels in bundle-sheath cells and stems and at very low levels in roots. This lower-level GUS expression was similar to the distribution of mRNA encoding the nonphotosynthetic form of the enzyme. We conclude that cis-acting DNA sequences controlling the expression of the C4 form in mesophyll cells and the chloroplast form in other cells and organs are co-located within the same 5′ region of the Pdk gene.

Identification of the Gene Encoding the Tryptophan Synthase β-Subunit from Chlamydomonas reinhardtii1

We report the isolation of a Chlamydomonas reinhardtii cDNA that encodes the β-subunit of tryptophan synthase (TSB). This cDNA was cloned by functional complementation of a trp-operon-deleted strain of Escherichia coli. Hybridization analysis indicated that the gene exists in a single copy. The predicted amino acid sequence showed the greatest identity to TSB polypeptides from other photosynthetic organisms. With the goal of identifying mutations in the gene encoding this enzyme, we isolated 11 recessive and 1 dominant single-gene mutation that conferred resistance to 5-fluoroindole. These mutations fell into three complementation groups, MAA2, MAA7, and TAR1. In vitro assays showed that mutations at each of these loci affected TSB activity. Restriction fragment-length polymorphism analysis suggested that MAA7 encodes TSB. MAA2 and TAR1 may act to regulate the activity of MAA7 or its protein product.

Intrinsic noise in gene regulatory networks

Cells are intrinsically noisy biochemical reactors: low reactant numbers can lead to significant statistical fluctuations in molecule numbers and reaction rates. Here we use an analytic model to investigate the emergent noise properties of genetic systems. We find for a single gene that noise is essentially determined at the translational level, and that the mean and variance of protein concentration can be independently controlled. The noise strength immediately following single gene induction is almost twice the final steady-state value. We find that fluctuations in the concentrations of a regulatory protein can propagate through a genetic cascade; translational noise control could explain the inefficient translation rates observed for genes encoding such regulatory proteins. For an autoregulatory protein, we demonstrate that negative feedback efficiently decreases system noise. The model can be used to predict the noise characteristics of networks of arbitrary connectivity. The general procedure is further illustrated for an autocatalytic protein and a bistable genetic switch. The analysis of intrinsic noise reveals biological roles of gene network structures and can lead to a deeper understanding of their evolutionary origin.

A single mutation that increases maize seed weight.

The maize endosperm-specific gene shrunken2 (Sh2) encodes the large subunit of the heterotetrameric starch synthetic enzyme adenosine diphosphoglucose pyrophosphorylase (AGP; EC 2.7.7.27). Here we exploit an in vivo, site-specific mutagenesis system to create short insertion mutations in a region of the gene known to be involved in the allosteric regulation of AGP. The site-specific mutagen is the transposable element dissociation (Ds). Approximately one-third (8 of 23) of the germinal revertants sequenced restored the wild-type sequence, whereas the remaining revertants contained insertions of 3 or 6 bp. All revertants retained the original reading frame 3' to the insertion site and involved the addition of tyrosine and/or serine. Each insertion revertant reduced total AGP activity and the amount of the SH2 protein. The revertant containing additional tyrosine and serine residues increased seed weight 11-18% without increasing or decreasing the percentage of starch. Other insertion revertants lacking an additional serine reduced seed weight. Reduced sensitivity to phosphate, a long-known inhibitor of AGP, was found in the high seed-weight revertant. This alteration is likely universally important since insertion of tyrosine and serine in the potato large subunit of AGP at the comparable position and expression in Escherichia coli also led to a phosphate-insensitive enzyme. These results show that single gene mutations giving rise to increased seed weight, and therefore perhaps yield, are clearly possible in a plant with a long history of intensive and successful breeding efforts.

Targeted replacement of the mycocerosic acid synthase gene in Mycobacterium bovis BCG produces a mutant that lacks mycosides.

A single gene (mas) encodes the multifunctional enzyme that catalyzes the synthesis of very long chain multiple methyl branched fatty acids called mycocerosic acids that are present only in slow-growing pathogenic mycobacteria and are thought to be important for pathogenesis. To achieve a targeted disruption of mas, an internal 2-kb segment of this gene was replaced with approximately the same size hygromycin-resistance gene (hyg), such that hyg was flanked by 4.7- and 1.4-kb segments of mas. Transformation of Mycobacterium bovis BCG with this construct in a plasmid that cannot replicate in mycobacteria yielded hygromycin-resistant transformants. Screening of 38 such transformants by PCR revealed several transformants representing homologous recombination with single crossover and one with double crossover. With primers representing the hyg termini and those representing the mycobacterial genome segments outside that used to make the transformation construct, the double-crossover mutant yielded PCR products expected from either side of hyg. Gene replacement was further confirmed by the absence of the vector and the 2-kb segment of mas replaced by hyg from the genome of the mutant. Thin-layer and radio-gas chromatographic analyses of the lipids derived from [1-14C]propionate showed that the mutant was incapable of synthesizing mycocerosic acids and mycosides. Thus, homologous recombination with double crossover was achieved in a slow-growing mycobacterium with an intron-containing RecA. The resulting mas-disrupted mutant should allow testing of the postulated roles of mycosides in pathogenesis.

Identification and characterization of a RING zinc finger gene (C-RZF) expressed in chicken embryo cells.

To identify changes in gene expression that occur in chicken embryo brain (CEB) cells as a consequence of their binding to the extracellular matrix molecule cytotactin/tenascin (CT/TN), a subtractive hybridization cloning strategy was employed. One of the cDNA clones identified was predicted to encode 381 amino acids and although it did not resemble any known sequences in the nucleic acid or protein data bases, it did contain the sequence motif for the cysteine-rich C3HC4 type of zinc finger, also known as a RING-finger. This sequence was therefore designated the chicken-RING zinc finger (C-RZF). In addition to the RING-finger, the C-RZF sequence also contained motifs for a leucine zipper, a nuclear localization signal, and a stretch of acidic amino acids similar to the activation domains of some transcription factors. Southern analysis suggested that C-RZF is encoded by a single gene. Northern and in situ hybridization analyses of E8 chicken embryo tissues indicated that expression of the C-RZF gene was restricted primarily to brain and heart. Western analysis of the nuclear and cytoplasmic fractions of chicken embryo heart cells and immunofluorescent staining of chicken embryo cardiocytes with anti-C-RZF antibodies demonstrated that the C-RZF protein was present in the nucleus. The data suggest that we have identified another member of the RING-finger family of proteins whose expression in CEB cells may be affected by CT/TN and whose nuclear localization and sequence motifs predict a DNA-binding function in the nucleus.