Iowa State Daily (February 8, 2012)

ContentsReaching new heightsCareer fair creates connectionsProfessor named to energy committeeVictory could stretch winning streak to threeShare love with Valentine's treatsIs Christianity really so different?Create homemade crafts to celebrate the holiday

Human intervertebral disc stiffness correlates better with the Otsu threshold computed from axial T2 map of its posterior annulus fibrosus than with clinical classifications

Degeneration of the intervertebral disc, sometimes associated with low back pain and abnormal spinal motions, represents a major health issue with high costs. A non-invasive degeneration assessment via qualitative or quantitative MRI (magnetic resonance imaging) is possible, yet, no relation between mechanical properties and T2 maps of the intervertebral disc (IVD) has been considered, albeit T2 relaxation time values quantify the degree of degeneration. Therefore, MRI scans and mechanical tests were performed on 14 human lumbar intervertebral segments freed from posterior elements and all soft tissues excluding the IVD. Degeneration was evaluated in each specimen using morphological criteria, qualitative T2 weighted images and quantitative axial T2 map data and stiffness was calculated from the load-deflection curves of in vitro compression, torsion, lateral bending and flexion/extension tests. In addition to mean T2, the OTSU threshold of T2 (TOTSU), a robust and automatic histogram-based method that computes the optimal threshold maximizing the distinction of two classes of values, was calculated for anterior, posterior, left and right regions of each annulus fibrosus (AF). While mean T2 and degeneration schemes were not related to the IVDs' mechanical properties, TOTSU computed in the posterior AF correlated significantly with those classifications as well as with all stiffness values. TOTSU should therefore be included in future degeneration grading schemes.

Robust Proximal Femur Segmentation in Conventional X-Ray Images via Random Forest Regression on Multi-resolution Gradient Features

In this paper, we propose a fully automatic, robust approach for segmenting proximal femur in conventional X-ray images. Our method is based on hierarchical landmark detection by random forest regression, where the detection results of 22 global landmarks are used to do the spatial normalization, and the detection results of the 59 local landmarks serve as the image cue for instantiation of a statistical shape model of the proximal femur. To detect landmarks in both levels, we use multi-resolution HoG (Histogram of Oriented Gradients) as features which can achieve better accuracy and robustness. The efficacy of the present method is demonstrated by experiments conducted on 150 clinical x-ray images. It was found that the present method could achieve an average point-to-curve error of 2.0 mm and that the present method was robust to low image contrast, noise and occlusions caused by implants.

Improved detection of amnestic MCI by means of discriminative vector quantization of single-trial cognitive ERP responses

Cognitive event-related potentials (ERPs) are widely employed in the study of dementive disorders. The morphology of averaged response is known to be under the influence of neurodegenerative processes and exploited for diagnostic purposes. This work is built over the idea that there is additional information in the dynamics of single-trial responses. We introduce a novel way to detect mild cognitive impairment (MCI) from the recordings of auditory ERP responses. Using single trial responses from a cohort of 25 amnestic MCI patients and a group of age-matched controls, we suggest a descriptor capable of encapsulating single-trial (ST) response dynamics for the benefit of early diagnosis. A customized vector quantization (VQ) scheme is first employed to summarize the overall set of ST-responses by means of a small-sized codebook of brain waves that is semantically organized. Each ST-response is then treated as a trajectory that can be encoded as a sequence of code vectors. A subject's set of responses is consequently represented as a histogram of activated code vectors. Discriminating MCI patients from healthy controls is based on the deduced response profiles and carried out by means of a standard machine learning procedure. The novel response representation was found to improve significantly MCI detection with respect to the standard alternative representation obtained via ensemble averaging (13% in terms of sensitivity and 6% in terms of specificity). Hence, the role of cognitive ERPs as biomarker for MCI can be enhanced by adopting the delicate description of our VQ scheme.

Watson–Crick and Sugar-Edge Base Pairing of Cytosine in the Gas Phase: UV and Infrared Spectra of Cytosine·2-Pyridone

