Si Lok

Increased exonic de novo mutation rate in individuals with schizophrenia

Schizophrenia is a severe psychiatric disorder that profoundly affects cognitive, behavioral and emotional processes. The wide spectrum of symptoms and clinical variability in schizophrenia suggest a complex genetic etiology, which is consistent with the numerous loci thus far identified by linkage, copy number variation and association studies. Although schizophrenia heritability may be as high as ∼80%, the genes responsible for much of this heritability remain to be identified. Here we sequenced the exomes of 14 schizophrenia probands and their parents. We identified 15 de novo mutations (DNMs) in eight probands, which is significantly more than expected considering the previously reported DNM rate. In addition, 4 of the 15 identified DNMs are nonsense mutations, which is more than what is expected by chance. Our study supports the notion that DNMs may account for some of the heritability reported for schizophrenia while providing a list of genes possibly involved in disease pathogenesis.

Recoding RNA editing of AZIN1 predisposes to hepatocellular carcinoma

A better understanding of human hepatocellular carcinoma (HCC) pathogenesis at the molecular level will facilitate the discovery of tumor-initiating events. Transcriptome sequencing revealed that adenosine-to-inosine (A→I) RNA editing of AZIN1 (encoding antizyme inhibitor 1) is increased in HCC specimens. A→I editing of AZIN1 transcripts, specifically regulated by ADAR1 (encoding adenosine deaminase acting on RNA-1), results in a serine-to-glycine substitution at residue 367 of AZIN1, located in β-strand 15 (β15) and predicted to cause a conformational change, induced a cytoplasmic-to-nuclear translocation and conferred gain-of-function phenotypes that were manifested by augmented tumor-initiating potential and more aggressive behavior. Compared with wild-type AZIN1 protein, the edited form has a stronger affinity to antizyme, and the resultant higher AZIN1 protein stability promotes cell proliferation through the neutralization of antizyme-mediated degradation of ornithine decarboxylase (ODC) and cyclin D1 (CCND1). Collectively, A→I RNA editing of AZIN1 may be a potential driver in the pathogenesis of human cancers, particularly HCC.

Complete Sequencing of pNDM-HK Encoding NDM-1 Carbapenemase from a Multidrug-Resistant Escherichia coli Strain Isolated in Hong Kong

Background The emergence of plasmid-mediated carbapenemases, such as NDM-1 in Enterobacteriaceae is a major public health issue. Since they mediate resistance to virtually all β-lactam antibiotics and there is often co-resistance to other antibiotic classes, the therapeutic options for infections caused by these organisms are very limited. Methodology We characterized the first NDM-1 producing E. coli isolate recovered in Hong Kong. The plasmid encoding the metallo-β-lactamase gene was sequenced. Principal Findings The plasmid, pNDM-HK readily transferred to E. coli J53 at high frequencies. It belongs to the broad host range IncL/M incompatibility group and is 88803 bp in size. Sequence alignment showed that pNDM-HK has a 55 kb backbone which shared 97% homology with pEL60 originating from the plant pathogen, Erwina amylovora in Lebanon and a 28.9 kb variable region. The plasmid backbone includes the mucAB genes mediating ultraviolet light resistance. The 28.9 kb region has a composite transposon-like structure which includes intact or truncated genes associated with resistance to β-lactams (bla TEM-1, bla NDM-1, Δbla DHA-1), aminoglycosides (aacC2, armA), sulphonamides (sul1) and macrolides (mel, mph2). It also harbors the following mobile elements: IS26, ISCR1, tnpU, tnpAcp2, tnpD, ΔtnpATn1 and insL. Certain blocks within the 28.9 kb variable region had homology with the corresponding sequences in the widely disseminated plasmids, pCTX-M3, pMUR050 and pKP048 originating from bacteria in Poland in 1996, in Spain in 2002 and in China in 2006, respectively. Significance The genetic support of NDM-1 gene suggests that it has evolved through complex pathways. The association with broad host range plasmid and multiple mobile genetic elements explain its observed horizontal mobility in multiple bacterial taxa.

Whole-genome bisulfite sequencing of multiple individuals reveals complementary roles of promoter and gene body methylation in transcriptional regulation

BackgroundDNA methylation is an important type of epigenetic modification involved in gene regulation. Although strong DNA methylation at promoters is widely recognized to be associated with transcriptional repression, many aspects of DNA methylation remain not fully understood, including the quantitative relationships between DNA methylation and expression levels, and the individual roles of promoter and gene body methylation.ResultsHere we present an integrated analysis of whole-genome bisulfite sequencing and RNA sequencing data from human samples and cell lines. We find that while promoter methylation inversely correlates with gene expression as generally observed, the repressive effect is clear only on genes with a very high DNA methylation level. By means of statistical modeling, we find that DNA methylation is indicative of the expression class of a gene in general, but gene body methylation is a better indicator than promoter methylation. These findings are general in that a model constructed from a sample or cell line could accurately fit the unseen data from another. We further find that promoter and gene body methylation have minimal redundancy, and either one is sufficient to signify low expression. Finally, we obtain increased modeling power by integrating histone modification data with the DNA methylation data, showing that neither type of information fully subsumes the other.ConclusionOur results suggest that DNA methylation outside promoters also plays critical roles in gene regulation. Future studies on gene regulatory mechanisms and disease-associated differential methylation should pay more attention to DNA methylation at gene bodies and other non-promoter regions.

