Airway inflammation accelerates pulmonary exacerbations in cystic fibrosis.

Airway inflammation underlies cystic fibrosis (CF) pulmonary exacerbations. In a prospective multicenter study of randomly selected, clinically stable adolescents and adults, we assessed relationships between 24 inflammation-associated molecules and the future occurrence of CF pulmonary exacerbation using proportional hazards models. We explored relationships for potential confounding or mediation by clinical factors and assessed sensitivities to treatments including CF transmembrane regulator (CFTR) protein synthesis modulators. Results from 114 participants, including seven on ivacaftor or lumacaftor-ivacaftor, representative of the US CF population during the study period, identified 10 biomarkers associated with future exacerbations mediated by percent predicted forced expiratory volume in 1 s. The findings were not sensitive to anti-inflammatory, antibiotic, and CFTR modulator treatments. The analyses suggest that combination treatments addressing RAGE-axis inflammation, protease-mediated injury, and oxidative stress might prevent pulmonary exacerbations. Our work may apply to other airway inflammatory diseases such as bronchiectasis and the acute respiratory distress syndrome.

SARS-CoV-2 innate effector associations and viral load in early nasopharyngeal infection.

COVID-19 causes severe disease with poor outcomes. We tested the hypothesis that early SARS-CoV-2 viral infection disrupts innate immune responses. These changes may be important for understanding subsequent clinical outcomes. We obtained residual nasopharyngeal swab samples from individuals who requested COVID-19 testing for symptoms at drive-through COVID-19 clinical testing sites operated by the University of Utah. We applied multiplex immunoassays, real-time polymerase chain reaction assays and quantitative proteomics to 20 virus-positive and 20 virus-negative samples. ACE-2 transcripts increased with infection (OR =17.4, 95% CI [CI] =4.78-63.8) and increasing viral N1 protein transcript load (OR =1.16, CI =1.10-1.23). Transcripts for two interferons (IFN) were elevated, IFN-λ1 (OR =71, CI =7.07-713) and IFN-λ2 (OR =40.2, CI =3.86-419), and closely associated with viral N1 transcripts (OR =1.35, CI =1.23-1.49 and OR =1.33 CI =1.20-1.47, respectively). Only transcripts for IP-10 were increased among systemic inflammatory cytokines that we examined (OR =131, CI =1.01-2620). We found widespread discrepancies between transcription and translation. IFN proteins were unchanged or decreased in infected samples (IFN-γ OR =0.90 CI =0.33-0.79, IFN-λ2,3 OR =0.60 CI =0.48-0.74) suggesting viral-induced shut-off of host antiviral protein responses. However, proteins for IP-10 (OR =3.74 CI =2.07-6.77) and several interferon-stimulated genes (ISG) increased with viral load (BST-1 OR =25.1, CI =3.33-188; IFIT1 OR =19.5, CI =4.25-89.2; IFIT3 OR =245, CI =15-4020; MX-1 OR =3.33, CI =1.44-7.70). Older age was associated with substantial modifications of some effects. Ambulatory symptomatic patients had an innate immune response with SARS-CoV-2 infection characterized by elevated IFN, proinflammatory cytokine and ISG transcripts, but there is evidence of a viral-induced host shut-off of antiviral responses. Our findings may characterize the disrupted immune landscape common in patients with early disease.

SARS-CoV-2 Innate Effector Associations and Viral Load in Early Nasopharyngeal Infection.

To examine innate immune responses in early SARS-CoV-2 infection that may change clinical outcomes, we compared nasopharyngeal swab data from 20 virus-positive and 20 virus-negative individuals. Multiple innate immune-related and ACE-2 transcripts increased with infection and were strongly associated with increasing viral load. We found widespread discrepancies between transcription and translation. Interferon proteins were unchanged or decreased in infected samples suggesting virally-induced shut-off of host anti-viral protein responses. However, IP-10 and several interferon-stimulated gene proteins increased with viral load. Older age was associated with modifications of some effects. Our findings may characterize the disrupted immune landscape of early disease.

