High-throughput and high-accuracy single-cell RNA isoform analysis using PacBio circular consensus sequencing

Although long-read single-cell RNA isoform sequencing (scISO-Seq) can reveal alternative RNA splicing in individual cells, it suffers from a low read throughput. Here, we introduce HIT-scISOseq, a method that removes most artifact cDNAs and concatenates multiple cDNAs for PacBio circular consensus sequencing (CCS) to achieve high-throughput and high-accuracy single-cell RNA isoform sequencing. HIT-scISOseq can yield >10 million high-accuracy long-reads in a single PacBio Sequel II SMRT Cell 8M. We also report the development of scISA-Tools that demultiplex HIT-scISOseq concatenated reads into single-cell cDNA reads with >99.99% accuracy and specificity. We apply HIT-scISOseq to characterize the transcriptomes of 3375 corneal limbus cells and reveal cell-type-specific isoform expression in them. HIT-scISOseq is a high-throughput, high-accuracy, technically accessible method and it can accelerate the burgeoning field of long-read single-cell transcriptomics.

Engineered T cells and their therapeutic applications in autoimmune diseases

Chimeric antigen receptor T cells (CAR-T cells) are engineered recombinant T cells, which were initially used to treat hematopoietic malignancies and are now widely used in the treatment of various diseases. Considering their intrinsic targeting efficiency, CAR-T cells show considerable potential in the treatment of autoimmune diseases. Furthermore, regulatory T cells (Treg), a subset of CD4 T cells exhibiting immunosuppressive functions, have attracted increasing attention regarding CAR-Treg cell production. In this review, we report on recent developments in preclinical and clinical studies on CAR-T cells in autoimmune diseases and provide an outlook on opportunities and challenges of CAR-T application in such diseases.

Forensic nanopore sequencing of STRs and SNPs using Verogen's ForenSeq DNA Signature Prep Kit and MinION

The MinION nanopore sequencing device (Oxford Nanopore Technologies, Oxford, UK) is the smallest commercially available sequencer and can be used outside of conventional laboratories. The use of the MinION for forensic applications, however, is hindered by the high error rate of nanopore sequencing. One approach to solving this problem is to identify forensic genetic markers that can consistently be typed correctly based on nanopore sequencing. In this pilot study, we explored the use of nanopore sequencing for single nucleotide polymorphism (SNP) and short tandem repeat (STR) profiling using Verogen’s (San Diego, CA, USA) ForenSeq DNA Signature Prep Kit. Thirty single-contributor samples and DNA standard material 2800 M were genotyped using the Illumina (San Diego, CA, USA) MiSeq FGx and MinION (with R9.4.1 flow cells) devices. With an optimized cutoff for allelic imbalance, all 94 identity-informative SNP loci could be genotyped reliably using the MinION device, with an overall accuracy of 99.958% (1 error among 2926 genotypes). STR typing was notably error prone, and its accuracy was locus dependent. We developed a custom-made bioinformatics workflow, and finally selected 13 autosomal STRs, 14 Y-STRs, and 4 X-STRs showing high consistency between nanopore and Illumina sequencing among the tested samples. These SNP and STR loci could be candidates for panel design for forensic analysis based on nanopore sequencing.

N6-Methyladenine DNA Modification in the Human Genome

DNA N⁶-methyladenine (6mA) modification is the most prevalent DNA modification in prokaryotes, but whether it exists in human cells and whether it plays a role in human diseases remain enigmatic. Here, we showed that 6mA is extensively present in the human genome, and we cataloged 881,240 6mA sites accounting for ∼0.051% of the total adenines. [G/C]AGG[C/T] was the most significantly associated motif with 6mA modification. 6mA sites were enriched in the coding regions and mark actively transcribed genes in human cells. DNA 6mA and N⁶-demethyladenine modification in the human genome were mediated by methyltransferase N6AMT1 and demethylase ALKBH1, respectively. The abundance of 6mA was significantly lower in cancers, accompanied by decreased N6AMT1 and increased ALKBH1 levels, and downregulation of 6mA modification levels promoted tumorigenesis. Collectively, our results demonstrate that DNA 6mA modification is extensively present in human cells and the decrease of genomic DNA 6mA promotes human tumorigenesis.

