Myofibroblasts are a main cell-type of collagen-producing cells during tissue fibrosis, but their origins remains controversial. While bone marrow-derived myofibroblasts in renal fibrosis has been reported, the cell origin and mechanisms regulating their transition into myofibroblasts remain undefined. In the present study, cell lineage tracing studies by adoptive transfer of GFP+ or dye-labelled macrophages identified that monocyte/macrophages from bone marrow can give rise to myofibroblasts via the process of macrophage-myofibroblast transition (MMT) in a mouse model of unilateral ureteric obstruction. The MMT cells were a major source of collagen-producing fibroblasts in the fibrosing kidney, accounting for more than 60% of α-SMA+ myofibroblasts. The MMT process occurred predominantly within M2-type macrophages and was regulated by TGF-β/Smad3 signalling as deletion of Smad3 in the bone marrow compartment of GFP+ chimeric mice prevented the M2 macrophage transition into the MMT cells and progressive renal fibrosis. In vitro studies in Smad3 null bone marrow macrophages also showed that Smad3 was required for TGF-β1-induced MMT and collagen production. In conclusion, we have demonstrated that bone marrow-derived fibroblasts originate from the monocyte/macrophage population via a process of MMT. This process contributes to progressive renal tissue fibrosis and is regulated by TGF-β/Smad3 signalling.

Acute pneumonia developed in a previously healthy man during the outbreak of severe acute respiratory syndrome (SARS) in southern China in March 2003. Antibiotic treatment was ineffective, and he died 8 days after illness onset. Human metapneumovirus was isolated from lung tissue. No other pathogen was found. Other etiologic agents should thus be sought in apparent SARS cases when coronavirus infection cannot be confirmed.

A key abstraction in representing proteomics knowledge is the notion of unique identifiers for individual entities (e.g. proteins) and the massive graph of relationships among them. These relationships are sometimes simple (e.g. synonyms) but are often more complex (e.g. one-to-many relationships in protein family membership).

We develop a statistical framework to study the relationship between chromatin features and gene expression. This can be used to predict gene expression of protein coding genes, as well as microRNAs. We demonstrate the prediction in a variety of contexts, focusing particularly on the modENCODE worm datasets. Moreover, our framework reveals the positional contribution around genes (upstream or downstream) of distinct chromatin features to the overall prediction of expression levels.

The role of oral microbiota in head and neck squamous cell carcinoma (HNSCC) is poorly understood. Here we sought to evaluate the association of the bacterial microbiome with host gene methylation and patient outcomes, and to explore its potential as a biomarker for early detection or intervention. Here we performed 16S rRNA gene amplicon sequencing in sixty-eight HNSCC patients across both tissue and oral rinse samples to identify oral bacteria with differential abundance between HNSCC and controls. A subset of thirty-one pairs of HNSCC tumor tissues and the adjacent normal tissues were characterized for host gene methylation profile using bisulfite capture sequencing. We observed significant enrichments of Fusobacterium and Peptostreptococcus in HNSCC tumor tissues when compared to the adjacent normal tissues, and in HNSCC oral rinses when compared to healthy subjects, while ten other bacterial genera were largely depleted. These HNSCC-related bacteria were discriminative for HNSCC and controls with area under the receiver operating curves (AUCs) of 0.84 and 0.86 in tissue and oral rinse samples, respectively. Moreover, Fusobacterium nucleatum abundance in HNSCC cases was strongly associated with non-smokers, lower tumor stage, lower rate of recurrence, and improved disease-specific survival. An integrative analysis identified that enrichment of F. nucleatum was associated with host gene promoter methylation, including hypermethylation of tumor suppressor genes LXN and SMARCA2, for which gene expressions were downregulated in the HNSCC cohort from The Cancer Genome Atlas. In conclusion, we identified a taxonomically defined microbial consortium associated with HNSCC that may have clinical potential regarding biomarkers for early detection or intervention. Host-microbe interactions between F. nucleatum enrichment and clinical outcomes or host gene methylation imply a potential role of F. nucleatum as a pro-inflammatory driver in initiating HNSCC without traditional risk factors, which warrants further investigation for the underlying mechanisms.

