Structural and genetic diversity in the secreted mucins MUC5AC and MUC5B.

The secreted mucins MUC5AC and MUC5B are large glycoproteins that play critical defensive roles in pathogen entrapment and mucociliary clearance. Their respective genes contain polymorphic and degenerate protein-coding variable number tandem repeats (VNTRs) that make the loci difficult to investigate with short reads. We characterize the structural diversity of MUC5AC and MUC5B by long-read sequencing and assembly of 206 human and 20 nonhuman primate (NHP) haplotypes. We find that human MUC5B is largely invariant (5,761-5,762 amino acids [aa]); however, seven haplotypes have expanded VNTRs (6,291-7,019 aa). In contrast, 30 allelic variants of MUC5AC encode 16 distinct proteins (5,249-6,325 aa) with cysteine-rich domain and VNTR copy-number variation. We group MUC5AC alleles into three phylogenetic clades: H1 (46%, ∼5,654 aa), H2 (33%, ∼5,742 aa), and H3 (7%, ∼6,325 aa). The two most common human MUC5AC variants are smaller than NHP gene models, suggesting a reduction in protein length during recent human evolution. Linkage disequilibrium and Tajima's D analyses reveal that East Asians carry exceptionally large blocks with an excess of rare variation (p < 0.05) at MUC5AC. To validate this result, we use Locityper for genotyping MUC5AC haplogroups in 2,600 unrelated samples from the 1000 Genomes Project. We observe a signature of positive selection in H1 among East Asians and a depletion of the likely ancestral haplogroup (H3). In Europeans, H3 alleles show an excess of common variation and deviate from Hardy-Weinberg equilibrium (p < 0.05), consistent with heterozygote advantage and balancing selection. This study provides a generalizable strategy to characterize complex protein-coding VNTRs for improved disease associations.

MagicalRsq-X: A cross-cohort transferable genotype imputation quality metric.

Since genotype imputation was introduced, researchers have been relying on the estimated imputation quality from imputation software to perform post-imputation quality control (QC). However, this quality estimate (denoted as Rsq) performs less well for lower-frequency variants. We recently published MagicalRsq, a machine-learning-based imputation quality calibration, which leverages additional typed markers from the same cohort and outperforms Rsq as a QC metric. In this work, we extended the original MagicalRsq to allow cross-cohort model training and named the new model MagicalRsq-X. We removed the cohort-specific estimated minor allele frequency and included linkage disequilibrium scores and recombination rates as additional features. Leveraging whole-genome sequencing data from TOPMed, specifically participants in the BioMe, JHS, WHI, and MESA studies, we performed comprehensive cross-cohort evaluations for predominantly European and African ancestral individuals based on their inferred global ancestry with the 1000 Genomes and Human Genome Diversity Project data as reference. Our results suggest MagicalRsq-X outperforms Rsq in almost every setting, with 7.3%-14.4% improvement in squared Pearson correlation with true R2, corresponding to 85-218 K variant gains. We further developed a metric to quantify the genetic distances of a target cohort relative to a reference cohort and showed that such metric largely explained the performance of MagicalRsq-X models. Finally, we found MagicalRsq-X saved up to 53 known genome-wide significant variants in one of the largest blood cell trait GWASs that would be missed using the original Rsq for QC. In conclusion, MagicalRsq-X shows superiority for post-imputation QC and benefits genetic studies by distinguishing well and poorly imputed lower-frequency variants.

Genetic variation in severe cystic fibrosis liver disease is associated with novel mechanisms for disease pathogenesis.

