NSD2 dimethylation at H3K36 promotes lung adenocarcinoma pathogenesis.

The etiological role of NSD2 enzymatic activity in solid tumors is unclear. Here we show that NSD2, via H3K36me2 catalysis, cooperates with oncogenic KRAS signaling to drive lung adenocarcinoma (LUAD) pathogenesis. In vivo expression of NSD2E1099K, a hyperactive variant detected in individuals with LUAD, rapidly accelerates malignant tumor progression while decreasing survival in KRAS-driven LUAD mouse models. Pathologic H3K36me2 generation by NSD2 amplifies transcriptional output of KRAS and several complementary oncogenic gene expression programs. We establish a versatile in vivo CRISPRi-based system to test gene functions in LUAD and find that NSD2 loss strongly attenuates tumor progression. NSD2 knockdown also blocks neoplastic growth of PDXs (patient-dervived xenografts) from primary LUAD. Finally, a treatment regimen combining NSD2 depletion with MEK1/2 inhibition causes nearly complete regression of LUAD tumors. Our work identifies NSD2 as a bona fide LUAD therapeutic target and suggests a pivotal epigenetic role of the NSD2-H3K36me2 axis in sustaining oncogenic signaling.

Comparative analysis of in-silico tools in identifying pathogenic variants in dominant inherited retinal diseases.

Inherited retinal diseases (IRDs) are a group of rare genetic eye conditions that cause blindness. Despite progress in identifying genes associated with IRDs, improvements are necessary for classifying rare autosomal dominant (AD) disorders. AD diseases are highly heterogenous, with causal variants being restricted to specific amino acid changes within certain protein domains, making AD conditions difficult to classify. Here, we aim to determine the top-performing in-silico tools for predicting the pathogenicity of AD IRD variants. We annotated variants from ClinVar and benchmarked 39 variant classifier tools on IRD genes, split by inheritance pattern. Using area-under-the-curve (AUC) analysis, we determined the top-performing tools and defined thresholds for variant pathogenicity. Top-performing tools were assessed using genome sequencing on a cohort of participants with IRDs of unknown etiology. MutScore achieved the highest accuracy within AD genes, yielding an AUC of 0.969. When filtering for AD gain-of-function and dominant negative variants, BayesDel had the highest accuracy with an AUC of 0.997. Five participants with variants in NR2E3, RHO, GUCA1A, and GUCY2D were confirmed to have dominantly inherited disease based on pedigree, phenotype, and segregation analysis. We identified two uncharacterized variants in GUCA1A (c.428T>A, p.Ile143Thr) and RHO (c.631C>G, p.His211Asp) in three participants. Our findings support using a multi-classifier approach comprised of new missense classifier tools to identify pathogenic variants in participants with AD IRDs. Our results provide a foundation for improved genetic diagnosis for people with IRDs.

Hepatitis B virus e antigen induces atypical metabolism and differentially regulates programmed cell deaths of macrophages.

Macrophages can undergo M1-like proinflammatory polarization with low oxidative phosphorylation (OXPHOS) and high glycolytic activities or M2-like anti-inflammatory polarization with the opposite metabolic activities. Here we show that M1-like macrophages induced by hepatitis B virus (HBV) display high OXPHOS and low glycolytic activities. This atypical metabolism induced by HBV attenuates the antiviral response of M1-like macrophages and is mediated by HBV e antigen (HBeAg), which induces death receptor 5 (DR5) via toll-like receptor 4 (TLR4) to induce death-associated protein 3 (DAP3). DAP3 then induces the expression of mitochondrial genes to promote OXPHOS. HBeAg also enhances the expression of glutaminases and increases the level of glutamate, which is converted to α-ketoglutarate, an important metabolic intermediate of the tricarboxylic acid cycle, to promote OXPHOS. The induction of DR5 by HBeAg leads to apoptosis of M1-like and M2-like macrophages, although HBeAg also induces pyroptosis of the former. These findings reveal novel activities of HBeAg, which can reprogram mitochondrial metabolism and trigger different programmed cell death responses of macrophages depending on their phenotypes to promote HBV persistence.

