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

Pseudogenes and host specialization in the emergent bacterial plant pathogen Xylella fastidiosa

Pseudogenes are regarded as "junk" DNA, representing vestigial functions no longer needed for fitness. Accordingly, a higher number of pseudogenes in a bacterial human pathogen was previously hypothesized to be a hallmark of host specialists. In this study, we tested this hypothesis on the emergent bacterial plant pathogen Xylella fastidiosa (Xf) to link pseudogene makeup and host specificity. Xf is an ideal subject for these studies by being a xylem-limited pathogen that underwent extensive genome reduction. Using natural host range data of 151 strains and Pseudofinder analysis on Xf whole genome sequences, we observed that Xf subsp. sandyi had the highest pseudogene content, followed by subsp. morus, while subsp. pauca, fastidiosa, and multiplex had the lowest. The first two subspecies are known to have a limited host range compared to the others, aligning with the hypothesis of a greater number of pseudogenes corresponding to narrower host range. Weed isolates are presumably host specialists because they had the highest pseudogene content. Using a thorough pseudogene map across genomes and empirical pathogenicity data on blueberries, we screened for genes potentially involved in blueberry specialization. Targets were identified by selecting sequences pseudogenized (i) in strains infecting hosts different from blueberry and (ii) only in blueberry asymptomatic strains. Six sequences were identified with a potential role in blueberry infection, including one that was common between the two criteria. Here, we generated hypotheses on host range and specificity of Xf strains that need to be tested experimentally to help understand this devastating plant pathogen.IMPORTANCEXylella fastidiosa is a highly destructive plant pathogen that infects hundreds of landscape and agriculturally important plant species mainly in Europe and the Americas. Nevertheless, the host range of specific genotypes and underlying mechanisms of host specificity remain unclear. These are important aspects to determine the potential risk of infection in specific areas depending on the genetic makeup of the pathogen population and hosts present. This study offers valuable insights into the role of pseudogenization in the genomes of different X. fastidiosa strains, linking it to host specialization. Despite the limited information available for the host range of different strains of this pathogen, this research proposes a relationship between the abundance of pseudogenes and host specificity. These findings are essential for predicting potential host shifts by this pathogen, aiding in the development of strategies to prevent its spread. Additionally, the identification of candidate genes putatively important for symptom development in blueberries offers targets for prevention and control efforts.

Retention and pseudogenization of aquaporin-10 in Rodentia

Vertebrates exhibit diversity in the presence and number of aquaporin (Aqp)-10 genes. In Rodentia, mice possess an Aqp10 pseudogene, whereas guinea pigs possess an intact Aqp10. However, Aqp10 retention and pseudogenization history in various rodent lineages remains unclear. Therefore, in this study, we aimed to investigate the molecular evolution of Aqp10 using the recent increasingly decoded rodent genome sequences. We analyzed Aqp10 in the genomes of 43 rodent species belonging to 14 families and found that Aqp10 was pseudogenized in 13 species of three families in the Myomorpha suborder. In contrast, a single intact Aqp10 was retained in the other 30 rodent species, with no Aqp10 pseudogene found in the Castorimorpha, Hystricomorpha, and Sciuromorpha suborders. Additionally, we investigated the tissue expression levels of aquaglyceroporin genes in guinea pigs and rats via reverse transcription-polymerase chain reaction and detected Aqp10 expression in the guinea pig intestines. Notably, none of the examined rat organs expressed Aqp10; however, Aqp7 was expressed in the rat intestines. In situ hybridization showed that guinea pig Aqp10 was expressed in the intestinal epithelial cells. Moreover, AQP10 was permeable to water, glycerol, urea, and boric acid in Xenopus oocytes. Overall, these results clarify the Aqp10 pseudogenization history in Rodentia and suggest guinea pigs as excellent small animal models to analyze the intestinal AQP10 functions.

Sex-stratified genome-wide meta-analysis identifies novel loci for cognitive decline in older adults

INTRODUCTION

Many complex traits and diseases show sex-specific biases in clinical presentation and prevalence.

METHODS

To understand sex-specific genetic architecture of cognitive decline across five cognitive domains (attention, memory, executive function, language, and visuospatial function) and global cognitive function, we performed sex-stratified genome-wide meta-analysis in 3021 older adults aged ≥ 65 years (female = 1545, male = 1476) from three prospective cohorts. Gene-based and gene-set enrichment analyses were conducted for each cognitive trait.

RESULTS

We identified a novel genome-wide significant (GWS) intergenic locus for decline of memory in males near RPS23P3 on chromosome 4 (rs6851574: minor allele frequency [MAF] = 0.39, Pmale = 4.10E-08, βmale = 0.19; Pinteraction = 3.76E-04). We also identified a subthreshold GWS locus for decline of executive function on chromosome 12 in females near the NDUFA12 gene, involved in oxidative phosphorylation (rs11107823: MAF = 0.12, Pfemale  = 9.35E-08, βfemale = 0.28; Pinteraction = 7.42E-06).

DISCUSSION

Sex-aware genetic studies can help in the identification of novel genetic loci and enhance sex-specific understanding of cognitive decline.

HIGHLIGHTS

Our genome-wide meta-analysis of single variants identified two new genetic associations, one in males and one in females. The novel male association was observed with the decline of memory in the intergenic region near the RPS23P3 gene on chromosome 4. This intergenic region has previously been implicated in several brain and cognition related traits, including anatomical brain aging, brain shape, and educational attainment. The novel female-specific association was observed with decline in executive function on chromosome 12 near the NDUFA12 gene, which is involved in oxidative phosphorylation. Sex-stratified analyses offer sex-specific understanding of biological pathways involved in cognitive aging.

Deciphering the Genetic Architecture of Parkinson's Disease in India

The genomic landscape of the Indian population, particularly for age-related disorders like Parkinson's disease (PD) remains underrepresented in global research. Genetic variability in PD has been studied predominantly in European populations, offering limited insights into its role within the Indian population. To address this gap, we conducted the first pan-India genomic survey of PD involving 4,806 cases and 6,364 controls, complemented by a meta-analysis integrating summary statistics from a multi-ancestry PD meta-analysis (N=611,485). We further leveraged RNA-sequencing data from lymphoblastoid cell lines of 731 individuals from the 1000 Genomes project to evaluate the expression of key loci across global populations. Our findings reveal a higher genetic burden of PD in the Indian population compared to Europeans, accounting for ∼30% of the previously unexplained heritability. Thirteen genome-wide significant loci were identified, including two novel loci, with an additional three loci uncovered through meta-analysis. Polygenic risk score analysis showed moderate transferability from European populations. Our results highlight the importance of genetic loci in immune function, lipid metabolism and SNCA aggregation in PD pathogenesis, with gene expression variability emphasizing population-specific differences. We also established South Asia's largest PD biobank, providing a foundation for patient-centric approaches to PD research and treatment in India.

Advancing the Indian cattle pangenome: characterizing non-reference sequences in Bos indicus

BACKGROUND

India harbors the world's largest cattle population, encompassing over 50 distinct Bos indicus breeds. This rich genetic diversity underscores the inadequacy of a single reference genome to fully capture the genomic landscape of Indian cattle. To comprehensively characterize the genomic variation within Bos indicus and, specifically, dairy breeds, we aim to identify non-reference sequences and construct a comprehensive pangenome.

RESULTS