A 472-SNP panel for pairwise kinship testing of second-degree relatives

Kinship testing based on genetic markers, as forensic short tandem repeats (STRs) and single nucleotide polymorphisms (SNPs), has valuable practical applications. Paternity and first-degree relationship can be accurately identified by current commonly-used forensic STRs and reported SNP markers. However, second-degree and more distant relationships remain challenging. Although ∼10⁵-10⁶ SNPs can be used to estimate relatedness of higher degrees, genome-wide genotyping and analysis may be impractical for forensic use. With rapid growth of human genome data sets, it is worthwhile to explore additional markers, especially SNPs, for kinship analysis. Here, we reported an autosomal SNP panel consisted of 342 SNP selected from >84 million SNPs and 131 SNPs from previous systems. We genotyped these SNPs in 136 Chinese individuals by multiplex amplicon Massively Parallel Sequencing, and performed pairwise gender-independent kinship testing. The specificity and sensitivity of these SNPs to distinguish second-degree relatives and the unrelated was 99.9% and 100%, respectively, compared with 53.7% and 99.9% of 19 commonly-used forensic STRs. Moreover, the specificity increased to 100% by the combined use of these STRs and SNPs. The 472-SNP panel could also greatly facilitate the discrimination among different relationships. We estimated that the power of ∼6.45 SNPs were equivalent to one forensic STR in the scenario of 2nd-degree relative pedigree. Altogether, we proposed a panel of 472 SNP markers for kinship analysis, which could be important supplementary of current forensic STRs to solve the problem of second-degree relative testing.

Neuroendocrine immunomodulation network dysfunction in SAMP8 mice and PrP-hAβPPswe/PS1ΔE9 mice: potential mechanism underlying cognitive impairment

Senescence-accelerated mouse prone 8 strain (SAMP8) and PrP-hAβPPswe/PS1ΔE9 (APP/PS1) mice are classic animal models of sporadic Alzheimer's disease and familial AD respectively. Our study showed that object recognition memory, spatial learning and memory, active and passive avoidance were deteriorated and neuroendocrine immunomodulation (NIM) network was imbalance in SAMP8 and APP/PS1 mice. SAMP8 and APP/PS1 mice had their own specific phenotype of cognition, neuroendocrine, immune and NIM molecular network. The endocrine hormone corticosterone, luteinizing hormone and follicle-stimulating hormone, chemotactic factor monocyte chemotactic protein-1, macrophage inflammatory protein-1β, regulated upon activation normal T cell expressed and secreted factor and eotaxin, pro-inflammatory factor interleukin-23, and the Th1 cell acting as cell immunity accounted for cognitive deficiencies in SAMP8 mice, while adrenocorticotropic hormone and gonadotropin-releasing hormone, colony stimulating factor granulocyte colony stimulating factor, and Th2 cell acting as humoral immunity in APP/PS1 mice. On the pathway level, chemokine signaling and T cell receptor signaling pathway played the key role in cognition impairments of two models, while cytokine-cytokine receptor interaction and natural killer cell mediated cytotoxicity were more important in cognitive deterioration of SAMP8 mice than APP/PS1 mice. This mechanisms of NIM network underlying cognitive impairment is significant for further understanding the pathogenesis of AD and can provide useful information for development of AD therapeutic drug.

Partners of patients with ulcerative colitis exhibit a biologically relevant dysbiosis in fecal microbial metacommunities

AIM

To investigate alterations in the fecal microbiome using 16S rRNA amplicon sequencing in couples in the same cohabitation environment.

METHODS

Fecal samples were collected from eight ulcerative colitis (UC) patients and their healthy partners at Lishui People’s Hospital, Zhejiang Province, China. DNA was extracted and the variable regions V3 and V4 of the 16S rRNA genes were PCR amplified using a two-step protocol. Clear reads were clustered into operational taxonomic units (OTUs) at the 97% sequence similarity level using UCLUST v1.2.22. The Wilcoxon rank-sum test (R v3.1.2) was used to compare inter-individual differences. Differences with a P value < 0.05 were considered statistically significant.

