
Citation: Kavita Singh, Andy Nong, Mark Feeley, Hing Man Chan. The use of biomonitoring equivalents for interpreting blood concentrations in population studies: a case for polychlorinated biphenyls[J]. AIMS Environmental Science, 2015, 2(1): 21-41. doi: 10.3934/environsci.2015.1.21
[1] | Chichia Chiu, Jui-Ling Yu . An optimal adaptive time-stepping scheme for solving reaction-diffusion-chemotaxis systems. Mathematical Biosciences and Engineering, 2007, 4(2): 187-203. doi: 10.3934/mbe.2007.4.187 |
[2] | Fabrizio Clarelli, Roberto Natalini . A pressure model of immune response to mycobacterium tuberculosis infection in several space dimensions. Mathematical Biosciences and Engineering, 2010, 7(2): 277-300. doi: 10.3934/mbe.2010.7.277 |
[3] | Chunlai Mu, Weirun Tao . Stabilization and pattern formation in chemotaxis models with acceleration and logistic source. Mathematical Biosciences and Engineering, 2023, 20(2): 2011-2038. doi: 10.3934/mbe.2023093 |
[4] | Yue Xing, Weihua Jiang, Xun Cao . Multi-stable and spatiotemporal staggered patterns in a predator-prey model with predator-taxis and delay. Mathematical Biosciences and Engineering, 2023, 20(10): 18413-18444. doi: 10.3934/mbe.2023818 |
[5] | Ali Moussaoui, Vitaly Volpert . The impact of immune cell interactions on virus quasi-species formation. Mathematical Biosciences and Engineering, 2024, 21(11): 7530-7553. doi: 10.3934/mbe.2024331 |
[6] | Fiona R. Macfarlane, Mark A. J. Chaplain, Tommaso Lorenzi . A hybrid discrete-continuum approach to model Turing pattern formation. Mathematical Biosciences and Engineering, 2020, 17(6): 7442-7479. doi: 10.3934/mbe.2020381 |
[7] | Tingfu Feng, Leyun Wu . Global dynamics and pattern formation for predator-prey system with density-dependent motion. Mathematical Biosciences and Engineering, 2023, 20(2): 2296-2320. doi: 10.3934/mbe.2023108 |
[8] | Marcello Delitala, Tommaso Lorenzi . Emergence of spatial patterns in a mathematical model for the co-culture dynamics of epithelial-like and mesenchymal-like cells. Mathematical Biosciences and Engineering, 2017, 14(1): 79-93. doi: 10.3934/mbe.2017006 |
[9] | Ali Moussaoui, Vitaly Volpert . The influence of immune cells on the existence of virus quasi-species. Mathematical Biosciences and Engineering, 2023, 20(9): 15942-15961. doi: 10.3934/mbe.2023710 |
[10] | Kalyan Manna, Malay Banerjee . Stability of Hopf-bifurcating limit cycles in a diffusion-driven prey-predator system with Allee effect and time delay. Mathematical Biosciences and Engineering, 2019, 16(4): 2411-2446. doi: 10.3934/mbe.2019121 |
Atrial fibrillation (AF) is a common cardiac arrhythmia, during which the regular and orderly electrical activity in the atrium is replaced by rapid and chaotic heartbeats. The ineffective contraction of the atrium in the post-diastolic phase during AF can induce thrombosis, while the shed thrombus can cause ischemic stroke when flowing to the brain. The presence of AF is independently associated with a 5-fold increased risk of stroke [1]. The left atrial appendage (LAA) is a finger- or stump-like extension of the left atrium (LA) with lobes that may harbor up to 90% of thrombi that occur in patients with AF [2]. The remodeling process associated with AF causes the LAA to function as a static pouch, predisposing to stagnation and thrombosis. It also should be noticed that the proportion of thrombus in LAA as the origin of cardiogenic emboli is as high as 100% for non-vascular AF patients [2].
Prevention strategies of ischemic stroke caused by AF include oral anticoagulation therapy (i.e., warfarin) and minimally invasive surgeries to exclude LAA [3,4,5]. Before anticoagulation therapy, the risk of stroke is stratified according to the CHA2DS2-VASc score which refers to age, sex, stroke history and so on. A higher score means a higher risk of stroke [6]. CHA2DS2-VASc score ≡2 means a patient at high risk and needs anticoagulation therapy. However, the scoring system still has some limitations. For example, a young male with AF whose CHA2DS2-VASc score is 0, still could develop ischemic stroke. Noticed that some clinical studies have shown that LA and LAA hemodynamic information which is difficult to obtain through current medical technology, can improve stroke risk stratification [6], thus computational fluid dynamics (CFD) is widely applied to capture blood flow characteristics [7,8,9]. The study of the LA wall motion indicates that LAA is ineffective and acts as a stasis blood reservoir during AF [9]. The blood evacuation rate in LAA decreases when AF occurs which might increase thrombosis risk [8]. Comparing the lack of active atrial contraction, the occurrence of high-frequency fibrillation might have a larger impact on the blood flow stagnation [10]. We find that only a few experimental and numerical research study the flow patterns in LA [11] and hemodynamic changes induced by different LAA morphology or AF [7,8,9,12]. At the same time, the direct simulation of thrombosis in the LA, especially LAA, has yet to be produced.
The purpose of this study was to simulate the thrombosis in LA using advanced computational modeling analyses. To this end, the patient-specific geometry was reconstructed and a thrombus model was combined with constant inlet pressure and realistic outflow. To the author*s knowledge, the present study is the first to use the CFD method to directly model thrombosis in patient-specific LA. The magnitudes of thrombosis variables, blood flow distribution and wall shear stress (WSS) related indices were analyzed to explore the impact of AF on the hemodynamics.
Computed tomography (CT) images of LA and four pulmonary veins (PVs) were obtained from an adult patient with non-valvular AF at Sir Run Run Shaw Hospital, Zhejiang University School of Medicine (Hangzhou, China). The study was conducted in accordance with the ethical standards of the Local Institution Review Board and with the 1964 Helsinki declaration and its later amendments. This project was approved by the Local Bioethics Committee. All patients have given their written informed consent before the study. The patient-specific LA geometry was constructed from the CT images by using Mimics 19.0 (Materialise, Belgium) (Figure 1-A). All the geometric boundaries of inlets and outlet were cropped to get a flat surface in GeoMagic Studio (GeoMagic Inc, USA) (Figure 1-B). The fluid domain was meshed using commercial software ANSYS-ICEM (ANSYS Inc, Canonsburg, USA), while five prism layers of unstructured tetrahedral meshes were created close to the atrial wall. The computational domain had approximately 2,420,000 elements, and the unstructured tetrahedral mesh had a maximum size of 7 mm. Finer meshes with more than 4,900,000 elements were created to assess the mesh sensitivity. The differences in peak wall shear stress and averaged wall shear stress between meshes were both less than 2%.
The blood was considered to be an incompressible and non-Newtonian fluid with a constant density of 1080 kg/m3 [13]. The LA wall was assumed to be rigid with no-slip conditions. The deformation of the LA should be considered for precise results in most heart flow simulations. However, the focus of the present study is on the application of a thrombus model in LA to simulate thrombosis. In addition, the effect of the atrial kick absence, a feature of the AF, on thrombosis is also explored. To build a fully developed blood flow environment before thrombus model was employed [13], Quemada model was used to simulate the shear-thinning characteristics of the blood [14] for two complete cardiac cycles (each cardiac cycle of 0.8 s duration including 37.5% ventricular diastole and 62.5% ventricular systole). The thrombus model was developed based on a hemodynamics-based model presented by Menichini, et al. [15], which simulated thrombosis by solving transport equations of flow residence time (RT), activated platelets (AP), resting platelets (RP), bound platelets (BP) and coagulation (C). The thrombus model considers the effect of the thrombus on the blood viscosity as well. Thrombus is identified through the local BP concentration. Details of the mathematical equations can be found in the study of Menichini, et al. [15], but there are some important changes described below.
The momentum of the flow is controlled by the simplified Navier-Stokes equation. We ignore the frictional force of thrombus on the fluid. For the transport of AP and RP, a non-normalized source item Si is applied which is the summary of $r_{1}$ and $r_{2}$. BP and C have been normalized by their initial values for avoiding numerical ill-conditioning. Besides, the judgment threshold value of the local time average shear strain rate is 1 s-1, which determines that the initial concentration of C is 100 or 0 nmol/L.
$ \rho \frac{\partial \boldsymbol{u}}{\partial t}+\rho \boldsymbol{u} \cdot \nabla \boldsymbol{u} = -\nabla p+\nabla \cdot\left(\mu\left(\nabla \boldsymbol{u}+(\nabla \boldsymbol{u})^{T}\right)\right) $ | (1) |
$ r_{1} = k_{1} \cdot \mathrm{AP} \cdot \mathrm{RP} \quad r_{2} = k_{2} \cdot \mathrm{RP} \cdot \mathrm{RRT} $ | (2) |
$k_{1}$ and $k_{2}$ are the kinetic constants [16]. AP and RP are the concentrations of the platelets respectively. RRT is the standardized RT which is normalized by cardiac cycle time.
The difficulty in obtaining patient-specific flowrate data in the outlet, mitral valve (MV), would be revealed when using medical measurement. The flowrate curves across the MV for the sinus rhythm (SR) cases were derived from the international regulation, which included E-wave and A-wave representing the first and the second phases, respectively [17]. The transmitral Doppler echocardiographic analysis of AF patients shows that A wave which is caused by the atrial kick in late diastole is absent in the flowrate curve [10,18]. The flowrate curves for AF cases were differed by dismissing the A wave [7]. For the four pulmonary veins inlets, constant pressure of 10 mmHg was adopted. Flowrate curves for SR and AF cases are shown in Figure 1-C, the peaks of two curves are nearly 46.9 ml/s. The maximum Reynolds number of all cases was found lower than 3000 when estimated on the basis of the peak flow rate and the MV diameter. Therefore, the blood flow was assumed as a laminar flow. Simulations were performed on ANSYS Workbench 16.1 (ANSYS Inc, USA). The time step in this transient study was 10 ms and all the cases were simulated for 20 cardiac cycles to get feasible results and the twentieth cycle data was post-processed.
Before applied to LA geometry, the thrombus model was tested in a backward-facing step (BFS) model (Figure 1-D). Quemada model with a 30% hematocrit was selected to simulate experimental blood [19], while in the LA model, the hematocrit of human blood was determined to be nearly 45% [14]. The boundary conditions were adopted from the experimental test, with 0.76 L/min at the inlets and zero pressure at the outlet [19]. Except for the bulk shear rate threshold, all other parameters were the same as in the LA model. The time step in this verified study was 0.05 s and the calculation time was 50 s. Verification results were compared with the experimental study [19] and the simulation results [15].
Figure 2 shows the thrombosis prediction results in the BFS model and they are compared with the published results. The thrombus initially appears in the downstream of the step then grows continuously. When the height of the thrombus reaches the step height (2.5 mm) and the maximum width is close to that of the step, thrombus growth almost stops in the radial direction. The thrombus growth turns slow after 14 s and almost stops growing when the total thrombus length reaches about 8.8 times the height of the step, with a maximum length of 22 mm. These results are consistent with the simulation data of Menichini, et al. [15] and match well with the experimental results [19] as we can use it for further study.
The velocity streamlines for two cases are illustrated in Figure 3. The initial streamlines seeds are averagely launched from PVs. At the first peak systole (0.4 s), the streamline distributions in the two cases are similar, the velocities in the inlets (PVs) and MV outlet are higher than in other areas of the LA. During the later peak systole (0.7 s), the blood flow gets highly disordered in the recirculating area, and the maximum velocity appears nearby the left inferior pulmonary vein (LIPV) rather than the PVs and MV in the AF case. Besides, there is little blood flowing through the MV in the AF case while blood can in the SR case. Mentioned that there is little blood flowing out from the LAA as well. Table 1 shows the proportions of blood that pass through each inlet boundaries during the last cardiac cycle. During the later peak systole for AF case, the flow ratios of the four inlets were not calculated because part of the inlet blood returned. In the SR case, the blood flow crossing the LIPV accounts for nearly 30% of the outflow while the minority of the blood (below 20%) goes through the right inferior pulmonary vein.
Inlet boundaries | 0.4 s | 0.7 s | ||
AF | SR | AF | SR | |
Left Superior Pulmonary Vein | 26.86% | 27.88% | - | 24.49% |
Left Inferior Pulmonary Vein | 30.27% | 30.56% | - | 29.02% |
Right Superior Pulmonary Vein | 26.14% | 25.35% | - | 27.20% |
Right Inferior Pulmonary Vein | 16.74% | 16.22% | - | 19.29% |
AF: atrial fibrillation; SR: sinus rhythm. |
The frictional force on the LA wall exerted by the blood flow is quite difficult to be obtained through current measurement technology but can be studied by the WSS. The distributions of WSS for two cases in peak systole phases are shown in Figure 4. Similar to the flow waveform, the distributions of WSS in two cases have little difference at the earlier peak systole. At the later peak systole, the maximum WSS values of two cases decrease greatly, especially in the AF case. Besides, the region around LIPV shows a bigger WSS than the other three inlets in the AF case (Figure 4-D), while areas around the four PVs have a similar distribution in the other three maps. The LAA has a small WSS value in all maps.
We choose BP as a thrombosis sign with a threshold of 200 nmol/L [20]. Figure 5 shows the thrombosis in the maximum cross-section of the LAA from 5.6 s with 1.6 s interval for AF case. The simulation data illustrates that thrombus forms since the seventh cardiac cycle. The thrombus first appears on the wall of the LAA, then grows inward and fills most space in the LAA at the end of the nineteenth cardiac cycle. The distributions and volume proportions (comparing with the total LA volume) of formed thrombus in the seventh and last cycles of two cases are shown in Figure 6. Thrombus only arises in the LAA with its proportions have increased tremendously for two cases from the eighth to the twentieth cycle. The difference of thrombosis ratio is clear between two cases at 16 s while AF case has a greater severe thrombosis with a 4.93% thrombus volume proportion.
The formation of BP is directly influenced by C, whose initial value is depending on the distribution of the average shear rate on the wall after the first two cardiac cycles. If the average shear rate is smaller than the judgment threshold, the initial value of C is 200 nmol/L [15]. In view of the importance of C initial value for outcomes, thrombosis in AF cases with judgment thresholds of 0.5 s-1 and 1.5 s-1 were calculated as well (Figure 7). For all cases with different thresholds, the thrombus is visible since the seventh cardiac cycle. In the end, the volume fraction of thrombus reaches 5.68% in the case with a threshold of 1.5 s-1, which is nearly twice that of the case with the judgment threshold of 0.5 s-1. Other than that, the relative thrombus fraction difference between cases with thresholds of 1.5 s-1 and 1 s-1 is 15.2%, which is smaller than that value between cases with thresholds of 1 s-1 and 0.5 s-1. Figure 7-D investigates the time-varying laws of thrombus proportions for three cases. The growth of thrombosis accelerates with the increase of judgment threshold and decelerates with the change of time. Furthermore, the difference of thrombus growth rates between cases with judgment thresholds of 1.5 s-1 and 1 s-1 is smaller than that between the latter and the case with a threshold of 0.5 s-1.
There are millions of patients suffering from non-valvular AF in China, with a 20% fatality rate and a 60% disability rate [21]. The management of AF in China is unproductive and accompanied by a high risk of mortality due to patients* low awareness and lack of proper treatment [22]. Clinical studies have shown that the usage of warfarin or new anticoagulants in AF patients with a high stroke risk can significantly reduce stroke events and improve patient outcomes [22]. Therefore, the risk assessment of thrombosis under AF conditions is crucial in the prevention of AF-related stroke. Hemodynamic characteristics of LA can improve stroke risk stratification [6], but some of them are difficult to perform high-precision measurement through non-invasive means in clinical work [7,23]. Thus we established a thrombus model to study the AF effect on the hemodynamics in the LA by numerical simulation and to build a fundament for further application of the CFD method to study thrombosis risk prediction [8].
In this study, we utilized the thrombus model to predict thrombosis in LA with the consideration of non-Newtonian properties of blood for investigating the hemodynamics difference between AF and SR cases to the same LA geometry while most previous studies just simplify blood as Newtonian flow [2,8,9,10,11,12]. The thrombus model was calibrated in the BFS model before applied to the LA model because of the absence of the thrombosis experiment in a spherical tube [19]. The results of the benchmark study were well matched with the experiment results.
Observing the streamlines in the second peak systole, in the AF case, there is little blood flowing through the MV due to the lack of atrial kick but toward the LIPV vice versa, while in the SR case, blood can pass through the MV. It will increase blood*s impact force on the wall which is consistent with the WSS distribution. The blood flow gets highly disordered in the recirculating area which means the formed thrombus is easier to separate from the wall and then move with the blood flow. Comparing to the other areas, the velocity and WSS are small in the LAA dramatically, indicating that the LAA region is susceptible to thrombosis and may be worthy of attention in the thrombosis risk assessment as previous studies have demonstrated that low WSS distributions identify regions more prone to thrombus formation [24].
There are many variables in the thrombus model, the distribution of BP predicts the thrombosis area directly while RT, C, AP, and RP are intermediate variables. C has the closest tie with the BP when compared to RT, AP and RP, as the formation of C decides the initiation of BP, in addition, BP will promote the formation of C as well. The impact of C on the thrombosis was also studied in this paper. Observation of BP above the judgment threshold shows that thrombosis mainly happens in the LAA which is consistent with the clinical statistics [2]. Thrombosis begins in the wall because of the tissue factor exposure and interaction of cells and fibrinogen there which are enhanced with a low shear rate [25]. The simulation results show that the thrombus distributions are distinguished between the different cardiac rhythms. In the AF case, thrombosis is more serious than the SR case with a ratio difference close to 18%. It also means that the impact of the occurrence of high-frequency fibrillation cannot be ignored in the simulation of AF [10]. C is an important variable that influences the initiation of thrombosis, and its initial value depends on the underlined threshold of the time-averaged shear stress rate from previous simulation results. From the sensitivity analysis of the threshold, we have found that a higher threshold means more LAA area has non zero initial C and a faster thrombosis growth rate. It can be deduced that the threshold we have chosen possessed certain rationality, noticing that the difference between thresholds of 1.5 s-1 and 1 s-1 is much smaller than that of cases with thresholds of 1 s-1 and 0.5 s-1.
Treating the LA wall as a rigid wall may be a key limitation of our simulation in this paper. In this study, we investigated the impact of the atrial kick absence by changing the outlet condition as previous researchers did [7]. For AF, there are more characters that will influence the hemodynamics in the LA, such as the high-frequency fibrillation of the LA wall. To the author*s knowledge, there are two ways to rebuild the deformation of LA, one is to set the atrial wall as the fibre-reinforced hyperplastic material and consider the fluid-structure interaction effect [9,26], while the other way is to load the deformation laws of the LA and the LAA, which can be derived from the MR or CT slice-images of the volunteers to the geometry directly [8,10,12,27]. We will consider the deformation of the LA and the LAA for the thrombosis prediction by using these two ways in our future study.
In addition, the boundary conditions are not strictly patient-specific due to the absence of clinical measurement data. The aim of the present study was the application of thrombus model to simulate the thrombosis in LA rather than the accurate thrombosis prediction in special patients. Our further study will focus on the latter strategy and PC-MRI measurement data would be used as boundary conditions for further improvement of the simulation.
LAA has stronger active contractility than LA, due to its thick pectinate muscle [28]. For healthy people, the contraction of LAA will empty the inside blood while AF patients vice versa [3,28]. The neglect of the influence of the LAA contraction may lead to thrombosis in the SR patient cases. More studies will be displayed on the LAA contraction and geometry in our future research to evaluate their impacts on thrombosis.
Thrombosis is a complex process involves the conversion of many substances [29]. In this paper, we only consider the effects of four substances. There is no doubt that other substances also have impacts on the viscosity property of the blood flow and thrombosis which should be further studied.
This is the first study that directly investigates the thrombosis in the LA with the discussion of the hemodynamics changes. The thrombus model used in this paper was calibrated in the BFS model and in good agreement with the experimental data. The hemodynamics of the AF cases and SR cases are compared, showing that the hemodynamics in LA is greatly affected by AF. Thrombosis observed in this study only occurs in the LAA which is consistent with clinical findings, thrombus first forms on the LAA wall then grows inward until filling up the inner space of LAA. AF can promote the degree of thrombosis. In the future, numerical simulation of the thrombosis and clinical statistics can be combined to improve the risk assessment mechanism of thrombosis and stroke.
All authors declare no conflicts of interest in this paper.
[1] | U.S. Environmental Protection Agency. Terminology Sevices. Available from: http://ofmpub.epa.gov/sor_internet/registry/termreg/home/overview/home.do. |
[2] | U.S. Environmental Protection Agency. Exposure Factors Handbook Glossary of Terms, 2011. Available from: http://www.epa.gov/ncea/efh/report.html. |
[3] | Agency for Toxic Substances and Disease Registry. Toxic Substances Portal. Minimal Risk Levels (MRLs), 2014. Available from: http://www.atsdr.cdc.gov/mrls/index.asp. |
[4] | Health Canada.Health-Based Tolerable Daily Intakes/Concentrations and Tumorigenic Doses/Concentrations for Priority Substances, 1996. Available from: http://www.hc-sc.gc.ca/ewh-semt/pubs/contaminants/hbct-jact/index-eng.php. |
[5] | Agency for Toxic Substances and Disease Registry (ATSDR) (2000) Toxicological profile for Polychlorinated Biphenyls (PCBs). Atlanta, GA: U.S. Department of Health and Human Services, Public Health Service. |
[6] | U.S. Environmental Protection Agency (EPA). Approaches for the Application of Physiologically Based Pharmacokinetic (PBPK) Models and Supporting Data in Risk Assessment. National Center for Environmental Assessment, Washington, DC; EPA/600/R-05/043F, 2006. Available from: National Technical Information Service, Springfield, VA, and online at http://epa.gov/ncea. |
[7] |
Hays SM, Becker RA, Leung HW, et al. (2007) Biomonitoring equivalents: a screening approach for interpreting biomonitoring results from a public health risk perspective. Regul Toxicol Pharmacol 47: 96-109. doi: 10.1016/j.yrtph.2006.08.004
![]() |
[8] |
Hays SM, Aylward LL, LaKind JS, et al. (2008) Guidelines for the derivation of biomonitoring equivalents: report from the biomonitoring equivalents expert workshop. Regul Toxicol Pharmacol 51: S4-15. doi: 10.1016/j.yrtph.2008.05.004
![]() |
[9] |
Krishnan K, Adamou T, Aylward LL, et al. (2011) Biomonitoring equivalents for 2,2',4,4',5-pentabromodiphenylether (PBDE-99). Regul Toxicol Pharmacol 60: 165-171. doi: 10.1016/j.yrtph.2011.03.011
![]() |
[10] |
Kirman CR, Aylward LL, Hays SM, et al. (2011) Biomonitoring equivalents for DDT/DDE. Regul Toxicol Pharmacol 60: 172-180. doi: 10.1016/j.yrtph.2011.03.012
![]() |
[11] |
Aylward LL, Barton HA, Hays SM. (2008) Biomonitoring Equivalents (BE) dossier for toluene (CAS No. 108-88-3). Regul Toxicol Pharmacol 51: S27-36. doi: 10.1016/j.yrtph.2008.05.009
![]() |
[12] |
Aylward LL, Hays SM, Gagné M, et al. (2010) Biomonitoring equivalents for hexachlorobenzene. Regul Toxicol Pharmacol 58: 25-32. doi: 10.1016/j.yrtph.2010.06.003
![]() |
[13] |
Aylward LL, Hays SM. (2008) Biomonitoring Equivalents (BE) dossier for 2,4-dichlorophenoxyacetic acid (2,4-D) (CAS No. 94-75-7). Regul Toxicol Pharmacol 51: S37-48. doi: 10.1016/j.yrtph.2008.05.006
![]() |
[14] |
Porta M, Zumeta E. (2002) Implementing the Stockholm Treaty on Persistent Organic Pollutants. Occup Environ Med 59: 651-653. doi: 10.1136/oem.59.10.651
![]() |
[15] |
Breivik K, Alcock R, Li YF, et al. (2004) Primary sources of selected POPs: regional and global scale emission inventories. Environ Pollut 128: 3-16. doi: 10.1016/j.envpol.2003.08.031
![]() |
[16] | Breivik K, Sweetman A, Pacyna JM, et al. (2002) Towards a global historical emission inventory for selected PCB congeners—a mass balance approach 1. Global production and consumption. Sci Total Environ 290: 181-198. |
[17] | U.S. Environmental Protection Agency (EPA) (2003) Non-dioxin-like PCBs: Effects and consideration in ecological risk assessment. Experimental Toxicology Division National Health and Environmental Effects Research Laboratory Office of Research and Development (NCEA-C-1340. ERASC-003. June). |
[18] | IARC monographs on the evaluation of carcinogenic risks to humans. Overall evaluations of carcinogenicity: an updating of IARC Monographs volumes 1 to 42. Supplement 7, 1987. Available from: http://monographs.iarc.fr/ENG/Monographs/suppl7/index.php. |
[19] | Van Oostdam J, Donaldson SG, Feeley M, et al. (2005) Human health implications of environmental contaminants in Arctic Canada: a review. Sci Total Environ 351-352: 165-246. |
[20] |
Laird BD, Goncharov AB, Chan HM. (2013) Body burden of metals and persistent organic pollutants among Inuit in the Canadian Arctic. Environ Int 59: 33-40. doi: 10.1016/j.envint.2013.05.010
![]() |
[21] | Agency for Toxic Substances and Disease Registry (ATSDR) (2011) Addendum to the Toxicological profile for Polychlorinated Biphenyls (PCBs). Atlanta, GA: U.S. Department of Health and Human Services, Public Health Service. |
[22] | Health Canada. Report on Human Biomonitoring of Environmental Chemicals in Canada, 2010. Available from: http://www.hc-sc.gc.ca/ewh-semt/pubs/contaminants/chms-ecms/section8-eng.php#n8_3. |
[23] | U.S. Environmental Protection Agency (EPA). Integrated Risk Information System (IRIS). Aroclor 1016 (CASRN 12674-11-2), 1993. Available from: http://www.epa.gov/iris/subst/0462.htm. |
[24] | Joint FAO/WHO Expert Committee on Food Additives. Meeting (67th : 2006 : Rome, Italy) Evaluation of certain food additives and contaminants : sixty-seventh report of the Joint FAO/WHO Expert Committee on Food Additives, 2007. Available from: http://www.who.int/ipcs/publications/jecfa/reports/trs940.pdf. |
[25] | Centers for Disease Control and Prevention. Polychlorinated Biphenyls (PCBs). DHHS (NIOSH) Publication Number 86-111, 1986. Available: http://www.cdc.gov/niosh/docs/86-111/. |
[26] | French Agency for Food, Environmental and Occupational Health and Safety (ANSES). What are the critical blood concentration levels for PCBs?, 2013. Available from: https://www.anses.fr/en/content/what-are-critical-blood-concentration-levels-pcbs. |
[27] | Health Canada. Part II: Health Canada Toxicological Reference Values (TRVs) and Chemical-Specific Factors, Version 2.0, 2010. |
[28] | European Food Safety Authority (EFSA). Opinion of the scientific panel on contaminants in the food chain on a request from the Commission related to the presence of non dioxin-like polychlorinated biphenyls (PCB) in feed and food, 2005. Available from: http://www.efsa.europa.eu/en/efsajournal/pub/284.htm. |
[29] | French Agency for Food, Environmental and Occupational Health and Safety (ANSES). Opinion of the French Food Safety Agency on the establishment of relevant maximum levels for non dioxin-like polychlorobiphenyls (NDL-PCB) in some foodstuffs, 2007. Available from: https://www.anses.fr/sites/default/files/documents/RCCP2006sa0305EN.pdf. |
[30] |
Weijs L, Roach AC, Yang RSH, et al. (2014) Lifetime PCB 153 bioaccumulation and pharmacokinetics in pilot whales: Bayesian population PBPK modeling and Markov chain Monte Carlo simulations. Chemosphere 94: 91-96. doi: 10.1016/j.chemosphere.2013.09.019
![]() |
[31] |
Weijs L, Covaci A, Yang RSH, et al. (2011) A non-invasive approach to study lifetime exposure and bioaccumulation of PCBs in protected marine mammals: PBPK modeling in harbor porpoises. Toxicol Appl Pharmacol 256: 136-145. doi: 10.1016/j.taap.2011.07.020
![]() |
[32] |
Weijs L, Yang RSH, Covaci A, et al. (2010) Physiologically based pharmacokinetic (PBPK) models for lifetime exposure to PCB 153 in male and female harbor porpoises (Phocoena phocoena): model development and evaluation. Environ Sci Technol 44: 7023-7030. doi: 10.1021/es101688h
![]() |
[33] |
Czub G, McLachlan MS. (2007) Influence of the temperature gradient in blubber on the bioaccumulation of persistent lipophilic organic chemicals in seals. Environ Toxicol Chem 26: 1600-1605. doi: 10.1897/06-548R.1
![]() |
[34] |
Klanjscek T, Nisbet RM, Caswell H, et al. (2007) A model for energetics and bioaccumulation in marine mammals with applications to the right whale. Ecol Appl 17: 2233-2250. doi: 10.1890/06-0426.1
![]() |
[35] |
Hickie B, Mackay D, De Koning J. (1999) Lifetime pharmacokinetic model for hydrophobic contaminants in marine mammals. Environ Toxicol Chem 18: 2622-2633. doi: 10.1002/etc.5620181132
![]() |
[36] |
Drouillard KG, Norstrom RJ, Fox GA, et al. (2003) Development and validation of a herring gull embryo toxicokinetic model for PCBs. Ecotoxicology 12: 55-68. doi: 10.1023/A:1022588913171
![]() |
[37] |
Drouillard KG, Fernie KJ, Smits JE, et al. (2001) Bioaccumulation and toxicokinetics of 42 polychlorinated biphenyl congeners in American kestrels (Falco sparverius). Environ Toxicol Chem 20: 2514-2522. doi: 10.1002/etc.5620201117
![]() |
[38] |
Sonne C, Gustavson K, Rigét FF, et al. (2009) Reproductive performance in East Greenland polar bears (Ursus maritimus) may be affected by organohalogen contaminants as shown by physiologically-based pharmacokinetic (PBPK) modelling. Chemosphere 77: 1558-1568. doi: 10.1016/j.chemosphere.2009.09.044
![]() |
[39] |
Lohitnavy M, Lu Y, Lohitnavy O, et al. (2008) A possible role of multidrug resistance-associated protein 2 (Mrp2) in hepatic excretion of PCB126, an environmental contaminant: PBPK/PD modeling. Toxicol Sci 104: 27-39. doi: 10.1093/toxsci/kfn026
![]() |
[40] |
Lee SK, Ou YC, Andersen ME, et al. (2007) A physiologically based pharmacokinetic model for lactational transfer of PCB 153 with or without PCB 126 in mice. Arch Toxicol 81: 101-111. doi: 10.1007/s00204-006-0130-0
![]() |
[41] |
Emond C, Charbonneau M, Krishnan K. (2005) Physiologically based modeling of the accumulation in plasma and tissue lipids of a mixture of PCB congeners in female Sprague-Dawley rats. J Toxicol Environ Health Part A 68: 1393-1412. doi: 10.1080/15287390590956551
![]() |
[42] |
Lee SK, Ou YC, Yang RSH. (2002) Comparison of pharmacokinetic interactions and physiologically based pharmacokinetic modeling of PCB 153 and PCB 126 in nonpregnant mice, lactating mice, and suckling pups. Toxicol Sci 65: 26-34. doi: 10.1093/toxsci/65.1.26
![]() |
[43] | Tuey D, Matthews B. (1977) Pharmacokinetics of 3,3',5,5'-tetrachlorobiphenyl in the male rat. Drug Metab Dispos 5: 444-450. |
[44] | Lutz R, Dedrick R, Matthews H, et al. (1977) A preliminary pharmacokinetic model for several chlorinated biphenyls in the rat. Drug Metab Dispos 5: 386-396. |
[45] | Nichols JW, Fitzsimmons PN, Whiteman FW, et al. (2004) A physiologically based toxicokinetic model for dietary uptake of hydrophobic organic compounds by fish I. Feeding studies with 2,2',5,5'-tetrachlorobiphenyl. Toxicol Sci 77: 206-218. |
[46] |
Karara A, Mcfarland V. (1992) A pharmacokinetic analysis of the uptake of polychlorinated biphenyls (PCBs) by golden shiners. Environ Toxicol Chem 11: 315-320. doi: 10.1002/etc.5620110305
![]() |
[47] |
Kulkarni MG, Karara AH. (1990) A pharmacokinetic model for the disposition of polychlorinated biphenyls ( PCBs ) in channel catfish. Aquat Toxicol 16: 141-150. doi: 10.1016/0166-445X(90)90083-2
![]() |
[48] | Pizza JC, Connor JMO. (1983) PCB dynamics in Hudson River striped bass II. Accumulation from dietary sources. Aquat Toxicol 3: 313-327. |
[49] |
Borchard RE, Hansen LG, Huber WG, et al. (1974) Pharmacokinetics of Aroclor 1254 components after intravenous administration to swine and sheep. Arch Environ Contam 2: 179-192. doi: 10.1007/BF01975471
![]() |
[50] |
Van Eijkeren JCH, Zeilmaker MJ, Kan CA, et al. (2006) A toxicokinetic model for the carry-over of dioxins and PCBs from feed and soil to eggs. Food Addit Contam 23: 509-517. doi: 10.1080/02652030500512045
![]() |
[51] |
Yu KN, Lam PKS, Cheung CCC, et al. (2002) Mathematical modeling of PCB bioaccumulation in Perna viridis. Mar Pollut Bull 45: 332-338. doi: 10.1016/S0025-326X(02)00099-1
![]() |
[52] | Verner M, Mcdougall R, Glynn A, et al. (2013) Is the relationship between prenatal exposure to PCB-153 and decreased birth weight attributable to pharmacokinetics? Environ Health Perspect 121: 1219-1225. |
[53] |
Ulaszewska MM, Ciffroy P, Tahraoui F, et al. (2012) Interpreting PCB levels in breast milk using a physiologically based pharmacokinetic model to reconstruct the dynamic exposure of Italian women. J Expo Sci Environ Epidemiol 22: 601-609. doi: 10.1038/jes.2012.36
![]() |
[54] |
Trnovec T, Dedík L, Jusko TA, et al. (2011) Assessment of exposure to PCB 153 from breast feeding and normal food intake in individual children using a system approach model. Chemosphere 85: 1687-1693. doi: 10.1016/j.chemosphere.2011.09.013
![]() |
[55] |
Redding LE, Sohn MD, McKone TE, et al. (2008) Population physiologically based pharmacokinetic modeling for the human lactational transfer of PCB-153 with consideration of worldwide human biomonitoring results. Environ Health Perspect 116: 1629-1634. doi: 10.1289/ehp.11519
![]() |
[56] |
Gascon M, Verner MA, Guxens M, et al. (2013) Evaluating the neurotoxic effects of lactational exposure to persistent organic pollutants (POPs) in Spanish children. Neurotoxicology 34: 9-15. doi: 10.1016/j.neuro.2012.10.006
![]() |
[57] |
Verner MA, Plusquellec P, Muckle G, et al. (2010) Alteration of infant attention and activity by polychlorinated biphenyls: unravelling critical windows of susceptibility using physiologically based pharmacokinetic modeling. Neurotoxicology 31: 424-431. doi: 10.1016/j.neuro.2010.05.011
![]() |
[58] |
Verner MA, Ayotte P, Muckle G, et al. (2009) A physiologically based pharmacokinetic model for the assessment of infant exposure to persistent organic pollutants in epidemiologic studies. Environ Health Perspect 117: 481-487. doi: 10.1289/ehp.0800047
![]() |
[59] |
Verner MA, Bachelet D, McDougall R, et al. (2011) A case study addressing the reliability of polychlorinated biphenyl levels measured at the time of breast cancer diagnosis in representing early-life exposure. Cancer Epidemiol Biomarkers Prev 20: 281-286. doi: 10.1158/1055-9965.EPI-10-0992
![]() |
[60] |
Verner MA, Charbonneau M, López-Carrillo L, et al. (2008) Physiologically based pharmacokinetic modeling of persistent organic pollutants for lifetime exposure assessment: a new tool in breast cancer epidemiologic studies. Environ Health Perspect 116: 886-892. doi: 10.1289/ehp.10917
![]() |
[61] |
Sonne C, Gustavson K, Rigét FF, et al. (2014) Physiologically based pharmacokinetic modeling of POPs in Greenlanders. Environ Int 64: 91-97. doi: 10.1016/j.envint.2013.12.006
![]() |
[62] |
Abass K, Huusko A, Nieminen P, et al. (2013) Estimation of health risk by using toxicokinetic modelling: a case study of polychlorinated biphenyl PCB153. J Hazard Mater 261: 1-10. doi: 10.1016/j.jhazmat.2013.07.011
![]() |
[63] |
Ayotte P, Dewailly E, Bruneau S, et al. (1995) Arctic air pollution and human health: what effects should be expected? Sci Total Environ 160-161: 529-537. doi: 10.1016/0048-9697(95)04387-G
![]() |
[64] |
Broding HC, Schettgen T, Göen T, et al. (2007) Development and verification of a toxicokinetic model of polychlorinated biphenyl elimination in persons working in a contaminated building. Chemosphere 68: 1427-1434. doi: 10.1016/j.chemosphere.2007.04.014
![]() |
[65] |
Arnot JA, Brown TN, Wania F. (2014) Estimating screening-level organic chemical half-lives in humans. Environ Sci Technol 48: 723-730. doi: 10.1021/es4029414
![]() |
[66] |
Tonnelier A, Coecke S, Zaldívar JM. (2012) Screening of chemicals for human bioaccumulative potential with a physiologically based toxicokinetic model. Arch Toxicol 86: 393-403. doi: 10.1007/s00204-011-0768-0
![]() |
[67] | Ritter R, Scheringer M, MacLeod M, et al. (2011) Intrinsic human elimination half-lives of polychlorinated biphenyls derived from the temporal evolution of cross-sectional biomonitoring data from the United Kingdom. Environ Health Perspect 119: 225-231. |
[68] |
Alcock RE, Sweetman AJ, Juan CY, et al. (2000) A generic model of human lifetime exposure to persistent organic contaminants: development and application to PCB-101. Environ Pollut 110: 253-265. doi: 10.1016/S0269-7491(99)00298-5
![]() |
[69] |
Aylward LL, Green E, Porta M, et al. (2014) Population variation in biomonitoring data for persistent organic pollutants (POPs): an examination of multiple population-based datasets for application to Australian pooled biomonitoring data. Environ Int 68: 127-138. doi: 10.1016/j.envint.2014.03.026
![]() |
[70] |
Becker K, Kaus S, Krause C, et al. (2002) German Environmental Survey 1998 (GerES III): Environmental pollutants in blood of the German population. Int J Hyg Environ Health 205: 297-308. doi: 10.1078/1438-4639-00155
![]() |
[71] |
Porta M, Gasull M, Puigdomènech E, et al. (2010) Distribution of blood concentrations of persistent organic pollutants in a representative sample of the population of Catalonia. Environ Int 36: 655-664. doi: 10.1016/j.envint.2010.04.013
![]() |
[72] |
Koppen G, Den Hond E, Nelen V, et al. (2009) Organochlorine and heavy metals in newborns: Results from the Flemish Environment and Health Survey (FLEHS 2002-2006). Environ Int 35: 1015-1022. doi: 10.1016/j.envint.2009.05.002
![]() |
[73] |
Schoeters G, Den Hond E, Colles A, et al. (2012) Concept of the Flemish human biomonitoring programme. Int J Hyg Environ Health 215: 102-108. doi: 10.1016/j.ijheh.2011.11.006
![]() |
[74] |
Schade G, Heinzow BU. (1998) Organochlorine pesticides and polychlorinated biphenyls in human milk of mothers living in northern Germany: Current extent of contamination , time trend from 1986 to 1997 and factors that influence the levels of contamination. Sci Total Environ 215: 31-39. doi: 10.1016/S0048-9697(98)00008-4
![]() |
[75] |
Ritter R, Scheringer M, MacLeod M, et al. (2009) A multi-individual pharmacokinetic model framework for interpreting time trends of persistent chemicals in human populations: application to a postban situation. Environ Health Perspect 117: 1280-1286. doi: 10.1289/ehp.0900648
![]() |
1. | Léon Masurel, Carlo Bianca, Annie Lemarchand, Space-velocity thermostatted kinetic theory model of tumor growth, 2021, 18, 1551-0018, 5525, 10.3934/mbe.2021279 | |
2. | Wenhai Shan, Pan Zheng, Global boundedness of the immune chemotaxis system with general kinetic functions, 2023, 30, 1021-9722, 10.1007/s00030-023-00840-4 | |
3. | Andrii Bomba, Serhii Baranovsky, Oksana Blavatska, Larysa Bachyshyna, Infectious disease model generalization based on diffuse perturbations under conditions of body's temperature reaction, 2022, 146, 00104825, 105561, 10.1016/j.compbiomed.2022.105561 | |
4. | Pan Zheng, Wenhai Shan, Global boundedness and stability analysis of the quasilinear immune chemotaxis system, 2023, 344, 00220396, 556, 10.1016/j.jde.2022.11.004 | |
5. | Carlo Bianca, A decade of thermostatted kinetic theory models for complex active matter living systems, 2024, 50, 15710645, 72, 10.1016/j.plrev.2024.06.015 | |
6. | Mengxin Chen, Spatiotemporal inhomogeneous pattern of a predator–prey model with delay and chemotaxis, 2023, 111, 0924-090X, 19527, 10.1007/s11071-023-08883-z | |
7. | Pan Zheng, Wenhai Shan, Guangyuan Liao, Stability Analysis of the Immune System Induced by Chemotaxis, 2023, 22, 1536-0040, 2527, 10.1137/22M1498310 | |
8. | Baranovsky S, Bomba A, Automated control of the immune response to the causants of infectious disease, 2023, 28, 27101673, 10, 10.15407/jai2023.03.010 | |
9. | Xiaoyan Gao, Liangying Miao, Global boundedness and asymptotic behavior of solutions to a three-dimensional immune chemotaxis system, 2024, 75, 0044-2275, 10.1007/s00033-024-02205-y | |
10. | Jianping Wang, Mingxin Wang, Global boundedness and finite time blow-up of solutions for a quasilinear chemotaxis-May–Nowak model, 2025, 35, 0218-2025, 1933, 10.1142/S0218202525500320 |
Inlet boundaries | 0.4 s | 0.7 s | ||
AF | SR | AF | SR | |
Left Superior Pulmonary Vein | 26.86% | 27.88% | - | 24.49% |
Left Inferior Pulmonary Vein | 30.27% | 30.56% | - | 29.02% |
Right Superior Pulmonary Vein | 26.14% | 25.35% | - | 27.20% |
Right Inferior Pulmonary Vein | 16.74% | 16.22% | - | 19.29% |
AF: atrial fibrillation; SR: sinus rhythm. |