Recommendations for Uniform Variant Calling of SARS-CoV-2 Genome Sequence across Bioinformatic Workflows.

Genomic sequencing of clinical samples to identify emerging variants of SARS-CoV-2 has been a key public health tool for curbing the spread of the virus. As a result, an unprecedented number of SARS-CoV-2 genomes were sequenced during the COVID-19 pandemic, which allowed for rapid identification of genetic variants, enabling the timely design and testing of therapies and deployment of new vaccine formulations to combat the new variants. However, despite the technological advances of deep sequencing, the analysis of the raw sequence data generated globally is neither standardized nor consistent, leading to vastly disparate sequences that may impact identification of variants. Here, we show that for both Illumina and Oxford Nanopore sequencing platforms, downstream bioinformatic protocols used by industry, government, and academic groups resulted in different virus sequences from same sample. These bioinformatic workflows produced consensus genomes with differences in single nucleotide polymorphisms, inclusion and exclusion of insertions, and/or deletions, despite using the same raw sequence as input datasets. Here, we compared and characterized such discrepancies and propose a specific suite of parameters and protocols that should be adopted across the field. Consistent results from bioinformatic workflows are fundamental to SARS-CoV-2 and future pathogen surveillance efforts, including pandemic preparation, to allow for a data-driven and timely public health response.

MYC overrides HIF-1α to regulate proliferating primary cell metabolism in hypoxia.

Hypoxia requires metabolic adaptations to sustain energetically demanding cellular activities. While the metabolic consequences of hypoxia have been studied extensively in cancer cell models, comparatively little is known about how primary cell metabolism responds to hypoxia. Thus, we developed metabolic flux models for human lung fibroblast and pulmonary artery smooth muscle cells proliferating in hypoxia. Unexpectedly, we found that hypoxia decreased glycolysis despite activation of hypoxia-inducible factor 1α (HIF-1α) and increased glycolytic enzyme expression. While HIF-1α activation in normoxia by prolyl hydroxylase (PHD) inhibition did increase glycolysis, hypoxia blocked this effect. Multi-omic profiling revealed distinct molecular responses to hypoxia and PHD inhibition, and suggested a critical role for MYC in modulating HIF-1α responses to hypoxia. Consistent with this hypothesis, MYC knockdown in hypoxia increased glycolysis and MYC over-expression in normoxia decreased glycolysis stimulated by PHD inhibition. These data suggest that MYC signaling in hypoxia uncouples an increase in HIF-dependent glycolytic gene transcription from glycolytic flux.

Towards increased accuracy and reproducibility in SARS-CoV-2 next generation sequence analysis for public health surveillance.

During the COVID-19 pandemic, SARS-CoV-2 surveillance efforts integrated genome sequencing of clinical samples to identify emergent viral variants and to support rapid experimental examination of genome-informed vaccine and therapeutic designs. Given the broad range of methods applied to generate new viral genomes, it is critical that consensus and variant calling tools yield consistent results across disparate pipelines. Here we examine the impact of sequencing technologies (Illumina and Oxford Nanopore) and 7 different downstream bioinformatic protocols on SARS-CoV-2 variant calling as part of the NIH Accelerating COVID-19 Therapeutic Interventions and Vaccines (ACTIV) Tracking Resistance and Coronavirus Evolution (TRACE) initiative, a public-private partnership established to address the COVID-19 outbreak. Our results indicate that bioinformatic workflows can yield consensus genomes with different single nucleotide polymorphisms, insertions, and/or deletions even when using the same raw sequence input datasets. We introduce the use of a specific suite of parameters and protocols that greatly improves the agreement among pipelines developed by diverse organizations. Such consistency among bioinformatic pipelines is fundamental to SARS-CoV-2 and future pathogen surveillance efforts. The application of analysis standards is necessary to more accurately document phylogenomic trends and support data-driven public health responses.

Expression of a Human Caveolin-1 Mutation in Mice Drives Inflammatory and Metabolic Defect-Associated Pulmonary Arterial Hypertension.

