Combining Partial True Discovery Guarantee Procedures.

Closed testing has recently been shown to be optimal for simultaneous true discovery proportion control. It is, however, challenging to construct true discovery guarantee procedures in such a way that it focuses power on some feature sets chosen by users based on their specific interest or expertise. We propose a procedure that allows users to target power on prespecified feature sets, that is, "focus sets." Still, the method also allows inference for feature sets chosen post hoc, that is, "nonfocus sets," for which we deduce a true discovery lower confidence bound by interpolation. Our procedure is built from partial true discovery guarantee procedures combined with Holm's procedure and is a conservative shortcut to the closed testing procedure. A simulation study confirms that the statistical power of our method is relatively high for focus sets, at the cost of power for nonfocus sets, as desired. In addition, we investigate its power property for sets with specific structures, for example, trees and directed acyclic graphs. We also compare our method with AdaFilter in the context of replicability analysis. The application of our method is illustrated with a gene ontology analysis in gene expression data.

Clarifying the biological and statistical assumptions of cross-sectional biological age predictors: an elaborate illustration using synthetic and real data.

BACKGROUND: There is divergence in the rate at which people age. The concept of biological age is postulated to capture this variability, and hence to better represent an individual's true global physiological state than chronological age. Biological age predictors are often generated based on cross-sectional data, using biochemical or molecular markers as predictor variables. It is assumed that the difference between chronological and predicted biological age is informative of one's chronological age-independent aging divergence ∆. METHODS: We investigated the statistical assumptions underlying the most popular cross-sectional biological age predictors, based on multiple linear regression, the Klemera-Doubal method or principal component analysis. We used synthetic and real data to illustrate the consequences if this assumption does not hold. RESULTS: The most popular cross-sectional biological age predictors all use the same strong underlying assumption, namely that a candidate marker of aging's association with chronological age is directly informative of its association with the aging rate ∆. We called this the identical-association assumption and proved that it is untestable in a cross-sectional setting. If this assumption does not hold, weights assigned to candidate markers of aging are uninformative, and no more signal may be captured than if markers would have been assigned weights at random. CONCLUSIONS: Cross-sectional methods for predicting biological age commonly use the untestable identical-association assumption, which previous literature in the field had never explicitly acknowledged. These methods have inherent limitations and may provide uninformative results, highlighting the importance of researchers exercising caution in the development and interpretation of cross-sectional biological age predictors.

RNA-sequencing to discover genes and signaling pathways associated with venous thromboembolism in glioblastoma patients: A case-control study.

BACKGROUND: Glioblastoma patients are at high risk of developing venous thromboembolism (VTE). Tumor-intrinsic features are considered to play a role, but the underlying pathophysiological mechanisms remain incompletely understood. OBJECTIVES: To identify tumor-expressed genes and signaling pathways that associate with glioblastoma-related VTE by using next generation RNA-sequencing (RNA-Seq). METHODS: The tumor gene expression profile of 23 glioblastoma patients with VTE and 23 glioblastoma patients without VTE was compared using an unpaired analysis. Ingenuity Pathway Analysis (IPA) core analysis was performed on the top 50 differentially expressed genes to explore associated functions and pathways. Based on full RNA-Seq data, molecular glioblastoma subtypes were determined by performing cluster analysis. RESULTS: Of the 19,327 genes, 1246 (6.4 %) were differentially expressed between glioblastoma patients with and without VTE (unadjusted P < 0.05). The most highly overexpressed gene was GLI1, a classical target gene in the Sonic Hedgehog (Shh) signaling pathway (log2 fold change: 3.7; unadjusted P < 0.0001, adjusted P = 0.219). In line, Shh signaling was among the top canonical pathways and processes associated with VTE. The proportion of patients with the proneural/neural glioblastoma subtype was higher among those with VTE than controls. CONCLUSION: Shh signaling may be involved in the development of glioblastoma-related VTE.

Histological Chorioamnionitis and Funisitis as New Risk Factors for Retinopathy of Prematurity: A Meta-analysis.