RESULTS

Fecal microbial communities were more similar among UC patients than their healthy partners (P = 0.024). UC individuals had a lower relative abundance of bacteria belonging to the Firmicutes, especially Blautia, Clostridium, Coprococcus and Roseburia (P < 0.05). Microbiota dysbiosis was detected in UC patients and their healthy partners. Relevant genera included Akkermansiam, Bacteroides, Escherichia, Lactobacillales, Klebsiella and Parabacteroides. The enriched pathways in fecal samples of UC patients were related to lipid and nucleotide metabolism. Additionally, the pathways involved in membrane transport and metabolism of cofactors and vitamins were more abundant in the healthy partners.

CONCLUSION

Our results suggested that the microbial composition might be affected in healthy partners cohabiting with UC patients, especially in terms of microbiota dysbiosis.

Exploring the efficacy of paternity and kinship testing based on single nucleotide polymorphisms

Short tandem repeats (STRs) are conventional genetic markers typically used for paternity and kinship testing. As supplementary markers of STRs, single nucleotide polymorphisms (SNPs) have less discrimination power but broader applicability to degraded samples. The rapid improvement of next-generation sequencing (NGS) and multiplex amplification technologies also make it possible now to simultaneously identify dozens or even hundreds of SNP loci in a single pool. However, few studies have been endeavored to kinship testing based on SNP loci. In this study, we genotyped 90 autosomal human identity SNP loci with NGS, and investigated their testing efficacies based on the likelihood ratio model in eight pedigree scenarios involving paternity, half/full-sibling, uncle/nephew, and first-cousin relationships. We found that these SNPs might be sufficient to discriminate paternity and full-sibling, but impractical for more distant relatives such as uncle and cousin. Furthermore, we conducted an in silico study to obtain the theoretical tendency of how testing efficacy varied with increasing number of SNP loci. For each testing battery in a given pedigree scenario, we obtained distributions of logarithmic likelihood ratio for both simulated relatives and unrelated controls. The proportion of the overlapping area between the two distributions was defined as a false testing level (FTL) to evaluate the testing efficacy. We estimated that 85, 127, 491, and 1,858 putative SNP loci were required to discriminate paternity, full-sibling, half-sibling/uncle-nephew, and first-cousin (FTL, 0.1%), respectively. To test a half-sibling or nephew, an additional uncle relative could be included to decrease the required number of putative SNP loci to ∼320 (FTL, 0.1%). As a systematic computation of paternity and kinship testing based only on SNPs, our results could be informative for further studies and applications on paternity and kinship testing using SNP loci.

Genetic diversity and evolutionary dynamics of Ebola virus in Sierra Leone

A novel Ebola virus (EBOV) first identified in March 2014 has infected more than 25,000 people in West Africa, resulting in more than 10,000 deaths. Preliminary analyses of genome sequences of 81 EBOV collected from March to June 2014 from Guinea and Sierra Leone suggest that the 2014 EBOV originated from an independent transmission event from its natural reservoir followed by sustained human-to-human infections. It has been reported that the EBOV genome variation might have an effect on the efficacy of sequence-based virus detection and candidate therapeutics. However, only limited viral information has been available since July 2014, when the outbreak entered a rapid growth phase. Here we describe 175 full-length EBOV genome sequences from five severely stricken districts in Sierra Leone from 28 September to 11 November 2014. We found that the 2014 EBOV has become more phylogenetically and genetically diverse from July to November 2014, characterized by the emergence of multiple novel lineages. The substitution rate for the 2014 EBOV was estimated to be 1.23 × 10(-3) substitutions per site per year (95% highest posterior density interval, 1.04 × 10(-3) to 1.41 × 10(-3) substitutions per site per year), approximating to that observed between previous EBOV outbreaks. The sharp increase in genetic diversity of the 2014 EBOV warrants extensive EBOV surveillance in Sierra Leone, Guinea and Liberia to better understand the viral evolution and transmission dynamics of the ongoing outbreak. These data will facilitate the international efforts to develop vaccines and therapeutics.

