Mining microsatellite markers from public expressed sequence tags databases for the study of threatened plants

Background Simple Sequence Repeats (SSRs) are widely used in population genetic studies but their classical development is costly and time-consuming. The ever-increasing available DNA datasets generated by high-throughput techniques offer an inexpensive alternative for SSRs discovery. Expressed Sequence Tags (ESTs) have been widely used as SSR source for plants of economic relevance but their application to non-model species is still modest. Methods Here, we explored the use of publicly available ESTs (GenBank at the National Center for Biotechnology Information-NCBI) for SSRs development in non-model plants, focusing on genera listed by the International Union for the Conservation of Nature (IUCN). We also search two model genera with fully annotated genomes for EST-SSRs, Arabidopsis and Oryza, and used them as controls for genome distribution analyses. Overall, we downloaded 16 031 555 sequences for 258 plant genera which were mined for SSRsand their primers with the help of QDD1. Genome distribution analyses in Oryza and Arabidopsis were done by blasting the sequences with SSR against the Oryza sativa and Arabidopsis thaliana reference genomes implemented in the Basal Local Alignment Tool (BLAST) of the NCBI website. Finally, we performed an empirical test to determine the performance of our EST-SSRs in a few individuals from four species of two eudicot genera, Trifolium and Centaurea. Results We explored a total of 14 498 726 EST sequences from the dbEST database (NCBI) in 257 plant genera from the IUCN Red List. We identify a very large number (17 102) of ready-to-test EST-SSRs in most plant genera (193) at no cost. Overall, dinucleotide and trinucleotide repeats were the prevalent types but the abundance of the various types of repeat differed between taxonomic groups. Control genomes revealed that trinucleotide repeats were mostly located in coding regions while dinucleotide repeats were largely associated with untranslated regions. Our results from the empirical test revealed considerable amplification success and transferability between congenerics. Conclusions The present work represents the first large-scale study developing SSRs by utilizing publicly accessible EST databases in threatened plants. Here we provide a very large number of ready-to-test EST-SSR (17 102) for 193 genera. The cross-species transferability suggests that the number of possible target species would be large. Since trinucleotide repeats are abundant and mainly linked to exons they might be useful in evolutionary and conservation studies. Altogether, our study highly supports the use of EST databases as an extremely affordable and fast alternative for SSR developing in threatened plants.

Isolation of an active gene and of two pseudogenes for mouse U7 small nuclear RNA

Three U7 RNA-related sequences were isolated from mouse genomic DNA libraries. Only one of the sequences completely matches the published mouse U7 RNA sequence, whereas the other two apparently represent pseudogenes. The matching sequence represents a functional gene, as it is expressed after microinjection into Xenopus laevis oocytes. Sequence variations of the conserved cis-acting 5' and 3' elements of U RNA genes may partly explain the low abundance of U7 RNA.

Comparative analysis of Tritrichomonas foetus (Riedmüller, 1928) cat genotype, T. foetus (Riedmüller, 1928) cattle genotype and Tritrichomonas suis (Davaine, 1875) at 10 DNA loci.

