Next Article in Journal
A Descriptive Study of the Carpal Joint of Healthy Donkeys Using Ultrasonography, Computed Tomography, and Magnetic Resonance Imaging
Next Article in Special Issue
Effects of Backfat Thickness on Oxidative Stress and Inflammation of Placenta in Large White Pigs
Previous Article in Journal
Changes in Neuropeptide Prohormone Genes among Cetartiodactyla Livestock and Wild Species Associated with Evolution and Domestication
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Veterinary Sciences
Volume 9
Issue 5
10.3390/vetsci9050248
vetsci-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Photoperiodic Modulation in Immune and Reproductive Systems in Japanese Quails (Coturnix japonica): A Morphometric Perspective

by
Khizar Hayat
1,2,3,*,
Ali Raza
4,†,
Aitzaz Anas
5,
Anas Sarwar Qureshi
2,
Sarmad Rehan
2,
Ameer Hamza Rabbani
6,
Hafiz Faseeh ur Rehman
3,
Abdul Ghaffar Qamar
7,
Tayyab Rehman
2,
Farah Deeba
4 and
Amber Salman
8
1
Jockey Club College of Veterinary Medicine and Life Sciences, City University of Hong Kong, Hong Kong 999077, China
2
Department of Anatomy, University of Agriculture, Faisalabad 38000, Pakistan
3
Department of Anatomy and Histology, University of Veterinary and Animal Sciences, Lahore 54000, Pakistan
4
Department of Clinical Medicine and Surgery, University of Agriculture, Faisalabad 38000, Pakistan
5
House Office, Independent Medical College, Faisalabad 38000, Pakistan
6
Department of Surgery, Cholistan University of Veterinary and Animal Sciences, Bahawalpur 63100, Pakistan
7
Department of Clinical Sciences, Riphah College of Veterinary Sciences, Lahore 54000, Pakistan
8
University Medical and Dental College, The University of Faisalabad, Faisalabad 38000, Pakistan
*
Author to whom correspondence should be addressed.
†
Current Address: Centre for Animal Science, Queensland Alliance for Agriculture and Food Innovation, University of Queensland, St. Lucia, QLD 4067, Australia.
Vet. Sci. 2022, 9(5), 248; https://0-doi-org.brum.beds.ac.uk/10.3390/vetsci9050248
Submission received: 13 April 2022 / Revised: 11 May 2022 / Accepted: 18 May 2022 / Published: 23 May 2022
(This article belongs to the Special Issue Advances in Animal Reproductive Physiology, Reproductive Endocrine and Obstetric Diseases)
Browse Figures
Versions Notes

Abstract

:
The present study was designed to elucidate a relationship between lymphoid organs and reproductive activity in male Japanese quails (Coturnix japonica) bred in a temperate region of Pakistan (30.3753° N, 69.3451° E) in response to photoperiodic changes. The research focused primarily on the relative morphological changes in primary (thymus and bursa of Fabricius) and secondary (spleen) lymphoid organs with respect to seasonal variations in the histomorphometry of testicular tissue. For this purpose, a comparable number of clinically healthy Japanese quails were exsanguinated during active (April–May), regressive (September–October) and inactive (January–February) reproductive phases. Following an extensive gross measurement of lymphoid and reproductive organs, a histomorphometric analysis was performed on sampled tissues by employing ImageJ® software. Blood was collected for hormonal and leukocytic analysis. One-way ANOVA was used for statistical comparison. Testes had the highest parenchymal development in the active phase (80.66 ± 21.22 µm) and the lowest in the inactive phase (27.80 ± 7.22 µm). Conversely, a percentage change was evident in the sizes of primary (bursa: 61.5%, thymus: 46.9%) and secondary (spleen: 23.9%) lymphoid organs during inactive and active reproductive phases. This study demonstrated that a physiological trade-off is imperative between immune and reproductive systems for optimum survivability and reproductive performance.

1. Introduction

The immune and reproductive systems are two major body systems that require extensive resource allocation therefore, trade-offs must be made between them during different physiological states. This phenomenon is particularly evident in seasonal breeders, such as birds [1]. Previous studies have demonstrated that androgens are primarily responsible for immune and reproductive changes in birds subsequent to neuroendocrine stimulation [2]. However, it is still unclear whether this phenomenon is limited to males only or whether estrogen could also play an important role in immune regression during a phase of higher reproductive activity [3,4].
Primary immune organs have been a focus of investigators as an integral part of the trade-off between immune and reproductive systems [5,6]. Despite being the major secondary immune organ, the spleen has been ignored altogether over the years. The spleen’s role is particularly important in birds because it acts as a pure lymphoid organ owing to a greater proportion of white pulp (containing aggregated lymphoid tissue) than red pulp [7,8]. Conversely, in mammals, red pulp comprises a greater portion of the spleen considering its auxiliary function in erythrocyte storage [9,10]. A comprehensive study was required to investigate if the spleen followed a similar trend as primary lymphoid organs. Most studies in this regard have focused on the physiological manifestations attributed to a compromised immune status. The authors in the current study employed the susceptibility toward parasitic infestation as an indirect indication of reduced immunogenic capabilities. It is widely accepted that blood and intestinal parasites could prove to be reasonable indicators of avian immune status. An increased level of parasitic infestations has been observed in birds during their active breeding season [11,12,13]. An anatomical basis of this physiological variation was analyzed among wild birds in a prior study [5]. However, limited literature has been published whereby both gross and histomorphometric parameters were studied concomitantly. Moreover, minimal information has been made available pertaining to the role of hypertrophied parenchyma relative to the organ stroma during active reproductive seasons.
Birds indigenous to tropical regions are known to undergo annual photoperiodic changes causing a reduced level of testosterone in plasma. This physiological change enhances immunity and depresses overall reproductive functioning during reproductively inactive season [14,15]. It has been reported that melatonin can enhance immune status of birds but at this time, the extent of this increase in melatonin countering the testosterone-mediated decrease in immunity during active reproductive season is not well-understood [16]. There has been some discussion on the number of reproductive phases in an annual breeding season [17]. Keeping in mind relatively recent studies, birds are thought to undergo two major transitions in one cycle with three distinct reproductive states, which are active, regressive and inactive [6]. The Japanese quail was selected for this study as a model bird because of its important role as an essential poultry bird for meat and egg production evolved over the last few decades [18,19].
Most prior studies relating to variations in the immune status of birds in relation to reproductive phases were conducted on wild birds [16,20,21,22]. Recent studies report a negative impact of stress on immune status causing its fairly mild regression [23,24]. Therefore, a study focused on investigating stress-associated immune compromises in commercially farmed birds is imperative to confirm the aforementioned trend in domesticated birds. The current study aimed to determine a seasonal relationship between immune and reproductive statuses of Japanese quails. To remove the impact of environmental factors and influences associated with breeding seasonality, captive birds were inducted in the present study. It was hypothesized that the immune and reproductive systems make a trade-off at a morphometric level during the annual breeding cycle. To observe this change in immune status, morphometric parameters of primary and secondary lymphoid organs were investigated in relation to a variation in testicular morphology, as a proxy for breeding seasonality in domestic birds. To support our observations, hematological and hormonal analyses were performed to correlate the morphometric changes of immune and reproductive organs with immune cells in blood and plasma testosterone levels.

