Vertebrate pseudogenes

Unveiling the functional nature of retrogenes in dinoflagellates

Retroposition is a gene duplication mechanism that uses RNA molecules as intermediaries to generate new gene copies. Dinoflagellates are proposed as an ideal model for exploring this process due to the tagging of retrogenes with DNA-encoded remnants of the dinoflagellate-specific splice-leader motif at their 5' end. We conducted a comprehensive search for retrogenes in dinoflagellate transcriptomes to uncover their functional nature and the processes underlying their redundancy. We obtained a high-confidence set of hypothetical functional retrogenes widespread through the dinoflagellate lineage. Through annotations and gene ontology enrichment analysis, we found that the functional diversity of retrogenes reflects the most prevalent and active processes during stress periods, particularly those involving post-translational modifications and cell signalling pathways. Additionally, the significant presence of retrogenes linked to specific biological processes involved in symbiosis and toxin production underscores the role of retrogenes in adaptation. The expression profile and codon composition similar to protein-coding genes confirm the operational status of retrogenes and strengthen the idea that retrogenes recapitulate parental gene expression and function. This study provides new evidence supporting widespread gene retroposition across dinoflagellates and highlights the functional link of retrogenes with the core activity of the cell.

Vitamin C: From Self-Sufficiency to Dietary Dependence in the Framework of Its Biological Functions and Medical Implications

Vitamin C is an organic compound biosynthesized in plants and most vertebrates. Since its discovery, the benefits of vitamin C use in the cure and prevention of various pathologies have been frequently reported, including its anti-oxidant, anti-inflammatory, anticoagulant, and immune modulatory properties. Vitamin C plays an important role in collagen synthesis and subsequent scurvy prevention. It is also required in vivo as a cofactor for enzymes involved in carnitine and catecholamine norepinephrine biosynthesis, peptide amidation, and tyrosine catabolism. Moreover, as an enzymatic cofactor, vitamin C is involved in processes of gene transcription and epigenetic regulation. The absence of the synthesis of L-gulono-1,4-lactone oxidase, a key enzyme in the pathway of vitamin C synthesis, is an inborn metabolism error in some fishes and several bird and mammalian species, including humans and non-human primates; it is caused by various changes in the structure of the original GULO gene, making these affected species dependent on external sources of vitamin C. The evolutionary cause of GULO gene pseudogenization remains controversial, as either dietary supplementation or neutral selection is evoked. An evolutionary improvement in the control of redox homeostasis was also considered, as potentially toxic H2O2 is generated as a byproduct in the vitamin C biosynthesis pathway. The inactivation of the GULO gene and the subsequent reliance on dietary vitamin C may have broader implications for aging and age-related diseases, as one of the most important actions of vitamin C is as an anti-oxidant. Therefore, an important aim for medical professionals regarding human and animal health should be establishing vitamin C homeostasis in species that are unable to synthesize it themselves, preventing pathologies such as cardiovascular diseases, cognitive decline, and even cancer.

Genome-scale evolution in local populations of wild chimpanzees

Analysis of genome-scale evolution has been difficult in large, endangered animals because opportunities to collect high-quality genetic samples are limited. There is a need for novel field-friendly, cost-effective genetic techniques. This study conducted an exome-wide analysis of a total of 42 chimpanzees (Pan troglodytes) across six African regions, providing insights into population discrimination techniques. Wild chimpanzee DNA was extracted noninvasively from collected fecal samples using the lysis-buffer storage method. To target genome-scale regions of host DNA, exome-capture sequencing was performed using cost-effective baits originally designed for humans (closely related to chimpanzees). Multivariate analysis effectively discriminated differences in local populations, aiding in the identification of samples' geographical origins. Exome-wide heterozygosity was negatively correlated significantly with genome-wide nonsynonymous-synonymous substitution ratios, suggesting that mutation loads exist at the local population level. Exome sequences revealed functional diversity and protein-coding gene divergence. Segregating pseudogenes were comprehensively annotated, with many being population-specific and others shared among populations. Focusing on multicopy chemosensory receptor genes, the segregating pseudogenes OR7D4 (an olfactory receptor) and TAS2R42 (a bitter taste receptor) were shared among western and eastern chimpanzees. Overall, our analytical framework offers ecological insights into chimpanzees and may be applicable to other organisms.

Molecular evolution of wound healing-related genes during cetacean secondary aquatic adaptation

The marine environment presents challenges for wound healing in cetaceans, despite their remarkable recovery abilities with minimal infections or complications. However, the molecular mechanism underlying this efficient wound healing remains underexplored. To better understand the molecular mechanisms behind wound healing in cetaceans, we investigated the evolutionary patterns of 37 wound healing-related genes in representative mammals. We found wound healing-related genes experience adaptive evolution in cetaceans: (1) Three extrinsic coagulation pathway-related genes-tissue factor (F3), coagulation factor VII (F7), and coagulation factor X (F10)-are subject to positive selection in cetaceans, which might promote efficient hemostasis after injury; positive selection in transforming growth factor-beta 2 (TGF-β2), transforming growth factor-beta 3 (TGF-β3), and platelet-derived growth factor D (PDGFD), which play immunological roles in wound healing, may help cetaceans enhance inflammatory response and tissue debridement. (2) Coagulation factor XII (F12) is the initiation factor in the intrinsic coagulation pathway. It had a premature stop codon mutation and was subjected to selective stress relaxation in cetaceans, suggesting that the early termination of F12 may help cetaceans avoid the risk of vascular blockage during diving. (3) Fibrinogen alpha chain (FGA) and FIII, which were detected to contain the specific amino acid substitutions in marine mammals, indicating similar evolutionary mechanisms might exist among marine mammals to maintain strong wound-healing ability. Thus, our research provides further impetus to study the evolution of the wound healing system in cetaceans and other marine mammals, extending knowledge of preventing coagulation disorder and atherosclerosis in humans.

Expansion and Functional Diversification of Long-Wavelength-Sensitive Opsin in Anabantoid Fishes

Gene duplication is one of the most important sources of novel genotypic diversity and the subsequent evolution of phenotypic diversity. Determining the evolutionary history and functional changes of duplicated genes is crucial for a comprehensive understanding of adaptive evolution. The evolutionary history of visual opsin genes is very dynamic, with repeated duplication events followed by sub- or neofunctionalization. While duplication of the green-sensitive opsins rh2 is common in teleost fish, fewer cases of multiple duplication events of the red-sensitive opsin lws are known. In this study, we investigate the visual opsin gene repertoire of the anabantoid fishes, focusing on the five lws opsin genes found in the genus Betta. We determine the evolutionary history of the lws opsin gene by taking advantage of whole-genome sequences of nine anabantoid species, including the newly assembled genome of Betta imbellis. Our results show that at least two independent duplications of lws occurred in the Betta lineage. The analysis of amino acid sequences of the lws paralogs of Betta revealed high levels of diversification in four of the seven transmembrane regions of the lws protein. Amino acid substitutions at two key-tuning sites are predicted to lead to differentiation of absorption maxima (λmax) between the paralogs within Betta. Finally, eye transcriptomics of B. splendens at different developmental stages revealed expression shifts between paralogs for all cone opsin classes. The lws genes are expressed according to their relative position in the lws opsin cluster throughout ontogeny. We conclude that temporal collinearity of lws expression might have facilitated subfunctionalization of lws in Betta and teleost opsins in general.

