Revealing the History and Mystery of RNA-Seq.

Advances in RNA-sequencing technologies have led to the development of intriguing experimental setups, a massive accumulation of data, and high demand for tools to analyze it. To answer this demand, computational scientists have developed a myriad of data analysis pipelines, but it is less often considered what the most appropriate one is. The RNA-sequencing data analysis pipeline can be divided into three major parts: data pre-processing, followed by the main and downstream analyses. Here, we present an overview of the tools used in both the bulk RNA-seq and at the single-cell level, with a particular focus on alternative splicing and active RNA synthesis analysis. A crucial part of data pre-processing is quality control, which defines the necessity of the next steps; adapter removal, trimming, and filtering. After pre-processing, the data are finally analyzed using a variety of tools: differential gene expression, alternative splicing, and assessment of active synthesis, the latter requiring dedicated sample preparation. In brief, we describe the commonly used tools in the sample preparation and analysis of RNA-seq data.

TRA2B Gene Splice Variant Linked to Seizures and Neurodevelopmental Delay: A Second Case Study.

In this study, we report a novel splice variant in the TRA2B gene identified in a patient presenting with seizures and neurodevelopmental delay. This paper represents the second investigation of pathogenic variants in the TRA2B gene in humans, reaffirming the conclusions of the initial study and underscoring the importance of this research. Comprehensive genetic testing, including whole genome sequencing, Sanger sequencing, and mRNA analysis, was performed on the proband and her parents. The proband harbored a de novo c.170+1G>A variant in the RS1 domain of Tra2ß, which was confirmed to be pathogenic through mRNA analysis, resulting in exon 2 deletion and a frameshift (p.Glu13Valfs*2). The clinical presentation of the patient was consistent with phenotypes described in one of the previous studies. These findings contribute to the dissemination and reinforcement of prior discoveries in the context of TRA2B-related syndrome and highlight the need for further investigation into the functional consequences and underlying pathogenic mechanisms associated with TRA2B mutations.

Transcriptome Analysis Reveals Multiple Genes and Complex Hormonal-Mediated Interactions with PEG during Adventitious Root Formation in Apple.

Adventitious root (AR) formation is a bottleneck for the mass propagation of apple rootstocks, and water stress severely restricts it. Different hormones and sugar signaling pathways in apple clones determine AR formation under water stress, but these are not entirely understood. To identify them, GL-3 stem cuttings were cultured on polyethylene glycol (PEG) treatment. The AR formation was dramatically decreased compared with the PEG-free control (CK) cuttings by increasing the endogenous contents of abscisic acid (ABA), zeatin riboside (ZR), and methyl jasmonate (JA-me) and reducing the indole-3-acetic acid (IAA) and gibberellic acid 3 (GA3) contents. We performed a transcriptomic analysis to identify the responses behind the phenotype. A total of 3204 differentially expressed genes (DEGs) were identified between CK and PEG, with 1702 upregulated and 1502 downregulated genes. Investigation revealed that approximately 312 DEGs were strongly enriched in hormone signaling, sugar metabolism, root development, and cell cycle-related pathways. Thus, they were selected for their possible involvement in adventitious rooting. However, the higher accumulation of ABA, ZR, and JA-me contents and the upregulation of their related genes, as well as the downregulation of sugar metabolism-related genes, lead to the inhibition of ARs. These results indicate that AR formation is a complicated biological process chiefly influenced by multiple hormonal signaling pathways and sugar metabolism. This is the first study to demonstrate how PEG inhibits AR formation in apple plants.

Improving clinical interpretation of five KRIT1 and PDCD10 intronic variants.

Cerebral cavernous malformation (CCM) is a vascular malformation of the central nervous system which may occur sporadically or segregate within families due to heterozygous variants in KRIT1/CCM1, MGC4607/CCM2 or PDCD10/CCM3. Intronic variants are not uncommon in familial CCM, but their clinical interpretation is often hampered by insufficient data supporting in silico predictions. Here, the mRNA analysis for two intronic unpublished variants (KRIT1 c.1147-7 T > G and PDCD10 c.395 + 2 T > G) and three previously published variants in KRIT1 but without data supporting their effects was carried out. This study demonstrated that all variants can induce a frameshift with the lack of residues located in the C-terminal regions and involved in protein-protein complex formation, which is essential for vascular homeostasis. These results support the introduction of mRNA analysis in the diagnostic pathway of familial CCM and expand the knowledge of abnormal splicing patterning in this disorder.

