Molecular features of colorectal cancer predisposition and progression

Both inherited genetic variations and somatically acquired mutations drive cancer development. The aim of this thesis was to gain insight into the molecular mechanisms underlying colorectal cancer (CRC) predisposition and tumor progression. Whereas one-third of CRC may develop in the context of hereditary predisposition, the known highly penetrant syndromes only explain a small fraction of all cases. Genome-wide association studies have shown that ten common single nucleotide polymorphisms (SNPs) modestly predispose to CRC. Our population-based sample series of around thousand CRC cases and healthy controls was genotyped for these SNPs. Tumors of heterozygous patients were analyzed for allelic imbalance, in an attempt to reveal the role of these SNPs in somatic tumor progression. The risk allele of rs6983267 at 8q24 was favored in the tumors significantly more often than the neutral allele, indicating that this germline variant is somatically selected for. No imbalance targeting the risk allele was observed in the remaining loci, suggesting that most of the low-penetrance CRC SNPs mainly play a role in the early stages of the neoplastic process. The ten SNPs were further analyzed in 788 CRC cases, 97 of which had a family history of CRC, to evaluate their combined contribution. A significant association appeared between the overall number of risk alleles and familial CRC and these ten SNPs seem to explain around 9% of the familial clustering of CRC. Finding more CRC susceptibility alleles may facilitate individualized risk prediction and cancer prevention in the future. Microsatellite instability (MSI), resulting from defective mismatch repair function, is a hallmark of Lynch syndrome and observed in a subset of all CRCs. Our aim was to identify microsatellite frameshift mutations that inactivate tumor suppressor genes in MSI CRCs. By sequencing microsatellite repeats of underexpressed genes we found six novel MSI target genes that were frequently mutated in 100 MSI CRCs: 51% in GLYR1, 47% in ABCC5, 43% in WDTC1, 33% in ROCK1, 30% in OR51E2, and 28% in TCEB3. Immunohistochemical staining of GLYR1 revealed defective protein expression in homozygously mutated tumors, providing further support for the loss of function hypothesis. Another mutation screening effort sought to identify MSI target genes with putative oncogenic functions. Microsatellites were similarly sequenced in genes that were overexpressed and, upon mutation, predicted to avoid nonsense-mediated mRNA decay. The mitotic checkpoint kinase TTK harbored protein-elongating mutations in 59% of MSI CRCs and the mutant protein was detected in heterozygous MSI CRC cells. No checkpoint dysregulation or defective protein localization was observable however, and the biological relevance of this mutation may hence be related to other mechanisms. In conclusion, these two large-scale and unbiased efforts identified frequently mutated genes that are likely to contribute to the development of this cancer type and may be utilized in developing diagnostic and therapeutic applications.

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.

The Transcriptional Regulation of Retrotransposon BARE

The purpose of this research project was to understand the steps of the retrotransposon BARE (BArley REtrotransposon) life cycle, from regulation of transcription to Virus-Like Particle (VLP) formation and ultimate integration back into the genome. Our study concentrates mainly on BARE1 transcriptional regulation because transcription is the crucial first step in the retrotransposon life cycle. The BARE element is a Class I LTR (Long Terminal Repeat) retrotransposon belonging to the Copia superfamily and was originally isolated in our research group. The LTR retrotransposons are transcribed from promoters in the LTRs and encode proteins for packaging of their transcripts, the reverse transcription of the transcripts into cDNA, and integration of the cDNA back into the genome. BARE1 is translated as a single polyprotein and cleaved into the capsid protein (GAG), integrase (IN), and reverse transcriptase-RNaseH (RT-RH) by the integral aspartic proteinase (AP). The BARE retrotransposon family comprises more than 104 copies in the barley (Hordeum vulgare) genome. The element is bound by long terminal repeats (LTRs, 1829 bp) containing promoters required for replication, signals for RNA processing, and motifs necessary for the integration of the cDNA. Members of the BARE1 subfamily are transcribed, translated, and form virus-like particles. Several basic questions concerning transcription are explored in the thesis: BARE1 transcription control, promoter choice in different barley tissues, start and termination sites for BARE transcripts, and BARE1 transcript polyadenylation (I). Polyadenylation is an important step during mRNA maturation, and determines its stability and translatability among other characteristics. Our work has found a novel way used by BARE1 to make extra GAG protein, which is critical for VLP formation. The discovery that BARE1 uses one RNA population for protein synthesis and another RNA population for making cDNA has established the most important step of the BARE1 life cycle (III). The relationship between BARE1 and BARE2 has been investigated. Besides BARE, we have examined the retrotransposon Cassandra (II), which uses a very different transcriptional mechanism and a fully parasitic life cycle. In general, this work is focused on BARE1 promoter activity, transcriptional regulation including differential promoter usage and RNA pools, extra GAG protein production and VLP formation. The results of this study give new insights into transcription regulation of LTR retrotransposons.

