The Spatio-Temporal Expression Profiles of Silkworm Pseudogenes Provide Valuable Insights into Their Biological Roles.

Background: Pseudogenes are sequences that have lost the ability to transcribe RNA molecules or encode truncated but possibly functional proteins. While they were once considered to be meaningless remnants of evolution, recent researches have shown that pseudogenes play important roles in various biological processes. However, the studies of pseudogenes in the silkworm, an important model organism, are limited and have focused on single or only a few specific genes. Objective: To fill these gaps, we present a systematic genome-wide studies of pseudogenes in the silkworm. Methods: We identified the pseudogenes in the silkworm using the silkworm genome assemblies, transcriptome, protein sequences from silkworm and its related species. Then we used transcriptome datasets from 832 RNA-seq analyses to construct spatio-temporal expression profiles for these pseudogenes. Additionally, we identified tissue-specifically expressed and differentially expressed pseudogenes to further understand their characteristics. Finally, the functional roles of pseudogenes as lncRNAs were systematically analyzed. Results: We identified a total of 4410 pseudogenes, which were grouped into 4 groups, including duplications (DUPs), unitary pseudogenes (Unitary), processed pseudogenes (retropseudogenes, RETs), and fragments (FRAGs). The most of pseudogenes in the domestic silkworm were generated before the divergence of wild and domestic silkworm, however, the domestication may also involve in the accumulation of pseudogenes. These pseudogenes were clearly divided into 2 cluster, a highly expressed and a lowly expressed, and the posterior silk gland was the tissue with the most tissue-specific pseudogenes (199), implying these pseudogenes may be involved in the development and function of silkgland. We identified 3299 lncRNAs in these pseudogenes, and the target genes of these lncRNAs in silkworm pseudogenes were enriched in the egg formation and olfactory function. Conclusions: This study replenishes the genome annotations for silkworm, provide valuable insights into the biological roles of pseudogenes. It will also contribute to our understanding of the complex gene regulatory networks in the silkworm and will potentially have implications for other organisms as well.

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.

Transforming 'Junk' DNA into Cancer Warriors: The Role of Pseudogenes in Hepatocellular Carcinoma.

In the dynamic landscape of hepatocellular carcinoma (HCC) or the liver cancer research, pseudogenes have emerged from the shadows of genetic obscurity to become central figures, significantly influencing the disease molecular development and clinical trajectory. This review explores a transformative shift in perspective, recognizing pseudogenes not as genetic remnants without function, but as critical regulators in the molecular underpinnings of HCC. Engaging in complex interactions such as microRNA sponging, gene expression modulation, and signaling pathway disruptions, pseudogenes orchestrate a part of the molecular complexity driving tumor genesis, progression, and drug resistance in the liver cancer. Their unique expression patterns in hepatoma tissues herald new opportunities for early HCC detection, offering insights into patient prognosis, and identifying novel targets for therapeutic intervention of this disease. Such advancements underscore the importance of pseudogenes in enriching our understanding and management of HCC, paving the way for more effective diagnostic strategies and targeted therapies in the ongoing battle against this challenging malignancy.

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.

The nature and distribution of putative non-functional alleles suggest only two independent events at the origins of Astyanax mexicanus cavefish populations.

BACKGROUND: Several studies suggested that cavefish populations of Astyanax mexicanus settled during the Late Pleistocene. This implies that the cavefish's most conspicuous phenotypic changes, blindness and depigmentation, and more cryptic characters important for cave life, evolved rapidly. RESULTS: Using the published genomes of 47 Astyanax cavefish from la Cueva de El Pachón, El Sótano de la Tinaja, La Cueva Chica and El Sótano de Molino, we searched for putative loss-of-function mutations in previously defined sets of genes, i.e., vision, circadian clock and pigmentation genes. Putative non-functional alleles for four vision genes were identified. Then, we searched genome-wide for putative non-functional alleles in these four cave populations. Among 512 genes with segregating putative non-functional alleles in cavefish that are absent in surface fish, we found an enrichment in visual perception genes. Among cavefish populations, different levels of shared putative non-functional alleles were found. Using a subset of 12 genes for which putative loss-of-function mutations were found, we extend the analysis of shared pseudogenes to 11 cave populations. Using a subset of six genes for which putative loss-of-function mutations were found in the El Sótano del Toro population, where extensive hybridization with surface fish occurs, we found a correlation between the level of eye regression and the amount of putative non-functional alleles. CONCLUSIONS: We confirm that very few putative non-functional alleles are present in a large set of vision genes, in accordance with the recent origin of Astyanax mexicanus cavefish. Furthermore, the genome-wide analysis indicates an enrichment of putative loss-of-function alleles in genes with vision-related GO-terms, suggesting that visual perception may be the function chiefly impacted by gene losses related to the shift from a surface to a cave environment. The geographic distribution of putative loss-of-function alleles newly suggests that cave populations from Sierra de Guatemala and Sierra de El Abra share a common origin, albeit followed by independent evolution for a long period. It also supports that populations from the Micos area have an independent origin. In El Sótano del Toro, the troglomorphic phenotype is maintained despite massive introgression of the surface genome.

