Genomic adaptations for arboreal locomotion in Asian flying treefrogs.

SignificanceTo adapt to arboreal lifestyles, treefrogs have evolved a suite of complex traits that support vertical movement and gliding, thus presenting a unique case for studying the genetic basis for traits causally linked to vertical niche expansion. Here, based on two de novo-assembled Asian treefrog genomes, we determined that genes involved in limb development and keratin cytoskeleton likely played a role in the evolution of their climbing systems. Behavioral and morphological evaluation and time-ordered gene coexpression network analysis revealed the developmental patterns and regulatory pathways of the webbed feet used for gliding in Rhacophorus kio.

Comparative genomic investigation of high-elevation adaptation in ectothermic snakes.

Several previous genomic studies have focused on adaptation to high elevations, but these investigations have been largely limited to endotherms. Snakes of the genus Thermophis are endemic to the Tibetan plateau and therefore present an opportunity to study high-elevation adaptations in ectotherms. Here, we report the de novo assembly of the genome of a Tibetan hot-spring snake (Thermophis baileyi) and then compare its genome to the genomes of the other two species of Thermophis, as well as to the genomes of two related species of snakes that occur at lower elevations. We identify 308 putative genes that appear to be under positive selection in Thermophis We also identified genes with shared amino acid replacements in the high-elevation hot-spring snakes compared with snakes and lizards that live at low elevations, including the genes for proteins involved in DNA damage repair (FEN1) and response to hypoxia (EPAS1). Functional assays of the FEN1 alleles reveal that the Thermophis allele is more stable under UV radiation than is the ancestral allele found in low-elevation lizards and snakes. Functional assays of EPAS1 alleles suggest that the Thermophis protein has lower transactivation activity than the low-elevation forms. Our analysis identifies some convergent genetic mechanisms in high-elevation adaptation between endotherms (based on studies of mammals) and ectotherms (based on our studies of Thermophis).

Comprehensive Analysis of lncRNA and mRNA Expression Profiles in Lung Cancer.

BACKGROUND: Lung cancer is one of the most common malignancies across the world. Long noncoding RNAs (lncRNAs) play an important role in the pathology. This study was to compare the lncRNA and mRNA of lung cancer. The aim of this study was to compare the lncRNA and mRNA expression profiles, their related biological functions, and pathways in lung cancer. METHODS: Human lung cancer mRNA expression datasets were searched and downloaded from NCBI-GEO. LncRNA expression profiles in lung cancer were detected by the Agilent Human LncRNA Microarray V3.0. Differential expression analysis was conducted between cases and controls in mRNA and lncRNA. The starBase web server v2.0 was used to decipher lncRNA-protein interactions. DAVID Bioinformatics Resources 6.7 was used to perform GO Biological Processes and KEGG pathway enrichment analysis of these dysregulated mRNA and lncRNA target genes. RESULTS: The study showed that differentially expressed lncRNA target genes included almost all of the differential expression genes from mRNA. Furthermore, these dysregulated lncRNAs reflected more comprehensive impairment in functional enrichment than dysregulated mRNAs. In addition, five top downregulated lncRNAs had more remarkable fold changes than top downregulated mRNAs, especially FENDRR. CONCLUSIONS: Amount of lncRNAs and mRNAs were differentially expressed in lung cancer. The degrees of difference in lncRNAs were more than mRNAs. It may provide valuable help for an effective strategy in diagnosis and prevention.

Co-delivery of photosensitizer and diclofenac through sequentially responsive bilirubin nanocarriers for combating hypoxic tumors.

