Rhesus Cytomegalovirus-encoded Fcγ-binding glycoproteins facilitate viral evasion from IgG-mediated humoral immunity

Human cytomegalovirus (HCMV) encodes four viral Fc-gamma receptors (vFcγRs) that counteract antibody-mediated activation in vitro , but their role in infection and pathogenesis is unknown. To examine the in vivo function of vFcγRs in animal hosts closely related to humans, we identified and characterized vFcγRs encoded by rhesus CMV (RhCMV). We demonstrate that Rh05, Rh152/151 and Rh173 represent the complete set of RhCMV vFcγRs, each displaying functional similarities to their respective HCMV orthologs with respect to antagonizing host FcγR activation in vitro . When RhCMV-naïve rhesus macaques were infected with vFcγR-deleted RhCMV, peak plasma viremia levels and anti-RhCMV antibody responses were comparable to wildtype infections. However, the duration of plasma viremia was significantly shortened in immunocompetent, but not in CD4+ T cell-depleted animals. Since vFcγRs were not required for superinfection, we conclude that vFcγRs delay control by virus-specific adaptive immune responses, particularly antibodies, during primary infection.

Multivalent cytomegalovirus glycoprotein B nucleoside modified mRNA vaccines did not demonstrate a greater antibody breadth

Human cytomegalovirus (HCMV) remains the most common congenital infection and infectious complication in immunocompromised patients. The most successful HCMV vaccine to date, an HCMV glycoprotein B (gB) subunit vaccine adjuvanted with MF59, achieved 50% efficacy against primary HCMV infection. A previous study demonstrated that gB/MF59 vaccinees were less frequently infected with HCMV gB genotype strains most similar to the vaccine strain than strains encoding genetically distinct gB genotypes, suggesting strain-specific immunity accounted for the limited efficacy. To determine whether vaccination with multiple HCMV gB genotypes could increase the breadth of anti-HCMV gB humoral and cellular responses, we immunized 18 female rabbits with monovalent (gB-1), bivalent (gB-1+gB-3), or pentavalent (gB-1+gB-2+gB-3+gB-4+gB-5) gB lipid nanoparticle-encapsulated nucleoside-modified RNA (mRNA-LNP) vaccines. The multivalent vaccine groups did not demonstrate a higher magnitude or breadth of the IgG response to the gB ectodomain or cell-associated gB compared to that of the monovalent vaccine. Also, the multivalent vaccines did not show an increase in the breadth of neutralization activity and antibody-dependent cellular phagocytosis against HCMV strains encoding distinct gB genotypes. Interestingly, peripheral blood mononuclear cell-derived gB-2-specific T-cell responses elicited by multivalent vaccines were of a higher magnitude compared to that of monovalent vaccinated animals against a vaccine-mismatched gB genotype at peak immunogenicity. Yet, no statistical differences were observed in T cell response against gB-3 and gB-5 variable regions among the three vaccine groups. Our data suggests that the inclusion of multivalent gB antigens is not an effective strategy to increase the breadth of anti-HCMV gB antibody and T cell responses. Understanding how to increase the HCMV vaccine protection breadth will be essential to improve the vaccine efficacy.

Multivariate analysis of FcR-mediated NK cell functions identifies unique clustering among humans and rhesus macaques

Rhesus macaques (RMs) are a common pre-clinical model used to test HIV vaccine efficacy and passive immunization strategies. Yet, it remains unclear to what extent the Fc-Fc receptor (FcR) interactions impacting antiviral activities of antibodies in RMs recapitulate those in humans. Here, we evaluated the FcR-related functionality of natural killer cells (NKs) from peripheral blood of uninfected humans and RMs to identify intra- and inter-species variation. NKs were screened for FcγRIIIa (human) and FcγRIII (RM) genotypes (FcγRIII(a)), receptor signaling, and antibody-dependent cellular cytotoxicity (ADCC), the latter mediated by a cocktail of monoclonal IgG1 antibodies with human or RM Fc. FcγRIII(a) genetic polymorphisms alone did not explain differences in NK effector functionality in either species cohort. Using the same parameters, hierarchical clustering separated each species into two clusters. Importantly, in principal components analyses, ADCC magnitude, NK contribution to ADCC, FcγRIII(a) cell-surface expression, and frequency of phosphorylated CD3ζ NK cells all contributed similarly to the first principal component within each species, demonstrating the importance of measuring multiple facets of NK cell function. Although ADCC potency was similar between species, we detected significant differences in frequencies of NK cells and pCD3ζ+ cells, level of cell-surface FcγRIII(a) expression, and NK-mediated ADCC (P<0.001), indicating that a combination of Fc-FcR parameters contribute to overall inter-species functional differences. These data strongly support the importance of multi-parameter analyses of Fc-FcR NK-mediated functions when evaluating efficacy of passive and active immunizations in pre- and clinical trials and identifying correlates of protection. The results also suggest that pre-screening animals for multiple FcR-mediated NK function would ensure even distribution of animals among treatment groups in future preclinical trials.

Immune Phenotype and Postoperative Complications After Elective Surgery

OBJECTIVES

To characterize and quantify accumulating immunologic alterations, pre and postoperatively in patients undergoing elective surgical procedures.

BACKGROUND

Elective surgery is an anticipatable, controlled human injury. Although the human response to injury is generally stereotyped, individual variability exists. This makes surgical outcomes less predictable, even after standardized procedures, and may provoke complications in patients unable to compensate for their injury. One potential source of variation is found in immune cell maturation, with phenotypic changes dependent on an individual's unique, lifelong response to environmental antigens.

METHODS

