Validity of biomarkers of early circulatory impairment to predict outcome: a retrospective analysis

Objectives: The definition of circulatory impairment in the premature infant is controversial. Current research suggests overdiagnosis and overtreatment. We aimed to analyse which biomarkers move clinicians to initiate cardiovascular treatment (CVT). The prognostic capacity for adverse outcome (death and/or moderate-severe brain damage by cranial ultrasound at term equivalent) of these biomarkers was evaluated. Study Design: Retrospective data analysis from preterm infants enrolled in a placebo-controlled trial on dobutamine for low superior vena cava (SVC) flow, who showed normal SVC flow within the first 24 h (not randomized). Five positive biomarkers were considered: MABP 4 mmol/L; BE < −9 mmol/L; SVC flow <51 ml/kg/min. Results: Ninety eight infants formed the study cohort. Thirty six received CVT (2–95 h). Logistic regression models adjusted for gestational age showed a positive association between CVT and the risk of death or moderate-severe abnormal cranial ultrasound at term equivalent [(OR 5.2, 95%CI: 1.8–15.1) p = 0.002]. MABP 4 mmol/L were the most prevalent biomarkers at start of treatment. Low BE, high serum lactate and low SVC flow at first echocardiography showed a trend toward being associated with adverse outcome, although not statistically significant. Conclusions: Low blood pressure and high lactate are the most prevalent biomarkers used for CVT prescription. Lactic acidosis and low SVC flow early after birth showed a trend toward being associated with adverse outcome. These findings support using a combination of biomarkers for inclusion in a placebo-controlled trial on CVT during transitional circulation

Albiglutide and cardiovascular outcomes in patients with type 2 diabetes and cardiovascular disease (Harmony Outcomes): a double-blind, randomised placebo-controlled trial

Background: Glucagon-like peptide 1 receptor agonists differ in chemical structure, duration of action, and in their effects on clinical outcomes. The cardiovascular effects of once-weekly albiglutide in type 2 diabetes are unknown. We aimed to determine the safety and efficacy of albiglutide in preventing cardiovascular death, myocardial infarction, or stroke. Methods: We did a double-blind, randomised, placebo-controlled trial in 610 sites across 28 countries. We randomly assigned patients aged 40 years and older with type 2 diabetes and cardiovascular disease (at a 1:1 ratio) to groups that either received a subcutaneous injection of albiglutide (30–50 mg, based on glycaemic response and tolerability) or of a matched volume of placebo once a week, in addition to their standard care. Investigators used an interactive voice or web response system to obtain treatment assignment, and patients and all study investigators were masked to their treatment allocation. We hypothesised that albiglutide would be non-inferior to placebo for the primary outcome of the first occurrence of cardiovascular death, myocardial infarction, or stroke, which was assessed in the intention-to-treat population. If non-inferiority was confirmed by an upper limit of the 95% CI for a hazard ratio of less than 1·30, closed testing for superiority was prespecified. This study is registered with ClinicalTrials.gov, number NCT02465515. Findings: Patients were screened between July 1, 2015, and Nov 24, 2016. 10 793 patients were screened and 9463 participants were enrolled and randomly assigned to groups: 4731 patients were assigned to receive albiglutide and 4732 patients to receive placebo. On Nov 8, 2017, it was determined that 611 primary endpoints and a median follow-up of at least 1·5 years had accrued, and participants returned for a final visit and discontinuation from study treatment; the last patient visit was on March 12, 2018. These 9463 patients, the intention-to-treat population, were evaluated for a median duration of 1·6 years and were assessed for the primary outcome. The primary composite outcome occurred in 338 (7%) of 4731 patients at an incidence rate of 4·6 events per 100 person-years in the albiglutide group and in 428 (9%) of 4732 patients at an incidence rate of 5·9 events per 100 person-years in the placebo group (hazard ratio 0·78, 95% CI 0·68–0·90), which indicated that albiglutide was superior to placebo (p&lt;0·0001 for non-inferiority; p=0·0006 for superiority). The incidence of acute pancreatitis (ten patients in the albiglutide group and seven patients in the placebo group), pancreatic cancer (six patients in the albiglutide group and five patients in the placebo group), medullary thyroid carcinoma (zero patients in both groups), and other serious adverse events did not differ between the two groups. There were three (&lt;1%) deaths in the placebo group that were assessed by investigators, who were masked to study drug assignment, to be treatment-related and two (&lt;1%) deaths in the albiglutide group. Interpretation: In patients with type 2 diabetes and cardiovascular disease, albiglutide was superior to placebo with respect to major adverse cardiovascular events. Evidence-based glucagon-like peptide 1 receptor agonists should therefore be considered as part of a comprehensive strategy to reduce the risk of cardiovascular events in patients with type 2 diabetes. Funding: GlaxoSmithKline

Ambient-noise tomography of the wider Vienna Basin region