OBJECTIVE: The role of placental inflammation in neonatal morbidities is underestimated due to lack of placental examination. This meta-analysis aims to assess the association between histological chorioamnionitis (HCA) with and without funisitis (FUN) and risk of retinopathy of prematurity (ROP). STUDY DESIGN: Forty-five studies reporting (unadjusted) data on HCA without FUN and HCA with FUN in neonates with ROP were included. Primary outcomes were any stage ROP and severe ROP. Potential confounders explored were gestational age (GA) at birth, birthweight, maternal steroid use, necrotizing enterocolitis, sepsis (suspected/proven) and mechanical ventilation duration. RESULTS: Neonates with HCA had increased risk for any stage ROP (odds ratio [OR] 1.8; 95% confidence interval [CI] 1.3-2.4) and severe ROP (OR 1.5; 95% CI 1.2-1.8) compared with neonates without HCA. The rates of any stage ROP (OR 1.8; 95% CI 1.4-2.2) and severe ROP (OR 1.4; 95% CI 1.1-1.6) were higher in neonates with FUN compared with neonates without FUN. Multivariate meta-regression analysis suggests that lower GA increases the effect size between FUN and severe ROP. CONCLUSION: This meta-analysis confirms that presence of HCA and FUN are risk factors for any stage ROP and severe ROP. Structured histological placental examination of HCA and FUN may be a tool to further refine the ROP risk profile. KEY POINTS: · This systematic review confirms that HCA is a risk factor for ROP.. · This meta-analysis reveals that FUN results in an even higher risk for developing ROP.. · Placental examination of HCA/FUN may be a tool to further refine the ROP risk profile..

Development and External Validation of a PET Radiomic Model for Prognostication of Head and Neck Cancer.

AIM: To build and externally validate an [18F]FDG PET radiomic model to predict overall survival in patients with head and neck squamous cell carcinoma (HNSCC). METHODS: Two multicentre datasets of patients with operable HNSCC treated with preoperative afatinib who underwent a baseline and evaluation [18F]FDG PET/CT scan were included (EORTC: n = 20, Unicancer: n = 34). Tumours were delineated, and radiomic features were extracted. Each cohort served once as a training and once as an external validation set for the prediction of overall survival. Supervised feature selection was performed using variable hunting with variable importance, selecting the top two features. A Cox proportional hazards regression model using selected radiomic features and clinical characteristics was fitted on the training dataset and validated in the external validation set. Model performances are expressed by the concordance index (C-index). RESULTS: In both models, the radiomic model surpassed the clinical model with validation C-indices of 0.69 and 0.79 vs. 0.60 and 0.67, respectively. The model that combined the radiomic features and clinical variables performed best, with validation C-indices of 0.71 and 0.82. CONCLUSION: Although assessed in two small but independent cohorts, an [18F]FDG-PET radiomic signature based on the evaluation scan seems promising for the prediction of overall survival for HNSSC treated with preoperative afatinib. The robustness and clinical applicability of this radiomic signature should be assessed in a larger cohort.

Cancer-Associated Fibroblasts Are Key Determinants of Cancer Cell Invasion in the Earliest Stage of Colorectal Cancer.

