The UCSC Genome Browser database: update 2010

Heterologous survey of 130 DNA transposons in human cells highlights their functional divergence and expands the genome engineering toolbox

Experimental studies on DNA transposable elements (TEs) have been limited in scale, leading to a lack of understanding of the factors influencing transposition activity, evolutionary dynamics, and application potential as genome engineering tools. We predicted 130 active DNA TEs from 102 metazoan genomes and evaluated their activity in human cells. We identified 40 active (integration-competent) TEs, surpassing the cumulative number (20) of TEs found previously. With this unified comparative data, we found that the Tc1/mariner superfamily exhibits elevated activity, potentially explaining their pervasive horizontal transfers. Further functional characterization of TEs revealed additional divergence in features such as insertion bias. Remarkably, in CAR-T therapy for hematological and solid tumors, Mariner2_AG (MAG), the most active DNA TE identified, largely outperformed two widely used vectors, the lentiviral vector and the TE-based vector SB100X. Overall, this study highlights the varied transposition features and evolutionary dynamics of DNA TEs and increases the TE toolbox diversity.

Unraveling the role of cloud computing in health care system and biomedical sciences

Cloud computing has emerged as a transformative force in healthcare and biomedical sciences, offering scalable, on-demand resources for managing vast amounts of data. This review explores the integration of cloud computing within these fields, highlighting its pivotal role in enhancing data management, security, and accessibility. We examine the application of cloud computing in various healthcare domains, including electronic medical records, telemedicine, and personalized patient care, as well as its impact on bioinformatics research, particularly in genomics, proteomics, and metabolomics. The review also addresses the challenges and ethical considerations associated with cloud-based healthcare solutions, such as data privacy and cybersecurity. By providing a comprehensive overview, we aim to assist readers in understanding the significance of cloud computing in modern medical applications and its potential to revolutionize both patient care and biomedical research.

Synthetic reversed sequences reveal default genomic states

Pervasive transcriptional activity is observed across diverse species. The genomes of extant organisms have undergone billions of years of evolution, making it unclear whether these genomic activities represent effects of selection or 'noise'1-4. Characterizing default genome states could help understand whether pervasive transcriptional activity has biological meaning. Here we addressed this question by introducing a synthetic 101-kb locus into the genomes of Saccharomyces cerevisiae and Mus musculus and characterizing genomic activity. The locus was designed by reversing but not complementing human HPRT1, including its flanking regions, thus retaining basic features of the natural sequence but ablating evolved coding or regulatory information. We observed widespread activity of both reversed and native HPRT1 loci in yeast, despite the lack of evolved yeast promoters. By contrast, the reversed locus displayed no activity at all in mouse embryonic stem cells, and instead exhibited repressive chromatin signatures. The repressive signature was alleviated in a locus variant lacking CpG dinucleotides; nevertheless, this variant was also transcriptionally inactive. These results show that synthetic genomic sequences that lack coding information are active in yeast, but inactive in mouse embryonic stem cells, consistent with a major difference in 'default genomic states' between these two divergent eukaryotic cell types, with implications for understanding pervasive transcription, horizontal transfer of genetic information and the birth of new genes.

Construction of single-cell cross-species chromatin accessibility landscapes with combinatorial-hybridization-based ATAC-seq

Despite recent advances in single-cell genomics, the lack of maps for single-cell candidate cis-regulatory elements (cCREs) in non-mammal species has limited our exploration of conserved regulatory programs across vertebrates and invertebrates. Here, we developed a combinatorial-hybridization-based method for single-cell assay for transposase-accessible chromatin using sequencing (scATAC-seq) named CH-ATAC-seq, enabling the construction of single-cell accessible chromatin landscapes for zebrafish, Drosophila, and earthworms (Eisenia andrei). By integrating scATAC censuses of humans, monkeys, and mice, we systematically identified 152 distinct main cell types and around 0.8 million cell-type-specific cCREs. Our analysis provided insights into the conservation of neural, muscle, and immune lineages across species, while epithelial cells exhibited a higher organ-origin heterogeneity. Additionally, a large-scale gene regulatory network (GRN) was constructed in four vertebrates by integrating scRNA-seq censuses. Overall, our study provides a valuable resource for comparative epigenomics, identifying the evolutionary conservation and divergence of gene regulation across different species.

A Guide to Genome-Wide Association Study Design for Diabetic Retinopathy

Diabetic retinopathy (DR) is the most common microvascular complication related to diabetes. There is evidence that genetics play an important role in DR pathogenesis, but the complexity of the disease makes genetic studies a challenge. This chapter is a practical overview of the basic steps for genome-wide association studies with respect to DR and its associated traits. Also described are approaches that can be adopted in future DR studies. This is intended to serve as a guide for beginners and to provide a framework for further in-depth analysis.

Folic Acid and Vitamin B12 Prevent Deleterious Effects of Rotenone on Object Novelty Recognition Memory and Kynu Expression in an Animal Model of Parkinson's Disease

Parkinson's disease (PD) is characterized by a range of motor signs, but cognitive dysfunction is also observed. Supplementation with folic acid and vitamin B12 is expected to prevent cognitive impairment. To test this in PD, we promoted a lesion within the substantia nigra pars compacta of rats using the neurotoxin rotenone. In the sequence, the animals were supplemented with folic acid and vitamin B12 for 14 consecutive days and subjected to the object recognition test. We observed an impairment in object recognition memory after rotenone administration, which was prevented by supplementation (p < 0.01). Supplementation may adjust gene expression through efficient DNA methylation. To verify this, we measured the expression and methylation of the kynureninase gene (Kynu), whose product metabolizes neurotoxic metabolites often accumulated in PD as kynurenine. Supplementation prevented the decrease in Kynu expression induced by rotenone in the substantia nigra (p < 0.05), corroborating the behavioral data. No differences were observed concerning the methylation analysis of two CpG sites in the Kynu promoter. Instead, we suggest that folic acid and vitamin B12 increased global DNA methylation, reduced the expression of Kynu inhibitors, maintained Kynu-dependent pathway homeostasis, and prevented the memory impairment induced by rotenone. Our study raises the possibility of adjuvant therapy for PD with folic acid and vitamin B12.

Domestic dog lineages reveal genetic drivers of behavioral diversification

Selective breeding of domestic dogs has generated diverse breeds often optimized for performing specialized tasks. Despite the heritability of breed-typical behavioral traits, identification of causal loci has proven challenging due to the complexity of canine population structure. We overcome longstanding difficulties in identifying genetic drivers of canine behavior by developing a framework for understanding relationships between breeds and the behaviors that define them, utilizing genetic data for over 4,000 domestic, semi-feral, and wild canids and behavioral survey data for over 46,000 dogs. We identify ten major canine genetic lineages and their behavioral correlates and show that breed diversification is predominantly driven by non-coding regulatory variation. We determine that lineage-associated genes converge in neurodevelopmental co-expression networks, identifying a sheepdog-associated enrichment for interrelated axon guidance functions. This work presents a scaffold for canine diversification that positions the domestic dog as an unparalleled system for revealing the genetic origins of behavioral diversity.

Pan-cancer surveys indicate cell cycle-related roles of primate-specific genes in tumors and embryonic cerebrum

BACKGROUND