We present a new 3-D shear-velocity model for the top 30 km of the crust in the wider Vienna Basin region based on surface waves extracted from ambient-noise cross-correlations. We use continuous seismic records of 63 broad-band stations of the AlpArray project to retrieve interstation Green’s functions from ambient-noise cross-correlations in the period range from 5 to 25 s. From these Green’s functions, we measure Rayleigh group traveltimes, utilizing all four components of the cross-correlation tensor, which are associated with Rayleigh waves (ZZ, RR, RZ and ZR), to exploit multiple measurements per station pair. A set of selection criteria is applied to ensure that we use high-quality recordings of fundamental Rayleigh modes. We regionalize the interstation group velocities in a 5 km × 5 km grid with an average path density of ∼20 paths per cell. From the resulting group-velocity maps, we extract local 1-D dispersion curves for each cell and invert all cells independently to retrieve the crustal shear-velocity structure of the study area. The resulting model provides a previously unachieved lateral resolution of seismic velocities in the region of ∼15 km. As major features, we image the Vienna Basin and Little Hungarian Plain as low-velocity anomalies, and the Bohemian Massif with high velocities. The edges of these features are marked with prominent velocity contrasts correlated with faults, such as the Alpine Front and Vienna Basin transfer fault system. The observed structures correlate well with surface geology, gravitational anomalies and the few known crystalline basement depths from boreholes. For depths larger than those reached by boreholes, the new model allows new insight into the complex structure of the Vienna Basin and surrounding areas, including deep low-velocity zones, which we image with previously unachieved detail. This model may be used in the future to interpret the deeper structures and tectonic evolution of the wider Vienna Basin region, evaluate natural resources, model wave propagation and improve earthquake locations, among others

Arrival angles of teleseismic fundamental mode Rayleigh waves across the AlpArray

The dense AlpArray network allows studying seismic wave propagation with high spatial resolution. Here we introduce an array approach to measure arrival angles of teleseismic Rayleigh waves. The approach combines the advantages of phase correlation as in the two-station method with array beamforming to obtain the phase-velocity vector. 20 earthquakes from the first two years of the AlpArray project are selected, and spatial patterns of arrival-angle deviations across the AlpArray are shown in maps, depending on period and earthquake location. The cause of these intriguing spatial patterns is discussed. A simple wave-propagation modelling example using an isolated anomaly and a Gaussian beam solution suggests that much of the complexity can be explained as a result of wave interference after passing a structural anomaly along the wave paths. This indicates that arrival-angle information constitutes useful additional information on the Earth structure, beyond what is currently used in inversions

Shear-wave velocity structure beneath the Dinarides from the inversion of Rayleigh-wave dispersion

Highlights • Rayleigh-wave phase velocity in the wider Dinarides region using the two-station method. • Uppermost mantle shear-wave velocity model of the Dinarides-Adriatic Sea region. • Velocity model reveals a robust high-velocity anomaly present under the whole Dinarides. • High-velocity anomaly reaches depth of 160 km in the northern Dinarides to more than 200 km under southern Dinarides. • New structural model incorporating delamination as one of the processes controlling the continental collision in the Dinarides. The interaction between the Adriatic microplate (Adria) and Eurasia is the main driving factor in the central Mediterranean tectonics. Their interplay has shaped the geodynamics of the whole region and formed several mountain belts including Alps, Dinarides and Apennines. Among these, Dinarides are the least investigated and little is known about the underlying geodynamic processes. There are numerous open questions about the current state of interaction between Adria and Eurasia under the Dinaric domain. One of the most interesting is the nature of lithospheric underthrusting of Adriatic plate, e.g. length of the slab or varying slab disposition along the orogen. Previous investigations have found a low-velocity zone in the uppermost mantle under the northern-central Dinarides which was interpreted as a slab gap. Conversely, several newer studies have indicated the presence of the continuous slab under the Dinarides with no trace of the low velocity zone. Thus, to investigate the Dinaric mantle structure further, we use regional-to-teleseismic surface-wave records from 98 seismic stations in the wider Dinarides region to create a 3D shear-wave velocity model. More precisely, a two-station method is used to extract Rayleigh-wave phase velocity while tomography and 1D inversion of the phase velocity are employed to map the depth dependent shear-wave velocity. Resulting velocity model reveals a robust high-velocity anomaly present under the whole Dinarides, reaching the depths of 160 km in the north to more than 200 km under southern Dinarides. These results do not agree with most of the previous investigations and show continuous underthrusting of the Adriatic lithosphere under Europe along the whole Dinaric region. The geometry of the down-going slab varies from the deeper slab in the north and south to the shallower underthrusting in the center. On-top of both north and south slabs there is a low-velocity wedge indicating lithospheric delamination which could explain the 200 km deep high-velocity body existing under the southern Dinarides

Crustal Thinning From Orogen to Back-Arc Basin: The Structure of the Pannonian Basin Region Revealed by P-to-S Converted Seismic Waves

