Pseudogenes: pseudo-functional or key regulators in health and disease?

Factors Influencing Gene Family Size Variation Among Related Species in a Plant Family, Solanaceae

Gene duplication and loss contribute to gene content differences as well as phenotypic divergence across species. However, the extent to which gene content varies among closely related plant species and the factors responsible for such variation remain unclear. Here, using the Solanaceae family as a model and Pfam domain families as a proxy for gene families, we investigated variation in gene family sizes across species and the likely factors contributing to the variation. We found that genes in highly variable families have high turnover rates and tend to be involved in processes that have diverged between Solanaceae species, whereas genes in low-variability families tend to have housekeeping roles. In addition, genes in high- and low-variability gene families tend to be duplicated by tandem and whole genome duplication, respectively. This finding together with the observation that genes duplicated by different mechanisms experience different selection pressures suggest that duplication mechanism impacts gene family turnover. We explored using pseudogene number as a proxy for gene loss but discovered that a substantial number of pseudogenes are actually products of pseudogene duplication, contrary to the expectation that most plant pseudogenes are remnants of once-functional duplicates. Our findings reveal complex relationships between variation in gene family size, gene functions, duplication mechanism, and evolutionary rate. The patterns of lineage-specific gene family expansion within the Solanaceae provide the foundation for a better understanding of the genetic basis underlying phenotypic diversity in this economically important family.

An APOO Pseudogene on Chromosome 5q Is Associated With Low-Density Lipoprotein Cholesterol Levels

BACKGROUND

Elevated levels of low-density lipoprotein cholesterol (LDL-C) are a major risk factor for cardiovascular disease via its contribution to the development and progression of atherosclerotic lesions. Although the genetic basis of LDL-C has been studied extensively, currently known genetic variants account for only ≈20% of the variation in LDL-C levels.

METHODS

Through an array-based association analysis in 1102 Amish subjects, we identified a variant strongly associated with LDL-C levels. Using a combination of genetic analyses, zebrafish models, and in vitro experiments, we sought to identify the causal gene driving this association.

RESULTS

We identified a founder haplotype associated with a 15 mg/dL increase in LDL-C on chromosome 5. After recombination mapping, the associated region contained 8 candidate genes. Using a zebrafish model to evaluate the relevance of these genes to cholesterol metabolism, we found that expression of the transcribed pseudogene, APOOP1, increased LDL-C and vascular plaque formation.

CONCLUSIONS

Based on these data, we propose that APOOP1 regulates levels of LDL-C in humans, thus identifying a novel mechanism of lipid homeostasis.

Single-nucleotide and copy-number variance related to severity of hypospadias

BACKGROUND

The genetic association of hypospadias-risk studies has been conducted in Caucasians, Chinese-Han populations and few in Indian populations. However, no comprehensive approach has been followed to assess genetic involvement in the severity of the disorder.

METHODS

The study evaluated to establish the correlation between genotyped single nucleotide and copy number variants (SNPs/CNVs) and severity of hypospadias by an association in a total 30 SNPs in genes related to sex hormone-biosynthesis and metabolism; embryonic-development and phospholipase-D-signalling pathways on 138 surgery-confirmed hypospadias-cases from North India (84 penile and 28 cases of penoscrotal-hypospadias as compared with 31 cases of glanular + coronal), and analyzed and identified CNVs in four familial cases (18 members) and three paired-sporadic cases (6 members) using array-based comparative-genomic-hybridization and validated in 32 hypospadias samples by TaqMan assay.

RESULTS

Based on odds ratio at 95% CI, Z Statistic and Significance Levels, STS gene-rs17268974 was associated with Penile-Hypospadias and 9-SNPs [seven-SNPs (rs5934740; rs5934842; rs5934913; rs6639811; rs3923341; rs17268974; rs5934937)] of STS gene; rs7562326-SRD5A2 and rs1877031-STARD3 were associated with penoscrotal-hypospadias. On aggregate analysis with p < 0.001, we identified homozygous-loss of Ch7:q34 (PRSS3P2, PRSS2). On validation in previously CNV-characterized and new (32 hypospadias cases), we identified PRSS3P2-loss in most of the grade 3 and 4 hypospadias. Hence, Grade 1 and 2 (coronal and granular) show no-PRSS3P2-loss and no-association with SNPs in STS; SRD5A2; STARD3-gene but Grade 3 and 4 (Penile and Penoscrotal) show PRSS3P2-loss accompanied with the association of SNPs in STS; SRD5A2; STARD3.

CONCLUSIONS

Hence, homozygous-loss of PRSS3P2 accompanied with the association of STS; SRD5A2; STARD3 may link to the severity of the disease.

Single nucleotide polymorphisms in the MYLKP1 pseudogene are associated with increased colon cancer risk in African Americans

INTRODUCTION

Pseudogenes are paralogues of functional genes historically viewed as defunct due to either the lack of regulatory elements or the presence of frameshift mutations. Recent evidence, however, suggests that pseudogenes may regulate gene expression, although the functional role of pseudogenes remains largely unknown. We previously reported that MYLKP1, the pseudogene of MYLK that encodes myosin light chain kinase (MLCK), is highly expressed in lung and colon cancer cell lines and tissues but not in normal lung or colon. The MYLKP1 promoter is minimally active in normal bronchial epithelial cells but highly active in lung adenocarcinoma cells. In this study, we further validate MYLKP1 as an oncogene via elucidation of the functional role of MYLKP1 genetic variants in colon cancer risk.

METHODS

Proliferation and migration assays were performed in MYLKP1-transfected colon and lung cancer cell lines (H441, A549) and commercially-available normal lung and colon cells. Fourteen MYLKP1 SNPs (MAFs >0.01) residing within the 4 kb MYLKP1 promoter region, the core 1.4 kb of MYLKP1 gene, and a 4 kb enhancer region were selected and genotyped in a colorectal cancer cohort. MYLKP1 SNP influences on activity of MYLKP1 promoter (2kb) was assessed by dual luciferase reporter assay.

RESULTS

Cancer cell lines, H441 and A549, exhibited increased MYLKP1 expression, increased MYLKP1 luciferase promoter activity, increased proliferation and migration. Genotyping studies identified two MYLKP1 SNPs (rs12490683; rs12497343) that significantly increase risk of colon cancer in African Americans compared to African American controls. Rs12490683 and rs12497343 further increase MYLKP1 promoter activity compared to the wild type MYLKP1 promoter.

CONCLUSION

MYLKP1 is a cancer-promoting pseudogene whose genetic variants differentially enhance cancer risk in African American populations.

Marfan Syndrome Variability: Investigation of the Roles of Sarcolipin and Calcium as Potential Transregulator of FBN1 Expression

Marfan syndrome (MFS) is an autosomal dominant connective tissue disorder that displays a great clinical variability. Previous work in our laboratory showed that fibrillin-1 (FBN1) messenger RNA (mRNA) expression is a surrogate endpoint for MFS severity. Therefore, an expression quantitative trait loci (eQTL) analysis was performed to identify trans-acting regulators of FBN1 expression, and a significant signal reached genome-wide significant threshold on chromosome 11. This signal delineated a region comprising one expressed gene, SLN (encoding sarcolipin), and a single pseudogene, SNX7-ps1 (CTD-2651C21.3). We first investigated the region and then looked for association between the genes in the region and FBN1 expression. For the first time, we showed that the SLN gene is weakly expressed in skin fibroblasts. There is no direct correlation between SLN and FBN1 gene expression. We showed that calcium influx modulates FBN1 gene expression. Finally, SLN gene expression is highly correlated to that of the neighboring SNX7-ps1. We were able to confirm the impact of calcium influx on FBN1 gene expression but we could not conclude regarding the role of sarcolipin and/or the eQTL locus in this regulation.

Striatal transcriptome of a mouse model of ADHD reveals a pattern of synaptic remodeling

