Plant Pattern recognition receptors: Exploring their evolution, diversification, and spatiotemporal regulation.

Plant genomes possess hundreds of candidate surface localized receptors capable of recognizing microbial components or modified-self molecules. Surface-localized pattern recognition receptors (PRRs) can recognize proteins, peptides, or structural microbial components as nonself, triggering complex signaling pathways leading to defense. PRRs possess diverse extracellular domains capable of recognizing epitopes, lipids, glycans and polysaccharides. Recent work highlights advances in our understanding of the diversity and evolution of PRRs recognizing pathogen components. We also discuss PRR functional diversification, pathogen strategies to evade detection, and the role of tissue and age-related resistance for effective plant defense.

UV-B Radiation Exhibited Tissue-Specific Regulation of Isoflavone Biosynthesis in Soybean Cell Suspension Cultures.

Isoflavones, a class of substances with high biological activity, are abundant in soybeans. This study investigated isoflavone biosynthesis in soybean cell suspension cultures under UV-B radiation. UV-B radiation enhanced the transcription level and activity of key enzymes involved in isoflavone synthesis in cell suspension cultures. As a result, the isoflavone contents significantly increased by 19.80% and 91.21% in hypocotyl and cotyledon suspension cultures compared with the control, respectively. Meanwhile, a significant difference was observed in the composition of isoflavones between soybean hypocotyl and cotyledon suspension cultures. Genistin was only detected in hypocotyl suspension cultures, whereas glycitin, daidzein, and genistein accumulated in cotyledon suspension cultures. Therefore, UV-B radiation exhibited tissue-specific regulation of isoflavone biosynthesis in soybean cell suspension cultures. The combination of suspension cultures and abiotic stress provides a novel technological approach to isoflavone accumulation.

Cold mediates maize root hair developmental plasticity via epidermis-specific transcriptomic responses.

Cold stress during early development limits maize (Zea mays L.) production in temperate zones. Low temperatures restrict root growth and reprogram gene expression. Here, we provide a systematic transcriptomic landscape of maize primary roots, their tissues, and cell types in response to cold stress. The epidermis exhibited a unique transcriptomic cold response, and genes involved in root hair formation were dynamically regulated in this cell type by cold. Consequently, activation of genes involved in root hair tip growth contributed to root hair recovery under moderate cold conditions. The maize root hair defective mutants roothair defective 5 (rth5) and roothair defective 6 (rth6) displayed enhanced cold tolerance with respect to primary root elongation. Furthermore, dehydration response element-binding protein 2.1 (dreb2.1) was the only member of the dreb subfamily of AP2/EREB transcription factor genes upregulated in primary root tissues and cell types but exclusively downregulated in root hairs upon cold stress. Plants overexpressing dreb2.1 significantly suppressed root hair elongation after moderate cold stress. Finally, the expression of rth3 was regulated by dreb2.1 under cold conditions, while rth6 transcription was regulated by dreb2.1 irrespective of the temperature regime. We demonstrated that dreb2.1 negatively regulates root hair plasticity at low temperatures by coordinating the expression of root hair defective genes in maize.

Prenatal stress induces sex- and tissue-specific alterations in insulin pathway of Wistar rats offspring.

BACKGROUND AND AIMS: Prenatal stress may lead to tissue and sex-specific cardiometabolic disorders in the offspring through imbalances in the insulin signaling pathway. Therefore, we aimed to determine the sex-specific adaptations of prenatal stress on the insulin signaling pathway of cardiac and hepatic tissue of adult offspring Wistar rats. METHODS: Wistar pregnant rats were divided into control and stress groups. Unpredictable stress protocol was performed from the 14th to the 21st day of pregnancy. After lactation, the dams were euthanized and blood was collected for corticosterone measurement and the offspring were separated into four groups according to sex and intervention (n=8/group). At 90 days old, the offspring were submitted to an oral glucose tolerance test (OGTT) and an insulin tolerance test (ITT). After euthanasia blood collection was used for biochemical analysis and the left ventricle and liver were used for protein expression and histological analysis. RESULTS: Stress increased maternal corticosterone levels, and in the offspring, decreased glucose concentration in both OGTT and ITT, reduced insulin receptor (Irß) and insulin receptor substrate-1 (IRS1) activation and reduced insulin receptor inhibition (PTP1B) in the liver of male offspring at 90 days old, without repercussions in cardiac tissue. Moreover, female offspring submitted to prenatal stress exhibited reduced fatty acid uptake, with lower hepatic CD36 expression, reduced high density lipoprotein (cHDL) and increased Castelli risk indexes I and II. CONCLUSIONS: Unpredictable prenatal stress evoked reduced insulin sensitivity and liver-specific impairment in insulin signaling activation in male while increasing markers of cardiovascular risk in females.

