Molecular and functional characterization of tumor-induced factor (TIF): Hamster homolog of CXCL3 (GROγ) displays tumor suppressive activity.

Previously our lab has created a mouse ovarian xenograft model of copy number variation (CNV)-mediated G protein-coupled receptor (GPCR) MAS-driven tumorigenesis, and RNA profiling identified a putative chemokine tumor-induced factor (Tif). Sequence analysis and chemotactic study suggested that Tif was likely to be a hamster homolog of human GROγ (CXCL3) [IJC 125 (2009): 1316-1327]. In the present study, we report the molecular and functional characterization of the Tif gene. Genomic study of CHO-K1 cells indicated that Tif gene consisted of 4 exons, characterized with an antisense B1 element which is embedded in the fourth exon. Two Tif transcripts were identified which shared identical sequences except that a string of 71-nt derived from the antisense B1 element was deficient in the shorter transcript. Of interests, B1-like RNA ladder was detected in xenografts. Functional studies showed that TIF induced chemotaxis and neovessel formation. Pharmacological studies suggested that TIF activated Gi-coupled CXCR2 and induced both calcium mobilization and ERK1/2 phosphorylation, and suppressed forskolin-stimulated cAMP accumulation. In addition, secreted matured TIF functioned as an autocrine factor and promoted anchorage-independent growth. Unexpectedly, TIF delayed the onset of tumor formation, possibly via suppressing proliferation of stromal fibroblasts. However, TIF did not exert any inhibitory effect on tumor growth. Potentially, TIF could be used for preventing cancer relapse.

Receptor tyrosine kinase-like orphan receptor serves as a potential target in cancer immunotherapy.

Receptor tyrosine kinase-like orphan receptor (ROR), consisting of ROR1 and ROR2, is a conserved family of receptor tyrosine kinase superfamily that plays crucial roles during embryonic development with limited expression in adult normal tissues. However, it is overexpressed in a range of hematological malignancies and solid tumors and functions in cellular processes including cell survival, polarity, and migration, serving as a potential target in cancer immunotherapy. This review summarizes the expression and structure of ROR in developmental morphogenesis and its function in cancers associated with Wnt5a signaling and highlights the cancer immunotherapy strategies targeting ROR.

CDR3 sequences in IgA nephropathy are shorter and exhibit reduced diversity.

Immunoglobulin (Ig) A nephropathy (IgAN) is the most common glomerulonephritis, which is characterized by the deposition of IgA antibody in the glomerulus. Systematic dissection of immune composition may contribute to a better understanding of the alternations in the immune system in IgAN. To this end, here we applied immune repertoire sequencing technology for parallel analysis of the complementary determining region 3 (CDR3) of all B cell receptors, including all five antibody subtypes (IgA, IgG, IgM, IgE and IgD), in 13 patients with IgAN and 7 healthy individuals. A significant decrease in CDR3 length was observed in the IgAN group. In particular, the JH6 family was significantly increased in IgAN. Amino acid usage was also altered in IgAN. Shannon, Simpson, Gini and Diversity 50 indices also revealed significant differences in the diversity of IgG, IgM and IgA antibodies as compared with controls. The proportions of IgA and IgG were increased, whereas IgM was decreased in IgAN. Moreover, a greater number of CDR3 sequences was shared between patients with IgAN. These findings suggest that the BCR immune repertoire is dramatically altered in IgAN, as characterized by shortened CDR3 length, as well as decreased overall diversity of CDR3.

An efficient and precise method for generating knockout cell lines based on CRISPR-Cas9 system.

Although the efficiency and versatility of CRISPR-Cas9 system has been greatly improved over conventional genome editing methods such as zinc finger or TALEN, it is still time-consuming and labor-intensive for screening knockout/knock-in cell clones due to differences of the targeted location or efficacies of guide RNAs (gRNAs). Here, we adapted a targeted knock-in strategy with CRISPR-Cas9 system and characterized the efficiency for generating single or double knockout cell lines. Specifically, a homology-arm based donor cassette consisting of genes encoding a fluorescence protein and antibiotic selection marker driven by a constitutive promoter was co-transfected with a gRNA expressing unit. Based on FACS sorting and antibiotic drug selection, positive cell clones were confirmed by genotyping and at the protein expression level. The results indicated that more than 70% of analyzed clones identified by cell sorting and selection were successfully targeted in both single and double knockout experiments. The procedure takes less than three weeks to obtain knockout cell lines. We believe that this methodology could be applicable and versatile in generating knockout cell clones with high efficiency in most cell lines.

