Evolution and function of developmentally dynamic pseudogenes in mammals

Single-cell transcriptional footprint for pseudogene SsCLEC9A is associated with antigen processing and presentation in Sus scrofa

The C-type lectin domain family 9 member A (CLEC9A) is widely recognized as the most critical receptor protein for cross presentation of dead cell associated antigens in animal dendritic cells (DCs). Surprisingly, we revealed for the first time that the sole CLEC9A (SsCLEC9A) in pigs becomes a pseudogene due to three causal mutations that occurred approximately 29.8-44.7 million years ago, challenging the significance of CLEC9A in immune cross-presentation across mammals. Interestingly, we found that SsCLEC9A can transcribe a mutated transcript encoding a truncated protein. Through fluorescence-activated cell sorting and single-cell RNA sequencing, we observed that SsCLEC9A mutant transcript is mainly expressed in DCs and correlated with the expression of its homolog CLEC7A. Further data showed that DCs with SsCLEC9A mutant transcripts exhibited reduced cellular interaction ability and downregulation of antigen presentation function, displaying the characteristics of mature DCs. In addition, introducing the conserved sequence of CLEC9A gene into FLT3L-induced bone marrow hematopoietic cells significantly increased the expression of genes involved in antigen processing and presentation. This study presents a natural mutation model of pseudogenes to understand its transcriptional adation, and provides a fundamental basis for rescuing SsCLEC9A to promote immunity in pigs in the future.

Human endogenous retrovirus-H long terminal repeat-associating 2: an emerging immune checkpoint for cancer immunotherapy

Human endogenous retrovirus-H long terminal repeat-associating 2 (HHLA2), a member of the B7 family of co-signaling molecules, is aberrantly expressed in various human cancers and has emerged as a promising target for cancer immunotherapy. It exhibits a unique structure and tissue distribution pattern compared to other B7 family members, where its expression is regulated by the complex physiological and tumor microenvironment. HHLA2 plays a crucial but contradictory role in immune modulation and is thereby associated with heterogeneous prognostic implications across different cancer types. It interacts with two distinct receptors: transmembrane and immunoglobulin domain-containing 2 (TMIGD2), which is predominantly expressed on naïve T and natural killer (NK) cells to deliver co-stimulatory signals to T cells and NK cells, and killer cell immunoglobulin-like receptor, three immunoglobulin domains, and long cytoplasmic tail (KIR3DL3), which is prevalent on terminally differentiated T and CD56dim CD16+ NK cells to transmit inhibitory signals. The expression dynamics of these receptors on immune cells contribute to the maintenance of immune response homeostasis. Therapeutic strategies targeting the HHLA2 immune checkpoint aim to selectively inhibit the immunosuppressive HHLA2-KIR3DL3 pathway while preserving the HHLA2-TMIGD2 signaling. Several anti-HHLA2 and anti-KIR3DL3 antibodies are currently under investigation in early clinical trials, building upon encouraging results observed in humanized mouse models. Notably, the nonoverlapping expression of HHLA2 and PD-L1 in tumors suggests potential synergistic benefits of combining HHLA2-KIR3DL3-targeted therapies with PD-1/PD-L1 blockade or anti-CTLA-4 to augment antitumor activity.

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.

Genetic variation in patent foramen ovale: a case-control genome-wide association study

Background

Patent foramen ovale (PFO) is a congenital defect between the atria, resulting in abnormal hemodynamics. We conducted a genome-wide association study (GWAS) to identify common genetic variants associated with PFO.

Methods

We performed a whole genome sequencing in a discovery cohort of 3,227 unrelated Chinese participants screened for PFO via contrast transthoracic echocardiography (cTTE). Single-nucleotide polymorphisms (SNPs) associated with PFO were further validated by Sanger sequencing and subsequently were evaluated in a validation cohort. Expression quantitative trait loci (eQTL) analysis was conducted using the GTEx database. Single-cell sequencing analyses with pseudotime trajectory modeling were employed to evaluate their expression in human fetal hearts.

Results

The case-control GWAS of discovery cohort ultimately included 517 cases and 517 demographically matched controls. Of the 7,040,407 variants assessed, we identified rs1227675732 (OR = 2.903; 95% CI, 1.961 to 4.297; p = 3.05 × 10-8), rs62206790 (OR = 2.780; 95% CI, 1.864 to 4.146; p = 2.02 × 10-7), rs879176184 (OR = 2.724; 95% CI, 1.822 to 4.073; p = 4.30 × 10-7) and rs13115019 (OR = 2.437; 95% CI, 1.702 to 3.488; p = 5.80 × 10-7) as high-risk variants for PFO, while rs57922961 (OR = 0.5081; 95% CI, 0.388 to 0.666; p = 6.82 × 10-7) was identified as protective variant. These variations were replicated in the validation cohort (111 cases and 152 controls). Single-cell sequencing showed that CNOT2, KCNMB4, MLLT10, IGBP1, and FRG1 were highly expressed with significant changes during heart development.