The Genome Organization of 5S rRNA Genes in the Model Organism Tribolium castaneum and Its Sibling Species Tribolium freemani

5S ribosomal DNAs (rDNAs) are arranged in tandem and are often under-represented in genome assemblies. In the present study, we performed a global and in-depth analysis of the 5S rDNAs in the model insect Tribolium castaneum and its closely related species Tribolium freemani. To accomplish this goal, we used our recently published genome assemblies based on Nanopore and PacBio long-read sequencing. Although these closely related species share the 5S rRNA gene sequence with high homology, they show a different organization of the 5S rDNA locus. Analysis of 5S rDNA arrays in T. castaneum revealed a typical tandemly repeated organization characterized by repeat units consisting of the 121 bp long 5S rRNA gene and the 71 bp long nontranscribed spacer (NTS). In contrast, T. freemani showed a much more complex organization of 5S rDNA arrays characterized by two patterns. The first is based on the association of 5S rRNA gene with arrays of a satellite DNA, representing the NTS sequence of the 5S rDNA genes in T. freemani. The second, more complex type is characterized by a somewhat less frequent occurrence of the 5S rRNA gene and its association with longer satellite DNA arrays that are regularly interrupted by Jockey-like retrotransposons. This organization, in which the ribosomal gene is associated with two completely different repetitive elements such as satellite DNAs and retrotransposons, suggests that the 5S rRNA gene, regardless of its crucial function in the genome, could be a subject of extremely dynamic genomic rearrangements.

Illuminating the function of the orphan transporter, SLC22A10, in humans and other primates

SLC22A10 is an orphan transporter with unknown substrates and function. The goal of this study is to elucidate its substrate specificity and functional characteristics. In contrast to orthologs from great apes, human SLC22A10, tagged with green fluorescent protein, is not expressed on the plasma membrane. Cells expressing great ape SLC22A10 orthologs exhibit significant accumulation of estradiol-17β-glucuronide, unlike those expressing human SLC22A10. Sequence alignments reveal a proline at position 220 in humans, which is a leucine in great apes. Replacing proline with leucine in SLC22A10-P220L restores plasma membrane localization and uptake function. Neanderthal and Denisovan genomes show proline at position 220, akin to modern humans, indicating functional loss during hominin evolution. Human SLC22A10 is a unitary pseudogene due to a fixed missense mutation, P220, while in great apes, its orthologs transport sex steroid conjugates. Characterizing SLC22A10 across species sheds light on its biological role, influencing organism development and steroid homeostasis.

The Ribosome Hypothesis: Decoding Mood Disorder Complexity

Several types of mood disorders lie along a continuum, with nebulous boundaries between them. Understanding the mechanisms that contribute to mood disorder complexity is critical for effective treatment. However, present treatments are largely centered around neurotransmission and receptor-based hypotheses, which, given the high instance of treatment resistance, fail to adequately explain the complexities of mood disorders. In this opinion piece, based on our recent results, we propose a ribosome hypothesis of mood disorders. We suggest that any hypothesis seeking to explain the diverse nature of mood disorders must incorporate infrastructure diversity that results in a wide range of effects. Ribosomes, with their mobility across neurites and complex composition, have the potential to become specialized during stress; thus, ribosome diversity and dysregulation are well suited to explaining mood disorder complexity. Here, we first establish a framework connecting ribosomes to the current state of knowledge associated with mood disorders. Then, we describe the potential mechanisms through which ribosomes could homeostatically regulate systems to manifest diverse mood disorder phenotypes and discuss approaches for substantiating the ribosome hypothesis. Investigating these mechanisms as therapeutic targets holds promise for transdiagnostic avenues targeting mood disorders.

Divergence time estimates for the hypoxia-inducible factor-1 alpha (HIF1α) reveal an ancient emergence of animals in low-oxygen environments

Unveiling the tempo and mode of animal evolution is necessary to understand the links between environmental changes and biological innovation. Although the earliest unambiguous metazoan fossils date to the late Ediacaran period, molecular clock estimates agree that the last common ancestor (LCA) of all extant animals emerged ~850 Ma, in the Tonian period, before the oldest evidence for widespread ocean oxygenation at ~635-560 Ma in the Ediacaran period. Metazoans are aerobic organisms, that is, they are dependent on oxygen to survive. In low-oxygen conditions, most animals have an evolutionarily conserved pathway for maintaining oxygen homeostasis that triggers physiological changes in gene expression via the hypoxia-inducible factor (HIFa). However, here we confirm the absence of the characteristic HIFa protein domain responsible for the oxygen sensing of HIFa in sponges and ctenophores, indicating the LCA of metazoans lacked the functional protein domain as well, and so could have maintained their transcription levels unaltered under the very low-oxygen concentrations of their environments. Using Bayesian relaxed molecular clock dating, we inferred that the ancestral gene lineage responsible for HIFa arose in the Mesoproterozoic Era, ~1273 Ma (Credibility Interval 957-1621 Ma), consistent with the idea that important genetic machinery associated with animals evolved much earlier than the LCA of animals. Our data suggest at least two duplication events in the evolutionary history of HIFa, which generated three vertebrate paralogs, products of the two successive whole-genome duplications that occurred in the vertebrate LCA. Overall, our results support the hypothesis of a pre-Tonian emergence of metazoans under low-oxygen conditions, and an increase in oxygen response elements during animal evolution.

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.

Interchromosomal Colocalization with Parental Genes Is Linked to the Function and Evolution of Mammalian Retrocopies

Retrocopies are gene duplicates arising from reverse transcription of mature mRNA transcripts and their insertion back into the genome. While long being regarded as processed pseudogenes, more and more functional retrocopies have been discovered. How the stripped-down retrocopies recover expression capability and become functional paralogs continually intrigues evolutionary biologists. Here, we investigated the function and evolution of retrocopies in the context of 3D genome organization. By mapping retrocopy-parent pairs onto sequencing-based and imaging-based chromatin contact maps in human and mouse cell lines and onto Hi-C interaction maps in 5 other mammals, we found that retrocopies and their parental genes show a higher-than-expected interchromosomal colocalization frequency. The spatial interactions between retrocopies and parental genes occur frequently at loci in active subcompartments and near nuclear speckles. Accordingly, colocalized retrocopies are more actively transcribed and translated and are more evolutionarily conserved than noncolocalized ones. The active transcription of colocalized retrocopies may result from their permissive epigenetic environment and shared regulatory elements with parental genes. Population genetic analysis of retroposed gene copy number variants in human populations revealed that retrocopy insertions are not entirely random in regard to interchromosomal interactions and that colocalized retroposed gene copy number variants are more likely to reach high frequencies, suggesting that both insertion bias and natural selection contribute to the colocalization of retrocopy-parent pairs. Further dissection implies that reduced selection efficacy, rather than positive selection, contributes to the elevated allele frequency of colocalized retroposed gene copy number variants. Overall, our results hint a role of interchromosomal colocalization in the "resurrection" of initially neutral retrocopies.

