Sex, tissue, and mitochondrial interactions modify the transcriptional response to rapamycin in Drosophila.

BACKGROUND: Many common diseases exhibit uncontrolled mTOR signaling, prompting considerable interest in the therapeutic potential of mTOR inhibitors, such as rapamycin, to treat a range of conditions, including cancer, aging-related pathologies, and neurological disorders. Despite encouraging preclinical results, the success of mTOR interventions in the clinic has been limited by off-target side effects and dose-limiting toxicities. Improving clinical efficacy and mitigating side effects require a better understanding of the influence of key clinical factors, such as sex, tissue, and genomic background, on the outcomes of mTOR-targeting therapies. RESULTS: We assayed gene expression with and without rapamycin exposure across three distinct body parts (head, thorax, abdomen) of D. melanogaster flies, bearing either their native melanogaster mitochondrial genome or the mitochondrial genome from a related species, D. simulans. The fully factorial RNA-seq study design revealed a large number of genes that responded to the rapamycin treatment in a sex-dependent and tissue-dependent manner, and relatively few genes with the transcriptional response to rapamycin affected by the mitochondrial background. Reanalysis of an earlier study confirmed that mitochondria can have a temporal influence on rapamycin response. CONCLUSIONS: We found significant and wide-ranging effects of sex and body part, alongside a subtle, potentially time-dependent, influence of mitochondria on the transcriptional response to rapamycin. Our findings suggest a number of pathways that could be crucial for predicting potential side effects of mTOR inhibition in a particular sex or tissue. Further studies of the temporal response to rapamycin are necessary to elucidate the effects of the mitochondrial background on mTOR and its inhibition.

The multi-tissue gene expression and physiological responses of water deprived Peromyscus eremicus.

The harsh and dry conditions of desert environments have resulted in genomic adaptations, allowing for desert organisms to withstand prolonged drought, extreme temperatures, and limited food resources. Here, we present a comprehensive exploration of gene expression across five tissues (kidney, liver, lung, gastrointestinal tract, and hypothalamus) and 19 phenotypic measurements to explore the whole-organism physiological and genomic response to water deprivation in the desert-adapted cactus mouse (Peromyscus eremicus). The findings encompass the identification of differentially expressed genes and correlative analysis between phenotypes and gene expression patterns across multiple tissues. Specifically, we found robust activation of the vasopressin renin-angiotensin-aldosterone system (RAAS) pathways, whose primary function is to manage water and solute balance. Animals reduced food intake during water deprivation, and upregulation of PCK1 highlights the adaptive response to reduced oral intake via its actions aimed at maintained serum glucose levels. Even with such responses to maintain water balance, hemoconcentration still occurred, prompting a protective downregulation of genes responsible for the production of clotting factors while simultaneously enhancing angiogenesis which is thought to maintain tissue perfusion. In this study, we elucidate the complex mechanisms involved in water balance in the desert-adapted cactus mouse, P. eremicus. By prioritizing a comprehensive analysis of whole-organism physiology and multi-tissue gene expression in a simulated desert environment, we describe the complex response of regulatory processes.

Development and tissue specific expression of RAPGEF1 (C3G) transcripts having exons encoding disordered segments with predicted regulatory function.

BACKGROUND: The ubiquitously expressed Guanine nucleotide exchange factor, RAPGEF1 (C3G), is essential for early development of mouse embryos. It functions to regulate gene expression and cytoskeletal reorganization, thereby controlling cell proliferation and differentiation. While multiple transcripts have been predicted, their expression in mouse tissues has not been investigated in detail. METHODS & RESULTS: Full length RAPGEF1 isoforms primarily arise due to splicing at two hotspots, one involving exon-3, and the other involving exons 12-14 incorporating amino acids immediately following the Crk binding region of the protein. These isoforms vary in expression across embryonic and adult organs. We detected the presence of unannotated, and unpredicted transcripts with incorporation of cassette exons in various combinations, specifically in the heart, brain, testis and skeletal muscle. Isoform switching was detected as myocytes in culture and mouse embryonic stem cells were differentiated to form myotubes, and embryoid bodies respectively. The cassette exons encode a serine-rich polypeptide chain, which is intrinsically disordered, and undergoes phosphorylation. In silico structural analysis using AlphaFold indicated that the presence of cassette exons alters intra-molecular interactions, important for regulating catalytic activity. LZerD based docking studies predicted that the isoforms with one or more cassette exons differ in interaction with their target GTPase, RAP1A. CONCLUSIONS: Our results demonstrate the expression of novel RAPGEF1 isoforms, and predict cassette exon inclusion as an additional means of regulating RAPGEF1 activity in various tissues and during differentiation.

