Submit or Track your Manuscript LOG-IN
Latest Blogs: https://researcherslinks.com/en/kahoot-login/ https://researcherslinks.com/en/blooket-login/ https://researcherslinks.com/en/comcast-login/ https://researcherslinks.com/en/gimkit-login/

Advances in Animal and Veterinary Sciences

AAVS_8_4_354-359

 

 

Research Article

 

Classification and Regression Tree (CRT) Analysis to Predict Body Weight of Potchefstroom Koekoek Laying Hens

 

Thobela Louis Tyasi1*, Kgotlelelo Maaposo Makgowo1, Kwena Mokoena1, Lebo Trudy Rashijane1, Madumetja Cyril Mathapo1, Lebogang William Danguru1, Kagisho Madikadike Molabe1, Paul Mogowe Bopape1, Nhlakanipho Divine Mathye1, Dannis Maluleke1, Busisiwe Gunya1, Masibonge Gxasheka2

1School of Agricultural & Environmental Sciences, Department of Agricultural Economics and Animal Production, University of Limpopo, Private Bag X1106, Sovenga 0727, Limpopo, South Africa; 2School of Agricultural & Environmental Sciences, Department of Plant Production, Soil Science & Agricultural Engineering, University of Limpopo, Private Bag X1106, Sovenga 0727, Limpopo, South Africa.

 

Abstract | Classification and regression tree analysis is a powerful statistical technique which helps to determine the most important variables in a particular dataset and helps to create a model. The study was conducted to identify linear body measurement traits (beak length, body length, keel length, chest circumference, toe length, body girth, shank length, back length, shank circumference and wing length) which could be employed in developing an effective prediction equation for body weight of Potchefstroom Koekoek laying hens. Eighty Potchefstroom Koekoek laying hens at twenty two weeks old were used. Pearson’s correlation together with classification and regression tree (CRT) methods were used for analysis. Descriptive statistics indicated that mean of body weight was 1.50 kg. Correlation findings revealed that body weight was positively significantly correlated (P < 0.05) with beak length (r = 0.23) and toe length (r = 0.21), respectively. CRT results demonstrated that beak length, wing length and back length play an important role in the body weight of Potchefstroom Koekoek laying hen chickens. This study suggests that body weight of laying hens could be estimated by some linear body measurement traits. The models established in the current study might be employed by chicken farmers when making selection during breeding to improve body weight. However, further studies need to be done to validate the use of classification and regression tree analysis in prediction of body weight from linear body measurement traits of chickens.

 

Keywords | Back length, Beak length, Body measurements, Correlation, Wing length

 

Received | November 21, 2019; Accepted | February 17, 2020; Published | March 20, 2020

*Correspondence | Thobela Louis Tyasi, School of Agricultural & Environmental Sciences, Department of Agricultural Economics and Animal Production, University of Limpopo, Private Bag X1106, Sovenga 0727, Limpopo, South Africa; Email: louis.tyasi@ul.ac.za

Citation | Tyasi TL, Makgowo KM, Mokoena K, Rashijane LT, Mathapo MC, Danguru LW, Molabe KM, Bopape PM, Mathye ND, Maluleke D, Gunya B, Gxasheka M (2020). Classification and regression tree (crt) analysis to predict body weight of potchefstroom koekoek laying hens. Adv. Anim. Vet. Sci. 8(4): 354-359.

DOI | http://dx.doi.org/10.17582/journal.aavs/2020/8.4.354.359

ISSN (Online) | 2307-8316; ISSN (Print) | 2309-3331

Copyright © 2020 Tyasi et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

 

INTRODUCTION

 

Potchefstroom Koekoek chicken genotype was developed from crossbreeding three chicken genotypes namely White Leghorn, Black Australorp and Bared Plymouth Rock (Tyasi et al., 2019; Mutibvu et al., 2019). This chicken genotype has gender-specific feather colour and patterns, which make it easy for sexing of day-old chicks (Heit, 2017). According to Dessie and Gatachew (2016), this chicken genotype is well acclimatized to the tropical regions which enable it to thrive in many of hot environment of the African countries, and is suitable for free-range production system and can also thrive and perform highly in small-scale rearing system.

 