Biased Diversity Metrics Revealed by Bacterial 16S Pyrotags Derived from Different Primer Sets

In recent years, PCR-based pyrosequencing of 16S rRNA genes has continuously increased our understanding of complex microbial communities in various environments of the Earth. However, there is always concern on the potential biases of diversity determination using different 16S rRNA gene primer sets and covered regions. Here, we first report how bacterial 16S rRNA gene pyrotags derived from a series of different primer sets resulted in the biased diversity metrics. In total, 14 types of pyrotags were obtained from two-end pyrosequencing of 7 amplicon pools generated by 7 primer sets paired by 1 of 4 forward primers (V1F, V3F, V5F, and V7F) and 1 of 4 reverse primers (V2R, V4R, V6R, and V9R), respectively. The results revealed that: i) the activated sludge exhibited a large bacterial diversity that represented a broad range of bacterial populations and served as a good sample in this methodology research; ii) diversity metrics highly depended on the selected primer sets and covered regions; iii) paired pyrotags obtained from two-end pyrosequencing of each short amplicon displayed different diversity metrics; iv) relative abundance analysis indicated the sequencing depth affected the determination of rare bacteria but not abundant bacteria; v) the primer set of V1F and V2R significantly underestimated the diversity of activated sludge; and vi) the primer set of V3F and V4R was highly recommended for future studies due to its advantages over other primer sets. All of these findings highlight the significance of this methodology research and offer a valuable reference for peer researchers working on microbial diversity determination.

Ammonia-oxidizing archaea and ammonia-oxidizing bacteria in six full-scale wastewater treatment bioreactors

In this study, dideoxy sequencing and 454 high-throughput sequencing were used to analyze diversities of the ammonia monooxygenase (amoA) genes and the 16S rRNA genes of ammonia-oxidizing archaea (AOA) and ammonia-oxidizing bacteria (AOB) in six municipal wastewater treatment plants. The results showed that AOB amoA genes were quite diverse in different wastewater treatment plants while the 16S rRNA genes were relatively conserved. Based on the observed complexity of amoA and 16S rRNA genes, most of the AOB can be assigned to the Nitrosomonas genus, with Nitrosomonas ureae, Nitrosomonas oligotropha, Nitrosomonas marina, and Nitrosomonas aestuarii being the four most dominant species. From the sequences of the AOA amoA genes, most AOA observed in this study belong to the CGI.1b group, i.e., the soil lineage. The AOB amoA and 16S rRNA genes were quantified by quantitative PCR and 454 high-throughput pyrosequencing, respectively. Although the results from the two approaches show some disconcordance, they both indicated that the abundance of AOB in activated sludge was very low.

Rab25 Is a Tumor Suppressor Gene with Antiangiogenic and Anti-Invasive Activities in Esophageal Squamous Cell Carcinoma

Esophageal squamous cell carcinoma (ESCC), the major histologic subtype of esophageal cancer, is a devastating disease characterized by distinctly high incidences and mortality rates. However, there remains limited understanding of molecular events leading to development and progression of the disease, which are of paramount importance to defining biomarkers for diagnosis, prognosis, and personalized treatment. By high-throughout transcriptome sequence profiling of nontumor and ESCC clinical samples, we identified a subset of significantly differentially expressed genes involved in integrin signaling. The Rab25 gene implicated in endocytic recycling of integrins was the only gene in this group significantly downregulated, and its downregulation was confirmed as a frequent event in a second larger cohort of ESCC tumor specimens by quantitative real-time PCR and immunohistochemical analyses. Reduced expression of Rab25 correlated with decreased overall survival and was also documented in ESCC cell lines compared with pooled normal tissues. Demethylation treatment and bisulfite genomic sequencing analyses revealed that downregulation of Rab25 expression in both ESCC cell lines and clinical samples was associated with promoter hypermethylation. Functional studies using lentiviral-based overexpression and suppression systems lent direct support of Rab25 to function as an important tumor suppressor with both anti-invasive and -angiogenic abilities, through a deregulated FAK-Raf-MEK1/2-ERK signaling pathway. Further characterization of Rab25 may provide a prognostic biomarker for ESCC outcome prediction and a novel therapeutic target in ESCC treatment.

Molecular Characterization of the 2011 Hong Kong Scarlet Fever Outbreak