Current genomic deep learning models display decreased performance in cell type-specific accessible regions.

BACKGROUND: A number of deep learning models have been developed to predict epigenetic features such as chromatin accessibility from DNA sequence. Model evaluations commonly report performance genome-wide; however, cis regulatory elements (CREs), which play critical roles in gene regulation, make up only a small fraction of the genome. Furthermore, cell type-specific CREs contain a large proportion of complex disease heritability. RESULTS: We evaluate genomic deep learning models in chromatin accessibility regions with varying degrees of cell type specificity. We assess two modeling directions in the field: general purpose models trained across thousands of outputs (cell types and epigenetic marks) and models tailored to specific tissues and tasks. We find that the accuracy of genomic deep learning models, including two state-of-the-art general purpose models-Enformer and Sei-varies across the genome and is reduced in cell type-specific accessible regions. Using accessibility models trained on cell types from specific tissues, we find that increasing model capacity to learn cell type-specific regulatory syntax-through single-task learning or high capacity multi-task models-can improve performance in cell type-specific accessible regions. We also observe that improving reference sequence predictions does not consistently improve variant effect predictions, indicating that novel strategies are needed to improve performance on variants. CONCLUSIONS: Our results provide a new perspective on the performance of genomic deep learning models, showing that performance varies across the genome and is particularly reduced in cell type-specific accessible regions. We also identify strategies to maximize performance in cell type-specific accessible regions.

Enhancing portability of trans-ancestral polygenic risk scores through tissue-specific functional genomic data integration.

Portability of trans-ancestral polygenic risk scores is often confounded by differences in linkage disequilibrium and genetic architecture between ancestries. Recent literature has shown that prioritizing GWAS SNPs with functional genomic evidence over strong association signals can improve model portability. We leveraged three RegulomeDB-derived functional regulatory annotations-SURF, TURF, and TLand-to construct polygenic risk models across a set of quantitative and binary traits highlighting functional mutations tagged by trait-associated tissue annotations. Tissue-specific prioritization by TURF and TLand provide a significant improvement in model accuracy over standard polygenic risk score (PRS) models across all traits. We developed the Trans-ancestral Iterative Tissue Refinement (TITR) algorithm to construct PRS models that prioritize functional mutations across multiple trait-implicated tissues. TITR-constructed PRS models show increased predictive accuracy over single tissue prioritization. This indicates our TITR approach captures a more comprehensive view of regulatory systems across implicated tissues that contribute to variance in trait expression.

Extracellular vesicles carry transcriptional 'dark matter' revealing tissue-specific information.