Identification of Wheat Glutamate Synthetase Gene Family and Expression Analysis under Nitrogen Stress.

Nitrogen (N), as the main component of biological macromolecules, maintains the basic process of plant growth and development. GOGAT, as a key enzyme in the N assimilation process, catalyzes α-ketoglutaric acid and glutamine to form glutamate. In this study, six GOGAT genes in wheat (Triticum aestivum L.) were identified and classified into two subfamilies, Fd-GOGAT (TaGOGAT2s) and NADH-GOGAT (TaGOGAT3s), according to the type of electron donor. Subcellular localization prediction showed that TaGOGAT3-D was localized in mitochondria and that the other five TaGOGATs were localized in chloroplasts. Via the analysis of promoter elements, many binding sites related to growth and development, hormone regulation and plant stress resistance regulations were found on the TaGOGAT promoters. The tissue-specificity expression analysis showed that TaGOGAT2s were mainly expressed in wheat leaves and flag leaves, while TaGOGAT3s were highly expressed in roots and leaves. The expression level of TaGOGATs and the enzyme activity of TaGOGAT3s in the leaves and roots of wheat seedlings were influenced by the treatment of N deficiency. This study conducted a systematic analysis of wheat GOGAT genes, providing a theoretical basis not only for the functional analysis of TaGOGATs, but also for the study of wheat nitrogen use efficiency (NUE).

GIGANTEA adjusts the response to shade at dusk by directly impinging on PHYTOCHROME INTERACTING FACTOR 7 function.

For plants adapted to bright light, a decrease in the amount of light received can be detrimental to their growth and survival. Consequently, in response to shade from surrounding vegetation, they initiate a suite of molecular and morphological changes known as the shade avoidance response through which stems and petioles elongate in search for light. Under sunlight-night cycles, the plant's responsiveness to shade varies across the day, being maximal at dusk time. While a role for the circadian clock in this regulation has long been proposed, mechanistic understanding of how it is achieved is incomplete. Here, we show that the clock component GIGANTEA (GI) directly interacts with the transcriptional regulator PHYTOCHROME INTERACTING FACTOR 7 (PIF7), a key player in the response to shade. GI represses PIF7 transcriptional activity and the expression of its target genes in response to shade, thereby fine-tuning the magnitude of the response to limiting light conditions. We find that under light/dark cycles, this function of GI is required to adequately modulate the gating of the response to shade at dusk. Importantly, we also show that this circuit primarily operates in epidermal cells, highlighting the relevance of tissue-specific clock-output connections for the regulation of plant development in resonance with the environment.

The cacao gene atlas: a transcriptome developmental atlas reveals highly tissue-specific and dynamically-regulated gene networks in Theobroma cacao L.

BACKGROUND: Theobroma cacao, the cocoa tree, is a tropical crop grown for its highly valuable cocoa solids and fat which are the basis of a 200-billion-dollar annual chocolate industry. However, the long generation time and difficulties associated with breeding a tropical tree crop have limited the progress of breeders to develop high-yielding disease-resistant varieties. Development of marker-assisted breeding methods for cacao requires discovery of genomic regions and specific alleles of genes encoding important traits of interest. To accelerate gene discovery, we developed a gene atlas composed of a large dataset of replicated transcriptomes with the long-term goal of progressing breeding towards developing high-yielding elite varieties of cacao. RESULTS: We describe the creation of the Cacao Transcriptome Atlas, its global characterization and define sets of genes co-regulated in highly organ- and temporally-specific manners. RNAs were extracted and transcriptomes sequenced from 123 different tissues and stages of development representing major organs and developmental stages of the cacao lifecycle. In addition, several experimental treatments and time courses were performed to measure gene expression in tissues responding to biotic and abiotic stressors. Samples were collected in replicates (3-5) to enable statistical analysis of gene expression levels for a total of 390 transcriptomes. To promote wide use of these data, all raw sequencing data, expression read mapping matrices, scripts, and other information used to create the resource are freely available online. We verified our atlas by analyzing the expression of genes with known functions and expression patterns in Arabidopsis (ACT7, LEA19, AGL16, TIP13, LHY, MYB2) and found their expression profiles to be generally similar between both species. We also successfully identified tissue-specific genes at two thresholds in many tissue types represented and a set of genes highly conserved across all tissues. CONCLUSION: The Cacao Gene Atlas consists of a gene expression browser with graphical user interface and open access to raw sequencing data files as well as the unnormalized and CPM normalized read count data mapped to several cacao genomes. The gene atlas is a publicly available resource to allow rapid mining of cacao gene expression profiles. We hope this resource will be used to help accelerate the discovery of important genes for key cacao traits such as disease resistance and contribute to the breeding of elite varieties to help farmers increase yields.