Despite the prevalence and high heritability of Attention-Deficit/Hyperactivity Disorder (ADHD), genetic etiology remains elusive. Clinical evidence points in part to reduced function of the striatum, but which specific genes are differentially expressed and how they sculpt striatal physiology to predispose ADHD are not well understood. As an exploratory tool, a polygenic mouse model of ADHD was recently developed through selective breeding for high home cage activity. Relative to the Control line, the High-Active line displays hyperactivity and motor impulsivity which are ameliorated with amphetamine. This study compared gene expression in the striatum between Control and High-Active mice to develop a coherent hypothesis for how genes might affect striatal physiology and predispose ADHD-like symptoms. To this end, striatal transcriptomes of High-Active and Control mice were analyzed after mice were treated with saline or amphetamines. The pseudogene Gm6180 for n-cofilin (Cfl1) displayed 20-fold higher expression in High-Active mice corresponding with reduced Cfl1 expression suggesting synaptic actin dysregulation. Latrophilin 3 (Lphn3), which is associated with ADHD in human populations and is involved in synapse structure, and its ligand fibronectin leucine rich transmembrane protein 3 (Flrt3), were downregulated in High-Active mice. Multiple genes were altered in High-Active mice in a manner predicted to downregulate the canonical Wnt pathway. A smaller and different set of genes including glyoxalase (Glo1) were differentially regulated in High-Active as compared to Control in response to amphetamine. Together, results suggest genes involved in excitatory synapse regulation and maintenance are downregulated in ADHD-like mice. Consistent with the molecular prediction, stereological analysis of the striatum from a separate set of mice processed for imunohistochemical detection of synaptophysin revealed approximately a 46% reduction in synaptophysin immunoreactivity in High-Active relative to Control. Results provide a new set of molecular targets related to synapse maintenance for the next generation of ADHD medicines.

Pseudogenes in gastric cancer pathogenesis: a review article

Cancer burden rises globally at an alarming pace. According to GLOBOCAN 2012, gastric cancer (GC) is regarded as the fifth most common malignancy in the world. Being twice as high in men as in women, GC is the third leading cause of cancer mortality in both sexes globally. Being labeled as 'junk DNA', pseudogenes were considered as nonfunctional 'trash', which contribute nothing to survival of the organism; therefore, a number of strategies have been developed to circumvent their accidental detection. Recent progresses have confirmed that pseudogenes can have broad and multifaceted spectrum of activities in human cancers in general and GC in particular. Furthermore, the mentioned functions are parental gene-dependent and/or -independent. Therefore, pseudogenes can be regarded as the emerging class of elaborate modulators of gene expression involved in pathogenesis of human cancers including gastric adenocarcinoma.

Landscape of the genome and host cell response of Mycobacterium shigaense reveals pathogenic features

A systems approach was used to explore the genome and transcriptome of Mycobacterium shigaense, a new opportunistic pathogen isolated from a patient with a skin infection, and the host response transcriptome was assessed using a macrophage infection model. The M. shigaense genome comprises 5,207,883 bp, with 67.2% G+C content and 5098 predicted coding genes. Evolutionarily, the bacterium belongs to a cluster in the phylogenetic tree along with three target opportunistic pathogenic strains, namely, M. avium, M. triplex and M. simiae. Potential virulence genes are indeed expressed by M. shigaense under culture conditions. Phenotypically, M. shigaense had similar infection and replication capacities in a macrophage model as the opportunistic species compared to M. tuberculosis. M. shigaense activated NF-κB, TNF, cytokines and chemokines in the host innate immune-related signaling pathways and elicited an early response shared with pathogenic bacilli except M. tuberculosis. M. shigaense upregulated specific host response genes such as TLR7, CCL4 and CXCL5. We performed an integrated and comparative analysis of M. shigaense. Multigroup comparison indicated certain differences with typical pathogenic bacilli in terms of gene features and the macrophage response.

Knockdown of SUMO1P3 represses tumor growth and invasion and enhances radiosensitivity in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is one of the most common types of malignant tumors with high recurrence and metastasis rates. Radiotherapy represents a major therapeutic option for HCC patients. However, the efficacy of radiotherapy has been limited due to the development of intrinsic radioresistance of the tumor cells. Small ubiquitin-like modifier 1 pseudogene 3 (SUMO1P3), one member of SUMO pseudogene family, is a novel identified lncRNA that was originally identified to be upregulated in gastric cancer. However, the detailed roles of SUMO1P3 in HCC development remain to be elucidated. Here, the expression of SUMO1P3 in HCC tissues and cells was examined by qRT-PCR. Cell proliferation, colony formation ability, invasion ability, apoptosis, and radiosensitivity were detected by MTT assay, colony formation assay, cell invasion assay, flow cytometry analysis, and survival fraction assay, respectively. We found that SUMO1P3 was significantly upregulated in HCC tissues and cells. Besides, SUMO1P3 was highly expressed in HCC patients with higher TNM stage. Furthermore, SUMO1P3 knockdown markedly suppressed cell proliferation, colony formation ability, and cell invasiveness, promoted apoptosis, and enhanced radiosensitivity of HCC cells. We concluded that the knockdown of SUMO1P3 repressed tumor growth, invasion, promoted apoptosis, and enhanced radiosensitivity in HCC, providing evidence that SUMO1P3 might be a potential novel biomarker and a therapeutic target for HCC.

The Complex Interaction of Mitochondrial Genetics and Mitochondrial Pathways in Psychiatric Disease

While accounting for only 2% of the body's weight, the brain utilizes up to 20% of the body's total energy. Not surprisingly, metabolic dysfunction and energy supply-and-demand mismatch have been implicated in a variety of neurological and psychiatric disorders. Mitochondria are responsible for providing the brain with most of its energetic demands, and the brain uses glucose as its exclusive energy source. Exploring the role of mitochondrial dysfunction in the etiology of psychiatric disease is a promising avenue to investigate further. Genetic analysis of mitochondrial activity is a cornerstone in understanding disease pathogenesis related to metabolic dysfunction. In concert with neuroimaging and pathological study, genetics provides an important bridge between biochemical findings and clinical correlates in psychiatric disease. Mitochondrial genetics has several unique aspects to its analysis, and corresponding special considerations. Here, we review the components of mitochondrial genetic analysis - nuclear DNA, mitochon-drial DNA, mitochondrial pathways, pseudogenes, nuclear-mitochondrial mismatch, and microRNAs - that could contribute to an observable clinical phenotype. Throughout, we highlight psychiatric diseases that can arise due to dysfunction in these processes, with a focus on schizophrenia and bipolar disorder.

The pseudogene-derived long non-coding RNA SFTA1P suppresses cell proliferation, migration, and invasion in gastric cancer

Pseudogenes were once regarded as transcriptionally inactive and without specific molecular function. However, current evidence shows that pseudogene-derived long non-coding RNAs (lncRNAs) may be crucial regulators of human cancer development, including gastric cancer (GC). In the present study, we report that a pseudogene-derived lncRNA named surfactant associated 1, pseudogene (SFTA1P), which is 693-nt long, was significantly down-regulated in GC tissues compared with that in the adjacent normal tissues. In addition, decreased SFTA1P expression was strongly correlated with advanced tumor lymph node metastasis (TNM) stage, larger tumor size, lymphatic metastasis, and poor prognosis of patients with GC. Moreover, gain-of-function experiments revealed that the overexpression of SFTA1P inhibits cell proliferation, migration, and invasion, thus verifying the tumor inhibitory role of SFTA1P in GC. Furthermore, we investigated the potential action mechanism of SFTA1P. Our results showed that down-regulation of SFTA1P may be associated with decreased TP53 expression. In summary, our work suggests that the pseudogene-derived lncRNA SFTA1P functions as a tumor suppressor in GC and thus may act as a potential diagnostic and therapeutic target of GC.

PTEN/PTENP1: 'Regulating the regulator of RTK-dependent PI3K/Akt signalling', new targets for cancer therapy

Regulation of the PI-3 kinase (PI3K)/Akt signalling pathway is essential for maintaining the integrity of fundamental cellular processes, cell growth, survival, death and metabolism, and dysregulation of this pathway is implicated in the development and progression of cancers. Receptor tyrosine kinases (RTKs) are major upstream regulators of PI3K/Akt signalling. The phosphatase and tensin homologue (PTEN), a well characterised tumour suppressor, is a prime antagonist of PI3K and therefore a negative regulator of this pathway. Loss or inactivation of PTEN, which occurs in many tumour types, leads to overactivation of RTK/PI3K/Akt signalling driving tumourigenesis. Cellular PTEN levels are tightly regulated by a number of transcriptional, post-transcriptional and post-translational regulatory mechanisms. Of particular interest, transcription of the PTEN pseudogene, PTENP1, produces sense and antisense transcripts that exhibit post-transcriptional and transcriptional modulation of PTEN expression respectively. These additional levels of regulatory complexity governing PTEN expression add to the overall intricacies of the regulation of RTK/PI-3 K/Akt signalling. This review will discuss the regulation of oncogenic PI3K signalling by PTEN (the regulator) with a focus on the modulatory effects of the sense and antisense transcripts of PTENP1 on PTEN expression, and will further explore the potential for new therapeutic opportunities in cancer treatment.