Considering that photodynamic therapy (PDT)-induced oxygen consumption and microvascular damage could exacerbate hypoxia to drive more glycolysis and angiogenesis, a novel approach to potentiate PDT and overcome the resistances of hypoxia is avidly needed. Herein, morpholine-modified PEGylated bilirubin was proposed to co-deliver chlorin e6, a photosensitizer, and diclofenac (Dc). In acidic milieu, the presence of morpholine could enable the nanocarriers to selectively accumulate in tumor cells, while PDT-generated reactive oxidative species (ROS) resulted in the collapse of bilirubin nanoparticles and rapid release of Dc. Combining with Dc showed a higher rate of apoptosis over PDT alone and simultaneously triggered a domino effect, including blocking the activity and expression of lactate dehydrogenase A (LDHA), interfering with lactate secretion, suppressing the activation of various angiogenic factors and thus obviating hypoxia-induced resistance-glycolysis and angiogenesis. In addition, inhibition of hypoxia-inducible factor-1α (HIF-1α) by Dc alleviated hypoxia-induced resistance. This study offered a sequentially responsive platform to achieve sufficient tumor enrichment, on-demand drug release and superior anti-tumor outcomes in vitro and in vivo.

Cholesterol crystals activate NLRP3 inflammasomes and promote gallstone formation by increasing mucin secretion.

Cholesterol crystals participate in cholesterol nucleation; however, the role of cholesterol crystals in gallstone development is unknown. Mucin secretion contributes to increased size of gallstones. Cholesterol crystals activate inflammasomes and participate in many sterile inflammation related human diseases. We investigated the role of cholesterol crystals and mucins in sterile inflammation and gallstone enlargement. We found that expression of mucin 5AC (MUC5AC), NOD-, LRR- and pyrin domain-containing protein 3 (NLRP3) and interleukin-1b (IL-1b) was increased significantly in tissues adjacent to gallstones. Experiments in vitro showed that cholesterol crystals promote MUC5AC secretion; they also increase expression of NLRP3, NLR family CARD domain-containing 4 (NLRC4), apoptosis-associated speck-like protein containing caspase recruitment domain (ASC), and cleaved caspase-1 in biliary epithelial cells. Inhibition of Inflammasomes by NLRP3, ASC or caspase-1 small interfering RNAs reduced MUC5AC secretion. Also, the IL-1 receptor antagonist, IL1RA, and caspase-1 inhibitor, Ac-YVAD, both inhibited MUC5AC secretion induced by cholesterol crystals. We found that inflammasome activation participates in cholesterol crystal induced mucin secretion and gallstone development.

Mangiferin, an anti-HIV-1 agent targeting protease and effective against resistant strains.

The anti-HIV-1 activity of mangiferin was evaluated. Mangiferin can inhibit HIV-1(Ⅲ)(B) induced syncytium formation at non-cytotoxic concentrations, with a 50% effective concentration (EC50) at 16.90 µM and a therapeutic index (TI) above 140. Mangiferin also showed good activities in other laboratory-derived strains, clinically isolated strains and resistant HIV-1 strains. Mechanism studies revealed that mangiferin might inhibit the HIV-1 protease, but is still effective against HIV peptidic protease inhibitor resistant strains. A combination of docking and pharmacophore methods clarified possible binding modes of mangiferin in the HIV-1 protease. The pharmacophore model of mangiferin consists of two hydrogen bond donors and two hydrogen bond acceptors. Compared to pharmacophore features found in commercially available drugs, three pharmacophoric elements matched well and one novel pharmacophore element was observed. Moreover, molecular docking analysis demonstrated that the pharmacophoric elements play important roles in binding HIV-1 protease. Mangiferin is a novel nonpeptidic protease inhibitor with an original structure that represents an effective drug development strategy for combating drug resistance.

High PKM2 expression is independently correlated with decreased overall survival in hepatocellular carcinoma.

Hepatocellular carcinoma (HCC) is one of the most lethal and malignant types of cancer that affects global human health. The present study aimed to investigate the effect of pyruvate kinase muscle isozyme M2 (PKM2) expression on the clinical features and prognosis of HCC. The present study employed univariate logistic regression to investigate the correlation between PKM2 expression and clinical features. Univariate and multivariate Cox regression analyses were performed to estimate the independent effect of PKM2 expression on survival status. The results revealed that patients in the high PKM2 group (≥11.25) exhibited significantly lower creatinine levels (P=0.043), higher fetoprotein levels (P<0.001), advanced stage (P<0.001) and higher grade (P=0.004) compared with patients with low PKM2 expression levels (<11.25). In addition, patients with high PKM2 expression exhibited poor prognosis compared with patients with low PKM2 expression. After correcting the covariates, PKM2 expression remains significantly associated with reduced overall survival (P<0.05). These findings suggested that PKM2 is an independent risk factor for HCC and provides valuable information for future studies on the pathogenesis of HCC and drug discovery.

