A Genome-Wide Investigation of Effects of Aberrant DNA Methylation on the Usage of Alternative Promoters in Hepatocellular Carcinoma

BackgroundThe alternative usage of promoters provides a way to regulate gene expression, has a significant influence on the transcriptome, and contributes to the cellular transformation of cancer. However, the function of alternative promoters (APs) in hepatocellular carcinoma (HCC) has not been systematically studied yet. In addition, the potential mechanism of regulation to the usage of APs remains unclear. DNA methylation, one of the most aberrant epigenetic modifications in cancers, is known to regulate transcriptional activity. Whether DNA methylation regulates the usage of APs needs to be explored. Here, we aim to investigate the effects of DNA methylation on usage of APs in HCC.MethodsPromoter activities were calculated based on RNA-seq data. Functional enrichment analysis was implemented to conduct GO terms. Correlation tests were used to detect the correlation between promoter activity and methylation status. The LASSO regression model was used to generate a diagnostic model. Kaplan–Meier analysis was used to compare the overall survival between high and low methylation groups. RNA-seq and whole-genome bisulfite sequencing (WGBS) in HCC samples were performed to validate the correlation of promoter activity and methylation.ResultsWe identified 855 APs in total, which could be well used to distinguish cancer from normal samples. The correlation of promoter activity and DNA methylation in APs was observed, and the APs with negative correlation were defined as methylation-regulated APs (mrAPs). Six mrAPs were identified to generate a diagnostic model with good performance (AUC = 0.97). Notably, the majority of mrAPs had CpG sites that could be used to predict clinical outcomes by methylation status. Finally, we verified 85.6% of promoter activity variation and 92.3% of methylation changes in our paired RNA-seq and WGBS samples, respectively. The negative correlation between promoter activity and methylation status was further confirmed in our HCC samples.ConclusionThe aberrant methylation status plays a critical role in the precision usage of APs in HCC, which sheds light on the mechanism of cancer development and provides a new insight into cancer screening and treatment.

Integrating RNA-Seq With GWAS Reveals a Novel SNP in Immune-Related HLA-DQB1 Gene Associated With Occupational Pulmonary Fibrosis Risk: A Multi-Stage Study

ObjectiveTo evaluate the association between single-nucleotide polymorphisms (SNPs) in RNA-seq identified mRNAs and silicosis susceptibility.MethodsA comprehensive RNA-seq was performed to screen for differently expressed mRNAs in the peripheral blood lymphocytes of eight subjects exposed to silica dust (four silicosis cases and four healthy controls). Following this, the SNPs located on the shortlisted mRNAs, which may affect silicosis susceptibility, were screened through silicosis-related genome-wide association studies (GWAS) (155 silicosis cases and 141 healthy controls), whereas functional expression quantitative trait locus (eQTL)-SNPs were identified using the GTEx database. Finally, the association between functional eQTL-SNPs and silicosis susceptibility (194 silicosis cases and 235 healthy controls) was validated.ResultsA total of 70 differentially expressed mRNAs (fold change > 2 or fold change < 0.5, P < 0.05) was obtained using RNA-seq. Furthermore, 476 SNPs located on the shortlisted mRNAs, which may affect silicosis susceptibility (P < 0.05) were obtained using GWAS, whereas subsequent six functional eQTL-SNPs were identified. The mutant A allele of rs9273410 in HLA-DQB1 indicated a potential increase in silicosis susceptibility in the validation stage (additive model: odds ratio (OR)= 1.31, 95% confidence interval (CI) = 0.99–1.74, P = 0.061), whereas the combination of GWAS and the validation results indicated that the mutant A allele of rs9273410 was associated with increased silicosis susceptibility (additive model: OR = 1.35, 95% CI =1.09–1.68, P = 0.006).ConclusionThe mutant A allele of rs9273410 was associated with increased silicosis susceptibility by modulating the expression of HLA-DQB1.

Vexitoxins: a novel class of conotoxin-like venom peptides from predatory gastropods of the genus Vexillum

Venoms of predatory marine cone snails (the family Conidae, order Neogastropoda) are intensely studied because of the broad range of biomedical applications of the neuropeptides that they contain, conotoxins. Meanwhile anatomy in some other neogastropod lineages strongly suggests that they have evolved similar venoms independently of cone snails, nevertheless their venom composition remains unstudied. Here we focus on the most diversified of these lineages, the genus Vexillum (the family Costellariidae). We have generated comprehensive multi-specimen, multi-tissue RNA-Seq data sets for three Vexillum species, and supported our findings in two species by proteomic profiling. We show that venoms of Vexillum are dominated by highly diversified short cysteine-rich peptides that in many aspects are very similar to conotoxins. Vexitoxins possess the same precursor organization, display overlapping cysteine frameworks and share several common post-translational modifications with conotoxins. Some vexitoxins show detectable sequence similarity to conotoxins, and are predicted to adopt similar domain conformations, including a pharmacologically relevant inhibitory cysteine-know motif (ICK). The tubular gL of Vexillum is a notably more recent evolutionary novelty than the conoidean venom gland. Thus, we hypothesize lower divergence between the toxin genes, and their somatic counterparts compared to that in conotoxins, and we find support for this hypothesis in the molecular evolution of the vexitoxin cluster V027. We use this example to discuss how future studies on vexitoxins can inform origin and evolution of conotoxins, and how they may help addressing standing questions in venom evolution.

