Instrumental variable method for the causal relationship between soluble receptor for advanced glycation end-products (sRAGE) and risk of pancreatic cancer

This thesis project is motivated by the potential problem of using observational data to draw inferences about a causal relationship in observational epidemiology research when controlled randomization is not applicable. Instrumental variable (IV) method is one of the statistical tools to overcome this problem. Mendelian randomization study uses genetic variants as IVs in genetic association study. In this thesis, the IV method, as well as standard logistic and linear regression models, is used to investigate the causal association between risk of pancreatic cancer and the circulating levels of soluble receptor for advanced glycation end-products (sRAGE). Higher levels of serum sRAGE were found to be associated with a lower risk of pancreatic cancer in a previous observational study (255 cases and 485 controls). However, such a novel association may be biased by unknown confounding factors. In a case-control study, we aimed to use the IV approach to confirm or refute this observation in a subset of study subjects for whom the genotyping data were available (178 cases and 177 controls). Two-stage IV method using generalized method of moments-structural mean models (GMM-SMM) was conducted and the relative risk (RR) was calculated. In the first stage analysis, we found that the single nucleotide polymorphism (SNP) rs2070600 of the receptor for advanced glycation end-products (AGER) gene meets all three general assumptions for a genetic IV in examining the causal association between sRAGE and risk of pancreatic cancer. The variant allele of SNP rs2070600 of the AGER gene was associated with lower levels of sRAGE, and it was neither associated with risk of pancreatic cancer, nor with the confounding factors. It was a potential strong IV (F statistic = 29.2). However, in the second stage analysis, the GMM-SMM model failed to converge due to non- concaveness probably because of the small sample size. Therefore, the IV analysis could not support the causality of the association between serum sRAGE levels and risk of pancreatic cancer. Nevertheless, these analyses suggest that rs2070600 was a potentially good genetic IV for testing the causality between the risk of pancreatic cancer and sRAGE levels. A larger sample size is required to conduct a credible IV analysis.^

Evolution of cis-regulatory elements in a repetitive sequence adjacent to a sea urchin aboral ectoderm -specific gene

The creation, preservation, and degeneration of cis-regulatory elements controlling developmental gene expression are fundamental genome-level evolutionary processes about which little is known. In this study, critical differences in cis-regulatory elements controlling the expression of the sea urchin aboral ectoderm-specific spec genes were identified and explored. In genomes of species within the Strongylocentrotidae family, multiple copies of a repetitive sequence element termed RSR were present, but RSRs were not detected in genomes of species outside Strongylocentrotidae. RSRs are invariably associated with spec genes, and in Strongylocentrotus purpuratus, the spec2a RSR functioned as a transcriptional enhancer displaying greater activity than RSRs from the spec1 or spec2c paralogs. Single base-pair differences at two cis-regulatory elements within the spec2a RSR greatly increased the binding affinities of four transcription factors: SpCCAAT-binding factor at one element and SpOtx, SpGoosecoid, and SpGATA-E at another. The cis-regulatory elements to which SpCCAAT-binding factor, SpOtx, SpGoosecoid, and SpGATA-E bound were recent evolutionary acquisitions that could act either to activate or repress transcription, depending on the cell type. These elements were found in the spec2a RSR ortholog in Strongylocentrotus pallidus but not in the RSR orthologs of Strongylocentrotus droebachiensis or Hemicentrotus pulcherrimus. These results indicate that spec genes exhibit a dynamic pattern of cis-regulatory element evolution while stabilizing selection preserves their aboral ectoderm expression domain. ^

Metagenomic Anlaysis of microsymbiont selection by the legume plant host

Rhizobium leguminosarum bv.viciae is able to establish nitrogen-fixing symbioses with legumes of the genera Pisum, Lens, Lathyrus and Vicia. Classic studies using trap plants (Laguerre et al., Young et al.) provided evidence that different plant hosts are able to select different rhizobial genotypes among those available in a given soil. However, these studies were necessarily limited by the paucity of relevant biodiversity markers. We have now reappraised this problem with the help of genomic tools. A well-characterized agricultural soil (INRA Bretennieres) was used as source of rhizobia. Plants of Pisum sativum, Lens culinaris, Vicia sativa and V. faba were used as traps. Isolates from 100 nodules were pooled, and DNA from each pool was sequenced (BGI-Hong Kong; Illumina Hiseq 2000, 500 bp PE libraries, 100 bp reads, 12 Mreads). Reads were quality filtered (FastQC, Trimmomatic), mapped against reference R. leguminosarum genomes (Bowtie2, Samtools), and visualized (IGV). An important fraction of the filtered reads were not recruited by reference genomes, suggesting that plant isolates contain genes that are not present in the reference genomes. For this study, we focused on three conserved genomic regions: 16S-23S rDNA, atpD and nodDABC, and a Single Nucleotide Polymorphism (SNP) analysis was carried out with meta / multigenomes from each plant. Although the level of polymorphism varied (lowest in the rRNA region), polymorphic sites could be identified that define the specific soil population vs. reference genomes. More importantly, a plant-specific SNP distribution was observed. This could be confirmed with many other regions extracted from the reference genomes (data not shown). Our results confirm at the genomic level previous observations regarding plant selection of specific genotypes. We expect that further, ongoing comparative studies on differential meta / multigenomic sequences will identify specific gene components of the plant-selected genotypes