SOX7 deficiency causes ventricular septal defects through its effects on endocardial-to-mesenchymal transition and the expression of Wnt4 and Bmp2.

SOX7 is a transcription factor-encoding gene located in a region on chromosome 8p23.1 that is recurrently deleted in individuals with ventricular septal defects (VSDs). We have previously shown that Sox7-/- embryos die of heart failure around E11.5. Here, we demonstrate that these embryos have hypocellular endocardial cushions with severely reduced numbers of mesenchymal cells. Ablation of Sox7 in the endocardium also resulted in hypocellular endocardial cushions, and we observed VSDs in rare E15.5 Sox7flox/-;Tie2-Cre and Sox7flox/flox;Tie2-Cre embryos that survived to E15.5. In atrioventricular explant studies, we showed that SOX7 deficiency leads to a severe reduction in endocardial-to-mesenchymal transition (EndMT). RNA-seq studies performed on E9.5 Sox7-/- heart tubes revealed severely reduced Wnt4 transcript levels. Wnt4 is expressed in the endocardium and promotes EndMT by acting in a paracrine manner to increase the expression of Bmp2 in the myocardium. Both WNT4 and BMP2 have been previously implicated in the development of VSDs in individuals with 46,XX sex reversal with dysgenesis of kidney, adrenals and lungs (SERKAL) syndrome and in individuals with short stature, facial dysmorphism and skeletal anomalies with or without cardiac anomalies 1 (SSFSC1) syndrome, respectively. We now show that Sox7 and Wnt4 interact genetically in the development of VSDs through their additive effects on endocardial cushion development with Sox7+/-;Wnt4+/- double heterozygous embryos having hypocellular endocardial cushions and perimembranous and muscular VSDs not seen in their Sox7+/- and Wnt4+/- littermates. These results provide additional evidence that SOX7, WNT4 and BMP2 function in the same pathway during mammalian septal development and that their deficiency can contribute to the development of VSDs in humans.

Systematic assessment of the contribution of structural variants to inherited retinal diseases.

Despite increasing success in determining genetic diagnosis for patients with inherited retinal diseases (IRDs), mutations in about 30% of the IRD cases remain unclear or unsettled after targeted gene panel or whole exome sequencing. In this study, we aimed to investigate the contributions of structural variants (SVs) to settling the molecular diagnosis of IRD with whole-genome sequencing (WGS). A cohort of 755 IRD patients whose pathogenic mutations remain undefined were subjected to WGS. Four SV calling algorithms including include MANTA, DELLY, LUMPY and CNVnator were used to detect SVs throughout the genome. All SVs identified by any one of these four algorithms were included for further analysis. AnnotSV was used to annotate these SVs. SVs that overlap with known IRD-associated genes were examined with sequencing coverage, junction reads and discordant read pairs. Polymerase Chain Reaction (PCR) followed by Sanger sequencing was used to further confirm the SVs and identify the breakpoints. Segregation of the candidate pathogenic alleles with the disease was performed when possible. A total of 16 candidate pathogenic SVs were identified in 16 families, including deletions and inversions, representing 2.1% of patients with previously unsolved IRDs. Autosomal dominant, autosomal recessive and X-linked inheritance of disease-causing SVs were observed in 12 different genes. Among these, SVs in CLN3, EYS and PRPF31 were found in multiple families. Our study suggests that the contribution of SVs detected by short-read WGS is about 0.25% of our IRD patient cohort and is significantly lower than that of single nucleotide changes and small insertions and deletions.

Ovol1/2 loss-induced epidermal defects elicit skin immune activation and alter global metabolism.

