Computational Techniques for Haplotype Inference and for Local Alignment Significance

This thesis which consists of an introduction and four peer-reviewed original publications studies the problems of haplotype inference (haplotyping) and local alignment significance. The problems studied here belong to the broad area of bioinformatics and computational biology. The presented solutions are computationally fast and accurate, which makes them practical in high-throughput sequence data analysis. Haplotype inference is a computational problem where the goal is to estimate haplotypes from a sample of genotypes as accurately as possible. This problem is important as the direct measurement of haplotypes is difficult, whereas the genotypes are easier to quantify. Haplotypes are the key-players when studying for example the genetic causes of diseases. In this thesis, three methods are presented for the haplotype inference problem referred to as HaploParser, HIT, and BACH. HaploParser is based on a combinatorial mosaic model and hierarchical parsing that together mimic recombinations and point-mutations in a biologically plausible way. In this mosaic model, the current population is assumed to be evolved from a small founder population. Thus, the haplotypes of the current population are recombinations of the (implicit) founder haplotypes with some point--mutations. HIT (Haplotype Inference Technique) uses a hidden Markov model for haplotypes and efficient algorithms are presented to learn this model from genotype data. The model structure of HIT is analogous to the mosaic model of HaploParser with founder haplotypes. Therefore, it can be seen as a probabilistic model of recombinations and point-mutations. BACH (Bayesian Context-based Haplotyping) utilizes a context tree weighting algorithm to efficiently sum over all variable-length Markov chains to evaluate the posterior probability of a haplotype configuration. Algorithms are presented that find haplotype configurations with high posterior probability. BACH is the most accurate method presented in this thesis and has comparable performance to the best available software for haplotype inference. Local alignment significance is a computational problem where one is interested in whether the local similarities in two sequences are due to the fact that the sequences are related or just by chance. Similarity of sequences is measured by their best local alignment score and from that, a p-value is computed. This p-value is the probability of picking two sequences from the null model that have as good or better best local alignment score. Local alignment significance is used routinely for example in homology searches. In this thesis, a general framework is sketched that allows one to compute a tight upper bound for the p-value of a local pairwise alignment score. Unlike the previous methods, the presented framework is not affeced by so-called edge-effects and can handle gaps (deletions and insertions) without troublesome sampling and curve fitting.

Postoperative soft tissue defects in the ankle area : The etiology and methods of reconstruction

The risk is obvious for soft tissue complications after operative treatment of the Achilles tendon, calcaneal bone or after ankle arthroplasty. Such complications after malleolar fractures are, however, seldom seen. The reason behind these complications is that the soft tissue in this region is tight and does not allow much tension to the wound area after surgery. Furthermore the area of operation may be damaged by swelling after the injury, or can be affected by peripheral vascular disease. While complications in this area are unavoidable, they can be diminished. This study attempts to highlight the possible predisposing factors leading to complications in these operations and on the other hand, to determine the solutions to solve soft tissue problems in this region. The study consists of five papers. The first article is a reprint on the soft tissue reconstruction of 25 patients after their complicated Achilles tendon surgeries were analysed. The second study reviews a series of 126 patients after having undergone an operative treatment of calcaneal bone fractures and analyses the complications and possible reasons behind them. The third part analyses a series of corrections of 35 soft tissue complications after calcaneal fracture operations. The fourth part reviews a series of 7 patients who had undergone complicated ankle arthroplasties. The last article presents a series of post operative lateral defects of the ankle treated with a less frequently used distally based peroneus brevis muscle flap and analyses the results. What can be conducted from these studies is that in general, the results after the correction of even severe soft tissue complications in the ankle region are good. For the small defects around the Achilles tendon, the local flaps are useful, but the larger defects are best treated with a free flap. We found that a long delay from trauma to surgery and a long operating time were predisposing factors that lead to soft tissue complications after operatively treated calcaneal bone fractures. The more severe the injury, the greater the risk for wound complication. Surprisingly, the long-term results after infected calcaneal osteosyntheses were acceptable and the calcaneal bone seems to tolerate chronic infections very well if the soft tissue is reconstructed successfully. Behind the complicated ankle arthroplasties, unexpectedly high number of cases experiencing arteriosclerosis of the lower extremity was found. These complications lead to ankle fusion but can be solved with a free flap if the vascularity is intact or can be reconstructed. For this reason a vascular examination of the lower extremity arteries of the patients going to ankle arthroplasty is strongly recommended. Moreover postoperative lateral malleolar wound infections which typically create lateral ankle defects can successfully be treated with a peroneus brevis muscle flap covered with a free skin graft.