2. Materials and Methods

2.1. Sample Collection

This study was conducted at the Department of Anatomy, University of Agriculture Faisalabad, Pakistan after approval from the Institutional Biosafety and Bioethics Committee (IBC, DGS No. 4141-44). Commercially available adult male birds (n = 50) were decapitated in each reproductive phase, namely active breeding season (Spring; April–May), regressive season (Autumn; September–October) and inactive season (Winter; January–February). To minimize the age variation among adult birds, feather analysis (primary and secondary coverts) was performed along with obtaining complete history (date of purchase of quail chicks) from the farmer [25]. After purchase, birds were transferred to the cages in animal house for 7 days for acclimatization. During this period, optimum environmental conditions such as light (5 Lux, 15.5 h per day), humidity (65%), temperature (25–27 °C) and ventilation were provided to minimize stress [26,27,28]. Feed and water were offered to the birds ad libitum. After acclimatization period, birds were decapitated, and organs were removed from the carcass via conventional postmortem technique. During slaughtering, blood sampling was also performed for lymphocytic and hormonal analysis. Organs were first washed in normal saline and then preserved in fixatives. Bursa of Fabricius, thymus and spleen were preserved in 10% aqueous solution of buffered formalin (10% NBF), while testes were preserved in Bouin’s solution for effective and rapid fixation [29]. Organs were kept in the fixatives for at least 72 h until further processing.

2.2. Hematological Analysis

Blood was collected in two separate vacutainers containing ethylenediaminetetraacetic acid (EDTA) and coagulation gel (Bolton®, MIDAS, Karachi, Pakistan). Lymphocyte percentage was counted after making a thin smear of uncoagulated blood under the microscope. EDTA-mixed blood was used for white blood cell enumeration with Natt–Herrick (N&H) method for total leukocyte count [30]. Blood was mixed with N&H solution in a micropipette and then poured in Neubar’s hemocytometer for analysis under the microscope at 400X magnification.

2.3. Gross Morphometry and Histological Sample Preparation

Gross morphometric values were measured before putting the organs in preservative fluids to avoid any variation from the in vivo state. Weight of the organs was measured via digital weighing balance, while length, width and diameter were measured using vernier calipers (VCs). The following formula was used to measure the relative weight of the organs in relation to the total weight of individual:
Relative Weight (RW) of organ: 100/Weight of bird × Absolute Weight (AW) of organ.
Zero error of VCs was removed prior to measuring the organs. After gross morphometric analysis, organs were immersed in fixative for 72 h. Conventional paraffin embedding technique was used for preparation of microscopic slides. Organs were dehydrated via increasing gradients of alcohol from 70% to 100%. Clearing of this alcohol was performed first using a mixture of xylene and alcohol (xylol) and later with 100% xylene. Paraffin infiltration was performed after melting of paraffin in hot air oven. Finally, paraffin embedding was conducted, and blocks were left overnight for further processing. Tissues were sectioned on manual rotary microtome (Microm HM315®, Rankin Biomed, Holly, MI, USA) at 4–6 µm thickness. After mounting of tissues on microscopic slides, a hotplate was used to dewax the tissue. Routine staining was performed using Hematoxylin and Eosin solution to stain the tissues on microscopic slides, and coverslips were mounted using DPX® (Sigma-Aldrich, Karachi, Pakistan) [29].

2.4. Histomorphometric Analysis

Compound light microscope (Anti-Mould SE®, Nikon, Tokyo, Japan) with camera was used to obtain the images for histomorphometric analysis of samples, which was performed manually using computerized software ImageJ (Research Services Branch, National Institute of Mental Health, Bethesda, MD, USA) [31]. Software was calibrated using a microscopic image of the stage micrometer. Multiple measurements were taken per slide to minimize human error. All measurements were saved, and mean values were calculated for graphical representations.

2.5. Hormonal Analysis

Serum was harvested from blood in gel vacutainers and stored at −20 °C for further analysis. Radioimmunoassay (RIA) was used for the quantification of serum testosterone levels using a commercially available kit (AccuBind ELISA Microwells Testosterone Kit®, Monobind, Lake Forest, CA, USA). Briefly, 10 µL of serum samples were taken in designated well in which 50 µL of testosterone enzyme reagent containing 0.7 mL of testosterone in steroid conjugate buffer. After 30 s, 50 µL of testosterone biotin solution was added to those wells, and the microplate was incubated at 27 °C for one hour. After decanting the microplate and drying with absorbent paper, it was washed with 350 µL of washing buffer three times. In each well, 100 µL of 3,3′,5,5′-tetramethylbenzidine substrate (TMB substrate) was added, and the microplate was incubated at 27 °C for 15 min. After incubation, 50 µL of Stop Solution was mixed, and the plate was shaken for 15–20 s. Finally, the microplate was read at 450 nm with reference filter set at 630 nm wavelength using a digital microplate reader (SHE3000®, Perfemed, Chemux Bioscience Inc., South San Francisco, CA, USA).

2.6. Statistical Analysis

Data were statistically analyzed using SPSS Version 21® (IBM, Armonk, NY, USA). One-way analysis of variance (ANOVA) was used with post hoc Tukey’s test for the analysis. Confidence interval (CI) was maintained at 95% (p < 0.05).