Skin epidermis constitutes the outer permeability barrier that protects the body from dehydration, heat loss, and myriad external assaults. Mechanisms that maintain barrier integrity in constantly challenged adult skin and how epidermal dysregulation shapes the local immune microenvironment and whole-body metabolism remain poorly understood. Here, we demonstrate that inducible and simultaneous ablation of transcription factor-encoding Ovol1 and Ovol2 in adult epidermis results in barrier dysregulation through impacting epithelial-mesenchymal plasticity and inflammatory gene expression. We find that aberrant skin immune activation then ensues, featuring Langerhans cell mobilization and T cell responses, and leading to elevated levels of secreted inflammatory factors in circulation. Finally, we identify failure to gain body weight and accumulate body fat as long-term consequences of epidermal-specific Ovol1/2 loss and show that these global metabolic changes along with the skin barrier/immune defects are partially rescued by immunosuppressant dexamethasone. Collectively, our study reveals key regulators of adult barrier maintenance and suggests a causal connection between epidermal dysregulation and whole-body metabolism that is in part mediated through aberrant immune activation.

DNA demethylase activity maintains intestinal cells in an undifferentiated state following loss of APC.

Although genome-wide hypomethylation is a hallmark of many cancers, roles for active DNA demethylation during tumorigenesis are unknown. Here, loss of the APC tumor suppressor gene causes upregulation of a DNA demethylase system and the concomitant hypomethylation of key intestinal cell fating genes. Notably, this hypomethylation maintained zebrafish intestinal cells in an undifferentiated state that was released upon knockdown of demethylase components. Mechanistically, the demethylase genes are directly activated by Pou5f1 and Cebpß and are indirectly repressed by retinoic acid, which antagonizes Pou5f1 and Cebpß. Apc mutants lack retinoic acid as a result of the transcriptional repression of retinol dehydrogenase l1 via a complex that includes Lef1, Groucho2, Ctbp1, Lsd1, and Corest. Our findings imply a model wherein APC controls intestinal cell fating through a switch in DNA methylation dynamics. Wild-type APC and retinoic acid downregulate demethylase components, thereby promoting DNA methylation of key genes and helping progenitors commit to differentiation.

TMEM16A inhibits autophagy and promotes the invasion of hypopharyngeal squamous cell carcinoma through mTOR pathway.

Previous studies have indicated that transmembrane protein 16A (TMEM16A) plays a crucial role in the pathogenesis and progression of various tumors by influencing multiple signaling pathways. However, the role of TMEM16A in regulating autophagy via the mammalian target of rapamycin (mTOR) pathway and its impact on the development of hypopharyngeal squamous cell carcinoma (HSCC) remain unclear. Immunohistochemistry and western blotting were used to assess the expression of TMEM16A in HSCC tissues and metastatic lymph nodes. Manipulation of TMEM16A expression levels was achieved in the FaDu cell line through overexpression or knockdown, followed by assessment of its biological effects using cell colony formation, wound healing, transwell, and invasion assays. Additionally, apoptosis and autophagy-related proteins, as well as autophagosome formation, were evaluated through western blotting, transmission electron microscopy, and immunofluorescence following TMEM16A knockdown or overexpression in FaDu cells. Our study revealed significantly elevated levels of TMEM16A in both HSCC tissues and metastatic lymph nodes compared to normal tissues. In vitro experiments demonstrated that silencing TMEM16A led to a notable suppression of HSCC cell proliferation, invasion, and migration. Furthermore, TMEM16A silencing effectively inhibited tumor growth in xenografted mice. Subsequent investigations indicated that knockdown of TMEM16A in HSCC cells could suppress mTOR activation, thereby triggering autophagic cell death by upregulating sequestosome-1 (SQSTM1/P62) and microtubule-associated protein light chain 3 II (LC3II). This study highlights the crucial role of TMEM16A in modulating autophagy in HSCC, suggesting its potential as a therapeutic target for the treatment of this malignancy.

A single cell RNA sequence atlas of the early Drosophila larval eye.

