Differences in molecular sampling and data processing explain variation among single-cell and single-nucleus RNA-seq experiments.

A mechanistic understanding of the biological and technical factors that impact transcript measurements is essential to designing and analyzing single-cell and single-nucleus RNA sequencing experiments. Nuclei contain the same pre-mRNA population as cells, but they contain a small subset of the mRNAs. Nonetheless, early studies argued that single-nucleus analysis yielded results comparable to cellular samples if pre-mRNA measurements were included. However, typical workflows do not distinguish between pre-mRNA and mRNA when estimating gene expression, and variation in their relative abundances across cell types has received limited attention. These gaps are especially important given that incorporating pre-mRNA has become commonplace for both assays, despite known gene length bias in pre-mRNA capture. Here, we reanalyze public data sets from mouse and human to describe the mechanisms and contrasting effects of mRNA and pre-mRNA sampling on gene expression and marker gene selection in single-cell and single-nucleus RNA-seq. We show that pre-mRNA levels vary considerably among cell types, which mediates the degree of gene length bias and limits the generalizability of a recently published normalization method intended to correct for this bias. As an alternative, we repurpose an existing post hoc gene length-based correction method from conventional RNA-seq gene set enrichment analysis. Finally, we show that inclusion of pre-mRNA in bioinformatic processing can impart a larger effect than assay choice itself, which is pivotal to the effective reuse of existing data. These analyses advance our understanding of the sources of variation in single-cell and single-nucleus RNA-seq experiments and provide useful guidance for future studies.

An Interleukin-25-Mediated Autoregulatory Circuit in Keratinocytes Plays a Pivotal Role in Psoriatic Skin Inflammation.

Psoriasis is a chronic autoinflammatory skin disease. Although interleukin-17, derived from lymphocytes, has been shown to be critical in psoriasis, the initiation and maintenance of chronic skin inflammation has not been well understood. IL-25 (also called IL-17E), another IL-17 family cytokine, is well known to regulate allergic responses and type 2 immunity. Here we have shown that IL-25, also highly expressed in the lesional skin of psoriasis patients, was regulated by IL-17 in murine skin of a imiquimod (IMQ)-induced psoriasis model. IL-25 injection induced skin inflammation, whereas germline or keratinocyte-specific deletion of IL-25 caused resistance to IMQ-induced psoriasis. Via IL-17RB expression in keratinocytes, IL-25 stimulated the proliferation of keratinocytes and induced the production of inflammatory cytokines and chemokines, via activation of the STAT3 transcription factor. Thus, our data demonstrate that an IL-17-induced autoregulatory circuit in keratinocytes is mediated by IL-25 and suggest that this circuit could be targeted in the treatment of psoriasis patients.

Co-inhibitory Molecule B7 Superfamily Member 1 Expressed by Tumor-Infiltrating Myeloid Cells Induces Dysfunction of Anti-tumor CD8+ T Cells.

The molecular mechanisms whereby CD8+ T cells become "exhausted" in the tumor microenvironment remain unclear. Programmed death ligand-1 (PD-L1) is upregulated on tumor cells and PD-1-PD-L1 blockade has significant efficacy in human tumors; however, most patients do not respond, suggesting additional mechanisms underlying T cell exhaustion. B7 superfamily member 1 (B7S1), also called B7-H4, B7x, or VTCN1, negatively regulates T cell activation. Here we show increased B7S1 expression on myeloid cells from human hepatocellular carcinoma correlated with CD8+ T cell dysfunction. B7S1 inhibition suppressed development of murine tumors. Putative B7S1 receptor was co-expressed with PD-1 but not T cell immunoglobulin and mucin-domain containing-3 (Tim-3) at an activated state of early tumor-infiltrating CD8+ T cells, and B7S1 promoted T cell exhaustion, possibly through Eomes overexpression. Combinatorial blockade of B7S1 and PD-1 synergistically enhanced anti-tumor immune responses. Collectively, B7S1 initiates dysfunction of tumor-infiltrating CD8+ T cells and may be targeted for cancer immunotherapy.

EMT activates exocytotic Rabs to coordinate invasion and immunosuppression in lung cancer.

