Statistical methods for biomedical signal analysis and processing

Statistical modelling and statistical learning theory are two powerful analytical frameworks for analyzing signals and developing efficient processing and classification algorithms. In this thesis, these frameworks are applied for modelling and processing biomedical signals in two different contexts: ultrasound medical imaging systems and primate neural activity analysis and modelling. In the context of ultrasound medical imaging, two main applications are explored: deconvolution of signals measured from a ultrasonic transducer and automatic image segmentation and classification of prostate ultrasound scans. In the former application a stochastic model of the radio frequency signal measured from a ultrasonic transducer is derived. This model is then employed for developing in a statistical framework a regularized deconvolution procedure, for enhancing signal resolution. In the latter application, different statistical models are used to characterize images of prostate tissues, extracting different features. These features are then uses to segment the images in region of interests by means of an automatic procedure based on a statistical model of the extracted features. Finally, machine learning techniques are used for automatic classification of the different region of interests. In the context of neural activity signals, an example of bio-inspired dynamical network was developed to help in studies of motor-related processes in the brain of primate monkeys. The presented model aims to mimic the abstract functionality of a cell population in 7a parietal region of primate monkeys, during the execution of learned behavioural tasks.

Studies of the atmospheric muon flux with the ANTARES detector

The thesis main topic is the determination of the vertical component of the atmospheric muon flux as a function of the sea depth at the ANTARES site. ANTARES is a Cherenkov neutrino telescope placed at 2500m depth in the Mediterranean Sea at 40 km from the southern cost of France. In order to retrieve back the physical flux from the experimental data a deconvolution algorithm has been perform which takes into consideration the trigger inefficiensies and the reconstruction errors on the zenith angle. The obtained results are in good agreement with other ANTARES indipendent analysis.

Analysis and modeling techniques for ultrasonic tissue characterization

This thesis introduces new processing techniques for computer-aided interpretation of ultrasound images with the purpose of supporting medical diagnostic. In terms of practical application, the goal of this work is the improvement of current prostate biopsy protocols by providing physicians with a visual map overlaid over ultrasound images marking regions potentially affected by disease. As far as analysis techniques are concerned, the main contributions of this work to the state-of-the-art is the introduction of deconvolution as a pre-processing step in the standard ultrasonic tissue characterization procedure to improve the diagnostic significance of ultrasonic features. This thesis also includes some innovations in ultrasound modeling, in particular the employment of a continuous-time autoregressive moving-average (CARMA) model for ultrasound signals, a new maximum-likelihood CARMA estimator based on exponential splines and the definition of CARMA parameters as new ultrasonic features able to capture scatterers concentration. Finally, concerning the clinical usefulness of the developed techniques, the main contribution of this research is showing, through a study based on medical ground truth, that a reduction in the number of sampled cores in standard prostate biopsy is possible, preserving the same diagnostic power of the current clinical protocol.

Statistical methods for analysis and processing of medical ultrasound: applications to segmentation and restoration

In this thesis two major topics inherent with medical ultrasound images are addressed: deconvolution and segmentation. In the first case a deconvolution algorithm is described allowing statistically consistent maximum a posteriori estimates of the tissue reflectivity to be restored. These estimates are proven to provide a reliable source of information for achieving an accurate characterization of biological tissues through the ultrasound echo. The second topic involves the definition of a semi automatic algorithm for myocardium segmentation in 2D echocardiographic images. The results show that the proposed method can reduce inter- and intra observer variability in myocardial contours delineation and is feasible and accurate even on clinical data.

Untersuchung zur Bestimmung quantitativer Parameter der Lungenventilation mittels CT und MRT