Functional implications of paralog genes in polyglutamine spinocerebellar ataxias

Polyglutamine (polyQ) spinocerebellar ataxias (SCAs) comprise a group of autosomal dominant neurodegenerative disorders caused by (CAG/CAA)n expansions. The elongated stretches of adjacent glutamines alter the conformation of the native proteins inducing neurotoxicity, and subsequent motor and neurological symptoms. Although the etiology and neuropathology of most polyQ SCAs have been extensively studied, only a limited selection of therapies is available. Previous studies on SCA1 demonstrated that ATXN1L, a human duplicated gene of the disease-associated ATXN1, alleviated neuropathology in mice models. Other SCA-associated genes have paralogs (i.e., copies at different chromosomal locations derived from duplication of the parental gene), but their functional relevance and potential role in disease pathogenesis remain unexplored. Here, we review the protein homology, expression pattern, and molecular functions of paralogs in seven polyQ dominant ataxias-SCA1, SCA2, MJD/SCA3, SCA6, SCA7, SCA17, and DRPLA. Besides ATXN1L, we highlight ATXN2L, ATXN3L, CACNA1B, ATXN7L1, ATXN7L2, TBPL2, and RERE as promising functional candidates to play a role in the neuropathology of the respective SCA, along with the parental gene. Although most of these duplicates lack the (CAG/CAA)n region, if functionally redundant, they may compensate for a partial loss-of-function or dysfunction of the wild-type genes in SCAs. We aim to draw attention to the hypothesis that paralogs of disease-associated genes may underlie the complex neuropathology of dominant ataxias and potentiate new therapeutic strategies.

Illuminating the Function of the Orphan Transporter, SLC22A10 in Humans and Other Primates

SLC22A10 is classified as an orphan transporter with unknown substrates and function. Here we describe the discovery of the substrate specificity and functional characteristics of SLC22A10. The human SLC22A10 tagged with green fluorescent protein was found to be absent from the plasma membrane, in contrast to the SLC22A10 orthologs found in great apes. Estradiol-17β-glucuronide accumulated in cells expressing great ape SLC22A10 orthologs (over 4-fold, p<0.001). In contrast, human SLC22A10 displayed no uptake function. Sequence alignments revealed two amino acid differences including a proline at position 220 of the human SLC22A10 and a leucine at the same position of great ape orthologs. Site-directed mutagenesis yielding the human SLC22A10-P220L produced a protein with excellent plasma membrane localization and associated uptake function. Neanderthal and Denisovan genomes show human-like sequences at proline 220 position, corroborating that SLC22A10 were rendered nonfunctional during hominin evolution after the divergence from the pan lineage (chimpanzees and bonobos). These findings demonstrate that human SLC22A10 is a unitary pseudogene and was inactivated by a missense mutation that is fixed in humans, whereas orthologs in great apes transport sex steroid conjugates.

Indications and management of preimplantation genetic testing for monogenic conditions: a committee opinion

This statement is offered to update and expand on the prior American Society for Reproductive Medicine preimplantation genetic testing (PGT) opinion, elucidate the current clinical and technical complexities specific to PGT for monogenic conditions, assist providers in supporting patient understanding of and access to this technology, and offer considerations for the development of future clinical and laboratory guidelines on PGT for monogenic conditions.

Micro Germline-Restricted Chromosome in Blue Tits: Evidence for Meiotic Functions

The germline-restricted chromosome (GRC) is likely present in all songbird species but differs widely in size and gene content. This extra chromosome has been described as either a microchromosome with only limited basic gene content or a macrochromosome with enriched gene functions related to female gonad and embryo development. Here, we assembled, annotated, and characterized the first micro-GRC in the blue tit (Cyanistes caeruleus) using high-fidelity long-read sequencing data. Although some genes on the blue tit GRC show signals of pseudogenization, others potentially have important functions, either currently or in the past. We highlight the GRC gene paralog BMP15, which is among the highest expressed GRC genes both in blue tits and in zebra finches (Taeniopygia guttata) and is known to play a role in oocyte and follicular maturation in other vertebrates. The GRC genes of the blue tit are further enriched for functions related to the synaptonemal complex. We found a similar functional enrichment when analyzing published data on GRC genes from two nightingale species (Luscinia spp.). We hypothesize that these genes play a role in maintaining standard maternal inheritance or in recombining maternal and paternal GRCs during potential episodes of biparental inheritance.

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.

Gapless genome assembly of azalea and multi-omics investigation into divergence between two species with distinct flower color

The genus Rhododendron (Ericaceae), with more than 1000 species highly diverse in flower color, is providing distinct ornamental values and a model system for flower color studies. Here, we investigated the divergence between two parental species with different flower color widely used for azalea breeding. Gapless genome assembly was generated for the yellow-flowered azalea, Rhododendron molle. Comparative genomics found recent proliferation of long terminal repeat retrotransposons (LTR-RTs), especially Gypsy, has resulted in a 125 Mb (19%) genome size increase in species-specific regions, and a significant amount of dispersed gene duplicates (13 402) and pseudogenes (17 437). Metabolomic assessment revealed that yellow flower coloration is attributed to the dynamic changes of carotenoids/flavonols biosynthesis and chlorophyll degradation. Time-ordered gene co-expression networks (TO-GCNs) and the comparison confirmed the metabolome and uncovered the specific gene regulatory changes underpinning the distinct flower pigmentation. B3 and ERF TFs were found dominating the gene regulation of carotenoids/flavonols characterized pigmentation in R. molle, while WRKY, ERF, WD40, C2H2, and NAC TFs collectively regulated the anthocyanins characterized pigmentation in the red-flowered R simsii. This study employed a multi-omics strategy in disentangling the complex divergence between two important azaleas and provided references for further functional genetics and molecular breeding.

Noncoding RNAs in Thyroid-Follicular-Cell-Derived Carcinomas

Simple Summary

Thyroid tumors represent the most common neoplastic pathology of the endocrine system. Mutations occurring in oncogenes and tumor suppressor genes are responsible for thyroid carcinogenesis; however, the complete mutational landscape characterizing these neoplasias has not been completely unveiled. It has been established that only the 2% of the human genome codes for proteins, suggesting that the vast majority of the genome has regulatory capabilities, which, if altered, could account for the onset of cancer. Hence, many scientific efforts are currently focused on the characterization of the heterogeneous class of noncoding RNAs, which represent an abundant part of the transcribed noncoding genome. In this review, we mainly focus on the involvement of microRNAs, long noncoding RNAs, and pseudogenes in thyroid cancer. The determination of the diagnosis, prognosis, and treatment of thyroid cancers based on the evaluation of the noncoding RNA network could allow the implementation of a more personalized approach to fighting these pathologies.

Abstract