3. Results

3.1. Anatomical Parameters

The absolute mass of the testes was significantly highest in the active reproductive phase (1.28 ± 0.16 g) as compared to the regressive (0.46 ± 0.12 g) and inactive phases (0.44 ± 0.14 g) with the latter ones differing insignificantly from each other. Interestingly, the relative weight of the testes differed significantly during all three seasons (1.74 ± 0.17%, 1.2 ± 0.12% and 0.55 ± 0.13%, respectively). Similar results were observed for the length and width of the testes where all three phases demonstrated significant changes from each other (Table 1).
Contrarily, the absolute mass of the bursa of Fabricius was highest in the inactive breeding phase (0.26 ± 0.05 g), while it was statistically indifferent during the regressive (0.19 ± 0.06 g) and active breeding seasons (0.18 ± 0.04 g). The relative weight was highest in the inactive breeding season (0.36 ± 0.04%) and differed significantly during all three seasons. The length and width of the bursa followed a similar trend as that of the absolute weight and relative weight, respectively (Figure 1).
The absolute weight and relative weight of the thymus followed an identical pattern where inactive, regressive and active reproductive phases differed significantly from each other (Table 1).
The spleen showed a similar trend where the absolute weight was statistically different during all three reproductive phases, i.e., active, regressive and inactive (Figure 1). The relative weight was highest in the inactive breeding phase (0.60 ± 0.07%); however, it differed insignificantly during the regressive and active seasons. The diameter of the spleen was significantly different during all three seasons (Table 1).

3.2. Histological Parameters

Histological examination of the testes showed that the diameter of seminiferous tubules (203.37 ± 65.22 µm) and the height of the spermatogenic epithelium (80.66 ± 21.22 µm) were significantly greater during the active reproductive season as compared to the regressive and inactive seasons. On the other hand, the luminal diameter of the seminiferous tubules was significantly different in the regressive season (114.4 ± 15.93 µm) as compared to the nonactive season (63.24 ± 19.11 µm). The difference was insignificant between regressive and active breeding phase (114.4 ± 15.93 µm and 122.71 ± 9.98 µm, respectively) (Figure 2). Contrarily, the thickness of the tunica albuginea was significantly higher in the inactive reproductive phase (191.10 ± 26.85 µm) as compared to the regressive and active phases (Table 2, Figure 3). A comparative analysis of histomorphometric parameters of the spleen demonstrated an opposite trend to that of the testes. The diameter of the splenic nodule (426.52 ± 29.47 µm) and thickness of capsule (82.19 ± 14.57 µm) were significantly higher in the nonactive reproductive phase, while they remained statistically non-significantly different during the regressive and active seasons (Figure 4). The histological analysis of the bursa showed that the diameter of the follicle and thickness of the capsule were significantly different during the three reproductive seasons. The diameter of the bursa follicle was highest in the nonactive reproductive season (508.52 ± 13.06 µm), while it was the lowest in the active reproductive season (195.5 ± 19.54 µm). Contrarily, the diameter of the bursa capsule was found to be highest in the active reproductive season (213.46 ± 23.49 µm), while it was the lowest in nonactive season (48.51 ± 7.94 µm). The height of the epithelium showed no significant difference (Figure 5).
The diameter of the thymic lobule differed significantly during the three reproductive seasons (p ≤ 0.05). The highest values were observed during the inactive reproductive phase (199.01 ± 7.03 µm), while the lowest value was observed during the active reproductive phase (105.63 ± 2.59 µm) with the regressive season (156.73 ± 5.47 µm) showing a gradual change in the diameter of the thymic lobule. In contrast, the difference in diameter of the thymic corpuscle was insignificantly different during all three reproductive seasons (Table 2).

3.3. Hematological Parameters

Immunological parameters such as Total Leukocyte Count (TLC) and Lymphocyte Count (LC) were significantly different during the three reproductive phases with the highest number being in the inactive breeding phase (189.3 ± 6.53 k/µL, 137.18 ± 4.92 k/µL) and the lowest number in the active breeding phase (Table 3).

3.4. Hormonal Analysis

The testosterone level was detected as being significantly higher during the active reproductive phase (2.45 ± 0.035 ng/mL) with a gradual reduction in the regressive phase, and the lowest values were seen in the nonactive reproductive phase (Table 3).

4. Discussion

