A novel 3-D segmentation algorithm for anatomic liver and tumor volume calculations for liver cancer treatment planning

Three-Dimensional (3-D) imaging is vital in computer-assisted surgical planning including minimal invasive surgery, targeted drug delivery, and tumor resection. Selective Internal Radiation Therapy (SIRT) is a liver directed radiation therapy for the treatment of liver cancer. Accurate calculation of anatomical liver and tumor volumes are essential for the determination of the tumor to normal liver ratio and for the calculation of the dose of Y-90 microspheres that will result in high concentration of the radiation in the tumor region as compared to nearby healthy tissue. Present manual techniques for segmentation of the liver from Computed Tomography (CT) tend to be tedious and greatly dependent on the skill of the technician/doctor performing the task. ^ This dissertation presents the development and implementation of a fully integrated algorithm for 3-D liver and tumor segmentation from tri-phase CT that yield highly accurate estimations of the respective volumes of the liver and tumor(s). The algorithm as designed requires minimal human intervention without compromising the accuracy of the segmentation results. Embedded within this algorithm is an effective method for extracting blood vessels that feed the tumor(s) in order to plan effectively the appropriate treatment. ^ Segmentation of the liver led to an accuracy in excess of 95% in estimating liver volumes in 20 datasets in comparison to the manual gold standard volumes. In a similar comparison, tumor segmentation exhibited an accuracy of 86% in estimating tumor(s) volume(s). Qualitative results of the blood vessel segmentation algorithm demonstrated the effectiveness of the algorithm in extracting and rendering the vasculature structure of the liver. Results of the parallel computing process, using a single workstation, showed a 78% gain. Also, statistical analysis carried out to determine if the manual initialization has any impact on the accuracy showed user initialization independence in the results. ^ The dissertation thus provides a complete 3-D solution towards liver cancer treatment planning with the opportunity to extract, visualize and quantify the needed statistics for liver cancer treatment. Since SIRT requires highly accurate calculation of the liver and tumor volumes, this new method provides an effective and computationally efficient process required of such challenging clinical requirements.^

Arc crust-magma interaction in the Andean Southern Volcanic Zone from thermobarometry, mineral composition, radiogenic isotope and rare earth element systematics of the Azufre-Planchon-Peteroa volcanic complex, Chile

The Andean Southern Volcanic Zone (SVZ) is a vast and complex continental arc that has been studied extensively to provide an understanding of arc-magma genesis, the origin and chemical evolution of the continental crust, and geochemical compositions of volcanic products. The present study focuses on distinguishing the magma/sub-arc crustal interaction of eruptive products from the Azufre-Planchon-Peteroa (APP 35°15'S) volcanic center and other major centers in the Central SVZ (CSVZ 37°S–42°S), Transitional SVZ (TSVZ 34.3–37.0°S), and Northern SVZ (NSVZ 33°S–34°30'S). New Hf and Nd isotopic and trace element data for SVZ centers are consistent with former studies that these magmas experienced variable depths of crystal fractionation, and that crustal assimilation is restricted to the lower crustal depths with an apparent role of garnet. Thermobarometric calculations applied to magma compositions constrain the depth of magma separation from mantle sources in all segments of the SVZ to(70-90 km). Magmatic separation at the APP complex occurs at an average depth of ~50 km which is confined to the mantle lithosphere and the base of the crust suggesting localized thermal abrasion both reservoirs. Thermobarometric calculations indicate that CSVZ primary magmas arise from a similar average depth of (~54 km) which confines magma separation to the asthenospheric mantle. The northwards along-arc Sr-Nd-Hf isotopic data and LREE enrichment accompanied with HREE depletion of SVZ mafic magmas correlates well with northward increasing crustal thickness and decreasing primary melt separation from mantle source regions indicating an increased involvement of lower crustal components in SVZ magma petrogenesis. ^ The study concludes that the development of mature subduction zones over millions of years of continuous magmatism requires that mafic arc derived melts stagnate at lower crustal levels due to density similarities and emplace at lower crustal depths. Basaltic underplating creates localized hot zone environments below major magmatic centers. These regions of high temperature/partial melting, and equilibration with underplated mafic rocks provides the mechanism that controls trace element and isotopic variability of primary magmas of the TSVZ and NSVZ from their baseline CSVZ-like precursors.^