Among the thyroid neoplasias originating from follicular cells, we can include well-differentiated carcinomas, papillary (PTC) and follicular (FTC) thyroid carcinomas, and the undifferentiated anaplastic (ATC) carcinomas. Several mutations in oncogenes and tumor suppressor genes have already been observed in these malignancies; however, we are still far from the comprehension of their full regulation-altered landscape. Even if only 2% of the human genome has the ability to code for proteins, most of the noncoding genome is transcribed, constituting the heterogeneous class of noncoding RNAs (ncRNAs), whose alterations are associated with the development of several human diseases, including cancer. Hence, many scientific efforts are currently focused on the elucidation of their biological role. In this review, we analyze the scientific literature regarding the involvement of microRNAs (miRNAs), long noncoding RNAs (lncRNAs), and pseudogenes in FTC, PTC, and ATC. Recent findings emphasized the role of lncRNAs in all steps of cancer progression. In particular, lncRNAs may control progression steps by regulating the expression of genes and miRNAs involved in cell proliferation, apoptosis, epithelial–mesenchymal transition, and metastatization. In conclusion, the determination of the diagnosis, prognosis, and treatment of cancer based on the evaluation of the ncRNA network could allow the implementation of a more personalized approach to fighting thyroid tumors.

Major and minor U small nuclear RNAs genes characterization in a neotropical fish genome: Chromosomal remodeling and repeat units dispersion in Parodontidae

In association with many proteins, small nuclear RNAs (snRNAs) organize the spliceosomes that play a significant role in processing precursor mRNAs during gene expression. According to snRNAs genic arrangements, two kinds of spliceosomes (major and minor) can be organized into eukaryotic cells. Although in situ localization of U1 and U2 snDNAs have been performed in fish karyotypes, studies with genomic characterization and functionality of U snRNAs integrated into chromosomal changes on Teleostei are still scarce. This study aimed to achieve a genomic characterization of the U snRNAs genes in Apareiodon sp. (2n = 54, ZZ/ZW), apply these data to recognize functional/defective copies, and map chromosomal changes involving snDNAs in Parodontidae species karyotype diversification. Nine snRNA multigene families (U1, U2, U4, U5, U6, U11, U12, U4atac and U6atac) arranged in putatively functional copies in the genome were analyzed. Proximal Sequence Elements (PSE) and TATA-box promoters occurrence, besides an entire transcribed region and conserved secondary structures, qualify them for spliceosome activity. In addition, several defective copies or pseudogenes were identified for the snRNAs that make up the major spliceosome. In situ localization of snDNAs in five species of Parodontidae demonstrated that U1, U2, and U4 snDNAs were involved in chromosomal location changes or units dispersion. The U snRNAs defective/pseudogenes units dispersion could be favored by the probable occurrence of active retrotransposition enzymes in the Apareiodon genome. The U2 and U4 snDNAs sites were involved in independent events in the differentiation of sex chromosomes among Parodontidae lineages. The study characterized U snRNA genes that compose major and minor spliceosomes in the Apareiodon sp. genome and proposes that their defective copies trigger chromosome differentiation and diversification events in Parodontidae.

Identification of multiple TAR DNA binding protein retropseudogene lineages during the evolution of primates

The TAR DNA Binding Protein (TARDBP) gene has become relevant after the discovery of its several pathogenic mutations. The lack of evolutionary history is in contrast to the amount of studies found in the literature. This study investigated the evolutionary dynamics associated with the retrotransposition of the TARDBP gene in primates. We identified novel retropseudogenes that likely originated in the ancestors of anthropoids, catarrhines, and lemuriformes, i.e. the strepsirrhine clade that inhabit Madagascar. We also found species-specific retropseudogenes in the Philippine tarsier, Bolivian squirrel monkey, capuchin monkey and vervet. The identification of a retropseudocopy of the TARDBP gene overlapping a lncRNA that is potentially expressed opens a new avenue to investigate TARDBP gene regulation, especially in the context of TARDBP associated pathologies.

Ankfn1-mutant vestibular defects require loss of both ancestral and derived paralogs for penetrance in zebrafish

How and to what degree gene duplication events create regulatory innovation, redundancy, or neofunctionalization remain important questions in animal evolution and comparative genetics. Ankfn1 genes are single copy in most invertebrates, partially duplicated in jawed vertebrates, and only the derived copy retained in most mammals. Null mutations in the single mouse homolog have vestibular and neurological abnormalities. Null mutation of the single Drosophila homolog is typically lethal with severe sensorimotor deficits in rare survivors. The functions and potential redundancy of paralogs in species with two copies are not known. Here, we define a vestibular role for Ankfn1 homologs in zebrafish based on the simultaneous disruption of each locus. Zebrafish with both paralogs disrupted showed vestibular defects and early lethality from swim bladder inflation failure. One intact copy at either locus was sufficient to prevent major phenotypes. Our results show that vertebrate Ankfn1 genes are required for vestibular-related functions, with at least partial redundancy between ancestral and derived paralogs.

Genome-wide identification of seven superoxide dismutase genes in the marine rotifer Brachionus rotundiformis and modulated expression and enzymatic activity in response to microplastics and nutritional status

Microplastics (MPs) pollution has attracted worldwide attention. Superoxide dismutase (SOD) is a sensitive indicator for assessing the toxic effects of MPs in aquatic organisms. However, few studies have been performed to identify all genes encoding SOD in aquatic invertebrates. Especially, effects of MPs on SOD activity and expression in aquatic organisms under starvation or a subsequent refeeding status are unclear. In the present study, all full-length genes encoding SOD were cloned and characterized from the marine rotifer Brachionus rotundiformis, including CuZnSOD1, CuZnSOD2, CuZnSOD3, CuZnSOD4, CuZnSOD5, MnSOD1, and MnSOD2. The CuZnSOD1, CuZnSOD2 and MnSOD2 are homologous to SODs from vertebrates and the other SOD proteins are rotifer-specific according to the results from the phylogenetic tree. The conserved signature sequences and binding sites of Cu²⁺, Zn²⁺and Mn²⁺ were also identified in the seven SOD proteins. Compared with feeding, starvation down-regulated SOD activity and mRNA expression of CuZnSOD2, CuZnSOD4, CuZnSOD5, MnSOD1 and MnSOD2 while refeeding maintained SOD activity comparable to the feeding level and up-regulated CuZnSOD5 and MnSOD2. Intake of MPs by B. rotundiformis was observed by examining fluorescence signals from the fluorescently-labeled microplastics under different nutritional status. Exposure to MPs reduced rotifer density and increased malondialdehyde (MDA) content and SOD activity in the rotifers under the refeeding condition, but did not affect these indicators under the feeding and starvation conditions. However, mRNA expression of some tested genes was responsive to MPs in the fed, starved and refed rotifers. The present study for the first time demonstrated a nutritional status-dependent effect of MPs on oxidative stress response, and provided more sensitive molecular biomarkers for assessing the toxicity of MPs using B. rotundiformis as a model animal.

Nutritional-status dependent effects of microplastics on activity and expression of alkaline phosphatase and alpha-amylase in Brachionus rotundiformis