Five representative genomes of prominent dairy breeds, including Gir, Kankrej, Tharparkar, Sahiwal, and Red Sindhi, were sequenced using 10X Genomics 'linked-read' technology. Assemblies generated from these linked-reads ranged from 2.70 Gb to 2.77 Gb, comparable to the Bos indicus Brahman reference genome. A pangenome of Bos indicus cattle was constructed by comparing the newly assembled genomes with the reference using alignment and graph-based methods, revealing 8 Mb and 17.7 Mb of novel sequence respectively. A confident set of 6,844 Non-reference Unique Insertions (NUIs) spanning 7.57 Mb was identified through both methods, representing the pangenome of Indian Bos indicus breeds. Comparative analysis with previously published pangenomes unveiled 2.8 Mb (37%) commonality with the Chinese indicine pangenome and only 1% commonality with the Bos taurus pangenome. Among these, 2,312 NUIs encompassing ~ 2 Mb, were commonly found in 98 samples of the 5 breeds and designated as Bos indicus Common Insertions (BICIs) in the population. Furthermore, 926 BICIs were identified within 682 protein-coding genes, 54 long non-coding RNAs (lncRNA), and 18 pseudogenes. These protein-coding genes were enriched for functions such as chemical synaptic transmission, cell junction organization, cell-cell adhesion, and cell morphogenesis. The protein-coding genes were found in various prominent quantitative trait locus (QTL) regions, suggesting potential roles of BICIs in traits related to milk production, reproduction, exterior, health, meat, and carcass. Notably, 63.21% of the bases within the BICIs call set contained interspersed repeats, predominantly Long Interspersed Nuclear Elements (LINEs). Additionally, 70.28% of BICIs are shared with other domesticated and wild species, highlighting their evolutionary significance.

CONCLUSIONS

This is the first report unveiling a robust set of NUIs defining the pangenome of Bos indicus breeds of India. The analyses contribute valuable insights into the genomic landscape of desi cattle breeds.

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.

Enhancing chemotherapy response prediction via matched colorectal tumor-organoid gene expression analysis and network-based biomarker selection

BACKGROUND

Colorectal cancer (CRC) presents significant challenges in chemotherapy response prediction due to its molecular heterogeneity. Current methods often fail to account for the complexity and variability inherent in individual tumors.

METHODS

We developed a novel approach using matched CRC tumor and organoid gene expression data. We applied Consensus Weighted Gene Co-expression Network Analysis (WGCNA) across three datasets: CRC tumors, matched organoids, and an independent organoid dataset with IC50 drug response values, to identify key gene modules and hub genes linked to chemotherapy response, particularly 5-fluorouracil (5-FU).

FINDINGS

Our integrative analysis identified significant gene modules and hub genes associated with CRC chemotherapy response. The predictive model built from these findings demonstrated superior accuracy over traditional methods when tested on independent datasets. The matched tumor-organoid data approach proved effective in capturing relevant biomarkers, enhancing prediction reliability.

INTERPRETATION

This study provides a robust framework for improving CRC chemotherapy response predictions by leveraging matched tumor and organoid gene expression data. Our approach addresses the limitations of previous methods, offering a promising strategy for personalized treatment planning in CRC. Future research should aim to validate these findings and explore the integration of more comprehensive drug response data.

FUNDING

This research was supported by US National Cancer Institute grant R37CA248289, and Sylvester Comprehensive Cancer Center. which receives funding from the National Cancer Institute award P30CA240139. This work was supported by National Institutes of Health (NIH) under the following grants: T32CA009501-31A1 and R37CA248289. This work was also supported by the MSK P30CA008748 grant.

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.

Polymorphic pseudogenes in the human genome - a comprehensive assessment

BACKGROUND

Over the past decade, variations of the coding portion of the human genome have become increasingly evident. In this study, we focus on polymorphic pseudogenes, a unique and relatively unexplored type of pseudogene whose inactivating mutations have not yet been fixed in the human genome at the global population level. Thus, polymorphic pseudogenes are characterized by the presence in the population of both coding alleles and non-coding alleles originating from Loss-of-Function (LoF) mutations. These alleles can be found both in heterozygosity and in homozygosity in different human populations and thus represent pseudogenes that have not yet been fixed in the population.

RESULTS

A methodical cross-population analysis of 232 polymorphic pseudogenes, including 35 new examples, reveals that human olfactory signalling, drug metabolism and immunity are among the systems most impacted by the variable presence of LoF variants at high frequencies. Within this dataset, a total of 179 genes presented polymorphic LoF variants in all analysed populations. Transcriptome and proteome analysis confirmed that although these genes may harbour LoF alleles, when the coding allele is present, the gene remains active and can play a functional role in various metabolic pathways, including drug/xenobiotic metabolism and immunity. The observation that many polymorphic pseudogenes are members of multigene families argues that genetic redundancy may play a key role in compensating for the inactivation of one paralogue.

CONCLUSIONS

The distribution, expression and integration of cellular/biological networks in relation to human polymorphic pseudogenes, provide novel insights into the architecture of the human genome and the dynamics of gene gain and loss with likely functional impact.

Androgen Receptor and Non-Coding RNAs' Interaction in Renal Cell Carcinoma

Renal cell carcinoma (RCC), the most prevalent among the urogenital cancers, accounts for around 3% of new cancer cases worldwide. Significantly, the incidence of RCC has doubled in developed world countries, ranking it as the sixth most common cancer in males, who represent two-thirds of RCC cases. Males with RCC exhibit a higher mortality rate and tend to develop a more aggressive form of the disease than females. Sex-related risk factors, including lifestyle and biological variations, explain this difference. The androgen receptor (AR) oncogenic signaling pathway has been extensively studied among the biological factors that affect RCC. Recent advancements in high-throughput RNA sequencing techniques have underscored the significant roles played by noncoding-RNAs (ncRNAs), previously dismissed as "junk". The oncogenic potential of AR is manifested through its dysregulation of the ncRNAs' availability and function, promoting RCC tumorigenesis. This review offers a summary of the most recent findings on the role and molecular mechanisms of the AR in dysregulating the ncRNAs that play a role in the progression of RCC and the possibility of utilizing ncRNAs to target AR as a potential therapeutic strategy.

Oxidative stress-induced gene expression changes in prostate epithelial cells in vitro reveal a robust signature of normal prostatic senescence and aging

Oxidative stress has long been postulated to play an essential role in aging mechanisms, and numerous forms of molecular damage associated with oxidative stress have been well documented. However, the extent to which changes in gene expression in direct response to oxidative stress are related to actual cellular aging, senescence, and age-related functional decline remains unclear. Here, we ask whether H2O2-induced oxidative stress and resulting gene expression alterations in prostate epithelial cells in vitro reveal gene regulatory changes typically observed in naturally aging prostate tissue and age-related prostate disease. While a broad range of significant changes observed in the expression of non-coding transcripts implicated in senescence-related responses, we also note an overrepresentation of gene-splicing events among differentially expressed protein-coding genes induced by H2O2. Additionally, the collective expression of these H2O2-induced DEGs is linked to age-related pathological dysfunction, with their protein products exhibiting a dense network of protein-protein interactions. In contrast, co-expression analysis of available gene expression data reveals a naturally occurring highly coordinated expression of H2O2-induced DEGs in normally aging prostate tissue. Furthermore, we find that oxidative stress-induced DEGs statistically overrepresent well-known senescence-related signatures. Our results show that oxidative stress-induced gene expression in prostate epithelial cells in vitro reveals gene regulatory changes typically observed in naturally aging prostate tissue and age-related prostate disease.

Pseudogene GSTM3P1 derived long non-coding RNA promotes ischemic acute kidney injury by target directed microRNA degradation of kidney-protective mir-668

