Synthesis of a multimodal molecular imaging probe based on a hyperbranched polymer architecture

Molecular imaging is utilised in modern medicine to aid in the diagnosis and treatment of disease by allowing its spatiotemporal state to be examined in vivo. This study focuses on the development of novel multimodal molecular imaging agents based on hyperbranched polymers that combine the complementary capabilities of optical fluorescence imaging and positron emission tomography-computed tomography (PET/CT) into one construct. RAFT-mediated polymerisation was used to prepare two hydrophilic hyperbranched polymers that were differentiated by their size and level of branching. The multiple functional end-groups facilitated covalent attachment of both near infrared fluorescent dyes for optical imaging, as well as a copper chelator allowing binding of 64Cu as a PET radio nuclei. In vivo multimodal imaging of mice using PET/CT and planar optical imaging was first used to assess the biodistribution of the polymeric materials and it was shown that the larger and more branched polymer had a significantly longer circulation time. The larger constructs were also shown to exhibit enhanced accumulation in solid tumours in a murine B16 melanoma model. Importantly, it was demonstrated that the PET modality gave rise to high sensitivity immediately after injection of the agent, while the optical modality facilitated extended longitudinal studies, thus highlighting how the complementary capabilities of the molecular imaging agents can be useful for studying various diseases, including cancer.

A mouse model for spondyloepiphyseal dysplasia congenita with secondary osteoarthritis due to a Col2a1 mutation

Progeny of mice treated with the mutagen N-ethyl-N-nitrosourea (ENU) revealed a mouse, designated Longpockets (Lpk), with short humeri, abnormal vertebrae, and disorganized growth plates, features consistent with spondyloepiphyseal dysplasia congenita (SEDC). The Lpk phenotype was inherited as an autosomal dominant trait. Lpk/+ mice were viable and fertile and Lpk/Lpk mice died perinatally. Lpk was mapped to chromosome 15 and mutational analysis of likely candidates from the interval revealed a Col2a1 missense Ser1386Pro mutation. Transient transfection of wild-type and Ser1386Pro mutant Col2a1 c-Myc constructs in COS-7 cells and CH8 chondrocytes demonstrated abnormal processing and endoplasmic reticulum retention of the mutant protein. Histology revealed growth plate disorganization in 14-day-old Lpk/+ mice and embryonic cartilage from Lpk/+ and Lpk/Lpk mice had reduced safranin-O and type-II collagen staining in the extracellular matrix. The wild-type and Lpk/+ embryos had vertical columns of proliferating chondrocytes, whereas those in Lpk/Lpk mice were perpendicular to the direction of bone growth. Electron microscopy of cartilage from 18.5 dpc wild-type, Lpk/+, and Lpk/Lpk embryos revealed fewer and less elaborate collagen fibrils in the mutants, with enlarged vacuoles in the endoplasmic reticulum that contained amorphous inclusions. Micro-computed tomography (CT) scans of 12-week-old Lpk/+ mice revealed them to have decreased bone mineral density, and total bone volume, with erosions and osteophytes at the joints. Thus, an ENU mouse model with a Ser1386Pro mutation of the Col2a1 C-propeptide domain that results in abnormal collagen processing and phenotypic features consistent with SEDC and secondary osteoarthritis has been established.

Impact of coronary tortuosity on coronary pressure: Numerical simulation study

Background: Coronary tortuosity (CT) is a common coronary angiographic finding. Whether CT leads to an apparent reduction in coronary pressure distal to the tortuous segment of the coronary artery is still unknown. The purpose of this study is to determine the impact of CT on coronary pressure distribution by numerical simulation. Methods: 21 idealized models were created to investigate the influence of coronary tortuosity angle (CTA) and coronary tortuosity number (CTN) on coronary pressure distribution. A 2D incompressible Newtonian flow was assumed and the computational simulation was performed using finite volume method. CTA of 30°, 60°, 90°, 120° and CTN of 0, 1, 2, 3, 4, 5 were discussed under both steady and pulsatile conditions, and the changes of outlet pressure and inlet velocity during the cardiac cycle were considered. Results: Coronary pressure distribution was affected both by CTA and CTN. We found that the pressure drop between the start and the end of the CT segment decreased with CTA, and the length of the CT segment also declined with CTA. An increase in CTN resulted in an increase in the pressure drop. Conclusions: Compared to no-CT, CT can results in more decrease of coronary blood pressure in dependence on the severity of tortuosity and severe CT may cause myocardial ischemia.

