Ensembl 2018

Investigating the relationship between depression and breast cancer: observational and genetic analyses

BACKGROUND

Both depression and breast cancer (BC) contribute to a substantial global burden of morbidity and mortality among women, and previous studies have observed a potential depression-BC link. We aimed to comprehensively characterize the phenotypic and genetic relationships between depression and BC.

METHODS

We first evaluated phenotypic association using longitudinal follow-up data from the UK Biobank (N = 250,294). We then investigated genetic relationships leveraging summary statistics from the hitherto largest genome-wide association study of European individuals conducted for depression (N = 500,199), BC (N = 247,173), and its subtypes based on the status of estrogen receptor (ER + : N = 175,475; ER - : N = 127,442).

RESULTS

Observational analysis suggested an increased hazard of BC in depression patients (HR = 1.10, 95%CIs = 0.95-1.26). A positive genetic correlation between depression and overall BC was observed ([Formula: see text] = 0.08, P = 3.00 × 10-4), consistent across ER + ([Formula: see text] = 0.06, P = 6.30 × 10-3) and ER - subtypes ([Formula: see text] = 0.08, P = 7.20 × 10-3). Several specific genomic regions showed evidence of local genetic correlation, including one locus at 9q31.2, and four loci at, or close, to 6p22.1. Cross-trait meta-analysis identified 17 pleiotropic loci shared between depression and BC. TWAS analysis revealed five shared genes. Bi-directional Mendelian randomization suggested risk of depression was causally associated with risk of overall BC (OR = 1.12, 95%Cis = 1.04-1.19), but risk of BC was not causally associated with risk of depression.

CONCLUSIONS

Our work demonstrates a shared genetic basis, pleiotropic loci, and a putative causal relationship between depression and BC, highlighting a biological link underlying the observed phenotypic relationship; these findings may provide important implications for future studies aimed reducing BC risk.

Association of HLA diversity with the risk of 25 cancers in the UK Biobank

BACKGROUND

The human leukocyte antigen (HLA) is a highly polymorphic region, and HLA diversity may play a role in presenting tumour-associated peptides and inducing immune responses. However, the effect of HLA diversity on cancers has not been fully assessed. We aimed to explore the role of HLA diversity on cancer development.

METHODS

A pan-cancer analysis was performed to evaluate the effect of HLA diversity, measured by HLA heterozygosity and HLA evolutionary divergence (HED), on the susceptibility of 25 cancers in the UK Biobank.

FINDINGS

We observed that the diversity of HLA class II locus was associated with a lower risk of lung cancer (OR_hetero = 0.94, 95% CI = 0.90-0.97, P = 1.29 × 10^-4) and head and neck cancer (OR_hetero = 0.91, 95% CI = 0.86-0.96, P = 1.56 × 10^-3). Besides, a lower risk of non-Hodgkin lymphoma was associated with an increased diversity of HLA class I (OR_hetero = 0.92, 95% CI = 0.87-0.98, P = 8.38 × 10^-3) and class II locus (OR_hetero = 0.89, 95% CI = 0.86-0.92, P = 1.65 × 10^-10). A lower risk of Hodgkin lymphoma was associated with the HLA class I diversity (OR_hetero = 0.85, 95% CI = 0.75-0.96, P = 0.011). The protective effect of HLA diversity was mainly observed in pathological subtypes with higher tumour mutation burden, such as lung squamous cell carcinoma (P = 9.39 × 10^-3) and diffuse large B cell lymphoma (P_{class I} = 4.12 × 10^-4; P_{class Ⅱ} = 4.71 × 10^-5), as well as the smoking subgroups of lung cancer (P = 7.45 × 10^-5) and head and neck cancer (P = 4.55 × 10^-3).

INTERPRETATION

We provided a systematic insight into the effect of HLA diversity on cancers, which might help to understand the etiological role of HLA on cancer development.

FUNDING

This study was supported by grants from the National Natural Science Foundation of China (82273705, 82003520); the Basic and Applied Basic Research Foundation of Guangdong Province, China (2021B1515420007); the Science and Technology Planning Project of Guangzhou, China (201804020094); Sino-Sweden Joint Research Programme (81861138006); the National Natural Science Foundation of China (81973131, 81903395, 81803319, 81802708).

Evaluating depressive symptoms, BDNF Val66Met, and APOE-ε4 as moderators of response to computerized cognitive training in heart failure

BACKGROUND

Depressive symptoms, brain-derived neurotrophic factor (BDNF) Val66Met, and apolipoprotein (APOE)-ε4 may moderate response to computerized cognitive training (CCT) interventions among patients with heart failure (HF).

OBJECTIVES

The purpose of this study was to examine moderators of intervention response to CCT over 8 months among patients with HF enrolled in a 3-arm randomized controlled trial. Outcomes were memory, serum BDNF, working memory, instrumental activities of daily living (IADLs), and health-related quality of life (HRQL).

METHODS

256 patients with HF were randomized to CCT, computerized crossword puzzles active control, and usual care control groups for 8 weeks. Data were collected at enrollment, baseline, 10 weeks, and 4 and 8 months. Mixed effects models were computed to evaluate moderators.

RESULTS

As previously reported, there were no statistically significant group by time effects in outcomes among the 3 groups over 8 months. Tests of moderation indicated that depressive symptoms and presence of BDNF Val66Met and APOE-ε4 were not statistically significant moderators of intervention response in outcomes of delayed recall memory, serum BDNF, working memory, IADLs, and HRQL. In post hoc analysis evaluating baseline global cognitive function, gender, age, and HF severity as moderators, no significant effects were found. HF severity was imbalanced among groups (P = .049) which may have influenced results.

CONCLUSIONS

Studies are needed to elucidate biological mechanisms of cognitive dysfunction in HF and test novel interventions to improve memory, serum BDNF, working memory, IADLs and HRQL. Patients may need to be stratified or randomized by HF severity within intervention trials.

Antigen-specific B cells direct T follicular-like helper cells into lymphoid follicles to mediate Mycobacterium tuberculosis control

Tuberculosis (TB), caused by Mycobacterium tuberculosis (Mtb), is a global cause of death. Granuloma-associated lymphoid tissue (GrALT) correlates with protection during TB, but the mechanisms of protection are not understood. During TB, the transcription factor IRF4 in T cells but not B cells is required for the generation of the TH1 and TH17 subsets of helper T cells and follicular helper T (TFH)-like cellular responses. A population of IRF4+ T cells coexpress the transcription factor BCL6 during Mtb infection, and deletion of Bcl6 (Bcl6fl/fl) in CD4+ T cells (CD4cre) resulted in reduction of TFH-like cells, impaired localization within GrALT and increased Mtb burden. In contrast, the absence of germinal center B cells, MHC class II expression on B cells, antibody-producing plasma cells or interleukin-10-expressing B cells, did not increase Mtb susceptibility. Indeed, antigen-specific B cells enhance cytokine production and strategically localize TFH-like cells within GrALT via interactions between programmed cell death 1 (PD-1) and its ligand PD-L1 and mediate Mtb control in both mice and macaques.

Progranulin deficiency results in sex-dependent alterations in microglia in response to demyelination

Heterozygous mutations in the granulin (GRN) gene, resulting in the haploinsufficiency of the progranulin (PGRN) protein, is a leading cause of frontotemporal lobar degeneration (FTLD). Complete loss of the PGRN protein causes neuronal ceroid lipofuscinosis (NCL), a lysosomal storage disorder. Polymorphisms in the GRN gene have also been associated with several other neurodegenerative diseases, including Alzheimer's disease (AD), and Parkinson's disease (PD). PGRN deficiency has been shown to cause myelination defects previously, but how PGRN regulates myelination is unknown. Here, we report that PGRN deficiency leads to a sex-dependent myelination defect with male mice showing more severe demyelination in response to cuprizone treatment. This is accompanied by exacerbated microglial proliferation and activation in the male PGRN-deficient mice. Interestingly, both male and female PGRN-deficient mice show sustained microglial activation after cuprizone removal and a defect in remyelination. Specific ablation of PGRN in microglia results in similar sex-dependent phenotypes, confirming a microglial function of PGRN. Lipid droplets accumulate in microglia specifically in male PGRN-deficient mice. RNA-seq analysis and mitochondrial function assays reveal key differences in oxidative phosphorylation in male versus female microglia under PGRN deficiency. A significant decrease in myelination and accumulation of myelin debris and lipid droplets in microglia were found in the corpus callosum regions of FTLD patients with GRN mutations. Taken together, our data support that PGRN deficiency leads to sex-dependent alterations in microglia with subsequent myelination defects.

Local genetic correlations exist among neurodegenerative and neuropsychiatric diseases