Epithelial-to-mesenchymal transition (EMT) underlies immunosuppression, drug resistance, and metastasis in epithelial malignancies. However, the way in which EMT orchestrates disparate biological processes remains unclear. Here, we identify an EMT-activated vesicular trafficking network that coordinates promigratory focal adhesion dynamics with an immunosuppressive secretory program in lung adenocarcinoma (LUAD). The EMT-activating transcription factor ZEB1 drives exocytotic vesicular trafficking by relieving Rab6A, Rab8A, and guanine nucleotide exchange factors from miR-148a-dependent silencing, thereby facilitating MMP14-dependent focal adhesion turnover in LUAD cells and autotaxin-mediated CD8+ T cell exhaustion, indicating that cell-intrinsic and extrinsic processes are linked through a microRNA that coordinates vesicular trafficking networks. Blockade of ZEB1-dependent secretion reactivates antitumor immunity and negates resistance to PD-L1 immune checkpoint blockade, an important clinical problem in LUAD. Thus, EMT activates exocytotic Rabs to drive a secretory program that promotes invasion and immunosuppression in LUAD.

The signaling suppressor CIS controls proallergic T cell development and allergic airway inflammation.

Transcription factors of the STAT family are critical in the cytokine-mediated functional differentiation of CD4(+) helper T cells. Signaling inhibitors of the SOCS family negatively regulate the activation of STAT proteins; however, their roles in the differentiation and function of helper T cells are not well understood. Here we found that the SOCS protein CIS, which was substantially induced by interleukin 4 (IL-4), negatively regulated the activation of STAT3, STAT5 and STAT6 in T cells. CIS-deficient mice spontaneously developed airway inflammation, and CIS deficiency in T cells led to greater susceptibility to experimental allergic asthma. CIS-deficient T cells showed enhanced differentiation into the TH2 and TH9 subsets of helper T cells. STAT5 and STAT6 regulated IL-9 expression by directly binding to the Il9 promoter. Our data thus demonstrate a critical role for CIS in controlling the proallergic generation of helper T cells.

N121T and N121S substitutions on the SARS-CoV-2 spike protein impact on serum neutralization.

The N121 site on the spike protein of SARS-CoV-2 is associated with heme and its metabolite, biliverdin, which can affect antibody binding. Both N121T and N121S substitutions have been observed in natural conditions and in a hamster model of dual infection with SARS-CoV-2 and Influenza A virus. Serum pseudotype neutralization assays against HIV-1 particles carrying wild-type, N121T, and N121S spikes with immune mouse and human sera revealed that N121T and N121S mutations had a greater impact on serum neutralization than biliverdin treatment. Although N121T and N121S substitutions are not currently major SARS-CoV-2 variants of concern, this study could provide fundamental information to prepare for potential future mutations at the N121 site of SARS-CoV-2.

A mouse xenograft long-term replication yields a SARS-CoV-2 Delta mutant with increased lethality.

We recently established a long-term SARS-CoV-2 infection model using lung-cancer xenograft mice and identified mutations that arose in the SARS-CoV-2 genome during long-term propagation. Here, we applied our model to the SARS-CoV-2 Delta variant, which has increased transmissibility and immune escape compared with ancestral SARS-CoV-2. We observed limited mutations in SARS-CoV-2 Delta during long-term propagation, including two predominant mutations: R682W in the spike protein and L330W in the nucleocapsid protein. We analyzed two representative isolates, Delta-10 and Delta-12, with both predominant mutations and some additional mutations. Delta-10 and Delta-12 showed lower replication capacity compared with SARS-CoV-2 Delta in cultured cells; however, Delta-12 was more lethal in K18-hACE2 mice compared with SARS-CoV-2 Delta and Delta-10. Mice infected with Delta-12 had higher viral titers, more severe histopathology in the lungs, higher chemokine expression, increased astrocyte and microglia activation, and extensive neutrophil infiltration in the brain. Brain tissue hemorrhage and mild vacuolation were also observed, suggesting that the high lethality of Delta-12 was associated with lung and brain pathology. Our long-term infection model can provide mutant viruses derived from SARS-CoV-2 Delta and knowledge about the possible contributions of emergent mutations to the properties of new variants.

Interleukin-17B Antagonizes Interleukin-25-Mediated Mucosal Inflammation.

