Genome-wide Screens Identify Lineage- and Tumor-Specific Genes Modulating MHC-I- and MHC-II-Restricted Immunosurveillance of Human Lymphomas.

Tumors frequently subvert major histocompatibility complex class I (MHC-I) peptide presentation to evade CD8+ T cell immunosurveillance, though how this is accomplished is not always well defined. To identify the global regulatory networks controlling antigen presentation, we employed genome-wide screening in human diffuse large B cell lymphomas (DLBCLs). This approach revealed dozens of genes that positively and negatively modulate MHC-I cell surface expression. Validated genes clustered in multiple pathways including cytokine signaling, mRNA processing, endosomal trafficking, and protein metabolism. Genes can exhibit lymphoma subtype- or tumor-specific MHC-I regulation, and a majority of primary DLBCL tumors displayed genetic alterations in multiple regulators. We established SUGT1 as a major positive regulator of both MHC-I and MHC-II cell surface expression. Further, pharmacological inhibition of two negative regulators of antigen presentation, EZH2 and thymidylate synthase, enhanced DLBCL MHC-I presentation. These and other genes represent potential targets for manipulating MHC-I immunosurveillance in cancers, infectious diseases, and autoimmunity.

RGS4 controls airway hyperresponsiveness through GAP-independent mechanisms.

Regulators of G protein signaling (RGS) proteins constrain G protein-coupled receptor (GPCR)-mediated and other responses throughout the body primarily, but not exclusively, through their GTPase-activating protein activity. Asthma is a highly prevalent condition characterized by airway hyper-responsiveness (AHR) to environmental stimuli resulting in part from amplified GPCR-mediated airway smooth muscle contraction. Rgs2 or Rgs5 gene deletion in mice enhances AHR and airway smooth muscle contraction, whereas RGS4 KO mice unexpectedly have decreased AHR because of increased production of the bronchodilator prostaglandin E2 (PGE2) by lung epithelial cells. Here, we found that knockin mice harboring Rgs4 alleles encoding a point mutation (N128A) that sharply curtails RGS4 GTPase-activating protein activity had increased AHR, reduced airway PGE2 levels, and augmented GPCR-induced bronchoconstriction compared with either RGS4 KO mice or WT controls. RGS4 interacted with the p85α subunit of PI3K and inhibited PI3K-dependent PGE2 secretion elicited by transforming growth factor beta in airway epithelial cells. Together, these findings suggest that RGS4 affects asthma severity in part by regulating the airway inflammatory milieu in a G protein-independent manner.

A Population Genomic Assessment of Three Decades of Evolution in a Natural Drosophila Population.

Population genetics seeks to illuminate the forces shaping genetic variation, often based on a single snapshot of genomic variation. However, utilizing multiple sampling times to study changes in allele frequencies can help clarify the relative roles of neutral and non-neutral forces on short time scales. This study compares whole-genome sequence variation of recently collected natural population samples of Drosophila melanogaster against a collection made approximately 35 years prior from the same locality-encompassing roughly 500 generations of evolution. The allele frequency changes between these time points would suggest a relatively small local effective population size on the order of 10,000, significantly smaller than the global effective population size of the species. Some loci display stronger allele frequency changes than would be expected anywhere in the genome under neutrality-most notably the tandem paralogs Cyp6a17 and Cyp6a23, which are impacted by structural variation associated with resistance to pyrethroid insecticides. We find a genome-wide excess of outliers for high genetic differentiation between old and new samples, but a larger number of adaptation targets may have affected SNP-level differentiation versus window differentiation. We also find evidence for strengthening latitudinal allele frequency clines: northern-associated alleles have increased in frequency by an average of nearly 2.5% at SNPs previously identified as clinal outliers, but no such pattern is observed at random SNPs. This project underscores the scientific potential of using multiple sampling time points to investigate how evolution operates in natural populations, by quantifying how genetic variation has changed over ecologically relevant timescales.

Type I IFNs facilitate innate immune control of the opportunistic bacteria Burkholderia cenocepacia in the macrophage cytosol.