Our study revealed that birds experienced the highest level of morphological development with regards to their reproductive organs during the active breeding phase. Similar annual changes have been reported in Coturnix japonica, Perdicula asiatica and Numidia meleagris [7,9,18]. The decrease in weight of the testes was significant as the season changed from active to regressive. However, the difference was relatively less during the inactive season. Conversely, the length and width of the testes decreased at an almost comparable rate during both physiological transitions. This variation in rate of decrease in the morphometric parameters of the testes was interesting to observe in terms of the functionality of the organ. Increased gonadal activity during the active reproductive season ensures that there would be maximum production of active spermatozoa, which can be used for a higher frequency of copulation [32]. The authors postulate that this variation in percent decline of weight during the two transitions was due to the presence of active spermatozoa in the seminiferous tubules during the active phase, which contributes toward the weight of the organ. An analysis of the photomicrographs shows that spermatozoa reach their minimum number during the regressive stage and disappear completely in the inactive season. After the regressive season, the number of spermatozoa in the seminiferous tubules decrease rapidly as the testicular size shrinks. The interstitial tissue has been shown to continue regressing at almost the same level which in turn causes a further decrease in the length and width of the organ [33,34].
This hypothesis was further tested and proven in the histomorphometric analysis of testicular tissue in which the height of the germinal epithelium (GE) was highest during the active breeding phase as opposed to the regressive and inactive ones. However, the decrease in the height of the spermatogenic epithelium (SE) was almost double when moving from the active to the regressive phase as compared to moving from the regressive to inactive phase (53% and 27%, respectively). This supports our initial hypothesis that due to meiotic activity, spermatogenic hyperplasia occurs in the testes, which adds to the weight of the organ exponentially. An increase in testicular parenchyma is directly related to an increase in the weight of the organ. These results were eerily similar to the ones reported in other avian species [5,35]. Moreover, we also observed that the diameter of the seminiferous tubules (STs), in keeping with the overall trend, differed slightly where the percent decrease in the first transition was relatively less than the second one, i.e., 25% and 39%, respectively. Due to the accumulation of free spermatozoa in the lumen of the STs and the hyperplasia of spermatogenic cells, the basement membrane of the STs may be stretched and even after the active season, it does not return to its normal stage immediately [36]. This leads to a relative slower decrease in the diameter of the STs as compared to the height of GE. The diameter of the SE lumen and thickness of the tunica albuginea followed the exact opposite trend where the active season showed the minimal values compared to the other two. It is understandable that due to an increase in the epithelial height, the lumen was decreased; however, as Figure 2 shows, the rate of increase in the luminal diameter was far greater during the first transition as compared to the second one. This corroborates the previously discussed phenomenon where a delay in regression/return of the basement membrane to its normal state was observed [33]. This theory needs further testing; therefore, the authors suggest an investigation employing special stains to differentiate between the types of connective tissue present in the basement membrane and their relative changes observed during different phases of the breeding season.
Morphometric parameters of the lymphoid organs showed a completely opposite trend, whereby the highest values were observed during the inactive breeding phase. The relative increase in the weight of the spleen and bursa was higher in the second transition as compared to the first one. Contrarily, an increase in the size of the thymus was relatively greater in the first seasonal change as compared to the second one. Similar trends were observed with other morphometric parameters of immune organs such as the length and width of the bursa and the diameter of spleen where the inactive breeding state showed the highest values as compared to the other two. These results agree with those observed in different species, such as willow tits [37]. Comparable results were reported for the splenic parameters of Indian tropical bush quails (Perdicula asiatica) as well [16]. In birds, the seasonal regression of the thymus is relatively more pronounced as compared to the bursa of Fabricius. During the active breeding season, birds can afford to reduce the thymic lymphoid tissue even beyond acceptable thresholds in certain dire situations. Consequently, an abrupt change in thymic weight is not in itself a measure of its activity, rather it is an organ’s attempt to return to its normal size after a significant decrease during the active breeding phase [38].
An organized immune system is present in specialized immune organs such as the bursa of Fabricius (bursal follicle), thymus (thymic lobule) and spleen (splenic nodule). Only one study has been reported regarding the histomorphometry of the thymus during different breeding seasons in any adult avian species [38]. Overall, the lymphoid tissue followed a similar trend in which the highest values were observed during the inactive season, and the lowest ones were in the active season. These results corroborated findings of a prior study conducted in Numidia meleagris [5,6]. However, unlike our present observations, the percent increase in their morphometric values varied significantly for all three organs. Both the thymic lobule and splenic nodule exhibited similar behavior with a steep decline in size (diameter) during the transition from the active to the regressive state (32% and 15%, respectively). However, a relatively lower rate was observed during the transition from the regressive to the inactive state (21% and 10%, respectively). This may correlate with the resource allocation from the reproductive system to the immune system. Immediately after the active breeding season, birds need to boost their immune system to cope with the climatic stress and increase their survivability for the next active season. To achieve this, a rapid reallocation of resources occurs, and thus we see a sharp increase in the lymphoid content. As the body moves from the regressive to the inactive season, other external factors, such as nestling stress, start cutting into the energy consumption thus causing the rate of increase in lymphatic aggregates to plumate [39]. Contradictory results were observed in the case of the bursa of Fabricius where a higher percentage change in bursal diameter was observed during the second transition as compared to the first one (46% and 28%, respectively). It may be related to the fact that lymphocyte production keeps on increasing in the bursa as it is the major primary lymphoid organ even in the inactive season to maintain the lymphocyte number in the body for higher chances of survivability. The same phenomenon can be observed in migratory birds where a significant decrease in the size of the spleen is observed during the migratory season [40]. However, this theory needs to be tested through a comparative analysis of the lymphocyte–lymphoblast population in the bursal follicles during all three phases. The capsule of the lymphoid organs is made up of connective tissue and is considered a stromal component as compared to the lymphoid aggregates, which are parenchymal in nature. The capsule thickness also varied in the three reproductive phases; however, this change was not identical to that of the parenchymal change during the same period, which leads us to believe that both vary independently from each other during the same period. Moreover, the change in the parenchymal components was much higher and more pronounced as compared to the stromal components, which supports our initial theory that it is the parenchymal content that shows major plasticity during different reproductive phases [41].
It was also interesting to compare the morphometric changes in primary and secondary lymphoid organs. As adult birds were used for this study, we were expecting to have a greater degree of morphometric variation (especially parenchymal changes) in the spleen as compared to the thymus and bursa. The reason behind this was that primary lymphoid organs tend to have a greater degree of functionality during the early years of life as compared to later ones [42]. However, this study revealed that the degree of variation was much higher in the case of primary lymphatic organs as compared to secondary ones. This is consistent with previous findings in which a variation was observed in the regression of primary immune organs in different avian species with some showing more regression compared to others [43]. The bursal and thymic lymphatic lobules doubled their diameter from the active season to inactive season, while the change in the splenic lobule was only around 75%. Contrarily, we also observed that the decrease in the lymphoid content of the thymus and bursa was much more pronounced during the active reproductive phase. We think this might be due to the difference in the types of cells present in these organs [44]. However, this needs further investigation where advanced techniques such as immunohistochemistry (IHC) can be used for the differentiation of these cells.
One of the factors that researchers explored while trying to analyze the immune status of wild birds is the extent of parasitic infestations (intestinal and blood infections) [11,12]. Previous studies have shown that wild birds face a higher level of parasitemia during the active breeding phase [6,16]. This conforms to the hypothesis that not only the lymphoid organs, but also the peripheral status of immunity undergoes major regression during higher reproductive activity. Similar results were seen in our study where the lymphocyte percentage and TLC were lowest in the active breeding season. Increased levels of leukocytes including lymphocytes during the inactive reproductive phase suggest that not only the initial responders of the immune response (neutrophils) but also the artillery (lymphocytes) is enhanced. These higher levels are required for proper bloodborne immune activity [45].
The results in this study confirm previous reports suggesting that the testosterone level changes with climatic conditions due to activation or regression of the gonads during the sexually active and inactive seasons, respectively [16,45,46]. The results reveal that the serum testosterone level was significantly higher during the active reproductive season as compared to the nonactive reproductive season. This further corroborated the initial hypothesis that an increase in testosterone was directly related to the lower levels of immune parameters.
Immunohistochemical (IHC) techniques could be employed in future endeavors to observe the differences in the ratios of different immune cells and the variation in the number of lymphocytes during active breeding seasons. There are some reports that state that this phenomenon is not readily observable in female birds. Therefore, future studies could also involve female birds to investigate the role of estrogen in controlling the morphometric values of immune organs.