The Drosophila eye has been an important model to understand principles of differentiation, proliferation, apoptosis and tissue morphogenesis. However, a single cell RNA sequence resource that captures gene expression dynamics from the initiation of differentiation to the specification of different cell types in the larval eye disc is lacking. Here, we report transcriptomic data from 13,000 cells that cover six developmental stages of the larval eye. Our data show cell clusters that correspond to all major cell types present in the eye disc ranging from the initiation of the morphogenetic furrow to the differentiation of each photoreceptor cell type as well as early cone cells. We identify dozens of cell type-specific genes whose function in different aspects of eye development have not been reported. These single cell data will greatly aid research groups studying different aspects of early eye development and will facilitate a deeper understanding of the larval eye as a model system.

Cwc27, associated with retinal degeneration, functions as a splicing factor in vivo.

Previous in vitro studies indicate that CWC27 functions as a splicing factor in the Bact spliceosome complex, interacting with CWC22 to form a landing platform for eIF4A3, a core component of the exon junction complex. However, the function of CWC27 as a splicing factor has not been validated in any in vivo systems. CWC27 variants have been shown to cause autosomal recessive retinal degeneration, in both syndromic and non-syndromic forms. The Cwc27K338fs/K338fs mouse model was shown to have significant retinal dysfunction and degeneration by 6 months of age. In this report, we have taken advantage of the Cwc27K338fs/K338fs mouse model to show that Cwc27 is involved in splicing in vivo in the context of the retina. Bulk RNA and single cell RNA-sequencing of the mouse retina showed that there were gene expression and splicing pattern changes, including alternative splice site usage and intron retention. Positive staining for CHOP suggests that ER stress may be activated in response to the splicing pattern changes and is a likely contributor to the disease mechanism. Our results provide the first evidence that CWC27 functions as a splicing factor in an in vivo context. The splicing defects and gene expression changes observed in the Cwc27K338fs/K338fs mouse retina provide insight to the potential disease mechanisms, paving the way for targeted therapeutic development.

Loss of salt tolerance during tomato domestication conferred by variation in a Na+ /K+ transporter.

Domestication has resulted in reduced salt tolerance in tomato. To identify the genetic components causing this deficiency, we performed a genome-wide association study (GWAS) for root Na+ /K+ ratio in a population consisting of 369 tomato accessions with large natural variations. The most significant variations associated with root Na+ /K+ ratio were identified within the gene SlHAK20 encoding a member of the clade IV HAK/KUP/KT transporters. We further found that SlHAK20 transports Na+ and K+ and regulates Na+ and K+ homeostasis under salt stress conditions. A variation in the coding sequence of SlHAK20 was found to be the causative variant associated with Na+ /K+ ratio and confer salt tolerance in tomato. Knockout mutations in tomato SlHAK20 and the rice homologous genes resulted in hypersensitivity to salt stress. Together, our study uncovered a previously unknown molecular mechanism of salt tolerance responsible for the deficiency in salt tolerance in cultivated tomato varieties. Our findings provide critical information for molecular breeding to improve salt tolerance in tomato and other crops.

Loss-of-function variants in UBAP1L cause autosomal recessive retinal degeneration.

PURPOSE: Inherited retinal diseases (IRDs) are a group of monogenic conditions that can lead to progressive blindness. Their missing heritability is still considerable, due in part to the presence of disease genes that await molecular identification. The purpose of this work was to identify novel genetic associations with IRDs. METHODS: Patients underwent a comprehensive ophthalmological evaluation using standard-of-care tests, such as detailed retinal imaging (macular optical coherence tomography and short-wavelength fundus autofluorescence) and electrophysiological testing. Exome and genome sequencing, as well as computer-assisted data analysis were used for genotyping and detection of DNA variants. A minigene-driven splicing assay was performed to validate the deleterious effects of 1 of such variants. RESULTS: We identified 8 unrelated families from Hungary, the United States, Israel, and The Netherlands with members presenting with a form of autosomal recessive and nonsyndromic retinal degeneration, predominantly described as rod-cone dystrophy but also including cases of cone/cone-rod dystrophy. Age of disease onset was very variable, with some patients experiencing first symptoms during their fourth decade of life or later. Myopia greater than 5 diopters was present in 5 of 7 cases with available refractive data, and retinal detachment was reported in 2 cases. All ascertained patients carried biallelic loss-of-function variants in UBAP1L (HGNC: 40028), a gene with unknown function and with homologies to UBAP1, encoding a protein involved in ubiquitin metabolism. One of these pathogenic variants, the intronic NM_001163692.2:c.910-7G>A substitution, was identified in 5 unrelated families. Minigene-driven splicing assays in HEK293T cells confirmed that this DNA change is responsible for the creation of a new acceptor splice site, resulting in aberrant splicing. CONCLUSION: We identified UBAP1L as a novel IRD gene. Although its function is currently unknown, UBAP1L is almost exclusively expressed in photoreceptors and the retinal pigment epithelium, hence possibly explaining the link between pathogenic variants in this gene and an ocular phenotype.

