Loss to gain: pseudogenes in microorganisms, focusing on eubacteria, and their biological significance.

Pseudogenes are defined as "non-functional" copies of corresponding parent genes. The cognition of pseudogenes continues to be refreshed through accumulating and updating research findings. Previous studies have predominantly focused on mammals, but pseudogenes have received relatively less attention in the field of microbiology. Given the increasing recognition on the importance of pseudogenes, in this review, we focus on several aspects of microorganism pseudogenes, including their classification and characteristics, their generation and fate, their identification, their abundance and distribution, their impact on virulence, their ability to recombine with functional genes, the extent to which some pseudogenes are transcribed and translated, and the relationship between pseudogenes and viruses. By summarizing and organizing the latest research progress, this review will provide a comprehensive perspective and improved understanding on pseudogenes in microorganisms. KEY POINTS: • Concept, classification and characteristics, identification and databases, content, and distribution of microbial pseudogenes are presented. • How pseudogenization contribute to pathogen virulence is highlighted. • Pseudogenes with potential functions in microorganisms are discussed.

The physical and evolutionary energy landscapes of devolved protein sequences corresponding to pseudogenes.

Protein evolution is guided by structural, functional, and dynamical constraints ensuring organismal viability. Pseudogenes are genomic sequences identified in many eukaryotes that lack translational activity due to sequence degradation and thus over time have undergone "devolution." Previously pseudogenized genes sometimes regain their protein-coding function, suggesting they may still encode robust folding energy landscapes despite multiple mutations. We study both the physical folding landscapes of protein sequences corresponding to human pseudogenes using the Associative Memory, Water Mediated, Structure and Energy Model, and the evolutionary energy landscapes obtained using direct coupling analysis (DCA) on their parent protein families. We found that generally mutations that have occurred in pseudogene sequences have disrupted their native global network of stabilizing residue interactions, making it harder for them to fold if they were translated. In some cases, however, energetic frustration has apparently decreased when the functional constraints were removed. We analyzed this unexpected situation for Cyclophilin A, Profilin-1, and Small Ubiquitin-like Modifier 2 Protein. Our analysis reveals that when such mutations in the pseudogene ultimately stabilize folding, at the same time, they likely alter the pseudogenes' former biological activity, as estimated by DCA. We localize most of these stabilizing mutations generally to normally frustrated regions required for binding to other partners.

Structural Features and Physiological Associations of Human 14-3-3ζ Pseudogenes.

There are about 14,000 pseudogenes that are mutated or truncated sequences resembling functional parent genes. About two-thirds of pseudogenes are processed, while others are duplicated. Although initially thought dead, emerging studies indicate they have functional and regulatory roles. We study 14-3-3ζ, an adaptor protein that regulates cytokine signaling and inflammatory diseases, including rheumatoid arthritis, cancer, and neurological disorders. To understand how 14-3-3ζ (gene symbol YWHAZ) performs diverse functions, we examined the human genome and identified nine YWHAZ pseudogenes spread across many chromosomes. Unlike the 32 kb exon-to-exon sequence in YWHAZ, all pseudogenes are much shorter and lack introns. Out of six, four YWHAZ exons are highly conserved, but the untranslated region (UTR) shows significant diversity. The putative amino acid sequence of pseudogenes is 78-97% homologous, resulting in striking structural similarities with the parent protein. The OMIM and Decipher database searches revealed chromosomal loci containing pseudogenes are associated with human diseases that overlap with the parent gene. To the best of our knowledge, this is the first report on pseudogenes of the 14-3-3 family protein and their implications for human health. This bioinformatics-based study introduces a new insight into the complexity of 14-3-3ζ's functions in biology.

Many purported pseudogenes in bacterial genomes are bona fide genes.

