Comparison between qPCR and RNA-seq reveals challenges of quantifying HLA expression.

Human leukocyte antigen (HLA) class I and II loci are essential elements of innate and acquired immunity. Their functions include antigen presentation to T cells leading to cellular and humoral immune responses, and modulation of NK cells. Their exceptional influence on disease outcome has now been made clear by genome-wide association studies. The exons encoding the peptide-binding groove have been the main focus for determining HLA effects on disease susceptibility/pathogenesis. However, HLA expression levels have also been implicated in disease outcome, adding another dimension to the extreme diversity of HLA that impacts variability in immune responses across individuals. To estimate HLA expression, immunogenetic studies traditionally rely on quantitative PCR (qPCR). Adoption of alternative high-throughput technologies such as RNA-seq has been hampered by technical issues due to the extreme polymorphism at HLA genes. Recently, however, multiple bioinformatic methods have been developed to accurately estimate HLA expression from RNA-seq data. This opens an exciting opportunity to quantify HLA expression in large datasets but also brings questions on whether RNA-seq results are comparable to those by qPCR. In this study, we analyze three classes of expression data for HLA class I genes for a matched set of individuals: (a) RNA-seq, (b) qPCR, and (c) cell surface HLA-C expression. We observed a moderate correlation between expression estimates from qPCR and RNA-seq for HLA-A, -B, and -C (0.2 ≤ rho ≤ 0.53). We discuss technical and biological factors which need to be accounted for when comparing quantifications for different molecular phenotypes or using different techniques.

Expression estimation and eQTL mapping for HLA genes with a personalized pipeline.

The HLA (Human Leukocyte Antigens) genes are well-documented targets of balancing selection, and variation at these loci is associated with many disease phenotypes. Variation in expression levels also influences disease susceptibility and resistance, but little information exists about the regulation and population-level patterns of expression. This results from the difficulty in mapping short reads originated from these highly polymorphic loci, and in accounting for the existence of several paralogues. We developed a computational pipeline to accurately estimate expression for HLA genes based on RNA-seq, improving both locus-level and allele-level estimates. First, reads are aligned to all known HLA sequences in order to infer HLA genotypes, then quantification of expression is carried out using a personalized index. We use simulations to show that expression estimates obtained in this way are not biased due to divergence from the reference genome. We applied our pipeline to the GEUVADIS dataset, and compared the quantifications to those obtained with reference transcriptome. Although the personalized pipeline recovers more reads, we found that using the reference transcriptome produces estimates similar to the personalized pipeline (r ≥ 0.87) with the exception of HLA-DQA1. We describe the impact of the HLA-personalized approach on downstream analyses for nine classical HLA loci (HLA-A, HLA-C, HLA-B, HLA-DRA, HLA-DRB1, HLA-DQA1, HLA-DQB1, HLA-DPA1, HLA-DPB1). Although the influence of the HLA-personalized approach is modest for eQTL mapping, the p-values and the causality of the eQTLs obtained are better than when the reference transcriptome is used. We investigate how the eQTLs we identified explain variation in expression among lineages of HLA alleles. Finally, we discuss possible causes underlying differences between expression estimates obtained using RNA-seq, antibody-based approaches and qPCR.

Assessing false paternity risk in simulated motherless cases from more than 20 000 real exclusion trios.

BACKGROUND: When the mother's DNA profile is not available for paternity testing, there is a smaller probability that a locus will exclude an alleged father. This study aims to evaluate the risk of potential false paternity inclusions in motherless cases. STUDY DESIGN AND METHODS: More than 20 000 duos were generated by removing the maternal genotypes from exclusion trios. After recalculating paternity in these duos, any found inclusions would be false. RESULTS: The use of an appropriate number of loci, mutation model, and mutation rates to analyze motherless paternity cases was robust against false inclusions. A single potential false inclusion was observed in a case wherein kinship plays a role. This result highlights the importance of testing the mother when available and of obtaining information on family circumstances for the proper handling of cases involving related individuals. CONCLUSION: The guidelines we used here were sufficient to avoid false inclusions in a data set of more than 20 000 motherless cases.

An immunogenetic view of COVID-19.

Meeting the challenges brought by the COVID-19 pandemic requires an interdisciplinary approach. In this context, integrating knowledge of immune function with an understanding of how genetic variation influences the nature of immunity is a key challenge. Immunogenetics can help explain the heterogeneity of susceptibility and protection to the viral infection and disease progression. Here, we review the knowledge developed so far, discussing fundamental genes for triggering the innate and adaptive immune responses associated with a viral infection, especially with the SARS-CoV-2 mechanisms. We emphasize the role of the HLA and KIR genes, discussing what has been uncovered about their role in COVID-19 and addressing methodological challenges of studying these genes. Finally, we comment on questions that arise when studying admixed populations, highlighting the case of Brazil. We argue that the interplay between immunology and an understanding of genetic associations can provide an important contribution to our knowledge of COVID-19.