Downregulated pseudogene CTNNAP1 promote tumor growth in human cancer by downregulating its cognate gene CTNNA1 expression

Accumulating evidence indicates that deregulation of cancer-associated pseudogene is involved in the pathogenesis of cancer. In the study, we demonstrated that pseudogene CTNNAP1, for the CTNNA1 gene, was dysregulated in colorectal cancer and the degree of dysregulation was remarkably associated with tumor node metastasis (TNM) stage (P<0.05). The mechanistic experiments revealed that pseudogene CTNNAP1 played a pivotal role in the regulation of its cognate gene CTNNA1 by competition for microRNA-141. Moreover, gain-of-function approaches showed that overexpression of CTNNAP1 or CTNNA1 significantly inhibited cell proliferation and tumor growth in vitro and in vivo by inducing G0/G1 cell cycle arrest. Our findings add a new regulatory circuit via competing endogenous RNA (ceRNA) cross-talk between pseudogene CTNNAP1 and its cognate gene CTNNA1, and provide new insights into potential diagnostic biomarker for monitoring human colorectal cancer.

Pseudogene BMI1P1 expression as a novel predictor for acute myeloid leukemia development and prognosis

The BMI1P1 levels of 144 de novo AML patients and 36 healthy donors were detected by real-time quantitative PCR (RQ-PCR). BMI1P1 was significantly down-regulated in AML compared with control (P < 0.001). A receiver operating characteristic (ROC) curve revealed that BMI1P1 expression could differentiate patients with AML from control subjects (AUC = 0.895, 95% CI: 0.835–0.954, P < 0.001). The percentage of blasts in bone marrow (BM) was significantly lower in BMI1P1 high-expressed group versus low-expressed group (P = 0.008). BMI1P1 high-expressed cases had significantly higher complete remission (CR) than BMI1P1 low-expressed cases (P = 0.023). Furthermore, Kaplan–Meier demonstrated that both whole AML cohort and non-M3-AML patients with low BMI1P1 expression showed shorter leukemia free survival (LFS, P = 0.002 and P = 0.01, respectively) and overall survival (OS, P < 0.001 and P = 0.011, respectively) than those with high BMI1P1 expression. Multivariate analysis also showed that BMI1P1 over-expression was an independent favorable prognostic factor for OS in both whole and non-M3 cohort of AML patients (HR = 0.462, 95% CI = 0.243–0.879, P = 0.019 and HR = 0.483, 95% CI = 0.254–0.919, P = 0.027). To further investigate the significance of BMI1P1 expression in the follow-up of AML patients, we monitored the BMI1P1 level in 26 de novo AML patients and found that the BMI1P1 level increased significantly from the initial diagnosis to post-CR (P < 0.001). These results indicated that BMI1P1 might contribute to the diagnosis of AML and the assessment of therapeutic effect.

The pseudogenes of barley

Pseudogenes have a reputation of being 'evolutionary relics' or 'junk DNA'. While they are well characterized in mammals, studies in more complex plant genomes have so far been hampered by the absence of reference genome sequences. Barley is one of the economically most important cereals and has a genome size of 5.1 Gb. With the first high-quality genome reference assembly available for a Triticeae crop, we conducted a whole-genome assessment of pseudogenes on the barley genome. We identified, characterized and classified 89 440 gene fragments and pseudogenes scattered along the chromosomes, with occasional hotspots and higher densities at the chromosome ends. Full-length pseudogenes (11 015) have preferentially retained their exon-intron structure. Retrotransposition of processed mRNAs only plays a marginal role in their creation. However, the distribution of retroposed pseudogenes reflects the Rabl configuration of barley chromosomes and thus hints at founding mechanisms. While parent genes related to the defense-response were found to be under-represented in cultivated barley, we detected several defense-related pseudogenes in wild barley accessions. The percentage of transcriptionally active pseudogenes is 7.2%, and these may potentially adopt new regulatory roles.The barley genome is rich in pseudogenes and small gene fragments mainly located towards chromosome tips or as tandemly repeated units. Our results indicate non-random duplication and pseudogenization preferences and improve our understanding of the dynamics of gene birth and death in large plant genomes and the mechanisms that lead to evolutionary innovations.

Long non-coding RNA MT1DP shunts the cellular defense to cytotoxicity through crosstalk with MT1H and RhoC in cadmium stress

Metallothioneins (MTs) are known to protect cells against oxidative stress, especially providing protection against cadmium (Cd) toxicity in hepatocytes. There are various gene variants and pseudogenes for MTs; however, there is little understanding on the functions of those non-coding MT members that are known to be expressed as long non-coding RNAs (lncRNAs) nowadays. Different from most protein-coding MT members, MT1DP was here found that remarkably induced to provoke cytotoxicity in hepatocytes in response to Cd treatment. MT1DP exerted such a pro-apoptotic function in Cd-treated hepatocytes through interacting with two partners: RhoC and MT1H. On one hand, MT1DP interacted with RhoC protein to increase the latter’s stability by preventing lysosome-dependent protein degradation. Therefore, upon Cd stress, MT1DP/RhoC complex was quickly reinforced to activate RhoC-CCN1/2-AKT signaling and potentiate Ca²⁺ influx, leading to enhanced Cd uptake and elevated Cd toxicity. On the other hand, MT1H, a protein-coding member of the MT family with little known function, was found to quickly respond to Cd exposure along with MT1DP. Mechanistically, MT1H and MT1DP were uncovered to mutually protect each other through a reciprocal ceRNA mechanism, building up a positive feedback loop to enforce MT1DP-conducted signaling upon Cd exposure. Moreover, MT1DP was found to contribute much more to the activation of RhoC-CCN1/2-AKT signaling than MT1H. Considered together, we here unveiled a mystery whether a pseudogene within the MT family, MT1DP, has actual biological functions in regulating Cd-induced cellular defense. Our findings unearthed an important role of pseudogene MT1DP in calibrating the cellular machinery to switch the cellular defense to cytotoxicity through crosslinking an interplay between its two partners, namely MT1H and RhoC, under cadmium stress.

Resistance to Plum Pox Virus (PPV) in apricot (Prunus armeniaca L.) is associated with down-regulation of two MATHd genes

BACKGROUND

Plum pox virus (PPV), causing Sharka disease, is one of the main limiting factors for Prunus production worldwide. In apricot (Prunus armeniaca L.) the major PPV resistance locus (PPVres), comprising ~ 196 kb, has been mapped to the upper part of linkage group 1. Within the PPVres, 68 genomic variants linked in coupling to PPV resistance were identified within 23 predicted transcripts according to peach genome annotation. Taking into account the predicted functions inferred from sequence homology, some members of a cluster of meprin and TRAF-C homology domain (MATHd)-containing genes were pointed as PPV resistance candidate genes.

RESULTS

Here, we have characterized the global apricot transcriptome response to PPV-D infection identifying six PPVres locus genes (ParP-1 to ParP-6) differentially expressed in resistant/susceptible cultivars. Two of them (ParP-3 and ParP-4), that encode MATHd proteins, appear clearly down-regulated in resistant cultivars, as confirmed by qRT-PCR. Concurrently, variant calling was performed using whole-genome sequencing data of 24 apricot cultivars (10 PPV-resistant and 14 PPV-susceptible) and 2 wild relatives (PPV-susceptible). ParP-3 and ParP-4, named as Prunus armeniaca PPVres MATHd-containing genes (ParPMC), are the only 2 genes having allelic variants linked in coupling to PPV resistance. ParPMC1 has 1 nsSNP, while ParPMC2 has 15 variants, including a 5-bp deletion within the second exon that produces a frameshift mutation. ParPMC1 and ParPMC2 are adjacent and highly homologous (87.5% identity) suggesting they are paralogs originated from a tandem duplication. Cultivars carrying the ParPMC2 resistant (mutated) allele show lack of expression in both ParPMC2 and especially ParPMC1.