BACKGROUND AND AIMS: It is not known why severe cystic fibrosis (CF) liver disease (CFLD) with portal hypertension occurs in only ~7% of people with CF. We aimed to identify genetic modifiers for severe CFLD to improve understanding of disease mechanisms. APPROACH AND RESULTS: Whole-genome sequencing was available in 4082 people with CF with pancreatic insufficiency (n = 516 with severe CFLD; n = 3566 without CFLD). We tested ~15.9 million single nucleotide polymorphisms (SNPs) for association with severe CFLD versus no-CFLD, using pre-modulator clinical phenotypes including (1) genetic variant ( SERPINA1 ; Z allele) previously associated with severe CFLD; (2) candidate SNPs (n = 205) associated with non-CF liver diseases; (3) genome-wide association study of common/rare SNPs; (4) transcriptome-wide association; and (5) gene-level and pathway analyses. The Z allele was significantly associated with severe CFLD ( p = 1.1 × 10 -4 ). No significant candidate SNPs were identified. A genome-wide association study identified genome-wide significant SNPs in 2 loci and 2 suggestive loci. These 4 loci contained genes [significant, PKD1 ( p = 8.05 × 10 -10 ) and FNBP1 ( p = 4.74 × 10 -9 ); suggestive, DUSP6 ( p = 1.51 × 10 -7 ) and ANKUB1 ( p = 4.69 × 10 -7 )] relevant to severe CFLD pathophysiology. The transcriptome-wide association identified 3 genes [ CXCR1 ( p = 1.01 × 10 -6 ) , AAMP ( p = 1.07 × 10 -6 ), and TRBV24 ( p = 1.23 × 10 -5 )] involved in hepatic inflammation and innate immunity. Gene-ranked analyses identified pathways enriched in genes linked to multiple liver pathologies. CONCLUSION: These results identify loci/genes associated with severe CFLD that point to disease mechanisms involving hepatic fibrosis, inflammation, innate immune function, vascular pathology, intracellular signaling, actin cytoskeleton and tight junction integrity and mechanisms of hepatic steatosis and insulin resistance. These discoveries will facilitate mechanistic studies and the development of therapeutics for severe CFLD.

Structural and genetic diversity in the secreted mucins, MUC5AC and MUC5B.

The secreted mucins MUC5AC and MUC5B play critical defensive roles in airway pathogen entrapment and mucociliary clearance by encoding large glycoproteins with variable number tandem repeats (VNTRs). These polymorphic and degenerate protein coding VNTRs make the loci difficult to investigate with short reads. We characterize the structural diversity of MUC5AC and MUC5B by long-read sequencing and assembly of 206 human and 20 nonhuman primate (NHP) haplotypes. We find that human MUC5B is largely invariant (5761-5762aa); however, seven haplotypes have expanded VNTRs (6291-7019aa). In contrast, 30 allelic variants of MUC5AC encode 16 distinct proteins (5249-6325aa) with cysteine-rich domain and VNTR copy number variation. We grouped MUC5AC alleles into three phylogenetic clades: H1 (46%, ~5654aa), H2 (33%, ~5742aa), and H3 (7%, ~6325aa). The two most common human MUC5AC variants are smaller than NHP gene models, suggesting a reduction in protein length during recent human evolution. Linkage disequilibrium (LD) and Tajima's D analyses reveal that East Asians carry exceptionally large MUC5AC LD blocks with an excess of rare variation (p<0.05). To validate this result, we used Locityper for genotyping MUC5AC haplogroups in 2,600 unrelated samples from the 1000 Genomes Project. We observed signatures of positive selection in H1 and H2 among East Asians and a depletion of the likely ancestral haplogroup (H3). In Africans and Europeans, H3 alleles show an excess of common variation and deviate from Hardy-Weinberg equilibrium, consistent with heterozygote advantage and balancing selection. This study provides a generalizable strategy to characterize complex protein coding VNTRs for improved disease associations.

The clinical utility of sequencing the entirety of CFTR.

BACKGROUND: Cystic fibrosis (CF) is caused by deleterious variants in each CFTR gene. We investigated the utility of whole-gene CFTR sequencing when fewer than two pathogenic or likely pathogenic (P/LP) variants were detected by conventional testing (sequencing of exons and flanking introns) of CFTR. METHODS: Individuals with features of CF and a CF-diagnostic sweat chloride concentration with zero or one P/LP variants identified by conventional testing enrolled in the CF Mutation Analysis Program (MAP) underwent whole-gene CFTR sequencing. Replication was performed on individuals enrolled in the CF Genome Project (CFGP), followed by phenotype review and interrogation of other genes. RESULTS: Whole-gene sequencing identified a second P/LP variant in 20/43 MAP enrollees (47 %) and 10/22 CFGP enrollees (45 %) who had one P/LP variant after conventional testing. No P/LP variants were detected when conventional testing was negative (MAP: n = 43; CFGP: n = 13). Genome-wide analysis was unable to find an alternative etiology in CFGP participants with fewer than two P/LP CFTR variants and CF could not be confirmed in 91 % following phenotype re-review. CONCLUSIONS: Whole-gene CFTR analysis is beneficial in individuals with one previously-identified P/LP variant and a CF-diagnostic sweat chloride. Negative conventional CFTR testing indicates that the phenotype should be re-evaluated.