RT-PCR analysis of mRNA revealed the splice-altering effect of rare intronic variants in monogenic disorders.

BACKGROUND: Variants perturbing the normal splicing of pre-mRNA can lead to human diseases. The splice-altering effect and eventual consequence on gene function was sometimes uncertain and hinders a definitive molecular diagnosis. METHODS: The impact of four rare intronic variants on splicing was analyzed through reverse transcription - polymerase chain reaction (RT-PCR) analysis of mRNA derived from the peripheral blood of patients. The results were compared with in-silico prediction. Potential implication on molecular diagnosis was discussed. RESULTS: Four rare intronic variants of SLC9A6, DLG3, GAA, and OCRL were identified in patients with suspected disorders, respectively. Although these four variants were all predicted to alter splicing by in-silico tools, RT-PCR analysis of mRNA derived from peripheral blood showed these variants affected splicing in different ways: c.899+3_899+6del of SLC9A6 resulted in one-exon skipping and an out-of-frame transcript; c.905-2A > G of DLG3 resulted in a mix of in-frame transcripts; c.1195-11T > A of GAA resulted in the in-frame insertion of nine nucleotides; c.723-2A > C of OCRL resulted in one-exon skipping and in-frame deletion of 102 nucleotides. The consequence revealed by mRNA analysis is essential for accurate interpretation of pathogenicity. CONCLUSION: Four intronic variants all caused aberrant mRNA splicing. For intronic variants with uncertain impact on splicing, mRNA analysis is helpful for ascertainment of alternative splicing and accurate interpretation of pathogenicity.

Integrated miRNA/mRNA Counter-Expression Analysis Highlights Oxidative Stress-Related Genes CCR7 and FOXO1 as Blood Markers of Coronary Arterial Disease.

Our interest in the mechanisms of atherosclerosis progression (ATHp) has led to the recent identification of 13 miRNAs and 1285 mRNAs whose expression was altered during ATHp. Here, we deepen the functional relationship among these 13 miRNAs and genes associated to oxidative stress, a crucial step in the onset and progression of vascular disease. We first compiled a list of genes associated to the response to oxidative stress (Oxstress genes) by performing a reverse Gene Ontology analysis (rGO, from the GO terms to the genes) with the GO terms GO0006979, GO1902882, GO1902883 and GO1902884, which included a total of 417 unique Oxstress genes. Next, we identified 108 putative targets of the 13 miRNAs among these unique Oxstress genes, which were validated by an integrated miRNA/mRNA counter-expression analysis with the 1285 mRNAs that yielded 14 genes, Map2k1, Mapk1, Mapk9, Dapk1, Atp2a2, Gata4, Fos, Egfr, Foxo1, Ccr7, Vkorc1l1, Rnf7, Kcnh3, and Mgat3. GO enrichment analysis and a protein-protein-interaction network analysis (PPI) identified most of the validated Oxstress transcripts as components of signaling pathways, highlighting a role for MAP signaling in ATHp. Lastly, expression of these Oxstress transcripts was measured in PBMCs from patients suffering severe coronary artery disease, a serious consequence of ATHp. This allowed the identification of FOXO1 and CCR7 as blood markers downregulated in CAD. These results are discussed in the context of the interaction of the Oxstress transcripts with the ATHp-associated miRNAs.

A novel de novo splicing mutation c.1444-2A>T in the TSC2 gene causes exon skipping and premature termination in a patient with tuberous sclerosis syndrome.