Nodes and biological processes identified on the basis of network analysis in the brain of the senescence accelerated mice as an Alzheimer's disease animal model

Harboring the behavioral and histopathological signatures of Alzheimer's disease (AD), senescence accelerated mouse-prone 8 (SAMP8) mice are currently considered a robust model for studying AD. However, the underlying mechanisms, prioritized pathways and genes in SAMP8 mice linked to AD remain unclear. In this study, we provide a biological interpretation of the molecular underpinnings of SAMP8 mice. Our results were derived from differentially expressed genes in the hippocampus and cerebral cortex of SAMP8 mice compared to age-matched SAMR1 mice at 2, 6, and 12 months of age using cDNA microarray analysis. On the basis of PPI, MetaCore and the co-expression network, we constructed a distinct genetic sub-network in the brains of SAMP8 mice. Next, we determined that the regulation of synaptic transmission and apoptosis were disrupted in the brains of SAMP8 mice. We found abnormal gene expression of RAF1, MAPT, PTGS2, CDKN2A, CAMK2A, NTRK2, AGER, ADRBK1, MCM3AP, and STUB1, which may have initiated the dysfunction of biological processes in the brains of SAMP8 mice. Specifically, we found microRNAs, including miR-20a, miR-17, miR-34a, miR-155, miR-18a, miR-22, miR-26a, miR-101, miR-106b, and miR-125b, that might regulate the expression of nodes in the sub-network. Taken together, these results provide new insights into the biological and genetic mechanisms of SAMP8 mice and add an important dimension to our understanding of the neuro-pathogenesis in SAMP8 mice from a systems perspective.

Cellular microRNA let-7c inhibits M1 protein expression of the H1N1 influenza A virus in infected human lung epithelial cells

The influenza virus (IV) triggers a series of signalling events inside host cells and induces complex cellular responses. Studies have suggested that host factors play an essential role in IV replication. MicroRNAs (miRNAs) represent a class of small non-coding RNAs that target mRNAs, triggering either translation repression or RNA degradation. Emerging research suggests that host-derived cellular miRNAs are involved in mediating the host–IV interaction. Using miRNA microarrays, we identified several miRNAs aberrantly expressed in IV-infected human lung epithelial cells (A549). Specifically, miR-let-7c was highly up-regulated in IV-infected A549 cells. PITA and miRanda database screening indicated that the let-7c seed sequence is a perfect complementary sequence match to the 3′ untranslated region (UTR) of viral gene M1 (+) cRNA, but not to PB2 and PA. As detected by a luciferase reporter system, let-7c directly targeted the 3′-UTR of M1 (+) cRNA, but not PB2 and PA. To experimentally identify the function of cellular let-7c, precursor let-7c was transfected into A549 cells. Let-7c down-regulated IV M1 expression at both the (+) cRNA and protein levels. Furthermore, transfection with a let-7c inhibitor enhanced the expression of M1. Therefore, let-7c may reduce IV replication by degrading M1 (+) cRNA. This is the first report indicating that cellular miRNA regulates IV replication through the degradation of viral gene (+) cRNA by matching the 3′-UTR of the viral cRNA. These findings suggest that let-7c plays a role in protecting host cells from the virus in addition to its known cellular functions.

[Screening for genes associated with cardiac fibrosis induced by aldosterone]

AIM

To investigate differently expressed genes associated with cardiac fibrosis induced independently by aldosterone.

METHODS

Fetal cardiac fibroblasts (FCFs)were isolated and cultured. Total RNA was extracted 8 hours after aldosterone administration. Then gene chips were used to screen these RNA samples. Some of candidate genes were confirmed by RT-PCR and Western blot.

RESULTS

Differently expressed 1519 genes were screened. Up-regulated genes were 714 while down-regulated genes were 805. The expression of CCL7, MMP-26 and IL31RA was tested by RT-PCR and western blot, the results is identical with those by gene chips.

CONCLUSION

Gene chip can efficiently single out differently expressed genes induced dependently by aldosterone in FCFs. CCL7, MMP-26 and IL31RA may be associated with cardiac fibrosis induced by aldosterone.