Despite having been extensively studied, it remains largely unclear why humans bear a particularly high risk of cancer. The antagonistic pleiotropy hypothesis predicts that primate-specific genes (PSGs) tend to promote tumorigenesis, while the molecular atavism hypothesis predicts that PSGs involved in tumors may represent recently derived duplicates of unicellular genes. However, these predictions have not been tested.

RESULTS

By taking advantage of pan-cancer genomic data, we find the upregulation of PSGs across 13 cancer types, which is facilitated by copy-number gain and promoter hypomethylation. Meta-analyses indicate that upregulated PSGs (uPSGs) tend to promote tumorigenesis and to play cell cycle-related roles. The cell cycle-related uPSGs predominantly represent derived duplicates of unicellular genes. We prioritize 15 uPSGs and perform an in-depth analysis of one unicellular gene-derived duplicate involved in the cell cycle, DDX11. Genome-wide screening data and knockdown experiments demonstrate that DDX11 is broadly essential across cancer cell lines. Importantly, non-neutral amino acid substitution patterns and increased expression indicate that DDX11 has been under positive selection. Finally, we find that cell cycle-related uPSGs are also preferentially upregulated in the highly proliferative embryonic cerebrum.

CONCLUSIONS

Consistent with the predictions of the atavism and antagonistic pleiotropy hypotheses, primate-specific genes, especially those PSGs derived from cell cycle-related genes that emerged in unicellular ancestors, contribute to the early proliferation of the human cerebrum at the cost of hitchhiking by similarly highly proliferative cancer cells.

An updated overview of experimental and computational approaches to identify non-canonical DNA/RNA structures with emphasis on G-quadruplexes and R-loops

Multiple types of non-canonical nucleic acid structures play essential roles in DNA recombination and replication, transcription, and genomic instability and have been associated with several human diseases. Thus, an increasing number of experimental and bioinformatics methods have been developed to identify these structures. To date, most reviews have focused on the features of non-canonical DNA/RNA structure formation, experimental approaches to mapping these structures, and the association of these structures with diseases. In addition, two reviews of computational algorithms for the prediction of non-canonical nucleic acid structures have been published. One of these reviews focused only on computational approaches for G4 detection until 2020. The other mainly summarized the computational tools for predicting cruciform, H-DNA and Z-DNA, in which the algorithms discussed were published before 2012. Since then, several experimental and computational methods have been developed. However, a systematic review including the conformation, sequencing mapping methods and computational prediction strategies for these structures has not yet been published. The purpose of this review is to provide an updated overview of conformation, current sequencing technologies and computational identification methods for non-canonical nucleic acid structures, as well as their strengths and weaknesses. We expect that this review will aid in understanding how these structures are characterised and how they contribute to related biological processes and diseases.

The Innovative Informatics Approaches of High-Throughput Technologies in Livestock: Spearheading the Sustainability and Resiliency of Agrigenomics Research

For more than a decade, next-generation sequencing (NGS) has been emerging as the mainstay of agrigenomics research. High-throughput technologies have made it feasible to facilitate research at the scale and cost required for using this data in livestock research. Scale frameworks of sequencing for agricultural and livestock improvement, management, and conservation are partly attributable to innovative informatics methodologies and advancements in sequencing practices. Genome-wide sequence-based investigations are often conducted worldwide, and several databases have been created to discover the connections between worldwide scientific accomplishments. Such studies are beginning to provide revolutionary insights into a new era of genomic prediction and selection capabilities of various domesticated livestock species. In this concise review, we provide selected examples of the current state of sequencing methods, many of which are already being used in animal genomic studies, and summarize the state of the positive attributes of genome-based research for cattle (Bos taurus), sheep (Ovis aries), pigs (Sus scrofa domesticus), horses (Equus caballus), chickens (Gallus gallus domesticus), and ducks (Anas platyrhyncos). This review also emphasizes the advantageous features of sequencing technologies in monitoring and detecting infectious zoonotic diseases. In the coming years, the continued advancement of sequencing technologies in livestock agrigenomics will significantly influence the sustained momentum toward regulatory approaches that encourage innovation to ensure continued access to a safe, abundant, and affordable food supplies for future generations.

Distinct phosphorylation states of mammalian CaMKIIβ control the induction and maintenance of sleep

The reduced sleep duration previously observed in Camk2b knockout mice revealed a role for Ca2+/calmodulin-dependent protein kinase II (CaMKII)β as a sleep-promoting kinase. However, the underlying mechanism by which CaMKIIβ supports sleep regulation is largely unknown. Here, we demonstrate that activation or inhibition of CaMKIIβ can increase or decrease sleep duration in mice by almost 2-fold, supporting the role of CaMKIIβ as a core sleep regulator in mammals. Importantly, we show that this sleep regulation depends on the kinase activity of CaMKIIβ. A CaMKIIβ mutant mimicking the constitutive-active (auto)phosphorylation state promotes the transition from awake state to sleep state, while mutants mimicking subsequent multisite (auto)phosphorylation states suppress the transition from sleep state to awake state. These results suggest that the phosphorylation states of CaMKIIβ differently control sleep induction and maintenance processes, leading us to propose a "phosphorylation hypothesis of sleep" for the molecular control of sleep in mammals.

Infection-induced membrane ruffling initiates danger and immune signaling via the mechanosensor PIEZO1

Microorganisms are generally sensed by receptors recognizing microbial molecules, which evoke changes in cellular activities and gene expression. Bacterial pathogens induce secretion of the danger signal ATP as an early alert response of intestinal epithelial cells, initiating overt inflammation. However, what triggers ATP secretion during infection is unclear. Here we show that the inherently mechanosensitive plasma membrane channel PIEZO1 acts as a sensor for bacterial entry. PIEZO1 is mechanically activated by invasion-induced membrane ruffles upstream of Ca²⁺ influx and ATP secretion. Mimicking mechanical stimuli of pathogen uptake with sterile beads equally elicits ATP secretion. Chemical or genetic PIEZO1 inactivation inhibits mechanically induced ATP secretion. Moreover, chemical or mechanical PIEZO1 activation evokes gene expression in immune and barrier pathways. Thus, mechanosensation of invasion-induced plasma membrane distortion initiates immune signaling upon infection, independently of detection of microbial molecules. Hence, PIEZO1-dependent detection of infection is driven by physical signals instead of chemical ligands.

eRNAs Identify Immune Microenvironment Patterns and Provide a Novel Prognostic Tool in Acute Myeloid Leukemia

Background: Enhancer RNAs (eRNAs) play an essential role in tumorigenesis as non-coding RNAs transcribed from enhancer regions. However, the landscape of eRNAs in acute myeloid leukemia (AML) and the potential roles of eRNAs in the tumor microenvironment (TME) remain unclear.

Method: Gene expression data collected from The Cancer Genome Atlas (TCGA) project were combined with Histone ChIP-seq so as to reveal the comprehensive landscape of eRNAs. Single-sample gene set enrichment analysis algorithm (ssGSEA) and ESTIMATE were employed to enumerate immune cell infiltration and tumor purity.

