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.

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.

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.

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.

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.

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.

AVIA v2.0: annotation, visualization and impact analysis of genomic variants and genes.

UNLABELLED: As sequencing becomes cheaper and more widely available, there is a greater need to quickly and effectively analyze large-scale genomic data. While the functionality of AVIA v1.0, whose implementation was based on ANNOVAR, was comparable with other annotation web servers, AVIA v2.0 represents an enhanced web-based server that extends genomic annotations to cell-specific transcripts and protein-level functional annotations. With AVIA's improved interface, users can better visualize their data, perform comprehensive searches and categorize both coding and non-coding variants. AVAILABILITY AND IMPLEMENTATION: AVIA is freely available through the web at http://avia.abcc.ncifcrf.gov. CONTACT: Hue.Vuong@fnlcr.nih.gov SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

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.

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.

Guanine holes are prominent targets for mutation in cancer and inherited disease.

Single base substitutions constitute the most frequent type of human gene mutation and are a leading cause of cancer and inherited disease. These alterations occur non-randomly in DNA, being strongly influenced by the local nucleotide sequence context. However, the molecular mechanisms underlying such sequence context-dependent mutagenesis are not fully understood. Using bioinformatics, computational and molecular modeling analyses, we have determined the frequencies of mutation at G • C bp in the context of all 64 5'-NGNN-3' motifs that contain the mutation at the second position. Twenty-four datasets were employed, comprising >530,000 somatic single base substitutions from 21 cancer genomes, >77,000 germline single-base substitutions causing or associated with human inherited disease and 16.7 million benign germline single-nucleotide variants. In several cancer types, the number of mutated motifs correlated both with the free energies of base stacking and the energies required for abstracting an electron from the target guanines (ionization potentials). Similar correlations were also evident for the pathological missense and nonsense germline mutations, but only when the target guanines were located on the non-transcribed DNA strand. Likewise, pathogenic splicing mutations predominantly affected positions in which a purine was located on the non-transcribed DNA strand. Novel candidate driver mutations and tissue-specific mutational patterns were also identified in the cancer datasets. We conclude that electron transfer reactions within the DNA molecule contribute to sequence context-dependent mutagenesis, involving both somatic driver and passenger mutations in cancer, as well as germline alterations causing or associated with inherited disease.

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.

Non-B DB v2.0: a database of predicted non-B DNA-forming motifs and its associated tools.

The non-B DB, available at http://nonb.abcc.ncifcrf.gov, catalogs predicted non-B DNA-forming sequence motifs, including Z-DNA, G-quadruplex, A-phased repeats, inverted repeats, mirror repeats, direct repeats and their corresponding subsets: cruciforms, triplexes and slipped structures, in several genomes. Version 2.0 of the database revises and re-implements the motif discovery algorithms to better align with accepted definitions and thresholds for motifs, expands the non-B DNA-forming motifs coverage by including short tandem repeats and adds key visualization tools to compare motif locations relative to other genomic annotations. Non-B DB v2.0 extends the ability for comparative genomics by including re-annotation of the five organisms reported in non-B DB v1.0, human, chimpanzee, dog, macaque and mouse, and adds seven additional organisms: orangutan, rat, cow, pig, horse, platypus and Arabidopsis thaliana. Additionally, the non-B DB v2.0 provides an overall improved graphical user interface and faster query performance.

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).

Non-B DB: a database of predicted non-B DNA-forming motifs in mammalian genomes.