A genomic perspective on HLA evolution.

Several decades of research have convincingly shown that classical human leukocyte antigen (HLA) loci bear signatures of natural selection. Despite this conclusion, many questions remain regarding the type of selective regime acting on these loci, the time frame at which selection acts, and the functional connections between genetic variability and natural selection. In this review, we argue that genomic datasets, in particular those generated by next-generation sequencing (NGS) at the population scale, are transforming our understanding of HLA evolution. We show that genomewide data can be used to perform robust and powerful tests for selection, capable of identifying both positive and balancing selection at HLA genes. Importantly, these tests have shown that natural selection can be identified at both recent and ancient timescales. We discuss how findings from genomewide association studies impact the evolutionary study of HLA genes, and how genomic data can be used to survey adaptive change involving interaction at multiple loci. We discuss the methodological developments which are necessary to correctly interpret genomic analyses involving the HLA region. These developments include adapting the NGS analysis framework so as to deal with the highly polymorphic HLA data, as well as developing tools and theory to search for signatures of selection, quantify differentiation, and measure admixture within the HLA region. Finally, we show that high throughput analysis of molecular phenotypes for HLA genes-namely transcription levels-is now a feasible approach and can add another dimension to the study of genetic variation.

Epidemiological aspects of syringomyelia in a 19-year old cohort of spinal cord injury patients.

OBJECTIVE: To determine the prevalence, clinical and radiological risk factors, and surgical management of post-traumatic syringomyelia (PTS) in a 19-year cohort study of Spinal Cord Injury (SCI) patients treated at a SCI rehabilitation center. METHODS: Retrospective study of SCI patients in whom PTS was radiologically confirmed between January 2000 and December 2018. Protocols for assessing signs and symptoms of PTS were applied prior to PTS diagnosis and treatment and later at neurosurgical and rehabilitation reviews. The variables analyzed were prevalence, demographic data, trauma event, clinical and radiological risk factors, location and size of the syrinx, and effectiveness of the surgical procedures. RESULTS: Over the 19-year period, review of 920 SCI patients revealed 85 patients who met the clinical and neuroradiological criteria for the diagnosis of PTS and who were prospectively followed. Road traffic accidents were the leading cause of injury (n = 58; 68.2%), syringomyelia was most commonly observed in the thoracic spine (n = 56; 65.9%), and upper extremity paresis was the most common indication for surgical treatment (n = 27; 45%). Surgical treatment was indicated in 48 patients and the operative procedures included 29 syringopleural shunts (60.4%), 17 adhesiolysis (35.4%), and two syringosubarachnoid shunts (4.1%). The prevalence of PTS was 9% and was higher in patients with ASIA impairment scale grade A injuries. Most patients with PTS (63/85, 74.1%) were treated surgically at the time of injury. There was a significant reduction both in the extent (p = 0.05) and largest area (p = 0.001) of the syrinx after surgical treatment. Reoperation rates were 47% and 37.9% for adhesiolysis and syringopleural shunting, respectively. CONCLUSION: Follow-up and routine clinical examination of SCI patients is critical for the diagnosis of PTS in patients with late neurological deterioration. Surgical treatment has a positive impact in reducing the size of the syrinx as seen on postoperative MRI.

High-throughput identification of functional regulatory SNPs in systemic lupus erythematosus.

Genome-wide association studies implicate multiple loci in risk for systemic lupus erythematosus (SLE), but few contain exonic variants, rendering systematic identification of non-coding variants essential to decoding SLE genetics. We utilized SNP-seq and bioinformatic enrichment to interrogate 2180 single-nucleotide polymorphisms (SNPs) from 87 SLE risk loci for potential binding of transcription factors and related proteins from B cells. 52 SNPs that passed initial screening were tested by electrophoretic mobility shift and luciferase reporter assays. To validate the approach, we studied rs2297550 in detail, finding that the risk allele enhanced binding to the transcription factor Ikaros (IKZF1), thereby modulating expression of IKBKE. Correspondingly, primary cells from genotyped healthy donors bearing the risk allele expressed higher levels of the interferon / NF-κB regulator IKKϵ. Together, these findings define a set of likely functional non-coding lupus risk variants and identify a new regulatory pathway involving rs2297550, Ikaros, and IKKϵ implicated by human genetics in risk for SLE.

