Exome sequencing and characterization of 49,960 individuals in the UK Biobank.

The UK Biobank is a prospective study of 502,543 individuals, combining extensive phenotypic and genotypic data with streamlined access for researchers around the world1. Here we describe the release of exome-sequence data for the first 49,960 study participants, revealing approximately 4 million coding variants (of which around 98.6% have a frequency of less than 1%). The data include 198,269 autosomal predicted loss-of-function (LOF) variants, a more than 14-fold increase compared to the imputed sequence. Nearly all genes (more than 97%) had at least one carrier with a LOF variant, and most genes (more than 69%) had at least ten carriers with a LOF variant. We illustrate the power of characterizing LOF variants in this population through association analyses across 1,730 phenotypes. In addition to replicating established associations, we found novel LOF variants with large effects on disease traits, including PIEZO1 on varicose veins, COL6A1 on corneal resistance, MEPE on bone density, and IQGAP2 and GMPR on blood cell traits. We further demonstrate the value of exome sequencing by surveying the prevalence of pathogenic variants of clinical importance, and show that 2% of this population has a medically actionable variant. Furthermore, we characterize the penetrance of cancer in carriers of pathogenic BRCA1 and BRCA2 variants. Exome sequences from the first 49,960 participants highlight the promise of genome sequencing in large population-based studies and are now accessible to the scientific community.

Exome sequencing of 20,791 cases of type 2 diabetes and 24,440 controls.

Protein-coding genetic variants that strongly affect disease risk can yield relevant clues to disease pathogenesis. Here we report exome-sequencing analyses of 20,791 individuals with type 2 diabetes (T2D) and 24,440 non-diabetic control participants from 5 ancestries. We identify gene-level associations of rare variants (with minor allele frequencies of less than 0.5%) in 4 genes at exome-wide significance, including a series of more than 30 SLC30A8 alleles that conveys protection against T2D, and in 12 gene sets, including those corresponding to T2D drug targets (P = 6.1 × 10-3) and candidate genes from knockout mice (P = 5.2 × 10-3). Within our study, the strongest T2D gene-level signals for rare variants explain at most 25% of the heritability of the strongest common single-variant signals, and the gene-level effect sizes of the rare variants that we observed in established T2D drug targets will require 75,000-185,000 sequenced cases to achieve exome-wide significance. We propose a method to interpret these modest rare-variant associations and to incorporate these associations into future target or gene prioritization efforts.

Pan-ancestry exome-wide association analyses of COVID-19 outcomes in 586,157 individuals.

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) causes coronavirus disease 2019 (COVID-19), a respiratory illness that can result in hospitalization or death. We used exome sequence data to investigate associations between rare genetic variants and seven COVID-19 outcomes in 586,157 individuals, including 20,952 with COVID-19. After accounting for multiple testing, we did not identify any clear associations with rare variants either exome wide or when specifically focusing on (1) 13 interferon pathway genes in which rare deleterious variants have been reported in individuals with severe COVID-19, (2) 281 genes located in susceptibility loci identified by the COVID-19 Host Genetics Initiative, or (3) 32 additional genes of immunologic relevance and/or therapeutic potential. Our analyses indicate there are no significant associations with rare protein-coding variants with detectable effect sizes at our current sample sizes. Analyses will be updated as additional data become available, and results are publicly available through the Regeneron Genetics Center COVID-19 Results Browser.

Skeletal Muscle Adiposity and Lung Function Trajectory in the Severe Asthma Research Program.

Rationale: Extrapulmonary manifestations of asthma, including fatty infiltration in tissues, may reflect systemic inflammation and influence lung function and disease severity. Objectives: To determine if skeletal muscle adiposity predicts lung function trajectory in asthma. Methods: Adult SARP III (Severe Asthma Research Program III) participants with baseline computed tomography imaging and longitudinal postbronchodilator FEV1% predicted (median follow-up 5 years [1,132 person-years]) were evaluated. The mean of left and right paraspinous muscle density (PSMD) at the 12th thoracic vertebral body was calculated (Hounsfield units [HU]). Lower PSMD reflects higher muscle adiposity. We derived PSMD reference ranges from healthy control subjects without asthma. A linear multivariable mixed-effects model was constructed to evaluate associations of baseline PSMD and lung function trajectory stratified by sex. Measurements and Main Results: Participants included 219 with asthma (67% women; mean [SD] body mass index, 32.3 [8.8] kg/m2) and 37 control subjects (51% women; mean [SD] body mass index, 26.3 [4.7] kg/m2). Participants with asthma had lower adjusted PSMD than control subjects (42.2 vs. 55.8 HU; P < 0.001). In adjusted models, PSMD predicted lung function trajectory in women with asthma (ß = -0.47 Δ slope per 10-HU decrease; P = 0.03) but not men (ß = 0.11 Δ slope per 10-HU decrease; P = 0.77). The highest PSMD tertile predicted a 2.9% improvement whereas the lowest tertile predicted a 1.8% decline in FEV1% predicted among women with asthma over 5 years. Conclusions: Participants with asthma have lower PSMD, reflecting greater muscle fat infiltration. Baseline PSMD predicted lung function decline among women with asthma but not men. These data support an important role of metabolic dysfunction in lung function decline.