Prospective multicenter randomized patient recruitment and sample collection to enable future measurements of sputum biomarkers of inflammation in an observational study of cystic fibrosis.

BACKGROUND: Biomarkers of inflammation predictive of cystic fibrosis (CF) disease outcomes would increase the power of clinical trials and contribute to better personalization of clinical assessments. A representative patient cohort would improve searching for believable, generalizable, reproducible and accurate biomarkers. METHODS: We recruited patients from Mountain West CF Consortium (MWCFC) care centers for prospective observational study of sputum biomarkers of inflammation. After informed consent, centers enrolled randomly selected patients with CF who were clinically stable sputum producers, 12 years of age and older, without previous organ transplantation. RESULTS: From December 8, 2014 through January 16, 2016, we enrolled 114 patients (53 male) with CF with continuing data collection. Baseline characteristics included mean age 27 years (SD = 12), 80% predicted forced expiratory volume in 1 s (SD = 23%), 1.0 prior year pulmonary exacerbations (SD = 1.2), home elevation 328 m (SD = 112) above sea level. Compared with other patients in the US CF Foundation Patient Registry (CFFPR) in 2014, MWCFC patients had similar distribution of sex, age, lung function, weight and rates of exacerbations, diabetes, pancreatic insufficiency, CF-related arthropathy and airway infections including methicillin-sensitive or -resistant Staphylococcus aureus, Pseudomonas aeruginosa, Burkholderia cepacia complex, fungal and non-tuberculous Mycobacteria infections. They received CF-specific treatments at similar frequencies. CONCLUSIONS: Randomly-selected, sputum-producing patients within the MWCFC represent sputum-producing patients in the CFFPR. They have similar characteristics, lung function and frequencies of pulmonary exacerbations, microbial infections and use of CF-specific treatments. These findings will plausibly make future interpretations of quantitative measurements of inflammatory biomarkers generalizable to sputum-producing patients in the CFFPR.

Dynamics of Bone Cell Interactions and Differential Responses to PTH and Antibody-Based Therapies.

We propose a mathematical model describing the dynamics of osteoblasts and osteoclasts in bone remodeling. The goal of this work is to develop an integrated modeling framework for bone remodeling and bone cell signaling dynamics that could be used to explore qualitatively combination treatments for osteoporosis in humans. The model has been calibrated using 57 checks from the literature. Specific global optimization methods based on qualitative objectives have been developed to perform the model calibration. We also added pharmacokinetics representations of three drugs to the model, which are teriparatide (PTH(1-34)), denosumab (a RANKL antibody) and romosozumab (a sclerostin antibody), achieving excellent goodness-of-fit of human clinical data. The model reproduces the paradoxical effects of PTH on the bone mass, where continuous administration of PTH results in bone loss but intermittent administration of PTH leads to bone gain, thus proposing an explanation of this phenomenon. We used the model to simulate different categories of osteoporosis. The main attributes of each disease are qualitatively well captured by the model, for example changes in bone turnover in the disease states. We explored dosing regimens for each disease based on the combination of denosumab and romosozumab, identifying adequate ratios and doses of both drugs for subpopulations of patients in function of categories of osteoporosis and the degree of severity of the disease.

Microbial Interactions in the Cystic Fibrosis Airway.

