Dieses Manuskript beschreibt ein Protokoll für die Neutronenradiographie und Computertomographie biologischer Proben unter Verwendung eines Hochfluss-Isotopenreaktors (HFIR) CG-1D-Beamlines zur Messung eines Metallimplantats in einem Rattenoberschenkelknochen, einer Mauslunge und einem krautigen Pflanzenwurzel-/Bodensystem.
Neutronen wurden in der Vergangenheit für eine breite Palette biologischer Anwendungen verwendet, bei denen Techniken wie Neutronenkleinwinkelstreuung, Neutronenspinecho, Beugung und inelastische Streuung zum Einsatz kamen. Im Gegensatz zu Neutronenstreutechniken, die Informationen im reziproken Raum erhalten, misst die auf Abschwächung basierende Neutronenbildgebung ein Signal im realen Raum, das in der Größenordnung von zehn Mikrometern aufgelöst ist. Das Prinzip der Neutronenbildgebung folgt dem Beer-Lambert-Gesetz und basiert auf der Messung der Neutronenabschwächung durch eine Probe. Eine stärkere Dämpfung weisen einige leichte Elemente (vor allem Wasserstoff) auf, die Hauptbestandteile biologischer Proben sind. Kontrastmittel wie Deuterium-, Gadolinium- oder Lithiumverbindungen können verwendet werden, um den Kontrast auf ähnliche Weise zu verbessern, wie dies in der medizinischen Bildgebung der Fall ist, einschließlich Techniken wie optischer Bildgebung, Magnetresonanztomographie, Röntgenstrahlen und Positronen-Emissions-Tomographie. Für biologische Systeme werden Neutronenradiographie und Computertomographie zunehmend eingesetzt, um die Komplexität des unterirdischen Wurzelnetzwerks von Pflanzen, seine Wechselwirkung mit Böden und die Dynamik des Wasserflusses in situ zu untersuchen. Darüber hinaus wurden Bemühungen, Kontrastdetails in Tierproben, wie Weichteilen und Knochen, zu verstehen, untersucht. Dieses Manuskript konzentriert sich auf die Fortschritte in der Neutronenbiobildgebung, wie z. B. Probenvorbereitung, Instrumentierung, Datenerfassungsstrategie und Datenanalyse mit der CG-1D-Neutronenbildgebungs-Beamline des Hochflussisotopenreaktors. Die genannten Fähigkeiten werden anhand einer Auswahl von Beispielen aus der Pflanzenphysiologie (krautiges Pflanzen-Wurzel-Boden-System) und biomedizinischen Anwendungen (Ratten-Femur und Mauslunge) veranschaulicht.
Das Prinzip der Neutronenradiographie (nR) beruht auf der Abschwächung von Neutronen durch die Materie, die sie durchqueren. Im Gegensatz zu Röntgenstrahlen, die von der Elektronenwolke eines Atoms gestreut werden, können Neutronen von seinem Atomkern absorbiert oder gestreut werden. Neutronen sind empfindlich gegenüber leichten Elementen wie Wasserstoff (H) und können daher zur Röntgenaufnahme biologischer Anwendungen wie Tier 1,2,3,4,5,6,7 oder menschlichem Gewebe 8,9 und unterirdischen Boden-/Wurzelsystemen 10,11,12,13,14 verwendet werden ,15. Die Neutronenbildgebung ist eine komplementäre Technik zur Röntgenbildgebung, die in der Lage ist, schwere Elemente zu detektieren16,17,18. Die dämpfungsbasierte nR wird durch die linearen Dämpfungskoeffizienten der Materialien innerhalb der Probe und durch die Dicke der Probe bestimmt, wie sie durch das Beer-Lambert-Gesetz beschrieben wird, das besagt, dass der transmittierte Strahl direkt proportional zur Materialmenge und der Weglänge durch das Material ist. Somit kann der Transmissionsgrad T wie folgt berechnet werden:
(1)
wobei I 0 bzw. I die einfallenden bzw. transmittierten Strahlintensitäten sind; μ und x sind der lineare Abschwächungskoeffizient bzw. die Dicke einer homogenen Probe. Der Dämpfungskoeffizient μ ist gegeben durch:
(2)
wobei σ der Neutronendämpfungsquerschnitt der Probe (sowohl Streuung als auch Absorption), ρ ihre Dichte, NA die Avogadro-Zahl und M ihre Molmasse ist.
Der Kontrast in der Radiographie biologischer Proben mit niederenergetischen Neutronen (d. h. Energien unter 0,5 eV) ist hauptsächlich auf eine Änderung der Dichte von H (für eine feste Probendicke) zurückzuführen. Dies liegt an der Wahrscheinlichkeit der Wechselwirkung eines Neutrons mit dem H-Kern, die größer ist als bei anderen Kernen, die in biologischen Proben vorhanden sind, und an der Tatsache, dass die Dichte des H-Atoms von größter Bedeutung ist, da es das am häufigsten vorkommende Atom in biologischen Proben ist.
Seit ihren Anfängen werden die nR- und Neutronen-Computertomographie (nCT) in großem Umfang für Materialien und technische Anwendungen eingesetzt 19,20,21,22,23. Die ersten Demonstrationsexperimente zur Neutronenempfindlichkeit gegenüber H in biologischen Proben begannen Mitte der 1950er Jahre24 mit den Messungen von Pflanzenproben. Die Arbeit wurde bis in die 1960er Jahre fortgesetzt, z. B. mit der Röntgenaufnahme einer menschlichen Brust25 oder einer Ratte26, in der der Einsatz von Kontrastmitteln wie Gadoliniumoxid (Gd2O3) erforscht wurde. Darüber hinaus wurde die Hypothese aufgestellt, dass der Kontrast zwischen menschlichem Tumorgewebe und normalem Gewebe auf einen lokalen Anstieg des H-Gehalts zurückzuführen ist. Während dieser ersten Versuche kam man zu dem Schluss, dass ein erhöhter Neutronenfluss und eine höhere räumliche Auflösung die Qualität von nR verbessern und wahrscheinlich seine Popularität als ergänzende Technik für industrielle oder biomedizinische Anwendungen erhöhen würden. Die jüngsten Studien umfassen nR- und nCT-Messungen an Krebsgewebeproben1 und Schnitten tierischer Organe 2,3,27 für biomedizinische und forensische Anwendungen.
Der High Flux Isotope Reactor (HFIR) befindet sich am Oak Ridge National Laboratory, Oak Ridge, Tennessee, und ist eine leistungsstarke Neutronenquelle, die Neutronen durch Spaltreaktion erzeugt. Diese Neutronen haben Energien in der Größenordnung von 2 MeV und werden im Reaktorbecken durch kinetische Reaktionen mit schwerem Wasser “gekühlt”, um Energien in der Größenordnung von 100-300 eV zu erreichen. Die Optimierung eines Neutronenexperiments, ob Streuung oder Bildgebung, beginnt mit dem Verständnis der Eigenschaften der Neutronenquelle und der Strahllinie wie ihrer Strahlintensität, Energieverteilung und der Wirkung des Hintergrunds (schnelle Neutronen, verzögerte Neutronen, Gammastrahlen). In der HFIR-Kälteleithalle, in der sich die abbildende Beamline befindet, werden die Neutronen durch kinetische Wechselwirkungen mit einem flüssigen H-Moderator weiter “gekühlt”. Sie werden dann in einem gekrümmten Leitsystem von der Sichtlinie der Quelle weg transportiert, wodurch schnelle Neutronen und Gamma-Verschmutzung eliminiert werden. Wie in Abbildung 1 dargestellt, ist die CG-1D-Neutronenabbildungsstrahllinie28,29 auf einer kalten Führung platziert, was bedeutet, dass der Neutronenenergiebereich von einigen meV bis zu einigen Dutzend eV variiert (in diesem Fall reicht die entsprechende nutzbare Neutronenwellenlänge von 0,8 bis 10 Å) mit einem Fluss im Bereich von 107 n/(cm2∙s) an der Probenposition. Ein motorisiertes Blenden-/Diffusorsystem definiert die Lochblendengeometrie des Bildgebungsinstruments. Neutronen legen eine Strecke von 6,59 m in einem mit Helium (He) gefüllten Flugrohr mit Aluminiumfenstern (Al) an jedem Ende zurück. Flugröhren werden verwendet, um Neutronen zu transportieren und gleichzeitig die Luftstreuung so zu begrenzen, dass der Verlust der Strahlintensität minimal ist. Für die in diesem Manuskript beschriebenen Messungen besteht der Diffusor aus einem 1 mm dicken 50 nm Aluminiumoxid (Al2O3) Nanopulver, das in einem Al-Behälter eingeschlossen ist. Der Diffusor reduziert die vom Neutronenleiter ausgehenden Strahlartefakte (die durch die Lochblendengeometrie einer abbildenden Strahllinie vergrößert werden), da sonst scharfe horizontale und vertikale Intensitätsschwankungen im Röntgenbild sichtbar sind und die Normalisierung der Daten schwierig wird. Für die hier dargestellten Experimente werden Neutronen mit einem 25 μm dicken Lithium-6-Fluorid/Zinksulfid-Phosphor (6LiF/ZnS:Ag) in Licht umgewandelt.
Die Kollimationsoptimierung hängt von der Position der Probe zum Detektor, der erforderlichen räumlichen Auflösung und der Erfassungszeit ab. Wenn die Probe einige Zentimeter vom Szintillator entfernt sitzt, führen hohe Kollimationen (L/D über 800, wobei L der Abstand von der Lochblende des Durchmessers D und des Detektors ist) zu einer besseren räumlichen Auflösung auf Kosten des Neutronenflusses. Eine niedrige Kollimation (L/D unter 800) ist für dynamische In-situ-Studien vorzuziehen, bei denen die Zeitauflösung die räumliche Auflösung überwiegt. Für die in diesem Manuskript beschriebenen Messungen betrug L/D und räumliche Auflösung etwa 355 bzw. 75 μm. Die zeitliche Auflösung variierte je nach Signal-Rausch-Verhältnis (SNR). Die Probe wurde so nah wie möglich am Szintillator positioniert, um geometrische Verzerrungen wie Unschärfen zu reduzieren. Translations- und Rotationstische stehen zur Verfügung, um die Probe in die Nähe der Detektoren zu bringen und eine Computertomographie (CT) durchzuführen. CG-1D bietet drei Arten von Detektoren: ein ladungsgekoppeltes Gerät (CCD) mit 2048 Pixeln x 2048 Pixeln mit einem Pixelabstand von 13,5 μm, einen wissenschaftlichen komplementären Metalloxid-Halbleiter (sCMOS)-Detektor mit 2560 Pixeln x 2160 Pixeln mit einem Pixelabstand von 6,5 μm und einen Micro-Channel-Plate-Detektor (MCP)30,31 mit 512 Pixeln x 512 Pixeln mit einer Pixelgröße von 55 μm. Gestreute Neutronen werden mit ~5 mm dickem Borkautschuk absorbiert, um den Detektorchip vor Neutronen zu schützen. Diese Absorption erzeugt Gammastrahlen, die durch Blei (Pb) gestoppt werden können, das zwischen dem Borkautschuk und dem Detektor platziert wird. Jeder Detektor ist für ein anderes Sichtfeld (FOV) sowie für räumliche und zeitliche Auflösungen optimiert. Für die Messungen des Oberschenkelknochens der Ratte und der Lunge der Maus wurde der CCD-Detektor aufgrund seines großen Sichtfelds (~ 7 cm x 7 cm) und seiner angemessenen räumlichen Auflösung von ca. 75 μm verwendet. Die nCT des Wurzel-Boden-Systems der Pflanze wurde mit dem sCMOS durchgeführt, da das Ziel darin bestand, nCTs so schnell wie möglich auf Kosten des FOV (das auf ~ 5 cm x 4,2 cm begrenzt war) zu erhalten. Die räumliche Auflösung hat also offensichtlich gelitten. In diesen Detektoren werden Neutronen zu Detektionszwecken entweder in Licht oder ein Alpha-Teilchen umgewandelt. Die Drehung der Probe um ihre vertikale Achse und die Aufnahme von Röntgenaufnahmen bei aufeinanderfolgenden Drehwinkeln ermöglicht die Aufnahme von nCT. Das 3-dimensionale volumetrisch gerenderte Modell der untersuchten Probe wird durch die Verwendung des hauseigenen iMARS3D-Python-basierten Jupyter-FBP-Notebooks (Filtered-Back-Projection), pyMBIR oder einer kommerziellen Software erhalten, die im Folgenden beschrieben wird.
Schließlich werden Neutronen, die nicht mit der Probe oder dem Detektor interagiert haben, in einer Strahlstoppposition ca. 1 m stromabwärts des Detektorsystems gesammelt, um das Hintergrundrauschen zu minimieren. Der Strahlanschlag CG-1D ist 0,75 m breit, 0,5 m hoch und 35 mm dick und besteht aus B4C in Epoxidharz. Der Strahlanschlag ist mit 10 mm 95% angereichertem Lithiumcarbonat (6 Li2CO3) in einem feuerfesten Epoxidharz verstärkt, wo der Neutronenstrahl auftrifft, mit einem Hohlraum, der mit 6Li, Blei (Pb) und Stahl ausgekleidet ist, um die hohe Rate der sekundären Gammastrahlen einzudämmen. Der Trägeranschlag wird direkt an der Stahlabschirmwand des Strahlrohrs befestigt. Ein Foto des CG-1D-Strahlrohrs ist in Abbildung 2 dargestellt.
Es wurden jeweils drei Rekonstruktionsprogramme verwendet, um die drei experimentellen Daten in 3D zu rekonstruieren. Die Rekonstruktion der Lungenprobe der Maus wurde mit Octopus32 durchgeführt, einer kommerziellen Rekonstruktionssoftware, die FBP verwendet. Die Octopus-Software befindet sich auf einem Server-PC und kann verwendet werden, um die an der Beamline gesammelten Daten zu rekonstruieren. Eine Rekonstruktionssoftware namens iMARS3D ist unter CG-1D verfügbar. Es basiert auf dem Open-Source-Code TomoPY33 mit zusätzlichen Funktionen wie automatischer Neigungskorrektur, Nachbearbeitungsfiltern usw. iMARS3D umfasst die Vorverarbeitung der Daten (Subtraktion des Hintergrunds und des Rauschens), das Zuschneiden, die Medianfilterung (zur Korrektur von Gammaschlägen und toten Pixeln), die automatische Korrektur von Schwankungen der Strahlintensität und die Korrektur der Probenneigung. Sobald Sinogramme erstellt sind, ist eine weitere Datenverarbeitung wie das Entfernen von Ringartefakten und das Glätten möglich. Die verschiedenen Schritte der Rekonstruktion werden auf dem Analyseserver gespeichert (und später in den freigegebenen Vorschlagsordner verschoben), während die endgültigen 2D-Schichten sofort im freigegebenen Vorschlagsordner gespeichert werden. Der Oberschenkelknochen der Ratte wurde mit iMARS3D rekonstruiert. Die Pflanzenwurzel-/Bodenprobe wurde durch Medianfilterung der Daten mit TomoPY vorverarbeitet, gefolgt von einer Neigungsachsenkorrektur mit der SciPy-Bibliothek von Python. Die Rekonstruktion wurde unter Verwendung eines intern entwickelten Python-Pakets mit dem Namen pyMBIR (erstellt unter Verwendung von Kernen aus der ASTRA-Toolbox34) durchgeführt, das eine Reihe von tomographischen Algorithmen vom Basis-FBP bis hin zu fortgeschrittenen modellbasierten iterativen Rekonstruktionstechniken35 implementiert, die qualitativ hochwertige Rekonstruktionen aus extrem spärlichen und verrauschten Neutronendatensätzen erhalten können. Alle gerenderten Volumina, die auf den oben genannten Rekonstruktionswerkzeugen basieren, werden im Dämpfungskontrast dargestellt. Die gesamte Visualisierung erfolgte mit dem kommerziellen Visualisierungs-, Segmentierungs- und Datenanalyse-Softwarepaket AMIRA36.
Ziel dieses Manuskripts ist es, das Verfahren der Neutronenbildgebung (nR und nCT) an der HFIR CG-1D-Beamline zu demonstrieren. Diese Studie veranschaulicht auch den aktuellen Stand der Technik der nR- und nCT-Fähigkeiten für biologische Proben, insbesondere für eine Mauslunge, einen Rattenknochen und pflanzliche Wurzel-/Bodensysteme. Die Mauslunge wurde ausgewählt, um die Komplementarität von Neutronen zur Messung des Lungengewebes zu veranschaulichen, während Röntgenstrahlen hauptsächlich empfindlich für Knochen sind. Die Knochenprobe, ein Oberschenkelknochen einer Ratte, war mit einem Titan (Ti)-Implantat versehen, was den Kontrast zwischen Knochen und Metall veranschaulichte und die Möglichkeit bot, die Grenzfläche zwischen Knochen und Metall zu sehen (die mit Röntgenstrahlen schwer zu messen ist, da Metalle sie stark abschwächen4). Schließlich veranschaulicht das Pflanzen-Wurzel-Wassersystem die dreidimensionale (3D) Fähigkeit von nCT, Wurzel-/Bodensysteme in situ zu messen. Darüber hinaus werden die Vor- und Nachteile der Verwendung von nR für biologische Proben aufgezeigt. Offensichtlich kann diese Methode sicher verwendet werden, um die Wasserdynamik in einem Pflanzen-Wurzel-System zu messen, kann aber aufgrund der mit der Strahlenbelastung verbundenen Risiken nicht als lebendes Tier- oder menschliches Bildgebungsverfahren betrachtet werden, so dass Studien entweder auf (tote) Mäuse oder pathologieähnliche Messungen beschränkt sind, bei denen beispielsweise eine Gewebeprobe von einem Patienten (Tier oder Mensch) entnommen und durch Fixierung präpariert wird, bevor sie in einem Neutronenstrahl gemessen wird.
Neutronenröntgen und CT biologischer Proben sind vielversprechende bildgebende Verfahren, die die Röntgenbildgebung oder die Magnetresonanztomographie ergänzen. Die entscheidenden Schritte bei der Durchführung eines Neutronenbildgebungsexperiments einer biologischen Probe beziehen sich auf ihre Präparation und ihren Einschluss an der Beamline. Die Optimierung eines Experiments wird durch die zu beantwortende wissenschaftliche Fragestellung vorangetrieben. Wenn die wissenschaftliche Fragestellung eine hohe räumliche Auflösung erfordert, um ein Phänomen zu beobachten, dann sind lange Aufnahmezeiten erforderlich, und der Nachteil der nCT (mit cm-großem Sichtfeld) besteht darin, dass es Stunden dauert, einen Scan durchzuführen. Dies ist vor allem auf den Unterschied im gesamten Neutronenfluss zurückzuführen, der an einem Reaktor im Vergleich zu einer Synchrotronquelle verfügbar ist, wo Röntgen-CT-Scans Sekunden bis Minuten für ein Sichtfeld von wenigen mm2 dauern können. Obwohl die Methode auf Ex-vivo-Gewebeproben angewendet werden kann, die von Tieren entnommen wurden, kann sie in vivo aufgrund des Strahlenexpositionsrisikos (z. B. Gammastrahlen, die von Neutronen erzeugt werden, und Neutronenwechselwirkungen mit den Atomen in der Probe) nicht auf lebende Tiere oder Menschen ausgedehnt werden. Es eignet sich jedoch gut für die Abbildung von Interaktionen zwischen Pflanzenwurzeln und Boden (Abbildung 7), wie z. B. der Dynamik der Wasseraufnahme.
Der Vorteil der schnellen nCT für die Pflanzendynamik ist die Empfindlichkeit gegenüber H im Wasser und das Fehlen von Strahlenschäden an der Pflanze, im Gegensatz zur Röntgen-CT. Darüber hinaus kann ein einzigartiger Kontrast durch die Verwendung von Neutronen in Knochen-Metall-Proben wie einem Ratten-Femur beobachtet werden, bei dem das Metall im Vergleich zum umgebenden Gewebe relativ transparent ist (Abbildung 5), wodurch möglicherweise Metallartefakte vermieden werden, die durch Röntgen-CT39 induziert werden. Tierische Gewebe, wie z. B. die Lunge der Maus (Abbildung 6), zeigen einen beeindruckenden Nachweis der Weichteilstruktur, da Neutronen empfindlich auf H reagieren, aber die räumliche Auflösung ist bei diesen Messungen etwas der limitierende Faktor. Der Kontrast wird durch die in biologischen Proben vorhandenen H-Atome gewährleistet19,39.
Mit den Fortschritten neuartiger Techniken wie der Neutronengitterinterferometrie und der Verbesserung der räumlichen Auflösung (einige Mikrometer wurden kürzlich berichtet42,43) könnte die Neutronenbildgebung noch neue Kontrastmechanismen für biologische Gewebe mit verbesserter räumlicher Auflösung bieten. Die Erforschung von Neutronen höherer Energie (um die Messung dicker Proben zu ermöglichen) verspricht auch die Möglichkeit, größere Ausschnitte eines tierischen Gewebes wie einer intakten Maus zu vermessen und bietet damit noch neue Möglichkeiten für die biomedizinische Forschung.
The authors have nothing to disclose.
Ein Teil dieser Forschung nutzte Ressourcen des Hochfluss-Isotopenreaktors, der vom ORNL betrieben und vom U.S. Department of Energy, Office of Science, User Facilities, im Rahmen des Vertrags DE-AC05-00OR22725 mit UT-Battelle, LLC gesponsert wurde. Ein Teil dieser Forschung wurde vom ORNL durch das Eugene Wigner Distinguished Staff Fellowship-Programm unterstützt. Diese Forschung wurde auch vom DOE Office of Science, Office of Biological and Environmental Research gesponsert. Die Oberschenkelproben von Ratten stammen aus Experimenten, die in Zusammenarbeit mit Dr. Rick Sumner am Rush University Medical Center mit Mitteln des NIH (R01AR066562) und des Orthopedic Research and Education Foundation-Smith and Nephew Award durchgeführt wurden. Das Team bedankt sich bei den HFIR-Support-Teams, die die Nutzung der Strahlrohre zur Neutronenstreuung ermöglichen.
Aluminum containers | custom | Made from aluminum plates or tubing (alternate is quartz), plant and mouse sample | |
Aluminum foil | Fisher | 01-213-100 | Mouse lung sample containment |
Deionized water or deuterium oxide | Water or D2O can be used to enhance contrast, plant sample | ||
Ethanol | Fisher | 04-355-223 | Mouse lung sample |
Gauze sponges | CardinalHealth | Fully submerged in phosphate-buffered saline (PBS) and used to wrap samples, rat femur sample | |
Growth chamber | Conviron | A1000 | Any growth chamber or greenhouse with controlled conditions would work, plant sample |
Laboratory balance | Weighing plant system can be used to measure actual water content in the soils, plant sample | ||
Pure silica sand | US Silica Co. | Flint#13 | Pure SiO2 provides low neutron attenuation compared to soils, plant sample |
Sprague-Dawley Rats | Harlan | Order Code: 002-US | Rat femur sample |
Titanium Rod | Goodfellow | TI007905 | Rat femur sample |