HIV-1 control in vivo is related to the number but not the fraction of infected cells with viral unspliced RNA.

In the absence of antiretroviral therapy (ART), a subset of individuals, termed HIV controllers, have levels of plasma viremia that are orders of magnitude lower than non-controllers (NC) who are at higher risk for HIV disease progression. In addition to having fewer infected cells resulting in fewer cells with HIV RNA, it is possible that lower levels of plasma viremia in controllers are due to a lower fraction of the infected cells having HIV-1 unspliced RNA (HIV usRNA) compared with NC. To directly test this possibility, we used sensitive and quantitative single-cell sequencing methods to compare the fraction of infected cells that contain one or more copies of HIV usRNA in peripheral blood mononuclear cells (PBMC) obtained from controllers and NC. The fraction of infected cells containing HIV usRNA did not differ between the two groups. Rather, the levels of viremia were strongly associated with the total number of infected cells that had HIV usRNA, as reported by others, with controllers having 34-fold fewer infected cells per million PBMC. These results reveal that viremic control is not associated with a lower fraction of proviruses expressing HIV usRNA, unlike what is reported for elite controllers, but is only related to having fewer infected cells overall, maybe reflecting greater immune clearance of infected cells. Our findings show that proviral silencing is not a key mechanism for viremic control and will help to refine strategies toward achieving HIV remission without ART.

Longitudinal Assessment of BNT162b2- and mRNA-1273-Induced Anti-SARS-CoV-2 Spike IgG Levels and Avidity Following Three Doses of Vaccination.

SARS-CoV-2 vaccination-induced protection against infection is likely to be affected by functional antibody features. To understand the kinetics of antibody responses in healthy individuals after primary series and third vaccine doses, sera from the recipients of the two licensed SARS-CoV-2 mRNA vaccines were assessed for circulating anti-SARS-CoV-2 spike IgG levels and avidity for up to 6 months post-primary series and 9 months after the third dose. Following primary series vaccination, anti-SARS-CoV-2 spike IgG levels declined from months 1 to 6, while avidity increased through month 6, irrespective of the vaccine received. The third dose of either vaccine increased anti-SARS-CoV-2 spike IgG levels and avidity and appeared to enhance antibody level persistence-generating a slower rate of decline in the 3 months following the third dose compared to the decline seen after the primary series alone. The third dose of both vaccines induced significant avidity increases 1 month after vaccination compared to the avidity response 6 months post-primary series vaccination (p ≤ 0.001). A significant difference in avidity responses between the two vaccines was observed 6 months post-third dose, where the BNT162b2 recipients had higher antibody avidity levels compared to the mRNA-1273 recipients (p = 0.020).

MoCoLo: a testing framework for motif co-localization.

Sequence-level data offers insights into biological processes through the interaction of two or more genomic features from the same or different molecular data types. Within motifs, this interaction is often explored via the co-occurrence of feature genomic tracks using fixed-segments or analytical tests that respectively require window size determination and risk of false positives from over-simplified models. Moreover, methods for robustly examining the co-localization of genomic features, and thereby understanding their spatial interaction, have been elusive. We present a new analytical method for examining feature interaction by introducing the notion of reciprocal co-occurrence, define statistics to estimate it and hypotheses to test for it. Our approach leverages conditional motif co-occurrence events between features to infer their co-localization. Using reverse conditional probabilities and introducing a novel simulation approach that retains motif properties (e.g. length, guanine-content), our method further accounts for potential confounders in testing. As a proof-of-concept, motif co-localization (MoCoLo) confirmed the co-occurrence of histone markers in a breast cancer cell line. As a novel analysis, MoCoLo identified significant co-localization of oxidative DNA damage within non-B DNA-forming regions that significantly differed between non-B DNA structures. Altogether, these findings demonstrate the potential utility of MoCoLo for testing spatial interactions between genomic features via their co-localization.

HIV infected CD4+ T cell clones are more stable than uninfected clones during long-term antiretroviral therapy.

Although combination antiretroviral therapy (ART) blocks HIV replication, it is not curative because infected CD4+ T cells that carry intact, infectious proviruses persist. Understanding the behavior of clones of infected T cells is important for understanding the stability of the reservoir; however, the stabilities of clones of infected T cells in persons on long-term ART are not well defined. We determined the relative stabilities of clones of infected and uninfected CD4+ T cells over time intervals of one to four years in three individuals who had been on ART for 9-19 years. The largest clones of uninfected T cells were larger than the largest clones of infected T cells. Clones of infected CD4+ T cells were more stable than clones of uninfected CD4+ T cells of a similar size. Individual clones of CD4+ T cells carrying intact, infectious proviruses can expand, contract, or remain stable over time.