Long non-coding RNAs (lncRNAs) are a group of epigenetic regulators that have been implicated in kidney diseases including acute kidney injury (AKI). However, very little is known about the specific lncRNAs involved in AKI and the mechanisms underlying their pathologic roles. Here, we report a new lncRNA derived from the pseudogene GSTM3P1, which mediates ischemic AKI by interacting with and promoting the degradation of mir-668, a kidney-protective microRNA. GSTM3P1 and its mouse orthologue Gstm2-ps1 were induced by hypoxia in cultured kidney proximal tubular cells. In mouse kidneys, Gstm2-ps1 was significantly upregulated in proximal tubules at an early stage of ischemic AKI. This transient induction of Gstm2-ps1 depends on G3BP1, a key component in stress granules. GSTM3P1 overexpression increased kidney proximal tubular apoptosis after ATP depletion, which was rescued by mir-668. Notably, kidney proximal tubule-specific knockout of Gstm2-ps1 protected mice from ischemic AKI, as evidenced by improved kidney function, diminished tubular damage and apoptosis, and reduced kidney injury biomarker (NGAL) induction. To test the therapeutic potential, Gstm2-ps1 siRNAs were introduced into cultured mouse proximal tubular cells or administered to mice. In cultured cells, Gstm2-ps1 knockdown suppressed ATP depletion-associated apoptosis. In mice, Gstm2-ps1 knockdown ameliorated ischemic AKI. Mechanistically, both GSTM3P1 and Gstm2-ps1 possessed mir-668 binding sites and downregulated the mature form of mir-668. Specifically, GSTM3P1 directly bound to mature mir-668 to induce its decay via target-directed microRNA degradation. Thus, our results identify GSTM3P1 as a novel lncRNA that promotes kidney tubular cell death in AKI by binding mir-668 to inducing its degradation.

Discovering the hidden function in fungal genomes

New molecular technologies have helped unveil previously unexplored facets of the genome beyond the canonical proteome, including microproteins and short ORFs, products of alternative splicing, regulatory non-coding RNAs, as well as transposable elements, cis-regulatory DNA, and other highly repetitive regions of DNA. In this Review, we highlight what is known about this 'hidden genome' within the fungal kingdom. Using well-established model systems as a contextual framework, we describe key elements of this hidden genome in diverse fungal species, and explore how these factors perform critical functions in regulating fungal metabolism, stress tolerance, and pathogenesis. Finally, we discuss new technologies that may be adapted to further characterize the hidden genome in fungi.

The Role of Long Non-Coding RNAs in Ovarian Cancer Cells

Among the most deadly malignancies that strike women worldwide, ovarian cancer is still one of the most common. The primary factor affecting a patient's survival is early lesion discovery. Unfortunately, because ovarian cancer is a sneaky illness that usually manifests as nonspecific symptoms only in advanced stages, its early detection and screening are challenging. A lot of research is being conducted on effective methods of diagnosing and treating ovarian cancer. Recently, non-coding RNAs (ncRNAs) have gained great popularity, which are considered to be the main regulators of many cellular processes, especially those occurring in cancer. LncRNAs are also being studied for their therapeutic use in the treatment of ovarian cancer and their use in diagnostics and as indicators of poor prognosis. In this article, we reviewed lncRNAs described in the literature that may play an important role in ovarian cancer.

Radiotherapy and breast cancer: finally, an lncRNA perspective on radiosensitivity and radioresistance

Radiotherapy (RT) serves as one of the key adjuvant treatments in management of breast cancer. Nevertheless, RT has two major problems: side effects and radioresistance. Given that patients respond differently to RT, it is imperative to understand the molecular mechanisms underlying these differences. Two-thirds of human genes do not encode proteins, as we have realized from genome-scale studies conducted after the advent of the genomic era; nevertheless, molecular understanding of breast cancer to date has been attained almost entirely based on protein-coding genes and their pathways. Long non-coding RNAs (lncRNAs) are a poorly understood but abundant class of human genes that yield functional non-protein-coding RNA transcripts. Here, we canvass the field to seek evidence for the hypothesis that lncRNAs contribute to radioresistance in breast cancer. RT-responsive lncRNAs ranging from "classical" lncRNAs discovered at the dawn of the post-genomic era (such as HOTAIR, NEAT1, and CCAT), to long intergenic lncRNAs such as LINC00511 and LINC02582, antisense lncRNAs such as AFAP-AS1 and FGD5-AS1, and pseudogene transcripts such as DUXAP8 were found during our screen of the literature. Radiation-related pathways modulated by these lncRNAs include DNA damage repair, cell cycle, cancer stem cells phenotype and apoptosis. Thus, providing a clear picture of these lncRNAs' underlying RT-relevant molecular mechanisms should help improve overall survival and optimize the best radiation dose for each individual patient. Moreover, in healthy humans, lncRNAs show greater natural expression variation than protein-coding genes, even across individuals, alluding to their exceptional potential for targeting in truly personalized, precision medicine.

The primate-specific presence of interferon regulatory factor-5 pseudogene 1

Interferon regulatory factor 5 (IRF5) is a key transcription factor in inflammatory and immune responses, with its dysregulation linked to autoimmune diseases. Using bioinformatic approaches, including Basic Local Alignment Search Tool (BLAST) for sequence similarity searches, BLAST-Like Alignment Tool (BLAT) for genome-wide alignments, and several phylogenetics software, such as Multiple Alignment using Fast Fourier Transform (MAFFT), for phylogenetic analyses, we characterized the structure, origin, and evolutionary history of the human IRF5 pseudogene 1 (IRF5P1). Our analyses reveal that IRF5P1 is a chimeric processed pseudogene containing sequences derived from multiple sources, including IRF5-like sequences from disparate organisms. We find that IRF5P1 is specific to higher primates, likely originating through an ancient retroviral integration event approximately 60 million years ago. Interestingly, IRF5P1 resides within the triple QxxK/R motif-containing (TRIQK) gene, and its antisense strand is predominantly expressed as part of the TRIQK pre-messenger RNA (mRNA). Analysis of publicly available RNA-seq data suggests potential expression of antisense IRF5P1 RNA. We hypothesize that this antisense RNA may regulate IRF5 expression through complementary binding to IRF5 mRNA, with human genetic variants potentially modulating this interaction. The conservation of IRF5P1 in the primate lineage suggests its positive effects on primate evolution and innate immunity. This study highlights the importance of investigating pseudogenes and their potential regulatory roles in shaping lineage-specific immune adaptations.

A Pilot Detection and Associate Study of Gene Presence-Absence Variation in Holstein Cattle

Presence-absence variations (PAVs) are important structural variations, wherein a genomic segment containing one or more genes is present in some individuals but absent in others. While PAVs have been extensively studied in plants, research in cattle remains limited. This study identified PAVs in 173 Holstein bulls using whole-genome sequencing data and assessed their associations with 46 economically important traits. Out of 28,772 cattle genes (from the longest transcripts), a total of 26,979 (93.77%) core genes were identified (present in all individuals), while variable genes included 928 softcore (present in 95-99% of individuals), 494 shell (present in 5-94%), and 371 cloud genes (present in <5%). Cloud genes were enriched in functions associated with hormonal and antimicrobial activities, while shell genes were enriched in immune functions. PAV-based genome-wide association studies identified associations between gene PAVs and 16 traits including milk, fat, and protein yields, as well as traits related to health and reproduction. Associations were found on multiple chromosomes, illustrating important associations on cattle chromosomes 7 and 15, involving olfactory receptor and immune-related genes, respectively. By examining the PAVs at the population level, the results of this research provided crucial insights into the genetic structures underlying the complex traits of Holstein cattle.

slc26a12-A novel member of the slc26 family, is located in tandem with slc26a2 in coelacanths, amphibians, reptiles, and birds

