In silico identification of conserved intercoding sequences in Leishmania genomes: Unraveling putative cis-regulatory elements

In silico analyses of Leishmania spp. genome data are a powerful resource to improve the understanding of these pathogens' biology. Trypanosomatids such as Leishmania spp. have their protein-coding genes grouped in long polycistronic units of functionally unrelated genes. The control of gene expression happens by a variety of posttranscriptional mechanisms. The high degree of synteny among Leishmania species is accompanied by highly conserved coding sequences (CDS) and poorly conserved intercoding untranslated sequences. To identify the elements involved in the control of gene expression, we conducted an in silico investigation to find conserved intercoding sequences (CICS) in the genomes of L major, L infantum, and L braziliensis. We used a combination of computational tools, such as Linux-Shell, PERL and R languages, BLAST, MSPcrunch, SSAKE, and Pred-A-Term algorithms to construct a pipeline which was able to: (i) search for conservation in target-regions, (ii) eliminate CICS redundancy and mask repeat elements, (iii) predict the mRNA's extremities, (iv) analyze the distribution of orthologous genes within the generated LeishCICS-clusters, (v) assign GO terms to the LeishCICS-clusters. and (vi) provide statistical support for the gene-enrichment annotation. We associated the LeishCICS-cluster data, generated at the end of the pipeline, with the expression profile oft. donovani genes during promastigote-amastigote differentiation, as previously evaluated by others (GEO accession: GSE21936). A Pearson's correlation coefficient greater than 0.5 was observed for 730 LeishCICS-clusters containing from 2 to 17 genes. The designed computational pipeline is a useful tool and its application identified potential regulatory cis elements and putative regulons in Leishmania. (C) 2012 Elsevier B.V. All rights reserved.

Trichoderma harzianum expressed sequence tags for identification of genes with putative roles in mycoparasitism against Fusarium solani

The plant pathogen Fusarium solani causes a disease root rot of common bean (Phaseolus vulgaris) resulting in great losses of yield in irrigated areas of the Southeast and Midwest regions of Brazil. Species of the genus Trichoderma have been used in the biological control of this pathogen as an alternative to chemical control. To gain new insights into the biocontrol mechanism used by Trichoderma harzianum against the phytopathogenic fungus, Fusarium solani, we performed a transcriptome analysis using expressed sequence tags (ESTs) and quantitative real-time PCR (RT-qPCR) approaches. A cDNA library from T. harzianum mycelium (isolate ALL42) grown on cell walls of F. solani (CWFS) was constructed and analyzed. A total of 2927 high quality sequences were selected from 3845 and 37.7% were identified as unique genes. The Gene Ontology analysis revealed that the majority of the annotated genes are involved in metabolic processes (80.9%), followed by cellular process (73.7%). We tested twenty genes that encode proteins with potential role in biological control. RT-qPCR analysis showed that none of these genes were expressed when T. harzianum was challenged with itself. These genes showed different patterns of expression during in vitro interaction between T. harzianum and F. solani. (C) 2012 Elsevier Inc. All rights reserved.

Towards the identification of a tumor suppressor gene on chromosome 3p12: Physical mapping and candidate gene screening

Renal cell carcinoma (RCC) is the most common malignant tumor of the kidney. Characterization of RCC tumors indicates that the most frequent genetic event associated with the initiation of tumor formation involves a loss of heterozygosity or cytogenetic aberration on the short arm of human chromosome 3. A tumor suppressor locus Nonpapillary Renal Carcinoma-1 (NRC-1, OMIM ID 604442) has been previously mapped to a 5–7 cM region on chromosome 3p12 and shown to induce rapid tumor cell death in vivo, as demonstrated by functional complementation experiments. ^ To identify the gene that accounts for the tumor suppressor activities of NRC-1, fine-scale physical mapping was conducted with a novel real-time quantitative PCR based method developed in this study. As a result, NRC-1 was mapped within a 4.6-Mb region defined by two unique sequences within UniGene clusters Hs.41407 and Hs.371835 (78,545Kb–83,172Kb in the NCBI build 31 physical map). The involvement of a putative tumor suppressor gene Robo1/Dutt1 was excluded as a candidate for NRC-1. Furthermore, a transcript map containing eleven candidate genes was established for the 4.6-Mb region. Analyses of gene expression patterns with real-time quantitative RT-PCR assays showed that one of the eleven candidate genes in the interval (TSGc28) is down-regulated in 15 out of 20 tumor samples compared with matched normal samples. Three exons of this gene have been identified by RACE experiments, although additional exon(s) seem to exist. Further gene characterization and functional studies are required to confirm the gene as a true tumor suppressor gene. ^ To study the cellular functions of NRC-1, gene expression profiles of three tumor suppressive microcell hybrids, each containing a functional copy of NRC-1, were compared with those of the corresponding parental tumor cell lines using 16K oligonucleotide microarrays. Differentially expressed genes were identified. Analyses based on the Gene Ontology showed that introduction of NRC-1 into tumor cell lines activates genes in multiple cellular pathways, including cell cycle, signal transduction, cytokines and stress response. NRC-1 is likely to induce cell growth arrest indirectly through WEE1. ^

