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.

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.

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.

A sequence-based variation map of 8.27 million SNPs in inbred mouse strains.

A dense map of genetic variation in the laboratory mouse genome will provide insights into the evolutionary history of the species and lead to an improved understanding of the relationship between inter-strain genotypic and phenotypic differences. Here we resequence the genomes of four wild-derived and eleven classical strains. We identify 8.27 million high-quality single nucleotide polymorphisms (SNPs) densely distributed across the genome, and determine the locations of the high (divergent subspecies ancestry) and low (common subspecies ancestry) SNP-rate intervals for every pairwise combination of classical strains. Using these data, we generate a genome-wide haplotype map containing 40,898 segments, each with an average of three distinct ancestral haplotypes. For the haplotypes in the classical strains that are unequivocally assigned ancestry, the genetic contributions of the Mus musculus subspecies--M. m. domesticus, M. m. musculus, M. m. castaneus and the hybrid M. m. molossinus--are 68%, 6%, 3% and 10%, respectively; the remaining 13% of haplotypes are of unknown ancestral origin. The considerable regional redundancy of the SNP data will facilitate imputation of the majority of these genotypes in less-densely typed classical inbred strains to provide a complete view of variation in additional strains.

A second generation human haplotype map of over 3.1 million SNPs.

We describe the Phase II HapMap, which characterizes over 3.1 million human single nucleotide polymorphisms (SNPs) genotyped in 270 individuals from four geographically diverse populations and includes 25-35% of common SNP variation in the populations surveyed. The map is estimated to capture untyped common variation with an average maximum r2 of between 0.9 and 0.96 depending on population. We demonstrate that the current generation of commercial genome-wide genotyping products captures common Phase II SNPs with an average maximum r2 of up to 0.8 in African and up to 0.95 in non-African populations, and that potential gains in power in association studies can be obtained through imputation. These data also reveal novel aspects of the structure of linkage disequilibrium. We show that 10-30% of pairs of individuals within a population share at least one region of extended genetic identity arising from recent ancestry and that up to 1% of all common variants are untaggable, primarily because they lie within recombination hotspots. We show that recombination rates vary systematically around genes and between genes of different function. Finally, we demonstrate increased differentiation at non-synonymous, compared to synonymous, SNPs, resulting from systematic differences in the strength or efficacy of natural selection between populations.

Genome-wide detection and characterization of positive selection in human populations.

With the advent of dense maps of human genetic variation, it is now possible to detect positive natural selection across the human genome. Here we report an analysis of over 3 million polymorphisms from the International HapMap Project Phase 2 (HapMap2). We used 'long-range haplotype' methods, which were developed to identify alleles segregating in a population that have undergone recent selection, and we also developed new methods that are based on cross-population comparisons to discover alleles that have swept to near-fixation within a population. The analysis reveals more than 300 strong candidate regions. Focusing on the strongest 22 regions, we develop a heuristic for scrutinizing these regions to identify candidate targets of selection. In a complementary analysis, we identify 26 non-synonymous, coding, single nucleotide polymorphisms showing regional evidence of positive selection. Examination of these candidates highlights three cases in which two genes in a common biological process have apparently undergone positive selection in the same population:LARGE and DMD, both related to infection by the Lassa virus, in West Africa;SLC24A5 and SLC45A2, both involved in skin pigmentation, in Europe; and EDAR and EDA2R, both involved in development of hair follicles, in Asia.

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.

Genomewide SNP variation reveals relationships among landraces and modern varieties of rice.

Rice, the primary source of dietary calories for half of humanity, is the first crop plant for which a high-quality reference genome sequence from a single variety was produced. We used resequencing microarrays to interrogate 100 Mb of the unique fraction of the reference genome for 20 diverse varieties and landraces that capture the impressive genotypic and phenotypic diversity of domesticated rice. Here, we report the distribution of 160,000 nonredundant SNPs. Introgression patterns of shared SNPs revealed the breeding history and relationships among the 20 varieties; some introgressed regions are associated with agronomic traits that mark major milestones in rice improvement. These comprehensive SNP data provide a foundation for deep exploration of rice diversity and gene-trait relationships and their use for future rice improvement.