Tissue-specific differences in Ca2+ sensitivity of the mitochondrial permeability transition pore (PTP). Experiments in male rat liver and heart.

Permeability transition pore (PTP) opening dissipates ion and electron gradients across the internal mitochondrial membrane (IMM), including excess Ca2+ in the mitochondrial matrix. After opening, immediate PTP closure must follow to prevent outer membrane disruption, loss of cytochrome c, and eventual apoptosis. Flickering, defined as the rapid alternative opening/closing of PTP, has been reported in heart, which undergoes frequent, large variations in Ca2+. In contrast, in tissues that undergo depolarization events less often, such as the liver, PTP would not need to be as dynamic and thus these tissues would not be as resistant to stress. To evaluate this idea, it was decided to follow the reversibility of the permeability transition (PT) in isolated murine mitochondria from two different tissues: the very dynamic heart, and the liver, which suffers depolarizations less frequently. It was observed that in heart mitochondria PT remained reversible for longer periods and at higher Ca2+ loads than in liver mitochondria. In all cases, Ca2+ uptake was inhibited by ruthenium red and PT was delayed by Cyclosporine A. Characterization of this phenomenon included measuring the rate of oxygen consumption, organelle swelling and Ca2+ uptake and retention. Results strongly suggest that there are tissue-specific differences in PTP physiology, as it resists many more Ca2+ additions before opening in a highly active organ such as the heart than in an organ that seldom suffers Ca2+ loading, such as the liver.

Navigating spatio-temporal microbiome dynamics: Environmental factors and trace elements shape the symbiont community of an invasive marine species.

The proliferation of marine invasive species is a mounting concern. While the role of microbial communities in invasive ascidian species is recognized, the role of seasonal shifts in microbiome composition remains largely unexplored. We sampled five individuals of the invasive ascidian Styela plicata quarterly from January 2020 to October 2021 in two harbours, examining gills, tunics, and surrounding water. By analysing Amplicon Sequence Variants (ASVs) and seawater trace elements, we found that compartment (seawater, tunic, or gills) was the primary differentiating factor, followed by harbour. Clear seasonal patterns were evident in seawater bacteria, less so in gills, and absent in tunics. We identified compartment-specific bacteria, as well as seasonal indicator ASVs and ASVs correlated with trace element concentrations. Among these bacteria, we found that Endozoicomonas, Hepatoplasma and Rhodobacteraceae species had reported functions which might be necessary for overcoming seasonality and trace element shifts. This study contributes to understanding microbiome dynamics in invasive holobiont systems, and the patterns found indicate a potential role in adaptation and invasiveness.

Understanding the multi-functionality and tissue-specificity of decellularized dental pulp matrix hydrogels for endodontic regeneration.