Effect of the Legume Plant Host on Rhizobial Soil Population Dynamics

Rhizobium leguminosarum bv viciae (Rlv) is a soil bacterium able to establish specific root-nodule symbioses with legumes of four different genera: Pisum, Vicia, Lens and Lathyrus. Rlv isolates from nodules of any of these legumes can nodulate any of them; however, it has been shown that plants select specific rhizobial genotypes from those present in the soil (1,2). We have previously shown this at the genomic level by following a population genomics approach. Pool genomic sequences from 100 isolates from each of four plant species: P. sativum, L. culinaris, V. faba and V. sativa, show different, specific profiles at the single nucleotide polymorphism (SNP) level for relevant genes. In this work, the extent of Rlv selection from a well-characterized soil population by different legume plant hosts: P. sativum, L. culinaris, V. faba and V. sativa, after a medium-term mesocosm study is described. Direct soil isolates from each of these mesocosm studies have been tested for specific rhizobial genes (glnII and fnrN) and symbiotic genes (nodC and nifH). Different populations were characterized further by Sanger sequencing of both the rpoB phylogenetic marker gene and the symbiotic genes nodC and nifH. The distribution and size of the rhizobial population for each legume host showed changes during the medium-term mesocosm study. Particularly, a non-symbiotic group of rhizobia was enriched by all four hosts, in contrast to the symbiotic rhizobia profile, which was specific for each legume plant host.

The Pangenome of the Plant-Selected Rhizobium leguminosarum bv viciae Sub-popuations

Rhizobium leguminosarum bv viciae (Rlv) is a bacterium able to establish effective symbioses with four different legume genera: Pisum, Lens, Lathyrus and Vicia. Classic studies using trap plants have previously shown that, given a choice, different plants prefer specific genotypes of rhizobia, which are adapted to the host (1, 2). In previous work we have performed a Pool-Seq analysis bases on pooled DNA samples from Rlv nodule isolates obtained from Pisum sativum, Lens culinaris, Vicia fava and V. sativa plants, used as rhizobial traps. This experiment allowed us to test the host preference hypothesis: different plant hosts select specific sub-populations of rhizobia from the available population present in a given soil. We have observed that plant-selected sub-populations are different at the single nucleotide polymorphism (SNP) level. We have selected individual isolates from each sub-population (9 fava-bean isolates, 14 pea isolates 9 vetch isolates and 9 lentil isolates) and sequenced their genomes at draft level (ca. 30x, 90 contigs). Genomic analyses have been carried out using J-species and CMG-Biotools. All the isolates had similar genome size (7.5 Mb) and number of genes (7,300). The resulting Average Nucleotide Identity (ANIm) tree showed that Rhizobium leguminosarum bv viciae is a highly diverse group. Each plant-selected subpopulation showed a closed pangenome and core genomes of similar size (11,500 and 4,800 genes, respectively). The addition of all four sub-population results in a larger, closed pangenome of 19,040 genes and a core genome of similar size (4,392 genes). Each sub-population contains a characteristic set of genes but no universal, plant-specific genes were found. The core genome obtained from all four sub-populations is probably a representative core genome for Rhizobium leguminosarum, given that the reference genome (Rhizobium leguminosarum bv. viciae strain 3841) contains most of the core genome. We have also analyzed the symbiotic cluster (nod), and different nod cluster genotypes were found in each sub-population. Supported by MINECO (Consolider-Ingenio 2010, MICROGEN Project, CSD2009-00006).

Mendel’s dwarfing gene: cDNAs from the Le alleles and function of the expressed proteins