Calculations on the end-group conformational barrier in carotenoids and a study of the bioavailability of xanthophyll esters in humans

The present study measures the increase in serum carotenoid concentration in 30 healthy individuals after supplementation with a low dose xanthophyll ester (3 and 6 mg of lutein equivalent/per day) when compared to a placebo. Serum levels of carotenoids were measured using HPLC and showed an increase in the concentration of lutein, zeaxanthin and four lutein metabolites proportional to dose. In order to further assess the importance of the end-group structure in carotenoids we have investigated the influence of the end-group type and functionality on the conformational energy barrier. We used the density functional method implemented on GAUSSIAN 98 to calculate the conformational energy curves for rotation of the P-ring or the E-ring relative to short polyene chains around the C6-C7 single bond. A large barrier is observed for the interconversion of conformers in the E-rings (8 kcal/mol) when compared to beta rings (2.3-3 kcal/mol).

(Table 1) Dissolved Hafnium and Neodymium in the South Atlantic measured on water bottle samples during POLARSTERN cruise ANT-XXIV/3

We present the first combined dissolved hafnium (Hf) and neodymium (Nd) concentrations and isotope compositions of deep water masses from the Atlantic sector of the Southern Ocean. Eight full depth profiles were analyzed for Hf and twelve for Nd. Hafnium concentrations are generally depleted in the upper few hundred meters ranging between 0.2 pmol/kg and 0.4 pmol/kg and increase to relatively constant values of around 0.6 pmol/kg in the deeper water column. At the stations north of the Polar Front (PF), Nd concentrations increase linearly from about 10 pmol/kg at depths of ~ 200 m to up to 31 pmol/kg close to the bottom indicating particle scavenging and release. Within the Weddell Gyre (WG), however, Nd concentrations are essentially constant at 25 pmol/kg at depths greater than ~ 1000 m. The distributions of both elements show a positive correlation with dissolved silicon implying a close linkage to diatom biogeochemistry. Hafnium essentially shows invariant isotope compositions with values averaging at epsilon-Hf = +4.6, whereas Nd isotopes mark distinct differences between water masses, such as modified North Atlantic Deep Water (NADW, epsilon-Nd = -11 to -10) and Antarctic Bottom Water (AABW, epsilon-Nd = -8.6 to -9.6), but also waters locally advected via the Agulhas Current can be identified by their unradiogenic Nd isotope compositions. Mixing calculations suggest that a small fraction of Nd is removed by particle scavenging during mixing of water masses north of the PF. Nevertheless, the Nd isotope composition has apparently not been significantly affected by uptake and release of Nd from particles, as indicated by mixing calculations. A mixing envelope of an approximated North Pacific and a North Atlantic end-member shows that Nd isotope and concentration patterns in the Lower Circumpolar Deep Water (LCDW) can be fully explained by ~ 30:70 percentage contributions of these respective end-members.

(Table 1, page 641) Hf isotope time series for the uppermost part of a ferromanganese crust from the San Pablo ferro-manganese pavement

The Hf isotope composition of seawater does not match that expected from dissolution of bulk continental crust. This mismatch is generally considered to be due to retention of unradiogenic Hf in resistant zircons during incomplete weathering of continental crust. During periods of intense glacial weathering, zircons should break down more efficiently, resulting in the release of highly unradiogenic Hf to the oceans. We test this hypothesis by comparing Nd and Hf isotope time series obtained from NW Atlantic ferromanganese crusts. Both isotope systems show a decrease associated with the onset of northern hemisphere glaciation. The observed changes display distinct trajectories in epsilon Nd- epsilon Hf space, which differ from previously reported arrays of bulk terrestrial material and seawater. Such patterns are consistent with the release of highly unradiogenic Hf from very old zircons, facilitated by enhanced mechanical weathering.

Neodymium and Strontium chemistry of Arctic Ocean sediments