High-throughput identification of functional regulatory SNPs in systemic lupus erythematosus.

Genome-wide association studies implicate multiple loci in risk for systemic lupus erythematosus (SLE), but few contain exonic variants, rendering systematic identification of non-coding variants essential to decoding SLE genetics. We utilized SNP-seq and bioinformatic enrichment to interrogate 2180 single-nucleotide polymorphisms (SNPs) from 87 SLE risk loci for potential binding of transcription factors and related proteins from B cells. 52 SNPs that passed initial screening were tested by electrophoretic mobility shift and luciferase reporter assays. To validate the approach, we studied rs2297550 in detail, finding that the risk allele enhanced binding to the transcription factor Ikaros (encoded by IKZF1), thereby modulating expression of IKBKE. Correspondingly, primary cells from genotyped healthy donors bearing the risk allele expressed higher levels of the interferon / NF-κB regulator IKKε. Together, these findings define a set of likely functional non-coding lupus risk variants and identify a regulatory pathway involving rs2297550, Ikaros, and IKKε implicated by human genetics in risk for SLE.

Rhinovirus infection of airway epithelial cells uncovers the non-ciliated subset as a likely driver of genetic susceptibility to childhood-onset asthma.

Asthma is a complex disease caused by genetic and environmental factors. Epidemiological studies have shown that in children, wheezing during rhinovirus infection (a cause of the common cold) is associated with asthma development during childhood. This has led scientists to hypothesize there could be a causal relationship between rhinovirus infection and asthma or that RV-induced wheezing identifies individuals at increased risk for asthma development. However, not all children who wheeze when they have a cold develop asthma. Genome-wide association studies (GWAS) have identified hundreds of genetic variants contributing to asthma susceptibility, with the vast majority of likely causal variants being non-coding. Integrative analyses with transcriptomic and epigenomic datasets have indicated that T cells drive asthma risk, which has been supported by mouse studies. However, the datasets ascertained in these integrative analyses lack airway epithelial cells. Furthermore, large-scale transcriptomic T cell studies have not identified the regulatory effects of most non-coding risk variants in asthma GWAS, indicating there could be additional cell types harboring these "missing regulatory effects". Given that airway epithelial cells are the first line of defense against rhinovirus, we hypothesized they could be mediators of genetic susceptibility to asthma. Here we integrate GWAS data with transcriptomic datasets of airway epithelial cells subject to stimuli that could induce activation states relevant to asthma. We demonstrate that epithelial cultures infected with rhinovirus significantly upregulate childhood-onset asthma-associated genes. We show that this upregulation occurs specifically in non-ciliated epithelial cells. This enrichment for genes in asthma risk loci, or 'asthma heritability enrichment' is also significant for epithelial genes upregulated with influenza infection, but not with SARS-CoV-2 infection or cytokine activation. Additionally, cells from patients with asthma showed a stronger heritability enrichment compared to cells from healthy individuals. Overall, our results suggest that rhinovirus infection is an environmental factor that interacts with genetic risk factors through non-ciliated airway epithelial cells to drive childhood-onset asthma.

Rhinovirus infection of airway epithelial cells uncovers the non-ciliated subset as a likely driver of genetic risk to childhood-onset asthma.

Asthma is a complex disease caused by genetic and environmental factors. Studies show that wheezing during rhinovirus infection correlates with childhood asthma development. Over 150 non-coding risk variants for asthma have been identified, many affecting gene regulation in T cells, but the effects of most risk variants remain unknown. We hypothesized that airway epithelial cells could also mediate genetic susceptibility to asthma given they are the first line of defense against respiratory viruses and allergens. We integrated genetic data with transcriptomics of airway epithelial cells subject to different stimuli. We demonstrate that rhinovirus infection significantly upregulates childhood-onset asthma-associated genes, particularly in non-ciliated cells. This enrichment is also observed with influenza infection but not with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) or cytokine activation. Overall, our results suggest that rhinovirus infection is an environmental factor that interacts with genetic risk factors through non-ciliated airway epithelial cells to drive childhood-onset asthma.

Non-coding autoimmune risk variant defines role for ICOS in T peripheral helper cell development.