Results: Most prognostic eRNAs were enriched in immune-related pathways. Two distinct immune microenvironment patterns, the immune-active subtype and the immune-resistant subtype, were identified in AML. We further developed an eRNA-derived score (E-score) that could quantify immune microenvironment patterns and predict the response to immune checkpoint inhibitor (ICI) treatment. Finally, we established a prognostic nomogram combining E-score and other clinical features, which showed great discriminative power in both the training set [Harrell’s concordance index (C index): 0.714 (0.651–0.777), p < 0.0001] and validation set [C index: 0.684 (0.614–0.755), p < 0.0001]. Calibration of the nomogram was also validated independently.

Conclusion: In this study, we systematically understood the roles of eRNAs in regulating TME diversity and complexity. Moreover, our E-score model provided the first predictive model for ICI treatment in AML.

Modeling uniquely human gene regulatory function via targeted humanization of the mouse genome

The evolution of uniquely human traits likely entailed changes in developmental gene regulation. Human Accelerated Regions (HARs), which include transcriptional enhancers harboring a significant excess of human-specific sequence changes, are leading candidates for driving gene regulatory modifications in human development. However, insight into whether HARs alter the level, distribution, and timing of endogenous gene expression remains limited. We examined the role of the HAR HACNS1 (HAR2) in human evolution by interrogating its molecular functions in a genetically humanized mouse model. We find that HACNS1 maintains its human-specific enhancer activity in the mouse embryo and modifies expression of Gbx2, which encodes a transcription factor, during limb development. Using single-cell RNA-sequencing, we demonstrate that Gbx2 is upregulated in the limb chondrogenic mesenchyme of HACNS1 homozygous embryos, supporting that HACNS1 alters gene expression in cell types involved in skeletal patterning. Our findings illustrate that humanized mouse models provide mechanistic insight into how HARs modified gene expression in human evolution.

Dosage sensitivity and exon shuffling shape the landscape of polymorphic duplicates in Drosophila and humans

Despite polymorphic duplicate genes' importance for the early stages of duplicate gene evolution, they are less studied than old gene duplicates. Two essential questions thus remain poorly addressed: how does dosage sensitivity, imposed by stoichiometry in protein complexes or by X chromosome dosage compensation, affect the emergence of complete duplicate genes? Do introns facilitate intergenic and intragenic chimaerism as predicted by the theory of exon shuffling? Here, we analysed new data for Drosophila and public data for humans, to characterize polymorphic duplicate genes with respect to dosage, exon-intron structures and allele frequencies. We found that complete duplicate genes are under dosage constraint induced by protein stoichiometry but potentially tolerated by X chromosome dosage compensation. We also found that in the intron-rich human genome, gene fusions and intragenic duplications extensively use intronic breakpoints generating in-frame proteins, in accordance with the theory of exon shuffling. Finally, we found that only a small proportion of complete or partial duplicates are at high frequencies, indicating the deleterious nature of dosage or gene structural changes. Altogether, we demonstrate how mechanistic factors including dosage sensitivity and exon-intron structure shape the short-term functional consequences of gene duplication.

Gene Regulation Analysis Reveals Perturbations of Autism Spectrum Disorder during Neural System Development

Autism spectrum disorder (ASD) is a neurodevelopmental disorder that impedes patients' cognition, social, speech and communication skills. ASD is highly heterogeneous with a variety of etiologies and clinical manifestations. The prevalence rate of ASD increased steadily in recent years. Presently, molecular mechanisms underlying ASD occurrence and development remain to be elucidated. Here, we integrated multi-layer genomics data to investigate the transcriptome and pathway dysregulations in ASD development. The RNA sequencing (RNA-seq) expression profiles of induced pluripotent stem cells (iPSCs), neural progenitor cells (NPCs) and neuron cells from ASD and normal samples were compared in our study. We found that substantially more genes were differentially expressed in the NPCs than the iPSCs. Consistently, gene set variation analysis revealed that the activity of the known ASD pathways in NPCs and neural cells were significantly different from the iPSCs, suggesting that ASD occurred at the early stage of neural system development. We further constructed comprehensive brain- and neural-specific regulatory networks by incorporating transcription factor (TF) and gene interactions with long 5 non-coding RNA(lncRNA) and protein interactions. We then overlaid the transcriptomes of different cell types on the regulatory networks to infer the regulatory cascades. The variations of the regulatory cascades between ASD and normal samples uncovered a set of novel disease-associated genes and gene interactions, particularly highlighting the functional roles of ELF3 and the interaction between STAT1 and lncRNA ELF3-AS 1 in the disease development. These new findings extend our understanding of ASD and offer putative new therapeutic targets for further studies.

Factors Influencing Green Purchase Intention: Moderating Role of Green Brand Knowledge

The current study aims to investigate the moderating effect of green brand knowledge (GBK) on the relationship of green brand positioning (GBP), attitude towards the green brand (ATGB), environmental concern (EC) and green purchase intention (GPI) in Pakistan. For this purpose, the data was collected from the individuals who were buying organic food by using purposive sampling, using cross-sectional research design and quantitative research approach. The Partial Least Square (PLS)-Structural Equation Modeling (SEM) technique results had shown that all the direct-effect relationships, namely, GBP, ATGB, EC variables have a positive and significant relationship with the GPI. While indirect-effect relationships have shown that the relationships of ATGB, EC and GPI are significantly moderated by GBK, which indicated that the effect of GBP, and EC toward GPI would be stronger when individuals have strong knowledge about green brands. In contrast, GBK is not significantly moderating the relationship between GBP and GPI. The empirical findings of this study fill a gap in the existing body of literature regarding the effects of GPI, ATGB and EC on green brands, as well as the moderating effect of GBK. As a result, this study provides insight into the topic, which has not been thoroughly investigated in earlier studies. Therefore, we consider that understanding this moderating effect is a positive contribution to the existing body of knowledge, which could help researchers explore this relationship in the future. This study could also help the owners and managers to know about the importance of these exogenous, and moderate variables to increase their customer’s green purchase intentions.

Hepatocellular adenoma in men: A nationwide assessment of pathology and correlation with clinical course

BACKGROUND & AIMS

Hepatocellular adenomas (HCA) rarely occur in males, and if so, are frequently associated with malignant transformation. Guidelines are based on small numbers of patients and advise resection of HCA in male patients, irrespective of size or subtype. This nationwide retrospective cohort study is the largest series of HCA in men correlating (immuno)histopathological and molecular findings with the clinical course.

METHODS

Dutch male patients with available histological slides with a (differential) diagnosis of HCA between 2000 and 2017 were identified through the Dutch Pathology Registry (PALGA). Histopathology and immunohistochemistry according to international guidelines were revised by two expert hepatopathologists. Next generation sequencing (NGS) was performed to confirm hepatocellular carcinoma (HCC) and/or subtype HCA. Final pathological diagnosis was correlated with recurrence, metastasis and death.

RESULTS

A total of 66 patients from 26 centres fulfilling the inclusion criteria with a mean (±SD) age of 45.0 ± 21.6 years were included. The diagnosis was changed after expert revision and NGS in 33 of the 66 patients (50%). After a median follow-up of 9.6 years, tumour-related mortality of patients with accessible clinical data was 1/18 (5.6%) in HCA, 5/14 (35.7%) in uncertain HCA/HCC and 4/9 (44.4%) in the HCC groups (P = .031). Four B-catenin mutated HCA were identified using NGS, which were not yet identified by immunohistochemistry and expert revision.