From eukaryotes to prokaryotes, all cells secrete extracellular vesicles (EVs) as part of their regular homeostasis, intercellular communication, and cargo disposal. Accumulating evidence suggests that small EVs carry functional small RNAs, potentially serving as extracellular messengers and liquid-biopsy markers. Yet, the complete transcriptomic landscape of EV-associated small RNAs during disease progression is poorly delineated due to critical limitations including the protocols used for sequencing, suboptimal alignment of short reads (20-50 nt), and uncharacterized genome annotations-often denoted as the 'dark matter' of the genome. In this study, we investigate the EV-associated small unannotated RNAs that arise from endogenous genes and are part of the genomic 'dark matter', which may play a key emerging role in regulating gene expression and translational mechanisms. To address this, we created a distinct small RNAseq dataset from human prostate cancer & benign tissues, and EVs derived from blood (pre- & post-prostatectomy), urine, and human prostate carcinoma epithelial cell line. We then developed an unsupervised data-based bioinformatic pipeline that recognizes biologically relevant transcriptional signals irrespective of their genomic annotation. Using this approach, we discovered distinct EV-RNA expression patterns emerging from the un-annotated genomic regions (UGRs) of the transcriptomes associated with tissue-specific phenotypes. We have named these novel EV-associated small RNAs as 'EV-UGRs' or "EV-dark matter". Here, we demonstrate that EV-UGR gene expressions are downregulated by ∼100 fold (FDR < 0.05) in the circulating serum EVs from aggressive prostate cancer subjects. Remarkably, these EV-UGRs expression signatures were regained (upregulated) after radical prostatectomy in the same follow-up patients. Finally, we developed a stem-loop RT-qPCR assay that validated prostate cancer-specific EV-UGRs for selective fluid-based diagnostics. Overall, using an unsupervised data driven approach, we investigate the 'dark matter' of EV-transcriptome and demonstrate that EV-UGRs carry tissue-specific Information that significantly alters pre- and post-prostatectomy in the prostate cancer patients. Although further validation in randomized clinical trials is required, this new class of EV-RNAs hold promise in liquid-biopsy by avoiding highly invasive biopsy procedures in prostate cancer.

Expression genome-wide association study identifies key regulatory variants enriched with metabolic and immune functions in four porcine tissues.

BACKGROUND: Integration of high throughput DNA genotyping and RNA-sequencing data enables the discovery of genomic regions that regulate gene expression, known as expression quantitative trait loci (eQTL). In pigs, efforts to date have been mainly focused on purebred lines for traits with commercial relevance as such growth and meat quality. However, little is known on genetic variants and mechanisms associated with the robustness of an animal, thus its overall health status. Here, the liver, lung, spleen, and muscle transcriptomes of 100 three-way crossbred female finishers were studied, with the aim of identifying novel eQTL regulatory regions and transcription factors (TFs) associated with regulation of porcine metabolism and health-related traits. RESULTS: An expression genome-wide association study with 535,896 genotypes and the expression of 12,680 genes in liver, 13,310 genes in lung, 12,650 genes in spleen, and 12,595 genes in muscle resulted in 4,293, 10,630, 4,533, and 6,871 eQTL regions for each of these tissues, respectively. Although only a small fraction of the eQTLs were annotated as cis-eQTLs, these presented a higher number of polymorphisms per region and significantly stronger associations with their target gene compared to trans-eQTLs. Between 20 and 115 eQTL hotspots were identified across the four tissues. Interestingly, these were all enriched for immune-related biological processes. In spleen, two TFs were identified: ERF and ZNF45, with key roles in regulation of gene expression. CONCLUSIONS: This study provides a comprehensive analysis with more than 26,000 eQTL regions identified that are now publicly available. The genomic regions and their variants were mostly associated with tissue-specific regulatory roles. However, some shared regions provide new insights into the complex regulation of genes and their interactions that are involved with important traits related to metabolism and immunity.

Integrating muti-omics data to identify tissue-specific DNA methylation biomarkers for cancer risk.

The relationship between tissue-specific DNA methylation and cancer risk remains inadequately elucidated. Leveraging resources from the Genotype-Tissue Expression consortium, here we develop genetic models to predict DNA methylation at CpG sites across the genome for seven tissues and apply these models to genome-wide association study data of corresponding cancers, namely breast, colorectal, renal cell, lung, ovarian, prostate, and testicular germ cell cancers. At Bonferroni-corrected P < 0.05, we identify 4248 CpGs that are significantly associated with cancer risk, of which 95.4% (4052) are specific to a particular cancer type. Notably, 92 CpGs within 55 putative novel loci retain significant associations with cancer risk after conditioning on proximal signals identified by genome-wide association studies. Integrative multi-omics analyses reveal 854 CpG-gene-cancer trios, suggesting that DNA methylation at 309 distinct CpGs might influence cancer risk through regulating the expression of 205 unique cis-genes. These findings substantially advance our understanding of the interplay between genetics, epigenetics, and gene expression in cancer etiology.

Deciphering tissue-specific expression profiles of mitochondrial genome-encoded tRNAs and rRNAs through transcriptomic profiling in buffalo.