5. Conclusions

The study revealed that Japanese quails follow a breeding cycle with varying degrees of reproductive activity such as active, regressive and inactive reproductive phases. This is evident both from the morphometric parameters and from the hormonal analysis. This change in reproduction is directly related to the immune status of birds. A trade-off between the reproductive status and immune statuses of birds must be made due to finite energy reserves in the body where an increase in activity of one system leads to a decrease in another and vice versa. Moreover, it was also observed that primary and secondary lymphoid organs do not show the same level of morphometric changes during this seasonal variation. A variation in morphometry was much more pronounced in primary lymphoid organs as compared to secondary ones. In addition, this variation is not limited to the morphometric parameters of organs only. It was also observed that the number of circulating immune cells (blood-borne) also follow a similar trend with the highest number being in the nonactive reproductive phase and vice versa.

Author Contributions

Conceptualization, K.H., A.S.Q. and S.R.; Formal analysis, K.H., A.H.R. and A.A.; Investigation, K.H., A.R., A.G.Q. and T.R.; Methodology, A.H.R.; Resources, A.S.Q., F.D., A.A. and A.S.; Supervision, A.S.Q., A.R. and S.R.; Visualization, A.H.R., H.F.u.R. and A.S.; Writing—original draft, K.H.; Writing—review and editing, A.S.Q., A.R. and S.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

The animal study protocol was approved by the Institutional Biosafety and Bioethics Committee (IBC) of the University of Agriculture, Faisalabad, Pakistan.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data generated in this study is already added in the tables of this article. If you need any further information, please feel free to contact authors.

Acknowledgments