The regulatory roles of small nucleolar RNAs within their host locus.

Small nucleolar RNAs (snoRNAs) are a class of conserved noncoding RNAs forming complexes with proteins to catalyse site-specific modifications on ribosomal RNA. Besides this canonical role, several snoRNAs are now known to regulate diverse levels of gene expression. While these functions are carried out in trans by mature snoRNAs, evidence has also been emerging of regulatory roles of snoRNAs in cis, either within their genomic locus or as longer transcription intermediates during their maturation. Herein, we review recent findings that snoRNAs can interact in cis with their intron to regulate the expression of their host gene. We also explore the ever-growing diversity of longer host-derived snoRNA extensions and their functional impact across the transcriptome. Finally, we discuss the role of snoRNA duplications into forging these new layers of snoRNA-mediated regulation, as well as their involvement in the genomic imprinting of their host locus.

Loss of function of metabolic traits in typhoidal Salmonella without apparent genome degradation.

Salmonella enterica serovar Typhi and Paratyphi A are the cause of typhoid and paratyphoid fever in humans, which are systemic life-threatening illnesses. Both serovars are exclusively adapted to the human host, where they can cause life-long persistent infection. A distinct feature of these serovars is the presence of a relatively high number of degraded coding sequences coding for metabolic pathways, most likely a consequence of their adaptation to a single host. As a result of convergent evolution, these serovars shared many of the degraded coding sequences although often affecting different genes in the same metabolic pathway. However, there are several coding sequences that appear intact in one serovar while clearly degraded in the other, suggesting differences in their metabolic capabilities. Here, we examined the functionality of metabolic pathways that appear intact in S. Typhi but that show clear signs of degradation in S. Paratyphi A. We found that, in all cases, the existence of single amino acid substitutions in S. Typhi metabolic enzymes, transporters, or transcription regulators resulted in the inactivation of these metabolic pathways. Thus, the inability of S. Typhi to metabolize Glucose-6-Phosphate or 3-phosphoglyceric acid is due to the silencing of the expression of the genes encoding the transporters for these compounds due to point mutations in the transcriptional regulatory proteins. In contrast, its inability to utilize glucarate or galactarate is due to the presence of point mutations in the transporter and enzymes necessary for the metabolism of these sugars. These studies provide additional support for the concept of adaptive convergent evolution of these two human-adapted S. enterica serovars and highlight a limitation of bioinformatic approaches to predict metabolic capabilities. IMPORTANCE: Salmonella enterica serovar Typhi and Paratyphi A are the cause of typhoid and paratyphoid fever in humans, which are systemic life-threatening illnesses. Both serovars can only infect the human host, where they can cause life-long persistent infection. Because of their adaptation to the human host, these bacterial pathogens have changed their metabolism, leading to the loss of their ability to utilize certain nutrients. In this study we examined the functionality of metabolic pathways that appear intact in S. Typhi but that show clear signs of degradation in S. Paratyphi A. We found that, in all cases, the existence of single amino acid substitutions in S. Typhi metabolic enzymes, transporters, or transcription regulators resulted in the inactivation of these metabolic pathways. These studies provide additional support for the concept of adaptive convergent evolution of these two human-adapted S. enterica serovars.

CpG Islands, Gene Expression and Pseudogenization: A Case for a Potential Trilogy.