Image-based midsole insert design and the material effects on heel plantar pressure distribution during simulated walking loads

The primary objective of this paper is to study the use of medical image-based finite element (FE) modelling in subjectspecific midsole design and optimisation for heel pressure reduction using a midsole plug under the calcaneus area (UCA). Plugs with different relative dimensions to the size of the calcaneus of the subject have been incorporated in the heel region of the midsole. The FE foot model was validated by comparing the numerically predicted plantar pressure with biomechanical tests conducted on the same subject. For each UCA midsole plug design, the effect of material properties and plug thicknesses on the plantar pressure distribution and peak pressure level during the heel strike phase of normal walking was systematically studied. The results showed that the UCA midsole insert could effectively modify the pressure distribution, and its effect is directly associated with the ratio of the plug dimension to the size of the calcaneus bone of the subject. A medium hardness plug with a size of 95% of the calcaneus has achieved the best performance for relieving the peak pressure in comparison with the pressure level for a solid midsole without a plug, whereas a smaller plug with a size of 65% of the calcaneus insert with a very soft material showed minimum beneficial effect for the pressure relief.

Computed wall stress may predict the growth of abdominal aortic aneurysm

Growth rate of abdominal aortic aneurysm (AAA) is thought to be an important indicator of the potential risk of rupture. Wall stress is also thought to be a trigger for its rupture. However, stress change during the expansion of an AAA is unclear. Forty-four patients with AAAs were included in this longitudinal follow-up study. They were assessed by serial abdominal ultrasonography and computerized tomography (CT) scans if a critical size was reached or a rapid expansion occurred. Patient-specific 3-dimensional AAA geometries were reconstructed from the follow-up CT images. Structural analysis was performed to calculate the wall stresses of the AAA models at both baseline and final visit. A non-linear large-strain finite element method was used to compute the wall stress distribution. The average growth rate was 0.66cm/year (range 0-1.32 cm/year). A significantly positive correlation between shoulder tress at baseline and growth rate was found (r=0.342; p=0.02). A higher shoulder stress is associated with a rapidly expanding AAA. Therefore, it may be useful for estimating the growth expansion of AAAs and further risk stratification of patients with AAAs.

CBCT image guidance in head and neck irradiation: The impact of daily and weekly imaging protocols

Purpose This study evaluated the impact of a daily and weekly image-guided radiotherapy protocols in reducing setup errors and setting of appropriate margins in head and neck cancer patients. Materials and methods Interfraction and systematic shifts for the hypothetical day 1–3 plus weekly imaging were extrapolated from daily imaging data from 31 patients (964 cone beam computed tomography (CBCT) scans). In addition, residual setup errors were calculated by taking the average shifts in each direction for each patient based on the first three shifts and were presumed to represent systematic setup error. The clinical target volume (CTV) to planning target volume (PTV) margins were calculated using van Herk formula and analysed for each protocol. Results The mean interfraction shifts for daily imaging were 0·8, 0·3 and 0·5 mm in the S-I (superior-inferior), L-R (left-right) and A-P (anterior-posterior) direction, respectively. On the other hand the mean shifts for day 1–3 plus weekly imaging were 0·9, 1·8 and 0·5 mm in the S-I, L-R and A-P direction, respectively. The mean day 1–3 residual shifts were 1·5, 2·1 and 0·7 mm in the S-I, L-R and A-P direction, respectively. No significant difference was found in the mean setup error for the daily and hypothetical day 1–3 plus weekly protocol. However, the calculated CTV to PTV margins for the daily interfraction imaging data were 1·6, 3·8 and 1·4 mm in the S-I, L-R and A-P directions, respectively. Hypothetical day 1–3 plus weekly resulted in CTV–PTV margins of 5, 4·2 and 5 mm in the S-I, L-R and A-P direction. Conclusions The results of this study show that a daily CBCT protocol reduces setup errors and allows setup margin reduction in head and neck radiotherapy compared to a weekly imaging protocol.