Tissue-nonspecific alkaline phosphatase (ALPL) and alpha-amylase (AMY) are essential in the immune and digestive systems, respectively. Microplastics (MPs) pose a risk to zooplankton which may be in a state of feeding, starvation, or subsequent refeeding. However, molecular characterization of both enzymes and the regulated mechanisms affected by nutritional statuses and MPs remain unclear in zooplankton. In the present study, four full-length genes encoding ALPL and two genes encoding AMY were cloned and characterized from an isolated marine rotifer, Brachionus rotundiformis, including alplA, alplB, alplC, alplD, amy2a, and amy2al. AMY activity and expression of amy2a and amy2al were reduced by starvation and recovered after refeeding compared with feeding. ALPL activity remained unchanged among different statuses, while alplA, alplB and alplD were down-regulated by starvation and refeeding compared with feeding. ALPL activity was not affected by exposure to 10, 100 and 1000 μg/L MPs in rotifers subjected to feeding, starvation and refeeding, whereas AMY activity was significantly enhanced by 1000 μg/L MPs in rotifers subjected to refeeding. Gene expression of the tested genes, except amy2a, was significantly responsive to MPs, especially in the feeding rotifers, depending on MPs concentrations and nutritional statuses. Two-way ANOVA confirmed that these changes were strongly associated with the interaction between MPs concentrations and nutritional statuses. The present study is the first to demonstrate a nutritional status-dependent impact of MPs on immune and digestive responses, and provides more sensitive molecular biomarkers for assessing MPs toxicity using the species as model animals.

Mining Proteomics Datasets to Uncover Functional Pseudogenes

Analysis of the proteome, combined with the human genome sequence, improved gene annotation and confirmed that some genes are actually expressed as proteins, including pseudogenes, alternative reading frames, and additional protein isoforms. Although sequencing a genome is a challenge in itself, the annotation of all genes encoding proteins is a bigger one. Here, we describe an in silico approach to identify evidence of pseudogene expression, as well as an experimental approach for the validation of the protein encoded by pseudogenes, including the steps necessary to quantify these proteins. This technique enables a comprehensive analysis for the expression of genes into proteins that were not suspected of existing.

The genetic informational network: how DNA conveys semantic information

The question of whether "genetic information" is a merely causal factor in development or can be made sense of semantically, in a way analogous to a language or other type of representation, has generated a long debate in the philosophy of biology. It is intimately connected with another intense debate, concerning the limits of genetic determinism. In this paper I argue that widespread attempts to draw analogies between genetic information and information contained in books, blueprints or computer programs, are fundamentally inadequate. In development, gene exons are the central part of an intricate and densely ramified semantic Genetic Informational Network. DNA in the entire genome is in a state of continuous positive and negative feedback with itself and with its 'environment', and is 'read' and acted upon by the cell in various alternative and complementary ways. The linear combinatorial coding relation between codons and amino acids is but one aspect of semantic genetic information, which is, when considered in its entirety, a far wider and richer concept.

Processed pseudogenes: A substrate for evolutionary innovation: Retrotransposition contributes to genome evolution by propagating pseudogene sequences with rich regulatory potential throughout the genome

Processed pseudogenes may serve as a genetic reservoir for evolutionary innovation. Here, we argue that through the activity of long interspersed element-1 retrotransposons, processed pseudogenes disperse coding and noncoding sequences rich with regulatory potential throughout the human genome. While these sequences may appear to be non-functional, a lack of contemporary function does not prohibit future development of biological activity. Here, we discuss the dynamic evolution of certain processed pseudogenes into coding and noncoding genes and regulatory elements, and their implication in wide-ranging biological and pathological processes. Also see the video abstract here: https://youtu.be/iUY_mteVoPI.

Systematic functional interrogation of human pseudogenes using CRISPRi

BACKGROUND

The human genome encodes over 14,000 pseudogenes that are evolutionary relics of protein-coding genes and commonly considered as nonfunctional. Emerging evidence suggests that some pseudogenes may exert important functions. However, to what extent human pseudogenes are functionally relevant remains unclear. There has been no large-scale characterization of pseudogene function because of technical challenges, including high sequence similarity between pseudogene and parent genes, and poor annotation of transcription start sites.

RESULTS

To overcome these technical obstacles, we develop an integrated computational pipeline to design the first genome-wide library of CRISPR interference (CRISPRi) single-guide RNAs (sgRNAs) that target human pseudogene promoter-proximal regions. We perform the first pseudogene-focused CRISPRi screen in luminal A breast cancer cells and reveal approximately 70 pseudogenes that affect breast cancer cell fitness. Among the top hits, we identify a cancer-testis unitary pseudogene, MGAT4EP, that is predominantly localized in the nucleus and interacts with FOXA1, a key regulator in luminal A breast cancer. By enhancing the promoter binding of FOXA1, MGAT4EP upregulates the expression of oncogenic transcription factor FOXM1. Integrative analyses of multi-omic data from the Cancer Genome Atlas (TCGA) reveal many unitary pseudogenes whose expressions are significantly dysregulated and/or associated with overall/relapse-free survival of patients in diverse cancer types.

CONCLUSIONS

Our study represents the first large-scale study characterizing pseudogene function. Our findings suggest the importance of nuclear function of unitary pseudogenes and underscore their underappreciated roles in human diseases. The functional genomic resources developed here will greatly facilitate the study of human pseudogene function.

Lack of CpG islands in human unitary pseudogenes and its implication

CpG islands (CGIs) are aggregation of CpG dinucleotides in the promoters of mammalian genes. These CGIs are present in almost all the housekeeping genes and some tissue-specific genes in the mammalian genome. Extensive research has been done on the prevalence and role of CGIs in protein-coding genes. However, little is known about CGIs in pseudogenes. In the current research project, we focused on CGIs in three main classes of pseudogenes e.g., duplicated pseudogenes (DPGs), processed pseudogenes (PPGs), and unitary pseudogenes (UPGs). We discovered a predominant absence of CGIs in the promoters of all three pseudogenes. We also compared the CGI profile of these pseudogenes with their parent genes and found that unitary pseudogenes (UPGs) differ from the DPGs and PPGs in the sense that in the latter, lack of CGIs is a consequential event while in UPGs, this lack of CGIs in their promoters is not a result of pseudogenization process. We also discussed the implication of the results obtained from this comparison. To our knowledge, this is the first-ever study highlighting this aspect of UPGs throwing new insights into the evolution of genome in general and especially in the context of pseudogenes.

The absence of thrombin-like activity in Bothrops erythromelas venom is due to the deletion of the snake venom thrombin-like enzyme gene

Snake venom thrombin-like enzymes (SVTLEs) are serine proteinases that clot fibrinogen. SVTLEs are distributed mainly in venoms from snakes of the Viperidae family, comprising venomous pit viper snakes. Bothrops snakes are distributed throughout Central and South American and are responsible for most venomous snakebites. Most Bothrops snakes display thrombin-like activity in their venoms, but it has been shown that some species do not present it. In this work, to understand SVTLE polymorphism in Bothrops snake venoms, we studied individual samples from two species of medical importance in Brazil: Bothrops jararaca, distributed in Southeastern Brazil, which displays coagulant activity on plasma and fibrinogen, and Bothrops erythromelas, found in Northeastern Brazil, which lacks direct fibrinogen coagulant activity but shows plasma coagulant activity. We tested the coagulant activity of venoms and the presence of SVTLE genes by a PCR approach. The SVTLE gene structure in B. jararaca is similar to the Bothrops atrox snake, comprising five exons. We could not amplify SVTLE sequences from B. erythromelas DNA, except for a partial pseudogene. These genes underwent a positive selection in some sites, leading to an amino acid sequence diversification, mostly in exon 2. The phylogenetic tree constructed using SVTLE coding sequences confirms that they are related to the chymotrypsin/kallikrein family. Interestingly, we found a B. jararaca specimen whose venom lacked thrombin-like activity, and its gene sequence was a pseudogene with SVTLE structure, presenting nonsense and frameshift mutations. Our results indicate an association of the lack of thrombin-like activity in B. jararaca and B. erythromelas venoms with mutations and deletions of snake venom thrombin-like enzyme genes.

