Abstract: FUS1 is a tumor suppressor gene located on human chromosome 3p21, and expression of Fus1 protein is highly regulated at various levels, leading to lost or greatly diminished tumor suppressor function in many lung cancers. Here we show that selected microRNAs (miRNA) interact with the 3'-untranslated region (3'UTR) of FUS1, leading to down-regulation of protein expression. Using computational methods, we first predicted that FUS1 is a target of three miRNAs, miR-93, miR-98, and miR-197, and then showed that exogenous overexpression of these miRNAs inhibited Fus1 protein expression. We then confirmed that the three miRNAs target the 3'UTR region of the FUS1 transcript and that individual deletion of the three miRNA target sites in the FUS1 3'UTR restores the expression level of Fus1 protein. We further found that miR-93 and miR-98 are expressed at higher levels in small-cell lung cancer cell lines (SCLC) than in non-small-cell lung cancer cell lines (NSCLC) and immortalized human bronchial epithelial cells (HBEC), and that miR-197 is expressed at higher levels in both SCLCs and NSCLCs than in HBECs. Finally, we found that elevated miR-93 and miR-197 expression is correlated with reduced Fus1 expression in NSCLC tumor specimens. These results suggest that the three miRNAs are negative regulators of Fus1 expression in lung cancers.
Abstract: Synthetic small duplex RNAs that are complementary to gene promoters can activate or inhibit target gene expression. The potency and robustness of gene modulation by these RNAs suggests that natural mechanisms may exist to facilitate recognition of sequences within gene promoters by endogenous small RNAs. Here, we describe computational methods for identifying potential miRNA target sites within gene promoters. These methods will facilitate investigations of whether miRNAs interact with sequences outside of 3'-untranslated regions and suggest new targets for the design of synthetic modulators of gene expression.
Abstract: Nonalcoholic fatty liver disease (NAFLD) is a burgeoning health problem of unknown etiology that varies in prevalence among ancestry groups. To identify genetic variants contributing to differences in hepatic fat content, we carried out a genome-wide association scan of nonsynonymous sequence variations (n = 9,229) in a population comprising Hispanic, African American and European American individuals. An allele in PNPLA3 (rs738409[G], encoding I148M) was strongly associated with increased hepatic fat levels (P = 5.9 x 10(-10)) and with hepatic inflammation (P = 3.7 x 10(-4)). The allele was most common in Hispanics, the group most susceptible to NAFLD; hepatic fat content was more than twofold higher in PNPLA3 rs738409[G] homozygotes than in noncarriers. Resequencing revealed another allele of PNPLA3 (rs6006460[T], encoding S453I) that was associated with lower hepatic fat content in African Americans, the group at lowest risk of NAFLD. Thus, variation in PNPLA3 contributes to ancestry-related and inter-individual differences in hepatic fat content and susceptibility to NAFLD.
Abstract: Coronary heart disease (CHD) is a major cause of death in Western countries. We used genome-wide association scanning to identify a 58-kilobase interval on chromosome 9p21 that was consistently associated with CHD in six independent samples (more than 23,000 participants) from four Caucasian populations. This interval, which is located near the CDKN2A and CDKN2B genes, contains no annotated genes and is not associated with established CHD risk factors such as plasma lipoproteins, hypertension, or diabetes. Homozygotes for the risk allele make up 20 to 25% of Caucasians and have a approximately 30 to 40% increased risk of CHD.
Abstract: The identification of DNA sequence variants underlying human complex phenotypes remains a significant challenge for several reasons: individual variants can have small phenotypic effects or low population frequencies, and multiple allelic variants may act in concert to affect a trait. We evaluated the combined effect of allelic variants in seven genes involved in high-density lipoprotein (HDL) metabolism, using forward stepwise regression. Analysis of all known common single-nucleotide polymorphisms (SNPs) in the seven candidate genes revealed four variants that were associated with incremental changes in HDL cholesterol levels in three independent samples. Conversely, analysis of 660 polymorphisms in eight genes that do not appear to be involved in HDL metabolism did not identify any associations with plasma HDL-cholesterol levels. These data indicate that several common SNPs act in concert to influence plasma levels of HDL cholesterol.
Abstract: An approach to understand quantitative traits was recently proposed based on the finding that nonsynonymous (NS) sequence variants in certain genes are preferentially enriched at one extreme of the population distribution. The NS variants, although individually rare, are cumulatively frequent and influence quantitative traits, such as plasma lipoprotein levels. Here, we use the NS variant technique to demonstrate that genetic variation in NPC1L1 contributes to variability in cholesterol absorption and plasma levels of low-density lipoproteins (LDLs). The ratio of plasma campesterol (a plant sterol) to lathosterol (a cholesterol precursor) was used to estimate relative cholesterol absorption in a population-based study. Nonsynonymous sequence variations in NPC1L1 were five times more common in low absorbers (n = 26 of 256) than in high absorbers (n = 5 of 256) (P < 0.001). The rare variants identified in low absorbers were found in 6% of 1,832 African-Americans and were associated with lower plasma levels of LDL cholesterol (LDL-C) (96 +/- 36 mg/dl vs. 105 +/- 36 mg/dl; P = 0.005). These data, together with prior findings, reveal a genetic architecture for LDL-C levels that does not conform to current models for quantitative traits and indicate that a significant fraction of genetic variance in LDL-C is due to multiple alleles with modest effects that are present at low frequencies in the population.
Abstract: Selected missense mutations in the proprotein convertase subtilisin/kexin type 9 serine protease gene (PCSK9) cause autosomal dominant hypercholesterolemia, whereas nonsense mutations in the same gene are associated with low plasma levels of low-density lipoprotein cholesterol (LDL-C). Here, DNA sequencing and chip-based oligonucleotide hybridization were used to determine whether other sequence variations in PCSK9 contribute to differences in LDL-C levels. The coding regions of PCSK9 were sequenced in the blacks and whites from the Dallas Heart Study (n=3,543) who had the lowest (<5th percentile) and highest (>95th percentile) plasma levels of LDL-C. Of the 17 missense variants identified, 3 (R46L, L253F, and A443T) were significantly and reproducibly associated with lower plasma levels of LDL-C (reductions ranging from 3.5% to 30%). None of the low-LDL-C variants were associated with increased hepatic triglyceride content, as measured by proton magnetic resonance spectroscopy. This finding is most consistent with the reduction in LDL-C being caused primarily by accelerating LDL clearance, rather than by reduced lipoprotein production. Association studies with 93 noncoding single-nucleotide polymorphisms (SNPs) at the PCSK9 locus identified 3 SNPs associated with modest differences in plasma LDL-C levels. Thus, a spectrum of sequence variations ranging in frequency (from 0.2% to 34%) and magnitude of effect (from a 3% increase to a 49% decrease) contribute to interindividual differences in LDL-C levels. These findings reveal that PCSK9 activity is a major determinant of plasma levels of LDL-C in humans and make it an attractive therapeutic target for LDL-C lowering.
Abstract: Secreted and cell surface proteins play important roles in cancer and are potential drug targets and tumor markers. Here, we describe a large-scale analysis of the genes encoding secreted and cell surface proteins in breast cancer. To identify these genes, we developed a novel signal sequence trap method called Escherichia coli ampicillin secretion trap (CAST). For CAST, we constructed a plasmid in which the signal sequence of beta-lactamase was deleted such that it does not confer ampicillin resistance. Eukaryotic cDNA libraries cloned into pCAST produced tens of thousands of ampicillin-resistant clones, 80% of which contained cDNA fragments encoding secreted and membrane spanning proteins. We identified 2,708 unique sequences from cDNA libraries made from surgical breast cancer specimens. We analyzed the expression of 1,287 of the 2,708 genes and found that 166 were overexpressed in breast cancers relative to normal breast tissues. Eighty-five percent of these genes had not been previously identified as markers of breast cancer. Twenty-three of the 166 genes (14%) were relatively tissue restricted, suggesting use as cancer-specific targets. We also identified several new markers of ovarian cancer. Our results indicate that CAST is a robust, rapid, and low cost method to identify cell surface and secreted proteins and is applicable to a variety of relevant biological questions.
Abstract: The low-density lipoprotein receptor (LDLR) prevents hypercholesterolemia and atherosclerosis by removing low-density lipoprotein (LDL) from circulation. Mutations in the genes encoding either LDLR or its ligand (APOB) cause severe hypercholesterolemia. Missense mutations in PCSK9, encoding a serine protease in the secretory pathway, also cause hypercholesterolemia. These mutations are probably gain-of-function mutations, as overexpression of PCSK9 in the liver of mice produces hypercholesterolemia by reducing LDLR number. To test whether loss-of-function mutations in PCSK9 have the opposite effect, we sequenced the coding region of PCSK9 in 128 subjects (50% African American) with low plasma levels of LDL and found two nonsense mutations (Y142X and C679X). These mutations were common in African Americans (combined frequency, 2%) but rare in European Americans (<0.1%) and were associated with a 40% reduction in plasma levels of LDL cholesterol. These data indicate that common sequence variations have large effects on plasma cholesterol levels in selected populations.
Abstract: We have cloned the chicken and mouse orthologues of the Caenorhabditis elegans heterochronic gene lin-41. During limb development, lin-41 is expressed in three phases over developmental time and most notably is associated with the developing autopod. Using chicken and mouse mutants and bead implantations, we report that lin-41 is genetically and biochemically downstream of both the Shh and Fgf signaling pathways. In C. elegans, it is proposed that lin-41 activity is temporally regulated by miRNAs (let-7 and lin-4) that bind to complementary sites in the lin-41 3'-untranslated region (UTR). Taking a bioinformatics approach, we also report the presence of potential miRNA binding sites in the 3'-UTR of chicken lin-41, including sites for the chicken orthologues of both C. elegans let-7 and lin-4. Finally, we show that these miRNAs and others are expressed in the chick limb consistent with the hypothesis that they regulate chicken Lin-41 activity in vivo.
Abstract: We describe a method for the generation of knowledge-based potentials and apply it to the observed torsional angles of known protein structures. The potential is derived using Bayesian reasoning, and is useful as a prior for further such reasoning in the presence of additional data. The potential takes the form of a probability density function, which is described by a small number of coefficients with the number of necessary coefficients determined by tests based on statistical significance and entropy. We demonstrate the methods in deriving one such potential corresponding to two dimensions, the Ramachandran plot. In contrast to traditional histogram-based methods, the function is continuous and differentiable. These properties allow us to use the function as a force term in the energy minimization of appropriately described structures. The method can easily be extended to other observable angles and higher dimensions, or to include sequence dependence and should find applications in structure determination and validation.
Abstract: Susceptibility to autoimmunity in B6.Sle1b mice is associated with extensive polymorphisms between two divergent haplotypes of the SLAM/CD2 family of genes. The B6.Sle1b-derived SLAM/CD2 family haplotype is found in many other laboratory mouse strains but only causes autoimmunity in the context of the C57Bl/6 (B6) genome. Phenotypic analyses have revealed variations in the structure and expression of several members of the SLAM/CD2 family in T and B lymphocytes from B6.Sle1b mice. T lymphocytes from B6.Sle1b mice have modified signaling responses to stimulation at 4-6 weeks of age. While autoimmunity may be mediated by a combination of genes in the SLAM/CD2 family cluster, the strongest candidate is Ly108, a specific isoform of which is constitutively upregulated in B6.Sle1b lymphocytes.
Abstract: Heritable variation in complex traits is generally considered to be conferred by common DNA sequence polymorphisms. We tested whether rare DNA sequence variants collectively contribute to variation in plasma levels of high density lipoprotein cholesterol (HDL-C). We sequenced three candidate genes (ABCA1, APOA1, and LCAT) that cause Mendelian forms of low HDL-C levels in individuals from a population-based study. Nonsynonymous sequence variants were significantly more common (16% versus 2%) in individuals with low HDL-C (<fifth percentile) than in those with high HDL-C (>95th percentile). Similar findings were obtained in an independent population, and biochemical studies indicated that most sequence variants in the low HDL-C group were functionally important. Thus, rare alleles with major phenotypic effects contribute significantly to low plasma HDL-C levels in the general population.
Abstract: ELXR (Exon Locator and Extractor for Resequencing) streamlines the process of determining exon/intron boundaries and designing PCR and sequencing primers for high-throughput resequencing of exons. We have pre-computed ELXR primer sets for all exons identified from the human, mouse, and rat mRNA reference sequence (RefSeq) public databases curated by the National Center for Biotechnology Information. The resulting exon-flanking PCR primer pairs have been compiled into a system called ELXRdb, which may be searched by keyword, gene name or RefSeq accession number.
Abstract: SNPCEQer II is a graphical user interface (GUI)-based application that integrates single nucleotide polymorphism (SNP) detection, SNP analysis and SNP editing in the Microsoft Windows (R) environment. SNPCEQer II detects SNPs in DNA sequences generated by the Beckman CEQ TM 2000 XL DNA analysis system. It provides tools to analyse SNPs by inspecting and comparing trace data (chromatograms) around putative SNPs with that of other related DNA sequences, and it can search for those SNPs in the National Center for Biotechnology Information (NCBI) databases. SNPCEQer II can determine the mutation type of a coding SNP and generate data for submission to the dbSNP database. The SNP report can be edited and printed, as can the chromatograms. SNPCEQer II is implemented in Visual C++.
Abstract: SNPCEQer identifies and reports SNPs in sequences obtained from the Beckman CEQ2000 DNA Analysis System. SNPCEQer aligns sequences obtained using CEQ2000 heterozygote detection analysis and reports discrepancies between individual sequences and the consensus sequence it generates from this set as SNPs when the individual base calls have high-quality values. SNPCEQer reported comparable numbers of SNPs to the UNIX-based PolyPhred (148 vs. 165, respectively) in regions amplified from eight genes. A total of 21 different SNPs was discovered. Each gene region was analyzed in 96-306 samples. SNPCEQer was designed to operate from Windows NT, making SNP detection more accessible to users without UNIX systems. SNPCEQer is available free of charge at http://innovation.swmed.edu.
Abstract: ARROGANT (ARRay OrGANizing Tool) is a software tool developed to facilitate the identification, annotation and comparison of large collections of genes or clones. The objective is to enable users to compile gene/clone collections from different databases, allowing them to design experiments and analyze the collections as well as associated experimental data efficiently. ARROGANT can relate different sequence identifiers to their common reference sequence using the UniGene database, allowing for the comparison of data from two different microarray experiments. ARROGANT has been successfully used to analyze microarray expression data for colon cancer, to compile genes potentially related to cardiac diseases for subsequent resequencing (to identify single nucleotide polymorphisms, SNPs), to design a new comprehensive human cDNA microarray for cancer, to combine and compare expression data generated by different microarrays and to provide annotation for genes on custom and Affymetrix chips.
Abstract: Searching for similarities between biological sequences is the principal means by which bioinformatics contributes to our understanding of biology. Of the various informatics tools developed to accomplish this task, the most widely used is BLAST, the basic local alignment search tool. This article discusses the principles, workings, applications and potential pitfalls of BLAST, focusing on the implementation developed at the National Center for Biotechnology Information.
Abstract: We have developed an algorithm that predicted 11,265 potentially polymorphic tandem repeats within transcribed sequences. We estimate that 22% (2,207/9,717) of the annotated clusters within UniGene contain at least one potentially polymorphic locus. Our predictions were tested by allelotyping a panel of approximately 30 individuals for 5% of these regions, confirming polymorphism for more than half the loci tested. Our study indicates that tandem-repeat polymorphisms in genes are more common than is generally believed. Approximately 8% of these loci are within coding sequences and, if polymorphic, would result in frameshifts. Our catalogue of putative polymorphic repeats within transcribed sequences comprises a large set of potentially phenotypic or disease-causing loci. In addition, from the anomalous character of the repetitive sequences within unannotated clusters, we also conclude that the UniGene cluster count substantially overestimates the number of genes in the human genome. We hypothesize that polymorphisms in repeated sequences occur with some baseline distribution, on the basis of repeat homogeneity, size, and sequence composition, and that deviations from that distribution are indicative of the nature of selection pressure at that locus. We find evidence of selective maintenance of the ability of some genes to respond very rapidly, perhaps even on intragenerational timescales, to fluctuating selective pressures.
Abstract: As the exponential growth of DNA sequence information in databases continues, the task of converting this deposited information into knowledge becomes more dependent on integrative sequence analysis and visualization tools. PANORAMA is an Internet-accessible software package that performs a variety of informatics analyses on a given DNA sequence and returns a visual and interactive representation of the results. Its design is modular, so that further sequence analysis tools can be integrated with minimal effort. The utility of PANORAMA is demonstrated in the analysis of 650 kb of human genomic DNA from chromosome region 3p21.3, a region of potential tumor suppressor genes involved in lung cancer, breast cancer, and other forms of cancer. PANORAMA aided in the discovery of genes and alternate splice forms of known exons, in the demarcation of intron-exon boundaries, and in the identification of promoter regions and polymorphisms, all of which contributed to a better understanding of the region. PANORAMA is available on the World Wide Web at http://atlas.swmed.edu.
Abstract: We have combined neutron solution scattering experiments with molecular dynamics simulation to isolate an excess experimental signal that is caused solely by N-acetyl-leucine-amide (NALA) correlations in aqueous solution. This excess signal contains information about how NALA molecule centers are correlated in water, and we show how these solute-solute correlations might be determined at dilute concentrations in the small angle region. We have tested qualitatively different pair distribution functions for NALA molecule centers-gas, cluster, and aqueous forms of gc(r)-and have found that the excess experimental signal is adequate enough to rule out gas and cluster pair distribution functions. The aqueous form of gc(r) that exhibits a solvent-separated minimum, and possibly longer-ranged correlations as well, is not only physically sound but reproduces the experimental data reasonably well. This work demonstrates that important information in the small angle region can be mined to resolve solute-solute correlations, their lengthscales, and thermodynamic consequences even at dilute concentrations. The hydration forces that operate on the microscopic scale of individual amino acid side chains, implied by the small angle scattering data, could have significant effects on the early stages of protein folding, on ligand binding, and on other intermolecular interactions.
Abstract: We use molecular dynamics to simulate recent neutron scattering experiments on aqueous solutions of N-acetyl-leucine-amide and N-acetyl-glutamine-amide, and break down the total scattering function into contributions from solute-solute, solute-water, water-water, and intramolecular correlations. We show that the shift of the main diffraction peak to smaller angle that is observed for leucine, but not for glutamine, is attributable primarily to alterations in water-water correlations relative to bulk. The perturbation of the water hydrogen-bonded network extends roughly two solvation layers from the hydrophobic side chain surface, and is characterized by a distribution of hydrogen bonded ring sizes that are more planar and are dominated by pentagons in particular than those near the hydrophilic side chain. The different structural organization of water near the hydrophobic solute that gives rise to the inward shift in the main neutron diffraction peak under ambient conditions may also provide insight into the same directional shift for pure liquid water as it is cooled and supercooled.
Abstract: Neutron scattering experiments are used to determine scattering profiles for aqueous solutions of hydrophobic and hydrophilic amino acid analogs. Solutions of hydrophobic solutes show a shift in the main diffraction peak to smaller angle as compared with pure water, whereas solutions of hydrophilic solutes do not. The same difference for solutions of hydrophobic and hydrophilic side chains is also predicted by molecular dynamics simulations. The neutron scattering curves of aqueous solutions of hydrophobic amino acids at room temperature are qualitatively similar to differences between the liquid molecular structure functions measured for ambient and supercooled water. The nonpolar solute-induced expansion of water structure reported here is also complementary to recent neutron experiments where compression of aqueous solvent structure has been observed at high salt concentration.
