ClinVar: improvements to accessing data.

ClinVar is a freely available, public archive of human genetic variants and interpretations of their relationships to diseases and other conditions, maintained at the National Institutes of Health (NIH). Submitted interpretations of variants are aggregated and made available on the ClinVar website (https://www.ncbi.nlm.nih.gov/clinvar/), and as downloadable files via FTP and through programmatic tools such as NCBI's E-utilities. The default view on the ClinVar website, the Variation page, was recently redesigned. The new layout includes several new sections that make it easier to find submitted data as well as summary data such as all diseases and citations reported for the variant. The new design also better represents more complex data such as haplotypes and genotypes, as well as variants that are in ClinVar as part of a haplotype or genotype but have no interpretation for the single variant. ClinVar's variant-centric XML had its production release in April 2019. The ClinVar website and E-utilities both have been updated to support the VCV (variation in ClinVar) accession numbers found in the variant-centric XML file. ClinVar's search engine has been fine-tuned for improved retrieval of search results.

ClinVar: improving access to variant interpretations and supporting evidence.

ClinVar (https://www.ncbi.nlm.nih.gov/clinvar/) is a freely available, public archive of human genetic variants and interpretations of their significance to disease, maintained at the National Institutes of Health. Interpretations of the clinical significance of variants are submitted by clinical testing laboratories, research laboratories, expert panels and other groups. ClinVar aggregates data by variant-disease pairs, and by variant (or set of variants). Data aggregated by variant are accessible on the website, in an improved set of variant call format files and as a new comprehensive XML report. ClinVar recently started accepting submissions that are focused primarily on providing phenotypic information for individuals who have had genetic testing. Submissions may come from clinical providers providing their own interpretation of the variant ('provider interpretation') or from groups such as patient registries that primarily provide phenotypic information from patients ('phenotyping only'). ClinVar continues to make improvements to its search and retrieval functions. Several new fields are now indexed for more precise searching, and filters allow the user to narrow down a large set of search results.

ClinVar: public archive of interpretations of clinically relevant variants.

ClinVar (https://www.ncbi.nlm.nih.gov/clinvar/) at the National Center for Biotechnology Information (NCBI) is a freely available archive for interpretations of clinical significance of variants for reported conditions. The database includes germline and somatic variants of any size, type or genomic location. Interpretations are submitted by clinical testing laboratories, research laboratories, locus-specific databases, OMIM®, GeneReviews™, UniProt, expert panels and practice guidelines. In NCBI's Variation submission portal, submitters upload batch submissions or use the Submission Wizard for single submissions. Each submitted interpretation is assigned an accession number prefixed with SCV. ClinVar staff review validation reports with data types such as HGVS (Human Genome Variation Society) expressions; however, clinical significance is reported directly from submitters. Interpretations are aggregated by variant-condition combination and assigned an accession number prefixed with RCV. Clinical significance is calculated for the aggregate record, indicating consensus or conflict in the submitted interpretations. ClinVar uses data standards, such as HGVS nomenclature for variants and MedGen identifiers for conditions. The data are available on the web as variant-specific views; the entire data set can be downloaded via ftp. Programmatic access for ClinVar records is available through NCBI's E-utilities. Future development includes providing a variant-centric XML archive and a web page for details of SCV submissions.

Enhancement of low-energy electron emission in 2D radioactive films.

High-energy radiation has been used for decades; however, the role of low-energy electrons created during irradiation has only recently begun to be appreciated. Low-energy electrons are the most important component of radiation damage in biological environments because they have subcellular ranges, interact destructively with chemical bonds, and are the most abundant product of ionizing particles in tissue. However, methods for generating them locally without external stimulation do not exist. Here, we synthesize one-atom-thick films of the radioactive isotope (125)I on gold that are stable under ambient conditions. Scanning tunnelling microscopy, supported by electronic structure simulations, allows us to directly observe nuclear transmutation of individual (125)I atoms into (125)Te, and explain the surprising stability of the 2D film as it underwent radioactive decay. The metal interface geometry induces a 600% amplification of low-energy electron emission (<10 eV; ref. ) compared with atomic (125)I. This enhancement of biologically active low-energy electrons might offer a new direction for highly targeted nanoparticle therapies.

Pulse energy measurement at the SXR instrument.

A gas monitor detector was implemented and characterized at the Soft X-ray Research (SXR) instrument to measure the average, absolute and pulse-resolved photon flux of the LCLS beam in the energy range between 280 and 2000 eV. The detector is placed after the monochromator and addresses the need to provide reliable absolute pulse energy as well as pulse-resolved measurements for the various experiments at this instrument. This detector provides a reliable non-invasive measurement for determining flux levels on the samples in the downstream experimental chamber and for optimizing signal levels of secondary detectors and for the essential need of data normalization. The design, integration into the instrument and operation are described, and examples of its performance are given.

Gene: a gene-centered information resource at NCBI.

The National Center for Biotechnology Information's (NCBI) Gene database (www.ncbi.nlm.nih.gov/gene) integrates gene-specific information from multiple data sources. NCBI Reference Sequence (RefSeq) genomes for viruses, prokaryotes and eukaryotes are the primary foundation for Gene records in that they form the critical association between sequence and a tracked gene upon which additional functional and descriptive content is anchored. Additional content is integrated based on the genomic location and RefSeq transcript and protein sequence data. The content of a Gene record represents the integration of curation and automated processing from RefSeq, collaborating model organism databases, consortia such as Gene Ontology, and other databases within NCBI. Records in Gene are assigned unique, tracked integers as identifiers. The content (citations, nomenclature, genomic location, gene products and their attributes, phenotypes, sequences, interactions, variation details, maps, expression, homologs, protein domains and external databases) is available via interactive browsing through NCBI's Entrez system, via NCBI's Entrez programming utilities (E-Utilities and Entrez Direct) and for bulk transfer by FTP.

RefSeq: an update on mammalian reference sequences.

The National Center for Biotechnology Information (NCBI) Reference Sequence (RefSeq) database is a collection of annotated genomic, transcript and protein sequence records derived from data in public sequence archives and from computation, curation and collaboration (http://www.ncbi.nlm.nih.gov/refseq/). We report here on growth of the mammalian and human subsets, changes to NCBI's eukaryotic annotation pipeline and modifications affecting transcript and protein records. Recent changes to NCBI's eukaryotic genome annotation pipeline provide higher throughput, and the addition of RNAseq data to the pipeline results in a significant expansion of the number of transcripts and novel exons annotated on mammalian RefSeq genomes. Recent annotation changes include reporting supporting evidence for transcript records, modification of exon feature annotation and the addition of a structured report of gene and sequence attributes of biological interest. We also describe a revised protein annotation policy for alternatively spliced transcripts with more divergent predicted proteins and we summarize the current status of the RefSeqGene project.

Current status and new features of the Consensus Coding Sequence database.

The Consensus Coding Sequence (CCDS) project (http://www.ncbi.nlm.nih.gov/CCDS/) is a collaborative effort to maintain a dataset of protein-coding regions that are identically annotated on the human and mouse reference genome assemblies by the National Center for Biotechnology Information (NCBI) and Ensembl genome annotation pipelines. Identical annotations that pass quality assurance tests are tracked with a stable identifier (CCDS ID). Members of the collaboration, who are from NCBI, the Wellcome Trust Sanger Institute and the University of California Santa Cruz, provide coordinated and continuous review of the dataset to ensure high-quality CCDS representations. We describe here the current status and recent growth in the CCDS dataset, as well as recent changes to the CCDS web and FTP sites. These changes include more explicit reporting about the NCBI and Ensembl annotation releases being compared, new search and display options, the addition of biologically descriptive information and our approach to representing genes for which support evidence is incomplete. We also present a summary of recent and future curation targets.

NCBI Reference Sequences (RefSeq): current status, new features and genome annotation policy.

The National Center for Biotechnology Information (NCBI) Reference Sequence (RefSeq) database is a collection of genomic, transcript and protein sequence records. These records are selected and curated from public sequence archives and represent a significant reduction in redundancy compared to the volume of data archived by the International Nucleotide Sequence Database Collaboration. The database includes over 16,00 organisms, 2.4 × 0(6) genomic records, 13 × 10(6) proteins and 2 × 10(6) RNA records spanning prokaryotes, eukaryotes and viruses (RefSeq release 49, September 2011). The RefSeq database is maintained by a combined approach of automated analyses, collaboration and manual curation to generate an up-to-date representation of the sequence, its features, names and cross-links to related sources of information. We report here on recent growth, the status of curating the human RefSeq data set, more extensive feature annotation and current policy for eukaryotic genome annotation via the NCBI annotation pipeline. More information about the resource is available online (see http://www.ncbi.nlm.nih.gov/RefSeq/).

Quantification of recombinant and platelet P2Y(1) receptors utilizing a [(125)I]-labeled high-affinity antagonist 2-iodo-N(6)-methyl-(N)-methanocarba-2'-deoxyadenosine-3',5'-bisphosphate ([(125)I]MRS2500).

The ADP-activated P2Y(1) receptor is broadly expressed and plays a crucial role in ADP-promoted platelet aggregation. We previously synthesized 2-iodo-N(6)-methyl-(N)-methanocarba-2'-deoxyadenosine-3',5'-bisphosphate (MRS2500), as a selective, high-affinity, competitive antagonist of this receptor. Here we report utilization of a trimethylstannyl precursor molecule for the multi-step radiochemical synthesis of a [(125)I]-labeled form of MRS2500. [(125)I]MRS2500 bound selectively to Sf9 insect cell membranes expressing the human P2Y(1) receptor but did not specifically bind to membranes isolated from empty vector-infected cells. Binding of [(125)I]MRS2500 to P2Y(1) receptors was saturable with a Kd of 1.2nM. Known agonists and antagonists of the P2Y(1) receptor inhibited [(125)I]MRS2500 binding to P2Y(1) receptor-expressing membranes with potencies in agreement with those previously observed in functional assays of this receptor. A high-affinity binding site for [(125)I]MRS2500 also was observed on intact human platelets (Kd=0.61nM) and mouse platelets (Kd=1.20nM) that exhibited the pharmacological selectivity of the P2Y(1) receptor. The densities of sites observed were 151 sites/platelet and 229 sites/platelet in human and mouse platelets, respectively. In contrast, specific binding was not observed in platelets isolated from P2Y(1) receptor(-/-) mice. Taken together, these data illustrate the synthesis and characterization of a novel P2Y(1) receptor radioligand and its utility for examining P2Y(1) receptors natively expressed on human and mouse platelets.

Functionalized congeners of P2Y1 receptor antagonists: 2-alkynyl (N)-methanocarba 2'-deoxyadenosine 3',5'-bisphosphate analogues and conjugation to a polyamidoamine (PAMAM) dendrimer carrier.

The P2Y(1) receptor is a prothrombotic G protein-coupled receptor (GPCR) activated by ADP. Preference for the North (N) ring conformation of the ribose moiety of adenine nucleotide 3',5'-bisphosphate antagonists of the P2Y(1) receptor was established by using a ring-constrained methanocarba (a bicyclo[3.1.0]hexane) ring as a ribose substitute. A series of covalently linkable N(6)-methyl-(N)-methanocarba-2'-deoxyadenosine-3',5'-bisphosphates containing extended 2-alkynyl chains was designed, and binding affinity at the human (h) P2Y(1) receptor determined. The chain of these functionalized congeners contained hydrophilic moieties, a reactive substituent, or biotin, linked via an amide. Variation of the chain length and position of an intermediate amide group revealed high affinity of carboxylic congener 8 (K(i) 23 nM) and extended amine congener 15 (K(i) 132 nM), both having a 2-(1-pentynoyl) group. A biotin conjugate 18 containing an extended epsilon-aminocaproyl spacer chain exhibited higher affinity than a shorter biotinylated analogue. Alternatively, click coupling of terminal alkynes of homologous 2-dialkynyl nucleotide derivatives to alkyl azido groups produced triazole derivatives that bound to the P2Y(1) receptor following deprotection of the bisphosphate groups. The preservation of receptor affinity of the functionalized congeners was consistent with new P2Y(1) receptor modeling and ligand docking. Attempted P2Y(1) antagonist conjugation to PAMAM dendrimer carriers by amide formation or palladium-catalyzed reaction between an alkyne on the dendrimer and a 2-iodopurine-derivatized nucleotide was unsuccessful. A dialkynyl intermediate containing the chain length favored in receptor binding was conjugated to an azide-derivatized dendrimer, and the conjugate inhibited ADP-promoted human platelet aggregation. This is the first example of attaching a strategically functionalized P2Y receptor antagonist to a PAMAM dendrimer to produce a multivalent conjugate exhibiting a desired biological effect, i.e., antithrombotic action.

Synthesis and pharmacological characterization of [(125)I]MRS5127, a high affinity, selective agonist radioligand for the A3 adenosine receptor.

A recently reported selective agonist of the human A(3) adenosine receptor (hA(3)AR), MRS5127 (1'R,2'R,3'S,4'R,5'S)-4'-[2-chloro-6-(3-iodobenzylamino)-purine]-2',3'-O-dihydroxy-bicyclo-[3.1.0]hexane, was radioiodinated and characterized pharmacologically. It contains a rigid bicyclic ring system in place of a 5'-truncated ribose moiety, and was selected for radiolabeling due to its nanomolar binding affinity at both human and rat A(3)ARs. The radioiodination of the N(6)-3-iodobenzyl substituent by iododestannylation of a 3-(trimethylstannyl)benzyl precursor was achieved in 73% yield, measured after purification by HPLC. [(125)I]MRS5127 bound to the human A(3)AR expressed in membranes of stably transfected HEK 293 cells. Specific binding was saturable, competitive, and followed a one-site binding model, with a K(d) value of 5.74+/-0.97nM. At a concentration equivalent to its K(d), non-specific binding comprised 27+/-2% of total binding. In kinetic studies, [(125)I]MRS5127 rapidly associated with the hA(3)AR (t(1/2)=0.514+/-0.014min), and the affinity calculated from association and dissociation rate constants was 3.50+/-1.46nM. The pharmacological profile of ligands in competition experiments with [(125)I]MRS5127 was consistent with the known structure-activity-relationship profile of the hA(3)AR. [(125)I]MRS5127 bound with similar high affinity (K(d), nM) to recombinant A(3)ARs from mouse (4.90+/-0.77), rabbit (2.53+/-0.11), and dog (3.35+/-0.54). For all of the species tested, MRS5127 exhibited A(3)AR agonist activity based on negative coupling to cAMP production. Thus, [(125)I]MRS5127 represents a new species-independent agonist radioligand for the A(3)AR. The major advantage of [(125)I]MRS5127 compared with previously used A(3)AR radioligands is its high affinity, low degree of non-specific binding, and improved A(3)AR selectivity.

The completion of the Mammalian Gene Collection (MGC).

Since its start, the Mammalian Gene Collection (MGC) has sought to provide at least one full-protein-coding sequence cDNA clone for every human and mouse gene with a RefSeq transcript, and at least 6200 rat genes. The MGC cloning effort initially relied on random expressed sequence tag screening of cDNA libraries. Here, we summarize our recent progress using directed RT-PCR cloning and DNA synthesis. The MGC now contains clones with the entire protein-coding sequence for 92% of human and 89% of mouse genes with curated RefSeq (NM-accession) transcripts, and for 97% of human and 96% of mouse genes with curated RefSeq transcripts that have one or more PubMed publications, in addition to clones for more than 6300 rat genes. These high-quality MGC clones and their sequences are accessible without restriction to researchers worldwide.

The consensus coding sequence (CCDS) project: Identifying a common protein-coding gene set for the human and mouse genomes.

Effective use of the human and mouse genomes requires reliable identification of genes and their products. Although multiple public resources provide annotation, different methods are used that can result in similar but not identical representation of genes, transcripts, and proteins. The collaborative consensus coding sequence (CCDS) project tracks identical protein annotations on the reference mouse and human genomes with a stable identifier (CCDS ID), and ensures that they are consistently represented on the NCBI, Ensembl, and UCSC Genome Browsers. Importantly, the project coordinates on manually reviewing inconsistent protein annotations between sites, as well as annotations for which new evidence suggests a revision is needed, to progressively converge on a complete protein-coding set for the human and mouse reference genomes, while maintaining a high standard of reliability and biological accuracy. To date, the project has identified 20,159 human and 17,707 mouse consensus coding regions from 17,052 human and 16,893 mouse genes. Three evaluation methods indicate that the entries in the CCDS set are highly likely to represent real proteins, more so than annotations from contributing groups not included in CCDS. The CCDS database thus centralizes the function of identifying well-supported, identically-annotated, protein-coding regions.

What everybody should know about the rat genome and its online resources.

It has been four years since the original publication of the draft sequence of the rat genome. Five groups are now working together to assemble, annotate and release an updated version of the rat genome. As the prevailing model for physiology, complex disease and pharmacological studies, there is an acute need for the rat's genomic resources to keep pace with the rat's prominence in the laboratory. In this commentary, we describe the current status of the rat genome sequence and the plans for its impending 'upgrade'. We then cover the key online resources providing access to the rat genome, including the new SNP views at Ensembl, the RefSeq and Genes databases at the US National Center for Biotechnology Information, Genome Browser at the University of California Santa Cruz and the disease portals for cardiovascular disease and obesity at the Rat Genome Database.

DNA sequence variation and selection of tag single-nucleotide polymorphisms at candidate genes for drought-stress response in Pinus taeda L.

Genetic association studies are rapidly becoming the experimental approach of choice to dissect complex traits, including tolerance to drought stress, which is the most common cause of mortality and yield losses in forest trees. Optimization of association mapping requires knowledge of the patterns of nucleotide diversity and linkage disequilibrium and the selection of suitable polymorphisms for genotyping. Moreover, standard neutrality tests applied to DNA sequence variation data can be used to select candidate genes or amino acid sites that are putatively under selection for association mapping. In this article, we study the pattern of polymorphism of 18 candidate genes for drought-stress response in Pinus taeda L., an important tree crop. Data analyses based on a set of 21 putatively neutral nuclear microsatellites did not show population genetic structure or genomewide departures from neutrality. Candidate genes had moderate average nucleotide diversity at silent sites (pi(sil) = 0.00853), varying 100-fold among single genes. The level of within-gene LD was low, with an average pairwise r2 of 0.30, decaying rapidly from approximately 0.50 to approximately 0.20 at 800 bp. No apparent LD among genes was found. A selective sweep may have occurred at the early-response-to-drought-3 (erd3) gene, although population expansion can also explain our results and evidence for selection was not conclusive. One other gene, ccoaomt-1, a methylating enzyme involved in lignification, showed dimorphism (i.e., two highly divergent haplotype lineages at equal frequency), which is commonly associated with the long-term action of balancing selection. Finally, a set of haplotype-tagging SNPs (htSNPs) was selected. Using htSNPs, a reduction of genotyping effort of approximately 30-40%, while sampling most common allelic variants, can be gained in our ongoing association studies for drought tolerance in pine.

Nucleotide diversity and linkage disequilibrium in loblolly pine.

Outbreeding species with large, stable population sizes, such as widely distributed conifers, are expected to harbor relatively more DNA sequence polymorphism. Under the neutral theory of molecular evolution, the expected heterozygosity is a function of the product 4N(e)mu, where N(e) is the effective population size and mu is the per-generation mutation rate, and the genomic scale of linkage disequilibrium is determined by 4N(e)r, where r is the per-generation recombination rate between adjacent sites. These parameters were estimated in the long-lived, outcrossing gymnosperm loblolly pine (Pinus taeda L.) from a survey of single nucleotide polymorphisms across approximately 18 kb of DNA distributed among 19 loci from a common set of 32 haploid genomes. Estimates of 4N(e)mu at silent and nonsynonymous sites were 0.00658 and 0.00108, respectively, and both were statistically heterogeneous among loci. By Tajima's D statistic, the site frequency spectrum of no locus was observed to deviate from that predicted by neutral theory. Substantial recombination in the history of the sampled alleles was observed and linkage disequilibrium declined within several kilobases. The composite likelihood estimate of 4N(e)r based on all two-site sample configurations equaled 0.00175. When geological dating, an assumed generation time (25 years), and an estimated divergence from Pinus pinaster Ait. are used, the effective population size of loblolly pine should be 5.6 x 10(5). The emerging narrow range of estimated silent site heterozygosities (relative to the vast range of population sizes) for humans, Drosophila, maize, and pine parallels the paradox described earlier for allozyme polymorphism and challenges simple equilibrium models of molecular evolution.

Comparative mapping in the Pinaceae.

A comparative genetic map was constructed between two important genera of the family Pinaceae. Ten homologous linkage groups in loblolly pine (Pinus taeda L.) and Douglas fir (Pseudotsuga menziesii [Mirb.] Franco) were identified using orthologous expressed sequence tag polymorphism (ESTP) and restriction fragment length polymorphism (RFLP) markers. The comparative mapping revealed extensive synteny and colinearity between genomes of the Pinaceae, consistent with the hypothesis of conservative chromosomal evolution in this important plant family. This study reports the first comparative map in forest trees at the family taxonomic level and establishes a framework for comparative genomics in Pinaceae.

Identification of quantitative trait loci influencing wood property traits in loblolly pine (Pinus taeda L.). III. QTL Verification and candidate gene mapping.

A long-term series of experiments to map QTL influencing wood property traits in loblolly pine has been completed. These experiments were designed to identify and subsequently verify QTL in multiple genetic backgrounds, environments, and growing seasons. Verification of QTL is necessary to substantiate a biological basis for observed marker-trait associations, to provide precise estimates of the magnitude of QTL effects, and to predict QTL expression at a given age or in a particular environment. Verification was based on the repeated detection of QTL among populations, as well as among multiple growing seasons for each population. Temporal stability of QTL was moderate, with approximately half being detected in multiple seasons. Fewer QTL were common to different populations, but the results are nonetheless encouraging for restricted applications of marker-assisted selection. QTL from larger populations accounted for less phenotypic variation than QTL detected in smaller populations, emphasizing the need for experiments employing much larger families. Additionally, 18 candidate genes related to lignin biosynthesis and cell wall structure were mapped genetically. Several candidate genes colocated with wood property QTL; however, these relationships must be verified in future experiments.

An arabinogalactan protein associated with secondary cell wall formation in differentiating xylem of loblolly pine.

Arabinogalactan proteins (AGPs) are abundant plant proteoglycans implicated in plant growth and development. Here, we report the genetic characterization, partial purification and immunolocalization of a classical AGP (PtaAGP6, accession number AF101785) in loblolly pine (Pinus taeda L.). A PtaAGP6 full-length cDNA clone was expressed in bacteria. PtaAGP6 resembles tomato LeAGP-1 and Arabidopsis AtAGP17-19 in that they all possess a subdomain composed of basic amino acids. The accessibility of this domain in the glycoprotein makes it possible to label the PtaAGP6 epitopes on the cell surface or in the cell wall with polyclonal antibodies raised against this subdomain. The antibodies recognize the peptide of the basic subdomain and bind to the intact protein molecule. A soluble protein-containing fraction was purified from the differentiating xylem of pine trees by using beta-glucosyl Yariv reagent (beta-glcY) and was recognized by antibodies against the basic subdomain. Immunolocalization studies showed that the PtaAGP6 epitopes are restricted to a file of cells that just precede secondary cell wall thickening, suggesting roles in xylem differentiation and wood formation. The location of apparent labeling of the PtaAGP6 epitopes is separated from the location of lignin deposition. Multiple single nucleotide polymorphisms (SNPs) were detected in EST variants. Denaturing HPLC analysis of PCR products suggests that PtaAGP6 is encoded by a single gene. Mobility variation in denaturing gel electrophoresis was used to map PtaAGP6 SNPs to a site on linkage group 5.