The interleukin-17 (IL-17) family of cytokines has emerged as a critical player in inflammatory diseases. Among them, IL-25 has been shown to be important in allergic inflammation and protection against parasitic infection. Here we have demonstrated that IL-17B, a poorly understood cytokine, functions to inhibit IL-25-driven inflammation. IL-17B and IL-25, both binding to the interleukin-17 receptor B (IL-17RB), were upregulated in their expression after acute colonic inflammation. Individual inhibition of these cytokines revealed opposing functions in colon inflammation: IL-25 was pathogenic but IL-17B was protective. Similarly opposing phenotypes were observed in Citrobacter rodentium infection and allergic asthma. Moreover, IL-25 was found to promote IL-6 production from colon epithelial cells, which was inhibited by IL-17B. Therefore, our data demonstrate that IL-17B is an anti-inflammatory cytokine in the IL-17 family.

Analysis of SARS-CoV-2 omicron mutations that emerged during long-term replication in a lung cancer xenograft mouse model.

SARS-CoV-2 Omicron has the largest number of mutations among all the known SARS-CoV-2 variants. The presence of these mutations might explain why Omicron is more infectious and vaccines have lower efficacy to Omicron than other variants, despite lower virulence of Omicron. We recently established a long-term in vivo replication model by infecting Calu-3 xenograft tumors in immunodeficient mice with parental SARS-CoV-2 and found that various mutations occurred majorly in the spike protein during extended replication. To investigate whether there are differences in the spectrum and frequency of mutations between parental SARS-CoV-2 and Omicron, we here applied this model to Omicron. At 30 days after infection, we found that the virus was present at high titers in the tumor tissues and had developed several rare sporadic mutations, mainly in ORF1ab with additional minor spike protein mutations. Many of the mutant isolates had higher replicative activity in Calu-3 cells compared with the original SARS-CoV-2 Omicron virus, suggesting that the novel mutations contributed to increased viral replication. Serial propagation of SARS-CoV-2 Omicron in cultured Calu-3 cells resulted in several rare sporadic mutations in various viral proteins with no mutations in the spike protein. Therefore, the genome of SARS-CoV-2 Omicron seems largely stable compared with that of the parental SARS-CoV-2 during extended replication in Calu-3 cells and xenograft model. The sporadic mutations and modified growth properties observed in Omicron might explain the emergence of Omicron sublineages. However, we cannot exclude the possibility of some differences in natural infection.

Medullins A-H, Sesquiterpenes from Perenniporia medulla-panis, and Their Cellular Signaling Mechanism in HDF Cells.

A chemical investigation of an ethyl acetate-soluble layer in the culture broth of Perenniporia medulla-panis resulted in the isolation of eight novel sesquiterpenes conjugated Gly (1), l-Val (2), l-Ala (3), l-Tyr (4), l-Thr (5), l-Ile (6), l-Leu (7), and l-Phe (8). Elucidation of their structures was performed through comprehensive spectroscopic analysis. The absolute configuration of the sesquiterpene skeleton was ascertained using modified Mosher's methods. The configurations of the amino acid units in compounds 2-8 were identified through acid hydrolysis followed by LC-MS analysis employing Marfey's method. Compounds 1-3 and 5-8 showed significant regulating effect on MAP kinase activity (p-ERK and p-JNK) in human diploid fibroblast (HDF) cells.

In vitro and in vivo suppression of SARS-CoV-2 replication by a modified, short, cell-penetrating peptide targeting the C-terminal domain of the viral spike protein.

Peptides are promising therapeutic agents for COVID-19 because of their specificity, easy synthesis, and ability to be fine-tuned. We previously demonstrated that a cell-permeable peptide corresponding to the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Spike C-terminal domain (CD) inhibits the interaction between viral spike and nucleocapsid proteins that results in SARS-CoV-2 replication in vitro. Here, we used docking studies to design R-t-Spike CD(D), a more potent short cell-penetrating peptide composed of all D-form amino acids and evaluated its inhibitory effect against the replication of SARS-CoV-2 S clade and other variants. R-t-Spike CD(D) was internalized into Vero cells and Calu-3 cells and suppressed the replication of SARS-CoV-2 S clade, delta variant, and omicron variant with higher potency than the original peptide. In hemizygous K18-hACE2 mice, intratracheal administration of R-t-Spike CD(D) effectively delivered the peptide to the trachea and lungs, whereas intranasal administration delivered the peptide mostly to the upper respiratory system and stomach, and a small amount to the lungs. Administration by either route reduced viral loads in mouse lungs and turbinates. Furthermore, intranasally administered R-t-Spike CD(D) mitigated pathological change in the lungs and increased the survival of mice after infection with the SARS-CoV-2 S clade or delta variant. Our data suggest that R-t-Spike CD(D) has potential as a therapeutic agent against SARS-CoV-2 infection.