Genetic correlation ([Formula: see text]) between traits can offer valuable insight into underlying shared biological mechanisms. Neurodegenerative diseases overlap neuropathologically and often manifest comorbid neuropsychiatric symptoms. However, global [Formula: see text] analyses show minimal [Formula: see text] among neurodegenerative and neuropsychiatric diseases. Importantly, local [Formula: see text] s can exist in the absence of global relationships. To investigate this possibility, we applied LAVA, a tool for local [Formula: see text] analysis, to genome-wide association studies of 3 neurodegenerative diseases (Alzheimer's disease, Lewy body dementia and Parkinson's disease) and 3 neuropsychiatric disorders (bipolar disorder, major depressive disorder and schizophrenia). We identified several local [Formula: see text] s missed in global analyses, including between (i) all 3 neurodegenerative diseases and schizophrenia and (ii) Alzheimer's and Parkinson's disease. For those local [Formula: see text] s identified in genomic regions containing disease-implicated genes, such as SNCA, CLU and APOE, incorporation of expression quantitative trait loci identified genes that may drive genetic overlaps between diseases. Collectively, we demonstrate that complex genetic relationships exist among neurodegenerative and neuropsychiatric diseases, highlighting putative pleiotropic genomic regions and genes. These findings imply sharing of pathogenic processes and the potential existence of common therapeutic targets.

ASAS-EGB: A statistical framework for estimating allele-specific alternative splicing events using transcriptome data

Ubiquitous Alternative splicing (AS) in eukaryotes greatly increases the biodiversity of proteins and involves in disease and cancer. Allele-specific AS studies can facilitate the identification of cis-acting elements because both alleles share the same cellular environment. Due to the limited information provided on the exons defined by AS events, we propose a statistical framework and algorithm ASAS-EGB for ASAS analysis using the gene transcriptome. The framework obtains exclusively compatible sets of gene isoforms supporting each event isoform, and uses both phased and non-phased SNPs within the exons on the gene isoforms for inference. Using this strategy, we have demonstrated ASAS-EGB can yield better ASAS inferential performance than using event isoforms. ASAS-EGB supports both single-end and paired-end RNA-seq data, and we have proved its robustness using RNA-seq replicates of individual NA12878. ASAS-EGB builds Bayesian models for ASAS analysis, and the MCMC method is used to solve the problem. With more detailed annotations for individual genomes and transcriptomes appearing in the future, the algorithm proposed by the paper can provide better support for these data to reveal the regulatory mechanisms of individual genomes.

A burst of genomic innovation at the origin of placental mammals mediated embryo implantation

The origin of embryo implantation in mammals ~148 million years ago was a dramatic shift in reproductive strategy, yet the molecular changes that established mammal implantation are largely unknown. Although progesterone receptor signalling predates the origin of mammals and is highly conserved in, and critical for, successful mammal pregnancy, it alone cannot explain the origin and subsequent diversity of implantation strategies throughout the placental mammal radiation. MiRNAs are known to be flexible and dynamic regulators with a well-established role in the pathophysiology of mammal placenta. We propose that a dynamic core microRNA (miRNA) network originated early in placental mammal evolution, responds to conserved mammal pregnancy cues (e.g. progesterone), and facilitates species-specific responses. Here we identify 13 miRNA gene families that arose at the origin of placental mammals and were subsequently retained in all descendent lineages. The expression of these miRNAs in response to early pregnancy molecules is regulated in a species-specific manner in endometrial epithelia of species with extreme implantation strategies (i.e. bovine and human). Furthermore, this set of miRNAs preferentially target proteins under positive selective pressure on the ancestral eutherian lineage. Discovery of this core embryo implantation toolkit and specifically adapted proteins helps explain the origin and evolution of implantation in mammals.

Single-Cell Genomics Reveals the Divergent Mitochondrial Genomes of Retaria (Foraminifera and Radiolaria)

Mitochondria originated from an ancient bacterial endosymbiont that underwent reductive evolution by gene loss and endosymbiont gene transfer to the nuclear genome. The diversity of mitochondrial genomes published to date has revealed that gene loss and transfer processes are ongoing in many lineages. Most well-studied eukaryotic lineages are represented in mitochondrial genome databases, except for the superphylum Retaria-the lineage comprising Foraminifera and Radiolaria. Using single-cell approaches, we determined two complete mitochondrial genomes of Foraminifera and two nearly complete mitochondrial genomes of radiolarians. We report the complete coding content of an additional 14 foram species. We show that foraminiferan and radiolarian mitochondrial genomes contain a nearly fully overlapping but reduced mitochondrial gene complement compared to other sequenced rhizarians. In contrast to animals and fungi, many protists encode a diverse set of proteins on their mitochondrial genomes, including several ribosomal genes; however, some aerobic eukaryotic lineages (euglenids, myzozoans, and chlamydomonas-like algae) have reduced mitochondrial gene content and lack all ribosomal genes. Similar to these reduced outliers, we show that retarian mitochondrial genomes lack ribosomal protein and tRNA genes, contain truncated and divergent small and large rRNA genes, and contain only 14 or 15 protein-coding genes, including nad1, -3, -4, -4L, -5, and -7, cob, cox1, -2, and -3, and atp1, -6, and -9, with forams and radiolarians additionally carrying nad2 and nad6, respectively. In radiolarian mitogenomes, a noncanonical genetic code was identified in which all three stop codons encode amino acids. Collectively, these results add to our understanding of mitochondrial genome evolution and fill in one of the last major gaps in mitochondrial sequence databases. IMPORTANCE We present the reduced mitochondrial genomes of Retaria, the rhizarian lineage comprising the phyla Foraminifera and Radiolaria. By applying single-cell genomic approaches, we found that foraminiferan and radiolarian mitochondrial genomes contain an overlapping but reduced mitochondrial gene complement compared to other sequenced rhizarians. An alternative genetic code was identified in radiolarian mitogenomes in which all three stop codons encode amino acids. Collectively, these results shed light on the divergent nature of the mitochondrial genomes from an ecologically important group, warranting further questions into the biological underpinnings of gene content variability and genetic code variation between mitochondrial genomes.

Effective extraction of polyribosomes exposes gene expression strategies in primary astrocytes

Regulation of mRNA translation in astrocytes gains a growing interest. However, until now, successful ribosome profiling of primary astrocytes has not been reported. Here, we optimized the standard 'polysome profiling' method and generated an effective protocol for polyribosome extraction, which enabled genome-wide assessment of mRNA translation dynamics along the process of astrocyte activation. Transcriptome (RNAseq) and translatome (Riboseq) data generated at 0, 24 and 48 h after cytokines treatment, revealed dynamic genome-wide changes in the expression level of ∼12 000 genes. The data clarify whether a change in protein synthesis rate results from a change in mRNA level or translation efficiency per se. It exhibit different expression strategies, based on changes in mRNA abundance and/or translation efficiency, which are specifically assigned to gene subsets depending on their function. Moreover, the study raises an important take-home message related to the possible presence of 'difficult to extract' polyribosome sub-groups, in all cell types, thus illuminating the impact of ribosomes extraction methodology on experiments addressing translation regulation.

Assessing Kidney Injury Induced by Mercuric Chloride in Guinea Pigs with In Vivo and In Vitro Experiments

Acute kidney injury, which is associated with high levels of morbidity and mortality, affects a significant number of individuals, and can be triggered by multiple factors, such as medications, exposure to toxic chemicals or other substances, disease, and trauma. Because the kidney is a critical organ, understanding and identifying early cellular or gene-level changes can provide a foundation for designing medical interventions. In our earlier work, we identified gene modules anchored to histopathology phenotypes associated with toxicant-induced liver and kidney injuries. Here, using in vivo and in vitro experiments, we assessed and validated these kidney injury-associated modules by analyzing gene expression data from the kidneys of male Hartley guinea pigs exposed to mercuric chloride. Using plasma creatinine levels and cell-viability assays as measures of the extent of renal dysfunction under in vivo and in vitro conditions, we performed an initial range-finding study to identify the appropriate doses and exposure times associated with mild and severe kidney injuries. We then monitored changes in kidney gene expression at the selected doses and time points post-toxicant exposure to characterize the mechanisms of kidney injury. Our injury module-based analysis revealed a dose-dependent activation of several phenotypic cellular processes associated with dilatation, necrosis, and fibrogenesis that were common across the experimental platforms and indicative of processes that initiate kidney damage. Furthermore, a comparison of activated injury modules between guinea pigs and rats indicated a strong correlation between the modules, highlighting their potential for cross-species translational studies.

NAD+ Metabolism Generates a Metabolic Vulnerability in Endocrine-Resistant Metastatic Breast Tumors in Females

Approximately 70% of human breast cancers express estrogen receptor-α (ERα), providing a potential target for endocrine therapy. However, 30% to 40% of patients with ER+ breast cancer still experience recurrence and metastasis, with a 5-year relative overall survival rate of 24%. In this study, we identified nicotinamide phosphoribosyltransferase (NAMPT), an important enzyme in nicotinamide adenine dinucleotide (NAD+) metabolism, to be increased in metastatic breast cancer (MBC) cells treated with fulvestrant (Fulv). We tested whether the blockade of NAD+ production via inhibition of NAMPT synergizes with standard-of-care therapies for ER+ MBC in vitro and in vivo. A synergistic effect was not observed when KPT-9274 was combined with palbociclib or tamoxifen or when Fulv was combined with other metabolic inhibitors. We show that NAMPT inhibitor KPT-9274 and Fulv works synergistically to reduce metastatic tumor burden. RNA-sequencing analysis showed that NAMPT inhibitor in combination with Fulv reversed the expression of gene sets associated with more aggressive tumor phenotype, and metabolomics analysis showed that NAMPT inhibition reduced the abundance of metabolites associated with several key tumor metabolic pathways. Targeting metabolic adaptations in endocrine-resistant MBC is a novel strategy, and alternative approaches aimed at improving the therapeutic response of metastatic ER+ tumors are needed. Our findings uncover the role of ERα-NAMPT crosstalk in MBC and the utility of NAMPT inhibition and antiestrogen combination therapy in reducing tumor burden and metastasis, potentially leading to new avenues of MBC treatment.