Characterization of species-specific repeated DNA sequences from B. nigra

The construction and characterization of two genome-specific recombinant DNA clones from B. nigra are described. Southern analysis showed that the two clones belong to a dispersed repeat family. They differ from each other in their length, distribution and sequence, though the average GC content is nearly the same (45%). These B genome-specific repeats have been used to analyse the phylogenetic relationships between cultivated and wild species of the family Brassicaceae.

Conformational Polymorphism in Telomeric Structures: Loop Orientation and Interloop Pairing in d(G4TnG4),

Sequence repeats constituting the telomeric regions of chromosomes are known to adopt a variety of unusual structures, consisting of a G tetraplex stem and short stretches of thymines or thymines and adenines forming loops over the stem. Detailed model building and molecular mechanics studies have been carried out for these telomeric sequences to elucidate different types of loop orientations and possible conformations of thymines in the loop. The model building studies indicate that a minimum of two thymines have to be interspersed between guanine stretches to form folded-back structures with loops across adjacent strands in a G tetraplex (both over the small as well as large groove), while the minimum number of thymines required to build a loop across the diagonal strands in a G tetraplex is three. For two repeat sequences, these hairpins, resulting from different types of folding, can dimerize in three distinct ways-i.e., with loops across adjacent strands and on same side, with loops across adjacent strands and on opposite sides, and with loops across diagonal strands and on opposite sides-to form hairpin dimer structures. Energy minimization studies indicate that all possible hairpin dimers have very similar total energy values, though different structures are stabilized by different types of interactions. When the two loops are on the same side, in the hairpin dimer structures of d(G(4)T(n)G(4)), the thymines form favorably stacked tetrads in the loop region and there is interloop hydrogen bonding involving two hydrogen bonds for each thymine-thymine pair. Our molecular mechanics calculations on various folded-back as well as parallel tetraplex structures of these telomeric sequences provide a theoretical rationale for the experimentally observed feature that the presence of intervening thymine stretches stabilizes folded-back structures, while isolated stretches of guanines adopt a parallel tetraplex structure

Chain folding and A:T pairing in human telomeric DNA: a model-building and molecular dynamics study

The various types of chain folding and possible intraloop as well as interloop base pairing in human telomeric DNA containing d(TTAG(3)) repeats have been investigated by model-building, molecular mechanics, and molecular dynamics techniques. Model-building and molecular mechanics studies indicate that it is possible to build a variety of energetically favorable folded-back structures with the two TTA loops on same side and the 5' end thymines in the two loops forming TATA tetrads involving a number of different intraloop as well as interloop A:T pairing schemes. In these folded-back structures, although both intraloop and interloop Watson-Crick pairing is feasible, no structure is possible with interloop Hoogsteen pairing. MD studies of representative structures indicate that the guanine-tetraplex stem is very rigid and, while the loop regions are relatively much more flexible, most of the hydrogen bonds remain intact throughout the 350-ps in vacuo simulation. The various possible TTA loop structures, although they are energetically similar, have characteristic inter proton distances, which could give rise to unique cross-peaks in two-dimensional nuclear Overhauser effect spectroscopy (NOESY) experiments. These folded-back structures with A:T pairings in the loop region help in rationalizing the data from chemical probing and other biochemical studies on human telomeric DNA.

The carboxy terminal WD domain of the pre-mRNA splicing factor Prp17p is critical for function

In Saccharomyces cerevisiae, Prp17p is required for the efficient completion of the second step of pre-mRNA splicing. The function and interacting factors for this protein have not been elucidated. We have performed a mutational analysis of yPrp17p to identify protein domains critical for function. A series of deletions were made throughout the region spanning the N-terminal 158 amino acids of the protein, which do not contain any identified structural motifs. The C-terminal portion (amino acids 160–455) contains a WD domain containing seven WD repeats. We determined that a minimal functional Prp17p consists of the WD domain and 40 amino acids N-terminal to it. We generated a three-dimensional model of the WD repeats in Prp17p based on the crystal structure of the [beta]-transducin WD domain. This model was used to identify potentially important amino acids for in vivo functional characterization. Through analysis of mutations in four different loops of Prp17p that lie between [beta] strands in the WD repeats, we have identified four amino acids, 235TETG238, that are critical for function. These amino acids are predicted to be surface exposed and may be involved in interactions that are important for splicing. Temperature-sensitive prp17 alleles with mutations of these four amino acids are defective for the second step of splicing and are synthetically lethal with a U5 snRNA loop I mutation, which is also required for the second step of splicing. These data reinforce the functional significance of this region within the WD domain of Prp17p in the second step of splicing.