Haploinsufficiency for SOX9, the master chondrogenesis transcription factor, can underlie campomelic dysplasia (CD), an autosomal dominant skeletal malformation syndrome, because heterozygous Sox9 null mice recapitulate the bent limb (campomelia) and some other phenotypes associated with CD. However, in vitro cell assays suggest haploinsufficiency may not apply for certain mutations, notably those that truncate the protein, but in these cases in vivo evidence is lacking and underlying mechanisms are unknown. Here, using conditional mouse mutants, we compared the impact of a heterozygous Sox9 null mutation (Sox9+/-) with the Sox9+/Y440X CD mutation that truncates the C-terminal transactivation domain but spares the DNA-binding domain. While some Sox9+/Y440X mice survived, all Sox9+/- mice died perinatally. However, the skeletal defects were more severe and IHH signaling in developing limb cartilage was significantly enhanced in Sox9+/Y440X compared with Sox9+/-. Activating Sox9Y440X specifically in the chondrocyte-osteoblast lineage caused milder campomelia, and revealed cell- and noncell autonomous mechanisms acting on chondrocyte differentiation and osteogenesis in the perichondrium. Transcriptome analyses of developing Sox9+/Y440X limbs revealed dysregulated expression of genes for the extracellular matrix, as well as changes consistent with aberrant WNT and HH signaling. SOX9Y440X failed to interact with β-catenin and was unable to suppress transactivation of Ihh in cell-based assays. We propose enhanced HH signaling in the adjacent perichondrium induces asymmetrically localized excessive perichondrial osteogenesis resulting in campomelia. Our study implicates combined haploinsufficiency/hypomorphic and dominant-negative actions of SOX9Y440X, cell-autonomous and noncell autonomous mechanisms, and dysregulated WNT and HH signaling, as the cause of human campomelia.

Whole genome duplication (WGD) has occurred in relatively few sexually reproducing invertebrates. Consequently, the WGD that occurred in the common ancestor of horseshoe crabs ~135 million years ago provides a rare opportunity to decipher the evolutionary consequences of a duplicated invertebrate genome. Here, we present a high-quality genome assembly for the mangrove horseshoe crab Carcinoscorpius rotundicauda (1.7 Gb, N50 = 90.2 Mb, with 89.8% sequences anchored to 16 pseudomolecules, 2n = 32), and a resequenced genome of the tri-spine horseshoe crab Tachypleus tridentatus (1.7 Gb, N50 = 109.7 Mb). Analyses of gene families, microRNAs, and synteny show that horseshoe crabs have undergone three rounds (3R) of WGD. Comparison of C. rotundicauda and T. tridentatus genomes from populations from several geographic locations further elucidates the diverse fates of both coding and noncoding genes. Together, the present study represents a cornerstone for improving our understanding of invertebrate WGD events on the evolutionary fates of genes and microRNAs, at both the individual and population level. We also provide improved genomic resources for horseshoe crabs, of applied value for breeding programs and conservation of this fascinating and unusual invertebrate lineage.

HOXA transcript at the distal tip (HOTTIP), a long noncoding RNA, is upregulated in pancreatic ductal adenocarcinoma (PDAC), but the HOTTIP-mediated oncogenic pathway is not fully understood. We identified canonical HOTTIP-HOXA13 targets, CYP26B1, CLIC5, CHI3L1 and UCP2-responsible for cell growth and cell invasion. Genome-wide analysis revealed that 38% of HOTTIP-regulated genes contain H3K4me3 and HOTTIP enrichment at their promoters, without HOXA13 binding. HOTTIP complexes with WDR5-MLL1 to trans-activate oncogenic proteins CYB5R2, SULT1A1, KIF26A, SLC1A4, and TSC22D1 by directly inducing H3K4me3 at their promoters. The WDR5, MLL1, and H3K4me3 levels at their promoters and their expression levels are sensitive to HOTTIP expression. These results indicate the importance of the noncanonical trans-acting HOTTIP-WDR5-MLL1 pathway in the HOTTIP regulatory mechanism by promoting oncogenic protein expression. Furthermore, HOTTIP is regulated by miR-497 in PDAC cells, but HOTTIP is negatively correlated with miR-497 levels in PDAC tissues. In conclusion, HOTTIP is upregulated in PDAC due to the loss of the inhibitory miR-497; HOTTIP promotes PDAC progression through the canonical HOTTIP-HOXA13 axis. A novel noncanonical trans-acting HOTTIP-WDR5-MLL1-H3K4me3 pathway is also delineated.