Pharmacophore-Map-Pick: A Method to Generate Pharmacophore Models for All Human GPCRs.

GPCR-based drug discovery is hindered by a lack of effective screening methods for most GPCRs that have neither ligands nor high-quality structures. With the aim to identify lead molecules for these GPCRs, we developed a new method called Pharmacophore-Map-Pick to generate pharmacophore models for all human GPCRs. The model of ADRB2 generated using this method not only predicts the binding mode of ADRB2-ligands correctly but also performs well in virtual screening. Findings also demonstrate that this method is powerful for generating high-quality pharmacophore models. The average enrichment for the pharmacophore models of the 15 targets in different GPCR families reached 15-fold at 0.5 % false-positive rate. Therefore, the pharmacophore models can be applied in virtual screening directly with no requirement for any ligand information or shape constraints. A total of 2386 pharmacophore models for 819 different GPCRs (99 % coverage (819/825)) were generated and are available at http://bsb.kiz.ac.cn/GPCRPMD.

In silico identification of anti-cancer compounds and plants from traditional Chinese medicine database.

There is a constant demand to develop new, effective, and affordable anti-cancer drugs. The traditional Chinese medicine (TCM) is a valuable and alternative resource for identifying novel anti-cancer agents. In this study, we aim to identify the anti-cancer compounds and plants from the TCM database by using cheminformatics. We first predicted 5278 anti-cancer compounds from TCM database. The top 346 compounds were highly potent active in the 60 cell lines test. Similarity analysis revealed that 75% of the 5278 compounds are highly similar to the approved anti-cancer drugs. Based on the predicted anti-cancer compounds, we identified 57 anti-cancer plants by activity enrichment. The identified plants are widely distributed in 46 genera and 28 families, which broadens the scope of the anti-cancer drug screening. Finally, we constructed a network of predicted anti-cancer plants and approved drugs based on the above results. The network highlighted the supportive role of the predicted plant in the development of anti-cancer drug and suggested different molecular anti-cancer mechanisms of the plants. Our study suggests that the predicted compounds and plants from TCM database offer an attractive starting point and a broader scope to mine for potential anti-cancer agents.

Metabolic modeling of common Escherichia coli strains in human gut microbiome.

The recent high-throughput sequencing has enabled the composition of Escherichia coli strains in the human microbial community to be profiled en masse. However, there are two challenges to address: (1) exploring the genetic differences between E. coli strains in human gut and (2) dynamic responses of E. coli to diverse stress conditions. As a result, we investigated the E. coli strains in human gut microbiome using deep sequencing data and reconstructed genome-wide metabolic networks for the three most common E. coli strains, including E. coli HS, UTI89, and CFT073. The metabolic models show obvious strain-specific characteristics, both in network contents and in behaviors. We predicted optimal biomass production for three models on four different carbon sources (acetate, ethanol, glucose, and succinate) and found that these stress-associated genes were involved in host-microbial interactions and increased in human obesity. Besides, it shows that the growth rates are similar among the models, but the flux distributions are different, even in E. coli core reactions. The correlations between human diabetes-associated metabolic reactions in the E. coli models were also predicted. The study provides a systems perspective on E. coli strains in human gut microbiome and will be helpful in integrating diverse data sources in the following study.

Interferon-lambda3 (IFN-λ3) and its cognate receptor subunits in tree shrews (Tupaia belangeri): genomic sequence retrieval, molecular identification and expression analysis.