rECHOmmend: An ECG-Based Machine Learning Approach for Identifying Patients at Increased Risk of Undiagnosed Structural Heart Disease Detectable by Echocardiography.

BACKGROUND: Timely diagnosis of structural heart disease improves patient outcomes, yet many remain underdiagnosed. While population screening with echocardiography is impractical, ECG-based prediction models can help target high-risk patients. We developed a novel ECG-based machine learning approach to predict multiple structural heart conditions, hypothesizing that a composite model would yield higher prevalence and positive predictive values to facilitate meaningful recommendations for echocardiography. METHODS: Using 2 232 130 ECGs linked to electronic health records and echocardiography reports from 484 765 adults between 1984 to 2021, we trained machine learning models to predict the presence or absence of any of 7 echocardiography-confirmed diseases within 1 year. This composite label included the following: moderate or severe valvular disease (aortic/mitral stenosis or regurgitation, tricuspid regurgitation), reduced ejection fraction <50%, or interventricular septal thickness >15 mm. We tested various combinations of input features (demographics, laboratory values, structured ECG data, ECG traces) and evaluated model performance using 5-fold cross-validation, multisite validation trained on 1 site and tested on 10 independent sites, and simulated retrospective deployment trained on pre-2010 data and deployed in 2010. RESULTS: Our composite rECHOmmend model used age, sex, and ECG traces and had a 0.91 area under the receiver operating characteristic curve and a 42% positive predictive value at 90% sensitivity, with a composite label prevalence of 17.9%. Individual disease models had area under the receiver operating characteristic curves from 0.86 to 0.93 and lower positive predictive values from 1% to 31%. Area under the receiver operating characteristic curves for models using different input features ranged from 0.80 to 0.93, increasing with additional features. Multisite validation showed similar results to cross-validation, with an aggregate area under the receiver operating characteristic curve of 0.91 across our independent test set of 10 clinical sites after training on a separate site. Our simulated retrospective deployment showed that for ECGs acquired in patients without preexisting structural heart disease in the year 2010, 11% were classified as high risk and 41% (4.5% of total patients) developed true echocardiography-confirmed disease within 1 year. CONCLUSIONS: An ECG-based machine learning model using a composite end point can identify a high-risk population for having undiagnosed, clinically significant structural heart disease while outperforming single-disease models and improving practical utility with higher positive predictive values. This approach can facilitate targeted screening with echocardiography to improve underdiagnosis of structural heart disease.

CT pectoralis muscle area is associated with DXA lean mass and correlates with emphysema progression in a tobacco-exposed cohort.

INTRODUCTION: Muscle loss is an important extrapulmonary manifestation of COPD. Dual energy X-ray absorptiometry (DXA) is the method of choice for body composition measurement but is not widely used for muscle mass evaluation. The pectoralis muscle area (PMA) is quantifiable by CT and predicts cross-sectional COPD-related morbidity. There are no studies that compare PMA with DXA measures or that evaluate longitudinal relationships between PMA and lung disease progression. METHODS: Participants from our longitudinal tobacco-exposed cohort had baseline and 6-year chest CT (n=259) and DXA (n=164) data. Emphysema was quantified by CT density histogram parenchymal scoring using the 15th percentile technique. Fat-free mass index (FFMI) and appendicular skeletal mass index (ASMI) were calculated from DXA measurements. Linear regression model relationships were reported using standardised coefficient (ß) with 95% CI. RESULTS: PMA was more strongly associated with DXA measures than with body mass index (BMI) in both cross-sectional (FFMI: ß=0.76 (95% CI 0.65 to 0.86), p<0.001; ASMI: ß=0.76 (95% CI 0.66 to 0.86), p<0.001; BMI: ß=0.36 (95% CI 0.25 to 0.47), p<0.001) and longitudinal (ΔFFMI: ß=0.43 (95% CI 0.28 to 0.57), p<0.001; ΔASMI: ß=0.42 (95% CI 0.27 to 0.57), p<0.001; ΔBMI: ß=0.34 (95% CI 0.22 to 0.46), p<0.001) models. Six-year change in PMA was associated with 6-year change in emphysema (ß=0.39 (95% CI 0.23 to 0.56), p<0.001) but not with 6-year change in airflow obstruction. CONCLUSIONS: PMA is an accessible measure of muscle mass and may serve as a useful clinical surrogate for assessing skeletal muscle loss in smokers. Decreased PMA correlated with emphysema progression but not lung function decline, suggesting a difference in the pathophysiology driving emphysema, airflow obstruction and comorbidity risk.