Delineating environmental envelopes to improve mapping of species distributions, via a hurdle model with CART &/or MaxEnT

Species distribution modelling (SDM) typically analyses species’ presence together with some form of absence information. Ideally absences comprise observations or are inferred from comprehensive sampling. When such information is not available, then pseudo-absences are often generated from the background locations within the study region of interest containing the presences, or else absence is implied through the comparison of presences to the whole study region, e.g. as is the case in Maximum Entropy (MaxEnt) or Poisson point process modelling. However, the choice of which absence information to include can be both challenging and highly influential on SDM predictions (e.g. Oksanen and Minchin, 2002). In practice, the use of pseudo- or implied absences often leads to an imbalance where absences far outnumber presences. This leaves analysis highly susceptible to ‘naughty-noughts’: absences that occur beyond the envelope of the species, which can exert strong influence on the model and its predictions (Austin and Meyers, 1996). Also known as ‘excess zeros’, naughty noughts can be estimated via an overall proportion in simple hurdle or mixture models (Martin et al., 2005). However, absences, especially those that occur beyond the species envelope, can often be more diverse than presences. Here we consider an extension to excess zero models. The two-staged approach first exploits the compartmentalisation provided by classification trees (CTs) (as in O’Leary, 2008) to identify multiple sources of naughty noughts and simultaneously delineate several species envelopes. Then SDMs can be fit separately within each envelope, and for this stage, we examine both CTs (as in Falk et al., 2014) and the popular MaxEnt (Elith et al., 2006). We introduce a wider range of model performance measures to improve treatment of naughty noughts in SDM. We retain an overall measure of model performance, the area under the curve (AUC) of the Receiver-Operating Curve (ROC), but focus on its constituent measures of false negative rate (FNR) and false positive rate (FPR), and how these relate to the threshold in the predicted probability of presence that delimits predicted presence from absence. We also propose error rates more relevant to users of predictions: false omission rate (FOR), the chance that a predicted absence corresponds to (and hence wastes) an observed presence, and the false discovery rate (FDR), reflecting those predicted (or potential) presences that correspond to absence. A high FDR may be desirable since it could help target future search efforts, whereas zero or low FOR is desirable since it indicates none of the (often valuable) presences have been ignored in the SDM. For illustration, we chose Bradypus variegatus, a species that has previously been published as an exemplar species for MaxEnt, proposed by Phillips et al. (2006). We used CTs to increasingly refine the species envelope, starting with the whole study region (E0), eliminating more and more potential naughty noughts (E1–E3). When combined with an SDM fit within the species envelope, the best CT SDM had similar AUC and FPR to the best MaxEnt SDM, but otherwise performed better. The FNR and FOR were greatly reduced, suggesting that CTs handle absences better. Interestingly, MaxEnt predictions showed low discriminatory performance, with the most common predicted probability of presence being in the same range (0.00-0.20) for both true absences and presences. In summary, this example shows that SDMs can be improved by introducing an initial hurdle to identify naughty noughts and partition the envelope before applying SDMs. This improvement was barely detectable via AUC and FPR yet visible in FOR, FNR, and the comparison of predicted probability of presence distribution for pres/absence.

Automatic ICD-10 classification of cancers from free-text death certificates