The parasitic protists in the genus Tritrichomonas cause significant disease in domestic cattle and cats. To assess the genetic diversity of feline and bovine isolates of Tritrichomonas foetus (Riedmüller, 1928) Wenrich and Emmerson, 1933, we used 10 different genetic regions, namely the protein coding genes of cysteine proteases 1, 2 and 4-9 (CP1, 2, 4-9) involved in the pathogenesis of the disease caused by the parasite. The cytosolic malate dehydrogenase 1 (MDH1) and internal transcribed spacer region 2 of the rDNA unit (ITS2) were included as additional markers. The gene sequences were compared with those of Tritrichomonas suis (Davaine, 1875) Morgan and Hawkins, 1948 and Tritrichomonas mobilensisCulberson et al., 1986. The study revealed 100% identity for all 10 genes among all feline isolates (=T. foetus cat genotype), 100% identity among all bovine isolates (=T. foetus cattle genotype) and a genetic distinctness of 1% between the cat and cattle genotypes of T. foetus. The cattle genotype of T. foetus was 100% identical to T. suis at nine loci (CP1, 2, 4-8, ITS2, MDH1). At CP9, three out of four T. suis isolates were identical to the T. foetus cattle genotype, while the T. suis isolate SUI-H3B sequence contained a single unique nucleotide substitution. Tritrichomonas mobilensis was 0.4% and 0.7% distinct from the cat and cattle genotypes of T. foetus, respectively. The genetic differences resulted in amino acid changes in the CP genes, most pronouncedly in CP2, potentially providing a platform for elucidation of genotype-specific host-pathogen interactions of T. foetus. On the basis of this data we judge T. suis and T. foetus to be subjective synonyms. For the first time, on objective nomenclatural grounds, the authority of T. suis is given to Davaine, 1875, rather than the commonly cited Gruby and Delafond, 1843. To maintain prevailing usage of T. foetus, we are suppressing the senior synomym T. suisDavaine, 1875 according to Article 23.9, because it has never been used as a valid name after 1899 and T. foetus is widely discussed as the cause of bovine trichomonosis. Thus bovine, feline and porcine isolates should all be given the name T. foetus. This promotes the stability of T. foetus for the veterinary and economically significant venereal parasite causing bovine trichomonosis.

Structures of four human pseudogenes for U7 small nuclear RNA.

Four U7 RNA-related sequences were isolated from a human genomic DNA library. None of the sequences completely match the published human U7 RNA sequence and all of them contain features typical of reverse-transcribed pseudogenes.

Histone-specific RNA 3' processing in nuclear extracts from mammalian cells.

The mature 3' ends of histone mRNAs are formed by endonucleolytic cleavage of longer precursor transcripts. This process occurs in the nucleus and can be regarded as the equivalent of the polyadenylation reaction involved in 3′-end-generation of all other mRNAs. A sea urchin H3 gene that failed to be properly processed in the Xenopus oocyte system proved particularly useful, because it allowed the identification of a processing component from sea urchins by a complementation assay. Nuclear extracts prepared from cells under various growth conditions have helped to reveal proliferation-dependent changes in the efficiency of histone RNA 3′ processing. RNA substrates for in vitro processing are best prepared by runoff transcription of specific DNA templates with bacterial or phage RNA polymerases. For this purpose, a restriction fragment containing the 3′-terminal region of a histone gene and including the conserved palindrome and spacer motifs is cloned into a polylinker sequence downstream of a strong promoter.

Cycling in the nucleus: regulation of RNA 3' processing and nuclear organization of replication-dependent histone genes.

The histones which pack new DNA during the S phase of animal cells are made from mRNAs that are cleaved at their 3' end but not polyadenylated. Some of the factors used in this reaction are unique to it while others are shared with the polyadenylation process that generates all other mRNAs. Recent work has begun to shed light on how the cell manages the assignment of these common components to the two 3' processing systems, and how it achieves their cell cycle-regulation and recruitment to the histone pre-mRNA. Moreover, recent and older findings reveal multiple connections between the nuclear organization of histone genes, their transcription and 3' end processing as well as the control of cell proliferation.

Development and characterization of novel microsatellite markers for Arion slug species

Seventeen polymorphic microsatellite markers were isolated and characterized in Arion vulgaris/lusitanicus, which belongs to the worst European slug pests with serious economic and ecological impact. These markers were tested on 23 individuals collected in a population in Switzerland. Numbers of alleles ranged from 2 to 14 per locus, observed and expected heterozygosities ranged from 0.174 to 0.87, and from 0.162 to 0.903, respectively. These loci were also successfully amplified and were polymorphic in the closely related species A. rufus and A. ater. These loci represent the first highly polymorphic nuclear markers described for A. vulgaris and pave the way for population genetics and molecular ecology research of the important Arion pest slugs.

Assessing the Phylogenetic Utility of DNA Barcoding Using the New Zealand Cicada Genus Kikihia