HDAC inhibitor and proteasome inhibitor induce cleavage and exosome-mediated secretion of HSP90 in mouse pluripotent stem cells.

Heat shock protein 90 (HSP90), one of the molecular chaperones, stabilizes several proteins necessary to maintain pluripotency of embryonic stem (ES) cells. Recently, we reported that HDAC inhibitors and proteasome inhibitors down-regulate HSP90 activity through HSP90 cleavage induced by reactive oxygen species (ROS) generation and caspase 10 activation in various cancer cells. In this study, we investigated HSP90 cleavage in mouse ES cells. HDAC inhibitors and proteasome inhibitors induced HSP90 cleavage in the mouse ES cell line R1, and the cleaved HSP90 was barely found in the cells and instead secreted out of the cells through the exosome. The HSP90 cleavage was associated with ROS generation and caspase 10 activation. In addition, HDAC inhibitor and proteasome inhibitor induced Fas expression, and the inhibition of caspase 8, a downstream molecule of Fas, blocked HSP90 cleavage. Therefore, HDAC inhibitor- and proteasome inhibitor-mediated HSP90 cleavage was induced by ROS generation and Fas expression. We observed similar results in mouse induced pluripotent stem (iPS) cells. Taken together, HSP90 cleavage was induced in mouse pluripotent cells similarly to cancer cells but differently regulated through Fas expression and exosomal secretion. These findings will be helpful in elucidating the regulation of HSP90 upon stress in pluripotent stem cells.

Epigenetic interplay between methylation and miRNA in bladder cancer: focus on isoform expression.

BACKGROUND: Various epigenetic factors are responsible for the non-genetic regulation on gene expression. The epigenetically dysregulated oncogenes or tumor suppressors by miRNA and/or DNA methylation are often observed in cancer cells. Each of these epigenetic regulators has been studied well in cancer progressions; however, their mutual regulatory relationship in cancer still remains unclear. In this study, we propose an integrative framework to systematically investigate epigenetic interactions between miRNA and methylation at the alternatively spliced mRNA level in bladder cancer. Each of these epigenetic regulators has been studied well in cancer progressions; however, their mutual regulatory relationship in cancer still remains unclear. RESULTS: The integrative analyses yielded 136 significant combinations (methylation, miRNA and isoform). Further, overall survival analysis on the 136 combinations based on methylation and miRNA, high and low expression groups resulted in 13 combinations associated with survival. Additionally, different interaction patterns were examined. CONCLUSIONS: Our study provides a higher resolution of molecular insight into the crosstalk between two epigenetic factors, DNA methylation and miRNA. Given the importance of epigenetic interactions and alternative splicing in cancer, it is timely to identify and understand the underlying mechanisms based on epigenetic markers and their interactions in cancer, leading to alternative splicing with primary functional impact.

Para sport translation of the IOC consensus on recording and reporting of data for injury and illness in sport.

In 2020, the IOC proposed a universal methodology for the recording and reporting of data for injury and illness in sport. Para sport is played by individuals with impairment, and they have a unique set of considerations not captured by these recommendations. Therefore, the aim of this addendum to IOC consensus statement was to guide the Para sport researcher through the complexities and nuances that should be taken into consideration when collecting, registering, reporting and interpreting data regarding Para athlete health. To develop this translation, experts in the field of Para sports medicine and epidemiology conducted a formal consensus development process, which began in March 2020 with the formation of a consensus group that worked over eight phases, incorporating three virtual consensus meetings to finalise the translation. This translation is consistent with the IOC consensus statement, yet provides more detailed Para athlete specific definitions and recommendations on study population, specifically, diagnostic and eligible impairment categorisation and recording of adaptive equipment, and defining and classifying health problems in the context of Para sport. Additionally, recommendations and Para athlete specific examples are described with regards to injury mechanism, mode of onset, injury and illness classification, duration, capturing and reporting exposure and risk. Finally, methods and considerations are provided to cater to the varied needs of athletes with impairment with respect to data collection tools. This harmonisation will allow the science to develop and facilitate a more accurate understanding of injury and illness patterns for tailoring evidence-informed prevention programmes and enabling better planning of medical services for Para sport events.

Alternative Splicing Regulation of Low-Frequency Genetic Variants in Exon 2 of TREM2 in Alzheimer's Disease by Splicing-Based Aggregation.