We enrolled 248 patients in a prospective trial facilitating comprehensive biospecimen and clinical data collection in patients scheduled to undergo elective surgery. Peripheral blood was collected preoperatively, and immediately on return to the postanesthesia care unit. Postoperative complications that occurred within 30 days after surgery were captured.

RESULTS

As this was an elective surgical cohort, outcomes were generally favorable. With a median follow-up of 6 months, the overall survival at 30 days was 100%. However, 20.5% of the cohort experienced a postoperative complication (infection, readmission, or system dysfunction). We identified substantial heterogeneity of immune senescence and terminal differentiation phenotypes in surgical patients. More importantly, phenotypes indicating increased T-cell maturation and senescence were associated with postoperative complications and were evident preoperatively.

CONCLUSIONS

The baseline immune repertoire may define an immune signature of resilience to surgical injury and help predict risk for surgical complications.

Prevalence and significance of clonal hematopoiesis of indeterminate potential in lung transplant recipients

BACKGROUND

Clonal hematopoiesis of indeterminate potential (CHIP), the age-related acquisition of somatic mutations that leads to an expanded blood cell clone, has been associated with development of a pro-inflammatory state. An enhanced or dysregulated inflammatory response may contribute to rejection after lung transplantation, however the prevalence of CHIP in lung recipients and influence of CHIP on allograft outcomes is unknown.

METHODS

We analyzed whole-exome sequencing data in 279 lung recipients to detect CHIP, defined by pre-specified somatic mutations in 74 genes known to promote clonal expansion of hematopoietic stem cells. We compared the burden of acute rejection (AR) over the first post-transplant year in lung recipients with vs. without CHIP using multivariable ordinal regression. Multivariate Cox proportional hazards models were used to assess the association between CHIP and CLAD-free survival. An exploratory analysis evaluated the association between the number of CHIP-associated variants and chronic lung allograft dysfunction (CLAD)-free survival.

RESULTS

We detected 64 CHIP-associated mutations in 45 individuals (15.7%), most commonly in TET2 (10.8%), DNMT3A (9.2%), and U2AF1 (9.2%). Patients with CHIP tended to be older but did not significantly differ from patients without CHIP in terms of race or native lung disease. Patients with CHIP did not have a higher incidence of AR over the first post-transplant year (p = 0.45) or a significantly increased risk of death or CLAD (adjusted HR 1.25, 95% CI 0.88-1.78). We did observe a significant association between the number of CHIP variants and CLAD-free survival, specifically patients with 2 or more CHIP-associated variants had an increased risk for death or CLAD (adjusted HR 3.79, 95% CI 1.98-7.27).

CONCLUSIONS

Lung recipients have a higher prevalence of CHIP and a larger variety of genes with CHIP-associated mutations compared with previous reports for the general population. CHIP did not increase the risk of AR, CLAD, or death in lung recipients.

Relationship of maternal cytomegalovirus-specific antibody responses and viral load to vertical transmission risk following primary maternal infection in a rhesus macaque model

Cytomegalovirus (CMV) is the most common congenital infection and cause of birth defects worldwide. Primary CMV infection during pregnancy leads to a higher frequency of congenital CMV (cCMV) than maternal re-infection, suggesting that maternal immunity confers partial protection. However, poorly understood immune correlates of protection against placental transmission contributes to the current lack of an approved vaccine to prevent cCMV. In this study, we characterized the kinetics of maternal plasma rhesus CMV (RhCMV) viral load (VL) and RhCMV-specific antibody binding and functional responses in a group of 12 immunocompetent dams with acute, primary RhCMV infection. We defined cCMV transmission as RhCMV detection in amniotic fluid (AF) by qPCR. We then leveraged a large group of past and current primary RhCMV infection studies in late-first/early-second trimester RhCMV-seronegative rhesus macaque dams, including immunocompetent (n = 15), CD4+ T cell-depleted with (n = 6) and without (n = 6) RhCMV-specific polyclonal IgG infusion before infection to evaluate differences between RhCMV AF-positive and AF-negative dams. During the first 3 weeks after infection, the magnitude of RhCMV VL in maternal plasma was higher in AF-positive dams in the combined cohort, while RhCMV glycoprotein B (gB)- and pentamer-specific binding IgG responses were lower magnitude compared to AF-negative dams. However, these observed differences were driven by the CD4+ T cell-depleted dams, as there were no differences in plasma VL or antibody responses between immunocompetent AF-positive vs AF-negative dams. Overall, these results suggest that levels of neither maternal plasma viremia nor humoral responses are associated with cCMV following primary maternal infection in healthy individuals. We speculate that other factors related to innate immunity are more important in this context as antibody responses to acute infection likely develop too late to influence vertical transmission. Yet, pre-existing CMV glycoprotein-specific and neutralizing IgG may provide protection against cCMV following primary maternal CMV infection even in high-risk, immunocompromised settings.

The pentameric complex is not required for vertical transmission of cytomegalovirus in seronegative pregnant rhesus macaques