Die zuverlässige Berechnung von quantitativen Parametern der Lungenventilation ist für ein Verständnis des Verhaltens der Lunge und insbesondere für die Diagnostik von Lungenerkrankungen von großer Bedeutung. Nur durch quantitative Parameter sind verlässliche und reproduzierbare diagnostische Aussagen über den Gesundheitszustand der Lunge möglich. Im Rahmen dieser Arbeit wurden neue quantitative Verfahren zur Erfassung der Lungenventilation basierend auf der dynamischen Computer- (CT) und Magnetresonanztomographie (MRT) entwickelt. Im ersten Teil dieser Arbeit wurde die Frage untersucht, ob das Aufblähen der Lunge in gesunden Schweinelungen und Lungen mit Akutem Lungenversagen (ARDS) durch einzelne, diskrete Zeitkonstanten beschrieben werden kann, oder ob kontinuierliche Verteilungen von Zeitkonstanten die Realität besser beschreiben. Hierzu wurden Serien dynamischer CT-Aufnahmen während definierter Beatmungsmanöver (Drucksprünge) aufgenommen und anschließend aus den Messdaten mittels inverser Laplace-Transformation die zugehörigen Verteilungen der Zeitkonstanten berechnet. Um die Qualität der Ergebnisse zu analysieren, wurde der Algorithmus im Rahmen von Simulationsrechnungen systematisch untersucht und anschließend in-vivo an gesunden und ARDS-Schweinelungen eingesetzt. Während in den gesunden Lungen mono- und biexponentielle Verteilungen bestimmt wurden, waren in den ARDS-Lungen Verteilungen um zwei dominante Zeitkonstanten notwendig, um die gemessenen Daten auf der Basis des verwendeten Modells verlässlich zu beschreiben. Es wurden sowohl diskrete als auch kontinuierliche Verteilungen gefunden. Die CT liefert Informationen über das solide Lungengewebe, während die MRT von hyperpolarisiertem 3He in der Lage ist, direkt das eingeatmete Gas abzubilden. Im zweiten Teil der Arbeit wurde zeitlich hochaufgelöst das Einströmen eines 3He-Bolus in die Lunge erfasst. Über eine Entfaltungsanalyse wurde anschließend das Einströmverhalten unter Idealbedingungen (unendlich kurzer 3He-Bolus), also die Gewebeantwortfunktion, berechnet und so eine Messtechnik-unabhängige Erfassung des Einströmens von 3He in die Lunge ermöglicht. Zentrale Fragestellung war hier, wie schnell das Gas in die Lunge einströmt. Im Rahmen von Simulationsrechnungen wurde das Verhalten eines Entfaltungsalgorithmus (basierend auf B-Spline Repräsentationen) systematisch analysiert. Zusätzlich wurde ein iteratives Entfaltungsverfahren eingesetzt. Aus zeitlich hochaufgelösten Messungen (7ms) an einer gesunden und einer ARDS-Schweinelunge konnte erstmals nachgewiesen werden, dass das Einströmen in-vivo in weniger als 0,1s geschieht. Die Ergebnisse zeigen Zeitkonstanten im Bereich von 4ms–50ms, wobei zwischen der gesunden Lungen und der ARDS-Lunge deutliche Unterschiede beobachtet wurden. Zusammenfassend ermöglichen daher die in dieser Arbeit vorgestellten Algorithmen eine objektivere Bestimmung quantitativer Parameter der Lungenventilation. Dies ist für die eindeutige Beschreibung ventilatorischer Vorgänge in der Lunge und somit für die Lungendiagnostik unerlässlich. Damit stehen quantitative Methoden für die Lungenfunktionsdiagnostik zur Verfügung, deren diagnostische Relevanz im Rahmen wissenschaftlicher und klinischer Studien untersucht werden kann.

Refined Estimation of Earthquake Source Parameters: Methods, Applications and Scaling Relationships

The objective of this work of thesis is the refined estimations of source parameters. To such a purpose we used two different approaches, one in the frequency domain and the other in the time domain. In frequency domain, we analyzed the P- and S-wave displacement spectra to estimate spectral parameters, that is corner frequencies and low frequency spectral amplitudes. We used a parametric modeling approach which is combined with a multi-step, non-linear inversion strategy and includes the correction for attenuation and site effects. The iterative multi-step procedure was applied to about 700 microearthquakes in the moment range 1011-1014 N•m and recorded at the dense, wide-dynamic range, seismic networks operating in Southern Apennines (Italy). The analysis of the source parameters is often complicated when we are not able to model the propagation accurately. In this case the empirical Green function approach is a very useful tool to study the seismic source properties. In fact the Empirical Green Functions (EGFs) consent to represent the contribution of propagation and site effects to signal without using approximate velocity models. An EGF is a recorded three-component set of time-histories of a small earthquake whose source mechanism and propagation path are similar to those of the master event. Thus, in time domain, the deconvolution method of Vallée (2004) was applied to calculate the source time functions (RSTFs) and to accurately estimate source size and rupture velocity. This technique was applied to 1) large event, that is Mw=6.3 2009 L’Aquila mainshock (Central Italy), 2) moderate events, that is cluster of earthquakes of 2009 L’Aquila sequence with moment magnitude ranging between 3 and 5.6, 3) small event, i.e. Mw=2.9 Laviano mainshock (Southern Italy).