Prolonged airway explant culture enables study of health, disease, and viral pathogenesis.

In vitro models play a major role in studying airway physiology and disease. However, the native lung's complex tissue architecture and non-epithelial cell lineages are not preserved in these models. Ex vivo tissue models could overcome in vitro limitations, but methods for long-term maintenance of ex vivo tissue has not been established. We describe methods to culture human large airway explants, small airway explants, and precision-cut lung slices for at least 14 days. Human airway explants recapitulate genotype-specific electrophysiology, characteristic epithelial, endothelial, stromal and immune cell populations, and model viral infection after 14 days in culture. These methods also maintain mouse, rabbit, and pig tracheal explants. Notably, intact airway tissue can be cryopreserved, thawed, and used to generate explants with recovery of function 14 days post-thaw. These studies highlight the broad applications of airway tissue explants and their use as translational intermediates between in vitro and in vivo studies.

Dynamic and prognostic proteomic associations with FEV1 decline in chronic obstructive pulmonary disease.

Rationale: Identification and validation of circulating biomarkers for lung function decline in COPD remains an unmet need. Objective: Identify prognostic and dynamic plasma protein biomarkers of COPD progression. Methods: We measured plasma proteins using SomaScan from two COPD-enriched cohorts, the Subpopulations and Intermediate Outcomes Measures in COPD Study (SPIROMICS) and Genetic Epidemiology of COPD (COPDGene), and one population-based cohort, Multi-Ethnic Study of Atherosclerosis (MESA) Lung. Using SPIROMICS as a discovery cohort, linear mixed models identified baseline proteins that predicted future change in FEV1 (prognostic model) and proteins whose expression changed with change in lung function (dynamic model). Findings were replicated in COPDGene and MESA-Lung. Using the COPD-enriched cohorts, Gene Set Enrichment Analysis (GSEA) identified proteins shared between COPDGene and SPIROMICS. Metascape identified significant associated pathways. Measurements and Main Results: The prognostic model found 7 significant proteins in common (p < 0.05) among all 3 cohorts. After applying false discovery rate (adjusted p < 0.2), leptin remained significant in all three cohorts and growth hormone receptor remained significant in the two COPD cohorts. Elevated baseline levels of leptin and growth hormone receptor were associated with slower rate of decline in FEV1. Twelve proteins were nominally but not FDR significant in the dynamic model and all were distinct from the prognostic model. Metascape identified several immune related pathways unique to prognostic and dynamic proteins. Conclusion: We identified leptin as the most reproducible COPD progression biomarker. The difference between prognostic and dynamic proteins suggests disease activity signatures may be different from prognosis signatures.

Design of the SPIROMICS Study of Early COPD Progression: SOURCE Study.

The biological mechanisms leading some tobacco-exposed individuals to develop early-stage chronic obstructive pulmonary disease (COPD) are poorly understood. This knowledge gap hampers development of disease-modifying agents for this prevalent condition. Accord-ingly, with National Heart, Lung and Blood Institute support, we initiated the SPIROMICS Study of Early COPD Progression (SOURCE), a multicenter observational cohort study of younger individuals with a history of cigarette smoking and thus at-risk for, or with, early-stage COPD. Our overall objectives are to identify those who will develop COPD earlier in life, characterize them thoroughly, and by contrasting them to those not developing COPD, define mechanisms of disease progression. SOURCE utilizes the established SPIROMICS clinical network. Its goal is to enroll n=649 participants, ages 30-55 years, all races/ethnicities, with ≥10 pack-years cigarette smoking, in either Global Initiative for Chronic Obstructive Lung Disease (GOLD) groups 0-2 or with Preserved Ratio Impaired Spirometry (PRISm); and an additional n=40 never-smoker controls. Participants undergo baseline and three-year follow-up visits, each including high-resolution computed tomography; respiratory oscillometry and spirometry (pre- and post-bronchodilator administration), exhaled breath condensate (baseline only); and extensive biospecimen collection, including sputum induction. Symptoms, interim healthcare utilization, and exacerbations are captured every six months via follow-up phone calls. An embedded bronchoscopy sub-study involving n=100 participants (including all never-smokers) will allow collection of lower airway samples for genetic, epigenetic, genomic, immunological, microbiome, mucin analyses, and basal cell culture. SOURCE should provide novel insights into the natural history of lung disease in younger individuals with a smoking history, and its biological basis.