BACKGROUND & AIMS: Improving clinical management of early stage colorectal cancers (T1CRCs) requires a better understanding of their underlying biology. Accumulating evidence shows that cancer-associated fibroblasts (CAFs) are important determinants of tumor progression in advanced colorectal cancer (CRC), but their role in the initial stages of CRC tumorigenesis is unknown. Therefore, we investigated the contribution of T1CAFs to early CRC progression. METHODS: Primary T1CAFs and patient-matched normal fibroblasts (NFs) were isolated from endoscopic biopsy specimens of histologically confirmed T1CRCs and normal mucosa, respectively. The impact of T1CAFs and NFs on tumor behavior was studied using 3-dimensional co-culture systems with primary T1CRC organoids and extracellular matrix (ECM) remodeling assays. Whole-transcriptome sequencing and gene silencing were used to pinpoint mediators of T1CAF functions. RESULTS: In 3-dimensional multicellular cultures, matrix invasion of T1CRC organoids was induced by T1CAFs, but not by matched NFs. Enhanced T1CRC invasion was accompanied by T1CAF-induced ECM remodeling and up-regulation of CD44 in epithelial cells. RNA sequencing of 10 NF-T1CAF pairs revealed 404 differentially expressed genes, with significant enrichment for ECM-related pathways in T1CAFs. Cathepsin H, a cysteine-type protease that was specifically up-regulated in T1CAFs but not in fibroblasts from premalignant lesions or advanced CRCs, was identified as a key factor driving matrix remodeling by T1CAFs. Finally, we showed high abundance of cathepsin H-expressing T1CAFs at the invasive front of primary T1CRC sections. CONCLUSIONS: Already in the earliest stage of CRC, cancer cell invasion is promoted by CAFs via direct interactions with epithelial cancer cells and stage-specific, cathepsin H-dependent ECM remodeling. RNA sequencing data of the 10 NF-T1CAF pairs can be found under GEO accession number GSE200660.

Molecular Landscape of Pelvic Organ Prolapse Provides Insights into Disease Etiology.

Pelvic organ prolapse (POP) represents a major health care burden in women, but its underlying pathophysiological mechanisms have not been elucidated. We first used a case-control design to perform an exome chip study in 526 women with POP and 960 control women to identify single nucleotide variants (SNVs) associated with the disease. We then integrated the functional interactions between the POP candidate proteins derived from the exome chip study and other POP candidate molecules into a molecular landscape. We found significant associations between POP and SNVs in 54 genes. The proteins encoded by 26 of these genes fit into the molecular landscape, together with 43 other POP candidate molecules. The POP landscape is located in and around epithelial cells and fibroblasts of the urogenital tract and harbors four interacting biological processes-epithelial-mesenchymal transition, immune response, modulation of the extracellular matrix, and fibroblast function-that are regulated by sex hormones and TGFB1. Our findings were corroborated by enrichment analyses of differential gene expression data from an independent POP cohort. Lastly, based on the landscape and using vaginal fibroblasts from women with POP, we predicted and showed that metformin alters gene expression in these fibroblasts in a beneficial direction. In conclusion, our integrated molecular landscape of POP provides insights into the biological processes underlying the disease and clues towards novel treatments.

Dynamic characteristics of skin reaction force in different body postures.

Mechanical stress influences scarring of a surgical wound. Several lines have been proposed for the best excision direction. It is unknown if these lines still apply when the body posture changes. The objective is to measure the skin reaction force in four directions and determine the direction of least force. Secondary objective is to determine if the reaction force varies in a different body posture. Skin reaction force was measured with the compressiometer in 30 participants on four different locations (forearm/upper arm/shoulder blade/lower back) in four directions (0°-45°-90°-135°) and two body postures. The direction of least skin reaction force changed with a different body posture and was significant for the forearm (p < 0.01) and shoulder blade (p = 0.05) The skin reaction force in all four direction changed significantly in a different body posture, except the 45° line in the upper arm and shoulder blade. Our results demonstrate that the skin reaction force in four directions in four locations varies with change in body posture. Focus should therefore not only lay on choosing the right direction, but also on managing skin tension postoperatively.

HLA expression as a risk factor for metastases of cutaneous squamous-cell carcinoma in organ- transplant recipients.