The early oceanographic history of the Arctic Ocean is important in regulating, and responding to, climatic changes. However, constraints on its oceanographic history preceding the Quaternary (the past 1.8 Myr) have become available only recently, because of the difficulties associated with obtaining continuous sediment records in such a hostile setting. Here, we use the neodymium isotope compositions of two sediment cores recovered near the North Pole to reconstruct over the past ~5 Myr the sources contributing to Arctic Intermediate Water, a water mass found today at depths of 200 to 1,500 m. We interpret high neodymium ratios for the period between 15 and 2 Myr ago, and for the glacial periods thereafter, as indicative of weathering input from the Siberian Putoranan basalts into the Arctic Ocean. Arctic Intermediate Water was then derived from brine formation in the Eurasian shelf regions, with only a limited contribution of intermediate water from the North Atlantic. In contrast, the modern circulation pattern, with relatively high contributions of North Atlantic Intermediate Water and negligible input from brine formation, exhibits low neodymium isotope ratios and is typical for the interglacial periods of the past 2 Myr. We suggest that changes in climatic conditions and the tectonic setting were responsible for switches between these two modes.

(Table 1) Sample locations, Hf and Nd isotopic compositions and concentrations of Atlantic seawater samples with hydrographic details

The first full water column hafnium isotopic compositions of Atlantic seawater have been obtained at seven locations from the Labrador Sea to the Drake Passage. Despite subpicomolar concentrations in seawater, a precision of the Hf isotopic measurements of <0.7 epsilon-Hf units was achieved. An overall epsilon-Hf range between -3.1 in the Labrador Sea and +4.4 in Antarctic bottom water was determined, the distribution of which broadly reflects continental weathering inputs. Within particular water column profiles, significant differences of up to 4 epsilon-Hf units occur. Combined with Nd isotope data of the same samples, it is evident that the Hf isotopic composition of seawater is too radiogenic for a given Nd isotopic composition and that the largest difference between expected and measured Hf isotopic compositions in seawater occurs near the oldest continental crust in the Labrador Sea. This corroborates the previous proposition, which was mainly based on ferromanganese crust data, that the Hf isotopic composition of seawater is controlled by incongruent weathering of continental crust and possibly, to some extent, by hydrothermal contributions. Hafnium concentrations in the ocean do not increase along the deep ocean conveyer indicating an oceanic residence time of only a few hundred years, which is significantly shorter than previously assumed. The Hf isotopic composition of past seawater can therefore serve as a proxy for short distance, basin scale mixing processes and the regime and intensity of nearby continental weathering processes.

Generalized atoms in molecules approach (GAIM approach)