An ECG-based machine learning model for predicting new-onset atrial fibrillation is superior to age and clinical features in identifying patients at high stroke risk.

BACKGROUND: Several large trials have employed age or clinical features to select patients for atrial fibrillation (AF) screening to reduce strokes. We hypothesized that a machine learning (ML) model trained to predict AF risk from 12­lead electrocardiogram (ECG) would be more efficient than criteria based on clinical variables in indicating a population for AF screening to potentially prevent AF-related stroke. METHODS: We retrospectively included all patients with clinical encounters in Geisinger without a prior history of AF. Incidence of AF within 1 year and AF-related strokes within 3 years of the encounter were identified. AF-related stroke was defined as a stroke where AF was diagnosed at the time of stroke or within a year after the stroke. The efficiency of five methods was evaluated for selecting a cohort for AF screening. The methods were selected from four clinical trials (mSToPS, GUARD-AF, SCREEN-AF and STROKESTOP) and the ECG-based ML model. We simulated patient selection for the five methods between the years 2011 and 2014 and evaluated outcomes for 1 year intervals between 2012 and 2015, resulting in a total of twenty 1-year periods. Patients were considered eligible if they met the criteria before the start of the given 1-year period or within that period. The primary outcomes were numbers needed to screen (NNS) for AF and AF-associated stroke. RESULTS: The clinical trial models indicated large proportions of the population with a prior ECG for AF screening (up to 31%), coinciding with NNS ranging from 14 to 18 for AF and 249-359 for AF-associated stroke. At comparable sensitivity, the ECG ML model indicated a modest number of patients for screening (14%) and had the highest efficiency in NNS for AF (7.3; up to 60% reduction) and AF-associated stroke (223; up to 38% reduction). CONCLUSIONS: An ECG-based ML risk prediction model is more efficient than contemporary AF-screening criteria based on age alone or age and clinical features at indicating a population for AF screening to potentially prevent AF-related strokes.

Deep Neural Networks Can Predict New-Onset Atrial Fibrillation From the 12-Lead ECG and Help Identify Those at Risk of Atrial Fibrillation-Related Stroke.

BACKGROUND: Atrial fibrillation (AF) is associated with substantial morbidity, especially when it goes undetected. If new-onset AF could be predicted, targeted screening could be used to find it early. We hypothesized that a deep neural network could predict new-onset AF from the resting 12-lead ECG and that this prediction may help identify those at risk of AF-related stroke. METHODS: We used 1.6 M resting 12-lead digital ECG traces from 430 000 patients collected from 1984 to 2019. Deep neural networks were trained to predict new-onset AF (within 1 year) in patients without a history of AF. Performance was evaluated using areas under the receiver operating characteristic curve and precision-recall curve. We performed an incidence-free survival analysis for a period of 30 years following the ECG stratified by model predictions. To simulate real-world deployment, we trained a separate model using all ECGs before 2010 and evaluated model performance on a test set of ECGs from 2010 through 2014 that were linked to our stroke registry. We identified the patients at risk for AF-related stroke among those predicted to be high risk for AF by the model at different prediction thresholds. RESULTS: The area under the receiver operating characteristic curve and area under the precision-recall curve were 0.85 and 0.22, respectively, for predicting new-onset AF within 1 year of an ECG. The hazard ratio for the predicted high- versus low-risk groups over a 30-year span was 7.2 (95% CI, 6.9-7.6). In a simulated deployment scenario, the model predicted new-onset AF at 1 year with a sensitivity of 69% and specificity of 81%. The number needed to screen to find 1 new case of AF was 9. This model predicted patients at high risk for new-onset AF in 62% of all patients who experienced an AF-related stroke within 3 years of the index ECG. CONCLUSIONS: Deep learning can predict new-onset AF from the 12-lead ECG in patients with no previous history of AF. This prediction may help identify patients at risk for AF-related strokes.

