Genome-wide analysis and experimentation of plant serine/threonine/tyrosine-specific protein kinases

Protein tyrosine phosphorylation plays an important role in cell growth, development and oncogenesis. No classical protein tyrosine kinase has hitherto been cloned from plants. Does protein tyrosine kinase exist in plants? To address this, we have performed a genomic survey of protein tyrosine kinase motifs in plants using the delineated tyrosine phosphorylation motifs from the animal system. The Arabidopsis thaliana genome encodes 57 different protein kinases that have tyrosine kinase motifs. Animal non-receptor tyrosine kinases, SRC, ABL, LYN, FES, SEK, KIN and RAS have structural relationship with putative plant tyrosine kinases. In an extended analysis, animal receptor and non-receptor kinases, Raf and Ras kinases, mixed lineage kinases and plant serine/threonine/tyrosine (STY) protein kinases, form a well-supported group sharing a common origin within the superfamily of STY kinases. We report that plants lack bona fide tyrosine kinases, which raise an intriguing possibility that tyrosine phosphorylation is carried out by dual-specificity STY protein kinases in plants. The distribution pattern of STY protein kinase families on Arabidopsis chromosomes indicates that this gene family is partly a consequence of duplication and reshuffling of the Arabidopsis genome and of the generation of tandem repeats. Genome-wide analysis is supported by the functional expression and characterization of At2g24360 and phosphoproteomics of Arabidopsis. Evidence for tyrosine phosphorylated proteins is provided by alkaline hydrolysis, anti-phosphotyrosine immunoblotting, phosphoamino acid analysis and peptide mass fingerprinting. These results report the first comprehensive survey of genome-wide and tyrosine phosphoproteome analysis of plant STY protein kinases.

A Unique Modification of the Eukaryotic Initiation Factor 5A Shows the Presence of the Complete Hypusine Pathway in Leishmania donovani

Deoxyhypusine hydroxylase (DOHH) catalyzes the final step in the post-translational synthesis of an unusual amino acid hypusine (N-(sic)-(4-amino-2-hydroxybutyl) lysine), which is present on only one cellular protein, eukaryotic initiation factor 5A (eIF5A). We present here the molecular and structural basis of the function of DOHH from the protozoan parasite, Leishmania donovani, which causes visceral leishmaniasis. The L. donovani DOHH gene is 981 bp and encodes a putative polypeptide of 326 amino acids. DOHH is a HEAT-repeat protein with eight tandem repeats of alpha-helical pairs. Four conserved histidine-glutamate sequences have been identified that may act as metal coordination sites. A similar to 42 kDa recombinant protein with a His-tag was obtained by heterologous expression of DOHH in Escherichia coli. Purified recombinant DOHH effectively catalyzed the hydroxylation of the intermediate, eIF5A-deoxyhypusine (eIF5A-Dhp), in vitro. L. donovani DOHH (LdDOHH) showed similar to 40.6% sequence identity with its human homolog. The alignment of L. donovani DOHH with the human homolog shows that there are two significant insertions in the former, corresponding to the alignment positions 159-162 (four amino acid residues) and 174-183 (ten amino acid residues) which are present in the variable loop connecting the N- and C-terminal halves of the protein, the latter being present near the substrate binding site. Deletion of the ten-amino-acid-long insertion decreased LdDOHH activity to 14% of the wild type recombinant LdDOHH. Metal chelators like ciclopirox olamine (CPX) and mimosine significantly inhibited the growth of L. donovani and DOHH activity in vitro. These inhibitors were more effective against the parasite enzyme than the human enzyme. This report, for the first time, confirms the presence of a complete hypusine pathway in a kinetoplastid unlike eubacteria and archaea. The structural differences between the L. donovani DOHH and the human homolog may be exploited for structure based design of selective inhibitors against the parasite.

An algorithm to find all palindromic sequences in proteins