Background: In 2012, mutations in Cav1 were found to be the driving mutation in several cases of heritable pulmonary arterial hypertension (PAH). These mutations replaced the last 21 amino acids of Cav1 with a novel 22-amino-acid sequence. Because previously only Cav1 knockouts had been studied in the context of PAH, examining the in vivo effects of this novel mutation holds promise for new understanding of the role of Cav1 in disease etiology. Methods: The new 22 amino acids created by the human mutation were knocked into the native mouse Cav1 locus. The mice underwent hemodynamic, energy balance, and inflammatory measurements, both at baseline and after being stressed with either a metabolic or an inflammatory challenge [low-dose lipopolysaccharide (LPS)]. To metabolically challenge the mice, they were injected with streptozotocin (STZ) and fed a high-fat diet for 12 weeks. Results: Very little mutant protein was found in vivo (roughly 2% of wild-type by mass spectrometry), probably because of degradation after failure to traffic from the endoplasmic reticulum. The homozygous mutants developed a mild, low-penetrance PAH similar to that described previously in knockouts, and neither baseline nor metabolic nor inflammatory stress resulted in pressures above normal in heterozygous animals. The homozygous mutants had increased lean mass and worsened oral glucose tolerance, as previously described in knockouts. Novel findings include the preservation of Cav2 and accessory proteins in the liver and the kidney, while they are lost with homozygous Cav1 mutation in the lungs. We also found that the homozygous mutants had a significantly lower tolerance to voluntary spontaneous exercise than the wild-type mice, with the heterozygous mice at an intermediate level. The mutants also had higher circulating monocytes, with both heterozygous and homozygous animals having higher pulmonary MCP1 and MCP5 proteins. The heterozygous animals also lost weight at an LPS challenge level at which the wild-type mice continued to gain weight. Conclusions: The Cav1 mutation identified in human patients in 2012 is molecularly similar to a knockout of Cav1. It results in not only metabolic deficiencies and mild pulmonary hypertension, as expected, but also an inflammatory phenotype and reduced spontaneous exercise.

A disease-associated frameshift mutation in caveolin-1 disrupts caveolae formation and function through introduction of a de novo ER retention signal.

Caveolin-1 (CAV1) is an essential component of caveolae and is implicated in numerous physiological processes. Recent studies have identified heterozygous mutations in the CAV1 gene in patients with pulmonary arterial hypertension (PAH), but the mechanisms by which these mutations impact caveolae assembly and contribute to disease remain unclear. To address this question, we examined the consequences of a familial PAH-associated frameshift mutation in CAV1, P158PfsX22, on caveolae assembly and function. We show that C-terminus of the CAV1 P158 protein contains a functional ER-retention signal that inhibits ER exit and caveolae formation and accelerates CAV1 turnover in Cav1-/- MEFs. Moreover, when coexpressed with wild-type (WT) CAV1 in Cav1-/- MEFs, CAV1-P158 functions as a dominant negative by partially disrupting WT CAV1 trafficking. In patient skin fibroblasts, CAV1 and caveolar accessory protein levels are reduced, fewer caveolae are observed, and CAV1 complexes exhibit biochemical abnormalities. Patient fibroblasts also exhibit decreased resistance to a hypo-osmotic challenge, suggesting the function of caveolae as membrane reservoir is compromised. We conclude that the P158PfsX22 frameshift introduces a gain of function that gives rise to a dominant negative form of CAV1, defining a new mechanism by which disease-associated mutations in CAV1 impair caveolae assembly.

Caveolin-1 is an aggresome-inducing protein.