Improved analysis of (e)CLIP data with RCRUNCH yields a compendium of RNA-binding protein binding sites and motifs

We present RCRUNCH, an end-to-end solution to CLIP data analysis for identification of binding sites and sequence specificity of RNA-binding proteins. RCRUNCH can analyze not only reads that map uniquely to the genome but also those that map to multiple genome locations or across splice boundaries and can consider various types of background in the estimation of read enrichment. By applying RCRUNCH to the eCLIP data from the ENCODE project, we have constructed a comprehensive and homogeneous resource of in-vivo-bound RBP sequence motifs. RCRUNCH automates the reproducible analysis of CLIP data, enabling studies of post-transcriptional control of gene expression.

Domain landscapes of somatic NF1 mutations in pheochromocytoma and paraganglioma

Pheochromocytoma and paraganglioma (PPGL), are rare neuroendocrine tumors arising from the adrenal medulla and extra-adrenal paraganglia, respectively. Up to about 60% are explained by germline or somatic mutations in one of the major known susceptibility genes e.g., inNF1,RET,VHL, SDHx,MAXandHRAS. Targeted Next Generation Sequencing was performed in 14 sporadic tumors using a panel including 26 susceptibility genes to characterize the mutation profile. A total of 6 germline and 8 somatic variants were identified. The most frequent somatic mutations were found in NF1(36%), four have not been reported earlier in PCC or PGL. Gene expression profile analysis showed that NF1 mutated tumors are classified into RTK3 subtype, cluster 2, with a high expression of genes associated with chromaffin cell differentiation, and into a RTK2 subtype, cluster 2, as well with overexpression of genes associated with cortisol biosynthesis. On the other hand, by analyzing the entire probe set on the array and TCGA data, ALDOC was found as the most significantly down regulated gene in NF1-mutated tumors compared to NF1-wild-type. Differential gene expression analysis showed a significant difference between Nt - and Ct-NF1 domains in mutated tumors probably engaging different cellular pathways. Notably, we had a metastatic PCC with a Ct-NF1 frameshift mutation and when performing protein docking analysis, Ct-NF1 showed an interaction with Nt-FAK suggesting their involvement in cell adhesion and cell growth. These results show that depending on the location of the NF1-mutation different pathways are activated in PPGLs. Further studies are required to clarify their clinical significance.

Characterization of genome-wide STR variation in 6487 human genomes

Short tandem repeats (STRs) are abundant and highly mutagenic in the human genome. Many STR loci have been associated with a range of human genetic disorders. However, most population-scale studies on STR variation in humans have focused on European ancestry cohorts or are limited by sequencing depth. Here, we depicted a comprehensive map of 366,013 polymorphic STRs (pSTRs) constructed from 6487 deeply sequenced genomes, comprising 3983 Chinese samples (~31.5x, NyuWa) and 2504 samples from the 1000 Genomes Project (~33.3x, 1KGP). We found that STR mutations were affected by motif length, chromosome context and epigenetic features. We identified 3273 and 1117 pSTRs whose repeat numbers were associated with gene expression and 3'UTR alternative polyadenylation, respectively. We also implemented population analysis, investigated population differentiated signatures, and genotyped 60 known disease-causing STRs. Overall, this study further extends the scale of STR variation in humans and propels our understanding of the semantics of STRs.

The targeted next-generation sequence revealed SMAD4, AKT1, and TP53 mutations from circulating cell-free DNA of breast cancer and its effect on protein structure - A computational approach

Breast cancer biomarkers that detect marginally advanced stages are still challenging. The detection of specific abnormalities, targeted therapy selection, prognosis, and monitoring of treatment effectiveness over time are all made possible by circulating free DNA (cfDNA) analysis. The proposed study will detect specific genetic abnormalities from the plasma cfDNA of a female breast cancer patient by sequencing a cancer-related gene panel (MGM455 - Oncotrack Ultima), including 56 theranostic genes (SNVs and small INDELs). Initially, we determined the pathogenicity of the observed mutations using PredictSNP, iStable, Align-GVGD, and ConSurf servers. As a next step, molecular dynamics (MD) was implemented to determine the functional significance of SMAD4 mutation (V465M). Lastly, the mutant gene relationships were examined using the Cytoscape plug-in GeneMANIA. Using ClueGO, we determined the gene's functional enrichment and integrative analysis. The structural characteristics of SMAD4 V465M protein by MD simulation analysis further demonstrated that the mutation was deleterious. The simulation showed that the native structure was more significantly altered by the SMAD4 (V465M) mutation. Our findings suggest that SMAD4 V465M mutation might be significantly associated with breast cancer, and other patient-found mutations (AKT1-E17K and TP53-R175H) are synergistically involved in the process of SMAD4 translocate to nuclease, which affects the target gene translation. Therefore, this combination of gene mutations could alter the TGF-β signaling pathway in BC. We further proposed that the SMAD4 protein loss may contribute to an aggressive phenotype by inhibiting the TGF-β signaling pathway. Thus, breast cancer's SMAD4 (V465M) mutation might increase their invasive and metastatic capabilities.Communicated by Ramaswamy H. Sarma.

11-Cl-BBQ, a select modulator of AhR-regulated transcription, suppresses lung cancer cell growth via activation of p53 and p27Kip1

Aryl hydrocarbon receptor (AhR) is a ligand-activated transcription factor and functions as a tumour suppressor in different cancer models. In the present study, we report detailed characterization of 11-chloro-7H-benzimidazo[2,1-a]benzo[de]iso-quinolin-7-one (11-Cl-BBQ) as a select modulator of AhR-regulated transcription (SMAhRT) with anti-cancer actions. Treatment of lung cancer cells with 11-Cl-BBQ induced potent and sustained AhR-dependent anti-proliferative effects by promoting G1 phase cell cycle arrest. Investigation of 11-Cl-BBQ-induced transcription in H460 cells with or without the AhR expression by RNA-sequencing revealed activation of p53 signalling. In addition, 11-Cl-BBQ suppressed multiple pathways involved in DNA replication and increased expression of cyclin-dependent kinase inhibitors, including p27Kip1 , in an AhR-dependent manner. CRISPR/Cas9 knockout of individual genes revealed the requirement for both p53 and p27Kip1 for the AhR-mediated anti-proliferative effects. Our results identify 11-Cl-BBQ as a potential lung cancer therapeutic, highlight the feasibility of targeting AhR and provide important mechanistic insights into AhR-mediated-anticancer actions.

Multiple ubiquitin E3 ligase genes antagonistically regulate chloroplast-associated protein degradation

The chloroplast is the most prominent member of a diverse group of plant organelles called the plastids, and it is characterized by its vital role in photosynthesis.1,2,3 Most of the ∼3,000 different proteins in chloroplasts are synthesized in the cytosol in precursor (preprotein) form, each with a cleavable transit peptide.4,5,6,7,8 Preproteins are imported via translocons in the outer and inner envelope membranes of the chloroplast, termed TOC and TIC, respectively.9,10,11,12,13 Discovery of the chloroplast-localized ubiquitin E3 ligase SUPPRESSOR OF PPI1 LOCUS1 (SP1) demonstrated that the nucleocytosolic ubiquitin-proteasome system (UPS) targets the TOC apparatus to dynamically control protein import and chloroplast biogenesis in response to developmental and environmental cues. The relevant UPS pathway is termed chloroplast-associated protein degradation (CHLORAD).14,15,16 Two homologs of SP1 exist, SP1-like1 (SPL1) and SPL2, but their roles have remained obscure. Here, we show that SP1 is ubiquitous in the Viridiplantae and that SPL2 and SPL1 appeared early during the evolution of the Viridiplantae and land plants, respectively. Through genetic and biochemical analysis, we reveal that SPL1 functions as a negative regulator of SP1, potentially by interfering with its ability to catalyze ubiquitination. In contrast, SPL2, the more distantly related SP1 homolog, displays partial functional redundancy with SP1. Both SPL1 and SPL2 modify the extent of leaf senescence, like SP1, but do so in diametrically opposite ways. Thus, SPL1 and SPL2 are bona fide CHLORAD system components with negative and positive regulatory functions that allow for nuanced control of this vital proteolytic pathway.

Comparative Transcriptomic Analysis Reveals the Functionally Segmented Intestine in Tunicate Ascidian