BACKGROUND: Mitochondria, essential for cellular energy production through oxidative phosphorylation (OXPHOS), integrate mt-DNA and nuclear-encoded genes. This cooperation extends to the mitochondrial translation machinery, involving crucial mtDNA-encoded RNAs: 22 tRNAs (mt-tRNAs) as adapters and two rRNAs (mt-rRNAs) for ribosomal assembly, enabling mitochondrial-encoded mRNA translation. Disruptions in mitochondrial gene expression can strongly impact energy generation and overall animal health. Our study investigates the tissue-specific expression patterns of mt-tRNAs and mt-rRNAs in buffalo. MATERIAL AND METHODS: To investigate the expression patterns of mt-tRNAs and mt-rRNAs in different tissues and gain a better understanding of tissue-specific variations, RNA-seq was performed on various tissues, such as the kidney, heart, brain, and ovary, from post-pubertal female buffaloes. Subsequently, we identified transcripts that were differentially expressed in various tissue comparisons. RESULTS: The findings reveal distinct expression patterns among specific mt-tRNA and mt-rRNA genes across various tissues, with some exhibiting significant upregulation and others demonstrating marked downregulation in specific tissue contexts. These identified variations reflect tissue-specific physiological roles, underscoring their significance in meeting the unique energy demands of each tissue. Notably, the brain exhibits the highest mtDNA copy numbers and an abundance of mitochondrial mRNAs of our earlier findings, potentially linked to the significant upregulation of mt-tRNAs in brain. This suggests a plausible association between mtDNA replication and the regulation of mtDNA gene expression. CONCLUSION: Overall, our study unveils the tissue-specific expression of mitochondrial-encoded non-coding RNAs in buffalo. As we proceed, our further investigations into tissue-specific mitochondrial proteomics and microRNA studies aim to elucidate the intricate mechanisms within mitochondria, contributing to tissue-specific mitochondrial attributes. This research holds promise to elucidate the critical role of mitochondria in animal health and disease.

Self-supervised learning for characterising histomorphological diversity and spatial RNA expression prediction across 23 human tissue types.

As vast histological archives are digitised, there is a pressing need to be able to associate specific tissue substructures and incident pathology to disease outcomes without arduous annotation. Here, we learn self-supervised representations using a Vision Transformer, trained on 1.7 M histology images across 23 healthy tissues in 838 donors from the Genotype Tissue Expression consortium (GTEx). Using these representations, we can automatically segment tissues into their constituent tissue substructures and pathology proportions across thousands of whole slide images, outperforming other self-supervised methods (43% increase in silhouette score). Additionally, we can detect and quantify histological pathologies present, such as arterial calcification (AUROC = 0.93) and identify missing calcification diagnoses. Finally, to link gene expression to tissue morphology, we introduce RNAPath, a set of models trained on 23 tissue types that can predict and spatially localise individual RNA expression levels directly from H&E histology (mean genes significantly regressed = 5156, FDR 1%). We validate RNAPath spatial predictions with matched ground truth immunohistochemistry for several well characterised control genes, recapitulating their known spatial specificity. Together, these results demonstrate how self-supervised machine learning when applied to vast histological archives allows researchers to answer questions about tissue pathology, its spatial organisation and the interplay between morphological tissue variability and gene expression.

Co-expression in tissue-specific gene networks links genes in cancer-susceptibility loci to known somatic driver genes.

BACKGROUND: The genetic background of cancer remains complex and challenging to integrate. Many somatic mutations within genes are known to cause and drive cancer, while genome-wide association studies (GWAS) of cancer have revealed many germline risk factors associated with cancer. However, the overlap between known somatic driver genes and positional candidate genes from GWAS loci is surprisingly small. We hypothesised that genes from multiple independent cancer GWAS loci should show tissue-specific co-regulation patterns that converge on cancer-specific driver genes. RESULTS: We studied recent well-powered GWAS of breast, prostate, colorectal and skin cancer by estimating co-expression between genes and subsequently prioritising genes that show significant co-expression with genes mapping within susceptibility loci from cancer GWAS. We observed that the prioritised genes were strongly enriched for cancer drivers defined by COSMIC, IntOGen and Dietlein et al. The enrichment of known cancer driver genes was most significant when using co-expression networks derived from non-cancer samples of the relevant tissue of origin. CONCLUSION: We show how genes within risk loci identified by cancer GWAS can be linked to known cancer driver genes through tissue-specific co-expression networks. This provides an important explanation for why seemingly unrelated sets of genes that harbour either germline risk factors or somatic mutations can eventually cause the same type of disease.