Not So Dead Genes-Retrocopies as Regulators of Their Disease-Related Progenitors and Hosts

Retroposition is RNA-based gene duplication leading to the creation of single exon nonfunctional copies. Nevertheless, over time, many of these duplicates acquire transcriptional capabilities. In human in most cases, these so-called retrogenes do not code for proteins but function as regulatory long noncoding RNAs (lncRNAs). The mechanisms by which they can regulate other genes include microRNA sponging, modulation of alternative splicing, epigenetic regulation and competition for stabilizing factors, among others. Here, we summarize recent findings related to lncRNAs originating from retrocopies that are involved in human diseases such as cancer and neurodegenerative, mental or cardiovascular disorders. Special attention is given to retrocopies that regulate their progenitors or host genes. Presented evidence from the literature and our bioinformatics analyses demonstrates that these retrocopies, often described as unimportant pseudogenes, are significant players in the cell’s molecular machinery.

Global gene methylation profiling of common warts caused by human papillomaviruses infection

Infection with the human papillomaviruses (HPV) often involves the epigenetic modification of the host genome. Despite its prevalence among the population, host genome methylation in HPV-induced warts is not clearly understood. In this study, genome-wide methylation profiling was carried out on paired healthy skin and wart samples in order to investigate the effects that benign HPV infection has on gene methylation status. To overcome this gap in knowledge, paired wart (n = 12) and normal skin (n = 12) samples were obtained from Arab males in order to perform DNA extraction and subsequent genome-wide methylation profiling on the Infinium Methylation EPIC Bead Chip microarray. Analysis of differential methylation revealed a clear pattern of discrimination between the wart and normal skin samples. In warts, the most differentially methylated (DM) genes included long non-coding RNAs (AC005884, AL049646.2, AC126121.2, AP001790.1, and AC107959.3), microRNAs (MIR374B, MIR596, MIR1255B1, MIR26B, and MIR196A2),snoRNAs (SNORD114-22, SNORD70, and SNORD114-31), pseudogenes (AC069366.1, RNU4ATAC11P, AC120057.1, NANOGP3, AC106038.2, TPT1P2, SDC4P, PKMP3, and VN2R3P), and protein-coding genes (AREG, GJB2, C12orf71, AC020909.2, S100A8, ZBED2, FABP7, and CYSLTR1). In addition, pathway analysis revealed that, among the most differentially methylated genes, STAT5A, RARA, MEF2D, MAP3K8, and THRA were the common regulators. It can be observed that HPV-induced warts involve a clear and unique epigenetic alteration to the host genome.

[Non-coding Natural Antisense RNA: Mechanisms of Action in the Regulation of Target Gene Expression and Its Clinical Implications]

Recent advances in high-throughput technologies have revealed that 75% of the human genome is transcribed to RNA, whereas only 3% of transcripts are translated into proteins. Consequently, many long non-coding RNAs (lncRNAs) have been identified, which has improved our understanding of the complexity of biological processes. LncRNAs comprise multiple classes of RNA transcripts that regulate the transcription, stability and translation of protein-coding genes in a genome. Natural antisense transcripts (NATs) form one such class, and the GENCODE v30 catalog contains 16193 lncRNA loci, of which 5611 are antisense loci. This review outlines our emerging understanding of lncRNAs, with a particular focus on how lncRNAs regulate gene expression using interferon-α1 (IFN-α1) mRNA and its antisense partner IFN-α1 antisense (as)RNA as an example. We have identified and characterized the asRNA that determines post-transcriptional IFN-α1 mRNA levels. IFN-α1 asRNA stabilizes IFN-α1 mRNA by cytoplasmic sense-antisense duplex formation, which may enhance the accessibility of an RNA stabilizer protein or decrease the affinity of an RNA decay factor for the RNA. IFN-α1 asRNA can also act as competing molecules in the competing endogenous (ce)RNA network with other members of the IFNA multigene family mRNAs/asRNAs, and other cellular mRNA transcripts. Furthermore, antisense oligoribonucleotides representing functional domains of IFN-α1 asRNA inhibit influenza virus proliferation in the respiratory tract of virus-infected animals. Thus, these findings support, at least in part, the rationale that dissecting the activity of NAT on gene expression regulation promises to reveal previously unanticipated biology, with potential to provide new therapeutic approaches to diseases.

Amalgamated cross-species transcriptomes reveal organ-specific propensity in gene expression evolution

The origins of multicellular physiology are tied to evolution of gene expression. Genes can shift expression as organisms evolve, but how ancestral expression influences altered descendant expression is not well understood. To examine this, we amalgamate 1,903 RNA-seq datasets from 182 research projects, including 6 organs in 21 vertebrate species. Quality control eliminates project-specific biases, and expression shifts are reconstructed using gene-family-wise phylogenetic Ornstein-Uhlenbeck models. Expression shifts following gene duplication result in more drastic changes in expression properties than shifts without gene duplication. The expression properties are tightly coupled with protein evolutionary rate, depending on whether and how gene duplication occurred. Fluxes in expression patterns among organs are nonrandom, forming modular connections that are reshaped by gene duplication. Thus, if expression shifts, ancestral expression in some organs induces a strong propensity for expression in particular organs in descendants. Regardless of whether the shifts are adaptive or not, this supports a major role for what might be termed preadaptive pathways of gene expression evolution.

PseudoChecker: an integrated online platform for gene inactivation inference

The rapid expansion of high-quality genome assemblies, exemplified by ongoing initiatives such as the Genome-10K and i5k, demands novel automated methods to approach comparative genomics. Of these, the study of inactivating mutations in the coding region of genes, or pseudogenization, as a source of evolutionary novelty is mostly overlooked. Thus, to address such evolutionary/genomic events, a systematic, accurate and computationally automated approach is required. Here, we present PseudoChecker, the first integrated online platform for gene inactivation inference. Unlike the few existing methods, our comparative genomics-based approach displays full automation, a built-in graphical user interface and a novel index, PseudoIndex, for an empirical evaluation of the gene coding status. As a multi-platform online service, PseudoChecker simplifies access and usability, allowing a fast identification of disruptive mutations. An analysis of 30 genes previously reported to be eroded in mammals, and 30 viable genes from the same lineages, demonstrated that PseudoChecker was able to correctly infer 97% of loss events and 95% of functional genes, confirming its reliability. PseudoChecker is freely available, without login required, at http://pseudochecker.ciimar.up.pt.

Comprehensive analysis of LDHAP5 pseudogene expression and potential pathogenesis in ovarian serous cystadenocarcinoma

Background

We aimed to identify differentially expressed pseudogenes and explore their potential functions in four types of common gynecological malignancies (e.g., cervical squamous cell carcinoma, ovarian serous cystadenocarcinoma, uterine corpus endometrial carcinoma, and uterine carcinosarcoma) using bioinformatics technology.