Conclusion

The identification of susceptible loci for PFO might provide insights into the pathogenesis of PFO and contribute to understanding heart development.

Clinical Trial Registration

https://www.chictr.org.cn/showproj.html?proj=40590, identifier ChiCTR1900024623.

Proteomics Can Rise to the Challenge of Pseudogenes' Coding Nature

Throughout the past decade, technological advances in genomics and transcriptomics have revealed pervasive translation throughout mammalian genomes. These putative proteins are usually excluded from proteomics analyses, as they are absent from common protein repositories. A sizable portion of these noncanonical proteins is translated from pseudogenes. Pseudogenes are commonly termed defective copies of coding genes unable to produce proteins. Here, we suggest that proteomics can help in their annotation. First, we define important terms and review specific examples underlining the caveats in pseudogene annotation and their coding potential. Then, we will discuss the challenges inherent to pseudogenes that have thus far rendered complex their confidence in omics data. Finally, we identify recent developments in experimental procedures, instrumentation, and computational methods in proteomics that put the field in a unique position to solve the pseudogene annotation conundrum.

Zebrafish: unraveling genetic complexity through duplicated genes

The zebrafish is an invaluable model organism for genetic, developmental, and disease research. Although its high conservation with humans is often cited as justification for its use, the zebrafish harbors oft-ignored genetic characteristics that may provide unique insights into gene structure and function. Zebrafish, along with other teleost fish, underwent an additional round of whole genome duplication after their split from tetrapods-resulting in an abundance of duplicated genes when compared to other vertebrates. These duplicated genes have evolved in distinct ways over the ensuing 350 million years. Thus, each gene within a duplicated gene pair has nuanced differences that create a unique identity. By investigating both members of the gene pair together, we can elucidate the mechanisms that underly protein structure and function and drive the complex interplay within biological systems, such as signal transduction cascades, genetic regulatory networks, and evolution of tissue and organ function. It is crucial to leverage such studies to explore these molecular dynamics, which could have far-reaching implications for both basic science and therapeutic development. Here, we will review the role of gene duplications and the existing models for gene divergence and retention following these events. We will also highlight examples within each of these models where studies comparing duplicated genes in the zebrafish have yielded key insights into protein structure, function, and regulation.

Network-based modelling reveals cell-type enriched patterns of non-coding RNA regulation during human skeletal muscle remodelling

A majority of human genes produce non-protein-coding RNA (ncRNA), and some have roles in development and disease. Neither ncRNA nor human skeletal muscle is ideally studied using short-read sequencing, so we used a customised RNA pipeline and network modelling to study cell-type specific ncRNA responses during muscle growth at scale. We completed five human resistance-training studies (n=144 subjects), identifying 61% who successfully accrued muscle-mass. We produced 288 transcriptome-wide profiles and found 110 ncRNAs linked to muscle growth in vivo, while a transcriptome-driven network model demonstrated interactions via a number of discrete functional pathways and single-cell types. This analysis included established hypertrophy-related ncRNAs, including CYTOR - which was leukocyte-associated (FDR = 4.9 ×10-7). Novel hypertrophy-linked ncRNAs included PPP1CB-DT (myofibril assembly genes, FDR = 8.15 × 10-8), and EEF1A1P24 and TMSB4XP8 (vascular remodelling and angiogenesis genes, FDR = 2.77 × 10-5). We also discovered that hypertrophy lncRNA MYREM shows a specific myonuclear expression pattern in vivo. Our multi-layered analyses established that single-cell-associated ncRNA are identifiable from bulk muscle transcriptomic data and that hypertrophy-linked ncRNA genes mediate their association with muscle growth via multiple cell types and a set of interacting pathways.

Network-based modelling reveals cell-type enriched patterns of non-coding RNA regulation during human skeletal muscle remodelling

A majority of human genes produce non-protein-coding RNA (ncRNA), and some have roles in development and disease. Neither ncRNA nor human skeletal muscle is ideally studied using short-read sequencing, so we used a customized RNA pipeline and network modelling to study cell-type specific ncRNA responses during muscle growth at scale. We completed five human resistance-training studies (n = 144 subjects), identifying 61% who successfully accrued muscle-mass. We produced 288 transcriptome-wide profiles and found 110 ncRNAs linked to muscle growth in vivo, while a transcriptome-driven network model demonstrated interactions via a number of discrete functional pathways and single-cell types. This analysis included established hypertrophy-related ncRNAs, including CYTOR-which was leukocyte-associated (false discovery rate [FDR] = 4.9 × 10-7). Novel hypertrophy-linked ncRNAs included PPP1CB-DT (myofibril assembly genes, FDR = 8.15 × 10-8), and EEF1A1P24 and TMSB4XP8 (vascular remodelling and angiogenesis genes, FDR = 2.77 × 10-5). We also discovered that hypertrophy lncRNA MYREM shows a specific myonuclear expression pattern in vivo. Our multi-layered analyses established that single-cell-associated ncRNA are identifiable from bulk muscle transcriptomic data and that hypertrophy-linked ncRNA genes mediate their association with muscle growth via multiple cell types and a set of interacting pathways.