10× single-cell sequencing revealed cellular composition heterogeneity in cardiac myxoma with malignant glandular properties.

Cardiac myxoma is the most common primary cardiac tumor in adults. The histogenesis and cellular composition of myxoma are still unclear. This study aims to reveal the role of myxoma cell components and their gene expression in tumor development. We obtained single living cells by enzymatic digestion of tissues from 4 cases of surgically resected cardiac myxoma. Of course, there was 1 case of glandular myxoma and 3 cases of nonglandular myxoma. Then, 10× single-cell sequencing was performed. We identified 12 types and 11 types of cell populations in glandular myxoma and nonglandular myxoma, respectively. Heterogeneous epithelial cells are the main components of glandular myxoma. The similarities and differences in T cells in both glandular and nonglandular myxoma were analyzed by KEGG and GO. The most important finding was that there was active communication between T cells and epithelial cells. These results clarify the possible tissue occurrence and heterogeneity of cardiac myxoma and provide a theoretical basis and guidance for clinical diagnosis and treatment.

A multicenter, randomized, double-blinded, placebo-controlled, phase â ¢ study evaluating the efficacy and safety of Xeligekimab (GR1501) in patients with moderate-to-severe plaque psoriasis.

BACKGROUND: Xeligekimab is a fully human monoclonal antibody that selectively neutralizes IL-17A and had shown potential efficacy in preliminary trials. OBJECTIVE: To evaluate the efficacy and safety of Xeligekimab in Chinese patients with moderate-to-severe psoriasis. METHODS: A total of 420 Chinese patients were randomized to 200 mg Xeligekimab every 2 weeks (n = 281) or placebo (n = 139) for the first 12 weeks, followed by extending the treatment schedule to GR1501 every 4 weeks for further 40 weeks. Efficacy was assessed by evaluating the Physician's Global Assessment (PGA) 0/1 and Psoriasis Area and Severity Index (PASI) 75/90/100 improvement. The safety profile was also evaluated. RESULTS: At week 12, The PASI 75/90/100 were achieved in 90.7%/74.4%/30.2%% patients in GR1501 group compared with 8.6%/1.4%/0% patients in placebo group, respectively. The PGA 0/1 were achieved in 74.4% patients of GR1501 group and 3.6% patients in placebo group, respectively. The PASI 75 and PGA 0/1 maintained until week 52. No unexpected adverse events were observed. CONCLUSION: Xeligekimab showed high efficacy and is well tolerated in Chinese patients with moderate-to-severe plaque psoriasis.

The dynamic role of nucleoprotein SHCBP1 in the cancer cell cycle and its potential as a synergistic target for DNA-damaging agents in cancer therapy.