BACKGROUND: Microbial genomes are largely comprised of protein coding sequences, yet some genomes contain many pseudogenes caused by frameshifts or internal stop codons. These pseudogenes are believed to result from gene degradation during evolution but could also be technical artifacts of genome sequencing or assembly. RESULTS: Using a combination of observational and experimental data, we show that many putative pseudogenes are attributable to errors that are incorporated into genomes during assembly. Within 126,564 publicly available genomes, we observed that nearly identical genomes often substantially differed in pseudogene counts. Causal inference implicated assembler, sequencing platform, and coverage as likely causative factors. Reassembly of genomes from raw reads confirmed that each variable affects the number of putative pseudogenes in an assembly. Furthermore, simulated sequencing reads corroborated our observations that the quality and quantity of raw data can significantly impact the number of pseudogenes in an assembler dependent fashion. The number of unexpected pseudogenes due to internal stops was highly correlated (R2 = 0.96) with average nucleotide identity to the ground truth genome, implying relative pseudogene counts can be used as a proxy for overall assembly correctness. Applying our method to assemblies in RefSeq resulted in rejection of 3.6% of assemblies due to significantly elevated pseudogene counts. Reassembly from real reads obtained from high coverage genomes showed considerable variability in spurious pseudogenes beyond that observed with simulated reads, reinforcing the finding that high coverage is necessary to mitigate assembly errors. CONCLUSIONS: Collectively, these results demonstrate that many pseudogenes in microbial genome assemblies are actually genes. Our results suggest that high read coverage is required for correct assembly and indicate an inflated number of pseudogenes due to internal stops is indicative of poor overall assembly quality.

Evidence That Aquaporin 11 (AQP11) in the Spiny Dogfish (Squalus acanthias) May Represent a Pseudogene.

Various attempts to amplify an AQP11 cDNA from tissues of the spiny dogfish (Squalus acanthias) were made. Two pairs of deoxy-inosine-containing degenerate primers were designed based on conserved amino acid sequences from an AQP11 alignment. These primers yielded some faint bands from gill cDNA that were sequenced. Blast searches with the sequences showed they were not AQP11. An elasmobranch AQP11 nucleotide sequence alignment was produced to identify conserved regions to make further degenerate primers. One primer pair produced a short 148 bp fragment showing particularly strong amplification in gill and intestine. It was sequenced and represented a piece of the AQP11 gene. However, as the fragment may have resulted from contaminating genomic DNA (in total RNA used to make cDNA), 5' and 3' RACE were performed to amplify the two ends of the putative cDNA. Furthermore, 5' and 3' RACE amplifications depend on the presence of a 5' cap nucleotide and a poly A tail, respectively on the putative AQP11 mRNA. Hence, successful amplification was only possible from cDNA and not genomic DNA. Nested RACE amplifications were performed using gill and intestinal RACE cDNA, but none of the DNA fragments sequenced were AQP11. Consequently, the spiny dogfish AQP11 gene may represent a pseudogene.

UBDP1 pseudogene and UBD network competitively bind miR­6072 to promote glioma progression.

Increasing evidence suggests that pseudogenes play crucial roles in various cancers, yet their functions and regulatory mechanisms in glioma pathogenesis remain enigmatic. In the present study, a novel pseudogene was identified, UBDP1, which is significantly upregulated in glioblastoma and positively correlated with the expression of its parent gene, UBD. Additionally, high levels of these paired genes are linked with a poor prognosis for patients. In the present study, clinical samples were collected followed by various analyses including microarray for long non­coding RNAs, reverse transcription­quantitative PCR, fluorescence in situ hybridization and western blotting. Cell lines were authenticated and cultured then subjected to various assays for proliferation, migration, and invasion to investigate the molecular mechanisms. Bioinformatic tools identified miRNA targets, and luciferase reporter assays validated these interactions. A tumor xenograft model in mice was used for in vivo studies. In vitro and in vivo studies have demonstrated that UBDP1, localized in the cytoplasm, functions as a tumor­promoting factor influencing cell proliferation, migration, invasion and tumor growth. Mechanistic investigations have indicated that UBDP1 exerts its oncogenic effects by decoying miR­6072 from UBD mRNA, thus forming a competitive endogenous RNA network, which results in the enhanced oncogenic activity of UBD. The present findings offered new insights into the role of pseudogenes in glioma progression, suggesting that targeting the UBDP1/miR­6072/UBD network may serve as a potential therapeutic strategy for glioma patients.