Nasopharyngeal carcinoma (NPC) is an Epstein-Barr virus (EBV)-associated epithelial malignancy. The high expression of BART-miRNAs (miR-BARTs) during latent EBV infection in NPC strongly supports their pathological importance in cancer progression. Recently, we found that several BART-miRNAs work co-operatively to modulate the DNA damage response (DDR) by reducing Ataxia-telangiectasia-mutated (ATM) activity. In this study, we further investigated the role of miR-BARTs on DDR. The immunohistochemical study showed that the DNA repair gene, BRCA1, is consistently down-regulated in primary NPCs. Using computer prediction programs and a series of reporter assays, we subsequently identified the negative regulatory role of BART2-3p, BART12, BART17-5p and BART19-3p in BRCA1 expression. The ectopic expression of these four miR-BARTs suppressed endogenous BRCA1 expression in EBV-negative epithelial cell lines, whereas BRCA1 expression was enhanced by repressing endogenous miR-BARTs activities in C666-1 cells. More importantly, suppressing BRCA1 expression in nasopharyngeal epithelial cell lines using miR-BART17-5p and miR-BART19-3p mimics reduced the DNA repair capability and increased the cell sensitivity to the DNA-damaging chemotherapeutic drugs, cisplatin and doxorubicin. Our findings suggest that miR-BARTs play a novel role in DDR and may facilitate the development of effective NPC therapies.

Epstein-Barr virus (EBV) has been recently found to generate novel circular RNAs (circRNAs) through backsplicing. However, comprehensive catalogs of EBV circRNAs in other cell lines and their functional characterization are still lacking. In this study, we have identified a list of putative EBV circRNAs in GM12878, an EBV-transformed lymphoblastoid cell line, with a significant majority encoded from the EBV latent genes. A novel EBV circRNA derived from the exon 5 of LMP-2 gene which exhibited highest prevalence, was further validated using RNase R assay and Sanger sequencing. This circRNA, which we term circLMP-2_e5, can be universally detected in a panel of EBV-positive cell lines modelling different latency programs. It ranges from lower expression in nasopharyngeal carcinoma (NPC) cells to higher expression in B cells, and is localized to both the cytoplasm and the nucleus. We provide evidence that circLMP-2_e5 is expressed concomitantly with its cognate linear LMP-2 RNA upon EBV lytic reactivation, and may be produced as a result of exon skipping, with its circularization possibly occurring without the involvement of cis elements in the short flanking introns. Furthermore, we show that circLMP-2_e5 is not involved in regulating cell proliferation, host innate immune response, its linear parental transcripts, or EBV lytic reactivation. Taken together, our study expands the current repertoire of putative EBV circRNAs, broadens our understanding of the biology of EBV circRNAs, and lays the foundation for further investigation of their function in the EBV life cycle and disease development.

Diabetic nephropathy (DN) is a major cause of end-stage renal disease, but treatment remains ineffective. C-reactive protein (CRP) is pathogenic in DN, which significantly correlated with dipeptidyl peptidase-4 (DPP4) expression in diabetic patients with unknown reason. Here, using our unique CRPtg-db/db mice, we observed human CRP markedly induced renal DPP4 associated with enhanced kidney injury compared with db/db mice. Interestingly, linagliptin, a US Food and Drug Administration (FDA)-approved specific DPP4 inhibitor, effectively blocked this CRP-driven DN in the CRPtg-db/db mice. Mechanistically, CRP evoked DPP4 in cultured renal tubular epithelial cells, where CD32b/nuclear factor κB (NF-κB) signaling markedly enriched p65 binding on the DPP4 promoter region to increase its transcription. Unexpectedly, we further discovered that CRP triggers dimerization of DPP4 with CD32b at protein level, forming a novel DPP4/CD32b/NF-κB signaling circuit for promoting CRP-mediated DN. More importantly, linagliptin effectively blocked the circuit, thereby inhibiting the CRP/CD32b/NF-κB-driven renal inflammation and fibrosis. Thus, DPP4 may represent a precise druggable target for CRP-driven DN.