Caveolin-1 (Cav1) drives the formation of flask-shaped membrane invaginations known as caveolae that participate in signaling, clathrin-independent endocytosis and mechanotransduction. Overexpression or mutations of Cav1 can lead to its mistrafficking, including its accumulation in a perinuclear compartment previously identified as the Golgi complex. Here, we show that in the case of overexpressed Cav1-GFP, this perinuclear compartment consists of cytoplasmic inclusion bodies generated in response to the accumulation of aggregates of misfolded proteins, known as aggresomes. Aggresomes containing Cav1-GFP are encased within vimentin cages, form in a microtubule-dependent manner, and are enriched in a number of key regulators of protein turnover, including ubiquitin, VCP/p97 and proteasomes. Interestingly, aggresome induction was cell-type dependent and was observed for many but not all Cav1 constructs tested. Furthermore, endogenous Cav1 accumulated in aggresomes formed in response to proteosomal inhibition. Our finding that Cav1 is both an aggresome-inducing and aggresome-localized protein provides new insights into how cells handle and respond to misfolded Cav1. They also raise the possibility that aggresome formation may contribute to some of reported phenotypes associated with overexpressed and/or mutant forms of Cav1.

Characterization of a caveolin-1 mutation associated with both pulmonary arterial hypertension and congenital generalized lipodystrophy.

Congenital generalized lipodystrophy (CGL) and pulmonary arterial hypertension (PAH) have recently been associated with mutations in the caveolin-1 ( CAV1 ) gene, which encodes the primary structural protein of caveolae. However, little is currently known about how these CAV1 mutations impact caveolae formation or contribute to the development of disease. Here, we identify a heterozygous F160X CAV1 mutation predicted to generate a C-terminally truncated mutant protein in a patient with both PAH and CGL using whole exome sequencing, and characterize the properties of CAV1 , caveolae-associated proteins and caveolae in skin fibroblasts isolated from the patient. We show that morphologically defined caveolae are present in patient fibroblasts and that they function in mechanoprotection. However, they exhibited several notable defects, including enhanced accessibility of the C-terminus of wild-type CAV1 in caveolae, reduced colocalization of cavin-1 with CAV1 and decreased stability of both 8S and 70S oligomeric CAV1 complexes that are necessary for caveolae formation. These results were verified independently in reconstituted CAV1 -/- mouse embryonic fibroblasts. These findings identify defects in caveolae that may serve as contributing factors to the development of PAH and CGL and broaden our knowledge of CAV1 mutations associated with human disease.

A little cooperation helps murine cytomegalovirus (MCMV) go a long way: MCMV co-infection rescues a chemokine salivary gland defect.

Cytomegaloviruses (CMVs) produce chemokines (vCXCLs) that have both sequence and functional homology to host chemokines. Assessment of vCXCL-1's role in CMV infection is limited to in vitro and in silico analysis due to CMVs species specificity. In this study, we used the murine CMV (MCMV) mouse model to evaluate the function of vCXCL-1 in vivo. Recombinant MCMVs expressing chimpanzee CMV vCXCL-1 (vCXCL-1CCMV) or host chemokine, mCXCL1, underwent primary dissemination to the popliteal lymph node, spleen and lung similar to the parental MCMV. However, neither of the recombinants expressing chemokines was recovered from the salivary gland (SG) at any time post-infection although viral DNA was detected. This implies that the virus does not grow in the SG or the overexpressed chemokine induces an immune response that leads to suppressed growth. Pointing to immune suppression of virus replication, recombinant viruses were isolated from the SG following infection of immune-ablated mice [i.e. SCID (severe combined immunodeficiency), NSG (non-obese diabetic SCID gamma) or cyclophosphamide treated]. Depletion of neutrophils or NK cells does not rescue the recovery of chemokine-expressing recombinants in the SG. Surprisingly we found that co-infection of parental virus and chemokine-expressing virus leads to the recovery of the recombinants in the SG. We suggest that parental virus reduces the levels of chemokine expression leading to a decrease in inflammatory monocytes and subsequent SG growth. Therefore, aberrant expression of the chemokines induces cells of the innate and adaptive immune system that curtail the growth and dissemination of the recombinants in the SG.

Assembly and Turnover of Caveolae: What Do We Really Know?