A palindrome is a set of characters that reads the same forwards and backwards. Since the discovery of palindromic peptide sequences two decades ago, little effort has been made to understand its structural, functional and evolutionary significance. Therefore, in view of this, an algorithm has been developed to identify all perfect palindromes (excluding the palindromic subset and tandem repeats) in a single protein sequence. The proposed algorithm does not impose any restriction on the number of residues to be given in the input sequence. This avant-garde algorithm will aid in the identification of palindromic peptide sequences of varying lengths in a single protein sequence.

Gene enrichment using antibodies to DNA/RNA hybrids : mapping the ribosomal DNA of slime mold and rat

A novel method for gene enrichment has been developed and applied to mapping the rRNA genes of two eucaryotic organisms. The method makes use of antibodies to DNA/RNA hybrids prepared by injecting rabbits with the synthetic hybrid poly(rA)•poly(dT). Antibodies which cross-react with non-hybrid nucleic acids were removed from the purified IgG fraction by adsorption on columns of DNA-Sepharose, oligo(dT)-cellulose, and poly(rA)-Sepharose. Subsequent purification of the specific DNA/RNA hybrid antibody was carried out on a column of oligo(dT)-cellulose to which poly(rA) was hybridized. Attachment of these antibodies to CNBr-activated Sepharose produced an affinity resin which specifically binds DNA/RNA hybrids.

In order to map the rDNA of the slime mold Dictyostelium discoideum, R-loops were formed using unsheared nuclear DNA and the 178 and 268 rRNAs of this organism. This mixture was passed through a column containing the affinity resin, and bound molecules containing R- loops were eluted by high salt. This purified rDN A was observed directly in the electron microscope. Evidence was obtained that there is a physical end to Dictyostelium rDN A molecules approximately 10 kilobase pairs (kbp) from the region which codes for the 268 rRNA. This finding is consistent with reports of other investigators that the rRNA genes exist as inverse repeats on extra-chromosomal molecules of DNA unattached to the remainder of the nuclear DNA in this organism.

The same general procedure was used to map the rRNA genes of the rat. Molecules of DNA which contained R-loops formed with the 188 and 288 rRNAs were enriched approximately 150- fold from total genomal rat DNA by two cycles of purification on the affinity column. Electron microscopic measurements of these molecules enabled the construction of an R-loop map of rat rDNA. Eleven of the observed molecules contained three or four R-loops or else two R-loops separated by a long spacer. These observations indicated that the rat rRNA genes are arranged as tandem repeats. The mean length of the repeating units was 37.2 kbp with a standard deviation of 1.3 kbp. These eleven molecules may represent repeating units of exactly the same length within the errors of the measurements, although a certain degree of length heterogeneity cannot be ruled out. If significantly shorter or longer repeating units exist, they are probably much less common than the 37.2 kbp unit.

The last section of the thesis describes the production of antibodies to non-histone chromosomal proteins which have been exposed to the ionic detergent sodium dodecyl sulfate (SDS). The presence of low concentrations of SDS did not seem to affect either production of antibodies or their general specificity. Also, a technique is described for the in situ immunofluorescent detection of protein antigens in polyacrylamide gels.

Associação entre polimorfismos dos genes IL1B, IL1RN e TNFA e a periodontite agressiva e a periodontite crônica severa