Circular RNAs (circRNAs) are abundantly expressed in cancer. Their resistance to exonucleases enables them to have potentially stable interactions with different types of biomolecules. Alternative splicing can create different circRNA isoforms that have different sequences and unequal interaction potentials. The study of circRNA function thus requires knowledge of complete circRNA sequences. Here we describe psirc, a method that can identify full-length circRNA isoforms and quantify their expression levels from RNA sequencing data. We confirm the effectiveness and computational efficiency of psirc using both simulated and actual experimental data. Applying psirc on transcriptome profiles from nasopharyngeal carcinoma and normal nasopharynx samples, we discover and validate circRNA isoforms differentially expressed between the two groups. Compared with the assumed circular isoforms derived from linear transcript annotations, some of the alternatively spliced circular isoforms have 100 times higher expression and contain substantially fewer microRNA response elements, showing the importance of quantifying full-length circRNA isoforms.

Previously, we reported that GPR30 activation by the receptor-specific, non-estrogenic ligand G-1 inhibited in vitro and in vivo growth of prostate cancer (PCa) cells via sustained Erk1/2 activation. Mechanism underlying the sustained Erk1/2 activation for PCa cell growth inhibition remains unclear. Here we report that G-1, through GPR30, activated Gαi1 proteins to sustain Erk1/2 activation but failed to activate adenylyl cyclase (AC) for cAMP production in PCa cells. The chemical-induced activation of AC-cAMP-PKA signaling attenuated Erk1/2 activity and blocked the cell growth inhibitory effects of G-1. Furthermore, PCa predominantly expressed Gαi1 proteins. Silencing of Gαi1 expression blocked the inhibitory effects of G-1 on PCa cell growth. By gene expression profiling, GPR30 activation by G-1 interfered expression of cell cycle regulators and machinery elements to modulate PCa cell growth and the RACGAP1 interactome to control metastatic properties. In this regard, we demonstrated that G-1 inhibited PCa cell migration and invasion with reduced formations of filopodia and stress fibers through a GPR30-dependent pathway. Taken together, our findings revealed the underlying mechanism for sustaining Erk1/2 activation upon GPR30 activation by G-1 in PCa cells and the GPR30-mediated pathways in controlling PCa cell growth and metastatic properties.

Protein interactions play significant roles in complex diseases. We analyzed peripheral blood mononuclear cells (PBMC) transcriptome using a multi-method strategy. We constructed a tissue-specific interactome (T2Di) and identified 420 molecular signatures associated with T2D-related comorbidity and symptoms, mainly implicated in inflammation, adipogenesis, protein phosphorylation and hormonal secretion. Apart from explaining the residual associations within the DIAbetes Genetics Replication And Meta-analysis (DIAGRAM) study, the T2Di signatures were enriched in pathogenic cell type-specific regulatory elements related to fetal development, immunity and expression quantitative trait loci (eQTL). The T2Di revealed a novel locus near a well-established GWAS loci AChE, in which SRRT interacts with JAZF1, a T2D-GWAS gene implicated in pancreatic function. The T2Di also included known anti-diabetic drug targets (e.g. PPARD, MAOB) and identified possible druggable targets (e.g. NCOR2, PDGFR). These T2Di signatures were validated by an independent computational method, and by expression data of pancreatic islet, muscle and liver with some of the signatures (CEBPB, SREBF1, MLST8, SRF, SRRT and SLC12A9) confirmed in PBMC from an independent cohort of 66 T2D and 66 control subjects. By combining prior knowledge and transcriptome analysis, we have constructed an interactome to explain the multi-layered regulatory pathways in T2D.

To find signature features shared by various ncRNA sub-types and characterize novel ncRNAs, we have developed a method, RNAfeature, to investigate >600 sets of genomic and epigenomic data with various evolutionary and biophysical scores. RNAfeature utilizes a fine-tuned intra-species wrapper algorithm that is followed by a novel feature selection strategy across species. It considers long distance effect of certain features (e.g. histone modification at the promoter region). We finally narrow down on 10 informative features (including sequences, structures, expression profiles and epigenetic signals). These features are complementary to each other and as a whole can accurately distinguish canonical ncRNAs from CDSs and UTRs (accuracies: >92% in human, mouse, worm and fly). Moreover, the feature pattern is conserved across multiple species. For instance, the supervised 10-feature model derived from animal species can predict ncRNAs in Arabidopsis (accuracy: 82%). Subsequently, we integrate the 10 features to define a set of noncoding potential scores, which can identify, evaluate and characterize novel noncoding RNAs. The score covers all transcribed regions (including unconserved ncRNAs), without requiring assembly of the full-length transcripts. Importantly, the noncoding potential allows us to identify and characterize potential functional domains with feature patterns similar to canonical ncRNAs (e.g. tRNA, snRNA, miRNA, etc) on ∼70% of human long ncRNAs (lncRNAs).