The mammalian immune system is constantly challenged by signals from both pathogenic and non-pathogenic microbes. Many of these non-pathogenic microbes have pathogenic potential if the immune system is compromised. The importance of type I interferons (IFNs) in orchestrating innate immune responses to pathogenic microbes has become clear in recent years. However, the control of opportunistic pathogens-and especially intracellular bacteria-by type I IFNs remains less appreciated. In this study, we use the opportunistic, Gram-negative bacterial pathogen Burkholderia cenocepacia (Bc) to show that type I IFNs are capable of limiting bacterial replication in macrophages, preventing illness in immunocompetent mice. Sustained type I IFN signaling through cytosolic receptors allows for increased expression of autophagy and linear ubiquitination mediators, which slows bacterial replication. Transcriptomic analyses and in vivo studies also show that LPS stimulation does not replicate the conditions of intracellular Gram-negative bacterial infection as it pertains to type I IFN stimulation or signaling. This study highlights the importance of type I IFNs in protection against opportunistic pathogens through innate immunity, without the need for damaging inflammatory responses.

Correction: Type I IFNs facilitate innate immune control of the opportunistic bacteria Burkholderia cenocepacia in the macrophage cytosol.

[This corrects the article DOI: 10.1371/journal.ppat.1009395.].

Markers of endothelial glycocalyx dysfunction in Clarkson disease.

BACKGROUND: Clarkson disease (monoclonal gammopathy-associated idiopathic systemic capillary leak syndrome, ISCLS) is a rare idiopathic condition marked by transient, relapsing-remitting episodes of systemic microvascular hyper-permeability, which liberates plasma fluid and macromolecules into the peripheral tissues. This pathology manifests clinically as the abrupt onset of hypotensive shock, hemoconcentration, and hypoalbuminemia. METHODS: We analysed endothelial glycocalyx (eGCX)-related markers in plasma from patients with ISCLS during acute disease flares and convalescence by ELISA and comprehensive proteomic profiling. We evaluated eGCX-related components and gene expression in cultured endothelial cells using RNA-sequencing, real-time PCR, and fluorescence staining. RESULTS: Serum levels of eGCX-related core components including hyaluronic acid (HA) and the core proteoglycan soluble syndecan-1 (sCD138) were elevated at baseline and during acute ISCLS flares. Serial measurements demonstrated that sCD138 levels peaked during the recovery (post-leak) phase of the illness. Proteomic analysis of matched acute and convalescent ISCLS plasma revealed increased abundance of eGCX-related proteins, including glypicans, thrombospondin-1 (TSP-1), and eGCX-degrading enzymes in acute compared to remission plasma. Abundance of endothelial cell damage markers did not differ in acute and baseline plasma. Expression of several eGCX-related genes and surface carbohydrate content in endothelial cells from patients with ISCLS did not differ significantly from that observed in healthy control cells. CONCLUSIONS: eGCX dysfunction, but not endothelial injury, may contribute to clinical symptoms of acute ISCLS. Serum levels of of eGCX components including sCD138 may be measured during acute episodes of ISCLS to monitor clinical status and therapeutic responses.

Decanalization of wing development accompanied the evolution of large wings in high-altitude Drosophila.

In higher organisms, the phenotypic impacts of potentially harmful or beneficial mutations are often modulated by complex developmental networks. Stabilizing selection may favor the evolution of developmental canalization--that is, robustness despite perturbation--to insulate development against environmental and genetic variability. In contrast, directional selection acts to alter the developmental process, possibly undermining the molecular mechanisms that buffer a trait's development, but this scenario has not been shown in nature. Here, we examined the developmental consequences of size increase in highland Ethiopian Drosophila melanogaster. Ethiopian inbred strains exhibited much higher frequencies of wing abnormalities than lowland populations, consistent with an elevated susceptibility to the genetic perturbation of inbreeding. We then used mutagenesis to test whether Ethiopian wing development is, indeed, decanalized. Ethiopian strains were far more susceptible to this genetic disruption of development, yielding 26 times more novel wing abnormalities than lowland strains in F2 males. Wing size and developmental perturbability cosegregated in the offspring of between-population crosses, suggesting that genes conferring size differences had undermined developmental buffering mechanisms. Our findings represent the first observation, to our knowledge, of morphological evolution associated with decanalization in the same tissue, underscoring the sensitivity of development to adaptive change.