DNA Barcoding (Hebert et al. 2003) has the potential to revolutionize the process of identifying and cataloguing biodiversity; however, significant controversy surrounds some of the proposed applications. In the seven years since DNA barcoding was introduced, the Web of Science records more than 600 studies that have weighed the pros and cons of this procedure. Unfortunately, the scientific community has been unable to come to any consensus on what threshold to use to differentiate species or even whether the barcoding region provides enough information to serve as an accurate species identification tool. The purpose of my thesis is to analyze mitochondrial DNA (mtDNA) barcoding’s potential to identify known species and provide a well-resolved phylogeny for the New Zealand cicada genus Kikihia. In order to do this, I created a phylogenetic tree for species in the genus Kikihia based solely on the barcoding region and compared it to a phylogeny previously created by Marshall et al. (2008) that benefits from information from other mtDNA and nuclear genes as well as species-specific song data. I determined how well the barcoding region delimits species that have been recognized based on morphology and song. In addition, I looked at the effect of sampling on the success of barcoding studies. I analyzed subsets of a larger, more densely sampled dataset for the Kikihia Muta Group to determine which aspects of my sampling strategy led to the most accurate identifications. Since DNA barcoding would by definition have problems in diagnosing hybrid individuals, I studied two species (K. “murihikua” and K. angusta) that are known to hybridize. Individuals that were not obvious hybrids (determined by morphology) were selected for the case study. Phylogenetic analysis of the barcoding region revealed insights into the reasons these two species could not be successfully differentiated using barcoding alone.

Elucidation of the role of 53BP1 in DNA double strand break (DSB) repair and tumor suppression

The protein p53 binding protein one (53BP1) was discovered in a yeast two-hybrid screen that used the DNA binding domain of p53 as bait. Cloning of full-length 53BP1 showed that this protein contains several protein domains which help make up the protein, which include two tandem BRCT domains and a amino-terminal serine/glutamine cluster domain (SCD). These are two protein domains are often seen in factors that are involved in the cellular response to DNA damage and control of cell cycle checkpoints and we hypothesize that 53BP1 is involved in the cellular response to DNA damage. In support of this hypothesis we observe that 53BP1 is phosphorylated and undergoes a dramatic nuclear re-localization in response to DNA damaging agents. 53BP1 also interacts with several factors that are important in the cellular response to DNA damage, such as the BRCA1 tumor suppressor, ATM and Rad3 related (ATR), and the phosphorylated version of the histone variant H2AX. Mice deficient in 53BP1 display increased sensitivity ionizing radiation (IR), a DNA damaging agent that introduces DNA double strand breaks (DSBs). In addition, 53BP1-deficient mice do not properly undergo the process of class switch recombination (CSR). We also observe that when a defect in 53BP1 is combined with a defect in p53; the resulting mice have an increased rate of formation of spontaneous tumors, notably the formation of B and T lineage lymphomas. The T lineage tumors arise by two distinct mechanisms: one driven by defects in cell cycle regulation and a second driven by defects in the ability to repair DNA DSBs. The B lineage tumors arise by the inability to repair DNA damage and over-expression of the oncogene c-myc. ^ With these observations, we conclude that not only does 53BP1 function in the cellular response to DNA damage, but it also works in concert with p53 to suppress tumor formation. ^

The molecular mechanisms of sensing nucleoside analogue-induced DNA damage