Solute carrier family 26 (Slc26) is a family of anion exchangers with 11 members in mammals (named Slc26a1-a11). Here, we identified a novel member of the slc26 family, slc26a12, located in tandem with slc26a2 in the genomes of several vertebrate lineages. BLAST and synteny analyses of various jawed vertebrate genome databases revealed that slc26a12 is present in coelacanths, amphibians, reptiles, and birds but not in cartilaginous fishes, lungfish, mammals, or ray-finned fishes. In some avian and reptilian lineages such as owls, penguins, egrets, and ducks, and most turtles examined, slc26a12 was lost or pseudogenized. Phylogenetic analysis showed that Slc26a12 formed an independent branch with the other Slc26 members and Slc26a12, Slc26a1 and Slc26a2 formed a single branch, suggesting that these three members formed a subfamily in Slc26. In jawless fish, hagfish have two genes homologous to slc26a2 and slc26a12, whereas lamprey has a single gene homologous to slc26a2. African clawed frogs express slc26a12 in larval gills, skin, and fins. These results show that slc26a12 was present at least before the separation of lobe-finned fish and tetrapods; the name slc26a12 is appropriate because the gene duplication occurred in the distant past.

Conservation of a Chromosome 8 Inversion and Exon Mutations Confirm Common Gulonolactone Oxidase Gene Evolution Among Primates, Including H. Neanderthalensis

Ascorbic acid functions as an antioxidant and facilitates other biochemical processes such as collagen triple helix formation, and iron uptake by cells. Animals which endogenously produce ascorbic acid have a functional gulonolactone oxidase gene (GULO); however, humans have a GULO pseudogene (GULOP) and depend on dietary ascorbic acid. In this study, the conservation of GULOP sequences in the primate haplorhini suborder were investigated and compared to the GULO sequences belonging to the primates strepsirrhini suborder. Phylogenetic analysis suggested that the conserved GULOP exons in the haplorhini primates experienced a high rate of mutations following the haplorhini/strepsirrhini divergence. This high mutation rate has decreased during the evolution of the haplorhini primates. Additionally, indels of the haplorhini GULOP sequences were conserved across the suborder. A separate analysis for GULO sequences and well-conserved GULOP sequences focusing on placental mammals identified an in-frame GULO sequence in the Brazilian guinea pig, and a potential GULOP sequence in the pika. Similar to haplorhini primates, the guinea pig and lagomorph species have experienced a high substitution rate when compared to the mammals used in this study. A shared synteny to examine the conservation of local genes near GULO/GULOP identified a conserved inversion around the GULO/GULOP locus between the haplorhini and strepsirrhini primates. Fischer's exact test did not support an association between GULOP and the chromosomal inversion. Mauve alignment showed that the inversion of the length of the syntenic block that the GULO/GULOP genes belonged to was variable. However, there were frequent rearrangements around ~ 2 million base pairs adjacent to GULOP involving the KIF13B and MSRA genes. These data may suggest that genes acquiring deleterious mutations in the coding sequence may respond to these deleterious mutations with rapid substitution rates.

Identifying the stability of housekeeping genes to be used for the quantitative real-time PCR normalization in retinal tissue of streptozotocin-induced diabetic rats

AIM

To investigate the stability of the seven housekeeping genes: beta-actin (ActB), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), 18s ribosomal unit 5 (18s), cyclophilin A (CycA), hypoxanthine-guanine phosphoribosyl transferase (HPRT), ribosomal protein large P0 (36B4) and terminal uridylyl transferase 1 (U6) in the diabetic retinal tissue of rat model.

METHODS

The expression of these seven genes in rat retinal tissues was determined using real-time quantitative reverse transcription polymerase chain reaction (RT-qPCR) in two groups; normal control rats and streptozotocin-induced diabetic rats. The stability analysis of gene expression was investigated using geNorm, NormFinder, BestKeeper, and comparative delta-Ct (ΔCt) algorithms.

RESULTS

The 36B4 gene was stably expressed in the retinal tissues of normal control animals; however, it was less stable in diabetic retinas. The 18s gene was expressed consistently in both normal control and diabetic rats' retinal tissue. That this gene was the best reference for data normalisation in RT-qPCR studies that used the retinal tissue of streptozotocin-induced diabetic rats. Furthermore, there was no ideal gene stably expressed for use in all experimental settings.

CONCLUSION

Identifying relevant genes is a need for achieving RT-qPCR validity and reliability and must be appropriately achieved based on a specific experimental setting.

Whole genome sequencing of mouse lines divergently selected for fatness (FLI) and leanness (FHI) revealed several genetic variants as candidates for novel obesity genes

BACKGROUND

Analysing genomes of animal model organisms is widely used for understanding the genetic basis of complex traits and diseases, such as obesity, for which only a few mouse models exist, however, without their lean counterparts.

OBJECTIVE

To analyse genetic differences in the unique mouse models of polygenic obesity (Fat line) and leanness (Lean line) originating from the same base population and established by divergent selection over more than 60 generations.

METHODS

Genetic variability was analysed using WGS. Variants were identified with GATK and annotated with Ensembl VEP. g.Profiler, WebGestalt, and KEGG were used for GO and pathway enrichment analysis. miRNA seed regions were obtained with miRPathDB 2.0, LncRRIsearch was used to predict targets of identified lncRNAs, and genes influencing adipose tissue amount were searched using the IMPC database.

RESULTS

WGS analysis revealed 6.3 million SNPs, 1.3 million were new. Thousands of potentially impactful SNPs were identified, including within 24 genes related to adipose tissue amount. SNP density was highest in pseudogenes and regulatory RNAs. The Lean line carries SNP rs248726381 in the seed region of mmu-miR-3086-3p, which may affect fatty acid metabolism. KEGG analysis showed deleterious missense variants in immune response and diabetes genes, with food perception pathways being most enriched. Gene prioritisation considering SNP GERP scores, variant consequences, and allele comparison with other mouse lines identified seven novel obesity candidate genes: 4930441H08Rik, Aff3, Fam237b, Gm36633, Pced1a, Tecrl, and Zfp536.

CONCLUSION

WGS revealed many genetic differences between the lines that accumulated over the selection period, including variants with potential negative impacts on gene function. Given the increasing availability of mouse strains and genetic polymorphism catalogues, the study is a valuable resource for researchers to study obesity.

Breaking paradigms: Long non-coding RNAs forming gene fusions with potential implications in cancer

Long non-coding RNAs (lncRNAs) are non-coding RNAs longer than 200 nucleotides with dynamic regulatory functions. They interact with a wide range of molecules such as DNA, RNA, and proteins to modulate diverse cellular functions through several mechanisms and, if deregulated, they can lead to cancer development and progression. Recently, it has been described that lncRNAs are susceptible to form gene fusions with mRNAs or other lncRNAs, breaking the paradigm of gene fusions consisting mainly of protein-coding genes. However, their biological significance in the tumor phenotype is still uncertain. Therefore, their recent identification opens a new line of research to study their biological role in tumorigenesis, and their potential as biomarkers with clinical relevance or as therapeutic targets. The present study aimed to review the lncRNA fusions identified so far and to know which of them have been associated with a potential function. We address the current challenges to deepen their study as well as the reasons why they represent a future therapeutic window in cancer.

The effect of psychoactive bacteria, Bifidobacterium longum Rosell®-175 and Lactobacillus rhamnosus JB-1, on brain proteome profiles in mice

RATIONALE

The gut microbiota may play an important role in the development and functioning of the mammalian central nervous system. The assumption of the experiment was to prove that the use of probiotic bacterial strains in the diet of mice modifies the expression of brain proteins involved in metabolic and immunological processes.

OBJECTIVES AND RESULTS

Albino Swiss mice were administered with Bifidobacterium longum Rosell®-175 or Lactobacillus rhamnosus JB-1 every 24 h for 28 days. Protein maps were prepared from hippocampal homogenates of euthanized mice. Selected proteins that were statistically significant were purified and concentrated and identified using MALDI-TOF mass spectrometry. Among the analysed samples, 13 proteins were identified. The mean volumes of calcyon, secreted frizzled-associated protein 3, and catalase in the hippocampus of mice from both experimental groups were statistically significantly higher than in the control group. In mice supplemented with Lactobacillus rhamnosus JB-1, a lower mean volume of fragrance binding protein 2, shadow of prion protein, and glycine receptor α4 subunit was observed compared to the control.