Comparative microsomal proteomics of a model lung cancer cell line NCI-H23 reveals distinct differences between molecular profiles of 3D and 2D cultured cells.

Lung cancer is the leading cause of cancer-related deaths in the USA and worldwide. Yet, about 95% of new drug candidates validated in preclinical phase eventually fail in clinical trials. Such a high attrition rate is attributed mostly to the inability of conventional two-dimensionally (2D) cultured cancer cells to mimic native three-dimensional (3D) growth of malignant cells in human tumors. To ascertain phenotypical differences between these two distinct culture conditions, we carried out a comparative proteomic analysis of a membrane fraction obtained from 3D- and 2D-cultured NSCLC model cell line NCI-H23. This analysis revealed a map of 1,166 (24%) protein species regulated in culture dependent manner, including differential regulation of a subset of cell surface-based CD molecules. We confirmed exclusive expression of CD99, CD146 and CD239 in 3D culture. Furthermore, label-free quantitation, targeting KRas proteoform-specific peptides, revealed upregulation of both wild type and monoallelic KRas4BG12C mutant at the surface of 3D cultured cells. In order to reduce the high attrition rate of new drug candidates, the results of this study strongly suggests exploiting base-line molecular profiling of a large number of patient-derived NSCLC cell lines grown in 2D and 3D culture, prior to actual drug candidate testing.

AVIA 3.0: interactive portal for genomic variant and sample level analysis.

SUMMARY: The Annotation, Visualization and Impact Analysis (AVIA) is a web application combining multiple features to annotate and visualize genomic variant data. Users can investigate functional significance of their genetic alterations across samples, genes and pathways. Version 3.0 of AVIA offers filtering options through interactive charts and by linking disease relevant data sources. Newly incorporated services include gene, variant and sample level reporting, literature and functional correlations among impacted genes, comparative analysis across samples and against data sources such as TCGA and ClinVar, and cohort building. Sample and data management is now feasible through the application, which allows greater flexibility with sharing, reannotating and organizing data. Most importantly, AVIA's utility stems from its convenience for allowing users to upload and explore results without any a priori knowledge or the need to install, update and maintain software or databases. Together, these enhancements strengthen AVIA as a comprehensive, user-driven variant analysis portal. AVAILABILITYAND IMPLEMENTATION: AVIA is accessible online at https://avia-abcc.ncifcrf.gov.

Electrostatically Driven Guanidinium Interaction Domains that Control Hydrogel-Mediated Protein Delivery In Vivo.

Protein biologics are an important class of drugs, but the necessity for frequent parenteral administration is a major limitation. Drug-delivery materials offer a potential solution, but protein-material adsorption can cause denaturation, which reduces their effectiveness. Here, we describe a new protein delivery platform that limits direct contact between globular protein domains and material matrix, yet from a single subcutaneous administration can be tuned for long-term drug release. The strategy utilizes complementary electrostatic interactions made between a suite of designed interaction domains (IDs), installed onto the terminus of a protein of interest, and a negatively charged self-assembled fibrillar hydrogel. These intermolecular interactions can be easily modulated by choice of ID to control material interaction and desorption energies, which allows regulation of protein release kinetics to fit desired release profiles. Molecular dynamics studies provided a molecular-level understanding of the mechanisms that govern release and identified optimal binding zones on the gel fibrils that facilitate strong ID-material interactions, which are crucial for sustained release of protein. This delivery platform can be easily loaded with cargo, is shear-thin syringe implantable, provides improved protein stability, is capable of a diverse range of in vitro release rates, and most importantly, can accomplish long-term control over in vivo protein delivery.

Clones of infected cells arise early in HIV-infected individuals.

In HIV-infected individuals on long-term antiretroviral therapy (ART), more than 40% of the infected cells are in clones. Although most HIV proviruses present in individuals on long-term ART are defective, including those in clonally expanded cells, there is increasing evidence that clones carrying replication-competent proviruses are common in patients on long-term ART and form part of the HIV reservoir that makes it impossible to cure HIV infection with current ART alone. Given the importance of clonal expansion in HIV persistence, we determined how soon after HIV acquisition infected clones can grow large enough to be detected (clones larger than ca. 1 × 105 cells). We studied 12 individuals sampled in early HIV infection (Fiebig stage III-V/VI) and 5 who were chronically infected. The recently infected individuals were started on ART at or near the time of diagnosis. We isolated more than 6,500 independent integration sites from peripheral blood mononuclear cells before ART was initiated and after 0.5-18 years of suppressive ART. Some infected clones could be detected approximately 4 weeks after HIV infection and some of these clones persisted for years. The results help to explain how the reservoir is established early and persists for years.

