Chuan-Yun Li

KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases

High-throughput experimental technologies often identify dozens to hundreds of genes related to, or changed in, a biological or pathological process. From these genes one wants to identify biological pathways that may be involved and diseases that may be implicated. Here, we report a web server, KOBAS 2.0, which annotates an input set of genes with putative pathways and disease relationships based on mapping to genes with known annotations. It allows for both ID mapping and cross-species sequence similarity mapping. It then performs statistical tests to identify statistically significantly enriched pathways and diseases. KOBAS 2.0 incorporates knowledge across 1327 species from 5 pathway databases (KEGG PATHWAY, PID, BioCyc, Reactome and Panther) and 5 human disease databases (OMIM, KEGG DISEASE, FunDO, GAD and NHGRI GWAS Catalog). KOBAS 2.0 can be accessed at http://kobas.cbi.pku.edu.cn.

Genes and (common) pathways underlying drug addiction

Drug addiction is a serious worldwide problem with strong genetic and environmental influences. Different technologies have revealed a variety of genes and pathways underlying addiction; however, each individual technology can be biased and incomplete. We integrated 2,343 items of evidence from peer-reviewed publications between 1976 and 2006 linking genes and chromosome regions to addiction by single-gene strategies, microrray, proteomics, or genetic studies. We identified 1,500 human addiction-related genes and developed KARG (http://karg.cbi.pku.edu.cn), the first molecular database for addiction-related genes with extensive annotations and a friendly Web interface. We then performed a meta-analysis of 396 genes that were supported by two or more independent items of evidence to identify 18 molecular pathways that were statistically significantly enriched, covering both upstream signaling events and downstream effects. Five molecular pathways significantly enriched for all four different types of addictive drugs were identified as common pathways which may underlie shared rewarding and addictive actions, including two new ones, GnRH signaling pathway and gap junction. We connected the common pathways into a hypothetical common molecular network for addiction. We observed that fast and slow positive feedback loops were interlinked through CAMKII, which may provide clues to explain some of the irreversible features of addiction.

Species differences in cannabinoid receptor 2 (CNR2 gene): identification of novel human and rodent CB2 isoforms, differential tissue expression, and regulation by cannabinoid receptor ligands

Cannabinoids, endocannabinoids and marijuana activate two well‐characterized cannabinoid receptors (CB‐Rs), CB1‐Rs and CB2‐Rs. The expression of CB1‐Rs in the brain and periphery has been well studied, but neuronal CB2‐Rs have received much less attention than CB1‐Rs. Many studies have now identified and characterized functional glial and neuronal CB2‐Rs in the central nervous system. However, many features of CB2‐R gene structure, regulation and variation remain poorly characterized in comparison with the CB1‐R. In this study, we report on the discovery of a novel human CB2 gene promoter transcribing testis (CB2A) isoform with starting exon located ca 45 kb upstream from the previously identified promoter transcribing the spleen isoform (CB2B). The 5′ exons of both CB2 isoforms are untranslated 5′UTRs and alternatively spliced to the major protein coding exon of the CB2 gene. CB2A is expressed higher in testis and brain than CB2B that is expressed higher in other peripheral tissues than CB2A. Species comparison found that the CB2 gene of human, rat and mouse genomes deviated in their gene structures and isoform expression patterns. mCB2A expression was increased significantly in the cerebellum of mice treated with the CB‐R mixed agonist, WIN55212‐2. These results provide much improved information about CB2 gene structure and its human and rodent variants that should be considered in developing CB2‐R‐based therapeutic agents.

Hominoid-Specific De Novo Protein-Coding Genes Originating from Long Non-Coding RNAs

Tinkering with pre-existing genes has long been known as a major way to create new genes. Recently, however, motherless protein-coding genes have been found to have emerged de novo from ancestral non-coding DNAs. How these genes originated is not well addressed to date. Here we identified 24 hominoid-specific de novo protein-coding genes with precise origination timing in vertebrate phylogeny. Strand-specific RNA–Seq analyses were performed in five rhesus macaque tissues (liver, prefrontal cortex, skeletal muscle, adipose, and testis), which were then integrated with public transcriptome data from human, chimpanzee, and rhesus macaque. On the basis of comparing the RNA expression profiles in the three species, we found that most of the hominoid-specific de novo protein-coding genes encoded polyadenylated non-coding RNAs in rhesus macaque or chimpanzee with a similar transcript structure and correlated tissue expression profile. According to the rule of parsimony, the majority of these hominoid-specific de novo protein-coding genes appear to have acquired a regulated transcript structure and expression profile before acquiring coding potential. Interestingly, although the expression profile was largely correlated, the coding genes in human often showed higher transcriptional abundance than their non-coding counterparts in rhesus macaque. The major findings we report in this manuscript are robust and insensitive to the parameters used in the identification and analysis of de novo genes. Our results suggest that at least a portion of long non-coding RNAs, especially those with active and regulated transcription, may serve as a birth pool for protein-coding genes, which are then further optimized at the transcriptional level.

A Human-Specific De Novo Protein-Coding Gene Associated with Human Brain Functions

To understand whether any human-specific new genes may be associated with human brain functions, we computationally screened the genetic vulnerable factors identified through Genome-Wide Association Studies and linkage analyses of nicotine addiction and found one human-specific de novo protein-coding gene, FLJ33706 (alternative gene symbol C20orf203). Cross-species analysis revealed interesting evolutionary paths of how this gene had originated from noncoding DNA sequences: insertion of repeat elements especially Alu contributed to the formation of the first coding exon and six standard splice junctions on the branch leading to humans and chimpanzees, and two subsequent substitutions in the human lineage escaped two stop codons and created an open reading frame of 194 amino acids. We experimentally verified FLJ33706s mRNA and protein expression in the brain. Real-Time PCR in multiple tissues demonstrated that FLJ33706 was most abundantly expressed in brain. Human polymorphism data suggested that FLJ33706 encodes a protein under purifying selection. A specifically designed antibody detected its protein expression across human cortex, cerebellum and midbrain. Immunohistochemistry study in normal human brain cortex revealed the localization of FLJ33706 protein in neurons. Elevated expressions of FLJ33706 were detected in Alzheimers brain samples, suggesting the role of this novel gene in human-specific pathogenesis of Alzheimers disease. FLJ33706 provided the strongest evidence so far that human-specific de novo genes can have protein-coding potential and differential protein expression, and be involved in human brain functions.

A nonsynonymous SNP in human cytosolic sialidase in a small Asian population results in reduced enzyme activity: potential link with severe adverse reactions to oseltamivir

The use of oseltamivir, widely stockpiled as one of the drugs for use in a possible avian influenza pandemic, has been reported to be associated with neuropsychiatric disorders and severe skin reactions, primarily in Japan. Here we identified a nonsynonymous SNP (single nucleotide polymorphism) in dbSNP database, R41Q, near the enzymatic active site of human cytosolic sialidase, a homologue of virus neuraminidase that is the target of oseltamivir. This SNP occurred in 9.29% of Asian population and none of European and African American population. Our structural analyses and Ki measurements using in vitro sialidase assays indicated that this SNP could increase the unintended binding affinity of human sialidase to oseltamivir carboxylate, the active form of oseltamivir, thus reducing sialidase activity. In addition, this SNP itself results in an enzyme with an intrinsically lower sialidase activity, as shown by its increased Km and decreased Vmax values. Theoretically administration of oseltamivir to people with this SNP might further reduce their sialidase activity. We note the similarity between the reported neuropsychiatric side effects of oseltamivir and the known symptoms of human sialidase-related disorders. We propose that this Asian-enriched sialidase variation caused by the SNP, likely in homozygous form, may be associated with certain severe adverse reactions to oseltamivir.

RNA Editome in Rhesus Macaque Shaped by Purifying Selection

Understanding of the RNA editing process has been broadened considerably by the next generation sequencing technology; however, several issues regarding this regulatory step remain unresolved – the strategies to accurately delineate the editome, the mechanism by which its profile is maintained, and its evolutionary and functional relevance. Here we report an accurate and quantitative profile of the RNA editome for rhesus macaque, a close relative of human. By combining genome and transcriptome sequencing of multiple tissues from the same animal, we identified 31,250 editing sites, of which 99.8% are A-to-G transitions. We verified 96.6% of editing sites in coding regions and 97.5% of randomly selected sites in non-coding regions, as well as the corresponding levels of editing by multiple independent means, demonstrating the feasibility of our experimental paradigm. Several lines of evidence supported the notion that the adenosine deamination is associated with the macaque editome – A-to-G editing sites were flanked by sequences with the attributes of ADAR substrates, and both the sequence context and the expression profile of ADARs are relevant factors in determining the quantitative variance of RNA editing across different sites and tissue types. In support of the functional relevance of some of these editing sites, substitution valley of decreased divergence was detected around the editing site, suggesting the evolutionary constraint in maintaining some of these editing substrates with their double-stranded structure. These findings thus complement the “continuous probing” model that postulates tinkering-based origination of a small proportion of functional editing sites. In conclusion, the macaque editome reported here highlights RNA editing as a widespread functional regulation in primate evolution, and provides an informative framework for further understanding RNA editing in human.

Hydrogen peroxide primes heart regeneration with a derepression mechanism

While the adult human heart has very limited regenerative potential, the adult zebrafish heart can fully regenerate after 20% ventricular resection. Although previous reports suggest that developmental signaling pathways such as FGF and PDGF are reused in adult heart regeneration, the underlying intracellular mechanisms remain largely unknown. Here we show that H2O2 acts as a novel epicardial and myocardial signal to prime the heart for regeneration in adult zebrafish. Live imaging of intact hearts revealed highly localized H2O2 (∼30 μM) production in the epicardium and adjacent compact myocardium at the resection site. Decreasing H2O2 formation with the Duox inhibitors diphenyleneiodonium (DPI) or apocynin, or scavenging H2O2 by catalase overexpression markedly impaired cardiac regeneration while exogenous H2O2 rescued the inhibitory effects of DPI on cardiac regeneration, indicating that H2O2 is an essential and sufficient signal in this process. Mechanistically, elevated H2O2 destabilized the redox-sensitive phosphatase Dusp6 and hence increased the phosphorylation of Erk1/2. The Dusp6 inhibitor BCI achieved similar pro-regenerative effects while transgenic overexpression of dusp6 impaired cardiac regeneration. H2O2 plays a dual role in recruiting immune cells and promoting heart regeneration through two relatively independent pathways. We conclude that H2O2 potentially generated from Duox/Nox2 promotes heart regeneration in zebrafish by unleashing MAP kinase signaling through a derepression mechanism involving Dusp6.

Genome wide association for substance dependence: convergent results from epidemiologic and research volunteer samples

BackgroundDependences on addictive substances are substantially-heritable complex disorders whose molecular genetic bases have been partially elucidated by studies that have largely focused on research volunteers, including those recruited in Baltimore. Maryland. Subjects recruited from the Baltimore site of the Epidemiological Catchment Area (ECA) study provide a potentially-useful comparison group for possible confounding features that might arise from selecting research volunteer samples of substance dependent and control individuals. We now report novel SNP (single nucleotide polymorphism) genome wide association (GWA) results for vulnerability to substance dependence in ECA participants, who were initially ascertained as members of a probability sample from Baltimore, and compare the results to those from ethnically-matched Baltimore research volunteers.ResultsWe identify substantial overlap between the home address zip codes reported by members of these two samples. We find overlapping clusters of SNPs whose allele frequencies differ with nominal significance between substance dependent vs control individuals in both samples. These overlapping clusters of nominally-positive SNPs identify 172 genes in ways that are never found by chance in Monte Carlo simulation studies. Comparison with data from human expressed sequence tags suggests that these genes are expressed in brain, especially in hippocampus and amygdala, to extents that are greater than chance.ConclusionThe convergent results from these probability sample and research volunteer sample datasets support prior genome wide association results. They fail to support the idea that large portions of the molecular genetic results for vulnerability to substance dependence derive from factors that are limited to research volunteers.

Species Differences in Cannabinoid Receptor 2 and Receptor Responses to Cocaine Self-Administration in Mice and Rats

The discovery of functional cannabinoid receptors 2 (CB2Rs) in brain suggests a potential new therapeutic target for neurological and psychiatric disorders. However, recent findings in experimental animals appear controversial. Here we report that there are significant species differences in CB2R mRNA splicing and expression, protein sequences, and receptor responses to CB2R ligands in mice and rats. Systemic administration of JWH133, a highly selective CB2R agonist, significantly and dose-dependently inhibited intravenous cocaine self-administration under a fixed ratio (FR) schedule of reinforcement in mice, but not in rats. However, under a progressive ratio (PR) schedule of reinforcement, JWH133 significantly increased breakpoint for cocaine self-administration in rats, but decreased it in mice. To explore the possible reasons for these conflicting findings, we examined CB2R gene expression and receptor structure in the brain. We found novel rat-specific CB2C and CB2D mRNA isoforms in addition to CB2A and CB2B mRNA isoforms. In situ hybridization RNAscope assays found higher levels of CB2R mRNA in different brain regions and cell types in mice than in rats. By comparing CB2R-encoding regions, we observed a premature stop codon in the mouse CB2R gene that truncated 13 amino-acid residues including a functional autophosphorylation site in the intracellular C-terminus. These findings suggest that species differences in the splicing and expression of CB2R genes and receptor structures may in part explain the different effects of CB2R-selective ligands on cocaine self-administration in mice and rats.