Congenital cytomegalovirus (cCMV) infection is the leading infectious cause of neonatal neurological impairment but essential virological determinants of transplacental CMV transmission remain unclear. The pentameric complex (PC), composed of five subunits, glycoproteins H (gH), gL, UL128, UL130, and UL131A, is essential for efficient entry into non-fibroblast cells in vitro . Based on this role in cell tropism, the PC is considered a possible target for CMV vaccines and immunotherapies to prevent cCMV. To determine the role of the PC in transplacental CMV transmission in a non-human primate model of cCMV, we constructed a PC-deficient rhesus CMV (RhCMV) by deleting the homologues of the HCMV PC subunits UL128 and UL130 and compared congenital transmission to PC-intact RhCMV in CD4+ T cell-depleted or immunocompetent RhCMV-seronegative, pregnant rhesus macaques (RM). Surprisingly, we found that the transplacental transmission rate was similar for PC-intact and PC-deleted RhCMV based on viral genomic DNA detection in amniotic fluid. Moreover, PC-deleted and PC-intact RhCMV acute infection led to similar peak maternal plasma viremia. However, there was less viral shedding in maternal urine and saliva and less viral dissemination in fetal tissues in the PC-deleted group. As expected, dams inoculated with PC-deleted RhCMV demonstrated lower plasma IgG binding to PC-intact RhCMV virions and soluble PC, as well as reduced neutralization of PC-dependent entry of the PC-intact RhCMV isolate UCD52 into epithelial cells. In contrast, binding to gH expressed on the cell surface and neutralization of entry into fibroblasts by the PC-intact RhCMV was higher for dams infected with PC-deleted RhCMV compared to those infected with PC-intact RhCMV. Our data demonstrates that the PC is dispensable for transplacental CMV infection in our non-human primate model.

One Sentence Summary

Congenital CMV transmission frequency in seronegative rhesus macaques is not affected by the deletion of the viral pentameric complex.

Relationship of maternal cytomegalovirus-specific antibody responses and viral load to vertical transmission risk following primary maternal infection in a rhesus macaque model

Cytomegalovirus (CMV) is the most common congenital infection and cause of birth defects worldwide. Primary CMV infection during pregnancy leads to a higher frequency of congenital CMV (cCMV) than maternal re-infection, suggesting that maternal immunity confers partial protection. However, poorly understood immune correlates of protection against placental transmission contributes to the current lack of an approved vaccine to prevent cCMV. In this study, we characterized the kinetics of maternal plasma rhesus CMV (RhCMV) viral load (VL) and RhCMV-specific antibody binding and functional responses in a group of 12 immunocompetent dams with acute, primary RhCMV infection. We defined cCMV transmission as RhCMV detection in amniotic fluid (AF) by qPCR. We then leveraged a large group of past and current primary RhCMV infection studies in late-first/early-second trimester RhCMV-seronegative rhesus macaque dams, including immunocompetent (n=15), CD4+ T cell-depleted with (n=6) and without (n=6) RhCMV-specific polyclonal IgG infusion before infection to evaluate differences between RhCMV AF-positive and AF-negative dams. During the first 3 weeks after infection, the magnitude of RhCMV VL in maternal plasma was higher in AF-positive dams in the combined cohort, while RhCMV glycoprotein B (gB)- and pentamer-specific binding IgG responses were lower magnitude compared to AF-negative dams. However, these observed differences were driven by the CD4+ T cell-depleted dams, as there were no differences in plasma VL or antibody responses between immunocompetent AF-positive vs AF-negative dams. Overall, these results suggest that levels of neither maternal plasma viremia nor humoral responses are associated with cCMV following primary maternal infection in healthy individuals. We speculate that other factors related to innate immunity are more important in this context as antibody responses to acute infection likely develop too late to influence vertical transmission. Yet, pre-existing CMV glycoprotein-specific and neutralizing IgG may provide protection against cCMV following primary maternal CMV infection even in high-risk, immunocompromised settings.

Late gene expression-deficient cytomegalovirus vectors elicit conventional T cells that do not protect against SIV

Rhesus cytomegalovirus-based (RhCMV-based) vaccine vectors induce immune responses that protect ~60% of rhesus macaques (RMs) from SIVmac239 challenge. This efficacy depends on induction of effector memory-based (EM-biased) CD8+ T cells recognizing SIV peptides presented by major histocompatibility complex-E (MHC-E) instead of MHC-Ia. The phenotype, durability, and efficacy of RhCMV/SIV-elicited cellular immune responses were maintained when vector spread was severely reduced by deleting the antihost intrinsic immunity factor phosphoprotein 71 (pp71). Here, we examined the impact of an even more stringent attenuation strategy on vector-induced immune protection against SIV. Fusion of the FK506-binding protein (FKBP) degradation domain to Rh108, the orthologue of the essential human CMV (HCMV) late gene transcription factor UL79, generated RhCMV/SIV vectors that conditionally replicate only when the FK506 analog Shield-1 is present. Despite lacking in vivo dissemination and reduced innate and B cell responses to vaccination, Rh108-deficient 68-1 RhCMV/SIV vectors elicited high-frequency, durable, EM-biased, SIV-specific T cell responses in RhCMV-seropositive RMs at doses of ≥ 1 × 106 PFU. Strikingly, elicited CD8+ T cells exclusively targeted MHC-Ia-restricted epitopes and failed to protect against SIVmac239 challenge. Thus, Rh108-dependent late gene expression is required for both induction of MHC-E-restricted T cells and protection against SIV.

Neuroimaging and immunological features of neurocognitive function related to substance use in people with HIV

This study sought to identify neuroimaging and immunological factors associated with substance use and that contribute to neurocognitive impairment (NCI) in people with HIV (PWH). We performed cross-sectional immunological phenotyping, neuroimaging, and neurocognitive testing on virally suppressed PWH in four substance groups: cocaine only users (COC), marijuana only users (MJ), dual users (Dual), and Non-users. Participants completed substance use assessments, multimodal MRI brain scan, neuropsychological testing, and blood and CSF sampling. We employed a two-stage analysis of 305 possible biomarkers of cognitive function associated with substance use. Feature reduction (Kruskal Wallis p-value < 0.05) identified 53 biomarkers associated with substance use (22 MRI and 31 immunological) for model inclusion along with clinical and demographic variables. We employed eXtreme Gradient Boosting (XGBoost) with these markers to predict cognitive function (global T-score). SHapley Additive exPlanations (SHAP) values were calculated to rank features for impact on model output and NCI. Participants were 110 PWH with sustained HIV viral suppression (33 MJ, 12 COC, 22 Dual, and 43 Non-users). The ten highest ranking biomarkers for predicting global T-score were 4 neuroimaging biomarkers including functional connectivity, gray matter volume, and white matter integrity; 5 soluble biomarkers (plasma glycine, alanine, lyso-phosphatidylcholine (lysoPC) aC17.0, hydroxy-sphingomyelin (SM.OH) C14.1, and phosphatidylcholinediacyl (PC aa) C28.1); and 1 clinical variable (nadir CD4 count). The results of our machine learning model suggest that substance use may indirectly contribute to NCI in PWH through both metabolomic and neuropathological mechanisms.