CONCLUSIONS

Expert revision with relevant immunohistochemistry may help the challenging but prognostically relevant distinction between HCA and well-differentiated HCC in male patients. NGS may be more important to subtype HCA than indicated in present guidelines.

DNA transposons mediate duplications via transposition-independent and -dependent mechanisms in metazoans

Despite long being considered as "junk", transposable elements (TEs) are now accepted as catalysts of evolution. One example is Mutator-like elements (MULEs, one type of terminal inverted repeat DNA TEs, or TIR TEs) capturing sequences as Pack-MULEs in plants. However, their origination mechanism remains perplexing, and whether TIR TEs mediate duplication in animals is almost unexplored. Here we identify 370 Pack-TIRs in 100 animal reference genomes and one Pack-TIR (Ssk-FB4) family in fly populations. We find that single-copy Pack-TIRs are mostly generated via transposition-independent gap filling, and multicopy Pack-TIRs are likely generated by transposition after replication fork switching. We show that a proportion of Pack-TIRs are transcribed and often form chimeras with hosts. We also find that Ssk-FB4s represent a young protein family, as supported by proteomics and signatures of positive selection. Thus, TIR TEs catalyze new gene structures and new genes in animals via both transposition-independent and -dependent mechanisms.

A systematic review of computational methods for predicting long noncoding RNAs

Accurately and rapidly distinguishing long noncoding RNAs (lncRNAs) from transcripts is prerequisite for exploring their biological functions. In recent years, many computational methods have been developed to predict lncRNAs from transcripts, but there is no systematic review on these computational methods. In this review, we introduce databases and features involved in the development of computational prediction models, and subsequently summarize existing state-of-the-art computational methods, including methods based on binary classifiers, deep learning and ensemble learning. However, a user-friendly way of employing existing state-of-the-art computational methods is in demand. Therefore, we develop a Python package ezLncPred, which provides a pragmatic command line implementation to utilize nine state-of-the-art lncRNA prediction methods. Finally, we discuss challenges of lncRNA prediction and future directions.

Vibrio cholerae cytotoxin MakA induces noncanonical autophagy resulting in the spatial inhibition of canonical autophagy

Autophagy plays an essential role in the defense against many microbial pathogens as a regulator of both innate and adaptive immunity. Some pathogens have evolved sophisticated mechanisms that promote their ability to evade or subvert host autophagy. Here, we describe a novel mechanism of autophagy modulation mediated by the recently discovered Vibrio cholerae cytotoxin, motility-associated killing factor A (MakA). pH-dependent endocytosis of MakA by host cells resulted in the formation of a cholesterol-rich endolysosomal membrane aggregate in the perinuclear region. Aggregate formation induced the noncanonical autophagy pathway driving unconventional LC3 (herein referring to MAP1LC3B) lipidation on endolysosomal membranes. Subsequent sequestration of the ATG12-ATG5-ATG16L1 E3-like enzyme complex, required for LC3 lipidation at the membranous aggregate, resulted in an inhibition of both canonical autophagy and autophagy-related processes, including the unconventional secretion of interleukin-1β (IL-1β). These findings identify a novel mechanism of host autophagy modulation and immune modulation employed by V. cholerae during bacterial infection.

Genome-wide analysis of pseudogenes reveals HBBP1's human-specific essentiality in erythropoiesis and implication in β-thalassemia

The human genome harbors 14,000 duplicated or retroposed pseudogenes. Given their functionality as regulatory RNAs and low conservation, we hypothesized that pseudogenes could shape human-specific phenotypes. To test this, we performed co-expression analyses and found that pseudogene exhibited tissue-specific expression, especially in the bone marrow. By incorporating genetic data, we identified a bone-marrow-specific duplicated pseudogene, HBBP1 (η-globin), which has been implicated in β-thalassemia. Extensive functional assays demonstrated that HBBP1 is essential for erythropoiesis by binding the RNA-binding protein (RBP), HNRNPA1, to upregulate TAL1, a key regulator of erythropoiesis. The HBBP1/TAL1 interaction contributes to a milder symptom in β-thalassemia patients. Comparative studies further indicated that the HBBP1/TAL1 interaction is human-specific. Genome-wide analyses showed that duplicated pseudogenes are often bound by RBPs and less commonly bound by microRNAs compared with retropseudogenes. Taken together, we not only demonstrate that pseudogenes can drive human evolution but also provide insights on their functional landscapes.

HRT Atlas v1.0 database: redefining human and mouse housekeeping genes and candidate reference transcripts by mining massive RNA-seq datasets

Housekeeping (HK) genes are constitutively expressed genes that are required for the maintenance of basic cellular functions. Despite their importance in the calibration of gene expression, as well as the understanding of many genomic and evolutionary features, important discrepancies have been observed in studies that previously identified these genes. Here, we present Housekeeping and Reference Transcript Atlas (HRT Atlas v1.0, www.housekeeping.unicamp.br) a web-based database which addresses some of the previously observed limitations in the identification of these genes, and offers a more accurate database of human and mouse HK genes and transcripts. The database was generated by mining massive human and mouse RNA-seq data sets, including 11 281 and 507 high-quality RNA-seq samples from 52 human non-disease tissues/cells and 14 healthy tissues/cells of C57BL/6 wild type mouse, respectively. User can visualize the expression and download lists of 2158 human HK transcripts from 2176 HK genes and 3024 mouse HK transcripts from 3277 mouse HK genes. HRT Atlas also offers the most stable and suitable tissue selective candidate reference transcripts for normalization of qPCR experiments. Specific primers and predicted modifiers of gene expression for some of these HK transcripts are also proposed. HRT Atlas has also been integrated with a regulatory elements resource from Epiregio server.

Many Different LINE-1 Retroelements Are Activated in Bladder Cancer

Human genomes contain about 100,000 LINE-1 (L1) retroelements, of which more than 100 are intact. L1s are normally tightly controlled by epigenetic mechanisms, which often fail in cancer. In bladder urothelial carcinoma (UC), particularly, L1s become DNA-hypomethylated, expressed and contribute to genomic instability and tumor growth. It is, however, unknown which individual L1s are activated. Following RNA-immunoprecipitation with a L1-specific antibody, third generation nanopore sequencing detected transcripts of 90 individual elements in the VM-Cub-1 UC line with high overall L1 expression. In total, 10 L1s accounted for >60% of the reads. Analysis of five specific L1s by RT-qPCR revealed generally increased expression in UC tissues and cell lines over normal controls, but variable expression among tumor cell lines from bladder, prostate and testicular cancer. Chromatin immunoprecipitation demonstrated active histone marks at L1 sequences with increased expression in VM-Cub-1, but not in a different UC cell line with low L1 expression. We conclude that many L1 elements are epigenetically activated in bladder cancer in a varied pattern. Our findings indicate that expression of individual L1s is highly heterogeneous between and among cancer types.

IncC conjugative plasmids and SXT/R391 elements repair double-strand breaks caused by CRISPR-Cas during conjugation

