Dynamic imaging using Motion-Compensated SmooThness Regularization on Manifolds (MoCo-SToRM)

We introduce an unsupervised motion-compensated reconstruction scheme for high-resolution free-breathing pulmonary MRI. We model the image frames in the time series as the deformed version of the 3D template image volume. We assume the deformation maps to be points on a smooth manifold in high-dimensional space. Specifically, we model the deformation map at each time instant as the output of a CNN-based generator that has the same weight for all time-frames, driven by a low-dimensional latent vector. The time series of latent vectors account for the dynamics in the dataset, including respiratory motion and bulk motion. The template image volume, the parameters of the generator, and the latent vectors are learned directly from the k-t space data in an unsupervised fashion. Our experimental results show improved reconstructions compared to state-of-the-art methods, especially in the context of bulk motion during the scans

Large-scale CFD simulations of the transitional and turbulent regime for the large human airways during rapid inhalation

The dynamics of unsteady flow in the human large airways during a rapid inhalation were investigated using highly detailed large-scale computational fluid dynamics on a subject-specific geometry. The simulations were performed to resolve all the spatial and temporal scales of the flow, thanks to the use of massive computational resources. A highly parallel finite element code was used, running on two supercomputers, solving the transient incompressible Navier–Stokes equations on unstructured meshes. Given that the finest mesh contained 350 million elements, the study sets a precedent for large-scale simulations of the respiratory system, proposing an analysis strategy for mean flow, fluctuations and wall shear stresses on a rapid and short inhalation (a so-called sniff). The geometry used encompasses the exterior face and the airways from the nasal cavity, through the trachea and up to the third lung bifurcation; it was derived from a contrast-enhanced computed tomography (CT) scan of a 48-year-old male. The transient inflow produces complex flows over a wide range of Reynolds numbers (Re). Thanks to the high fidelity simulations, many features involving the flow transition were observed, with the level of turbulence clearly higher in the throat than in the nose. Spectral analysis revealed turbulent characteristics persisting downstream of the glottis, and were captured even with a medium mesh resolution. However a fine mesh resolution was found necessary in the nasal cavity to observe transitional features. This work indicates the potential of large-scale simulations to further understanding of airway physiological mechanics, which is essential to guide clinical diagnosis; better understanding of the flow also has implications for the design of interventions such as aerosol drug delivery.We acknowledge PRACE for awarding us access to resource FERMI based in Italy at Bologna hosted by Cineca. This work was financially supported by the PRACE project Pra04 693 (2011050693 to the Fourth PRACE regular call). The second author gratefully acknowledges support from project ‘MatComPhys’ under the European Research Executive Agency FP7-PEOPLE-2011- IEF framework. The third author was supported by the Engineering and Physical Sciences Research Council [grant number EP/ M506345/1].Peer ReviewedPostprint (author's final draft

Multiple novel prostate cancer susceptibility signals identified by fine-mapping of known risk loci among Europeans

Genome-wide association studies (GWAS) have identified numerous common prostate cancer (PrCa) susceptibility loci. We have fine-mapped 64 GWAS regions known at the conclusion of the iCOGS study using large-scale genotyping and imputation in 25 723 PrCa cases and 26 274 controls of European ancestry. We detected evidence for multiple independent signals at 16 regions, 12 of which contained additional newly identified significant associations. A single signal comprising a spectrum of correlated variation was observed at 39 regions; 35 of which are now described by a novel more significantly associated lead SNP, while the originally reported variant remained as the lead SNP only in 4 regions. We also confirmed two association signals in Europeans that had been previously reported only in East-Asian GWAS. Based on statistical evidence and linkage disequilibrium (LD) structure, we have curated and narrowed down the list of the most likely candidate causal variants for each region. Functional annotation using data from ENCODE filtered for PrCa cell lines and eQTL analysis demonstrated significant enrichment for overlap with bio-features within this set. By incorporating the novel risk variants identified here alongside the refined data for existing association signals, we estimate that these loci now explain ∼38.9% of the familial relative risk of PrCa, an 8.9% improvement over the previously reported GWAS tag SNPs. This suggests that a significant fraction of the heritability of PrCa may have been hidden during the discovery phase of GWAS, in particular due to the presence of multiple independent signals within the same regio

Large-scale CFD simulations of the transitional and turbulent regime for the large human airways during rapid inhalation

The dynamics of unsteady flow in the human large airways during a rapid inhalation were investigated using highly detailed large-scale computational fluid dynamics on a subject-specific geometry. The simulations were performed to resolve all the spatial and temporal scales of the flow, thanks to the use of massive computational resources. A highly parallel finite element code was used, running on two supercomputers, solving the transient incompressible Navier–Stokes equations on unstructured meshes. Given that the finest mesh contained 350 million elements, the study sets a precedent for large-scale simulations of the respiratory system, proposing an analysis strategy for mean flow, fluctuations and wall shear stresses on a rapid and short inhalation (a so-called sniff). The geometry used encompasses the exterior face and the airways from the nasal cavity, through the trachea and up to the third lung bifurcation; it was derived from a contrast-enhanced computed tomography (CT) scan of a 48-year-old male. The transient inflow produces complex flows over a wide range of Reynolds numbers (Re). Thanks to the high fidelity simulations, many features involving the flow transition were observed, with the level of turbulence clearly higher in the throat than in the nose. Spectral analysis revealed turbulent characteristics persisting downstream of the glottis, and were captured even with a medium mesh resolution. However a fine mesh resolution was found necessary in the nasal cavity to observe transitional features. This work indicates the potential of large-scale simulations to further understanding of airway physiological mechanics, which is essential to guide clinical diagnosis; better understanding of the flow also has implications for the design of interventions such as aerosol drug delivery.We acknowledge PRACE for awarding us access to resource FERMI based in Italy at Bologna hosted by Cineca. This work was financially supported by the PRACE project Pra04 693 (2011050693 to the Fourth PRACE regular call). The second author gratefully acknowledges support from project ‘MatComPhys’ under the European Research Executive Agency FP7-PEOPLE-2011- IEF framework. The third author was supported by the Engineering and Physical Sciences Research Council [grant number EP/ M506345/1].Peer Reviewe

Dynamics of Airflow in a Short Inhalation

During a rapid inhalation, such as a sniff, the flowin the airways accelerates and decays quickly. The consequences for flow development and convective transport of an inhaled gas were investigated in a subject geometry extending from the nose to the bronchi. The progress of flow transition and the advance of an inhaled non-absorbed gas were determined using highly resolved simulations of a sniff 0.5 s long, 1 l s21 peak flow, 364 ml inhaled volume. In the nose, the distribution of airflow evolved through three phases: (i) an initial transient of about 50 ms, roughly the filling time for a nasal volume, (ii) quasi-equilibrium over themajority of the inhalation, and (iii) a terminating phase. Flow transition commenced in the supraglottic region within 20 ms, resulting in largeamplitude fluctuations persisting throughout the inhalation; in the nose, fluctuations that arose nearer peak flow were of much reduced intensity and diminished in the flow decay phase. Measures of gas concentration showed non-uniform build-up and wash-out of the inhaled gas in the nose. At the carina, the form of the temporal concentration profile reflected both shear dispersion and airway filling defects owing to recirculation regions.This research was supported by EPSRC Doctoral Training Award EP/P505550/1.Electronic supplementary material is available at http://dx.doi.org/10.1098/rsif.2014.0880 or via http://rsif.royalsocietypublishing.org.Peer Reviewe

Pre- and post-operative visualization of neonatal esophageal atresia/tracheoesophageal fistula via magnetic resonance imaging

Esophageal atresia (EA) is a relatively uncommon congenital anomaly, often observed in conjunction with tracheoesophageal fistula (TEF). Surgical repair in neonates typically takes place with little information about the pre-existing EA/TEF structure because there are currently no acceptable tools for evaluating EA/TEF anatomy prior to repair; chest x-ray radiograph does not identify malformation sub-type or gap length, while x-ray computed tomography (CT) demonstrates an unacceptably high exposure to ionizing radiation. There is a need for safe imaging methods to evaluate pre-operative EA/TEF anatomy, which would add value in surgical planning; this need may be met with high-resolution structural MRI. We report three cases of Type-C EA/TEF in neonates. Patients were imaged prior to surgical repair using high-resolution ultrashort echo time (UTE) magnetic resonance imaging (MRI) to visualize tracheoesophageal anatomy and allow for informed surgical planning and risk management. One of the three patients was imaged post-repair to evaluate surgical efficacy and evolution of the tracheoesophageal anatomy

Multiple loci on 8q24 associated with prostate cancer susceptibility

Previous studies have identified multiple loci on 8q24 associated with prostate cancer risk. We performed a comprehensive analysis of SNP associations across 8q24 by genotyping tag SNPs in 5,504 prostate cancer cases and 5,834 controls. We confirmed associations at three previously reported loci and identified additional loci in two other linkage disequilibrium blocks (rs1006908: per-allele OR = 0.87, P = 7.9 x 10(-8); rs620861: OR = 0.90, P = 4.8 x 10(-8)). Eight SNPs in five linkage disequilibrium blocks were independently associated with prostate cancer susceptibility

Identification of seven new prostate cancer susceptibility loci through a genome-wide association study

Prostate cancer (PrCa) is the most frequently diagnosed male cancer in developed countries. To identify common PrCa susceptibility alleles, we have previously conducted a genome-wide association study in which 541, 129 SNPs were genotyped in 1,854 PrCa cases with clinically detected disease and 1,894 controls. We have now evaluated promising associations in a second stage, in which we genotyped 43,671 SNPs in 3,650 PrCa cases and 3,940 controls, and a third stage, involving an additional 16,229 cases and 14,821 controls from 21 studies. In addition to previously identified loci, we identified a further seven new prostate cancer susceptibility loci on chromosomes 2, 4, 8, 11, and 22 (P=1.6×10−8 to P=2.7×10−33)