We present the results of P-to-S receiver function analysis to improve the 3D image of the sedimentary layer, the upper crust, and lower crust in the Pannonian Basin area. The Pannonian Basin hosts deep sedimentary depocentres superimposed on a complex basement structure and it is surrounded by mountain belts. We processed waveforms from 221 three-component broadband seismological stations. As a result of the dense station coverage, we were able to achieve so far unprecedented spatial resolution in determining the velocity structure of the crust. We applied a three-fold quality control process; the first two being applied to the observed waveforms and the third to the calculated radial receiver functions. This work is the first comprehensive receiver function study of the entire region. To prepare the inversions, we performed station-wise H-Vp/Vs grid search, as well as Common Conversion Point migration. Our main focus was then the S-wave velocity structure of the area, which we determined by the Neighborhood Algorithm inversion method at each station, where data were sub-divided into back-azimuthal bundles based on similar Ps delay times. The 1D, nonlinear inversions provided the depth of the discontinuities, shear-wave velocities and Vp/Vs ratios of each layer per bundle, and we calculated uncertainty values for each of these parameters. We then developed a 3D interpolation method based on natural neighbor interpolation to obtain the 3D crustal structure from the local inversion results. We present the sedimentary thickness map, the first Conrad depth map and an improved, detailed Moho map, as well as the first upper and lower crustal thickness maps obtained from receiver function analysis. The velocity jump across the Conrad discontinuity is estimated at less than 0.2 km/s over most of the investigated area. We also compare the new Moho map from our approach to simple grid search results and prior knowledge from other techniques. Our Moho depth map presents local variations in the investigated area: the crust-mantle boundary is at 20–26 km beneath the sedimentary basins, while it is situated deeper below the Apuseni Mountains, Transdanubian and North Hungarian Ranges (28–33 km), and it is the deepest beneath the Eastern Alps and the Southern Carpathians (40–45 km). These values reflect well the Neogene evolution of the region, such as crustal thinning of the Pannonian Basin and orogenic thickening in the neighboring mountain belts

Shape Optimization Using the Finite Element Method on Multiple Meshes with Nitsche Coupling

An important step in shape optimization with partial differential equation constraints is to adapt the geometry during each optimization iteration. Common strategies are to employ mesh deformation or remeshing, where one or the other typically lacks robustness or is computationally expensive. This paper proposes a different approach, in which the computational domain is represented by multiple, independent nonmatching meshes. The individual meshes can move independently, hence mesh deformation or remeshing is entirely avoided if the geometry can be parameterized by a combination of rigid motions and scaling. For general geometry changes, we present a deformation scheme tailored to nonmatching meshes. This deformation scheme is robust because the nonmatching mesh interfaces are free to move and computationally cheap because the scheme is applied only on a subset of the meshes. To solve the state and corresponding adjoint equations we use the multimesh finite element method. This method weakly enforces continuity over the nonmatching mesh interfaces by using Nitsche and additional stability terms. To obtain the shape derivatives we analyze both the strong formulation (Hadamard formulation) and weak formulation (method of mappings). We demonstrate the capabilities of our approach on the optimal placement of heat emitting wires in a cable to minimize the chance of overheating, the drag minimization in Stokes flow, and the orientation of nine objects in Stokes flow.publishedVersio

Code for simulations in Chapter 2 of ResearchBrief

Contains the code for the simulations in Chapter 2 of ResearchBrief "Mesh dependence in PDE-constrained optimisation: An application in tidal turbine array layouts

Variational data assimilation for transient blood flow simulations - Cerebral aneurysms as an illustrative example

Several cardiovascular diseases are caused from localised abnormal blood flow such as in the case of stenosis or aneurysms. Prevailing theories propose that the development is caused by abnormal wall shear stress in focused areas. Computational fluid mechanics have arisen as a promising tool for a more precise and quantitative analysis, in particular because the anatomy is often readily available even by standard imaging techniques such as magnetic resonance and computed tomography angiography. However, computational fluid mechanics rely on accurate initial and boundary conditions, which are difficult to obtain. In this paper, we address the problem of recovering high‐resolution information from noisy and low‐resolution physical measurements of blood flow (for example, from phase‐contrast magnetic resonance imaging [PC‐MRI]) using variational data assimilation based on a transient Navier‐Stokes model. Numerical experiments are performed in both 3D (2D space and time) and 4D (3D space and time) and with pulsatile flow relevant for physiological flow in cerebral aneurysms. The results demonstrate that, with suitable regularisation, the model accurately reconstructs flow, even in the presence of significant noise

Apparent diffusion coefficient estimates based on 24 hours tracer movement support glymphatic transport in human cerebral cortex

The recently proposed glymphatic system suggests that bulk flow is important for clearing waste from the brain, and as such may underlie the development of e.g. Alzheimer’s disease. The glymphatic hypothesis is still controversial and several biomechanical modeling studies at the micro-level have questioned the system and its assumptions. In contrast, at the macro-level, there are many experimental findings in support of bulk flow. Here, we will investigate to what extent the CSF tracer distributions seen in novel magnetic resonance imaging (MRI) investigations over hours and days are suggestive of bulk flow as an additional component to diffusion. In order to include the complex geometry of the brain, the heterogeneous CSF flow around the brain, and the transport over the time-scale of days, we employed the methods of partial differential constrained optimization to identify the apparent diffusion coefficient (ADC) that would correspond best to the MRI findings. We found that the computed ADC in the cortical grey matter was 5–26% larger than the ADC estimated with DTI, which suggests that diffusion may not be the only mechanism governing transport