Identification of putative 1p36 tumor suppressor genes involved in oligodendroglioma development

Oligodendrogliomas are primary neoplasms of the central nervous system (CNS). One of the most common and characteristic chromosomal abnormalities observed in oligodendroglioma is allelic loss of 1p (Reifenberger et al., 1994; Bello et al., 1995). Since 1p loss has been reported for both well-differentiated and anaplastic oligodendroglioma, it is believed to occur early in tumor development (Bello et al., 1995). This allelic loss also has clinical significance, for oligodendroglioma patients with 1p loss generally respond significantly better to combination chemotherapy and have longer average survival than do oligodendroglioma patients without 1p loss (Cairncross et al., 1998). To date, no genes on 1p have been implicated as essential to the development or treatment response of oligodendroglioma. In order to localize and/or identify a gene involved in oligodendroglioma development, I tested 170 oligodendrogliomas for deletions of 1p and tested 26 tumors for differential expression of genes in the region of 1p36. Evidence obtained from these methods implicated two genes, SHREW1 and the gene encoding DNA fragmentation factor beta (DFFB). The function for the SHREW1 locus is currently not well known, but preliminary data suggests that it a novel member of adherens junctions. The DFFB gene is an enhancer for apoptosis. Thus, both SHREW1 and DFFB may be candidates for an oligodendroglioma tumor suppressor. Mutational analysis of both genes did not uncover any mutations. Future studies will evaluate other mechanisms that may be responsible for inactivation of these genes in oligodendrogliomas. ^

Identification of a novel vertebrate circadian clock-regulated gene encoding the protein nocturnin

Photoreceptors of the Xenopus laevis retina are the site of a circadian clock. As part of a differential display screen for rhythmic gene products in this system, we have identified a photoreceptor-specific mRNA expressed in peak abundance at night. cDNA cloning revealed an open reading frame encoding a putative 388 amino acid protein that we have named “nocturnin” (for night-factor). This protein has strong sequence similarity to the C-terminal domain of the yeast transcription factor, CCR4, as well as a leucine zipper-like dimerization motif. Nocturnin mRNA levels exhibit a high amplitude circadian rhythm and nuclear run-on analysis indicates that it is controlled by the retinal circadian clock at the level of transcription. Our observations suggest that nocturnin may function through protein–protein interaction either as a component of the circadian clock or as a downstream effector of clock function.

Aldehyde dehydrogenase class 3 expression: Identification of a cornea-preferred gene promoter in transgenic mice

Aldehyde dehydrogenase class 3 (ALDH3) constitutes 20–40% of the total water-soluble proteins in the mammalian cornea. Here, we show by Northern blot analysis that ALDH3 expression in the mouse is at least 500-fold higher in the cornea than in any other tissue examined, with very low levels of expression detected in the stomach, urinary bladder, ocular lens, and lung. Histochemical localization reveals that this exceptional level of expression in the mouse cornea occurs in the anterior epithelial cells and that little ALDH3 is present in the keratocytes or corneal endothelial cells. A 13-kbp mouse ALDH3 promoter fragment containing >12 kbp of the 5′ flanking sequence, the 40-bp untranslated first exon, and 29 bp of intron 1 directed cat reporter gene expression to tissues that express the endogenous ALDH3 gene, except that transgene promoter activity was higher in the stomach and bladder than in the cornea. By contrast, when driven by a 4.4-kbp mouse ALDH3 promoter fragment [1,050-bp 5′ flanking region, exon 1, intron 1 (3.4 kbp), and 7 bp of exon 2] expression of the cat reporter gene was confined to the corneal epithelial cells, except for very low levels in the liver, effectively reproducing the corneal expression pattern of the endogenous ALDH3 gene. These results indicate that tissue-specific expression of ALDH3 is determined by positive and negative elements in the 5′ flanking region of the gene and suggests putative silencers located in intron 1. We demonstrate regulatory sequences capable of directing cornea-specific gene expression, affording the opportunity for genetic engineering in this transparent tissue.

Phylogeny of genes for secretion NTPases: Identification of the widespread tadA subfamily and development of a diagnostic key for gene classification