Clinical Studies on Epidural and Spinal Postoperative Analgesia : With Special Reference to Continuous Techniques of Administration and Adjuvant Drugs

Continuous epidural analgesia (CEA) and continuous spinal postoperative analgesia (CSPA) provided by a mixture of local anaesthetic and opioid are widely used for postoperative pain relief. E.g., with the introduction of so-called microcatheters, CSPA found its way particularly in orthopaedic surgery. These techniques, however, may be associated with dose-dependent side-effects as hypotension, weakness in the legs, and nausea and vomiting. At times, they may fail to offer sufficient analgesia, e.g., because of a misplaced catheter. The correct position of an epidural catheter might be confirmed by the supposedly easy and reliable epidural stimulation test (EST). The aims of this thesis were to determine a) whether the efficacy, tolerability, and reliability of CEA might be improved by adding the α2-adrenergic agonists adrenaline and clonidine to CEA, and by the repeated use of EST during CEA; and, b) the feasibility of CSPA given through a microcatheter after vascular surgery. Studies I IV were double-blinded, randomized, and controlled trials; Study V was of a diagnostic, prospective nature. Patients underwent arterial bypass surgery of the legs (I, n=50; IV, n=46), total knee arthroplasty (II, n=70; III, n=72), and abdominal surgery or thoracotomy (V, n=30). Postoperative lumbar CEA consisted of regular mixtures of ropivacaine and fentanyl either without or with adrenaline (2 µg/ml (I) and 4 µg/ml (II)) and clonidine (2 µg/ml (III)). CSPA (IV) was given through a microcatheter (28G) and contained either ropivacaine (max. 2 mg/h) or a mixture of ropivacaine (max. 1 mg/h) and morphine (max. 8 µg/h). Epidural catheter tip position (V) was evaluated both by EST at the moment of catheter placement and several times during CEA, and by epidurography as reference diagnostic test. CEA and CSPA were administered for 24 or 48 h. Study parameters included pain scores assessed with a visual analogue scale, requirements of rescue pain medication, vital signs, and side-effects. Adrenaline (I and II) had no beneficial influence as regards the efficacy or tolerability of CEA. The total amounts of epidurally-infused drugs were even increased in the adrenaline group in Study II (p=0.02, RM ANOVA). Clonidine (III) augmented pain relief with lowered amounts of epidurally infused drugs (p=0.01, RM ANOVA) and reduced need for rescue oxycodone given i.m. (p=0.027, MW-U; median difference 3 mg (95% CI 0 7 mg)). Clonidine did not contribute to sedation and its influence on haemodynamics was minimal. CSPA (IV) provided satisfactory pain relief with only limited blockade of the legs (no inter-group differences). EST (V) was often related to technical problems and difficulties of interpretation, e.g., it failed to identify the four patients whose catheters were outside the spinal canal already at the time of catheter placement. As adjuvants to lumbar CEA, clonidine only slightly improved pain relief, while adrenaline did not provide any benefit. The role of EST applied at the time of epidural catheter placement or repeatedly during CEA remains open. The microcatheter CSPA technique appeared effective and reliable, but needs to be compared to routine CEA after peripheral arterial bypass surgery.

Modern computed tomography techniques in stroke management

Technological development of fast multi-sectional, helical computed tomography (CT) scanners has allowed computed tomography perfusion (CTp) and angiography (CTA) in evaluating acute ischemic stroke. This study focuses on new multidetector computed tomography techniques, namely whole-brain and first-pass CT perfusion plus CTA of carotid arteries. Whole-brain CTp data is acquired during slow infusion of contrast material to achieve constant contrast concentration in the cerebral vasculature. From these data quantitative maps are constructed of perfused cerebral blood volume (pCBV). The probability curve of cerebral infarction as a function of normalized pCBV was determined in patients with acute ischemic stroke. Normalized pCBV, expressed as a percentage of contralateral normal brain pCBV, was determined in the infarction core and in regions just inside and outside the boundary between infarcted and noninfarcted brain. Corresponding probabilities of infarction were 0.99, 0.96, and 0.11, R² was 0.73, and differences in perfusion between core and inner and outer bands were highly significant. Thus a probability of infarction curve can help predict the likelihood of infarction as a function of percentage normalized pCBV. First-pass CT perfusion is based on continuous cine imaging over a selected brain area during a bolus injection of contrast. During its first passage, contrast material compartmentalizes in the intravascular space, resulting in transient tissue enhancement. Functional maps such as cerebral blood flow (CBF), and volume (CBV), and mean transit time (MTT) are then constructed. We compared the effects of three different iodine concentrations (300, 350, or 400 mg/mL) on peak enhancement of normal brain tissue and artery and vein, stratified by region-of-interest (ROI) location, in 102 patients within 3 hours of stroke onset. A monotonic increasing peak opacification was evident at all ROI locations, suggesting that CTp evaluation of patients with acute stroke is best performed with the highest available concentration of contrast agent. In another study we investigated whether lesion volumes on CBV, CBF, and MTT maps within 3 hours of stroke onset predict final infarct volume, and whether all these parameters are needed for triage to intravenous recombinant tissue plasminogen activator (IV-rtPA). The effect of IV-rtPA on the affected brain by measuring salvaged tissue volume in patients receiving IV-rtPA and in controls was investigated also. CBV lesion volume did not necessarily represent dead tissue. MTT lesion volume alone can serve to identify the upper size limit of the abnormally perfused brain, and those with IV-rtPA salvaged more brain than did controls. Carotid CTA was compared with carotid DSA in grading of stenosis in patients with stroke symptoms. In CTA, the grade of stenosis was determined by means of axial source and maximum intensity projection (MIP) images as well as a semiautomatic vessel analysis. CTA provides an adequate, less invasive alternative to conventional DSA, although tending to underestimate clinically relevant grades of stenosis.