BACKGROUND: The promoters of mammalian genes contain clusters of CG dinucleotides known as CpG islands. Most mammalian housekeeping genes predominantly contain CpG islands (CGIs), facilitating gene transcription. Numerous studies have explored the physiological implications of the relationship between CGIs and gene expression. However, the evolutionary implications of this relationship remain largely unexplored. Pseudogenes, in contrast, are genomic remnants that have lost their function over evolutionary time. METHODS: In our current research, we employed comparative genomic techniques to demonstrate a correlation between the absence of gene expression due to a lack of CGIs in the gene promoters and pseudogenization. RESULTS: We showed that there is a significant enrichment of tissue-specific genes in the functional orthologs of pseudogenes. We also found a significant correlation between the lack of CGIs and enriched tissue specificity in these functional orthologs of pseudogenes. CONCLUSIONS: We inferred that perhaps tissue-specific genes are more prone to the process of pseudogenization. In this way, because of their impact on gene expression, CGIs may affect the fate of a gene. To our knowledge, this is the first study to propose a connection between CGIs, gene expression, and the pseudogenization process and discuss the evolutionary implications of this potential trilogy.

Molecular Evidence for Relaxed Selection on the Enamel Genes of Toothed Whales (Odontoceti) with Degenerative Enamel Phenotypes.

Different species of toothed whales (Odontoceti) exhibit a variety of tooth forms and enamel types. Some odontocetes have highly prismatic enamel with Hunter-Schreger bands, whereas enamel is vestigial or entirely lacking in other species. Different tooth forms and enamel types are associated with alternate feeding strategies that range from biting and grasping prey with teeth in most oceanic and river dolphins to the suction feeding of softer prey items without the use of teeth in many beaked whales. At the molecular level, previous studies have documented inactivating mutations in the enamel-specific genes of some odontocete species that lack complex enamel. At a broader scale, however, it is unclear whether enamel complexity across the full diversity of extant Odontoceti correlates with the relative strength of purifying selection on enamel-specific genes. Here, we employ sequence alignments for seven enamel-specific genes (ACP4, AMBN, AMELX, AMTN, ENAM, KLK4, MMP20) in 62 odontocete species that are representative of all extant families. The sequences for 33 odontocete species were obtained from databases, and sequences for the remaining 29 species were newly generated for this study. We screened these alignments for inactivating mutations (e.g., frameshift indels) and provide a comprehensive catalog of these mutations in species with one or more inactivated enamel genes. Inactivating mutations are rare in Delphinidae (oceanic dolphins) and Platanistidae/Inioidea (river dolphins) that have higher enamel complexity scores. By contrast, mutations are much more numerous in clades such as Monodontidae (narwhal, beluga), Ziphiidae (beaked whales), Physeteroidea (sperm whales), and Phocoenidae (porpoises) that are characterized by simpler enamel or even enamelless teeth. Further, several higher-level taxa (e.g., Hyperoodon, Kogiidae, Monodontidae) possess shared inactivating mutations in one or more enamel genes, which suggests loss of function of these genes in the common ancestor of each clade. We also performed selection (dN/dS) analyses on a concatenation of these genes and used linear regression and Spearman's rank-order correlation to test for correlations between enamel complexity and two different measures of selection intensity (# of inactivating mutations per million years, dN/dS values). Selection analyses revealed that relaxed purifying selection is especially prominent in physeteroids, monodontids, and phocoenids. Linear regressions and correlation analyses revealed a strong negative correlation between selective pressure (dN/dS values) and enamel complexity. Stronger purifying selection (low dN/dS) is found on branches with more complex enamel and weaker purifying selection (higher dN/dS) occurs on branches with less complex enamel or enamelless teeth. As odontocetes diversified into a variety of feeding modes, in particular, the suction capture of prey, a reduced reliance on the dentition for prey capture resulted in the relaxed selection of genes that are critical to enamel development.

Stage II oesophageal carcinoma: peril in disguise associated with cellular reprogramming and oncogenesis regulated by pseudogenes.

