More stories

  • in

    Carcass traits and meat quality of goats fed with cactus pear (Opuntia ficus-indica Mill) silage subjected to an intermittent water supply

    Morphometric measurements are subjective and used to assess the carcass development and quantitatively measure the muscular distribution in the carcass with estimates of its conformation. In the present study there were not significative differences observed for these parameters or for carcass compactness index (CCI), inferring that the use of cactus pear silage as well as intermittent water supply combined or alone did not alter animal growth and/or carcass conformation, maintaining the muscle pattern achieved by the control diet (usual) and demonstrating body and carcass uniformity. Since animals used in this study were homogeneous and had similar age and body performance, as indicated by the carcass morphometric measurements and by the difference between the empty carcass and hot carcass weights, which resulted in the sum of head + limb with an average of 8.2 ± 0.13 kg between treatments, giving an idea that the animals were similar in chronological age, since the allometric growth of the body occurs from the extremities to the interior of the body.The significant difference between treatments with inclusion of cactus pear silage for hot carcass yield (HCY) and cold carcass yield (CCY) may be related to the weight of the full gastrointestinal tract, which showed higher values for animals fed with a higher proportion of Tifton 85 grass hay in the diet (0% CPS). Increasing the NDF content of the diet reduces the passage rate of digesta, and the emptying of the gastrointestinal tract (GT) that cause a distension of the rumen-reticulum and increase the weight of the gastrointestinal tract, resulting in lower HCY and consequently lower CCY. While the diets with inclusion of CPS increase NFC content, such as pectin, which have higher rates of rumen degradability and, higher rates of passage7,8,9.Measurements and evaluations carried out on the carcass, such as the carcass compactness index and loin eye area (LEA), are parameters that quantitatively measure the muscle distribution in the carcass, an edible part of greater financial return, which indicates the conformation of these animals3, while the body condition score (BCS) and the measure C, which are highly correlated, measure the distribution of fat on the carcass, giving an idea of the carcass finish, in which the higher these variables, the greater the proportion of fat that allows for less water loss due to carcass cooling10. These variables in the present study were also not influenced by the levels of cactus pear silage and water restrictions, presenting an overall mean of 0.17 kg/cm, 7.68 cm, 2.42 points and 0.7 mm respectively, and consequently did not influence the losses due to cooling, which presented an average loss of 1.48%.The main cuts of the goat carcass are the neck, leg, shoulder, loin, and rib. Their economic values differ, and their proportions become an important index to evaluate the carcass quality9. The cuts of greatest importance and commercial values are the leg and the loin, called noble cuts because they present greater yield and muscle tenderness, being interesting that they present a good proportion in the carcass, for providing greater edible tissue content, mainly muscle.Carcasses with similar weight tend to have equivalent proportions of cuts, as they exhibit isogonic growth. As the cold carcass weight (CCW) and the conformation of the animals were similar, with similar morphometric measurements, they had a direct relationship in the absence of an effect on commercial cuts.The commercial value of the carcass, whether through carcass yield and/or the proportions of the cuts, is also linked to tissue composition, thus the dissection of the leg represents an estimate of measuring the tissue composition of the carcass, in which is sought a greater proportion of muscle, intermediate proportion of fat and less bone in carcasses11. In this way, diets with cactus pear silage and the different levels of intermittent water supply resulted in the constancy in the amount of muscle, fat, and bone in legs of goats. The similarity in muscle proportion is related to the lack of effects on slaughter weight and CCW, as the weight of muscles is highly correlated to carcass weight. The average muscle yield was above 60% in all treatments, confirming that the animals showed good efficiency to the diets and adapted well to the water supply levels. Although the diets with cactus silage had high amounts of metabolizable energy (ME) and no difference in DM intake, the energy input was similar that not influencing carcass weights and carcass compactness index. That is, it did not influence muscle deposition in the carcass, probably due to synchronicity of energy and protein.As for the weight and proportion of bone tissue, it is believed that because this is a tissue with early development in relation to muscle and fat2, diets in the final stages of growth (average of 8 months) would hardly change their participation in the tissue composition, where the relationship of this tissue with the others is usually only increased when there are changes in the proportion of muscle and/or fat.Water restriction, as long as it is moderate and acute, mainly affects the loss of body water and not tissues, which does not cause deleterious effects on animal productivity and growth.The muscle:fat ratio indicates the state of leg fattening, while the muscle:bone ratio estimates the carcass muscularity, both being attributes of quality3. The similarity previously reported in the weight of fat, bone and muscle corroborates that these relationships also do not have differences. The same occurs for the leg muscularity index (LMI), due to the weight of the five muscles used to determine the index and the length of the femur which had been similar between the animals.Nevertheless, when considering fat as a percentage of participation in leg weight, it is possible to observe that the intermittency in water supply in both intervals (24 and 48 h) reduced the proportion of fat in the leg. Although in this research, the water supply levels did not affect the daily intake of dry matter from animals, with average intake of 650.67 g/kg DM, ranging from 599 to 682 g/kg DM between treatments7, during days of water deprivation, fat mobilization for energy availability may occur, possibly offsetting water stress and influencing not only feed intake, on these days of deprivation but also affecting energy metabolism, which results in the mobilization of energy reserves2.When the physicochemical composition of the meat was evaluated, it was observed that the diets and water supply levels probably did not affect the reserves of muscle glycogen during the pre-slaughter management as can be seen through pHinitial and pHfinal. The pHinitial right after slaughter should be close to neutrality, as well as in the live animal, indicating that the animal did not suffer from stress during the pre-slaughter period. The pHfinal, on the other hand, is expected to show a considerable variation, between 5.55 and 6.2 for goat meat; and due be inversely proportional to the concentration of muscle glycogen at the time of slaughter, that is, a more intense expenditure of glycogen stores results in less lactic acid production and higher pHfinal10,12,13. In this research, the pHfinal had an average of 5.74, a pH higher than the isoelectric point of muscle proteins (5.2–5.3). This result is favorable, since it is above the neutral charge and presenting an excessive negative charge that provides the repulsion of filaments, which allows water molecules to bind and improve the organoleptic characteristics of the meat, through succulence and texture of meat13 evaluated by cooking loss, moisture, and shear force, principally. The cooking loss (CL), moisture and shear force (SF) were within the values recommended (20–35% CL, moisture above 70% and SF up to 44.13 Newton (N) for goat meat) to classify the meat as soft and tender14. Statistically, interactions were found between the supply of silage and intermittent water supply, in which goats on a diet without cactus pear silage and without intermittent water supply showed higher values of cooking losses and shear force.Higher concentrations of collagen content and/or greater activities of calpastatin (which inhibit the action of calpains), as well as larger fascicles and greater number of fibers present in each muscle fascicle, as was visually observed in the meat of the animals in this research, can lead to reductions in meat tenderness15. Because goat carcasses are generally small, with low marbling degree and a thin layer of subcutaneous fat, there is rapid heat dissipation at the beginning of the post-mortem period, which can lead to cold shortening, muscle hardening, and less tender meats16.pHfinal of the meat has a high correlation with color parameters (L*—lightness, a*—redness, b*—yellowness and Chroma), as the pHfinal can affect the reaction of myoglobin to oxymyoglobin. The b* index in meat, on the other hand, may be related to the concentration of fat and/or the presence of carotenoids in the diet which can be affected by forage preservation processes, such as silage and hay, which significantly reduces by up to 80% carotenoids levels13. It is believed that the carotenoid concentrations in the diet of this study were similar between treatments and consequently in values of b* of meat. Values of a* and Chroma directly depend on the content and state of the heme pigments in the muscle, due to the chemical state of iron (Fe), playing an important role in meat color10. These parameters showed no significant difference between treatments, however, higher values of a* and Chroma in meat are desired, as a result of the increase in oxymyoglobin and decrease in metmyoglobin that provides the meat’s “bloom”. According to Dawson et al.17, the minimum critical value for meat luminosity (L*) is 34. Lower values of L are related to elevating pHfinal, which results in the high concentration of metmyoglobin, making the meat darker, which causes rejection by consumers for associating dark meat to as old meat.The meat’s presentation and more precisely its color is an important factor that can influence a consumer’s purchase decision, as it gives us the idea of freshness and meat’ quality. The L* and a* color parameters are the most representative for these characteristics18. Although in our research it did not have a significant effect on the color parameters, we can indicate that the meat obtained in this research would be well accepted by consumers, because Hopkins19 suggests that consumers will consider meat color acceptable when the L* value is equal to or exceeds 34, and a* value below 19 or equal to or exceeds 9.5 according to Khliji et al.18. In the present study, all values for L* remained above this aforementioned threshold and the values of a* remained within these values which suggests that meats from all diets and water supply levels had an acceptable color for consumers.When evaluating the chemical composition of meat, no significant differences were observed between treatments, except for the ash content, that remained above the average values found in the literature, which is 0.99–1.10%16. It is believed that because cactus pear is a rich source of Ca, Mg, K and with increasing level of cactus pear silage in the diet31, these minerals were consumed in larger amounts, which could have resulted in a higher proportion of minerals in the meat of animals that received 42% cactus pear silage.The lipid fatty acid profile in meat has a major impact on sensory properties and nutritional quality, influencing acceptance and health for consumers20,21. Intermittent water supply, cactus pear silage, and interaction between water supply and cactus pear silage did not influence most fatty acids present in the Longissimus lumborum muscle of the animals under study, except only a few saturated fatty acids e.g. docosanoic acid (C22:0), tricosanoic acid (C23:0), BCFA, anteiso-tridecanoic acid (C13:0 anteiso) and anteiso-pentadecanoic acid (C15:0 anteiso).Biohydrogenation of ruminal bacteria results in a circumstantial variety of fatty acids (FA), which will be absorbed in the intestine and later incorporated into the meat of goats. In addition to the diet and the biohydrogenation, the meat lipid profile can vary due to de novo synthesis, desaturation, duration of the feeding period and differences in pathways of various FA by the animal organism22.A high concentration of saturated fatty acids present in meat is not desirable, as there is evidence that saturated fatty acids, mainly C16:0, as well as myristic (C14:0) and lauric (C12:0) increase the blood cholesterol and low-density lipoproteins (LDL) concentration, due to interferences with hepatic LDL receptors23, however, in the studied treatments, there were no significant differences for these fatty acids. On the other hand, C18:0 has no impact on cholesterol levels, due to being poorly digested and easily desaturated to C18:1 by Δ9-desaturase24, present in the cell endoplasmic reticulum. This fatty acid is not harmful to health and is considered the only desirable SFA. As the levels of C18:0 in diets tend to be minimal, their main origin is the biohydrogenation of PUFA and de novo syntheses in diets with a high energy pattern25.In addition to carrying out the biohydrogenation process, ruminal bacteria synthesize a series of FA, mainly those of odd and branched chain, that comprise mainly the lipids of the bacterial membrane26,27, to maintain membrane fluidity. Linear odd-chains fatty acids are formed when propionyl-CoA, instead of acetyl-CoA, is used as a de novo synthesis initiator25. On the other hand, iso and anteiso FA are synthesized by the precursors branched-chain amino acids (valine, leucine, and isoleucine) and their corresponding branched- short-chain carboxylic acids (isobutyric, isovaleric and 2-methyl butyric acids)28.There is an increasing interest to study odd-and branched-chain fatty acids (OBCFAs) from animal products, mainly in milk due to its higher concentration compared to meat. Researchers reported that several OBCFAs have potential health benefits in humans29 as improved gut health30 and presenting anti-cancer activity31, as well as improve the sensory characteristics of the meat, providing a greater sensation of tenderness and juiciness, because BCFA content are associated with a less consistent fat in meat from lambs due to its lower melting point and its chain structure32.The FAs profile in the ruminal bacteria is largely composed by OBCFAs (C15:0; anteiso C15:0; iso C15:0; C17:0; iso C17:0; C17:1 and anteiso C17:0) in the bacteria membrane lipids24. Thus, the higher concentration of OBCFAs might be the result of the difference in the rumen bacterial populations induced by variation in the dietary carbohydrate, that is, a higher concentration of cellulolytic bacteria in relation to amylolytic bacteria, due to the high neutral detergent fiber (NDF) content in the diet with 0% cactus forage silage. It is also known that amylolytic bacteria produce more linear odd chain and anteiso FAs than iso FAs, whereas cellulolytic bacteria produce more iso FAs28,32. As the Tifton 85 grass hay-based diet had the highest neutral detergent fiber corrected for ash and protein (NDFap) and starch content (highest % of ground corn), the meat of those animals had higher concentrations of anteiso C15:0 and anteiso C13:0 compared to animals fed diets with the inclusion of cactus pear silage, also influencing the total sum of branched chain fatty acids.Although levels of intermittent water supply have generated punctual changes in tricosanoic acid (C23:0) SFA, the same was not observed for MUFA and PUFA, due to changes in the rumen environment, promoted by water restrictions, which were not sufficient to circumstantially modify biohydrogenation, resulting in similarities in concentrations of unsaturated fatty acids in goat meat.The animals subjected to 24 h of intermittent water supply (IWS) presented the highest concentration of C23:0 in relation to other treatments, which is interesting because it is involved in the synthesis of ceramide and reduces the risk of diabetes in humans33.The cactus pear has high non-fibrous carbohydrate (NFC) content (mainly pectin), having 59.5% high and medium rumen degradation carbohydrates which provide a higher production rate and removal of short-chain fatty acids and changes in rumen bacterial populations34. The inclusion of CPS resulted in a higher passage rate of digesta, affected biohydrogenation, and resulted in the escape of intermediate fatty acids isomers that are absorbed in the small intestine. Consequently, there was changing composition of fatty acids in the muscle of these animals, with a significant effect being observed only in the cis-13 C18:1. Furthermore, diets with high proportions of cactus pear silage (CPS), such as 42% CPS diet, can decrease ruminal pH and affect the final stages of biohydrogenation, resulting in the escape of intermediate fatty acids isomers, that are absorbed in the small intestine, which can explain the similarity of the C20:1 in 42% CPS diet from the Tifton hay-based diet, with differences between goat meat from 21% CPS diet and Tifton hay-based diet.Oleic acid (c9-C18:1) was the MUFA with the highest participation in the lipid profile of goat meat, which is interesting because it has a hypocholesterolemic effect, being a desirable fatty acid (DFA) for not reducing the serum high density lipoproteins (HDL) levels and thus prevent cardiovascular disease by reducing LDL levels35. The high concentrations of c9-C18:1 in ruminant meat come from the food intake, the effect of biohydrogenation, and mainly of the high activity of Δ9-desaturase, necessary for animal biosynthesis through desaturation of C18:0 to c9-C18:127. This fatty acid in the lipid profile of red meat varies between 30 and 43%36, confirming that the meat in the present study had a good concentration of this fatty acid.Much of unsaturated fatty acids, which have 18 carbons or 16 carbons, are largely converted to C18:0 and C16:0 through biohydrogenation, and when this process is not 100% completed, in addition to the PUFA that pass through this process intact, some product intermediates are formed, reaching the duodenum and are absorbed by the animal, in which significant amounts of cis and trans-monounsaturated, such as vaccenic fatty acid (t11-C18:1), reach the duodenum and are absorbed, later composing the muscle tissue22.The literature indicates that the precursor of conjugated linoleic acid (CLA) in the meat of animals is trans vaccenic acid (t11-C18:1), so the enzyme ∆9-desaturase, besides acting in the conversion of stearic into oleic fatty acid, also converts the trans-vaccenic acid to its corresponding CLA isomer, c9t11-C18:236. This pathway is more expressive in the mammary gland, and as the concentration of vaccenic acid (t11-C18:1) was not different, the concentration of CLA was not affected by the supply of silage and intermittent water supply, in the same way, that there are also no differences in the activity of ∆9-desaturase. Nevertheless, it is worth noting that in the human adipose tissue there is also the presence of ∆9-desaturase, and therefore, increased intake of vaccenic fatty acid could have the same beneficial effects associated with the intake of CLA, where the dietary vaccenic fatty acid shows 19–30% conversion rate37.Tifton hay is a natural source of n-3 fatty acids, mainly C18:3 n-3 with up to 20% participation in the lipid profile2, allowing a certain part of these PUFAs to be absorbed and increased in the tissue muscle, with 10 to 30% PUFAs in the diet generally escaping from biohydrogenation.Linoleic fatty acid (c9c12 C18:2) and α-linolenic acid (C18:3 n-3) are essential fatty acids for humans, that serve as precursors of the n-3 and n-6 pathways, distinct families, but synthesized by some of the same enzymes (∆4-desaturase, ∆5-desaturase, and ∆6-desaturase)25. Arachidonic fatty acid (C20:4 n-6) comes from elongation and desaturation of linoleic acid, where its concentrations, even close to that of its precursor, may indicate that there was a high activity of ∆6-desaturase (desaturation to γ-linolenic), elongase (elongation of γ-linolenic to dihomo-gamma-linolenic) and ∆5-desaturase. This fatty acid was influenced by the diets, presenting lower concentrations in the meat of animals fed the 42% cactus pear silage when compared to the Tifton hay diet (0% cactus pear silage).A higher concentration of long-chain PUFA n-3, docosahexaenoic (C22:6 n-3), was observed in the muscle of animals fed on Tifton hay. This was probably due to the high concentration of C18:3 n-3, precursor of C22:6 n-3, that the hay presents in relation to the cactus pear silage.The ratios and proportions of fatty acids are used to determine nutritional and nutraceutical values of the product or diet, and mainly, to indicate the cholesterolemic potential4. It is interesting that the n-6/n-3 ratio is low due to the pro-inflammatory properties of n-6; it is recommended to decrease its intake to assist in disease prevention38, while n-3 fatty acids are anti-inflammatory, antithrombotic, antiarrhythmic and reduce blood lipids, with vasodilating properties, being interesting that they present a higher proportion24. n-6 fatty acids tend to have a higher percentage in meat, and this directly influences the formation of n-3 isomers, since linoleic acid, when in excess, can reduce the synthesis of linolenic acid metabolites. The percentage of FA in one group can interfere with the metabolism of the other, reducing its incorporation into tissue lipids and altering its general biological effects38. Therefore, it is not recommended that the n-6/n-3 ratio be kept above 5 or 639, demonstrating that the averages of the current research remained acceptable.In relation to atherogenicity index (AI) and thrombogenicity index (TI), Ulbricht and Southgate39 proposed that sheep meat should have values of up to 1.0 and 1.58, respectively, and the lower the values for these indices in the lipid fraction, the greater the prevention of early stages of cardiovascular diseases. In the present study, the general averages observed were 0.29 for the AI, and 0.81 for the TI, although there were no significant differences, all treatments are within the recommended range, despite having been used as comparative standard to sheep, due to the absence of the proposed standard for goat meat.The h:H ratio did not differ for diets and water supply levels, but had an average of 1.90, below the reference value for meat products, which is 2.0. Values above 2.0 are recommended and favorable40, as it indicates a higher proportion of hypocholesterolemic fatty acids, that are beneficial to human health.The ∆9-desaturase enzyme that acts on both the mammary gland and adipose tissue, responsible for the transformation of SFA into unsaturated fatty acids (UFA), as well as in the endogenous conversion of CLA37 did not differ between treatments. On the other hand, the elongase showed less activity. Probably there was a greater “de novo” synthesis which resulted in a greater accumulation of palmitic fatty acid, and a reduction in the activity of the elongase enzyme.The crossbred goats demonstrated to present efficient mechanisms for adapting to water restrictions, especially when receiving feed with higher water content, such as cactus pear silage, being able to replace Tifton hay with 42% cactus pear silage in the diet for goats in confinement without negatively affecting the carcass traits and meat quality. Because, although these animals have shown some differences in the indices of tenderness and juiciness of their meats, however, all presented values of juiciness and tenderness compatible with meat extremely appreciated by the consumer market, and even goat meat showing some fatty acids with different concentrations induced by the supply of silage and water intermittence, the final lipid profile was appropriate to the health of consumers, observed by the absence of differences in the total concentrations of PUFA and in the main nutraceutical parameters (DFA, n-6/n-3; h:H; AI and TI).These results are relevant, indicating that goat feedlots in regions with low water availability may adopt strategies of lesser demand for drinking water and considerable concentrations of cactus pear silage in the diet, can reduce production costs without considerably affecting the product to be marketed, and therefore, provide higher profitability of the system. More

  • in

    Genetic and ecological drivers of molt in a migratory bird

    Stefansson, S. O., Björnsson, B. T., Ebbesson, L. O. E. & McCormick, S. D. Smoltification. In Fish Larval Physiology (eds Finn, R. N. & Kapoor, B. G.) 639–681 (CRC Press, 2020).Chapter 

    Google Scholar 
    Kaleka, A. S., Kaur, N. & Bali, G. K. Larval development and molting. In Edible Insects (ed. Mikkola, H.) 17 (IntechOpen, 2019).
    Google Scholar 
    Butler, L. K. & Rohwer, V. G. Feathers and molt. in Ornithology: Foundation, Analysis, and Application (eds Morrison, M. L. et al.) 242–270 (JHU Press, 2018).
    Google Scholar 
    Swaddle, J. P., Witter, M. S., Cuthill, I. C., Budden, A. & McCowen, P. Plumage condition affects flight performance in common starlings: Implications for developmental homeostasis, abrasion and moult. J. Avian Biol. 27, 103–111 (1996).Article 

    Google Scholar 
    Norris, D. R., Marra, P. P., Montgomerie, R., Kyser, T. K. & Ratcliffe, L. M. Reproductive effort, molting latitude, and feather color in a migratory songbird. Science 306, 2249–2250 (2004).Article 
    ADS 
    CAS 

    Google Scholar 
    Delhey, K., Peters, A. & Kempenaers, B. Cosmetic coloration in birds: Occurrence, function, and evolution. Am. Nat. 169, S145–S158 (2007).Article 

    Google Scholar 
    Tomotani, B. M. & Muijres, F. T. A songbird compensates for wing molt during escape flights by reducing the molt gap and increasing angle of attack. J. Exp. Biol. 222, 195396 (2019).Article 

    Google Scholar 
    Galván, I., Negro, J. J., Rodriguez, A. & Carrascal, L. M. On showy dwarfs and sober giants: Body size as a constraint for the evolution of bird plumage colouration. Acta Ornithol. 48, 65–80 (2013).Article 

    Google Scholar 
    Speakman, J. R. & Król, E. Maximal heat dissipation capacity and hyperthermia risk: Neglected key factors in the ecology of endotherms. J. Anim. Ecol. 79, 726–746 (2010).
    Google Scholar 
    Wolf, B. O. & Walsberg, G. E. The role of the plumage in heat transfer processes of birds. Am. Zool. 40, 575–584 (2000).
    Google Scholar 
    Berthold, P. & Querner, U. Genetic basis of moult, wing length, and body weight in a migratory bird species, Sylvia atricapilla. Experientia 38, 801–802 (1982).Article 

    Google Scholar 
    Gwinner, E., Neusser, V., Engl, D., Schmidl, D. & Bals, L. Haltung, Zucht und Eiaufzucht afrikanischer und europäischer Schwarzkehlchen Saxicola torquata. Gefied. Welt 111, 118–120 (1987).
    Google Scholar 
    Berthold, P. & Querner, U. Microevolutionary aspects of bird migration based on experimental results. Isr. J. Ecol. Evol. 41, 377–385 (1995).
    Google Scholar 
    Helm, B. & Gwinner, E. Timing of postjuvenal molt in African (Saxicola torquata axillaris) and European (Saxicola torquata rubicola) stonechats: Effects of genetic and environmental factors. Auk 116, 589–603 (1999).Article 

    Google Scholar 
    Helm, B. & Gwinner, E. Timing of molt as a buffer in the avian annual cycle. Acta Zool. Sin. 52, 703–706 (2006).
    Google Scholar 
    Rohwer, S., Ricklefs, R. E., Rohwer, V. G. & Copple, M. M. Allometry of the duration of flight feather molt in birds. PLoS Biol. 7, e1000132 (2009).Article 

    Google Scholar 
    Jenni, L. & Winkler, R. The Biology of Moult in Birds (Bloomsbury Publishing, 2020).
    Google Scholar 
    Tonra, C. M. & Reudink, M. W. Expanding the traditional definition of molt-migration. Auk Ornithol. Adv. 135, 1123–1132 (2018).
    Google Scholar 
    Rohwer, S., Butler, L. K., Froehlich, D. R., Greenberg, R. & Marra, P. P. Ecology and demography of east–west differences in molt scheduling of Neotropical migrant passerines. Birds Two Worlds Ecol. Evol. Migr. (R. Greenb. PP Marra, Eds.). Johns Hopkins Univ. Press. Balt. Maryl., 87–105 (2005).Bensch, S., Åkesson, S. & Irwin, D. E. The use of AFLP to find an informative SNP: Genetic differences across a migratory divide in willow warblers. Mol. Ecol. 11, 2359–2366 (2002).Article 
    CAS 

    Google Scholar 
    Ruegg, K. Genetic, morphological, and ecological characterization of a hybrid zone that spans a migratory divide. Evol. Int. J. Org. Evol. 62, 452–466 (2008).Article 

    Google Scholar 
    Delmore, K. E., Fox, J. W. & Irwin, D. E. Dramatic intraspecific differences in migratory routes, stopover sites and wintering areas, revealed using light-level geolocators. Proc. R. Soc. B Biol. Sci. 279, 4582–4589 (2012).Article 

    Google Scholar 
    Delmore, K. E. et al. Individual variability and versatility in an eco-evolutionary model of avian migration. Proc. R. Soc. B 287, 20201339 (2020).Article 

    Google Scholar 
    Procházka, P. et al. Across a migratory divide: divergent migration directions and non-breeding grounds of Eurasian reed warblers revealed by geolocators and stable isotopes. J. Avian Biol. 49, 012516 (2018).Article 

    Google Scholar 
    Bensch, S., Grahn, M., Müller, N., Gay, L. & Åkesson, S. Genetic, morphological, and feather isotope variation of migratory willow warblers show gradual divergence in a ring. Mol. Ecol. 18, 3087–3096 (2009).Article 

    Google Scholar 
    Rohwer, S. & Irwin, D. E. Molt, orientation, and avian speciation. Auk 128, 419–425 (2011).Article 

    Google Scholar 
    Pageau, C., Sonnleitner, J., Tonra, C. M., Shaikh, M. & Reudink, M. W. Evolution of winter molting strategies in European and North American migratory passerines. Ecol. Evol. 11, 13247–13258 (2021).Article 

    Google Scholar 
    Butler, L. K., Rohwer, S. & Rogers, M. Prebasic molt and molt-related movements in Ash-throated Flycatchers. Condor 108, 647–660 (2006).Article 

    Google Scholar 
    Barry, J. H., Butler, L. K., Rohwer, S. & Rohwer, V. G. Documenting molt-migration in Western Kingbird (Tyrannus verticalis) using two measures of collecting effort. Auk 126, 260–267 (2009).Article 

    Google Scholar 
    Hobson, K. A. & Wassenaar, L. I. Linking breeding and wintering grounds of neotropical migrant songbirds using stable hydrogen isotopic analysis of feathers. Oecologia 109, 142–148 (1996).Article 
    ADS 
    CAS 

    Google Scholar 
    Hobson, K. A. & Wassenaar, L. I. Tracking Animal Migration with Stable Isotopes (Academic Press, 2018).
    Google Scholar 
    Rubenstein, D. R. & Hobson, K. A. From birds to butterflies: Animal movement patterns and stable isotopes. Trends Ecol. Evol. 19, 256–263 (2004).Article 

    Google Scholar 
    Bearhop, S. et al. Assortative mating as a mechanism for rapid evolution of a migratory divide. Science 310, 502–504 (2005).Article 
    ADS 
    CAS 

    Google Scholar 
    Eppig, J. T. et al. The mouse genome database (MGD): Comprehensive resource for genetics and genomics of the laboratory mouse. Nucleic Acids Res. 40, D881–D886 (2012).Article 
    CAS 

    Google Scholar 
    Contina, A., Bridge, E. S. & Kelly, J. F. Exploring novel candidate genes from the mouse genome informatics database: Potential implications for avian migration research. Integr. Zool. 11, 240 (2016).Article 

    Google Scholar 
    Yang, J. et al. Common SNPs explain a large proportion of the heritability for human height. Nat. Genet. 42, 565–569 (2010).Article 
    CAS 

    Google Scholar 
    Thompson, C. W. Is the Painted Bunting actually two species? Problems determining species limits between allopatric populations. Condor 93, 987–1000 (1991).Article 

    Google Scholar 
    Contina, A., Bridge, E. S., Seavy, N. E., Duckles, J. M. & Kelly, J. F. Using geologgers to investigate bimodal isotope patterns in Painted Buntings (Passerina ciris). Auk 130, 265 (2013).Article 

    Google Scholar 
    Besozzi, E., Chew, B., Allen, D. C. & Contina, A. Stable isotope analysis of an aberrant Painted Bunting (Passerina ciris) feather suggests post-molt movements. Wilson J. Ornithol. 133, 151 (2021).Article 

    Google Scholar 
    Sharp, A. et al. Spatial and Temporal Scale-Dependence of the Strength of Migratory Connectivity in a North American Passerine. https://assets.researchsquare.com/files/rs-1483049/v1/72236b63-952d-4870-89e7-461056b8625b.pdf?c=1648893558 (2022).Pyle, P. et al. Temporal, spatial, and annual variation in the occurrence of molt-migrant passerines in the Mexican monsoon region. Condor 111, 583–590 (2009).Article 

    Google Scholar 
    Bridge, E. S., Fudickar, A. M., Kelly, J. F., Contina, A. & Rohwer, S. Causes of bimodal stable isotope signatures in the feathers of a molt-migrant songbird. Can. J. Zool. 89, 951 (2011).Article 
    CAS 

    Google Scholar 
    Seutin, G., White, B. N. & Boag, P. T. Preservation of avian blood and tissue samples for DNA analyses. Can. J. Zool. 69, 82–90 (1991).Article 
    CAS 

    Google Scholar 
    Ali, O. A. et al. RAD capture rapture: Flexible and efficient sequence-based genotyping. Genetics 202, 389–400 (2016).Article 
    CAS 

    Google Scholar 
    Contina, A. et al. Characterization of SNP markers for the Painted Bunting (Passerina ciris) and their relevance in population differentiation and genome evolution studies. Conserv. Genet. Resour. 11, 5–10 (2019).Article 
    ADS 

    Google Scholar 
    Catchen, J., Hohenlohe, P. A., Bassham, S., Amores, A. & Cresko, W. A. Stacks: An analysis tool set for population genomics. Mol. Ecol. 22, 3124–3140 (2013).Article 

    Google Scholar 
    Parker, P., Li, B., Li, H. & Wang, J. The genome of Darwin’s Finch (Geospiza fortis). Gigascience 10, 100040 (2012).
    Google Scholar 
    Langmead, B. & Salzberg, S. L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 9, 357–359 (2012).Article 
    CAS 

    Google Scholar 
    Van der Auwera, G. A. et al. From FastQ data to high-confidence variant calls: The genome analysis toolkit best practices pipeline. Curr. Protoc. Bioinform. 43, 1–33 (2013).
    Google Scholar 
    McKenna, A. et al. The genome analysis toolkit: A MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 20, 1297–1303 (2010).Article 
    CAS 

    Google Scholar 
    Danecek, P. et al. The variant call format and VCFtools. Bioinformatics 27, 2156–2158 (2011).Article 
    CAS 

    Google Scholar 
    Anderson, E. genoscapeRtools: Tools for Building Migratory Bird Genoscapes (2019).Purcell, S. et al. PLINK: A tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 81, 559–575 (2007).Article 
    CAS 

    Google Scholar 
    Alexander, D. H., Novembre, J. & Lange, K. Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 19, 1655–1664 (2009).Article 
    CAS 

    Google Scholar 
    Alexander, D. H. & Lange, K. Enhancements to the ADMIXTURE algorithm for individual ancestry estimation. BMC Bioinform. 12, 246 (2011).Article 

    Google Scholar 
    Francis, R. M. pophelper: An R package and web app to analyse and visualize population structure. Mol. Ecol. Resour. 17, 27–32 (2017).Article 
    CAS 

    Google Scholar 
    Chew, B., Kelly, J. & Contina, A. Stable isotopes in avian research: a step by step protocol to feather sample preparation for stable isotope analysis of carbon (δ13C), nitrogen (δ15N), and hydrogen (δ2H). Version 1.1. https://doi.org/10.17504/protocols.io.z2uf8ew (2019).Wassenaar, L. I. & Hobson, K. A. Comparative equilibration and online technique for determination of non-exchangeable hydrogen of keratins for use in animal migration studies. Isotopes Environ. Health Stud. 39(3), 211–217 (2003).Article 
    CAS 

    Google Scholar 
    Bowen, G. J., Wassenaar, L. I. & Hobson, K. A. Global application of stable hydrogen and oxygen isotopes to wildlife forensics. Oecologia 143, 337–348 (2005).Article 
    ADS 

    Google Scholar 
    R Core Team: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, 2021).Wassenaar, L. I. & Hobson, K. A. Stable-hydrogen isotope heterogeneity in keratinous materials: Mass spectrometry and migratory wildlife tissue subsampling strategies. Rapid Commun. Mass Spectrom. 20, 2505–2510 (2006).Article 
    ADS 
    CAS 

    Google Scholar 
    Zhou, X., Carbonetto, P. & Stephens, M. Polygenic modeling with Bayesian sparse linear mixed models. PLoS Genet. 9, e1003264 (2013).Article 
    CAS 

    Google Scholar 
    Guan, Y. & Stephens, M. Bayesian variable selection regression for genome-wide association studies and other large-scale problems. Ann. Appl. Stat. 5, 455 (2011).Article 
    MathSciNet 
    MATH 

    Google Scholar 
    Marchini, J., Cardon, L. R., Phillips, M. S. & Donnelly, P. The effects of human population structure on large genetic association studies. Nat. Genet. 36, 512–517 (2004).Article 
    CAS 

    Google Scholar 
    Browning, S. R. & Browning, B. L. Rapid and accurate haplotype phasing and missing-data inference for whole-genome association studies by use of localized haplotype clustering. Am. J. Hum. Genet. 81, 1084–1097 (2007).Article 
    CAS 

    Google Scholar 
    Chaves, J. A. et al. Genomic variation at the tips of the adaptive radiation of Darwin’s finches. Mol. Ecol. 25, 5282–5295 (2016).Article 
    CAS 

    Google Scholar 
    Zhang, Y.-W. et al. mrMLM v4.0.2: An R platform for multi-locus genome-wide association studies. Genom. Proteom. Bioinform. 18, 481–487 (2020).Article 

    Google Scholar 
    Grabherr, M. G. et al. Genome-wide synteny through highly sensitive sequence alignment: Satsuma. Bioinformatics 26, 1145–1151 (2010).Article 
    CAS 

    Google Scholar 
    Quinlan, A. R. & Hall, I. M. BEDTools: A flexible suite of utilities for comparing genomic features. Bioinformatics 26, 841–842 (2010).Article 
    CAS 

    Google Scholar 
    Ellis, N., Smith, S. J. & Pitcher, C. R. Gradient forests: Calculating importance gradients on physical predictors. Ecology 93, 156–168 (2012).Article 

    Google Scholar 
    Hijmans, R. J., Cameron, S. E., Parra, J. L., Jones, P. G. & Jarvis, A. Very high resolution interpolated climate surfaces for global land areas. Int. J. Climatol. 25, 1965–1978 (2005).Article 

    Google Scholar 
    Anderson, E. C. snps2assays: Prepare SNP Assay Orders from ddRAD or RAD Loci (2015).Li, H. & Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25, 1754–1760 (2009).Article 
    CAS 

    Google Scholar 
    Patterson, N., Price, A. L. & Reich, D. Population structure and eigenanalysis. PLoS Genet. 2, e190 (2006).Article 

    Google Scholar 
    Price, A. L. et al. Principal components analysis corrects for stratification in genome-wide association studies. Nat. Genet. 38, 904–909 (2006).Article 
    CAS 

    Google Scholar 
    Ruegg, K. et al. Ecological genomics predicts climate vulnerability in an endangered southwestern songbird. Ecol. Lett. 21, 1085–1096 (2018).Article 

    Google Scholar 
    Bay, R. A. et al. Genomic signals of selection predict climate-driven population declines in a migratory bird. Science 359, 83–86 (2018).Article 
    ADS 
    CAS 

    Google Scholar 
    Hedenström, A. Adaptations to migration in birds: Behavioural strategies, morphology and scaling effects. Philos. Trans. R. Soc. B Biol. Sci. 363, 287–299 (2008).Article 

    Google Scholar 
    Buehler, D. M. & Piersma, T. Travelling on a budget: Predictions and ecological evidence for bottlenecks in the annual cycle of long-distance migrants. Philos. Trans. R. Soc. B Biol. Sci. 363, 247–266 (2008).Article 

    Google Scholar 
    Schieltz, P. C. & Murphy, M. E. The contribution of insulation changes to the energy cost of avian molt. Can. J. Zool. 75, 396–400 (1997).Article 

    Google Scholar 
    Carling, M. D. & Thomassen, H. A. The role of environmental heterogeneity in maintaining reproductive isolation between hybridizing Passerina (Aves: Cardinalidae) buntings. Int. J. Ecol. 2012, 1–11 (2012).Article 

    Google Scholar 
    Irwin, D. E. Incipient ring speciation revealed by a migratory divide. Mol. Ecol. 18, 2923–2925 (2009).Article 

    Google Scholar 
    Thomas, D. W., Blondel, J., Perret, P., Lambrechts, M. M. & Speakman, J. R. Energetic and fitness costs of mismatching resource supply and demand in seasonally breeding birds. Science 291, 2598–2600 (2001).Article 
    ADS 
    CAS 

    Google Scholar 
    Rohwer, V. G., Rohwer, S. & Ortiz-Ramirez, M. F. Molt biology of resident and migrant birds of the monsoon region of west Mexico. Ornitol. Neotrop. 20, 565–584 (2009).
    Google Scholar 
    Bensch, S., Andersson, T. & Åkesson, S. Morphological and molecular variation across a migratory divide in willow warblers, Phylloscopus trochilus. Evolution 53, 1925–1935 (1999).Article 

    Google Scholar 
    Turbek, S. P., Scordato, E. S. C. & Safran, R. J. The role of seasonal migration in population divergence and reproductive isolation. Trends Ecol. Evol. 33, 164–175 (2018).Article 

    Google Scholar 
    Scordato, E. S. C. et al. Migratory divides coincide with reproductive barriers across replicated avian hybrid zones above the Tibetan Plateau. Ecol. Lett. 23, 231–241 (2020).Article 

    Google Scholar 
    Battey, C. J. et al. A migratory divide in the Painted Bunting (Passerina ciris). Am. Nat. 191, 259–268 (2018).Article 
    CAS 

    Google Scholar 
    Contina, A. et al. Genetic structure of the Painted Bunting and its implications for conservation of migratory populations. Ibis 161, 372 (2019).Article 

    Google Scholar 
    Butler, L. K. The grass is always greener: Do monsoon rains matter for molt of the Vermilion Flycatcher (Pyrocephalus rubinus)? Auk 130, 297–307 (2013).Article 

    Google Scholar 
    Turbek, S. P. et al. A migratory divide spanning two continents is associated with genomic and ecological divergence. Evolution 76, 722 (2022).Article 

    Google Scholar 
    Dietz, M. W., Daan, S. & Masman, D. Energy requirements for molt in the kestrel Falco tinnunculus. Physiol. Zool. 65, 1217–1235 (1992).Article 

    Google Scholar 
    Vézina, F., Gustowska, A., Jalvingh, K. M., Chastel, O. & Piersma, T. Hormonal correlates and thermoregulatory consequences of molting on metabolic rate in a northerly wintering shorebird. Physiol. Biochem. Zool. 82, 129–142 (2009).Article 

    Google Scholar 
    Bazzi, G. et al. Candidate genes have sex-specific effects on timing of spring migration and moult speed in a long-distance migratory bird. Curr. Zool. 63, 479–486 (2017).CAS 

    Google Scholar 
    Busby, L. et al. Sonic hedgehog specifies flight feather positional information in avian wings. Development 147, 188821 (2020).Article 

    Google Scholar 
    Eichberger, T. et al. GLI2-specific transcriptional activation of the bone morphogenetic protein/Activin antagonist Follistatin in human epidermal cells. J. Biol. Chem. 283, 12426–12437 (2008).Article 
    CAS 

    Google Scholar 
    Matzuk, M. M. et al. Multiple defects and perinatal death in mice deficient in follistatin. Nature 374, 360–363 (1995).Article 
    ADS 
    CAS 

    Google Scholar 
    Patel, K., Makarenkova, H. & Jung, H.-S. The role of long range, local and direct signalling molecules during chick feather bud development involving the BMPs, follistatin and the Eph receptor tyrosine kinase Eph-A4. Mech. Dev. 86, 51–62 (1999).Article 
    CAS 

    Google Scholar 
    Nakamura, M. et al. Control of pelage hair follicle development and cycling by complex interactions between follistatin and activin. FASEB J. 17, 1–22 (2003).Article 
    MathSciNet 

    Google Scholar 
    Pays, L., Charvet, I., Hemming, F. J. & Saxod, R. Close link between cutaneous nerve pattern development and feather morphogenesis demonstrated by experimental production of neo-apteria and ectopic feathers: Implication of chondroitin sulphate proteoglycans and other matrix molecules. Anat. Embryol. 195, 457–466 (1997).Article 
    CAS 

    Google Scholar 
    Pyle, P., Saracco, J. F. & DeSante, D. F. Evidence of widespread movements from breeding to molting grounds by North American landbirds. Auk Ornithol. Adv. 135, 506–520 (2018).
    Google Scholar 
    De Mita, S. et al. Detecting selection along environmental gradients: Analysis of eight methods and their effectiveness for outbreeding and selfing populations. Mol. Ecol. 22, 1383–1399 (2013).Article 

    Google Scholar 
    Lotterhos, K. E. & Whitlock, M. C. Evaluation of demographic history and neutral parameterization on the performance of FST outlier tests. Mol. Ecol. 23, 2178–2192 (2014).Article 

    Google Scholar 
    Frichot, E., Schoville, S. D., de Villemereuil, P., Gaggiotti, O. E. & François, O. Detecting adaptive evolution based on association with ecological gradients: Orientation matters!. Heredity (Edinb.) 115, 22–28 (2015).Article 
    CAS 

    Google Scholar 
    Trivedi, A. K. et al. Temperature alters the hypothalamic transcription of photoperiod responsive genes in induction of seasonal response in migratory redheaded buntings. Mol. Cell. Endocrinol. 493, 110454 (2019).Article 
    CAS 

    Google Scholar  More

  • in

    Tracking microbes in extreme environments

    In 2008, I was investigating the methane bubbling up on the beaches and in shallow waters of Mocha Island, off the coast of central Chile. I became intrigued by how microorganisms could thrive in methane-rich areas and changed my research focus from marine biology to extreme environments. I wanted to understand how methane acts as a source of energy and carbon for microbes.Since then, I have explored a number of bizarre environments. In 2010, I went in a submarine down to 200 metres in the Black Sea, one of the world’s largest anoxic water bodies. There, I found mats of filamentous bacteria that survive on sulfur compounds.In 2017, I studied the microbes in Canada’s tailing ponds, artificial lakes of water, sand and clay waste that are left behind after petroleum extraction. And I sampled the microorganisms living in 100 °C Antarctic hot springs in 2022.I came home to Chile in 2018 and began collaborating with an international team researching the geomicrobiology of thermal features, including hot springs, geysers and volcanoes. After travelling with the group to Argentina’s active volcanic region, I got funding to explore the microbial communities that exist beneath hydrothermal vents in southern Chile, where the oceanic crust is subducting beneath the continental plate.In this image, I am in the Atacama Desert in South America, the driest non-polar desert on the planet. I am measuring 80–100 °C steam released from a fumarole containing yellow sulfur, which crystallizes at its opening as the vapour cools. I also sampled sub-surface microbes that are flushed out with the fluids. We’ll sequence their DNA to assess the microbial communities and their biological interactions.My goal is to learn more about subsurface microbes in extreme environments. I want to understand how microbial forces shaped the planet and how these communities might shift in the future with climate change. More

  • in

    Reply to: Plant traits alone are good predictors of ecosystem properties when used carefully

    Plant Ecology and Nature Conservation Group, Wageningen University, Wageningen, the NetherlandsFons van der Plas & Liesje MommerSystematic Botany and Functional Biodiversity, Life Science, Leipzig University, Leipzig, GermanyThomas Schröder-Georgi, Alexandra Weigelt, Kathryn Barry & Christian WirthGerman Centre for Integrative Biodiversity Research Halle-Jena-Leipzig, Leipzig, GermanyAlexandra Weigelt, Kathryn Barry, Adriana Alzate, Nico Eisenhauer, Anke Hildebrandt, Christiane Roscher & Christian WirthTerrestrial Ecology Research Group, School of Life Sciences Weihenstephan, Technical University of Munich, Munich, GermanySebastian Meyer & Wolfgang WeisserAquaculture and Fisheries Group, Wageningen University and Research Centre, Wageningen, the NetherlandsAdriana AlzateAgroécologie, AgroSup Dijon, Institut National de la Recherche Agronomique, Université de Bourgogne, Université de Bourgogne Franche-Comté, Dijon, FranceRomain L. BarnardEidgenössische Technische Hochschule Zürich, Zurich, SwitzerlandNina BuchmannDepartment of Experimental Plant Ecology, Institute for Water and Wetland Research, Radboud University Nijmegen, Nijmegen, the NetherlandsHans de KroonInstitute of Ecology and Evolution, University Jena, Jena, GermanyAnne Ebeling & Winfried VoigtInstitute of Biology, Leipzig University, Leipzig, GermanyNico EisenhauerHumboldt-Universität zu Berlin, Berlin, GermanyChristof EngelsInstitute of Plant Sciences, University of Bern, Bern, SwitzerlandMarkus FischerMax Planck Institute for Biogeochemistry, Jena, GermanyGerd Gleixner, Ernst-Detlef Schulze & Christian WirthHelmholtz Centre for Environmental Research, Leipzig, GermanyAnke HildebrandtFriedrich Schiller University Jena, Jena, GermanyAnke HildebrandtGeoecology, University of Tübingen, Tübingen, GermanyEva Koller-France & Yvonne OelmannInstitute of Geography and Geoecology, Karlsruhe Institute of Technology, Karlsruhe, GermanySophia Leimer & Wolfgang WilckeEcotron Européen de Montpellier, Centre National de la Recherche Scientifique, Montferrier-sur-Lez, FranceAlexandru MilcuCentre d’Ecologie Fonctionnelle et Evolutive, Unité Mixte de Recherche 5175 (Centre National de la Recherche Scientifique-Université de Montpellier-Université Paul-Valéry Montpellier-Ecole Pratique des Hautes Etudes), Montpellier, FranceAlexandru MilcuDepartment of Evolutionary Biology and Environmental Studies, University of Zurich, Zurich, SwitzerlandPascal A. NiklausUFZ, Helmholtz Centre for Environmental Research, Department Physiological Diversity, Leipzig, GermanyChristiane RoscherInstitute of Landscape Ecology, University of Münster, Münster, GermanyChristoph ScherberCentre for Biodiversity Monitoring, Zoological Research Museum Alexander Koenig, Bonn, GermanyChristoph ScherberGeobotany, Faculty of Biology, University of Freiburg, Freiburg, GermanyMichael Scherer-LorenzenCentre of Biodiversity and Sustainable Land Use, University of Göttingen, Göttingen, GermanyStefan ScheuJ.F. Blumenbach Institute of Zoology and Anthropology, Animal Ecology, University of Göttingen, Göttingen, GermanyStefan ScheuDepartment of Geography, University of Zurich, Zurich, SwitzerlandBernhard SchmidInstitute of Ecology, College of Urban and Environmental Sciences, Peking University, Beijing, ChinaBernhard SchmidLeuphana University Lüneburg, Institute of Ecology, Lüneburg, GermanyVicky TempertonAgroecology, Department of Crop Sciences, University of Göttingen, Göttingen, GermanyTeja TscharntkeF.v.d.P. wrote the initial draft of the manuscript. T.S.-G., A.W., K.B., S.M., A.A., R.L.B., N.B., H.d.K., A.E., N.E., C.E., M.F., G.G., A.H., E.K.-F., S.L., A.M., L.M., P.A.N., Y.O., C.R., C.S., M.S.-L., S.S., B.S., E.-D.S., V.T., T.T., W.V., W. Weisser, W. Wilcke and C.W. helped edit the manuscript. More

  • in

    Modelling the impact of non-pharmaceutical interventions on the spread of COVID-19 in Saudi Arabia

    Here, we examine infection dynamics in the four regions of focus to learn more about how various control interventions performed in each region. Since the first documented cases emerged in these regions, the virus was able to spread freely across much of the first and second phases with a gradual increase in the control interventions. In the four regions of Makkah, Madinah, Eastern, and Riyadh, cases peaked on 12th May (6397 cases; 95% CI 5960–9697), 15th May (1967 cases; 95% CI 1625–2308), 23rd June (10367; 95% CI 8948–11785) and 11th June (11273 cases; 95% CI 11068–12491) respectively according to the fitted model shown in Fig. 2. As the epidemic progressed, more measures were adopted to contain the disease, and the disease’s infectiousness sharply decreased after the third period. There are a few factors responsible for the sudden declining trend: first our model is dependent on official data on cases that have been documented, and these data will only ever reflect a portion of the overall number of cases. Second, different regions developed different testing strategies, and some locations altered their approach to testing during the course of the time period that was investigated. It is possible that the beginning of the Hajj term (period 4) was a contributing factor in the decrease in the number of documented cases. Additionally, previous to this time period, the government indicated that it would be increasing the size of its local testing in order to detect new cases. It is possible that the efficacy of interventions would be reduced if increases are found to be occurring during the falling phase of an epidemic. This may result in measures being kept in place for a longer period of time than they would have been had more accurate data been provided.Figure 2With the use of the Delay Rejection Adaptive Metropolis method, the relevant parameters were estimated for each of the four areas of interest by fitting the data from 13th March until 25th September.Full size imageWe estimate the effective reproduction number (R_t) as an indicator of SARS-CoV-2 transmission before and after the interventions. Figure 3 depicts the dramatic shift in the rate of SARS-CoV-2 transmission as a result of decreased social contact and other control measures. At the beginning of the pandemic, (R_t) for SARS-CoV-2 in Saudi regions was between 4 and 6 as illustrated in Tables 6 and 7. In other words, on average each case spread to between four and six others. Considering that each new generation of SARS-CoV-2 cases occurs every five days, it is evident that this pandemic was rapidly expanding out of control. Moreover, we assumed that the transmission rate and the documented infection rate did not change during the first two periods since interventions were carried out gradually until a complete lockdown took place. As more measures were introduced, the spread of the disease began to decrease. Therefore, our data were based on the weekly reported number of documented SARS-CoV-2 cases broken down by region. As a result, it became clear that the reliability of the (R_t) value was relatively high for transmission.Figure 3Distribution of Rt estimates derived from 10000 MCMC samples for Makkah, Madinah, Eastern, and Riyadh, respectively. The black dot in the centre of each violin plot denotes the median, the thick bar in the plot denotes the interquartile range, and the thin bar in the plot denotes the lowest and maximum values. The mean and the credible interval for 95%, which is shown in parentheses, are labelled below or above, respectively.Full size imageThe effects of the events and interventions on the dynamics of SARS-CoV-2 in the regions of interest are considered. First, if the controls remained in phase four in Makkah, our model projects that the total number of documented cases would increase to 81047 (95% CI 79421–82672). In Al-Madinah, the cumulative number of documented cases would have increased to 22997 (95% CI 19578–26415). The number of cumulative documented cases may have reached 80520 (95% CI 78335–82704) if controls stayed steady in the Eastern region at the level they were at in phase four. If the pattern shown during the fourth period is taken into account, we estimate that there would have been 67150 (95% CI 63731–70568) documented infections in the Riyadh region. Figure 4 illustrates these findings.Figure 4The relevant parameters were estimated for each of the four regions of interest by first fitting the data of each region, and then predicting using the parameters from period 4. This was done with each of the four regions of interest separately.Full size imageWe now explore the impact of controls remaining in place at the same level as that implemented in phase three. In that case, the number of documented cases in Makkah would have increased to 116641 (95% CI 105015–128266). Similarly, the total number of documented cases in Al-Madinah would have increased to 53877 (95% CI 50458-57295) if the outbreak had been allowed to continue at the same level. If the controls had remained unchanged from how they were in phase three in the Eastern region, the total number of documented cases would have been 310459 (95% CI 298362–334981). Finally, in Riyadh this would have resulted in 665241 documented cases (95% CI 651822 to 678659). Figure 5 highlights these findings.Figure 5For each of the four areas of interest, the relevant parameters were estimated by first fitting the data of each region and then predicting using the parameters from period 3. This was carried out for each of the four areas of interest independently.Full size imageWe now investigate the impact of second-period controls remaining in place. In that case, the number of documented cases would increase to 1236642 (95% 1218314–1251626), 442865 (95% CI 439446–456283), 454031 (95% CI 441846–466215), and 2322624 (95% CI 1919206-3026042) in the regions of Makkah, Madinah, Eastern, and Riyadh, respectively (see Fig. 6).Figure 6For each of the four regions of interest, the relevant variables were identified by first fitting the data from each area and then making predictions using the parameters from period 2. Each of the four regions of interest was done separately.Full size imageThe efficacy of NPIs is dependent on when they are adopted, with earlier adoption resulting in greater success in lowering transmission rates of infectious diseases. In the early stages of COVID-19, Saudi regions made the decision to gradually implement measures in order to understand the severity of the disease and reduce the economic and social costs of lockdowns, as well as the political costs. In Fig. 6, if the government were to rely on the interventions of the second phase, then the number of cases of infection would considerably rise owing to the ineffectiveness of the measures. In the third period as in Fig. 5, the government made it possible to relax some of the control measures, but it is ultimately up to each area to decide whether they will maintain the same level of control or whether they will increase or decrease it. In comparison to the control measures carried out during the third and fourth periods, this led to significantly improved outcomes. The reason that these time periods were chosen is that there was no stiffening of the NPI response in most Saudi regions during the first two periods and control interventions were improved later on.Significant undetected infections resulted in the fast spread of new coronaviruses (SARS-CoV-2) which is illustrated in Fig. 7. The proportion of undocumented infections, including asymptomatic cases and undocumented symptomatic individuals who did not seek medical treatment or be tested for mild symptoms, was greater than that of Wuhan at the onset of the pandemic26, which may be a result of the following factors: first, the medical configuration was not optimal and public awareness was limited during the onset of the pandemic while the undocumented rate progressively increased; Second, contact tracing procedures employed in Saudi regions may have become overwhelmed if the number of early-stage cases in Saudi regions rises substantially. The discrepancy between the predicted proportions of asymptomatic (undocumented) cases may be attributable to the difficulty in the un-identifiability of parameters in epidemiological models. There were a substantial number of asymptomatic infected individuals with high infectivity in Saudi regions, where the epidemic situation escalated rapidly. Our research emphasises the frequency of asymptomatic SARS-CoV-2 cases and their role in transmission in order to increase people’s knowledge of asymptomatic cases and to serve as a guide for the prevention and control of SARS-CoV-2.Table 3 Estimated transmission rate in Saudi Regions.Full size tableTable 4 Estimated ascertainable infection rate.Full size tableIn this model, we fitted dynamic transmission rates because of varied preventable measures by the Saudi government at the level of the country or region. After a series of actions taken by the government, regions and cities went into lockdown, resulting in a decrease in the transmission rate as in Table 3. Before the interventions were introduced, in the first two periods of our study, we assumed the transmission rate did not change since individual and community responses had not effectively taken place. After severe interventions were implemented, the transmission rates were allowed to vary in later periods and reduced gradually due to the control measures that reduced the spread of disease27. Estimates of documented infection rates are presented in Table 4. Our model estimates show the documented infection rate has continued to decrease in the last two periods. Thus, the parameters we fit across periods are a measure of how effective the lockdown was in bringing down the documented infection rate28.Risk of resurgenceThe risk of resurgence in Saudi Arabia’s four regions has been examined in this section after the relaxation of intervention measures. There will be a rise in disease activity if control measures are relaxed without taking into account increases in the number of cases being detected, isolated, and/or traced. We predict the first week of no new cases of infection and the week when all current infections in Saudi Arabia will be eradicated.In the Makkah region, had the trend continued into the fourth period, the number of documented infections would have dropped to zero on average by the 6th September (23rd August to 27th September), and all infections would have been eradicated by the 26th of October (7th October to 14th November). On the 28th June, the number of weekly active infections (including presymptomatic, symptomatic, and asymptomatic cases) reached its highest point of 230,230 (95% CI 226811–234364), and on 8th September, that number dropped to 44023 (95% CI 40604–47441).Therefore, the number of documented infections would have reached zero in Al-Madinah region on average on 6th November (23rd October to 22nd November), and all infections would have been eliminated by 1st December (27th November to 14th December). On 23rd June, weekly active infections (including presymptomatic, symptomatic, and asymptomatic cases) peaked at 130,134 (95% CI 126715–133552) and then declined to 60023 (95% CI 58604–63441) on 25th September.If the trend had continued as it did in the fourth period in the Eastern region, the average number of documented infections would have reached zero on 2nd November (from 23rd October to 18th November), and the total eradication of infections would have happened on 1st December (26th November to 22nd December). The number of weekly active infections (including presymptomatic, symptomatic, and asymptomatic cases) peaked at 65000 (95% CI 61581–68418) during the week of July 23rd and subsequently decreased to 800 (95% CI 765–834) on 8th of September.Lastly, the model predicted that the number of weekly active infections in the Riyadh region (including presymptomatic, symptomatic, and asymptomatic infections) peaked on 28th June at 562332 (95% CI 513379–619542) and then decreased to 188215 (95% CI 174796–191633) on 18th September. On average, we expected that the number of documented infections would have decreased to zero on 18th October (7th October to 14th November) and that the total number of infections would have been eliminated on 1st December if the trend continued as it did in the fourth period (20th November to 23rd December). Figure 7 illustrates these findings. We found that if control measures were lifted 30 days following the first day of zero documented cases.Figure 7The estimated number of infected cases that were active (presymptomatic, symptomatic, and asymptomatic) during the research period in the areas of Makkah, Madinah, Eastern, and Riyadh respectively.Full size imageThe probability of resurgence, which we define as the number of active documented cases greater than 100 could be as high as 0.96 in Eastern, 0.95 in Madinah, 0.97 in Makkah, and 0.96 in Riyadh. If we adopt more stringent conditions of lifting controls after observing no confirmed cases for a continuous period of 30 days, the probability of resurgence decreases to 0.31, 0.28, 0.30, and 0.30, with probable resurgence occurring on 13th February, 7th February, 2nd January, and 8th January for Eastern, Makkah, Madinah and Riyadh, respectively (Fig. 8). Despite the use of a simplified model, these results emphasize the hazards of ignoring undetermined occurrences when modifying intervention techniques.Figure 8Figure demonstrating the effect of relaxing all control measures in all four regions 30 days following the first day without confirmed cases.Full size imageSensitivity analysisFor the purpose of testing the robustness of our research results, we conducted a series of sensitivity analyses by varying the durations of the latent and infectious periods, the ratio of transmissibility in asymptomatic (undocumented) cases to symptomatic (documented) cases, and the initial documented infection rate. We conduct eight sensitivity analyses (S1 to S8) within each model for each region of Saudi Arabia to assess the robustness of our model results. For instance, the sensitivity analysis performed for S1 was based on the changes of the latent period and pre-symptomatic infectious period, respectively, and other parameters remain the same. These modifications were carried out with the help of reference15,29, and the same approaches were used for the other parts of the sensitivity analysis, which is summarised in Table 5.Table 5 Description of essential model parameters that were not fitted in the MCMC, where (D_e) refers to the latent period, (D_p) refers to the pre-symptomatic infectious period, (D_i) refers to the symptomatic infectious period, (gamma _0) refers to the initial ascertain rate and (alpha) refers to the ratio of the transmission rate for P and A to I.Full size tableIn particular, for (S1), we raised the incubation period to 7 days (upper 95% CI based on ref15) and the pre-symptomatic infectious period to 3 days (upper 95% CI based on ref29). Therefore we set (D_e = 4) and (D_p=3), and modified (E_0) and (P_0) as needed. The transmissibility of the undocumented cases was assumed to be 0.46 (lower 95 % CI according to ref.31) of the infection cases for (S2); for (S3), the transmissibility of the asymptomatic (undocumented) cases was assumed to be 0.62 (upper 95 % CI according to ref31). We assumed that in (S4), the initial documented infection rate was (gamma _0) = 0.14 (lower 95 % CI according to ref13) and adjusted (A_0), (P_0) and E(0) accordingly. Similarly for (S5) we assumed the initial documented infection rate was (gamma _0) = 0.42 (upper 95 % CI according to ref13) and adjusted (P_0), (A_0), and (E_0) accordingly. In (S6) we set the variables (D_ e=3) and (D_p=1.1), and altered the values of (P_0) and (E_ 0) as necessary in accordance with13. In (S7) we assumed that the transmission rate of asymptomatic (undocumented) cases was half that of documented cases by setting 0.5. Finally, in (S8) we assumed that the infectious period ((D_i)) was double that of symptomatic cases by setting 6 days. Both (S7) and (S8) were based on30. The results of our sensitivity analysis are summarised in Tables 6 and 7. We note that the variation in the model predictions of (R_t) varies from setting to setting. However, these variations appear to be fairly small, proposing the robustness of the results to the specification of associated values in fairly realistic ranges13,13. Our sensitivity analysis provides information about the importance of each parameter to the model representing the transmission of SARS-CoV-2. An increase (or decrease) in parameter values, while other parameters’ values remain the same, contributes to an increase (or decrease) in effective reproduction numbers. For example, an increase in infectious period would result in a higher effective reproduction number at the beginning of the epidemic and a longer time required to clear all infections in Saudi regions32. Our sensitivity analysis indicates that almost all model parameters may have an important role in spreading this virus among susceptible people. In particular, the contact rate from person-to-person and the transition rate of asymptomatic (undetected) individuals play a significant role in disease spread. Our important findings, of a significant decrease in (R_t) after interventions and the existence of a substantial number of presymptomatic and asymptomatic cases, were found to be robust. This highlights that Saudi authorities should pay attention to intervention strategies in the event of a resurgence of cases and quarantining those who were in contact with active cases can effectively reduce the disease33. In Tables 6 and 7 we show the estimated effective reproduction number (R_t) associated with 95% CIs obtained from those eight sensitivity analyses for all four regions and all five time periods.Table 6 Sensitivity analysis of the effective reproduction number (R_t) for Eastern and Madinah.Full size tableTable 7 Sensitivity analysis of estimated effective reproduction number (R_t) for Makkah and Riyadh.Full size table More

  • in

    Landscape management strategies for multifunctionality and social equity

    The Global Assessment Report on Biodiversity and Ecosystem Services: Summary for Policy-Makers (IPBES, 2019)DeFries, R. & Nagendra, H. Ecosystem management as a wicked problem. Science 356, 265–270 (2017).Article 
    CAS 

    Google Scholar 
    Turkelboom, F. et al. When we cannot have it all: ecosystem services trade-offs in the context of spatial planning. Ecosyst. Serv. 29, 566–578 (2018).Article 

    Google Scholar 
    Lee, H. & Lautenbach, S. A quantitative review of relationships between ecosystem services. Ecol. Indic. 66, 340–351 (2016).Article 

    Google Scholar 
    Bennett, E. M., Peterson, G. D. & Gordon, L. J. Understanding relationships among multiple ecosystem services. Ecol. Lett. 12, 1394–1404 (2009).Article 

    Google Scholar 
    Goldstein, J. H. et al. Integrating ecosystem-service tradeoffs into land-use decisions. Proc. Natl Acad. Sci. USA 109, 7565–7570 (2012).Article 
    CAS 

    Google Scholar 
    Vallet, A., Locatelli, B. & Pramova, E. Ecosystem Services and Social Equity: Who Controls, Who Benefits and Who Loses? (CIFOR, 2020); https://doi.org/10.17528/cifor/007849Neyret, M. et al. Assessing the impact of grassland management on landscape multifunctionality. Ecosyst. Serv. 52, 101366 (2021).Linders, T. E. W. et al. Stakeholder priorities determine the impact of an alien tree invasion on ecosystem multifunctionality. People Nat. 3, 658–672 (2021).Article 

    Google Scholar 
    Herzig, A., Ausseil, A.-G. & Dymond, J. in Ecosystem Services in New Zealand—Conditions and Trends (ed. Dymond, J. R.) 511–523 (Manaaki Whenua Press, 2014).Chan, K. M. A., Shaw, M. R., Cameron, D. R., Underwood, E. C. & Daily, G. C. Conservation planning for ecosystem services. PLoS Biol. 4, e379 (2006).Article 

    Google Scholar 
    Pennington, D. N. et al. Cost-effective land use planning: optimizing land use and land management patterns to maximize social benefits. Ecol. Econ. 139, 75–90 (2017).Article 

    Google Scholar 
    Hölting, L. et al. Including stakeholders’ perspectives on ecosystem services in multifunctionality assessments. Ecosyst. People 16, 354–368 (2020).Article 

    Google Scholar 
    Plieninger, T. et al. Exploring futures of ecosystem services in cultural landscapes through participatory scenario development in the Swabian Alb, Germany. Ecol. Soc. 18, 39 (2013).Article 

    Google Scholar 
    Tasser, E., Schirpke, U., Zoderer, B. M. & Tappeiner, U. Towards an integrative assessment of land-use type values from the perspective of ecosystem services. Ecosyst. Serv. 42, 101082 (2020).Article 

    Google Scholar 
    Sayer, J. et al. Ten principles for a landscape approach to reconciling agriculture, conservation, and other competing land uses. Proc. Natl Acad. Sci. USA 110, 8349–8356 (2013).Article 
    CAS 

    Google Scholar 
    Vallet, A. et al. Linking equity, power, and stakeholders: roles in relation to ecosystem services. Ecol. Soc. 24, 14 (2019).Article 

    Google Scholar 
    Allan, E. et al. Land use intensification alters ecosystem multifunctionality via loss of biodiversity and changes to functional composition. Ecol. Lett. 18, 834–843 (2015).Article 

    Google Scholar 
    Hector, A. & Bagchi, R. Biodiversity and ecosystem multifunctionality. Nature 448, 188–190 (2007).Article 
    CAS 

    Google Scholar 
    Manning, P. et al. Redefining ecosystem multifunctionality. Nat. Ecol. Evol. 2, 427–436 (2018).Article 

    Google Scholar 
    Raudsepp-Hearne, C., Peterson, G. D. & Bennett, E. M. Ecosystem service bundles for analyzing tradeoffs in diverse landscapes. Proc. Natl Acad. Sci. USA 107, 5242–5247 (2010).Article 
    CAS 

    Google Scholar 
    Daniel, T. C. et al. Contributions of cultural services to the ecosystem services agenda. Proc. Natl Acad. Sci. USA 109, 8812–8819 (2012).Article 
    CAS 

    Google Scholar 
    Gunton, R. M. et al. Beyond ecosystem services: valuing the invaluable. Trends Ecol. Evol. 32, 249–257 (2017).Article 

    Google Scholar 
    Peter, S., Le Provost, G., Mehring, M., Müller, T. & Manning, P. Cultural worldviews consistently explain bundles of ecosystem service prioritisation across rural Germany. People Nat. 4, 218–230 (2022).Article 

    Google Scholar 
    Haines-Young, R. & Potschin, M. in Ecosystem Ecology (eds Raffaelli, D. G. & Frid, C. L. J.) 110–139 (Cambridge Univ. Press, 2010).Fischer, M. et al. Implementing large-scale and long-term functional biodiversity research: the Biodiversity Exploratories. Basic Appl. Ecol. 11, 473–485 (2010).Article 

    Google Scholar 
    Wilson, E. O. Half-Earth: Our Planet’s Fight for Life (Norton, 2017).Bastin, J.-F. et al. The global tree restoration potential. Science 365, 76–79 (2019).Article 
    CAS 

    Google Scholar 
    Clapp, J. & Moseley, W. G. This food crisis is different: COVID-19 and the fragility of the neoliberal food security order. J. Peasant Stud. 47, 1393–1417 (2020).Article 

    Google Scholar 
    Kirwan, J. & Maye, D. Food security framings within the UK and the integration of local food systems. J. Rural Stud. 29, 91–100 (2013).Article 

    Google Scholar 
    Ellis, E. C. To conserve nature in the Anthropocene, half Earth is not nearly enough. One Earth 1, 163–167 (2019).Article 

    Google Scholar 
    Boetzl, F. A. et al. A multitaxa assessment of the effectiveness of agri-environmental schemes for biodiversity management. Proc. Natl Acad. Sci. USA 118, e2016038118 (2021).Tyllianakis, E. & Martin-Ortega, J. Agri-environmental schemes for biodiversity and environmental protection: how we are not yet ‘hitting the right keys’. Land Use Policy 109, 105620 (2021).Article 

    Google Scholar 
    Arroyo-Rodríguez, V. et al. Designing optimal human-modified landscapes for forest biodiversity conservation. Ecol. Lett. 23, 1404–1420 (2020).Article 

    Google Scholar 
    Gilroy, J. J. et al. Cheap carbon and biodiversity co-benefits from forest regeneration in a hotspot of endemism. Nat. Clim. Change 4, 503–507 (2014).Article 

    Google Scholar 
    Lindenmayer, D. B. et al. Avoiding bio-perversity from carbon sequestration solutions: avoiding bio-perversity in carbon markets. Conserv. Lett. 5, 28–36 (2012).Article 

    Google Scholar 
    Stoll-Kleemann, S. & O’Riordan, T. in The Encyclopedia of the Anthropocene Vol. 3 (eds DellaSala, D. A. & Goldstein, M. I.) 347–353 (Elsevier, 2018).Schaich, H., Bieling, C. & Plieninger, T. Linking ecosystem services with cultural landscape research. GAIA 19, 269–277 (2010).Article 

    Google Scholar 
    O’Connor, L. M. J. et al. Balancing conservation priorities for nature and for people in Europe. Science 372, 856–860 (2021).Article 

    Google Scholar 
    Büscher, B. et al. Half-Earth or Whole Earth? Radical ideas for conservation, and their implications. Oryx 51, 407–410 (2017).Article 

    Google Scholar 
    van der Plas, F. et al. Towards the development of general rules describing landscape heterogeneity–multifunctionality relationships. J. Appl. Ecol. 56, 168–179 (2019).Article 

    Google Scholar 
    Almeida, I., Rösch, C. & Saha, S. Converting monospecific into mixed forests: stakeholders’ views on ecosystem services in the Black Forest Region. Ecol. Soc. 26, 28 (2021).Meyer, M. A. & Früh-Müller, A. Patterns and drivers of recent agricultural land-use change in southern Germany. Land Use Policy 99, 104959 (2020).Article 

    Google Scholar 
    Kastner, T. et al. Global agricultural trade and land system sustainability: implications for ecosystem carbon storage, biodiversity, and human nutrition. One Earth 4, 1425–1443 (2021).Rasmussen, L. V. et al. Social–ecological outcomes of agricultural intensification. Nat. Sustain. 1, 275–282 (2018).Article 

    Google Scholar 
    Lindborg, R. et al. How spatial scale shapes the generation and management of multiple ecosystem services. Ecosphere 8, e01741 (2017).Article 

    Google Scholar 
    Duarte, G. T., Santos, P. M., Cornelissen, T. G., Ribeiro, M. C. & Paglia, A. P. The effects of landscape patterns on ecosystem services: meta-analyses of landscape services. Landsc. Ecol. 33, 1247–1257 (2018).Article 

    Google Scholar 
    Le Provost, G. et al. The supply of multiple ecosystem services requires biodiversity across spatial scales. Nat. Ecol. Evol. https://doi.org/10.1038/s41559-022-01918-5 (2022).Martin, D. A. et al. Land-use trajectories for sustainable land system transformations: identifying leverage points in a global biodiversity hotspot. Proc. Natl Acad. Sci. USA 119, e2107747119 (2022).Article 
    CAS 

    Google Scholar 
    Seabloom, E. W., Borer, E. T. & Tilman, D. Grassland ecosystem recovery after soil disturbance depends on nutrient supply rate. Ecol. Lett. 23, 1756–1765 (2020).Article 

    Google Scholar 
    Messinger, J. & Winterbottom, B. African forest landscape restoration initiative (AFR100): restoring 100 million hectares of degraded and deforested land in Africa. Nat. Faune 30, 14–17 (2016).
    Google Scholar 
    Whittingham, M. J. The future of agri-environment schemes: biodiversity gains and ecosystem service delivery? J. Appl. Ecol. 48, 509–513 (2011).Article 

    Google Scholar 
    Le Clec’h, S. et al. Assessment of spatial variability of multiple ecosystem services in grasslands of different intensities. J. Environ. Manage. 251, 109372 (2019).Article 

    Google Scholar 
    Forschungsethische Grundsätze und Prüfverfahren in den Sozial‐ und Wirtschaftswissenschaften Output 9, Berufungsperiode 5 (German Data Forum, 2017).Strukturdaten Reutlingen—Statistisches Bundesamt (Bundeswahlleiter, 2020); https://www.bundeswahlleiter.de/europawahlen/2019/strukturdaten/bund-99/land-8/kreis-8415.htmlStrukturdaten Uckermark—Statistisches Bundesamt (Bundeswahlleiter, 2020); https://www.bundeswahlleiter.de/europawahlen/2019/strukturdaten/bund-99/land-12/kreis-12073.htmlStrukturdaten Unstrut-Hainich-Kreis—Statistisches Bundesamt (Bundeswahlleiter, 2020); https://www.bundeswahlleiter.de/europawahlen/2019/strukturdaten/bund-99/land-16/kreis-16064.htmlBlüthgen, N. et al. A quantitative index of land-use intensity in grasslands: integrating mowing, grazing and fertilization. Basic Appl. Ecol. 13, 207–220 (2012).Article 

    Google Scholar 
    Ostrowski, A., Lorenzen, K., Petzold, E. & Schindler, S. Land use intensity index (LUI) calculation tool of the Biodiversity Exploratories project for grassland survey data from three different regions in Germany since 2006, BEXIS 2 module. Zenodo https://doi.org/10.5281/zenodo.3865579 (2020).Schall, P. et al. The impact of even‐aged and uneven‐aged forest management on regional biodiversity of multiple taxa in European beech forests. J. Ecol. 55, 267–278 (2018).Statistisches Jahrbuch über Ernährung, Landwirtschaft und Forsten der Bundesrepublik Deutschland Vol. 63 (Bundesministerium für Ernährung und Landwirtschaft, 2019).Simons, N. K. & Weisser, W. W. Agricultural intensification without biodiversity loss is possible in grassland landscapes. Nat. Ecol. Evol. https://doi.org/10.1038/s41559-017-0227-2 (2017).Zinke, O. Heupreise steigen: Futter für die Bauern knapp und teuer. Agrarheute https://www.agrarheute.com/markt/futtermittel/heupreise-steigen-futter-fuer-bauern-knapp-teuer-571946 (2020).Bois de Chez Nous (Lignum, 2021); https://www.lignum.ch/files/images/Downloads_francais/Shop/20010_Bois_de_chez_nous.pdfGerman Timber Company—Internationaler Holzhandel (German Timber Company, 2021); https://www.germantimber.company/. Accessed 2021-11-24Holzeinschlag nach Holzartengruppen, Holzsorten, ausgewählten Besitzarten (Statistisches Bundesamt, 2022); https://www.destatis.de/DE/Themen/Branchen-Unternehmen/Landwirtschaft-Forstwirtschaft-Fischerei/Wald-Holz/Tabellen/holzeinschlag-deutschland.htmlJahresjagdstrecke Bundesrepublik Deutschland, 2019–2020 (Deutsche Jagdverband, 2020); https://www.jagdverband.de/sites/default/files/2021-01/2021-01_Infografik_Jahresjagdstrecke_Bundesrepublik_Deutschland_2019_2020.jpgHeinze, E. et al. Habitat use of large ungulates in northeastern Germany in relation to forest management. For. Ecol. Manage. 261, 288–296 (2011).Article 

    Google Scholar 
    Conant, R. T., Cerri, C. E. P., Osborne, B. B. & Paustian, K. Grassland management impacts on soil carbon stocks: a new synthesis. Ecol. Appl. 27, 662–668 (2017).Article 

    Google Scholar 
    Hermes, J., Albert, C. & von Haaren, C. Mapping and Assessing Local Recreation as a Cultural Ecosystem Service in Germany. UVP-Report https://doi.org/10.17442/uvp-report.034.08 (2020).Hermes, J., Albert, C. & von Haaren, C. Assessing the aesthetic quality of landscapes in Germany. Ecosyst. Serv. 31, 296–307 (2018).Article 

    Google Scholar 
    Ehrhart, S. & Schraml, U. Perception and evaluation of natural forest dynamics. Allg. Forst Jagdztg. 185, 166–183 (2014).
    Google Scholar 
    Villanueva-Rivera, L. J. & Pijanowski, B. C. soundecology: Soundscape ecology. R package version 1.3.3 (2018).Meyer, S., Wesche, K., Krause, B. & Leuschner, C. Dramatic losses of specialist arable plants in central Germany since the 1950s/60s—a cross-regional analysis. Divers. Distrib. 19, 1175–1187 (2013).Article 

    Google Scholar 
    Sasaki, K., Hotes, S., Kadoya, T., Yoshioka, A. & Wolters, V. Landscape associations of farmland bird diversity in Germany and Japan. Glob. Ecol. Conserv. 21, e00891 (2020).Article 

    Google Scholar 
    Peña, L., Casado-Arzuaga, I. & Onaindia, M. Mapping recreation supply and demand using an ecological and a social evaluation approach. Ecosyst. Serv. 13, 108–118 (2015).Article 

    Google Scholar 
    Schägner, J. P., Brander, L., Paracchini, M.-L., Hartje, V. & Maes, J. Mapping recreational ecosystem services and its values across Europe: a combination of GIS and meta-analysis. In European Association of Environmental and Resource Economists 22nd Annual Conference (2016).R Core Team. R: A Language and Environment for Statistical Computing v.4.2.1 (R Foundation for Statistical Computing, 2022).Rust Programming Language https://www.rust-lang.org/ v 1.44Le Provost, G. et al. Contrasting responses of above- and belowground diversity to multiple components of land-use intensity. Nat. Commun. 12, 3918 (2021).Article 

    Google Scholar 
    Gini, C. On the measurement of concentration and variability of characters (English translation from Italian by Fulvio de Santis in 2005). Metron 63, 1–38 (1914). More

  • in

    Plastic plumage colouration in response to experimental humidity supports Gloger’s rule

    West-Eberhard, M. J. Developmental Plasticity and Evolution (Oxford University Press, 2003).Book 

    Google Scholar 
    Piersma, T. & Van Gils, J. A. The Flexible Phenotype: A Body-Centred Integration of Ecology, Physiology, and Behaviour (Oxford University Press, 2011).
    Google Scholar 
    Piersma, T. & Drent, J. Phenotypic flexibility and the evolution of organismal design. Trends Ecol. Evol. 18, 228–233 (2003).Article 

    Google Scholar 
    Tabari, H. Climate change impact on flood and extreme precipitation increases with water availability. Sci. Rep. 10, 1–10 (2020).
    Google Scholar 
    Beck, H. E. et al. Present and future Köppen–Geiger climate classification maps at 1-km resolution. Sci. Data 5, 180214 (2018).Article 

    Google Scholar 
    Bogert, C. M. Thermoregulation in reptiles, a factor in evolution. Evolution 3, 195–211 (1949).Article 
    CAS 

    Google Scholar 
    Rensch, B. Das Prinzip geographischer Rassenkreise und das Problem der Artbildung (Gebrueder Borntraeger, 1929).
    Google Scholar 
    Clusella Trullas, S., van Wyk, J. H. & Spotila, J. R. Thermal melanism in ectotherms. J. Therm. Biol. 32, 235–245 (2007).Article 

    Google Scholar 
    Delhey, K. A review of Gloger’s rule, an ecogeographical rule of colour: Definitions, interpretations and evidence. Biol. Rev. 94, 1294–1316 (2019).
    Google Scholar 
    Stuart-Fox, D., Newton, E. & Clusella-Trullas, S. Thermal consequences of colour and near-infrared reflectance. Philos. Trans. R. Soc. B: Biol. Sci. 372, 20160345 (2017).Article 

    Google Scholar 
    Friedman, N. R. & Remês, V. Ecogeographical gradients in plumage coloration among Australasian songbird clades. Glob. Ecol. Biogeogr. 26, 261–274 (2017).Article 

    Google Scholar 
    Delhey, K. Darker where cold and wet: Australian birds follow their own version of Gloger’s rule. Ecography 41, 673–683 (2018).Article 

    Google Scholar 
    Galván, I., Rodríguez-Martínez, S. & Carrascal, L. M. Dark pigmentation limits thermal niche position in birds. Funct. Ecol. 32, 1531–1540 (2018).Article 

    Google Scholar 
    Medina, I. et al. Reflection of near-infrared light confers thermal protection in birds. Nat. Commun 9, 3610 (2018).Article 
    ADS 

    Google Scholar 
    Aldrich, J. W. & James, F. C. Ecogeographic variation in the American Robin (Turdus migratorius). Auk 108, 230–249 (1991).
    Google Scholar 
    Morales, H. E. et al. Neutral and selective drivers of colour evolution in a widespread Australian passerine. J. Biogeogr. 44, 522–536 (2017).Article 

    Google Scholar 
    Griffith, S. C., Owens, I. P. & Burke, T. Environmental determination of a sexually selected trait. Nature 400, 358–360 (1999).Article 
    ADS 
    CAS 

    Google Scholar 
    Fargallo, J. A., Laaksonen, T., Korpimäki, E. & Wakamatsu, K. A melanin-based trait reflects environmental growth conditions of nestling male Eurasian kestrels. Evol. Ecol. 21, 157–171 (2007).Article 

    Google Scholar 
    Fargallo, J. A., Martínez, F., Wakamatsu, K., Serrano, D. & Blanco, G. Sex-dependent expression and fitness consequences of sunlight derived color phenotypes. Am. Nat. 191, 726–743 (2018).Article 

    Google Scholar 
    Beebe, W. Geographic variation in birds, with especial reference to the effects of humidity. Zoologica 1, 3–41 (1907).
    Google Scholar 
    Bieber, H. Fellverdunklung beim hauskaninchen nach kälteeinwirkung. Zeitschrift für Säugetierkunde 38, 33–38 (1972).
    Google Scholar 
    Johnston, R. F. & Selander, R. K. House sparrows: Rapid evolution of races in North America. Science 144, 548–550 (1964).Article 
    ADS 
    CAS 

    Google Scholar 
    Galván, I., Wakamatsu, K. & Alonso-Álvarez, C. Black bib size is associated with feather content of pheomelanin in male house sparrows. Pigment Cell Melanoma Res. 27, 1159–1161 (2014).Article 

    Google Scholar 
    Endler, J. A. On the measurement and classification of colour in studies of animal colour patterns. Biol. J. Linn. Soc. 41, 315–352 (1990).Article 

    Google Scholar 
    Montgomerie, R. Analyzing colors. In Bird Colouration I. Mechanisms and Measurements (eds Hill, E. G. & McGraw, K. J.) (Harvard University Press, 2006).
    Google Scholar 
    McGraw, K. J., Dale, J. & Mackillop, E. A. Social environment during molt and the expression of melanin-based plumage pigmentation in male house sparrows (Passer domesticus). Behav. Ecol. Sociobiol. 53, 116–122 (2003).Article 

    Google Scholar 
    Lessells, C. M. & Boag, P. T. Unrepeatable repeatabilities a common mistake. Auk 104, 116–121 (1987).Article 

    Google Scholar 
    Anderson, T. R. Biology of the Ubiquitous House Sparrow (Oxford University Press, 2006).Book 

    Google Scholar 
    Gelman, A. & Hill, J. Data Analysis Using Regression and Multilevel/Hierarchical Models (Cambridge University Press, 2006).Book 

    Google Scholar 
    Nakagawa, S., Ockendon, N., Gillespie, D. O., Hatchwell, B. J. & Burke, T. Assessing the function of house sparrows’ bib size using a flexible meta-analysis method. Behav. Ecol. 18, 831–840 (2007).Article 

    Google Scholar 
    Hill, G. E. & McGraw, K. J. Bird Coloration, Volume I: Mechanisms and Measurements (Harvard University Press, 2006).Book 

    Google Scholar 
    D’Alba, L. & Shawkey, M. D. Melanosomes: Biogenesis, properties, and evolution of an ancient organelle. Physiol. Rev. 99, 1–19 (2018).Article 

    Google Scholar 
    Delhey, K., Burger, C., Fiedler, W. & Peters, A. Seasonal changes in colour: A comparison of structural, melanin- and carotenoid-based plumage colours. PLoS ONE 5, e11582 (2010).Article 
    ADS 

    Google Scholar 
    Galván, I., Mousseau, T. A. & Møller, A. P. Bird population declines due to radiation exposure at Chernobyl are stronger in species with pheomelanin-based coloration. Oecologia 165, 827–835 (2011).Article 
    ADS 

    Google Scholar 
    Meunier, J., Pinto, S. F., Burri, R. & Roulin, A. Eumelanin-based coloration and fitness parameters in birds: A meta-analysis. Behav. Ecol. Sociobiol. 65, 559–567 (2011).Article 

    Google Scholar 
    Roulin, A., Almasi, B., Meichtry-Stier, K. S. & Jenni, L. Eumelanin- and pheomelanin-based colour advertise resistance to oxidative stress in opposite ways. J. Evol. Biol. 24, 2241–2247 (2011).Article 
    CAS 

    Google Scholar 
    Gasparini, J. et al. Strength and cost of an induced immune response are associated with a heritable melanin-based colour trait in female tawny owls. J. Anim. Ecol. 78, 608–616 (2009).Article 

    Google Scholar 
    Fargallo, J. A. et al. Sex-specific phenotypic integration: Endocrine profiles, coloration, and behavior in fledgling boobies. Behav. Ecol. 25, 76–87 (2013).Article 

    Google Scholar 
    Wittkopp, P. J. & Beldade, P. Development and evolution of insect pigmentation: Genetic mechanisms and the potential consequences of pleiotropy. Semin. Cell Dev. Biol. 20, 65–71 (2009).Article 
    CAS 

    Google Scholar 
    Hubbard, J. K., Uy, J. A. C., Hauber, M. E., Hoekstra, H. E. & Safran, R. J. Vertebrate pigmentation: From underlying genes to adaptive function. Trends Genet. 26, 231–239 (2010).Article 
    CAS 

    Google Scholar 
    McKinnon, J. S. & Pierotti, M. E. Colour polymorphism and correlated characters: Genetic mechanisms and evolution. Mol. Ecol. 19, 5101–5125 (2010).Article 

    Google Scholar 
    Poston, J. P., Hasselquist, D., Stewart, I. R. & Westneat, D. F. Dietary amino acids influence plumage traits and immune responses of male house sparrows, Passer domesticus, but not as expected. Anim. Behav. 70, 1171–1181 (2005).Article 

    Google Scholar 
    McGraw, K. J. Dietary mineral content influences the expression of melanin-based ornamental coloration. Behav. Ecol. 18, 137–142 (2007).Article 

    Google Scholar 
    Fargallo, J. A., Martínez-Padilla, J., Toledano-Díaz, A., Santiago-Moreno, J. & Dávila, J. A. Sex and testosterone effects on growth, immunity and melanin coloration of nestling Eurasian kestrels. J. Anim. Ecol. 76, 201–209 (2007).Article 

    Google Scholar 
    Fitze, P. S. & Richner, H. Differential effects of a parasite on ornamental structures based on melanins and carotenoids. Behav. Ecol. 13, 401–407 (2002).Article 

    Google Scholar 
    Roulin, A., Altwegg, R., Jensen, H., Steinsland, I. & Schaub, M. Sex-dependent selection on an autosomal melanic female ornament promotes the evolution of sex ratio bias. Ecol. Lett. 13, 616–626 (2010).Article 

    Google Scholar 
    Sharma, A. Effect of ambient humidity on UV/visible photodegradation of melanin thin films. Photochem. Photobiol. 86, 852–855 (2010).Article 
    CAS 

    Google Scholar 
    Burtt, E. H. The adaptiveness of animal colors. Bioscience 31, 723–729 (1981).Article 

    Google Scholar 
    Heppner, F. The metabolic significance of differential absorption of radiant energy by black and white birds. Condor 72, 50–59 (1970).Article 

    Google Scholar 
    Clusella-Trullas, S., Terblanche, J. S., Blackburn, T. M. & Chown, S. L. Testing the thermal melanism hypothesis: A macrophysiological approach. Funct. Ecol. 22, 232–238 (2008).Article 

    Google Scholar 
    Zink, R. M. & Remsen, J. V. Evolutionary processes and patterns of geographic variation in birds. Curr. Ornithol. 4, 1–69 (1986).
    Google Scholar 
    Burtt, E. H. & Ichida, J. M. Gloger’s rule, feather-degrading bacteria, and color variation among song sparrows. Condor 106, 681–686 (2004).Article 

    Google Scholar 
    Ruiz-De-Castaneda, R., Burtt, E. H. Jr., Gonzalez-Braojos, S. & Moreno, J. Bacterial degradability of an intrafeather unmelanized ornament: A role for feather-degrading bacteria in sexual selection?. Biol. J. Linn. Soc. 105, 409–419 (2012).Article 

    Google Scholar 
    Goldstein, G. et al. Bacterial degradation of black and white feathers. Auk 121, 656–659 (2004).Article 

    Google Scholar 
    Ducrest, A. L., Keller, L. & Roulin, A. Pleiotropy in the melanocortin system, coloration and behavioural syndromes. Trends Ecol. Evol. 23, 502–510 (2008).Article 

    Google Scholar 
    Kim, S. Y., Fargallo, J. A., Vergara, P. & Martínez-Padilla, J. Multivariate heredity of melanin-based coloration, body mass and immunity. Heredity 111, 139–146 (2013).Article 
    CAS 

    Google Scholar 
    Horrocks, N. P. C. et al. Environmental proxies of antigen exposure explain variation in immune investment better than indices of pace of life. Oecologia 177, 281–290 (2015).Article 
    ADS 

    Google Scholar 
    McLean, N., Van Der Jeugd, H. P. & van de Pol, M. High intra-specific variation in avian body condition responses to climate limits generalisation across species. PLoS ONE 13, e0192401 (2018).Article 

    Google Scholar 
    Gardner, J. L. et al. Spatial variation in avian bill size is associated with humidity in summer among Australian passerines. Clim. Change Responses 3, 11 (2016).Article 

    Google Scholar 
    Gerson, A. R. et al. Flight at low ambient humidity increases protein catabolism in migratory birds. Science 333, 1434–1436 (2011).Article 
    ADS 
    CAS 

    Google Scholar  More