Interactions in the airway ecology of cystic fibrosis may alter organism persistence and clinical outcomes. Better understanding of such interactions could guide clinical decisions. We used generalized estimating equations to fit logistic regression models to longitudinal 2-year patient cohorts in the Cystic Fibrosis Foundation Patient Registry, 2003 to 2011, in order to study associations between the airway organisms present in each calendar year and their presence in the subsequent year. Models were adjusted for clinical characteristics and multiple observations per patient. Adjusted models were tested for sensitivity to cystic fibrosis-specific treatments. The study included 28,042 patients aged 6 years and older from 257 accredited U.S. care centers and affiliates. These patients had produced sputum specimens for at least two consecutive years that were cultured for methicillin-sensitive Staphylococcus aureus, methicillin-resistant S. aureus, Pseudomonas aeruginosa, Burkholderia cepacia complex, Stenotrophomonas maltophilia, Achromobacter xylosoxidans, and Candida and Aspergillus species. We analyzed 99.8% of 538,458 sputum cultures from the patients during the study period. Methicillin-sensitive S. aureus was negatively associated with subsequent Paeruginosa. Paeruginosa was negatively associated with subsequent B. cepacia complex, Axylosoxidans, and Smaltophilia. Bcepacia complex was negatively associated with the future presence of all bacteria studied, as well as with that of Aspergillus species. Paeruginosa, B. cepacia complex, and S. maltophilia were each reciprocally and positively associated with Aspergillus species. Independently of patient characteristics, the organisms studied interact and alter the outcomes of treatment decisions, sometimes in unexpected ways. By inhibiting P. aeruginosa, methicillin-sensitive S. aureus may delay lung disease progression. Paeruginosa and B. cepacia complex may inhibit other organisms by decreasing airway biodiversity, potentially worsening lung disease.

Statistical science: a grammar for research.

I greatly appreciate the invitation to give this lecture with its century long history. The title is a warning that the lecture is rather discursive and not highly focused and technical. The theme is simple. That statistical thinking provides a unifying set of general ideas and specific methods relevant whenever appreciable natural variation is present. To be most fruitful these ideas should merge seamlessly with subject-matter considerations. By contrast, there is sometimes a temptation to regard formal statistical analysis as a ritual to be added after the serious work has been done, a ritual to satisfy convention, referees, and regulatory agencies. I want implicitly to refute that idea.

The design of empirical studies: towards a unified view.

A broad review is given of the general principles underlying study design with emphasis on applications in medical and epidemiological contexts. The main theme of the paper is that, while the distinction between interventionist studies, that is experiments, and purely observational ones is important, there are many common threads. A wide range of specific applications are used in outline to illustrate the discussion.

Effect of genome and environment on metabolic and inflammatory profiles.

Twin and family studies have established the contribution of genetic factors to variation in metabolic, hematologic and immunological parameters. The majority of these studies analyzed single or combined traits into pre-defined syndromes. In the present study, we explore an alternative multivariate approach in which a broad range of metabolic, hematologic, and immunological traits are analyzed simultaneously to determine the resemblance of monozygotic (MZ) twin pairs, twin-spouse pairs and unrelated, non-cohabiting individuals. A total of 517 participants from the Netherlands Twin Register, including 210 MZ twin pairs and 64 twin-spouse pairs, took part in the study. Data were collected on body composition, blood pressure, heart rate, and multiple biomarkers assessed in fasting blood samples, including lipid levels, glucose, insulin, liver enzymes, hematological measurements and cytokine levels. For all 51 measured traits, pair-wise Pearson correlations, correcting for family relatedness, were calculated across all the individuals in the cohort. Hierarchical clustering techniques were applied to group the measured traits into sub-clusters based on similarity. Sub-clusters were observed among metabolic traits and among inflammatory markers. We defined a phenotypic profile as the collection of all the traits measured for a given individual. Average within-pair similarity of phenotypic profiles was determined for the groups of MZ twin pairs, spouse pairs and pairs of unrelated individuals. The average similarity across the full phenotypic profile was higher for MZ twin pairs than for spouse pairs, and lowest for pairs of unrelated individuals. Cohabiting MZ twins were more similar in their phenotypic profile compared to MZ twins who no longer lived together. The correspondence in the phenotypic profile is therefore determined to a large degree by familial, mostly genetic, factors, while household factors contribute to a lesser degree to profile similarity.

Exploring relationships between personality and anatomy performance.