Nucleoside analogues are antimetabolites effective in the treatment of a wide variety of solid tumors and hematological malignancies. Upon being metabolized to their active triphosphate form, these agents are incorporated into DNA during replication or excision repair synthesis. Because DNA polymerases have a greatly decreased affinity for primers terminated by most nucleoside analogues, their incorporation causes stalling of replication forks. The molecular mechanisms that recognize blocked replication may contribute to drug resistance but have not yet been elucidated. Here, several molecules involved in sensing nucleoside analogue-induced stalled replication forks have been identified and examined for their contribution to drug resistance. ^ The phosphorylation of the DNA damage sensor, H2AX, was characterized in response to nucleoside analogues and found to be dependent on both time and drug concentration. This response was most evident in the S-phase fraction and was associated with an inhibition of DNA synthesis, S-phase accumulation, and activation of the S-phase checkpoint pathway (Chk1-Cdc25A-Cdk2). Exposure of the Chk1 inhibitor, 7-hydroxystaurosporine (UCN-01), to cultures previously treated with nucleoside analogues caused increased apoptosis, clonogenic death, and a further log-order increase in H2AX phosphorylation, suggesting enhanced DNA damage. Ataxia-telangiectasia mutated (ATM) has been identified as a key DNA damage signaling kinase for initiating cell cycle arrest, DNA repair, and apoptosis while the Mre11-Rad50-Nbs1 (MRN) complex is known for its functions in double-strand break repair. Activated ATM and the MRN complex formed distinct nuclear foci that colocalized with phosphorylated H2AX after inhibition of DNA synthesis by the nucleoside analogues, gemcitabine, ara-C, and troxacitabine. Since double-strand breaks were undetectable, this response was likely due to stalling of replication forks. A similar DNA damage response was observed in human lymphocytes after exposure to ionizing radiation and in acute myelogenous leukemia blasts during therapy with the ara-C prodrug, CP-4055. Deficiencies in ATM, Mre11, and Rad50 led to a two- to five-fold increase in gemcitabine sensitivity, suggesting that these molecules contribute to drug resistance. Based on these results, a model is proposed for the sensing of nucleoside analogue-induced stalled replication forks that includes H2AX, ATM, and the Mre11-Rad50-Nbs1 complex. ^

Clinical relevance of polymorphic DNA copy number variation (CNV) in African American women with stage I-II breast cancer

Objectives. The chief goal of this study was to analyze copy number variation (CNV) in breast cancer tumors from 25 African American women with early stage breast cancer (BC) using molecular inversion probes (MIP) in order to: (1) compare the degree of CNV in tumors compared to normal lymph nodes, and (2) determine whether gains and/or losses of genes in specific chromosomes differ between pathologic subtypes of breast cancer defined by known prognostic markers, (3) determine whether gains/losses in CN are associated with known oncogenes or tumor suppressor genes, and (4) determine whether increased gains/losses in CN for specific chromosomes were associated with differences in breast cancer recurrence. ^ Methods. Twenty to 37 nanograms of DNA extracted from 25 formalin-fixed paraffin embedded (FFPE) tumor samples and matched normal lymph nodes were added to individual tubes. Oligonucleotide probes with recognition sequences at each terminus were hybridized with a genomic target sequence to form a circular structure. Probes are released from genomic DNA obtained from FFPE samples, and those which have been correctly "circularized" in the proper allele/nucleotide reaction combination are amplified using polymerase chain reaction (PCR) primers. Amplicons were fluorescently labeled and the tag sequences released from the genome homology regions by treatment with uracil-N-glycosylase to cleave the probe at the site where uracils are present, and detected using a complementary tag array developed by Affymetrix. ^ Results. Analysis of CN gains and losses from tumors and normal tissues showed marked differences in tumors with numerous chromosomes affected. Similar changes were not observed in normal lymph nodes. When tumors were stratified into four groups based on expression or lack of expression of the estrogen receptor and HER2/neu, distinct patterns of CNV for different chromosomes were observed. Gains or losses in CN for specific chromosomes correlated with amplifications/deletions of particular oncogenes or tumor suppressor genes (i.e. such as found on chromosome 17) known to be associated with aggressive tumor phenotype and poor prognosis. There was a trend for increases in CN observed for chromosome 17 to correlate inversely with time to recurrence of BC (p=0.14 for trend). CNV was also observed for chromosomes 5, 8, 10, 11, and 16, which are known sites for several breast cancer susceptibility alleles. ^ Conclusions. This study is the first to validate the MIP technique, to correlate differences in gene expression with known prognostic tumor markers, and to correlate significant increases/decreases in CN with known tumor markers associated with prognosis. The results of this study may have far reaching public health implications towards identifying new high-risk groups based on genomic differences in CNP, both with respect to prognosis and response to therapy, and to eventually identify new therapeutic targets for prevention and treatment of this disease. ^