Quantitative CT Characteristics of Cluster Phenotypes in the Severe Asthma Research Program Cohorts.

Background Clustering key clinical characteristics of participants in the Severe Asthma Research Program (SARP), a large, multicenter prospective observational study of patients with asthma and healthy controls, has led to the identification of novel asthma phenotypes. Purpose To determine whether quantitative CT (qCT) could help distinguish between clinical asthma phenotypes. Materials and Methods A retrospective cross-sectional analysis was conducted with the use of qCT images (maximal bronchodilation at total lung capacity [TLC], or inspiration, and functional residual capacity [FRC], or expiration) from the cluster phenotypes of SARP participants (cluster 1: minimal disease; cluster 2: mild, reversible; cluster 3: obese asthma; cluster 4: severe, reversible; cluster 5: severe, irreversible) enrolled between September 2001 and December 2015. Airway morphometry was performed along standard paths (RB1, RB4, RB10, LB1, and LB10). Corresponding voxels from TLC and FRC images were mapped with use of deformable image registration to characterize disease probability maps (DPMs) of functional small airway disease (fSAD), voxel-level volume changes (Jacobian), and isotropy (anisotropic deformation index [ADI]). The association between cluster assignment and qCT measures was evaluated using linear mixed models. Results A total of 455 participants were evaluated with cluster assignments and CT (mean age ± SD, 42.1 years ± 14.7; 270 women). Airway morphometry had limited ability to help discern between clusters. DPM fSAD was highest in cluster 5 (cluster 1 in SARP III: 19.0% ± 20.6; cluster 2: 18.9% ± 13.3; cluster 3: 24.9% ± 13.1; cluster 4: 24.1% ± 8.4; cluster 5: 38.8% ± 14.4; P < .001). Lower whole-lung Jacobian and ADI values were associated with greater cluster severity. Compared to cluster 1, cluster 5 lung expansion was 31% smaller (Jacobian in SARP III cohort: 2.31 ± 0.6 vs 1.61 ± 0.3, respectively, P < .001) and 34% more isotropic (ADI in SARP III cohort: 0.40 ± 0.1 vs 0.61 ± 0.2, P < .001). Within-lung Jacobian and ADI SDs decreased as severity worsened (Jacobian SD in SARP III cohort: 0.90 ± 0.4 for cluster 1; 0.79 ± 0.3 for cluster 2; 0.62 ± 0.2 for cluster 3; 0.63 ± 0.2 for cluster 4; and 0.41 ± 0.2 for cluster 5; P < .001). Conclusion Quantitative CT assessments of the degree and intraindividual regional variability of lung expansion distinguished between well-established clinical phenotypes among participants with asthma from the Severe Asthma Research Program study. © RSNA, 2022 Online supplemental material is available for this article. See also the editorial by Verschakelen in this issue.

Super U-Net: a modularized generalizable architecture.

Objective: To develop and validate a novel convolutional neural network (CNN) termed "Super U-Net" for medical image segmentation. Methods: Super U-Net integrates a dynamic receptive field module and a fusion upsampling module into the classical U-Net architecture. The model was developed and tested to segment retinal vessels, gastrointestinal (GI) polyps, skin lesions on several image types (i.e., fundus images, endoscopic images, dermoscopic images). We also trained and tested the traditional U-Net architecture, seven U-Net variants, and two non-U-Net segmentation architectures. K-fold cross-validation was used to evaluate performance. The performance metrics included Dice similarity coefficient (DSC), accuracy, positive predictive value (PPV), and sensitivity. Results: Super U-Net achieved average DSCs of 0.808±0.0210, 0.752±0.019, 0.804±0.239, and 0.877±0.135 for segmenting retinal vessels, pediatric retinal vessels, GI polyps, and skin lesions, respectively. The Super U-net consistently outperformed U-Net, seven U-Net variants, and two non-U-Net segmentation architectures (p < 0.05). Conclusion: Dynamic receptive fields and fusion upsampling can significantly improve image segmentation performance.