While keto-amino cytosine is the dominant species in aqueous solution, spectroscopic studies in molecular beams and in noble gas matrices show that other cytosine tautomers prevail in apolar environments. Each of these offers two or three H-bonding sites (Watson–Crick, wobble, sugar-edge). The mass- and isomer-specific S1 ← S0 vibronic spectra of cytosine·2-pyridone (Cyt·2PY) and 1-methylcytosine·2PY are measured using UV laser resonant two-photon ionization (R2PI), UV/UV depletion, and IR depletion spectroscopy. The UV spectra of the Watson–Crick and sugar-edge isomers of Cyt·2PY are separated using UV/UV spectral hole-burning. Five different isomers of Cyt·2PY are observed in a supersonic beam. We show that the Watson–Crick and sugar-edge dimers of keto-amino cytosine with 2PY are the most abundant in the beam, although keto-amino-cytosine is only the third most abundant tautomer in the gas phase. We identify the different isomers by combining three different diagnostic tools: (1) methylation of the cytosine N1–H group prevents formation of both the sugar-edge and wobble isomers and gives the Watson–Crick isomer exclusively. (2) The calculated ground state binding and dissociation energies, relative gas-phase abundances, excitation and the ionization energies are in agreement with the assignment of the dominant Cyt·2PY isomers to the Watson–Crick and sugar-edge complexes of keto-amino cytosine. (3) The comparison of calculated ground state vibrational frequencies to the experimental IR spectra in the carbonyl stretch and NH/OH/CH stretch ranges strengthen this identification.

Characterization of microcirculation in multiple sclerosis lesions by dynamic texture parameter analysis (DTPA)

OBJECTIVE Texture analysis is an alternative method to quantitatively assess MR-images. In this study, we introduce dynamic texture parameter analysis (DTPA), a novel technique to investigate the temporal evolution of texture parameters using dynamic susceptibility contrast enhanced (DSCE) imaging. Here, we aim to introduce the method and its application on enhancing lesions (EL), non-enhancing lesions (NEL) and normal appearing white matter (NAWM) in multiple sclerosis (MS). METHODS We investigated 18 patients with MS and clinical isolated syndrome (CIS), according to the 2010 McDonald's criteria using DSCE imaging at different field strengths (1.5 and 3 Tesla). Tissues of interest (TOIs) were defined within 27 EL, 29 NEL and 37 NAWM areas after normalization and eight histogram-based texture parameter maps (TPMs) were computed. TPMs quantify the heterogeneity of the TOI. For every TOI, the average, variance, skewness, kurtosis and variance-of-the-variance statistical parameters were calculated. These TOI parameters were further analyzed using one-way ANOVA followed by multiple Wilcoxon sum rank testing corrected for multiple comparisons. RESULTS Tissue- and time-dependent differences were observed in the dynamics of computed texture parameters. Sixteen parameters discriminated between EL, NEL and NAWM (pAVG = 0.0005). Significant differences in the DTPA texture maps were found during inflow (52 parameters), outflow (40 parameters) and reperfusion (62 parameters). The strongest discriminators among the TPMs were observed in the variance-related parameters, while skewness and kurtosis TPMs were in general less sensitive to detect differences between the tissues. CONCLUSION DTPA of DSCE image time series revealed characteristic time responses for ELs, NELs and NAWM. This may be further used for a refined quantitative grading of MS lesions during their evolution from acute to chronic state. DTPA discriminates lesions beyond features of enhancement or T2-hypersignal, on a numeric scale allowing for a more subtle grading of MS-lesions.

Hydrostatic intestinal edema induced signaling pathways: potential role of mechanical forces.

BACKGROUND: Hydrostatic intestinal edema initiates a signal transduction cascade that results in smooth muscle contractile dysfunction. Given the rapid and concurrent alterations in the mechanical properties of edematous intestine observed with the development of edema, we hypothesize that mechanical forces may serve as a stimulus for the activation of certain signaling cascades. We sought to examine whether isolated similar magnitude mechanical forces induced the same signal transduction cascades associated with edema. METHODS: The distal intestine from adult male Sprague Dawley rats was stretched longitudinally for 2 h to 123% its original length, which correlates with the interstitial stress found with edema. We compared wet-to-dry ratios, myeloperoxidase activity, nuclear signal transduction and activator of transcription (STAT)-3 and nuclear factor (NF)-kappa B DNA binding, STAT-3 phosphorylation, myosin light chain phosphorylation, baseline and maximally stimulated intestinal contractile strength, and inducible nitric oxide synthase (iNOS) and sodium hydrogen exchanger 1-3 messenger RNA (mRNA) in stretched and adjacent control segments of intestine. RESULTS: Mechanical stretch did not induce intestinal edema or an increase in myeloperoxidase activity. Nuclear STAT-3 DNA binding, STAT-3 phosphorylation, and nuclear NF-kappa B DNA binding were significantly increased in stretched seromuscular samples. Increased expression of sodium hydrogen exchanger 1 was found but not an increase in iNOS expression. Myosin light chain phosphorylation was significantly decreased in stretched intestine as was baseline and maximally stimulated intestinal contractile strength. CONCLUSION: Intestinal stretch, in the absence of edema/inflammatory/ischemic changes, leads to the activation of signaling pathways known to be altered in intestinal edema. Edema may initiate a mechanotransductive cascade that is responsible for the subsequent activation of various signaling cascades known to induce contractile dysfunction.