Parallel Evolution of Cold Tolerance within Drosophila melanogaster.

Drosophila melanogaster originated in tropical Africa before expanding into strikingly different temperate climates in Eurasia and beyond. Here, we find elevated cold tolerance in three distinct geographic regions: beyond the well-studied non-African case, we show that populations from the highlands of Ethiopia and South Africa have significantly increased cold tolerance as well. We observe greater cold tolerance in outbred versus inbred flies, but only in populations with higher inversion frequencies. Each cold-adapted population shows lower inversion frequencies than a closely-related warm-adapted population, suggesting that inversion frequencies may decrease with altitude in addition to latitude. Using the FST-based "Population Branch Excess" statistic (PBE), we found only limited evidence for parallel genetic differentiation at the scale of ∼4 kb windows, specifically between Ethiopian and South African cold-adapted populations. And yet, when we looked for single nucleotide polymorphisms (SNPs) with codirectional frequency change in two or three cold-adapted populations, strong genomic enrichments were observed from all comparisons. These findings could reflect an important role for selection on standing genetic variation leading to "soft sweeps". One SNP showed sufficient codirectional frequency change in all cold-adapted populations to achieve experiment-wide significance: an intronic variant in the synaptic gene Prosap. Another codirectional outlier SNP, at senseless-2, had a strong association with our cold trait measurements, but in the opposite direction as predicted. More generally, proteins involved in neurotransmission were enriched as potential targets of parallel adaptation. The ability to study cold tolerance evolution in a parallel framework will enhance this classic study system for climate adaptation.

A Thousand Fly Genomes: An Expanded Drosophila Genome Nexus.

The Drosophila Genome Nexus is a population genomic resource that provides D. melanogaster genomes from multiple sources. To facilitate comparisons across data sets, genomes are aligned using a common reference alignment pipeline which involves two rounds of mapping. Regions of residual heterozygosity, identity-by-descent, and recent population admixture are annotated to enable data filtering based on the user's needs. Here, we present a significant expansion of the Drosophila Genome Nexus, which brings the current data object to a total of 1,121 wild-derived genomes. New additions include 305 previously unpublished genomes from inbred lines representing six population samples in Egypt, Ethiopia, France, and South Africa, along with another 193 genomes added from recently-published data sets. We also provide an aligned D. simulans genome to facilitate divergence comparisons. This improved resource will broaden the range of population genomic questions that can addressed from multi-population allele frequencies and haplotypes in this model species. The larger set of genomes will also enhance the discovery of functionally relevant natural variation that exists within and between populations.

C-C motif chemokine receptor 2 and 7 synergistically control inflammatory monocyte recruitment but the infecting virus dictates monocyte function in the brain.

Inflammatory monocytes (iMO) are recruited from the bone marrow to the brain during viral encephalitis. C-C motif chemokine receptor (CCR) 2 deficiency substantially reduces iMO recruitment for most, but not all encephalitic viruses. Here we show CCR7 acts synergistically with CCR2 to control this process. Following Herpes simplex virus type-1 (HSV-1), or La Crosse virus (LACV) infection, we find iMO proportions are reduced by approximately half in either Ccr2 or Ccr7 knockout mice compared to control mice. However, Ccr2/Ccr7 double knockouts eliminate iMO recruitment following infection with either virus, indicating these receptors together control iMO recruitment. We also find that LACV induces a more robust iMO recruitment than HSV-1. However, unlike iMOs in HSV-1 infection, LACV-recruited iMOs do not influence neurological disease development. LACV-induced iMOs have higher expression of proinflammatory and proapoptotic but reduced mitotic, phagocytic and phagolysosomal transcripts compared to HSV-1-induced iMOs. Thus, virus-specific activation of iMOs affects their recruitment, activation, and function.

Intrinsic endothelial hyperresponsiveness to inflammatory mediators drives acute episodes in models of Clarkson disease.