There is increasing recognition in medicine of the importance of noncognitive factors, including personality, for performance, and for good medical practice. The personality domain of conscientiousness is a well-established predictor of performance in workplace and academic settings. This study investigates the relationships between the "Big Five" personality domains, the facets of conscientiousness and performance in a practical anatomy examination. First- and second-year undergraduate medical students (n = 85) completed a paper-based questionnaire, which included a 50-item measure of the Big Five personality domains (neuroticism, extraversion, openness to experience, agreeableness, and conscientiousness) and a 60-item measure of the six conscientiousness facets (orderliness, dutifulness, achievement-striving, self-discipline, self-efficacy, and cautiousness) from the International Personality Item Pool (IPIP). In addition, routinely-collected academic performance scores from the end of semester anatomy practical examinations (spotters) were obtained. Anatomy examination performance correlated moderately with conscientiousness (r = 0.24, P = 0.03). Of the six facets of conscientiousness, a positive relationship was observed between anatomy examination performance and achievement striving (r = 0.22, P = 0.05). In conclusion, this study found that performance in an anatomy examination was related to higher levels of conscientiousness and, more specifically, to higher levels of achievement striving. The results have implications for selection and assessment in medicine.

Aging as accelerated accumulation of somatic variants: whole-genome sequencing of centenarian and middle-aged monozygotic twin pairs.

It has been postulated that aging is the consequence of an accelerated accumulation of somatic DNA mutations and that subsequent errors in the primary structure of proteins ultimately reach levels sufficient to affect organismal functions. The technical limitations of detecting somatic changes and the lack of insight about the minimum level of erroneous proteins to cause an error catastrophe hampered any firm conclusions on these theories. In this study, we sequenced the whole genome of DNA in whole blood of two pairs of monozygotic (MZ) twins, 40 and 100 years old, by two independent next-generation sequencing (NGS) platforms (Illumina and Complete Genomics). Potentially discordant single-base substitutions supported by both platforms were validated extensively by Sanger, Roche 454, and Ion Torrent sequencing. We demonstrate that the genomes of the two twin pairs are germ-line identical between co-twins, and that the genomes of the 100-year-old MZ twins are discerned by eight confirmed somatic single-base substitutions, five of which are within introns. Putative somatic variation between the 40-year-old twins was not confirmed in the validation phase. We conclude from this systematic effort that by using two independent NGS platforms, somatic single nucleotide substitutions can be detected, and that a century of life did not result in a large number of detectable somatic mutations in blood. The low number of somatic variants observed by using two NGS platforms might provide a framework for detecting disease-related somatic variants in phenotypically discordant MZ twins.

Assessing the phenotypic effects in the general population of rare variants in genes for a dominant Mendelian form of diabetes.

Genome sequencing can identify individuals in the general population who harbor rare coding variants in genes for Mendelian disorders and who may consequently have increased disease risk. Previous studies of rare variants in phenotypically extreme individuals display ascertainment bias and may demonstrate inflated effect-size estimates. We sequenced seven genes for maturity-onset diabetes of the young (MODY) in well-phenotyped population samples (n = 4,003). We filtered rare variants according to two prediction criteria for disease-causing mutations: reported previously in MODY or satisfying stringent de novo thresholds (rare, conserved and protein damaging). Approximately 1.5% and 0.5% of randomly selected individuals from the Framingham and Jackson Heart Studies, respectively, carry variants from these two classes. However, the vast majority of carriers remain euglycemic through middle age. Accurate estimates of variant effect sizes from population-based sequencing are needed to avoid falsely predicting a substantial fraction of individuals as being at risk for MODY or other Mendelian diseases.

B cell exchange across the blood-brain barrier in multiple sclerosis.

In multiple sclerosis (MS) pathogenic B cells likely act on both sides of the blood-brain barrier (BBB). However, it is unclear whether antigen-experienced B cells are shared between the CNS and the peripheral blood (PB) compartments. We applied deep repertoire sequencing of IgG heavy chain variable region genes (IgG-VH) in paired cerebrospinal fluid and PB samples from patients with MS and other neurological diseases to identify related B cells that are common to both compartments. For the first time to our knowledge, we found that a restricted pool of clonally related B cells participated in robust bidirectional exchange across the BBB. Some clusters of related IgG-VH appeared to have undergone active diversification primarily in the CNS, while others have undergone active diversification in the periphery or in both compartments in parallel. B cells are strong candidates for autoimmune effector cells in MS, and these findings suggest that CNS-directed autoimmunity may be triggered and supported on both sides of the BBB. These data also provide a powerful approach to identify and monitor B cells in the PB that correspond to clonally amplified populations in the CNS in MS and other inflammatory states.