ANKRD1, the gene encoding cardiac ankyrin repeat protein, is a novel dilated cardiomyopathy gene.

OBJECTIVES: We evaluated ankyrin repeat domain 1 (ANKRD1), the gene encoding cardiac ankyrin repeat protein (CARP), as a novel candidate gene for dilated cardiomyopathy (DCM) through mutation analysis of a cohort of familial or idiopathic DCM patients, based on the hypothesis that inherited dysfunction of mechanical stretch-based signaling is present in a subset of DCM patients. BACKGROUND: CARP, a transcription coinhibitor, is a member of the titin-N2A mechanosensory complex and translocates to the nucleus in response to stretch. It is up-regulated in cardiac failure and hypertrophy and represses expression of sarcomeric proteins. Its overexpression results in contractile dysfunction. METHODS: In all, 208 DCM patients were screened for mutations/variants in the coding region of ANKRD1 using polymerase chain reaction, denaturing high-performance liquid chromatography, and direct deoxyribonucleic acid sequencing. In vitro functional analyses of the mutation were performed using yeast 2-hybrid assays and investigating the effect on stretch-mediated gene expression in myoblastoid cell lines using quantitative real-time reverse transcription-polymerase chain reaction. RESULTS: Three missense heterozygous ANKRD1 mutations (P105S, V107L, and M184I) were identified in 4 DCM patients. The M184I mutation results in loss of CARP binding with Talin 1 and FHL2, and the P105S mutation in loss of Talin 1 binding. Intracellular localization of mutant CARP proteins is not altered. The mutations result in differential stretch-induced gene expression compared with wild-type CARP. CONCLUSIONS: ANKRD1 is a novel DCM gene, with mutations present in 1.9% of DCM patients. The ANKRD1 mutations may cause DCM as a result of disruption of the normal cardiac stretch-based signaling.

Assessment of Collimator Jaw Optimization in Reducing Normal Tissue Irradiation with Intensity Modulated Radiation Therapy

Purpose: To evaluate normal tissue dose reduction in step-and-shoot intensity-modulated radiation therapy (IMRT) on the Varian 2100 platform by tracking the multileaf collimator (MLC) apertures with the accelerator jaws. Methods: Clinical radiation treatment plans for 10 thoracic, 3 pediatric and 3 head and neck patients were converted to plans with the jaws tracking each segment’s MLC apertures. Each segment was then renormalized to account for the change in collimator scatter to obtain target coverage within 1% of that in the original plan. The new plans were compared to the original plans in a commercial radiation treatment planning system (TPS). Reduction in normal tissue dose was evaluated in the new plan by using the parameters V5, V10, and V20 in the cumulative dose-volume histogram for the following structures: total lung minus GTV (gross target volume), heart, esophagus, spinal cord, liver, parotids, and brainstem. In order to validate the accuracy of our beam model, MLC transmission measurements were made and compared to those predicted by the TPS. Results: The greatest change between the original plan and new plan occurred at lower dose levels. The reduction in V20 was never more than 6.3% and was typically less than 1% for all patients. The reduction in V5 was 16.7% maximum and was typically less than 3% for all patients. The variation in normal tissue dose reduction was not predictable, and we found no clear parameters that indicated which patients would benefit most from jaw tracking. Our TPS model of MLC transmission agreed with measurements with absolute transmission differences of less than 0.1 % and thus uncertainties in the model did not contribute significantly to the uncertainty in the dose determination. Conclusion: The amount of dose reduction achieved by collimating the jaws around each MLC aperture in step-and-shoot IMRT does not appear to be clinically significant.