Clarkson disease, or monoclonal gammopathy-associated idiopathic systemic capillary leak syndrome (ISCLS), is a rare, relapsing-remitting disorder featuring the abrupt extravasation of fluids and proteins into peripheral tissues, which in turn leads to hypotensive shock, severe hemoconcentration, and hypoalbuminemia. The specific leakage factor(s) and pathways in ISCLS are unknown, and there is no effective treatment for acute flares. Here, we characterize an autonomous vascular endothelial defect in ISCLS that was recapitulated in patient-derived endothelial cells (ECs) in culture and in a mouse model of disease. ISCLS-derived ECs were functionally hyperresponsive to permeability-inducing factors like VEGF and histamine, in part due to increased endothelial nitric oxide synthase (eNOS) activity. eNOS blockade by administration of N(γ)-nitro-l-arginine methyl ester (l-NAME) ameliorated vascular leakage in an SJL/J mouse model of ISCLS induced by histamine or VEGF challenge. eNOS mislocalization and decreased protein phosphatase 2A (PP2A) expression may contribute to eNOS hyperactivation in ISCLS-derived ECs. Our findings provide mechanistic insights into microvascular barrier dysfunction in ISCLS and highlight a potential therapeutic approach.

Lack of association between classical HLA genes and asymptomatic SARS-CoV-2 infection.

Human genetic studies of critical COVID-19 pneumonia have revealed the essential role of type I interferon-dependent innate immunity to SARS-CoV-2 infection. Conversely, an association between the HLA-B∗15:01 allele and asymptomatic SARS-CoV-2 infection in unvaccinated individuals was recently reported, suggesting a contribution of pre-existing T cell-dependent adaptive immunity. We report a lack of association of classical HLA alleles, including HLA-B∗15:01, with pre-omicron asymptomatic SARS-CoV-2 infection in unvaccinated participants in a prospective population-based study in the United States (191 asymptomatic vs. 945 symptomatic COVID-19 cases). Moreover, we found no such association in the international COVID Human Genetic Effort cohort (206 asymptomatic vs. 574 mild or moderate COVID-19 cases and 1,625 severe or critical COVID-19 cases). Finally, in the Human Challenge Characterisation study, the three HLA-B∗15:01 individuals infected with SARS-CoV-2 developed symptoms. As with other acute primary infections studied, no classical HLA alleles favoring an asymptomatic course of SARS-CoV-2 infection were identified.

A Simple, Rapid, and Effective Heparinase Protocol to Enable Nucleic Acid Study from Frozen Heparinized Plasma.

Cell-free RNAs (cfRNAs) are promising analytes as non-invasive biomarkers and have even greater potential if tied in with metabolomics. Plasma is an optimal source for cfRNAs but is often derived from a variety of anticoagulants. Plasma obtained in heparin is suitable for metabolomics but is difficult to utilize for qPCR-based downstream analysis. In the present study, we aimed to develop a simple, time-efficient, and cost-effective heparinase protocol, followed by library preparation and sequencing of human plasma cfRNAs drawn and stored in heparin at -80 °C for several years. Blood was collected in CPT™ sodium heparin tubes from patients with chronic HCV infection (NCT02400216) at the National Institutes of Health (NIH) Clinical Center. Plasma cfRNAs were treated with heparinase I and used for library preparation and next-generation sequencing (NGS). Heparinase treatment maintained RNA integrity and allowed for successful library preparation for all the study subjects even with 7 ng of cfRNAs as starting material. The classification report derived from Pavian R package v1.2.0 showed no artificial reads. The abundance of chordate over microbial reads suggests no addition of experimental error through heparinase I treatment. We report a novel and practical approach to heparinase treatment for human plasma collected and frozen in sodium heparin for several years. This is an effective demonstration of utilizing heparin plasma for NGS and downstream transcriptomic research, which could then be integrated with metabolomics from the same samples, maximizing efficiency and minimizing blood draws.

Genetically diverse mouse models of SARS-CoV-2 infection reproduce clinical variation in type I interferon and cytokine responses in COVID-19.