Analysis of synonymous codon usage bias in 09H1N1

A novel subtype of influenza A virus 09H1N1 has rapidly spread across the world. Evolutionary analyses of this virus have revealed that 09H1N1 is a triple reassortant of segments from swine, avian and human influenza viruses. In this study, we investigated factors shaping the codon usage bias of 09H1N1 and carried out cluster analysis of 60 strains of influenza A virus from different subtypes based on their codon usage bias. We discovered that more preferentially used codons of 09H1N1 are A-ended or U-ended, and the intra-genomic codon usage bias of 09H1N1 is quite low. Base composition constraint, dinucleotide biases and translational selection are the main factors influencing the codon usage bias of 09H1N1. At the genome level, we find that the codon usage bias of 09H1N1 is similar to H1N1 (A/swine/Kansas/77778/2007H1N1), H9N2 from Asia, H1N2 from Asia and North America and H3N2 from North America. Our results provide insight for understanding the processes governing evolution, regulation of gene expression, and revealing the evolution of 09H1N1.

[A genetic optimization designing method for microorganism detection genechip probe based on genetic algorithm]

A new automatic selection approach of microorganism specific fragment combination is presented in this paper. Genetic algorithm is used to search optimal solution on the basis of classification ability of SNP combination, which is evaluated by the rough set theory. Other related experimental parameters are also been incorporated. Experimental results show that the method can find the best SNP combination pattern efficiently and accurately, which implies that it is a reliable approach to the genechip probe design.

Protocatechuic aldehyde inhibits hepatitis B virus replication both in vitro and in vivo

Natural compounds provide a large reservoir of potentially active anti-hepatitis B virus (HBV) agents. We examined the direct effects of protocatechuic aldehyde (PA; derived from the Chinese herb, Salvia miltiorrhiza) on HBV replication in HepG2 2.2.15 cell line and duck hepatitis B virus (DHBV) replication in ducklings in vivo. The extracellular HBV DNA, hepatitis B e antigen (HBeAg) and hepatitis B surface antigen (HBsAg) concentrations in cell culture medium were determined by quantitative real-time PCR and ELISA, respectively. DHBV in duck serum was analyzed by dot blot. PA appeared to downregulate the secretion of HBsAg and HBeAg as well as the release of HBV DNA from HepG2 2.2.15 in a dose- and time-dependent manner at concentrations between 24 and 48 microg/mL. PA (25, 50, or 100 mg/kg, intraperitoneally, twice daily) also reduced viremia in DHBV-infected ducks. We provide the first evidence that PA, a novel anti-HBV substance derived from traditional Chinese herb S. miltiorrhiza, can efficiently inhibits HBV replication in HepG2 2.2.15 cell line in vitro and inhibit DHBV replication in ducks in vivo. PA therefore warrants further investigation as a potential therapeutic agent for HBV infections.

[Sequence analysis of the complete genome of bocavirus WLL-1]

Human bocavirus, which was firstly discovered in 2005, is a new human parvovirus associated with lower respiratory tract infection in children. In this study, a human bocavirus, named WLL-1 isolate, was identified in Wenlin County, Zhejiang Province. The genome of bocavirus WLL-1 has been sequenced and analyzed. Phylogenetic analyses showed that WLL-1 shares 99% homology with other bocaviruses recently reported, but also has some special variations.

Antisense oligonucleotides targeted against asialoglycoprotein receptor 1 block human hepatitis B virus replication