Identification and characterization of a novel CXC chemokine in xenograft tumor induced by mas-overexpressing cells.

Overexpressions of G protein-coupled receptor (GPCR) with elevated downstream signaling events have been reported in various tumors. However, the cellular mechanism that GPCR overexpression leads to tumor formation is largely unknown. The orphan GPCR mas was originally isolated from a human epidermoid carcinoma. In vivo studies of mas-overexpressing cells suggested that xenograft tumor formation was positively correlated with the levels of mas expression. Histochemical analysis indicated that xenograft tumor consisted of mas-transfected and stromal cells. Biochemical analyses revealed that cells overexpressing mas exhibited significantly increased anchorage-independent growth, whereas there was no significant difference in cell proliferation in comparison with empty vector-transfected control cells. Expression profiling using mRNA differential display and Northern analysis indicated an elevated expression of GRO and a novel CXC chemokines, tumor-induced factor (TIF), in mas-transfected cells and xenograft tumor. Bacterially expressed recombinant TIF was found to act as a neutrophil chemoattractant in a chemotactic assay. These results suggest that mas overexpression enables anchorage-independent growth of transformed cells, and interplays of secreted chemokines with stromal cells modulate xenograft tumor formation. Importantly, a novel CXC chemokine, TIF, was identified in the xenograft tumor tissues.

Integrated sequencing of exome and mRNA of large-sized single cells.

Current approaches of single cell DNA-RNA integrated sequencing are difficult to call SNPs, because a large amount of DNA and RNA is lost during DNA-RNA separation. Here, we performed simultaneous single-cell exome and transcriptome sequencing on individual mouse oocytes. Using microinjection, we kept the nuclei intact to avoid DNA loss, while retaining the cytoplasm inside the cell membrane, to maximize the amount of DNA and RNA captured from the single cell. We then conducted exome-sequencing on the isolated nuclei and mRNA-sequencing on the enucleated cytoplasm. For single oocytes, exome-seq can cover up to 92% of exome region with an average sequencing depth of 10+, while mRNA-sequencing reveals more than 10,000 expressed genes in enucleated cytoplasm, with similar performance for intact oocytes. This approach provides unprecedented opportunities to study DNA-RNA regulation, such as RNA editing at single nucleotide level in oocytes. In future, this method can also be applied to other large cells, including neurons, large dendritic cells and large tumour cells for integrated exome and transcriptome sequencing.

Long range haplotyping of paired-homologous chromosomes by single-chromosome sequencing of a single cell.

The longest possible haplotype is chromosome haplotype that is a set of co-inherited alleles occurred on a single strand chromosome inherited from one parent. Standard whole-genome shotgun sequencing technologies are limited by the inability to independently study the haplotype of homologous chromosomes due to the short-reads sequencing strategy and disturbance of homologue chromosomes. Here, we investigated several types of chromosomal abnormalities by a dilution-based method to separate an intact copy of homologous chromosome from human metaphase cells, and then single chromosomes were independently amplified by whole-genome amplification methods, converted into barcoded sequencing libraries, and sequenced in multiplexed pools by Illumina sequencers. We analyzed single chromosome derived from single metaphase cells of one patient with balanced chromosomal translocation t(3;5)(q24;q13), one patient with (47, XXY) karyotype and one with (47, XY, 21+) Down syndrome. We determined the translocation region of chromosomes in patient with t(3;5)(q24;q13) balanced chromosomal translocation by shallow whole-genome sequencing, which is helpful to pinpoint the chromosomal break point. We showed that SCS can physically separate and independently sequence three copies of chromosome 21 of Down syndrome patient. SCS has potential applications in personal genomics, single-cell genomics, and clinical diagnosis, particularly in revealing chromosomal level of genetic diseases.

Statistical analysis of mutant allele frequency level of circulating cell-free DNA and blood cells in healthy individuals.