Multi-omic characterization of allele-specific regulatory variation in hybrid pigs.

Hybrid mapping is a powerful approach to efficiently identify and characterize genes regulated through mechanisms in cis. In this study, using reciprocal crosses of the phenotypically divergent Duroc and Lulai pig breeds, we perform a comprehensive multi-omic characterization of regulatory variation across the brain, liver, muscle, and placenta through four developmental stages. We produce one of the largest multi-omic datasets in pigs to date, including 16 whole genome sequenced individuals, as well as 48 whole genome bisulfite sequencing, 168 ATAC-Seq and 168 RNA-Seq samples. We develop a read count-based method to reliably assess allele-specific methylation, chromatin accessibility, and RNA expression. We show that tissue specificity was much stronger than developmental stage specificity in all of DNA methylation, chromatin accessibility, and gene expression. We identify 573 genes showing allele specific expression, including those influenced by parent-of-origin as well as allele genotype effects. We integrate methylation, chromatin accessibility, and gene expression data to show that allele specific expression can be explained in great part by allele specific methylation and/or chromatin accessibility. This study provides a comprehensive characterization of regulatory variation across multiple tissues and developmental stages in pigs.

Tissue of origin detection for cancer tumor using low-depth cfDNA samples through combination of tumor-specific methylation atlas and genome-wide methylation density in graph convolutional neural networks.

BACKGROUND: Cell free DNA (cfDNA)-based assays hold great potential in detecting early cancer signals yet determining the tissue-of-origin (TOO) for cancer signals remains a challenging task. Here, we investigated the contribution of a methylation atlas to TOO detection in low depth cfDNA samples. METHODS: We constructed a tumor-specific methylation atlas (TSMA) using whole-genome bisulfite sequencing (WGBS) data from five types of tumor tissues (breast, colorectal, gastric, liver and lung cancer) and paired white blood cells (WBC). TSMA was used with a non-negative least square matrix factorization (NNLS) deconvolution algorithm to identify the abundance of tumor tissue types in a WGBS sample. We showed that TSMA worked well with tumor tissue but struggled with cfDNA samples due to the overwhelming amount of WBC-derived DNA. To construct a model for TOO, we adopted the multi-modal strategy and used as inputs the combination of deconvolution scores from TSMA with other features of cfDNA. RESULTS: Our final model comprised of a graph convolutional neural network using deconvolution scores and genome-wide methylation density features, which achieved an accuracy of 69% in a held-out validation dataset of 239 low-depth cfDNA samples. CONCLUSIONS: In conclusion, we have demonstrated that our TSMA in combination with other cfDNA features can improve TOO detection in low-depth cfDNA samples.

Partitioning and aggregating cross-tissue and tissue-specific genetic effects to identify gene-trait associations.

TWAS have shown great promise in extending GWAS loci to a functional understanding of disease mechanisms. In an effort to fully unleash the TWAS and GWAS information, we propose MTWAS, a statistical framework that partitions and aggregates cross-tissue and tissue-specific genetic effects in identifying gene-trait associations. We introduce a non-parametric imputation strategy to augment the inaccessible tissues, accommodating complex interactions and non-linear expression data structures across various tissues. We further classify eQTLs into cross-tissue eQTLs and tissue-specific eQTLs via a stepwise procedure based on the extended Bayesian information criterion, which is consistent under high-dimensional settings. We show that MTWAS significantly improves the prediction accuracy across all 47 tissues of the GTEx dataset, compared with other single-tissue and multi-tissue methods, such as PrediXcan, TIGAR, and UTMOST. Applying MTWAS to the DICE and OneK1K datasets with bulk and single-cell RNA sequencing data on immune cell types showcases consistent improvements in prediction accuracy. MTWAS also identifies more predictable genes, and the improvement can be replicated with independent studies. We apply MTWAS to 84 UK Biobank GWAS studies, which provides insights into disease etiology.