Type III IFNs (IFN-λs) constitute a new subfamily with antiviral activities by signaling through a unique receptor complex composed of IFN-λs receptor 1 (IFNλR1) and interleukin-10 receptor 2 (IL10R2). As tree shrews (Tupaia belangeri) have shown susceptiblility to several human viruses, they are a potentially important model for analyzing viral infection. However, little is known about their IFN-λs system. We used the tree shrew genome to retrieve IFN-λs and their receptor contig sequences by BLASTN and BLASTZ algorithms, and GenScan was used to scan transcripts from the putative contig sequences. RT-PCR and bioinformatic methods were then used to clone and characterize the IFN-λs system. Due to its highest identity with human IFN-λ3, we opted to define one intact IFN-λ gene, tsIFN-λ3, as well as its two receptor subunits, tsIFNλR1 and tsIL10R2. Additionally, our results showed that tsIFN-λ3 contained many features conserved in IFN-λ3 genes from other mammals, including conserved signal peptide cleavage and glycosylation sites, and several residues responsible for binding to the type III IFNR. We also found six transcript variants in the receptors: three in tsIFNλR1, wherein different extracellular regions exist in three transmembrane proteins, resulting in different affinities with IFN-λs; and three more variants in tsIL10R2, encoding one transmembrane and two soluble proteins. Based on tissue distribution in the liver, heart, brain, lung, intestine, kidney, spleen, and stomach, we found that IFN-λs receptor complex was expressed in a variety of organs although the expression level differed markedly between them. As the first study to find transcript variants in IL-10R2, our study offers novel insights that may have important implications for the role of IFN-λs in tree shrews' susceptibility with a variety of human viruses, bolstering the arguments for using tree shrews as an animal model in the study of human viral infections.

[Analysis of the molecular characteristics and cloning of full-length coding sequence of interleukin-2 in tree shrews].

While the tree shrew (Tupaia belangeri chinensis) is an excellent animal model for studying the mechanisms of human diseases, but few studies examine interleukin-2 (IL-2), an important immune factor in disease model evaluation. In this study, a 465 bp of the full-length IL-2 cDNA encoding sequence was cloned from the RNA of tree shrew spleen lymphocytes, which were then cultivated and stimulated with ConA (concanavalin). Clustal W 2.0 was used to compare and analyze the sequence and molecular characteristics, and establish the similarity of the overall structure of IL-2 between tree shrews and other mammals. The homology of the IL-2 nucleotide sequence between tree shrews and humans was 93%, and the amino acid homology was 80%. The phylogenetic tree results, derived through the Neighbour-Joining method using MEGA5.0, indicated a close genetic relationship between tree shrews, Homo sapiens, and Macaca mulatta. The three-dimensional structure analysis showed that the surface charges in most regions of tree shrew IL-2 were similar to between tree shrews and humans; however, the N-glycosylation sites and local structures were different, which may affect antibody binding. These results provide a fundamental basis for the future study of IL-2 monoclonal antibody in tree shrews, thereby improving their utility as a model.

[Genome-wide prediction of interferon family members of tree shrew and their molecular characteristics analysis].

Interferons (IFNs) represent proteins with antiviral activities that are secreted from cells in response to a variety of stimuli. In addition to antiviral, antibacterial and anti-parasitic host-defense functions they are now also recognized as crucial regulators of cell proliferation, differentiation, survival and death as well as activators of specialized cell functions particularly in the immune system and play important roles in infectious and inflammatory diseases, autoimmunity and cancer. Tree shrews (Tupaia belangeri) were found to be susceptible to several human viruses and therefore are widely regarded as good models for analyzing mechanism of human diseases. In this report, we have forecasted the interferon family members of tree shrew from its genome mainly using the methods like Blast (whole genome shotgun sequence) and gene prediction. Our data show that tree shrew interferon system includes: type I IFN: α (five subtypes), ß, ω, κ, epsilon, δ; type II IFN: γ; type III IFN: λ1, λ2/3. Furthermore, the predicted structures of α and λ have similar character with those of other mammals. However, there are some differences in cysteine position and N-glycosylation numbers between human and Tree shrew IFNs. These results provide fundamental basis for further molecular cloning and function analysis of tree shrew IFNs in future.

[Cloning of full-length coding sequence of tree shrew CD4 and prediction of its molecular characteristics].