The authors would like to thank Michael Doube (Jockey Club College of Veterinary Medicine and Life Sciences, City University of Hong Kong, Hong Kong SAR) for final editing and proofreading of the manuscript. Moreover, we would like to acknowledge the faculty members in the Department of Anatomy, at the University of Agriculture, Faisalabad. Pakistan for their valuable input and contribution.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Williams, T.D. Physiological Adaptations for Breeding in Birds; Princeton University Press: Princeton, NJ, USA, 2012; ISBN 978-1-4008-4279-7. [Google Scholar]
  2. Ruiz, M.; French, S.S.; Demas, G.E.; Martins, E.P. Food supplementation and testosterone interact to influence reproductive behavior and immune function in Sceloporus graciosus. Horm. Behav. 2010, 57, 134–139. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Hasselquist, D. Comparative immunoecology in birds: Hypotheses and tests. J. Ornithol. 2007, 148, 571–582. [Google Scholar] [CrossRef]
  4. Koutsos, E.A.; Klasing, K.C. Chapter 17—Factors Modulating the Avian Immune System. In Avian Immunology, 2nd ed.; Schat, K.A., Kaspers, B., Kaiser, P., Eds.; Academic Press: Boston, MA, USA, 2014; pp. 299–313. ISBN 978-0-12-396965-1. [Google Scholar]
  5. Akbar, Z.; Qureshi, A.S. Effects of Seasonal Variation in Different Reproductive Phases on the Cellular Response of Bursa and Testes in Japanese Quail (Coturnix japonica). Pak. Vet. J. 2012, 32, 525–529. [Google Scholar]
  6. Ali, M.Z.; Qureshi, A.S.; Rehan, S.; Akbar, S.Z.; Manzoor, A. Seasonal variations in histomorphology of testes and bursa, immune parameters and serum testosterone concentration in male guinea fowl (Numida meleagris). Pak. Vet. J. 2015, 35, 88–92. [Google Scholar]
  7. John, J.L. The avian spleen: A neglected organ. Q. Rev. Biol. 1994, 69, 327–351. [Google Scholar] [CrossRef]
  8. Smith, K.G.; Hunt, J.L. On the use of spleen mass as a measure of avian immune system strength. Oecologia 2004, 138, 28–31. [Google Scholar] [CrossRef]
  9. Brendolan, A.; Rosado, M.M.; Carsetti, R.; Selleri, L.; Dear, T.N. Development and function of the mammalian spleen. Bioessays 2007, 29, 166–177. [Google Scholar] [CrossRef]
  10. Pabst, R. Plasticity and heterogeneity of lymphoid organs: What are the criteria to call a lymphoid organ primary, secondary or tertiary? Immunol. Lett. 2007, 112, 1–8. [Google Scholar] [CrossRef]
  11. Møller, A.P.; Erritzøe, J.; Saino, N. Seasonal Changes in Immune Response and Parasite Impact on Hosts. Am. Nat. 2003, 161, 657–671. [Google Scholar] [CrossRef] [Green Version]
  12. Dunn, J.C.; Goodman, S.J.; Benton, T.G.; Hamer, K.C. Avian blood parasite infection during the non-breeding season: An overlooked issue in declining populations? BMC Ecol. 2013, 13, 30. [Google Scholar] [CrossRef] [Green Version]
  13. Cornelius, J.; Zylberberg, M.; Breuner, C.; Gleiss, A.C.; Hahn, T. Assessing the role of reproduction and stress in the spring emergence of haematozoan parasites in birds. J. Exp. Biol. 2014, 217, 841–849. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Majewski, P.; Markowska, M.; Pawlak, J.; Piesiewicz, A.; Turkowska, E.; Skwarlo-Sonta, K. Pineal gland and melatonin: Impact on the seasonality of immune defence in mammals and birds. Adv. Neuroimmune Biol. 2012, 3, 95–108. [Google Scholar] [CrossRef]
  15. Lopes, L.E.; Teixeira, J.P.G.; Meireles, R.C.; Bastos, D.S.S.; de Oliveira, L.L.; Solar, R. High seasonal variation of plasma testosterone levels for a tropical grassland bird resembles patterns of temperate birds. Physiol. Biochem. Zool. 2021, 94, 143–151. [Google Scholar] [CrossRef] [PubMed]
  16. Singh, S.S.; Haldar, C. Melatonin prevents testosterone-induced suppression of immune parameters and splenocyte proliferation in Indian tropical jungle bush quail, Perdicula asiatica. Gen. Comp. Endocrinol. 2005, 141, 226–232. [Google Scholar] [CrossRef]
  17. Aire, T.A. Morphological changes in the efferent ducts during the main phases of the reproductive cycle of birds. J. Morphol. 2002, 253, 64–75. [Google Scholar] [CrossRef]
  18. Nasar, A.; Rahman, A.; Hoque, N.; Talukder, A.K.; Das, Z.C. A survey of Japanese quail (Coturnix coturnix japonica) farming in selected areas of Bangladesh. Vet. World 2016, 9, 940–947. [Google Scholar] [CrossRef] [Green Version]
  19. Lukanov, H. Domestic quail (Coturnix japonica domestica), is there such farm animal? Worlds Poult. Sci. J. 2019, 75, 547–558. [Google Scholar] [CrossRef]
  20. Singh, S.S.; Haldar, C. Peripheral melatonin modulates seasonal immunity and reproduction of Indian tropical male bird Perdicula asiatica. Comp. Biochem. Physiol. Part. A Mol. Integr. Physiol. 2007, 146, 446–450. [Google Scholar] [CrossRef]
  21. Hau, M. Timing of Breeding in Variable Environments: Tropical Birds as Model Systems. Horm. Behav. 2001, 40, 281–290. [Google Scholar] [CrossRef]
  22. Pap, P.L.; Czirják, G.Á.; Vágási, C.I.; Barta, Z.; Hasselquist, D. Sexual dimorphism in immune function changes during the annual cycle in house sparrows. Naturwissenschaften 2010, 97, 891–901. [Google Scholar] [CrossRef]
  23. Kuhlman, J.R.; Martin, L.B. Captivity affects immune redistribution to skin in a wild bird. Funct. Ecol. 2010, 24, 830–837. [Google Scholar] [CrossRef]
  24. Martin, L.B.; Brace, A.J.; Urban, A.; Coon, C.A.C.; Liebl, A.L. Does immune suppression during stress occur to promote physical performance? J. Exp. Biol. 2012, 215, 4097–4103. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Leopold, A.S. Age Determination in Quail. J. Wildl. Manag. 1939, 3, 261. [Google Scholar] [CrossRef]
  26. Pereira, D.F.; Batista, E.D.; Sanches, F.T.; Gabriel Filho, L.R.; Bueno, L.G. Diferenças comportamentais de poedeiras em diferentes ambientes térmicos. Energ. Na Agric. 2015, 30, 33–40. [Google Scholar] [CrossRef] [Green Version]
  27. Raharjo, S.; Rahayu, E.S.; Purnomo, S.H. Factors affecting quail egg production under the changing climate at Kulonprogo Regency, Indonesia. IOP Conf. Ser. Earth Environ. Sci. 2018, 200, 012012. [Google Scholar] [CrossRef] [Green Version]
  28. Molino, A.B.; Garcia, E.A.; Santos, G.C.; Vieira Filho, J.A.; Baldo, G.A.A.; Almeida Paz, I.C.L. Photostimulation of Japanese quail. Poult. Sci. 2015, 94, 156–161. [Google Scholar] [CrossRef] [PubMed]
  29. Bancroft, J.D.; Gamble, M. Theory and Practice of Histological Techniques; Elsevier Health Sciences: Amsterdam, The Netherlands, 2008; ISBN 0-443-10279-1. [Google Scholar]
  30. Carisch, L.; Stirn, M.; Hatt, J.M.; Federer, K.; Hofmann-Lehmann, R.; Riond, B. White blood cell count in birds: Evaluation of a commercially available method. BMC Vet. Res. 2019, 15, 93. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Schneider, C.A.; Rasband, W.S.; Eliceiri, K.W. NIH Image to ImageJ: 25 years of image analysis. Nat. Methods 2012, 9, 671–675. [Google Scholar] [CrossRef]
  32. Kriegsfeld, L.J.; Ubuka, T.; Bentley, G.E.; Tsutsui, K. Seasonal control of gonadotropin-inhibitory hormone (GnIH) in birds and mammals. Front. Neuroendocr. 2015, 37, 65–75. [Google Scholar] [CrossRef] [Green Version]
  33. Aire, T.A.; Ozegbe, P.C. The testicular capsule and peritubular tissue of birds: Morphometry, histology, ultrastructure and immunohistochemistry. J. Anat. 2007, 210, 731–740. [Google Scholar] [CrossRef]
  34. Leska, A.; Kiezun, J.; Kaminska, B.; Dusza, L. Seasonal changes in the expression of the androgen receptor in the testes of the domestic goose (Anser anser f. domestica). Gen. Comp. Endocrinol. 2012, 179, 63–70. [Google Scholar] [CrossRef] [PubMed]
  35. Dixit, A.S.; Singh, N.S. Environmental control of seasonal reproduction in the wild and captive Eurasian Tree Sparrow (Passer montanus) with respect to variations in gonadal mass, histology, and sex steroids. Can. J. Zool. 2013, 91, 302–312. [Google Scholar] [CrossRef]
  36. Tae, H.; Jang, B.; Ahn, D.; Choi, E.; Kang, H.; Kim, N.; Lee, J.; Park, S.; Yang, H.; Kim, I. Morphometric studies on the testis of Korean ring-necked pheasant (Phasianus colchicus karpowi) during the breeding and non-breeding seasons. Vet. Res. Commun. 2005, 29, 629–643. [Google Scholar] [CrossRef] [PubMed]
  37. Silverin, B.; Fänge, R.; Viebke, P.-A.; Westin, J. Seasonal Changes in Mass and Histology of the Spleen in Willow Tits Parus montanus. J. Avian Biol. 1999, 30, 255–262. [Google Scholar] [CrossRef]
  38. Höhn, E.O. Seasonal recrudescence of the thymus in adult birds. Can. J. Biochem. Physiol. 1956, 34, 90–101. [Google Scholar] [CrossRef] [PubMed]
  39. Węgrzyn, E. Resource allocation between growth and endothermy allows rapid nestling development at low feeding rates in a species under high nest predation. J. Avian Biol. 2013, 44, 383–389. [Google Scholar] [CrossRef]
  40. Deerenberg, C.; Biebach, H.; Bauchinger, U. Spleen size variation during long-distance migration in the garden warbler Sylvia borin. Avian Sci. 2002, 2, 217–226. [Google Scholar]
  41. Thompson, H.L.; Smithey, M.J.; Surh, C.D.; Nikolich-Žugich, J. Functional and homeostatic impact of age-related changes in lymph node stroma. Front. Immunol. 2017, 8, 706. [Google Scholar] [CrossRef] [Green Version]
  42. Ciriaco, E.; Píñera, P.P.; Díaz-Esnal, B.; Laurà, R. Age-related changes in the avian primary lymphoid organs (thymus and bursa of Fabricius). Microsc. Res. Tech. 2003, 62, 482–487. [Google Scholar] [CrossRef]
  43. Mešťanová, V.; Varga, I. Morphological view on the evolution of the immunity and lymphoid organs of vertebrates, focused on thymus. Biologia 2016, 71, 1080–1097. [Google Scholar] [CrossRef]
  44. Kannan, T.A.; Ramesh, G.; Ushakumari, S.; Vairamuthu, S.; Sivakumar, M. Age related ultrastructural changes in lymphoid organs of Nandanam layer chicken (Gallus domesticus). J. Entomol. Zool. Stud. 2018, 6, 1417–1421. [Google Scholar]
  45. Bentley, G.E.; Demas, G.E.; Nelson, R.J.; Ball, G.F. Melatonin, immunity and cost of reproductive state in male European starlings. Proc. R. Soc. Lond. Ser. B Biol. Sci. 1998, 265, 1191–1195. [Google Scholar] [CrossRef] [PubMed]
  46. Abdul-Rahman, I.I.; Obese, F.Y.; Jeffcoate, I.A. Developmental changes in the histological structure of the testes, and testosterone profiles in male guinea fowls (Numida meleagris). Theriogenology 2017, 101, 114–122. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Vetsci 09 00248 g001 550