Inflammation in response to severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) infection drives severity of coronavirus disease 2019 (COVID-19) and is influenced by host genetics. To understand mechanisms of inflammation, animal models that reflect genetic diversity and clinical outcomes observed in humans are needed. We report a mouse panel comprising the genetically diverse Collaborative Cross (CC) founder strains crossed to human ACE2 transgenic mice (K18-hACE2) that confers susceptibility to SARS-CoV-2. Infection of CC x K18- hACE2 resulted in a spectrum of survival, viral replication kinetics, and immune profiles. Importantly, in contrast to the K18-hACE2 model, early type I interferon (IFN-I) and regulated proinflammatory responses were required for control of SARS-CoV-2 replication in PWK x K18-hACE2 mice that were highly resistant to disease. Thus, virus dynamics and inflammation observed in COVID-19 can be modeled in diverse mouse strains that provide a genetically tractable platform for understanding anti-coronavirus immunity.

Genetically diverse mouse models of SARS-CoV-2 infection reproduce clinical variation in type I interferon and cytokine responses in COVID-19.

Inflammation in response to severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) infection drives severity of coronavirus disease 2019 (COVID-19) and is influenced by host genetics. To understand mechanisms of inflammation, animal models that reflect genetic diversity and clinical outcomes observed in humans are needed. We report a mouse panel comprising the genetically diverse Collaborative Cross (CC) founder strains crossed to human ACE2 transgenic mice (K18-hACE2) that confers susceptibility to SARS-CoV-2. Infection of CC x K18-hACE2 resulted in a spectrum of survival, viral replication kinetics, and immune profiles. Importantly, in contrast to the K18-hACE2 model, early type I interferon (IFN-I) and regulated proinflammatory responses were required for control of SARS-CoV-2 replication in PWK x K18-hACE2 mice that were highly resistant to disease. Thus, virus dynamics and inflammation observed in COVID-19 can be modeled in diverse mouse strains that provide a genetically tractable platform for understanding anti-coronavirus immunity.

Lack of association between HLA and asymptomatic SARS-CoV-2 infection.

Human genetic studies of critical COVID-19 pneumonia have revealed the essential role of type I interferon-dependent innate immunity to SARS-CoV-2 infection. Conversely, an association between the HLA-B*15:01 allele and asymptomatic SARS-CoV-2 infection in unvaccinated individuals was recently reported, suggesting a contribution of pre-existing T cell-dependent adaptive immunity. We report a lack of association of classical HLA alleles, including HLA-B*15:01, with pre-omicron asymptomatic SARS-CoV-2 infection in unvaccinated participants in a prospective population-based study in the US (191 asymptomatic vs. 945 symptomatic COVID-19 cases). Moreover, we found no such association in the international COVID Human Genetic Effort cohort (206 asymptomatic vs. 574 mild or moderate COVID-19 cases and 1,625 severe or critical COVID-19 cases). Finally, in the Human Challenge Characterisation study, the three HLA-B*15:01 individuals infected with SARS-CoV-2 developed symptoms. As with other acute primary infections, no classical HLA alleles favoring an asymptomatic course of SARS-CoV-2 infection were identified. These findings suggest that memory T-cell immunity to seasonal coronaviruses does not strongly influence the outcome of SARS-CoV-2 infection in unvaccinated individuals.

Hepatitis E virus genotype 3 in wild rats, United States.

The role of rodents in the epidemiology of zoonotic hepatitis E virus (HEV) infection has been a subject of considerable debate. Seroprevalence studies suggest widespread HEV infection in commensal Rattus spp. rats, but experimental transmission has been largely unsuccessful and recovery of zoonotic genotype 3 HEV RNA from wild Rattus spp. rats has never been confirmed. We surveyed R. rattus and R. norvegicus rats from across the United States and several international populations by using a hemi-nested reverse transcription PCR approach. We isolated HEV RNA in liver tissues from 35 of 446 rats examined. All but 1 of these isolates was relegated to the zoonotic HEV genotype 3, and the remaining sequence represented the recently discovered rat genotype from the United States and Germany. HEV-positive rats were detected in urban and remote localities. Genetic analyses suggest all HEV genotype 3 isolates obtained from wild Rattus spp. rats were closely related.

Invasion facilitates hybridization with introgression in the Rattus rattus species complex.