Leveraging population-based exome screening to impact clinical care: The evolution of variant assessment in the Geisinger MyCode research project.

Exome and genome sequencing are increasingly utilized in research studies and clinical care and can provide clinically relevant information beyond the initial intent for sequencing, including medically actionable secondary findings. Despite ongoing debate about sharing this information with patients and participants, a growing number of clinical laboratories and research programs routinely report secondary findings that increase the risk for selected diseases. Recently, there has been a push to maximize the potential benefit of this practice by implementing proactive genomic screening at the population level irrespective of medical history, but the feasibility of deploying population-scale proactive genomic screening requires scaling key elements of the genomic data evaluation process. Herein, we describe the motivation, development, and implementation of a population-scale variant-first screening pipeline combining bioinformatics-based filtering with a manual review process to screen for clinically relevant findings in research exomes generated through the DiscovEHR collaboration within Geisinger's MyCode® research project. Consistent with other studies, this pipeline yields a screen-positive detection rate between 2.1 and 2.6% (depending on inclusion of those with prior indication-based testing) in 130,048 adult MyCode patient-participants screened for clinically relevant findings in 60 genes. Our variant-first pipeline affords cost and time savings by filtering out negative cases, thereby avoiding analysis of each exome one-by-one, as typically employed in the diagnostic setting. While research is still needed to fully appreciate the benefits of population genomic screening, MyCode provides the first demonstration of a program at scale to help shape how population genomic screening is integrated into routine clinical care.

PheWAS and Beyond: The Landscape of Associations with Medical Diagnoses and Clinical Measures across 38,662 Individuals from Geisinger.

Most phenome-wide association studies (PheWASs) to date have used a small to moderate number of SNPs for association with phenotypic data. We performed a large-scale single-cohort PheWAS, using electronic health record (EHR)-derived case-control status for 541 diagnoses using International Classification of Disease version 9 (ICD-9) codes and 25 median clinical laboratory measures. We calculated associations between these diagnoses and traits with ∼630,000 common frequency SNPs with minor allele frequency > 0.01 for 38,662 individuals. In this landscape PheWAS, we explored results within diseases and traits, comparing results to those previously reported in genome-wide association studies (GWASs), as well as previously published PheWASs. We further leveraged the context of functional impact from protein-coding to regulatory regions, providing a deeper interpretation of these associations. The comprehensive nature of this PheWAS allows for novel hypothesis generation, the identification of phenotypes for further study for future phenotypic algorithm development, and identification of cross-phenotype associations.

Profiling and Leveraging Relatedness in a Precision Medicine Cohort of 92,455 Exomes.

Large-scale human genetics studies are ascertaining increasing proportions of populations as they continue growing in both number and scale. As a result, the amount of cryptic relatedness within these study cohorts is growing rapidly and has significant implications on downstream analyses. We demonstrate this growth empirically among the first 92,455 exomes from the DiscovEHR cohort and, via a custom simulation framework we developed called SimProgeny, show that these measures are in line with expectations given the underlying population and ascertainment approach. For example, within DiscovEHR we identified ∼66,000 close (first- and second-degree) relationships, involving 55.6% of study participants. Our simulation results project that >70% of the cohort will be involved in these close relationships, given that DiscovEHR scales to 250,000 recruited individuals. We reconstructed 12,574 pedigrees by using these relationships (including 2,192 nuclear families) and leveraged them for multiple applications. The pedigrees substantially improved the phasing accuracy of 20,947 rare, deleterious compound heterozygous mutations. Reconstructed nuclear families were critical for identifying 3,415 de novo mutations in ∼1,783 genes. Finally, we demonstrate the segregation of known and suspected disease-causing mutations, including a tandem duplication that occurs in LDLR and causes familial hypercholesterolemia, through reconstructed pedigrees. In summary, this work highlights the prevalence of cryptic relatedness expected among large healthcare population-genomic studies and demonstrates several analyses that are uniquely enabled by large amounts of cryptic relatedness.

Elevated plasma level of Pentraxin 3 is associated with emphysema and mortality in smokers.

