Genome Capture Sequencing Selectively Enriches Bacterial DNA and Enables Genome-Wide Measurement of Intrastrain Genetic Diversity in Human Infections.

Within-host evolution produces genetic diversity in bacterial strains that cause chronic human infections. However, the lack of facile methods to measure bacterial allelic variation in clinical samples has limited understanding of intrastrain diversity's effects on disease. Here, we report a new method termed genome capture sequencing (GenCap-Seq) in which users inexpensively make hybridization probes from genomic DNA or PCR amplicons to selectively enrich and sequence targeted bacterial DNA from clinical samples containing abundant human or nontarget bacterial DNA. GenCap-Seq enables accurate measurement of allele frequencies over targeted regions and is scalable from specific genes to entire genomes, including the strain-specific accessory genome. The method is effective with samples in which target DNA is rare and inhibitory and DNA-degrading substances are abundant, including human sputum and feces. In proof-of-principle experiments, we used GenCap-Seq to investigate the responses of diversified Pseudomonas aeruginosa populations chronically infecting the lungs of people with cystic fibrosis to in vivo antibiotic exposure, and we found that treatment consistently reduced intrastrain genomic diversity. In addition, analysis of gene-level allele frequency changes suggested that some genes without conventional resistance functions may be important for bacterial fitness during in vivo antibiotic exposure. GenCap-Seq's ability to scalably enrich targeted bacterial DNA from complex samples will enable studies on the effects of intrastrain and intraspecies diversity in human infectious disease. IMPORTANCE Genetic diversity evolves in bacterial strains during human infections and could affect disease manifestations and treatment resistance. However, the extent of diversity present in vivo and its changes over time are difficult to measure by conventional methods. We developed a novel approach, GenCap-Seq, to enrich microbial DNA from complex human samples like sputum and feces for genome-wide measurements of bacterial allelic diversity. The approach is inexpensive, scalable to encompass entire targeted genomes, and works in the presence of abundant untargeted nucleic acids and inhibiting substances. We used GenCap-Seq to investigate in vivo responses of diversified bacterial strains to antibiotic treatment. This method will enable new ideas about the effects of intrastrain diversity on human infections to be tested.

Club cell secretory protein and lung function in children with cystic fibrosis.

BACKGROUND: Club cell secretory protein (CC16) exerts anti-inflammatory functions in lung disease. We sought to determine the relation of serum CC16 deficits and genetic variants that control serum CC16 to lung function among children with cystic fibrosis (CF). METHODS: We used longitudinal data from CF children (EPIC Study) with no positive cultures for Pseudomonas aeruginosa prior to enrollment. Circulating levels of CC16 and an inflammatory score (generated from CRP, SAA, calprotectin, G-CSF) were compared between participants with the lowest and highest FEV1 levels in adolescence (LLF and HLF groups, respectively; N = 130-per-group). Single nucleotide variants (SNVs) in the SCGB1A1, EHF-APIP loci were tested for association with circulating CC16 and with decline of FEV1 and FEV1/FVC% predicted levels between ages 7-16 using mixed models. RESULTS: Compared with the HLF group, the LLF group had lower levels of CC16 (geometric means: 8.2 vs 6.5 ng/ml, respectively; p = 0.0002) and higher levels of the normalized inflammatory score (-0.21 vs 0.21, p = 0.0007). Participants in the lowest CC16 and highest inflammation tertile had the highest odds for having LLF (p<0.0001 for comparison with participants in the highest CC16 and lowest inflammation tertile). Among seven SNVs associated with circulating CC16, the top SNV rs3741240 was associated with decline of FEV1/FVC and, marginally, FEV1 (p = 0.003 and 0.025, respectively; N = 611 participants, 20,801 lung function observations). CONCLUSIONS: Serum CC16 deficits are strongly associated with severity of CF lung disease and their effects are additive with systemic inflammation. The rs3741240 A allele is associated with low circulating CC16 and, possibly, accelerated lung function decline in CF.

An exome-wide study of renal operational tolerance.