Tuberous sclerosis complex (TSC) syndrome is a neurocutaneous syndrome that affects the brain, skin, and kidneys that has an adverse impact on the patient's health and quality of life. There have been several recent advances that elucidate the genetic complex of this disorder that will help understand the basic neurobiology of this disorder. We report a Tunisian patient with clinical manifestations of TSC syndrome. We investigated the causative molecular defect in this patient using PCR followed by direct sequencing. Subsequently, in silico studies and mRNA analysis were performed to study the pathogenicity of the new variation found in the TSC2. Bioinformatics tools predicted that the novel mutation c.1444-2A>T have pathogenic effects on splicing machinery. RT-PCR followed by sequencing revealed that the mutation c.1444-2A>T generates two aberrant transcripts. The first, with exon 15 skipping, is responsible for the loss of 52 amino acids, which causes the production of an aberrant protein isoform. The second, with the inclusion of 122 nucleotides of intron 14, is responsible for the creation of new premature termination codons (TGA), which causes the production of a truncated TSC2 protein. This study highlighted the clinical features of a Tunisian patient with TSC syndrome and revealed a splicing mutation c.1444-2A>T within intron 14 of TSC2 gene, which is present for the first time using Sanger sequencing approach, as a disease-causing mutation in a Tunisian patient with TSC syndrome.

Deep intronic variant c.5999-277G>A of F8 gene may be a hot spot mutation for mild hemophilia A patients without mutation in exonic DNA.

BACKGROUND: In 10%-18% of mild-type hemophilia A (HA) patients, mutations cannot be found by routine DNA analysis. OBJECTIVE: We aimed to identify the genetic defects by mRNA analysis of F8 gene in mild HA patients without mutation in exonic DNA. PATIENTS AND METHODS: From 2006 to 2016, we identified F8 exon mutations in 39 of 49 mild HA patients using routine genetic testing. We then evaluated the 10 remaining patients from six unrelated families without exonic DNA mutation by performing cDNA sequence analysis. RESULTS: Nine of the 10 (90%) patients were confirmed to have F8 gene mutation. Eight patients from four unrelated families were notably found to have presence of an aberrant 675-bp fragment. Sequencing of this fragment showed that there were two separate new alternative splicing exons of 35 bp and 55 bp within intron 18, which formed a 90-bp insertion between exon 18 and exon 19 (E18ins90bpE19) in the mRNA. Based on direct sequencing, this alternative splicing transcript appears to have resulted from deep intronic variant c.5999-277G>A of intron 18. CONCLUSIONS: Our study suggests that deep intronic variant of c.5999-277G>A may be a hot spot mutation for mild hemophilia patients without mutation in exonic DNA.

CTNS mRNA molecular analysis revealed a novel mutation in a child with infantile nephropathic cystinosis: a case report.

BACKGROUND: Cystinosis is an autosomal recessive lysosomal storage disorder characterized by accumulation of cystine in lysosomes throughout the body. Cystinosis is caused by mutations in the CTNS gene that encodes the lysosomal cystine carrier protein cystinosin. CTNS mutations result in either complete absence or reduced cystine transporting function of the protein. The diagnosis of nephropathic cystinosis is generally based on measuring leukocyte cystine level, demonstration of corneal cystine crystals by the slit lamp examination and confirmed by genetic analysis of the CTNS gene. CASE PRESENTATION: A boy born to consanguineous Caucasian parents had the characteristic clinical features of the infantile nephropathic cystinosis including renal Fanconi syndrome (polydipsia/polyuria, metabolic acidosis, hypokalemia, hypophosphatemia, low molecular weight proteinuria, glycosuria, cystine crystals in the cornea) and elevated WBC cystine levels. Initially we performed RFLP analysis of the common in the Northern European population 57-kb deletion of proband's DNA, then a direct Sanger sequencing which revealed no mutations in the coding part of the CTNS gene. To confirm the diagnosis we performed RT-PCR analysis of total RNA obtained from patient-derived fibroblasts in combination with cDNA sequencing. This revealed the skipping of exon 4 and exon 5 in the CTNS in our patient. Therefore, we detected a novel 9-kb homozygous deletion in the CTNS gene at genomic DNA level, spanning region from intron 3 to intron 5. In order to identify the inheritance pattern of the deletion we analyzed DNA of proband's mother and father. Both parents were found to be heterozygous carriers of the CTNS mutation. CONCLUSIONS: Analysis of CTNS gene transcript allowed to identify a large homozygous deletion in the patient with infantile nephropathic cystinosis. Mutational detection at RNA level may be an efficient tool to establish the genetic defect in some cystinosis patients.

