Database resources of the National Center for Biotechnology Information

The National Center for Biotechnology Information (NCBI) provides online information resources for biology, including the GenBank® nucleic acid sequence database and the PubMed® database of citations and abstracts published in life science journals. NCBI provides search and retrieval operations for most of these data from 35 distinct databases. The E-utilities serve as the programming interface for most of these databases. Resources receiving significant updates in the past year include PubMed, PMC, Bookshelf, SciENcv, the NIH Comparative Genomics Resource (CGR), NCBI Virus, SRA, RefSeq, foreign contamination screening tools, Taxonomy, iCn3D, ClinVar, GTR, MedGen, dbSNP, ALFA, ClinicalTrials.gov, Pathogen Detection, antimicrobial resistance resources, and PubChem. These resources can be accessed through the NCBI home page at https://www.ncbi.nlm.nih.gov.

The status of the human gene catalogue

Scientists have been trying to identify every gene in the human genome since the initial draft was published in 2001. In the years since, much progress has been made in identifying protein-coding genes, currently estimated to number fewer than 20,000, with an ever-expanding number of distinct protein-coding isoforms. Here we review the status of the human gene catalogue and the efforts to complete it in recent years. Beside the ongoing annotation of protein-coding genes, their isoforms and pseudogenes, the invention of high-throughput RNA sequencing and other technological breakthroughs have led to a rapid growth in the number of reported non-coding RNA genes. For most of these non-coding RNAs, the functional relevance is currently unclear; we look at recent advances that offer paths forward to identifying their functions and towards eventually completing the human gene catalogue. Finally, we examine the need for a universal annotation standard that includes all medically significant genes and maintains their relationships with different reference genomes for the use of the human gene catalogue in clinical settings.

The status of the human gene catalogue

Scientists have been trying to identify all of the genes in the human genome since the initial draft of the genome was published in 2001. Over the intervening years, much progress has been made in identifying protein-coding genes, and the estimated number has shrunk to fewer than 20,000, although the number of distinct protein-coding isoforms has expanded dramatically. The invention of high-throughput RNA sequencing and other technological breakthroughs have led to an explosion in the number of reported non-coding RNA genes, although most of them do not yet have any known function. A combination of recent advances offers a path forward to identifying these functions and towards eventually completing the human gene catalogue. However, much work remains to be done before we have a universal annotation standard that includes all medically significant genes, maintains their relationships with different reference genomes, and describes clinically relevant genetic variants.

Database resources of the National Center for Biotechnology Information in 2023

The National Center for Biotechnology Information (NCBI) provides online information resources for biology, including the GenBank® nucleic acid sequence database and the PubMed® database of citations and abstracts published in life science journals. NCBI provides search and retrieval operations for most of these data from 35 distinct databases. The E-utilities serve as the programming interface for most of these databases. New resources include the Comparative Genome Resource (CGR) and the BLAST ClusteredNR database. Resources receiving significant updates in the past year include PubMed, PMC, Bookshelf, IgBLAST, GDV, RefSeq, NCBI Virus, GenBank type assemblies, iCn3D, ClinVar, GTR, dbGaP, ALFA, ClinicalTrials.gov, Pathogen Detection, antimicrobial resistance resources, and PubChem. These resources can be accessed through the NCBI home page at https://www.ncbi.nlm.nih.gov.

A joint NCBI and EMBL-EBI transcript set for clinical genomics and research

Comprehensive genome annotation is essential to understand the impact of clinically relevant variants. However, the absence of a standard for clinical reporting and browser display complicates the process of consistent interpretation and reporting. To address these challenges, Ensembl/GENCODE1 and RefSeq2 launched a joint initiative, the Matched Annotation from NCBI and EMBL-EBI (MANE) collaboration, to converge on human gene and transcript annotation and to jointly define a high-value set of transcripts and corresponding proteins. Here, we describe the MANE transcript sets for use as universal standards for variant reporting and browser display. The MANE Select set identifies a representative transcript for each human protein-coding gene, whereas the MANE Plus Clinical set provides additional transcripts at loci where the Select transcripts alone are not sufficient to report all currently known clinical variants. Each MANE transcript represents an exact match between the exonic sequences of an Ensembl/GENCODE transcript and its counterpart in RefSeq such that the identifiers can be used synonymously. We have now released MANE Select transcripts for 97% of human protein-coding genes, including all American College of Medical Genetics and Genomics Secondary Findings list v3.0 (ref. 3) genes. MANE transcripts are accessible from major genome browsers and key resources. Widespread adoption of these transcript sets will increase the consistency of reporting, facilitate the exchange of data regardless of the annotation source and help to streamline clinical interpretation.