The vertebrate intestinal system consists of separate segments that remarkably differ in morphology and function. However, the origin of intestinal segmentation remains unclear. In this study, we investigated the segmentation of the intestine in a tunicate ascidian species, Ciona savignyi, by performing RNA sequencing. The gene expression profiles showed that the whole intestine was separated into three segments. Digestion, ion transport and signal transduction, and immune-related pathway genes were enriched in the proximal, middle, and distal parts of the intestine, respectively, implying that digestion, absorption, and immune function appear to be regional specializations in the ascidian intestine. We further performed a multi-species comparison analysis and found that the Ciona intestine showed a similar gene expression pattern to vertebrates, indicating tunicates and vertebrates might share the conserved intestinal functions. Intriguingly, vertebrate pancreatic homologous genes were expressed in the digestive segment of the Ciona intestine, suggesting that the proximal intestine might play the part of pancreatic functions in C. savignyi. Our results demonstrate that the tunicate intestine can be functionally separated into three distinct segments, which are comparable to the corresponding regions of the vertebrate intestinal system, offering insights into the functional evolution of the digestive system in chordates.

Integrative comparative analysis of avian chromosome evolution by in-silico mapping of the gene ontology of homologous synteny blocks and evolutionary breakpoint regions

Avian chromosomes undergo more intra- than interchromosomal rearrangements, which either induce or are associated with genome variations among birds. Evolving from a common ancestor with a karyotype not dissimilar from modern chicken, two evolutionary elements characterize evolutionary change: homologous synteny blocks (HSBs) constitute common conserved parts at the sequence level, while evolutionary breakpoint regions (EBRs) occur between HSBs, defining the points where rearrangement occurred. Understanding the link between the structural organization and functionality of HSBs and EBRs provides insight into the mechanistic basis of chromosomal change. Previously, we identified gene ontology (GO) terms associated with both; however, here we revisit our analyses in light of newly developed bioinformatic algorithms and the chicken genome assembly galGal6. We aligned genomes available for six birds and one lizard species, identifying 630 HSBs and 19 EBRs. We demonstrate that HSBs hold vast functionality expressed by GO terms that have been largely conserved through evolution. Particularly, we found that genes within microchromosomal HSBs had specific functionalities relevant to neurons, RNA, cellular transport and embryonic development, and other associations. Our findings suggest that microchromosomes may have conserved throughout evolution due to the specificity of GO terms within their HSBs. The detected EBRs included those found in the genome of the anole lizard, meaning they were shared by all saurian descendants, with others being unique to avian lineages. Our estimate of gene richness in HSBs supported the fact that microchromosomes contain twice as many genes as macrochromosomes.

Natural selection and convergent evolution of the HOX gene family in Carnivora

HOX genes play a central role in the development and regulation of limb patterns. For mammals in the order Carnivora, limbs have evolved in different forms, and there are interesting cases of phenotypic convergence, such as the pseudothumb of the giant and red pandas, and the flippers or specialized limbs of the pinnipeds and sea otter. However, the molecular bases of limb development remain largely unclear. Here, we studied the molecular evolution of the HOX9 ~ 13 genes of 14 representative species in Carnivora and explored the molecular evolution of other HOX genes. We found that only one limb development gene, HOXC10, underwent convergent evolution between giant and red pandas and was thus an important candidate gene related to the development of pseudothumbs. No signals of amino acid convergence and natural selection were found in HOX9 ~ 13 genes between pinnipeds and sea otter, but there was evidence of positive selection and rapid evolution in four pinniped species. Overall, few HOX genes evolve via natural selection or convergent evolution, and these could be important candidate genes for further functional validation. Our findings provide insights into potential molecular mechanisms of the development of specialized pseudothumbs and flippers (or specialized limbs).

Rani ZZ, Sulaiman SA, Chow YP, Abdullah N, Ahmad N, Ismail N, Abdul Jalal N, Kamaruddin MA, Saperi AA, Jamal R. Hypercholesterolemia in the Malaysian Cohort Participants: Genetic and Non-Genetic Risk Factors

Hypercholesterolemia was prevalent in 44.9% of The Malaysian Cohort participants, of which 51% were Malay. This study aimed to identify the variants involved in hypercholesterolemia among Malays and to determine the association between genetic and non-genetic risk factors. This nested case–control study included 25 Malay participants with the highest low-density lipoprotein cholesterol (LDL-C, >4.9 mmol/L) and total cholesterol (TC, >7.5 mmol/L) and 25 participants with the lowest LDL-C/TC. Genomic DNA was extracted, and whole-exome sequencing was performed using the Ion ProtonTM system. All variants were annotated, filtered, and cross-referenced against publicly available databases. Forty-five selected variants were genotyped in 677 TMC Malay participants using the MassARRAY® System. The association between genetic and non-genetic risk factors was determined using logistic regression analysis. Age, fasting blood glucose, tobacco use, and family history of hyperlipidemia were significantly associated with hypercholesterolemia. Participants with the novel OSBPL7 (oxysterol-binding protein-like 7) c.651_652del variant had 17 times higher odds for hypercholesterolemia. Type 2 diabetes patients on medication and those with PCSK9 (proprotein convertase subtilisin/kexin type 9) rs151193009 had low odds for hypercholesterolemia. Genetic predisposition can interact with non-genetic factors to increase hypercholesterolemia risk in Malaysian Malays.

8-Oxoguanine DNA Glycosylase 1 Upregulation as a Risk Factor for Obesity and Colorectal Cancer

DNA damage has been extensively studied as a potentially helpful tool in assessing and preventing cancer, having been widely associated with the deregulation of DNA damage repair (DDR) genes and with an increased risk of cancer. Adipose tissue and tumoral cells engage in a reciprocal interaction to establish an inflammatory microenvironment that enhances cancer growth by modifying epigenetic and gene expression patterns. Here, we hypothesize that 8-oxoguanine DNA glycosylase 1 (OGG1)—a DNA repair enzyme—may represent an attractive target that connects colorectal cancer (CRC) and obesity. In order to understand the mechanisms underlying the development of CRC and obesity, the expression and methylation of DDR genes were analyzed in visceral adipose tissue from CRC and healthy participants. Gene expression analysis revealed an upregulation of OGG1 expression in CRC participants (p < 0.005) and a downregulation of OGG1 in normal-weight healthy patients (p < 0.05). Interestingly, the methylation analysis showed the hypermethylation of OGG1 in CRC patients (p < 0.05). Moreover, expression patterns of OGG1 were found to be regulated by vitamin D and inflammatory genes. In general, our results showed evidence that OGG1 can regulate CRC risk through obesity and may act as a biomarker for CRC.

The acetabular labrum tissue shows unique transcriptome signatures compared to cartilage and responds to combined cyclic compression and surface shearing

The labrum is a fibrocartilaginous ring surrounding the acetabulum. Loss of labrum function contributes to the degeneration of the hip joint, leading to osteoarthritis. Successful labrum restoration requires profound knowledge about the tissue being replaced. The aim of this study was to characterize the transcriptome and the mechanobiological function of the labrum. RNA-seq was performed to compare the transcriptome of bovine labrum against articular cartilage tissue. Differential expression and gene ontology (GO) term pathway analysis were applied using the SUSHI framework. Bovine labrum explants were cultured for 5 days with / without mechanical loading and targeted gene expression was analyzed by real time quantitative polymerase chain reaction. More than 6'000 genes were significantly differentially expressed in the labrum compared to cartilage. Up- and downregulated genes were associated with the GO term extracellular matrix organization. The study established an extracellular matrix gene expression profile of healthy labrum tissue and identified significantly upregulated extracellular matrix related genes compared to cartilage tissue. Mechanical loading significantly upregulated aggrecan (ACAN), cartilage oligomeric matrix protein (COMP), fibronectin (FN1) and proteoglycan 4 (PRG4). MMP1/3/9 and IL6, which were upregulated by an inflammatory stimulus (IL-1b), were statistically unaffected by the loading, although IL6 was upregulated in each donor immediately after the loading. Unique ECM related features may guide the development of labrum tissue-engineering solutions. Despite the transcriptome differences between labrum and cartilage tissue, gene expression response to mechanical loading showed similarities with previously reported responses in cartilage, indicating a preserved tissue adaptation mechanism to mechanical loading. Running title: Acetabular Labrum Mechanobiology.

Higher-order modular regulation of the human proteome

Operons are transcriptional modules that allow bacteria to adapt to environmental changes by coordinately expressing the relevant set of genes. In humans, biological pathways and their regulation are more complex. If and how human cells coordinate the expression of entire biological processes is unclear. Here, we capture 31 higher-order co-regulation modules, which we term progulons, by help of supervised machine-learning on proteomics data. Progulons consist of dozens to hundreds of proteins that together mediate core cellular functions. They are not restricted to physical interactions or co-localisation. Progulon abundance changes are primarily controlled at the level of protein synthesis and degradation. Implemented as a web app at www.proteomehd.net/progulonFinder, our approach enables the targeted search for progulons of specific cellular processes. We use it to identify a DNA replication progulon and reveal multiple new replication factors, validated by extensive phenotyping of siRNA-induced knockdowns. Progulons provide a new entry point into the molecular understanding of biological processes.

Engineering an inhibitor-resistant human CSF1R variant for microglia replacement