Characterization of Plasma Immunoglobulin G Responses in Elite Neutralizers of Human Cytomegalovirus

BACKGROUND

Human cytomegalovirus (HCMV) is the most common infectious complication of organ transplantation and cause of birth defects worldwide. There are limited therapeutic options and no licensed vaccine to prevent HCMV infection or disease. To inform development of HCMV antibody-based interventions, a previous study identified individuals with potent and broad plasma HCMV-neutralizing activity, termed elite neutralizers (ENs), from a cohort of HCMV-seropositive (SP) blood donors. However, the specificities and functions of plasma antibodies associated with EN status remained undefined.

METHODS

We sought to determine the plasma antibody specificities, breadth, and Fc-mediated antibody effector functions associated with the most potent HCMV-neutralizing responses in plasma from ENs (n = 25) relative to that from SP donors (n = 19). We measured antibody binding against various HCMV strains and glycoprotein targets and evaluated Fc-mediated effector functions, antibody-dependent cellular cytotoxicity (ADCC), and antibody-dependent cellular phagocytosis (ADCP).

RESULTS

We demonstrate that ENs have elevated immunoglobulin G binding responses against multiple viral glycoproteins, relative to SP donors. Our study also revealed potent HCMV-specific antibody-dependent cellular cytotoxicity and antibody-dependent cellular phagocytosis activity of plasma from ENs.

CONCLUSIONS

We conclude that antibody responses against multiple glycoprotein specificities may be needed to achieve potent plasma neutralization and that potently HCMV elite-neutralizing plasma antibodies can also mediate polyfunctional responses.

Self-assembling peptide nanofiber HIV vaccine elicits robust vaccine-induced antibody functions and modulates Fc glycosylation

To develop vaccines for certain key global pathogens such as HIV, it is crucial to elicit both neutralizing and non-neutralizing Fc-mediated effector antibody functions. Clinical evidence indicates that non-neutralizing antibody functions including antibody-dependent cellular cytotoxicity (ADCC) and antibody-dependent cellular phagocytosis (ADCP) contribute to protection against several pathogens. In this study, we demonstrated that conjugation of HIV Envelope (Env) antigen gp120 to a self-assembling nanofiber material named Q11 induced antibodies with higher breadth and functionality when compared to soluble gp120. Immunization with Q11-conjugated gp120 vaccine (gp120-Q11) demonstrated higher tier 1 neutralization, ADCP, and ADCC as compared to soluble gp120. Moreover, Q11 conjugation altered the Fc N-glycosylation profile of antigen-specific antibodies, leading to a phenotype associated with increased ADCC in animals immunized with gp120-Q11. Thus, this nanomaterial vaccine strategy can enhance non-neutralizing antibody functions possibly through modulation of immunoglobulin G Fc N-glycosylation.

Association between germline variants and somatic mutations in colorectal cancer

Colorectal cancer (CRC) is a heterogeneous disease with evidence of distinct tumor types that develop through different somatically altered pathways. To better understand the impact of the host genome on somatically mutated genes and pathways, we assessed associations of germline variations with somatic events via two complementary approaches. We first analyzed the association between individual germline genetic variants and the presence of non-silent somatic mutations in genes in 1375 CRC cases with genome-wide SNPs data and a tumor sequencing panel targeting 205 genes. In the second analysis, we tested if germline variants located within previously identified regions of somatic allelic imbalance were associated with overall CRC risk using summary statistics from a recent large scale GWAS (n≃125 k CRC cases and controls). The first analysis revealed that a variant (rs78963230) located within a CNA region associated with TLR3 was also associated with a non-silent mutation within gene FBXW7. In the secondary analysis, the variant rs2302274 located in CDX1/PDGFRB frequently gained/lost in colorectal tumors was associated with overall CRC risk (OR = 0.96, p = 7.50e-7). In summary, we demonstrate that an integrative analysis of somatic and germline variation can lead to new insights about CRC.

Genetic Regulation of DNA Methylation Yields Novel Discoveries in GWAS of Colorectal Cancer

BACKGROUND

Colorectal cancer has a strong epigenetic component that is accompanied by frequent DNA methylation (DNAm) alterations in addition to heritable genetic risk. It is of interest to understand the interrelationship of germline genetics, DNAm, and colorectal cancer risk.

METHODS

We performed a genome-wide methylation quantitative trait locus (meQTL) analysis in 1,355 people, assessing the pairwise associations between genetic variants and lymphocytes methylation data. In addition, we used penalized regression with cis-genetic variants ± 1 Mb of methylation to identify genome-wide heritable DNAm. We evaluated the association of genetically predicted methylation with colorectal cancer risk based on genome-wide association studies (GWAS) of over 125,000 cases and controls using the multivariate sMiST as well as univariately via examination of marginal association with colorectal cancer risk.

RESULTS

Of the 142 known colorectal cancer GWAS loci, 47 were identified as meQTLs. We identified four novel colorectal cancer-associated loci (NID2, ATXN10, KLHDC10, and CEP41) that reside over 1 Mb outside of known colorectal cancer loci and 10 secondary signals within 1 Mb of known loci.