In addition to containing highly dynamic nanoscale domains, the plasma membranes of many cell types are decorated with caveolae, flask-shaped domains enriched in the structural protein caveolin-1 (Cav1). The importance of caveolae in numerous cellular functions and processes has become well-recognized, and recent years have seen dramatic advances in our understanding of how caveolae assemble and the mechanisms control the turnover of Cav1. At the same time, work from our lab and others have revealed that commonly utilized strategies such as overexpression and tagging of Cav1 have unexpectedly complex consequences on the trafficking and fate of Cav1. Here, we discuss the implications of these findings for current models of caveolae biogenesis and Cav1 turnover. In addition, we discuss how disease-associated mutants of Cav1 impact caveolae assembly and outline open questions in this still-emerging area.

Diminished reovirus capsid stability alters disease pathogenesis and littermate transmission.

Reovirus is a nonenveloped mammalian virus that provides a useful model system for studies of viral infections in the young. Following internalization into host cells, the outermost capsid of reovirus virions is removed by endosomal cathepsin proteases. Determinants of capsid disassembly kinetics reside in the viral σ3 protein. However, the contribution of capsid stability to reovirus-induced disease is unknown. In this study, we found that mice inoculated intramuscularly with a serotype 3 reovirus containing σ3-Y354H, a mutation that reduces viral capsid stability, succumbed at a higher rate than those infected with wild-type virus. At early times after inoculation, σ3-Y354H virus reached higher titers than wild-type virus at several sites within the host. Animals inoculated perorally with a serotype 1 reassortant reovirus containing σ3-Y354H developed exaggerated myocarditis accompanied by elaboration of pro-inflammatory cytokines. Surprisingly, unchallenged littermates of mice infected with σ3-Y354H virus displayed higher titers in the intestine, heart, and brain than littermates of mice inoculated with wild-type virus. Together, these findings suggest that diminished capsid stability enhances reovirus replication, dissemination, lethality, and host-to-host spread, establishing a new virulence determinant for nonenveloped viruses.

Microtubule motors power plasma membrane tubulation in clathrin-independent endocytosis.

How the plasma membrane is bent to accommodate clathrin-independent endocytosis remains uncertain. Recent studies suggest Shiga and cholera toxin induce membrane curvature required for their uptake into clathrin-independent carriers by binding and cross-linking multiple copies of their glycosphingolipid receptors on the plasma membrane. But it remains unclear if toxin-induced sphingolipid crosslinking provides sufficient mechanical force for deforming the plasma membrane, or if host cell factors also contribute to this process. To test this, we imaged the uptake of cholera toxin B-subunit into surface-derived tubular invaginations. We found that cholera toxin mutants that bind to only one glycosphingolipid receptor accumulated in tubules, and that toxin binding was entirely dispensable for membrane tubulations to form. Unexpectedly, the driving force for tubule extension was supplied by the combination of microtubules, dynein and dynactin, thus defining a novel mechanism for generating membrane curvature during clathrin-independent endocytosis.

B cell responses to HIV antigen are a potent correlate of viremia in HIV-1 infection and improve with PD-1 blockade.

Infection with Human Immunodeficiency Virus Type 1 (HIV-1) induces defects of both cellular and humoral immune responses. Impaired CD4+ T cell help and B cell dysfunction may partially explain the low frequency of broadly neutralizing antibodies in HIV-infected individuals. To understand the extent of B cell dysfunction during HIV infection, we assessed the level of B cell activation at baseline and after stimulation with a variety of antigens. Increased levels of viremia were associated with higher baseline expression of the activation marker CD86 on B cells and with decreased ability of B cells to increase expression of CD86 after in vitro stimulation with inactivated HIV-1. In a series of cell isolation experiments B cell responses to antigen were enhanced in the presence of autologous CD4+ T cells. HIV infected individuals had a higher frequency of PD-1 expression on B cells compared to HIV- subjects and PD-1 blockade improved B cell responsiveness to HIV antigen, suggesting that inhibitory molecule expression during HIV-1 infection may contribute to some of the observed B cell defects. Our findings demonstrate that during chronic HIV infection, B cells are activated and lose full capacity to respond to antigen, but suppression of inhibitory pressures as well as a robust CD4+ T cell response may help preserve B cell function.