BACKGROUND: Malignant tumours seriously threaten human life and health, and effective treatments for cancer are still being explored. The ability of SHC SH2 domain-binding protein 1 (SHCBP1) to induce cell cycle disturbance and inhibit tumour growth has been increasingly studied, but its dynamic role in the tumour cell cycle and corresponding effects leading to mitotic catastrophe and DNA damage have rarely been studied. RESULTS: In this paper, we found that the nucleoprotein SHCBP1 exhibits dynamic spatiotemporal expression during the tumour cell cycle, and SHCBP1 knockdown slowed cell cycle progression by inducing spindle disorder, as reflected by premature mitotic entry and multipolar spindle formation. This dysfunction was caused by G2/M checkpoint impairment mediated by downregulated WEE1 kinase and NEK7 (a member of the mammalian NIMA-related kinase family) expression and upregulated centromere/kinetochore protein Zeste White 10 (ZW10) expression. Moreover, both in vivo and in vitro experiments confirmed the significant inhibitory effects of SHCBP1 knockdown on tumour growth. Based on these findings, SHCBP1 knockdown in combination with low-dose DNA-damaging agents had synergistic tumouricidal effects on tumour cells. In response to this treatment, tumour cells were forced into the mitotic phase with considerable unrepaired DNA lesions, inducing mitotic catastrophe. These synergistic effects were attributed not only to the abrogation of the G2/M checkpoint and disrupted spindle function but also to the impairment of the DNA damage repair system, as demonstrated by mass spectrometry-based proteomic and western blotting analyses. Consistently, patients with low SHCBP1 expression in tumour tissue were more sensitive to radiotherapy. However, SHCBP1 knockdown combined with tubulin-toxic drugs weakened the killing effect of the drugs on tumour cells, which may guide the choice of chemotherapeutic agents in clinical practice. CONCLUSION: In summary, we elucidated the role of the nucleoprotein SHCBP1 in tumour cell cycle progression and described a novel mechanism by which SHCBP1 regulates tumour progression and through which targeting SHCBP1 increases sensitivity to DNA-damaging agent therapy, indicating its potential as a cancer treatment.

Children with severe neurological symptoms associated with SARS-CoV-2 infection during Omicron pandemic in China.

BACKGROUND: To analyze the clinical characteristics and outcomes of children with severe neurological symptoms associated with the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) infection during the Omicron pandemic in China. METHODS: This study used a questionnaire to obtain data from pediatric intensive care unit (PICU) centers in seven tertiary hospitals in Northeast China from December 1, 2022, to January 31, 2023. RESULTS: A total of 255 patients were confirmed to have SARS-CoV-2 infection, and 45 patients (17.65 %) were included in this study. Of these, seven (15.6%) patients died, and the median time from admission to death was 35 h (IQR, 14-120 h). Twenty (52.6%) survivors experienced neurological sequelae. Patients with platelet counts lower than 100 × 109/L had a higher incidence of complications such as multiple organ dysfunction, mechanical ventilation rate, and mortality. Cranial magnetic resonance imaging (MRI) always reveals cerebral tissue edema, with some severe lesions forming a softening site. CONCLUSION: Children infected with SARS-CoV-2 often exhibit severe neurological symptoms, and in some cases, they may rapidly develop malignant cerebral edema or herniation, leading to a fatal outcome. An early decrease in platelet count may associated with an unfavorable prognosis. IMPACT: Since early December 2022, China has gradually adjusted its prevention and control policy of SARS-CoV-2; Omicron outbreaks have occurred in some areas for a relatively short period. Due to the differences in ethnicity, endemic strains and vaccination status, there was a little difference from what has been reported about children with SARS-CoV-2 infection with severe neurological symptoms in abroad. This is the first multicenter clinical study in children with nervous system involvement after acute SARS-CoV-2 infection in China, and helpful for pediatricians to have a more comprehensive understanding of the clinical symptoms and prognosis of such disease.

Rapid and selective RRS determination of ferrocyanide with a nanogold surface molecularly imprinted polymethacrylic acid probe.

In this work, a nanogold surface molecularly imprinted polymer spectral probe (AuNP@MIP) for selectively identifying ferrocyanide was prepared under microwave irradiation using nanogold as the core, ferrocyanide as the template ion, methacrylic acid as the monomer, and ethylene glycol dimethacrylate as the cross-linking agent. AuNP@MIP was found to produce a resonance Rayleigh scattering (RRS) peak at 370 nm. When potassium ferrocyanide (K4Fe(CN)6) was present, a AuNP@MIP-Fe(CN)6 complex was formed, producing RRS-energy transfer (RRS-ET). With an increase in ferrocyanide concentration within a certain range, the RRS intensity at 370 nm decreased linearly, and the detection range was 0.02-0.40 µmol L-1, with a detection limit as low as 0.006 µmol L-1 ferrocyanide. This new method has the advantages of simplicity, rapidity, and selectivity when applied for the determination of K4Fe(CN)6 in table salt.