Lung Transplantation and Survival in Children with Cystic Fibrosis.

BACKGROUND: The effects of lung transplantation on the survival and quality of life in children with cystic fibrosis are uncertain. METHODS: We used data from the U.S. Cystic Fibrosis Foundation Patient Registry and from the Organ Procurement and Transplantation Network to identify children with cystic fibrosis who were on the waiting list for lung transplantation during the period from 1992 through 2002. We performed proportional-hazards survival modeling, using multiple clinically relevant covariates that were available before the children were on the waiting list and the interactions of these covariates with lung transplantation as a time-dependent covariate. The data were insufficient in quality and quantity for a retrospective quality-of-life analysis. RESULTS: A total of 248 of the 514 children on the waiting list underwent lung transplantation in the United States during the period from 1992 through 2002. Proportional-hazards modeling identified four variables besides transplantation that were associated with changes in survival. Burkholderia cepacia infection was associated with a trend toward decreased survival, regardless of whether the patient underwent transplantation. A diagnosis of diabetes before the patient was placed on the waiting list decreased survival while the patient was on the waiting list but did not decrease survival after transplantation, whereas older age did not affect waiting-list survival but decreased post-transplantation survival. Staphylococcus aureus infection increased waiting-list survival but decreased post-transplantation survival. Using age, diabetes status, and S. aureus infection status as covariates, we estimated the effect of transplantation on survival for each patient group, expressed as a hazard factor of less than 1 for a benefit and more than 1 for a risk of harm. Five patients had a significant estimated benefit, 283 patients had a significant risk of harm, 102 patients had an insignificant benefit, and 124 patients had an insignificant risk of harm associated with lung transplantation. CONCLUSIONS: Our analyses estimated clearly improved survival for only 5 of 514 patients on the waiting list for lung transplantation. Prolongation of life by means of lung transplantation should not be expected in children with cystic fibrosis. A prospective, randomized trial is needed to clarify whether and when patients derive a survival and quality-of-life benefit from lung transplantation.

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.

A high-density association screen of 155 ion transport genes for involvement with common migraine.

The clinical overlap between monogenic Familial Hemiplegic Migraine (FHM) and common migraine subtypes, and the fact that all three FHM genes are involved in the transport of ions, suggest that ion transport genes may underlie susceptibility to common forms of migraine. To test this leading hypothesis, we examined common variation in 155 ion transport genes using 5257 single nucleotide polymorphisms (SNPs) in a Finnish sample of 841 unrelated migraine with aura cases and 884 unrelated non-migraine controls. The top signals were then tested for replication in four independent migraine case-control samples from the Netherlands, Germany and Australia, totalling 2835 unrelated migraine cases and 2740 unrelated controls. SNPs within 12 genes (KCNB2, KCNQ3, CLIC5, ATP2C2, CACNA1E, CACNB2, KCNE2, KCNK12, KCNK2, KCNS3, SCN5A and SCN9A) with promising nominal association (0.00041 < P < 0.005) in the Finnish sample were selected for replication. Although no variant remained significant after adjusting for multiple testing nor produced consistent evidence for association across all cohorts, a significant epistatic interaction between KCNB2 SNP rs1431656 (chromosome 8q13.3) and CACNB2 SNP rs7076100 (chromosome 10p12.33) (pointwise P = 0.00002; global P = 0.02) was observed in the Finnish case-control sample. We conclude that common variants of moderate effect size in ion transport genes do not play a major role in susceptibility to common migraine within these European populations, although there is some evidence for epistatic interaction between potassium and calcium channel genes, KCNB2 and CACNB2. Multiple rare variants or trans-regulatory elements of these genes are not ruled out.

A genomewide association study of skin pigmentation in a South Asian population.