BACKGROUND: Solid organ-transplant recipients (SOTR) have an increased risk of cutaneous squamous-cell carcinoma (cSCC), metastasis and death from cSCC. In immunocompetent patients with mucosal SCC, downregulation of HLA class I is associated with poor prognosis. Since the degree of HLA expression on tumor cells could play a role in immunogenicity and pathophysiology of cSCC metastasis, we hypothesized that decreased HLA expression is associated with an increased risk of metastasis. METHODS: We compared HLA expression between primary metastasized cSCCs, their metastases, and non-metastasized cSCCs from the same patients. Samples were stained for HLA-A, HLA-B/-C and quantified by calculating the difference in immunoreactivity score (IRS) of the primary cSCC compared with all non-metastasized cSCCs. RESULTS: The mean IRS score for HLA-B/C expression was 2.07 point higher in metastasized compared to non-metastasized cSCCs (p = 0.065, 95 % CI -0.18-4.32). 83.3 % of the primary metastasized cSCCs had an IRS score of 4 or higher, compared to 42.9 % in non-metastasized cSCCs. Moderately to poorly differentiated cSCCs had more HLA class I expression compared to well-differentiated cSCCs. CONCLUSION: Contrary to immunocompetent patients, HLA-B/C expression tends to be upregulated in metastasized cSCC compared to non-metastasized cSCC in SOTR, suggesting that different tumor escape mechanisms play a role in SOTR compared to immunocompetent patients.

Closed testing with Globaltest, with application in metabolomics.

The Globaltest is a powerful test for the global null hypothesis that there is no association between a group of features and a response of interest, which is popular in pathway testing in metabolomics. Evaluating multiple feature sets, however, requires multiple testing correction. In this paper, we propose a multiple testing method, based on closed testing, specifically designed for the Globaltest. The proposed method controls the familywise error rate simultaneously over all possible feature sets, and therefore allows post hoc inference, that is, the researcher may choose feature sets of interest after seeing the data without jeopardizing error control. To circumvent the exponential computation time of closed testing, we derive a novel shortcut that allows exact closed testing to be performed on the scale of metabolomics data. An R package ctgt is available on comprehensive R archive network for the implementation of the shortcut procedure, with applications on several real metabolomics data examples.

Cerebrospinal Fluid and Plasma Amine Profiles in Interictal Migraine.

OBJECTIVE: Impaired amine metabolism has been associated with the etiology of migraine, that is, why patients continue to get migraine attacks. However, evidence from cerebrospinal fluid (CSF) is lacking. Here, we evaluated individual amine levels, global amine profiles, and amine pathways in CSF and plasma of interictal migraine patients and healthy controls. METHODS: CSF and plasma were sampled between 8:30 am and 1:00 pm, randomly and interchangeably over the time span to avoid any diurnal and seasonal influences, from healthy volunteers and interictal migraine patients, matched for age, sex, and sampling time. The study was approved by the local medical ethics committee. Individual amines (n = 31), global amine profiles, and specific amine pathways were analyzed using a validated ultraperformance liquid chromatography mass spectrometry platform. RESULTS: We analyzed n = 99 participants with migraine with aura, n = 98 with migraine without aura, and n = 96 healthy volunteers. Univariate analysis with Bonferroni correction indicated that CSF L-arginine was reduced in migraine with aura (10.4%, p < 0.001) and without aura (5.0%, p = 0.03). False discovery rate-corrected CSF L-phenylalanine was also lower in migraine with aura (6.9%, p = 0.011) and without aura (8.1%, p = 0.001), p = 0.088 after Bonferroni correction. Multivariate analysis revealed that CSF global amine profiles were similar for both types of migraine (p = 0.64), but distinct from controls (p = 0.009). Global profile analyses were similar in plasma. The strongest associated pathways with migraine were related to L-arginine metabolism. INTERPRETATION: L-Arginine was decreased in the CSF (but not in plasma) of interictal patients with migraine with or without aura, and associated pathways were altered. This suggests that dysfunction of nitric oxide signaling is involved in susceptibility to getting migraine attacks. ANN NEUROL 2023;93:715-728.

NetTDP: permutation-based true discovery proportions for differential co-expression network analysis.