Inferences on earthquake kinematic properties from data inversion: two different approaches

During my PhD, starting from the original formulations proposed by Bertrand et al., 2000 and Emolo & Zollo 2005, I developed inversion methods and applied then at different earthquakes. In particular large efforts have been devoted to the study of the model resolution and to the estimation of the model parameter errors. To study the source kinematic characteristics of the Christchurch earthquake we performed a joint inversion of strong-motion, GPS and InSAR data using a non-linear inversion method. Considering the complexity highlighted by superficial deformation data, we adopted a fault model consisting of two partially overlapping segments, with dimensions 15x11 and 7x7 km2, having different faulting styles. This two-fault model allows to better reconstruct the complex shape of the superficial deformation data. The total seismic moment resulting from the joint inversion is 3.0x1025 dyne.cm (Mw = 6.2) with an average rupture velocity of 2.0 km/s. Errors associated with the kinematic model have been estimated of around 20-30 %. The 2009 Aquila sequence was characterized by an intense aftershocks sequence that lasted several months. In this study we applied an inversion method that assumes as data the apparent Source Time Functions (aSTFs), to a Mw 4.0 aftershock of the Aquila sequence. The estimation of aSTFs was obtained using the deconvolution method proposed by Vallée et al., 2004. The inversion results show a heterogeneous slip distribution, characterized by two main slip patches located NW of the hypocenter, and a variable rupture velocity distribution (mean value of 2.5 km/s), showing a rupture front acceleration in between the two high slip zones. Errors of about 20% characterize the final estimated parameters.

Etablierung transgener BY-2-Zelllinien zur In-vivo-Lokalisierung pflanzlicher Cytoskelettproteine

Bei den Pflanzen sind viele Fragen bezüglich der Organisation und Regulation des bei der Zellteilung und differenzierung wichtigen Auf-, Ab- und Umbaus des Mikrotubuli-Netzwerkes noch immer offen, insbesondere was die Rolle des γ-Tubulins betrifft. Ziel der vorliegenden Arbeit war die Etablierung von BY-2 Modell-Zelllinien (Nicotiana), die verschiedene mit fluoreszierenden Proteinen (FP) markierte Elemente des Cytoskeletts exprimieren, um eine fluoreszenzmikroskopische Detektion in vivo zu ermöglichen.rnAls Grundlage für alle weiteren Versuche wurde eine zuverlässige Methode zur A. tumefaciens vermittelten stabilen Transfektion von BY-2 Zellen erarbeitet. Für die Expression von FP-markierten Cytoskelettproteinen, wurden entsprechende Fusionskonstrukte kloniert und via A. tumefaciens in BY-2 Zellen transferiert. So gelang zunächst die Herstellung transgener Zelllinien, die GFP-markiertes α- bzw. γ-Tubulin exprimierten. Diese sollten später als Basis für die Untersuchung des dynamischen Mikrotubuli-Netzwerkes bzw. dessen Regulation dienen. In beiden Zelllinien standen die Konstrukte zunächst unter Kontrolle eines doppelten 35S-Promotors, was zu einer starken, konstitutiven Expression der Transgene führte. Fluoreszenzmikroskopisch konnten Strukturen, an deren Aufbau Mikrotubuli beteiligt sind, detektiert werden. Aufgrund einer starken Hintergrundfluoreszenz, vermutlich bedingt durch die konstitutive Überexpression, war die Darstellung feinerer Bereiche, wie sie im Cytoskelett häufig auftreten, jedoch äußerst schwierig. Deshalb wurde eine schwächere bzw. adäquate Expressionsrate angestrebt. rnPhysiologische Expressionsraten sollten vor allem durch den endogenen γ-Tubulin-Promotor ermöglicht werden. Da die entsprechende Sequenz noch unbekannt war, wurde sie zunächst bestimmt und in ein passendes Konstrukt integriert. Fluoreszenzmikroskopische Untersuchungen der resultierenden Zelllinie ließen auf eine stark reduzierte Expressionsrate schließen. Tatsächlich war die Detektion von Cytoskelettstrukturen, wenn überhaupt, erst bei deutlich längeren Belichtungszeiten möglich. Bedingt durch die langen Belichtungszeiten wurde die Dokumentation durch eine latente pflanzentypische Autofluoreszenz der Zellen erschwert. Auch wenn hier keine detailreicheren Aufnahmen der Cytoskelettstrukturen möglich waren, ist die Zellkultur für weiterführende Untersuchungen, z.B. in Studien bezüglich des zeitlichen Expressionsmusters des γ-Tubulins, potentiell geeignet. Der Einsatz eines sensibleren Mikroskopsystems ist allerdings erforderlich. rnUm klären zu können, inwieweit γ-Tubulin mit den Mikrotubuli co-lokalisiert, wurden Zelllinien benötigt, bei denen die entsprechenden Elemente unterschiedlich markiert waren. Zu diesem Zweck wurde der Einsatz von RFP-markiertem Tubulin getestet. Eine deutliche Überexpression von RFP alleine war möglich. Trotz mehrfacher Wiederholung der Versuche war aber keine Expression von RFP-markiertem α-Tubulin in BY-2 Zellen zur Visualisierung der Mikrotubuli detektierbar. Die DNA-Sequenzen waren im Genom nachweisbar, eine Transkription jedoch nicht. Möglicherweise spielten hier gene silencing Effekte eine Rolle. Das verwendete RFP (TagRFP) und GFP stammten aus unterschiedlichen Organismen, aus einer Seeanemone bzw. einer Qualle. Eine Lösung könnte der Austausch des TagRFP durch ein Quallen-Derivat, das in einer von grün unterscheidbaren Farbe fluoresziert, bringen. Da bereits BY-2 Zelllinien vorliegen, die GFP-markiertes α- bzw. γ-Tubulin exprimieren, sollte es, nach Klonieren eines entsprechenden Konstruktes, zeitnah möglich sein, eine doppelt transfizierte Zelllinie herzustellen.