Mutations in WDR62, encoding a centrosomal and nuclear protein, in Indian primary microcephaly families with cortical malformations

Primary microcephaly is an autosomal recessive disorder characterized by smaller than normal brain size and mental retardation. It is genetically heterogeneous with seven loci: MCPH1-MCPH7. We have previously reported genetic analysis of 35 families, including the identification of the MCPH7 gene STIL. Of the 35 families, three families showed linkage to the MCPH2 locus. Recent whole-exome sequencing studies have shown that the WDR62 gene, located in the MCPH2 candidate region, is mutated in patients with severe brain malformations. We therefore sequenced the WDR62 gene in our MCPH2 families and identified two novel homozygous protein truncating mutations in two families. Affected individuals in the two families had pachygyria, microlissencephaly, band heterotopias, gyral thickening, and dysplastic cortex. Using immunofluorescence study, we showed that, as with other MCPH proteins, WDR62 localizes to centrosomes in A549, HepG2, and HaCaT cells. In addition, WDR62 was also localized to nucleoli. Bioinformatics analysis predicted two overlapping nuclear localization signals and multiple WD-40 repeats in WDR62. Two other groups have also recently identified WDR62 mutations in MCPH2 families. Our results therefore add further evidence that WDR62 is the MCPH2 gene. The present findings will be helpful in genetic diagnosis of patients linked to the MCPH2 locus.

A Monoclonal Antibody against Human Notch1 Ligand-Binding Domain Depletes Subpopulation of Putative Breast Cancer Stem-like Cells

Overexpression of Notch receptors and ligands has been associated with various cancers and developmental disorders, making Notch a potential therapeutic target. Here, we report characterization of Notch1 monoclonal antibodies (mAb) with therapeutic potential. The mAbs generated against epidermal growth factor (EGF) repeats 11 to 15 inhibited binding of Jagged1 and Delta-like4 and consequently, signaling in a dose-dependent manner, the antibodies against EGF repeats 11 to 12 being more effective than those against repeats 13 to 15. These data emphasize the role of EGF repeats 11 to 12 in ligand binding. One of the mAbs, 602.101, which specifically recognizes Notch1, inhibited ligand-dependent expression of downstream target genes of Notch such as HES-1, HES-5, and HEY-L in the breast cancer cell line MDA-MB-231. The mAb also decreased cell proliferation and induced apoptotic cell death. Furthermore, exposure to this antibody reduced CD44(Hi)/CD24(Low) subpopulation in MDA-MB-231 cells, suggesting a decrease in the cancer stem-like cell subpopulation. This was confirmed by showing that exposure to the antibody decreased the primary, secondary, and tertiary mammosphere formation efficiency of the cells. Interestingly, effect of the antibody on the putative stem-like cells appeared to be irreversible, because the mammosphere-forming efficiency could not be salvaged even after antibody removal during the secondary sphere formation. The antibody also modulated expression of genes associated with stemness and epithelial-mesenchymal transition. Thus, targeting individual Notch receptors by specific mAbs is a potential therapeutic strategy to reduce the potential breast cancer stem-like cell subpopulation. Mol Cancer Ther; 11(1); 77-86. (C) 2011 AACR.

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.

As(3)p: A Fast Algorithm to Search Structurally Similar Protiens

Protein structure comparison is essential for understanding various aspects of protein structure, function and evolution. It can be used to explore the structural diversity and evolutionary patterns of protein families. In view of the above, a new algorithm is proposed which performs faster protein structure comparison using the peptide backbone torsional angles. It is fast, robust, computationally less expensive and efficient in finding structural similarities between two different protein structures and is also capable of identifying structural repeats within the same protein molecule.

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.

RepEx: Repeat extractor for biological sequences

Genomic sequences are far from being random but are made up of systematically ordered and information rich patterns. These repeated sequence patterns have been vastly utilized for their fundamental importance in understanding the genome function and organization. To this end, a comprehensive toolkit, RepEx, has been developed which extracts repeat (inverted, everted and mirror) patterns from the given genome sequence(s) without any constraints. The toolkit can also be used to fetch the inverted repeats present in the protein sequence (s). Further, it is capable of extracting exact and degenerate repeats with a user defined spacer intervals. It is remarkably more precise and sensitive when compared to the existing tools. An example with comprehensive case studies and a performance evaluation of the proposed toolkit has been presented to authenticate its efficiency and accuracy. (C) 2013 Elsevier Inc. All rights reserved.

Differential proteomic and genomic profiling of mouse striatal cell model of Huntington's disease and control; probable implications to the disease biology