This thesis involves two parts. The first is a new-proposed theoretical approach called generalized atoms in molecules (GAIM). The second is a computational study on the deamination reaction of adenine with OH⁻/nH₂O (n=0, 1, 2, 3) and 3H₂O. The GAIM approach aims to solve the energy of each atom variationally in the first step and then to build the energy of a molecule from each atom. Thus the energy of a diatomic molecule (A-B) is formulated as a sum of its atomic energies, EA and EB. Each of these atomic energies is expressed as, EA = Hᴬ + Vₑₑᴬᴬ + 1/2Vₑₑᴬ<>ᴮ EB = Hᴮ + Vₑₑᴮᴮ + 1/2Vₑₑᴬ<>ᴮ where; Hᴬ and Hᴮ are the kinetic and nuclear attraction energy of electrons of atoms A and B, respectively; Vₑₑᴬᴬ and Vₑₑᴮᴮ are the interaction energy between the electrons on atoms A and B, respectively; and Vₑₑᴬ<>ᴮ is the interaction energy between the electrons of atom A with the electrons of atom B. The energy of the molecule is then minimized subject to the following constraint, |ρA(r)dr + |ρB(r)dr = N where ρA(r) and ρB(r) are the electron densities of atoms A and B, respectively, and N is the number of electrons. The initial testing of the performance of GAIM was done through calculating dissociation curves for H₂, LiH, Li₂, BH, HF, HCl, N₂, F₂, and Cl₂. The numerical results show that GAIM performs very well with H₂, LiH, Li₂, BH, HF, and HCl. GAIM shows convergence problems with N₂, F₂, and Cl₂ due to difficulties in reordering the degenerate atomic orbitals Pₓ, Py, and Pz in N, F, and Cl atoms. Further work for the development of GAIM is required. Deamination of adenine results in one of several forms of premutagenic lesions occurring in DNA. In this thesis, mechanisms for the deamination reaction of adenine with OH⁻/nH₂O, (n = 0, 1, 2, 3) and 3H₂O were investigated. HF/6-31G(d), B3LYP/6-31G(d), MP2/6-31G(d), and B3LYP/6-31+G(d) levels of theory were employed to optimize all the geometries. Energies were calculated at the G3MP2B3 and CBS-QB3 levels of theory. The effect of solvent (water) was computed using the polarizable continuum model (PCM). Intrinsic reaction coordinate (IRC) calculations were performed for all transition states. Five pathways were investigated for the deamination reaction of adenine with OH⁻/nH₂O and 3H₂O. The first four pathways (A-D) begin with by deprotonation at the amino group of adenine by OH⁻, while pathway E is initiated by tautomerization of adenine. For all pathways, the next two steps involve the formation of a tetrahedral intermediate followed by dissociation to yield products via a 1,3-hydrogen shift. Deamination with a single OH⁻ has a high activation barrier (190 kJ mol⁻¹ using G3MP2B3 level) for the rate-determining step. Addition of one water molecule reduces this barrier by 68 kJ mol⁻¹ calculated at G3MP2B3 level. Adding more water molecules decreases the overall activation energy of the reaction, but the effect becomes smaller with each additional water molecule. The most plausible mechanism is pathway E, the deamination reaction of adenine with 3H₂O, which has an overall G3MP2B3 activation energy of 139 and 137 kJ mol⁻¹ in the gas phase and PCM, respectively. This barrier is lower than that for the deamination with OH⁻/3H₂O by 6 and 2 kJ mol⁻¹ in the gas phase and PCM, respectively.

Radiogenic Hf isotopic composition of continental eolian dust in Pacific Ocean sediments

The inorganic silicate fraction extracted from bulk pelagic sediments from the North Pacific Ocean is eolian dust. It monitors the composition of continental crust exposed to erosion in Asia. 176Lu/177Hf ratios of modern dust are subchondritic between 0.011 and 0.016 but slightly elevated with respect to immature sediments. Modern dust samples display a large range in Hf isotopic composition (IC), -4.70 < epsilon-Hf < +16.45, which encompasses that observed for the time series of DSDP cores 885/886 and piston core LL44-GPC3 extending back to the late Cretaceous. Hafnium and neodymium isotopic results are consistent with a dominantly binary mixture of dust contributed from island arc volcanic material and dust from central Asia. The Hf-Nd isotopic correlation for all modern dust samples, epsilon-Hf= =0.78 epsilon-Nd = +5.66 (n =22, R**2 =0.79), is flatter than those reported so far for terrestrial reservoirs. Moreover, the variability in epsilon-Hf of Asian dust exceeds that predicted on the basis of corresponding epsilon-Nd values (34.76 epsilon-Hf < +2.5; -10.96< epsilon-Nd <-10.1). This is attributed to: (1) the fixing of an important unradiogenic fraction of Hf in zircons, balanced by radiogenic Hf that is mobile in the erosional cycle, (2) the elevated Lu/Hf ratio in chemical sediments which, given time, results in a Hf signature that is radiogenic compared with Hf expected from its corresponding Nd isotopic components, and (3) the possibility that diagenetic resetting of marine sediments may incorporate a significant radiogenic Hf component into diagenetically grown minerals such as illite. Together, these processes may explain the variability and more radiogenic character of Hf isotopes when compared to the Nd isotopic signatures of Asian dust. The Hf-Nd isotope time series of eolian dust are consistent with the results of modern dust except two samples that have extremely radiogenic Hf for their Nd (epsilon-Hf =+8.6 and +10.3, epsilon-Nd =39.5 and 39.8). These data may point to a source contribution of dust unresolved by Nd and Pb isotopes. The Hf IC of eolian dust input to the oceans may be more variable and more radiogenic than previously anticipated. The Hf signature of Pacific seawater, however, has varied little over the past 20 Myr, especially across the drastic increase of eolian dust flux from Asia around 3.5 Ma. Therefore, continental contributions to seawater Hf appear to be riverine rather than eolian. Current predictions regarding the relative proportions of source components to seawater Hf must account for the presence of a variable and radiogenic continental component. Data on the IC and flux of river-dissolved Hf to the oceans are urgently required to better estimate contributions to seawater Hf. This then would permit the use of Hf isotopes as a monitor of past changes in erosion.