INTRODUCTION: Pseudogenes have been implicated for their role in regulating cellular differentiation and organismal development. However, their role in promoting cancer-associated differentiation has not been well-studied. This study explores the tumour landscape of oesophageal carcinoma to identify pseudogenes that may regulate events of differentiation to promote oncogenic transformation. MATERIALS AND METHOD: De-regulated differentiation-associated pseudogenes were identified using DeSeq2 followed by 'InteractiVenn' analysis to identify their expression pattern. Gene expression dependent and independent enrichment analyses were performed with GSEA and ShinyGO, respectively, followed by quantification of cellular reprogramming, extent of differentiation and pleiotropy using three unique metrics. Stage-specific gene regulatory networks using Bayesian Network Splitting Average were generated, followed by network topology analysis. MEME, STREME and Tomtom were employed to identify transcription factors and miRNAs that play a regulatory role downstream of pseudogenes to initiate cellular reprogramming and further promote oncogenic transformation. The patient samples were stratified based on the expression pattern of pseudogenes, followed by GSEA, mutation analysis and survival analysis using GSEA, MAF and 'survminer', respectively. RESULTS: Pseudogenes display a unique stage-wise expression pattern that characterizes stage II (SII) ESCA with a high rate of cellular reprogramming, degree of differentiation and pleiotropy. Gene regulatory network and associated topology indicate high robustness, thus validating high pleiotropy observed for SII. Pseudogene-regulated expression of SOX2, FEV, PRRX1 and TFAP2A in SII may modulate cellular reprogramming and promote oncogenesis. Additionally, patient stratification-based mutational analysis in SII signifies APOBEC3A (A3A) as a potential hallmark of homeostatic mutational events of reprogrammed cells which in addition to de-regulated APOBEC3G leads to distinct events of hypermutations. Further enrichment analysis for both cohorts revealed the critical role of combinatorial expression of pseudogenes in cellular reprogramming. Finally, survival analysis reveals distinct genes that promote poor prognosis in SII ESCA and patient-stratified cohorts, thus providing valuable prognostic bio-markers along with markers of differentiation and oncogenesis for distinct landscapes of pseudogene expression. CONCLUSION: Pseudogenes associated with the events of differentiation potentially aid in the initiation of cellular reprogramming to facilitate oncogenic transformation, especially during SII ESCA. Despite a better overall survival of SII, patient stratification reveals combinatorial de-regulation of pseudogenes as a notable marker for a high degree of cellular differentiation with a unique mutational landscape.

The development of a next-generation sequencing panel targeting cannabinoid synthase genes to distinguish between marijuana and hemp.

Hemp and marijuana, both derived from Cannabis sativa L. (C. sativa), are subject to divergent legal regulations due to their different Δ9-tetrahydrocannabinol (Δ9-THC) contents. Cannabinoid synthase genes are considered the key enzymes that determine the chemical composition or chemotype of a particular cultivar. However, existing methods for crop type differentiation based on previous synthase gene theories have limitations in terms of precision and specificity, and a wider range of cannabis varieties must be considered when examining cannabis-based genetic markers. A custom next-generation sequencing (NGS) panel was developed targeting all synthase genes, including Δ9-THC acid synthase, cannabidiolic acid synthase, and cannabichromenic acid synthase, as well as the pseudogenes across diverse C. sativa samples, spanning reference hemp and marijuana, commercial hemp derivatives, and seized marijuana extracts. Interpretation of NGS data revealed a relationship between genotypes and underlying chemotypes, with the principal component analysis indicating a clear distinction between hemp and marijuana clusters. This differentiation was attributed to variations in both synthase genes and pseudogene variants. Finally, this study proposes a genetic cannabis classification method using a differentiation flow chart with novel synthase markers. The flow chart successfully differentiated hemp from marijuana with a 1.3% error rate (n = 147).

Pseudogenomic insights into the evolution of Mycobacterium ulcerans.

BACKGROUND: Buruli ulcer (BU) disease, caused by Mycobacterium ulcerans (MU), and characterized by necrotic ulcers is still a health problem in Africa and Australia. The genome of the bacterium has several pseudogenes due to recent evolutionary events and environmental pressures. Pseudogenes are genetic elements regarded as nonessential in bacteria, however, they are less studied due to limited available tools to provide understanding of their evolution and roles in MU pathogenicity. RESULTS: This study developed a bioinformatic pipeline to profile the pseudogenomes of sequenced MU clinical isolates from different countries. One hundred and seventy-two MU genomes analyzed revealed that pseudogenomes of African strains corresponded to the two African lineages 1 and 2. Pseudogenomes were lineage and location specific and African lineage 1 was further divided into A and B. Lineage 2 had less relaxation in positive selection than lineage 1 which may signify different evolutionary points. Based on the Gil-Latorre model, African MU strains may be in the latter stages of evolutionary adaption and are adapting to an environment rich in metabolic resources with a lower temperature and decreased UV radiation. The environment fosters oxidative metabolism and MU may be less reliant on some secondary metabolites. In-house pseudogenomes from Ghana and Cote d'Ivoire were different from other African strains, however, they were identified as African strains. CONCLUSION: Our bioinformatic pipeline provides pseudogenomic insights to complement other whole genome analyses, providing a better view of the evolution of the genome of MU and suggest an adaptation model which is important in understanding transmission. MU pseudogene profiles vary based on lineage and country, and an apparent reduction in insertion sequences used for the detection of MU which may adversely affect the sensitivity of diagnosis.