Fine-mapping and functional studies implicate rs117701653, a non-coding single nucleotide polymorphism in the CD28/CTLA4/ICOS locus, as a risk variant for rheumatoid arthritis and type 1 diabetes. Here, using DNA pulldown, mass spectrometry, genome editing and eQTL analysis, we establish that the disease-associated risk allele is functional, reducing affinity for the inhibitory chromosomal regulator SMCHD1 to enhance expression of inducible T-cell costimulator (ICOS) in memory CD4+ T cells from healthy donors. Higher ICOS expression is paralleled by an increase in circulating T peripheral helper (Tph) cells and, in rheumatoid arthritis patients, of blood and joint fluid Tph cells as well as circulating plasmablasts. Correspondingly, ICOS ligation and carriage of the rs117701653 risk allele accelerate T cell differentiation into CXCR5-PD-1high Tph cells producing IL-21 and CXCL13. Thus, mechanistic dissection of a functional non-coding variant in human autoimmunity discloses a previously undefined pathway through which ICOS regulates Tph development and abundance.

FIKElectricity: A Electricity Consumption Dataset from Three Restaurant Kitchens in Portugal.

Industrial Kitchens are very energy-intensive businesses, consuming between five and seven times more energy per square meter than other commercial spaces like office spaces and retail stores. Still, very little research has been carried out on improving the energy efficiency of this industry. This paper presents the FIKElectricity dataset, a collection of electricity data from three Portuguese restaurant kitchens during their daily operation. The duration of the datasets spans three to four weeks in each industrial kitchen and comprises aggregated and individual appliance consumption, collected at 1 Hz and [Formula: see text] Hz, respectively. The public release of FIKElectricity is expected to draw more attention from the research community to this overlooked industrial sector. The data collection and post-processing methods are thoroughly described in this paper, as well as the dataset organization. Examples showing the dataset's quality and instructions for its reuse are also provided.

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.

Carbapenem-Resistant Pseudomonas aeruginosa Spondylodiscitis Treated with Ceftazidime-Avibactam: A Case Report with Literature Review.

Pyogenic spondylodiscitis (PS) is a highly morbid and potentially fatal bacterial infection with an increasing incidence in recent decades. Its diagnosis and treatment are challenging, especially with the expansion of multidrug- or extensively drug-resistant bacteria. We report a rare case of PS caused by carbapenem-resistant Pseudomonas aeruginosa (CRPA) that was treated with ceftazidime-avibactam (C/A). The choice of C/A therapy was based on the patient's bacterial sensitivity profile and intolerance to the initial therapeutic regimen (polymyxin B and meropenem). The total antimicrobial treatment time was seven weeks. The evolution of the clinical course met the cure criteria, which was characterized by remission of signs and symptoms, normalization of inflammatory markers, and radiological improvement over 18 months of clinical follow-up. This is a rare case of CRPA spondylodiscitis that responded to C/A treatment.

Identification of a regulatory pathway governing TRAF1 via an arthritis-associated non-coding variant.

TRAF1/C5 was among the first loci shown to confer risk for inflammatory arthritis in the absence of an associated coding variant, but its genetic mechanism remains undefined. Using Immunochip data from 3,939 patients with juvenile idiopathic arthritis (JIA) and 14,412 control individuals, we identified 132 plausible common non-coding variants, reduced serially by single-nucleotide polymorphism sequencing (SNP-seq), electrophoretic mobility shift, and luciferase studies to the single variant rs7034653 in the third intron of TRAF1. Genetically manipulated experimental cells and primary monocytes from genotyped donors establish that the risk G allele reduces binding of Fos-related antigen 2 (FRA2), encoded by FOSL2, resulting in reduced TRAF1 expression and enhanced tumor necrosis factor (TNF) production. Conditioning on this JIA variant eliminated attributable risk for rheumatoid arthritis, implicating a mechanism shared across the arthritis spectrum. These findings reveal that rs7034653, FRA2, and TRAF1 mediate a pathway through which a non-coding functional variant drives risk of inflammatory arthritis in children and adults.

Mapping the dynamic genetic regulatory architecture of HLA genes at single-cell resolution.

The human leukocyte antigen (HLA) locus plays a critical role in complex traits spanning autoimmune and infectious diseases, transplantation and cancer. While coding variation in HLA genes has been extensively documented, regulatory genetic variation modulating HLA expression levels has not been comprehensively investigated. Here we mapped expression quantitative trait loci (eQTLs) for classical HLA genes across 1,073 individuals and 1,131,414 single cells from three tissues. To mitigate technical confounding, we developed scHLApers, a pipeline to accurately quantify single-cell HLA expression using personalized reference genomes. We identified cell-type-specific cis-eQTLs for every classical HLA gene. Modeling eQTLs at single-cell resolution revealed that many eQTL effects are dynamic across cell states even within a cell type. HLA-DQ genes exhibit particularly cell-state-dependent effects within myeloid, B and T cells. For example, a T cell HLA-DQA1 eQTL ( rs3104371 ) is strongest in cytotoxic cells. Dynamic HLA regulation may underlie important interindividual variability in immune responses.