Integrated miRNA-mRNA analysis reveals regulatory pathways underlying the curly fleece trait in Chinese tan sheep.

BACKGROUND: Tan sheep is an indigenous Chinese breed well known for its beautiful curly fleece. One prominent breed characteristic of this sheep breed is that the degree of curliness differs markedly between lambs and adults, but the molecular mechanisms regulating the shift are still not well understood. In this study, we identified 49 differentially expressed (DE) microRNAs (miRNAs) between Tan sheep at the two stages through miRNA-seq, and combined the data with that in our earlier Suppression Subtractive Hybridization cDNA (SSH) library study to elucidate the mechanisms underlying curly fleece formation. RESULTS: Thirty-six potential miRNA-mRNA target pairs were identified using computational methods, including 25 DE miRNAs and 10 DE genes involved in the MAPK signaling pathway, steroid biosynthesis and metabolic pathways. With the differential expressions between lambs and adults confirmed by qRT-PCR, some miRNAs were already annotated in the genome, but some were novel miRNAs. Inhibition of KRT83 expression by miR-432 was confirmed by both gene knockdown with siRNA and overexpression, which was consistent with the miRNAs and targets prediction results. CONCLUSION: Our study represents the comprehensive analysis of mRNA and miRNA in Tan sheep and offers detailed insight into the development of curly fleece as well as the potential mechanisms controlling curly hair formation in humans.

Functional mRNA analysis reveals aberrant splicing caused by novel intronic mutation in WDR45 in NBIA patient.

WDR45 gene-associated neurodegeneration with brain iron accumulation (NBIA), referred to as beta-propeller protein-associated neurodegeneration (BPAN), is a rare disorder that presents with a very nonspecific clinical phenotype in children constituting global developmental delay. This case report illustrates the power of a combination of trio exome sequencing, in silico splicing analysis, and mRNA analysis to provide sufficient evidence for pathogenicity of a relatively intronic variant in WDR45, and in so doing, find a genetic diagnosis for a 6-year-old patient with developmental delay and seizures, a diagnosis which may otherwise have only been found once the characteristic MRI patterns of the disease became more obvious in young adulthood.

Human single follicle growth in vitro from cryopreserved ovarian tissue after slow freezing or vitrification.