CONCLUSIONS

Leveraging information of DNAm regulation into genetic association of colorectal cancer risk reveals novel pathways in colorectal cancer tumorigenesis. Our summary statistics-based framework sMiST provides a powerful approach by combining information from the effect through methylation and residual direct effects of the meQTLs on disease risk. Further validation and functional follow-up of these novel pathways are needed.

IMPACT

Using genotype, DNAm, and GWAS, we identified four new colorectal cancer risk loci. We studied the landscape of genetic regulation of DNAm via single-SNP and multi-SNP meQTL analyses.

HIV envelope antigen valency on peptide nanofibers modulates antibody magnitude and binding breadth

A major challenge in developing an effective vaccine against HIV-1 is the genetic diversity of its viral envelope. Because of the broad range of sequences exhibited by HIV-1 strains, protective antibodies must be able to bind and neutralize a widely mutated viral envelope protein. No vaccine has yet been designed which induces broadly neutralizing or protective immune responses against HIV in humans. Nanomaterial-based vaccines have shown the ability to generate antibody and cellular immune responses of increased breadth and neutralization potency. Thus, we have developed supramolecular nanofiber-based immunogens bearing the HIV gp120 envelope glycoprotein. These immunogens generated antibody responses that had increased magnitude and binding breadth compared to soluble gp120. By varying gp120 density on nanofibers, we determined that increased antigen valency was associated with increased antibody magnitude and germinal center responses. This study presents a proof-of-concept for a nanofiber vaccine platform generating broad, high binding antibody responses against the HIV-1 envelope glycoprotein.

Genetic architectures of proximal and distal colorectal cancer are partly distinct

OBJECTIVE

An understanding of the etiologic heterogeneity of colorectal cancer (CRC) is critical for improving precision prevention, including individualized screening recommendations and the discovery of novel drug targets and repurposable drug candidates for chemoprevention. Known differences in molecular characteristics and environmental risk factors among tumors arising in different locations of the colorectum suggest partly distinct mechanisms of carcinogenesis. The extent to which the contribution of inherited genetic risk factors for CRC differs by anatomical subsite of the primary tumor has not been examined.

DESIGN

To identify new anatomical subsite-specific risk loci, we performed genome-wide association study (GWAS) meta-analyses including data of 48 214 CRC cases and 64 159 controls of European ancestry. We characterised effect heterogeneity at CRC risk loci using multinomial modelling.

RESULTS

We identified 13 loci that reached genome-wide significance (p<5×10-8) and that were not reported by previous GWASs for overall CRC risk. Multiple lines of evidence support candidate genes at several of these loci. We detected substantial heterogeneity between anatomical subsites. Just over half (61) of 109 known and new risk variants showed no evidence for heterogeneity. In contrast, 22 variants showed association with distal CRC (including rectal cancer), but no evidence for association or an attenuated association with proximal CRC. For two loci, there was strong evidence for effects confined to proximal colon cancer.

CONCLUSION

Genetic architectures of proximal and distal CRC are partly distinct. Studies of risk factors and mechanisms of carcinogenesis, and precision prevention strategies should take into consideration the anatomical subsite of the tumour.

Large-Scale Hypothesis Testing for Causal Mediation Effects with Applications in Genome-wide Epigenetic Studies

In genome-wide epigenetic studies, it is of great scientific interest to assess whether the effect of an exposure on a clinical outcome is mediated through DNA methylations. However, statistical inference for causal mediation effects is challenged by the fact that one needs to test a large number of composite null hypotheses across the whole epigenome. Two popular tests, the Wald-type Sobel's test and the joint significant test using the traditional null distribution are underpowered and thus can miss important scientific discoveries. In this paper, we show that the null distribution of Sobel's test is not the standard normal distribution and the null distribution of the joint significant test is not uniform under the composite null of no mediation effect, especially in finite samples and under the singular point null case that the exposure has no effect on the mediator and the mediator has no effect on the outcome. Our results explain why these two tests are underpowered, and more importantly motivate us to develop a more powerful Divide-Aggregate Composite-null Test (DACT) for the composite null hypothesis of no mediation effect by leveraging epigenome-wide data. We adopted Efron's empirical null framework for assessing statistical significance of the DACT test. We showed analytically that the proposed DACT method had improved power, and could well control type I error rate. Our extensive simulation studies showed that, in finite samples, the DACT method properly controlled the type I error rate and outperformed Sobel's test and the joint significance test for detecting mediation effects. We applied the DACT method to the US Department of Veterans Affairs Normative Aging Study, an ongoing prospective cohort study which included men who were aged 21 to 80 years at entry. We identified multiple DNA methylation CpG sites that might mediate the effect of smoking on lung function with effect sizes ranging from -0.18 to -0.79 and false discovery rate controlled at level 0.05, including the CpG sites in the genes AHRR and F2RL3. Our sensitivity analysis found small residual correlations (less than 0.01) of the error terms between the outcome and mediator regressions, suggesting that our results are robust to unmeasured confounding factors.

Reduced plasmablast frequency is associated with seronegative myasthenia gravis

BACKGROUND

The immunopathology of autoimmune seronegative myasthenia gravis (SN MG) is poorly understood. Our objective was to determine immune profiles associated with a diagnosis of SN MG.

METHODS

We performed high-dimensional flow cytometry on blood samples from SN MG patients (N = 68), healthy controls (N = 46), and acetylcholine receptor antibody (AChR+) MG patients (N = 27). We compared 12 immune cell subsets in SN MG to controls using logistic modeling via a discovery-replication design. An exploratory analysis fit a multinomial model comparing AChR+ MG and controls to SN MG.