Objective Death certificates provide an invaluable source for cancer mortality statistics; however, this value can only be realised if accurate, quantitative data can be extracted from certificates – an aim hampered by both the volume and variable nature of certificates written in natural language. This paper proposes an automatic classification system for identifying cancer related causes of death from death certificates. Methods Detailed features, including terms, n-grams and SNOMED CT concepts were extracted from a collection of 447,336 death certificates. These features were used to train Support Vector Machine classifiers (one classifier for each cancer type). The classifiers were deployed in a cascaded architecture: the first level identified the presence of cancer (i.e., binary cancer/nocancer) and the second level identified the type of cancer (according to the ICD-10 classification system). A held-out test set was used to evaluate the effectiveness of the classifiers according to precision, recall and F-measure. In addition, detailed feature analysis was performed to reveal the characteristics of a successful cancer classification model. Results The system was highly effective at identifying cancer as the underlying cause of death (F-measure 0.94). The system was also effective at determining the type of cancer for common cancers (F-measure 0.7). Rare cancers, for which there was little training data, were difficult to classify accurately (F-measure 0.12). Factors influencing performance were the amount of training data and certain ambiguous cancers (e.g., those in the stomach region). The feature analysis revealed a combination of features were important for cancer type classification, with SNOMED CT concept and oncology specific morphology features proving the most valuable. Conclusion The system proposed in this study provides automatic identification and characterisation of cancers from large collections of free-text death certificates. This allows organisations such as Cancer Registries to monitor and report on cancer mortality in a timely and accurate manner. In addition, the methods and findings are generally applicable beyond cancer classification and to other sources of medical text besides death certificates.

Solutions for archiving data in long-term studies: A reply to Whitlock et al.

In our recent paper [1], we discussed some potential undesirable consequences of public data archiving (PDA) with specific reference to long-term studies and proposed solutions to manage these issues. We reaffirm our commitment to data sharing and collaboration, both of which have been common and fruitful practices supported for many decades by researchers involved in long-term studies. We acknowledge the potential benefits of PDA (e.g., [2]), but believe that several potential negative consequences for science have been underestimated [1] (see also 3 and 4). The objective of our recent paper [1] was to define practices to simultaneously maximize the benefits and minimize the potential unwanted consequences of PDA.

G-protein beta3 subunit gene (GNB3) variant in causation of essential hypertension

-Essential hypertensives display enhanced signal transduction through pertussis toxin-sensitive G proteins. The T allele of a C825T variant in exon 10 of the G protein beta3 subunit gene (GNB3) induces formation of a splice variant (Gbeta3-s) with enhanced activity. The T allele of GNB3 was shown recently to be associated with hypertension in unselected German patients (frequency=0.31 versus 0.25 in control). To confirm and extend this finding in a different setting, we performed an association study in Australian white hypertensives. This involved an extensively examined cohort of 110 hypertensives, each of whom were the offspring of 2 hypertensive parents, and 189 normotensives whose parents were both normotensive beyond age 50 years. Genotyping was performed by polymerase chain reaction and digestion with BseDI, which either cut (C allele) or did not cut (T allele) the 268-bp polymerase chain reaction product. T allele frequency in the hypertensive group was 0.43 compared with 0.25 in the normotensive group (chi2=22; P=0.00002; odds ratio=2.3; 95% CI=1.7 to 3.3). The T allele tracked with higher pretreatment blood pressure: diastolic=105+/-7, 109+/-16, and 128+/-28 mm Hg (mean+/-SD) for CC, CT, and TT, respectively (P=0.001 by 1-way ANOVA). Blood pressures were higher in female hypertensives with a T allele (P=0.006 for systolic and 0.0003 for diastolic by ANOVA) than they were in male hypertensives. In conclusion, the present study of a group with strong family history supports a role for a genetically determined, physiologically active splice variant of the G protein beta3 subunit gene in the causation of essential hypertension.

The inferior whorl for detecting diabetic peripheral neuropathy using corneal confocal microscopy