BACKGROUND: Pentraxin 3 (PTX3) influences innate immunity and inflammation, host defence, the complement cascade and angiogenesis. PTX3 expression in lung and blood of subjects with tobacco exposure, and its potential relationship with disease pattern and clinical outcome are poorly understood. METHODS: Using independent platforms and cohorts, we identified associations of PTX3 gene expression in lung tissue and plasma from current and former tobacco smokers (with and without chronic obstructive pulmonary disease, COPD) to disease phenotypes including quantitative CT determined emphysema, lung function, symptoms and survival. Two putative regulatory variants of the PTX3 gene were examined for association with COPD manifestations. The relationship between plasma PTX3 and hyaluronic acid levels was further examined. RESULTS: PTX3 gene expression in lung tissue was directly correlated with emphysema severity (p<0.0001). Circulating levels of PTX3 were inversely correlated with FEV1 (p=0.006), and positively associated with emphysema severity (p=0.004) and mortality (p=0.008). Two PTX3 gene regulatory variants were associated with a lower risk for emphysema and expiratory airflow obstruction, and plasma levels of PTX3 and hyaluronic acid were related. CONCLUSIONS: These data show strong and overlapping associations of lung and blood PTX3 levels, and PTX3 regulatory gene variants, with the severity of airflow obstruction, emphysema and mortality among smokers. These findings have potential implications regarding the pathogenesis of smoking-related lung diseases and warrant further exploration for the use of PTX3 as a predictive biomarker.

Gene coexpression networks reveal novel molecular endotypes in alpha-1 antitrypsin deficiency.

BACKGROUND: Alpha-1 antitrypsin deficiency (AATD) is a genetic condition that causes early onset pulmonary emphysema and airways obstruction. The complete mechanisms via which AATD causes lung disease are not fully understood. To improve our understanding of the pathogenesis of AATD, we investigated gene expression profiles of bronchoalveolar lavage (BAL) and peripheral blood mononuclear cells (PBMCs) in AATD individuals. METHODS: We performed RNA-Seq on RNA extracted from matched BAL and PBMC samples isolated from 89 subjects enrolled in the Genomic Research in Alpha-1 Antitrypsin Deficiency and Sarcoidosis (GRADS) study. Subjects were stratified by genotype and augmentation therapy. Supervised and unsupervised differential gene expression analyses were performed using Weighted Gene Co-expression Network Analysis (WGCNA) to identify gene profiles associated with subjects' clinical variables. The genes in the most significant WGCNA module were used to cluster AATD individuals. Gene validation was performed by NanoString nCounter Gene Expression Assay. RESULT: We observed modest effects of AATD genotype and augmentation therapy on gene expression. When WGCNA was applied to BAL transcriptome, one gene module, ME31 (2312 genes), correlated with the highest number of clinical variables and was functionally enriched with numerous immune T-lymphocyte related pathways. This gene module identified two distinct clusters of AATD individuals with different disease severity and distinct PBMC gene expression patterns. CONCLUSIONS: We successfully identified novel clusters of AATD individuals where severity correlated with increased immune response independent of individuals' genotype and augmentation therapy. These findings may suggest the presence of previously unrecognised disease endotypes in AATD that associate with T-lymphocyte immunity and disease severity.

Transcriptomics of bronchoalveolar lavage cells identifies new molecular endotypes of sarcoidosis.

BACKGROUND: Sarcoidosis is a multisystem granulomatous disease of unknown origin with a variable and often unpredictable course and pattern of organ involvement. In this study we sought to identify specific bronchoalveolar lavage (BAL) cell gene expression patterns indicative of distinct disease phenotypic traits. METHODS: RNA sequencing by Ion Torrent Proton was performed on BAL cells obtained from 215 well-characterised patients with pulmonary sarcoidosis enrolled in the multicentre Genomic Research in Alpha-1 Antitrypsin Deficiency and Sarcoidosis (GRADS) study. Weighted gene co-expression network analysis and nonparametric statistics were used to analyse genome-wide BAL transcriptome. Validation of results was performed using a microarray expression dataset of an independent sarcoidosis cohort (Freiburg, Germany; n=50). RESULTS: Our supervised analysis found associations between distinct transcriptional programmes and major pulmonary phenotypic manifestations of sarcoidosis including T-helper type 1 (Th1) and Th17 pathways associated with hilar lymphadenopathy, transforming growth factor-ß1 (TGFB1) and mechanistic target of rapamycin (MTOR) signalling with parenchymal involvement, and interleukin (IL)-7 and IL-2 with airway involvement. Our unsupervised analysis revealed gene modules that uncovered four potential sarcoidosis endotypes including hilar lymphadenopathy with increased acute T-cell immune response; extraocular organ involvement with PI3K activation pathways; chronic and multiorgan disease with increased immune response pathways; and multiorgan involvement, with increased IL-1 and IL-18 immune and inflammatory responses. We validated the occurrence of these endotypes using gene expression, pulmonary function tests and cell differentials from Freiburg. CONCLUSION: Taken together, our results identify BAL gene expression programmes that characterise major pulmonary sarcoidosis phenotypes and suggest the presence of distinct disease molecular endotypes.