Background: Renal operational tolerance is a rare and beneficial state of prolonged renal allograft function in the absence of immunosuppression. The underlying mechanisms are unknown. We hypothesized that tolerance might be driven by inherited protein coding genetic variants with large effect, at least in some patients. Methods: We set up a European survey of over 218,000 renal transplant recipients and collected DNAs from 40 transplant recipients who maintained good allograft function without immunosuppression for at least 1 year. We performed an exome-wide association study comparing the distribution of moderate to high impact variants in 36 tolerant patients, selected for genetic homogeneity using principal component analysis, and 192 controls, using an optimal sequence-kernel association test adjusted for small samples. Results: We identified rare variants of HOMER2 (3/36, FDR 0.0387), IQCH (5/36, FDR 0.0362), and LCN2 (3/36, FDR 0.102) in 10 tolerant patients vs. 0 controls. One patient carried a variant in both HOMER2 and LCN2. Furthermore, the three genes showed an identical variant in two patients each. The three genes are expressed at the primary cilium, a key structure in immune responses. Conclusion: Rare protein coding variants are associated with operational tolerance in a sizable portion of patients. Our findings have important implications for a better understanding of immune tolerance in transplantation and other fields of medicine.ClinicalTrials.gov, identifier: NCT05124444.

Inflammasome Genetic Variants, Macrophage Function, and Clinical Outcomes in Cystic Fibrosis.