Huntington's disease (HD) is an autosomal dominant disorder of central nervous system caused by expansion of CAG repeats in exon1 of the huntingtin gene (Htt). Among various dysfunctions originated from the mutation in Htt gene, transcriptional deregulation has been considered to be one of the most important abnormalities. Large numbers of investigations identified altered expressions of genes in brains of HD patients and many models of HD. In this study we employed 2D SDS-PAGE/MALDI-MS coupled with 2D-DIGE and real-time PCR experiments of an array of genes focused to HD pathway to determine altered protein and gene expressions in STHdh(Q111)/Hdh(Q111) cells, a cell model of HD and compared with STHdh(Q7)/Hdh(Q7) cells, its wild type counterpart. We annotated 76 proteins from these cells and observed differential expressions of 31 proteins (by 2D-DIGE) involved in processes like unfolded protein binding, negative regulation of neuron apoptosis, response to superoxides etc. Our PCR array experiments identified altered expressions of 47 genes. Altogether significant alteration of 77 genes/proteins could be identified in this HD cell line with potential relevance to HD biology. Biological significance: In this study we intended to find out differential proteomic and genomic profiles in HD condition. We used the STHdh cells, a cellular model for HD and control. These are mouse striatal neuronal cell lines harboring 7 and 111 knock -in CAG repeats in their two alleles. The 111Q containing cell line (STHdh(Q111)/Hdh(Q111)) mimics diseased condition, whereas the 7Q containing ones (STHdh(Q7)/Hdh(Q7)), serves as the proper control cell line. Proteomic experiments were performed earlier to obtain differential expressions of proteins in R6/2 mice models, Hdh(Q) knock -in mice and in plasma and CSF from HD patients. However, no earlier report on proteomic alterations in these two HD cell lines and control was available in literature. It was, therefore, an important objective to find out differential expressions of proteins in these two cell lines. In this study, we annotated 76 proteins from STHdh(Q7)/Hdh(Q7) and STHdh(Q111)/Hdh(Q111) cells using 2D-gel/mass spectrometry. Next, by performing 2D-DIGE, we observed differential expressions of 31 proteins (16 upregulated and 15 downregulated) between these two cell lines. We also performed customized qRT-PCR array focused to HD pathway and found differential expressions of 47 genes (8 gene exptessions increased and 39 genes were decreased significantly). A total of 77 genes/proteins (Htt downregulated in both the studies) were found to be significantly altered from both the experimental paradigms. We validated the differential expressions of Vim, Hypk, Ran, Dstn, Hspa5 and Sod2 either by qRT-PCR or Western blot analysis or both. Out of these 77, similar trends in alteration of 19 out of 31 and 38 out of 47 proteins/genes were reported in earlier studies. Thus our study confirmed earlier observations on differential gene/protein expressions in HD and are really useful. Additionally, we observed differential expression of some novel genes/proteins. One of this was Hypk, a Htt-interacting chaperone protein with the ability to solubilize mHtt aggregated structures in cell lines. We propose that downregulation of Hypk in STHdh-Qm (Q111)/Hdh(Q111) has a causal effect towards HD pathogenesis. Thus the novel findings from our study need further research and might be helpful to understand the molecular mechanism behind HD pathogenesis. (C) 2015 Elsevier B.V. All rights reserved.

Repeat-Associated Fission Yeast-Like Regional Centromeres in the Ascomycetous Budding Yeast Candida tropicalis

The centromere, on which kinetochore proteins assemble, ensures precise chromosome segregation. Centromeres are largely specified by the histone H3 variant CENP-A (also known as Cse4 in yeasts). Structurally, centromere DNA sequences are highly diverse in nature. However, the evolutionary consequence of these structural diversities on de novo CENP-A chromatin formation remains elusive. Here, we report the identification of centromeres, as the binding sites of four evolutionarily conserved kinetochore proteins, in the human pathogenic budding yeast Candida tropicalis. Each of the seven centromeres comprises a 2 to 5 kb non-repetitive mid core flanked by 2 to 5 kb inverted repeats. The repeat-associated centromeres of C. tropicalis all share a high degree of sequence conservation with each other and are strikingly diverged from the unique and mostly non-repetitive centromeres of related Candida species-Candida albicans, Candida dubliniensis, and Candida lusitaniae. Using a plasmid-based assay, we further demonstrate that pericentric inverted repeats and the underlying DNA sequence provide a structural determinant in CENP-A recruitment in C. tropicalis, as opposed to epigenetically regulated CENP-A loading at centromeres in C. albicans. Thus, the centromere structure and its influence on de novo CENP-A recruitment has been significantly rewired in closely related Candida species. Strikingly, the centromere structural properties along with role of pericentric repeats in de novo CENP-A loading in C. tropicalis are more reminiscent to those of the distantly related fission yeast Schizosaccharomyces pombe. Taken together, we demonstrate, for the first time, fission yeast-like repeat-associated centromeres in an ascomycetous budding yeast.