The effect of DNA methylation in carcinogenesis

CpG island methylation within single gene promoters can silence expression of associated genes. We first extended these studies to bidirectional gene pairs controlled by single promoters. We showed that hypermethylation of bidirectional promoter-associated CpG island silences gene pairs (WNT9A/CD558500, CTDSPL/BC040563, and KCNK15/BF 195580) simultaneously. Hypomethylation of these promoters by 5-aza-2'-deoxycytidine treatment reactivated or enhanced gene expression bidirectionally. These results were further confirmed by luciferase assays. Methylation of WNT9A/CD558500 and CTDSPL/BC040563 promoters occurs frequently in primary colon cancers and acute lymphoid leukemia, respectively. ^ Next we sought to understand the origins of hypermethylation in cancer. CpG islands associated with tumor suppressor genes are normally free from methylation, but can be hypermethylated in cancer. It remains poorly understood how these genes are protected from methylation in normal tissues. In our studies, we aimed to determine if cis-acting elements in these genes are responsible for this protection, using the tumor suppressor gene p16 as a model. We found that Alu repeats located both upstream and downstream of the p16 promoter become hypermethylated with age. In colon cancer samples, the methylation level is particularly high, and the promoter can also be affected. Therefore, the protection in the promoter against methylation spreading could fail during tumorigenesis. This methylation pattern in p16 was also observed in cell lines of different tissue origins, and their methylation levels were found to be inversely correlated with that of active histone modification markers (H3K4-3me and H3K9-Ac). To identify the mechanism of protection against methylation spreading, we constructed serial deletions of the p16 protected region and used silencing of a neomycin reporter gene to evaluate the protective effects of these fragments. A 126 bp element was identified within the region which exerts bidirectional protection against DNA methylation, independently of its transcriptional activity. The protective strength of this element is comparable to that of the HS4 insulator. During long-term culture, the presence of this element significantly slowed methylation spreading. In conclusion, we have found that an element located in the p16 promoter is responsible for protection against DNA methylation spreading in normal tissues. The failure of protective cis-elements may be a general feature of tumor-suppressor gene silencing during tumorigenesis. ^

Regulation of squamous cell differentiation in human head and neck carcinoma cells by retinoids

Retinoids have been found to be effective in the prevention of premalignant lesions and second primary cancers in the upper aerodigestive tract. Further development of retinoids for prevention and therapy of head and neck squamous cell carcinoma (HNSCC) requires a better understanding of their mechanism of action on the growth and differentiation of such cells. I have chosen to employ cultured HNSCC cell lines as a model system for investigating the mechanism underlying the effects of retinoids. These cells are useful because all-trans retinoic acid (ATRA) inhibits their proliferation. Furthermore, two HNSCC cell lines were found to express three squamous differentiation (SqD) markers characteristic of normal keratinocytes and ATRA suppressed the expression of these markers as reported for normal keratinocytes. It is thought that nuclear retinoic acid receptors (RARs and RXRs), which act as DNA-binding transcription modulating factors, mediate the effects of retinoids on the growth and differentiation of normal and tumor cells. I found that all four cell lines examined expressed RAR-$\alpha ,$ RAR-$\tau ,$ and RXR-$\alpha$ and three of four expressed RAR-$\beta .$ ATRA treatment increased the level of RAR-$\alpha ,$ -$\beta ,$ and -$\tau$ in four cell lines. Two HNSCC cell lines that exhibited a progressive increase in the expression of SqD markers during growth in culture also showed a concurrent decrease in RAR-$\beta$ level. Moreover, increasing concentrations of RA suppressed the SqD marker while inducing RAR-$\beta$ mRNA. Several synthetic retinoids which exhibit a preference for binding to specific nuclear RARs showed a differential ability to inhibit cell proliferation, transactivate transcription of the reporter genes (CAT and luciferase) from the RA response element (RARE) of the RAR-$\beta$ gene, and induce RAR-$\beta$ expression. Those retinoids that were effective inducers of RAR-$\beta$ also suppressed SqD effectively, indicating an inverse relationship exists between the expression of RAR-$\beta$ and SqD. This inverse relationship suggests a role for RAR-$\beta$ in the suppression of SqD. ^