RESULTS

An increase in CD19⁺ CD20^- CD38^hi plasmablast frequencies was associated with lower odds of being a SN MG case in both the discovery and replication analyses (discovery P-value = .0003, replication P-value = .0021). Interleukin (IL) -21 producing helper T cell frequencies were associated with a diagnosis of AChR+ MG (P = .004).

CONCLUSIONS

Reduced plasmablast frequencies are strongly associated with a SN MG diagnosis and may be a useful diagnostic biomarker in the future.

Intake of Dietary Fruit, Vegetables, and Fiber and Risk of Colorectal Cancer According to Molecular Subtypes: A Pooled Analysis of 9 Studies

Protective associations of fruits, vegetables, and fiber intake with colorectal cancer risk have been shown in many, but not all epidemiologic studies. One possible reason for study heterogeneity is that dietary factors may have distinct effects by colorectal cancer molecular subtypes. Here, we investigate the association of fruit, vegetables, and fiber intake with four well-established colorectal cancer molecular subtypes separately and in combination. Nine observational studies including 9,592 cases with molecular subtypes for microsatellite instability (MSI), CpG island methylator phenotype (CIMP), and somatic mutations in BRAF and KRAS genes, and 7,869 controls were analyzed. Both case-only logistic regression analyses and polytomous logistic regression analyses (with one control set and multiple case groups) were used. Higher fruit intake was associated with a trend toward decreased risk of BRAF-mutated tumors [OR 4th vs. 1st quartile = 0.82 (95% confidence interval, 0.65-1.04)] but not BRAF-wildtype tumors [1.09 (0.97-1.22); P difference as shown in case-only analysis = 0.02]. This difference was observed in case-control studies and not in cohort studies. Compared with controls, higher fiber intake showed negative association with colorectal cancer risk for cases with microsatellite stable/MSI-low, CIMP-negative, BRAF-wildtype, and KRAS-wildtype tumors (P trend range from 0.03 to 3.4e-03), which is consistent with the traditional adenoma-colorectal cancer pathway. These negative associations were stronger compared with MSI-high, CIMP-positive, BRAF-mutated, or KRAS-mutated tumors, but the differences were not statistically significant. These inverse associations for fruit and fiber intake may explain, in part, inconsistent findings between fruit or fiber intake and colorectal cancer risk that have previously been reported. SIGNIFICANCE: These analyses by colorectal cancer molecular subtypes potentially explain the inconsistent findings between dietary fruit or fiber intake and overall colorectal cancer risk that have previously been reported.

A general framework for functionally informed set-based analysis: Application to a large-scale colorectal cancer study

Genome-wide association studies (GWAS) have successfully identified tens of thousands of genetic variants associated with various phenotypes, but together they explain only a fraction of heritability, suggesting many variants have yet to be discovered. Recently it has been recognized that incorporating functional information of genetic variants can improve power for identifying novel loci. For example, S-PrediXcan and TWAS tested the association of predicted gene expression with phenotypes based on GWAS summary statistics by leveraging the information on genetic regulation of gene expression and found many novel loci. However, as genetic variants may have effects on more than one gene and through different mechanisms, these methods likely only capture part of the total effects of these variants. In this paper, we propose a summary statistics-based mixed effects score test (sMiST) that tests for the total effect of both the effect of the mediator by imputing genetically predicted gene expression, like S-PrediXcan and TWAS, and the direct effects of individual variants. It allows for multiple functional annotations and multiple genetically predicted mediators. It can also perform conditional association analysis while adjusting for other genetic variants (e.g., known loci for the phenotype). Extensive simulation and real data analyses demonstrate that sMiST yields p-values that agree well with those obtained from individual level data but with substantively improved computational speed. Importantly, a broad application of sMiST to GWAS is possible, as only summary statistics of genetic variant associations are required. We apply sMiST to a large-scale GWAS of colorectal cancer using summary statistics from ∼120, 000 study participants and gene expression data from the Genotype-Tissue Expression (GTEx) project. We identify several novel and secondary independent genetic loci.

Epigenome-wide association analysis of daytime sleepiness in the Multi-Ethnic Study of Atherosclerosis reveals African-American-specific associations

STUDY OBJECTIVES

Daytime sleepiness is a consequence of inadequate sleep, sleep-wake control disorder, or other medical conditions. Population variability in prevalence of daytime sleepiness is likely due to genetic and biological factors as well as social and environmental influences. DNA methylation (DNAm) potentially influences multiple health outcomes. Here, we explored the association between DNAm and daytime sleepiness quantified by the Epworth Sleepiness Scale (ESS).

METHODS

We performed multi-ethnic and ethnic-specific epigenome-wide association studies for DNAm and ESS in the Multi-Ethnic Study of Atherosclerosis (MESA; n = 619) and the Cardiovascular Health Study (n = 483), with cross-study replication and meta-analysis. Genetic variants near ESS-associated DNAm were analyzed for methylation quantitative trait loci and followed with replication of genotype-sleepiness associations in the UK Biobank.

RESULTS

In MESA only, we detected four DNAm-ESS associations: one across all race/ethnic groups; three in African-Americans (AA) only. Two of the MESA AA associations, in genes KCTD5 and RXRA, nominally replicated in CHS (p-value < 0.05). In the AA meta-analysis, we detected 14 DNAm-ESS associations (FDR q-value < 0.05, top association p-value = 4.26 × 10-8). Three DNAm sites mapped to genes (CPLX3, GFAP, and C7orf50) with biological relevance. We also found evidence for associations with DNAm sites in RAI1, a gene associated with sleep and circadian phenotypes. UK Biobank follow-up analyses detected SNPs in RAI1, RXRA, and CPLX3 with nominal sleepiness associations.