CONCLUSIONS

Accordingly, we hypothesize that ParPMC2 is a pseudogene that mediates down-regulation of its functional paralog ParPMC1 by silencing. As a whole, results strongly support ParPMC1 and/or ParPMC2 as host susceptibility genes required for PPV infection which silencing may confer PPV resistance trait. This finding may facilitate resistance breeding by marker-assisted selection and pave the way for gene edition approaches in Prunus.

Over-expression of oncigenic pesudogene DUXAP10 promotes cell proliferation and invasion by regulating LATS1 and β-catenin in gastric cancer

BACKGROUND

Recently, the pesudogenes have emerged as critical regulators in human cancers tumorigenesis and progression, and been identified as a key revelation in post-genomic biology. However, the expression pattern, biological function and mechanisms responsible for these molecules in human gastric cancer (GC) are not fully understood.

METHODS

In this study, we globally assessed the transcriptomic differences of pesudogenes in gastric cancer using publicly available microarray data. DUXAP10 expression levels in GC tissues and cells was detected using quantitative real-time PCR (qPCR). DUXAP10 siRNAs and over-expression vector were transfected into GC cells to down-regulate or up-regulate DUXAP10 expression. Loss- and gain-of function assays were performed to investigate the role of DUXAP10 in GC cells cell proliferation, and invasion. RIP, RNA pulldown, and ChIP assays were used to determine the mechanism of DUXAP10's regulation of underlying targets.

RESULTS

The pesudogene DUXAP10 is the only pseudogene that significantly over-expressed in all four GEO datasets, and frequently over-expressed in many other cancers including Liver Hepatocellular carcinoma, Bladder cancer, and Esophageal Cancer. High DUXAP10 expression is associated with GC patients poor prognosis, and knockdown of DUXAP10 significantly inhibits cells proliferation, migration and invasion in GC. Mechanistic investigation shows that DUXAP10 can interact with PRC2 and LSD1 to repress LATS1 expression at transcriptional level, and bind with HuR to maintain the stability of β-catenin mRNA and increase its protein levels at post-transcriptional level.

CONCLUSIONS

Overall, our findings illuminate how increased DUXAP10 confers an oncogenic function in GC development and progression that may serve as a candidate prognostic biomarker and target for clinical management of GC.

Genome-Wide Identification and Comparative Analysis of the 3-Hydroxy-3-methylglutaryl Coenzyme A Reductase (HMGR) Gene Family in Gossypium

Terpenes are the largest and most diverse class of secondary metabolites in plants and play a very important role in plant adaptation to environment. 3-Hydroxy-3-methylglutaryl coenzyme A reductase (HMGR) is a rate-limiting enzyme in the process of terpene biosynthesis in the cytosol. Previous study found the HMGR genes underwent gene expansion in Gossypium raimondii, but the characteristics and evolution of the HMGR gene family in Gossypium genus are unclear. In this study, genome-wide identification and comparative study of HMGR gene family were carried out in three Gossypium species with genome sequences, i.e., G. raimondii, Gossypium arboreum, and Gossypium hirsutum. In total, nine, nine and 18 HMGR genes were identified in G. raimondii, G. arboreum, and G. hirsutum, respectively. The results indicated that the HMGR genes underwent gene expansion and a unique gene cluster containing four HMGR genes was found in all the three Gossypium species. The phylogenetic analysis suggested that the expansion of HMGR genes had occurred in their common ancestor. There was a pseudogene that had a 10-bp deletion resulting in a frameshift mutation and could not be translated into functional proteins in G. arboreum and the A-subgenome of G. hirsutum. The expression profiles of the two pseudogenes showed that they had tissue-specific expression. Additionally, the expression pattern of the pseudogene in the A-subgenome of G. hirsutum was similar to its paralogous gene in the D-subgenome of G. hirsutum. Our results provide useful information for understanding cytosolic terpene biosynthesis in Gossypium species.

Widespread modulation of gene expression by copy number variation in skeletal muscle

Copy number variation (CNV) is a frequently observed deviation from the diploid state due to duplication or deletion of genomic regions. Although intensively analyzed for association with diseases and production traits, the specific mechanisms and extent by which such variations affect the phenotype are incompletely understood. We present an integrative study on CNV and genome-wide gene expression in Brazilian Bos indicus cattle. We analyzed CNVs inferred from SNP-chip data for effects on gene expression measured with RNA-seq in skeletal muscle samples of 183 steers. Local effects, where expression changes coincided with CNVs in the respective genes, were restricted to immune genes. Distal effects were attributable to several high-impact CNVs that modulated remote expression in an orchestrated and intertwined fashion. These CNVs were located in the vicinity of major skeletal muscle pathway regulators and associated genes were enriched for proteolysis, autophagy, and muscle structure development. From association analysis between CNVs and several meat quality and production traits, we found CNV-associated expression effects to also manifest at the phenotype level. Based on genome sequences of the population founders, we further demonstrate that CNVs with impact on expression and phenotype are passed on from one generation to another.

Sex differences in the late first trimester human placenta transcriptome

Background

Development of the placenta during the late first trimester is critical to ensure normal growth and development of the fetus. Developmental differences in this window such as sex-specific variation are implicated in later placental disease states, yet gene expression at this time is poorly understood.

Methods

RNA-sequencing was performed to characterize the transcriptome of 39 first trimester human placentas using chorionic villi following genetic testing (17 females, 22 males). Gene enrichment analysis was performed to find enriched canonical pathways and gene ontologies in the first trimester. DESeq2 was used to find sexually dimorphic gene expression. Patient demographics were analyzed for sex differences in fetal weight at time of chorionic villus sampling and birth.

Results

RNA-sequencing analyses detected 14,250 expressed genes, with chromosome 19 contributing the greatest proportion (973/2852, 34.1% of chromosome 19 genes) and Y chromosome contributing the least (16/568, 2.8%). Several placenta-enriched genes as well as histone-coding genes were identified to be unique to the first trimester and common to both sexes. Further, we identified 58 genes with significantly different expression between males and females: 25 X-linked, 15 Y-linked, and 18 autosomal genes. Genes that escape X inactivation were highly represented (59.1%) among X-linked genes upregulated in females. Many genes differentially expressed by sex consisted of X/Y gene pairs, suggesting that dosage compensation plays a role in sex differences. These X/Y pairs had roles in parallel, ancient canonical pathways important for eukaryotic cell growth and survival: chromatin modification, transcription, splicing, and translation.

Conclusions

This study is the first characterization of the late first trimester placenta transcriptome, highlighting similarities and differences among the sexes in ongoing human pregnancies resulting in live births. Sexual dimorphism may contribute to pregnancy outcomes, including fetal growth and birth weight, which was seen in our cohort, with males significantly heavier than females at birth. This transcriptome provides a basis for development of early diagnostic tests of placental function that can indicate overall pregnancy heath, fetal-maternal health, and long-term adult health.

Electronic supplementary material

The online version of this article (10.1186/s13293-018-0165-y) contains supplementary material, which is available to authorized users.

Improved ontology for eukaryotic single-exon coding sequences in biological databases

Efficient extraction of knowledge from biological data requires the development of structured vocabularies to unambiguously define biological terms. This paper proposes descriptions and definitions to disambiguate the term 'single-exon gene'. Eukaryotic Single-Exon Genes (SEGs) have been defined as genes that do not have introns in their protein coding sequences. They have been studied not only to determine their origin and evolution but also because their expression has been linked to several types of human cancer and neurological/developmental disorders and many exhibit tissue-specific transcription. Unfortunately, the term 'SEGs' is rife with ambiguity, leading to biological misinterpretations. In the classic definition, no distinction is made between SEGs that harbor introns in their untranslated regions (UTRs) versus those without. This distinction is important to make because the presence of introns in UTRs affects transcriptional regulation and post-transcriptional processing of the mRNA. In addition, recent whole-transcriptome shotgun sequencing has led to the discovery of many examples of single-exon mRNAs that arise from alternative splicing of multi-exon genes, these single-exon isoforms are being confused with SEGs despite their clearly different origin. The increasing expansion of RNA-seq datasets makes it imperative to distinguish the different SEG types before annotation errors become indelibly propagated in biological databases. This paper develops a structured vocabulary for their disambiguation, allowing a major reassessment of their evolutionary trajectories, regulation, RNA processing and transport, and provides the opportunity to improve the detection of gene associations with disorders including cancers, neurological and developmental diseases.