Macromolecular transport systems in bacteria currently are classified by function and sequence comparisons into five basic types. In this classification system, type II and type IV secretion systems both possess members of a superfamily of genes for putative NTP hydrolase (NTPase) proteins that are strikingly similar in structure, function, and sequence. These include VirB11, TrbB, TraG, GspE, PilB, PilT, and ComG1. The predicted protein product of tadA, a recently discovered gene required for tenacious adherence of Actinobacillus actinomycetemcomitans, also has significant sequence similarity to members of this superfamily and to several unclassified and uncharacterized gene products of both Archaea and Bacteria. To understand the relationship of tadA and tadA-like genes to those encoding the putative NTPases of type II/IV secretion, we used a phylogenetic approach to obtain a genealogy of 148 NTPase genes and reconstruct a scenario of gene superfamily evolution. In this phylogeny, clear distinctions can be made between type II and type IV families and their constituent subfamilies. In addition, the subgroup containing tadA constitutes a novel and extremely widespread subfamily of the family encompassing all putative NTPases of type IV secretion systems. We report diagnostic amino acid residue positions for each major monophyletic family and subfamily in the phylogenetic tree, and we propose an easy method for precisely classifying and naming putative NTPase genes based on phylogeny. This molecular key-based method can be applied to other gene superfamilies and represents a valuable tool for genome analysis.

Identification of a Saccharomyces cerevisiae Gene that Is Required for G1 Arrest in Response to the Lipid Oxidation Product Linoleic Acid Hydroperoxide*

Reactive oxygen species cause damage to all of the major cellular constituents, including peroxidation of lipids. Previous studies have revealed that oxidative stress, including exposure to oxidation products, affects the progression of cells through the cell division cycle. This study examined the effect of linoleic acid hydroperoxide, a lipid peroxidation product, on the yeast cell cycle. Treatment with this peroxide led to accumulation of unbudded cells in asynchronous populations, together with a budding and replication delay in synchronous ones. This observed modulation of G1 progression could be distinguished from the lethal effects of the treatment and may have been due to a checkpoint mechanism, analogous to that known to be involved in effecting cell cycle arrest in response to DNA damage. By examining several mutants sensitive to linoleic acid hydroperoxide, the YNL099c open reading frame was found to be required for the arrest. This gene (designated OCA1) encodes a putative protein tyrosine phosphatase of previously unknown function. Cells lacking OCA1 did not accumulate in G1 on treatment with linoleic acid hydroperoxide, nor did they show a budding, replication, or Start delay in synchronous cultures. Although not essential for adaptation or immediate cellular survival, OCA1 was required for growth in the presence of linoleic acid hydroperoxide, thus indicating that it may function in linking growth, stress responses, and the cell cycle. Identification of OCA1 establishes cell cycle arrest as an actively regulated response to oxidative stress and will enable further elucidation of oxidative stress-responsive signaling pathways in yeast.

Identification of the promoter of the mouse obese gene.

Primer extension and RACE (rapid amplification of cDNA ends) assays were used to identify and sequence the 5' terminus of mouse ob mRNA. This sequence was used to obtain a recombinant bacteriophage containing the first exon of the encoding gene. DNA sequence analysis of the region immediately upstream of the first exon of the mouse ob gene revealed DNA sequences corresponding to presumptive cis-regulatory elements. A canonical TATA box was observed 30-34 base pairs upstream from the start site of transcription and a putative binding site for members of the C/EBP family of transcription factors was identified immediately upstream from the TATA box. Nuclear extracts prepared from primary adipocytes contained a DNA binding activity capable of avid and specific interaction with the putative C/EBP response element; antibodies to C/EBP alpha neutralized the DNA binding activity present in adipocyte nuclear extracts. When linked to a firefly luciferase reporter and transfected into primary adipocytes, the presumptive promoter of the mouse ob gene facilitated luciferase expression. When transfected into HepG2 cells, which lack C/EBP alpha, the mouse ob promoter was only weakly active. Supplementation of C/EBP alpha by cotransfection with a C/EBP alpha expression vector markedly stimulated luciferase expression. Finally, an ob promoter variant mutated at the C/EBP response element was inactive in both primary adipocytes and HepG2 cells. These observations provide evidence for identification of a functional promoter capable of directing expression of the mouse ob gene.

Identification of putative noncoding RNAs among the RIKEN mouse full-length cDNA collection

With the sequencing and annotation of genomes and transcriptomes of several eukaryotes, the importance of noncoding RNA (ncRNA)-RNA molecules that are not translated to protein products-has become more evident. A subclass of ncRNA transcripts are encoded by highly regulated, multi-exon, transcriptional units, are processed like typical protein-coding mRNAs and are increasingly implicated in regulation of many cellular functions in eukaryotes. This study describes the identification of candidate functional ncRNAs from among the RIKEN mouse full-length cDNA collection, which contains 60,770 sequences, by using a systematic computational filtering approach. We initially searched for previously reported ncRNAs and found nine murine ncRNAs and homologs of several previously described nonmouse ncRNAs. Through our computational approach to filter artifact-free clones that lack protein coding potential, we extracted 4280 transcripts as the largest-candidate set. Many clones in the set had EST hits, potential CpG islands surrounding the transcription start sites, and homologies with the human genome. This implies that many candidates are indeed transcribed in a regulated manner. Our results demonstrate that ncRNAs are a major functional subclass of processed transcripts in mammals.