Although the capability of DNA to form a variety of non-canonical (non-B) structures has long been recognized, the overall significance of these alternate conformations in biology has only recently become accepted en masse. In order to provide access to genome-wide locations of these classes of predicted structures, we have developed non-B DB, a database integrating annotations and analysis of non-B DNA-forming sequence motifs. The database provides the most complete list of alternative DNA structure predictions available, including Z-DNA motifs, quadruplex-forming motifs, inverted repeats, mirror repeats and direct repeats and their associated subsets of cruciforms, triplex and slipped structures, respectively. The database also contains motifs predicted to form static DNA bends, short tandem repeats and homo(purine•pyrimidine) tracts that have been associated with disease. The database has been built using the latest releases of the human, chimp, dog, macaque and mouse genomes, so that the results can be compared directly with other data sources. In order to make the data interpretable in a genomic context, features such as genes, single-nucleotide polymorphisms and repetitive elements (SINE, LINE, etc.) have also been incorporated. The database is accessed through query pages that produce results with links to the UCSC browser and a GBrowse-based genomic viewer. It is freely accessible at http://nonb.abcc.ncifcrf.gov.

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.

Optimized method for computing (18)O/(16)O ratios of differentially stable-isotope labeled peptides in the context of postdigestion (18)O exchange/labeling.

Differential (18)O/(16)O stable isotope labeling of peptides that relies on enzyme-catalyzed oxygen exchange at their carboxyl termini in the presence of H(2)(18)O has been widely used for relative quantitation of peptides/proteins. The role of tryptic proteolysis in bottom-up shotgun proteomics and low reagent costs have made trypsin-catalyzed (18)O postdigestion exchange a convenient and affordable stable isotope labeling approach. However, it is known that trypsin-catalyzed (18)O exchange at the carboxyl terminus is in many instances inhomogeneous/incomplete. The extent of the (18)O exchange/incorporation fluctuates from peptide to peptide mostly due to variable enzyme-substrate affinity. Thus, accurate calculation and interpretation of peptide ratios are analytically complicated and in some regard deficient. Therefore, a computational approach capable of improved measurement of actual (18)O incorporation for each differentially labeled peptide pair is needed. In this regard, we have developed an algorithmic method that relies on the trapezoidal rule to integrate peak intensities of all detected isotopic species across a particular peptide ion over the retention time, which fits the isotopic manifold to Poisson distributions. Optimal values for manifold fitting were calculated and then (18)O/(16)O ratios derived via evolutionary programming. The algorithm is tested using trypsin-catalyzed (18)O postdigestion exchange to differentially label bovine serum albumin (BSA) at a priori determined ratios. Both accuracy and precision are improved utilizing this rigorous mathematical approach. We further demonstrate the effectiveness of this method to accurately calculate (18)O/(16)O ratios in a large scale proteomic quantitation of detergent resistant membrane microdomains (DRMMs) isolated from cells expressing wild-type HIV-1 Gag and its nonmyristylated mutant.

Label-free global serum proteomic profiling reveals novel celecoxib-modulated proteins in familial adenomatous polyposis patients.

Celecoxib, a selective inhibitor of cyclooxygenase-2 (Cox-2), was efficacious in clinical prevention trials of patients with familial adenomatous polyposis (FAP) and sporadic colorectal cancer. To identify as yet poorly defined molecular determinants of celecoxib efficacy, a multidimensional serum fractionation approach was used coupled with nanospray tandem mass spectrometry to perform label-free global proteomic profiling of serum samples from the FAP/celecoxib prevention trial. Subsequently, the application of an algorithm for large-scale biomarker discovery on comparative serum proteomic profiles of pre- and post-celecoxib treatment samples identified 83 potentially celecoxib-responsive proteins from various cellular compartments, biological processes and molecular functions. Celecoxib modulation of some of these proteins was confirmed in serum samples of FAP patients and colorectal cancer cell lines by Western blotting. Thus, using a shotgun procedure to rapidly identify important celecoxib-modulated proteins, this pilot study has uncovered novel systemic changes some of which are highly relevant for carcinogenesis and vascular biology. Validation of selected markers, especially those involved in key signaling networks and those considered molecular indicators of cardiovascular pathology, in larger celecoxib clinical trials is expected to provide insights into the molecular mechanisms of celecoxib and the efficacy/toxicity issues related to its use as a chemopreventive agent.