Improving the accuracy of radiation pneumonitis dose response models

The prognosis for lung cancer patients remains poor. Five year survival rates have been reported to be 15%. Studies have shown that dose escalation to the tumor can lead to better local control and subsequently better overall survival. However, dose to lung tumor is limited by normal tissue toxicity. The most prevalent thoracic toxicity is radiation pneumonitis. In order to determine a safe dose that can be delivered to the healthy lung, researchers have turned to mathematical models predicting the rate of radiation pneumonitis. However, these models rely on simple metrics based on the dose-volume histogram and are not yet accurate enough to be used for dose escalation trials. The purpose of this work was to improve the fit of predictive risk models for radiation pneumonitis and to show the dosimetric benefit of using the models to guide patient treatment planning. The study was divided into 3 specific aims. The first two specifics aims were focused on improving the fit of the predictive model. In Specific Aim 1 we incorporated information about the spatial location of the lung dose distribution into a predictive model. In Specific Aim 2 we incorporated ventilation-based functional information into a predictive pneumonitis model. In the third specific aim a proof of principle virtual simulation was performed where a model-determined limit was used to scale the prescription dose. The data showed that for our patient cohort, the fit of the model to the data was not improved by incorporating spatial information. Although we were not able to achieve a significant improvement in model fit using pre-treatment ventilation, we show some promising results indicating that ventilation imaging can provide useful information about lung function in lung cancer patients. The virtual simulation trial demonstrated that using a personalized lung dose limit derived from a predictive model will result in a different prescription than what was achieved with the clinically used plan; thus demonstrating the utility of a normal tissue toxicity model in personalizing the prescription dose.

Molecular characterization of a gene required for entry into S phase of the mitotic cell cycle in the yeast {\it Saccharomyces cerevisiae\/}

A strain of Saccaromyces cerevisiae (SC3B) with a temperature sensitive defect in the synthesis of DNA has been isolated. This defect is due to a single recessive mutation in a gene named INS1 required for the initiation of S phase. Arrested cells carrying the ins1$\sp{ts}$ allele are defective in the completion of G1 to S phase transition events including SPB duplication or separation, initiation of DNA synthesis, normal control of budding, and bud neck stability. The mutation and a gene which complements the mutation were mapped to chromosome IV. The complementing gene was proved to be the wild type allele of the temperature sensitive mutation by genetic linkage of an integrated clone. A very low abundance 4.2 kb RNA message was observed in the strain SC3B which increased greatly in this strain transformed with a multiple copy plasmid carrying the complementing clone. The wild type gene was sequenced and found to encode a 1268 amino acid protein of with a molecular weight of 142,655 Daltons. Computer assisted searches for similar DNA sequences revealed no significant homology matches. However, searches for protein sequence homology revealed a protein (the DIS3 gene product of S. pombe) with a similar sequence over a 534 amino acid stretch to the predicted INS1 gene product. A later search revealed a near identical sequence for a gene (SRK1) also isolated from S. cerevisiae. ^

Regulation of Spec genes in aboral ectoderm cells of sea urchin {\it Strongylocentrotus purpuratus}