Noninvasive techniques in assessing coronary artery disease

Conventional invasive coronary angiography is the clinical gold standard for detecting of coronary artery stenoses. Noninvasive multidetector computed tomography (MDCT) in combination with retrospective ECG gating has recently been shown to permit visualization of the coronary artery lumen and detection of coronary artery stenoses. Single photon emission tomography (SPECT) perfusion imaging has been considered the reference method for evaluation of nonviable myocardium, but magnetic resonance imaging (MRI) can accurately depict structure, function, effusion, and myocardial viability, with an overall capacity unmatched by any other single imaging modality. Magnetocardiography (MCG) provides noninvasively information about myocardial excitation propagation and repolarization without the use of electrodes. This evolving technique may be considered the magnetic equivalent to electrocardiography. The aim of the present series of studies was to evaluate changes in the myocardium assessed with SPECT and MRI caused by coronary artery disease, examine the capability of multidetector computed tomography coronary angiography (MDCT-CA) to detect significant stenoses in the coronary arteries, and MCG to assess remote myocardial infarctions. Our study showed that in severe, progressing coronary artery disease laser treatment does not improve global left ventricular function or myocardial perfusion, but it does preserve systolic wall thickening in fixed defects (scar). It also prevents changes from ischemic myocardial regions to scar. The MCG repolarization variables are informative in remote myocardial infarction, and may perform as well as the conventional QRS criteria in detection of healed myocardial infarction. These STT abnormalities are more pronounced in patients with Q-wave infarction than in patients with non-Q-wave infarctions. MDCT-CA had a sensitivity of 82%, a specificity of 94%, a positive predictive value of 79%, and a negative predictive value of 95% for stenoses over 50% in the main coronary arteries as compared with conventional coronary angiography in patients with known coronary artery disease. Left ventricular wall dysfunction, perfusion defects, and infarctions were detected in 50-78% of sectors assigned to calcifications or stenoses, but also in sectors supplied by normally perfused coronary arteries. Our study showed a low sensitivity (sensitivity 63%) in detecting obstructive coronary artery disease assessed by MDCT in patients with severe aortic stenosis. Massive calcifications complicated correct assessment of the lumen of coronary arteries.

Microvascular reconstruction in extremity soft tissue sarcoma surgery

Soft tissue sarcomas (STS) are rare tumors of soft tissue occurring most frequently in the extremities. Modern treatment of extremity STS is based on limb-sparing surgery combined with radiotherapy. To prevent local recurrence, a healthy tissue margin of 2.5 cm around the resected tumor is required. This results in large defects of soft tissue and bone, necessitating the use of reconstructive surgery to achieve wound closure. When local or pedicled soft tissue flaps are unavailable, reconstruction with free flaps is used. Free flaps are elevated at a distant site, and have their blood flow restored at the recipient site through microvascular anastomosis. When limb-sparing surgery is made impossible, amputation is the only option. Proximal amputation such as forequarter amputation (FQA) causes considerable morbidity, but is nevertheless warranted for carefully selected patients for cure or palliation. 116 patients treated in 1985 - 2006 were included in the study. Of these, 93 patients treated with limb-sparing surgery and microvascular reconstructive surgery after resection of extremity STS. 25 patients who underwent FQA were also included. Patients were identified and their medical records retrospectively reviewed. In all, 105 free flap procedures were performed for 103 patients. A total of 95 curatively treated STS patients were included in survival analysis. The latissimus dorsi, used in 56% of cases, was the most frequently used free flap. Free flap success rate was 96%. There were 9% microvascular anastomosis complications and 15% wound complications. For curatively treated STS patients, local recurrence-free survival at 5 years was 73.1%, metastasis-free survival 58.3%, and overall disease-specific survival 68.9%. Functional results were good, with 75% of patients regaining normal or near-normal function after lower extremity, and 55% after upper extremity STS resection. Among curatively treated forequarter amputees, 5-year disease-free survival was 44%. In the palliatively treated group median time until disease death was 14 months. Microvascular reconstruction after extremity soft tissue sarcoma resection is safe and reliable, and produces well-healing wounds allowing early oncological treatment. Oncological outcome after these procedures is comparable to that of other extremity sarcoma patients. Functional results are generally good. Forequarter amputation is a useful treatment option for soft tissue tumors of the shoulder girdle and proximal upper extremity. When free flap coverage of extended forequarter amputation is required, the preferable flap is a fillet flap from the amputated extremity. Acceptable oncological outcome is achieved for curatively treated FQA patients. In the palliatively treated patient considerable periods of increased quality of life can be achieved.