Consensus coding sequence (CCDS) database: a standardized set of human and mouse protein-coding regions supported by expert curation

The Consensus Coding Sequence (CCDS) project provides a dataset of protein-coding regions that are identically annotated on the human and mouse reference genome assembly in genome annotations produced independently by NCBI and the Ensembl group at EMBL-EBI. This dataset is the product of an international collaboration that includes NCBI, Ensembl, HUGO Gene Nomenclature Committee, Mouse Genome Informatics and University of California, Santa Cruz. Identically annotated coding regions, which are generated using an automated pipeline and pass multiple quality assurance checks, are assigned a stable and tracked identifier (CCDS ID). Additionally, coordinated manual review by expert curators from the CCDS collaboration helps in maintaining the integrity and high quality of the dataset. The CCDS data are available through an interactive web page (https://www.ncbi.nlm.nih.gov/CCDS/CcdsBrowse.cgi) and an FTP site (ftp://ftp.ncbi.nlm.nih.gov/pub/CCDS/). In this paper, we outline the ongoing work, growth and stability of the CCDS dataset and provide updates on new collaboration members and new features added to the CCDS user interface. We also present expert curation scenarios, with specific examples highlighting the importance of an accurate reference genome assembly and the crucial role played by input from the research community.

XX Disorder of Sex Development is associated with an insertion on chromosome 9 and downregulation of RSPO1 in dogs (Canis lupus familiaris)

Remarkable progress has been achieved in understanding the mechanisms controlling sex determination, yet the cause for many Disorders of Sex Development (DSD) remains unknown. Of particular interest is a rare XX DSD subtype in which individuals are negative for SRY, the testis determining factor on the Y chromosome, yet develop testes or ovotestes, and both of these phenotypes occur in the same family. This is a naturally occurring disorder in humans (Homo sapiens) and dogs (C. familiaris). Phenotypes in the canine XX DSD model are strikingly similar to those of the human XX DSD subtype. The purposes of this study were to identify 1) a variant associated with XX DSD in the canine model and 2) gene expression alterations in canine embryonic gonads that could be informative to causation. Using a genome wide association study (GWAS) and whole genome sequencing (WGS), we identified a variant on C. familiaris autosome 9 (CFA9) that is associated with XX DSD in the canine model and in affected purebred dogs. This is the first marker identified for inherited canine XX DSD. It lies upstream of SOX9 within the canine ortholog for the human disorder, which resides on 17q24. Inheritance of this variant indicates that XX DSD is a complex trait in which breed genetic background affects penetrance. Furthermore, the homozygous variant genotype is associated with embryonic lethality in at least one breed. Our analysis of gene expression studies (RNA-seq and PRO-seq) in embryonic gonads at risk of XX DSD from the canine model identified significant RSPO1 downregulation in comparison to XX controls, without significant upregulation of SOX9 or other known testis pathway genes. Based on these data, a novel mechanism is proposed in which molecular lesions acting upstream of RSPO1 induce epigenomic gonadal mosaicism.

A curated catalog of canine and equine keratin genes