Hematopoietic stem cell transplantation (HSCT) can replace endogenous microglia with circulation-derived macrophages but has high mortality. To mitigate the risks of HSCT and expand the potential for microglia replacement, we engineered an inhibitor-resistant CSF1R that enables robust microglia replacement. A glycine to alanine substitution at position 795 of human CSF1R (G795A) confers resistance to multiple CSF1R inhibitors, including PLX3397 and PLX5622. Biochemical and cell-based assays show no discernable gain or loss of function. G795A- but not wildtype-CSF1R expressing macrophages efficiently engraft the brain of PLX3397-treated mice and persist after cessation of inhibitor treatment. To gauge translational potential, we CRISPR engineered human-induced pluripotent stem cell-derived microglia (iMG) to express G795A. Xenotransplantation studies demonstrate that G795A-iMG exhibit nearly identical gene expression to wildtype iMG, respond to inflammatory stimuli, and progressively expand in the presence of PLX3397, replacing endogenous microglia to fully occupy the brain. In sum, we engineered a human CSF1R variant that enables nontoxic, cell type, and tissue-specific replacement of microglia.

Exploring microbial functional biodiversity at the protein family level-From metagenomic sequence reads to annotated protein clusters

Metagenomics has enabled accessing the genetic repertoire of natural microbial communities. Metagenome shotgun sequencing has become the method of choice for studying and classifying microorganisms from various environments. To this end, several methods have been developed to process and analyze the sequence data from raw reads to end-products such as predicted protein sequences or families. In this article, we provide a thorough review to simplify such processes and discuss the alternative methodologies that can be followed in order to explore biodiversity at the protein family level. We provide details for analysis tools and we comment on their scalability as well as their advantages and disadvantages. Finally, we report the available data repositories and recommend various approaches for protein family annotation related to phylogenetic distribution, structure prediction and metadata enrichment.

Zebra-Sphinx: Modeling Sphingolipidoses in Zebrafish

Sphingolipidoses are inborn errors of metabolism due to the pathogenic mutation of genes that encode for lysosomal enzymes, transporters, or enzyme cofactors that participate in the sphingolipid catabolism. They represent a subgroup of lysosomal storage diseases characterized by the gradual lysosomal accumulation of the substrate(s) of the defective proteins. The clinical presentation of patients affected by sphingolipid storage disorders ranges from a mild progression for some juvenile- or adult-onset forms to severe/fatal infantile forms. Despite significant therapeutic achievements, novel strategies are required at basic, clinical, and translational levels to improve patient outcomes. On these bases, the development of in vivo models is crucial for a better understanding of the pathogenesis of sphingolipidoses and for the development of efficacious therapeutic strategies. The teleost zebrafish (Danio rerio) has emerged as a useful platform to model several human genetic diseases owing to the high grade of genome conservation between human and zebrafish, combined with precise genome editing and the ease of manipulation. In addition, lipidomic studies have allowed the identification in zebrafish of all of the main classes of lipids present in mammals, supporting the possibility to model diseases of the lipidic metabolism in this animal species with the advantage of using mammalian lipid databases for data processing. This review highlights the use of zebrafish as an innovative model system to gain novel insights into the pathogenesis of sphingolipidoses, with possible implications for the identification of more efficacious therapeutic approaches.

Vorasidenib and ivosidenib in IDH1-mutant low-grade glioma: a randomized, perioperative phase 1 trial

Vorasidenib and ivosidenib inhibit mutant forms of isocitrate dehydrogenase (mIDH) and have shown preliminary clinical activity against mIDH glioma. We evaluated both agents in a perioperative phase 1 trial to explore the mechanism of action in recurrent low-grade glioma (IGG) and select a molecule for phase 3 testing. Primary end-point was concentration of D-2-hydroxyglutarate (2-HG), the metabolic product of mIDH enzymes, measured in tumor tissue from 49 patients with mIDH1-R132H nonenhancing gliomas following randomized treatment with vorasidenib (50 mg or 10 mg once daily, q.d.), ivosidenib (500 mg q.d. or 250 mg twice daily) or no treatment before surgery. Tumor 2-HG concentrations were reduced by 92.6% (95% credible interval (CrI), 76.1-97.6) and 91.1% (95% CrI, 72.0-97.0) in patients treated with vorasidenib 50 mg q.d. and ivosidenib 500 mg q.d., respectively. Both agents were well tolerated and follow-up is ongoing. In exploratory analyses, 2-HG reduction was associated with increased DNA 5-hydroxymethylcytosine, reversal of 'proneural' and 'stemness' gene expression signatures, decreased tumor cell proliferation and immune cell activation. Vorasidenib, which showed brain penetrance and more consistent 2-HG suppression than ivosidenib, was advanced to phase 3 testing in patients with mIDH LGGs. Funded by Agios Pharmaceuticals, Inc. and Servier Pharmaceuticals LLC; ClinicalTrials.gov number NCT03343197.

Profiling TREM2 expression in amyotrophic lateral sclerosis

BACKGROUND AND OBJECTIVES

There is growing evidence of the contribution of neuroinflammation, and in particular microglia, in the pathogenesis of amyotrophic lateral sclerosis (ALS). TREM2 gene plays a crucial role in shaping microglia in neurodegenerative conditions. To deepen the understanding of TREM2 in ALS and investigate the performance of TREM2 as a biomarker, we profiled TREM2 expression levels in spinal cord, cerebrospinal fluid and blood of patients with sporadic ALS. We also wanted to investigate whether the combined measurement of sTREM2 in fluids could improve the diagnostic yield of total and phosphorylated TDP-43 levels.

METHODS

We performed a case-control study to profile overall and transcript-specific TREM2 mRNA levels by RT-qPCR and protein expression levels by Western-blot in postmortem specimens of spinal cord from ALS patients and controls. In parallel, we measured soluble TREM2 (sTREM2) protein levels and full length and phosphorylated TDP-43 (tTDP-43 and pTDP-43) by ELISA in CSF and serum from ALS patients vs healthy controls. Patients were prospectively recruited from an ALS unit of a tertiary hospital and fulfilled El Escorial revised criteria. After bivariate analysis, a logistic regression model was developed to identify adjusted estimates of the association of sTREM2 levels in CSF and serum with ALS status.

RESULTS

Overall and transcript-specific TREM2 mRNA were upregulated in the spinal cord of ALS patients (n = 21) compared to controls (n = 19). Similar changes were observed in TREM2 protein levels (p < 0.01) in spinal cord of ALS patients vs healthy controls. We also detected significantly higher sTREM2 levels in CSF (p-value < 0.01) of ALS patients (n = 46) vs controls (n = 46) and serum (p-value < 0.001) of ALS patients (n = 100) vs controls (n = 100). In a logistic regression model, both CSF and serum sTREM2 remained independently associated with ALS status with OR = 3.41 (CI 95 %=1.34-8.66) (p-value < 0.05) and OR = 3.38 (CI 95 %: 1.86-6.16) (p-value < 0.001), respectively. We also observed that pTDP-43 levels in CSF is an independent predictor of ALS (p-value < 0.05).

CONCLUSIONS

Our results support the role of TREM2 in ALS pathophysiology and demonstrates that the three TREM2 transcripts are deregulated in ALS in postmortem human specimens of spinal cord. We hypothesise about the possible influence of systemic-peripheral inflammation in the disease. Finally, we conclude that pTDP-43 levels in CSF could be a biomarker of ALS, and sTREM2 measurement in CSF and blood emerge as potential non-invasive biomarker in ALS.

Gβγ subunits colocalize with RNA polymerase II and regulate transcription in cardiac fibroblasts

Gβγ subunits mediate many different signaling processes in various compartments of the cell, including the nucleus. To gain insight into the functions of nuclear Gβγ signaling, we investigated the functional role of Gβγ signaling in the regulation of GPCR-mediated gene expression in primary rat neonatal cardiac fibroblasts. We identified a novel, negative, regulatory role for the Gβ₁γ dimer in the fibrotic response. Depletion of Gβ₁ led to derepression of the fibrotic response at the mRNA and protein levels under basal conditions and an enhanced fibrotic response after sustained stimulation of the angiotensin II type I receptor. Our genome-wide chromatin immunoprecipitation experiments revealed that Gβ₁ colocalized and interacted with RNA polymerase II on fibrotic genes in an angiotensin II-dependent manner. Additionally, blocking transcription with inhibitors of Cdk9 prevented association of Gβγ with transcription complexes. Together, our findings suggest that Gβ₁γ is a novel transcriptional regulator of the fibrotic response that may act to restrict fibrosis to conditions of sustained fibrotic signaling. Our work expands the role for Gβγ signaling in cardiac fibrosis and may have broad implications for the role of nuclear Gβγ signaling in other cell types.

mtR_find: A Parallel Processing Tool to Identify and Annotate RNAs Derived from the Mitochondrial Genome

RNAs originating from mitochondrial genomes are abundant in transcriptomic datasets produced by high-throughput sequencing technologies, primarily in short-read outputs. Specific features of mitochondrial small RNAs (mt-sRNAs), such as non-templated additions, presence of length variants, sequence variants, and other modifications, necessitate the need for the development of an appropriate tool for their effective identification and annotation. We have developed mtR_find, a tool to detect and annotate mitochondrial RNAs, including mt-sRNAs and mitochondria-derived long non-coding RNAs (mt-lncRNA). mtR_find uses a novel method to compute the count of RNA sequences from adapter-trimmed reads. When analyzing the published datasets with mtR_find, we identified mt-sRNAs significantly associated with the health conditions, such as hepatocellular carcinoma and obesity, and we discovered novel mt-sRNAs. Furthermore, we identified mt-lncRNAs in early development in mice. These examples show the immediate impact of miR_find in extracting a novel biological information from the existing sequencing datasets. For benchmarking, the tool has been tested on a simulated dataset and the results were concordant. For accurate annotation of mitochondria-derived RNA, particularly mt-sRNA, we developed an appropriate nomenclature. mtR_find encompasses the mt-ncRNA transcriptomes in unpreceded resolution and simplicity, allowing re-analysis of the existing transcriptomic databases and the use of mt-ncRNAs as diagnostic or prognostic markers in the field of medicine.