CONCLUSIONS

We identified methylation sites in multiple genes possibly implicated in daytime sleepiness. Most significant DNAm-ESS associations were specific to AA. Future work is needed to identify mechanisms driving ancestry-specific methylation effects.

Associations of variants In the hexokinase 1 and interleukin 18 receptor regions with oxyhemoglobin saturation during sleep

Sleep disordered breathing (SDB)-related overnight hypoxemia is associated with cardiometabolic disease and other comorbidities. Understanding the genetic bases for variations in nocturnal hypoxemia may help understand mechanisms influencing oxygenation and SDB-related mortality. We conducted genome-wide association tests across 10 cohorts and 4 populations to identify genetic variants associated with three correlated measures of overnight oxyhemoglobin saturation: average and minimum oxyhemoglobin saturation during sleep and the percent of sleep with oxyhemoglobin saturation under 90%. The discovery sample consisted of 8,326 individuals. Variants with p < 1 × 10(-6) were analyzed in a replication group of 14,410 individuals. We identified 3 significantly associated regions, including 2 regions in multi-ethnic analyses (2q12, 10q22). SNPs in the 2q12 region associated with minimum SpO2 (rs78136548 p = 2.70 × 10(-10)). SNPs at 10q22 were associated with all three traits including average SpO2 (rs72805692 p = 4.58 × 10(-8)). SNPs in both regions were associated in over 20,000 individuals and are supported by prior associations or functional evidence. Four additional significant regions were detected in secondary sex-stratified and combined discovery and replication analyses, including a region overlapping Reelin, a known marker of respiratory complex neurons.These are the first genome-wide significant findings reported for oxyhemoglobin saturation during sleep, a phenotype of high clinical interest. Our replicated associations with HK1 and IL18R1 suggest that variants in inflammatory pathways, such as the biologically-plausible NLRP3 inflammasome, may contribute to nocturnal hypoxemia.

Transcriptome-wide association studies accounting for colocalization using Egger regression

Integrating genome-wide association (GWAS) and expression quantitative trait locus (eQTL) data into transcriptome-wide association studies (TWAS) based on predicted expression can boost power to detect novel disease loci or pinpoint the susceptibility gene at a known disease locus. However, it is often the case that multiple eQTL genes colocalize at disease loci, making the identification of the true susceptibility gene challenging, due to confounding through linkage disequilibrium (LD). To distinguish between true susceptibility genes (where the genetic effect on phenotype is mediated through expression) and colocalization due to LD, we examine an extension of the Mendelian randomization (MR) egger regression method that allows for LD while only requiring summary association data for both GWAS and eQTL. We derive the standard TWAS approach in the context of MR and show in simulations that the standard TWAS does not control type I error for causal gene identification when eQTLs have pleiotropic or LD-confounded effects on disease. In contrast, LD-aware MR-Egger (LDA MR-Egger) regression can control type I error in this case while attaining similar power as other methods in situations where these provide valid tests. However, when the direct effects of genetic variants on traits are correlated with the eQTL associations, all of the methods we examined including LDA MR-Egger regression can have inflated type I error. We illustrate these methods by integrating gene expression within a recent large-scale breast cancer GWAS to provide guidance on susceptibility gene identification.

Testing for the indirect effect under the null for genome-wide mediation analyses

Mediation analysis helps researchers assess whether part or all of an exposure's effect on an outcome is due to an intermediate variable. The indirect effect can help in designing interventions on the mediator as opposed to the exposure and better understanding the outcome's mechanisms. Mediation analysis has seen increased use in genome-wide epidemiological studies to test for an exposure of interest being mediated through a genomic measure such as gene expression or DNA methylation (DNAm). Testing for the indirect effect is challenged by the fact that the null hypothesis is composite. We examined the performance of commonly used mediation testing methods for the indirect effect in genome-wide mediation studies. When there is no association between the exposure and the mediator and no association between the mediator and the outcome, we show that these common tests are overly conservative. This is a case that will arise frequently in genome-wide mediation studies. Caution is hence needed when applying the commonly used mediation tests in genome-wide mediation studies. We evaluated the performance of these methods using simulation studies, and performed an epigenome-wide mediation association study in the Normative Aging Study, analyzing DNAm as a mediator of the effect of pack-years on FEV₁ .

Reclassification of genetic-based risk predictions as GWAS data accumulate

Background

Disease risk assessments based on common genetic variation have gained widespread attention and use in recent years. The clinical utility of genetic risk profiles depends on the number and effect size of identified loci, and how stable the predicted risks are as additional loci are discovered. Changes in risk classification for individuals over time would undermine the validity of common genetic variation for risk prediction. In this analysis, we quantified reclassification of genetic risk based on past and anticipated future GWAS data.

Methods

We identified disease-associated SNPs via the NHGRI GWAS catalog and recent large scale genome-wide association study (GWAS). We calculated the genomic risk for a simulated cohort of 100,000 individuals based on a multiplicative odds ratio model using cumulative GWAS-identified SNPs at four time points: 2007, 2009, 2011, and 2013. Individuals were classified as Higher Risk (population adjusted odds >2), Average Risk (between 0.5 and 2), and Lower Risk (<0.5) for each time point and we compared classifications between time points for breast cancer (BrCa), prostate cancer (PrCa), diabetes mellitus type 2 (T2D), and cardiovascular heart disease (CHD). We estimated future reclassification using the anticipated number of undiscovered SNPs.

Results