Age Effects Aggressive Behavior: RNA-Seq Analysis in Cattle with Implications for Studying Neoteny Under Domestication

The reactive type of aggression is regulated mostly by the brain’s prefrontal cortex; however, the molecular changes underlying aggressiveness in adults have not been fully characterized. We used an RNA-seq approach to investigate differential gene expression in the prefrontal cortex of bovines from the aggressive Lidia breed at different ages: young three-year old and adult four-year-old bulls. A total of 50 up and 193 down-regulated genes in the adult group were identified. Furthermore, a cross-species comparative analysis retrieved 29 genes in common with previous studies on aggressive behaviors, representing an above-chance overlap with the differentially expressed genes in adult bulls. We detected changes in the regulation of networks such as synaptogenesis, involved in maintenance and refinement of synapses, and the glutamate receptor pathway, which acts as excitatory driver in aggressive responses. The reduced reactive aggression typical of domestication has been proposed to form part of a retention of juvenile traits as adults (neoteny).

Genome-Wide Identification and Analysis of Cell Cycle Genes in Birch

Research Highlights: This study identified the cell cycle genes in birch that likely play important roles during the plant’s growth and development. This analysis provides a basis for understanding the regulatory mechanism of various cell cycles in Betula pendula Roth. Background and Objectives: The cell cycle factors not only influence cell cycles progression together, but also regulate accretion, division, and differentiation of cells, and then regulate growth and development of the plant. In this study, we identified the putative cell cycle genes in the B. pendula genome, based on the annotated cell cycle genes in Arabidopsis thaliana (L.) Heynh. It can be used as a basis for further functional research. Materials and Methods: RNA-seq technology was used to determine the transcription abundance of all cell cycle genes in xylem, roots, leaves, and floral tissues. Results: We identified 59 cell cycle gene models in the genome of B. pendula, with 17 highly expression genes among them. These genes were BpCDKA.1, BpCDKB1.1, BpCDKB2.1, BpCKS1.2, BpCYCB1.1, BpCYCB1.2, BpCYCB2.1, BpCYCD3.1, BpCYCD3.5, BpDEL1, BpDpa2, BpE2Fa, BpE2Fb, BpKRP1, BpKRP2, BpRb1, and BpWEE1. Conclusions: By combining phylogenetic analysis and tissue-specific expression data, we identified 17 core cell cycle genes in the Betulapendula genome.

QTL Mapping Low-Temperature Germination Ability in the Maize IBM Syn10 DH Population

Chilling injury poses a serious threat to seed emergence of spring-sowing maize in China, which has become one of the main climatic limiting factors affecting maize production in China. It is of great significance to mine the key genes controlling low-temperature tolerance during seed germination and study their functions for breeding new maize varieties with strong low-temperature tolerance during germination. In this study, 176 lines of the intermated B73 × Mo17 (IBM) Syn10 doubled haploid (DH) population, which comprised 6618 bin markers, were used for QTL analysis of low-temperature germination ability. The results showed significant differences in germination related traits under optimum-temperature condition (25 °C) and low-temperature condition (10 °C) between two parental lines. In total, 13 QTLs were detected on all chromosomes, except for chromosome 5, 7, 10. Among them, seven QTLs formed five QTL clusters on chromosomes 1, 2, 3, 4, and 9 under the low-temperature condition, which suggested that there may be some genes regulating multiple germination traits at the same time. A total of 39 candidate genes were extracted from five QTL clusters based on the maize GDB under the low-temperature condition. To further screen candidate genes controlling low-temperature germination, RNA-Seq, in which RNA was extracted from the germination seeds of B73 and Mo17 at 10 °C, was conducted, and three B73 upregulated genes and five Mo17 upregulated genes were found by combined analysis of RNA-Seq and QTL located genes. Additionally, the variations of Zm00001d027976 (GLABRA2), Zm00001d007311 (bHLH transcription factor), and Zm00001d053703 (bZIP transcription factor) were found by comparison of amino sequence between B73 and Mo17. This study will provide a theoretical basis for marker-assisted breeding and lay a foundation for further revealing molecular mechanism of low-temperature germination tolerance in maize.