High-throughput single-molecule long-read RNA sequencing analysis of tissue-specific genes and isoforms in lettuce (Lactuca sativa L.).

Lettuce is one of the most widely cultivated and consumed dicotyledonous vegetables globally. Despite the availability of its reference genome sequence, lettuce gene annotation remains incomplete, impeding comprehensive research and the broad application of genomic resources. Long-read RNA isoform sequencing (Iso-Seq) offers substantial advantages for analyzing RNA alternative splicing and aiding gene annotation, yet it faces throughput limitations. We present the HIT-ISOseq method tailored for bulk sample analysis, significantly enhancing RNA sequencing throughput on the PacBio platform by concatenating cDNA. Here we show, HIT-ISOseq generates 3-4 cDNA molecules per CCS read in lettuce, yielding 15.7 million long reads per PacBio Sequel II SMRT Cell 8 M. We validate its effectiveness in analyzing six lettuce tissue samples, including roots, stems, and leaves, revealing tissue-specific gene expression patterns and RNA isoforms. Leveraging diverse tissue long-read RNA sequencing, we refine the transcript annotation of the lettuce reference genome, expanding its GO and KEGG annotation repertoire. Collectively, this study serves as a foundational reference for genome annotation and the analysis of multi-sample isoform expression, utilizing high-throughput long-read transcriptome sequencing.

An integrative multi-context Mendelian randomization method for identifying risk genes across human tissues.

Mendelian randomization (MR) provides valuable assessments of the causal effect of exposure on outcome, yet the application of conventional MR methods for mapping risk genes encounters new challenges. One of the issues is the limited availability of expression quantitative trait loci (eQTLs) as instrumental variables (IVs), hampering the estimation of sparse causal effects. Additionally, the often context- or tissue-specific eQTL effects challenge the MR assumption of consistent IV effects across eQTL and GWAS data. To address these challenges, we propose a multi-context multivariable integrative MR framework, mintMR, for mapping expression and molecular traits as joint exposures. It models the effects of molecular exposures across multiple tissues in each gene region, while simultaneously estimating across multiple gene regions. It uses eQTLs with consistent effects across more than one tissue type as IVs, improving IV consistency. A major innovation of mintMR involves employing multi-view learning methods to collectively model latent indicators of disease relevance across multiple tissues, molecular traits, and gene regions. The multi-view learning captures the major patterns of disease relevance and uses these patterns to update the estimated tissue relevance probabilities. The proposed mintMR iterates between performing a multi-tissue MR for each gene region and joint learning the disease-relevant tissue probabilities across gene regions, improving the estimation of sparse effects across genes. We apply mintMR to evaluate the causal effects of gene expression and DNA methylation for 35 complex traits using multi-tissue QTLs as IVs. The proposed mintMR controls genome-wide inflation and offers insights into disease mechanisms.

Single-cell chromatin accessibility and transposable element landscapes reveal shared features of tissue-residing immune cells.

Tissue adaptation is required for regulatory T (Treg) cell function within organs. Whether this program shares aspects with other tissue-localized immune populations is unclear. Here, we analyzed single-cell chromatin accessibility data, including the transposable element (TE) landscape of CD45+ immune cells from colon, skin, adipose tissue, and spleen. We identified features of organ-specific tissue adaptation across different immune cells. Focusing on tissue Treg cells, we found conservation of the Treg tissue adaptation program in other tissue-localized immune cells, such as amphiregulin-producing T helper (Th)17 cells. Accessible TEs can act as regulatory elements, but their contribution to tissue adaptation is not understood. TE landscape analysis revealed an enrichment of specific transcription factor binding motifs in TE regions within accessible chromatin peaks. TEs, specifically from the LTR family, were located in enhancer regions and associated with tissue adaptation. These findings broaden our understanding of immune tissue residency and provide an important step toward organ-specific immune interventions.

HNF4A and HNF1A exhibit tissue specific target gene regulation in pancreatic beta cells and hepatocytes.