Bacteria have evolved defence mechanisms against bacteriophages. Restriction-modification systems provide innate immunity by degrading invading DNAs that lack proper methylation. CRISPR-Cas systems provide adaptive immunity by sampling the genome of past invaders and cutting the DNA of closely related DNA molecules. These barriers also restrict horizontal gene transfer mediated by conjugative plasmids. IncC conjugative plasmids are important contributors to the global dissemination of multidrug resistance among pathogenic bacteria infecting animals and humans. Here, we show that IncC conjugative plasmids are highly resilient to host defence systems during entry into a new host by conjugation. Using a TnSeq strategy, we uncover a conserved operon containing five genes (vcrx089-vcrx093) that confer a novel host defence evasion (hde) phenotype. We show that vcrx089-vcrx090 promote resistance against type I restriction-modification, whereas vcrx091-vcxr093 promote CRISPR-Cas evasion by repairing double-strand DNA breaks via recombination between short sequence repeats. vcrx091, vcrx092 and vcrx093 encode a single-strand binding protein, and a single-strand annealing recombinase and double-strand exonuclease related to Redβ and λExo of bacteriophage λ, respectively. Homologous genes of the integrative and conjugative element R391 also provide CRISPR-Cas evasion. Hence, the conserved hde operon considerably broadens the host range of large families of mobile elements spreading multidrug resistance.

Parallel Genetic Origin of Foot Feathering in Birds

Understanding the genetic basis of similar phenotypes shared between lineages is a long-lasting research interest. Even though animal evolution offers many examples of parallelism, for many phenotypes little is known about the underlying genes and mutations. We here use a combination of whole-genome sequencing, expression analyses, and comparative genomics to study the parallel genetic origin of ptilopody (Pti) in chicken. Ptilopody (or foot feathering) is a polygenic trait that can be observed in domesticated and wild avian species and is characterized by the partial or complete development of feathers on the ankle and feet. In domesticated birds, ptilopody is easily selected to fixation, though extensive variation in the type and level of feather development is often observed. By means of a genome-wide association analysis, we identified two genomic regions associated with ptilopody. At one of the loci, we identified a 17-kb deletion affecting PITX1 expression, a gene known to encode a transcription regulator of hindlimb identity and development. Similarly to pigeon, at the second loci, we observed ectopic expression of TBX5, a gene involved in forelimb identity and a key determinant of foot feather development. We also observed that the trait evolved only once as foot-feathered birds share the same haplotype upstream TBX5. Our findings indicate that in chicken and pigeon ptilopody is determined by the same set of genes that affect similar molecular pathways. Our study confirms that ptilopody has evolved through parallel evolution in chicken and pigeon.

Identification of a missense ARSA mutation in metachromatic leukodystrophy and its potential pathogenic mechanism

Background

Metachromatic leukodystrophy (MLD) is a rare inherited lysosomal disorder caused by mutations in ARSA. The biological processes of MLD disease caused by candidate pathogenic mutations in the ARSA gene remain unclear.

Methods

We used whole‐exome sequencing (WES) and Sanger sequencing to identify the pathogenic mutation in a Chinese family. Literature review and protein three‐dimensional structure prediction were performed to analyze the potential pathogenesis of the identified mutations. Overexpression cell models of wild‐type and mutated ARSA genes were constructed. The accumulated sulfatides and expression profiles in the cell models were detected, and a series of bioinformatics analyses were carried out to compare the biological changes caused by the candidate pathogenic mutations.

Results

We identified an ARSA c.925G>A homozygous mutation from a Chinese late‐infantile MLD patient, the first report of this mutation in East Asia. The literature and protein structure analysis indicated that three types of mutations at c.925G (c.925G>A, c.925G>T, c.925G>C) were pathogenic. The overexpression of wild‐type or mutated ARSA genes influenced the accumulation of sulfatides. The co‐expression modules in the mutated cell models were constructed by genes related to calcium signaling and vesicle transport.

Conclusion

Our results identified a pathogenic mutation, ARSA homozygosity c.925G>A, from a Chinese MLD family. The pathogenic mechanism of the ARSA mutation in MLD was identified, which may suggest new approaches to diagnosis and treatment.

Site-1 protease ablation in the osterix-lineage in mice results in bone marrow neutrophilia and hematopoietic stem cell alterations

Site-1 protease (S1P) ablation in the osterix-lineage in mice drastically reduces bone development and downregulates bone marrow-derived skeletal stem cells. Here we show that these mice also suffer from spina bifida occulta with a characteristic lack of bone fusion in the posterior neural arches. Molecular analysis of bone marrow-derived non-red blood cell cells, via single-cell RNA-Seq and protein mass spectrometry, demonstrate that these mice have a much-altered bone marrow with a significant increase in neutrophils and Ly6C-expressing leukocytes. The molecular composition of bone marrow neutrophils is also different as they express more and additional members of the stefin A (Stfa) family of proteins. In vitro, recombinant Stfa1 and Stfa2 proteins have the ability to drastically inhibit osteogenic differentiation of bone marrow stromal cells, with no effect on adipogenic differentiation. FACS analysis of hematopoietic stem cells show that despite a decrease in hematopoietic stem cells, S1P ablation results in an increased production of granulocyte-macrophage progenitors, the precursors to neutrophils. These observations indicate that S1P has a role in the lineage specification of hematopoietic stem cells and/or their progenitors for development of a normal hematopoietic niche. Our study designates a fundamental requirement of S1P for maintaining a balanced regenerative capacity of the bone marrow niche.

OncoBase: a platform for decoding regulatory somatic mutations in human cancers