Neural functions in cancer: Data analyses and database construction

Recent studies have revealed that neural functions are involved in possibly every aspect of a cancer development, serving as bridges connecting microenvironmental stressors, activities of intracellular subsystems, and cell survival. Elucidation of the functional roles played by the neural system could provide the missing links in developing a systems-level understanding of cancer biology. However, the existing information is highly fragmented and scattered across the literature and internet databases, making it difficult for cancer researchers to use. We have conducted computational analyses of transcriptomic data of cancer tissues in TCGA and tissues of healthy organs in GTEx, aiming to demonstrate how the functional roles by the neural genes could be derived and what non-neural functions they are associated with, across different stages of 26 cancer types. Several novel discoveries are made, including i) the expressions of certain neural genes can predict the prognosis of a cancer patient; ii) cancer metastasis tends to involve specific neural functions; iii) cancers of low survival rates involve more neural interactions than those with high survival rates; iv) more malignant cancers involve more complex neural functions; and v) neural functions are probably induced to alleviate stresses and help the associated cancer cells to survive. A database, called NGC, is developed for organizing such derived neural functions and associations, along with gene expressions and functional annotations collected from public databases, aiming to provide an integrated and publicly available information resource to enable cancer researchers to take full advantage of the relevant information in their research, facilitated by tools provided by NGC.

LncCat: An ORF attention model to identify LncRNA based on ensemble learning strategy and fused sequence information

Background

Long non-coding RNA (lncRNA) is one of the most essential forms of transcripts, playing crucial regulatory roles in the development of cancers and diseases without protein-coding ability. It was assumed that short ORFs (sORFs) in lncRNA were weak to translate proteins. However, recent research has shown that sORFs can encode peptides, which increases the difficulty to identify lncRNA. Therefore, identifying lncRNAs with sORFs facilitates finding novel regulatory factors.

Results

In this paper, we propose LncCat for identifying lncRNA based on category boosting (CatBoost) and ORF-attention features. LncCat combines five types of features to encode transcript sequences and employs CatBoost to build a prediction model. In addition, the visualization comparison reveals that the ORF-attention features between lncRNAs and protein-coding transcripts are significantly distinct. The comparison results show that LncCat outperforms competing methods on several benchmark datasets. For Matthew's Correlation Coefficient (MCC), LncCat achieves 0.9503, 0.9219, 0.8591, 0.8672, and 0.9047 on the human, mouse, zebrafish, wheat, and chicken datasets, with improvements ranging from 1.90% to 7.82%, 1.49-17.63%, 6.11-21.50%, 3.02-51.64% and 5.35-26.90%, respectively. Moreover, LncCat dramatically improves the MCC by at least 11.90%, 12.96% and 42.61% on sORF test datasets of human, mouse, and zebrafish, respectively.

Conclusions

Experiments indicate that LncCat performs better both on long ORF and sORF datasets, and ORF-attention features show positive effects on predicting lncRNA. In brief, LncCat is a reliable method for identifying lncRNA. Additionally, a user-friendly web server is developed for academics at http://cczubio.top/lnccat.

Data resources and computational methods for lncRNA-disease association prediction

Increasing interest has been attracted in deciphering the potential disease pathogenesis through lncRNA-disease association (LDA) prediction, regarding to the diverse functional roles of lncRNAs in genome regulation. Whilst, computational models and algorithms benefit systematic biology research, even facilitate the classical biological experimental procedures. In this review, we introduce representative diseases associated with lncRNAs, such as cancers, cardiovascular diseases, and neurological diseases. Current publicly available resources related to lncRNAs and diseases have also been included. Furthermore, all of the 64 computational methods for LDA prediction have been divided into 5 groups, including machine learning-based methods, network propagation-based methods, matrix factorization- and completion-based methods, deep learning-based methods, and graph neural network-based methods. The common evaluation methods and metrics in LDA prediction have also been discussed. Finally, the challenges and future trends in LDA prediction have been discussed. Recent advances in LDA prediction approaches have been summarized in the GitHub repository at https://github.com/sheng-n/lncRNA-disease-methods.

UXT attenuates the CGAS-STING1 signaling by targeting STING1 for autophagic degradation

STING1 (stimulator of interferon response cGAMP interactor 1), the pivotal adaptor protein of CGAS (cyclic GMP-AMP synthase)-STING1 signaling, is critical for type I IFN production of innate immunity. However, excessive or prolonged activation of STING1 is associated with autoinflammatory and autoimmune diseases. Thus, preventing STING1 from over-activation is important to maintain immune homeostasis. Here, we reported that UXT (ubiquitously expressed prefoldin like chaperone), a small chaperone-like protein, was essential to prevent the excessive activation of STING1-mediated type I IFN signaling through autophagic degradation of STING1 via SQSTM1 (sequestosome 1). Upon DNA mimics or cyclic GMP-AMP (cGAMP) stimulation, UXT specifically interacted with STING1 and promoted STING1 degradation through selective macroautophagy/autophagy. Moreover, UXT was required for more efficient autophagic degradation of STING1 by facilitating the interaction of SQSTM1 and STING1. The in vivo role of UXT in attenuating the CGAS-STING1 signaling was further confirmed in the mouse model of DNA-virus infection and the TMPD (2,6,10,14-tetramethylpentadecane)-induced murine lupus model. Intriguingly, the expression of UXT was consistently impaired and exhibited a remarkable inverse correlation with type I IFN signature in the leukocytes and PBMCs (peripheral blood mononuclear cells) of several large SLE (systemic lupus erythematosus) cohorts. Importantly, the replenishment of UXT effectively suppressed the production of IFNs and ISGs in the PBMCs of SLE patients. Taken together, our study reveals a novel regulatory role of UXT in autophagic degradation of STING1 to maintain immune homeostasis. UXT might be a potential therapeutic target for alleviating aberrant type I IFNs in autoimmune diseasesAbbreviations: 3-MA: 3-methyladenine; BMDMs: bone marrow-derived macrophages; cGAMP: cyclic GMP-AMP; CGAS: cyclic gmp-amp synthase; cKO: conditional knockout; CXCL10: C-X-C motif chemokine ligand 10; GAPDH: glyceraldehyde-3-phosphate dehydrogenase; HSV-1: herpes simplex virus type 1; HTDNA: herring testes DNA; IFIT1: interferon induced protein with tetratricopeptide repeats 1; IFNA4: interferon alpha 4; IFNB: interferon beta; IRF3: interferon regulatory factor 3; ISD: interferon stimulatory DNA; ISGs: IFN-stimulated genes; MAP1LC3B/LC3B: microtubule associated protein 1 light chain 3 beta; MEFs: mouse embryonic fibroblasts; RNA-seq: RNA sequencing; PBMCs: peripheral blood mononuclear cells; RSAD2: radical S-adenosyl methionine domain containing 2; SLE: systemic lupus erythematosus; SQSTM1: sequestosome 1; STING1: stimulator of interferon response cGAMP interactor 1; TBK1: TANK binding kinase 1; TMPD: 2,6,10,14-tetramethylpentadecane; UXT: ubiquitously expressed prefoldin like chaperone.

Transcriptome annotation of 17 porcine tissues using nanopore sequencing technology

The annotation of animal genomes plays an important role in elucidating molecular mechanisms behind the genetic control of economically important traits. Here, we employed long-read sequencing technology, Oxford Nanopore Technology, to annotate the pig transcriptome across 17 tissues from two Yorkshire littermate pigs. More than 9.8 million reads were obtained from a single flow cell, and 69 781 unique transcripts at 50 108 loci were identified. Of these transcripts, 16 255 were found to be novel isoforms, and 22 344 were found at loci that were novel and unannotated in the Ensembl (release 102) and NCBI (release 106) annotations. Novel transcripts were mostly expressed in cerebellum, followed by lung, liver, spleen, and hypothalamus. By comparing the unannotated transcripts to existing databases, there were 21 285 (95.3%) transcripts matched to the NT database (v5) and 13 676 (61.2%) matched to the NR database (v5). Moreover, there were 4324 (19.4%) transcripts matched to the SwissProt database (v5), corresponding to 11 356 proteins. Tissue-specific gene expression analyses showed that 9749 transcripts were highly tissue-specific, and cerebellum contained the most tissue-specific transcripts. As the same samples were used for the annotation of cis-regulatory elements in the pig genome, the transcriptome annotation generated by this study provides an additional and complementary annotation resource for the Functional Annotation of Animal Genomes effort to comprehensively annotate the pig genome.

miRinGO: Prediction of Biological Processes Indirectly Targeted by Human microRNAs