Sputum biomarkers and the prediction of clinical outcomes in patients with cystic fibrosis.

Lung function, acute pulmonary exacerbations (APE), and weight are the best clinical predictors of survival in cystic fibrosis (CF); however, underlying mechanisms are incompletely understood. Biomarkers of current disease state predictive of future outcomes might identify mechanisms and provide treatment targets, trial endpoints and objective clinical monitoring tools. Such CF-specific biomarkers have previously been elusive. Using observational and validation cohorts comprising 97 non-transplanted consecutively-recruited adult CF patients at the Intermountain Adult CF Center, University of Utah, we identified biomarkers informative of current disease and predictive of future clinical outcomes. Patients represented the majority of sputum producers. They were recruited March 2004-April 2007 and followed through May 2011. Sputum biomarker concentrations were measured and clinical outcomes meticulously recorded for a median 5.9 (interquartile range 5.0 to 6.6) years to study associations between biomarkers and future APE and time-to-lung transplantation or death. After multivariate modeling, only high mobility group box-1 protein (HMGB-1, mean=5.84 [log ng/ml], standard deviation [SD] =1.75) predicted time-to-first APE (hazard ratio [HR] per log-unit HMGB-1=1.56, p-value=0.005), number of future APE within 5 years (0.338 APE per log-unit HMGB-1, p<0.001 by quasi-Poisson regression) and time-to-lung transplantation or death (HR=1.59, p=0.02). At APE onset, sputum granulocyte macrophage colony stimulating factor (GM-CSF, mean 4.8 [log pg/ml], SD=1.26) was significantly associated with APE-associated declines in lung function (-10.8 FEV(1)% points per log-unit GM-CSF, p<0.001 by linear regression). Evaluation of validation cohorts produced similar results that passed tests of mutual consistency. In CF sputum, high HMGB-1 predicts incidence and recurrence of APE and survival, plausibly because it mediates long-term airway inflammation. High APE-associated GM-CSF identifies patients with large acute declines in FEV(1)%, possibly providing a laboratory-based objective decision-support tool for determination of an APE diagnosis. These biomarkers are potential CF reporting tools and treatment targets for slowing long-term progression and reducing short-term severity.

Naive antibody gene-segment frequencies are heritable and unaltered by chronic lymphocyte ablation.

A diverse antibody repertoire is essential for an effective adaptive immune response to novel molecular surfaces. Although past studies have observed common patterns of V-segment use, as well as variation in V-segment use between individuals, the relative contributions to variance from genetics, disease, age, and environment have remained unclear. Using high-throughput sequence analysis of monozygotic twins, we show that variation in naive V(H) and D(H) segment use is strongly determined by an individual's germ-line genetic background. The inherited segment-use profiles are resilient to differential environmental exposure, disease processes, and chronic lymphocyte depletion therapy. Signatures of the inherited profiles were observed in class switched germ-line use of each individual. However, despite heritable segment use, the rearranged complementarity-determining region-H3 repertoires remained highly specific to the individual. As it has been previously demonstrated that certain V-segments exhibit biased representation in autoimmunity, lymphoma, and viral infection, we anticipate our findings may provide a unique mechanism for stratifying individual risk profiles in specific diseases.

Synthetic antibodies designed on natural sequence landscapes.

We present a method for synthetic antibody library generation that combines the use of high-throughput immune repertoire analysis and a novel synthetic technology. The library design recapitulates positional amino acid frequencies observed in natural antibody repertoires. V-segment diversity in four heavy (V(H)) and two kappa (V(κ)) germlines was introduced based on the analysis of somatically hypermutated donor-derived repertoires. Complementarity-determining region 3 length and amino acid designs were based on aggregate frequencies of all V(H) and V(κ) sequences in the data set. The designed libraries were constructed through an adaptation of a novel gene synthesis technology that enables precise positional control of amino acid composition and incorporation frequencies. High-throughput pyrosequencing was used to monitor the fidelity of construction and characterize genetic diversity in the final 3.6×10(10) transformants. The library exhibited Fab expression superior to currently reported synthetic approaches of equivalent diversity, with greater than 93% of clones observed to successfully display both a correctly folded heavy chain and a correctly folded light chain. Genetic diversity in the library was high, with 95% of 7.0×10(5) clones sequenced observed only once. The obtained library diversity explores a comparable sequence space as the donor-derived natural repertoire and, at the same time, is able to access novel recombined diversity due to lack of segmental linkage. The successful isolation of low- and subnanomolar-affinity antibodies against a diverse panel of receptors, growth factors, enzymes, antigens from infectious reagents, and peptides confirms the functional viability of the design strategy.