The major gibberellin (GA) controlling stem elongation in pea (Pisum sativum L.) is GA1, which is formed from GA20 by 3β-hydroxylation. This step, which limits GA1 biosynthesis in pea, is controlled by the Le locus, one of the original Mendelian loci. Mutations in this locus result in dwarfism. We have isolated cDNAs encoding a GA 3β-hydroxylase from lines of pea carrying the Le, le, le-3, and led alleles. The cDNA sequences from le and le-3 each contain a base substitution resulting in single amino acid changes relative to the sequence from Le. The cDNA sequence from led, a mutant derived from an le line, contains both the le “mutation” and a single-base deletion, which causes a shift in reading frame and presumably a null mutation. cDNAs from each line were expressed in Escherichia coli. The expression product for the clone from Le converted GA9 to GA4, and GA20 to GA1, with Km values of 1.5 μM and 13 μM, respectively. The amino acid substitution in the clone from le increased Km for GA9 100-fold and reduced conversion of GA20 to almost nil. Expression products from le and le-3 possessed similar levels of 3β-hydroxylase activity, and the expression product from led was inactive. Our results suggest that the 3β-hydroxylase cDNA is encoded by Le. Le transcript is expressed in roots, shoots, and cotyledons of germinating pea seedlings, in internodes and leaves of established seedlings, and in developing seeds.

β-Catenin mutations in cell lines established from human colorectal cancers

β-catenin has functions as both an adhesion and a signaling molecule. Disruption of these functions through mutations of the β-catenin gene (CTNNB1) may be important in the development of colorectal tumors. We examined the entire coding sequence of β-catenin by reverse transcriptase–PCR (RT-PCR) and direct sequencing of 23 human colorectal cancer cell lines from 21 patients. In two cell lines, there was apparent instability of the β-catenin mRNA. Five different mutations (26%) were found in the remaining 21cell lines (from 19 patients). A three-base deletion (codon 45) was identified in the cell line HCT 116, whereas cell lines SW 48, HCA 46, CACO 2, and Colo 201 each contained single-base missense mutations (codons 33, 183, 245, and 287, respectively). All 23 cell lines had full-length β-catenin protein that was detectable by Western blotting and that coprecipitated with E-cadherin. In three of the cell lines with CTNNB1 mutations, complexes of β-catenin with α-catenin and APC were detectable. In SW48 and HCA 46, however, we did not detect complexes of β-catenin protein with α-catenin and APC, respectively. These results show that selection of CTNNB1 mutations occurs in up to 26% of colorectal cancers from which cell lines are derived. In these cases, mutation selection is probably for altered β-catenin function, which may significantly alter intracellular signaling and intercellular adhesion and may serve as a complement to APC mutations in the early stages of tumorigenesis.

Horizontal gene transfer and mutation: Ngrol genes in the genome of Nicotiana glauca

Ngrol genes (NgrolB, NgrolC, NgORF13, and NgORF14) that are similar in sequence to genes in the left transferred DNA (TL-DNA) of Agrobacterium rhizogenes have been found in the genome of untransformed plants of Nicotiana glauca. It has been suggested that a bacterial infection resulted in transformation of Ngrol genes early in the evolution of the genus Nicotiana. Although the corresponding four rol genes in TL-DNA provoked hairy-root syndrome in plants, present-day N. glauca and plants transformed with Ngrol genes did not exhibit this phenotype. Sequenced complementation analysis revealed that the NgrolB gene did not induce adventitious roots because it contained two point mutations. Single-base site-directed mutagenesis at these two positions restored the capacity for root induction to the NgrolB gene. When the NgrolB, with these two base substitutions, was positioned under the control of the cauliflower mosaic virus 35S promoter (P35S), transgenic tobacco plants exhibited morphological abnormalities that were not observed in P35s-RirolB plants. In contrast, the activity of the NgrolC gene may have been conserved after an ancient infection by bacteria. Discussed is the effect of the horizontal gene transfer of the Ngrol genes and mutations in the NgrolB gene on the phenotype of ancient plants during the evolution of N. glauca.

In silico detection of control signals: mRNA 3′-end-processing sequences in diverse species

We have investigated mRNA 3′-end-processing signals in each of six eukaryotic species (yeast, rice, arabidopsis, fruitfly, mouse, and human) through the analysis of more than 20,000 3′-expressed sequence tags. The use and conservation of the canonical AAUAAA element vary widely among the six species and are especially weak in plants and yeast. Even in the animal species, the AAUAAA signal does not appear to be as universal as indicated by previous studies. The abundance of single-base variants of AAUAAA correlates with their measured processing efficiencies. As found previously, the plant polyadenylation signals are more similar to those of yeast than to those of animals, with both common content and arrangement of the signal elements. In all species examined, the complete polyadenylation signal appears to consist of an aggregate of multiple elements. In light of these and previous results, we present a broadened concept of 3′-end-processing signals in which no single exact sequence element is universally required for processing. Rather, the total efficiency is a function of all elements and, importantly, an inefficient word in one element can be compensated for by strong words in other elements. These complex patterns indicate that effective tools to identify 3′-end-processing signals will require more than consensus sequence identification.