Pathological roles of miRNAs and pseudogene-derived lncRNAs in human cancers, and their comparison as prognosis/diagnosis biomarkers.

This review examines and compares the diagnostic and prognostic capabilities of miRNAs and lncRNAs derived from pseudogenes in cancer patients. Additionally, it delves into their roles in cancer pathogenesis. Both miRNAs and pseudogene-derived lncRNAs have undergone thorough investigation as remarkably sensitive and specific cancer biomarkers, offering significant potential for cancer detection and monitoring. . Extensive research is essential to gain a complete understanding of the precise roles these non-coding RNAs play in cancer, allowing the development of novel targeted therapies and biomarkers for improved cancer detection and treatment approaches.

Construction and analysis of pseudogene-related ceRNA network in breast cancer.

Breast cancer (BC) is one of the leading causes of cancer-related deaths in women. The present study explored the potential role of pseudogenes in BC via construction and analysis of a competing endogenous RNA (ceRNA) network through a three-step process. First, we screened differentially expressed genes in nine BC datasets. Then the gene-pseudogenes pairs (nine hub genes) were selected according to the functional enrichment and correlation analysis. Second, the candidate hub genes and interacting miRNAs were used to construct the ceRNA network. Further analysis of the ceRNA network revealed a crucial ceRNA module with two genes-pseudogene pairs and two miRNAs. The in-depth analysis identified the GBP1/hsa-miR-30d-5p/GBP1P1 axis as a potential tumorigenic axis in BC patients. In the third step, the GBP1/hsa-miR-30d-5p/GBP1P1 axis expression level was assessed in 40 tumor/normal BC patients and MCF-7 cell lines. The expression of GBP1 and GBP1P1 was significantly higher in the tumor compared to the normal tissue. However, the expression of hsa-miR-30d-5p was lower in tumor samples. Then, we introduced the GBP1P1 pseudogene into the MCF-7 cell line to evaluate its effect on GBP1 and hsa-miR-30d-5p expression. As expected, the GBP1 level increased while the hsa-miR-30d-5p level decreased in the GBP1P1-overexprsssing cell line. In addition, the oncogenic properties of MCF-7 (cell viability, clonogenicity, and migration) were improved after GBP1P1 overexpression. In conclusion, we report a ceRNA network that may provide new insight into the role of pseudogenes in BC development.

Prediction of Tumor Microenvironment Characteristics and Treatment Response in Lung Squamous Cell Carcinoma by Pseudogene OR7E47P-related Immune Genes.

OBJECTIVE: Pseudogenes are initially regarded as nonfunctional genomic sequences, but some pseudogenes regulate tumor initiation and progression by interacting with other genes to modulate their transcriptional activities. Olfactory receptor family 7 subfamily E member 47 pseudogene (OR7E47P) is expressed broadly in lung tissues and has been identified as a positive regulator in the tumor microenvironment (TME) of lung adenocarcinoma (LUAD). This study aimed to elucidate the correlation between OR7E47P and tumor immunity in lung squamous cell carcinoma (LUSC). METHODS: Clinical and molecular information from The Cancer Genome Atlas (TCGA) LUSC cohort was used to identify OR7E47P-related immune genes (ORIGs) by weighted gene correlation network analysis (WGCNA). Based on the ORIGs, 2 OR7E47P clusters were identified using non-negative matrix factorization (NMF) clustering, and the stability of the clustering was tested by an extreme gradient boosting classifier (XGBoost). LASSO-Cox and stepwise regressions were applied to further select prognostic ORIGs and to construct a predictive model (ORPScore) for immunotherapy. The Botling cohorts and 8 immunotherapy cohorts (the Samstein, Braun, Jung, Gide, IMvigor210, Lauss, Van Allen, and Cho cohorts) were included as independent validation cohorts. RESULTS: OR7E47P expression was positively correlated with immune cell infiltration and enrichment of immune-related pathways in LUSC. A total of 57 ORIGs were identified to classify the patients into 2 OR7E47P clusters (Cluster 1 and Cluster 2) with distinct immune, mutation, and stromal programs. Compared to Cluster 1, Cluster 2 had more infiltration by immune and stromal cells, lower mutation rates of driver genes, and higher expression of immune-related proteins. The clustering performed well in the internal and 5 external validation cohorts. Based on the 7 ORIGs (HOPX, STX2, WFS, DUSP22, SLFN13, GGCT, and CCSER2), the ORPScore was constructed to predict the prognosis and the treatment response. In addition, the ORPScore was a better prognostic factor and correlated positively with the immunotherapeutic response in cancer patients. The area under the curve values ranged from 0.584 to 0.805 in the 6 independent immunotherapy cohorts. CONCLUSION: Our study suggests a significant correlation between OR7E47P and TME modulation in LUSC. ORIGs can be applied to molecularly stratify patients, and the ORPScore may serve as a biomarker for clinical decision-making regarding individualized prognostication and immunotherapy.