A scarlet fever outbreak occurred in Hong Kong in 2011. The majority of cases resulted in the isolation of Streptococcus pyogenes emm12 with multiple antibiotic resistances. Phylogenetic analysis of 22 emm12 scarlet fever outbreak isolates, 7 temporally and geographically matched emm12 non-scarlet fever isolates, and 18 emm12 strains isolated during 2005-2010 indicated the outbreak was multiclonal. Genome sequencing of 2 nonclonal scarlet fever isolates (HKU16 and HKU30), coupled with diagnostic polymerase chain reaction assays, identified 2 mobile genetic elements distributed across the major lineages: a 64.9-kb integrative and conjugative element encoding tetracycline and macrolide resistance and a 46.4-kb prophage encoding superantigens SSA and SpeC and the DNase Spd1. Phenotypic comparison of HKU16 and HKU30 with the S. pyogenes M1T1 strain 5448 revealed that HKU16 displays increased adherence to HEp-2 human epithelial cells, whereas HKU16, HKU30, and 5448 exhibit equivalent resistance to neutrophils and virulence in a humanized plasminogen murine model. However, in contrast to M1T1, the virulence of HKU16 and HKU30 was not associated with covRS mutation. The multiclonal nature of the emm12 scarlet fever isolates suggests that factors such as mobile genetic elements, environmental factors, and host immune status may have contributed to the 2011 scarlet fever outbreak.

Genomic Sequencing and Comparative Analysis of Epstein-Barr Virus Genome Isolated from Primary Nasopharyngeal Carcinoma Biopsy

Whether certain Epstein-Barr virus (EBV) strains are associated with pathogenesis of nasopharyngeal carcinoma (NPC) is still an unresolved question. In the present study, EBV genome contained in a primary NPC tumor biopsy was amplified by Polymerase Chain Reaction (PCR), and sequenced using next-generation (Illumina) and conventional dideoxy-DNA sequencing. The EBV genome, designated HKNPC1 (Genbank accession number JQ009376) is a type 1 EBV of approximately 171.5 kb. The virus appears to be a uniform strain in line with accepted monoclonal nature of EBV in NPC but is heterogeneous at 172 nucleotide positions. Phylogenetic analysis with the four published EBV strains, B95-8, AG876, GD1, and GD2, indicated HKNPC1 was more closely related to the Chinese NPC patient-derived strains, GD1 and GD2. HKNPC1 contains 1,589 single nucleotide variations (SNVs) and 132 insertions or deletions (indels) in comparison to the reference EBV sequence (accession number NC007605). When compared to AG876, a strain derived from Ghanaian Burkitts lymphoma, we found 322 SNVs, of which 76 were non-synonymous SNVs and were shared amongst the Chinese GD1, GD2 and HKNPC1 isolates. We observed 88 non-synonymous SNVs shared only by HKNPC1 and GD2, the only other NPC tumor-derived strain reported thus far. Non-synonymous SNVs were mainly found in the latent, tegument and glycoprotein genes. The same point mutations were found in glycoprotein (BLLF1 and BALF4) genes of GD1, GD2 and HKNPC1 strains and might affect cell type specific binding. Variations in LMP1 and EBNA3B epitopes and mutations in Cp (11404 C>T) and Qp (50134 G>C) found in GD1, GD2 and HKNPC1 could potentially affect CD8+ T cell recognition and latent gene expression pattern in NPC, respectively. In conclusion, we showed that whole genome sequencing of EBV in NPC may facilitate discovery of previously unknown variations of pathogenic significance.

Identification of PTK6, via RNA Sequencing Analysis, as a Suppressor of Esophageal Squamous Cell Carcinoma

BACKGROUND & AIMS Esophageal squamous cell carcinoma (ESCC) is the most commonly observed histologic subtype of esophageal cancer. ESCC is believed to develop via accumulation of numerous genetic alterations, including inactivation of tumor suppressor genes and activation of oncogenes. We searched for transcripts that were altered in human ESCC samples compared with nontumor tissues. METHODS We performed integrative transcriptome sequencing (RNA-Seq) analysis using ESCC samples from 3 patients and adjacent nontumor tissues to identify transcripts that were altered in ESCC tissue. We performed molecular and functional studies of the transcripts identified and investigated the mechanisms of alteration. RESULTS We identified protein tyrosine kinase 6 (PTK6) as a transcript that was significantly down-regulated in ESCC tissues and cell lines compared with nontumor tissues or immortalized normal esophageal cell lines. The promoter of the PTK6 gene was inactivated in ESCC tissues at least in part via hypermethylation and histone deacetylation. Knockdown of PTK6 in KYSE30 ESCC cells using small hairpin RNAs increased their ability to form foci, migrate, and invade extracellular matrix in culture and form tumors in nude mice. Overexpression of PTK6 in these cells reduced their proliferation in culture and tumor formation in mice. PTK6 reduced phosphorylation of Akt and glycogen synthase kinase (GSK)3β, leading to activation of β-catenin. CONCLUSIONS PTK6 was identified as a transcript that is down-regulated in human ESCC tissues via epigenetic modification at the PTK6 locus. Its product appears to regulate cell proliferation by reducing phosphorylation of Akt and GSK3β, leading to activation of β-catenin. Reduced levels of PTK6 promote growth of xenograft tumors in mice; it might be developed as a marker of ESCC.