A periodontite é um processo inflamatório crônico de origem bacteriana mediado por citocinas, em especial, interleucina-1 (IL1) e fator de necrose tumoral (TNFα). Polimorfismos genéticos de IL1 e TNFA têm sido associados com a variação de expressão dessas proteínas, o que poderia justificar as diferenças interindividuais de manifestação da doença. O objetivo do presente estudo foi investigar possíveis associações entre os genes IL1B, IL1RN e TNFA e a suscetibilidade à periodontite agressiva e à periodontite crônica severa. Foram selecionados 145 pacientes do Estado do Rio de Janeiro, 43 com periodontite agressiva (PAgr) (33,1 4,8 anos), 52 com periodontite crônica severa (PCr) (50,6 5,8 anos) e 50 controles (40,1 7,8 anos). Os DNAs genômicos dos integrantes dos grupos PAgr, PCr e controle foram obtidos através da coleta de células epiteliais bucais raspadas da parte interna da bochecha com cotonete. Os SNPs IL1B -511C>T, IL1B +3954C>T e TNFA -1031T>C foram analisados pela técnica de PCR-RFLP, utilizando as enzimas de restrição Ava I Taq I e Bpi I, respectivamente. O polimorfismo de número variável de repetições in tandem (VNTR) no intron 2 do gene IL1RN foi feita pela análise direta dos amplicons. Todos os polimorfismos foram analisados por eletroforese em gel de poliacrilamida 8%. As frequências alélica e genotípica do polimorfismo IL1B +3954C>T no grupo PCr foram significativamente diferentes das observadas no grupo controle (p=0,003 e p=0,041, respectivamente). A freqüência do alelo A2 do polimorfismo IL1RN VNTR intron2 no grupo PAgr foi significativamente maior do que no grupo controle (p=0,035). Não houve associação entre os polimorfismos IL1B -511C>T e TNFA -1031T>C e as periodontites agressiva e crônica. A presença dos alelos 2 nos genótipos combinados de IL1RN VNTR intron2 e IL1B +3954C>T no grupo PCr foi significativamente maior quando comparada ao grupo controle (p=0,045). Entretanto, não se observou associação entre as combinações genotípicas IL1B -511C>T / IL1B +3954C>T e IL1RN VNTR / IL1B -511C>T e a predisposição à doença periodontal. De acordo com os nossos resultados podemos sugerir que, para a população estudada, o polimorfismo IL1B +3954C>T interfere no desenvolvimento da periodontite crônica, enquanto a presença do alelo A2 do polimorfismo IL1RN VNTR intron2 pode ser considerado como indicador de risco para a periodontite agressiva. O presente estudo também nos permite sugerir que a ausência de homozigose dos alelos 1 nos genótipos combinados de IL1RN VNTR intron2 e IL1B +3954C>T pode representar maior suscetibilidade à periodontite crônica severa.

Cloning of bradykinin precursor cDNAs from skin of Bombina maxima reveals novel bombinakinin M antagonists and a bradykinin potential peptide

Bombinakinin M (DLPKINRKGP-bradykinin) is a bradykinin-related peptide purified from skin secretions of the frog Bombina maxima. As previously reported, its biosynthesis is characterized by a tandem repeats with various copy numbers of the peptide and sometimes co-expressed with other structure-function distinguishable peptides. At present study, two novel cDNAs encoding bombinakinin M and its variants were cloned from a cDNA library from the skin of the frog. The encoded two precursor proteins are common in that each contains three repeats of a novel 16-amino acid peptide unit and one copy of kinestatin at their N- and C-terminal parts, respectively. They differ in that the first precursor contains two copies of bombinakinin M and the second one contains one copy of a novel bombinakinin M variant. Bombinakinin M was found to elicit concentration-dependent contractile effects on guinea pig ileum, with an EC50 value of 4 nM that is four times higher than that of bradykinin (1 nM). Interestingly, the synthetic peptide (DYTIRTRLH-amide), as deduced from the 16-amino acid peptide repeats in the newly cloned cDNAs, possessed weak inhibitory activity on the contractile effects of bombinakinin M, but not on that of bradykinin. Furthermore, the newly identified bombinakinin M variant (DLSKMSFLHG-Ile(1)-bradykinin), did not show contractile activity on guinea pig ileum, but showed potentiation effect on the myotropic activity of bradykinin. In a molar raito of 1:58, it augmented the activity of bradykinin up to two-fold. (C) 2004 Elsevier B.V. All rights reserved.

MicroRNAs and other small RNAs enriched in the Arabidopsis RNA-dependent RNA polymerase-2 mutant.

