Functional genomics implicates natural killer cells as potential key drivers in the pathogenesis of ankylosing spondylitis.

Objective: Multiple lines of evidence indicate that ankylosing spondylitis (AS) is a lymphocyte-driven disease. However, which lymphocyte populations are critical in AS pathogenesis is not known. In this study, we aimed to identify the key cell types mediating the genetic risk in AS using an unbiased integrative functional genomics approach. Methods: We integrated GWAS data with epigenomic and transcriptomic datasets of immune cells in healthy humans. To quantify enrichment of cell type-specific open chromatin regions or gene expression in AS risk loci, we used three published methods which have identified cell types for other diseases. Additionally, we performed co-localization analyses between GWAS risk loci and genetic variants associated with gene expression (eQTL) to find putative target genes of AS risk variants. Results: Natural killer (NK) cell-specific open chromatin regions are significantly enriched in heritability for AS, compared to other immune cell types such as T cells, B cells, and monocytes. This finding was consistent between two AS GWAS. Using RNA-seq data, we validated that genes in AS risk loci are enriched in NK cell-specific gene expression. Expression levels of AS-associated genes, such as RUNX3, TBX21, TNFRSF1A, and NPEPPS, were found to be highest in NK cells compared to five T cell subsets. Using the human Space-Time Gut Cell Atlas we found significant upregulation of AS-associated genes predominantly in NK cells. Co-localization analysis revealed four AS risk loci affecting regulation of candidate target genes in NK cells: two known loci, ERAP1 and TNFRSF1A, and two under-studied loci, ENTR1 (aka SDCCAG3) and B3GNT2. Conclusion: Our results point to NK cells as potential key drivers in the development of AS and highlight four putative target genes for functional follow-up in NK cells.

History of tuberculosis disease is associated with genetic regulatory variation in Peruvians.

A quarter of humanity is estimated to be latently infected with Mycobacterium tuberculosis (Mtb) with a 5-10% risk of developing tuberculosis (TB) disease. Variability in responses to Mtb infection could be due to host or pathogen heterogeneity. Here, we focused on host genetic variation in a Peruvian population and its associations with gene regulation in monocyte-derived macrophages and dendritic cells (DCs). We recruited former household contacts of TB patients who previously progressed to TB (cases, n=63) or did not progress to TB (controls, n=63). Transcriptomic profiling of monocyte-derived dendritic cells (DCs) and macrophages measured the impact of genetic variants on gene expression by identifying expression quantitative trait loci (eQTL). We identified 330 and 257 eQTL genes in DCs and macrophages (False Discovery Rate (FDR) < 0.05), respectively. Five genes in DCs showed interaction between eQTL variants and TB progression status. The top eQTL interaction for a protein-coding gene was with FAH, the gene encoding fumarylacetoacetate hydrolase, which mediates the last step in mammalian tyrosine catabolism. FAH expression was associated with genetic regulatory variation in cases but not controls. Using public transcriptomic and epigenomic data of Mtb-infected monocyte-derived dendritic cells, we found that Mtb infection results in FAH downregulation and DNA methylation changes in the locus. Overall, this study demonstrates effects of genetic variation on gene expression levels that are dependent on history of infectious disease and highlights a candidate pathogenic mechanism through pathogen-response genes. Furthermore, our results point to tyrosine metabolism and related candidate TB progression pathways for further investigation.

A hotspot for posttranslational modifications on the androgen receptor dimer interface drives pathology and anti-androgen resistance.

Mutations of the androgen receptor (AR) associated with prostate cancer and androgen insensitivity syndrome may profoundly influence its structure, protein interaction network, and binding to chromatin, resulting in altered transcription signatures and drug responses. Current structural information fails to explain the effect of pathological mutations on AR structure-function relationship. Here, we have thoroughly studied the effects of selected mutations that span the complete dimer interface of AR ligand-binding domain (AR-LBD) using x-ray crystallography in combination with in vitro, in silico, and cell-based assays. We show that these variants alter AR-dependent transcription and responses to anti-androgens by inducing a previously undescribed allosteric switch in the AR-LBD that increases exposure of a major methylation target, Arg761. We also corroborate the relevance of residues Arg761 and Tyr764 for AR dimerization and function. Together, our results reveal allosteric coupling of AR dimerization and posttranslational modifications as a disease mechanism with implications for precision medicine.