MicroRNAs (miRNAs) are small non-coding RNAs that are known for their role in the post-transcriptional regulation of target genes. Typically, their functions are predicted by first identifying their target genes and then finding biological processes enriched in these targets. Current tools for miRNA functional analysis use only genes with physical binding sites as their targets and exclude other genes that are indirectly targeted transcriptionally through transcription factors. Here, we introduce a method to predict gene ontology (GO) annotations indirectly targeted by microRNAs. The proposed method resulted in better performance in predicting known miRNA-GO term associations compared to the canonical approach. To facilitate miRNA GO enrichment analysis, we developed an R Shiny application, miRinGO, that is freely available online at GitHub.

IRF4 International Consortium. A multimorphic mutation in IRF4 causes human autosomal dominant combined immunodeficiency

Interferon regulatory factor 4 (IRF4) is a transcription factor (TF) and key regulator of immune cell development and function. We report a recurrent heterozygous mutation in IRF4, p.T95R, causing an autosomal dominant combined immunodeficiency (CID) in seven patients from six unrelated families. The patients exhibited profound susceptibility to opportunistic infections, notably Pneumocystis jirovecii, and presented with agammaglobulinemia. Patients' B cells showed impaired maturation, decreased immunoglobulin isotype switching, and defective plasma cell differentiation, whereas their T cells contained reduced TH17 and TFH populations and exhibited decreased cytokine production. A knock-in mouse model of heterozygous T95R showed a severe defect in antibody production both at the steady state and after immunization with different types of antigens, consistent with the CID observed in these patients. The IRF4T95R variant maps to the TF's DNA binding domain, alters its canonical DNA binding specificities, and results in a simultaneous multimorphic combination of loss, gain, and new functions for IRF4. IRF4T95R behaved as a gain-of-function hypermorph by binding to DNA with higher affinity than IRF4WT. Despite this increased affinity for DNA, the transcriptional activity on IRF4 canonical genes was reduced, showcasing a hypomorphic activity of IRF4T95R. Simultaneously, IRF4T95R functions as a neomorph by binding to noncanonical DNA sites to alter the gene expression profile, including the transcription of genes exclusively induced by IRF4T95R but not by IRF4WT. This previously undescribed multimorphic IRF4 pathophysiology disrupts normal lymphocyte biology, causing human disease.

The SPOC domain is a phosphoserine binding module that bridges transcription machinery with co- and post-transcriptional regulators

The heptad repeats of the C-terminal domain (CTD) of RNA polymerase II (Pol II) are extensively modified throughout the transcription cycle. The CTD coordinates RNA synthesis and processing by recruiting transcription regulators as well as RNA capping, splicing and 3'end processing factors. The SPOC domain of PHF3 was recently identified as a CTD reader domain specifically binding to phosphorylated serine-2 residues in adjacent CTD repeats. Here, we establish the SPOC domains of the human proteins DIDO, SHARP (also known as SPEN) and RBM15 as phosphoserine binding modules that can act as CTD readers but also recognize other phosphorylated binding partners. We report the crystal structure of SHARP SPOC in complex with CTD and identify the molecular determinants for its specific binding to phosphorylated serine-5. PHF3 and DIDO SPOC domains preferentially interact with the Pol II elongation complex, while RBM15 and SHARP SPOC domains engage with writers and readers of m6A, the most abundant RNA modification. RBM15 positively regulates m6A levels and mRNA stability in a SPOC-dependent manner, while SHARP SPOC is essential for its localization to inactive X-chromosomes. Our findings suggest that the SPOC domain is a major interface between the transcription machinery and regulators of transcription and co-transcriptional processes.

Cardiomyopathy prevalence exceeds 30% in individuals with TTN variants and early atrial fibrillation

PURPOSE

TTN truncating variants (TTNtvs) represent the largest known genetic cause of dilated cardiomyopathies (DCMs), however their penetrance for DCM in general populations is low. More broadly, patients with cardiomyopathies (CMs) often exhibit other cardiac conditions, such as atrial fibrillation (Afib), which has also been linked to TTNtvs. This retrospective analysis aims to characterize the relationship between different cardiac conditions in those with TTNtvs and identify individuals with the highest risk of DCM.

METHODS

In this work we leverage longitudinal electronic health record and exome sequencing data from approximately 450,000 individuals in 2 health systems to statistically confirm and pinpoint the genetic footprint of TTNtv-related diagnoses aside from CM, such as Afib, and determine whether vetting additional significantly associated phenotypes better stratifies CM risk across those with TTNtvs. We focused on TTNtvs in exons with a percentage spliced in >90% (hiPSI TTNtvs), a representation of constitutive cardiac expression.

RESULTS

When controlling for CM and Afib, other cardiac conditions retained only nominal association with TTNtvs. A sliding window analysis of TTNtvs across the locus confirms that the association is specific to hiPSI exons for both CM and Afib, with no meaningful associations in percent spliced in ≤90% exons (loPSI TTNtvs). The combination of hiPSI TTNtv status and early Afib diagnosis (before age 60) found a subset of TTNtv individuals at high risk for CM. The prevalence of CM in this subset was 33%, a rate that was 3.5 fold higher than that in individuals with hiPSI TTNtvs (9% prevalence), 5-fold higher than that in individuals without TTNtvs with early Afib (6% prevalence), and 80-fold higher than that in the general population.

CONCLUSION

Our retrospective analyses revealed that those with hiPSI TTNtvs and early Afib (∼1/2900) have a high prevalence of CM (33%), far exceeding that in other individuals with TTNtvs and in those without TTNtvs with an early Afib diagnosis. These results show that combining phenotypic information along with genomic population screening can identify patients at higher risk for progressing to symptomatic heart failure.

Predicting Modifiers of Genotype-Phenotype Correlations in Craniofacial Development

Most human birth defects are phenotypically variable even when they share a common genetic basis. Our understanding of the mechanisms of this variation is limited, but they are thought to be due to complex gene-environment interactions. Loss of the transcription factor Gata3 associates with the highly variable human birth defects HDR syndrome and microsomia, and can lead to disruption of the neural crest-derived facial skeleton. We have demonstrated that zebrafish gata3 mutants model the variability seen in humans, with genetic background and candidate pathways modifying the resulting phenotype. In this study, we sought to use an unbiased bioinformatic approach to identify environmental modifiers of gata3 mutant craniofacial phenotypes. The LINCs L1000 dataset identifies chemicals that generate differential gene expression that either positively or negatively correlates with an input gene list. These chemicals are predicted to worsen or lessen the mutant phenotype, respectively. We performed RNA-seq on neural crest cells isolated from zebrafish across control, Gata3 loss-of-function, and Gata3 rescue groups. Differential expression analyses revealed 551 potential targets of gata3. We queried the LINCs database with the 100 most upregulated and 100 most downregulated genes. We tested the top eight available chemicals predicted to worsen the mutant phenotype and the top eight predicted to lessen the phenotype. Of these, we found that vinblastine, a microtubule inhibitor, and clofibric acid, a PPAR-alpha agonist, did indeed worsen the gata3 phenotype. The Topoisomerase II and RNA-pol II inhibitors daunorubicin and triptolide, respectively, lessened the phenotype. GO analysis identified Wnt signaling and RNA polymerase function as being enriched in our RNA-seq data, consistent with the mechanism of action of some of the chemicals. Our study illustrates multiple potential pathways for Gata3 function, and demonstrates a systematic, unbiased process to identify modifiers of genotype-phenotype correlations.

Dysregulation of PRMT5 in chronic lymphocytic leukemia promotes progression with high risk of Richter's transformation

Richter's Transformation (RT) is a poorly understood and fatal progression of chronic lymphocytic leukemia (CLL) manifesting histologically as diffuse large B-cell lymphoma. Protein arginine methyltransferase 5 (PRMT5) is implicated in lymphomagenesis, but its role in CLL or RT progression is unknown. We demonstrate herein that tumors uniformly overexpress PRMT5 in patients with progression to RT. Furthermore, mice with B-specific overexpression of hPRMT5 develop a B-lymphoid expansion with increased risk of death, and Eµ-PRMT5/TCL1 double transgenic mice develop a highly aggressive disease with transformation that histologically resembles RT; where large-scale transcriptional profiling identifies oncogenic pathways mediating PRMT5-driven disease progression. Lastly, we report the development of a SAM-competitive PRMT5 inhibitor, PRT382, with exclusive selectivity and optimal in vitro and in vivo activity compared to available PRMT5 inhibitors. Taken together, the discovery that PRMT5 drives oncogenic pathways promoting RT provides a compelling rationale for clinical investigation of PRMT5 inhibitors such as PRT382 in aggressive CLL/RT cases.

Cell type-specific histone acetylation profiling of Alzheimer's disease subjects and integration with genetics