Novel Candidate Genes Involved in an Initial Stage of White Striping Development in Broiler Chickens

White striping (WS) is a myopathy characterized by the appearance of white stripes parallel to the muscle fibers in the breast of broiler chickens, composed of adipose and connective tissues. This condition causes economic losses and, although common, its etiology remains poorly understood. Hence, the objective was to identify genes and biological mechanisms involved in the early stages of WS using a paternal broiler line that grows slightly slower than commercial ones, at 35 days of age, through the RNA sequencing of the pectoralis major muscle. Thirty genes were differentially expressed between normal and WS-affected chickens, with 23 upregulated and 7 downregulated in the affected broilers. Of these, 14 genes are novel candidates for WS and are implicated in biological processes related to muscle development (CEPBD, DUSP8, METTL21EP, NELL2, and UBE3D), lipid metabolism (PDK4, DDIT4, FKBP5, DGAT2, LIPG, TDH, and RGCC), and collagen (COL4A5 and COL4A6). Genes related to changes in muscle fiber type and the processes of apoptosis, autophagy, proliferation, and differentiation are possibly involved with the initial stage of WS development. In contrast, the genes linked to lipid metabolism and collagen may have their expression altered due to the progression of the myopathy.

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.

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

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

Many purported pseudogenes in bacterial genomes are bona fide genes

BACKGROUND

Microbial genomes are largely comprised of protein coding sequences, yet some genomes contain many pseudogenes caused by frameshifts or internal stop codons. These pseudogenes are believed to result from gene degradation during evolution but could also be technical artifacts of genome sequencing or assembly.

RESULTS

Using a combination of observational and experimental data, we show that many putative pseudogenes are attributable to errors that are incorporated into genomes during assembly. Within 126,564 publicly available genomes, we observed that nearly identical genomes often substantially differed in pseudogene counts. Causal inference implicated assembler, sequencing platform, and coverage as likely causative factors. Reassembly of genomes from raw reads confirmed that each variable affects the number of putative pseudogenes in an assembly. Furthermore, simulated sequencing reads corroborated our observations that the quality and quantity of raw data can significantly impact the number of pseudogenes in an assembler dependent fashion. The number of unexpected pseudogenes due to internal stops was highly correlated (R2 = 0.96) with average nucleotide identity to the ground truth genome, implying relative pseudogene counts can be used as a proxy for overall assembly correctness. Applying our method to assemblies in RefSeq resulted in rejection of 3.6% of assemblies due to significantly elevated pseudogene counts. Reassembly from real reads obtained from high coverage genomes showed considerable variability in spurious pseudogenes beyond that observed with simulated reads, reinforcing the finding that high coverage is necessary to mitigate assembly errors.

CONCLUSIONS

Collectively, these results demonstrate that many pseudogenes in microbial genome assemblies are actually genes. Our results suggest that high read coverage is required for correct assembly and indicate an inflated number of pseudogenes due to internal stops is indicative of poor overall assembly quality.

EndoQuad: a comprehensive genome-wide experimentally validated endogenous G-quadruplex database