The Arabidopsis genome contains a highly complex and abundant population of small RNAs, and many of the endogenous siRNAs are dependent on RNA-Dependent RNA Polymerase 2 (RDR2) for their biogenesis. By analyzing an rdr2 loss-of-function mutant using two different parallel sequencing technologies, MPSS and 454, we characterized the complement of miRNAs expressed in Arabidopsis inflorescence to considerable depth. Nearly all known miRNAs were enriched in this mutant and we identified 13 new miRNAs, all of which were relatively low abundance and constitute new families. Trans-acting siRNAs (ta-siRNAs) were even more highly enriched. Computational and gel blot analyses suggested that the minimal number of miRNAs in Arabidopsis is approximately 155. The size profile of small RNAs in rdr2 reflected enrichment of 21-nt miRNAs and other classes of siRNAs like ta-siRNAs, and a significant reduction in 24-nt heterochromatic siRNAs. Other classes of small RNAs were found to be RDR2-independent, particularly those derived from long inverted repeats and a subset of tandem repeats. The small RNA populations in other Arabidopsis small RNA biogenesis mutants were also examined; a dcl2/3/4 triple mutant showed a similar pattern to rdr2, whereas dcl1-7 and rdr6 showed reductions in miRNAs and ta-siRNAs consistent with their activities in the biogenesis of these types of small RNAs. Deep sequencing of mutants provides a genetic approach for the dissection and characterization of diverse small RNA populations and the identification of low abundance miRNAs.

Attention deficit/hyperactivity disorder children with a 7-repeat allele of the dopamine receptor D4 gene have extreme behavior but normal performance on critical neuropsychological tests of attention

An association of the dopamine receptor D4 (DRD4) gene located on chromosome 11p15.5 and attention deficit/hyperactivity disorder (ADHD) has been demonstrated and replicated by multiple investigators. A specific allele [the 7-repeat of a 48-bp variable number of tandem repeats (VNTR) in exon 3] has been proposed as an etiological factor in attentional deficits manifested in some children diagnosed with this disorder. In the current study, we evaluated ADHD subgroups defined by the presence or absence of the 7-repeat allele of the DRD4 gene, using neuropsychological tests with reaction time measures designed to probe attentional networks with neuroanatomical foci in D4-rich brain regions. Despite the same severity of symptoms on parent and teacher ratings for the ADHD subgroups, the average reaction times of the 7-present subgroup showed normal speed and variability of response whereas the average reaction times of the 7-absent subgroup showed the expected abnormalities (slow and variable responses). This was opposite the primary prediction of the study. The 7-present subgroup seemed to be free of some of the neuropsychological abnormalities thought to characterize ADHD.

The sequence and de novo assembly of the giant panda genome

Using next-generation sequencing technology alone, we have successfully generated and assembled a draft sequence of the giant panda genome. The assembled contigs (2.25 gigabases (Gb)) cover approximately 94% of the whole genome, and the remaining gaps (0.05 Gb) seem to contain carnivore-specific repeats and tandem repeats. Comparisons with the dog and human showed that the panda genome has a lower divergence rate. The assessment of panda genes potentially underlying some of its unique traits indicated that its bamboo diet might be more dependent on its gut microbiome than its own genetic composition. We also identified more than 2.7 million heterozygous single nucleotide polymorphisms in the diploid genome. Our data and analyses provide a foundation for promoting mammalian genetic research, and demonstrate the feasibility for using next-generation sequencing technologies for accurate, cost-effective and rapid de novo assembly of large eukaryotic genomes.

Genome-wide identification and characterization of cytochrome P450 monooxygenase genes in the ciliate Tetrahymena thermophila

Background: Cytochrome P450 monooxygenases play key roles in the metabolism of a wide variety of substrates and they are closely associated with endocellular physiological processes or detoxification metabolism under environmental exposure. To date, however, none has been systematically characterized in the phylum Ciliophora. T. thermophila possess many advantages as a eukaryotic model organism and it exhibits rapid and sensitive responses to xenobiotics, making it an ideal model system to study the evolutionary and functional diversity of the P450 monooxygenase gene family. Results: A total of 44 putative functional cytochrome P450 genes were identified and could be classified into 13 families and 21 sub-families according to standard nomenclature. The characteristics of both the conserved intron-exon organization and scaffold localization of tandem repeats within each P450 family clade suggested that the enlargement of T. thermophila P450 families probably resulted from recent separate small duplication events. Gene expression patterns of all T. thermophila P450s during three important cell physiological stages (vegetative growth, starvation and conjugation) were analyzed based on EST and microarray data, and three main categories of expression patterns were postulated. Evolutionary analysis including codon usage preference, sit-especific selection and gene-expression evolution patterns were investigated and the results indicated remarkable divergences among the T. thermophila P450 genes. Conclusion: The characterization, expression and evolutionary analysis of T. thermophila P450 monooxygenase genes in the current study provides useful information for understanding the characteristics and diversities of the P450 genes in the Ciliophora, and provides the baseline for functional analyses of individual P450 isoforms in this model ciliate species.