STUDY QUESTION: What is the effect of human ovarian tissue cryopreservation on single follicular development in vitro? SUMMARY ANSWER: Vitrification had a greater negative effect on growth and gene expression of human ovarian follicles when compared with fresh follicles. WHAT IS KNOWN ALREADY: For human ovarian cortex cryopreservation, the conventional option is slow freezing while more recently vitrification has been demonstrated to maintain good quality and function of ovarian tissues. STUDY DESIGN, SIZE, DURATION: Ovarian tissues were collected from 11 patients. For every patient, the ovarian cortex was divided into three samples: Fresh, slow-rate freezing (Slow) and vitrification (Vit). Tissue histology was performed and follicles were isolated for single-cell mRNA analysis and in vitro culture (IVC) in 1% alginate for 8 days. PARTICIPANTS/MATERIALS, SETTING, METHODS: Follicle morphology was assessed with hematoxylin-eosin analysis. Follicles were individually embedded in alginate (1% w/v) and cultured in vitro for 8 days. Follicle survival and growth were assessed by microscopy. Follicle viability was observed after Calcein-AM and ethidium homodimer-I (Ca-AM/EthD-I) staining. Expression of genes, including GDF9 (growth differentiation factor 9), BMP15 (bone morphogenetic protein 15) and ZP3 (zona pellucida glycoprotein 3) in oocytes and AMH (anti-Mullerian hormone), FSHR (FSH receptor), CYP11A (cholesterol side-chain cleavage cytochrome P450) and STAR (steroidogenic acute regulatory protein) in GCs, was evaluated by single-cell mRNA analysis. MAIN RESULTS AND THE ROLE OF CHANCE: A total of 129 follicles were separated from ovarian cortex (Fresh n = 44; Slow n = 40; Vit n = 45). The percentage of damaged oocytes and granulosa cells was significantly higher in both the Slow and Vit groups, as compared with Fresh control (P< 0.05). The growth of follicles in vitro was significantly delayed in the Vit group compared with the Fresh group (P< 0.05). Both slow freezing (P< 0.05) and vitrification (P< 0.05) down-regulated the mRNA levels of ZP3 and CYP11A compared with Fresh group, while there was no significant difference between the Slow and Vit groups (P> 0.05). Vitrification also down-regulates AMH mRNA levels compared with Fresh group (P< 0.05). LIMITATIONS, REASONS FOR CAUTION: Only short-term IVC studies (8 days) are reported. Further study should be performed to examine and improve follicular development in a long-term culture system after cryopreservation. WIDER IMPLICATIONS OF THE FINDINGS: This is the first comparison of gene expression and growth of single human ovarian follicles in vitro after either slow freezing or vitrification. With the decreased gene expression and growth during IVC, damage by cryopreservation still exists and needs to be minimized during the long-term IVC of follicles in the future for eventual clinical application. STUDY FUNDING/COMPETING INTERESTS: This work was supported by the National Natural Science Foundation of China (31230047, 81571386, 81471508, 31429004 and 81501247), National Natural Science Foundation of Beijing (7142166) and Mega-projects of Science Research for the 12th five-year plan (2012ba132b05). There are no conflicts of interest to declare.

A large cohort of ß(+)-thalassemia in Thailand: molecular, hematological and diagnostic considerations.

We report the molecular and hematological characteristics associated with a large cohort of ß(+)-thalassemia in Thailand. Study was done on 21,068 unrelated subjects referred to our center in northeast Thailand for hemoglobinopathies investigation. Among 21,068 subjects, 2637 (12.5%) were found to carry ß-thalassemia. Of these 2637 cases, 705 (26.7%) carried ß(+)-thalassemia with eight different mutations including 6 promoter mutations; NT-28 (A-G), NT-31 (A-G), NT-50 (G-A), NT-86 (C-G), NT-87 (C-A) and NT-90 (C-T) and two missense mutations; Hb Malay (codon 19; AAC-AGC) and Hb Dhonburi (codon 126; GTG-GGG). Hematological features of carriers with these ß(+)-thalassemia (n=528) were compared with those with ß(0)-thalassemia (n=309). Data for Hb E-ß(+)-thalassemia (n=177) were also presented along with Hb E-ß(0)-thalassemia in our series (n=94). All patients with Hb E-ß(+)-thalassemia were associated with mild thalassemia intermedia phenotypes. Most of the ß(+)-thalassemia carriers had elevated Hb A2 and mild hypochromic microcytosis, some demonstrated borderline MCV and MCH values which, could compromise carrier screening. Analysis of α/ß-globin mRNA ratio in representative cases with normal, Hb E trait, ß(+)-thalassemia trait, Hb Dhonburi trait and ß(0)-thalassemia trait demonstrated the average values of 1.1, 1.7, 2.1, 1.7 and 3.1, respectively which is helpful in identification and differentiation of the cases.

Characterization of an apparently synonymous F5 mutation causing aberrant splicing and factor V deficiency.