Dental pulp is the only soft tissue in the tooth which plays a crucial role in maintaining intrinsic multi-functional behaviors of the dentin-pulp complex. Nevertheless, the restoration of fully functional pulps after pulpitis or pulp necrosis, termed endodontic regeneration, remained a major challenge for decades. Therefore, a bioactive and in-situ injectable biomaterial is highly desired for tissue-engineered pulp regeneration. Herein, a decellularized matrix hydrogel derived from porcine dental pulps (pDDPM-G) was prepared and characterized through systematic comparison against the porcine decellularized nerve matrix hydrogel (pDNM-G). The pDDPM-G not only exhibited superior capabilities in facilitating multi-directional differentiation of dental pulp stem cells (DPSCs) during 3D culture, but also promoted regeneration of pulp-like tissues after DPSCs encapsulation and transplantation. Further comparative proteomic and transcriptome analyses revealed the differential compositions and potential mechanisms that endow the pDDPM-G with highly tissue-specific properties. Finally, it was realized that the abundant tenascin C (TNC) in pDDPM served as key factor responsible for the activation of Notch signaling cascades and promoted DPSCs odontoblastic differentiation. Overall, it is believed that pDDPM-G is a sort of multi-functional and tissue-specific hydrogel-based material that holds great promise in endodontic regeneration and clinical translation. STATEMENT OF SIGNIFICANCE: Functional hydrogel-based biomaterials are highly desirable for endodontic regeneration treatments. Decellularized extracellular matrix (dECM) preserves most extracellular matrix components of its native tissue, exhibiting unique advantages in promoting tissue regeneration and functional restoration. In this study, we prepared a porcine dental pulp-derived dECM hydrogel (pDDPM-G), which exhibited superior performance in promoting odontogenesis, angiogenesis, and neurogenesis of the regenerating pulp-like tissue, further showed its tissue-specificity compared to the peripheral nerve-derived dECM hydrogel. In-depth proteomic and transcriptomic analyses revealed that the activation of tenascin C-Notch axis played an important role in facilitating odontogenic regeneration. This biomaterial-based study validated the great potential of the dental pulp-specific pDDPM-G for clinical applications, and provides a springboard for research strategies in ECM-related regenerative medicine.

Molecular adaptations in response to exercise training are associated with tissue-specific transcriptomic and epigenomic signatures.

Regular exercise has many physical and brain health benefits, yet the molecular mechanisms mediating exercise effects across tissues remain poorly understood. Here we analyzed 400 high-quality DNA methylation, ATAC-seq, and RNA-seq datasets from eight tissues from control and endurance exercise-trained (EET) rats. Integration of baseline datasets mapped the gene location dependence of epigenetic control features and identified differing regulatory landscapes in each tissue. The transcriptional responses to 8 weeks of EET showed little overlap across tissues and predominantly comprised tissue-type enriched genes. We identified sex differences in the transcriptomic and epigenomic changes induced by EET. However, the sex-biased gene responses were linked to shared signaling pathways. We found that many G protein-coupled receptor-encoding genes are regulated by EET, suggesting a role for these receptors in mediating the molecular adaptations to training across tissues. Our findings provide new insights into the mechanisms underlying EET-induced health benefits across organs.

Comparative transcriptomic analysis delineates adaptation strategies of Rana kukunoris toward cold stress on the Qinghai-Tibet Plateau.

BACKGROUND: Cold hardiness is fundamental for amphibians to survive during the extremely cold winter on the Qinghai-Tibet plateau. Exploring the gene regulation mechanism of freezing-tolerant Rana kukunoris could help us to understand how the frogs survive in winter. RESULTS: Transcriptome of liver and muscle of R. kukunoris collected in hibernation and spring were assisted by single molecule real-time (SMRT) sequencing technology. A total of 10,062 unigenes of R. kukunoris were obtained, and 9,924 coding sequences (CDS) were successfully annotated. Our examination of the mRNA response to whole body freezing and recover in the frogs revealed key genes concerning underlying antifreeze proteins and cryoprotectants (glucose and urea). Functional pathway analyses revealed differential regulated pathways of ribosome, energy supply, and protein metabolism which displayed a freeze-induced response and damage recover. Genes related to energy supply in the muscle of winter frogs were up-regulated compared with the muscle of spring frogs. The liver of hibernating frogs maintained modest levels of protein synthesis in the winter. In contrast, the liver underwent intensive high levels of protein synthesis and lipid catabolism to produce substantial quantity of fresh proteins and energy in spring. Differences between hibernation and spring were smaller than that between tissues, yet the physiological traits of hibernation were nevertheless passed down to active state in spring. CONCLUSIONS: Based on our comparative transcriptomic analyses, we revealed the likely adaptive mechanisms of R. kukunoris. Ultimately, our study expands genetic resources for the freezing-tolerant frogs.

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

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

Enriched atlas of lncRNA and protein-coding genes for the GRCg7b chicken assembly and its functional annotation across 47 tissues.