Direct repeat sequences in the Streptomyces chitinase-63 promoter direct both glucose repression and chitin induction

The chi63 promoter directs glucose-sensitive, chitin-dependent transcription of a gene involved in the utilization of chitin as carbon source. Analysis of 5′ and 3′ deletions of the promoter region revealed that a 350-bp segment is sufficient for wild-type levels of expression and regulation. The analysis of single base changes throughout the promoter region, introduced by random and site-directed mutagenesis, identified several sequences to be important for activity and regulation. Single base changes at −10, −12, −32, −33, −35, and −37 upstream of the transcription start site resulted in loss of activity from the promoter, suggesting that bases in these positions are important for RNA polymerase interaction. The sequences centered around −10 (TATTCT) and −35 (TTGACC) in this promoter are, in fact, prototypical of eubacterial promoters. Overlapping the RNA polymerase binding site is a perfect 12-bp direct repeat sequence. Some base changes within this direct repeat resulted in constitutive expression, suggesting that this sequence is an operator for negative regulation. Other base changes resulted in loss of glucose repression while retaining the requirement for chitin induction, suggesting that this sequence is also involved in glucose repression. The fact that cis-acting mutations resulted in glucose resistance but not inducer independence rules out the possibility that glucose repression acts exclusively by inducer exclusion. The fact that mutations that affect glucose repression and chitin induction fall within the same direct repeat sequence module suggests that the direct repeat sequence facilitates both chitin induction and glucose repression.

Structural basis of DNA bending and oriented heterodimer binding by the basic leucine zipper domains of Fos and Jun

Interactions among transcription factors that bind to separate sequence elements require bending of the intervening DNA and juxtaposition of interacting molecular surfaces in an appropriate orientation. Here, we examine the effects of single amino acid substitutions adjacent to the basic regions of Fos and Jun as well as changes in sequences flanking the AP-1 site on DNA bending. Substitution of charged amino acid residues at positions adjacent to the basic DNA-binding domains of Fos and Jun altered DNA bending. The change in DNA bending was directly proportional to the change in net charge for all heterodimeric combinations between these proteins. Fos and Jun induced distinct DNA bends at different binding sites. Exchange of a single base pair outside of the region contacted in the x-ray crystal structure altered DNA bending. Substitution of base pairs flanking the AP-1 site had converse effects on the opposite directions of DNA bending induced by homodimers and heterodimers. These results suggest that Fos and Jun induce DNA bending in part through electrostatic interactions between amino acid residues adjacent to the basic region and base pairs flanking the AP-1 site. DNA bending by Fos and Jun at inverted binding sites indicated that heterodimers bind to the AP-1 site in a preferred orientation. Mutation of a conserved arginine within the basic regions of Fos and transversion of the central C:G base pair in the AP-1 site to G:C had complementary effects on the orientation of heterodimer binding and DNA bending. The conformational variability of the Fos–Jun–AP-1 complex may contribute to its functional versatility at different promoters.

In vitro selection for altered divalent metal specificity in the RNase P RNA

The ribozyme RNase P absolutely requires divalent metal ions for catalytic function. Multiple Mg2+ ions contribute to the optimal catalytic efficiency of RNase P, and it is likely that the tertiary structure of the ribozyme forms a specific metal-binding pocket for these ions within the active-site. To identify base moieties that contribute to catalytic metal-binding sites, we have used in vitro selection to isolate variants of the Escherichia coli RNase P RNA with altered specificities for divalent metal. RNase P RNA variants with increased activity in Ca2+ were enriched over 18 generations of selection for catalysis in the presence of Ca2+, which is normally disfavored relative to Mg2+. Although a wide spectrum of mutations was found in the generation-18 clones, only a single point mutation was common to all clones: a cytosine-to-uracil transition at position 70 (E. coli numbering) of RNase P. Analysis of the C70U point mutant in a wild-type background confirmed that the identity of the base at position 70 is the sole determinant of Ca2+ selectivity. It is noteworthy that C70 lies within the phylogenetically well conserved J3/4-P4-J2/4 region, previously implicated in Mg2+ binding. Our finding that a single base change is sufficient to alter the metal preference of RNase P is further evidence that the J3/4-P4-J2/4 domain forms a portion of the ribozyme’s active site.