Coagulation factor V (FV) deficiency is a rare autosomal recessive bleeding disorder. We investigated a patient with severe FV deficiency (FV:C < 3%) and moderate bleeding symptoms. Thrombin generation experiments showed residual FV expression in the patient's plasma, which was quantified as 0.7 ± 0.3% by a sensitive prothrombinase-based assay. F5 gene sequencing identified a novel missense mutation in exon 4 (c.578G>C, p.Cys193Ser), predicting the abolition of a conserved disulphide bridge, and an apparently synonymous variant in exon 8 (c.1281C>G). The observation that half of the patient's F5 mRNA lacked the last 18 nucleotides of exon 8 prompted us to re-evaluate the c.1281C>G variant for its possible effects on splicing. Bioinformatics sequence analysis predicted that this transversion would activate a cryptic donor splice site and abolish an exonic splicing enhancer. Characterization in a F5 minigene model confirmed that the c.1281C>G variant was responsible for the patient's splicing defect, which could be partially corrected by a mutation-specific morpholino antisense oligonucleotide. The aberrantly spliced F5 mRNA, whose stability was similar to that of the normal mRNA, encoded a putative FV mutant lacking amino acids 427-432. Expression in COS-1 cells indicated that the mutant protein is poorly secreted and not functional. In conclusion, the c.1281C>G mutation, which was predicted to be translationally silent and hence neutral, causes FV deficiency by impairing pre-mRNA splicing. This finding underscores the importance of cDNA analysis for the correct assessment of exonic mutations.

Splicing defects caused by exonic mutations in PKD1 as a new mechanism of pathogenesis in autosomal dominant polycystic kidney disease.

The correct splicing of precursor-mRNA depends on the actual splice sites plus exonic and intronic regulatory elements recognized by the splicing machinery. Surprisingly, an increasing number of examples reveal that exonic mutations disrupt the binding of splicing factors to these sequences or generate new splice sites or regulatory elements, causing disease. This contradicts the general assumption that missense mutations disrupt protein function and that synonymous mutations are merely polymorphisms. Autosomal dominant polycystic kidney disease (ADPKD) is a common inherited disorder caused mainly by mutations in the PKD1 gene. Recently, we analyzed a substantial number of PKD1 missense or synonymous mutations to further characterize their consequences on pre-mRNA splicing. Our results showed that one missense and 2 synonymous mutations induce significant defects in pre-mRNA splicing. Thus, it appears that aberrant splicing as a result of exonic mutations is a previously unrecognized cause of ADPKD.

A high-throughput protocol for message RNA quantification using RNA dot-blots.

This study develops a method to rapidly measure the relative abundance of mRNA in total RNA samples using a dot-blotting technique and biotin-labeled detection probes that recognize the polyadenylate tail on mRNA. We demonstrate the effectiveness of this technique by determining the relative total amounts of mRNA in three tissues of turtles (Trachemys scripta elegans) exposed to normoxic versus anoxic conditions. The data emphasize the usefulness of the method for the simple and rapid analysis of relative total mRNA levels for a variety of comparison purposes.

MicroRNA expression profile of gastric cancer stem cells in the MKN-45 cancer cell line.

MicroRNAs (miRNAs) play important roles in post-transcriptional gene silencing of target messenger RNAs, which are involved in virtually all biological processes. Previously, we have demonstrated that spheroid body-forming cells from the MKN-45 cancer cell line possessed gastric cancer stem cell (CSC) properties. In this study, we aimed to determine the miRNA profile of the gastric CSCs and to explore the role of miRNAs in gastric CSCs. Human miRNA microarrays, which contain probes specific for 1887 human miRNAs were used to determine the expression profiles of the gastric CSCs. A total of 182 miRNAs with a more than 2-fold change were identified to be differentially expressed between the spheroid body-forming cells and the parental cells. Of these miRNAs, 9 miRNAs were over-expressed in the spheroid body-forming cells, while the other 173 miRNAs were all under-expressed, indicating that the role of most miRNAs in human gastric CSCs may be tumor suppressors. The results of microarray analysis were validated by quantitative real-time polymerase chain reaction, and the consistence rate is 70% (7 out of 10). The target genes for the validated miRNAs were predicted by using three online software programs, miRanda, PicTar, and TargetScan. Most of the potential targets of the miRNAs were relevant to the regulation of actin cytoskeleton, focal adhesion, extracellular matrix-receptor interaction, and the pathways in cancer. Especially, several genes are associated with some pivotal signaling pathways of the 'stem cell genes'. Evaluating the characteristic miRNAs of the gastric CSCs may be a new method for studying gastric cancer and developing therapeutic strategies, which aimed at eradicating the subpopulation of CSCs in gastric cancer.