CONCLUSION

The psychobiotics Bifidobacterium longum Rosell®-175 and Lactobacillus rhamnosus JB-1enhances expression of proteins involved in the activation and maturation of nerve cells, as well as myelination and homeostatic regulation of neurogenesis in mice. The tested psychobiotics cause a decrease in the expression of proteins associated with CNS development and in synaptic transmission, thereby reducing the capacity for communication between nerve cells. The results of the study indicate that psychobiotic bacteria can be used in auxiliary treatment of neurological disorders.

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.

Blood Transcriptomics Identifies Multiple Gene Expression Pathways Associated with the Clinical Efficacy of Hymenoptera Venom Immunotherapy

Allergen-specific venom immunotherapy (VIT) is a well-established therapy for Hymenoptera venom allergy (HVA). However, the precise mechanism underlying its clinical effect remains uncertain. Our study aimed to identify the molecular mechanisms associated with VIT efficiency. We prospectively included 19 patients with HVA undergoing VIT (sampled before the beginning of VIT, after reaching the maintenance dose, one year after finishing VIT, and after a sting challenge) and 9 healthy controls. RNA sequencing of whole blood was performed on an Illumina sequencing platform. Longitudinal transcriptomic profiling revealed the importance of the inhibition of the NFκB pathway and the downregulation of DUX4 transcripts for the early protection and induction of tolerance after finishing VIT. Furthermore, successful treatment was associated with inhibiting Th2, Th17, and macrophage alternative signalling pathways in synergy with the inhibition of the PPAR pathway and further silencing of the Th2 response. The immune system became activated when reaching the maintenance dose and was suppressed after finishing VIT. Finally, successful VIT restores the immune system's balance to a state similar to that of healthy individuals. Our results underline the important role of the inhibition of four pathways in the clinical effect of VIT: Th2, Th17, NFκB, and macrophage signalling. Two biomarkers specific for successful VIT, regardless of the time of sampling, were C4BPA and RPS10-NUDT3 and should be further tested as potential biomarkers.

A human stomach cell type transcriptome atlas

BACKGROUND

The identification of cell type-specific genes and their modification under different conditions is central to our understanding of human health and disease. The stomach, a hollow organ in the upper gastrointestinal tract, provides an acidic environment that contributes to microbial defence and facilitates the activity of secreted digestive enzymes to process food and nutrients into chyme. In contrast to other sections of the gastrointestinal tract, detailed descriptions of cell type gene enrichment profiles in the stomach are absent from the major single-cell sequencing-based atlases.

RESULTS

Here, we use an integrative correlation analysis method to predict human stomach cell type transcriptome signatures using unfractionated stomach RNAseq data from 359 individuals. We profile parietal, chief, gastric mucous, gastric enteroendocrine, mitotic, endothelial, fibroblast, macrophage, neutrophil, T-cell, and plasma cells, identifying over 1600 cell type-enriched genes.

CONCLUSIONS

We uncover the cell type expression profile of several non-coding genes strongly associated with the progression of gastric cancer and, using a sex-based subset analysis, uncover a panel of male-only chief cell-enriched genes. This study provides a roadmap to further understand human stomach biology.

Four classic "de novo" genes all have plausible homologs and likely evolved from retro-duplicated or pseudogenic sequences

Despite being previously regarded as extremely unlikely, the idea that entirely novel protein-coding genes can emerge from non-coding sequences has gradually become accepted over the past two decades. Examples of "de novo origination", resulting in lineage-specific "orphan" genes, lacking coding orthologs, are now produced every year. However, many are likely cases of duplicates that are difficult to recognize. Here, I re-examine the claims and show that four very well-known examples of genes alleged to have emerged completely "from scratch"- FLJ33706 in humans, Goddard in fruit flies, BSC4 in baker's yeast and AFGP2 in codfish-may have plausible evolutionary ancestors in pre-existing genes. The first two are likely highly diverged retrogenes coding for regulatory proteins that have been misidentified as orphans. The antifreeze glycoprotein, moreover, may not have evolved from repetitive non-genic sequences but, as in several other related cases, from an apolipoprotein that could have become pseudogenized before later being reactivated. These findings detract from various claims made about de novo gene birth and show there has been a tendency not to invest the necessary effort in searching for homologs outside of a very limited syntenic or phylostratigraphic methodology. A robust approach is used for improving detection that draws upon similarities, not just in terms of statistical sequence analysis, but also relating to biochemistry and function, to obviate notable failures to identify homologs.

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.

Blastomeres of 8-cell mouse embryos differ in their ability to generate embryonic stem cells and produce lines with different transcriptional signatures

Embryonic stem cell (ESC) derivation from single blastomeres of 8-cell mouse embryos results in lower derivation rates than that from whole blastocysts, raising a biological question about the developmental potential of sister blastomeres. We aimed to assess the ability of 8-cell blastomeres to produce epiblast cells and ESC lines after isolation, and the properties of the resulting lines. Our results revealed unequal competence among sister blastomeres to produce ESC lines. At least half of the blastomeres possess a lower potential to generate ESCs, although culture conditions and blastomeres plasticity can redirect their non-pluripotent fate towards the epiblast lineage, allowing us to generate up to seven lines from the same embryo. Lines originated from the same embryo segregated into two groups according to their transcriptional signatures. While the expression of genes related to pluripotency and development was higher in one group, no differences were found in their trilineage differentiation ability. These results may help to improve our understanding of the ESC derivation process from single blastomeres and cell fate determination in the preimplantation mouse embryos.

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.

CCIVR2 facilitates comprehensive identification of both overlapping and non-overlapping antisense transcripts within specified regions

Pairs of sense and antisense transcriptions that are adjacent at their 5' and 3' regions are called divergent and convergent transcription, respectively. However, the structural properties of divergent/convergent transcription in different species or RNA biotypes are poorly characterized. Here, we developed CCIVR2, a program that facilitates identification of both overlapping and non-overlapping antisense transcripts produced from divergent/convergent transcription whose transcription start sites (TSS) or transcript end sites (TES) are located within a specified region. We used CCIVR2 to analyze antisense transcripts starting around the sense TSS (from divergent transcription) or ending around the sense TES (from convergent transcription) in 11 different species and found species- and RNA biotype-specific features of divergent/convergent transcription. Furthermore, we confirmed that CCIVR2 enables the identification of multiple sense/antisense transcript pairs from divergent transcription, including those with known functions in processes such as embryonic stem cell differentiation and TGFβ stimulation. CCIVR2 is therefore a valuable bioinformatics tool that facilitates the characterization of divergent/convergent transcription in different species and aids the identification of functional sense/antisense transcript pairs from divergent transcription in specified biological processes.

RNA-RNA competitive interactions: a molecular civil war ruling cell physiology and diseases

The idea that proteins are the main determining factors in the functioning of cells and organisms, and their dysfunctions are the first cause of pathologies, has been predominant in biology and biomedicine until recently. This protein-centered view was too simplistic and failed to explain the physiological and pathological complexity of the cell. About 80% of the human genome is dynamically and pervasively transcribed, mostly as non-protein-coding RNAs (ncRNAs), which competitively interact with each other and with coding RNAs generating a complex RNA network regulating RNA processing, stability, and translation and, accordingly, fine-tuning the gene expression of the cells. Qualitative and quantitative dysregulations of RNA-RNA interaction networks are strongly involved in the onset and progression of many pathologies, including cancers and degenerative diseases. This review will summarize the RNA species involved in the competitive endogenous RNA network, their mechanisms of action, and involvement in pathological phenotypes. Moreover, it will give an overview of the most advanced experimental and computational methods to dissect and rebuild RNA networks.

A pan-cancer analysis of pituitary tumor-transforming 3, pseudogene

BACKGROUND