Identification of independent association signals and putative functional variants for breast cancer risk through fine-scale mapping of the 12p11 locus

BACKGROUND: Multiple recent genome-wide association studies (GWAS) have identified a single nucleotide polymorphism (SNP), rs10771399, at 12p11 that is associated with breast cancer risk. METHOD: We performed a fine-scale mapping study of a 700 kb region including 441 genotyped and more than 1300 imputed genetic variants in 48,155 cases and 43,612 controls of European descent, 6269 cases and 6624 controls of East Asian descent and 1116 cases and 932 controls of African descent in the Breast Cancer Association Consortium (BCAC; http://bcac.ccge.medschl.cam.ac.uk/ ), and in 15,252 BRCA1 mutation carriers in the Consortium of Investigators of Modifiers of BRCA1/2 (CIMBA). Stepwise regression analyses were performed to identify independent association signals. Data from the Encyclopedia of DNA Elements project (ENCODE) and the Cancer Genome Atlas (TCGA) were used for functional annotation. RESULTS: Analysis of data from European descendants found evidence for four independent association signals at 12p11, represented by rs7297051 (odds ratio (OR) = 1.09, 95 % confidence interval (CI) = 1.06-1.12; P = 3 × 10(-9)), rs805510 (OR = 1.08, 95 % CI = 1.04-1.12, P = 2 × 10(-5)), and rs1871152 (OR = 1.04, 95 % CI = 1.02-1.06; P = 2 × 10(-4)) identified in the general populations, and rs113824616 (P = 7 × 10(-5)) identified in the meta-analysis of BCAC ER-negative cases and BRCA1 mutation carriers. SNPs rs7297051, rs805510 and rs113824616 were also associated with breast cancer risk at P < 0.05 in East Asians, but none of the associations were statistically significant in African descendants. Multiple candidate functional variants are located in putative enhancer sequences. Chromatin interaction data suggested that PTHLH was the likely target gene of these enhancers. Of the six variants with the strongest evidence of potential functionality, rs11049453 was statistically significantly associated with the expression of PTHLH and its nearby gene CCDC91 at P < 0.05. CONCLUSION: This study identified four independent association signals at 12p11 and revealed potentially functional variants, providing additional insights into the underlying biological mechanism(s) for the association observed between variants at 12p11 and breast cancer risk

Identification of a novel prostate cancer susceptibility variant in the KLK3 gene transcript

Abstract Genome-wide association studies (GWAS) have identified more than 30 prostate cancer (PrCa) susceptibility loci. One of these (rs2735839) is located close to a plausible candidate susceptibility gene, KLK3, which encodes prostate-specific antigen (PSA). PSA is widely used as a biomarker for PrCa detection and disease monitoring. To refine the association between PrCa and variants in this region, we used genotyping data from a two-stage GWAS using samples from the UK and Australia, and the Cancer Genetic Markers of Susceptibility (CGEMS) study. Genotypes were imputed for 197 and 312 single nucleotide polymorphisms (SNPs) from HapMap2 and the 1000 Genome Project, respectively. The most significant association with PrCa was with a previously unidentified SNP, rs17632542 (combined P = 3.9 × 10−22). This association was confirmed by direct genotyping in three stages of the UK/Australian GWAS, involving 10,405 cases and 10,681 controls (combined P = 1.9 × 10−34). rs17632542 is also shown to be associated with PSA levels and it is a non-synonymous coding SNP (Ile179Thr) in KLK3. Using molecular dynamic simulation, we showed evidence that this variant has the potential to introduce alterations in the protein or affect RNA splicing. We propose that rs17632542 may directly influence PrCa risk.

Sequence and identification of the barley yellow dwarf virus coat protein gene

The nucleotide sequence of the coat protein gene of barley yellow dwarf virus (BYDV, PAV serotype) was determined, and the amino acid sequence was deduced. The open reading frame, encoding a protein of relative molecular mass (Mr) 22,047, was confirmed as the coat protein gene by comparison with amino acid sequences of tryptic peptides derived from dissociated virions. In addition, a fragment of this gene expressed in Escherichia coli produced a product which was recognized by antibodies prepared against purified BYDV virions. An overlapping reading frame encoding an Mr 17,147 protein is contained completely within the coat protein gene. © 1988.