Chronic hepatitis B virus (HBV) infection is a major worldwide public health problem. Better therapeutics and treatment strategies are urgently needed because of ineffective clinical treatment. Our previous study showed that asialoglycoprotein receptor 1 (ASGPR1) was upregulated by HBV but downregulated by lamivudine in HepG2.2.15 cells. It has also been reported that ASGPR is a candidate receptor for HBV attachment to hepatocytes. Therefore, as a major subunit of ASGPR, ASGPR1, might be a potential target for anti-HBV drugs. To validate this hypothesis, antisense oligonucleiotides (ASODNs) were used to downregulate ASGPR1 level in HepG2.2.15 cells. By using the MFOLD web server and BLAST searches, five ASODNs theoretically targeting ASGPR1 were selected. After 72 h post-transfection, HBV-DNA level in cell medium were examined by real-time polymerase chain reaction (PCR). Hepatitis B surface antigen (HBsAg) and Hepatitis B e antigen (HBeAg) were detected using enzyme-linked immunosorbent assay (ELISA). ASGPR1 mRNA and protein level were measured by semi-quantitative reverse transcriptase (RT)-PCR and Western blot analysis respectively. The results showed that ASODN2 significantly downregulated ASGPR1 level. It also reduced HBV-DNA, HBsAg and HBeAg level in cell medium as observed with lamivudine. In contrast, the sense sequence and scrambled sequence of ASODN2 had no effect on ASGPR1 and HBV markers in HepG2.2.15 cells. This indicated that ASODN2 could specifically reduce HBV replication in vitro. Additionally, cell proliferation and apoptosis assay suggested that downregulation of ASGPR1 did not affect cell viability. We, therefore, proposed that ASODNs targeted against ASGPR1 could block HBV replication without the influence of other changes, and ASGPR1 could be targeted for anti-HBV drug development.

Differentially expressed cellular genes following HBV: potential targets of anti-HBV drugs?

The aim of the study was to screen for cellular genes that are differentially expressed following hepatitis B virus (HBV) infection, in an attempt to identify potential targets of anti-HBV drugs. An oligonucleotide microarray containing 231 virus-infection-associated genes was prepared. Differential gene expression in HepG2.2.15 cells compared to control with HepG2 cells was analysed by this in-house microarray. The change in gene expression in HepG2.2.15 cells treated by lamivudine on days 4 and 8 after exposure was also studied. Semi-quantitative reverse transcriptase-polymerase chain reaction (RT-PCR) was used to comfirm the differentially expressed genes induced by HBV and lamivudine. There were 31 upregulated and four downregulated genes in HepG2.2.15 cells compared with the HepG2 control cells. Eleven genes were consistently altered by lamivudine at both time points. Of the 31 genes that were upregulated in HepG2.2.15 cells, there were seven genes which were downregulated by lamivudine. Of the four downregulated genes, there was one gene which was upregulated by lamivudine. Of the differentially expressed genes induced by HBV and lamivudine, the expression of five genes was confirmed by semi-quantitative RT-PCR. These results shed new light on the effects of HBV and lamivudine on cellular gene expression. Differentially expressed genes induced by HBV and lamivudine could potentially become new anti-HBV drug targets in novel therapies.

Effective siRNA targets screening for human telomerase reverse transcriptase

AIM: To study the inhibitory effects of siRNAs targeting different hTERT sequences and to screen the effective siRNA sequence.

METHODS: Five double-stranded siRNAs targeting coding and non-coding regions of hTERT gene were designed and synthesized by T7 transcription system in vitro. siRNA4 sequence was screened by full length gene targeting technique and the rest of the siRNA sequences were selected randomly. After being purified by ethanol precipitation, the siRNAs were transfected to the human hepatocellular carcinoma cell (HepG2) by Lipofectamine^TM2000. At 48-72 h after siRNAs transfection, MTT assay, RT-PCR and Western-blot were applied to evaluate the effects of siRNAs on cell growth, mRNA and protein expression level of hTERT gene, respectively.

RESULTS: Compared to the control cells, the cells treated with the five double-stranded siRNAs exhibited different degrees of inhibition of cell proliferation in a dose-dependent manner. siRNA2 and siRNA4, exhibited obvious effects of inhibiting hTERT mRNA and protein expression in HepG2 cells.

CONCLUSION: siRNAs targeting different hTERT sequences have significantly various inhibitory effects on hTERT gene expression. The siRNA sequence screened by full length gene targeting technique has comparable inhibitory effect with the rest siRNA sequences screened by random selection, suggesting that siRNAs and antisense oligonucleic acids may have the same effective target sites. Compared with chemical synthesis method, synthesizing double-stranded siRNA by T7 transcription system in vitro is a rapid, simple, and inexpensive method suitable for screening high-effect siRNA targeting site for specific gene.