SIGNIFICANCE: Prevention and treatment of Buruli ulcer is still a problem but large whole genome datasets on M. ulcerans are readily available. However, genomic studies fail to thoroughly investigate pseudogenes to probe evolutionary changes in the bacteria, and this can be attributed to the lack of bioinformatic tools. This work studied pseudogenes in Mycobacterium ulcerans (MU) to understand its adapted niche and evolutionary differences across African strains. Our results posit an MU niche-adapted model important in understanding transmission. Also, MU pseudogene profiles vary based on lineage and country, suggesting their influence on pseudogenization patterns in the genome. We further identify a reduction in insertion sequences that are used for the detection of the bacteria which may affect the sensitivity of diagnosis.

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.

Numerous insertions of mitochondrial DNA in the genome of the northern mole vole, Ellobius talpinus.

BACKGROUND: Ellobius talpinus is a subterranean rodent representing an attractive model in population ecology studies due to its highly special lifestyle and sociality. In such studies, mitochondrial DNA (mtDNA) is widely used. However, if nuclear copies of mtDNA, aka NUMTs, are present, they may co-amplify with the target mtDNA fragment, generating misleading results. The aim of this study was to determine whether NUMTs are present in E. talpinus. METHODS AND RESULTS: PCR amplification of the putative mtDNA CytB-D-loop fragment using 'universal' primers from 56 E. talpinus samples produced multiple double peaks in 90% of the sequencing chromatograms. To reveal NUMTs, molecular cloning and sequencing of PCR products of three specimens was conducted, followed by phylogenetic analysis. The pseudogene nature of three out of the seven detected haplotypes was confirmed by their basal positions in relation to other Ellobius haplotypes in the phylogenetic tree. Additionally, 'haplotype B' was basal in relation to other E. talpinus haplotypes and found present in very distant sampling sites. BLASTN search revealed 195 NUMTs in the E. talpinus nuclear genome, including fragments of all four PCR amplified pseudogenes. Although the majority of the NUMTs studied were short, the entire mtDNA had copies in the nuclear genome. The most numerous NUMTs were found for rrnL, COXI, and D-loop. CONCLUSIONS: Numerous NUMTs are present in E. talpinus and can be difficult to discriminate against mtDNA sequences. Thus, in future population or phylogenetic studies in E. talpinus, the possibility of cryptic NUMTs amplification should always be taken into account.

Phylogeny, Ecology, and Gene Families Covariation Shaped the Olfactory Subgenome of Rodents.

Olfactory receptor (OR) genes represent the largest multigenic family in mammalian genomes and encode proteins that bind environmental odorant molecules. The OR repertoire is extremely variable among species and is subject to many gene duplications and losses, which have been linked to ecological adaptations in mammals. Although they have been studied on a broad taxonomic scale (i.e., placental), finer sampling has rarely been explored in order to better capture the mechanisms that drove the evolution of the OR repertoire. Among placental mammals, rodents are well-suited for this task, as they exhibit diverse life history traits, and genomic data are available for most major families and a diverse array of lifestyles. In this study, 53 rodent published genomes were mined for their OR subgenomes. We retrieved more than 85,000 functional and pseudogene OR sequences that were subsequently classified into phylogenetic clusters. Copy number variation among rodents is similar to that of other mammals. Using our OR counts along with comparative phylogenetic approaches, we demonstrated that ecological niches such as diet, period of activity, and a fossorial lifestyle strongly impacted the proportion of OR pseudogenes. Within the OR subgenome, phylogenetic inertia was the main factor explaining the relative variations of the 13 OR gene families. However, a striking exception was a convergent 10-fold expansion of the OR family 14 among the phylogenetically divergent subterranean mole-rat lineages belonging to Bathyergidae and Spalacidae families. This study illustrates how the diversity of the OR repertoire has evolved among rodents, both shaped by selective forces stemming from species life history traits and neutral evolution along the rodent phylogeny.