Existing methods for differential network analysis could only infer whether two networks of interest have differences between two groups of samples, but could not quantify and localize network differences. In this work, a novel method, permutation-based Network True Discovery Proportions (NetTDP), is proposed to quantify the number of edges (correlations) or nodes (genes) for which the co-expression networks are different. In the NetTDP method, we propose an edge-level statistic and a node-level statistic, and detect true discoveries of edges and nodes in the sense of differential co-expression network, respectively, by the permutation-based sumSome method. Furthermore, the NetTDP method could further localize the differences by inferring the TDPs for edge or gene subsets of interest, which can be selected post hoc. Our NetTDP method allows inference on data-driven modules or biology-driven gene sets, and remains valid even when these sub-networks are optimized using the same data. Experimental results on both simulation data sets and five real data sets show the effectiveness of the proposed method in inferring the quantification and localization of differential co-expression networks. The R code is available at https://github.com/LiminLi-xjtu/NetTDP.

Immune Determinants of Viral Clearance in Hospitalised COVID-19 Patients: Reduced Circulating Naïve CD4+ T Cell Counts Correspond with Delayed Viral Clearance.

Virus-specific cellular and humoral responses are major determinants for protection from critical illness after SARS-CoV-2 infection. However, the magnitude of the contribution of each of the components to viral clearance remains unclear. Here, we studied the timing of viral clearance in relation to 122 immune parameters in 102 hospitalised patients with moderate and severe COVID-19 in a longitudinal design. Delayed viral clearance was associated with more severe disease and was associated with higher levels of SARS-CoV-2-specific (neutralising) antibodies over time, increased numbers of neutrophils, monocytes, basophils, and a range of pro-inflammatory cyto-/chemokines illustrating ongoing, partially Th2 dominating, immune activation. In contrast, early viral clearance and less critical illness correlated with the peak of neutralising antibodies, higher levels of CD4 T cells, and in particular naïve CD4+ T cells, suggesting their role in early control of SARS-CoV-2 possibly by proving appropriate B cell help. Higher counts of naïve CD4+ T cells also correlated with lower levels of MIF, IL-9, and TNF-beta, suggesting an indirect role in averting prolonged virus-induced tissue damage. Collectively, our data show that naïve CD4+ T cell play a critical role in rapid viral T cell control, obviating aberrant antibody and cytokine profiles and disease deterioration. These data may help in guiding risk stratification for severe COVID-19.

Comparison of trends in Clostridioides difficile infections in hospitalised patients during the first and second waves of the COVID-19 pandemic: A retrospective sentinel surveillance study.

Background: During the COVID-19 pandemic, several factors, such as improved hand hygiene, social distancing, and restricted hospital referral, may have had an influence on the epidemiology of Clostridioides difficile infections (CDI). Methods: The annual CDI incidence rate of nine hospitals participating in the Dutch sentinel CI surveillance with complete data was compared between 2020 and the previous five surveillance years. Trends in characteristics of hospitalised CDI patients in 21-24 participating hospitals were compared between the first (March 13-May 12, 2020) or second Dutch COVID-19 wave (September 17, 2020-January 1, 2021) and the same calendar periods in 2015 through 2019. All analyses were adjusted for trend changes over time. Findings: The annual CDI incidence rate in 2020 was lower compared to previous years. During the second wave, the percentage of CDI patients with severe CDI was higher compared to earlier (25·8% in 2020 vs 17·9% in 2015-2019 (RR 1·6; 95%CI 1·1-2·3)). After adjustment for delayed C. difficile diagnostics (≥8 days from start symptoms), the increase disappeared. Delayed C. difficile diagnostics was indeed more common during the second wave (RR 1·7; 95%CI 1·1-2·6), but only for community-onset CDI (CO-CDI). Interpretation: This study shows that a higher percentage of severe CDI cases was observed during the second COVID-19 wave. This may partially be caused by delayed diagnostics, potentially due to decreased visits to a physician or restricted hospital referral for CO-CDI patients. Funding: Dutch ministry of Health.

A prospective matched case-control study on the genomic epidemiology of colistin-resistant Enterobacterales from Dutch patients.