Comprehensive Identification of Mitochondrial Pseudogenes (NUMTs) in the Human Telomere-to-Telomere Reference Genome.

Practices related to mitochondrial research have long been hindered by the presence of mitochondrial pseudogenes within the nuclear genome (NUMTs). Even though partially assembled human reference genomes like hg38 have included NUMTs compilation, the exhaustive NUMTs within the only complete reference genome (T2T-CHR13) remain unknown. Here, we comprehensively identified the fixed NUMTs within the reference genome using human pan-mitogenome (HPMT) from GeneBank. The inclusion of HPMT serves the purpose of establishing an authentic mitochondrial DNA (mtDNA) mutational spectrum for the identification of NUMTs, distinguishing it from the polymorphic variations found in NUMTs. Using HPMT, we identified approximately 10% of additional NUMTs in three human reference genomes under stricter thresholds. And we also observed an approximate 6% increase in NUMTs in T2T-CHR13 compared to hg38, including NUMTs on the short arms of chromosomes 13, 14, and 15 that were not assembled previously. Furthermore, alignments based on 20-mer from mtDNA suggested the presence of more mtDNA-like short segments within the nuclear genome, which should be avoided for short amplicon or cell free mtDNA detection. Finally, through the assay of transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) on cell lines before and after mtDNA elimination, we concluded that NUMTs have a minimal impact on bulk ATAC-seq, even though 16% of sequencing data originated from mtDNA.

Tracking the Diversity and Chromosomal Distribution of the Olfactory Receptor Gene Repertoires of Three Anurans Species.

Olfaction is a crucial capability for most vertebrates and is realized through olfactory receptors in the nasal cavity. The enormous diversity of olfactory receptors has been created by gene duplication, following a birth-and-death model of evolution. The olfactory receptor genes of the amphibians have received relatively little attention up to now, although recent studies have increased the number of species for which data are available. This study analyzed the diversity and chromosomal distribution of the OR genes of three anuran species (Engystomops pustulosus, Bufo bufo and Hymenochirus boettgeri). The OR genes were identified through searches for homologies, and sequence filtering and alignment using bioinformatic tools and scripts. A high diversity of OR genes was found in all three species, ranging from 917 in B. bufo to 1194 in H. boettgeri, and a total of 2076 OR genes in E. pustulosus. Six OR groups were recognized using an evolutionary gene tree analysis. While E. pustulosus has one of the highest numbers of genes of the gamma group (which detect airborne odorants) yet recorded in an anuran, B. bufo presented the smallest number of pseudogene sequences ever identified, with no pseudogenes in either the beta or epsilon groups. Although H. boettgeri shares many morphological adaptations for an aquatic lifestyle with Xenopus, and presented a similar number of genes related to the detection of water-soluble odorants, it had comparatively far fewer genes related to the detection of airborne odorants. This study is the first to describe the complete OR repertoire of the three study species and represents an important contribution to the understanding of the evolution and function of the sense of smell in vertebrates.

Unravelling the due importance of pseudogenes and their resurrection in plants.