Rescue of dystrophin expression in mdx mouse muscle by RNA/DNA oligonucleotides

Chimeric RNA/DNA oligonucleotides (“chimeraplasts”) have been shown to induce single base alterations in genomic DNA both in vitro and in vivo. The mdx mouse strain has a point mutation in the dystrophin gene, the consequence of which is a muscular dystrophy resulting from deficiency of the dystrophin protein in skeletal muscle. To test the feasibility of chimeraplast-mediated gene therapy for muscular dystrophies, we used a chimeraplast (designated “MDX1”) designed to correct the point mutation in the dystrophin gene in mdx mice. After direct injection of MDX1 into muscles of mdx mice, immunohistochemical analysis revealed dystrophin-positive fibers clustered around the injection site. Two weeks after single injections into tibialis anterior muscles, the maximum number of dystrophin-positive fibers (approximately 30) in any muscle represented 1–2% of the total number of fibers in that muscle. Ten weeks after single injections, the range of the number of dystrophin-positive fibers was similar to that seen after 2 wk, suggesting that the expression was stable, as would be predicted for a gene-conversion event. Staining with exon-specific antibodies showed that none of these were “revertant fibers.” Furthermore, dystrophin from MDX1-injected muscles was full length by immunoblot analysis. No dystrophin was detectable by immunohistochemical or immunoblot analysis after control chimeraplast injections. Finally, reverse transcription–PCR analysis demonstrated the presence of transcripts with the wild-type dystrophin sequence only in mdx muscles injected with MDX1 chimeraplasts. These results provide the foundation for further studies of chimeraplast-mediated gene therapy as a therapeutic approach to muscular dystrophies and other genetic disorders of muscle.

Mutations in a NIMA-related kinase gene, Nek1, cause pleiotropic effects including a progressive polycystic kidney disease in mice

We previously have described a mouse model for polycystic kidney disease (PKD) caused by either of two mutations, kat or kat2J, that map to the same locus on chromosome 8. The homozygous mutant animals have a latent onset, slowly progressing form of PKD with renal pathology similar to the human autosomal-dominant PKD. In addition, the mutant animals show pleiotropic effects that include facial dysmorphism, dwarfing, male sterility, anemia, and cystic choroid plexus. We previously fine-mapped the kat2J mutation to a genetic distance of 0.28 ± 0.12 centimorgan between D8Mit128 and D8Mit129. To identify the underlying molecular defect in this locus, we constructed an integrated genetic and physical map of the critical region surrounding the kat2J mutation. Cloning and expression analysis of the transcribed sequences from this region identified Nek1, a NIMA (never in mitosis A)-related kinase as a candidate gene. Further analysis of the Nek1 gene from both kat/kat and kat2J/kat2J mutant animals identified a partial internal deletion and a single-base insertion as the molecular basis for these mutations. The complex pleiotropic phenotypes seen in the homozygous mutant animals suggest that the NEK1 protein participates in different signaling pathways to regulate diverse cellular processes. Our findings identify a previously unsuspected role for Nek1 in the kidney and open a new avenue for studying cystogenesis and identifying possible modes of therapy.

Generation of longer cDNA fragments from serial analysis of gene expression tags for gene identification

We have developed a technique called the generation of longer cDNA fragments from serial analysis of gene expression (SAGE) tags for gene identification (GLGI), to convert SAGE tags of 10 bases into their corresponding 3′ cDNA fragments covering hundred bases. A primer containing the 10-base SAGE tag is used as the sense primer, and a single base anchored oligo(dT) primer is used as an antisense primer in PCR, together with Pfu DNA polymerase. By using this approach, a cDNA fragment extending from the SAGE tag toward the 3′ end of the corresponding sequence can be generated. Application of the GLGI technique can solve two critical issues in applying the SAGE technique: one is that a longer fragment corresponding to a SAGE tag, which has no match in databases, can be generated for further studies; the other is that the specific fragment corresponding to a SAGE tag can be identified from multiple sequences that match the same SAGE tag. The development of the GLGI method provides several potential applications. First, it provides a strategy for even wider application of the SAGE technique for quantitative analysis of global gene expression. Second, a combined application of SAGE/GLGI can be used to complete the catalogue of the expressed genes in human and in other eukaryotic species. Third, it can be used to identify the 3′ cDNA sequence from any exon within a gene. It can also be used to confirm the reality of exons predicted by bioinformatic tools in genomic sequences. Fourth, a combined application of SAGE/GLGI can be applied to define the 3′ boundary of expressed genes in the genomic sequences in human and in other eukaryotic genomes.