Exclusion of an Exotic Top Quark with -4/3 Electric Charge Using Soft Lepton Tagging

We present a measurement of the electric charge of the top quark using $\ppbar$ collisions corresponding to an integrated luminosity of 2.7~fb$^{-1}$ at the CDF II detector. We reconstruct $\ttbar$ events in the lepton+jets final state and use kinematic information to determine which $b$-jet is associated with the leptonically- or hadronically-decaying $t$-quark. Soft lepton taggers are used to determine the $b$-jet flavor. Along with the charge of the $W$ boson decay lepton, this information permits the reconstruction of the top quark's electric charge. Out of 45 reconstructed events with $2.4\pm0.8$ expected background events, 29 are reconstructed as $\ttbar$ with the standard model $+$2/3 charge, whereas 16 are reconstructed as $\ttbar$ with an exotic $-4/3$ charge. This is consistent with the standard model and excludes the exotic scenario at 95\% confidence level. This is the strongest exclusion of the exotic charge scenario and the first to use soft leptons for this purpose.

Multivariate Techniques for Identifying Diffractive Interactions at the LHC

Close to one half of the LHC events are expected to be due to elastic or inelastic diffractive scattering. Still, predictions based on extrapolations of experimental data at lower energies differ by large factors in estimating the relative rate of diffractive event categories at the LHC energies. By identifying diffractive events, detailed studies on proton structure can be carried out. The combined forward physics objects: rapidity gaps, forward multiplicity and transverse energy flows can be used to efficiently classify proton-proton collisions. Data samples recorded by the forward detectors, with a simple extension, will allow first estimates of the single diffractive (SD), double diffractive (DD), central diffractive (CD), and non-diffractive (ND) cross sections. The approach, which uses the measurement of inelastic activity in forward and central detector systems, is complementary to the detection and measurement of leading beam-like protons. In this investigation, three different multivariate analysis approaches are assessed in classifying forward physics processes at the LHC. It is shown that with gene expression programming, neural networks and support vector machines, diffraction can be efficiently identified within a large sample of simulated proton-proton scattering events. The event characteristics are visualized by using the self-organizing map algorithm.

Measurement of the tt̅ production cross section in pp̅ collisions at √s=1.96 TeV using soft electron b-tagging

We present a measurement of the top quark pair production cross section in ppbar collisions at sqrt(s)=1.96 TeV using a data sample corresponding to 1.7/fb of integrated luminosity collected with the Collider Detector at Fermilab. We reconstruct ttbar events in the lepton+jets channel. The dominant background is the production of W bosons in association with multiple jets. To suppress this background, we identify electrons from the semileptonic decay of heavy-flavor jets. We measure a production cross section of 7.8 +/- 2.4 (stat) +/- 1.6 (syst) +/- 0.5 (lumi) pb. This is the first measurement of the top pair production cross section with soft electron tags in Run II of the Tevatron.

Measurement of the tt̅ production cross section in 2 fb-1 of pp̅ collisions at √s=1.96 TeV using lepton plus jets events with soft muon b tagging

We present a measurement of the tt̅ production cross section in pp̅ collisions at √s=1.96  TeV using events containing a high transverse momentum electron or muon, three or more jets, and missing transverse energy. Events consistent with tt̅ decay are found by identifying jets containing candidate heavy-flavor semileptonic decays to muons. The measurement uses a CDF run II data sample corresponding to 2  fb-1 of integrated luminosity. Based on 248 candidate events with three or more jets and an expected background of 79.5±5.3 events, we measure a production cross section of 9.1±1.6  pb.