The Spec genes of the sea urchin Stronylocentrotus purpuratus serves as an excellent model for studying cell type-specific gene expression during early embryogenesis. The Spec1/Spec2 genes encode cytosolic calcium-binding proteins related to the calmodulin/troponin C/myosin light chain superfamily. Members of the Spec gene family are activated shortly after the sixth cleavage as the lineage-specific founder cells giving rise to aboral ectoderm are established, and the accumulation of the Spec mRNAs is limited exclusively to aboral ectoderm cell lineages. In this dissertation, the transcriptional regulation of the Spec genes was studied. Sequence comparisons of the Spec gene 5$\sp\prime$ flanking regions showed that a DNA block of approximately 800 bp from the 3$\sp\prime$ end of the first exon to the 5$\sp\prime$ end of a repetitive DNA element, termed RSR, was highly conserved. In Spec2a, the conserved region was a continuous stretch of DNA, but in Spec1 and Spec2c, DNA insertions interrupt the conserved sequence block and alter the relative placement of the RSR element and other 5$\sp\prime$ flanking DNA. Thus, drastic rearrangements have occurred within the putative control regions of the Spec genes. In vivo expression experiments using the sea urchin embryo gene-transfer system showed that while the 5$\sp\prime$ flanking regions of all three Spec genes conferred proper temporal activation to the reporter CAT gene, only the Spec2a 5$\sp\prime$ flanking region could restrict lacZ gene expression to aboral ectoderm cells. However, the Spec2a conserved region alone was not sufficient to confer proper spatial expression, suggesting that negative spatial elements are also associated with the proper activation of Spec2a. A major positive regulatory region, defined as the RSR enhancer, was identified between base pairs $-$631 and $-$443 on Spec2a. The RSR enhancer was essential for maximal activity and conferred preferential aboral ectoderm expression to a lacZ reporter gene. DNaseI footprinting and band-shift analysis of the RSR enhancer revealed multiple DNA-elements. One of the elements, an A/T-rich sequence called the A/T palindrome was studied in detail. This element binds a single 45-kDa nuclear protein, the A/T palindrome binding protein (A/TBP), whose DNA-binding specificity suggests a possible relationship with the bicoid-class homeodomain proteins. Mutated A/T palindromes are incapable of binding the 45-kDa protein and lower promoter activity by 8-fold. DNA-binding activity for A/TBP is low in unfertilized eggs, increases by the 16-cell stage and continues rising in blastulae. These data suggest that A/TBP plays a major role in the activation of the Spec2a gene in aboral ectoderm cells. ^

A dosimetric study of radionuclide therapy for bone marrow ablation

In a phase I clinical trial, six multiple myeloma patients, who were non-responsive to conventional therapy and were scheduled for bone marrow transplantation, received Holmium-166 ($\sp{166}$Ho) labeled to a bone seeking agent, DOTMP (1,4,7,10-tetraazacyclododecane-1,4,7,10-tetramethylene-phosphonic acid), for the purpose of bone marrow ablation. The specific aims of my research within this protocol were to evaluate the toxicity and efficacy of $\sp{166}$Ho DOTMP by quantifying the in vivo pharmacokinetics and radiation dosimetry, and by correlating these results to the biologic response observed. The reproducibility of pharmacokinetics from multiple injections of $\sp{166}$Ho DOTMP administered to these myeloma patients was demonstrated from both blood and whole body retention. The skeletal concentration of $\sp{166}$Ho DOTMP was heterogenous in all six patients: high in the ribs, pelvis, and lumbar vertebrae regions, and relatively low in the femurs, arms, and head.^ A novel technique was developed to calculate the radiation dose to the bone marrow in each skeletal ROI, and was applied to all six $\sp{166}$Ho DOTMP patients. Radiation dose estimates for the bone marrow calculated using the standard MIRD "S" factors were compared with the average values derived from the heterogenous distribution of activity in the skeleton (i.e., the regional technique). The results from the two techniques were significantly different; the average of the dose estimates from the regional technique were typically 30% greater. Furthermore, the regional technique provided a range of radiation doses for the entire marrow volume, while the MIRD "S" factors only provided a single value. Dose volume histogram analysis of data from the regional technique indicated a range of dose estimates that varied by a factor of 10 between the high dose and low dose regions. Finally, the observed clinical response of cells and abnormal proteins measured in bone marrow aspirates and peripheral blood samples were compared with radiation dose estimates for the bone marrow calculated from the standard and regional technique. The results showed the regional technique values correlated more closely to several clinical response parameters. (Abstract shortened by UMI.) ^

Phosphorylation in the regulation of muscle-specific transcription