Automated quantification of COVID-19 severity and progression using chest CT images.

OBJECTIVE: To develop and test computer software to detect, quantify, and monitor progression of pneumonia associated with COVID-19 using chest CT scans. METHODS: One hundred twenty chest CT scans from subjects with lung infiltrates were used for training deep learning algorithms to segment lung regions and vessels. Seventy-two serial scans from 24 COVID-19 subjects were used to develop and test algorithms to detect and quantify the presence and progression of infiltrates associated with COVID-19. The algorithm included (1) automated lung boundary and vessel segmentation, (2) registration of the lung boundary between serial scans, (3) computerized identification of the pneumonitis regions, and (4) assessment of disease progression. Agreement between radiologist manually delineated regions and computer-detected regions was assessed using the Dice coefficient. Serial scans were registered and used to generate a heatmap visualizing the change between scans. Two radiologists, using a five-point Likert scale, subjectively rated heatmap accuracy in representing progression. RESULTS: There was strong agreement between computer detection and the manual delineation of pneumonic regions with a Dice coefficient of 81% (CI 76-86%). In detecting large pneumonia regions (> 200 mm3), the algorithm had a sensitivity of 95% (CI 94-97%) and specificity of 84% (CI 81-86%). Radiologists rated 95% (CI 72 to 99) of heatmaps at least "acceptable" for representing disease progression. CONCLUSION: The preliminary results suggested the feasibility of using computer software to detect and quantify pneumonic regions associated with COVID-19 and to generate heatmaps that can be used to visualize and assess progression. KEY POINTS: • Both computer vision and deep learning technology were used to develop computer software to quantify the presence and progression of pneumonia associated with COVID-19 depicted on CT images. • The computer software was tested using both quantitative experiments and subjective assessment. • The computer software has the potential to assist in the detection of the pneumonic regions, monitor disease progression, and assess treatment efficacy related to COVID-19.

Growth and tolerance of healthy, term infants fed lower protein extensively hydrolyzed or amino acid-based formula: double-blind, randomized, controlled trial.

BACKGROUND: Optimal protein level in hypoallergenic infant formulas is an area of ongoing investigation. The aim was to evaluate growth of healthy term infants who received extensively hydrolyzed (EH) or amino acid (AA)-based formulas with reduced protein. METHODS: In this prospective, multi-center, double-blind, controlled, parallel group study, infants were randomized to receive a marketed EH casein infant formula at 2.8 g protein/100 kcal (Control) or one of two investigational formulas: EH casein formula at 2.4 g protein/100 kcal (EHF) or AA-based formula at 2.4 g total protein equivalents/100 kcal (AAF). Control and EHF each had 2 × 107 CFU Lactobacillus rhamnosus GG/100 kcal. Anthropometrics were measured and recall of formula intake, tolerance, and stool characteristics was collected at 14, 30, 60, 90, 120 days of age. Primary outcome was weight growth rate (g/day) between 14 and 120 days of age (analyzed by ANOVA). Medically confirmed adverse events were recorded throughout the study. RESULTS: No group differences in weight or length growth rate from 14 to 120 days were detected. With the exception of significant differences at several study time points for males, no group differences were detected in mean head circumference growth rates. However, mean achieved weight, length, and head circumference demonstrated normal growth throughout the study period. No group differences in achieved weight or length (males and females) and head circumference (females) were detected and means were within the WHO growth 25th and 75th percentiles from 14 to 120 days of age. With the exception of Day 90, there were no statistically significant group differences in achieved head circumference for males; means remained between the WHO 50th and 75th percentiles for growth at Days 14, 30, and 60 and continued along the 75th percentile through Day 120. No differences in study discontinuation due to formula were detected. The number of participants for whom at least one adverse event was reported was similar among groups. CONCLUSIONS: This study demonstrated hypoallergenic infant formulas at 2.4 g protein/100 kcal were safe, well-tolerated, and associated with appropriate growth in healthy term infants from 14 to 120 days of age. TRIAL REGISTRATION: ClinicalTrials.gov, ClinicalTrials.gov Identifier: NCT01354366 . Registered 13 May 2011.