Whole-exome and whole-genome sequencing have revealed millions of somatic mutations associated with different human cancers, and the vast majority of them are located outside of coding sequences, making it challenging to directly interpret their functional effects. With the rapid advances in high-throughput sequencing technologies, genome-scale long-range chromatin interactions were detected, and distal target genes of regulatory elements were determined using three-dimensional (3D) chromatin looping. Herein, we present OncoBase (http://www.oncobase.biols.ac.cn/), an integrated database for annotating 81 385 242 somatic mutations in 68 cancer types from more than 120 cancer projects by exploring their roles in distal interactions between target genes and regulatory elements. OncoBase integrates local chromatin signatures, 3D chromatin interactions in different cell types and reconstruction of enhancer-target networks using state-of-the-art algorithms. It employs informative visualization tools to display the integrated local and 3D chromatin signatures and effects of somatic mutations on regulatory elements. Enhancer-promoter interactions estimated from chromatin interactions are integrated into a network diffusion system that quantitatively prioritizes somatic mutations and target genes from a large pool. Thus, OncoBase is a useful resource for the functional annotation of regulatory noncoding regions and systematically benchmarking the regulatory effects of embedded noncoding somatic mutations in human carcinogenesis.

ENdb: a manually curated database of experimentally supported enhancers for human and mouse

Enhancers are a class of cis-regulatory elements that can increase gene transcription by forming loops in intergenic regions, introns and exons. Enhancers, as well as their associated target genes, and transcription factors (TFs) that bind to them, are highly associated with human disease and biological processes. Although some enhancer databases have been published, most only focus on enhancers identified by high-throughput experimental techniques. Therefore, it is highly desirable to construct a comprehensive resource of manually curated enhancers and their related information based on low-throughput experimental evidences. Here, we established a comprehensive manually-curated enhancer database for human and mouse, which provides a resource for experimentally supported enhancers, and to annotate the detailed information of enhancers. The current release of ENdb documents 737 experimentally validated enhancers and their related information, including 384 target genes, 263 TFs, 110 diseases and 153 functions in human and mouse. Moreover, the enhancer-related information was supported by experimental evidences, such as RNAi, in vitro knockdown, western blotting, qRT-PCR, luciferase reporter assay, chromatin conformation capture (3C) and chromosome conformation capture-on-chip (4C) assays. ENdb provides a user-friendly interface to query, browse and visualize the detailed information of enhancers. The database is available at http://www.licpathway.net/ENdb.

A genome-wide analysis of DNA methylation identifies a novel association signal for Lp(a) concentrations in the LPA promoter

Lipoprotein(a) [Lp(a)] is a major cardiovascular risk factor, which is largely genetically determined by one major gene locus, the LPA gene. Many aspects of the transcriptional regulation of LPA are poorly understood and the role of epigenetics has not been addressed yet. Therefore, we conducted an epigenome-wide analysis of DNA methylation on Lp(a) levels in two population-based studies (total n = 2208). We identified a CpG site in the LPA promoter which was significantly associated with Lp(a) concentrations. Surprisingly, the identified CpG site was found to overlap the SNP rs76735376. We genotyped this SNP de-novo in three studies (total n = 7512). The minor allele of rs76735376 (1.1% minor allele frequency) was associated with increased Lp(a) values (p = 1.01e-59) and explained 3.5% of the variation of Lp(a). Statistical mediation analysis showed that the effect on Lp(a) is rather originating from the base change itself and is not mediated by DNA methylation levels. This finding is supported by eQTL data from 208 liver tissue samples from the GTEx project, which shows a significant association of the rs76735376 minor allele with increased LPA expression. To evaluate, whether the association signal at rs76735376 may actually be derived from a stronger eQTL signal in LD with this SNP, eQTL association results of all correlated SNPs (r2≥0.1) were integrated with genetic association results. This analysis pinpointed to rs10455872 as the potential trigger of the effect of rs76735376. Furthermore, both SNPs coincide with short apo(a) isoforms. Adjusting for both, rs10455872 and the apo(a) isoforms diminished the effect size of rs76735376 to 5.38 mg/dL (p = 0.0463). This indicates that the effect of rs76735376 can be explained by both an independent effect of the SNP and a strong correlation with rs10455872 and apo(a) isoforms.

Targeted Sequencing of Tubo-ovarian and Peritoneal High-grade Serous Carcinoma With Wild-type p53 Immunostaining Pattern

BACKGROUND/AIM

We aimed to demonstrate the use of next-generation sequencing (NGS) to confirm the presence of tumor protein 53 (TP53) mutations in tubo-ovarian and peritoneal high-grade serous carcinoma (HGSC) with a wild-type p53 immunostaining pattern and investigate whether the TP53 mutational status is altered by chemotherapy.

MATERIALS AND METHODS

A commercial NGS panel comprising 171 genes was used to analyze the genetic profiles of 15 HGSC samples. Paired specimens obtained before and after chemotherapy were available for four patients.

RESULTS

All examined samples exhibited TP53 mutations. For all the patients who underwent neoadjuvant or postoperative adjuvant chemotherapy, TP53 mutations identified in samples obtained after chemotherapy were the same as those detected in pre-chemotherapeutic samples.

CONCLUSION

HGSCs exhibit TP53 mutations even though a subset of HGSCs displayed a wild-type p53 immunostaining pattern. Chemotherapy does not affect the TP53 mutational status in HGSC.

Clinicopathologic and Molecular Characteristics of Mesonephric Adenocarcinoma Arising From the Uterine Body

Mesonephric adenocarcinoma (MNAC) is a rare tumor of the female genital tract mainly occurring in the uterine cervix. To date, only a few cases of MNAC arising from of the uterine body (UB-MNAC) have been reported. The clinicopathologic and molecular characteristics of UB-MNAC remain unknown. In this study, we investigated the clinical, histopathologic, immunohistochemical, and genetic features of UB-MNAC. In total, 11 cases were included. Six patients developed metastatic disease, most commonly in lungs (5/6). Histopathologically, UB-MNAC was characterized by an admixture of tubular, glandular, papillary, retiform, glomeruloid, sex cord-like, and comedonecrosis-like architectural patterns. Three adverse pathologic characteristics, including advanced International Federation of Gynecology and Obstetrics stage, high mitotic activity, and presence of lymphovascular the invasion, were independent factors predicting the development of metastasis. All cases were positive for GATA-binding protein 3 and paired box 2 expression and showed wild-type p53, patchy p16, and preserved PTEN expression, as indicated by immunohistochemistry. Next-generation sequencing using 12 samples (11 primary tumors and 1 metastatic tumor) revealed 42 single nucleotide variations in 16 genes, mostly in KRAS (10/12) and ARID1A (9/12). Copy number variation was found in 16 genomic regions, and consisted of 57 gains and 10 losses, with 1q gain (11/12) being the most prevalent. In conclusion, UB-MNAC displays an aggressive biological behavior, with a tendency to metastasize to the lungs. Adverse pathologic characteristics reflect the aggressive nature of UB-MNAC. Distinct molecular features of UB-MNAC include frequent somatic mutations of KRAS and ARID1A and gain of 1q.

Bioinformatics analysis of the prognosis and biological significance of HMGB1, HMGB2, and HMGB3 in gastric cancer

High mobility group box (HMGB) consists primarily of HMGB1, HMGB2, and HMGB3 proteins. Although abnormal HMGB expression is associated with various tumors, the relationship with gastric cancer (GC) remains unclear. In this study, HMGB1, HMGB2, and HMGB3 expression was analyzed using the Oncomine and TCGA databases. Correlations between HMGB1, HMGB2, and HMGB3 and clinicopathological factors were analyzed. cBioPortal was used to analyze HMGB1, HMGB2, and HMGB3 genetic alterations and its gene regulation network in GC tissue. HMGB1, HMGB2, and HMGB3 expression was higher in tumor tissues than in normal tissues, especially in GC. High HMGB1, HMGB2, and HMGB3 expression may predict a poor prognosis among patients with GC (hazard ratios [HR] = 1.90; 95% confidence interval [CI]: [1.30-2.78]) and human digestive system neoplasm (HR = 1.85; 95% CI [1.64-2.10]). These findings suggest that HMGB1, HMGB2, and HMGB3 may be useful prognostic indicators for patients with GC.

Senataxin homologue Sen1 is required for efficient termination of RNA polymerase III transcription

R-loop disassembly by the human helicase Senataxin contributes to genome integrity and to proper transcription termination at a subset of RNA polymerase II genes. Whether Senataxin also contributes to transcription termination at other classes of genes has remained unclear. Here, we show that Sen1, one of two fission yeast homologues of Senataxin, promotes efficient termination of RNA polymerase III (RNAP3) transcription in vivo. In the absence of Sen1, RNAP3 accumulates downstream of RNAP3-transcribed genes and produces long exosome-sensitive 3'-extended transcripts. Importantly, neither of these defects was affected by the removal of R-loops. The finding that Sen1 acts as an ancillary factor for RNAP3 transcription termination in vivo challenges the pre-existing view that RNAP3 terminates transcription autonomously. We propose that Sen1 is a cofactor for transcription termination that has been co-opted by different RNA polymerases in the course of evolution.

Pergola-web: a web server for the visualization and analysis of longitudinal behavioral data using repurposed genomics tools and standards

We present a new web application to query and visualize time-series behavioral data: the Pergola web-server. This server provides a user-friendly interface for exploring longitudinal behavioral data taking advantage of the Pergola Python library. Using the server, users can process the data applying some basic operations, such as binning or grouping, while formatting the data into existing genomic formats. Thanks to this repurposing of genomics standards, the application automatically renders an interactive data visualization based on sophisticated genome visualization tools. Our tool allows behavioral scientists to share, display and navigate complex behavioral data comprising multiple individuals and multiple data types, in a scalable and flexible manner. A download option allows for further analysis using genomic tools. The server can be a great resource for the field in a time where behavioral science is entering a data-intensive cycle thanks to high-throughput behavioral phenotyping platforms. Pergola is publicly available at http://pergola.crg.eu/.

High levels of PIWI-interacting RNAs are present in the small RNA landscape of prostate epithelium from vitamin D clinical trial specimens

BACKGROUND

Vitamin D, a hormone that acts through the nuclear vitamin D receptor (VDR), upregulates antitumorigenic microRNA in prostate epithelium. This may contribute to the lower levels of aggressive prostate cancer (PCa) observed in patients with high serum vitamin D. The small noncoding RNA (ncRNA) landscape includes many other RNA species that remain uncharacterized in prostate epithelium and their potential regulation by vitamin D is unknown.

METHODS

Laser capture microdissection (LCM) followed by small-RNA sequencing was used to identify ncRNAs in the prostate epithelium of tissues from a vitamin D-supplementation trial. VDR chromatin immunoprecipitation-sequencing was performed to identify vitamin D genomic targets in primary prostate epithelial cells.

RESULTS

Isolation of epithelium by LCM increased sample homogeneity and captured more diversity in ncRNA species compared with publicly available small-RNA sequencing data from benign whole prostate. An abundance of PIWI-interacting RNAs (piRNAs) was detected in normal prostate epithelium. The obligate binding partners of piRNAs, PIWI-like (PIWIL) proteins, were also detected in prostate epithelium. High prostatic vitamin D levels were associated with increased expression of piRNAs. VDR binding sites were located near several ncRNA biogenesis genes and genes regulating translation and differentiation.

CONCLUSIONS

Benign prostate epithelium expresses both piRNA and PIWIL proteins, suggesting that these small ncRNA may serve an unknown function in the prostate. Vitamin D may increase the expression of prostatic piRNAs. VDR binding sites in primary prostate epithelial cells are consistent with its reported antitumorigenic functions and a role in ncRNA biogenesis.

Isolation and strategies of novel tetranucleotide microsatellites with polymorphisms from different chromosomes of the rhesus monkey (Macaca mulatta)

A total of 45 tetranucleotide chromosome-specific microsatellite markers with polymorphism were developed successfully based on three reference rhesus monkey genomes and on In-silico PCR prescreening. The polymorphic information content (PIC) values of 45 polymorphic microsatellite loci ranged from 0.487 to 0.879, with an average of 0.715, which were proven to be moderate to highly polymorphic. We detected 315 alleles on 45 microsatellite loci in 24 Rhesus monkeys. The number of alleles ranged from 3 to 15 and the mean number of alleles was 7 for each locus. Accordingly, the observed and expected heterozygosities obtained were between 0.417 and 1.0 and between 0.550 and 0.908, with an average value of 0.736 and 0.767, respectively. Genetic information demonstrated that 10 loci significantly deviated from Hardy-Weinberg equilibrium (P < 0.05). All 45 primers were not significant with regard to linkage disequilibrium (P > 0.001). Pearson correlation indicated that the PIC value exhibited a significant negative correlation with the loci number (r = - 0.741, P = 0.022), whereas the positive correlation with the number of the samples (r = 0.847, P = 0.070) was not significant. This may be attributed to the presence of random particularities within the loci. The T test of the sample groups indicated that the PIC difference was not significant when the number of samples was set at 10 and/or ≥ 15 (P = 0.7472 ~ 0.8564). These polymorphic and valuable microsatellite loci will facilitate further conservation genetics studies for rhesus monkeys and can be further applied to develop novel genetic markers for other species.

eDiVA-Classification and prioritization of pathogenic variants for clinical diagnostics

Mendelian diseases have shown to be an and efficient model for connecting genotypes to phenotypes and for elucidating the function of genes. Whole‐exome sequencing (WES) accelerated the study of rare Mendelian diseases in families, allowing for directly pinpointing rare causal mutations in genic regions without the need for linkage analysis. However, the low diagnostic rates of 20–30% reported for multiple WES disease studies point to the need for improved variant pathogenicity classification and causal variant prioritization methods. Here, we present the exome Disease Variant Analysis (eDiVA; http://ediva.crg.eu), an automated computational framework for identification of causal genetic variants (coding/splicing single‐nucleotide variants and small insertions and deletions) for rare diseases using WES of families or parent–child trios. eDiVA combines next‐generation sequencing data analysis, comprehensive functional annotation, and causal variant prioritization optimized for familial genetic disease studies. eDiVA features a machine learning‐based variant pathogenicity predictor combining various genomic and evolutionary signatures. Clinical information, such as disease phenotype or mode of inheritance, is incorporated to improve the precision of the prioritization algorithm. Benchmarking against state‐of‐the‐art competitors demonstrates that eDiVA consistently performed as a good or better than existing approach in terms of detection rate and precision. Moreover, we applied eDiVA to several familial disease cases to demonstrate its clinical applicability.

Dark Proteome Database: Studies on Dark Proteins

The dark proteome, as we define it, is the part of the proteome where 3D structure has not been observed either by homology modeling or by experimental characterization in the protein universe. From the 550.116 proteins available in Swiss-Prot (as of July 2016), 43.2% of the eukarya universe and 49.2% of the virus universe are part of the dark proteome. In bacteria and archaea, the percentage of the dark proteome presence is significantly less, at 12.6% and 13.3% respectively. In this work, we present a necessary step to complete the dark proteome picture by introducing the map of the dark proteome in the human and in other model organisms of special importance to mankind. The most significant result is that around 40% to 50% of the proteome of these organisms are still in the dark, where the higher percentages belong to higher eukaryotes (mouse and human organisms). Due to the amount of darkness present in the human organism being more than 50%, deeper studies were made, including the identification of ‘dark’ genes that are responsible for the production of so-called dark proteins, as well as the identification of the ‘dark’ tissues where dark proteins are over represented, namely, the heart, cervical mucosa, and natural killer cells. This is a step forward in the direction of gaining a deeper knowledge of the human dark proteome.

GenTree, an integrated resource for analyzing the evolution and function of primate-specific coding genes

The origination of new genes contributes to phenotypic evolution in humans. Two major challenges in the study of new genes are the inference of gene ages and annotation of their protein-coding potential. To tackle these challenges, we created GenTree, an integrated online database that compiles age inferences from three major methods together with functional genomic data for new genes. Genome-wide comparison of the age inference methods revealed that the synteny-based pipeline (SBP) is most suited for recently duplicated genes, whereas the protein-family–based methods are useful for ancient genes. For SBP-dated primate-specific protein-coding genes (PSGs), we performed manual evaluation based on published PSG lists and showed that SBP generated a conservative data set of PSGs by masking less reliable syntenic regions. After assessing the coding potential based on evolutionary constraint and peptide evidence from proteomic data, we curated a list of 254 PSGs with different levels of protein evidence. This list also includes 41 candidate misannotated pseudogenes that encode primate-specific short proteins. Coexpression analysis showed that PSGs are preferentially recruited into organs with rapidly evolving pathways such as spermatogenesis, immune response, mother–fetus interaction, and brain development. For brain development, primate-specific KRAB zinc-finger proteins (KZNFs) are specifically up-regulated in the mid-fetal stage, which may have contributed to the evolution of this critical stage. Altogether, hundreds of PSGs are either recruited to processes under strong selection pressure or to processes supporting an evolving novel organ.

New insights into the pathogenicity of non-synonymous variants through multi-level analysis

Precise classification of non-synonymous single nucleotide variants (SNVs) is a fundamental goal of clinical genetics. Next-generation sequencing technology is effective for establishing the basis of genetic diseases. However, identification of variants that are causal for genetic diseases remains a challenge. We analyzed human non-synonymous SNVs from a multilevel perspective to characterize pathogenicity. We showed that computational tools, though each having its own strength and weakness, tend to be overly dependent on the degree of conservation. For the mutations at non-degenerate sites, the amino acid sites of pathogenic substitutions show a distinct distribution in the classes of protein domains compared with the sites of benign substitutions. Overlooked disease susceptibility of genes explains in part the failures of computational tools. The more pathogenic sites observed, the more likely the gene is expressed in a high abundance or in a high tissue-specific manner, and have a high node degree of protein-protein interaction. The destroyed functions due to some false-negative mutations may arise because of a reprieve from the epigenetic repressed state which shouldn't happen in multiple biological conditions, instead of the defective protein. Our work adds more to our knowledge of non-synonymous SNVs' pathogenicity, thus will benefit the field of clinical genetics.

Design, optimization and validation of genes commonly used in expression studies on DMH/AOM rat colon carcinogenesis model

Colorectal cancer (CRC), also known as colon cancer, is the third most common form of cancer worldwide in men and the second in women and is characterized by several genetic alterations, among them the expression of several genes. 1,2-dimethylhydrazine (DMH) and its metabolite azoxymethane (AOM) are procarcinogens commonly used to induce colon cancer in rats (DMH/AOM rat model). This rat model has been used to study changes in mRNA expression in genes involved in this pathological condition. However, a lack of proper detailed PCR primer design in the literature limits the reproducibility of the published data. The present study aims to design, optimize and validate the qPCR, in accordance with the MIQE (Minimum Information for Publication of Quantitative Real-Time PCR Experiments) guidelines, for seventeen genes commonly used in the DMH/AOM rat model of CRC (Apc, Aurka, Bax, Bcl2, β-catenin, Ccnd1, Cdkn1a, Cox2, Gsk3beta, IL-33, iNOs, Nrf2, p53, RelA, Smad4, Tnfα and Vegfa) and two reference genes (Actb or β-actin and B2m). The specificity of all primer pairs was empirically validated on agarose gel, and furthermore, the melting curve inspection was checked as was their efficiency (%) ranging from 90 to 110 with a correlation coefficient of r² > 0.980. Finally, a pilot study was performed to compare the robustness of two candidate reference genes.

Differential genome-wide profiling of alternative polyadenylation sites in nasopharyngeal carcinoma by high-throughput sequencing

Background

Alternative polyadenylation (APA) is a widespread phenomenon in the posttranscriptional regulation of gene expression that generates mRNAs with alternative 3′-untranslated regions (3’UTRs). APA contributes to the pathogenesis of various diseases, including cancer. However, the potential role of APA in the development of nasopharyngeal carcinoma (NPC) remains largely unknown.

Methods

A strategy of sequencing APA sites (SAPAS) based on second-generation sequencing technology was carried out to explore the global patterns of APA sites and identify genes with tandem 3’UTRs in samples from 6 NPC and 6 normal nasopharyngeal epithelial tissue (NNET). Sequencing results were then validated using quantitative RT-PCR in a larger cohort of 16 NPC and 16 NNET samples.

Results

The sequencing data showed that the use of tandem APA sites was prevalent in NPC, and numerous genes with APA-switching events were discovered. In total, we identified 195 genes with significant differences in the tandem 3’UTR length between NPC and NNET; including 119 genes switching to distal poly (A) sites and 76 genes switching to proximal poly (A) sites. Several gene ontology (GO) terms were enriched in the list of genes with switched APA sites, including regulation of cell migration, macromolecule catabolic process, protein catabolic process, proteolysis, small conjugating protein ligase activity, and ubiquitin-protein ligase activity.

Conclusions

APA site-switching events are prevalent in NPC. APA-mediated regulation of gene expression may play an important role in the development of NPC, and more detailed studies targeting genes with APA-switching events may contribute to the development of novel future therapeutic strategies for NPC.

Electronic supplementary material

The online version of this article (10.1186/s12929-018-0477-6) contains supplementary material, which is available to authorized users.

Estimating the effects of transcription factors binding and histone modifications on gene expression levels in human cells

Transcription factors and histone modifications are vital for the regulation of gene expression. Hence, to estimate the effects of transcription factors binding and histone modifications on gene expression, we construct a statistical model for the genome-wide 15 transcription factors binding data, 10 histone modifications profiles and DNase-I hypersensitivity data in three mammalian. Remarkably, our results show POLR2A and H3K36me3 can highly and consistently predict gene expression in three cell lines. And H3K4me3, H3K27me3 and H3K9ac are more reliable predictors than other histone modifications in human embryonic stem cells. Moreover, genome-wide statistical redundancies exist within and between transcription factors and histone modifications, and these phenomena may be caused by the regulation mechanism. In further study, we find that even though transcription factors and histone modifications offer similar effects on expression levels of genome-wide genes, the effects of transcription factors and histone modifications on predictive abilities are different for genes in independent biological processes.

LSCC SNP variant regulates SOX2 modulation of VDAC3

Lung squamous cell carcinoma (LSCC) is a genomically complex malignancy with no effective treatments. Recent studies have found a large number of DNA alterations such as SOX2 amplification in LSCC patients. As a stem cell transcription factor, SOX2 is important for the maintenance of pluripotent cells and may play a role in cancer. To study the downstream mechanisms of SOX2, we employed expression quantitative trait loci (eQTLs) technology to investigate how the presence of SOX2 affects the expression of target genes. We discovered unique eQTLs, such as rs798827-VDAC3 (FDR p-value = 0.0034), that are only found in SOX2-active patients but not in SOX2-inactive patients. SNP rs798827 is within strong linkage disequilibrium (r2 = 1) to rs58163073, where rs58163073 [T] allele increases the binding affinity of SOX2 and allele [TA] decreases it. In our analysis, SOX2 silencing downregulates VDAC3 in two LSCC cell lines. Chromatin conformation capturing data indicates that this SNP is located within the same Topologically Associating Domain (TAD) of VDAC3, further suggesting SOX2's role in the regulation of VDAC3 through the binding of rs58163073. By first subgrouping patients based on SOX2 activity, we made more relevant eQTL discoveries and our analysis can be applied to other diseases.