Classification and regression tree (CRT) is a data mining algorithm commonly used where the dependent variable is categorical which help to identify the class within which a dependent variable fall into, and also used where the dependent variable is continuous which help to predict the value (Mohamad et al., 2012). According to Aytekin et al. (2018), CRT is a statistical technique which is applicable in different forms of data including continuous, ordinal and nominal variables. Based on our knowledge, there is limited information about the estimation of body weight from linear body measurement traits using regression tree data mining algorithms. Body weight and body measurement traits are important characteristics to chicken breeders and processors (Ojo et al., 2014). Body weight plays a vital function in discovering some economic traits in farm livestock (Ojedapo, 2013). Linear body measurement traits have been indicated as the good predictor of body weight in different chicken genotypes including Chinese indigenous Dagu chicken genotype (Tyasi et al., 2017), Uganda indigenous chickens (Semakula et al., 2011), Senegal indigenous chickens (Gueye et al., 1998), two commercial meat-type (Ross and Anak Titan) chicken genotype (Ajayi et al., 2008) and Cobb broiler chicken genotype (Latshaw and Bishop 2001). Ukwu et al. (2014) indicated that more studies need to be conducted for prediction of body weight of some chicken genotypes using linear body measurement traits.

 

The objectives of the present study were: 1) to determine the relationship between body weight (BW) and linear body measurement traits which are beak length (BKL), body length (BDL), keel length (KL), chest circumference (CC), toe length (TL), body girth (BG), shank length (SL), back length (BCL), shank circumference (SC) and wing length (WL) of Potchefstroom koekoek using correlation analysis, 2) to develop a model for predicting body weight from linear body measurement traits using classification and regression tree data mining algorithm. This study will help Potchefstroom koekoek chicken farmers in selecting the best linear body measurement trait during breeding to improve body weight.

 

MATERIALS AND METHODS

 

Study Area

The study was carried out at the University of Limpopo experimental farm, South Africa. The location, ambient temperatures, latitudes, longitude and rainfall of the study area are the same as described by Alabi et al. (2012).

 

Experimental Birds and Management

Potchefstroom Koekoek chicken genotype was used for the current study. Chicks for this study were hatched at the university farm and reared intensively in a closed space on deep litter since day old until they are pullets. At the age of twenty two weeks, the chickens were separated per sex and eighty (80) laying hens were randomly selected from the population for the experiment. The layers were fed with 16.10% crude protein diet 11.97 MJkg/DM. The diet comprised of ingredients like: maize (64%), maize gluten meal (11.67%), soya Hipro (4.37%), fish meal (5%), full fat soya (4.91%), Di sodium phosphate (1.33%), L-lysine (0.20%), CaCO3 (8.17), DL-methionine (0.20%) and vitamin trace element premix (0.15%) was performed following the national research council (NRC) feeding standard. Water and feed were given to the chickens without restriction.

 

Data Collection

Body weight and linear body measurements were collected as described by Ige et al. (2016) and Tyasi et al. (2018). Briefly, body weight was measured with a weighing scale in kilograms (kg). The following linear body measurement traits were measured once per bird using a measuring tape and ruler in centimetres (cm): beak length (BKL), body length (BDL), keel length (KL), chest circumference (CC), toe length (TL), body girth (BG), shank length (SL), back length (BCL), shank circumference (SC) and wing length (WL). Linear body measurement traits were measured as follows: WL was measured as the length from the humorous-coracoid junction to the distal tip of the phalange digits; KL was recorded as the length between the cranial and the caudal terminals of the keel bone; SL was recorded as the length of tarso-metatarsus from the hock joint to the metatarsal pad; CC was recorded as the circumference of the chest; BKL was measured as the length of the beak; BDL was measured as the distance between the base of the neck and the cloaca; TL was measured as the length of the third toe; BG was taken when a measuring tape is looped around the region of the breast under the wing. BCL was measured from the nadir of the curve to the base of the tail, SC was taken as the width of the shank. All the measurements were collected with the same person to avoid individual variations on measuring.

 

Classification and Regression Tree (Crt) Method