The complexities of a genome are underpinned to the vast expanses of the intergenic region, which constitutes ∼97-98% of the genome. This region is essentially composed of what is colloquially referred to as the "junk DNA" and is composed of various elements like transposons, repeats, pseudogenes, etc. The latter have long been considered as dead elements merely contributing to transcriptional noise in the genome. Many studies now describe the previously unknown regulatory functions of these genes. Recent advances in the Next-generation sequencing (NGS) technologies have allowed unprecedented access to these regions. With the availability of whole genome sequences of more than 788 different plant species in past 20 years, genome annotation has become feasible like never before. Different bioinformatic pipelines are available for the identification of pseudogenes. However, still little is known about their biological functions. The functional validation of these genes remains challenging and research in this area is still in infancy, particularly in plants. CRISPR/Cas-based genome editing could provide solutions to understand the biological roles of these genes by allowing creation of precise edits within these genes. The possibility of pseudogene reactivation or resurrection as has been demonstrated in a few studies might open new avenues of genetic manipulation to yield a desirable phenotype. This review aims at comprehensively summarizing the progress made with regards to the identification of pseudogenes and understanding their biological functions in plants.

Human VDAC pseudogenes: an emerging role for VDAC1P8 pseudogene in acute myeloid leukemia.

BACKGROUND: Voltage-dependent anion selective channels (VDACs) are the most abundant mitochondrial outer membrane proteins, encoded in mammals by three genes, VDAC1, 2 and 3, mostly ubiquitously expressed. As 'mitochondrial gatekeepers', VDACs control organelle and cell metabolism and are involved in many diseases. Despite the presence of numerous VDAC pseudogenes in the human genome, their significance and possible role in VDAC protein expression has not yet been considered. RESULTS: We investigated the relevance of processed pseudogenes of human VDAC genes, both in physiological and in pathological contexts. Using high-throughput tools and querying many genomic and transcriptomic databases, we show that some VDAC pseudogenes are transcribed in specific tissues and pathological contexts. The obtained experimental data confirm an association of the VDAC1P8 pseudogene with acute myeloid leukemia (AML). CONCLUSIONS: Our in-silico comparative analysis between the VDAC1 gene and its VDAC1P8 pseudogene, together with experimental data produced in AML cellular models, indicate a specific over-expression of the VDAC1P8 pseudogene in AML, correlated with a downregulation of the parental VDAC1 gene.

A Prognostic Model of Pseudogenes in Acute Myeloid Leukemia.

BACKGROUND: Acute myeloid leukemia (AML) is a hematologic malignancy characterized by the abnormal proliferation of myeloid hematopoietic cells and it is urgently needed to develop new molecular biomarkers to predict clinical outcomes and improve therapeutic effects. METHODS: The differentially expressed genes were identified by comparing TCGA with GETx data. Univariate LASSO and multivariate cox regression analysis were performed to identify prognosis-associated pseudogenes. Based on the overall survival of related pseudogenes, we used them to construct a prognostic model for AML patients. Moreover, we built the pseudogenes-miRNA-mRNA ceRNA networks and explored their involved biological functions and pathways via GO and KEGG enrichment analysis. RESULTS: Seven prognosis-associated pseudogenes were identified, including CCDC150P1, DPY19L1P1, FTH1P8, GTF2IP4, HLA-K, NAPSB, and PDCD6IPP2. The risk model based on these 7 pseudogenes could accurately predict the 1-year, 3-year, and 5-year survival rates. The GO and KEGG enrichment analyses demonstrated that these prognosis-associated pseudogenes were significantly enriched in cell cycle, myeloid leukocyte differentiation, regulation of hemopoiesis, and other critical cancer-related biological functions and pathways. We systematically and comprehensively analyzed the prognostic role of pseudogenes in AML. CONCLUSIONS: The prognostic model of pseudogenes we identified is an independent predictor of overall survival in AML and could be used as biomarker for AML treatment.

Inferring pseudogene-MiRNA associations based on an ensemble learning framework with similarity kernel fusion.