From teeth to pad: tooth loss and development of keratinous structures in sirenians.

Sirenians are a well-known example of morphological adaptation to a shallow-water grazing diet characterized by a modified feeding apparatus and orofacial morphology. Such adaptations were accompanied by an anterior tooth reduction associated with the development of keratinized pads, the evolution of which remains elusive. Among sirenians, the recently extinct Steller's sea cow represents a special case for being completely toothless. Here, we used µ-CT scans of sirenian crania to understand how motor-sensor systems associated with tooth innervation responded to innovations such as keratinized pads and continuous dental replacement. In addition, we surveyed nine genes associated with dental reduction for signatures of loss of function. Our results reveal how patterns of innervation changed with modifications of the dental formula, especially continuous replacement in manatees. Both our morphological and genomic data show that dental development was not completely lost in the edentulous Steller's sea cows. By tracing the phylogenetic history of tooth innervation, we illustrate the role of development in promoting the innervation of keratinized pads, similar to the secondary use of dental canals for innervating neomorphic keratinized structures in other tetrapod groups.

ORFans in Mitochondrial Genomes of Marine Polychaete Polydora.

Most characterized metazoan mitochondrial genomes are compact and encode a small set of proteins that are essential for oxidative phosphorylation, as well as rRNA and tRNA for their expression. However, in rare cases, invertebrate taxa have additional open reading frames (ORFs) in their mtDNA sequences. Here, we sequenced and analyzed the mitochondrial genome of a polychaete worm, Polydora cf. ciliata, part of whose life cycle takes place in low-oxygen conditions. In the mitogenome, we found three "ORFan" regions (544, 1,060, and 427 bp) that have no resemblance to any standard metazoan mtDNA gene but lack stop codons in one of the reading frames. Similar regions are found in the mitochondrial genomes of three other Polydora species and Bocardiella hamata. All five species share the same gene order in their mitogenomes, which differ from that of other known Spionidae mitogenomes. By analyzing the ORFan sequences, we found that they are under purifying selection pressure and contain conservative regions. The codon adaptation indices (CAIs) of the ORFan genes were in the same range of values as the CAI of conventional protein-coding genes in corresponding mitochondrial genomes. The analysis of the P. cf. ciliata mitochondrial transcriptome showed that ORFan-544, ORFan-427, and a portion of the ORFan-1060 are transcribed. Together, this suggests that ORFan-544 and ORFan-427 encode functional proteins. It is likely that the ORFans originated when the Polydora/Bocardiella species complex separated from the rest of the Spionidae, and this event coincided with massive gene rearrangements in their mitochondrial genomes and tRNA-Met duplication.

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.

Duplications of Human Longevity-Associated Genes Across Placental Mammals.

Natural selection has shaped a wide range of lifespans across mammals, with a few long-lived species showing negligible signs of ageing. Approaches used to elucidate the genetic mechanisms underlying mammalian longevity usually involve phylogenetic selection tests on candidate genes, detections of convergent amino acid changes in long-lived lineages, analyses of differential gene expression between age cohorts or species, and measurements of age-related epigenetic changes. However, the link between gene duplication and evolution of mammalian longevity has not been widely investigated. Here, we explored the association between gene duplication and mammalian lifespan by analyzing 287 human longevity-associated genes across 37 placental mammals. We estimated that the expansion rate of these genes is eight times higher than their contraction rate across these 37 species. Using phylogenetic approaches, we identified 43 genes whose duplication levels are significantly correlated with longevity quotients (False Discovery Rate (FDR) < 0.05). In particular, the strong correlation observed for four genes (CREBBP, PIK3R1, HELLS, FOXM1) appears to be driven mainly by their high duplication levels in two ageing extremists, the naked mole rat (Heterocephalus glaber) and the greater mouse-eared bat (Myotis myotis). Further sequence and expression analyses suggest that the gene PIK3R1 may have undergone a convergent duplication event, whereby the similar region of its coding sequence was independently duplicated multiple times in both of these long-lived species. Collectively, this study identified several candidate genes whose duplications may underlie the extreme longevity in mammals, and highlighted the potential role of gene duplication in the evolution of mammalian long lifespans.