Gene atlases for livestock are steadily improving thanks to new genome assemblies and new expression data improving the gene annotation. However, gene content varies across databases due to differences in RNA sequencing data and bioinformatics pipelines, especially for long non-coding RNAs (lncRNAs) which have higher tissue and developmental specificity and are harder to consistently identify compared to protein coding genes (PCGs). As done previously in 2020 for chicken assemblies galgal5 and GRCg6a, we provide a new gene atlas, lncRNA-enriched, for the latest GRCg7b chicken assembly, integrating "NCBI RefSeq", "EMBL-EBI Ensembl/GENCODE" reference annotations and other resources such as FAANG and NONCODE. As a result, the number of PCGs increases from 18,022 (RefSeq) and 17,007 (Ensembl) to 24,102, and that of lncRNAs from 5789 (RefSeq) and 11,944 (Ensembl) to 44,428. Using 1400 public RNA-seq transcriptome representing 47 tissues, we provided expression evidence for 35,257 (79%) lncRNAs and 22,468 (93%) PCGs, supporting the relevance of this atlas. Further characterization including tissue-specificity, sex-differential expression and gene configurations are provided. We also identified conserved miRNA-hosting genes with human counterparts, suggesting common function. The annotated atlas is available at gega.sigenae.org.

ApCPS2 contributes to medicinal diterpenoid biosynthesis and defense against insect herbivore in Andrographis paniculata.

Kalmegh (Andrographis paniculata) spatiotemporally produces medicinally-important ent-labdane-related diterpenoids (ent-LRDs); andrographolide (AD), 14-deoxy-11,12-didehydroandrographolide (DDAD), neoandrographolide (NAD). ApCPS1 and ApCPS2, the ent-copalyl pyrophosphate (ent-CPP)-producing class II diterpene synthases (diTPSs) were identified, but their contributions to ent-CPP precursor supply for ent-LRD biosynthesis were not well understood. Here, we characterized ApCPS4, an additional ent-CPP-forming diTPS. Further, we elucidated in planta function of the ent-CPP-producing diTPSs (ApCPS1,2,4) by integrating transcript-metabolite co-profiles, biochemical analysis and gene functional characterization. ApCPS1,2,4 localized to the plastids, where diterpenoid biosynthesis occurs in plants, but ApCPS1,2,4 transcript expression patterns and ent-LRD contents revealed a strong correlation of ApCPS2 expression and ent-LRD accumulation in kalmegh. ApCPS1,2,4 upstream sequences differentially activated ß-glucuronidase (GUS) in Arabidopsis and transiently-transformed kalmegh. Similar to higher expression of ApCPS1 in kalmegh stem, ApCPS1 upstream sequence activated GUS in stem/hypocotyl of Arabidopsis and kalmegh. However, ApCPS2,4 upstream sequences weakly activated GUS expression in Arabidopsis, which was not well correlated with ApCPS2,4 transcript expression in kalmegh tissues. Whereas, ApCPS2,4 upstream sequences could activate GUS expression at a considerable level in kalmegh leaf and roots/calyx, respectively, suggesting the involvement of transcriptional regulator(s) of ApCPS2,4 that might participate in kalmegh-specific diterpenoid pathway. Interestingly, ApCPS2-silenced kalmegh showed a drastic reduction in AD, DDAD and NAD contents and compromised defense against insect herbivore Spodoptera litura. However, ent-LRD contents and herbivore defense in ApCPS1 or ApCPS4-silenced plants remained largely unaltered. Overall, these results suggested an important role of ApCPS2 in producing ent-CPP for medicinal ent-LRD biosynthesis and defense against insect herbivore.

Epigenetic analyses in forensic medicine: future and challenges.

The possibility of using epigenetics in forensic investigation has gradually risen over the last few years. Epigenetic changes with their dynamic nature can either be inherited or accumulated throughout a lifetime and be reversible, prompting investigation of their use across various fields. In forensic sciences, multiple applications have been proposed, such as the discrimination of monozygotic twins, identifying the source of a biological trace left at a crime scene, age prediction, determination of body fluids and tissues, human behavior association, wound healing progression, and determination of the post-mortem interval (PMI). Despite all these applications, not all the studies considered the impact of PMI and post-sampling effects on the epigenetic modifications and the tissue-specificity of the epigenetic marks.This review aims to highlight the substantial forensic significance that epigenetics could support in various forensic investigations. First, basic concepts in epigenetics, describing the main epigenetic modifications and their functions, in particular, DNA methylation, histone modifications, and non-coding RNA, with a particular focus on forensic applications, were covered. For each epigenetic marker, post-mortem stability and tissue-specificity, factors that should be carefully considered in the study of epigenetic biomarkers in the forensic context, have been discussed. The advantages and limitations of using post-mortem tissues have been also addressed, proposing directions for these innovative strategies to analyze forensic specimens.

