|Computer Modeling in Engineering & Sciences|
Reducing the Range of Cancer Risk on BI-RADS 4 Subcategories via Mathematical Modelling
1Near East University, Faculty of Arts and Sciences, Nicosia, 99138, North Cyprus
2Near East University, Mathematics Research Center, Nicosia, 99138, North Cyprus
3Near East University, Faculty of Medicine, Nicosia, 99138, North Cyprus
*Corresponding Author: Evren Hincal. Email: email@example.com
Received: 14 October 2021; Accepted: 07 February 2022
Abstract: Breast Imaging Reporting and Data System, also known as BI-RADS is a universal system used by radiologists and doctors. It constructs a comprehensive language for the diagnosis of breast cancer. BI-RADS 4 category has a wide range of cancer risk since it is divided into 3 categories. Mathematical models play an important role in the diagnosis and treatment of cancer. In this study, data of 42 BI-RADS 4 patients taken from the Center for Breast Health, Near East University Hospital is utilized. Regarding the analysis, a mathematical model is constructed by dividing the population into 4 compartments. Sensitivity analysis is applied to the parameters with the desired outcome of a reduced range of cancer risk. Numerical simulations of the parameters are demonstrated. The results of the model have revealed that an increase in the lactation rate and early menopause have a negative correlation with the chance of being diagnosed with BI-RADS 4 whereas a positive correlation increase in age, the palpable mass, and family history is distinctive. Furthermore, the negative effects of smoking and late menopause on BI-RADS 4C diagnosis are vehemently outlined. Consequently, the model showed that the percentages of parameters play an important role in the diagnosis of BI-RADS 4 subcategories. All things considered, with the assistance of the most effective parameters, the range of cancer risks in BI-RADS 4 subcategories will decrease.
Keywords: Mathematical model; breast cancer; BI-RADS 4; sensitivity analysis; reducing cancer risk; suspicious malignancy; ordinary differential equations
Cancer is a serious disease that occurs as a result of an uncontrolled proliferation and growth of cells in any organ or tissue of the body. There is a requirement of building new cells in the event of cell death or damage. When the order of this formation process fails, damaged or abnormal cells can grow. The growth or proliferation of these kind of cells may result in tumors (which are sometimes called as lumps of tissues) [1,2]. Tumors can be benign or malignant. Malignant tumors are known as cancerous tumors and they have the capability to proliferate, invade the surrounding structures and make metastatic-colonies. Benign tumors do not possess these kind of properties .
According to the World Health Organization (WHO), cancer is a persuasive cause of death, vehemently prevalent in the world. In 2020, approximately 10 million people passed away because of cancer .
Cancer may begin in almost any part of the body and it is named dependent on the tissue in which it occurs. There are more than 200 types of identified cancers. The most common and mortal types of cancer include lung, stomach, liver, colon, and breast cancer. Female breast cancer is the most common cancer type with an estimated 2,3 million new cases (11.7%) around the globe . It has the highest mortality rate in contrast to the other types of cancer for women [6,7] and the statistics of W.H.O. revealed that 685000 people out of 2,3 million diagnoses passed away in 2020 . In fact, breast cancer is the primary cause of mortality; especially in post-menopausal women accounting 23% of all cancer related deaths .
Breast cancer generally originates from the epithelial cells of the breast ducts and glands. The most common clinical presentation is a newly formed breast lump. However, bloody nipple discharge, skin dimpling or ulceration, abnormal change in the breast size, redness or skin edema may be a sign or symptom of breast cancer. In spite of this fact, sometimes there may be a radiological abnormality during the screening but no clinical finding [9,10].
There are risk factors associated with breast cancer. To illustrate, these include age, female gender, genetic mutations, family history with ovarian or breast cancer, obesity, personal history with other types of cancer, etc. Dense breast on mammography is also a risk factor and, in this case, the diagnostic capability of mammography is relatively insufficient [9,10].
Breastfeeding, late first menstrual period and premature menopause are some of the factors that may reduce the risk of breast cancer [11,12]. As in other diseases, early diagnosis is the inevitable key factor for the treatment of breast cancer. After the age of 40, a screening mammography in every 1–2 year is recommended while after 50 s, screening mammography should be performed once a year.
In asymptomatic women, a screening mammography should be performed in order to increase the chance of successful treatment. On the other hand, in symptomatic women, the frequency of diagnostic mammography performed is higher due to the abnormality in the breast. Diagnostic mammography enables doctors to determine the exact size, location, and the shape of the lump or other pathological finding, exactly .
American College of Radiology (ACR) developed a categorization to reduce and group the variation of radiologists’ findings and descriptions for the diagnosis. Breast Imaging Reporting and Data System is given to the name of this categorization, which is also referred as BI-RADS, with initials of the name. BI-RADS system has 7 categories, starting from 0 till 6. The categorization of the BI-RADS is determined according to the obtained data from the patient’s ultrasound, mammography, etc. In the case of Category 0, further examinations are necessary in order to identify the obtained lump . BI-RADS 1 is the category defined as negative findings and Category 2 is read as benign structures with 100%. In other words, in Category 2, some instances of evaluated harmless lumps are demonstrated. Category 3 is considered as benign with 98%, another control is generally suggested after 6 months. Biopsy is recommended for Category 4 since the lesions are suspicious with an enormous range, 2% to 95%. Category 5 has the probability of 95% malignancy and Category 6 is the class of proven malignancy, i.e., the patient is 100% breast cancer .
BI-RADS 4 category has a wide range of cancer probability which may portend contradictions for the doctors while diagnosis. Due to this reason, it is divided into 3 subcategories as BI-RADS 4A, BI-RADS 4B, and BI-RADS 4C. Category 4A has a relatively lower (2%–10%) malignancy rate; malignancy of Category 4B is moderate with 10%–50% and Category 4C has high suspicion of malignancy with 50%–95%. The categorization of findings may arise controversies between radiologists since their interpretations can be different .
Mathematical modelling has been widely used in the diagnosis or treatment of diseases for many years. Preventing the spread or occurrence of the tumor/cancer cells via control structures is one of the leading aims of mathematical modelling in health sciences. After identifying the problem, the population that is going to be examined should be separated into sensible number of compartments. Subsequently, with the help of parameters and variables, relationship between the compartments is constructed by means of differential equations . The kind of differential equations (ODE, PDE, FDE, etc.) may differ and it depends on the type or structure of the disease and field of mathematicians.
For the effect of any disease, obtaining the basic reproduction number, denoted by , with mathematical methods is a remarkable way. One of the most popular methods for obtaining of diseases is the Next Generation Matrix Method. With this method, formulation of is evaluated by using matrices that are generated from the model . The results of different kinds of mathematical models in health sciences provide doctors and patients new perspectives and ideas for the control, diagnosis and treatment of diseases. Thus, introducing new parameters or variables for the models provides further but more importantly, proved results about the diseases.
In , a mathematical model is constructed for tumor growth while  modelled growth of tumor with single irradiation. Fractional mathematical modelling is used in  for the breast cancer and in  for the system of cardiac conduction. Jajarmi et al.  represents a fractional system about the dynamics of immunogenic tumors. Treatments like drug intervention, chemotherapy and immunotherapy are modelled in the studies [23–25].
In this paper, BI-RADS subcategories in breast radiological analyses is studied. The study includes 42 patients of Near East University Hospital, Center for Breast Health. Patients are diagnosed with BI-RADS 4 and divided into subcategories according to the findings. The principal focus of the study is to determine the effect of introduced parameters on the BI-RADS 4 subcategories and to narrow the range of cancer risk in the subcategories of BI-RADS 4.
The paper is organized as follows: First of all, the information about the used data and methods are demonstrated in Section 2. Then, the constructed model and its analyses are examined in Section 3. Afterwards, numerical simulations of the constructed model are presented in Section 4. Lastly, in Section 5, the conclusions and discussions that are obtained are evaluated.
In this section, data and methods that are used during the study are explained precisely.
In this study, the data of 42 patients that are obtained in the Center for Breast Health, Near East University Hospital are used. The parameters are determined according to the obtained data.
With the help of collected data, the population is divided into 4 compartments; susceptible individuals , BI-RADS 4A diagnosed individuals , BI-RADS 4B diagnosed individuals , and BI-RADS 4C diagnosed individuals . For obtaining the change in compartments at time , a mathematical model is proposed with Ordinary Differential Equations (ODEs). In other words, ODEs are constructed according to the relationship between compartments, by means of parameters. Sensitivity analysis is applied to the parameters, separately, to see their effects on compartments. The idea of sensitivity analysis is that it shows how a small change in one or more parameters can affect the disease to spread or die out . A program is created via MatLab in order to support the results of sensitivity analysis.
In this model, values are obtained by using the next generation matrix method and it represents number of expected new BI-RADS 4 diagnoses. Matrices are constructed by using the inputs and outputs of equations. As in infectious diseases, values should be less than 1 to say that there will be no new diagnosis for BI-RADS 4.
In this section, the constructed model is given and the existence of the solution is proved. Afterwards, its analysis is given.
The population, denoted by , is divided into 4 compartments; susceptible individuals, denoted by , BI-RADS 4A diagnosed individuals, denoted by , BI-RADS 4B diagnosed individuals, denoted by , and BI-RADS 4C diagnosed individuals, denoted by . That is, , at time . Model is constructed as follows:
The description of variables and parameters are given in Tables 1 and 2. In the construction of the model, it is assumed that transition between the subcategories process step by step. That is, a member of BI-RADS 4A cannot be a member of BI-RADS 4C directly; it should be a member of BI-RADS 4B. Also, it is assumed that increase in smoking causes any patient to be diagnosed as BI-RADS 4B and BI-RADS 4C because of the idea that smoking may increase the breast cancer risk in long term . In the study, an early menopause is considered as going through menopause before the age of 40 while a late menopause means going through menopause after the age of 55.
Theorem 3.1. Let be the solution of the constructed system with the initial conditions and . Then, the following set:
is positive, invariant, and all of the solutions in stay in with respect to the constructed system .
Proof. Addition of all of the terms that are on the right side of the system gives
From the equality, it is clear that . Applying integration to the both sides with respect to yields
for some arbitrary constant . With the use of Rota and Birkhoff for the above differential inequality, it can be obtained that as tends to infinity, , . Consequently, the solutions of the given system enter the region Thus, it is certain that the model is feasible by means of biology and it is enough to consider the dynamics on the model in .
The Disease Free Equilibrium Point, denoted by , is the point where the disease dies out. In this case, is the point where BI-RADS 4 category does not exist in the population and hence no one is diagnosed as BI-RADS 4. In order to reach , value(s) of the disease should be less than 1 . of this model is unique and it obtained as
It is obvious that attracts the region so that
In this model, 3 different are obtained since there are 3 different categories of disease. For the use of Next Generation Matrix Method, the following matrices are constructed. denotes the matrix of new diagnoses of BI-RADS 4 subcategories while the matrix is written with the rest of the system.
Dominant eigenvalues of the matrix multiplication are the formulas of BI-RADS 4 subcategories.
are the values of BI-RADS 4A, BI-RADS 4B, and BI-RADS 4C, respectively.
Theorem 3.2.1. Disease Free Equilibrium, , is globally asymptotically stable whenever , , and .
Proof. Consider the Lyapunov function
where . It is clear that the constructed function is always positive and it equals to at the point . So, it is enough to show that .
It is obvious that . For the rest of the terms
if and hence
Similarly, and if and , respectively. Thus, is globally asymptotically stable if , , and .
Table 3 shows the importance of high lactation rate and early menopause for the BI-RADS 4A categorization. On the other hand, increase in the smoking rate causes the patients of BI-RADS 4A to be diagnosed as BI-RADS 4B or even higher category. The most effective parameter is the family history. Increase in the members of family history causes increase in .
According to the Table 4, likewise BI-RADS 4A, increase in the rate of family history causes an increase in , i.e., in BI-RADS 4B. However, increase in the smoking rate causes patients to be diagnosed as BI-RADS 4C. Also, higher lactation rate and earlier menopause causes a decrease in BI-RADS 4B category.
In Table 5, as in BI-RADS 4A and BI-RADS 4B, increase in the number of members diagnosed as breast cancer before in family causes a serious increase in . For BI-RADS 4C as well, higher lactation rate and earlier menopause means lower BI-RADS 4C diagnoses.
In this section, numerical simulations for the distribution of BI-RADS 4 subcategories and sensitivity analysis are shown.
Fig. 1 shows that in time, BI-RADS 4A category is expected to extinct and hence BI-RADS 4 category is going to be separated into 2 categories. This may lead a decrease in the range of cancer probability for BI-RADS 4.
Here the idea of sensitivity analysis is to increase the values of parameters (from the obtained data) by 10% and see the trend of BI-RADS 4A, BI-RADS 4B, and BI-RADS 4C at time .
Figs. 2–4 show the effect of parameter on the subcategories of BI-RADS 4. According to the figures, age causes an increase in each of the BI-RADS 4 subcategories. In other words, as people get older, the chance of being diagnosed as BI-RADS 4 increases.
Figs. 5–7 present the effect of parameter on the subcategories of BI-RADS 4. As expected, increase in the palpable mass raises the probability of BI-RADS 4 diagnosis for individuals. It is clear from the figures that increase in the palpable mass mostly increases the risk of BI-RADS 4C diagnosis.
In Figs. 8–10, the effect and hence the importance of parameters and on the subcategories BI-RADS 4A, BI-RADS 4B, and BIRADS 4C are presented, respectively. The figures showed that increase in the lactation rate might lead a meaningful decrease in BI-RADS 4 subcategories. Especially for the categories BI-RADS 4B and BI-RADS 4C, there is a significant positive effect in the case of active lactation.
In the Figs. 11–13, the effect of parameters and on the subcategories BI-RADS 4A, BI-RADS 4B, and BIRADS 4C are given, respectively. An early menopause causes a decrease in the subcategories according to the figures. For all of the subcategories, a meaningful change exists in this case.
Effect of the parameters and , i.e., the family history on the subcategories BI-RADS 4A, BI-RADS 4B, and BIRADS 4C are presented in the Figs. 14–16, respectively. Family history causes an enormous increase in the diagnosis of BI-RADS 4 subcategories as can be seen from the figures.
Fig. 17 is obtained to see and show the effect of smoking on the subcategory BI-RADS 4C. Increase in the smoking rate leads an increase in the BI-RADS 4C subcategory, which is expected to happen.
The study is constructed with the aim of reducing the huge range of cancer risk for BI-RADS 4 subcategories. In the analysis of the model, disease-free equilibrium point is found with the proof of its existence. Furthermore, it has shown that it is globally asymptotically stable under some conditions. Consequently, it is possible to have a population without BI-RADS 4 category at .
Simulation of the BI-RADS 4 subcategories, given in Fig. 1, shows that Category of BI-RADS 4 will be divided into 2 categories in the coming years. Hence, it is expected to have more significant properties for the categories so that the diagnosis of patients can be divided into 2 categories which may lead convenience for the radiologists and doctors.
Sensitivity analysis applied to the parameters of the model in Subsection 3 aims to identify the effects of parameters on the BI-RADS 4 subcategories. According to the Figs. 2–4, age has a negative effect on all of the BI-RADS 4 subcategories since it causes an increase in patients. That is, increase in age increases the chance of diagnosing as BI-RADS 4 for people. As it is expected, increase in the palpable mass leads to an increase for each BI-RADS 4 subcategories. The results and the pattern of diagnosis with the increase are obtained and shown in Figs. 5–7 for each BI-RADS 4 subcategory, separately.
An increase in the lactation rate has a positive effect on all BI-RADS 4 subcategories, because it causes an expressive decrease in diagnosis. The results of an increase in the percentage of lactation rate are given in Figs. 8–10 for BI-RADS 4 subcategories, respectively. Especially in the subcategories BI-RADS 4B and BI-RADS 4C, the lactation rate can be one of the most important parameters. This can be interpreted as giving birth decreases the risk of being breast cancer, substantially.
Early menopause is also an important parameter so that women have fewer chance to be diagnosed as BI-RADS 4 (especially BI-RADS 4B and BI-RADS 4C) if they go through menopause earlier. The graph of this result is given in Figs. 11–13.
The numerical results of the model showed that family history is another important parameter for the BI-RADS 4. An increase in the members with cancer history develops an inclination towards BI-RADS 4 diagnosis. This result is shown in Figs. 14–16 for each subcategory, individually. Thus, it can be interpreted as; regular health controls are important and necessary for a person with family history. That is, the breast cancer screenings should be more frequent and maybe detailed or alternative screening modalities like breast magnetic resonance imaging should be used.
In Fig. 17, the effect of increase in smoking in BI-RADS 4C is given. According to the results, smoking plays more important role than some parameters like age or palpable mass in BI-RADS 4C. It causes an increase in BI-RADS 4C which may result in the possession of a tendency towards breast cancer. The option of being a non-smoker is the mere precaution that can be taken against cancer.
The negative effects of late menopause on BI-RADS 4C is shown in Fig. 18. Going through menopause in old ages increases the risk of BI-RADS 4C diagnosis according to the obtained results.
In conclusion, the results showed that being non-smoker, increasing the lactation rate and regular controls (especially for the people with family history) are some of the significant parameters that can be controlled by people themselves. With this way and with the use of screening methods, people can reduce the risk of being cancer since they enable early diagnosis and so treatment.
In addition, the model revealed that BI-RADS 4 subcategories will be divided into 2 compartments in time. This may lead to a meaningful decrease in the range of cancer risk in future with the control of studied parameters. For further studies, a non-standard finite difference scheme can be introduced into systems as in  and .
Funding Statement: The authors received no specific funding for this study.
Conflicts of Interest: The authors declare that they have no conflicts of interest to report regarding the present study.
|This work is licensed under a Creative Commons Attribution 4.0 International License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.|