Genetic variants in the MRPS30 region and postmenopausal breast cancer risk.

BACKGROUND: Genome-wide association studies have identified several genomic regions that are associated with breast cancer risk, but these provide an explanation for only a small fraction of familial breast cancer aggregation. Genotype by environment interactions may contribute further to such explanation, and may help to refine the genomic regions of interest. METHODS: We examined genotypes for 4,988 SNPs, selected from recent genome-wide studies, and four randomized hormonal and dietary interventions among 2,166 women who developed invasive breast cancer during the intervention phase of the Women's Health Initiative (WHI) clinical trial (1993 to 2005), and one-to-one matched controls. These SNPs derive from 3,224 genomic regions having pairwise squared correlation (r2) between adjacent regions less than 0.2. Breast cancer and SNP associations were identified using a test statistic that combined evidence of overall association with evidence for SNPs by intervention interaction. RESULTS: The combined 'main effect' and interaction test led to a focus on two genomic regions, the fibroblast growth factor receptor two (FGFR2) and the mitochondrial ribosomal protein S30 (MRPS30) regions. The ranking of SNPs by significance level, based on this combined test, was rather different from that based on the main effect alone, and drew attention to the vicinities of rs3750817 in FGFR2 and rs7705343 in MRPS30. Specifically, rs7705343 was included with several FGFR2 SNPs in a group of SNPs having an estimated false discovery rate < 0.05. In further analyses, there were suggestions (nominal P < 0.05) that hormonal and dietary intervention hazard ratios varied with the number of minor alleles of rs7705343. CONCLUSIONS: Genotype by environment interaction information may help to define genomic regions relevant to disease risk. Combined main effect and intervention interaction analyses raise novel hypotheses concerning the MRPS30 genomic region and the effects of hormonal and dietary exposures on postmenopausal breast cancer risk.

Immune reconstitution and risk of Kaposi sarcoma and non-Hodgkin lymphoma in HIV-infected adults.

OBJECTIVE: Given the well documented occurrence of immune reconstitution inflammatory syndrome (IRIS) in HIV-infected patients who recently started combination antiretroviral therapy (cART), we examined whether cART initiation increased the risk of Kaposi sarcoma and non-Hodgkin lymphoma (NHL) using data from the Concerted Action on SeroConversion to AIDS and Death in Europe (CASCADE) collaboration. DESIGN: A nested matched case-control study design was used to assess the effects of individual CD4 cell trajectories and exposure to cART close to the time of cancer diagnosis. METHODS: Cases were patients diagnosed with either cancer during follow-up with a minimum of two consecutive CD4 cell readings within the year preceding diagnosis. For each case, up to 10 controls, matched by sex and cohort, were selected by random sampling. Changes in CD4 cell count, calculated by simple and piecewise linear regression, and recent exposure to cART were compared within matched case-control sets using conditional logistic regression. RESULTS: Using data on 689 cases and 4588 controls, we found that an initially low and decreasing CD4 cell count during the year prior to cancer diagnosis is predictive of both Kaposi sarcoma and NHL. Most of this cancer risk is explained by the immunodeficiency characteristic of the period before cART initiation; however, an increased cancer risk was seen in patients who initiated cART in the previous 3 months (odds ratio 2.31; 95% confidence interval 1.33, 4.00). CONCLUSION: Although IRIS may transiently increase the risk of Kaposi sarcoma or NHL in HIV-infected patients, the timely initiation of cART remains the best strategy to avoid the development of these malignancies.