PURPOSE: In vivo corneal confocal microscopy (CCM) is increasingly used as a surrogate endpoint in studies of diabetic polyneuropathy (DPN). However, it is not clear whether imaging the central cornea provides optimal diagnostic utility for DPN. Therefore, we compared nerve morphology in the central cornea and the inferior whorl, a more distal and densely innervated area located inferior and nasal to the central cornea. METHODS: A total of 53 subjects with type 1/type 2 diabetes and 15 age-matched control subjects underwent detailed assessment of neuropathic symptoms (NPS), deficits (neuropathy disability score [NDS]), quantitative sensory testing (vibration perception threshold [VPT], cold and warm threshold [CT/WT], and cold- and heat-induced pain [CIP/HIP]), and electrophysiology (sural and peroneal nerve conduction velocity [SSNCV/PMNCV], and sural and peroneal nerve amplitude [SSNA/PMNA]) to diagnose patients with (DPN+) and without (DPN-) neuropathy. Corneal nerve fiber density (CNFD) and length (CNFL) in the central cornea, and inferior whorl length (IWL) were quantified. RESULTS: Comparing control subjects to DPN- and DPN+ patients, there was a significant increase in NDS (0 vs. 2.6 ± 2.3 vs. 3.3 ± 2.7, P < 0.01), VPT (V; 5.4 ± 3.0 vs. 10.6 ± 10.3 vs. 17.7 ± 11.8, P < 0.01), WT (°C; 37.7 ± 3.5 vs. 39.1 ± 5.1 vs. 41.7 ± 4.7, P < 0.05), and a significant decrease in SSNCV (m/s; 50.2 ± 5.4 vs. 48.4 ± 5.0 vs. 39.5 ± 10.6, P < 0.05), CNFD (fibers/mm2; 37.8 ± 4.9 vs. 29.7 ± 7.7 vs. 27.1 ± 9.9, P < 0.01), CNFL (mm/mm2; 27.5 ± 3.6 vs. 24.4 ± 7.8 vs. 20.7 ± 7.1, P < 0.01), and IWL (mm/mm2; 35.1 ± 6.5 vs. 26.2 ± 10.5 vs. 23.6 ± 11.4, P < 0.05). For the diagnosis of DPN, CNFD, CNFL, and IWL achieved an area under the curve (AUC) of 0.75, 0.74, and 0.70, respectively, and a combination of IWL-CNFD achieved an AUC of 0.76. CONCLUSIONS: The parameters of CNFD, CNFL, and IWL have a comparable ability to diagnose patients with DPN. However, IWL detects an abnormality even in patients without DPN. Combining IWL with CNFD may improve the diagnostic performance of CCM.

Gravity-induced coronal plane joint moments in adolescent idiopathic scoliosis

Background Adolescent Idiopathic Scoliosis is the most common type of spinal deformity, and whilst the risk of progression appears to be biomechanically mediated (larger deformities are more likely to progress), the detailed biomechanical mechanisms driving progression are not well understood. Gravitational forces in the upright position are the primary sustained loads experienced by the spine. In scoliosis they are asymmetrical, generating moments about the spinal joints which may promote asymmetrical growth and deformity progression. Using 3D imaging modalities to estimate segmental torso masses allows the gravitational loading on the scoliotic spine to be determined. The resulting distribution of joint moments aids understanding of the mechanics of scoliosis progression. Methods Existing low-dose CT scans were used to estimate torso segment masses and joint moments for 20 female scoliosis patients. Intervertebral joint moments at each vertebral level were found by summing the moments of each of the torso segment masses above the required joint. Results The patients’ mean age was 15.3 years (SD 2.3; range 11.9 – 22.3 years); mean thoracic major Cobb angle 52° (SD 5.9°; range 42°-63°) and mean weight 57.5 kg (SD 11.5 kg; range 41 – 84.7 kg). Joint moments of up to 7 Nm were estimated at the apical level. No significant correlation was found between the patients’ major Cobb angles and apical joint moments. Conclusions Patients with larger Cobb angles do not necessarily have higher joint moments, and curve shape is an important determinant of joint moment distribution. These findings may help to explain the variations in progression between individual patients. This study suggests that substantial corrective forces are required of either internal instrumentation or orthoses to effectively counter the gravity-induced moments acting to deform the spinal joints of idiopathic scoliosis patients.

Sonography of paediatric appendicitis and its secondary signs