Fulvestrant (ICI-182,780) has recently been shown to effectively suppress prostate cancer cell growth in vitro and in vivo. But it is unclear whether microRNAs play a role in regulating oncogene expression in fulvestrant-treated prostate cancer. Here, this study reports hsa-miR-765 as the first fulvestrant-driven, ERβ-regulated miRNA exhibiting significant tumor suppressor activities like fulvestrant, against prostate cancer cell growth via blockage of cell-cycle progression at the G2/M transition, and cell migration and invasion possibly via reduction of filopodia/intense stress-fiber formation. Fulvestrant was shown to upregulate hsa-miR-765 expression through recruitment of ERβ to the 5'-regulatory-region of hsa-miR-765. HMGA1, an oncogenic protein in prostate cancer, was identified as a downstream target of hsa-miR-765 and fulvestrant in cell-based experiments and a clinical study. Both the antiestrogen and the hsa-miR-765 mimic suppressed HMGA1 protein expression. In a neo-adjuvant study, levels of hsa-miR-765 were increased and HMGA1 expression was almost completely lost in prostate cancer specimens from patients treated with a single dose (250 mg) of fulvestrant 28 days before prostatectomy. These findings reveal a novel fulvestrant signaling cascade involving ERβ-mediated transcriptional upregulation of hsa-miR-765 that suppresses HMGA1 protein expression as part of the mechanism underlying the tumor suppressor action of fulvestrant in prostate cancer.

Comprehensive whole-genome structural variation detection is challenging with current approaches. With diploid cells as DNA source and the presence of numerous repetitive elements, short-read DNA sequencing cannot be used to detect structural variation efficiently. In this report, we show that genome mapping with long, fluorescently labeled DNA molecules imaged on nanochannel arrays can be used for whole-genome structural variation detection without sequencing. While whole-genome haplotyping is not achieved, local phasing (across >150-kb regions) is routine, as molecules from the parental chromosomes are examined separately. In one experiment, we generated genome maps from a trio from the 1000 Genomes Project, compared the maps against that derived from the reference human genome, and identified structural variations that are >5 kb in size. We find that these individuals have many more structural variants than those published, including some with the potential of disrupting gene function or regulation.

Severe acute respiratory syndrome (SARS) is a recently emerged infection from a novel coronavirus (CoV). Apart from fever and respiratory complications, gastrointestinal symptoms are frequently observed in patients with SARS but the significance remains undetermined. Herein, we describe the clinical, pathologic, and virologic features of the intestinal involvement of this new viral infection.

Very often genome-wide data analysis requires the interoperation of multiple databases and analytic tools. A large number of genome databases and bioinformatics applications are available through the web, but it is difficult to automate interoperation because: 1) the platforms on which the applications run are heterogeneous, 2) their web interface is not machine-friendly, 3) they use a non-standard format for data input and output, 4) they do not exploit standards to define application interface and message exchange, and 5) existing protocols for remote messaging are often not firewall-friendly. To overcome these issues, web services have emerged as a standard XML-based model for message exchange between heterogeneous applications. Web services engines have been developed to manage the configuration and execution of a web services workflow.

Epstein-Barr virus (EBV) is a herpesvirus associated with nasopharyngeal carcinoma (NPC), gastric carcinoma (GC), and other malignancies. EBV is the first human virus found to express microRNAs (miRNAs), the functions of which remain largely unknown. We report on the regulation of a cellular protein named p53 up-regulated modulator of apoptosis (PUMA) by an EBV miRNA known as miR-BART5, which is abundantly expressed in NPC and EBV-GC cells. Modulation of PUMA expression by miR-BART5 and anti-miR-BART5 oligonucleotide was demonstrated in EBV-positive cells. In addition, PUMA was found to be significantly underexpressed in approximately 60% of human NPC tissues. Although expression of miR-BART5 rendered NPC and EBV-GC cells less sensitive to proapoptotic agents, apoptosis can be triggered by depleting miR-BART5 or inducing the expression of PUMA. Collectively, our findings suggest that EBV encodes an miRNA to facilitate the establishment of latent infection by promoting host cell survival.