We profile genome-wide histone 3 lysine 27 acetylation (H3K27ac) of 3 major brain cell types from hippocampus and dorsolateral prefrontal cortex (dlPFC) of subjects with and without Alzheimer's Disease (AD). We confirm that single nucleotide polymorphisms (SNPs) associated with late onset AD (LOAD) show a strong tendency to reside in microglia-specific gene regulatory elements. Despite this significant colocalization, we find that microglia harbor more acetylation changes associated with age than with amyloid-β (Aβ) load. In contrast, we detect that an oligodendrocyte-enriched glial (OEG) population contains the majority of differentially acetylated peaks associated with Aβ load. These differential peaks reside near both early onset risk genes (APP, PSEN1, PSEN2) and late onset AD risk loci (including BIN1, PICALM, CLU, ADAM10, ADAMTS4, SORL1, FERMT2), Aβ processing genes (BACE1), as well as genes involved in myelinating and oligodendrocyte development processes. Interestingly, a number of LOAD risk loci associated with differentially acetylated risk genes contain H3K27ac peaks that are specifically enriched in OEG. These findings implicate oligodendrocyte gene regulation as a potential mechanism by which early onset and late onset risk genes mediate their effects, and highlight the deregulation of myelinating processes in AD. More broadly, our dataset serves as a resource for the study of functional effects of genetic variants and cell type specific gene regulation in AD.

Single-cell transcriptomes and T cell receptors of vaccine-expanded apolipoprotein B-specific T cells

Atherosclerotic cardiovascular diseases are the major cause of death worldwide. CD4 T cells responding to Apolipoprotein B (ApoB), the core protein of most lipoproteins, have been identified as critical disease modulators. In healthy individuals, ApoB-reactive (ApoB+) CD4 T cells are mostly regulatory T cells (Tregs), which exert anti-inflammatory effects. Yet, they may obtain pro-inflammatory features and thus become proatherogenic. Evidence from animal studies suggests that vaccination against certain major histocompatibility complex (MHC) II-binding ApoB peptides induces an expansion of ApoB+ Tregs and thus confers atheroprotection. To date, in-depth phenotyping of vaccine-expanded ApoB+ T cells has not yet been performed. To this end, we vaccinated C57BL/6J mice with the ApoB-peptide P6 (ApoB978-993 TGAYSNASSTESASY) and performed single-cell RNA sequencing of tetramer-sorted P6+ T cells. P6+ cells were clonally expanded (one major, two minor clones) and formed a transcriptional cluster distinct from clusters mainly containing non-expanded P6+ and P6- cells. Transcriptomic profiling revealed that most expanded P6+ cells had a strong Treg signature and highly expressed genes mediating suppressive functions. Yet, some expanded P6+ cells only had a residual Treg signature and expressed genes related to T helper 1 (TH1) cells, which are proatherogenic. Modeling the T cell receptor (TCR) and P6:MHC-II interaction showed that only three amino acid residues in the α and β chain contact the P6 peptide in the MHC-II groove and thus determine the specificity of this TCR to P6. Our data begin to reveal the vaccination-induced response to an ApoB epitope.

Bi-allelic variants in NAE1 cause intellectual disability, ischiopubic hypoplasia, stress-mediated lymphopenia and neurodegeneration

Neddylation has been implicated in various cellular pathways and in the pathophysiology of numerous diseases. We identified four individuals with bi-allelic variants in NAE1, which encodes the neddylation E1 enzyme. Pathogenicity was supported by decreased NAE1 abundance and overlapping clinical and cellular phenotypes. To delineate how cellular consequences of NAE1 deficiency would lead to the clinical phenotype, we focused primarily on the rarest phenotypic features, based on the assumption that these would best reflect the pathophysiology at stake. Two of the rarest features, neuronal loss and lymphopenia worsening during infections, suggest that NAE1 is required during cellular stress caused by infections to protect against cell death. In support, we found that stressing the proteasome system with MG132-requiring upregulation of neddylation to restore proteasomal function and proteasomal stress-led to increased cell death in fibroblasts of individuals with NAE1 genetic variants. Additionally, we found decreased lymphocyte counts after CD3/CD28 stimulation and decreased NF-κB translocation in individuals with NAE1 variants. The rarest phenotypic feature-delayed closure of the ischiopubic rami-correlated with significant downregulation of RUN2X and SOX9 expression in transcriptomic data of fibroblasts. Both genes are involved in the pathophysiology of ischiopubic hypoplasia. Thus, we show that NAE1 plays a major role in (skeletal) development and cellular homeostasis during stress. Our approach suggests that a focus on rare phenotypic features is able to provide significant pathophysiological insights in diseases caused by mutations in genes with pleiotropic effects.

Pathway-based integration of multi-omics data reveals lipidomics alterations validated in an Alzheimer's disease mouse model and risk loci carriers

Alzheimer's disease (AD) is a highly prevalent neurodegenerative disorder. Despite increasing evidence of the importance of metabolic dysregulation in AD, the underlying metabolic changes that may impact amyloid plaque formation are not understood, particularly for late-onset AD. This study analyzed genome-wide association studies (GWAS), transcriptomics, and proteomics data obtained from several data repositories to obtain differentially expressed (DE) multi-omics elements in mouse models of AD. We characterized the metabolic modulation in these data sets using gene ontology, transcription factor, pathway, and cell-type enrichment analyses. A predicted lipid signature was extracted from genome-scale metabolic networks (GSMN) and subsequently validated in a lipidomic data set derived from cortical tissue of ABCA-7 null mice, a mouse model of one of the genes associated with late-onset AD. Moreover, a metabolome-wide association study (MWAS) was performed to further characterize the association between dysregulated lipid metabolism in human blood serum and genes associated with AD risk. We found 203 DE transcripts, 164 DE proteins, and 58 DE GWAS-derived mouse orthologs associated with significantly enriched metabolic biological processes. Lipid and bioenergetic metabolic pathways were significantly over-represented across the AD multi-omics data sets. Microglia and astrocytes were significantly enriched in the lipid-predominant AD-metabolic transcriptome. We also extracted a predicted lipid signature that was validated and robustly modeled class separation in the ABCA7 mice cortical lipidome, with 11 of these lipid species exhibiting statistically significant modulations. MWAS revealed 298 AD single nucleotide polymorphisms-metabolite associations, of which 70% corresponded to lipid classes. These results support the importance of lipid metabolism dysregulation in AD and highlight the suitability of mapping AD multi-omics data into GSMNs to identify metabolic alterations.

Genome-wide identification of Arabidopsis non-AUG-initiated upstream ORFs with evolutionarily conserved regulatory sequences that control protein expression levels

This study identified four novel regulatory non-AUG-initiated upstream ORFs (uORFs) with evolutionarily conserved sequences in Arabidopsis and elucidated the mechanism by which a non-AUG-initiated uORF promotes main ORF translation. Upstream open reading frames (uORFs) are short ORFs found in the 5'-untranslated regions (5'-UTRs) of eukaryotic transcripts and can influence the translation of protein-coding main ORFs (mORFs). Recent genome-wide ribosome profiling studies have revealed that hundreds or thousands of uORFs initiate translation at non-AUG start codons. However, the physiological significance of these non-AUG uORFs has so far been demonstrated for only a few of them. In this study, to identify physiologically important regulatory non-AUG uORFs in Arabidopsis, we took an approach that combined bioinformatics and experimental analysis. Since physiologically important non-AUG uORFs are likely to be conserved across species, we first searched the Arabidopsis genome for non-AUG-initiated uORFs with evolutionarily conserved sequences. Then, we examined the effects of the conserved non-AUG uORFs on the expression of the downstream mORFs using transient expression assays. As a result, three inhibitory and one promotive non-AUG uORFs were identified. Among the inhibitory non-AUG uORFs, two exerted repressive effects on mORF expression in an amino acid sequence-dependent manner. These two non-AUG uORFs are likely to encode regulatory peptides that cause ribosome stalling, thereby enhancing their repressive effects. In contrast, one of the identified regulatory non-AUG uORFs promoted mORF expression by alleviating the inhibitory effect of a downstream AUG-initiated uORF. These findings provide insights into the mechanisms that enable non-AUG uORFs to play regulatory roles despite their low translation initiation efficiencies.

Novel homozygous leptin receptor mutation in an infant with monogenic obesity

Monogenic obesity can be caused by a mutation in one of the single genes involved in hunger and satiety. The most common mutations affect melanocortin 4 (MC4) followed by the leptin gene and its receptor. Leptin receptor (LEPR) gene mutation is an extremely rare endocrine disease characterized by early-onset obesity, hyperphagia in addition to pituitary hormone deficiency, and metabolic abnormalities. We report the case of a 12-month-old male infant born of a non-consanguineous marriage. He presented to us with rapid weight gain from 2 months of age along with hyperphagia. Biochemistry revealed a deranged lipid profile, elevated transaminases, and markedly raised serum leptin levels. On genetic analysis, a novel mutation was detected, which was a homozygous variation In exon 12 of the LEPR gene (chr1:g.65608901G>A) that resulted in the synonymous amino acid change of lysine at codon 584 proximal to donor splice site (p.Lys584). The in silico prediction of the variant was 'damaging' by MutationTaster2. The mutation was classified as a 'variant of uncertain significance' due to a lack of published literature and had to be correlated carefully with the clinical symptoms. It was recommended to do Sanger sequencing of the parents and other family members. However, due to financial constraints, the family could not afford the same. At the time of writing, funds were being arranged for procuring setmelanotide, which is a novel and effective therapy for monogenic obesity due to LepR mutation.