Cell-free DNA (cfDNA) in plasma has emerged as a potential important biomarker in clinical diagnostics, particularly in cancer. However, somatic mutations are also commonly found in healthy individuals, which interfere with the effectiveness for cancer diagnostics. This study examined the background somatic mutations in white blood cells (WBC) and cfDNA in healthy controls based on sequencing data from 821 non-cancer individuals and several cancer samples with the aim of understanding the patterns of mutations detected in cfDNA. We determined the mutation allele frequencies in both WBC and cfDNA using a panel of 50 cancer-associated genes that covers 20 K-nucleotide region and ultra-deep sequencing with average depth >40000-fold. Our results showed that most of the mutations in cfDNA were highly correlated to WBC. We also observed that the NPM1 gene was the most frequently mutated gene in both WBC and cfDNA. Our study highlighted the importance of sequencing both cfDNA and WBC to improve the sensitivity and accuracy for calling cancer-related mutations from circulating tumour DNA, and shedded light on developing a strategy for early cancer diagnosis by cfDNA sequencing.

PGWD: Integrating Personal Genome for Warfarin Dosing.

Warfarin is a drug normally used in the prevention of thrombosis and the formation of blood clots. The dosage of warfarin is strongly affected by genetic variants of CYP2C9 and VKORC1 genes. Current technologies for detecting the variants of these genes are mainly based on real-time PCR. In recent years, due to the rapidly dropping cost of whole genome sequencing and genotyping, more and more people get their whole genome sequenced or genotyped. However, current software for warfarin dosing prediction is based on low-throughput genetic information from either real-time PCR or melting curve methods. There is no bioinformatics tool available that can take the high-throughput genome sequencing data as input and determine the accurate dosage of warfarin. Here, we present PGWD, a web tool that analyzes personal genome sequencing data and integrates with clinical information for warfarin dosing.

T Cell Repertoire Diversity Is Decreased in Type 1 Diabetes Patients.

Type 1 diabetes mellitus (T1D) is an immune-mediated disease. The autoreactive T cells in T1D patients attack and destroy their own pancreatic cells. In order to systematically investigate the potential autoreactive T cell receptors (TCRs), we used a high-throughput immune repertoire sequencing technique to profile the spectrum of TCRs in individual T1D patients and controls. We sequenced the T cell repertoire of nine T1D patients, four type 2 diabetes (T2D) patients, and six nondiabetic controls. The diversity of the T cell repertoire in T1D patients was significantly decreased in comparison with T2D patients (P=7.0E-08 for CD4+ T cells, P=1.4E-04 for CD8+ T cells) and nondiabetic controls (P=2.7E-09 for CD4+ T cells, P=7.6E-06 for CD8+ T cells). Moreover, T1D patients had significantly more highly-expanded T cell clones than T2D patients (P=5.2E-06 for CD4+ T cells, P=1.9E-07 for CD8+ T cells) and nondiabetic controls (P=1.7E-07 for CD4+ T cells, P=3.3E-03 for CD8+ T cells). Furthermore, we identified a group of highly-expanded T cell receptor clones that are shared by more than two T1D patients. Although further validation in larger cohorts is needed, our data suggest that T cell receptor diversity measurements may become a valuable tool in investigating diabetes, such as using the diversity as an index to distinguish different types of diabetes.

PGWD: integrating personal genome for warfarin dosing.

Warfarin is a drug normally used in the prevention of thrombosis and the formation of blood clots. The dosage of warfarin is strongly affected by genetic variants of CYP2C9 and VKORC1 genes. Current technologies for detecting the variants of these genes are mainly based on real-time PCR. In recent years, due to the rapidly dropping cost of whole genome sequencing and genotyping, more and more people get their whole genome sequenced or genotyped. However, current software for warfarin dosing prediction is based on low-throughput genetic information either from real-time PCR, or melting curve methods. There is no bioinformatics tool available that can take the high-throughput genome sequencing data as input and determine the accurate dosage of warfarin. Here, we present PGWD, a web tool that analyzes personal genome sequencing data, and integrates with clinical information for warfarin dosing.

Virtual Pharmacist: A Platform for Pharmacogenomics.