Direct molecular dissection of tumor parenchyma from tumor stroma in tumor xenograft using mass spectrometry-based glycoproteomics.

The most widely used cancer animal model is the human-murine tumor xenograft. Unbiased molecular dissection of tumor parenchyma versus stroma in human-murine xenografts is critical for elucidating dysregulated protein networks/pathways and developing therapeutics that may target these two functionally codependent compartments. Although antibody-reliant technologies (e.g., immunohistochemistry, imaging mass cytometry) are capable of distinguishing tumor-proper versus stromal proteins, the breadth or extent of targets is limited. Here, we report an antibody-free targeted cross-species glycoproteomic (TCSG) approach that enables direct dissection of human tumor parenchyma from murine tumor stroma at the molecular/protein level in tumor xenografts at a selectivity rate presently unattainable by other means. This approach was used to segment/dissect and obtain the protein complement phenotype of the tumor stroma and parenchyma of the metastatic human lung adenocarcinoma A549 xenograft, with no need for tissue microdissection prior to mass-spectrometry analysis. An extensive molecular map of the tumor proper and the associated microenvironment was generated along with the top functional N-glycosylated protein networks enriched in each compartment. Importantly, immunohistochemistry-based cross-validation of selected parenchymal and stromal targets applied on human tissue samples of lung adenocarcinoma and normal adjacent tissue is indicative of a noteworthy translational capacity for this unique approach that may facilitate identifications of novel targets for next generation antibody therapies and development of real time preclinical tumor models.

Proviruses with identical sequences comprise a large fraction of the replication-competent HIV reservoir.

The major obstacle to curing HIV infection is the persistence of cells with intact proviruses that can produce replication-competent virus. This HIV reservoir is believed to exist primarily in CD4+ T-cells and is stable despite years of suppressive antiretroviral therapy. A potential mechanism for HIV persistence is clonal expansion of infected cells, but how often such clones carry replication-competent proviruses has been controversial. Here, we used single-genome sequencing to probe for identical HIV sequence matches among viruses recovered in different viral outgrowth cultures and between the sequences of outgrowth viruses and proviral or intracellular HIV RNA sequences in uncultured blood mononuclear cells from eight donors on suppressive ART with diverse proviral populations. All eight donors had viral outgrowth virus that was fully susceptible to their current ART drug regimen. Six of eight donors studied had identical near full-length HIV RNA sequences recovered from different viral outgrowth cultures, and one of the two remaining donors had identical partial viral sequence matches between outgrowth virus and intracellular HIV RNA. These findings provide evidence that clonal expansion of HIV-infected cells is an important mechanism of reservoir persistence that should be targeted to cure HIV infection.

Somatic Genomics and Clinical Features of Lung Adenocarcinoma: A Retrospective Study.

BACKGROUND: Lung adenocarcinoma (LUAD) is the most common histologic subtype of lung cancer and has a high risk of distant metastasis at every disease stage. We aimed to characterize the genomic landscape of LUAD and identify mutation signatures associated with tumor progression. METHODS AND FINDINGS: We performed an integrative genomic analysis, incorporating whole exome sequencing (WES), determination of DNA copy number and DNA methylation, and transcriptome sequencing for 101 LUAD samples from the Environment And Genetics in Lung cancer Etiology (EAGLE) study. We detected driver genes by testing whether the nonsynonymous mutation rate was significantly higher than the background mutation rate and replicated our findings in public datasets with 724 samples. We performed subclonality analysis for mutations based on mutant allele data and copy number alteration data. We also tested the association between mutation signatures and clinical outcomes, including distant metastasis, survival, and tumor grade. We identified and replicated two novel candidate driver genes, POU class 4 homeobox 2 (POU4F2) (mutated in 9 [8.9%] samples) and ZKSCAN1 (mutated in 6 [5.9%] samples), and characterized their major deleterious mutations. ZKSCAN1 was part of a mutually exclusive gene set that included the RTK/RAS/RAF pathway genes BRAF, EGFR, KRAS, MET, and NF1, indicating an important driver role for this gene. Moreover, we observed strong associations between methylation in specific genomic regions and somatic mutation patterns. In the tumor evolution analysis, four driver genes had a significantly lower fraction of subclonal mutations (FSM), including TP53 (p = 0.007), KEAP1 (p = 0.012), STK11 (p = 0.0076), and EGFR (p = 0.0078), suggesting a tumor initiation role for these genes. Subclonal mutations were significantly enriched in APOBEC-related signatures (p < 2.5×10-50). The total number of somatic mutations (p = 0.0039) and the fraction of transitions (p = 5.5×10-4) were associated with increased risk of distant metastasis. Our study's limitations include a small number of LUAD patients for subgroup analyses and a single-sample design for investigation of subclonality. CONCLUSIONS: These data provide a genomic characterization of LUAD pathogenesis and progression. The distinct clonal and subclonal mutation signatures suggest possible diverse carcinogenesis pathways for endogenous and exogenous exposures, and may serve as a foundation for more effective treatments for this lethal disease. LUAD's high heterogeneity emphasizes the need to further study this tumor type and to associate genomic findings with clinical outcomes.