Although evidence regarding pituitary tumor-transforming 3, pseudogene (PTTG3P) involvement in human cancers has been acquired via human and animal model-based molecular studies, there is a lack of pan-cancer analysis of this gene in human tumors.

METHODS

Tumor-causing effects of PTTG3P in 24 human tumors were explored using The Cancer Genome Atlas (TCGA) datasets from different bioinformatics databases and applying in silico tools such as The University of ALabama at Birmingham CANcer (UALCAN), Human Protein Atlas (HPA), Kaplan Meier (KM) plotter, cBioPortal, Search Tool for the Retrieval of Interacting Genes/Proteins (STRING), Cytoscape, Database for Annotation, Visualization, and Integrated Discovery (DAVID), Tumor IMmune Estimation Resource (TIMER), and Comparative Toxicogenomics Database (CTD). Then, via in vitro experiments, including RNA sequencing (RNA-seq) and targeted bisulfite sequencing (bisulfite-seq), expression and promoter methylation levels of PTTG3P were verified in cell lines.

RESULTS

The PTTG3P expression was overexpressed across 23 malignancies and its overexpression was further found significantly effecting the overall survival (OS) durations of the esophageal carcinoma (ESCA) and head and neck cancer (HNSC) patients. This important information helps us to understand that PTTG3P plays a significant role in the development and progression of ESCA and HNSC. As for PTTG3P functional mechanisms, this gene along with its other binding partners was significantly concentrated in "Oocyte meiosis", "Cell cycle", "Ubiquitin mediated proteolysis", and "Progesterone-mediated oocyte maturation". Moreover, ESCA and HNSC tissues having the higher expression of PTTG3P were found to have lower promoter methylation levels of PTTG3P and higher CD8+ T immune cells level. Additionally, PTTG3P expression-regulatory drugs were also explored in the current manuscript for designing appropriate treatment strategies for ESCA and HNSC with respect to PTTG3P expression.

CONCLUSION

Our pan-cancer based findings provided a comprehensive account of the oncogenic role and utilization of PTTG3P as a novel molecular biomarker of ESCA and HNSC.

Pseudogenes in Cancer: State of the Art

Pseudogenes are duplicates of protein-coding genes that have accumulated multiple detrimental alterations, rendering them unable to produce the protein they encode. Initially disregarded as "junk DNA" due to their perceived lack of functionality, research on their biological roles has been hindered by this assumption. Nevertheless, recent focus has shifted towards these molecules due to their abnormal expression in cancer phenotypes. In this review, our objective is to provide a thorough overview of the current understanding of pseudogene formation, the mechanisms governing their expression, and the roles they may play in promoting tumorigenesis.

Genome-wide association study identifies multiple HLA loci for sarcoidosis susceptibility

Sarcoidosis is a complex systemic disease. Our study aimed to (1) identify novel alleles associated with sarcoidosis susceptibility; (2) provide an in-depth evaluation of HLA alleles and sarcoidosis susceptibility and (3) integrate genetic and transcription data to identify risk loci that may more directly impact disease pathogenesis. We report a genome-wide association study of 1335 sarcoidosis cases and 1264 controls of European descent (EA) and investigate associated alleles in a study of African Americans (AA: 1487 cases and 1504 controls). The EA and AA cohort was recruited from multiple United States sites. HLA alleles were imputed and tested for association with sarcoidosis susceptibility. Expression quantitative locus and colocalization analysis were performed using a subset of subjects with transcriptome data. Forty-nine SNPs in the HLA region in HLA-DRA, -DRB9, -DRB5, -DQA1 and BRD2 genes were significantly associated with sarcoidosis susceptibility in EA, rs3129888 was also a risk variant for sarcoidosis in AA. Classical HLA alleles DRB1*0101, DQA1*0101 and DQB1*0501, which are highly correlated, were also associated with sarcoidosis. rs3135287 near HLA-DRA was associated with HLA-DRA expression in peripheral blood mononuclear cells and bronchoalveolar lavage from subjects and lung tissue and whole blood from GTEx. We identified six novel SNPs (out of the seven SNPs representing the 49 significant SNPs) and nine HLA alleles associated with sarcoidosis susceptibility in the largest EA population. We also replicated our findings in an AA population. Our study reiterates the potential role of antigen recognition and/or presentation HLA class II genes in sarcoidosis pathogenesis.

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

BACKGROUND

Voltage-dependent anion selective channels (VDACs) are the most abundant mitochondrial outer membrane proteins, encoded in mammals by three genes, VDAC1, 2 and 3, mostly ubiquitously expressed. As 'mitochondrial gatekeepers', VDACs control organelle and cell metabolism and are involved in many diseases. Despite the presence of numerous VDAC pseudogenes in the human genome, their significance and possible role in VDAC protein expression has not yet been considered.

RESULTS

We investigated the relevance of processed pseudogenes of human VDAC genes, both in physiological and in pathological contexts. Using high-throughput tools and querying many genomic and transcriptomic databases, we show that some VDAC pseudogenes are transcribed in specific tissues and pathological contexts. The obtained experimental data confirm an association of the VDAC1P8 pseudogene with acute myeloid leukemia (AML).

CONCLUSIONS

Our in-silico comparative analysis between the VDAC1 gene and its VDAC1P8 pseudogene, together with experimental data produced in AML cellular models, indicate a specific over-expression of the VDAC1P8 pseudogene in AML, correlated with a downregulation of the parental VDAC1 gene.

c-Myc/microRNA-17-92 Axis Phase-Dependently Regulates PTEN and p21 Expression via ceRNA during Reprogramming to Mouse Pluripotent Stem Cells

Induced pluripotent stem cells (iPSCs) are promising cell sources for regenerative medicine and disease modeling. iPSCs are commonly established by introducing the defined reprogramming factors Oct4, Sox2, Klf4, and c-Myc. However, iPSC reprogramming efficiency remains low. Although recent studies have identified microRNAs that contribute to efficient reprogramming, the underlying molecular mechanisms are not completely understood. miR-17-92 is highly expressed in embryonic stem cells and may play an important role in regulating stem cell properties. Therefore, we examined the role of miR-17-92 in the induction of mouse iPSC production. c-Myc-mediated miR-17-92 upregulation increased reprogramming efficiency, whereas CRISPR/Cas9-based deletion of the miR-17-92 cluster decreased reprogramming efficiency. A combination of in silico and microarray analyses revealed that Pten and cyclin-dependent kinase inhibitor 1 (known as p21) are common target genes of miR-17 and miR-20a, which are transcribed from the miR-17-92 cluster. Moreover, miR-17-92 downregulated p21 in the early phase and PTEN in the mid-to-late phase of reprogramming. These downregulations were perturbed by introducing the 3' UTR of PTEN and p21, respectively, suggesting that PTEN and p21 mRNAs are competing endogenous RNAs (ceRNA) against miR-17-92. Collectively, we propose that the c-Myc-mediated expression of miR-17-92 is involved in iPSC reprogramming through the phase-dependent inhibition of PTEN and p21 in a ceRNA manner, thus elucidating an underlying mechanism of iPSC reprogramming.

The non-coding competing endogenous RNAs in acute myeloid leukemia: biological and clinical implications