Materials and methods

We identified up-regulated and down-regulated pseudogenes and built a pseudogene-miRNA-mRNA regulatory network through public datasets to explore their potential functions in carcinogenesis and cancer prognosis.

Results

Among the 63 up-regulated pseudogenes identified, LDHAP5 demonstrated the greatest potential as a candidate pseudogene due to its significant association with poor overall survival in ovarian serous cystadenocarcinoma. KEGG pathway analysis revealed that LDHAP5 showed significant enrichment in MicroRNAs in cancer, Pathway in cancer and PI3K-AKT signaling pathway. Further analysis revealed that EGFR was the potential target mRNA of LDHAP5, which may play an important role in ovarian serous cystadenocarcinoma.

Conclusions

LDHAP5 was associated with the occurrence and prognosis of ovarian serous cystadenocarcinoma, and thus shows potential as a novel therapeutic target against such cancer.

The Zmat2 gene in non-mammalian vertebrates: Organizational simplicity within a divergent locus in fish

ZMAT2 is among the least-studied of mammalian proteins and genes, even though it is the ortholog of Snu23, a protein involved in pre-mRNA splicing in yeast. Here we have used data from genomic and gene expression repositories to examine the Zmat2 gene and locus in 8 terrestrial vertebrates, 10 ray-finned fish, and 1 lobe-finned fish representing > 500 million years of evolutionary diversification. The analyses revealed that vertebrate Zmat2 genes are similar to their mammalian counterparts, as in 16/19 species studied they contain 6 exons, and in 18/19 encode a single conserved protein. However, unlike in mammals, no Zmat2 pseudogenes were identified in these vertebrates, although an expressed Zmat2 paralog was characterized in flycatcher that resembled a DNA copy of a processed and retro-transposed mRNA, and thus could be a proto-pseudogene captured during its evolutionary journey from active to inert. The Zmat2 locus in terrestrial vertebrates, and in spotted gar and coelacanth, also shares additional genes with its mammalian counterparts, including Histidyl-tRNA synthetase (Hars), Hars2, and others, but these are absent from the Zmat2 locus in teleost fish, in which Stem-loop-binding protein 2 (Slbp2) and Lymphocyte cytosolic protein 2a (Lcp2a) are present instead. Taken together, these observations argue that a recognizable Zmat2 was present in the earliest vertebrate ancestors, and postulate that during chromosomal tetraploidization and subsequent re-diploidization during modern teleost evolution, the duplicated Zmat2 gene was retained and the original lost. This study also highlights how information from genomic resources can be leveraged to reveal new biologically significant insights.

Complex Analysis of Retroposed Genes' Contribution to Human Genome, Proteome and Transcriptome

Gene duplication is a major driver of organismal evolution. One of the main mechanisms of gene duplications is retroposition, a process in which mRNA is first transcribed into DNA and then reintegrated into the genome. Most gene retrocopies are depleted of the regulatory regions. Nevertheless, examples of functional retrogenes are rapidly increasing. These functions come from the gain of new spatio-temporal expression patterns, imposed by the content of the genomic sequence surrounding inserted cDNA and/or by selectively advantageous mutations, which may lead to the switch from protein coding to regulatory RNA. As recent studies have shown, these genes may lead to new protein domain formation through fusion with other genes, new regulatory RNAs or other regulatory elements. We utilized existing data from high-throughput technologies to create a complex description of retrogenes functionality. Our analysis led to the identification of human retroposed genes that substantially contributed to transcriptome and proteome. These retrocopies demonstrated the potential to encode proteins or short peptides, act as cis- and trans- Natural Antisense Transcripts (NATs), regulate their progenitors’ expression by competing for the same microRNAs, and provide a sequence to lncRNA and novel exons to existing protein-coding genes. Our study also revealed that retrocopies, similarly to retrotransposons, may act as recombination hot spots. To our best knowledge this is the first complex analysis of these functions of retrocopies.

Zmat2 in mammals: conservation and diversification among genes and Pseudogenes

BACKGROUND

Recent advances in genetics and genomics present unique opportunities for enhancing our understanding of mammalian biology and evolution through detailed multi-species comparative analysis of gene organization and expression. Yet, of the more than 20,000 protein coding genes found in mammalian genomes, fewer than 10% have been examined in any detail. Here we elucidate the power of data available in publicly-accessible genomic and genetic resources by querying them to evaluate Zmat2, a minimally studied gene whose human ortholog has been implicated in spliceosome function and in keratinocyte differentiation.

RESULTS

We find extensive conservation in coding regions and overall structure of Zmat2 in 18 mammals representing 13 orders and spanning ~ 165 million years of evolutionary development, and in their encoded proteins. We identify a tandem duplication in the Zmat2 gene and locus in opossum, but not in other monotremes, marsupials, or other mammals, indicating that this event occurred subsequent to the divergence of these species from one another. We also define a collection of Zmat2 pseudogenes in half of the mammals studied, and suggest based on phylogenetic analysis that they each arose independently in the recent evolutionary past.

CONCLUSIONS

Mammalian Zmat2 genes and ZMAT2 proteins illustrate conservation of structure and sequence, along with the development and diversification of pseudogenes in a large fraction of species. Collectively, these observations also illustrate how the focused identification and interpretation of data found in public genomic and gene expression resources can be leveraged to reveal new insights of potentially high biological significance.

The genomic context of retrocopies increases their chance of functional relevancy in mammals

Described as "junk" DNA, pseudogenes are dead structures of previously active genes present in genomes. Pseudogenes are categorized into two main classes: processed pseudogenes, formed through retrotransposition, and non-processed pseudogenes, typically originated from gene decay following duplication events. The term "processed pseudogene" has changed to "retrocopy" since they are likely to evolve new functional roles and became a retrogene. Here, we surveyed 38,080 retrocopies from chimpanzee, dog, human, mouse, and rat genomes to assess their potential adaptive value. The retrocopies inserted in the same chromosome of the parental gene have higher chances of remain potentially "active" (absence of premature stop codons and frameshifts) (~26.1%), while those placed into a different chromosome have a twofold decrease chance of continuing potentially "active" (~7.52%). The genomic context of their placement seems associated with their expression. Retrocopies placed in intragenic regions and the same sense of the "host" gene have higher chances of being expressed relative to other genomic contexts. The proximity of retrocopies to their parental gene is associated with a lower decay rate, and their location likely influence their expression. Thus, despite their unclear role, retrocopies are probably involved in adaptive processes. Our results evidence natural selection acting in retrocopies.

Re-recognition of pseudogenes: From molecular to clinical applications