A new SOS-DFPT approximation for NMR shielding calculations : the Loc.3 correction applied to the catalytic mechanism of Serine Proteases

Thèse numérisée par la Direction des bibliothèques de l'Université de Montréal.

Path Optimization in Free Energy Calculations

Free energy calculations are a computational method for determining thermodynamic quantities, such as free energies of binding, via simulation.

Currently, due to computational and algorithmic limitations, free energy calculations are limited in scope.

In this work, we propose two methods for improving the efficiency of free energy calculations.

First, we expand the state space of alchemical intermediates, and show that this expansion enables us to calculate free energies along lower variance paths.

We use Q-learning, a reinforcement learning technique, to discover and optimize paths at low computational cost.

Second, we reduce the cost of sampling along a given path by using sequential Monte Carlo samplers.

We develop a new free energy estimator, pCrooks (pairwise Crooks), a variant on the Crooks fluctuation theorem (CFT), which enables decomposition of the variance of the free energy estimate for discrete paths, while retaining beneficial characteristics of CFT.

Combining these two advancements, we show that for some test models, optimal expanded-space paths have a nearly 80% reduction in variance relative to the standard path.

Additionally, our free energy estimator converges at a more consistent rate and on average 1.8 times faster when we enable path searching, even when the cost of path discovery and refinement is considered.

A new SOS-DFPT approximation for NMR shielding calculations : the Loc.3 correction applied to the catalytic mechanism of Serine Proteases

Thèse numérisée par la Direction des bibliothèques de l'Université de Montréal.

Alkenone stable carbon isotopes, carbonate stable carbon and oxygen isotopes, planktic foraminifera boron isotopes and atmospheric CO2 calculations and box model results from Ras il-Pellegrin section, Malta

The development of a permanent, stable ice sheet in East Antarctica happened during the middle Miocene, about 14 million years (Myr) ago. The middle Miocene therefore represents one of the distinct phases of rapid change in the transition from the "greenhouse" of the early Eocene to the "icehouse" of the present day. Carbonate carbon isotope records of the period immediately following the main stage of ice sheet development reveal a major perturbation in the carbon system, represented by the positive d13C excursion known as carbon maximum 6 ("M6"), which has traditionally been interpreted as reflecting increased burial of organic matter and atmospheric pCO2 drawdown. More recently, it has been suggested that the d13C excursion records a negative feedback resulting from the reduction of silicate weathering and an increase in atmospheric pCO2. Here we present high-resolution multi-proxy (alkenone carbon and foraminiferal boron isotope) records of atmospheric carbon dioxide and sea surface temperature across CM6. Similar to previously published records spanning this interval, our records document a world of generally low (~300 ppm) atmospheric pCO2 at a time generally accepted to be much warmer than today. Crucially, they also reveal a pCO2 decrease with associated cooling, which demonstrates that the carbon burial hypothesis for CM6 is feasible and could have acted as a positive feedback on global cooling.