The contents of this dissertation include studies on the mechanisms by which FGF and growth factor down-stream kinases inactivate myogenin; characterization of myogenin phosphorylation and its role in regulation of myogenin activity; analysis the C-terminal transcriptional activation domain of myogenin; studies on the nuclear localization of myogenin and characterization of proteins that interact with PKC.^ Activation of muscle transcription by the MyoD family requires their heterodimerization with ubiquitous bHLH proteins such as the E2A gene products E12 and E47. I have shown that dimerization with E2A products potentiates phosphorylation of myogenin at serine 43 in its amino-terminus and serine 170 in the carboxyl-terminal transcription activation domains. Mutations of these sites resulted in enhanced transcriptional activity of myogenin, suggesting that their phosphorylation diminishes myogenin's transcriptional activity. Consistent with the role of phosphorylation at serine 170, analysis of the carboxyl-terminal transcriptional activation domain by deletion has revealed a stretch of residues from 157 to 170 which functions as a negative element for myogenin activity.^ In addition to inducing phosphorylation of myogenin, E12 also localizes myogenin to the nucleus. The DNA binding and dimerization mutants of myogenin show various deficiencies in nuclear localization. Cotransfection of E12 with the DNA binding mutants, but not a dimerization mutant, greatly enhances their nuclear binding. These data suggest that the nuclear localization signal is located in the DNA binding region and myogenin can also be nuclear localized by virtue of dimerizing with a nuclear protein.^ FGF is one of the most potent inhibitors of myogenesis and activates many down-stream pathways to exert its functions. One of these pathway is the MAP kinase pathway. Studies have shown that Raf-1 and Erk-1 kinase inactivate transactivation by myogenin and E proteins independent of DNA binding. The other is the PKC pathway. In transfected cells, FGF induces phosphorylation of thr-87 that maps to the previously identified PKC sites in the DNA binding domain of myogenin. Myogenin mutant T-N87 could resist the inhibition directed to the bHLH domain by FGF, suggesting that FGF inactivates myogenin by inducing phosphorylation of this site. In C2 myotubes, where FGF receptors are lost, the phosphatase inhibitor, okadaic acid, and phorbal ester PdBu, can also induce the phosphorylation of thr-87. This result supports the previous observation and suggests that in myotubes, other mechanisms, such as innervation, may inactivate myogenin through PKC induced phosphorylation.^ Many functions of PKC have been well documented, yet, little is known about the activators or effectors of PKC or proteins that mediate PKC nuclear localizations. Identification of PKC binding proteins will help to understand the molecular mechanism of PKC function. Two proteins that interact with the C kinase (PICKS) have been characterized, PICK-1 and PICK-2. PICK1 interacts with two conserved regions in the catalytic domain of PKC. It is localized to the perinuclear region and is phosphorylated in response to PKC activation. PICK2 is a novel protein with homology to the heat shock protein family. It interacts extensively with the catalytic domain of PKC and is localized in the cytoplasm in a punctate pattern. PICK1 and PICK2 may play important roles in mediating the actions of PKC. ^

A case of EDTA-dependent pseudothrombocytopenia: simple recognition of an underdiagnosed and misleading phenomenon

1. BMC Clin Pathol. 2014 May 1;14:19. doi: 10.1186/1472-6890-14-19. eCollection 2014. A case of EDTA-dependent pseudothrombocytopenia: simple recognition of an underdiagnosed and misleading phenomenon. Nagler M, Keller P, Siegrist D, Alberio L. Author information: Department of Hematology and Central Hematology Laboratory, Inselspital University Hospital and University of Berne, CH-3010 Berne, Switzerland. BACKGROUND: EDTA-dependent pseudothrombocytopenia (EDTA-PTCP) is a common laboratory phenomenon with a prevalence ranging from 0.1-2% in hospitalized patients to 15-17% in outpatients evaluated for isolated thrombocytopenia. Despite its harmlessness, EDTA-PTCP frequently leads to time-consuming, costly and even invasive diagnostic investigations. EDTA-PTCP is often overlooked because blood smears are not evaluated visually in routine practice and histograms as well as warning flags of hematology analyzers are not interpreted correctly. Nonetheless, EDTA-PTCP may be diagnosed easily even by general practitioners without any experiences in blood film examinations. This is the first report illustrating the typical patterns of a platelet (PLT) and white blood cell (WBC) histograms of hematology analyzers. CASE PRESENTATION: A 37-year-old female patient of Caucasian origin was referred with suspected acute leukemia and the crew of the emergency unit arranged extensive investigations for work-up. However, examination of EDTA blood sample revealed atypical lymphocytes and an isolated thrombocytopenia together with typical patterns of WBC and PLT histograms: a serrated curve of the platelet histogram and a peculiar peak on the left side of the WBC histogram. EDTA-PTCP was confirmed by a normal platelet count when examining citrated blood. CONCLUSION: Awareness of typical PLT and WBC patterns may alert to the presence of EDTA-PTCP in routine laboratory practice helping to avoid unnecessary investigations and over-treatment. PMCID: PMC4012027 PMID: 24808761 [PubMed]