Figure 1. Graphical representation of gross morphometric values of organs (testes, bursa of Fabricius, thymus and spleen) in different phases of a reproductive season. AW: Average weight, RW: Relative weight, L: Length, W: Width and D: Diameter. All values are represented as means ± SEM, and parameters with different superscripts (a,b,c) are significantly different from each other (p < 0.05).
Figure 1. Graphical representation of gross morphometric values of organs (testes, bursa of Fabricius, thymus and spleen) in different phases of a reproductive season. AW: Average weight, RW: Relative weight, L: Length, W: Width and D: Diameter. All values are represented as means ± SEM, and parameters with different superscripts (a,b,c) are significantly different from each other (p < 0.05).
Vetsci 09 00248 g001
Vetsci 09 00248 g002a 550Vetsci 09 00248 g002b 550
Figure 2. Boxplot representing histomorphometric parameters in immune organs during active (A), regressive (B) and inactive breeding phases (C). All values are represented with their median, and upper and lower quartiles are shown in the above figure. Overlapping whiskers between different seasons show a relatively lower level of significance as compared to the ones in the which the difference is greater. Parameters with different superscripts (a,b,c) are significantly different from each other (p < 0.05).
Figure 2. Boxplot representing histomorphometric parameters in immune organs during active (A), regressive (B) and inactive breeding phases (C). All values are represented with their median, and upper and lower quartiles are shown in the above figure. Overlapping whiskers between different seasons show a relatively lower level of significance as compared to the ones in the which the difference is greater. Parameters with different superscripts (a,b,c) are significantly different from each other (p < 0.05).
Vetsci 09 00248 g002aVetsci 09 00248 g002b
Vetsci 09 00248 g003 550
Figure 3. Photomicrographs of testes of Japanese quails during different breeding phases of annual reproductive cycle. (a) Active breeding phase, (b) Regressive breeding phase and (c) Inactive breeding phase. Curved black arrows demonstrate the measurement of diameter; green arrows show the lumen of seminiferous tubules, and red arrows show height of spermatogenic epithelium. (H&E 200X).
Figure 3. Photomicrographs of testes of Japanese quails during different breeding phases of annual reproductive cycle. (a) Active breeding phase, (b) Regressive breeding phase and (c) Inactive breeding phase. Curved black arrows demonstrate the measurement of diameter; green arrows show the lumen of seminiferous tubules, and red arrows show height of spermatogenic epithelium. (H&E 200X).
Vetsci 09 00248 g003
Vetsci 09 00248 g004 550
Figure 4. Photomicrographs of spleen of Japanese quails during different breeding phases of annual reproductive cycle. (a) Active breeding phase, (b) Regressive breeding phase and (c) Inactive breeding phase. Red circles show the splenic nodule (present in white pulp) of spleen (H&E 200X).
Figure 4. Photomicrographs of spleen of Japanese quails during different breeding phases of annual reproductive cycle. (a) Active breeding phase, (b) Regressive breeding phase and (c) Inactive breeding phase. Red circles show the splenic nodule (present in white pulp) of spleen (H&E 200X).
Vetsci 09 00248 g004
Vetsci 09 00248 g005 550
Figure 5. Photomicrographs of bursa of Fabricius of Japanese quails during different breeding phases of annual reproductive cycle. (a) Active breeding phase, (b) Regressive breeding phase and (c) Inactive breeding phase. Red circles show the bursal nodule; black arrows show the epithelium of bursa, and red arrows show bursal capsule (H&E 200X).
Figure 5. Photomicrographs of bursa of Fabricius of Japanese quails during different breeding phases of annual reproductive cycle. (a) Active breeding phase, (b) Regressive breeding phase and (c) Inactive breeding phase. Red circles show the bursal nodule; black arrows show the epithelium of bursa, and red arrows show bursal capsule (H&E 200X).
Vetsci 09 00248 g005
Table
Table 1. Gross anatomical measurements of different organs of Japanese quails (Coturnix japonica) during different reproductive phases.
Table 1. Gross anatomical measurements of different organs of Japanese quails (Coturnix japonica) during different reproductive phases.
ParameterAnnual Reproductive Season
Active Season
(April–May)		Regressive Season
(September–October)		Nonactive Season
(January–February)
Testes
Absolute Weight (g)1.28 ± 0.16 a0.46 ± 0.12 b0.44 ± 0.14 b
Relative Weight (%)1.74 ± 0.17 a1.2 ± 0.12 b0.55 ± 0.13 c
Length (cm)2.07 ± 0.12 a1.51 ± 0.15 b0.71 ± 0.1 c
Width (cm)1.66 ± 0.13 a1.05 ± 0.02 b0.57 ± 0.07 c
Bursa of Fabricius
Absolute Weight (g)0.18 ± 0.04 b0.19 ± 0.06 b0.26 ± 0.05 a
Relative Weight (%)0.25 ± 0.07 c0.27 ± 0.035 b0.36 ± 0.04 a
Length (cm)0.98 ± 0.05 b1.13 ± 0.02 b1.45 ± 0.09 a
Width (cm)0.26 ± 0.04 c0.49 ± 0.05 b0.71 ± 0.07 a
Thymus
Absolute Weight (g)0.11 ± 0.06 c0.23 ± 0.04 b0.36 ± 0.07 a
Relative Weight (%)0.16 ± 0.01 c0.29 ± 0.03 b0.45 ± 0.08 a
Spleen
Absolute Weight (g)0.19 ± 0.03 c0.27 ± 0.05 b0.54 ± 0.04 a
Relative Weight (%)0.24 ± 0.08 b0.31 ± 0.07 b0.60 ± 0.07 a
Diameter (cm)0.78 ± 0.02 c1.1 ± 0.06 b1.65 ± 0.07 a
All values are represented as mean ± SD, whereby different superscripts (a,b,c) are used to identify a statistical significance of p < 0.05.
Table
Table 2. Histomorphometric analysis of Japanese quails (Coturnix japonica) during different reproductive phases.
Table 2. Histomorphometric analysis of Japanese quails (Coturnix japonica) during different reproductive phases.
ParameterAnnual Reproductive Season
Active Phase
(April–May)		Regressive Phase
(September–October)		Nonactive Phase
(January–February)
Testes
Diameter of seminiferous tubules (µm)203.37 ± 65.22 a151.71 ± 28.22 b91.04 ± 11.41 b
Height of spermatogenic epithelium (µm)80.66 ± 21.22 a37.06 ± 9.67 b27.80 ± 7.22 c
Lumen of seminiferous tubules (µm)122.71 ± 9.98 b114.4 ± 15.93 b63.24 ± 19.11 a
Thickness of tunica albuginea (µm)86.81 ± 17.75 b105.29 ± 22.15 b191.10 ± 26.85 a
Bursa of Fabricius
Diameter of follicles (µm)195.5 ± 58.54 c271.13 ± 69.39 b508.52 ± 91.06 a
Thickness of capsule (µm)213.46 ± 78.49 a101.42 ± 37.08 b48.51 ± 23.94 c
Epithelial height (µm)4.57 ± 1.09 a3.61 ± 0.73 a2.53 ± 0.4 a
Thymus
Mean diameter of thymic lobule (µm)105.63 ± 2.59 c156.73 ± 5.47 b199.01 ± 7.03 a
Diameter of thymic corpuscle (µm)12.11 ± 1.34 a15.16 ± 2.71 a18.79 ± 3.25 a
Spleen
Diameter of splenic nodule (µm)324.58 ± 33.34 b383.74 ± 28.86 b426.52 ± 29.47 a
Thickness of capsule (µm)46.68 ± 11.96 b57.74 ± 12.03 b82.19 ± 14.57 a
All values are represented as mean ± SD, whereby different superscripts (a,b,c) are used to identify a statistical significance of p < 0.05.
Table
Table 3. Immune and hormonal analysis of Japanese quails (Coturnix japonica) during different reproductive phases.
Table 3. Immune and hormonal analysis of Japanese quails (Coturnix japonica) during different reproductive phases.
ParameterAnnual Reproductive Season
Active Phase
(April–May)		Regressive Phase
(September–October)		Nonactive Phase
(January–February)
Immunological Analysis
Total Leukocyte Count (×103/µL)109.5 ± 1.95 a157.25 ± 2.37 b189.3 ± 6.53 c
Lymphocyte Count (×103/µL)79.34 ± 2.79 a113.94 ± 3.16 b137.18 ± 4.92 c
Hormonal Analysis
Serum Testosterone (ng/mL)2.45 ± 0.035 a0.60 ± 0.024 b0.27 ± 0.011 c
All values are represented as mean ± SD, whereby different superscripts (a,b,c) are used to identify a statistical significance of p < 0.05.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hayat, K.; Raza, A.; Anas, A.; Qureshi, A.S.; Rehan, S.; Rabbani, A.H.; Rehman, H.F.u.; Qamar, A.G.; Rehman, T.; Deeba, F.; et al. Photoperiodic Modulation in Immune and Reproductive Systems in Japanese Quails (Coturnix japonica): A Morphometric Perspective. Vet. Sci. 2022, 9, 248. https://0-doi-org.brum.beds.ac.uk/10.3390/vetsci9050248

AMA Style

Hayat K, Raza A, Anas A, Qureshi AS, Rehan S, Rabbani AH, Rehman HFu, Qamar AG, Rehman T, Deeba F, et al. Photoperiodic Modulation in Immune and Reproductive Systems in Japanese Quails (Coturnix japonica): A Morphometric Perspective. Veterinary Sciences. 2022; 9(5):248. https://0-doi-org.brum.beds.ac.uk/10.3390/vetsci9050248

Chicago/Turabian Style

Hayat, Khizar, Ali Raza, Aitzaz Anas, Anas Sarwar Qureshi, Sarmad Rehan, Ameer Hamza Rabbani, Hafiz Faseeh ur Rehman, Abdul Ghaffar Qamar, Tayyab Rehman, Farah Deeba, and et al. 2022. "Photoperiodic Modulation in Immune and Reproductive Systems in Japanese Quails (Coturnix japonica): A Morphometric Perspective" Veterinary Sciences 9, no. 5: 248. https://0-doi-org.brum.beds.ac.uk/10.3390/vetsci9050248

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop