Genomic copy number variants: evidence for association with antibody response to anthrax vaccine adsorbed

Bioinformatics analysis of an immunotherapy responsiveness-related gene signature in predicting lung adenocarcinoma prognosis

Background

Immune therapy has become first-line treatment option for patients with lung cancer, but some patients respond poorly to immune therapy, especially among patients with lung adenocarcinoma (LUAD). Novel tools are needed to screen potential responders to immune therapy in LUAD patients, to better predict the prognosis and guide clinical decision-making. Although many efforts have been made to predict the responsiveness of LUAD patients, the results were limited. During the era of immunotherapy, this study attempts to construct a novel prognostic model for LUAD by utilizing differentially expressed genes (DEGs) among patients with differential immune therapy responses.

Methods

Transcriptome data of 598 patients with LUAD were downloaded from The Cancer Genome Atlas (TCGA) database, which included 539 tumor samples and 59 normal control samples, with a mean follow-up time of 29.69 months (63.1% of patients remained alive by the end of follow-up). Other data sources including three datasets from the Gene Expression Omnibus (GEO) database were analyzed, and the DEGs between immunotherapy responders and nonresponders were identified and screened. Univariate Cox regression analysis was applied with the TCGA cohort as the training set and GSE72094 cohort as the validation set, and least absolute shrinkage and selection operator (LASSO) Cox regression were applied in the prognostic-related genes which fulfilled the filter criteria to establish a prognostic formula, which was then tested with time-dependent receiver operating characteristic (ROC) analysis. Enriched pathways of the prognostic-related genes were analyzed with Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses, and tumor immune microenvironment (TIME), tumor mutational burden, and drug sensitivity tests were completed with appropriate packages in R (The R Foundation of Statistical Computing). Finally, a nomogram incorporating the prognostic formula was established.

Results

A total of 1,636 DEGs were identified, 1,163 prognostic-related DEGs were extracted, and 34 DEGs were selected and incorporated into the immunotherapy responsiveness-related risk score (IRRS) formula. The IRRS formula had good performance in predicting the overall prognoses in patients with LUAD and had excellent performance in prognosis prediction in all LUAD subgroups. Moreover, the IRRS formula could predict anticancer drug sensitivity and immunotherapy responsiveness in patients with LUAD. Mechanistically, immune microenvironments varied profoundly between the two IRRS groups; the most significantly varied pathway between the high-IRRS and low-IRRS groups was ribonucleoprotein complex biogenesis, which correlated closely with the TP53 and TTN mutation burdens. In addition, we established a nomogram incorporating the IRRS, age, sex, clinical stage, T-stage, N-stage, and M-stage as predictors that could predict the prognoses of 1-year, 3-year, and 5-year survival in patients with LUAD, with an area under curve (AUC) of 0.718, 0.702, and 0.68, respectively.

Conclusions

The model we established in the present study could predict the prognosis of LUAD patients, help to identify patients with good responses to anticancer drugs and immunotherapy, and serve as a valuable tool to guide clinical decision-making.

MAPT Locus in Parkinson's Disease Patients of Ashkenazi Origin: A Stratified Analysis

Introduction: MAPT locus is associated with Parkinson's disease (PD), which is located within a large inversion region of high linkage disequilibrium (LD). We aimed to determine whether the H2-haplotype protective effect and its effect size depends on the GBA1 or LRRK2 risk allele carrier status, and to further characterize genetic alterations that might contribute to its effect. Methods: LD analysis was performed using whole-genome sequencing data of 202 unrelated Ashkenazi Jewish (AJ) PDs. A haplotype-divergent variant was genotyped in a cohort of 1200 consecutively recruited AJ-PDs. The odd ratios were calculated using AJ-non-neuro cases from the gnomAD database as the controls in an un-stratified and a stratified manner according to the mutation carrier status, and the effect on the Age at Motor Symptom Onset (AMSO) was examined. Expression and splicing quantitative trait locus (eQTL and sQTL) analyses were carried out using brain tissues from a database. Results: The H2 haplotype exhibited significant association with PD protection, with a similar effect size in GBA1 carriers, LRRK2-G2019S carriers, and non-carriers (OR = 0.77, 0.69, and 0.82, respectively), and there was no effect on AMSO. The LD interval was narrowed to approximately 1.2 Mb. The H2 haplotype carried potential variants in candidate genes (MAPT and SPPL2C); structural deletions and segmental duplication (KANSL1); and variants affecting gene expression and intron excision ratio in brain tissues (LRRC37A/2). Conclusions: Our results demonstrate that H2 is associated with PD and its protective effect is not influenced by the GBA1/LRRK2 risk allele carrier status. This effect may be genetically complex, resulting from different levels of variations such as missense mutations in relevant genes, structural variations, epigenetic modifications, and RNA expression changes, which may operate independently or in synergy.

17q21.31 sub-haplotypes underlying H1-associated risk for Parkinson's disease are associated with LRRC37A/2 expression in astrocytes

BACKGROUND

Parkinson's disease (PD) is genetically associated with the H1 haplotype of the MAPT 17q.21.31 locus, although the causal gene and variants underlying this association have not been identified.

METHODS

To better understand the genetic contribution of this region to PD and to identify novel mechanisms conferring risk for the disease, we fine-mapped the 17q21.31 locus by constructing discrete haplotype blocks from genetic data. We used digital PCR to assess copy number variation associated with PD-associated blocks, and used human brain postmortem RNA-seq data to identify candidate genes that were then further investigated using in vitro models and human brain tissue.

RESULTS

We identified three novel H1 sub-haplotype blocks across the 17q21.31 locus associated with PD risk. Protective sub-haplotypes were associated with increased LRRC37A/2 copy number and expression in human brain tissue. We found that LRRC37A/2 is a membrane-associated protein that plays a role in cellular migration, chemotaxis and astroglial inflammation. In human substantia nigra, LRRC37A/2 was primarily expressed in astrocytes, interacted directly with soluble α-synuclein, and co-localized with Lewy bodies in PD brain tissue.

CONCLUSION

These data indicate that a novel candidate gene, LRRC37A/2, contributes to the association between the 17q21.31 locus and PD via its interaction with α-synuclein and its effects on astrocytic function and inflammatory response. These data are the first to associate the genetic association at the 17q21.31 locus with PD pathology, and highlight the importance of variation at the 17q21.31 locus in the regulation of multiple genes other than MAPT and KANSL1, as well as its relevance to non-neuronal cell types.

Epigenome-wide association study and network analysis for IgA Nephropathy from CD19+ B-cell in Chinese Population

IgA nephropathy (IgAN) is the most common primary glomerular disease in China and worldwide. The proliferation of B cells is known to be associated with both risk and prognosis of IgAN, but the epigenetic mechanism underlying this association is unknown. In this study we carried out the first Epigenome-wide Association Study (EWAS) by using the latest Infinium Methylation EPIC BeadChip on 184 B cell-specific samples (92 case/control pairs) for Chinese IgAN population. After rigorous data normalization and residual batch effect correction, linear mixed effect model was used to detect methylation CpG sites associated with IgAN adjusting for age, gender and smoking. False discovery rate (FDR) less than 10% was used to account for multiple testing. Weighted gene co-methylation networks were generated to identify gene modules highly correlated with IgAN. A permutation test was performed to account for the potential effect of overfitting. After adjusting clinical covariates and potential technical batch effects, three CpGs corresponding to PCDH17, TERT, WDR82 genes and three in the intergenic regions passed the genome-wide significant threshold. Methylation network analysis identified an additional IgAN associated gene module, containing 72 significant CpGs including GALNT6, IQSEC1, CDC16 and SYS1, involved in the pathway related to tubular atrophy/interstitial fibrosis of IgAN. These results suggested important DNA methylation and gene targets in CD19⁺ B cells for the pathogenesis of IgAN.

Recent developments in the understanding and use of anthrax vaccine adsorbed: achieving more with less

Anthrax Vaccine Adsorbed (AVA, BioThrax™) is the only Food and Drug Administration (FDA) approved vaccine for the prevention of anthrax in humans. Recent improvements in pre-exposure prophylaxis (PrEP) use of AVA include intramuscular (IM) administration and simplification of the priming series to three doses over 6 months. Administration IM markedly reduced the frequency, severity and duration of injection site reactions. Refinement of animal models for inhalation anthrax, identification of immune correlates of protection and cross-species modeling have created opportunities for reductions in the PrEP booster schedule and were pivotal in FDA approval of a post-exposure prophylaxis (PEP) indication. Clinical and nonclinical studies of accelerated PEP schedules and divided doses may provide prospects for shortening the PEP antimicrobial treatment period. These data may assist in determining feasibility of expanded coverage in a large-scale emergency when vaccine demand may exceed availability. Enhancements to the AVA formulation may broaden the vaccine's PEP application.

A Highly Polymorphic Copy Number Variant in the NSF Gene is Associated with Cocaine Dependence

Cocaine dependence is a complex psychiatric disorder involving both genetic and environmental factors. Several neurotransmitter systems mediate cocaine’s effects, dependence and relapse, being the components of the neurotransmitter release machinery good candidates for the disorder. Previously, we identified a risk haplotype for cocaine dependence in the NSF gene, encoding the protein N-Ethylmaleimide-Sensitive Factor essential for synaptic vesicle turnover. Here we examined the possible contribution to cocaine dependence of a large copy number variant (CNV) that encompasses part of the NSF gene. We performed a case-control association study in a discovery sample (359 cases and 356 controls) and identified an association between cocaine dependence and the CNV (P = 0.013), that was confirmed in the replication sample (508 cases and 569 controls, P = 7.1e-03) and in a pooled analysis (P = 1.8e-04), with an over-representation of low number of copies in cases. Subsequently, we studied the functional impact of the CNV on gene expression and found that the levels of two NSF transcripts were significantly increased in peripheral blood mononuclear cells (PBMC) along with the number of copies of the CNV. These results, together with a previous study from our group, support the role of NSF in the susceptibility to cocaine dependence.

The personal touch: strategies toward personalized vaccines and predicting immune responses to them

The impact of vaccines on public health and wellbeing has been profound. Smallpox has been eradicated, polio is nearing eradication, and multiple diseases have been eliminated from certain areas of the world. Unfortunately, we now face diseases such as hepatitis C, malaria or tuberculosis, as well as new and re-emerging pathogens for which we lack effective vaccines. Empirical approaches to vaccine development have been successful in the past, but may not be up to the current infectious disease challenges facing us. New, directed approaches to vaccine design, development, and testing need to be developed. Ideally these approaches will capitalize on cutting-edge technologies, advanced analytical and modeling strategies, and up-to-date knowledge of both pathogen and host. These approaches will pay particular attention to the causes of inter-individual variation in vaccine response in order to develop new vaccines tailored to the unique needs of individuals and communities within the population.