Routinely reported ejection fraction and mortality in clinical practice: where does the nadir of risk lie?

AIMS: We investigated the relationship between clinically assessed left ventricular ejection fraction (LVEF) and survival in a large, heterogeneous clinical cohort. METHODS AND RESULTS: Physician-reported LVEF on 403 977 echocardiograms from 203 135 patients were linked to all-cause mortality using electronic health records (1998-2018) from US regional healthcare system. Cox proportional hazards regression was used for analyses while adjusting for many patient characteristics including age, sex, and relevant comorbidities. A dataset including 45 531 echocardiograms and 35 976 patients from New Zealand was used to provide independent validation of analyses. During follow-up of the US cohort, 46 258 (23%) patients who had undergone 108 578 (27%) echocardiograms died. Overall, adjusted hazard ratios (HR) for mortality showed a u-shaped relationship for LVEF with a nadir of risk at an LVEF of 60-65%, a HR of 1.71 [95% confidence interval (CI) 1.64-1.77] when ≥70% and a HR of 1.73 (95% CI 1.66-1.80) at LVEF of 35-40%. Similar relationships with a nadir at 60-65% were observed in the validation dataset as well as for each age group and both sexes. The results were similar after further adjustments for conditions associated with an elevated LVEF, including mitral regurgitation, increased wall thickness, and anaemia and when restricted to patients reported to have heart failure at the time of the echocardiogram. CONCLUSION: Deviation of LVEF from 60% to 65% is associated with poorer survival regardless of age, sex, or other relevant comorbidities such as heart failure. These results may herald the recognition of a new phenotype characterized by supra-normal LVEF.

Genomics-First Evaluation of Heart Disease Associated With Titin-Truncating Variants.

BACKGROUND: Truncating variants in the Titin gene (TTNtvs) are common in individuals with idiopathic dilated cardiomyopathy (DCM). However, a comprehensive genomics-first evaluation of the impact of TTNtvs in different clinical contexts, and the evaluation of modifiers such as genetic ancestry, has not been performed. METHODS: We reviewed whole exome sequence data for >71 000 individuals (61 040 from the Geisinger MyCode Community Health Initiative (2007 to present) and 10 273 from the PennMedicine BioBank (2013 to present) to identify anyone with TTNtvs. We further selected individuals with TTNtvs in exons highly expressed in the heart (proportion spliced in [PSI] >0.9). Using linked electronic health records, we evaluated associations of TTNtvs with diagnoses and quantitative echocardiographic measures, including subanalyses for individuals with and without DCM diagnoses. We also reviewed data from the Jackson Heart Study to validate specific analyses for individuals of African ancestry. RESULTS: Identified with a TTNtv in a highly expressed exon (hiPSI) were 1.2% individuals in PennMedicine BioBank and 0.6% at Geisinger. The presence of a hiPSI TTNtv was associated with increased odds of DCM in individuals of European ancestry (odds ratio [95% CI]: 18.7 [9.1-39.4] {PennMedicine BioBank} and 10.8 [7.0-16.0] {Geisinger}). hiPSI TTNtvs were not associated with DCM in individuals of African ancestry, despite a high DCM prevalence (odds ratio, 1.8 [0.2-13.7]; P=0.57). Among 244 individuals of European ancestry with DCM in PennMedicine BioBank, hiPSI TTNtv carriers had lower left ventricular ejection fraction (ß=-12%, P=3×10-7), and increased left ventricular diameter (ß=0.65 cm, P=9×10-3). In the Geisinger cohort, hiPSI TTNtv carriers without a cardiomyopathy diagnosis had more atrial fibrillation (odds ratio, 2.4 [1.6-3.6]) and heart failure (odds ratio, 3.8 [2.4-6.0]), and lower left ventricular ejection fraction (ß=-3.4%, P=1×10-7). CONCLUSIONS: Individuals of European ancestry with hiPSI TTNtv have an abnormal cardiac phenotype characterized by lower left ventricular ejection fraction, irrespective of the clinical manifestation of cardiomyopathy. Associations with arrhythmias, including atrial fibrillation, were observed even when controlling for cardiomyopathy diagnosis. In contrast, no association between hiPSI TTNtvs and DCM was discerned among individuals of African ancestry. Given these findings, clinical identification of hiPSI TTNtv carriers may alter clinical management strategies.