Cost-effectiveness of natriuretic peptide-based screening and collaborative care: a report from the STOP-HF (St Vincent's Screening TO Prevent Heart Failure) study

AIMS: Prevention of cardiovascular disease and heart failure (HF) in a cost-effective manner is a public health goal. This work aims to assess the cost-effectiveness of the St Vincent's Screening TO Prevent Heart Failure (STOP-HF) intervention.

METHODS AND RESULTS: This is a substudy of 1054 participants with cardiovascular risk factors [median age 65.8 years, interquartile range (IQR) 57.8:72.4, with 4.3 years, IQR 3.4:5.2, follow-up]. Annual natriuretic peptide-based screening was performed, with collaborative cardiovascular care between specialist physicians and general practitioners provided to patients with BNP levels >50 pg/mL. Analysis of cost per case prevented and cost-effectiveness per quality-adjusted life year (QALY) gained was performed. The primary clinical endpoint of LV dysfunction (LVD) with or without HF was reduced in intervention patients [odds ratio (OR) 0.60; 95% confidence interval (CI) 0.38-0.94; P = 0.026]. There were 157 deaths and/or emergency hospitalizations for major adverse cardiac events (MACE) in the control group vs. 102 in the intervention group (OR 0.68; 95% CI 0.49-0.93; P = 0.01). The cost per case of LVD/HF prevented was €9683 (sensitivity range -€843 to €20 210), whereas the cost per MACE prevented was €3471 (sensitivity range -€302 to €7245). Cardiovascular hospitalization savings offset increased outpatient and primary care costs. The cost per QALY gain was €1104 and the intervention has an 88% probability of being cost-effective at a willingness to pay threshold of €30 000.

CONCLUSION: Among patients with cardiovascular risk factors, natriuretic peptide-based screening and collaborative care reduced LVD, HF, and MACE, and has a high probability of being cost-effective.

TRIAL REGISTRATION: NCT00921960.

Natriuretic peptide-based screening and collaborative care for heart failure:the STOP-HF randomized trial

IMPORTANCE: Prevention strategies for heart failure are needed.

OBJECTIVE: To determine the efficacy of a screening program using brain-type natriuretic peptide (BNP) and collaborative care in an at-risk population in reducing newly diagnosed heart failure and prevalence of significant left ventricular (LV) systolic and/or diastolic dysfunction.

DESIGN, SETTING, AND PARTICIPANTS: The St Vincent's Screening to Prevent Heart Failure Study, a parallel-group randomized trial involving 1374 participants with cardiovascular risk factors (mean age, 64.8 [SD, 10.2] years) recruited from 39 primary care practices in Ireland between January 2005 and December 2009 and followed up until December 2011 (mean follow-up, 4.2 [SD, 1.2] years).

INTERVENTION: Patients were randomly assigned to receive usual primary care (control condition; n=677) or screening with BNP testing (n=697). Intervention-group participants with BNP levels of 50 pg/mL or higher underwent echocardiography and collaborative care between their primary care physician and specialist cardiovascular service.

MAIN OUTCOMES AND MEASURES: The primary end point was prevalence of asymptomatic LV dysfunction with or without newly diagnosed heart failure. Secondary end points included emergency hospitalization for arrhythmia, transient ischemic attack, stroke, myocardial infarction, peripheral or pulmonary thrombosis/embolus, or heart failure.

RESULTS: A total of 263 patients (41.6%) in the intervention group had at least 1 BNP reading of 50 pg/mL or higher. The intervention group underwent more cardiovascular investigations (control, 496 per 1000 patient-years vs intervention, 850 per 1000 patient-years; incidence rate ratio, 1.71; 95% CI, 1.61-1.83; P<.001) and received more renin-angiotensin-aldosterone system-based therapy at follow-up (control, 49.6%; intervention, 56.5%; P=.01). The primary end point of LV dysfunction with or without heart failure was met in 59 (8.7%) of 677 in the control group and 37 (5.3%) of 697 in the intervention group (odds ratio [OR], 0.55; 95% CI, 0.37-0.82; P = .003). Asymptomatic LV dysfunction was found in 45 (6.6%) of 677 control-group patients and 30 (4.3%) of 697 intervention-group patients (OR, 0.57; 95% CI, 0.37-0.88; P = .01). Heart failure occurred in 14 (2.1%) of 677 control-group patients and 7 (1.0%) of 697 intervention-group patients (OR, 0.48; 95% CI, 0.20-1.20; P = .12). The incidence rates of emergency hospitalization for major cardiovascular events were 40.4 per 1000 patient-years in the control group vs 22.3 per 1000 patient-years in the intervention group (incidence rate ratio, 0.60; 95% CI, 0.45-0.81; P = .002).

CONCLUSION AND RELEVANCE: Among patients at risk of heart failure, BNP-based screening and collaborative care reduced the combined rates of LV systolic dysfunction, diastolic dysfunction, and heart failure.

TRIAL REGISTRATION: clinicaltrials.gov Identifier: NCT00921960.