The effects of DNA cytosine methylation on H -ras promoter structure and function

The effect of DNA cytosine methylation on H-ras promoter activity was assessed using a transient expression system employing the plasmid H-rasCAT (NaeI H-ras promoter linked to the chloramphenicol acetyltransferase (CAT) gene). This 551 bp promoter is 80% GC rich, enriched with 168 CpG dinucleotides, and contains six functional GC box elements which represent major DNA methylation target sites. Prokaryotic methyltransferases HhaI (CGm$\sp5$CG) and HpaII (Cm$\sp5$CGG) alone or in combination with a human placental methyltransferase (HP MTase) were used to introduce methyl groups at different CpG sites within the promoter. To test for functional promoter activity, the methylated plasmids were introduced into CV-1 cells and CAT activity assessed 48 h post-transfection. Methylation at specific HhaI and HpaII sites reduced CAT expression by 70%, whereas more extensive methylation at generalized CpG sites with HP MTase inactivated the promoter $>$95%. The inhibition of H-ras promoter activity was not attributable to methylation-induced differences in DNA uptake or stability in the cell, topological form of the plasmid, or methylation effects in nonpromoter regions. We also observed that DNA cytosine methylation of a 360 bp promoter fragment by HP MTase induced a local change in DNA conformation. Using three independent methodologies (nitrocellulose filter binding assays, gel mobility shifts, and Southwestern blots), we determined that this change in promoter conformation affected the interaction of nuclear proteins with cis-regulatory sequences residing in the promoter region. The results provide evidence to suggest that DNA methylation may regulate gene expression by inducing changes in local promoter conformation which in turn alters the interactions between DNA and protein factors required for transcription. The results provide supportive evidence for the hypothesis of Cedar and Riggs, who postulated that DNA methylation may regulate gene expression by altering the binding affinities of proteins for DNA. ^

Characterization of the MLC box: A conserved {\it cis\/}-acting DNA element present in the chicken MLC1s/1c gene

The slow/cardiac alkali myosin light chain (MLC1s/1c) is a member of a multigene family whose protein products are essential for activation of the myosin ATPase. In the adult, the MLC1s/1c isoform is expressed in both cardiac and slow-twitch skeletal muscles, while it is expressed by all skeletal muscles during development.^ To elucidate the molecular mechanisms that underlie the transcriptional regulation of MLC1s/1c gene expression, the immediate 5$\sp\prime$ flanking region of the gene was isolated and shown to be capable of directing reporter gene expression. Analysis of this region revealed a 110 bp muscle-specific enhancer that includes a myocyte-specific enhancer-binding factor 2 (MEF-2) site, E-boxes, which are potential binding sites for the basic-helix-loop-helix proteins such as MyoD, and a MLC box. The focus of the thesis was to identify the role of the MLC box in expression of the MLC1s/1c gene.^ The MLC box is a member of the family of CArG box containing cis-acting DNA elements. Mutagenesis showed that the MLC box is necessary, but not sufficient, for the expression of a reporter gene linked to the 5$\sp\prime$ flanking region of the MLC1s/1c gene. Linker scanner and site-directed mutagenesis identified a number of potential sites within the 110 bp muscle-specific enhancer that may cooperate with the MLC box. These are the MEF-2 site, the E-box site, and a 10 bp element located upstream of the MEF-2 site that does not have sequence similarity with any known cis-acting element. The MLC box is capable of binding to factors present in muscle nuclear extracts, as well as to human recombinant serum response factor (SRF). Binding of SRF to the MLC box was correlated with the ability of the 5$\sp\prime$ flanking region of the MLC1s/1c gene to drive reporter gene expression. Results suggest a model in which binding of SRF to the MLC box activates expression of the MLC1s/1c gene while binding of the factors present in the nuclear extracts suppresses the expression of the gene. (Abstract shortened with permission of author.) ^