AncestryPainter: A Graphic Program for Displaying Ancestry Composition of Populations and Individuals

Ancestry composition of populations and individuals has been extensively investigated in recent years due to advances in the genotyping and sequencing technologies. As the number of populations and individuals used for ancestry inference increases remarkably, say more than 100 populations or 1000 individuals, it is usually challenging to present the ancestry composition in a traditional way using a rectangular graph. To address this issue, we developed a program, AncestryPainter, which can illustrate the ancestry composition of populations and individuals with a rounded and nice-looking graph to save space. Individuals are depicted as length-fixed bars partitioned into colored segments representing different ancestries, and the population of interest can be highlighted as a pie chart in the center of the circle plot. In addition, AncestryPainter can also be applied to display personal ancestry in a way similar to that for displaying population ancestry. AncestryPainter is publicly available at http://www.picb.ac.cn/PGG/resource.php.     

CircPlant:An Integrated Tool for circRNA Detection and Functional Prediction in Plants

The recent discovery of circular RNAs (circRNAs) and characterization of their functional roles have opened a new avenue for understanding the biology of genomes. circRNAs have been implicated to play important roles in a variety of biological processes, but their precise functions remain largely elusive. Currently, a few approaches are available for novel circRNA prediction, but almost all these methods are intended for animal genomes. Considering that the major differences between the organization of plant and mammal genomes cannot be neglected, a plant-specific method is needed to enhance the validity of plant circRNA identification. In this study, we present CircPlant, an integrated tool for the exploration of plant circRNAs, potentially acting as competing endogenous RNAs (ceRNAs), and their potential functions. With the incorporation of several unique plant-specific criteria, CircPlant can accurately detect plant circRNAs from high-throughput RNA-seq data. Based on comparison tests on simulated and real RNA-seq datasets from Arabidopsis thaliana and Oryza sativa, we show that CircPlant outperforms all evaluated competing tools in both accuracy and efficiency. CircPlant is freely available at http://bis.zju.edu.cn/circplant.     

TSNAdb: A Database for Tumor-specific Neoantigens from Immunogenomics Data Analysis

Tumor-specific neoantigens have attracted much attention since they can be used as biomarkers to predict therapeutic effects of immune checkpoint blockade therapy and as potential targets for cancer immunotherapy. In this study, we developed a comprehensive tumor-specific neoantigen database (TSNAdb v1.0), based on pan-cancer immunogenomic analyses of somatic mutation data and human leukocyte antigen (HLA) allele information for 16 tumor types with 7748 tumor samples from The Cancer Genome Atlas (TCGA) and The Cancer Immunome Atlas (TCIA). We predicted binding affinities between mutant/wild-type peptides and HLA class I molecules by NetMHCpan v2.8/v4.0, and presented detailed information of 3,707,562/1,146,961 potential neoantigens generated by somatic mutations of all tumor samples. Moreover, we employed recurrent mutations in combination with highly frequent HLA alleles to predict potential shared neoantigens across tumor patients, which would facilitate the discovery of putative targets for neoantigen-based cancer immunotherapy. TSNAdb is freely available at http://biopharm.zju.edu.cn/tsnadb.     

IRESbase:A Comprehensive Database of Experimentally Validated Internal Ribosome Entry Sites

Internal ribosome entry sites (IRESs) are functional RNA elements that can directly recruit ribosomes to an internal position of the mRNA in a cap-independent manner to initiate translation. Recently, IRES elements have attracted much attention for their critical roles in various processes including translation initiation of a new type of RNA, circular RNA (circRNA), with no 5′ cap to support classical cap-dependent translation. Thus, an integrative data resource of IRES elements with experimental evidence will be useful for further studies. In this study, we present IRESbase, a comprehensive database of IRESs, by curating the experimentally validated functional minimal IRES elements from literature and annotating their host linear and circular RNAs. The current version of IRESbase contains 1328 IRESs, including 774 eukaryotic IRESs and 554 viral IRESs from 11 eukaryotic organisms and 198 viruses, respectively. As IRESbase collects only IRES of minimal length with functional evidence, the median length of IRESs in IRESbase is 174 nucleotides. By mapping IRESs to human circRNAs and long non-coding RNAs (lncRNAs), 2191 circRNAs and 168 lncRNAs were found to contain at least one entire or partial IRES sequence. IRESbase is available at http://reprod.njmu.edu.cn/cgi-bin/iresbase/index.php.     