TREM2 is among the most well-known Alzheimer's disease (AD) risk genes; however, the functional roles of its AD-associated variants remain to be elucidated, and most known risk alleles are low-frequency variants whose investigation is challenging. Here, we utilized a splicing-guided aggregation method in which multiple low-frequency TREM2 variants were bundled together to investigate the functional impact of those variants on alternative splicing in AD. We analyzed whole genome sequencing (WGS) and RNA-seq data generated from cognitively normal elderly controls (CN) and AD patients in two independent cohorts, representing three regions in the frontal lobe of the human brain: the dorsolateral prefrontal cortex (CN = 213 and AD = 376), frontal pole (CN = 72 and AD = 175), and inferior frontal (CN = 63 and AD = 157). We observed an exon skipping event in the second exon of TREM2, with that exon tending to be more frequently skipped (p = 0.0012) in individuals having at least one low-frequency variant that caused loss-of-function for a splicing regulatory element. In addition, genes differentially expressed between AD patients with high vs. low skipping of the second exon (i.e., loss of a TREM2 functional domain) were significantly enriched in immune-related pathways. Our splicing-guided aggregation method thus provides new insight into the regulation of alternative splicing of the second exon of TREM2 by low-frequency variants and could be a useful tool for further exploring the potential molecular mechanisms of multiple, disease-associated, low-frequency variants.

Topoisomerase IIIß Deficiency Induces Neuro-Behavioral Changes and Brain Connectivity Alterations in Mice.

Topoisomerase IIIß (Top3ß), the only dual-activity topoisomerase in mammals that can change topology of both DNA and RNA, is known to be associated with neurodevelopment and mental dysfunction in humans. However, there is no report showing clear associations of Top3ß with neuropsychiatric phenotypes in mice. Here, we investigated the effect of Top3ß on neuro-behavior using newly generated Top3ß deficient (Top3ß-/-) mice. We found that Top3ß-/- mice showed decreased anxiety and depression-like behaviors. The lack of Top3ß was also associated with changes in circadian rhythm. In addition, a clear expression of Top3ß was demonstrated in the central nervous system of mice. Positron emission tomography/computed tomography (PET/CT) analysis revealed significantly altered connectivity between many brain regions in Top3ß-/- mice, including the connectivity between the olfactory bulb and the cerebellum, the connectivity between the amygdala and the olfactory bulb, and the connectivity between the globus pallidus and the optic nerve. These connectivity alterations in brain regions are known to be linked to neurodevelopmental as well as psychiatric and behavioral disorders in humans. Therefore, we conclude that Top3ß is essential for normal brain function and behavior in mice and that Top3ß could be an interesting target to study neuropsychiatric disorders in humans.

Personal Health Information Inference Using Machine Learning on RNA Expression Data from Patients With Cancer: Algorithm Validation Study.

BACKGROUND: As the need for sharing genomic data grows, privacy issues and concerns, such as the ethics surrounding data sharing and disclosure of personal information, are raised. OBJECTIVE: The main purpose of this study was to verify whether genomic data is sufficient to predict a patient's personal information. METHODS: RNA expression data and matched patient personal information were collected from 9538 patients in The Cancer Genome Atlas program. Five personal information variables (age, gender, race, cancer type, and cancer stage) were recorded for each patient. Four different machine learning algorithms (support vector machine, decision tree, random forest, and artificial neural network) were used to determine whether a patient's personal information could be accurately predicted from RNA expression data. Performance measurement of the prediction models was based on the accuracy and area under the receiver operating characteristic curve. We selected five cancer types (breast carcinoma, kidney renal clear cell carcinoma, head and neck squamous cell carcinoma, low-grade glioma, and lung adenocarcinoma) with large samples sizes to verify whether predictive accuracy would differ between them. We also validated the efficacy of our four machine learning models in analyzing normal samples from 593 cancer patients. RESULTS: In most samples, personal information with high genetic relevance, such as gender and cancer type, could be predicted from RNA expression data alone. The prediction accuracies for gender and cancer type, which were the best models, were 0.93-0.99 and 0.78-0.94, respectively. Other aspects of personal information, such as age, race, and cancer stage, were difficult to predict from RNA expression data, with accuracies ranging from 0.0026-0.29, 0.76-0.96, and 0.45-0.79, respectively. Among the tested machine learning methods, the highest predictive accuracy was obtained using the support vector machine algorithm (mean accuracy 0.77), while the lowest accuracy was obtained using the random forest method (mean accuracy 0.65). Gender and race were predicted more accurately than other variables in the samples. On average, the accuracy of cancer stage prediction ranged between 0.71-0.67, while the age prediction accuracy ranged between 0.18-0.23 for the five cancer types. CONCLUSIONS: We attempted to predict patient information using RNA expression data. We found that some identifiers could be predicted, but most others could not. This study showed that personal information available from RNA expression data is limited and this information cannot be used to identify specific patients.