This data mining algorithm is a binary decision tree. This method was proposed by Breiman et al. (1984) and it is very simple and easy to visualize. In the present study, CRT was applied to predict the dependent variable (body weight) from independent variables (beak length, body length, keel length, chest circumference, toe length, body girth, shank length, back length, shank circumference and wing length. The CRT was conducted as described by Yakubu (2012); Oguntunji (2017) and Eyduran et al. (2019). Briefly, 10 fold cross-validation an error estimation method known to be the most acceptable method of estimation was employed to provide estimates of the future prediction error for each node and explained the variation observed in the target variable was estimated as follows:

 

S2x = (1 – S2e) x 100

S2e = risk value/S2y

Where:

S2x = explained variation, S2e = unexplained variation and S2y = variance of the root node (standard deviation of the root node)2.

 

Statistical Analysis

Data was analysed by Statistical Analysis System (SAS, 2009) for descriptive statistics and R-software was used for computing of correlation between variables and constructing heat map of correlation coefficients. Statistical Package for Social Sciences (IBM SPSS.25) software was used for computing of classification and regression tree data mining algorithm.

 

RESULTS

 

Descriptive Statistics

Boxplot (Figure 1) is presented the summary of body weight distribution in Potchefstroom Koekoek laying hens. The summary showed the minimum, 25th percentile (first quartile), median, 75th percentile (third quartile) and maximum values of body weight. Potchefstroom Koekoek chicken genotype boxplot indicated that the values of minimum (>1.0kg), first quartile, median (>1.40kg), third quartile (1.60kg) and maximum (>1.60kg).

 

 

Table 1 shows the summary of body weight and linear body measurement traits including beak length, body length, keel length, chest circumference, toe length, body girth, shank length, back length, shank circumference and wing length. Body weight of Potchefstroom Koekoek was noted as 1.50 ± 0.02 kg while for wing length, back length beak length, shank length, shank circumference, chest circumference, keel length, body girth, toe length and body length were noted as 4.70 ± 0.71, 0.49 ± 0.01, 0.22 ± 0.01, 8.17 ± 0.16, 4.69 ± 0.09, 26.56 ± 0.55, 11.22 ± 0.17, 36.83 ± 0.56, 6.28 ± 0.11 and 33.15 ± 0.38 cm, respectively.

 

Table 1: Descriptive statistics of body weight and linear body measurement traits of Potchefstroom Koekoek hens.

 

Trait Mean ± SE CV Minimum Maximum
BW (kg) 1.50 ± 0.02 13.7 1.10 1.90
WL (cm) 4.70 ± 0.71 134 0.93 22.5
BCL (cm) 0.49 ± 0.01 25.0 0.30 0.79
BKL (cm) 0.22 ± 0.01 29.2 0.13 0.34
SL (cm) 8.17 ± 0.16 17.0 4.20 9.90
SC (cm) 4.69 ± 0.09 17.4 3.90 9.10
CC (cm) 26.6 ± 0.55 18.7 18.6 34.7
KL (cm) 11.2 ± 0.17 13.2 8.29 14.4
BG (cm) 36.8 ± 0.56 13.5 30.1 47.0
TL (cm) 6.28 ± 0.11 15.3 4.74 8.30
BDL (cm) 33.2 ± 0.38 10.4 27.2 38.6

 

Standard error of means: SE, Coefficient of variation: CV, Body weight: BW, Beak length: BKL, Body length: BDL, Keel length: KL, Chest circumference: CC, Toe length: TL, Body girth: BG, Shank length: SL, Back length: BCL, Shank circumference: SC and Wing length: WL.

 

Phenotypic Correlation

Phenotypic correlation between body weight and linear body measurement traits are presented in Figure 2. Pearson correlation results revealed that only two linear body measurement traits (toe length, r = 0.21* and beak length, r = 0.23*) which were positive significant correlation with body weight. Our results showed wing length had a negative highly significant correlation with back length (r = -0.56**), chest circumference (r = -0.29**) and keel length (r = -0.27**), while beak length had a positive highly significant correlation with keel length (r = 0.56**), toe length (r = 0.42**) and body length (r = 0.39**). Shank length had a positive highly significant correlation with chest circumference (r = 0.33**) but negative highly significant correlation with shank circumference (r = -0.33**), Chest circumference had a positive highly significant correlation with keel length (r = 0.47**), body girth (r = 0.34**), toe length (r = 0.46**), and body length (r = 0.56**). Keel length had a positive highly significant correlation with body girth (r = 0.28**), toe length (r = 0.46**), and body length (r = 0.59**). Body girth had a positive highly significant correlation with toe length (r = 0.55**) and body length (r = 0.41**). Toe length had a positive highly significant correlation with body length (r = 0.38**). Correlation results also indicated that body weight had positive non-significant correlation with wing length (r = 0.16ns), body girth (r = 0.01), chest circumference (r = 0.14), keel length (r = 0.15) and body length (r = 0.17) nevertheless, negative non-significant correlation with shank circumference (r = -0.03), back length (r = -0.08) and shank length (r = -0.05), respectively.


 

Classification and Regression Tree Method

The classification and regression tree (CRT) model in Potchefstroom Koekoek chicken genotype (Figure 3) with the body weight as a dependent variable and linear body measurement traits as the independent variables. This model contained a total of eight (8) nodes on which five of them were terminal nodes (node 1, 5, 6, 7 and 8), respectively. Node 0 is the root node containing the descriptive statistics of body weight (mean = 1.498 kg, standard deviation = 0.204 and n = 80). Node 0 on the basis of wing length was divided into node 1 (≤ 1.145 cm) and node 2 (˃ 1.145 cm), respectively. Node 2 was divided on the basis of beak length into node 3 ((≤ 0.200) and node 4 (˃ 0.200 cm). Node 3 on the basis of wing length was divided into node 5 (≤ 2.365 cm) and node 6 (˃2.365 cm). The last node to be divided in this model was node 4 which was divided on the basis of back length into node 7 (≤0.565 cm) and node 8 (˃0.565 cm), respectively. In all the terminal nodes, node 8 appeared to be the best node as it was recorded the highest predicted mean (1.830 kg) than node 1 (1.325 kg), node 5 (1.375 kg), node 6 (1.544 kg), and node 7 (1.527 kg), respectively. This model showed that node 6 had the lowest variance [(0.030)2 = 0.0009] and the variance of the root node or dependent variable (body weight) was S2y = (0.204)2 = 0.042. The unexplained variation in the body weight was S2e = risk value ÷ S2y = 0.018 ÷ 0.042 = 0.429 and the variation in the model was explained as S2y = 1 – S2e = 1- 0.042 = 0.57, respectively.

 

 

DISCUSSION

 

Body weight measurement is vital for any selection and breeding programme, feeding, vaccination and drug dosage in animal production (Olawumi and Farinnako, 2017). We firstly explored the distribution of Potchefstroom Koekoek laying hen’s body weight by using boxplot. The boxplot summary suggests that our data was well distributed without outliers. The data was summarized for all the studied traits and results indicated body weight of Potchefstroom Koekoek laying hens had 1.50kg. Our descriptive statistics were lower than of Tyasi et al. (2017) in Chinese Dagu male indigenous chickens, and Latshaw and Bishop (2001) in Cobb broiler chicken genotype, but higher than of Yakubu and Salako (2009) in Nigerian indigenous chickens. The variations might be due to the breed differences. The present study was conducted to investigate the association between body weight and linear body measurement traits of Potchefstroom Koekoek laying hens. The results demonstrate that by increasing toe length and beak length in chicken laying hen the body weight also increase. Therefore, toe length and beak length might be the selection criteria during breeding to improve body weight in Potchefstroom Koekoek chicken laying hens. These findings suggest that body weight and linear body measurement traits of Potchefstroom Koekoek laying hens need to be improved genetically. These results demonstrated that body weight does not depend on shank circumference, back length and shank length. Our findings are in accordance with findings of Tyasi et al. (2019). Other studies (Ajayi et al., 2008; Semakula et al., 2011) also found similar results between body weight and linear body measurement traits of chicken genotype. However, Egena et al. (2014) showed dissimilar correlation results between body weight and linear body measurement traits of Nigerian indigenous chicken genotype. The variation might be due to different chicken genotypes. Correlation results suggest that there is a relationship between body measurement traits of Potchefstroom Koekoek laying hens. However, the findings are not demonstrating which traits might be used to estimate the body weight. Hence, regression analysis was performed to predict the equations for the estimation of body weight using linear body measurement traits which were positively significant correlated with body weight. Our findings suggest that wing length, beak length and back length explained 57% of variation in body weight of Potchefstroom Koekoek chicken. The model developed from the current study showed that wing length had the highest significant role in body weight followed by beak length, respectively. It was also showed that back length also played a remarkable role in body weight of Potchefstroom Koekoek laying hen chicken. Our findings are supported by Dzungwe et al. (2018) who indicated that wing length is one of the linear body measurement traits which had a significant role in the body weight of chickens. Similar results found by Ojo et al. (2014) who also indicated that wing length plays an important in body weight of Japanese quail. However, Tyasi et al. (2017) reported that shank circumference and shank length played a vital role in the body weight of Chinese Dagu male chicken genotype. The variation with the current study might be due to sex differences.

 

In conclusion, the findings suggest that there is a relationship between body measurements. However, only beak length and toe length were positive significantly correlated with body weight. Classification and regression tree results suggest that beak length, wing length and back length had a role on body weight of Potchefstroom Koekoek laying hen chickens. Our findings might help chicken farmers, researchers and extension officers to recognize the association between body weight and linear body measurement traits in Potchefstroom Koekoek laying hen chickens. This study might also help farmers during selection to improve body weight in breeding. Further studies are needed to investigate the relationship between body weight and linear body measurement traits in larger sample of the same or different chicken breeds.

 

ACKNOWLEDGEMENTS

 

Authors wish to express their appreciation to the Department of Agricultural Economics and Animal Production, the University of Limpopo for financial support and to the University of Limpopo Experimental Farmworkers for their support during data collection.

 

CONFLICTS OF INTEREST

 

Authors for this work declare no conflict of interest.

 

AUTHOR CONTRIBUTIONS

 

TLT designed the experiment, analysed the data and wrote the manuscript. KMM, KM, TLR, MCM, LWD, KMM, PMB, NDM and DM performed fieldwork and wrote the manuscript. KM and KMM reviewed the manuscript. TLT, BG and MG read, edited and approved the final manuscript.

 

REFERENCES

 

  • Ajayi FO, Ejiofor O, Ironkwe MO (2008). Estimation of Body Weight from Linear Body Measurements in Two Commercial Meat-Type Chicken. Global J. Agric. Sci. 7(1): 57-59. http://dx.doi.org/10.4314/gjass.v7i1.2361.
  • Alabi OJ, Egena SSA, Ng’ambi JW, Norris D (2012). Comparative study of the three indigenous chicken breeds of South Africa: Body weight and linear body measurements. Agric. J. 7: 220–225. https://doi.org/10.3923/aj.2012.220.225
  • Aytekin C, Ni X, Cricri F, Aksu E (2018). Clustering and Unsupervised Anomaly Detection with l2Normalized Deep Auto-Encoder Representations. International Joint Conference on Neural Networks (IJCNN), Rio de Janeiro 1-6. https://doi.org/10.1109/IJCNN.2018.8489068
  • Breiman L, Friedman JH, Olshen RA, Stone CI (1984). Classification and regression trees. Belmont, Calif.: Wadsworth.
  • Dessie T, Gatachew F (2016). The Potchefstroom koekoek breed. International livestock research institute.DOI: 10.18805/ijar.B-944. https://doi.org/10.18805/ijar.B-944
  • Dzungwe JT, Gwaza DS, Egahi JO (2018). Statistical Modeling of Body Weight and Body Linear Measurements of the French Broiler Guinea Fowl in the Humid Tropics of Nigeria. Poult. Fish Wildlife Sci. 6: 197.
  • Egena SSA, Ijaiya AT, Ogah DM, Aya VE (2014). Principal component analysis of body measurements in a population of indigenous Nigerian chickens raised under extensive management system. Slovak J. Anim. Sci. 47:77–87.
  • Eyduran E, Celik S, Keskin I, Sevgenler H, Iqbal F (2019). Classification and regression tree (cart) algorithm (siniflama ve regresyon ağaci algoritmasi). Department of Animal Sciences and Business Administration, Igdir Universitesi, Turkey. https://www.researchgate.net/publication/335432905_CLASSIFICATION_AND_REGRESSION_TREE_CART_ALGORITHM_SINIFLAMA_VE_REGRESYON_AGACI_ALGORITMASI.
  • Guèye EF, Ndiaye A, Branckaert RDS (1998). Prediction of body weight on the basis of body measurements in mature indigenous chickens in Senegal. Lives. Res. Rural Dev. 10(3).
  • Heit HL (2017). A comparative study on the growth rates between outbred and inbred chickens. Int. Livest. Res. Institute.
  • Ige AO, Mudasiru IT, Rafiu BR (2016). Effect of Genotype on Growth Traits Characteristics of Two Commercial Broiler Chickens in a Derived Savannah Zone of Nigeria. Int. J. Res. Studies Agric. Sci. 2: 26–32. https://doi.org/10.20431/2454-6224.0206004
  • Latshaw JD, Bishop BL (2001). Estimating body weight and body composition of chickens by using non-invasive measurements. Poult. Sci. 80:868–873. https://doi.org/10.1093/ps/80.7.868
  • Mohammad MT, Rafeeq M, Bajwa MA, Awan MA, Abbas F, Waheed A, Bukhariand FA, Akhtar P (2012). Prediction of body weight from body measurements using regression tree (rt) method for indigenous sheep breeds in Baluchistan, Pakistan. J. Anim. Plant Sci. 22(1).
  • Mutibvu T, Chimonyo M, Halimani TE (2019). Effect of strain, sex and rearing among carcass and fat yield of Naked Neck, Ovambo and Potchefstroom Koekoek chickens. Indian J. Anim. Res. 2: 50–53. https://doi.org/10.18805/ijar.B-944
  • Oguntunji AO (2017). Regression tree analysis for predicting body weight of Nigerian muscovy duck (Cairina moschata). Genetika. 49 (20):743–753. https://doi.org/10.2298/GENSR1702743O
  • Ojedapo LO (2013). Evaluation of body weight and other linear body parameters of Marshall Broiler for Repeatability estimates. Int. J. Appl. Agric. Res. 9: 175–181.
  • Ojo V, Fayeye TR, Ayorindeamd KL, Olojede H (2014). Relationship between Body Weight and Linear Body Measurements in Japanese quail. J. Sci. Res. 6: 175–183. https://doi.org/10.3329/jsr.v6i1.16368
  • Olawumi S, Farinnako A (2017). Evaluation of the Relationship between Body Weight and Linear Measurements in West African Dwarf Goat as Influenced by Sex and Agro-vegetational Zone. Sci. Int. 5(2): 63–67. https://doi.org/10.17311/sciintl.2017.63.67
  • Semakula J, Lusembo P, Kugonza DR, Mutetikka D, Ssennyonjo J, Mwesigwa M (2011). Estimation of live body weight using zoometrical measurements for improved marketing of indigenous chicken in the Lake Victoria basin of Uganda. Liv. Res. Rural Dev. 23(8).
  • Tyasi TL, Norris D, Ng’ambi JW, Mabelebele M (2019). Combining abilities and heterosis of body weight in a diallel cross from three South African indigenous chicken genotypes. Appl. Ecol. Environmen. Res. 17(4): 9717–9723 https://doi.org/10.15666/aeer/1704_97179723.
  • Tyasi TL, Qin N, Jing Y, Mu F, Zhu H, Liu D, Yuan S, Xu R (2017). Assessment of relationship between body weight and body measurement traits of indigenous Chinese Dagu chickens using path analysis. Indian J. Anim. Res. 51(3):588–593.
  • Tyasi TL, Xu R, Chen X, Niu X, Qin N, Sun X, Zhang F, Zhu H (2018). Prediction of carcass weight from body measurement traits of Chinese indigenous Dagu male chickens using path coefficient analysis. Indian J. Anim. Sci. 88: 744–748.
  • Ukwu HO, Okoro VMO, Nosine RJ (2014). Statistical modeling of body weight and linear body measurements in Nigerian indigenous chicken. J. Agric. Vet. Sci. 7: 27–31. https://doi.org/10.9790/2380-07152730
  • Yakubu A (2012). Application of regression tree methodology in predicting the body weight of Uda sheep. Anim. Sci. Biotech. 45 (2): 484–490.
  • Yakubu A, Salako AE (2009). Path coefficient analysis of body weight and morphological traits of Nigerian indigenous chickens. Egyptian Poult. Sci. 29 (3): 837-850.
  •  

     

     

     

    Advances in Animal and Veterinary Sciences

    September

    Vol. 12, Iss. 9, pp. 1622-1845

    Featuring

    Click here for more

    Subscribe Today

    Receive free updates on new articles, opportunities and benefits


    Subscribe Unsubscribe