Identification and expression profiling of microRNAs in the brain, liver and gonads of marine medaka (Oryzias melastigma) and in response to hypoxia.

The marine medaka (Oryzias melastigma) has been increasingly used as a fish model for detecting environmental stresses and chemical contaminants in the marine environment. Recent mammalian studies have shown that environmental stresses can alter the expression profiles of microRNAs (miRNAs), leading to transgenerational effects. Here, we use high-throughput Illumina RNA sequencing (RNA-Seq) for miRNA transcriptome analysis of brain, liver, and gonads from sexually mature male and female marine medaka. A total of 128,883,806 filtered sequence reads were generated from six small RNA libraries, identifying a total of 2,125,663 non-redundant sequences. These sequences were aligned and annotated to known animal miRNAs (miRBase) using the BLAST method. A total of 223 distinct miRNA types were identified, with the greatest number expressed in brain tissue. Our data suggested that 55 miRNA types from 34 families are common to all tested tissues, while some of the miRNAs are tissue-enriched or sex-enriched. Quantitative real-time PCR analysis further demonstrated that let-7a, miR-122, and miR-9-3p were downregulated in hypoxic female medaka, while miR-2184 was specifically upregulated in the testis of hypoxic male fish. This is the first study to identify miRNAs in O. melastigma using small RNA deep sequencing technology. Because miRNA expression is highly conserved between marine medaka and other vertebrates, marine medaka may serve as a good model for studies on the functional roles of miRNAs in hypoxia stress response and signaling in marine fish.

Targeted next-generation sequencing on Hirschsprung disease: a pilot study exploits DNA pooling.

To adopt an efficient approach of identifying rare variants possibly related to Hirschsprung disease (HSCR), a pilot study was set up to evaluate the performance of a newly designed protocol for next generation targeted resquencing. In total, 20 Chinese HSCR patients and 20 Chinese sex-matched individuals with no HSCR were included, for which coding sequences (CDS) of 62 genes known to be in signaling pathways relevant to enteric nervous system development were selected for capture and sequencing. Blood DNAs from eight pools of five cases or controls were enriched by PCR-based RainDance technology (RDT) and then sequenced on a 454 FLX platform. As technical validation, five patients from case Pool-3 were also independently enriched by RDT, indexed with barcode and sequenced with sufficient coverage. Assessment for CDS single nucleotide variants showed DNA pooling performed well (specificity/sensitivity at 98.4%/83.7%) at the common variant level; but relatively worse (specificity/sensitivity at 65.5%/61.3%) at the rare variant level. Further Sanger sequencing only validated five out of 12 rare damaging variants likely involved in HSCR. Hence more improvement at variant detection and sequencing technology is needed to realize the potential of DNA pooling for large-scale resequencing projects.

Identification of microRNA-like RNAs in mycelial and yeast phases of the thermal dimorphic fungus Penicillium marneffei.

BACKGROUND: Penicillium marneffei is the most important thermal dimorphic fungus causing systemic mycosis in China and Southeast Asia. While miRNAs are increasingly recognized for their roles in post-transcriptional regulation of gene expression in animals and plants, miRNAs in fungi were less well studied and their potential roles in fungal dimorphism were largely unknown. Based on P. marneffei genome sequence, we hypothesize that miRNA-like RNAs (milRNAs) may be expressed in the dimorphic fungus. METHODOLOGY/PRINCIPAL FINDINGS: We attempted to identify milRNAs in P. marneffei in both mycelial and yeast phase using high-throughput sequencing technology. Small RNAs were more abundantly expressed in mycelial than yeast phase. Sequence analysis revealed 24 potential milRNA candidates, including 17 candidates in mycelial and seven in yeast phase. Two genes, dcl-1 and dcl-2, encoding putative Dicer-like proteins and the gene, qde-2, encoding Argonaute-like protein, were identified in P. marneffei. Phylogenetic analysis showed that dcl-2 of P. marneffei was more closely related to the homologues in other thermal dimorphic pathogenic fungi than to Penicillium chrysogenum and Aspergillus spp., suggesting the co-evolution of dcl-2 among the thermal dimorphic fungi. Moreover, dcl-2 demonstrated higher mRNA expression levels in mycelial than yeast phase by 7 folds (P<0.001). Northern blot analysis confirmed the expression of two milRNAs, PM-milR-M1 and PM-milR-M2, only in mycelial phase. Using dcl-1(KO), dcl-2(KO), dcl(DKO) and qde-2(KO) deletion mutants, we showed that the biogenesis of both milRNAs were dependent on dcl-2 but not dcl-1 or qde-2. The mRNA expression levels of three predicted targets of PM-milR-M1 were upregulated in knockdown strain PM-milR-M1 (KD), supporting regulatory function of milRNAs. CONCLUSIONS/SIGNIFICANCE: Our findings provided the first evidence for differential expression of milRNAs in different growth phases of thermal dimorphic fungi and shed light on the evolution of fungal proteins involved in milRNA biogenesis and possible role of post-transcriptional control in governing thermal dimorphism.

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.

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.

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.

Dissemination of pHK01-like incompatibility group IncFII plasmids encoding CTX-M-14 in Escherichia coli from human and animal sources.

Few studies have compared CTX-M encoding plasmids identified in different ecological sources. This study aimed to analyze and compare the molecular epidemiology of plasmids encoding CTX-M-14 among strains from humans and animals. The CTX-M-14 encoding plasmids in 160 Escherichia coli isolates from animal faecal (14 pigs, 16 chickens, 12 cats, 8 cattle, 5 dogs and 3 rodents), human faecal (45 adults and 20 children) and human urine (37 adults) sources in 2002-2010 were characterized by molecular methods. The replicon types of the CTX-M-14 encoding plasmids were IncFII (n=61), I1-Iγ (n=24), other F types (n=23), B/O (n=10), K (n=6), N (n=3), A/C (n=1), HI1 (n=1), HI2 (n=1) and nontypeable (n=30). The genetic environment, ISEcp1 -bla(CTX-M-14) - IS903 was found in 89.7% (52/58), 87.7% (57/65) and 86.5% (32/37) of the animal faecal, human faecal and human urine isolates, respectively. Subtyping of the 61 IncFII incompatibility group plasmids by replicon sequence typing, plasmid PCR-restriction fragment length polymorphism and marker genes (yac, malB, eitA/eitC and parB/A) profiles showed that 31% (18/58), 30.6% (20/65) and 37.8% (14/37) of the plasmids originating from animal faecal, human faecal and human urine isolates, respectively, were pHK01-like. These 52 pHK01-like plasmids originated from diverse human (20 faecal isolates from 2002, 2007 to 2008, 14 urinary isolates from 2004) and animal (all faecal, 1 cattle, 1 chicken, 5 pigs, 9 cats, 1 dog, 1 rodent from 2008 to 2010) sources. In conclusion, this study highlights the importance of the IncFII group, pHK01-like plasmids in the dissemination of CTX-M-14 among isolates from diverse sources.

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.

Complete sequencing of the FII plasmid pHK01, encoding CTX-M-14, and molecular analysis of its variants among Escherichia coli from Hong Kong.

OBJECTIVES: We characterized plasmids encoding CTX-M-14 ß-lactamase originating from Escherichia coli isolates recovered from patients with uncomplicated cystitis or individuals with faecal colonization in Hong Kong from 2002 to 2004. METHODS: Plasmids carrying CTX-M-14 were studied by conjugation, replicon typing, S1 nuclease-PFGE and plasmid PCR-restriction fragment length polymorphism (RFLP). The complete sequence of pHK01, a 70 kb plasmid encoding CTX-M-14 from an E. coli strain, was determined and the results compared with reference plasmids and aligned with GenBank data. RESULTS: The blaCTX-M-14 plasmids could be transferred in 23 of 44 E. coli strains tested. Among the 23 transconjugants, the replicon types of the CTX-M-14-encoding plasmid were FII (n=13), I1-Iγ (n=4), F1B (n=2), FII and I1-Iγ (n=1), K (80 kb, n=1) and undetermined (n=2). Plasmid pHK01 (FII replicon) shares a high degree of homology with R100 except mainly for a 11 kb variable region containing blaCTX-M-14 (with an upstream ISEcp1 and a downstream truncated IS903), an iron transport system, an outer membrane protein (malB, maltoporin) and a putative toxin-antitoxin plasmid stability system (yacABC). It was highly related to blaCTX-M-14 (pKF3-70) and blaCTX-M-24 (pEG356) plasmids reported from mainland China in 2006 and Vietnam in 2007, respectively. Subtyping by a plasmid PCR-RFLP scheme showed that 10 of the 13 FII plasmids originating from isolates collected by multiple laboratories exhibited either identical or highly similar profiles. CONCLUSIONS: This study showed that narrow host-range FII plasmids play important roles in the dissemination of CTX-M-14. FII plasmids closely related to pHK01 have disseminated widely in the Hong Kong community.

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.