Involvement of miRNAs and Pseudogenes in Cancer

Our understanding of cancer pathways has been changed by the determination of noncoding transcripts in the human genome in recent years. miRNAs and pseudogenes are key players of the noncoding transcripts from the genome, and alteration of their expression levels provides clues for significant biomarkers in pathogenesis of diseases. Especially, miRNAs and pseudogenes have both oncogenic and tumor-suppressive roles in each step of cancer tumorigenesis. In this current study, association between oncogenes and miRNAs-pseudogenes was reviewed and determined in human cancer by the CellMiner web-tool.

Network analysis of pseudogene-gene relationships: from pseudogene evolution to their functional potentials

Pseudogenes are fossil relatives of genes. Pseudogenes have long been thought of as "junk DNAs", since they do not code proteins in normal tissues. Although most of the human pseudogenes do not have noticeable functions, ∼20% of them exhibit transcriptional activity. There has been evidence showing that some pseudogenes adopted functions as lncRNAs and work as regulators of gene expression. Furthermore, pseudogenes can even be "reactivated" in some conditions, such as cancer initiation. Some pseudogenes are transcribed in specific cancer types, and some are even translated into proteins as observed in several cancer cell lines. All the above have shown that pseudogenes could have functional roles or potentials in the genome. Evaluating the relationships between pseudogenes and their gene counterparts could help us reveal the evolutionary path of pseudogenes and associate pseudogenes with functional potentials. It also provides an insight into the regulatory networks involving pseudogenes with transcriptional and even translational activities.In this study, we develop a novel approach integrating graph analysis, sequence alignment and functional analysis to evaluate pseudogene-gene relationships, and apply it to human gene homologs and pseudogenes. We generated a comprehensive set of 445 pseudogene-gene (PGG) families from the original 3,281 gene families (13.56%). Of these 438 (98.4% PGG, 13.3% total) were non-trivial (containing more than one pseudogene). Each PGG family contains multiple genes and pseudogenes with high sequence similarity. For each family, we generate a sequence alignment network and phylogenetic trees recapitulating the evolutionary paths. We find evidence supporting the evolution history of olfactory family (both genes and pseudogenes) in human, which also supports the validity of our analysis method. Next, we evaluate these networks in respect to the gene ontology from which we identify functions enriched in these pseudogene-gene families and infer functional impact of pseudogenes involved in the networks. This demonstrates the application of our PGG network database in the study of pseudogene function in disease context.

Human Accelerated Regions and Other Human-Specific Sequence Variations in the Context of Evolution and Their Relevance for Brain Development

The review discusses, in a format of a timeline, the studies of different types of genetic variants, present in Homo sapiens, but absent in all other primate, mammalian, or vertebrate species, tested so far. The main characteristic of these variants is that they are found in regions of high evolutionary conservation. These sequence variations include single nucleotide substitutions (called human accelerated regions), deletions, and segmental duplications. The rationale for finding such variations in the human genome is that they could be responsible for traits, specific to our species, of which the human brain is the most remarkable. As became obvious, the vast majority of human-specific single nucleotide substitutions are found in noncoding, likely regulatory regions. A number of genes, associated with these human-specific alleles, often through novel enhancer activity, were in fact shown to be implicated in human-specific development of certain brain areas, including the prefrontal cortex. Human-specific deletions may remove regulatory sequences, such as enhancers. Segmental duplications, because of their large size, create new coding sequences, like new functional paralogs. Further functional study of these variants will shed light on evolution of our species, as well as on the etiology of neurodevelopmental disorders.

Inferring the evolution of the major histocompatibility complex of wild pigs and peccaries using hybridisation DNA capture-based sequencing

The major histocompatibility complex (MHC) is a key genomic model region for understanding the evolution of gene families and the co-evolution between host and pathogen. To date, MHC studies have mostly focused on species from major vertebrate lineages. The evolution of MHC classical (Ia) and non-classical (Ib) genes in pigs has attracted interest because of their antigen presentation roles as part of the adaptive immune system. The pig family Suidae comprises over 18 extant species (mostly wild), but only the domestic pig has been extensively sequenced and annotated. To address this, we used a DNA-capture approach, with probes designed from the domestic pig genome, to generate MHC data for 11 wild species of pigs and their closest living family, Tayassuidae. The approach showed good efficiency for wild pigs (~80% reads mapped, ~87× coverage), compared to tayassuids (~12% reads mapped, ~4× coverage). We retrieved 145 MHC loci across both families. Phylogenetic analyses show that the class Ia and Ib genes underwent multiple duplications and diversifications before suids and tayassuids diverged from their common ancestor. The histocompatibility genes mostly form orthologous groups and there is genetic differentiation for most of these genes between Eurasian and sub-Saharan African wild pigs. Tests of selection showed that the peptide-binding region of class Ib genes was under positive selection. These findings contribute to better understanding of the evolutionary history of the MHC, specifically, the class I genes, and provide useful data for investigating the immune response of wild populations against pathogens.

Competing endogenous RNA networks in human cancer: hypothesis, validation, and perspectives

Non-coding RNAs represent a majority of the human transcriptome. However, less is known about the functions and regulatory mechanisms of most non-coding species. Moreover, little is known about the potential non-coding functions of coding RNAs. The competing endogenous RNAs (ceRNAs) hypothesis is proposed recently. This hypothesis describes potential communication networks among all transcript RNA species mediated by miRNAs and miRNA-recognizing elements (MREs) within RNA transcripts. Here we review the evolution of the ceRNA hypothesis, summarize the validation experiments and discusses the significance and perspectives of this hypothesis in human cancer.

Transcription of a B chromosome CAP-G pseudogene does not influence normal Condensin Complex genes in a grasshopper

Parasitic B chromosomes invade and persist in natural populations through several mechanisms for transmission advantage (drive). They may contain gene-derived sequences which, in some cases, are actively transcribed. A further interesting question is whether B-derived transcripts become functional products. In the grasshopper Eyprepocnemis plorans, one of the gene-derived sequences located on the B chromosome shows homology with the gene coding for the CAP-G subunit of condensin I. We show here, by means of fluorescent in situ hybridization coupled with tyramide signal amplification (FISH-TSA), that this gene is located in the distal region of the B24 chromosome variant. The DNA sequence located in the B chromosome is a pseudogenic version of the CAP-G gene (B-CAP-G). In two Spanish populations, we found active transcription of B-CAP-G, but it did not influence the expression of CAP-D2 and CAP-D3 genes coding for corresponding condensin I and II subunits, respectively. Our results indicate that the transcriptional regulation of the B-CAP-G pseudogene is uncoupled from the standard regulation of the genes that constitute the condensin complex, and suggest that some of the B chromosome known effects may be related with its gene content and transcriptional activity, thus opening new exciting avenues for research.

HMGA1-pseudogenes and cancer

Pseudogenes are DNA sequences with high homology to the corresponding functional gene, but, because of the accumulation of various mutations, they have lost their initial functions to code for proteins. Consequently, pseudogenes have been considered until few years ago dysfunctional relatives of the corresponding ancestral genes, and then useless in the course of genome evolution. However, several studies have recently established that pseudogenes are owners of key biological functions. Indeed, some pseudogenes control the expression of functional genes by competitively binding to the miRNAs, some of them generate small interference RNAs to negatively modulate the expression of functional genes, and some of them even encode functional mutated proteins. Here, we concentrate our attention on the pseudogenes of the HMGA1 gene, that codes for the HMGA1a and HMGA1b proteins having a critical role in development and cancer progression. In this review, we analyze the family of HMGA1 pseudogenes through three aspects: classification, characterization, and their possible function and involvement in cancer.

A Matter of Taste: Lineage-Specific Loss of Function of Taste Receptor Genes in Vertebrates