Acute myeloid leukemia (AML) is a hematologic carcinoma that has seen a considerable improvement in patient prognosis because of genetic diagnostics and molecularly-targeted therapies. Nevertheless, recurrence and drug resistance remain significant obstacles to leukemia treatment. It is critical to investigate the underlying molecular mechanisms and find solutions. Non-coding RNAs (ncRNAs), such as microRNAs (miRNAs), circular RNAs, long non-coding RNAs, and pseudogenes, have been found to be crucial components in driving cancer. The competing endogenous RNA (ceRNA) mechanism has expanded the complexity of miRNA-mediated gene regulation. A great deal of literature has shown that ncRNAs are essential to the biological functions of the ceRNA network (ceRNET). NcRNAs can compete for the same miRNA response elements to influence miRNA-target RNA interactions. Recent evidence suggests that ceRNA might be a potential biomarker and therapeutic strategy. So far, however, there have been no comprehensive studies on ceRNET about AML. What is not yet clear is the clinical application of ceRNA in AML. This study attempts to summarize the development of research on the related ceRNAs in AML and the roles of ncRNAs in ceRNET. We also briefly describe the mechanisms of ceRNA and ceRNET. What's more significant is that we explore the clinical value of ceRNAs to provide accurate diagnostic and prognostic biomarkers as well as therapeutic targets. Finally, limitations and prospects are considered.

The X-factor in ART: does the use of assisted reproductive technologies influence DNA methylation on the X chromosome?

BACKGROUND

Assisted reproductive technologies (ART) may perturb DNA methylation (DNAm) in early embryonic development. Although a handful of epigenome-wide association studies of ART have been published, none have investigated CpGs on the X chromosome. To bridge this knowledge gap, we leveraged one of the largest collections of mother-father-newborn trios of ART and non-ART (natural) conceptions to date to investigate sex-specific DNAm differences on the X chromosome. The discovery cohort consisted of 982 ART and 963 non-ART trios from the Norwegian Mother, Father, and Child Cohort Study (MoBa). To verify our results from the MoBa cohort, we used an external cohort of 149 ART and 58 non-ART neonates from the Australian 'Clinical review of the Health of adults conceived following Assisted Reproductive Technologies' (CHART) study. The Illumina EPIC array was used to measure DNAm in both datasets. In the MoBa cohort, we performed a set of X-chromosome-wide association studies ('XWASs' hereafter) to search for sex-specific DNAm differences between ART and non-ART newborns. We tested several models to investigate the influence of various confounders, including parental DNAm. We also searched for differentially methylated regions (DMRs) and regions of co-methylation flanking the most significant CpGs. Additionally, we ran an analogous model to our main model on the external CHART dataset.

RESULTS

In the MoBa cohort, we found more differentially methylated CpGs and DMRs in girls than boys. Most of the associations persisted after controlling for parental DNAm and other confounders. Many of the significant CpGs and DMRs were in gene-promoter regions, and several of the genes linked to these CpGs are expressed in tissues relevant for both ART and sex (testis, placenta, and fallopian tube). We found no support for parental DNAm-dependent features as an explanation for the observed associations in the newborns. The most significant CpG in the boys-only analysis was in UBE2DNL, which is expressed in testes but with unknown function. The most significant CpGs in the girls-only analysis were in EIF2S3 and AMOT. These three loci also displayed differential DNAm in the CHART cohort.

CONCLUSIONS

Genes that co-localized with the significant CpGs and DMRs associated with ART are implicated in several key biological processes (e.g., neurodevelopment) and disorders (e.g., intellectual disability and autism). These connections are particularly compelling in light of previous findings indicating that neurodevelopmental outcomes differ in ART-conceived children compared to those naturally conceived.

Astrocytes express aberrant immunoglobulins as putative gatekeeper of astrocytes to neuronal progenitor conversion

Using multi-omics analyses including RNAseq, RT-PCR, RACE-PCR, and shotgun proteomic with enrichment strategies, we demonstrated that newborn rat astrocytes produce neural immunoglobulin constant and variable heavy chains as well as light chains. However, their edification is different from the ones found in B cells and they resemble aberrant immunoglobulins observed in several cancers. Moreover, the complete enzymatic V(D)J recombination complex has also been identified in astrocytes. In addition, the constant heavy chain is also present in adult rat astrocytes, whereas in primary astrocytes from human fetus we identified constant and variable kappa chains as well as the substitution lambda chains known to be involved in pre-B cells. To gather insights into the function of these neural IgGs, CRISPR-Cas9 of IgG2B constant heavy chain encoding gene (Igh6), IgG2B overexpression, proximal labeling of rat astrocytes IgG2B and targets identification through 2D gels were performed. In Igh6 KO astrocytes, overrepresentation of factors involved in hematopoietic cells, neural stem cells, and the regulation of neuritogenesis have been identified. Moreover, overexpression of IgG2B in astrocytes induces the CRTC1-CREB-BDNF signaling pathway known to be involved in gliogenesis, whereas Igh6 KO triggers the BMP/YAP1/TEAD3 pathway activated in astrocytes dedifferentiation into neural progenitors. Proximal labeling experiments revealed that IgG2B is N-glycosylated by the OST complex, addressed to vesicle membranes containing the ATPase complex, and behaves partially like CD98hc through its association with LAT1. These experiments also suggest that proximal IgG2B-LAT1 interaction occurs concomitantly with MACO-1 and C2CD2L, at the heart of a potentially novel cell signaling platform. Finally, we demonstrated that these chains are synthesized individually and associated to recognize specific targets. Indeed, intermediate filaments Eif4a2 and Pdia6 involved in astrocyte fate constitute targets for these neural IgGs. Taken together, we hypothese that neural aberrant IgG chains may act as gatekeepers of astrocytes' fate.

Specialized Ribosomes in Health and Disease

Ribosomal heterogeneity exists within cells and between different cell types, at specific developmental stages, and occurs in response to environmental stimuli. Mounting evidence supports the existence of specialized ribosomes, or specific changes to the ribosome that regulate the translation of a specific group of transcripts. These alterations have been shown to affect the affinity of ribosomes for certain mRNAs or change the cotranslational folding of nascent polypeptides at the exit tunnel. The identification of specialized ribosomes requires evidence of the incorporation of different ribosomal proteins or of modifications to rRNA and/or protein that lead(s) to physiologically relevant changes in translation. In this review, we summarize ribosomal heterogeneity and specialization in mammals and discuss their relevance to several human diseases.

Polystyrene nanoplastics affect transcriptomic and epigenomic signatures of human fibroblasts and derived induced pluripotent stem cells: Implications for human health

Plastic pollution is increasing at an alarming rate yet the impact of this pollution on human health is poorly understood. Because human induced pluripotent stem cells (hiPSC) are frequently derived from dermal fibroblasts, these cells offer a powerful platform for the identification of molecular biomarkers of environmental pollution in human cells. Here, we describe a novel proof-of-concept for deriving hiPSC from human dermal fibroblasts deliberately exposed to polystyrene (PS) nanoplastic particles; unexposed hiPSC served as controls. In parallel, unexposed hiPSC were exposed to low and high concentrations of PS nanoparticles. Transcriptomic and epigenomic signatures of all fibroblasts and hiPSCs were defined using RNA-seq and whole genome methyl-seq, respectively. Both PS-treated fibroblasts and derived hiPSC showed alterations in expression of ESRRB and HNF1A genes and circuits involved in the pluripotency of stem cells, as well as in pathways involved in cancer, inflammatory disorders, gluconeogenesis, carbohydrate metabolism, innate immunity, and dopaminergic synapse. Similarly, the expression levels of identified key transcriptional and DNA methylation changes (DNMT3A, ESSRB, FAM133CP, HNF1A, SEPTIN7P8, and TTC34) were significantly affected in both PS-exposed fibroblasts and hiPSC. This study illustrates the power of human cellular models of environmental pollution to narrow down and prioritize the list of candidate molecular biomarkers of environmental pollution. This knowledge will facilitate the deciphering of the origins of environmental diseases.

Na+/Cl- cotransporter 2 is not fish-specific and is widely found in amphibians, non-avian reptiles, and select mammals