Cloning, characterization, and gene expression analysis of a novel cadmium metallothionein gene in Tetrahymena pigmentosa

A novel cadmium-inducible metallothionein (MT) gene (Tpig-MT1) was cloned and sequenced from the ciliate Tetrahymena pigmentosa. The number of deduced amino acids is 118. The polypeptide possesses CCC and CC clusters characteristic of typical Tetrahymena Cd-inducible MTs. The structure of Tpig-MT1 is different from the reported Cd-MT in T. pyriformis, T. thermophila and T. pigmentosa. Tpig-MT1 contains two intragenic tandem repeats with 72.9% identity described as Tpig-MT1 (repeat A1) and Tpig-MT1 (repeat A2). The transcriptional response of Tpig-MT1 gene to different heavy metals (Cd, Cu, Zn, Hg, Pb) and oxidative stress (H2O2) was measured using real-time quantitative PCR. The results showed that the gene was quickly induced (1 h) by the five heavy metals and the order of expression level was Hg>Pb>Cd>Cu>Zn. The induction effect of H2O2 was 5-fold after about 15 min, but soon decreased to a non-significant level (30 min). The genetic diversity of Tetrahymena MT genes is discussed in relation to the unique structure of the Tpig-MT1 gene and other reported Cd-MT isoforms. (C) 2008 Elsevier B.V. All rights reserved.

Cloning and characterization of interferon stimulated genes Viperin and ISG15, and their promoters from snakehead Channa argus

By suppression subtractive hybridization, rapid amplification of cDNA ends and gene walking methods, interferon stimulated genes (ISGs), Viperin and ISG15, and their promoters have been cloned and characterized from snakehead Channa argus. The Viperin cDNA was found to be 1474 nt and contain an open reading frame (ORF) of 1059 nt that translates into a putative peptide of 352 amino acid (aa). The putative peptide of Viperin shows high identity to that in teleosts and mammals except for the N-terminal 70 aa. The ISG15 cDNA was found to be 758 nt and contain an ORF of 468 nt that translates into a putative peptide of 155 aa. The putative peptide of ISG15 is composed of two tandem repeats of ubiquitin-like (UBL) domains, and a canonical conjugation motif (LRGG) at C-terminal. Viperin and ISG15 promoter regions were characterized by the presence of interferon stimulating response elements (ISRE) and gamma-IFN activation sites (GAS). ISRE is a feature of IFN-induced gene promoter and partially overlaps interferon regulatory factor (IRF) 1 and IRF2 recognition sites. GAS is responsible for the gamma-IFN mediated transcription. One conserved site for NF-kappa B was found in the promoter region of Viperin. This is the first report of conservative binding motif for NF-kappa B in accordance with the consensus sequence (GGGRN-NYYCC) among teleost ISG promoters. Moreover, there were also TATA, CAAT and Sp1 transcription factor sites in Viperin and ISG15 promoters. In 5' untranslated region (UTR), snakehead ISG15 gene contains a single intron, which differs from Viperin gene. The transcripts of Vipeirn and ISG15 mRNA were mainly expressed in head kidney, posterior kidney, spleen and gill. The expression levels in liver were found to increase obviously in response to induction by IFN-inducer poly I : C.

The complete mitochondrial genome sequence of the cutlassfish Trichiurus japonicus (Perciformes: Trichiuridae) : Genome characterization and phylogenetic considerations