Background: Colistin is a last-resort treatment option for infections with multidrug-resistant Gram-negative bacteria. However, colistin resistance is increasing. Methods: A six-month prospective matched case-control study was performed in which 22 Dutch laboratories with 32 associated hospitals participated. Laboratories were invited to send a maximum of five colistin-resistant Escherichia coli or Klebsiella pneumoniae (COLR-EK) isolates and five colistin-susceptible isolates (COLS-EK) to the reference laboratory, matched for patient location, material of origin and bacterial species. Epidemiological/clinical data were collected and included in the analysis. Characteristics of COLR-EK/COLS-EK isolates were compared using logistic regression with correction for variables used for matching. Forty-six ColR-EK/ColS-EK pairs were analysed by next-generation sequencing (NGS) for whole-genome multi-locus sequence typing and identification of resistance genes, including mcr genes. To identify chromosomal mutations potentially leading to colistin resistance, NGS reads were mapped against gene sequences of pmrAB, phoPQ, mgrB and crrB. Results: In total, 72 COLR-EK/COLS-EK pairs (75% E. coli and 25% K. pneumoniae) were included. Twenty-one percent of COLR-EK patients had received colistin, in contrast to 3% of COLS-EK patients (OR > 2.9). Of COLR-EK isolates, five contained mcr-1 and two mcr-9. One isolate lost mcr-9 after repeated sub-culturing, but retained colistin resistance. Among 46 sequenced COLR-EK isolates, genetic diversity was large and 19 (41.3%) isolates had chromosomal mutations potentially associated with colistin resistance. Conclusions: Colistin resistance is present but uncommon in the Netherlands and caused by the mcr gene in a minority of COLR-EK isolates. There is a need for surveillance of colistin resistance using appropriate susceptibility testing methods.

Study Protocol PROMETHEUS: Prospective Multicenter Study to Evaluate the Correlation Between Safety Margin and Local Recurrence After Thermal Ablation Using Image Co-registration in Patients with Hepatocellular Carcinoma.

PURPOSE: The primary objective is to determine the minimal ablation margin required to achieve a local recurrence rate of < 10% in patients with hepatocellular carcinoma undergoing thermal ablation. Secondary objectives are to analyze the correlation between ablation margins and local recurrence and to assess efficacy. MATERIALS AND METHODS: This study is a prospective, multicenter, non-experimental, non-comparative, open-label study. Patients > 18 years with Barcelona Clinic Liver Cancer stage 0/A hepatocellular carcinoma (or B with a maximum of two lesions < 5 cm each) are eligible. Patients will undergo dual-phase contrast-enhanced computed tomography directly before and after ablation. Ablation margins will be quantitatively assessed using co-registration software, blinding assessors (i.e. two experienced radiologists) for outcome. Presence and location of recurrence are evaluated independently on follow-up scans by two other experienced radiologists, blinded for the quantitative margin analysis. A sample size of 189 tumors (~ 145 patients) is required to show with 80% power that the risk of local recurrence is confidently below 10%. A two-sided binomial z-test will be used to test the null hypothesis that the local recurrence rate is ≥ 10% for patients with a minimal ablation margin ≥ 2 mm. Logistic regression will be used to find the relationship between minimal ablation margins and local recurrence. Kaplan-Meier estimates are used to assess local and overall recurrence, disease-free and overall survival. DISCUSSION: It is expected that this study will result in a clear understanding of the correlation between ablation margins and local recurrence. Using co-registration software in future patients undergoing ablation for hepatocellular carcinoma may improve intraprocedural evaluation of technical success. Trial registration The Netherlands Trial Register (NL9713), https://www.trialregister.nl/trial/9713 .

The prognostic value of extracranial vascular characteristics on procedural duration and revascularization success in endovascularly treated acute ischemic stroke patients.