Accumulating evidence shows that pseudogenes can function as microRNAs (miRNAs) sponges and regulate gene expression. Mining potential interactions between pseudogenes and miRNAs will facilitate the clinical diagnosis and treatment of complex diseases. However, identifying their interactions through biological experiments is time-consuming and labor intensive. In this study, an ensemble learning framework with similarity kernel fusion is proposed to predict pseudogene-miRNA associations, named ELPMA. First, four pseudogene similarity profiles and five miRNA similarity profiles are measured based on the biological and topology properties. Subsequently, similarity kernel fusion method is used to integrate the similarity profiles. Then, the feature representation for pseudogenes and miRNAs is obtained by combining the pseudogene-pseudogene similarities, miRNA-miRNA similarities. Lastly, individual learners are performed on each training subset, and the soft voting is used to yield final decision based on the prediction results of individual learners. The k-fold cross validation is implemented to evaluate the prediction performance of ELPMA method. Besides, case studies are conducted on three investigated pseudogenes to validate the predict performance of ELPMA method for predicting pseudogene-miRNA interactions. Therefore, all experiment results show that ELPMA model is a feasible and effective tool to predict interactions between pseudogenes and miRNAs.

A Pathogenic Variant Reclassified to the Pseudogene PMS2P1 in a Patient with Suspected Hereditary Cancer.

The PMS2 gene is involved in DNA repair by the mismatch repair pathway. Deficiencies in this mechanism have been associated with Lynch Syndrome (LS), which is characterized by a high risk for colorectal, endometrial, ovarian, breast, and other cancers. Germinal pathogenic variants of PMS2 are associated with up to 5% of all cases of LS. The prevalence is overestimated for the existence of multiple homologous pseudogenes. We report the case of a 44-year-old woman diagnosed with breast cancer at 34 years without a relevant cancer family history. The presence of pathogenic variant NM_000535.7:c.1A > T, (p.Met1Leu) in PMS2 was determined by next-generation sequencing analysis with a panel of 322 cancer-associated genes and confirmed by capillary sequencing in the patient. The variant was determined in six family members (brothers, sisters, and a son) and seven non-cancerous unrelated individuals. Analysis of the amplified region showed high homology of PMS2 with five of its pseudogenes. We determined that the variant is associated with the PMS2P1 pseudogene following sequence alignment analysis. We propose considering the variant c.1A > T, (p.Met1Leu) in PMS2 for reclassification as not hereditary cancer-related, given the impact on the diagnosis and treatment of cancer patients and families carrying this variant.

Functional Characterization of a Phf8 Processed Pseudogene in the Mouse Genome.

Most pseudogenes are generated when an RNA transcript is reverse-transcribed and integrated into the genome at a new location. Pseudogenes are often considered as an imperfect and silent copy of a functional gene because of the accumulation of numerous mutations in their sequence. Here we report the presence of Pfh8-ps, a Phf8 retrotransposed pseudogene in the mouse genome, which has no disruptions in its coding sequence. We show that this pseudogene is mainly transcribed in testis and can produce a PHF8-PS protein in vivo. As the PHF8-PS protein has a well-conserved JmjC domain, we characterized its enzymatic activity and show that PHF8-PS does not have the intrinsic capability to demethylate H3K9me2 in vitro compared to the parental PHF8 protein. Surprisingly, PHF8-PS does not localize in the nucleus like PHF8, but rather is mostly located at the cytoplasm. Finally, our proteomic analysis of PHF8-PS-associated proteins revealed that PHF8-PS interacts not only with mitochondrial proteins, but also with prefoldin subunits (PFDN proteins) that deliver unfolded proteins to the cytosolic chaperonin complex implicated in the folding of cytosolic proteins. Together, our findings highlighted PHF8-PS as a new pseudogene-derived protein with distinct molecular functions from PHF8.

The pseudogene DUXAP10 contributes to gefitinib resistance in NSCLC by repressing OAS2 expression.