Comparative proteomics of a model MCF10A-KRasG12V cell line reveals a distinct molecular signature of the KRasG12V cell surface.

Oncogenic Ras mutants play a major role in the etiology of most aggressive and deadly carcinomas in humans. In spite of continuous efforts, effective pharmacological treatments targeting oncogenic Ras isoforms have not been developed. Cell-surface proteins represent top therapeutic targets primarily due to their accessibility and susceptibility to different modes of cancer therapy. To expand the treatment options of cancers driven by oncogenic Ras, new targets need to be identified and characterized at the surface of cancer cells expressing oncogenic Ras mutants. Here, we describe a mass spectrometry-based method for molecular profiling of the cell surface using KRasG12V transfected MCF10A (MCF10A-KRasG12V) as a model cell line of constitutively activated KRas and native MCF10A cells transduced with an empty vector (EV) as control. An extensive molecular map of the KRas surface was achieved by applying, in parallel, targeted hydrazide-based cell-surface capturing technology and global shotgun membrane proteomics to identify the proteins on the KRasG12V surface. This method allowed for integrated proteomic analysis that identified more than 500 cell-surface proteins found unique or upregulated on the surface of MCF10A-KRasG12V cells. Multistep bioinformatic processing was employed to elucidate and prioritize targets for cross-validation. Scanning electron microscopy and phenotypic cancer cell assays revealed changes at the cell surface consistent with malignant epithelial-to-mesenchymal transformation secondary to KRasG12V activation. Taken together, this dataset significantly expands the map of the KRasG12V surface and uncovers potential targets involved primarily in cell motility, cellular protrusion formation, and metastasis.

Whole exome sequencing in families at high risk for Hodgkin lymphoma: identification of a predisposing mutation in the KDR gene.

Hodgkin lymphoma shows strong familial aggregation but no major susceptibility genes have been identified to date. The goal of this study was to identify high-penetrance variants using whole exome sequencing in 17 Hodgkin lymphoma prone families with three or more affected cases or obligate carriers (69 individuals), followed by targeted sequencing in an additional 48 smaller HL families (80 individuals). Alignment and variant calling were performed using standard methods. Dominantly segregating, rare, coding or potentially functional variants were further prioritized based on predicted deleteriousness, conservation, and potential importance in lymphoid malignancy pathways. We selected 23 genes for targeted sequencing. Only the p.A1065T variant in KDR (kinase insert domain receptor) also known as VEGFR2 (vascular endothelial growth factor receptor 2) was replicated in two independent Hodgkin lymphoma families. KDR is a type III receptor tyrosine kinase, the main mediator of vascular endothelial growth factor induced proliferation, survival, and migration. Its activity is associated with several diseases including lymphoma. Functional experiments have shown that p.A1065T, located in the activation loop, can promote constitutive autophosphorylation on tyrosine in the absence of vascular endothelial growth factor and that the kinase activity was abrogated after exposure to kinase inhibitors. A few other promising mutations were identified but appear to be "private". In conclusion, in the largest sequenced cohort of Hodgkin lymphoma families to date, we identified a causal mutation in the KDR gene. While independent validation is needed, this mutation may increase downstream tumor cell proliferation activity and might be a candidate for targeted therapy.

Clonally expanded CD4+ T cells can produce infectious HIV-1 in vivo.

Reservoirs of infectious HIV-1 persist despite years of combination antiretroviral therapy and make curing HIV-1 infections a major challenge. Most of the proviral DNA resides in CD4(+)T cells. Some of these CD4(+)T cells are clonally expanded; most of the proviruses are defective. It is not known if any of the clonally expanded cells carry replication-competent proviruses. We report that a highly expanded CD4(+) T-cell clone contains an intact provirus. The highly expanded clone produced infectious virus that was detected as persistent plasma viremia during cART in an HIV-1-infected patient who had squamous cell cancer. Cells containing the intact provirus were widely distributed and significantly enriched in cancer metastases. These results show that clonally expanded CD4(+)T cells can be a reservoir of infectious HIV-1.