Biological invasions result in novel species interactions, which can have significant evolutionary impacts on both native and invading taxa. One evolutionary concern with invasions is hybridization among lineages that were previously isolated, but make secondary contact in their invaded range(s). Black rats, consisting of several morphologically very similar but genetically distinct taxa that collectively have invaded six continents, are arguably the most successful mammalian invaders on the planet. We used mitochondrial cytochrome b sequences, two nuclear gene sequences (Atp5a1 and DHFR) and nine microsatellite loci to examine the distribution of three invasive black rat lineages (Rattus tanezumi, Rattus rattus I and R. rattus IV) in the United States and Asia and to determine the extent of hybridization among these taxa. Our analyses revealed two mitochondrial lineages that have spread to multiple continents, including a previously undiscovered population of R. tanezumi in the south-eastern United States, whereas the third lineage (R. rattus IV) appears to be confined to Southeast Asia. Analyses of nuclear DNA (both sequences and microsatellites) suggested significant hybridization is occurring among R. tanezumi and R. rattus I in the United States and also suggest hybridization between R. tanezumi and R. rattus IV in Asia, although further sampling of the latter species pair in Asia is required. Furthermore, microsatellite analyses suggest unidirectional introgression from both R. rattus I and R. rattus IV into R. tanezumi. Within the United States, introgression appears to be occurring to such a pronounced extent that we were unable to detect any nuclear genetic signal for R. tanezumi, and a similar pattern was detected in Asia.

Gene Regulatory Evolution in Cold-Adapted Fly Populations Neutralizes Plasticity and May Undermine Genetic Canalization.

The relationships between adaptive evolution, phenotypic plasticity, and canalization remain incompletely understood. Theoretical and empirical studies have made conflicting arguments on whether adaptive evolution may enhance or oppose the plastic response. Gene regulatory traits offer excellent potential to study the relationship between plasticity and adaptation, and they can now be studied at the transcriptomic level. Here, we take advantage of three closely related pairs of natural populations of Drosophila melanogaster from contrasting thermal environments that reflect three separate instances of cold tolerance evolution. We measure the transcriptome-wide plasticity in gene expression levels and alternative splicing (intron usage) between warm and cold laboratory environments. We find that suspected adaptive changes in both gene expression and alternative splicing tend to neutralize the ancestral plastic response. Further, we investigate the hypothesis that adaptive evolution can lead to decanalization of selected gene regulatory traits. We find strong evidence that suspected adaptive gene expression (but not splicing) changes in cold-adapted populations are more vulnerable to the genetic perturbation of inbreeding than putatively neutral changes. We find some evidence that these patterns may reflect a loss of genetic canalization accompanying adaptation, although other processes including hitchhiking recessive deleterious variants may contribute as well. Our findings augment our understanding of genetic and environmental effects on gene regulation in the context of adaptive evolution.

The evolution of larger size in high-altitude Drosophila melanogaster has a variable genetic architecture.

Important uncertainties persist regarding the genetic architecture of adaptive trait evolution in natural populations, including the number of genetic variants involved, whether they are drawn from standing genetic variation, and whether directional selection drives them to complete fixation. Here, we take advantage of a unique natural population of Drosophila melanogaster from the Ethiopian highlands, which has evolved larger body size than any other known population of this species. We apply a bulk segregant quantitative trait locus mapping approach to 4 unique crosses between highland Ethiopian and lowland Zambian populations for both thorax length and wing length. Results indicated a persistently variable genetic basis for these evolved traits (with largely distinct sets of quantitative trait loci for each cross), and at least a moderately polygenic architecture with relatively strong effects present. We complemented these mapping experiments with population genetic analyses of quantitative trait locus regions and gene ontology enrichment analysis, generating strong hypotheses for specific genes and functional processes that may have contributed to these adaptive trait changes. Finally, we find that the genetic architectures indicated by our quantitative trait locus mapping results for size traits mirror those from similar experiments on other recently evolved traits in this species. Collectively, these studies suggest a recurring pattern of polygenic adaptation in this species, in which causative variants do not approach fixation and moderately strong effect loci are present.