Solute carrier 12 (Slc12) is a family of electroneutral cation-coupled chloride (Cl-) cotransporters. Na+/K+/2Cl- (Nkcc) and Na+/Cl- cotransporters (Ncc) belong to the Nkcc/Ncc subfamily. Human and mouse possess one gene for the Na+/Cl- cotransporter (ncc gene: slc12a3), whereas teleost fishes possess multiple ncc genes, slc12a3 (ncc1) and slc12a10 (ncc2), in addition to their species-specific paralogs. Amphibians and squamates have two ncc genes: slc12a3 (ncc1) and ncc3. However, the evolutionary relationship between slc12a10 and ncc3 remains unresolved, and the presence of slc12a10 (ncc2) in mammals has not been clarified. Synteny and phylogenetic analyses of vertebrate genome databases showed that ncc3 is the ortholog of slc12a10, and slc12a10 is present in most ray-finned fishes, coelacanths, amphibians, reptiles, and a few mammals (e.g., platypus and horse) but pseudogenized or deleted in birds, most mammals, and some ray-finned fishes (pufferfishes). This shows that slc12a10 is widely present among bony vertebrates and pseudogenized or deleted independently in multiple lineages. Notably, as compared with some fish that show varied slc12a10 tissue expression profile, spotted gar, African clawed frog, red-eared slider turtle, and horse express slc12a10 in the ovaries or premature gonads. In horse tissues, an unexpectedly large number of splicing variants for Slc12a10 have been cloned, many of which encode truncated forms of Slc12a10, suggesting that the functional constraints of horse slc12a10 are weakened, which may be in the process of becoming a pseudogene. Our results elaborate on the evolution of Nkcc/Ncc subfamily of Slc12 in vertebrates.NEW & NOTEWORTHY slc12a10 is not a fish-specific gene and is present in a few mammals (e.g., platypus and horse), non-avian reptiles, amphibians, but was pseudogenized or deleted in most mammals (e.g., human, mouse, cat, cow, and rhinoceros), birds, and some ray-finned fishes (pufferfishes).

A comprehensive survey of C. elegans argonaute proteins reveals organism-wide gene regulatory networks and functions

Argonaute (AGO) proteins associate with small RNAs to direct their effector function on complementary transcripts. The nematode Caenorhabditis elegans contains an expanded family of 19 functional AGO proteins, many of which have not been fully characterized. In this work, we systematically analyzed every C. elegans AGO using CRISPR-Cas9 genome editing to introduce GFP::3xFLAG tags. We have characterized the expression patterns of each AGO throughout development, identified small RNA binding complements, and determined the effects of ago loss on small RNA populations and developmental phenotypes. Our analysis indicates stratification of subsets of AGOs into distinct regulatory modules, and integration of our data led us to uncover novel stress-induced fertility and pathogen response phenotypes due to ago loss.

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

BACKGROUND

Preeclampsia (PE) is one of the most serious pregnancy complications with unknown pathogenesis. Emerging evidence has demonstrated that Fms-related tyrosine kinase 1 (FLT1) is highly involved in PE development. As a pseudogene of FLT1, FLT1P1 increased in PE samples. However, its functions remain largely unknown.

METHODS AND RESULTS

In this study, co-expression analysis was performed to identify the potential target genes of FTL1P1. Then chromatin isolation using RNA purification (ChIRP) method was employed to explore the interactomes of FLT1P1, including interacting with DNA fragments and proteins. We found that in PE samples, both FLT1P1 and FLT1 were highly expressed and closely correlated. ChIRP-protein data revealed that FLT1P1 interacts with translation- and transcription-related proteins, including 4 transcription factors (TFs). ChIRP-DNA analysis revealed that FLT1P1 preferentially interacted with DNA fragments downstream of transcription start sites (TSSs). Functional analysis of its interacting genes revealed that they were enriched in transcriptional regulation and apoptosis-related pathways. Twenty-six TFs, including CREB1 and SRF, were extracted from the potential FLT1P1-interacting gene sets and were potential targets of FLT1P1. CREB1 could bind to FLT1 promoter, and was negatively correlated with FLT1 at the expression level, making it a potential regulator of FLT1.

CONCLUSIONS

Our study extensively investigated the interactome profiles of FLT1P1, especially the prompter region of TF gene CREB1, and revealed the potential molecular regulatory mechanisms of FLT1 expression in PE samples. Our results provide a novel view of PE pathogenesis, and suggest that FLT1P1 could serve as a potential therapeutic target in PE diagnosis and treatment.

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.

Protein-Coding Region Derived Small RNA in Exosomes from Influenza A Virus-Infected Cells

Exosomes may function as multifactorial mediators of cell-to-cell communication, playing crucial roles in both physiological and pathological processes. Exosomes released from virus-infected cells may contain RNA and proteins facilitating infection spread. The purpose of our study was to analyze how the small RNA content of exosomes is affected by infection with the influenza A virus (IAV). Exosomes were isolated by ultracentrifugation after hemadsorption of virions and their small RNA content was identified using high-throughput sequencing. As compared to mock-infected controls, 856 RNA transcripts were significantly differentially expressed in exosomes from IAV-infected cells, including fragments of 458 protein-coding (pcRNA), 336 small, 28 long intergenic non-coding RNA transcripts, and 33 pseudogene transcripts. Upregulated pcRNA species corresponded mainly to proteins associated with translation and antiviral response, and the most upregulated among them were RSAD2, CCDC141 and IFIT2. Downregulated pcRNA species corresponded to proteins associated with the cell cycle and DNA packaging. Analysis of differentially expressed pseudogenes showed that in most cases, an increase in the transcription level of pseudogenes was correlated with an increase in their parental genes. Although the role of exosome RNA in IAV infection remains undefined, the biological processes identified based on the corresponding proteins may indicate the roles of some of its parts in IAV replication.

Single-Cell Transcriptome Analysis Reveals Paraspeckles Expression in Osteosarcoma Tissues

Nuclear paraspeckles are subnuclear bodies contracted by nuclear-enriched abundant transcript 1 (NEAT1) long non-coding RNA, localised in the interchromatin space of mammalian cell nuclei. Paraspeckles have been critically involved in tumour progression, metastasis and chemoresistance. To this date, there are limited findings to suggest that paraspeckles, NEAT1 and heterogeneous nuclear ribonucleoproteins (hnRNPs) directly or indirectly play roles in osteosarcoma progression. Herein, we analysed NEAT1, paraspeckle proteins (SFPQ, PSPC1 and NONO) and hnRNP members (HNRNPK, HNRNPM, HNRNPR and HNRNPD) gene expression in 6 osteosarcoma tumour tissues using the single-cell RNA-sequencing method. The normalised data highlighted that the paraspeckles transcripts were highly abundant in osteoblastic OS cells, except NEAT1, which was highly expressed in myeloid cell 1 and 2 subpopulations.

An Epigenetic LINE-1-Based Mechanism in Cancer

In the last fifty years, large efforts have been deployed in basic research, clinical oncology, and clinical trials, yielding an enormous amount of information regarding the molecular mechanisms of cancer and the design of effective therapies. The knowledge that has accumulated underpins the complexity, multifactoriality, and heterogeneity of cancer, disclosing novel landscapes in cancer biology with a key role of genome plasticity. Here, we propose that cancer onset and progression are determined by a stress-responsive epigenetic mechanism, resulting from the convergence of upregulation of LINE-1 (long interspersed nuclear element 1), the largest family of human retrotransposons, genome damage, nuclear lamina fragmentation, chromatin remodeling, genome reprogramming, and autophagy activation. The upregulated expression of LINE-1 retrotransposons and their protein products plays a key role in these processes, yielding an increased plasticity of the nuclear architecture with the ensuing reprogramming of global gene expression, including the reactivation of embryonic transcription profiles. Cancer phenotypes would thus emerge as a consequence of the unscheduled reactivation of embryonic gene expression patterns in an inappropriate context, triggering de-differentiation and aberrant proliferation in differentiated cells. Depending on the intensity of the stressing stimuli and the level of LINE-1 response, diverse degrees of malignity would be generated.