Introduction: Vascular anatomy might affect endovascular treatment success in acute ischemic stroke patients with large vessel occlusion. We investigated the prognostic value of extracranial vascular characteristics on procedural time and revascularization success in patients with large vessel occlusion in the anterior cerebral circulation. Patients and methods: We included 828 patients endovascularly treated within 6.5 hours of symptom onset from the Dutch MR CLEAN-Registry. We evaluated aortic arch configuration, stenosis and tortuosity of supra-aortic arteries, and internal carotid arteries (ICAs) on pre-intervention CTA. We constructed logistic prediction models for outcome variables procedural duration (≥60 minutes) and non-successful revascularization (extended thrombolysis in cerebral infarction (eTICI) of 0-2A) using baseline characteristics and assessed the effect of extracranial vascular characteristics on model performance. Results: Cervical ICA tortuosity and stenosis ≥99% improved prediction of long procedural duration compared with baseline characteristics from area under the curve of 0.61 (95% CI: 0.57-0.65) to 0.66 (95% CI: 0.62-0.70) (P < 0.001). Cervical ICA tortuosity was significantly associated with non-successful recanalization. Prediction of non-successful revascularization did not improve after including aortic arch elongation, acute take-off angle, aortic variant, origin stenosis of supra-aortic arteries, and cervical ICA tortuosity, with an area under the curve of 0.63 (95% CI: 0.59-0.67) compared with 0.59 (95% CI: 0.55-0.63) (P = 0.11). Conclusion: Extracranial vascular characteristics have additional prognostic value for procedural duration, but not for revascularization success, compared with baseline characteristics. Performance of both prediction models is limited in patients treated for large vessel occlusion.

Simultaneous confidence intervals for ranks with application to ranking institutions.

When a ranking of institutions such as medical centers or universities is based on a numerical measure of performance provided with a standard error, confidence intervals (CIs) should be calculated to assess the uncertainty of these ranks. We present a novel method based on Tukey's honest significant difference test to construct simultaneous CIs for the true ranks. When all the true performances are equal, the probability of coverage of our method attains the nominal level. In case the true performance measures have no exact ties, our method is conservative. For this situation, we propose a rescaling method to the nominal level that results in shorter CIs while keeping control of the simultaneous coverage. We also show that a similar rescaling can be applied to correct a recently proposed Monte-Carlo based method, which is anticonservative. After rescaling, the two methods perform very similarly. However, the rescaling of the Monte-Carlo based method is computationally much more demanding and becomes infeasible when the number of institutions is larger than 30-50. We discuss another recently proposed method similar to ours based on simultaneous CIs for the true performance. We show that our method provides uniformly shorter CIs for the same confidence level. We illustrate the superiority of our new methods with a data analysis for travel time to work in the United States and on rankings of 64 hospitals in the Netherlands.

Echo planar imaging-induced errors in intracardiac 4D flow MRI quantification.

PURPOSE: To assess errors associated with EPI-accelerated intracardiac 4D flow MRI (4DEPI) with EPI factor 5, compared with non-EPI gradient echo (4DGRE). METHODS: Three 3T MRI experiments were performed comparing 4DEPI to 4DGRE: steady flow through straight tubes, pulsatile flow in a left-ventricle phantom, and intracardiac flow in 10 healthy volunteers. For each experiment, 4DEPI was repeated with readout and blip phase-encoding gradient in different orientations, parallel or perpendicular to the flow direction. In vitro flow rates were compared with timed volumetric collection. In the left-ventricle phantom and in vivo, voxel-based speed and spatio-temporal median speed were compared between sequences, as well as mitral and aortic transvalvular net forward volume. RESULTS: In steady-flow phantoms, the flow rate error was largest (12%) for high velocity (>2 m/s) with 4DEPI readout gradient parallel to the flow. Voxel-based speed and median speed in the left-ventricle phantom were ≤5.5% different between sequences. In vivo, mean net forward volume inconsistency was largest (6.4 ± 8.5%) for 4DEPI with nonblip phase-encoding gradient parallel to the main flow. The difference in median speed for 4DEPI versus 4DGRE was largest (9%) when the 4DEPI readout gradient was parallel to the flow. CONCLUSIONS: Velocity and flow rate are inaccurate for 4DEPI with EPI factor 5 when flow is parallel to the readout or blip phase-encoding gradient. However, mean differences in flow rate, voxel-based speed, and spatio-temporal median speed were acceptable (≤10%) when comparing 4DEPI to 4DGRE for intracardiac flow in healthy volunteers.