Mapping the dynamic genetic regulatory architecture of HLA genes at single-cell resolution.

The human leukocyte antigen (HLA) locus plays a critical role in complex traits spanning autoimmune and infectious diseases, transplantation, and cancer. While coding variation in HLA genes has been extensively documented, regulatory genetic variation modulating HLA expression levels has not been comprehensively investigated. Here, we mapped expression quantitative trait loci (eQTLs) for classical HLA genes across 1,073 individuals and 1,131,414 single cells from three tissues, using personalized reference genomes to mitigate technical confounding. We identified cell-type-specific cis-eQTLs for every classical HLA gene. Modeling eQTLs at single-cell resolution revealed that many eQTL effects are dynamic across cell states even within a cell type. HLA-DQ genes exhibit particularly cell-state-dependent effects within myeloid, B, and T cells. Dynamic HLA regulation may underlie important interindividual variability in immune responses.

Comparison of DGGE and immunohistochemistry in the detection of TP53 variants in a Brazilian sample of sporadic breast tumors.

The presence of TP53 gene mutations in breast cancer has been associated with worse prognosis. These mutations interfere with the ability of the p53 protein, a transcription factor, to regulate the expression of target genes. Unlike the wild-type protein, which is rapidly degraded in cells, mutated forms have increased half-life and accumulate in tumor cells. Immunohistochemistry (IHC) is widely used in Brazil in the determination of breast cancer patients' prognosis. However, this technique is not able to detect many altered forms of the p53 protein (false-negative results) and readily detects the accumulation of wild-type p53 (false-positive results) that is associated with non-tumoral processes. For these reasons, we have set out to compare the efficiency of IHC with a molecular technique that detects gene variations at the DNA level in the evaluation of Brazilian patients with sporadic breast cancer. We have used denaturing gradient gel electrophoresis (DGGE) to study the TP53 status in 45 tumors, finding 26 allelic variants, most of them located in exon 4. Comparing the two techniques, IHC showed a false-negative rate of 64% and a false-positive rate of 50%. These results confirm the inability of IHC to correctly detect TP53 status, reason because it should not be routinely used to establish prognosis of breast cancer patients in Brazilian Pathology Laboratories. We recommend the utilization of a screening method, such as DGGE, followed by sequencing of altered exonic fragments to correctly detect TP53 gene variants and establish the prognosis of breast cancer patients.

Genetic diversity and statistical parameters of 15 autosomal STR loci in the Pomeranian subpopulation of Espirito Santo State, Brazil.

Allelic frequencies and other population data analysis are reported for the 15 autosomal Short Tandem Repeats (STR) loci included in the PowerPlex(®)16 kit (CSF1PO, D13S317, D16S539, D18S51, D21S11, D3S1358, D5S818, D7S820, D8S1179, FGA, Penta D, Penta E, TPOX, TH01 and vWA) in Pomeranian's descendants from the Espirito Santo State (ES), Brazil, third largest population of Pomeranian's descendants in the world. They chose the mountain region of the state for their preferred geographic location, and they have a very peculiar lifestyle with a selective mating behavior which has maintained their characteristics as a relatively pure subpopulation. Blood samples were obtained from 82 unrelated volunteers from 11 different cities of Espirito Santo State, where there are the Pomeranian's descendants. All 15 loci analyzed showed Power of Discrimination (PD) values > 0.75. Except the TPOX locus, all analyzed loci were at Hardy-Weinberg equilibrium. This subpopulation has not yet been characterized for STR allelic frequencies used for forensic and genetic identification studies.

HLApers: HLA Typing and Quantification of Expression with Personalized Index.

The plethora of RNA-seq data which have been generated in the recent years constitutes an attractive resource to investigate HLA variation and its relationship with normal and disease phenotypes, such as cancer. However, next generation sequencing (NGS) brings new challenges to HLA analysis because of the mapping bias introduced by aligning short reads originated from polymorphic genes to a single reference genome. Here we describe HLApers, a pipeline which adapts widely used tools for analysis of standard RNA-seq data to infer HLA genotypes and estimate expression. By generating reliable expression estimates for each HLA allele that an individual carries, HLApers allows a better understanding of the relationship between HLA alleles and phenotypes manifested by an individual.