Genetic effects of sequence-conserved enhancer-like elements on human complex traits.

BACKGROUND: The vast majority of findings from human genome-wide association studies (GWAS) map to non-coding sequences, complicating their mechanistic interpretations and clinical translations. Non-coding sequences that are evolutionarily conserved and biochemically active could offer clues to the mechanisms underpinning GWAS discoveries. However, genetic effects of such sequences have not been systematically examined across a wide range of human tissues and traits, hampering progress to fully understand regulatory causes of human complex traits. RESULTS: Here we develop a simple yet effective strategy to identify functional elements exhibiting high levels of human-mouse sequence conservation and enhancer-like biochemical activity, which scales well to 313 epigenomic datasets across 106 human tissues and cell types. Combined with 468 GWAS of European (EUR) and East Asian (EAS) ancestries, these elements show tissue-specific enrichments of heritability and causal variants for many traits, which are significantly stronger than enrichments based on enhancers without sequence conservation. These elements also help prioritize candidate genes that are functionally relevant to body mass index (BMI) and schizophrenia but were not reported in previous GWAS with large sample sizes. CONCLUSIONS: Our findings provide a comprehensive assessment of how sequence-conserved enhancer-like elements affect complex traits in diverse tissues and demonstrate a generalizable strategy of integrating evolutionary and biochemical data to elucidate human disease genetics.

Recent advances in understanding microRNA function and regulation in C. elegans.

MicroRNAs (miRNAs) were first discovered in C. elegans as essential post-transcriptional regulators of gene expression. Since their initial discovery, miRNAs have been implicated in numerous areas of physiology and disease in all animals examined. In recent years, the C. elegans model continues to contribute important advances to all areas of miRNA research. Technological advances in tissue-specific miRNA profiling and genome editing have driven breakthroughs in understanding biological functions of miRNAs, mechanism of miRNA action, and regulation of miRNAs. In this review, we highlight these new C. elegans findings from the past five to seven years.

Comparative RNA-Seq Analysis Revealed Tissue-Specific Splicing Variations during the Generation of the PDX Model.

Tissue-specific gene expression generates fundamental differences in the function of each tissue and affects the characteristics of the tumors that are created as a result. However, it is unclear how much the tissue specificity is conserved during grafting of the primary tumor into an immune-compromised mouse model. Here, we performed a comparative RNA-seq analysis of four different primary-patient derived xenograft (PDX) tumors. The analysis revealed a conserved RNA biotype distribution of primary-PDX pairs, as revealed by previous works. Interestingly, we detected significant changes in the splicing pattern of PDX, which was mainly comprised of skipped exons. This was confirmed by splicing variant-specific RT-PCR analysis. On the other hand, the correlation analysis for the tissue-specific genes indicated overall strong positive correlations between the primary and PDX tumor pairs, with the exception of gastric cancer cases, which showed an inverse correlation. These data propose a tissue-specific change in splicing events during PDX formation as a variable factor that affects primary-PDX integrity.

Comparison of Alternative Splicing Landscapes Revealed by Long-Read Sequencing in Hepatocyte-Derived HepG2 and Huh7 Cultured Cells and Human Liver Tissue.

The long-read RNA sequencing developed by Oxford Nanopore Technologies provides a direct quantification of transcript isoforms, thereby making it possible to present alternative splicing (AS) profiles as arrays of single splice variants with different abundances. Additionally, AS profiles can be presented as arrays of genes characterized by the degree of alternative splicing (the DAS-the number of detected splice variants per gene). Here, we successfully utilized the DAS to reveal biological pathways influenced by the alterations in AS in human liver tissue and the hepatocyte-derived malignant cell lines HepG2 and Huh7, thus employing the mathematical algorithm of gene set enrichment analysis. Furthermore, analysis of the AS profiles as abundances of single splice variants by using the graded tissue specificity index τ provided the selection of the groups of genes expressing particular splice variants specifically in liver tissue, HepG2 cells, and Huh7 cells. The majority of these splice variants were translated into proteins products and appeal to be in focus regarding further insights into the mechanisms underlying cell malignization. The used metrics are intrinsically suitable for transcriptome-wide AS profiling using long-read sequencing.