Atypical absence seizures and gene variants: A gene-based review of etiology, electro-clinical features, and associated epilepsy syndrome.

Atypical absence seizures are generalized non-convulsive seizures that often occur in children with cognitive impairment. They are common in refractory epilepsy and have been recognized as one of the hallmarks of developmental epileptic encephalopathies. Notably, pathogenic variants associated with AAS, such as GABRG2, GABRG3, SLC6A1, CACNB4, SCN8A, and SYNGAP1, are also linked to developmental epileptic encephalopathies. Atypical absences differ from typical absences in that they are frequently drug-resistant and the prognosis is dependent on the etiology or related epileptic syndromes. To improve clinicians' understanding of atypical absences and provide novel perspectives for clinical treatment, we have reviewed the electro-clinical characteristics, etiologies, treatment, and prognosis of atypical absences, with a focus on the etiology of advancements in gene variants, shedding light on potential avenues for improved clinical management.

Integrated analysis of whole genome and transcriptome sequencing uncovers genetic differences between Zi goose and Xianghai flying goose.

Zi goose is a famous indigenous breed originating from northeast China with high annual egg production. Xianghai flying goose is a composite breed and is bred by crosses of the wild swan goose and the Zi goose. Our previous study revealed significant differences in muscle fiber characteristics between the two populations. Here, we aimed to reveal the underlying genetic basis of the above phenotype differences through whole-genome and transcriptome analysis. A total of 20 blood samples (10 Zi geese and 10 Xianghai flying geese) were used for whole genome sequencing, and eight breast muscle tissue samples (four Zi geese and four Xianghai flying geese) were used for RNA sequencing. Using the FST and XP-EHH analysis, some highly differentiated genome regions annotated with egg production (RORB, WNT4, BMPR1B) and breast muscle development (WNT7B) between the two populations were detected. RNA-sequencing analysis revealed differentially expressed genes related to muscle development (IGF1, PAX7). Moreover, several genes were detected by both genome and transcriptome analysis, and some of them were reported to be associated with muscle growth (SLIT2, PREX1) and intramuscular fat (COL6A1). These findings will help researchers better understand the genetic basis related to egg production and muscle development in geese.

3'aQTL-atlas: an atlas of 3'UTR alternative polyadenylation quantitative trait loci across human normal tissues.

Genome-wide association studies (GWAS) have identified thousands of non-coding single-nucleotide polymorphisms (SNPs) associated with human traits and diseases. However, functional interpretation of these SNPs remains a significant challenge. Our recent study established the concept of 3' untranslated region (3'UTR) alternative polyadenylation (APA) quantitative trait loci (3'aQTLs), which can be used to interpret ∼16.1% of GWAS SNPs and are distinct from gene expression QTLs and splicing QTLs. Despite the growing interest in 3'aQTLs, there is no comprehensive database for users to search and visualize them across human normal tissues. In the 3'aQTL-atlas (https://wlcb.oit.uci.edu/3aQTLatlas), we provide a comprehensive list of 3'aQTLs containing ∼1.49 million SNPs associated with APA of target genes, based on 15,201 RNA-seq samples across 49 human Genotype-Tissue Expression (GTEx v8) tissues isolated from 838 individuals. The 3'aQTL-atlas provides a ∼2-fold increase in sample size compared with our published study. It also includes 3'aQTL searches by Gene/SNP across tissues, a 3'aQTL genome browser, 3'aQTL boxplots, and GWAS-3'aQTL colocalization event visualization. The 3'aQTL-atlas aims to establish APA as an emerging molecular phenotype to explain a large fraction of GWAS risk SNPs, leading to significant novel insights into the genetic basis of APA and APA-linked susceptibility genes in human traits and diseases.