HNF4A and HNF1A encode transcription factors that are important for the development and function of the pancreas and liver. Mutations in both genes have been directly linked to Maturity Onset Diabetes of the Young (MODY) and type 2 diabetes (T2D) risk. To better define the pleiotropic gene regulatory roles of HNF4A and HNF1A, we generated a comprehensive genome-wide map of their binding targets in pancreatic and hepatic cells using ChIP-Seq. HNF4A was found to bind and regulate known (ACY3, HAAO, HNF1A, MAP3K11) and previously unidentified (ABCD3, CDKN2AIP, USH1C, VIL1) loci in a tissue-dependent manner. Functional follow-up highlighted a potential role for HAAO and USH1C as regulators of beta cell function. Unlike the loss-of-function HNF4A/MODY1 variant I271fs, the T2D-associated HNF4A variant (rs1800961) was found to activate AKAP1, GAD2 and HOPX gene expression, potentially due to changes in DNA-binding affinity. We also found HNF1A to bind to and regulate GPR39 expression in beta cells. Overall, our studies provide a rich resource for uncovering downstream molecular targets of HNF4A and HNF1A that may contribute to beta cell or hepatic cell (dys)function, and set up a framework for gene discovery and functional validation.

Tissue-specific and endogenous protein labeling with split fluorescent proteins.

The ability to label proteins by fusion with genetically encoded fluorescent proteins is a powerful tool for understanding dynamic biological processes. However, current approaches for expressing fluorescent protein fusions possess drawbacks, especially at the whole organism level. Expression by transgenesis risks potential overexpression artifacts while fluorescent protein insertion at endogenous loci is technically difficult and, more importantly, does not allow for tissue-specific study of broadly expressed proteins. To overcome these limitations, we have adopted the split fluorescent protein system mNeonGreen21-10/11 (split-mNG2) to achieve tissue-specific and endogenous protein labeling in zebrafish. In our approach, mNG21-10 is expressed under a tissue-specific promoter using standard transgenesis while mNG211 is inserted into protein-coding genes of interest using CRISPR/Cas-directed gene editing. Each mNG2 fragment on its own is not fluorescent, but when co-expressed the fragments self-assemble into a fluorescent complex. Here, we report successful use of split-mNG2 to achieve differential labeling of the cytoskeleton genes tubb4b and krt8 in various tissues. We also demonstrate that by anchoring the mNG21-10 component to specific cellular compartments, the split-mNG2 system can be used to manipulate protein localization. Our approach should be broadly useful for a wide range of applications.

Long-read RNA sequencing identifies region- and sex-specific C57BL/6J mouse brain mRNA isoform expression and usage.

Alternative splicing (AS) contributes to the biological heterogeneity between species, sexes, tissues, and cell types. Many diseases are either caused by alterations in AS or by alterations to AS. Therefore, measuring AS accurately and efficiently is critical for assessing molecular phenotypes, including those associated with disease. Long-read sequencing enables more accurate quantification of differentially spliced isoform expression than short-read sequencing approaches, and third-generation platforms facilitate high-throughput experiments. To assess differences in AS across the cerebellum, cortex, hippocampus, and striatum by sex, we generated and analyzed Oxford Nanopore Technologies (ONT) long-read RNA sequencing (lrRNA-Seq) C57BL/6J mouse brain cDNA libraries. From > 85 million reads that passed quality control metrics, we calculated differential gene expression (DGE), differential transcript expression (DTE), and differential transcript usage (DTU) across brain regions and by sex. We found significant DGE, DTE, and DTU across brain regions and that the cerebellum had the most differences compared to the other three regions. Additionally, we found region-specific differential splicing between sexes, with the most sex differences in DTU in the cortex and no DTU in the hippocampus. We also report on two distinct patterns of sex DTU we observed, sex-divergent and sex-specific, that could potentially help explain sex differences in the prevalence and prognosis of various neurological and psychiatric disorders in future studies. Finally, we built a Shiny web application for researchers to explore the data further. Our study provides a resource for the community; it underscores the importance of AS in biological heterogeneity and the utility of long-read sequencing to better understand AS in the brain.