Pseudogenes were initially regarded as "nonfunctional" genomic elements that did not have protein-coding abilities due to several endogenous inactivating mutations. Although pseudogenes are widely expressed in prokaryotes and eukaryotes, for decades, they have been largely ignored and classified as gene "junk" or "relics". With the widespread availability of high-throughput sequencing analysis, especially omics technologies, knowledge concerning pseudogenes has substantially increased. Pseudogenes are evolutionarily conserved and derive primarily from a mutation or retrotransposon, conferring the pseudogene with a "gene repository" role to store and expand genetic information. In contrast to previous notions, pseudogenes have a variety of functions at the DNA, RNA and protein levels for broadly participating in gene regulation to influence the development and progression of certain diseases, especially cancer. Indeed, some pseudogenes have been proven to encode proteins, strongly contradicting their "trash" identification, and have been confirmed to have tissue-specific and disease subtype-specific expression, indicating their own value in disease diagnosis. Moreover, pseudogenes have been correlated with the life expectancy of patients and exhibit great potential for future use in disease treatment, suggesting that they are promising biomarkers and therapeutic targets for clinical applications. In this review, we summarize the natural properties, functions, disease involvement and clinical value of pseudogenes. Although our knowledge of pseudogenes remains nascent, this field deserves more attention and deeper exploration.

PseudoFuN: Deriving functional potentials of pseudogenes from integrative relationships with genes and microRNAs across 32 cancers

Background

Long thought “relics” of evolution, not until recently have pseudogenes been of medical interest regarding regulation in cancer. Often, these regulatory roles are a direct by-product of their close sequence homology to protein-coding genes. Novel pseudogene-gene (PGG) functional associations can be identified through the integration of biomedical data, such as sequence homology, functional pathways, gene expression, pseudogene expression, and microRNA expression. However, not all of the information has been integrated, and almost all previous pseudogene studies relied on 1:1 pseudogene–parent gene relationships without leveraging other homologous genes/pseudogenes.

Results

We produce PGG families that expand beyond the current 1:1 paradigm. First, we construct expansive PGG databases by (i) CUDAlign graphics processing unit (GPU) accelerated local alignment of all pseudogenes to gene families (totaling 1.6 billion individual local alignments and >40,000 GPU hours) and (ii) BLAST-based assignment of pseudogenes to gene families. Second, we create an open-source web application (PseudoFuN [Pseudogene Functional Networks]) to search for integrative functional relationships of sequence homology, microRNA expression, gene expression, pseudogene expression, and gene ontology. We produce four “flavors” of CUDAlign-based databases (>462,000,000 PGG pairwise alignments and 133,770 PGG families) that can be queried and downloaded using PseudoFuN. These databases are consistent with previous 1:1 PGG annotation and also are much more powerful including millions of de novo PGG associations. For example, we find multiple known (e.g., miR-20a-PTEN-PTENP1) and novel (e.g., miR-375-SOX15-PPP4R1L) microRNA-gene-pseudogene associations in prostate cancer. PseudoFuN provides a “one stop shop” for identifying and visualizing thousands of potential regulatory relationships related to pseudogenes in The Cancer Genome Atlas cancers.

Conclusions

Thousands of new PGG associations can be explored in the context of microRNA-gene-pseudogene co-expression and differential expression with a simple-to-use online tool by bioinformaticians and oncologists alike.

Emerging Role of USP8, HMGA, and Non-Coding RNAs in Pituitary Tumorigenesis

Two novel molecular mechanisms with a driver role in pituitary tumorigenesis have been recently identified. They are (a) mutations in the Ubiquitin-Specific Protease 8 (USP8) gene in corticotroph tumors and (b) overexpression of the HMGA1 and HMGA2 genes in most of the pituitary tumors. Moreover, deregulated expression of the non-coding RNAs has been very frequently observed in this neoplasia. The aim of this review is to better elucidate the role, the mechanisms, and the possible clinical impact of these novel alterations in the development of pituitary neoplasia.

dreamBase: DNA modification, RNA regulation and protein binding of expressed pseudogenes in human health and disease

Although thousands of pseudogenes have been annotated in the human genome, their transcriptional regulation, expression profiles and functional mechanisms are largely unknown. In this study, we developed dreamBase (http://rna.sysu.edu.cn/dreamBase) to facilitate the investigation of DNA modification, RNA regulation and protein binding of potential expressed pseudogenes from multidimensional high-throughput sequencing data. Based on ∼5500 ChIP-seq and DNase-seq datasets, we identified genome-wide binding profiles of various transcription-associated factors around pseudogene loci. By integrating ∼18 000 RNA-seq data, we analysed the expression profiles of pseudogenes and explored their co-expression patterns with their parent genes in 32 cancers and 31 normal tissues. By combining microRNA binding sites, we demonstrated complex post-transcriptional regulation networks involving 275 microRNAs and 1201 pseudogenes. We generated ceRNA networks to illustrate the crosstalk between pseudogenes and their parent genes through competitive binding of microRNAs. In addition, we studied transcriptome-wide interactions between RNA binding proteins (RBPs) and pseudogenes based on 458 CLIP-seq datasets. In conjunction with epitranscriptome sequencing data, we also mapped 1039 RNA modification sites onto 635 pseudogenes. This database will provide insights into the transcriptional regulation, expression, functions and mechanisms of pseudogenes as well as their roles in biological processes and diseases.

Comprehensive elaboration of database resources utilized in next-generation sequencing-based tumor somatic mutation detection

The rapid evolution of next-generation sequencing (NGS)-based tumor genomic profile detection and the emergence of molecularly targeted therapies have enabled precision oncology. In NGS-based analysis, various types of databases have been developed to perform different functions. However, many problems still exist when using these public databases. Therefore, it is important to better understand the characteristics and limitations of each database and have them complement each other to provide useful clinical evidence for NGS testing. In this review, we elaborate on the important role of databases and their concrete applications in NGS-based somatic mutation detection. We introduce the typically used databases for sequence alignment, variant filtration, and variant interpretation, and compare the differences between the databases with similar functions. Subsequently, we determine the limitations of each database and provide the corresponding solutions. Furthermore, we present an overview diagram to clearly illustrate the database used in the entire NGS-based somatic mutation detection pipeline.

Karyotypic Diversity and Evolutionary Trends in Neotropical Electric Fish of the Genus Gymnotus (Gymnotiformes: Gymnotidae)

Electric fish of the order Gymnotiformes are endemic to the Neotropical region, and their highest diversity is observed in the Amazon region. The family Gymnotidae, which consists of the genera Electrophorus and Gymnotus, is a natural group and is located at the base of the phylogeny of the order. Gymnotus is a widely distributed and specious genus with high karyotypic diversity, especially concerning to the diploid number and the locations of repetitive sequences. Our karyotyping results in five species of the family Gymnotidae (Gymnotus ucamara, Gymnotus cf. stenoleucus, Gymnotus cf. pedanopterus, Gymnotus mamiraua, and Gymnotus carapo "Maranhão") corroborate the proposal of plasticity of the diploid number in this group. Moreover, in this study, we propose that the 5S ribosomal DNA (rDNA) sequences were species-specific markers that act as a potential biogeographical marker for the genus. Besides, the sequence's location, particularly in G. mamiraua from Central Amazon, shows a close relationship with 5S of the Gymnotus species, with 54 chromosomes, from the Paraná-Paraguay basin in the Center-South of Brazil. Considering that the ancestral diploid number for Gymnotidae is 52 chromosomes, we also suggest that the trend in the family is toward a decrease in the chromosome number. However, the carapo clade stands out in this regard, with an increase and a decrease in chromosome number; this pattern may be reinforced with the ecologic behaviors and the geodispersal patterns of this clade.