We present Virtual Pharmacist, a web-based platform that takes common types of high-throughput data, namely microarray SNP genotyping data, FASTQ and Variant Call Format (VCF) files as inputs, and reports potential drug responses in terms of efficacy, dosage and toxicity at one glance. Batch submission facilitates multivariate analysis or data mining of targeted groups. Individual analysis consists of a report that is readily comprehensible to patients and practioners who have basic knowledge in pharmacology, a table that summarizes variants and potential affected drug response according to the US Food and Drug Administration pharmacogenomic biomarker labeled drug list and PharmGKB, and visualization of a gene-drug-target network. Group analysis provides the distribution of the variants and potential affected drug response of a target group, a sample-gene variant count table, and a sample-drug count table. Our analysis of genomes from the 1000 Genome Project underlines the potentially differential drug responses among different human populations. Even within the same population, the findings from Watson's genome highlight the importance of personalized medicine. Virtual Pharmacist can be accessed freely at http://www.sustc-genome.org.cn/vp or installed as a local web server. The codes and documentation are available at the GitHub repository (https://github.com/VirtualPharmacist/vp). Administrators can download the source codes to customize access settings for further development.

Comprehensive analysis of the T-cell receptor beta chain gene in rhesus monkey by high throughput sequencing.

Profiling immune repertoires by high throughput sequencing enhances our understanding of immune system complexity and immune-related diseases in humans. Previously, cloning and Sanger sequencing identified limited numbers of T cell receptor (TCR) nucleotide sequences in rhesus monkeys, thus their full immune repertoire is unknown. We applied multiplex PCR and Illumina high throughput sequencing to study the TCRß of rhesus monkeys. We identified 1.26 million TCRß sequences corresponding to 643,570 unique TCRß sequences and 270,557 unique complementarity-determining region 3 (CDR3) gene sequences. Precise measurements of CDR3 length distribution, CDR3 amino acid distribution, length distribution of N nucleotide of junctional region, and TCRV and TCRJ gene usage preferences were performed. A comprehensive profile of rhesus monkey immune repertoire might aid human infectious disease studies using rhesus monkeys.

Quantitative assessment of single-cell whole genome amplification methods for detecting copy number variation using hippocampal neurons.

Single-cell genomic analysis has grown rapidly in recent years and finds widespread applications in various fields of biology, including cancer biology, development, immunology, pre-implantation genetic diagnosis, and neurobiology. To date, the amplification bias, amplification uniformity and reproducibility of the three major single cell whole genome amplification methods (GenomePlex WGA4, MDA and MALBAC) have not been systematically investigated using mammalian cells. In this study, we amplified genomic DNA from individual hippocampal neurons using three single-cell DNA amplification methods, and sequenced them at shallow depth. We then systematically evaluated the GC-bias, reproducibility, and copy number variations among individual neurons. Our results showed that single-cell genome sequencing results obtained from the MALBAC and WGA4 methods are highly reproducible and have a high success rate. The MALBAC displays significant biases towards high GC content. We then attempted to correct the GC bias issue by developing a bioinformatics pipeline, which allows us to call CNVs in single cell sequencing data, and chromosome level and sub-chromosomal level CNVs among individual neurons can be detected. We also proposed a metric to determine the CNV detection limits. Overall, MALBAC and WGA4 have better performance than MDA in detecting CNVs.

Determining common insertion sites based on retroviral insertion distribution across tumors.

A CIS (common insertion site) indicates a genome region that is hit more frequently by retroviral insertions than expected by chance. Such a region is strongly related to cancer gene loci, which leads to the detection of cancer genes. An algorithm for detecting CISs should satisfy the following: (1) it does not require any prior knowledge of underlying insertion distribution; (2) it can resolve the insertion biases caused by hotspots; (3) it can detect CISs of any biological width; (4) it can identify noises resulting from statistic mistakes and non-CIS insertions; and (5) it can identify the widths of CISs as accurately as possible. We develop a method to resolve these difficulties. We verify a region's significance from two perspectives: distribution width and distribution depth. The former indicates how many insertions in a region while the latter evaluates the insertion distribution across the tumors in a region. We compare our method with kernel density estimation and sliding window on the simulated data, showing that our method not only identifies cancer-related insertions effectively, but also filters noises correctly. The experiments on the real data show that taking insertion distribution into account can highlight significant CISs. We detect 53 novel CISs, some of which have been proven correct by the biological literature.