Vertebrates can perceive at least five different taste qualities, each of which is thought to have a specific role in the evolution of different species. The avoidance of potentially poisonous foods, which are generally bitter or sour tasting, and the search for more nutritious ones, those with high-fat and high-sugar content, are two of the most well-known examples. The study of taste genes encoding receptors that recognize ligands triggering taste sensations has helped to reconstruct several evolutionary adaptations to dietary changes. In addition, an increasing number of studies have focused on pseudogenes, genomic DNA sequences that have traditionally been considered defunct relatives of functional genes mostly because of the presence of deleterious mutations interrupting their open reading frames. The study of taste receptor pseudogenes has helped to shed light on how the evolutionary history of taste in vertebrates has been the result of a succession of gene gain and loss processes. This dynamic role in evolution has been explained by the "less-is-more" hypothesis, suggesting gene loss as a mechanism of evolutionary change in response to a dietary shift. This mini-review aims at depicting the major lineage-specific loss of function of taste receptor genes in vertebrates, stressing their evolutionary importance and recapitulating signatures of natural selection and their correlations with food habits.

Identification of potential cancer-related pseudogenes in lung adenocarcinoma based on ceRNA hypothesis

Pseudogenes are initially regarded as non-functional genomic fossils resulted from inactivating gene mutations during evolution. Far from being silent, pseudogenes are proved to regulate the expression of protein-coding genes through function as microRNA sponge in vivo. The aim of our study was to propose an integrative systems biology approach to identify disease pseudogenes base on competitive endogenous RNA (ceRNA) hypothesis. Here, we applied our method to lung adenocarcinoma (LUAD) RNASeq data from TCGA and identified 33 candidate pseudogenes. We described the characteristics of the candidate pseudogenes and performed functional enrichment. Through analyzing neighboring genes we found these pseudogenes were surrounded by tumor genes and may involve in tumor pathway. Furthermore, the DNA methylation analysis indicated that 21 pseudogenes co-methylated with their competitive mRNAs. In the co-methylated network, we discovered 6 differentially expressed pseudogenes, which we termed potential LUAD-associated pseudogenes. We further revealed that the 3 ceRNA triples (miR-21-5p-NKAPP1-PRDM11, miR-29c-3p-MSTO2P-EZH2 and miR-29c-3p-RPLP0P2-EZH2), whose high risk groups were associated with the poor prognosis of LUAD, may be considered as potential prognostic signatures. Moreover, by integrating target information of microRNA we also provided a new perspective for the discovery of potential small molecule drugs. This work may facilitate cancer research and serve as the basis for future efforts to understand the role of pseudogenes, develop novel biomarkers and improve knowledge of tumor biology.

Computational Prediction of Position Effects of Apparently Balanced Human Chromosomal Rearrangements

Interpretation of variants of uncertain significance, especially chromosomal rearrangements in non-coding regions of the human genome, remains one of the biggest challenges in modern molecular diagnosis. To improve our understanding and interpretation of such variants, we used high-resolution three-dimensional chromosomal structural data and transcriptional regulatory information to predict position effects and their association with pathogenic phenotypes in 17 subjects with apparently balanced chromosomal abnormalities. We found that the rearrangements predict disruption of long-range chromatin interactions between several enhancers and genes whose annotated clinical features are strongly associated with the subjects' phenotypes. We confirm gene-expression changes for a couple of candidate genes to exemplify the utility of our analysis of position effect. These results highlight the important interplay between chromosomal structure and disease and demonstrate the need to utilize chromatin conformational data for the prediction of position effects in the clinical interpretation of non-coding chromosomal rearrangements.

Neuronal-expressed microRNA-targeted pseudogenes compete with coding genes in the human brain

MicroRNAs orchestrate brain functioning via interaction with microRNA recognition elements (MRE) on target transcripts. However, the global impact of potential competition on the microRNA pool between coding and non-coding brain transcripts that share MREs with them remains unexplored. Here we report that non-coding pseudogene transcripts carrying MREs (PSG^+MRE) often show duplicated origin, evolutionary conservation and higher expression in human temporal lobe neurons than comparable duplicated MRE-deficient pseudogenes (PSG^-MRE). PSG^+MRE participate in neuronal RNA-induced silencing complexes (RISC), indicating functional involvement. Furthermore, downregulation cell culture experiments validated bidirectional co-regulation of PSG^+MRE with MRE-sharing coding transcripts, frequently not their mother genes, and with targeted microRNAs; also, PSG^+MRE single-nucleotide polymorphisms associated with schizophrenia, bipolar disorder and autism, suggesting interaction with mental diseases. Our findings indicate functional roles of duplicated PSG^+MRE in brain development and cognition, supporting physiological impact of the reciprocal co-regulation of PSG^+MRE with MRE-sharing coding transcripts in human brain neurons.

A comprehensive genomic pan-cancer classification using The Cancer Genome Atlas gene expression data

BACKGROUND

The Cancer Genome Atlas (TCGA) has generated comprehensive molecular profiles. We aim to identify a set of genes whose expression patterns can distinguish diverse tumor types. Those features may serve as biomarkers for tumor diagnosis and drug development.

METHODS

Using RNA-seq expression data, we undertook a pan-cancer classification of 9,096 TCGA tumor samples representing 31 tumor types. We randomly assigned 75% of samples into training and 25% into testing, proportionally allocating samples from each tumor type.

RESULTS

We could correctly classify more than 90% of the test set samples. Accuracies were high for all but three of the 31 tumor types, in particular, for READ (rectum adenocarcinoma) which was largely indistinguishable from COAD (colon adenocarcinoma). We also carried out pan-cancer classification, separately for males and females, on 23 sex non-specific tumor types (those unrelated to reproductive organs). Results from these gender-specific analyses largely recapitulated results when gender was ignored. Remarkably, more than 80% of the 100 most discriminative genes selected from each gender separately overlapped. Genes that were differentially expressed between genders included BNC1, FAT2, FOXA1, and HOXA11. FOXA1 has been shown to play a role for sexual dimorphism in liver cancer. The differentially discriminative genes we identified might be important for the gender differences in tumor incidence and survival.

CONCLUSIONS

We were able to identify many sets of 20 genes that could correctly classify more than 90% of the samples from 31 different tumor types using TCGA RNA-seq data. This accuracy is remarkable given the number of the tumor types and the total number of samples involved. We achieved similar results when we analyzed 23 non-sex-specific tumor types separately for males and females. We regard the frequency with which a gene appeared in those sets as measuring its importance for tumor classification. One third of the 50 most frequently appearing genes were pseudogenes; the degree of enrichment may be indicative of their importance in tumor classification. Lastly, we identified a few genes that might play a role in sexual dimorphism in certain cancers.

Down-regulation of pseudogene Vimentin 2p is associated with poor outcome in de novo acute myeloid leukemia

OBJECTIVES

This study was intended to investigate the expression status of Vimentin 2p (VIM 2p), a pseudogene of Vimentin, and further analyze its clinical significance in AML patients.

METHODS

Real-time quantitative PCR (RQ-PCR) was employed to explore the expression status of VIM 2p in 128 patients with de novo AML and 36 healthy controls.

RESULTS

The expression level of VIM 2p was significantly decreased compared with healthy controls (P< 0.001). The patients with low VIM 2p expression were identified in 93 of 128 (73%) of AML patients. No significant differences could be observed in sex, age, blood parameters, FAB/WHO subtypes, karyotype risks and ten gene mutations (FLT3-ITD, NPM1, C-KIT, IDH1/IDH2, DNMT3 A, C/EBPA, N/K-RAS and U2AF1) between VIM 2p low-expressed and high-expressed patients (P> 0.05). Patients with low VIM 2p expression had significantly shorter overall survival (OS) than those with high VIM 2p expression in whole AML cases (median 7 vs. 13 months, respectively, P= 0.032), besides cytogenetically normal AML (CN-AML) and non-M3 AML cohort (P= 0.042 and 0.045, respectively).

CONCLUSIONS

These findings indicate that VIM 2p down-regulation is a common event in AML and may be associated with poor clinical outcome.

Xp11.22 deletions encompassing CENPVL1, CENPVL2, MAGED1 and GSPT2 as a cause of syndromic X-linked intellectual disability