Cystic fibrosis (CF) is characterized by chronic airway infection, inflammation, and tissue damage that lead to progressive respiratory failure. NLRP3 and NLRC4 are cytoplasmic pattern recognition receptors that activate the inflammasome, initiating a caspase-1-mediated response. We hypothesized that gain-of-function inflammasome responses are associated with worse outcomes in children with CF. We genotyped nonsynonymous variants in NLRP3 and the NLRC4 pathway from individuals in the EPIC (Early Pseudomonas Infection Control) Observational Study cohort and tested for association with CF outcomes. We generated knockouts of NLRP3 and NLRC4 in human macrophage-like cells and rescued knockouts with wild-type or variant forms of NLRP3 and NLRC4. We identified a SNP in NLRP3, p.(Q705K), that was associated with a higher rate of P. aeruginosa colonization (N = 609; P = 0.01; hazard ratio, 2.3 [Cox model]) and worsened lung function over time as measured by forced expiratory volume in 1 second (N = 445; P = 0.001 [generalized estimating equation]). We identified a SNP in NLRC4, p.(A929S), that was associated with a lower rate of P. aeruginosa colonization as part of a composite of rare variants (N = 405; P = 0.045; hazard ratio, 0.68 [Cox model]) and that was individually associated with protection from lung function decline (P < 0.001 [generalized estimating equation]). Rescue of the NLRP3 knockout with the p.(Q705K) variant produced significantly more IL-1ß in response to NLRP3 stimulation than rescue with the wild type (P = 0.020 [Student's t test]). We identified a subset of children with CF at higher risk of early lung disease progression. Knowledge of these genetic modifiers could guide therapies targeting inflammasome pathways.

Characterization of TP53 mutations in Pap test DNA of women with and without serous ovarian carcinoma.

OBJECTIVE: Pap tests hold promise as a molecular diagnostic for serous ovarian cancer, but previous studies reported limited sensitivity. Furthermore, the presence of somatic mutations in normal tissue is increasingly recognized as a challenge to the specificity of mutation-based cancer diagnostics. We applied an ultra-deep sequencing method with the goal of improving sensitivity and characterizing the landscape of low-frequency somatic TP53 mutations in Pap tests. METHODS: We used CRISPR-DS to deeply sequence (mean Duplex depth ~3000×) the TP53 gene in 30 Pap tests from 21 women without cancer and 9 women with serous ovarian carcinoma with known TP53 driver mutations. Mutations were annotated and compared to those in the TP53 cancer database. RESULTS: The tumor-derived mutation was identified in 3 of 8 Pap tests from women with ovarian cancer and intact tubes. In addition, 221 low-frequency (≲0.001) exonic TP53 mutations were identified in Pap tests from women with ovarian cancer (94 mutations) and without ovarian cancer (127 mutations). Many of these mutations resembled TP53 mutations found in cancer: they impaired protein activity, were predicted to be pathogenic, and clustered in exons 5 to 8 and hotspot codons. Cancer-like mutations were identified in all women but at higher frequency in women with ovarian cancer. CONCLUSIONS: Pap tests have low sensitivity for ovarian cancer detection and carry abundant low-frequency TP53 mutations. These mutations are more frequently pathogenic in women with ovarian cancer. Determining whether low-frequency TP53 mutations in normal gynecologic tissues are associated with an increased cancer risk warrants further study.

Ultra-Sensitive TP53 Sequencing for Cancer Detection Reveals Progressive Clonal Selection in Normal Tissue over a Century of Human Lifespan.

High-accuracy next-generation DNA sequencing promises a paradigm shift in early cancer detection by enabling the identification of mutant cancer molecules in minimally invasive body fluid samples. We demonstrate 80% sensitivity for ovarian cancer detection using ultra-accurate Duplex Sequencing to identify TP53 mutations in uterine lavage. However, in addition to tumor DNA, we also detect low-frequency TP53 mutations in nearly all lavages from women with and without cancer. These mutations increase with age and share the selection traits of clonal TP53 mutations commonly found in human tumors. We show that low-frequency TP53 mutations exist in multiple healthy tissues, from newborn to centenarian, and progressively increase in abundance and pathogenicity with older age across tissue types. Our results illustrate that subclonal cancer evolutionary processes are a ubiquitous part of normal human aging, and great care must be taken to distinguish tumor-derived from age-associated mutations in high-sensitivity clinical cancer diagnostics.

Exonic sequencing identifies TLR1 genetic variation associated with mortality in Thais with melioidosis.

Melioidosis, an infectious disease caused by the bacterium Burkholderia pseudomallei, is a common cause of sepsis in Southeast Asia. We investigated whether novel TLR1 coding variants are associated with outcome in Thai patients with melioidosis. We performed exonic sequencing on a discovery set of patients with extreme phenotypes (mild vs. severe) of bacteremic melioidosis. We analysed the association of missense variants in TLR1 with severe melioidosis in a by-gene analysis. We then genotyped key variants and tested the association with death in two additional sets of melioidosis patients. Using a by-gene analysis, TLR1 was associated with severe bacteremic melioidosis (P = 0.016). One of the eight TLR1 variants identified, rs76600635, a common variant in East Asians, was associated with in-hospital mortality in a replication set of melioidosis patients (adjusted odds ratio 1.71, 95% CI 1.01-2.88, P = 0.04.) In a validation set of patients, the point estimate of effect of the association of rs76600635 with 28-day mortality was similar but not statistically significant (adjusted odds ratio 1.81, 95% CI 0.96-3.44, P = 0.07). Restricting the validation set analysis to patients recruited in a comparable fashion to the discovery and replication sets, rs76600635 was significantly associated with 28-day mortality (adjusted odds ratio 3.88, 95% CI 1.43-10.56, P = 0.01). Exonic sequencing identifies TLR1 as a gene associated with a severe phenotype of bacteremic melioidosis. The TLR1 variant rs76600635, common in East Asian populations, may be associated with poor outcomes from melioidosis. This variant has not been previously associated with outcomes in sepsis and requires further study.

Correction: Whole genome sequencing of extreme phenotypes identifies variants in CD101 and UBE2V1 associated with increased risk of sexually acquired HIV-1.

[This corrects the article DOI: 10.1371/journal.ppat.1006703.].

Mitochondrial DNA Mutations are Associated with Ulcerative Colitis Preneoplasia but Tend to be Negatively Selected in Cancer.

The role of mitochondrial DNA (mtDNA) mutations in cancer remains controversial. Ulcerative colitis is an inflammatory bowel disease that increases the risk of colorectal cancer and involves mitochondrial dysfunction, making it an ideal model to study the role of mtDNA in tumorigenesis. Our goal was to comprehensively characterize mtDNA mutations in ulcerative colitis tumorigenesis using Duplex Sequencing, an ultra-accurate next-generation sequencing method. We analyzed 46 colon biopsies from non-ulcerative colitis control patients and ulcerative colitis patients with and without cancer, including biopsies at all stages of dysplastic progression. mtDNA was sequenced at a median depth of 1,364x. Mutations were classified by mutant allele frequency: clonal > 0.95, subclonal 0.01-0.95, and very low frequency (VLF) < 0.01. We identified 208 clonal and subclonal mutations and 56,764 VLF mutations. Mutations were randomly distributed across the mitochondrial genome. Clonal and subclonal mutations increased in number and pathogenicity in early dysplasia, but decreased in number and pathogenicity in cancer. Most clonal, subclonal, and VLF mutations were C>T transitions in the heavy strand of mtDNA, which likely arise from DNA replication errors. A subset of VLF mutations were C>A transversions, which are probably due to oxidative damage. VLF transitions and indels were less abundant in the non-D-loop region and decreased with progression. Our results indicate that mtDNA mutations are frequent in ulcerative colitis preneoplasia but negatively selected in cancers. IMPLICATIONS: While mtDNA mutations might contribute to early ulcerative colitis tumorigenesis, they appear to be selected against in cancer, suggesting that functional mitochondria might be required for malignant transformation in ulcerative colitis.

Cirrhosis in Hemochromatosis: Independent Risk Factors in 368 HFE p.C282Y Homozygotes.

INTRODUCTION AND AIM: We sought to identify independent risk factors for cirrhosis in HFE p.C282Y homozygotes in a cross-sectional study. MATERIAL AND METHODS: We evaluated 368 p.C282Y homozygotes who underwent liver biopsy and compared characteristics of those with and without cirrhosis. We performed multivariable logistic regression on cirrhosis with: age; sex; race/ethnicity; diabetes; blood pints/units donated voluntarily; erythrocyte pints/units received; iron supplement use; alcohol intake, g/d; body mass index, kg/m2; swollen/tender 2nd/3rd metacarpophalangeal joints; elevated alanine aminotransferase; elevated aspartate aminotransferase; steatosis/fatty liver; iron removed by phlebotomy, g; and GNPAT p.D519G positivity. RESULTS: Mean age of 368 participants (73.6% men) was 47 ± 13 (standard deviation) y. Cirrhosis was diagnosed in 86 participants (23.4%). Participants with cirrhosis had significantly greater mean age, proportion of men, diabetes prevalence, mean daily alcohol intake, prevalence of swollen/ tender 2nd/3rd metacarpophalangeal joints, mean serum ferritin, elevated alanine aminotransferase, elevated aspartate aminotransferase, and mean iron removed; and significantly fewer mean blood pints/units donated. GNPAT p.D519G positivity was detected in 82 of 188 participants (43.6%). In a multivariable model for cirrhosis, there were four significant positive associations: age (10-y intervals) (odds ratio 2.2 [95% confidence interval 1.5, 3.3]); diabetes (3.3; [1.1, 9.7]); alcohol intake (14 g alcohol drinks/d) (1.5 [1.2, 1.8]); and iron removed, g (1.3 [1.2, 1.4]). There was no statistical evidence of two-way interactions between these variables. CONCLUSION: In conclusion, cirrhosis in HFE p.C282Y homozygotes is significantly associated with age, diabetes, daily alcohol intake, and iron removed by phlebotomy, taking into account the effect of other variables.

Variation in Cilia Protein Genes and Progression of Lung Disease in Cystic Fibrosis.

RATIONALE: Cystic fibrosis, like primary ciliary dyskinesia, is an autosomal recessive disorder characterized by abnormal mucociliary clearance and obstructive lung disease. We hypothesized that genes underlying the development or function of cilia may modify lung disease severity in persons with cystic fibrosis. OBJECTIVES: To test this hypothesis, we compared variants in 93 candidate genes in both upper and lower tertiles of lung function in a large cohort of children and adults with cystic fibrosis with those of a population control dataset. METHODS: Variants within candidate genes were tested for association using the SKAT-O test, comparing cystic fibrosis cases defined by poor (n = 127) or preserved (n = 127) lung function with population controls (n = 3,269 or 3,148, respectively). Associated variants were then tested for association with related phenotypes in independent datasets. RESULTS: Variants in DNAH14 and DNAAF3 were associated with poor lung function in cystic fibrosis, whereas variants in DNAH14 and DNAH6 were associated with preserved lung function in cystic fibrosis. Associations between DNAH14 and lung function were replicated in disease-related phenotypes characterized by obstructive lung disease in adults. CONCLUSIONS: Genetic variants within DNAH6, DNAH14, and DNAAF3 are associated with variation in lung function among persons with cystic fibrosis.

Whole genome sequencing of extreme phenotypes identifies variants in CD101 and UBE2V1 associated with increased risk of sexually acquired HIV-1.

Host genetic variation modifying HIV-1 acquisition risk can inform development of HIV-1 prevention strategies. However, associations between rare or intermediate-frequency variants and HIV-1 acquisition are not well studied. We tested for the association between variation in genic regions and extreme HIV-1 acquisition phenotypes in 100 sub-Saharan Africans with whole genome sequencing data. Missense variants in immunoglobulin-like regions of CD101 and, among women, one missense/5' UTR variant in UBE2V1, were associated with increased HIV-1 acquisition risk (p = 1.9x10-4 and p = 3.7x10-3, respectively, for replication). Both of these genes are known to impact host inflammatory pathways. Effect sizes increased with exposure to HIV-1 after adjusting for the independent effect of increasing exposure on acquisition risk. TRIAL REGISTRATION: ClinicalTrials.gov NCT00194519; NCT00557245.

The Epithelial Sodium Channel Is a Modifier of the Long-Term Nonprogressive Phenotype Associated with F508del CFTR Mutations.

Cystic fibrosis (CF) remains the most lethal genetic disease in the Caucasian population. However, there is great variability in clinical phenotypes and survival times, even among patients harboring the same genotype. We identified five patients with CF and a homozygous F508del mutation in the CFTR gene who were in their fifth or sixth decade of life and had shown minimal changes in lung function over a longitudinal period of more than 20 years. Because of the rarity of this long-term nonprogressive phenotype, we hypothesized these individuals may carry rare genetic variants in modifier genes that ameliorate disease severity. Individuals at the extremes of survival time and lung-function trajectory underwent whole-exome sequencing, and the sequencing data were filtered to include rare missense, stopgain, indel, and splicing variants present with a mean allele frequency of <0.2% in general population databases. Epithelial sodium channel (ENaC) mutants were generated via site-directed mutagenesis and expressed for Xenopus oocyte assays. Four of the five individuals carried extremely rare or never reported variants in the SCNN1D and SCNN1B genes of the ENaC. Separately, an independently enriched rare variant in SCNN1D was identified in the Exome Variant Server database associated with a milder pulmonary disease phenotype. Functional analysis using Xenopus oocytes revealed that two of the three variants in δ-ENaC encoded by SCNN1D exhibited hypomorphic channel activity. Our data suggest a potential role for δ-ENaC in controlling sodium reabsorption in the airways, and advance the plausibility of ENaC as a therapeutic target in CF.

Homozygosity for the WRN Helicase-Inactivating Variant, R834C, does not confer a Werner syndrome clinical phenotype.

Loss-of-function mutations in the WRN helicase gene cause Werner syndrome- a progeroid syndrome with an elevated risk of cancer and other age-associated diseases. Large numbers of single nucleotide polymorphisms have been identified in WRN. We report here the organismal, cellular, and molecular phenotypes of variant rs3087425 (c. 2500C > T) that results in an arginine to cysteine substitution at residue 834 (R834C) and up to 90% reduction of WRN helicase activity. This variant is present at a high (5%) frequency in Mexico, where we identified 153 heterozygous and three homozygous individuals among 3,130 genotyped subjects. Family studies of probands identified ten additional TT homozygotes. Biochemical analysis of WRN protein purified from TT lymphoblast cell lines confirmed that the R834C substitution strongly and selectively reduces WRN helicase, but not exonuclease activity. Replication track analyses showed reduced replication fork progression in some homozygous cells following DNA replication stress. Among the thirteen TT homozygotes, we identified a previously unreported and statistically significant gender bias in favor of males (p = 0.0016), but none of the clinical findings associated with Werner syndrome. Our results indicate that WRN helicase activity alone is not rate-limiting for the development of clinical WS.

GNPAT p.D519G is independently associated with markedly increased iron stores in HFE p.C282Y homozygotes.

BACKGROUND: GNPAT p.D519G positivity is significantly increased in HFE p.C282Y homozygotes with markedly increased iron stores. We sought to determine associations of p.D519G and iron-related variables with iron stores in p.C282Y homozygotes. METHODS: We defined markedly increased iron stores as serum ferritin >2247pmol/L (>1000µg/L) and either hepatic iron >236µmol/g dry weight or iron >10g by induction phlebotomy (men and women). We defined normal or mildly elevated iron stores as serum ferritin <674.1pmol/L (<300µg/L) or either age≥40y with iron ≤2.5g iron by induction phlebotomy or age≥50y with ≤3.0g iron by induction phlebotomy (men only). We compared participant subgroups using univariate methods. Using multivariable logistic regression, we evaluated associations of markedly increased iron stores with these variables: age; iron supplement use (dichotomous); whole blood units donated; erythrocyte units received as transfusion; daily alcohol consumption, g; and p.D519G positivity (heterozygosity or homozygosity). RESULTS: The mean age of 56 participants (94.6% men) was 55±10 (SD) y; 41 had markedly increased iron stores. Prevalences of swollen/tender 2nd/3rd metacarpophalangeal joints and elevated aspartate or alanine aminotransferase were significantly greater in participants with markedly increased iron stores. Only participants with markedly increased iron stores had cirrhosis. In multivariable analyses, p.D519G positivity was the only exposure variable significantly associated with markedly increased iron stores (odds ratio 9.9, 95% CI [1.6, 60.3], p=0.0126). CONCLUSIONS: GNPAT p.D519G is strongly associated with markedly increased iron stores in p.C282Y homozygotes after correction for age, iron-related variables, and alcohol consumption.

Ultra-deep sequencing detects ovarian cancer cells in peritoneal fluid and reveals somatic TP53 mutations in noncancerous tissues.

Current sequencing methods are error-prone, which precludes the identification of low frequency mutations for early cancer detection. Duplex sequencing is a sequencing technology that decreases errors by scoring mutations present only in both strands of DNA. Our aim was to determine whether duplex sequencing could detect extremely rare cancer cells present in peritoneal fluid from women with high-grade serous ovarian carcinomas (HGSOCs). These aggressive cancers are typically diagnosed at a late stage and are characterized by TP53 mutations and peritoneal dissemination. We used duplex sequencing to analyze TP53 mutations in 17 peritoneal fluid samples from women with HGSOC and 20 from women without cancer. The tumor TP53 mutation was detected in 94% (16/17) of peritoneal fluid samples from women with HGSOC (frequency as low as 1 mutant per 24,736 normal genomes). Additionally, we detected extremely low frequency TP53 mutations (median mutant fraction 1/13,139) in peritoneal fluid from nearly all patients with and without cancer (35/37). These mutations were mostly deleterious, clustered in hotspots, increased with age, and were more abundant in women with cancer than in controls. The total burden of TP53 mutations in peritoneal fluid distinguished cancers from controls with 82% sensitivity (14/17) and 90% specificity (18/20). Age-associated, low frequency TP53 mutations were also found in 100% of peripheral blood samples from 15 women with and without ovarian cancer (none with hematologic disorder). Our results demonstrate the ability of duplex sequencing to detect rare cancer cells and provide evidence of widespread, low frequency, age-associated somatic TP53 mutation in noncancerous tissue.

Mitochondrial-targeted catalase is good for the old mouse proteome, but not for the young: 'reverse' antagonistic pleiotropy?

Reactive oxygen species (ROS) are highly reactive oxygen-containing molecules associated with aging and a broad spectrum of pathologies. We have previously shown that transgenic expression of the antioxidant enzyme catalase targeted to the mitochondria (mCAT) in mice reduces ROS, attenuates age-related disease, and increases lifespan. However, it has been increasingly recognized that ROS also has beneficial roles in signaling, hormesis, stress response, and immunity. We therefore hypothesized that mCAT might be beneficial only when ROS approaches pathological levels in older age and might not be advantageous at a younger age when basal ROS is low. We analyzed abundance and turnover of the global proteome in hearts and livers of young (4 month) and old (20 month) mCAT and wild-type (WT) mice. In old hearts and livers of WT mice, protein half-lives were reduced compared to young, while in mCAT mice the reverse was observed; the longest half-lives were seen in old mCAT mice and the shortest in young mCAT. Protein abundance of old mCAT hearts recapitulated a more youthful proteomic expression profile (P-value < 0.01). However, young mCAT mice partially phenocopied the older wild-type proteome (P-value < 0.01). Age strongly interacts with mCAT, consistent with antagonistic pleiotropy in the reverse of the typical direction. These findings underscore the contrasting roles of ROS in young vs. old mice and indicate the need for better understanding of the interaction between dose and age in assessing the efficacy of therapeutic interventions in aging, including mitochondrial antioxidants.