Gefitinib, an epidermal growth factor receptor-tyrosine kinase inhibitor (EGFR-TKI),is the currently recommended first-line therapy for advanced EGFR-mutant lung cancer, and understanding the mechanism of resistance is the key to formulating therapeutic strategies for EGFR-TKIs. In this study, we evaluate the expression patterns and potential biological functions of the pseudogene DUXAP10 in gefitinib resistance. We find that pseudogene DUXAP10 expression is significantly upregulated in NSCLC gefitinib-resistant cells and tissues. Gain and loss of function assays reveal that knockdown of DUXAP10 by siRNA reverses gefitinib resistance both in vitro and in vivo. Furthermore, DUXAP10 interacts with the histone methyltransferase enhancer of zeste homolog 2 (EZH2) to repress the expression of 2',5'-oligoadenylate synthetase (OAS2). Overall, our study highlights the pivotal role of DUXAP10 in gefitinib resistance, and the DUXAP10/EZH2/OAS2 axis might be a promising therapeutic target to overcome acquired gefitinib resistance in NSCLC.

Global DNA and protein interactomes of FLT1P1 (Fms-related tyrosine kinase 1 pseudogene 1) revealed its molecular regulatory functions associated with preeclampsia.

BACKGROUND: Preeclampsia (PE) is one of the most serious pregnancy complications with unknown pathogenesis. Emerging evidence has demonstrated that Fms-related tyrosine kinase 1 (FLT1) is highly involved in PE development. As a pseudogene of FLT1, FLT1P1 increased in PE samples. However, its functions remain largely unknown. METHODS AND RESULTS: In this study, co-expression analysis was performed to identify the potential target genes of FTL1P1. Then chromatin isolation using RNA purification (ChIRP) method was employed to explore the interactomes of FLT1P1, including interacting with DNA fragments and proteins. We found that in PE samples, both FLT1P1 and FLT1 were highly expressed and closely correlated. ChIRP-protein data revealed that FLT1P1 interacts with translation- and transcription-related proteins, including 4 transcription factors (TFs). ChIRP-DNA analysis revealed that FLT1P1 preferentially interacted with DNA fragments downstream of transcription start sites (TSSs). Functional analysis of its interacting genes revealed that they were enriched in transcriptional regulation and apoptosis-related pathways. Twenty-six TFs, including CREB1 and SRF, were extracted from the potential FLT1P1-interacting gene sets and were potential targets of FLT1P1. CREB1 could bind to FLT1 promoter, and was negatively correlated with FLT1 at the expression level, making it a potential regulator of FLT1. CONCLUSIONS: Our study extensively investigated the interactome profiles of FLT1P1, especially the prompter region of TF gene CREB1, and revealed the potential molecular regulatory mechanisms of FLT1 expression in PE samples. Our results provide a novel view of PE pathogenesis, and suggest that FLT1P1 could serve as a potential therapeutic target in PE diagnosis and treatment.

Differences in pseudogene evolution contributed to the contrasting flavors of turnip and Chiifu, two Brassica rapa subspecies.

Pseudogenes are important resources for investigation of genome evolution and genomic diversity because they are nonfunctional but have regulatory effects that influence plant adaptation and diversification. However, few systematic comparative analyses of pseudogenes in closely related species have been conducted. Here, we present a turnip (Brassica rapa ssp. rapa) genome sequence and characterize pseudogenes among diploid Brassica species/subspecies. The results revealed that the number of pseudogenes was greatest in Brassica oleracea (CC genome), followed by B. rapa (AA genome) and then Brassica nigra (BB genome), implying that pseudogene differences emerged after species differentiation. In Brassica AA genomes, pseudogenes were distributed asymmetrically on chromosomes because of numerous chromosomal insertions/rearrangements, which contributed to the diversity among subspecies. Pseudogene differences among subspecies were reflected in the flavor-related glucosinolate (GSL) pathway. Specifically, turnip had the highest content of pungent substances, probably because of expansion of the methylthioalkylmalate synthase-encoding gene family in turnips; these genes were converted into pseudogenes in B. rapa ssp. pekinensis (Chiifu). RNA interference-based silencing of the gene encoding 2-oxoglutarate-dependent dioxygenase 2, which is also associated with flavor and anticancer substances in the GSL pathway, resulted in increased abundance of anticancer compounds and decreased pungency of turnip and Chiifu. These findings revealed that pseudogene differences between turnip and Chiifu influenced the evolution of flavor-associated GSL metabolism-related genes, ultimately resulting in the different flavors of turnip and Chiifu.