Dynamic characterisation of four nine-story large-panel R.C. buildings in Bishkek (Kyrgyzstan): A comparison between experimental ambient vibration analysis and numerical finite element modeling

With the outlook of improving seismic vulnerability assessment for the city of Bishkek (Kyrgyzstan), the global dynamic behaviour of four nine-storey r.c. large-panel buildings in elastic regime is studied. The four buildings were built during the Soviet era within a serial production system. Since they all belong to the same series, they have very similar geometries both in plan and in height. Firstly, ambient vibration measurements are performed in the four buildings. The data analysis composed of discrete Fourier transform, modal analysis (frequency domain decomposition) and deconvolution interferometry, yields the modal characteristics and an estimate of the linear impulse response function for the structures of the four buildings. Then, finite element models are set up for all four buildings and the results of the numerical modal analysis are compared with the experimental ones. The numerical models are finally calibrated considering the first three global modes and their results match the experimental ones with an error of less then 20%.

Calibration of a portal imaging device for high-precision dosimetry: a Monte Carlo study

Today electronic portal imaging devices (EPID's) are used primarily to verify patient positioning. They have, however, also the potential as 2D-dosimeters and could be used as such for transit dosimetry or dose reconstruction. It has been proven that such devices, especially liquid filled ionization chambers, have a stable dose response relationship which can be described in terms of the physical properties of the EPID and the pulsed linac radiation. For absolute dosimetry however, an accurate method of calibration to an absolute dose is needed. In this work, we concentrate on calibration against dose in a homogeneous water phantom. Using a Monte Carlo model of the detector we calculated dose spread kernels in units of absolute dose per incident energy fluence and compared them to calculated dose spread kernels in water at different depths. The energy of the incident pencil beams varied between 0.5 and 18 MeV. At the depth of dose maximum in water for a 6 MV beam (1.5 cm) and for a 18 MV beam (3.0 cm) we observed large absolute differences between water and detector dose above an incident energy of 4 MeV but only small relative differences in the most frequent energy range of the beam energy spectra. It is shown that for a 6 MV beam the absolute reference dose measured at 1.5 cm water depth differs from the absolute detector dose by 3.8%. At depth 1.2 cm in water, however, the relative dose differences are almost constant between 2 and 6 MeV. The effects of changes in the energy spectrum of the beam on the dose responses in water and in the detector are also investigated. We show that differences larger than 2% can occur for different beam qualities of the incident photon beam behind water slabs of different thicknesses. It is therefore concluded that for high-precision dosimetry such effects have to be taken into account. Nevertheless, the precise information about the dose response of the detector provided in this Monte Carlo study forms the basis of extracting directly the basic radiometric quantities photon fluence and photon energy fluence from the detector's signal using a deconvolution algorithm. The results are therefore promising for future application in absolute transit dosimetry and absolute dose reconstruction.

Oral heroin in opioid-dependent patients: pharmacokinetic comparison of immediate and extended release tablets

In diacetylmorphine prescription programs for heavily dependent addicts, diacetylmorphine is usually administered intravenously, but this may not be possible due to venosclerosis or when heroin abuse had occurred via non-intravenous routes. Since up to 25% of patients administer diacetylmorphine orally, we characterised morphine absorption after single oral doses of immediate and extended release diacetylmorphine in 8 opioid addicts. Plasma concentrations were determined by liquid chromatography-mass spectrometry. Non-compartmental methods and deconvolution were applied for data analysis. Mean (+/-S.D.) immediate and extended release doses were 719+/-297 and 956+/-404 mg, with high absolute morphine bioavailabilities of 56-61%, respectively. Immediate release diacetylmorphine caused rapid morphine absorption, peaking at 10-15 min. Morphine absorption was considerably slower and more sustained for extended release diacetylmorphine, with only approximately 30% of maximal immediate release absorption being reached after 10 min and maintained for 3-4h, with no relevant food interaction. The relative extended to immediate release bioavailability was calculated to be 86% by non-compartmental analysis and 93% by deconvolution analysis. Thus, immediate and extended release diacetylmorphine produce the intended morphine exposures. Both are suitable for substitution treatments. Similar doses can be applied if used in combination or sequentially.

Application of MFBD Algorithms to Image Reconstruction Under Anisoplanatic Conditions

All optical systems that operate in or through the atmosphere suffer from turbulence induced image blur. Both military and civilian surveillance, gun-sighting, and target identification systems are interested in terrestrial imaging over very long horizontal paths, but atmospheric turbulence can blur the resulting images beyond usefulness. My dissertation explores the performance of a multi-frame-blind-deconvolution technique applied under anisoplanatic conditions for both Gaussian and Poisson noise model assumptions. The technique is evaluated for use in reconstructing images of scenes corrupted by turbulence in long horizontal-path imaging scenarios and compared to other speckle imaging techniques. Performance is evaluated via the reconstruction of a common object from three sets of simulated turbulence degraded imagery representing low, moderate and severe turbulence conditions. Each set consisted of 1000 simulated, turbulence degraded images. The MSE performance of the estimator is evaluated as a function of the number of images, and the number of Zernike polynomial terms used to characterize the point spread function. I will compare the mean-square-error (MSE) performance of speckle imaging methods and a maximum-likelihood, multi-frame blind deconvolution (MFBD) method applied to long-path horizontal imaging scenarios. Both methods are used to reconstruct a scene from simulated imagery featuring anisoplanatic turbulence induced aberrations. This comparison is performed over three sets of 1000 simulated images each for low, moderate and severe turbulence-induced image degradation. The comparison shows that speckle-imaging techniques reduce the MSE 46 percent, 42 percent and 47 percent on average for low, moderate, and severe cases, respectively using 15 input frames under daytime conditions and moderate frame rates. Similarly, the MFBD method provides, 40 percent, 29 percent, and 36 percent improvements in MSE on average under the same conditions. The comparison is repeated under low light conditions (less than 100 photons per pixel) where improvements of 39 percent, 29 percent and 27 percent are available using speckle imaging methods and 25 input frames and 38 percent, 34 percent and 33 percent respectively for the MFBD method and 150 input frames. The MFBD estimator is applied to three sets of field data and the results presented. Finally, a combined Bispectrum-MFBD Hybrid estimator is proposed and investigated. This technique consistently provides a lower MSE and smaller variance in the estimate under all three simulated turbulence conditions.

Fenofibrate metabolism in the cynomolgus monkey using ultraperformance liquid chromatography-quadrupole time-of-flight mass spectrometry-based metabolomics

Fenofibrate, widely used for the treatment of dyslipidemia, activates the nuclear receptor, peroxisome proliferator-activated receptor alpha. However, liver toxicity, including liver cancer, occurs in rodents treated with fibrate drugs. Marked species differences occur in response to fibrate drugs, especially between rodents and humans, the latter of which are resistant to fibrate-induced cancer. Fenofibrate metabolism, which also shows species differences, has not been fully determined in humans and surrogate primates. In the present study, the metabolism of fenofibrate was investigated in cynomolgus monkeys by ultraperformance liquid chromatography-quadrupole time-of-flight mass spectrometry (UPLC-QTOFMS)-based metabolomics. Urine samples were collected before and after oral doses of fenofibrate. The samples were analyzed in both positive-ion and negative-ion modes by UPLC-QTOFMS, and after data deconvolution, the resulting data matrices were subjected to multivariate data analysis. Pattern recognition was performed on the retention time, mass/charge ratio, and other metabolite-related variables. Synthesized or purchased authentic compounds were used for metabolite identification and structure elucidation by liquid chromatographytandem mass spectrometry. Several metabolites were identified, including fenofibric acid, reduced fenofibric acid, fenofibric acid ester glucuronide, reduced fenofibric acid ester glucuronide, and compound X. Another two metabolites (compound B and compound AR), not previously reported in other species, were characterized in cynomolgus monkeys. More importantly, previously unknown metabolites, fenofibric acid taurine conjugate and reduced fenofibric acid taurine conjugate were identified, revealing a previously unrecognized conjugation pathway for fenofibrate.

CSF from Parkinson disease patients differentially affects cultured microglia and astrocytes.

BACKGROUND: Excessive and abnormal accumulation of alpha-synuclein (α-synuclein) is a factor contributing to pathogenic cell death in Parkinson's disease. The purpose of this study, based on earlier observations of Parkinson's disease cerebrospinal fluid (PD-CSF) initiated cell death, was to determine the effects of CSF from PD patients on the functionally different microglia and astrocyte glial cell lines. Microglia cells from human glioblastoma and astrocytes from fetal brain tissue were cultured, grown to confluence, treated with fixed concentrations of PD-CSF, non-PD disease control CSF, or control no-CSF medium, then photographed and fluorescently probed for α-synuclein content by deconvolution fluorescence microscopy. Outcome measures included manually counted cell growth patterns from day 1-8; α-synuclein density and distribution by antibody tagged 3D model stacked deconvoluted fluorescent imaging. RESULTS: After PD-CSF treatment, microglia growth was reduced extensively, and a non-confluent pattern with morphological changes developed, that was not evident in disease control CSF and no-CSF treated cultures. Astrocyte growth rates were similarly reduced by exposure to PD-CSF, but morphological changes were not consistently noted. PD-CSF treated microglia showed a significant increase in α-synuclein content by day 4 compared to other treatments (p ≤ 0.02). In microglia only, α-synuclein aggregated and redistributed to peri-nuclear locations. CONCLUSIONS: Cultured microglia and astrocytes are differentially affected by PD-CSF exposure compared to non-PD-CSF controls. PD-CSF dramatically impacts microglia cell growth, morphology, and α-synuclein deposition compared to astrocytes, supporting the hypothesis of cell specific susceptibility to PD-CSF toxicity.

Upregulation of defensins in burn sheep small intestine.

OBJECTIVE: The aim of this study was to visualize and localize the sheep antimicrobials, beta-defensins 1, 2, and 3, (SBD-1, SBD-2, SBD-3), sheep neutrophil defensin alpha (SNP-1), and the cathelicidin LL-37 in sheep small intestine after burn injury, our hypothesis being that these compounds would be upregulated in an effort to overcome a compromised endothelial lining. Response to burn injury includes the release of proinflammatory cytokines and systemic immune suppression that, if untreated, can progress to multiple organ failure and death, so protective mechanisms have to be initiated and implemented. METHODS: Tissue sections were probed with antibodies to the antimicrobials and then visualized with fluorescently labeled secondary antibodies and subjected to fluorescence deconvolution microscopy and image reconstruction. RESULTS: In both the sham and burn samples, all the aforementioned antimicrobials were seen in each of the layers of small intestine, the highest concentration being localized to the epithelium. SBD-2, SBD-3, and SNP-1 were upregulated in both enterocytes and Paneth cells, while SNP-1 and LL-37 showed increases in both the inner circular and outer longitudinal muscle layers of the muscularis externa following burn injury. Each of the defensins, except SBD-1, was also seen in between the muscle layers of the externa and while burn caused slight increases of SBD-2, SBD-3, and SNP-1 in this location, LL-37 content was significantly decreased. CONCLUSION: That while each of these human antimicrobials is present in multiple layers of sheep small intestine, SBD-2, SBD-3, SNP-1, and LL-37 are upregulated in the specific layers of the small intestine.