HeteroMeth: A Database of Cell-to-cell Heterogeneity in DNA Methylation

DNA methylation is an important epigenetic mark that plays a vital role in gene expression and cell differentiation. The average DNA methylation level among a group of cells has been extensively documented. However, the cell-to-cell heterogeneity in DNA methylation, which reflects the differentiation of epigenetic status among cells, remains less investigated. Here we established a gold standard of the cell-to-cell heterogeneity in DNA methylation based on single-cell bisulfite sequencing (BS-seq) data. With that, we optimized a computational pipeline for estimating the heterogeneity in DNA methylation from bulk BS-seq data. We further built HeteroMeth, a database for searching, browsing, visualizing, and downloading the data for heterogeneity in DNA methylation for a total of 141 samples in humans, mice, Arabidopsis, and rice. Three genes are used as examples to illustrate the power of HeteroMeth in the identification of unique features in DNA methylation. The optimization of the computational strategy and the construction of the database in this study complement the recent experimental attempts on single-cell DNA methylomes and will facilitate the understanding of epigenetic mechanisms underlying cell differentiation and embryonic development. HeteroMeth is publicly available at http://qianlab.genetics.ac.cn/HeteroMeth.     

A Comprehensive Molecular Interaction Map for Rheumatoid Arthritis

Background

Computational biology contributes to a variety of areas related to life sciences and, due to the growing impact of translational medicine - the scientific approach to medicine in tight relation with basic science -, it is becoming an important player in clinical-related areas. In this study, we use computation methods in order to improve our understanding of the complex interactions that occur between molecules related to Rheumatoid Arthritis (RA).

Methodology

Due to the complexity of the disease and the numerous molecular players involved, we devised a method to construct a systemic network of interactions of the processes ongoing in patients affected by RA. The network is based on high-throughput data, refined semi-automatically with carefully curated literature-based information. This global network has then been topologically analysed, as a whole and tissue-specifically, in order to translate the experimental molecular connections into topological motifs meaningful in the identification of tissue-specific markers and targets in the diagnosis, and possibly in the therapy, of RA.

Significance

We find that some nodes in the network that prove to be topologically important, in particular AKT2, IL6, MAPK1 and TP53, are also known to be associated with drugs used for the treatment of RA. Importantly, based on topological consideration, we are also able to suggest CRKL as a novel potentially relevant molecule for the diagnosis or treatment of RA. This type of finding proves the potential of in silico analyses able to produce highly refined hypotheses, based on vast experimental data, to be tested further and more efficiently. As research on RA is ongoing, the present map is in fieri, despite being -at the moment- a reflection of the state of the art. For this reason we make the network freely available in the standardised and easily exportable .xml CellDesigner format at ‘www.picb.ac.cn/ClinicalGenomicNTW/temp.html’ and ‘www.celldesigner.org’.     

DEBKS: A Tool to Detect Differentially Expressed Circular RNAs

Circular RNAs (circRNAs) are involved in various biological processes and disease pathogenesis. However, only a small number of functional circRNAs have been identified among hundreds of thousands of circRNA species, partly because most current methods are based on circular junction counts and overlook the fact that a circRNA is formed from the host gene by back-splicing (BS). To distinguish the expression difference originating from BS or the host gene, we present differentially expressed back-splicing (DEBKS), a software program to streamline the discovery of differential BS events between two rRNA-depleted RNA sequencing (RNA-seq) sample groups. By applying to real and simulated data and employing RT-qPCR for validation, we demonstrate that DEBKS is efficient and accurate in detecting circRNAs with differential BS events between paired and unpaired sample groups. DEBKS is available at https://github.com/yangence/DEBKS as open-source software.     

Circular RNA expression profiles and features in human tissues: a study using RNA-seq data

Background