Mitochondrial genome sequence and structure analysis has become a powerful tool for studying molecular evolution and phylogenetic relationships. To understand the systematic status of Trichiurus japonicus in suborder Scombroidei, we determined the complete mitochondrial genome (mitogenome) sequence using the long-polymerase chain reaction (long-PCR) and shotgun sequencing method. The entire mitogenome is 16,796 by in length and has three unusual features, including (1) the absence of tRNA(Pro) gene, (2) the possibly nonfunctional light-strand replication origin (O-L) showing a shorter loop in secondary structure and no conserved motif (5'-GCCGG-3'), (3) two sets of the tandem repeats at the 5' and 3' ends of the control region. The three features seem common for Trichiurus mitogenomes, as we have confirmed them in other three T. japonicus individuals and in T nanhaiensis. Phylogenetic analysis does not support the monophyly of Trichiuridae, which is against the morphological result. T. japonicus is most closely related to those species of family Scombridae; they in turn have a sister relationship with Perciformes members including suborders Acanthuroidei, Caproidei, Notothenioidei, Zoarcoidei, Trachinoidei, and some species of Labroidei, based on the current dataset of complete mitogenome. T japonicus together with T. brevis, T lepturus and Aphanopus carbo form a clade distinct from Lepidopus caudatus in terms of the complete Cyt b sequences. T. japonicus mitogenome, as the first discovered complete mitogenome of Trichiuridae, should provide important information on both genomics and phylogenetics of Trichiuridae. (C) 2009 Elsevier B.V. All rights reserved.

A centromeric satellite sequence in the Pacific oyster (Crassostrea gigas Thunberg) identified by fluorescence in situ hybridization

A highly repetitive satellite sequence was previously identified in the Pacific oyster Crassostrea gigas Thunberg. The sequence has 168 bp per unit, present in tandem repeats, and accounts for 1% to 4% of the genome. We studied the chromosomal location of this satellite sequence by fluorescence in situ hybridization (FISH), A probe was made by polymerase chain reaction and incorporation of digoxigenin-11-dUTP. Hybridization was detected with fluorescein-labeled antidigoxigenin antibodies. FISH signals were located at centromeric regions of 7 pairs of the Pacific oyster chromosomes. No interstitial site was found. Signals were strong and consistent on chromosomes 1, 2, 4, and 7, but weak or variable oil chromosomes 5, 8, and 10. No signal was observed on chromosomes 3, 6, and 9. Our results showed that this sequence is clearly a centromeric satellite, disputing its previous assignment to the telomeric and submetacentric regions of 2 chromosomes. No signal was detected in the American oyster (Crassostrea virginica Gmelin).

Predictive Modeling of Loss-of-Heterozygosity Events in Yeast

During mitotic cell cycles, DNA experiences many types of endogenous and exogenous damaging agents that could potentially cause double strand breaks (DSB). In S. cerevisiae, DSBs are primarily repaired by mitotic recombination and as a result, could lead to loss-of-heterozygosity (LOH). Genetic recombination can happen in both meiosis and mitosis. While genome-wide distribution of meiotic recombination events has been intensively studied, mitotic recombination events have not been mapped unbiasedly throughout the genome until recently. Methods for selecting mitotic crossovers and mapping the positions of crossovers have recently been developed in our lab. Our current approach uses a diploid yeast strain that is heterozygous for about 55,000 SNPs, and employs SNP-Microarrays to map LOH events throughout the genome. These methods allow us to examine selected crossovers and unselected mitotic recombination events (crossover, noncrossover and BIR) at about 1 kb resolution across the genome. Using this method, we generated maps of spontaneous and UV-induced LOH events. In this study, we explore machine learning and variable selection techniques to build a predictive model for where the LOH events occur in the genome.

Randomly from the yeast genome, we simulated control tracts resembling the LOH tracts in terms of tract lengths and locations with respect to single-nucleotide-polymorphism positions. We then extracted roughly 1,100 features such as base compositions, histone modifications, presence of tandem repeats etc. and train classifiers to distinguish control tracts and LOH tracts. We found interesting features of good predictive values. We also found that with the current repertoire of features, the prediction is generally better for spontaneous LOH events than UV-induced LOH events.