The Genetic Basis of the First Patient with Wiedemann-Rautenstrauch Syndrome in the Russian Federation.

Bi-allelic pathogenic variations within POLR3A have been associated with a spectrum of hereditary disorders. Among these, a less frequently observed condition is Wiedemann-Rautenstrauch syndrome (WRS), also known as neonatal progeroid syndrome. This syndrome typically manifests neonatally and is characterized by growth retardation, evident generalized lipodystrophy with distinctively localized fat accumulations, sparse scalp hair, and atypical facial features. Our objective was to elucidate the underlying molecular mechanisms of Wiedemann-Rautenstrauch syndrome (WRS). In this study, we present a clinical case of a 7-year-old female patient diagnosed with WRS. Utilizing whole-exome sequencing (WES), we identified a novel missense variant c.3677T>C (p.Leu1226Pro) in the POLR3A gene (NM_007055.4) alongside two cis intronic variants c.1909+22G>A and c.3337-11T>C. Via the analysis of mRNA derived from fibroblasts, we reconfirmed the splicing-affecting nature of the c.3337-11T>C variant. Furthermore, our investigation led to the reclassification of the c.3677T>C (p.Leu1226Pro) variant as a likely pathogenic variant. Therefore, this is the first case demonstrating the molecular genetics of a patient with Wiedemann-Rautenstrauch syndrome from the Russian Federation. A limited number of clinical cases have been documented until this moment; therefore, broadening the linkage between phenotype and molecular changes in the POLR3A gene will significantly contribute to the comprehensive understanding of the molecular basis of POLR3A-related disorders.

Large-Scale Single-Nucleus RNA Sequencing Compatible with Complex Archived Samples.

Transcriptome profiling at single-cell resolution allows us to identify and assess functional cell types and cellular states, including those within degenerating ocular tissues in retinitis pigmentosa. The technology is particularly valuable when studying tissues with high cellular heterogeneity, or when specific cell types are of interest. In this chapter, we introduce a detailed protocol of a medium-throughput single-nucleus RNA sequencing technique that utilizes frozen tissue as input sample. This protocol can be executed by any researcher with basic training in molecular biology techniques. With this protocol, a single experimenter can easily process two samples per day up to cDNA amplification, and library preparations can be done in batches of 8. Routinely we can obtain ~20 K nuclei per eye from 3 to 4 library preparations.

The co-inhibitory receptor TIGIT regulates NK cell function and is upregulated in human intrahepatic CD56bright NK cells.

The crosstalk between NK cells and their surrounding environment is enabled through activating and inhibitory receptors, which tightly control NK cell activity. The co-inhibitory receptor TIGIT decreases NK cell cytotoxicity and is involved in NK cell exhaustion, but has also been associated with liver regeneration, highlighting that the contribution of human intrahepatic CD56bright NK cells in regulating tissue homeostasis remains incompletely understood. A targeted single-cell mRNA analysis revealed distinct transcriptional differences between matched human peripheral blood and intrahepatic CD56bright NK cells. Multiparameter flow cytometry identified a cluster of intrahepatic NK cells with overlapping high expression of CD56, CD69, CXCR6, TIGIT and CD96. Intrahepatic CD56bright NK cells also expressed significantly higher protein surface levels of TIGIT, and significantly lower levels of DNAM-1 compared to matched peripheral blood CD56bright NK cells. TIGIT+ CD56bright NK cells showed diminished degranulation and TNF-α production following stimulation. Co-incubation of peripheral blood CD56bright NK cells with human hepatoma cells or primary human hepatocyte organoids resulted in migration of NK cells into hepatocyte organoids and upregulation of TIGIT and downregulation of DNAM-1 expression, in line with the phenotype of intrahepatic CD56bright NK cells. Intrahepatic CD56bright NK cells represent a transcriptionally, phenotypically, and functionally distinct population of NK cells that expresses higher levels of TIGIT and lower levels of DNAM-1 than matched peripheral blood CD56bright NK cells. Increased expression of inhibitory receptors by NK cells within the liver environment can contribute to tissue homeostasis and reduction of liver inflammation.