Whole exome sequencing reveals a C-terminal germline variant in CEBPA-associated acute myeloid leukemia: 45-year follow up of a large family.

Familial acute myeloid leukemia is rare and linked to germline mutations in RUNX1, GATA2 or CCAAT/enhancer binding protein-α (CEBPA). We re-evaluated a large family with acute myeloid leukemia originally seen at NIH in 1969. We used whole exome sequencing to study this family, and conducted in silico bioinformatics analysis, protein structural modeling and laboratory experiments to assess the impact of the identified CEBPA Q311P mutation. Unlike most previously identified germline mutations in CEBPA, which were N-terminal frameshift mutations, we identified a novel Q311P variant that was located in the C-terminal bZip domain of C/EBPα. Protein structural modeling suggested that the Q311P mutation alters the ability of the CEBPA dimer to bind DNA. Electrophoretic mobility shift assays showed that the Q311P mu-tant had attenuated binding to DNA, as predicted by the protein modeling. Consistent with these findings, we found that the Q311P mutation has reduced transactivation, consistent with a loss-of-function mutation. From 45 years of follow up, we observed incomplete penetrance (46%) of CEBPA Q311P. This study of a large multi-generational pedigree reveals that a germline mutation in the C-terminal bZip domain can alter the ability of C/EBP-α to bind DNA and reduces transactivation, leading to acute myeloid leukemia.

LEDGF/p75 interacts with mRNA splicing factors and targets HIV-1 integration to highly spliced genes.

The host chromatin-binding factor LEDGF/p75 interacts with HIV-1 integrase and directs integration to active transcription units. To understand how LEDGF/p75 recognizes transcription units, we sequenced 1 million HIV-1 integration sites isolated from cultured HEK293T cells. Analysis of integration sites showed that cancer genes were preferentially targeted, raising concerns about using lentivirus vectors for gene therapy. Additional analysis led to the discovery that introns and alternative splicing contributed significantly to integration site selection. These correlations were independent of transcription levels, size of transcription units, and length of the introns. Multivariate analysis with five parameters previously found to predict integration sites showed that intron density is the strongest predictor of integration density in transcription units. Analysis of previously published HIV-1 integration site data showed that integration density in transcription units in mouse embryonic fibroblasts also correlated strongly with intron number, and this correlation was absent in cells lacking LEDGF. Affinity purification showed that LEDGF/p75 is associated with a number of splicing factors, and RNA sequencing (RNA-seq) analysis of HEK293T cells lacking LEDGF/p75 or the LEDGF/p75 integrase-binding domain (IBD) showed that LEDGF/p75 contributes to splicing patterns in half of the transcription units that have alternative isoforms. Thus, LEDGF/p75 interacts with splicing factors, contributes to exon choice, and directs HIV-1 integration to transcription units that are highly spliced.

The somatic autosomal mutation matrix in cancer genomes.

DNA damage in somatic cells originates from both environmental and endogenous sources, giving rise to mutations through multiple mechanisms. When these mutations affect the function of critical genes, cancer may ensue. Although identifying genomic subsets of mutated genes may inform therapeutic options, a systematic survey of tumor mutational spectra is required to improve our understanding of the underlying mechanisms of mutagenesis involved in cancer etiology. Recent studies have presented genome-wide sets of somatic mutations as a 96-element vector, a procedure that only captures the immediate neighbors of the mutated nucleotide. Herein, we present a 32 × 12 mutation matrix that captures the nucleotide pattern two nucleotides upstream and downstream of the mutation. A somatic autosomal mutation matrix (SAMM) was constructed from tumor-specific mutations derived from each of 909 individual cancer genomes harboring a total of 10,681,843 single-base substitutions. In addition, mechanistic template mutation matrices (MTMMs) representing oxidative DNA damage, ultraviolet-induced DNA damage, (5m)CpG deamination, and APOBEC-mediated cytosine mutation, are presented. MTMMs were mapped to the individual tumor SAMMs to determine the maximum contribution of each mutational mechanism to the overall mutation pattern. A Manhattan distance across all SAMM elements between any two tumor genomes was used to determine their relative distance. Employing this metric, 89.5% of all tumor genomes were found to have a nearest neighbor from the same tissue of origin. When a distance-dependent 6-nearest neighbor classifier was used, 10.4% of the SAMMs had an Undetermined tissue of origin, and 92.2% of the remaining SAMMs were assigned to the correct tissue of origin. [corrected]. Thus, although tumors from different tissues may have similar mutation patterns, their SAMMs often display signatures that are characteristic of specific tissues.