Circular RNA (circRNA) is one type of noncoding RNA that forms a covalently closed continuous loop. Similar to long noncoding RNA (lncRNA), circRNA can act as microRNA (miRNA) ‘sponges’ to regulate gene expression, and its abnormal expression is related to diseases such as atherosclerosis, nervous system disorders and cancer. So far, there have been no systematic studies on circRNA abundance and expression profiles in human adult and fetal tissues.

Results

We explored circRNA expression profiles using RNA-seq data for six adult and fetal normal tissues (colon, heart, kidney, liver, lung, and stomach) and four gland normal tissues (adrenal gland, mammary gland, pancreas, and thyroid gland). A total of 8120, 25,933 and 14,433 circRNAs were detected by at least two supporting junction reads in adult, fetal and gland tissues, respectively. Among them, 3092, 14,241 and 6879 circRNAs were novel when compared to the published results. In each adult tissue type, we found at least 1000 circRNAs, among which 36.97–50.04% were tissue-specific. We reported 33 circRNAs that were ubiquitously expressed in all the adult tissues we examined. To further explore the potential “housekeeping” function of these circRNAs, we constructed a circRNA-miRNA-mRNA regulatory network containing 17 circRNAs, 22 miRNAs and 90 mRNAs. Furthermore, we found that both the abundance and the relative expression level of circRNAs were higher in fetal tissue than adult tissue. The number of circRNAs in gland tissues, especially in mammary gland (9665 circRNA candidates), was higher than that of other adult tissues (1160–3777).

Conclusions

We systematically investigated circRNA expression in a variety of human adult and fetal tissues. Our observation of different expression level of circRNAs in adult and fetal tissues suggested that circRNAs might play their role in a tissue-specific and development-specific fashion. Analysis of circRNA-miRNA-mRNA network provided potential targets of circRNAs. High expression level of circRNAs in mammary gland might be attributed to the rich innervation.

     

Similarities in mRNA expression of peripheral blood drug metabolizing enzymes and cancer marker genes with biopsy samples of head and neck cancer patients

Purpose: To develop peripheral blood mRNA expression profiles of drug metabolizing enzymes (DMEs) as a surrogate to monitor tobacco induced head and neck squamous cell carcinoma (HNSCC), attempts were made to investigate (i) similarities in alterations with the cancer marker genes in biopsy samples and (ii) if alterations similar to that seen in biopsy samples are reflected in peripheral blood.

Methods: Total RNA from eight soft gingival tissues and eight biopsy samples of HNSCC patients and total DNA and RNA from blood of healthy controls (n?=?150) and HNSCC patients (n?=?150) was processed for expression and genotyping studies. Blood from patients receiving chemo-radiotherapy was processed for follow-up study.

Results: qRT-PCR revealed significant increase in mRNA expression of DMEs in biopsy and blood samples of HNSCC patients when compared to controls. Similar alterations were observed in cancer marker genes in these samples. Patients with variant genotypes of DMEs showed greater magnitude of alterations in mRNA expression when compared to wild type controls. Responders of chemo-radiotherapy showed significant decline in induction of mRNA expression of DMEs and cancer marker genes

Conclusions: The data suggest that peripheral blood expression profiles could be used to monitor tobacco-induced HNSCC as well as the treatment response.     

iKcr-PseEns: Identify lysine crotonylation sites in histone proteins with pseudo components and ensemble classifier

Lysine crotonylation (Kcr) is an evolution-conserved histone posttranslational modification (PTM), occurring in both human somatic and mouse male germ cell genomes. It is important for male germ cell differentiation. Information of Kcr sites in proteins is very useful for both basic research and drug development. But it is time-consuming and expensive to determine them by experiments alone. Here, we report a novel predictor called iKcr-PseEns that is established by incorporating five tiers of amino acid pairwise couplings into the general pseudo amino acid composition. It has been observed via rigorous cross-validations that the new predictor's sensitivity (Sn), specificity (Sp), accuracy (Acc), and stability (MCC) are 90.53%, 95.27%, 94.49%, and 0.826, respectively. For the convenience of most experimental scientists, a user-friendly web-server for iKcr-PseEns has been established at http://www.jci-bioinfo.cn/iKcr-PseEns, by which users can easily obtain their desired results without the need to go through the complicated mathematical equations involved.     