Sonographic diagnosis of appendicitis in children is an important clinical tool, often obviating the need for potentially harmful ionising radiation from computed tomography (CT) scans and unnecessary appendectomies. Established criteria do not commonly account for the sonographic secondary signs of acute appendicitis as an adjunct or corollary to an identifiably inflamed appendix. If one of, or combinations of these secondary signs are a reliable positive and/or negative indicator of the condition, diagnostic accuracy may be improved. This will be of particular importance in cases where the appendix cannot be easily identified, possibly providing referring clinicians with a less equivocal diagnosis. Acute appendicitis (AA) is the most common emergency presentation requiring surgical intervention among both adults and children. During 2010-11 in Australia 25000 appendicectomies were performed on adults and children, more than double the number of the next most common surgical procedure [1]. Ultrasound has been commonly used to diagnose AA since the 1980s, however the best imaging modality or combination of modalities to accurately and cost-effectively diagnose the condition is still debated. A study by Puylaert advocated ultrasound in all presentations [2], whereas others suggested it only as a first line modality [3–5]. Conversely, York et al state that it is not appropriate as it delays treatment [6]. CT has been shown to more accurately diagnose AA than ultrasound, however its inherent radiation risks warrant cautionary use in children [7]. Improved accuracy in the diagnosis of suspected AA using ultrasound would enable surgeons to make a decision without the need to expose children to the potentially harmful effects of CT. Secondary signs of appendicitis are well established [8], although research into their predictive values has only recently been undertaken [9,10] indicating their potential diagnostic benefit in the absence of an identifiable appendix. The purpose of this review is to examine the history of appendiceal sonography, established sonographic criteria, paediatric specific techniques and the predictive value of secondary signs.

Ratios of T-cell immune-effectors with tumour associated macrophages and PD-1/PD-L1 axis immune-checkpoint molecules, as prognosticators in diffuse large B cell lymphoma: A population-based study

Background Risk-stratification of diffuse large B-cell lymphoma (DLBCL) requires identification of patients with disease that is not cured despite initial R-CHOP. Although the prognostic importance of the tumour microenvironment (TME) is established, the optimal strategy to quantify it is unknown. Methods The relationship between immune-effector and inhibitory (checkpoint) genes was assessed by NanoString™ in 252 paraffin-embedded DLBCL tissues. A model to quantify net anti-tumoural immunity as an outcome predictor was tested in 158 R-CHOP treated patients, and validated in tissue/blood from two independent R-CHOP treated cohorts of 233 and 140 patients respectively. Findings T and NK-cell immune-effector molecule expression correlated with tumour associated macrophage and PD-1/PD-L1 axis markers consistent with malignant B-cells triggering a dynamic checkpoint response to adapt to and evade immune-surveillance. A tree-based survival model was performed to test if immune-effector to checkpoint ratios were prognostic. The CD4*CD8:(CD163/CD68)*PD-L1 ratio was better able to stratify overall survival than any single or combination of immune markers, distinguishing groups with disparate 4-year survivals (92% versus 47%). The immune ratio was independent of and added to the revised international prognostic index (R-IPI) and cell-of-origin (COO). Tissue findings were validated in 233 DLBCL R-CHOP treated patients. Furthermore, within the blood of 140 R-CHOP treated patients immune-effector:checkpoint ratios were associated with differential interim-PET/CT+ve/-ve expression.

Highly efficient spin-conversion effect leading to energy up-converted electroluminescence in singlet fission photovoltaics

Free charge generation in donor-acceptor (D-A) based organic photovoltaic diodes (OPV) progresses through formation of charge-transfer (CT) and charge-separated (CS) states and excitation decay to the triplet level is considered as a terminal loss. On the other hand a direct excitation decay to the triplet state is beneficial for multiexciton harvesting in singlet fission photovoltaics (SF-PV) and the formation of CT-state is considered as a limiting factor for multiple triplet harvesting. These two extremes when present in a D-A system are expected to provide important insights into the mechanism of free charge generation and spin-character of bimolecular recombination in OPVs. Herein, we present the complete cycle of events linked to spin conversion in the model OPV system of rubrene/C60. By tracking the spectral evolution of photocurrent generation at short-circuit and close to open-circuit conditions we are able to capture spectral changes to photocurrent that reveal the triplet character of CT-state. Furthermore, we unveil an energy up-conversion effect that sets in as a consequence of triplet population build-up where triplet-triplet annihilation (TTA) process effectively regenerates the singlet excitation. This detailed balance is shown to enable a rare event of photon emission just above the open-circuit voltage (VOC) in OPVs.