D3G Release 23.02 ================== D3G provides a set of genomic data to support the development of oligonucleotide therapeutics. Genome sequences and gene models are compiled for seven species, human, crab-eating macaque, rhesus macaque, common marmoset, mouse, rat, and rabbit. Expression data are supplemented for two species of crab-eating macaque and common marmoset. What's new ----------- * Human genome sequence is now hg38, patch 14 (updated from patch 13) * RefSeq genes were updated for human and mouse * Rhesus macaque (Macaca mulatta) and rabbit (Oryctolagus cuniculus) are additionally supported. * An experimental track of human long read RNA-seq is included in the genome browser visualization Data use / embargo ------------------- Here we open our original data sets to support development of oligonucleotide therapeutics under [CC-BY](https://creativecommons.org/licenses/by/4.0/). We encourage anyone (for example, in academic institute, commercial company, regulatory agency, or others) to use the data sets freely for development or assessment of drugs. Meanwhile, we request users to respect the embargo on the publication of genome-wide analysis based on this data set, as we are currently preparing manuscripts to provide genome-scale analysis of the non-human primates with complete description of the experimental and computational details with raw data. Exceptions to the policy are for analyses on a couple of locus, gene families, and oligonucleotide sequences, rather than comprehensive large-scale analysis. Data source ------------ * Human (Homo sapiens) - Genome assembly: GRCh38/hg38.p14 (GCA_000001405.29) - Gene models - RefSeq: NCBI Homo sapiens Release 110 (2022-04-12) - Gencode: V42 * Crab-eating macaque (Cynomolgus macaque) - Genome assembly: macFasRKS1912v2 (GCF_012559485.2) - Gene models - RefSeq: NCBI Macaca fascicularis Annotation Release 102 (direct download from NCBI) - TRaC: 22.02 - built from CAGE and RNA-seq profiles (data from PRJDB9546 and others), and MinION sequencing of full-length cDNAs prepared by cap-trapping. - Expression: CAGE profiles obtained from the above * Rhesus macaque (Macaca mulatta) - Genome assembly: rheMac10 (GCF_003339765.1) - Gene models - RefSeq: Macaca mulatta Annotation Release 103 (2019-06-26) * Common marmoset (Callithrix jacchus) - Genome assembly: calJac4 (GCF_009663435.1) - Gene models - RefSeq: NCBI Callithrix jacchus Annotation Release 105 (2020-07-11) - TRaC: 22.02 - built from CAGE profiles and RNA-seq profiles (data from PRJDB9547 and others), and MinION sequencing of full-length cDNAs prepared by cap-trapping. - Expression: CAGE profiles obtained from the above * Mouse (Mus musculus) - Genome assembly: GRCm39/mm39 (GCA_000001635.9) - Gene models - RefSeq: NCBI Mus musculus Annotation Release 109 (2022-11-25) - Gencode: VM31 * Rat (Rattus norvegicus) - Genome assembly: mRatBN7.2/rn7 (GCA_000001895.4) - Gene models - RefSeq: NCBI Rattus norvegicus Annotation Release 108 (2021-01-21) * Rabbit (Oryctolagus cuniculus) - Genome assembly: UM_NZW_1.0chr, a modified version of UM_NZW_1.0 (GCF_009806435.1) where chromosomal names (such as chr1, chr2, etc) are used as sequence names. - Gene models - RefSeq: NCBI Oryctolagus cuniculus Annotation Release 103 (2022-11-08) The genome assembly of crab-eating macaque (GCF_012559485.2) is a polished version of the one previously described [^1]. The data files, genome assemblies and gene models, were directly downloaded from NCBI and compiled for D3G. The genome assembly and gene models of rabbit were also directly downloaded from NCBI, and compiled for D3G with a modification of using chromosomal names (such as chr1, chr2, etc) as sequence names. Expression data for crab-eating macaque and common marmoset are based on our own CAGE data profiles. The other data sets, reference genome [^2] and gene models (RefSeq [^3] and Gencode [^4]) for human, mouse, rhesus macaque, and rat) were obtained from [the UCSC Genome Browser Database](https://hgdownload.soe.ucsc.edu/downloads.html) in Feb 2023. [^1]: Jayakumar V, Nishimura O, Kadota M, Hirose N, Sano H, Murakawa Y, Yamamoto Y, Nakaya M, Tsukiyama T, Seita Y, Nakamura S, Kawai J, Sasaki E, Ema M, Kuraku S, Kawaji H, Sakakibara Y. Chromosomal-scale de novo genome assemblies of Cynomolgus Macaque and Common Marmoset. Sci Data. 2021;8:159. doi: https://dx.doi.org/10.1038/s41597-021-00935-6 [^2]: Church DM, Schneider VA, Graves T, et al. Modernizing reference genome assemblies. PLoS Biol. 2011;9(7):e1001091. doi:10.1371/journal.pbio.1001091 [^3]: O'Leary NA, Wright MW, Brister JR, et al. Reference sequence (RefSeq) database at NCBI: current status, taxonomic expansion, and functional annotation. Nucleic Acids Res. 2016;44(D1):D733-D745. doi:10.1093/nar/gkv1189 [^4]: Harrow J, Frankish A, Gonzalez JM, et al. GENCODE: the reference human genome annotation for The ENCODE Project. Genome Res. 2012;22(9):1760-1774. doi:10.1101/gr.135350.111 Data files ----------- ``` |-- md5sum.txt |-- hg38 | |-- gene | | |-- hg38_ncbiRefSeqCurated.bed12.gz | | |-- hg38_ncbiRefSeqCuratedProteinCoding.bed12.gz | | |-- hg38_ncbiRefSeqCuratedProteinCoding_prespliced.fa.gz | | |-- hg38_ncbiRefSeqCuratedProteinCoding_spliced.fa.gz | | |-- hg38_ncbiRefSeqCurated_prespliced.fa.gz | | |-- hg38_ncbiRefSeqCurated_spliced.fa.gz | | |-- hg38_ncbiRefSeqPredicted.bed12.gz | | `-- hg38_wgEncodeGencodeCompV42.bed12.gz | `-- genome | `-- hg38.p14.fa.gz |-- macFasRKS1912v2 | |-- expression | | |-- macFasRKS1912v2_ncbiRefseqPromoterExpCounts.txt.gz | | `-- macFasRKS1912v2_ncbiRefseqPromoterExpCpm.txt.gz | |-- gene | | |-- macFasRKS1912v2_Trac2202.bed12.gz | | |-- macFasRKS1912v2_ncbiRefSeq.bed12.gz | | |-- macFasRKS1912v2_ncbiRefSeq_prespliced.fa.gz | | `-- macFasRKS1912v2_ncbiRefSeq_spliced.fa.gz | `-- genome | `-- macFasRKS1912v2.fa.gz |-- rheMac10 | |-- gene | | |-- rheMac10_ncbiRefSeq.bed12.gz | | |-- rheMac10_ncbiRefSeq_prespliced.fa.gz | | `-- rheMac10_ncbiRefSeq_spliced.fa.gz | `-- genome | `-- rheMac10.fa.gz |-- calJac4 | |-- expression | | |-- calJac4_ncbiRefseqPromoterExpCounts.txt.gz | | `-- calJac4_ncbiRefseqPromoterExpCpm.txt.gz | |-- gene | | |-- calJac4_Trac2202.bed12.gz | | |-- calJac4_ncbiRefSeq.bed12.gz | | |-- calJac4_ncbiRefSeq_prespliced.fa.gz | | `-- calJac4_ncbiRefSeq_spliced.fa.gz | `-- genome | `-- calJac4.fa.gz |-- mm39 | |-- gene | | |-- mm39_ncbiRefSeqCurated.bed12.gz | | |-- mm39_ncbiRefSeqCuratedProteinCoding.bed12.gz | | |-- mm39_ncbiRefSeqCuratedProteinCoding_prespliced.fa.gz | | |-- mm39_ncbiRefSeqCuratedProteinCoding_spliced.fa.gz | | |-- mm39_ncbiRefSeqCurated_prespliced.fa.gz | | |-- mm39_ncbiRefSeqCurated_spliced.fa.gz | | |-- mm39_ncbiRefSeqPredicted.bed12.gz | | `-- mm39_wgEncodeGencodeCompVM31.bed12.gz | `-- genome | `-- mm39.fa.gz `-- rn7 | |-- gene | | |-- rn7_ncbiRefSeqCurated.bed12.gz | | |-- rn7_ncbiRefSeqCuratedProteinCoding.bed12.gz | | |-- rn7_ncbiRefSeqCuratedProteinCoding_prespliced.fa.gz | | |-- rn7_ncbiRefSeqCuratedProteinCoding_spliced.fa.gz | | |-- rn7_ncbiRefSeqCurated_prespliced.fa.gz | | |-- rn7_ncbiRefSeqCurated_spliced.fa.gz | | `-- rn7_ncbiRefSeqPredicted.bed12.gz | `-- genome | `-- rn7.fa.gz |-- UM_NZW_1.0chr |-- gene | |-- UM_NZW_1.0chr_ncbiRefSeq.bed12.gz | |-- UM_NZW_1.0chr_ncbiRefSeq_prespliced.fa.gz | `-- UM_NZW_1.0chr_ncbiRefSeq_spliced.fa.gz `-- genome `-- UM_NZW_1.0chr.fa.gz ``` Files with `.fa.gz` suffix in their names contain nucleotide sequences in FASTA format, ones with `.bed12.gz` suffix contain exon-intron structure of gene models in [BED12 format (BED format with 12 columns)](https://genome.ucsc.edu/FAQ/FAQformat.html#format1), and ones with `.matrix.txt.gz` suffix contain tab-delimited text of expression intensities. Nucleotide sequences for gene models are compiled by using their genomic coordinates and the genome sequences, only for a selected subset of the gene models. Files containing `prespliced` in their names represent immature transcripts before splicing, the same to pre-mRNA for protein coding transcripts (the term "pre-spliced" is used to be compatible with long noncoding RNAs). Note that the protein coding transcripts here include ones on both nuclear and mitochondrial DNAs. It means that protein coding transcripts compiled based on RefSeq [^2] includes `YP_` (for human) or `NP_` (for mouse) entries in addition to `NM_` ones. Identifier (ID) for RefSeq gene models consists of `REFSEQ_TRANSCRIPT_ID|GENE_SYMBOL;ENTREZ_GENE_ID`, as seen in the example of `NM_000454.4|SOD1;6647` in the `hg38NcbiRefSeqCurated.bed12.gz` file. In FASTA files of nucleotide sequences, other information such as genomic coordinates are concatenated after `|`. As for expression data , files containing `ExpCounts` in their names represent counts of 5'-ends of CAGE read alignments with the genome sequences, and the ones with `ExpCpm` represent expression intensities where the read counts were normalized by CPM, counts per million with RLE (relative log expression) method [^5]. The files which names include 'ncbiRefseqPromoter' indicate promoter level expression data (according to CAGE peaks) of RefSeq gene models (not CAGE peaks). Experimental data and their computational processing ----------------------------------------------------- We generated the following data to construct 5'-end complete gene models by a novel approach (TRaC, Transcript models based on RNA-seq, CAGE, and long read sequencing by Oxford Nanopore sequencing). * RNA-seq * CAGE, based on ssCAGE protocol [^6] * Long read (Oxford Nanopore) sequencing of full-length cDNA prepared by cap-trapping [^5]: Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol 2010; 11:R106 doi: http://dx.doi.org/10.1186/gb-2010-11-10-r106 [^6]: Morioka MS, Kawaji H, Nishiyori-Sueki H, et al. Cap Analysis of Gene Expression (CAGE): A Quantitative and Genome-Wide Assay of Transcription Start Sites. Methods Mol Biol. 2020;2120:277-301. doi:10.1007/978-1-0716-0327-7_20 Leadership and contact ---------------------- The original data is made in a collaboration among [RIKEN PMI](https://www.riken.jp/research/labs/pmi/), [RIKEN IMS](https://www.riken.jp/research/labs/ims/), [Shiga University of Medical Science](https://www.shiga-med.ac.jp/), [CIEA (Central Institute for Experimental Animals)](https://www.ciea.or.jp/), [Keio University](https://www.keio.ac.jp/), [RIKEN BDR](https://www.riken.jp/research/labs/bdr/), [DBCLS (Database Center for Life Science)](http://dbcls.rois.ac.jp/), [NIBIOHN (National Institutes of Biomedical Innovation, Health and Nutrition)](https://www.nibiohn.go.jp/), and [TMIMS (Tokyo Metropolitan Institute of Medical Science)](http://www.igakuken.or.jp/). Please contact us via e-mail below for any questions, comments, suggestions, or collaboration: d3g@ml.riken.jp How to cite ------------ Please refer our database like this: D3G: Database for Drug Development based on Genome and RNA sequences, https://d3g.riken.jp, 2023 Acknowledgement ---------------- We thank to [AMED (Japan Agency for Medical Research and Development)](https://www.amed.go.jp/), [NIHS (National Institute of Health Sciences)](http://www.nihs.go.jp/), [JPMA (Japan Pharmaceutical Manufacturers Association)](http://www.jpma.or.jp/), and [the FANTOM consortium](https://fantom.gsc.riken.jp/) for relevant advices and fruitful discussion. The experiments and the database is financially supported by AMED under Grant Number JP17kk0305008 and JP20kk0305013.