By searching a clinical database of over 60,000 individuals referred for array-based CNV analyses and online resources, we identified four males from three families with intellectual disability, developmental delay, hypotonia, joint hypermobility and relative macrocephaly who carried small, overlapping deletions of Xp11.22. The maximum region of overlap between their deletions spanned ~430 kb and included two pseudogenes, CENPVL1 and CENPVL2, whose functions are not known, and two protein coding genes-the G1 to S phase transition 2 gene (GSPT2) and the MAGE family member D1 gene (MAGED1). Deletions of this ~430 kb region have not been previously implicated in human disease. Duplications of GSPT2 have been documented in individuals with intellectual disability, but the phenotypic consequences of a loss of GSPT2 function have not been elucidated in humans or mouse models. Changes in MAGED1 have not been associated with intellectual disability in humans, but loss of MAGED1 function is associated with neurocognitive and neurobehavioral phenotypes in mice. In all cases, the Xp11.22 deletion was inherited from an unaffected mother. Studies performed on DNA from one of these mothers did not show evidence of skewed X-inactivation. These results suggest that deletions of an ~430 kb region on chromosome Xp11.22 that encompass CENPVL1, CENPVL2, GSPT2 and MAGED1 cause a distinct X-linked syndrome characterized by intellectual disability, developmental delay, hypotonia, joint hypermobility and relative macrocephaly. Loss of GSPT2 and/or MAGED1 function may contribute to the intellectual disability and developmental delay seen in males with these deletions.

The pseudogene-derived long noncoding RNA SFTA1P is down-regulated and suppresses cell migration and invasion in lung adenocarcinoma

Pseudogenes were once considered to be genomic fossils without biological function. Interestingly, recent evidence showed that a lot of pseudogenes are transcribed in human cancers, and their alterations contribute to multiple cancer development and progression. It is apparent that many pseudogenes transcribe noncoding RNAs and contribute to the role noncoding genome plays in human cancers. On this basis, some pseudogene transcripts are currently ranked among the classes of long noncoding RNAs. In this study, we identified a new pseudogene-derived long noncoding RNA termed SFTA1P by analyzing the microarray data of non-small cell lung cancer from Gene Expression Omnibus datasets. We found that SFTA1P expression was significantly decreased in non-small cell lung cancer tissues compared with normal tissues in non-small cell lung cancer microarray data. Moreover, decreased SFTA1P expression is only correlated with lung adenocarcinoma patients' poor survival time but not with lung squamous cell carcinoma patients' survival. In addition, gain-of-function studies including growth curves, migration, invasion assays, and in vivo studies were performed to verify the tumor suppressor role of SFTA1P in non-small cell lung cancer. Finally, the potential underlying pathways involved in SFTA1P were investigated by analyzing the SFTA1P-correlated genes in The Cancer Genome Atlas lung adenocarcinoma and normal tissues RNA sequencing data. Taken together, these findings demonstrate that pseudogene-derived long noncoding RNA SFTA1P exerts the tumor suppressor functions in human lung adenocarcinoma. Our investigation reveals the novel roles of pseudogene in lung adenocarcinoma, which may serve as a new target for lung adenocarcinoma diagnosis and therapy.

The origin of B chromosomes in yellow-necked mice (Apodemus flavicollis)-Break rules but keep playing the game

B chromosomes (Bs) are known for more than hundred years but their origin, structure and pattern of evolution are not well understood. In the past few years new methodological approaches, involving isolation of Bs followed by whole DNA amplification, DNA probe generation, and fluorescent in situ hybridization (FISH) or the B chromosome DNA sequencing, has allowed detailed analysis of their origin and molecular structure in different species. In this study we explored the origin of Bs in the yellow-necked wood mouse, Apodemus flavicollis, using generation of microdissected DNA probes followed by FISH on metaphase chromosomes. Bs of A. flavicollis were successfully isolated and DNA was used as the template for B-specific probes for the first time. We revealed homology of DNA derived from the analyzed B chromosomes to the pericentromeric region (PR) of sex chromosomes and subtelomeric region of two pairs of small autosomes, but lower homology to the rest of the Y chromosome. Moreover, all analysed Bs had the same structure regardless of their number per individual or the great geographic distance between examined populations from the Balkan Peninsula (Serbia) and Eastern Europe (south region of Russia and central Belarus). Therefore, it was suggested that B chromosomes in A. flavicollis have a unique common origin from the PR of sex chromosomes, and/or similar evolutionary pattern.

Mechanism of Deletion Removing All Dystrophin Exons in a Canine Model for DMD Implicates Concerted Evolution of X Chromosome Pseudogenes

Duchenne muscular dystrophy (DMD) is a lethal, X-linked, muscle-wasting disorder caused by mutations in the large, 2.4-Mb dystrophin gene. The majority of DMD-causing mutations are sporadic, multi-exon, frameshifting deletions, with the potential for variable immunological tolerance to the dystrophin protein from patient to patient. While systemic gene therapy holds promise in the treatment of DMD, immune responses to vectors and transgenes must first be rigorously evaluated in informative preclinical models to ensure patient safety. A widely used canine model for DMD, golden retriever muscular dystrophy, expresses detectable amounts of near full-length dystrophin due to alternative splicing around an intronic point mutation, thereby confounding the interpretation of immune responses to dystrophin-derived gene therapies. Here we characterize a naturally occurring deletion in a dystrophin-null canine, the German shorthaired pointer. The deletion spans 5.6 Mb of the X chromosome and encompasses all coding exons of the DMD and TMEM47 genes. The sequences surrounding the deletion breakpoints are virtually identical, suggesting that the deletion occurred through a homologous recombination event. Interestingly, the deletion breakpoints are within loci that are syntenically conserved among mammals, yet the high homology among this subset of ferritin-like loci is unique to the canine genome, suggesting lineage-specific concerted evolution of these atypical sequence elements.

The role of long non-coding RNAs in nasopharyngeal carcinoma: As systemic review

Recent development of cutting edge research found that long noncoding RNAs (lncRNAs) plays important roles in carcinogenesis and progression. In Southeast Asia and North Africa, nasopharyngeal carcinoma (NPC) is the most common aggressive squamous cell carcinoma. Nasopharyngeal carcinoma is most frequently occurring in males. However, nasopharyngeal carcinoma is caused by a combination of several factors as viral, environmental factors, and heredity. Till now, the potential pathway or mechanism of NPC is not well known. In our present review, we strongly emphasized on long noncoding RNAs (lncRNAs) and its significant role in nasopharyngeal carcinoma. It has been showed that lncRNAs regulate the development and progression of different types of cancers, including NPC. In addition, it has been found that chromatin organization, transcriptional and post-transcriptional events are regulated by lncRNAs. Our present review summarizes the roles of lncRNAs in nasopharyngeal carcinoma and provides an overview of the feasibility of lncRNAs as diagnosis, prognosis and potential treatment for NPC patients.

Next-Generation Sequencing of Protein-Coding and Long Non-protein-Coding RNAs in Two Types of Exosomes Derived from Human Whole Saliva

Exosomes are small extracellular vesicles containing microRNAs and mRNAs that are produced by various types of cells. We previously used ultrafiltration and size-exclusion chromatography to isolate two types of human salivary exosomes (exosomes I, II) that are different in size and proteomes. We showed that salivary exosomes contain large repertoires of small RNAs. However, precise information regarding long RNAs in salivary exosomes has not been fully determined. In this study, we investigated the compositions of protein-coding RNAs (pcRNAs) and long non-protein-coding RNAs (lncRNAs) of exosome I, exosome II and whole saliva (WS) by next-generation sequencing technology. Although 11% of all RNAs were commonly detected among the three samples, the compositions of reads mapping to known RNAs were similar. The most abundant pcRNA is ribosomal RNA protein, and pcRNAs of some salivary proteins such as S100 calcium-binding protein A8 (protein S100-A8) were present in salivary exosomes. Interestingly, lncRNAs of pseudogenes (presumably, processed pseudogenes) were abundant in exosome I, exosome II and WS. Translationally controlled tumor protein gene, which plays an important role in cell proliferation, cell death and immune responses, was highly expressed as pcRNA and pseudogenes in salivary exosomes. Our results show that salivary exosomes contain various types of RNAs such as pseudogenes and small RNAs, and may mediate intercellular communication by transferring these RNAs to target cells as gene expression regulators.

Araport11: a complete reannotation of the Arabidopsis thaliana reference genome