HBx protein of hepatitis B virus (HBV) can form complex with mitochondrial HSP60 and HSP70

HBx, a transcriptional transactivating protein of hepatitis B virus (HBV), is required for viral infection and has been implicated in virus-mediated liver oncogenesis. However, the molecular mechanism for its influence on cell remains largely unknown. It was proved that HBx need the help of host cell proteins to exert its function by binding to them. During purifying of GSTX (fusion protein of GST and HBx) expressed in E. coli, we found that it can bind specifically with GrpE (HSP60) and DnaK (HSP70) of E. coli while GST cannot. Using GST pull-down, two-dimensional gel electrophoresis and mass spectrum, we found that GSTX can also bind to human mitochondrial HSP60 and HSP70, which are homologues of GrpE and DnaK. These interactions between HBx and mitochondrial HSP60 and HSP70 are supported by the result of co-immunoprecipitation experiment. It means that HBx can form complex with E. coli and human HSP60 and HSP70. The implication of HBx, HSP60 and HSP70 complex in molecular mechanism of virus infection is discussed.

HBVPathDB: A database of HBV infection-related molecular interaction network

AIM: To describe molecules or genes interaction between hepatitis B viruses (HBV) and host, for understanding how virus’ and host’s genes and molecules are networked to form a biological system and for perceiving mechanism of HBV infection.

METHODS: The knowledge of HBV infection-related reactions was organized into various kinds of pathways with carefully drawn graphs in HBVPathDB. Pathway information is stored with relational database management system (DBMS), which is currently the most efficient way to manage large amounts of data and query is implemented with powerful Structured Query Language (SQL). The search engine is written using Personal Home Page (PHP) with SQL embedded and web retrieval interface is developed for searching with Hypertext Markup Language (HTML).

RESULTS: We present the first version of HBVPathDB, which is a HBV infection-related molecular interaction network database composed of 306 pathways with 1050 molecules involved. With carefully drawn graphs, pathway information stored in HBVPathDB can be browsed in an intuitive way. We develop an easy-to-use interface for flexible accesses to the details of database. Convenient software is implemented to query and browse the pathway information of HBVPathDB. Four search page layout options-category search, gene search, description search, unitized search-are supported by the search engine of the database. The database is freely available at http://www.bio-inf.net/HBVPathDB/HBV/.

CONCLUSION: The conventional perspective HBVPathDB have already contained a considerable amount of pathway information with HBV infection related, which is suitable for in-depth analysis of molecular interaction network of virus and host. HBVPathDB integrates pathway data-sets with convenient software for query, browsing, visualization, that provides users more opportunity to identify regulatory key molecules as potential drug targets and to explore the possible mechanism of HBV infection based on gene expression datasets.

Screening of specific antigens for SARS clinical diagnosis using a protein microarray

In this study several SARS-CoV structural proteins and fragments were expressed in E. coli as GST or TRX fusion proteins. They were fabricated on a microarray and tested with sera from SARS patients. Antigenic screening indicated that recombinant GST-N2 fusion protein, the carboxy-terminus 213aa-423aa of N protein, was strongest positive and weakest non-specific compared with others. An indirect antibody ELISA method was developed and clinical positive and negative sera for their antibodies against GST-N2 fusion protein were assayed. 311 out of the 442 sera from clinical SARS inpatients, as well as 229 out of 302 sera from convalescent patients gave positive reactivities; positive rates were 70.4% and 75.8% respectively. Sera from a total of 2726 non-SARS patients and healthy individuals were tested and the false positive rate was only 0.07%. When the sensitivity control sample was diluted 1 : 64, it yielded OD values above the cutoff value. Reported data showed that this was a relatively high degree of sensitivity and specificity for SARS-CoV antibody testing. The data indicate that GST-N2 fusion protein, which was screened by protein microarray, may be a valuable diagnostic antigen for the development of serological assays for SARS. In addition, protein microarray assay presents a higher positive rate and sensitivity (86.1% and 1 : 200) compared with the traditional ELISA screening method, and could provide a rapid, parallel and high-throughput antigen screening platform.