G-quadruplexes (G4s) are non-canonical four-stranded structures and are emerging as novel genetic regulatory elements. However, a comprehensive genomic annotation of endogenous G4s (eG4s) and systematic characterization of their regulatory network are still lacking, posing major challenges for eG4 research. Here, we present EndoQuad (https://EndoQuad.chenzxlab.cn/) to address these pressing issues by integrating high-throughput experimental data. First, based on high-quality genome-wide eG4s mapping datasets (human: 1181; mouse: 24; chicken: 2) generated by G4 ChIP-seq/CUT&Tag, we generate a reference set of genome-wide eG4s. Our multi-omics analyses show that most eG4s are identified in one or a few cell types. The eG4s with higher occurrences across samples are more structurally stable, evolutionarily conserved, enriched in promoter regions, mark highly expressed genes and associate with complex regulatory programs, demonstrating higher confidence level for further experiments. Finally, we integrate millions of functional genomic variants and prioritize eG4s with regulatory functions in disease and cancer contexts. These efforts have culminated in the comprehensive and interactive database of experimentally validated DNA eG4s. As such, EndoQuad enables users to easily access, download and repurpose these data for their own research. EndoQuad will become a one-stop resource for eG4 research and lay the foundation for future functional studies.

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.

Study of the cross-talk between Fasciola hepatica juveniles and the intestinal epithelial cells of the host by transcriptomics in an in vitro model

Fasciolosis is a globally widespread trematodiasis with a major economic and veterinary impact. Therefore, this disease is responsible for millions of dollars in losses to the livestock industry, and also constitutes an emerging human health problem in endemic areas. The ubiquitous nature of Fasciola hepatica, the main causative agent, is one of the key factors for the success of fasciolosis. Accordingly, this parasite is able to subsist in a wide variety of ecosystems and hosts, thanks to the development of a plethora of strategies for adaption and immune evasion. Fasciolosis comprises a growing concern due to its high prevalence rates, together with the emergence of strains of the parasite resistant to the treatment of choice (triclabendazole). These facts highlight the importance of developing novel control measures which allow for an effective protection against the disease before F. hepatica settles in a niche inaccessible to the immune system. However, knowledge about the initial phases of the infection, including the migration mechanisms of the parasite and the early innate host response, is still scarce. Recently, our group developed an in vitro host-parasite interaction model that allowed the early events to be unveiled after the first contact between the both actors. This occurs shortly upon ingestion of F. hepatica metacercariae and the emergence of the newly excysted juveniles (FhNEJ) in the host duodenum. Here, we present a transcriptomic analysis of such model using an approach based on RNA sequencing (RNA-Seq), which reveals changes in gene expression related to proteolysis and uptake of metabolites in FhNEJ. Additionally, contact with the parasite triggered changes in host intestinal cells related to pseudogenes expression and host defence mechanisms, including immune response, among others. In sum, these results provide a better understanding of the early stages of fasciolosis at molecular level, and a pool of targets that could be used in future therapeutic strategies against the disease.

Integrating massive RNA-seq data to elucidate transcriptome dynamics in Drosophila melanogaster

The volume of ribonucleic acid (RNA)-seq data has increased exponentially, providing numerous new insights into various biological processes. However, due to significant practical challenges, such as data heterogeneity, it is still difficult to ensure the quality of these data when integrated. Although some quality control methods have been developed, sample consistency is rarely considered and these methods are susceptible to artificial factors. Here, we developed MassiveQC, an unsupervised machine learning-based approach, to automatically download and filter large-scale high-throughput data. In addition to the read quality used in other tools, MassiveQC also uses the alignment and expression quality as model features. Meanwhile, it is user-friendly since the cutoff is generated from self-reporting and is applicable to multimodal data. To explore its value, we applied MassiveQC to Drosophila RNA-seq data and generated a comprehensive transcriptome atlas across 28 tissues from embryogenesis to adulthood. We systematically characterized fly gene expression dynamics and found that genes with high expression dynamics were likely to be evolutionarily young and expressed at late developmental stages, exhibiting high nonsynonymous substitution rates and low phenotypic severity, and they were involved in simple regulatory programs. We also discovered that human and Drosophila had strong positive correlations in gene expression in orthologous organs, revealing the great potential of the Drosophila system for studying human development and disease.

Evolution and function of developmentally dynamic pseudogenes in mammals

BACKGROUND

Pseudogenes are excellent markers for genome evolution, which are emerging as crucial regulators of development and disease, especially cancer. However, systematic functional characterization and evolution of pseudogenes remain largely unexplored.

RESULTS

To systematically characterize pseudogenes, we date the origin of human and mouse pseudogenes across vertebrates and observe a burst of pseudogene gain in these two lineages. Based on a hybrid sequencing dataset combining full-length PacBio sequencing, sample-matched Illumina sequencing, and public time-course transcriptome data, we observe that abundant mammalian pseudogenes could be transcribed, which contribute to the establishment of organ identity. Our analyses reveal that developmentally dynamic pseudogenes are evolutionarily conserved and show an increasing weight during development. Besides, they are involved in complex transcriptional and post-transcriptional modulation, exhibiting the signatures of functional enrichment. Coding potential evaluation suggests that 19% of human pseudogenes could be translated, thus serving as a new way for protein innovation. Moreover, pseudogenes carry disease-associated SNPs and conduce to cancer transcriptome perturbation.

CONCLUSIONS

Our discovery reveals an unexpectedly high abundance of mammalian pseudogenes that can be transcribed and translated, and these pseudogenes represent a novel regulatory layer. Our study also prioritizes developmentally dynamic pseudogenes with signatures of functional enrichment and provides a hybrid sequencing dataset for further unraveling their biological mechanisms in organ development and carcinogenesis in the future.