The flowering plant Arabidopsis thaliana is a dicot model organism for research in many aspects of plant biology. A comprehensive annotation of its genome paves the way for understanding the functions and activities of all types of transcripts, including mRNA, the various classes of non-coding RNA, and small RNA. The TAIR10 annotation update had a profound impact on Arabidopsis research but was released more than 5 years ago. Maintaining the accuracy of the annotation continues to be a prerequisite for future progress. Using an integrative annotation pipeline, we assembled tissue-specific RNA-Seq libraries from 113 datasets and constructed 48 359 transcript models of protein-coding genes in eleven tissues. In addition, we annotated various classes of non-coding RNA including microRNA, long intergenic RNA, small nucleolar RNA, natural antisense transcript, small nuclear RNA, and small RNA using published datasets and in-house analytic results. Altogether, we identified 635 novel protein-coding genes, 508 novel transcribed regions, 5178 non-coding RNAs, and 35 846 small RNA loci that were formerly unannotated. Analysis of the splicing events and RNA-Seq based expression profiles revealed the landscapes of gene structures, untranslated regions, and splicing activities to be more intricate than previously appreciated. Furthermore, we present 692 uniformly expressed housekeeping genes, 43% of whose human orthologs are also housekeeping genes. This updated Arabidopsis genome annotation with a substantially increased resolution of gene models will not only further our understanding of the biological processes of this plant model but also of other species.

Pseudogene PTENP1 Functions as a Competing Endogenous RNA (ceRNA) to Regulate PTEN Expression by Sponging miR-499-5p

Increasing evidence has shown that pseudogenes can widely regulate gene expression. However, little is known about the specific role of PTENP1 and miR-499-5p in insulin resistance. The relative transcription level of PTENP1 was examined in db/db mice and high fat diet (HFD)-fed mice by real-time PCR. To explore the effect of PTENP1 on insulin resistance, adenovirus overexpressing or inhibiting vectors were injected through the tail vein. Bioinformatics predictions and a luciferase reporter assay were used to explore the interaction between PTENP1 and miR-499-5p. The relative transcription level of PTENP1 was largely enhanced in db/db mice and HFD-fed mice. Furthermore, the overexpression of PTENP1 resulted in impaired Akt/GSK activation as well as glycogen synthesis, while PTENP1 inhibition led to the improved activation of Akt/GSK and enhanced glycogen contents. More importantly, PTENP1 could directly bind miR-499-5p, thereby becoming a sink for miR-499-5p. PTENP1 overexpression results in the impairment of the insulin-signaling pathway and may function as a competing endogenous RNA for miR-499-5p, thereby contributing to insulin resistance.

Transcriptome and long noncoding RNA sequencing of three extracellular vesicle subtypes released from the human colon cancer LIM1863 cell line

Previously we reported that LIM1863 colorectal cancer (CRC) cells secrete three distinct extracellular vesicle subtypes – two subpopulations of exosomes (apical EpCAM-Exos and basolateral A33-Exos) and shed microvesicles (sMVs) – with distinct protein and miRNA signatures. Here, we extend our omics approach to understand the fundamental role of LIM1863-derived EVs by performing a comprehensive analysis of their mRNAs and long non-coding RNAs (lncRNAs) using RNA-Seq. We show that 2,389 mRNAs, 317 pseudogene transcripts, 1,028 lncRNAs and 206 short non-coding RNAs selectively distributed to (i.e., are enriched in) LIM1863 EVs, relative to the parent cell. An Ensembl/UniProtKB analysis revealed 1,937 mRNAs encode canonical proteins, 348 isoforms (including splice-variant proteins), and 119 ‘missing proteins’ (i.e., annotated in Ensembl but not UniProtKB). Further dissection of our protein/RNA data revealed that 6/151 observed RNA binding proteins have the potential to interact with ~75% of EV-enriched RNAs. Intriguingly, the co-existence of U1 and U2 ribonucleoproteins and their cognate snRNAs in LIM1863 EVs suggests a possible association of CRC EVs with recipient cell splicing events. Our data reveal several potential lncRNA CRC biomarkers and novel splicing/fusion genes that, collectively, will advance our understanding of EV biology in CRC and accelerate the development of EV-based diagnostics and therapeutics.

The correlation between expression profiles measured in single cells and in traditional bulk samples

Reverse transcription quantitative PCR (RT-qPCR) is already an established tool for mRNA detection and quantification. Since recently, this technique has been successfully employed for gene expression analyses, and also in individual cells (single cell RT-qPCR). Although the advantages of single cell measurements have been proven several times, a study correlating the expression measured on single cells, and in bulk samples consisting of a large number of cells, has been missing. Here, we collected a large data set to explore the relation between gene expression measured in single cells and in bulk samples, reflected by qPCR Cq values. We measured the expression of 95 genes in 12 bulk samples, each containing thousands of astrocytes, and also in 693 individual astrocytes. Combining the data, we described the relation between Cq values measured in bulk samples with either the percentage of the single cells that express the given genes, or the average expression of the genes across the single cells. We show that data obtained with single cell RT-qPCR are fully consistent with measurements in bulk samples. Our results further provide a base for quality control in single cell expression profiling, and bring new insights into the biological process of cellular expression.

Genomic organization and evolution of olfactory receptors and trace amine-associated receptors in channel catfish, Ictalurus punctatus

BACKGROUND

Channel catfish (Ictalurus punctatus) live in turbid waters with limited visibility to chase prey within a certain distance. This can be compensated through detecting specific water-soluble substances by the olfactory receptors (ORs) and trace amine associated receptors (TAARs) expressed on the olfactory epithelium.

METHODS

We identified the OR and TAAR repertoires in channel catfish, and characterized the genomic organizations of these two gene families by data mining available genomic resources.

RESULTS

A total of 47 putative OR genes and 36 putative TAAR genes were identified in the channel catfish genome, including 27 functional OR genes and 28 functional TAAR genes. Phylogenetic and orthogroup analyses were conducted to illustrate the evolutionary dynamics of the vertebrate ORs and TAARs. Collinear analysis revealed the presence of two conserved orthologous blocks that contain OR genes between the catfish genome and zebrafish genome. The complete loss of a conserved motif in fish OR family H may contribute to the divergence of family H from other families. The dN/dS analysis indicated that the highest degree of selection pressure was imposed on TAAR subfamily 14 among all fish ORs and TAARs.

CONCLUSIONS

The present study provides understanding of the evolutionary dynamics of the two gene families (OR and TAAR) associated with olfaction in channel catfish.

GENERAL SIGNIFICANCE

This is the first systematic study of ORs and TAARs in catfish, which could provide valuable genomic resources for further investigation of olfactory mechanisms in teleost fish.

A Genomic Analysis of Factors Driving lincRNA Diversification: Lessons from Plants

Transcriptomic analyses from across eukaryotes indicate that most of the genome is transcribed at some point in the developmental trajectory of an organism. One class of these transcripts is termed long intergenic noncoding RNAs (lincRNAs). Recently, attention has focused on understanding the evolutionary dynamics of lincRNAs, particularly their conservation within genomes. Here, we take a comparative genomic and phylogenetic approach to uncover factors influencing lincRNA emergence and persistence in the plant family Brassicaceae, to which Arabidopsis thaliana belongs. We searched 10 genomes across the family for evidence of > 5000 lincRNA loci from A. thaliana. From loci conserved in the genomes of multiple species, we built alignments and inferred phylogeny. We then used gene tree/species tree reconciliation to examine the duplication history and timing of emergence of these loci. Emergence of lincRNA loci appears to be linked to local duplication events, but, surprisingly, not whole genome duplication events (WGD), or transposable elements. Interestingly, WGD events are associated with the loss of loci for species having undergone relatively recent polyploidy. Lastly, we identify 1180 loci of the 6480 previously annotated A. thaliana lincRNAs (18%) with elevated levels of conservation. These conserved lincRNAs show higher expression, and are enriched for stress-responsiveness and cis-regulatory motifs known as conserved noncoding sequences (CNSs). These data highlight potential functional pathways and suggest that CNSs may regulate neighboring genes at both the genomic and transcriptomic level. In sum, we provide insight into processes that may influence lincRNA diversification by providing an evolutionary context for previously annotated lincRNAs.