Risk reclassification occurred for all four phenotypes from 2007 to 2013. During the most recent interval (2011-2013), the degree of risk reclassification ranged from 16.3 % for CHD to 24.4 % for PrCa. Many individuals classified as Higher Risk at earlier time points were subsequently reclassified into a lower risk category. From 2011 to 2013, the degree of such downward risk reclassification ranged from 24.9 % for T2D to 55 % for CHD. The percent of individuals classified as Higher Risk increased as more SNPs were discovered, ranging from an increase of 5 % for CHD to 9 % for PrCa from 2007 to 2013. Reclassification continued to occur when we modeled the discovery of anticipated SNPs based on doubling current sample size.

Conclusion

Risk estimates from common genetic variation show large reclassification rates. Identifying disease-associated SNPs facilitates the clinically relevant task of identifying higher-risk individuals. However, the large amount of reclassification that we demonstrated in individuals initially classified as Higher Risk but later as Average Risk or Lower Risk, suggests that caution is currently warranted in basing clinical decisions on common genetic variation for many complex diseases.

Electronic supplementary material

The online version of this article (doi:10.1186/s13073-016-0272-5) contains supplementary material, which is available to authorized users.

Genome-wide association study identifies multiple susceptibility loci for pancreatic cancer

We performed a multistage genome-wide association study including 7,683 individuals with pancreatic cancer and 14,397 controls of European descent. Four new loci reached genome-wide significance: rs6971499 at 7q32.3 (LINC-PINT, per-allele odds ratio (OR) = 0.79, 95% confidence interval (CI) 0.74-0.84, P = 3.0 × 10(-12)), rs7190458 at 16q23.1 (BCAR1/CTRB1/CTRB2, OR = 1.46, 95% CI 1.30-1.65, P = 1.1 × 10(-10)), rs9581943 at 13q12.2 (PDX1, OR = 1.15, 95% CI 1.10-1.20, P = 2.4 × 10(-9)) and rs16986825 at 22q12.1 (ZNRF3, OR = 1.18, 95% CI 1.12-1.25, P = 1.2 × 10(-8)). We identified an independent signal in exon 2 of TERT at the established region 5p15.33 (rs2736098, OR = 0.80, 95% CI 0.76-0.85, P = 9.8 × 10(-14)). We also identified a locus at 8q24.21 (rs1561927, P = 1.3 × 10(-7)) that approached genome-wide significance located 455 kb telomeric of PVT1. Our study identified multiple new susceptibility alleles for pancreatic cancer that are worthy of follow-up studies.

A noninvasive radiotelemetry system to monitor heart rate for assessing stress responses of bovines

A noninvasive radiotelemetry system was developed to monitor heart rates of cows and to view and analyze data. The system was validated by comparing heart rate data of two restrained heifers collected simultaneously using telemetric and direct electrocardiogram measurements and by acquiring data over 72 h from two dry cows housed in an experimental handling facility consisting of a free-stall pen, a holding pen, a pass-through stall, and a second holding pen. Telemetric and direct measurements in response to pharmacological elevation of heart rates were essentially identical. For cows in the experimental facility, peristimulus-time histograms indexed to standing or lying showed that average heart rates for cows increased 4.0 +/- 1.4 beats/min after cows stood and decreased 4.8 +/- 1.0 beats/min after cows lay. Similarly, the average heart rate for the cow naive to the facility increased from 60 to 86 beats/min and remained elevated for 6.3 min when heart rate was indexed to maximal heart rate within +/- 3 min of entry into the pass-through stall. Heart rate for the naive cow increased consistently from around 60 to over 160 beats/min during repeated agonistic encounters between animals. Heart rate for the other cow was not affected by the encounters. These results show clearly that heart rate can be used to monitor animal anxiety.

Carrying as colic "therapy": a randomized controlled trial

In healthy infants, crying behavior is reduced significantly by "supplemental" carrying; that is, increased carrying throughout the day in addition to that which occurs during feeding and in response to crying. To determine whether the recommendation to increase carrying would be effective as a therapy for colic, 66 mothers of infants 4 weeks of age or less who came to their pediatricians with complaints of crying problems ("colic") were randomized to receive standard pediatric advice (standard group) or standard advice plus the recommendation to increase supplemental carrying by 50% (supplemental group). Overall, the supplemental group carried their infants 6.1 hours/d throughout the intervention period, an increase of 2.2 hours/d (56%) more than that provided by the standard group. Despite this significant increase in carrying, there was no difference between groups in the duration or frequency of crying, fussing, or cry/fuss at any time throughout the intervention period. When the greatest treatment effect was expected at 6 weeks, the supplemental group infants cried only 3 minutes less per day (95% confidence interval: 37 minutes less to 32 minutes more per day). We conclude that, compared with standard pediatric advice to be "responsive," supplemental carrying does not reduce crying and fussing behavior further in infants who have colic. In marked contrast to healthy infants, this apparent resistance to increased carrying may indicate an important difference in state regulation and control in infants with colic.

A simple method of drying virus on inanimate objects for virucidal testing

A simple method for drying virus on inanimate objects (cover slips) under vacuum in the cold is described. Following this procedure virus maintains high titers (10(6-7)) for periods of 1-3 wk at -70 degrees C depending on the virus. For virucidal assay of disinfectants, cover slips are exposed to medium simulating the disinfectant (virus control) or disinfectant in an upright position in an Ultra-Vu cuvette. Cover slips are readily removed and placed in tissue culture medium for dilution of virus and determination of virus titer. Cytotoxicity of disinfectant is determined by exposing cover slip without virus to disinfectant, then placing it in medium, diluting the medium and incubating with the indicator cells. The use of this technique results in high titers of virus on cover slips, which are inanimate objects requiring minimal manipulation. The titration of virus or cytotoxicity in microplates is cell, medium, serum, and labware economical.