Keratins represent a large protein family with essential structural and functional roles in epithelial cells of skin, hair follicles, and other organs. During evolution the genes encoding keratins have undergone multiple rounds of duplication and humans have two clusters with a total of 55 functional keratin genes in their genomes. Due to the high similarity between different keratin paralogs and species-specific differences in gene content, the currently available keratin gene annotation in species with draft genome assemblies such as dog and horse is still imperfect. We compared the National Center for Biotechnology Information (NCBI) (dog annotation release 103, horse annotation release 101) and Ensembl (release 87) gene predictions for the canine and equine keratin gene clusters to RNA-seq data that were generated from adult skin of five dogs and two horses and from adult hair follicle tissue of one dog. Taking into consideration the knowledge on the conserved exon/intron structure of keratin genes, we annotated 61 putatively functional keratin genes in both the dog and horse, respectively. Subsequently, curators in the RefSeq group at NCBI reviewed their annotation of keratin genes in the dog and horse genomes (Annotation Release 104 and Annotation Release 102, respectively) and updated annotation and gene nomenclature of several keratin genes. The updates are now available in the NCBI Gene database (https://www.ncbi.nlm.nih.gov/gene).

Reference sequence (RefSeq) database at NCBI: current status, taxonomic expansion, and functional annotation

The RefSeq project at the National Center for Biotechnology Information (NCBI) maintains and curates a publicly available database of annotated genomic, transcript, and protein sequence records (http://www.ncbi.nlm.nih.gov/refseq/). The RefSeq project leverages the data submitted to the International Nucleotide Sequence Database Collaboration (INSDC) against a combination of computation, manual curation, and collaboration to produce a standard set of stable, non-redundant reference sequences. The RefSeq project augments these reference sequences with current knowledge including publications, functional features and informative nomenclature. The database currently represents sequences from more than 55,000 organisms (>4800 viruses, >40,000 prokaryotes and >10,000 eukaryotes; RefSeq release 71), ranging from a single record to complete genomes. This paper summarizes the current status of the viral, prokaryotic, and eukaryotic branches of the RefSeq project, reports on improvements to data access and details efforts to further expand the taxonomic representation of the collection. We also highlight diverse functional curation initiatives that support multiple uses of RefSeq data including taxonomic validation, genome annotation, comparative genomics, and clinical testing. We summarize our approach to utilizing available RNA-Seq and other data types in our manual curation process for vertebrate, plant, and other species, and describe a new direction for prokaryotic genomes and protein name management.

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.

Sequence variations in equine candidate genes For XX and XY inherited disorders of sexual development

Inherited disorders of sexual development (DSD) cause sterility and infertility in horses. Mutations causing such disorders have been identified in other mammals, but there is little information on the molecular causes in horses. While the equine genome sequence has made it possible to identify candidate genes, additional tools are needed to routinely screen them for causative mutations. In this study, we designed a screening panel of polymerase chain reaction primer pairs for 15 equine genes. These are the candidate genes for testicular or ovotesticular XX DSD and XY DSD, the latter of which includes gonadal dysgenesis, androgen insensitivity syndrome (AIS), persistent Mullerian duct syndrome and isolated cryptorchidism. Six horses with testicular or ovotesticular XX DSD and controls were screened. In addition, candidate genes for androgen insensitivity syndrome, persistent Mullerian duct syndrome and isolated cryptorchidism were screened in normal horses. While no sequence variants were uniquely associated with XX DSD, the 38 sequence variants identified can serve as intragenic markers in genome-wide association studies or linkage studies to hasten mutation identification in equine XX DSD and XY DSD.

A molecular diagnostic test for persistent Müllerian duct syndrome in miniature schnauzer dogs

In persistent Müllerian duct syndrome (PMDS), Müllerian ducts fail to regress in males during sexual differentiation. In the canine miniature schnauzer model, PMDS is caused by a C to T transition in exon 3 of the Müllerian inhibiting substance type II receptor (MISRII), which introduces a DdeI restriction site. Here we report a molecular diagnostic test for PMDS in the miniature schnauzer to identify affected dogs and carriers. As our test results suggest that the mutation is identical by descent in affected dogs of this breed, the test could be used to eliminate this mutation from the miniature schnauzer breed worldwide.

A single base pair mutation encoding a premature stop codon in the MIS type II receptor is responsible for canine persistent Müllerian duct syndrome

Müllerian inhibiting substance (MIS), a secreted glycoprotein in the transforming growth factor-beta family of growth factors, mediates regression of the Müllerian ducts during embryonic sex differentiation in males. In persistent Müllerian duct syndrome (PMDS), rather than undergoing involution, the Müllerian ducts persist in males, giving rise to the uterus, fallopian tubes, and upper vagina. Genetic defects in MIS or its receptor (MISRII) have been identified in patients with PMDS. The phenotype in the canine model of PMDS derived from the miniature schnauzer breed is strikingly similar to that of human patients. In this model, PMDS is inherited as a sex-limited autosomal recessive trait. Previous studies indicated that a defect in the MIS receptor or its downstream signaling pathway was likely to be causative of the canine syndrome. In this study, the canine PMDS phenotype and clinical sequelae are described in detail. Affected and unaffected members of this pedigree are genotyped, identifying a single base pair substitution in MISRII that introduces a stop codon in exon 3. The homozygous mutation terminates translation at 80 amino acids, eliminating much of the extracellular domain and the entire transmembrane and intracellular signaling domains. Findings in this model could enable insights to be garnered from correlation of detailed clinical descriptions with molecular defects, which are not otherwise possible in the human syndrome.

Mutations in the RSPO1 coding region are not the main cause of canine SRY-negative XX sex reversal in several breeds

This report details a case of SRY-negative XX sex reversal in a mixed breed dog and surveys affected dogs of several breeds for mutations in RSPO1 coding regions. Genomic DNA from the mixed breed case was evaluated for mutations in candidate genes. Sequencing identified a homozygous G to A transition in RSPO1 exon 4 that changes a highly conserved amino acid codon in the thrombospondin domain. The possibility that this was a single nucleotide polymorphism (SNP) could not be excluded by genotyping family members. Therefore, the coding region of RSPO1 was sequenced in a survey of affected dogs, which identified a T to C transition (exon 3) in some, the above G to A transition (exon 4) in others, and no change in the remaining affected dogs. Genotypes at these base pair positions were not uniquely associated with the affected phenotype in any breed, indicating the identified transitions are most likely SNPs, not causative mutations for this canine disorder. However, the possibility that polymorphisms play a modifier role, such as changing threshold or severity of phenotypic expression in a mixed breed dog, cannot be excluded. This study emphasizes the importance of canine pedigree, breed, and population studies in evaluating candidate mutations.

Linkage to CFA29 Detected in a Genome-Wide Linkage Screen of a Canine Pedigree Segregating Sry-Negative XX Sex Reversal

Canine Sry-negative XX sex reversal is a disorder of gonadal development wherein individuals having a female karyotype develop testes or ovotestes. In this study, linkage mapping was undertaken in a pedigree derived from one proven carrier American cocker spaniel founder male and beagle females. All affected dogs in the analysis were XX true hermaphrodites and confirmed to be Sry negative by polymerase chain reaction. A genome-wide linkage screen conducted using 245 microsatellite markers revealed highest LOD score of 3.4 (marker CPH9) on CFA29. Fine mapping with additional microsatellites in the region containing CPH9 localized the Sry-negative XX sex reversal locus to a 5.4-Mb candidate region between markers CPH9 and FH3003 (LOD score 3.15). Insignificant LOD scores were found at genome-wide screen or fine mapping markers that were within 10 Mb of 45 potential candidate genes reported to have a role in mammalian sex determination or differentiation. Together, these results suggest that a novel locus on CFA29 may be responsible for sex reversal in this pedigree.

Exclusion of Candidate Genes for Canine SRY-Negative XX Sex Reversal

In mammals, the Y-linked SRY gene is normally responsible for testis induction, yet testis development can occur in the absence of Y-linked genes, including SRY. The canine model of SRY-negative XX sex reversal could lead to the discovery of novel genes in the mammalian sex determination pathway. The autosomal genes causing testis induction in this disorder in dogs, humans, pigs, and horses are presently unknown. In goats, a large deletion is responsible for sex reversal linked to the polled (hornless) phenotype. However, this region has been excluded as being causative of the canine disorder, as have WT1 and DMRT1 in more recent studies. The purpose of this study was to determine whether microsatellite marker alleles near or within five candidate genes (GATA4, FOG2, LHX1, SF1, SOX9) are associated with the affected phenotype in a pedigree of canine SRY-negative XX sex reversal. Primer sequences flanking nucleotide repeats were designed within genomic sequences of canine candidate gene homologues. Fluorescence-labeled polymorphic markers were used to screen a subset of the multigenerational pedigree, and marker alleles were determined by software. Our results indicate that the mutation causing canine SRY-negative XX sex reversal in this pedigree is unlikely to be located in regions containing these candidates.

Exclusion of Lhx9 as a Candidate Gene for Sry-Negative XX Sex Reversal in the American Cocker Spaniel Model

XX sex reversal is known in 17 breeds of dogs. In the American cocker spaniel, it segregates as an autosomal recessive trait, and the affected animals lack the testis determining Sry gene. In the search for an autosomal gene that causes this trait, we considered the possibility of Lhx9, a gene encoding LIM homeobox containing transcription factor 9, as a candidate gene. An American cocker spaniel pedigree showing Sry-negative XX sex reversal phenotype was genotyped with an intronic Lhx9 microsatellite marker. Segregation of the Lhx9 marker in the pedigree indicated that a mutation in canine Lhx9 is not likely to be the cause of Sry-negative XX sex reversal. In addition, using the recently available 7.6X canine genomic sequence, we report the location and genomic organization of canine Lhx9.