We have conducted a multistage genomewide association study, using 1,620,742 single-nucleotide polymorphisms to systematically investigate the genetic factors influencing intrinsic skin pigmentation in a population of South Asian descent. Polymorphisms in three genes--SLC24A5, TYR, and SLC45A2--yielded highly significant replicated associations with skin-reflectance measurements, an indirect measure of melanin content in the skin. The associations detected in these three genes, in an additive manner, collectively account for a large fraction of the natural variation of skin pigmentation in a South Asian population. Our study is the first to interrogate polymorphisms across the genome, to find genetic determinants of the natural variation of skin pigmentation within a human population.

Lung transplantation and survival in children with cystic fibrosis.

BACKGROUND: The effects of lung transplantation on the survival and quality of life in children with cystic fibrosis are uncertain. METHODS: We used data from the U.S. Cystic Fibrosis Foundation Patient Registry and from the Organ Procurement and Transplantation Network to identify children with cystic fibrosis who were on the waiting list for lung transplantation during the period from 1992 through 2002. We performed proportional-hazards survival modeling, using multiple clinically relevant covariates that were available before the children were on the waiting list and the interactions of these covariates with lung transplantation as a time-dependent covariate. The data were insufficient in quality and quantity for a retrospective quality-of-life analysis. RESULTS: A total of 248 of the 514 children on the waiting list underwent lung transplantation in the United States during the period from 1992 through 2002. Proportional-hazards modeling identified four variables besides transplantation that were associated with changes in survival. Burkholderia cepacia infection was associated with a trend toward decreased survival, regardless of whether the patient underwent transplantation [corrected]. A diagnosis of diabetes before the patient was placed on the waiting list decreased survival while the patient was on the waiting list but did not decrease survival after transplantation, whereas older age did not affect waiting-list survival but decreased post-transplantation survival. Staphylococcus aureus infection increased waiting-list survival but decreased post-transplantation survival. Using age, diabetes status, and S. aureus infection status as covariates, we estimated the effect of transplantation on survival for each patient group, expressed as a hazard factor of less than 1 for a benefit and more than 1 for a risk of harm. Five patients had a significant estimated benefit, 283 patients had a significant risk of harm, 102 patients had an insignificant benefit, and 124 patients had an insignificant risk of harm associated with lung transplantation [corrected]. CONCLUSIONS: Our analyses estimated clearly improved survival for only 5 of 514 patients on the waiting list for lung transplantation. Prolongation of life by means of lung transplantation should not be expected in children with cystic fibrosis. A prospective, randomized trial is needed to clarify whether and when patients derive a survival and quality-of-life benefit from lung transplantation.

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.

Allele-specific KRT1 expression is a complex trait.

The differential expression of alleles occurs commonly in humans and is likely an important genetic factor underlying heritable differences in phenotypic traits. Understanding the molecular basis of allelic expression differences is thus an important challenge. Although many genes have been shown to display differential allelic expression, this is the first study to examine in detail the cumulative effects of multiple cis-regulatory polymorphisms responsible for allele-specific expression differences. We have used a variety of experimental approaches to identify and characterize cis-regulatory polymorphisms responsible for the extreme allele-specific expression differences of keratin-1 (KRT1) in human white blood cells. The combined data from our analyses provide strong evidence that the KRT1 allelic expression differences result from the haplotypic combinations and interactions of five cis-regulatory single nucleotide polymorphisms (SNPs) whose alleles differ in their affinity to bind transcription factors and modulate KRT1 promoter activity. Two of these cis-regulatory SNPs bind transcriptional activators with the alleles on the high-expressing KRT1 haplotype pattern having a higher affinity than the alleles on the low-expressing haplotype pattern. In contrast, the other three cis-regulatory SNPs bind transcriptional inhibitors with the alleles on the low-expressing haplotype pattern having a higher affinity than the alleles on the high-expressing haplotype pattern. Our study provides important new insights into the degree of complexity that the cis-regulatory sequences responsible for allele-specific transcriptional regulation have. These data suggest that allelic expression differences result from the cumulative contribution of multiple DNA sequence polymorphisms, with each having a small effect, and that allele-specific expression can thus be viewed as a complex trait.