基层医疗机构人才队伍建设的探讨

??????? 在医改中硬件是基础,软件是根本,基层医疗机构人才队伍建设问题至关重要,本文就如何吸引毕业生下沉到基层、如何提升基层现存医疗队伍的技术水平提出建议,并为如何实现2020年培养30万名全科医生的总体目标,提出利用社会融资方法培养农村全科医生的构想。     

Histone 3 modifications and blood pressure in the Beijing Truck Driver Air Pollution Study

Context: Histone modifications regulate gene expression; dysregulation has been linked with cardiovascular diseases. Associations between histone modification levels and blood pressure in humans are unclear.

Objective: We examine the relationship between global histone concentrations and various markers of blood pressure.

Materials and methods: Using the Beijing Truck Driver Air Pollution Study, we investigated global peripheral white blood cell histone modifications (H3K9ac, H3K9me3, H3K27me3, and H3K36me3) associations with pre- and post-work measurements of systolic (SBP) and diastolic (DBP) blood pressure, mean arterial pressure (MAP), and pulse pressure (PP) using multivariable mixed-effect models.

Results: H3K9ac was negatively associated with pre-work SBP and MAP; H3K9me3 was negatively associated with pre-work SBP, DBP, and MAP; and H3K27me3 was negatively associated with pre-work SBP. Among office workers, H3K9me3 was negatively associated with pre-work SBP, DBP, and MAP. Among truck drivers, H3K9ac and H3K27me were negatively associated with pre-work SBP, and H3K27me3 was positively associated with post-work PP.

Discussion and conclusion: Epigenome-wide H3K9ac, H3K9me3, and H3K27me3 were negatively associated with multiple pre-work blood pressure measures. These associations substantially changed during the day, suggesting an influence of daily activities. Blood-based histone modification biomarkers are potential candidates for studies requiring estimations of morning/pre-work blood pressure.     

CA9 as a biomarker in preoperative biopsy of small solid renal masses for diagnosis of clear cell renal cell carcinoma

Purpose: The purpose of our study was to test the utility of CA9 expression in preoperative biopsy samples to identify the ccRCC.

Materials and methods: A total of 55 patients with a small solid renal mass (≤4?cm) entered into this study. The immunohistochemical staining (51 samples) and RT-PCR (33 samples) were performed to detect CA9 expression.

Results: For immunohistochemistry detection, CA9 was positive in 31/34 of biopsy samples of ccRCCs. CA9 was also positive in five suspected diagnosis of ccRCC. For RT-PCR detection, CA9 was positive in 25/25 biopsy samples of ccRCCs. The diagnostic accuracy of CA9 expression for ccRCC was 100%. RT-PCR was performed in four biopsies where immunohistochemistry could not be performed because of limited tissue materials or necrosis. Two ccRCCs with a negative staining by immunohistochemistry were CA9 positive by RT-PCR.

Conclusions: CA9 may be potentially useful biomarker in help making a diagnosis of ccRCC in the biopsy of renal mass. RT-PCR might be a preferred method to immunohistochemistry for the detection of CA9 in renal biopsy samples.     

The Genome Sequence Archive Family: Toward Explosive Data Growth and Diverse Data Types

The Genome Sequence Archive (GSA) is a data repository for archiving raw sequence data, which provides data storage and sharing services for worldwide scientific communities. Considering explosive data growth with diverse data types, here we present the GSA family by expanding into a set of resources for raw data archive with different purposes, namely, GSA (https://ngdc.cncb.ac.cn/gsa/), GSA for Human (GSA-Human, https://ngdc.cncb.ac.cn/gsa-human/), and Open Archive for Miscellaneous Data (OMIX, https://ngdc.cncb.ac.cn/omix/). Compared with the 2017 version, GSA has been significantly updated in data model, online functionalities, and web interfaces. GSA-Human, as a new partner of GSA, is a data repository specialized in human genetics-related data with controlled access and security. OMIX, as a critical complement to the two resources mentioned above, is an open archive for miscellaneous data. Together, all these resources form a family of resources dedicated to archiving explosive data with diverse types, accepting data submissions from all over the world, and providing free open access to all publicly available data in support of worldwide research activities.