The tree shrews, as an ideal animal model receiving extensive attentions to human disease research, demands essential research tools, in particular cellular markers and monoclonal antibodies for immunological studies. In this paper, a 1 365 bp of the full-length CD4 cDNA encoding sequence was cloned from total RNA in peripheral blood of tree shrews, the sequence completes two unknown fragment gaps of tree shrews predicted CD4 cDNA in the GenBank database, and its molecular characteristics were analyzed compared with other mammals by using biology software such as Clustal W2.0 and so forth. The results showed that the extracellular and intracellular domains of tree shrews CD4 amino acid sequence are conserved. The tree shrews CD4 amino acid sequence showed a close genetic relationship with Homo sapiens and Macaca mulatta. Most regions of the tree shrews CD4 molecule surface showed positive charges as humans. However, compared with CD4 extracellular domain D1 of human, CD4 D1 surface of tree shrews showed more negative charges, and more two N-glycosylation sites, which may affect antibody binding. This study provides a theoretical basis for the preparation and functional studies of CD4 monoclonal antibody.

[An extension strategy of Discovery Studio 2.0 for non-bonded interaction energy automatic calculation at the residue level].

Non-bonded interaction forces play crucial roles in molecular recognition and binding in biological systems. However, it is difficult for traditional methods to automatically calculate and batch the non-bonded energy at the residue level. In recent years, many studies have focused on non-bonded interactions and developed tools to calculate and analyze such interactions. In this study, we present a highly automated approach for the calculation of non-bonded energy. Our strategy invoked protocols relevant to non-bonded interactions within Discovery Studio 2.0 (DS2.0, Accelrys Inc.) bottom module using Perl script, and determined the direct command line operation of calculating non-bonded interaction energy batches without accessing the graphical interface of DS. This approach extended the DS2.0 module and was applied to a recent study of complex structure analysis.

A novel enzyme-linked immunosorbent assay for screening HIV-1 fusion inhibitors targeting HIV-1 Gp41 core structure.

The gp41 subunit of the human immunodeficiency virus type 1 (HIV-1) envelope glycoprotein mediates the fusion of viral and host cell membranes. As the HIV-1 enters the host cells, the 2 helical regions, HR1 and HR2, in the ectodomain of gp41 can form a 6-helix bundle, which brings the viral and target cell membranes to close proximity and serves as an attractive target for developing HIV-1 fusion inhibitors. Now, there are several cell- and molecule-based assays to identify potential HIV-1 fusion inhibitors targeting gp41. However, these assays cannot be used universally because they are time-consuming, inconvenient, and expensive. In the present study, the authors expressed and purified GST-HR121 and C43-30a proteins that were derived from the HIV-1 gp41 ectodomain region. GST-HR121 has a function similar to the HR1 peptide of gp41, whereas C43-30a is an HR2-derived peptide that added 50 amino acid residues (aa) in the N-terminal of C43. Further research found they could interact with each other, and a potential HIV-1 fusion inhibitor could inhibit this interaction. On the basis of this fact, a novel, rapid, and economic enzyme-linked immunosorbent assay was established, which can be developed for high-throughput screening of HIV-1 fusion inhibitors.

[Cloning of full-length coding sequence of tree shrew CD3E and prediction of its molecular characteristics].

The use of tree shrews (Tupaia belangeri) in human disease studies demands essential research tools, in particular cellular markers and their monoclonal antibodies for immunological studies. Here we cloned the full-length cDNAs encoding CD3E from total RNA of the spleen, liver and peripheral blood of tree shrews and analyzed their structural characteristics in comparison with other mammals by Discovery Studio software. The results showed that the open reading frame sequence of tree shrew CD3E was 582 bp, encoding 194 amino acids. The overall structure of tree shrew CD3E protein was similar to its counterparts of other mammals, intracellular and transmembrane domain highly conserved. However, detailed analysis revealed two potential glycosylation sites and different surface charges in the extracellular domain. Availability of the entire open-reading-frame and related sequence information would therefore facilitate the preparation of monoclonal antibodies against tree shrew CD3 and further studies for its function.