High incidence of injuries at the Pyeongchang 2018 Paralympic Winter Games: a prospective cohort study of 6804 athlete days.

OBJECTIVE: To describe the epidemiology of sports injury at the Pyeongchang 2018 Paralympic Winter Games. METHODS: 567 athletes from 49 countries were monitored daily for 12 days over the Pyeongchang 2018 Paralympic Winter Games (6804 athlete days). Injury data were obtained daily from teams with their own medical support (41 teams and 557 athletes) and teams without their own medical support (8 teams and 10 athletes) through two electronic data capturing systems. RESULTS: 112 of 567 athletes (19.8%) reported a total of 142 injuries, with an injury incidence rate (IR) of 20.9 per 1000 athlete days (95% CI 17.4 to 25.0). The highest IR was reported for para snowboard (IR of 40.5 per 1000 athlete days [95% CI 28.5 to 57.5]; p<0.02), particularly in the lower limb and head/face/neck anatomical areas. Across all sports at the Games, acute traumatic injuries (IR of 16.2 per 1000 athlete days [95% CI 13.2 to 19.8]) and injuries to the shoulder/arm/elbow complex (IR of 5.7 per 1000 athlete days [95% CI 4.2 to 7.8]) were most common. However, most injuries (78.9%) did not require time loss. CONCLUSION: The new Paralympic Winter Games sport of Para snowboard requires attention to implement actions that will reduce injury risk. The shoulder was the most injured single joint-a consistent finding in elite para sport.

In Situ Thermal Atomic Layer Etching for Sub-5 nm InGaAs Multigate MOSFETs.

Thermal atomic layer etching (ALE) was demonstrated on ternary III-V compound semiconductors. In particular, thermal ALE on InGaAs and InAlAs was achieved with sequential, self-limiting fluorination and ligand-exchange reactions using hydrogen fluoride (HF) as the fluorination reactant and dimethylaluminum chloride (DMAC) as the ligand-exchange reactant. Thermal ALE was investigated on planar surfaces and three-dimensional nanostructures. The measured radial etch rates on In0.53Ga0.47As and In0.52Al0.48As vertical nanowires (VNWs) at 300 °C were 0.24 and 0.62 Å/cycle, respectively. An optimized thermal ALE process did not increase the surface roughness after 200 cycles. The etching process also displayed selectivity and orientation dependence. This new thermal ALE process in combination with in situ atomic layer deposition (ALD) was used to fabricate InGaAs gate-all-around structures with minimum width down to 3 nm. The in situ ALE-ALD process produced a sharp vertical MOS interface. Finally, the merits of thermal ALE were demonstrated in the fabrication of n-channel InGaAs FinFETs with record ON-state and OFF-state transistor performance. On the basis of this transistor demonstration, thermal ALE shows great promise for high-volume device manufacturing.

Alternative Splicing Regulation of an Alzheimer's Risk Variant in CLU.

Clusterin (CLU) is one of the risk genes most associated with late onset Alzheimer's disease (AD), and several genetic variants in CLU are associated with AD risk. However, the functional role of known AD risk genetic variants in CLU has been little explored. We investigated the effect of an AD risk variant (rs7982) in the 5th exon of CLU on alternative splicing by using an integrative approach of brain-tissue-based RNA-Seq and whole genome sequencing data from Accelerating Medicines Partnership-Alzheimer's Disease (AMP-AD). RNA-Seq data were generated from three regions in the temporal lobe of the brain-the temporal cortex, superior temporal gyrus, and parahippocampal gyrus. The rs7982 was significantly associated with intron retention (IR) of the 5th exon of CLU; as the number of alternative alleles (G) increased, the IR rates decreased more significantly in females than in males. Our results suggest a sex-dependent role of rs7982 in AD pathogenesis via splicing regulation.