SARS-CoV-2-specific T cell responses in patients with multisystem inflammatory syndrome in children.

Multisystem inflammatory syndrome in children (MIS-C) is a severe complication of SARS-CoV-2 infections that occurs in the pediatric population. We sought to characterize T cell responses in MIS-C compared to COVID-19 and pediatric hyperinflammatory syndromes. MIS-C was distinct from COVID-19 and hyperinflammatory syndromes due to an expansion of T cells expressing TRBV11-2 that was not associated with HLA genotype. Children diagnosed with MIS-C, but who were negative for SARS-CoV-2 by PCR and serology, did not display Vß skewing. There was no difference in the proportion of T cells that became activated after stimulation with SARS-CoV-2 peptides in children with MIS-C compared to convalescent COVID-19. The frequency of SARS-CoV-2-specific TCRs and the antigens recognized by these TCRs were comparable in MIS-C and COVID-19. Expansion of Vß11-2+ T cells was a specific biomarker of MIS-C patients with laboratory confirmed SARS-CoV-2 infections. Children with MIS-C had robust antigen-specific T cell responses to SARS-CoV-2.

Aedes aegypti Strain Subjected to Long-Term Exposure to Bacillus thuringiensis svar. israelensis Larvicides Displays an Altered Transcriptional Response to Zika Virus Infection.

Bacillus thuringiensis svar. israelensis (Bti) larvicides are effective in controlling Aedes aegypti; however, the effects of long-term exposure need to be properly evaluated. We established an Ae. aegypti strain that has been treated with Bti for 30 generations (RecBti) and is still susceptible to Bti, but females exhibited increased susceptibility to Zika virus (ZIKV). This study compared the RecBti strain to a reference strain regarding: first, the relative transcription of selected immune genes in ZIKV-challenged females (F30) with increased susceptibility detected in a previous study; then, the whole transcriptomic profile using unchallenged females (F35). Among the genes compared by RT-qPCR in the ZIKV-infected and uninfected females from RecBti (F30) and the reference strain, hop, domeless, relish 1, defensin A, cecropin D, and gambicin showed a trend of repression in RecBti infected females. The transcriptome of RecBti (F35) unchallenged females, compared with a reference strain by RNA-seq, showed a similar profile and only 59 differentially expressed genes were found among 9202 genes analyzed. Our dataset showed that the long-term Bti exposure of the RecBti strain was associated with an alteration of the expression of genes potentially involved in the response to ZIKV infection in challenged females, which is an important feature found under this condition.

Whole Genome Sequencing Analysis of Bacillus thuringiensis GR007 Reveals Multiple Pesticidal Protein Genes.

Bacillus thuringiensis (Bt) are soil ubiquitous bacteria. They produce a great variability of insecticidal proteins, where certain of these toxins are used worldwide for pest control. Through their adaptation to diverse ecosystems, certain Bt strains have acquired genetic mobile elements by horizontal transfer, harboring genes that encode for different virulent factors and pesticidal proteins (PP). Genomic characterization of Bt strains provides a valuable source of PP with potential biotechnological applications for pest control. In this work, we have sequenced the complete genome of the bacterium Bt GR007 strain that is toxic to Spodoptera frugiperda and Manduca sexta larvae. Four replicons (one circular chromosome and three megaplasmids) were identified. The two largest megaplasmids (pGR340 and pGR157) contain multiple genes that codify for pesticidal proteins: 10 cry genes (cry1Ab, cry1Bb, cry1Da, cry1Fb, cry1Hb, cry1Id, cry1Ja, cry1Ka, cry1Nb, and cry2Ad), two vip genes (vip3Af and vip3Ag), two binary toxin genes (vpa2Ac and vpb1Ca), five genes that codify for insecticidal toxin components (Tc's), and a truncated cry1Bd-like gene. In addition, genes that codify for several virulent factors were also found in this strain. Proteomic analysis of the parasporal crystals of GR007 revealed that they are composed of eight Cry proteins. Further cloning of these genes for their individual expression in Bt acrystalliferous strain, by means of their own intrinsic promoter showed expression of seven Cry proteins. These proteins display differential toxicity against M. sexta and S. frugiperda larvae, where Cry1Bb showed to be the most active protein against S. frugiperda larvae and Cry1Ka the most active protein against M. sexta larvae.