More stories

  • in

    Analysis on ecological characteristics of Mississippian coral reefs in Langping, Guangxi

    Notwithstanding constraints on the amount of hard data, according to our integrated analysis, the developmental environment and ecology of reef communities have an important impact on the appearance of reefs.Analysis of environmental conditions for reef developmentSettings of reef developmentThe F/F extinction event in Late Devonian caused the complete recession of the reef-building communities based on stromatoporoid-coral assemblages7,17. The Carboniferous is generally considered to be a sub-optimal period for the development of framed reefs. After the biological mass extinction, microorganisms and algae rebuilt new reef-building ecosystems18,19. Some short-term biological frame reefs developed with low diversity, limited reef-building organisms, small sizes, and restricted distribution20. Harsh climate and marine conditions occurred in the Mississippian, including extensive marine hypoxia, repeated glacial and interglacial climate changes, and frequent changes of sea level and seawater surface temperature, potentially hindering the recovery of Early Carboniferous metazoan reefs7,21.Metazoans gradually began to participate in reef building in Early Viséan. A large number of biogenic structures formed by corals and bryozoans began to appear, including a small number of sponge reefs/mounds in the middle and late stage of Viséan. The richness and biodiversity of the Mississippian post-zoobenthic reefs flourished in the late Viséan during which corals, bryozoans, sponges, calcareous microorganisms, and some calcareous algae became the main builders3 and large-scale reefs could also be seen in some areas although most of the Viséan metazoan reefs were tabular or laminar. Thus, the metazoan skeletal reefs in the middle to late Viséan were considered to have been resurrected due to relatively warm climatic conditions and higher sea levels after a period of complete disintegration at the end of the Devonian and recession at the beginning of the Carboniferous7.Consequently, the coral reefs in the study area were the products of shallow benthic communities thriving in relatively favourable conditions of Late Viséan-Serpukhovian, which was common for reef development at that time7. Thus, it is expected that more synchronous reefs would to be identified in southern China, or even in the study area in the future.Paleogeography of reefsLangping is located in Dian-Qian-Gui Basin22 regionally (Fig. 2), in the eastern end of Tethys tectonic domain and at the interjunction of Tethys and Pacific structure globally. The Carboniferous Dian-Qian-Gui Basin is adjacent to the Tethys Basin. During the Early Carboniferous, the continent of Gondwana was close to the equator but was separated from the northern continent by the Tethys, where the tropical currents flowed freely from east to west. The benthic warm-water organisms were distributed widely with high abundance and diversity on both sides of the shallow shelf of Tethys.Figure 2Paleogeographic map of southern China in Viséan-Serpukhovian (modified from Feng23, Yao8, and Maillet24). This figure was obtained from articles by Feng23, Yao8 and Maillet24 respectively. The author modified the picture with CorelDRAW (version 2022, and the URL link: https://www.coreldraw.com/cn/). QG Qian-Gui Basin, DQGX Dian-Qian-Gui-Xiang platform.Full size imageViséan-Serpukhovian ecosystems experienced dramatic climate changes and widespread glaciation25. However, the Viséan was also a key layer for a variety of biological structures, with abundant coral reefs and a high diversity of shallow benthic communities, peaking in the late Viséan. Newly discovered post-faunal reefs in Tianlin were mainly formed in the late Viséan-Serpukhovian period, which coincided with frequent sea level fluctuations and possible glacial changes. It seems counterintuitive that tropical coral communities developed during glacial period. However, recent studies suggest that the persistent warm ocean currents on the platform helped some coral species survive from Carboniferous glacial events24. While other areas of symbiotic reefs were poorly developed, Tianlin may provide an ecological sanctuary for corals associated with ocean currents26.Sedimentation of reef developmentAccording to the regional geological structure, the slope model for Langping paleocarbonate platform was obviously different from that of steep slope platform margin, which could be directly affected by waves. Langping palaeo-platform could be regarded as one of the small blocks (block fault barrier) separated from a large platform (continental margin sea basin)13. The relative positions of these blocks were crucial for the emergence and growth of reefs.In situ development of mud-crystalline tuffs and muddy tuffs with weak hydrodynamic conditions is common in the Langping area, and evaporites are poorly developed. There were patch reefs and reef layers in different sizes in the wide intraclast beach, where obviously developed reef beach complexes were rare. The fragments of carbonate base broken by storm in the clastic beach haven’t been observed. The study area is considered to be gentle-slope open platform27,28 based on sedimentary characteristics. It suggests that the study area was far away from the margin of steep-slope platform that directly affected by waves, and more consistent with less energetic internal environments of gentle-slope platform.On the vast platform of Langping gentle slop, deep water lead to low water energy. While in the coastal area, the water energy was relatively strong, thus coarse-grained bioclastic beach and a small amount of point reef could be developed. The beaches were irregular-shaped due to long term transportation and reformation effects of waves and water flow, showing low and gentle slope angles. Dispersed reef-beach complexes at the platform margin slightly impacted inner-platform seawater and the water flows smoothly29.Therefore, it can be assumed that Langping reefs developed along the intertidal shallows of the terrace. The seawater around Langping carbonate platform in Late Viséan-Serpukhovian was relatively shallow while the water flow was strong. Remains of crinoids, brachiopods, a few foraminifera, and solitary corals were likely broken by strong currents, and deposited in situ with a small amount of gravels and lime-mud (Fig. 3). The clastic beach was unstable, suggesting large-scale wave-resistant structures could not be formed quickly30 due to insufficient cohesive and consolidating organisms. In addition, the circumferential impact of water in extensive terraces leads to mud-lime deposition, which is detrimental to most benthic organisms. However, bondstone was more likely to be formed by some binding algaes in the platform (Fig. 4). Therefore, neither the surrounding or the inner region of the platform could provide favourable living conditions for coral reefs to develop over for a long term. The gently sloping terrace environment of Lanping resulted in significant differences in growth size, wave resistance and reef-building capacity between corals in the study area and those on the edge of the steeply sloping terrace.Figure 3Clastic beaches in the Langping. Various clastic beaches developed in the study area. Diverse composition, fragmentation degree and sorting of the clast indicate different water conditions of formation. This figure is modified by the author from field photos with CorelDRAW (version 2022, and the URL link: https://www.coreldraw.com/cn/).Full size imageFigure 4Algal bondstone in the Langping. Bondstone formed by various algaes living in still water. Morphology of bondstone correlates water environment and deposition of mud. Vast algal bondstones indicate deep water and high deposition rate. This figure is modified by the author from field photos with CorelDRAW (version 2022, and the URL link: https://www.coreldraw.com/cn/).Full size imageAt the same time, the warm climate of the late Viséan-Serpukhovian, the good circulation of seawater around the Langping platform, and the abundant supply of oxygen and nutrients were a series of favourable conditions that facilitated the growth of reef-building corals, which led to uplifts being formed on clastic beach, including patch reefs and reef layers with certain sizes. These uplifts impeded waves and provided a protected nearshore environment, though they were much smaller than those developed at the steep-slope platform margin. The inhabitants on the beach could not resist strong waves. These rises were therefore known as reef-beach complexes and could only persist where waves and currents were mild28. They were essentially different from the framework coral reefs which developed on steep-slope platform margin that reflected changed hydrodynamic conditions, nutrient sources, reef sizes, and growth rates.Another potentially favourable factor in the study area could be the deeper water area in the gentle-slope sedimentary environment, which could provide more stable conditions and reduce the damaging effects of global glacial events and large scale sea level fluctuations on reef-building communities25. The frequent fault activities in Dian-Qian-Gui Basin caused the rise and fall of equivalent sea level. More influence of sea-level fluctuations and hydrodynamic conditions would be exerted on Langping platform due to its small size. Furthermore, reef growth promoted by reef-building communities would be frequently disturbed. The sediments displaying evidence of multicycle sedimentation, different components, and diversely fragile clasts in the study area provided direct evidence of frequently changing environment.Alternatively, the sedimentary environment of Langping platform provided conditions favorable for reef-building communities to develop and reefs to grow rapidly. These factors directly or indirectly determined the ecology of reef-building communities and the general appearance of reef development in the study area.Overall, the environmental factor is the primary factor affecting the overall development trend of reefs.Inferred ecological characteristics of reef communitiesResponse of reef-building corals to hydrodynamic conditionsHydrodynamic conditions are very important factors for reef development, which directly determine the abundance and distribution of each reef-building population and are key factors influencing sedimentation and reef growth, and was particularly evident at Langping. Evidence from the fossils suggested the reef-building corals were also changed in response (Table 1). The hydrodynamic condition changes during the development of reefs are inferred based on analysis of the vertical sediments and microfacies changes of coral reefs in the study area31. How these ancient reef-building corals adapted to hydrodynamic conditions was reconstructed combining the evolution of reef-building communities with the study.Table.1 General situation of reef-building coral population in Langping.Full size tableThe Xiadong coral reef started with colonization and expansion of Diphyphyllum on the bioclastic hard substrates32,33. They grew vertically into upright clusters (Fig. 5A) and were insensitive to more sediment in a relatively calm, turbid water environment34. The relatively dense clumped Siphonodendron and massive Lithostrotion (Fig. 5B) were better suited to the turbulent water environment, becoming dominant over time, with Diphyphyllum subordinate with the continuous increase of the water energy, as indicated by the characteristics of sediment particles from fine to coarse. After flourishing for a period, the Siphonodendron–Lithostrotion assemblage eventually waned, likely due to the failure to adapt to the increasing hydrodynamic conditions. Diphyphyllum had persisted combined with Syringopora, to maintain the growth of the reef. However, this assemblage subsequently declined as a result of strong hydrodynamic conditions and finally died out in response to continuous falling of sea level. Consequently, the reefs stopped developing.Figure 5Sketch of coral cluster with upright growing morphology. Most reef-building corals in Langping grow vertically into cluster colonies. This type of morphology is very favourable for corals to get more living space and is important to reef-building. (A) Cluster coral individuals grow uprightly with certain distance between each other. (B) Polygonal columnar coral individuals grow closely to resist strong water flow. This figure is made by the author with CorelDRAW (version 2022, and the URL link: https://www.coreldraw.com/cn/).Full size imageThe Longjiangdong multi-layer reef was composed of three relatively independent, flat reef layers, suggesting three distinct periods of reef development. Diverse species were identified in the reef, with colonial coral Diphyphyllum contributing greatly to reef growth. Diphyphyllum clusters colonized in patchy form on substrates composed of bioclasts or lithic gravels (Fig. 6A). The first reef-building process was brief, ending under high-energy water conditions after a period of growing (Fig. 6B). Subsequently the hydrodynamic conditions became weaker and favorable. Then Diphyphyllum once again flourished. Diphyphyllum clumps in the unit grew closely together in strong currents, with larger and more sparse individuals than in the lower units. A relatively low energy environment was formed between the Diphyphyllum clusters (Fig. 6C). Subsequently, Diphyphyllum could only grow in a limited area of suitability due to the disturbance of high-energy water brought about by short-term sea-level rise and fall. Afterwards, the environment became more favourable and Diphyphyllum expanded rapidly. As a result, the upper unit of Longjiangdong coral reef was formed, in which Diphyphyllum individuals were slightly larger than those in the first two units. Finally, because the kinetic energy of the water continued to weaken, the plaster deposition forced the whole coral reef to stop growing (Fig. 6D).Figure 6Micrographs of sediments in different positions of reef. (A) Calcareous bioclastic limestone, with biological particles accounting for about 70% of the debris. Abundant and diverse organisms indicate a medium-energy environment of the subtidal zone. Samples were taken from the bioclastic beach at the reef base. (B) Slightly larger bioclastics but lower biologic content than that in (A) suggest an increasing water energy. (C) Various bioclastic particles account for about 80% of the clastic particles contained in the calcareous bio-granular rock. The obviously small benthos indicate a low-energy environment in the subtidal zone barriered by the Diphyphyllum clusters. (D) Bioclastic grainstone is mainly composed of marl, with fine clastic particles (about 30%) and bedding. Low biomass indicates a low-energy environment of the subtidal zone. This figure is modified by the author with CorelDRAW (version 2022, and the URL link: https://www.coreldraw.com/cn/). Meaning of the letters in the figure: C crinoids, BF brachiopods, F foraminifera, B bryozoan, P pelletoid, MF mollusk shell fragment.Full size imageLongjiangdong patch reef started to develop in a relatively deep water environment. Diphyphyllum initially colonized and expanded in favorable conditions with the increase of water energy. Then the reef-builders transitioned from a single coral species to an assemblage of Diphyphyllum–Caninia–Lithostrotionella. These three coral species grew independently and contribute almost equally to the structure of the reef. However, the structure and function of the coral community were not yet stable enough. It was easily influenced by the weakening hydrodynamics and the increasing sedimentation, resulting in only small patch reefs.The Xinzhai layer reef was initialized by colonization and expansion of Lonsdaleia on bioclastic beach. Large coral clusters were formed in the presence of turbulent water. With the weakening of hydrodynamic conditions, an unknown branchlike organism and Antheria communities continued to develop separately in this area. Slender branchlike organisms expanded rapidly in these low-energy water environments until they were replaced by some individual corals as hydrodynamic energy increased. Each builder was short-lived in this layer reef, departing from the reef just at the beginning of colonization and expansion, due to rapidly changed hydrodynamic conditions.The evolution of reef-building corals in these four reefs indicated that both the coral assemblages and coral individuals would constantly adapt to the changing hydrodynamic conditions in Langping as sea level rose and fell. Although this was a reactive adjustment of coral populations in response to long-term environmental impacts, it was clearly positive for the building and development of coral reefs.Impact of disturbance on reef communitiesDisturbance is a relatively discontinuous event, which is ubiquitous in nature. It may indirectly affect the composition and population structure of reef communities by changing the environmental conditions, thus affect the structure and function of reef communities, even the evolution of the reef35. The major disturbances evident in these Mississippian framework reefs were associated with frequent changes of water flow, and drastic changes of climate and weather. These seem to be most obvious in the Langping platform due to its small size, with more frequent environmental influence evident on the reef communities in the study area.The most direct effect of disturbance events on coral reefs is the disruption of continuously evolving reef communities, which is common in coral reef studies. After the interruption caused by disturbances, some communities gradually recover due to the absence of continuous disturbance, or the dominant biota may be substituted by invading communities. The winner after interruption is decided by random factors to a large extent, in a ‘Competitive lottery’36. The conditions for the emergence of ‘Competitive lottery’ also include the need for species in a community to have similar abilities to invade discontinuities and to tolerate environmental conditions.Certainly, low-intensity disturbance does not necessarily produce discontinuity, but medium-intensity disturbance without discontinuity could directly impact on community species diversity. According to the ‘Moderate disturbance hypothesis’, moderate disturbance is conducive to a higher level of community diversity37. In environmental conditions with moderate intensity of disturbance, most species will not disappear entirely. The dominant pioneer species will also be restrained by disturbance to a certain extent, so large number of species can coexist, attaining the highest diversity35.The reef-builders in Langping are diverse compared with the Late Carboniferous reefs in Ziyun County10, which also developed in Dian-Qian-Gui Basin. More than 4 reef-building corals are identified in Xiadong reef, while 4 and 3 are in Xinzhai layer reef, Longjiangdong patch reef respectively. These reef-building corals, mostly Diphyphyllum, Lithostrotion, Siphonodendron and Lonsdaleia, were distributed irregularly in the reefs. Their ecological niche and function were likely similar and none of them was obviously dominant in the community (Fig. 7). This is in line with ‘Competitive lottery’ theory and the ‘Moderate disturbance hypothesis’.Figure 7Different species occupied the discontinuity surface irregularly. (A) Different reef-builders colonized and grew on the same hard substrate. (B) and (C) show detailed morphology of colony corals of (A). (D) Colony corals and a large number of individual corals grew together in a limited area, indicating equal colonization on the newly formed discontinuity surface. This figure is modified by the author from field photos with CorelDRAW (version 2022, and the URL link: https://www.coreldraw.com/cn/).Full size imageThe stability of a classical reef ecosystem includes the ability to withstand external disturbances and the ability to return to its original state once the disturbance is removed37,38. It is generally accepted that communities with high diversity are always more stable although ecosystem stability is not absolutely correlated to biodiversity35.There have been no reef-building corals with strong resistance and rapid recovery ability in the communities in Langping. None of these corals succeeded in developing into dominant species that can build reef shelves, which made the reefs in Langping mostly appear in the form of small patch reefs or reef layers. However, formation of the large reef in Xiadong Village, patch reef in Longjiangdong Village, and layer reef in Xinzhai Village were all related to their relatively high diversity of reef-building corals. Compared with the situation where only one reef-building organism dominated the Bianping large coral reef, Wengdao large phylloid algal reef and Ivanovia cf. manchurica patch reefs in Ziyun County10, Guizhou province, the different coral assemblages in Langping area could effectively adapt to changing hydrodynamic conditions and maintain reef growth.Species diversity increased by disturbance stabilized the ecosystemas shown during the construction of coral reefs in Langping.Effects of non-reef-builders on reef-building coralsBesides reef-building corals, there were a large number of reef-dwellers and off-reef organisms in the study area. Reef-dwellers referred to the species that didn’t directly contribute to reef growth in the community, mainly including various benthos and algaes39. Off-reef organisms are not part of the reef-building community, but also play an important role in participating in energy flow and providing organic matter to the reef ecosystem40.Common reef-dwelling organisms include crinoids, brachiopods, gastropods, various algae, foraminifera, bryozoans and individual corals. Crinoids were overwhelmingly dominant in numbers in the reef samples studied here.Carboniferous echinodermata in Guangxi Province reached its peak in Middle-to-Late Mississippian. In terms of amount and distribution, thick limestone with echinodermata debris in the carbonate platform were often dominated by crinoids41,42,43,44,45,46. The large number of crinoids in Langping excluded other metazoans and restricted the development of benthic reef-builders in Late Viséan-Serpukhovian in Langping, leading to poorly developed reef-building communities.Microorganisms and algaes had limited success in stablishing on the moving clastic beach in frequently disturbed water. There has not been obvious evidence of extensive “algal turf” in the coastal area of Langping platform. Only a few corals bonded by algal mats were observed47 (Fig. 8). In addition to their significant contribution to primary productivity, macroalgae were considered to play an important role in two aspects of coral reef ecosystems. One was to promote reef construction by its own binding and consolidation48,49. The other was to create a good condition for zoobenthos larvae to dwell and develop, thereby improving species diversity50. The limited productivity of algae in Langping constrained coral reef trophic inputs, which may then have limited populations of dependent metazoans. As a result, algaes and other metazoans were unable to achieve a variety of reef-building patterns, such as bonding, bounding, entanglement51,52. The reef framework in the study area was not stable in the presence of strong water flow, and the biological communities could not deal with frequent environmental changes, which were directly related to poor development of calcareous algae.Figure 8Micrographs of microbes and algaes. (A) Encrustations (indicated by black arrows) with distinct thickness around coral clusters formed by microbe and algal mats through bonding mud. The encrustations were formed before the clastic deposition (indicated by white arrows), showing the corals were living then. Microbes and algaes inside of the dense coral clusters had little impact on corals. (B) Single polarized micrograph showed clear and smooth boundaries of coral individuals without encrustation or drilling hole made by microbes or algaes. Few corals surrounded by bonding algaes could be observed in Langping, indicating that algaes were poorly developed between coral clusters. This figure is modified by the author from field photos with CorelDRAW (version 2022, and the URL link: https://www.coreldraw.com/cn/).Full size imageInfluence of coral morphology on reef developmentThe accumulation of reef structure had obvious impact on communities. Large reef structures could support abundance and diverse biota by modifying local environments and creating diverse conditions. Consequently the reef-building communities thrived between disturbances, stabilizing reef construction. In terms of large reef, the framework-building corals would play a key role in reef construction regardless of which kind of patterns was adopted. Therefore, reef-building corals with large size, rapid growth vertically, and strong resistance would become the biggest contributors to reef frame construction.The main reef-building corals in Langping were composed of Diphyphyllum, Lithostrotion, Siphonodendron, and Lonsdaleia, etc., being the dominant builders. These corals were similar in morphology such as cluster colony, thick and strong skeleton, and densely packed individuals (Fig. 9), which enabled them to resist water flow. At the same time, the upright colonies were adaptable to relatively calm water, being insensitive to mud deposition. The ecological characteristics of the Langping corals matched the gently sloping environment, the deep water environment and the rapidly changing energy of the currents. These cluster corals were able to colonize hard substrates and expand rapidly, thus altering the surrounding environment. The visible carbonate uplifts were formed with a large amount of benthos grouped into reef-building communities. These distinct uplifts constructed by coral clusters in different water conditions are composed of coral reefs of different sizes and appearance in the study area.Figure 9Main reef-building corals in the study area. (A) Diphyphyllum, (B) Lithostrotion, (C) Siphonodendron, (D) Lonsdaleia. (A) Rapidly grew clusters of main reef-building corals. The strong individuals are packed tightly when growing to support each other. This figure is modified by the author from field photos with CorelDRAW (version 2022, and the URL link: https://www.coreldraw.com/cn/).Full size imageThe complex and diverse local environments formed by large coral reefs can significantly increase benthic populations and improve reef species diversity. As a result, the nutrient flow in the community becomes complicated, and nutrients could be recycled effectively by reducing loss caused by water flow. Therefore, the overall productivity of large coral reef communities was always high. Complex trophic structure satisfied most of the benthos in the community with sufficient nutrients and inorganic salts.The morphology of reef-building corals in Langping enabled them to become predominant species in various water environments, which promoted the continued domed growth of coral reefs and facilitates the development of reef-building communities that form a variety of reefs. It suggests that the morphology of reef-building corals was a key prerequisite for reef development.In conclusion, coral reef communities are always constrained and influenced by environmental conditions. However, the ecology of the inhabitants is also an important factor in the formation of coral reefs. More

  • in

    The effects of visitors and social isolation from a peer on the behavior of a mixed-species pair of captive gibbons

    Kazarov, E. The Role of Zoos in Creating a Conservation Ethic in Visitors. SIT Digital Collections (2022). at https://digitalcollections.sit.edu/isp_collection/584.Hosey, G. How does the zoo environment affect the behaviour of captive primates?. Appl. Anim. Behav. Sci. 90, 107–129 (2005).
    Google Scholar 
    Morgan, K. & Tromborg, C. Sources of stress in captivity. Appl. Anim. Behav. Sci. 102, 262–302 (2007).
    Google Scholar 
    Sherwen, S. & Hemsworth, P. The visitor effect on zoo animals: Implications and opportunities for zoo animal welfare. Animals 9, 366 (2019).PubMed Central 

    Google Scholar 
    Chamove, A., Hosey, G. & Schaetzel, P. Visitors excite primates in zoos. Zoo Biol. 7, 359–369 (1988).
    Google Scholar 
    Tetley, C. L. & O’Hara, S. J. Ratings of animal personality as a tool for improving the breeding, management and welfare of zoo mammals. Anim. Welf. UFAW J. 21(4), 463 (2012).CAS 

    Google Scholar 
    Stoinski, T. S., Jaicks, H. F. & Drayton, L. A. Visitor effects on the behavior of captive western lowland gorillas: The importance of individual differences in examining welfare. Zoo Biol. 31(5), 586–599 (2012).PubMed 

    Google Scholar 
    Queiroz, M. B. & Young, R. J. The different physical and behavioural characteristics of zoo mammals that influence their response to visitors. Animals 8(8), 139 (2018).PubMed Central 

    Google Scholar 
    Fanson, K. V. & Wielebnowski, N. C. Effect of housing and husbandry practices on adrenocortical activity in captive Canada lynx (Lynx canadensis). Anim. Welf. 22, 159–165 (2013).CAS 

    Google Scholar 
    Pirovino, M. et al. Fecal glucocorticoid measurements and their relation to rearing, behavior, and environmental factors in the population of pileated gibbons (Hylobates pileatus) held in European zoos. Int. J. Primatol. 32(5), 1161–1178 (2011).
    Google Scholar 
    Williams, I., Hoppitt, W. & Grant, R. The effect of auditory enrichment, rearing method and social environment on the behavior of zoo-housed psittacines (Aves: Psittaciformes); implications for welfare. Appl. Anim. Behav. Sci. 186, 85–92 (2017).
    Google Scholar 
    Fernandez, E., Tamborski, M., Pickens, S. & Timberlake, W. Animal–visitor interactions in the modern zoo: Conflicts and interventions. Appl. Anim. Behav. Sci. 120, 1–8 (2009).
    Google Scholar 
    Hosey, G. & Skyner, L. Self-injurious behavior in zoo primates. Int. J. Primatol. 28, 1431–1437 (2007).
    Google Scholar 
    Mallapur, A., Sinha, A. & Waran, N. Influence of visitor presence on the behaviour of captive lion-tailed macaques (Macaca silenus) housed in Indian zoos. Appl. Anim. Behav. Sci. 94, 341–352 (2005).
    Google Scholar 
    Davey, G. Visitors’ Effects on the Welfare of Animals in the Zoo: A Review. J. Appl. Anim. Welf. Sci. 10, 169–183 (2007).CAS 
    PubMed 

    Google Scholar 
    Jones, H., McGregor, P., Farmer, H. & Baker, K. The influence of visitor interaction on the behavior of captive crowned lemurs (Eulemur coronatus) and implications for welfare. Zoo Biol. 35, 222–227 (2016).CAS 
    PubMed 

    Google Scholar 
    Cook, S. & Hosey, G. R. Interaction sequences between chimpanzees and human visitors at the zoo. Zoo Biol. 14(5), 431–440 (1995).
    Google Scholar 
    Baker, K. C. Benefits of positive human interaction for socially-housed chimpanzees. Anim. Welf. (South Mimms, Engl.nd) 13(2), 239 (2004).CAS 

    Google Scholar 
    Carder, G. & Semple, S. Visitor effects on anxiety in two captive groups of western lowland gorillas. Appl. Anim. Behav. Sci. 115, 211–220 (2008).
    Google Scholar 
    Wood, W. Interactions among environmental enrichment, viewing crowds, and zoo chimpanzees (Pantroglodytes). Zoo Biol. 17, 211–230 (1998).
    Google Scholar 
    Todd, P., Macdonald, C. & Coleman, D. Visitor-associated variation in captive Diana monkey (Cercopithecus diana diana) behaviour. Appl. Anim. Behav. Sci. 107, 162–165 (2007).
    Google Scholar 
    Davis, N., Schaffner, C. & Smith, T. Evidence that zoo visitors influence HPA activity in spider monkeys (Ateles geoffroyii rufiventris). Appl. Anim. Behav. Sci. 90, 131–141 (2005).
    Google Scholar 
    Sherwen, S. L. et al. Effects of visual contact with zoo visitors on black-capped capuchin welfare. Appl. Anim. Behav. Sci. 167, 65–73 (2015).
    Google Scholar 
    Choo, Y., Todd, P. & Li, D. Visitor effects on zoo orangutans in two novel, naturalistic enclosures. Appl. Anim. Behav. Sci. 133, 78–86 (2011).
    Google Scholar 
    Sherwen, S., Magrath, M., Butler, K., Phillips, C. & Hemsworth, P. A multi-enclosure study investigating the behavioural response of meerkats to zoo visitors. Appl. Anim. Behav. Sci. 156, 70–77 (2014).
    Google Scholar 
    Hosey, G. & Druck, P. The influence of zoo visitors on the behaviour of captive primates. Appl. Anim. Behav. Sci. 18, 19–29 (1987).
    Google Scholar 
    Mitchell, G. et al. More on the ‘influence’of zoo visitors on the behaviour of captive primates. Appl. Anim. Behav. Sci. 35(2), 189–198 (1992).
    Google Scholar 
    Sellinger, R. & Ha, J. The effects of visitor density and intensity on the behavior of two captive jaguars (Panthera onca). J. Appl. Anim. Welfare Sci. 8, 233–244 (2005).CAS 

    Google Scholar 
    Azevedo, C., Lima, M., Silva, V., Young, R. & Rodrigues, M. Visitor Influence on the Behavior of Captive Greater Rheas (Rhea americana, Rheidae Aves). J. Appl. Anim. Welfare Sci. 15, 113–125 (2012).
    Google Scholar 
    Das Gupta, M., Das, A., Sumy, M. C. & Islam, M. M. An explorative study on visitor’s behaviour and their effect on the behaviour of primates at Chittagong zoo. Bangladesh J. Vet. Anim. Sci. 5(2), 24–32 (2017).
    Google Scholar 
    Hemsworth, P. Human–animal interactions in livestock production. Appl. Anim. Behav. Sci. 81, 185–198 (2003).
    Google Scholar 
    Stoinski, T., Czekala, N., Lukas, K. & Maple, T. Urinary androgen and corticoid levels in captive, male Western lowland gorillas (Gorilla g. gorilla): Age- and social group-related differences. Am. J. Primatol. 56, 73–87 (2002).CAS 
    PubMed 

    Google Scholar 
    Stoinski, T., Lukas, K., Kuhar, C. & Maple, T. Factors influencing the formation and maintenance of all-male gorilla groups in captivity. Zoo Biol. 23, 189–203 (2004).
    Google Scholar 
    Olsson, I. & Westlund, K. More than numbers matter: The effect of social factors on behaviour and welfare of laboratory rodents and non-human primates. Appl. Anim. Behav. Sci. 103, 229–254 (2007).
    Google Scholar 
    Martin, J. E. Early life experiences: Activity levels and abnormal behaviours in resocialised chimpanzees. Anim Welf. 11(4), 419–436 (2002).CAS 

    Google Scholar 
    Birkett, L. P. & Newton-Fisher, N. E. How abnormal is the behaviour of captive, zoo-living chimpanzees?. PLoS ONE 6(6), e20101 (2011).ADS 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Ballen, C., Shine, R. & Olsson, M. Effects of early social isolation on the behaviour and performance of juvenile lizards Chamaeleo calyptratus. Anim. Behav. 88, 1–6 (2014).
    Google Scholar 
    Coe, C., Mendoza, S., Smotherman, W. & Levine, S. Mother-infant attachment in the squirrel monkey: Adrenal response to separation. Behav. Biol. 22, 256–263 (1978).CAS 
    PubMed 

    Google Scholar 
    Mendoza, S., Smotherman, W., Miner, M., Kaplan, J. & Levine, S. Pituitary-adrenal response to separation in mother and infant squirrel monkeys. Dev. Psychobiol. 11, 169–175 (1978).CAS 
    PubMed 

    Google Scholar 
    Gilbert, M. & Baker, K. Social buffering in adult male rhesus macaques (Macaca mulatta): Effects of stressful events in single vs. pair housing. J. Med. Primatol. 40, 71–78 (2010).PubMed 
    PubMed Central 

    Google Scholar 
    Schapiro, S. Effects of social manipulations and environmental enrichment on behavior and cell-mediated immune responses in rhesus macaques. Pharmacol. Biochem. Behav. 73, 271–278 (2002).CAS 
    PubMed 

    Google Scholar 
    Chen, W. et al. Effects of social isolation and re-socialization on cognition and ADAR1 (p110) expression in mice. PeerJ 4, e2306 (2016).PubMed 
    PubMed Central 

    Google Scholar 
    Glatston, A., Geilvoet-Soeteman, E., Hora-Pecek, E. & Van Hooff, J. The influence of the zoo environment on social behavior of groups of cotton-topped tamarins Saguinus oedipus oedipus. Zoo Biol. 3, 241–253 (1984).
    Google Scholar 
    Mitchell, G. et al. Effects of visitors and cage changes on the behaviors of mangabeys. Zoo Biol. 10, 417–423 (1991).
    Google Scholar 
    Geissmann, T. & Orgeldinger, M. The relationship between duet songs and pair bonds in siamangs Hylobates syndactylus. Anim. Behav. 60, 805–809 (2000).CAS 
    PubMed 

    Google Scholar 
    Palombit, R. Pair bonds in monogamous apes: A comparison of the siamang hylobates syndactylus and the white-handed gibbon hylobates lar. Behaviour 133, 321–356 (1996).
    Google Scholar 
    Rutberg, A. The evolution of monogamy in primates. J. Theor. Biol. 104, 93–112 (1983).ADS 
    CAS 
    PubMed 

    Google Scholar 
    Giorgi, A., Montebovi, G., Vitale, A. & Alleva, E. A behavioural case study of early social isolation of a subadult white-handed gibbon (Hylobates lar). Folia Primatol. 89, 287–294 (2018).
    Google Scholar 
    Skynner, L. A., Amory, J. R. & Hosey, G. The effect of visitors on the self-injurious behaviour of a male pileated gibbon (Hylobates pileatus). Zool. Garten 74(1), 38–41 (2004).
    Google Scholar 
    Smith, K. & Kuhar, C. Siamangs (Hylobates syndactylus) and white-cheeked gibbons (Hylobates leucogenys) show few behavioral differences related to zoo attendance. J. Appl. Anim. Welfare Sci. 13, 154–163 (2010).CAS 

    Google Scholar 
    Lukas, K. E. et al. Longitudinal study of delayed reproductive success in a pair of white-cheeked gibbons (Hylobates leucogenys). Zoo Biol. 21, 413–434 (2002).
    Google Scholar 
    Cooke, C. & Schillaci, M. Behavioral responses to the zoo environment by white handed gibbons. Appl. Anim. Behav. Sci. 106, 125–133 (2007).
    Google Scholar 
    Mootnick, A. & Baker, E. Masturbation in captiveHylobates (gibbons). Zoo Biol. 13, 345–353 (1994).
    Google Scholar 
    Geissmann, T. Reassessment of age of sexual maturity in gibbons (hylobates spp.). American Journal of Primatology 23, 11–22 (1991).Altmann, J. Observational study of behavior: Sampling methods. Behaviour 49(3–4), 227–266 (1974).CAS 
    PubMed 

    Google Scholar 
    Pomerantz, O. & Terkel, J. Effects of positive reinforcement training techniques on the psychological welfare of zoo-housed chimpanzees (Pan troglodytes). Am. J. Primatol. 71, 687–695 (2009).PubMed 

    Google Scholar 
    Orgeldinger, M. Protective and territorial behavior in captive siamangs (Hylobates syndactylus). Zoo Biol. 16, 309–325 (1997).
    Google Scholar 
    Fox, J. et al. Package ‘car’. Vienna: R Foundation for Statistical Computing, 16 https://cran.uni-muenster.de/web/packages/car/car.pdf (2012).Magnusson, A., Skaug, H., Nielsen, A., Berg, C., Kristensen, K., Maechler, M., van Bentham, K., Bolker, B., Brooks, M. & Brooks, M. M. Package ‘glmmtmb’. R Package Version 0.2. 0 (2017).Hartig, F., & Hartig, M. F. Package ‘DHARMa’. Vienna, Austria: R Development Core Team (2017).Troisi, A. Displacement activities as a behavioral measure of stress in nonhuman primates and human subjects. Stress 5, 47–54 (2002).PubMed 

    Google Scholar 
    Baker, K. & Aureli, F. Behavioural indicators of anxiety: An empirical test in chimpanzees. Behaviour 134, 1031–1050 (1997).
    Google Scholar 
    Vick, S. J. & Paukner, A. Variation and context of yawns in captive chimpanzees (Pan troglodytes). Am. J. Primatol. Off. J. Am. Soc. Primatol. 72(3), 262–269 (2010).
    Google Scholar 
    Norscia, I. & Palagi, E. When play is a family business: Adult play, hierarchy, and possible stress reduction in common marmosets. Primates 52, 101–104 (2010).PubMed 

    Google Scholar 
    Held, S. & Špinka, M. Animal play and animal welfare. Anim. Behav. 81, 891–899 (2011).
    Google Scholar 
    Davey, G. Visitor behavior in zoos: A review. Anthrozoös 19, 143–157 (2006).
    Google Scholar 
    Nimon, A. & Dalziel, F. Cross-species interaction and communication: a study method applied to captive siamang (Hylobates syndactylus) and long-billed corella (Cacatua tenuirostris) contacts with humans. Appl. Anim. Behav. Sci. 33, 261–272 (1992).
    Google Scholar 
    Suomi, S. Early determinants of behaviour: Evidence from primate studies. Br. Med. Bull. 53, 170–184 (1997).CAS 
    PubMed 

    Google Scholar 
    Anderson, J. & Chamove, A. Self-aggression and social aggression in laboratory-reared macaques. J. Abnorm. Psychol. 89, 539–550 (1980).CAS 
    PubMed 

    Google Scholar 
    Mallapur, A. & Choudhury, B. Behavioral abnormalities in captive nonhuman primates. J. Appl. Anim. Welfare Sci. 6, 275–284 (2003).CAS 

    Google Scholar 
    Barlow, C., Caldwell, C. & Lee, P. Individual differences and response to visitors in zoo-housed diana monkeys (Cercopithecus diana diana). Cabdirect.org (2022). at https://www.cabdirect.org/cabdirect/abstract/20123180753.Gartner, M. & Weiss, A. Studying primate personality in zoos: Implications for the management, welfare and conservation of great apes. International Zoo Yearbook 52, 79–91 (2018).
    Google Scholar 
    Mitchell, G., Raymond, E., Ruppenthal, G. & Harlow, H. Long-term effects of total social isolation upon behavior of rhesus monkeys. Psychol. Rep. 18, 567–580 (1966).
    Google Scholar 
    Martín, O., Vinyoles, D., García-Galea, E. & Maté, C. Improving the welfare of a zoo-housed male drill (Mandrillus leucophaeus poensis) aggressive toward visitors. J. Appl. Anim. Welfare Sci. 19, 323–334 (2016).
    Google Scholar 
    Ross, S., Melber, L., Gillespie, K. & Lukas, K. The impact of a modern, naturalistic exhibit design on visitor behavior: A cross-facility comparison. Visitor Stud. 15, 3–15 (2012).
    Google Scholar 
    Quadros, S., Goulart, V., Passos, L., Vecci, M. & Young, R. Zoo visitor effect on mammal behaviour: Does noise matter?. Appl. Anim. Behav. Sci. 156, 78–84 (2014).
    Google Scholar 
    Bonnie, K., Ang, M. & Ross, S. Effects of crowd size on exhibit use by and behavior of chimpanzees (Pan troglodytes) and Western lowland gorillas (Gorilla gorilla) at a zoo. Appl. Anim. Behav. Sci. 178, 102–110 (2016).
    Google Scholar  More

  • in

    Effect of marigold (Tagetes erecta L.) on soil microbial communities in continuously cropped tobacco fields

    Chen, X. L. et al. Effects of Meloidogyne incognitaon the fungal community in tobaccorhizosphere. Rev. Bras. Cienc. Solo. 46, e0210127 (2022).
    Google Scholar 
    Zhang, S. X. et al. Research progresses on continuous cropping obstacles of tobacco. Soils 47(5), 823–829 (2015).CAS 

    Google Scholar 
    Luo, J. Y. et al. Effects of soil salinity onrhizosphere soil microbes in transgenic Bt cotton fields. J. Integr. Agric. 16, 1624–1633 (2017).CAS 

    Google Scholar 
    Chaparro, J. M. et al. Manipulating the soil microbiome to increase soil health and plant fertility. Biol. Fertil. Soils 48, 489–499 (2012).
    Google Scholar 
    Newton, A., Begg, G. & Swanston, J. Deployment of diversity for enhanced crop function. Ann. Appl. Biol. 154, 309–322 (2009).
    Google Scholar 
    Li, X. G. et al. Effects of intercropping with Atractylodeslancea and application of bio-organic fertiliser on soil invertebrates, disease control and peanut productivity in continuouspeanut cropping field in subtropical China. Agrofor. Syst. 88, 41–52 (2014).
    Google Scholar 
    Ahmed, W. et al. Ralstonia solanacearum, a deadly pathogen: Revisiting the bacterial wilt biocontrol practices in tobacco and other Solanaceae. Rhizosphere 21, 100479 (2022).
    Google Scholar 
    Gómez-Rodrıguez, O., Zavaleta-Mejıa, E., Gonzalez-Hernandez, V., Livera-Munoz, M. & Cárdenas-Soriano, E. Allelopathyand microclimatic modification of intercropping with marigold on tomato early blight disease development. Field Crops Res. 83, 27–34 (2003).
    Google Scholar 
    Weidenhamer, J. D., Montgomery, T. M., Cipollini, D. F., Weston, P. A. & Mohney, B. K. Plandensity and rhizosphere chemistry: Does marigold root exudate composition respond to intra-and interspecific competition?. J. Chem. Ecol. 45(5–6), 525–533 (2019).CAS 
    PubMed 

    Google Scholar 
    Ploeg, A. T. Effects of selected marigold varieties on root-knot nematodes and tomato and melon yields. Plant Dis. 86(5), 505–508 (2002).PubMed 

    Google Scholar 
    Hooks, C. R., Wang, K. H., Ploeg, A. & McSorley, R. Using marigold (Tagetes spp.) as a cover crop to protect crops fromplant-parasitic nematodes. Appl. Soil Ecol. 46, 307–320 (2010).
    Google Scholar 
    Li, W., Xu, J., Chen, H. & Qi, Y. Phytochemicals and their biological activities of plants in tagetes l.-sciencedirect. Chin. Herbal Med. 4(2), 103–117 (2012).
    Google Scholar 
    Weidenhamer, J. D., Mohney, B. K., Shihada, N. & Rupasinghe, M. Spatial and temporal dynamics of root exudation: How important is heterogeneity in allelopathic interactions?. J. Chem. Ecol. 40(8), 940–952 (2014).CAS 
    PubMed 

    Google Scholar 
    Marotti, I. et al. Thiophene occurrence in different tagetes species: Agricultural biomasses as sources ofbiocidal substances. J. Sci. Food Agric. 90(7), 1210–1217 (2010).CAS 
    PubMed 

    Google Scholar 
    Barto, E. K. et al. The fungal fastlane: Common mycorrhizal networks extendbioactive zones of allelochemicals in soils. PLoS ONE 6, e27195 (2011).ADS 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Evenhuis, A., Korthals, G. & Molendijk, L. Tagetes patula as an effective catch crop forlong-term control of Pratylenchus penetrans. Nematology 6, 877–881 (2004).
    Google Scholar 
    Wu, W. T. et al. Effects of marigold-tobacco rotation on soil nematode community composition. Southwest China J. Agric. Sci. 32(2), 342–348 (2019).
    Google Scholar 
    Reynolds, L. B., Potter, J. W. & Ball-Coelho, B. R. Crop rotation with sp. is an alternative to chemical fumigation for control of root-lesion nematodes. Agron. J. 92(5), 957–966 (2000).
    Google Scholar 
    El-Hamawi, M., Youssef, M. & Zawam, H. S. Management of Meloidogyne incognita, the root-knot nematode, on soybean asaffected by marigold and sea ambrosia (damsisa) plants. J. Pest Sci. 77, 95–98 (2004).
    Google Scholar 
    Kumar, N., Krishnappa, K., Reddy, B., Ravichandra, N. & Karuna, K. Intercropping for the management of root-knotnematode, Meloidogyne incognitain vegetable-based cropping systems. Indian J. Nematol. 35, 46–49 (2005).
    Google Scholar 
    Zhang, J. et al. Crop rotation with marigold promotes soil bacterial structure to assist in mitigating clubroot Incidence in Chinese Cabbage. Plants 11(17), 2295 (2022).CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Xia, T. Y. et al. Microbial diversity of tobacco rhizospheresoil in different growth stages of marigold-tobacco intercropping system. Southwest China J. Agric. Sci. 31(4), 680–686 (2018).
    Google Scholar 
    Wei, H. Y. et al. Effects of marigold diversified cropping with angelica on fungal community in soils. Plant Prot. 41(5), 69–74 (2015).MathSciNet 
    CAS 

    Google Scholar 
    Li, Y. et al. Intercropping with marigold promotes soil health and microbialstructure to assist in mitigating tobacco bacterial wilt. J. Plant Pathol. 102, 731–742 (2020).
    Google Scholar 
    Caporaso, J. G. et al. Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proc. Natl. Acad. Sci. 108, 4516–4522 (2011).ADS 
    CAS 
    PubMed 

    Google Scholar 
    Edgar, R. C. UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 10, 996–998 (2013).CAS 
    PubMed 

    Google Scholar 
    Wang, Q., Garrity, G. M., Tiedje, J. M. & Cole, J. R. Naive bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. MicroBiol. 73, 5261–5267 (2007).ADS 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Altschul, S. F., Gish, W., Miller, W., Myers, E. W. & Lipman, D. J. Basic local alignment search tool. J. Mol. Biol. 215, 403–410 (1990).CAS 
    PubMed 

    Google Scholar 
    Irikiin, I. et al. Rhizobacterial community-level, sole carbon source utilization pattern aff ects the delay in the bacterial wilt of tomato grown in rhizobacterial community model system. Appl. Soil Ecol. 34(1), 27–32 (2006).
    Google Scholar 
    Wu, M. N. et al. Soil fungistasis and its relations to soil microbial composition and diversity: A case study of a series of soils with different fungistasis. J. Environ. Sci. 20(7), 871–877 (2008).CAS 

    Google Scholar 
    Mendes, L. W. et al. Soil-Borne microbiome: Linking diversity to function. Microb. Ecol. 70(1), 255–265 (2015).CAS 
    PubMed 

    Google Scholar 
    Jaiswal, A. K. et al. Linking the belowground microbial composition, diversity and activity to soilborne disease suppression and growth promotion of tomato amended with biochar. Sci. Rep. 7, 44382 (2017).ADS 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Raaijmakers, J. M. & Mazzola, M. Soil immune responses soil microbiomes may be harnessed for plant health. Science 352, 1392–1393 (2016).ADS 
    CAS 
    PubMed 

    Google Scholar 
    Kušlienė, G., Rasmussen, J., Kuzyakov, Y. & Eriksen, J. Medium-term response of microbial community to rhizodeposits of white clover and ryegrass and tracing of active processes induced by 13C and 15N labelled exudates. Soil Biol. Biochem. 76, 22–33 (2014).
    Google Scholar 
    Mohammadi, K. Soil microbial activity and biomass as influenced by tillage and fertilization in wheat production. Am.-Eurasian J. Agric. Environ. Sci. 10, 330–337 (2011).
    Google Scholar 
    Wang, G. H. et al. Research progress of Acidobacteria ecology in soils. Biotechnol. Bull. 32(2), 14–20 (2016).
    Google Scholar 
    Wei, H., Wang, L., Hassan, M. & Xie, B. Succession of the functional microbial communities and the metabolic functions in maize straw composting process. Bioresour. Technol. 256, 333–341 (2018).CAS 
    PubMed 

    Google Scholar 
    Wang, Y., Liu, L., Yang, J., Duan, Y. & Zhao, Z. The diversity of microbial community and function varied in response to different agricultural residues composting. Sci. Total Environ. 715, 136983 (2020).ADS 
    CAS 
    PubMed 

    Google Scholar 
    Glass, N. L., Schmoll, M., Cate, J. H. & Coradetti, S. Plant cell wall deconstruction by ascomycete fungi. Annu. Rev. Microbiol. 67, 477–498 (2013).CAS 
    PubMed 

    Google Scholar 
    Li, Y. et al. Linking soil fungal community structure and function to soil organic carbon chemical composition in intensively managed subtropical bamboo forests. Soil Biol. Biochem. 107, 19–31 (2017).CAS 

    Google Scholar 
    Martins, L. F., Kolling, D., Camassola, M., Dillon, A. J. & Ramos, L. P. Comparison of Penicillium echinulatumand Trichoderma reeseicellulases in relation to their activity against various cellulosic substrates. Bioresour. Technol. 99, 1417–1424 (2008).CAS 
    PubMed 

    Google Scholar  More

  • in

    Tidal effects on periodical variations in the occurrence of singing humpback whales in coastal waters of Chichijima Island, Ogasawara, Japan

    Morrison, M. A., Francis, M. P., Hartill, B. W. & Parkinson, D. M. Diurnal and tidal variation in the abundance of the fish fauna of a temperate tidal mudflat. Estuar. Coast. Shelf Sci. 54, 793–807 (2002).Article 
    ADS 

    Google Scholar 
    Ribeiro, J. et al. Seasonal, tidal and diurnal changes in fish assemblages in the Ria Formosa lagoon (Portugal). Estuar. Coast. Shelf Sci. 67, 461–474 (2006).Article 
    ADS 

    Google Scholar 
    Takemura, A., Rahman, M. S. & Park, Y. J. External and internal controls of lunar-related reproductive rhythms in fishes. J. Fish Biol. 76, 7–26 (2010).Article 
    CAS 
    PubMed 

    Google Scholar 
    Mendes, S., Turrell, W., Lütkebohle, T. & Thompson, P. Influence of the tidal cycle and a tidal intrusion front on the spatio-temporal distribution of coastal bottlenose dolphins. Mar. Ecol. Prog. Ser. 239, 221–229 (2002).Article 
    ADS 

    Google Scholar 
    Johnston, D. W., Thorne, L. H. & Read, A. J. Fin whales Balaenoptera physalus and minke whales Balaenoptera acutorostrata exploit a tidally driven island wake ecosystem in the Bay of Fundy. Mar. Ecol. Prog. Ser. 305, 287–295 (2005).Article 
    ADS 

    Google Scholar 
    Ichikawa, K. et al. Dugong (Dugong dugon) vocalization patterns recorded by automatic underwater sound monitoring systems. J. Acoust. Soc. Am. 119, 3726–3733 (2006).Article 
    ADS 
    PubMed 

    Google Scholar 
    Akamatsu, T. et al. Seasonal and diurnal presence of finless porpoises at a corridor to the ocean from their habitat. Mar. Biol. 157, 1879–1887 (2010).Article 

    Google Scholar 
    Li, S. et al. Seasonal, lunar and tidal influences on habitat use of indo-pacific humpback dolphins in Beibu gulf, China. Zool. Stud. https://doi.org/10.6620/ZS.2018.57-01 (2018).Article 
    PubMed 
    PubMed Central 

    Google Scholar 
    Zamon, J. E. Seal predation on salmon and forage fish schools as a function of tidal currents in the San Juan Islands, Washington, USA. Fish. Oceanogr. 10, 353–366 (2001).Article 

    Google Scholar 
    Van Parijs, S. M., Hastie, G. D. & Thompson, P. M. Geographical variation in temporal and spatial vocalization patterns of male harbour seals in the mating season. Anim. Behav. 58, 1231–1239 (1999).Article 
    PubMed 

    Google Scholar 
    Bortolotto, G. A., Danilewicz, D., Hammond, P. S., Thomas, L. & Zerbini, A. N. Whale distribution in a breeding area: Spatial models of habitat use and abundance of western South Atlantic humpback whales. Mar. Ecol. Prog. Ser. 585, 213–227 (2017).Article 
    ADS 

    Google Scholar 
    Johnson, J. H. & Wolman, A. A. The humpback whale, Megaptera novaeangliae. Mar. Fish. Rev. 46, 30–37 (1984).
    Google Scholar 
    Kobayashi, N. et al. Spatial distribution and habitat use patterns of humpback whales in Okinawa, Japan. Mammal Study 41, 207–214 (2016).Article 

    Google Scholar 
    Mori, K., Sata, F., Yamaguchi, M., Suganuma, H. & Ueyanagi, S. Distribution, migration and local movements of humpback whale (Megaptera novaeangliae) in the adjacent waters of the Ogasawara (Bonin) Islands Japan. J. Fac. Mar. Sci. Technol. Tokai Univ. 45, 197–213 (1998).
    Google Scholar 
    Rasmussen, K., Calambokidis, J. & Steiger, G. H. Distribution and migratory destinations of humpback whales off the Pacific coast of Central America during the boreal winters of 1996–2003. Mar. Mammal Sci. 28, 1–13 (2012).Article 

    Google Scholar 
    Calambokidis, J. et al. SPLASH: structure of populations, levels of abuncance and status of humpback whales in the North Pacific. Final report for Contract AB133F-03-RP-00078, to U.S. Dept. of Comm. Western Administrative Center, Seattle, WA. https://cascadiaresearch.org/files/SPLASH-contract-Report-May08.pdf (2008).Hill, M. et al. Found: A missing breeding ground for endangered western North Pacific humpback whales in the Mariana Archipelago. Endanger. Species Res. 41, 91–103 (2020).Article 

    Google Scholar 
    Payne, R. S. & McVay, S. Songs of humpback whales. Science 173, 585–597 (1971).Article 
    ADS 
    CAS 
    PubMed 

    Google Scholar 
    Winn, H. E. & Winn, L. The song of the humpback whale Megaptera novaeangliae in the West Indies. Mar. Biol. 47, 97–114 (1978).Article 

    Google Scholar 
    Tyack, P. Interactions between singing Hawaiian humpback whales and conspecifics nearby. Behav. Ecol. Sociobiol. 8, 105–116 (1981).Article 

    Google Scholar 
    Herman, L. M. The multiple functions of male song within the humpback whale (Megaptera novaeangliae) mating system: Review, evaluation, and synthesis. Biol. Rev. 92, 1795–1818 (2017).Article 
    PubMed 

    Google Scholar 
    Au, W. W. L., Mobley, J., Burgess, W. C., Lammers, M. O. & Nachtigall, P. E. Seasonal and diurnal trends of chorusing humpback whales wintering in waters off western Maui. Mar. Mammal Sci. 16, 530–544 (2000).Article 

    Google Scholar 
    Cerchio, S., Collins, T., Strindberg, S., Bennett, C. & Rosenbaum, H. Humpback whale singing activity off northern Angola: An indication of the migratory cycle, breeding habitat and impact of seismic surveys on singer number in Breeding. Int. Whal. Comm. P. SC/62/SH12 (2010).Kobayashi, N., Okabe, H., Higashi, N., Miyahara, H. & Uchida, S. Diel patterns in singing activity of humpback whales in a winter breeding area in Okinawan (Ryukyuan) waters. Mar. Mammal Sci. 37, 982–992 (2021).Article 

    Google Scholar 
    Munger, L. M., Lammers, M. O., Fisher-Pool, P. & Wong, K. Humpback whale (Megaptera novaeangliae) song occurrence at American Samoa in long-term passive acoustic recordings, 2008–2009. J. Acoust. Soc. Am. 132, 2265–2272 (2012).Article 
    ADS 
    PubMed 

    Google Scholar 
    Barlow, D. R., Fournet, M. & Sharpe, F. Incorporating tides into the acoustic ecology of humpback whales. Mar. Mammal Sci. 35, 234–251 (2019).Article 

    Google Scholar 
    Chenoweth, E., Gabriele, C. & Hill, D. Tidal influences on humpback whale habitat selection near headlands. Mar. Ecol. Prog. Ser. 423, 279–289 (2011).Article 
    ADS 

    Google Scholar 
    Sousa-Lima, R. S., Clark, C. W. & Road, S. W. Modeling the effect of boat traffic on singing activity of humpback whales (Megaptera novaeangliae) in the abrolhos national marine park, Brazil. Can. Acoust 36, 174–181 (2008).
    Google Scholar 
    Cerchio, S., Strindberg, S., Collins, T., Bennett, C. & Rosenbaum, H. Seismic surveys negatively affect humpback whale singing activity off Northern Angola. PLoS ONE 9, e86464. https://doi.org/10.1371/journal.pone.0086464 (2014).Article 
    ADS 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Darling, J. D. & Mori, K. Recent observations of humpback whales (Megaptera novaeangliae) in Japanese waters off Ogasawara and Okinawa. Can. J. Zool. 71, 325–333 (1993).Article 

    Google Scholar 
    Calambokidis, J. et al. Movements and population structure of humpback whales in the North Pacific. Mar. Mammal Sci. 17, 769–794 (2001).Article 

    Google Scholar 
    Wessel, P., Smith, W. H. F., Scharroo, R., Luis, J. & Wobbe, F. Generic mapping tools: Improved version released. Eos Trans. Am. Geophys. Union 94, 409–410 (2013).Article 
    ADS 

    Google Scholar 
    Helweg, D. A. & Herman, L. M. Diurnal patterns of behaviour and group membership of humpback whales (Megaptera novaeangliae) wintering in Hawaiian waters. Ethology 98, 298–311 (1994).Article 

    Google Scholar 
    Darling, J. D. & Berube, M. Interactions of singing humpback whales with other males. Mar. Mammal Sci. 17, 570–584 (2001).Article 

    Google Scholar 
    Whitlow, W. L. et al. Acoustic properties of humpback whale songs. J. Acoust. Soc. Am. 120, 1103–1110 (2006).Article 

    Google Scholar 
    Japan Coast Guard. Sailing Directions for South and East Coasts of Honshu. (1981).Tsujii, K. et al. Change in singing behavior of humpback whales caused by shipping noise. PLoS ONE 13, e0204112. https://doi.org/10.1371/journal.pone.0204112 (2018).Article 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Ryan, J. P. et al. Humpback whale song occurrence reflects ecosystem variability in feeding and migratory habitat of the northeast Pacific. PLoS ONE 14, e0222456. https://doi.org/10.1371/journal.pone.0222456 (2019).Article 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    R Core Team. R: A Language and Environment for Statistical Computing. 4.0.0 version. R Foundation for Statistical Computing, Vienna, Austria. http://www.R-project.org/ (2020).Wood, S.N. Generalized Additive Models: An Introduction with R 2nd edn, (Chapman and Hall/CRC, 2017). More

  • in

    Numerical analysis of the relationship between mixing regime, nutrient status, and climatic variables in Lake Biwa

    Model validationBased on the time-series validations of water temperature and DO concentration, model accuracy improved gradually, despite several discrepancies at the beginning of the simulation (Supplementary Fig. S1). The model is primarily driven by a set of boundary data, including wind speed, solar radiation, and precipitation data24,25. From this perspective, more high-quality boundary data promotes better numerical reproducibility. However, meteorological data collection was challenging due to the early observation equipment limitations and low observational accuracy compared to current data. The temporal inconsistency of accuracy in observational data has been eliminated to a large extent by fitting a regression curve24. Spatial resolution is the other issue. Possessing spatially constant values for all boundary conditions complicates the numerical reproducibility of variations on finer scales.The relationship between turnovers and the curve shape of water temperature versus DO concentration is theoretically sound27,28. In the last stage of stratification in the lake, water temperature and DO concentration near the bottom are more likely to slightly increase due to thermal diffusion and DO supplies from the upper water. If a turnover occurs, the whole column of water is mixed strongly (Supplementary Fig. S3). Bottom water temperature decreases due to surface water cooling, and DO concentration increases, due to surface water replenishment and increased oxygen solubility. If the turnover fails, only the partial column of water is mixed, causing a delay in the timing of deep-water renewal (Supplementary Fig. S3). However, the upper water in later months, like that in March, has been rapidly warmed, resulting in an increase in the bottom water temperature. For example, in 2007 and 2016, the simulated water temperature and DO concentration fluctuated within a limited range in February and then skyrocketed in March, after mixing with the warmed surface water (blue points in Supplementary Fig. S4). On the other hand, explicit definitions of turnover timing are challenging. The threshold used to judge turnover timing is reliable because the results matched the observation. The turnover timing varied by 36 days in Lake Biwa during the simulation period, which is comparable to that observed in other lakes, such as approximately 21 days in Heiligensee, Germany over a 17-year timespan29, 16 days in Lake Washington over a 40-year timespan30, and 28 days in Blelham Tarn over a 41-year timespan31.Variables affecting the mixing regimeDetermining variables that affect the mixing regime is essential to improve understanding and enable future projections16,17,18. Air temperature, wind speed, cloud cover, precipitation, water density, and lake transparency are all potential variables. We, here, compared the above variables to the turnover timing in Lake Biwa. The meteorological inputs in this study provided data for air temperature, wind speed, cloud cover, and precipitation. Water density and particulate organic carbon (POC) concentration representing lake transparency were the model’s outputs. The annual averages and cold season (November–April) values of the above variables were calculated over the simulation period (Supplementary Fig. S6). Annual averages illustrate general long-term warming trends18, while cold season values particularly determine the timing of turnover17. However, in Lake Biwa, air temperature during the cold season fluctuated greatly compared to the annual averages. A random forest analysis17 has been conducted between the turnover timing and the above two variable sets (cold season values versus annual averages) in Lake Biwa, and the cold season values better explained the turnover timing (35.39% versus 18.48%). The results agree with the conclusion drawn from the previous sensitivity tests, which indicated the relative importance of air temperature and solar radiation during winter based on 40 scenarios32.The importance of variables was estimated based on the random forest analysis using the cold season data (Fig. 4a). Wind speed dominates the timing of turnover, which is consistent with the previous studies17,25. The POC concentration, the difference in water density between the surface and bottom, and cloud cover have moderate effects on the timing of turnover. However, air temperature is less important, which is contrary to the turnover mechanism17,24,32. A re-confirmation was conducted of the relationship between turnover timing and air temperature (Fig. 4b and Supplementary Fig. S7). The cool air generally encourages an early turnover, albeit with several anomalies. The turnover timing between 1976 and 1990 remained constant independent of climate change, and the period coincidently had a substantial nutrient fluctuation (Fig. 3). As a result, it is essential to investigate the nutrient status further.Figure 4Analysis results of the relationship between potential variables and turnover timing: (a) the importance of variables importance using a random forest analysis, and (b) the relationship between the cold season air temperature and the timing of turnover. Variable importance is calculated using the percentage increase in mean square error (MSE) and the increase in node purity. Higher values illustrate the greater importance of the variable. Variables include air temperature (AT), precipitation (pptn.), cloud cover (CC), the difference in density (DD), POC, and wind speed (WS).Full size imageLake nutrient concentrationsBecause phosphorus is the limiting nutrient in Lake Biwa and DIP concentrations can be effectively limited by regulating external loadings as practiced (Fig. 3), DIP concentrations become the focus of this discussion for nutrient status. However, the DIP concentrations disproportionately responded to the external loadings of total phosphorus (TP) in Lake Biwa. Although external TP loading itself fails to determine lake phosphorus concentrations due to the hydrodynamics of lakes33, Lake Biwa exhibited insignificant changes in the inflow rate or the retention time (and see an example of the surface flow in Supplementary Fig. S8). Therefore, it can be assumed that the hydraulic loading remained constant, and the input nutrient concentrations were proportionate to the external nutrient loadings in Lake Biwa. This finding contradicts a recent meta-analysis that highlighted a deterministic relationship between input nutrient concentrations and lake nutrient concentrations, based on steady-state mass balance models6. The possible reason is the dynamics of the lake’s ecosystem22, which have been considered in this study. For example, the surface DIP concentrations were almost nonexistent regardless of the external TP loadings in Lake Biwa, supporting that phosphorus is the limiting nutrient in Lake Biwa34,35. The low DIP concentrations at the surface may be caused by the rapid recycling of phosphorus because the amount of phosphorus available for phytoplankton is easily affected by the feedback mechanism between phytoplankton photosynthesis and the phosphorus released from the water35,36.Hypoxia and strategiesThe variations in DO concentration are the public’s top concern as it relates to hypoxia, a key indicator of water quality. Lake bottom, among all water depths, is more sensitive to small changes in oxygen conditions12. In Lake Biwa, the annual minimum DO concentrations ranged from 2 to 5.5 mg/L over the last 60 years (Supplementary Fig. S9). The decrease in DO concentrations in the early period, typically till the 1980s, was mainly caused by nutrient enrichments (Fig. 3). The nutrient enrichment-induced heavy eutrophication eventually accelerates the rate of DO depletion2. After eutrophication was controlled in the 1980s, climate change became the dominant stressor23. There remains much uncertainty surrounding the relationship between climatic variables-related turnover timing and hypoxia in Lake Biwa12. We, therefore, first investigate the relationship between hypoxia and turnover timing, and then concentrate on nutrients to alleviate hypoxia.Although the relationship between turnover timing and DO concentrations is quite weak (R2 = 0.10), there is a general decrease in DO concentrations with increasing turnover timing (Fig. 5a). On the other hand, a linear relationship has been found between DIP concentrations and DO concentrations, with an R2 of 0.67 (Fig. 5b). The slope of –0.841 μgP/mgDO means an increase in DIP concentrations by approximately 0.841 μgP/L causes a decrease in DO concentrations by 1 mg/L. Note that the simulation results were compared over the whole period, and eutrophication-induced hypoxia differs theoretically from climate-induced hypoxia. Additional testing has been conducted to distinguish the effects of two stressors (eutrophication- and climate-induced hypoxia; Supplementary Fig. S10). Before 1980 when eutrophication progressed, the annual minimum DO concentrations and the DIP concentrations had a stronger linear relationship (R2 = 0.89). Although waste-water treatment has improved conditions in the lake, climate change induced alteration of turnover timing may adversely influence water quality. However, the relationship weakened dramatically with an R2 of 0.10 after 1980, when climate change dominated hypoxia. The lower R2 value indicates that climate-related hypoxia is more complex as concluded previously37,38. The two possibilities are as follows. First, there can be a legacy of hypoxia related to eutrophication. The DO recovery at the bottom of Lake Biwa was complicated by the low DO concentration in 1980 and the delayed timing of turnover; similar phenomena have been observed in the Lake of Zurich22. Second, ecosystem dynamics could help explain the difficulty in predicting hypoxia at the bottom. Phytoplankton fully exploits phosphorus at the surface, as explained above, then the death and sinking of the surface phytoplankton are accompanied by the sedimentation of phosphorus to the bottom as modeled. Bacteria break down the sinking phytoplankton, releasing phosphorus and consuming DO in the process. Additional DO consumption lowers the bottom DO concentration, which in turn encourages phosphorus release from the sediment in a low DO environment22. Such unfavorable feedback between DIP and DO concentrations are strengthened by prolonged stratification and eventually accelerates the development of hypoxia. However, future research is necessary because this numerical model simplified the relationship between water and sediment. The sinking of organic carbon into sediment is integrated in the model, and due to the decomposition of organic carbon in the sediment, nutrients are released into and oxygen is depleted in the water. Despite that, the trends between DO and DIP concentrations stay the same under climate change (Fig. 5b), and thus controlling lake phosphorus is beneficial to the Lake Biwa hypoxia.Figure 5The linear regression results of the relationship: (a) between turnover timing and annual minimum concentration of DO, (b) between the annual minimum concentration of DO and annual average concentration of DIP. The simulation results at the monitoring station were used for analysis.Full size image More

  • in

    Genomic and ecological evidence shed light on the recent demographic history of two related invasive insects

    Gandhi, K. J. K. & Herms, D. A. Direct and indirect effects of alien insect herbivores on ecological processes and interactions in forests of eastern North America. Biol. Invasions 12, 389–405 (2010).
    Google Scholar 
    Desurmont, G. A. et al. Alien interference: disruption of infochemical networks by invasive insect herbivores. Plant. Cell Environ. 37, 1854–1865 (2014).PubMed 

    Google Scholar 
    Kenis, M. et al. Ecological effects of invasive alien insects. Biol. Invasions 11, 21–45 (2009).
    Google Scholar 
    Paini, D. R. et al. Global threat to agriculture from invasive species. Proc. Natl. Acad. Sci. 113, 7575–7579 (2016).ADS 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Bradshaw, C. J. A. et al. Massive yet grossly underestimated global costs of invasive insects. Nat. Commun. 7, 1–8 (2016).
    Google Scholar 
    Sherpa, S. et al. Unravelling the invasion history of the Asian tiger mosquito in Europe. Mol. Ecol. 28, 2360–2377 (2019).PubMed 

    Google Scholar 
    Sherpa, S. et al. Landscape does matter: Disentangling founder effects from natural and human-aided post-introduction dispersal during an ongoing biological invasion. J. Anim. Ecol. 89, 2027–2042 (2020).PubMed 

    Google Scholar 
    Sherpa, S. & Després, L. The evolutionary dynamics of biological invasions: A multi‐approach perspective. Evol. Appl. (2021).North, H. L., McGaughran, A. & Jiggins, C. Insights into invasive species from whole-genome resequencing. Mol. Ecol. (2021).Ma, L. et al. Rapid and strong population genetic differentiation and genomic signatures of climatic adaptation in an invasive mealybug. Divers. Distrib. 26, 610–622 (2020).
    Google Scholar 
    Ortego, J., Céspedes, V., Millán, A. & Green, A. J. Genomic data support multiple introductions and explosive demographic expansions in a highly invasive aquatic insect. Mol. Ecol. 30, 4189–4203 (2021).PubMed 

    Google Scholar 
    Varone, L., Logarzo, G., Briano, J., Hight, S. & Carpenter, J. Cactoblastis cactorum (Berg) (Lepidoptera: Pyralidae) use of Opuntia host species in Argentina. Biol. Invasions 16, 2367–2380 (2014).
    Google Scholar 
    Singer, M. C., Ng, D. & Moore, R. A. Genetic variation in oviposition preference between butterfly populations. J. Insect Behav. 4, 531–535 (1991).
    Google Scholar 
    Forister, M. L. Oviposition preference and larval performance within a diverging lineage of lycaenid butterflies. Ecol. Entomol. 29, 264–272 (2004).
    Google Scholar 
    Wiklund, C. The concept of oligophagy and the natural habitats and host plants of Papilio machaon L. Fennoscandia. Insect Syst. Evol. 5, 151–160 (1974).
    Google Scholar 
    Courtney, S. P. & Forsberg, J. Host use by two pierid butterflies varies with host density. Funct. Ecol. 2, 67–75 (1988).
    Google Scholar 
    Franklin, J. Species distribution models in conservation biogeography: developments and challenges. Divers. Distrib. 19, 1217–1223 (2013).
    Google Scholar 
    Peterson, A. et al. Ecological niches and geographic distributions. Monographs in Population Biology vol. 49 (2011).Alvarado-Serrano, D. F. & Knowles, L. L. Ecological niche models in phylogeographic studies: Applications, advances and precautions. Mol. Ecol. Resour. 14, 233–248 (2014).PubMed 

    Google Scholar 
    Carrera-Martínez, R., Aponte-Díaz, L. A., Ruiz-Arocho, J., Lorenzo-Ramos, A. & Jenkins, D. A. The effects of the invasive Harrisia cactus mealybug (Hypogeococcus sp.) and exotic lianas (Jasminum fluminense) on Puerto Rican native cacti survival and reproduction. Biol. Invasions 21, 3269–3284 (2019).
    Google Scholar 
    Acevedo-Rodríguez, P. & Strong, M. T. Catalogue of seed plants of the West Indies. Smithson. Contrib. to Bot. 98, 1–1192 (2012).
    Google Scholar 
    Carrera-Martínez, R., Aponte-Díaz, L., Ruiz-Arocho, J. & Jenkins, D. A. Symptomatology of infestation by Hypogeococcus pungens: Contrasts between host species. Haseltonia 2015, 14–18 (2015).
    Google Scholar 
    Aponte-Díaz, L., Ruiz-Arocho, J., Carrera-Martínez, R. & Ee, B. Contrasting effects of the invasive Hypogeococcus sp. (Hemiptera: Pseudococcidae) infestation on seed germination of Pilosocereus royenii (Cactaceae), a Puerto Rican native cactus. Caribb. J. Sci. 50, 212–218 (2020).
    Google Scholar 
    California Department of Food and Agriculture. Harrisia Cactus Mealybug | Hypogeococcus pungens | Pest rating proposals and final ratings. https://blogs.cdfa.ca.gov/Section3162/?p=5881 (2018).Poveda-Martínez, D. et al. Species complex diversification by host plant use in an herbivorous insect: The source of Puerto Rican cactus mealybug pest and implications for biological control. Ecol. Evol. 10, 10463–10480 (2020).PubMed 
    PubMed Central 

    Google Scholar 
    Segarra-Carmona, A. E., Ramírez-Lluch, A., Cabrera-Asencio, I. & Jiménez-López, A. N. First report of a new invasive mealybug, the Harrisia cactus mealybug Hypogeococcus pungens (Hemiptera: Pseudococcidae). J. Agric. Univ. Puerto Rico 94, 183–187 (2010).
    Google Scholar 
    Poveda-Martínez, D. et al. Untangling the Hypogeococcus pungens species complex (Hemiptera: Pseudococcidae) for Argentina, Australia, and Puerto Rico based on host plant associations and genetic evidence. PLoS ONE 14, e0220366 (2019).PubMed 
    PubMed Central 

    Google Scholar 
    McKenzie, H. L. Mealybugs of California. (Univ of California Press, 1967).Hamon, A. B. A cactus mealybug, Hypogeococcus festerianus (Lizer y Trelles). Florida (Homoptera Coccoidea Pseudococcidae). Entomol. Circ. Div. Plant Ind. Florida Dep. Agric. Consum. Serv. 263, 2 (1984).Hodges, A. & Hodges, G. Hypogeococcus pungens Granara de Willink (Insecta: Hemiptera: Pseudococcidae), a mealybug. EDIS 2009, (2009).Halbert, S. Entomology section. Triology 35, 2–4 (1996).
    Google Scholar 
    Aguirre, M. B. et al. Analysis of biological traits of Anagyrus cachamai and Anagyrus lapachosus to assess their potential as biological control candidate agents against Harrisia cactus mealybug pest in Puerto Rico. Biocontrol 64, 539–551 (2019).CAS 

    Google Scholar 
    Aguirre, M. B. et al. Influence of competition and intraguild predation between two candidate biocontrol parasitoids on their potential impact against Harrisia cactus mealybug, Hypogeococcus sp. (Hemiptera: Pseudococcidae). Sci. Rep. 11, 13377 (2021).ADS 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Eaton, D. A. R. & Overcast, I. ipyrad: Interactive assembly and analysis of RADseq datasets. Bioinformatics 36, 2592–2594 (2020).CAS 
    PubMed 

    Google Scholar 
    Poveda-Martínez, D., Salinas, N., Aguirre, M. B., Sánchez-Restrepo, A. F. & Hight, S., Diaz-Soltero, H. Dataset generated in Genomic and ecological evidence shed light on the recent demographic history of two related invasive insects. https://doi.org/10.6084/m9.figshare.15167082.v2 (2022).Frichot, E., Mathieu, F., Trouillon, T., Bouchard, G. & François, O. Fast and efficient estimation of individual ancestry coefficients. Genetics 196, 973–983 (2014).PubMed 
    PubMed Central 

    Google Scholar 
    Excoffier, L., Dupanloup, I., Huerta-Sánchez, E., Sousa, V. C. & Foll, M. Robust demographic inference from genomic and SNP data. PLoS Genet. 9, e1003905 (2013).PubMed 
    PubMed Central 

    Google Scholar 
    Gattepaille, L. M., Jakobsson, M. & Blum, M. G. B. Inferring population size changes with sequence and SNP data: Lessons from human bottlenecks. Heredity (Edinb). 110, 409–419 (2013).CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Born‐Schmidt, G. et al. The implementation of the mexican strategy on invasive species: How far have we come? Invasive Alien Species Obs. Issues from Around World 4, 153–164 (2021).
    Google Scholar 
    McFadyen, R. E. & Tomley, A. J. Preliminary indications of success in the biological control of Harrisia cactus (Eriocereus martinii Lab.) in Queensland. In Proceedings of the First Conference of the Council of Australian Weed Science Societies held at National Science Centre, Parkville, Victoria, Australia, 12–14 April 1978 108–112 (Council of Australian Weed Science Societies, 1978).McFadyen, R. E. & Tomley, A. J. The successful biological control of Harrisia cactus (Eriocereus martinii) in Queensland. In Proceedings of the Sixth Australian Weeds Conference, Volume 1, City of Gold Coast, Queensland, Australia, 13–18 September, 1981 139–143 (Queensland Weed Society, 1981).Paterson, I. D. et al. Biological control of Cactaceae in South Africa. African Entomol. 19, 230–246 (2011).
    Google Scholar 
    Sutton, G. F., Klein, H. & Paterson, I. D. Evaluating the efficacy of Hypogeococcus sp. as a biological control agent of the cactaceous weed Cereus jamacaru in South Africa. Biocontrol 63, 493–503 (2018).
    Google Scholar 
    Paterson, I. D. et al. Biological control of Cactaceae in South Africa. African Entomol. 29, 713–734 (2021).
    Google Scholar 
    McFadyen, R. E. Harrisia (Eriocereus) martinii (Labour.) Britton—Harrisia cactus Acanthocereus tetragonus (L.) Hummelink—sword pear. (ed. Julien, M., McFadyen, R., & Cullen, J.), Biological control of weeds in Australia 274– 281. (CSIRO Publishing, 2012).Julien, M. H. & Griffiths, M. Biological control of weeds: A world catalogue of agents and their target weeds. (Cab International, 1998).Houston, W. A. & Elder, R. Biocontrol of Harrisia cactus Harrisia martinii by the mealybug Hypogeococcus festerianus (Hemiptera: Pseudococcidae) in salt-influenced habitats in Australia. Austral Entomol. 58, 696–703 (2019).
    Google Scholar 
    Hofmeister, N., Werner, S. & Lovette, I. Environmental correlates of genetic variation in the invasive European starling in North America. Mol. Ecol. 30, 1251–1263 (2021).PubMed 

    Google Scholar 
    Driscoe, A. L. et al. Host plant associations and geography interact to shape diversification in a specialist insect herbivore. Mol. Ecol. 28, 4197–4211 (2019).CAS 
    PubMed 

    Google Scholar 
    Vidal, M. C., Quinn, T. W., Stireman, J. O. 3rd., Tinghitella, R. M. & Murphy, S. M. Geography is more important than host plant use for the population genetic structure of a generalist insect herbivore. Mol. Ecol. 28, 4317–4334 (2019).PubMed 

    Google Scholar 
    Poveda-Martínez, D. et al. Spatial and host related genomic variation in partially sympatric cactophagous moth species. Mol. Ecol. 31, 356–371 (2021).PubMed 

    Google Scholar 
    Cao, L., Wei, S., Hoffmann, A. A., Wen, J. & Chen, M. Rapid genetic structuring of populations of the invasive fall webworm in relation to spatial expansion and control campaigns. Divers. Distrib. 22, 1276–1287 (2016).
    Google Scholar 
    Sih, A. et al. Predator–prey naïveté, antipredator behavior, and the ecology of predator invasions. Oikos 119, 610–621 (2010).
    Google Scholar 
    Yang, Q.-Q. et al. Introgressive hybridization between two non-native apple snails in China: Widespread hybridization and homogenization in egg morphology. Pest Manag. Sci. 76, 4231–4239 (2020).CAS 
    PubMed 

    Google Scholar 
    Cordeiro, E. M. G. et al. Hybridization and introgression between Helicoverpa armigera and H zea: An adaptational bridge. BMC Evol. Biol. 20, 61 (2020).CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Pardo-Diaz, C. et al. Adaptive introgression across species boundaries in Heliconius butterflies. PLOS Genet. 8, e1002752 (2012).CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Caltagirone, L. E. Landmark examples in classical biological control. Annu. Rev. Entomol. 26, 213–232 (1981).
    Google Scholar 
    Goldson, S. L., Phillips, C. B. & Barlow, N. D. The value of parasitoids in biological control. New Zeal. J. Zool. 21, 91–96 (1994).
    Google Scholar 
    Wang, Z., Liu, Y., Shi, M., Huang, J. & Chen, X. Parasitoid wasps as effective biological control agents. J. Integr. Agric. 18, 705–715 (2019).
    Google Scholar 
    Miller, G., & Lugo. A. E. Guide to the ecological systems of Puerto Rico. IITF-GTR-35 (2009).Bolger, A. M., Lohse, M. & Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120 (2014).CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Andrews, S. FastQC: A Quality control tool for high throughput sequence data. (2010).Ewels, P., Magnusson, M., Lundin, S. & Käller, M. MultiQC: summarize analysis results for multiple tools and samples in a single report. Bioinformatics 32, 3047–3048 (2016).CAS 
    PubMed 
    PubMed Central 

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

    Google Scholar 
    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).CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Foll, M. & Gaggiotti, O. A genome-scan method to identify selected loci appropriate for both dominant and codominant markers: A Bayesian perspective. Genetics 180, 977–993 (2008).PubMed 
    PubMed Central 

    Google Scholar 
    Linck, E. & Battey, C. J. Minor allele frequency thresholds strongly affect population structure inference with genomic data sets. Mol. Ecol. Resour. 19, 639–647 (2019).CAS 
    PubMed 

    Google Scholar 
    Keenan, K., McGinnity, P., Cross, T. F., Crozier, W. W. & Prodöhl, P. A. diveRsity: An R package for the estimation and exploration of population genetics parameters and their associated errors. Methods Ecol. Evol. 4, 782–788 (2013).
    Google Scholar 
    Goudet, J. Hierfstat, a package for R to compute and test hierarchical F-statistics. Mol. Ecol. Notes 5, 184–186 (2005).
    Google Scholar 
    Kamvar, Z. N., Tabima, J. F. & Grünwald, N. J. Poppr: An R package for genetic analysis of populations with clonal, partially clonal, and/or sexual reproduction. PeerJ 2, e281 (2014).PubMed 
    PubMed Central 

    Google Scholar 
    Tajima, F. The effect of change in population size on DNA polymorphism. Genetics 123, 597–601 (1989).CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Frichot, E. & François, O. LEA: An R package for landscape and ecological association studies. Methods Ecol. Evol. 6, 925–929 (2015).
    Google Scholar 
    Pembleton, L. W., Cogan, N. O. I. & Forster, J. W. St AMPP: An R package for calculation of genetic differentiation and structure of mixed-ploidy level populations. Mol. Ecol. Resour. 13, 946–952 (2013).CAS 
    PubMed 

    Google Scholar 
    Cockerham, C. C. Drift and mutation with a finite number of allelic states. Proc. Natl. Acad. Sci. 81, 530–534 (1984).ADS 
    CAS 
    PubMed 
    PubMed Central 
    MATH 

    Google Scholar 
    Lynch, M. & Conery, J. S. The origins of genome complexity. Science 302, 1401–1404 (2003).ADS 
    CAS 
    PubMed 

    Google Scholar 
    Rozas, J. et al. DnaSP 6: DNA sequence polymorphism analysis of large data sets. Mol. Biol. Evol. 34, 3299–3302 (2017).CAS 
    PubMed 

    Google Scholar 
    Keightley, P. D., Ness, R. W., Halligan, D. L. & Haddrill, P. R. Estimation of the spontaneous mutation rate per nucleotide site in a Drosophila melanogaster full-sib family. Genetics 196, 313–320 (2014).CAS 
    PubMed 

    Google Scholar 
    R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria (2021).Neu, C. W., Byers, C. R. & Peek, J. M. A technique for analysis of utilization-availability data. J. Wildl. Manage. 38, 541–545 (1974).
    Google Scholar 
    Soberón, J. & Peterson, A. Interpretation of models of fundamental ecological niches and species’ distributional areas. Biodivers. Informatics 2, 1-10 (2005).Jorge, S. & Miguel, N. Niches and distributional areas: Concepts, methods, and assumptions. Proc. Natl. Acad. Sci. 106, 19644–19650 (2009).
    Google Scholar 
    Phillips, S. J., Anderson, R. P. & Schapire, R. E. Maximum entropy modeling of species geographic distributions. Ecol. Modell. 190, 231–259 (2006).
    Google Scholar 
    Cobos, M. E., Peterson, A., Barve, N. & Osorio-Olvera, L. Kuenm: An R package for detailed development of ecological niche models using Maxent. PeerJ 7, e6281 (2019).PubMed 
    PubMed Central 

    Google Scholar 
    Fick, S. E. & Hijmans, R. J. WorldClim 2: New 1-km spatial resolution climate surfaces for global land areas. Int. J. Climatol. 37, 4302–4315 (2017).
    Google Scholar 
    Title, P. O. & Bemmels, J. B. ENVIREM: An expanded set of bioclimatic and topographic variables increases flexibility and improves performance of ecological niche modeling. Ecography (Cop.) 41, 291–307 (2018).
    Google Scholar 
    Warren, B. H. et al. Evaluating alternative explanations for an association of extinction risk and evolutionary uniqueness in multiple insular lineages. Evolution 72, 2005–2024 (2018).PubMed 

    Google Scholar 
    Warren, D. L., Glor, R. E. & Turelli, M. Environmental niche equivalency versus conservatism: Quantitative approaches to niche evolution. Evol. Int. J. Org. Evol. 62, 2868–2883 (2008).
    Google Scholar 
    Schoener, T. W. The anolis lizards of Bimini: Resource partitioning in a complex fauna. Ecology 49, 704–726 (1968).
    Google Scholar 
    Van der Vaart, A. W. Asymptotic Statistics (UK Cam, 1998).MATH 

    Google Scholar  More

  • in

    The study of aggression and affiliation motifs in bottlenose dolphins’ social networks

    Subjects and facilityWe observed two groups of Atlantic bottlenose dolphins (six different individuals in total) housed at the marine zoo “Marineland Mallorca”. One of the groups was composed of four individuals (G1) and the other was constituted by five individuals (G2). The two adult males and one of the females were the same in both groups (Table 1). Group composition changed due to the transfer of individuals to another pool of the zoo and due to the arrival of new individuals from another aquatic park.Table 1 Age, sex, group, and identification number in the network of the subject dolphins. M male, F female.Full size tableThe dolphins were kept in three outdoor interconnecting pools: the main performance pool (1.6 million liters of water), a medical pool (37.8 thousand liters of water) and a small pool (636.8 thousand liters of water). During the observational periods, the dolphins had free access to all the pools. Underwater viewing at the main and the small pool was available through the transparent walls around the rim of the pools.Ethics statementThis study was approved by the UIB Committee of Research Ethics and Marineland Mallorca. This research was conducted in compliance with the standards of the European Association of Zoos and Aquaria (EAZA). All subjects tested in this study were housed in Marineland Mallorca following the Directive 1999/22/EC on the keeping of animals in zoos. This study was strictly non-invasive and did not affect the welfare of dolphins.Behavioral observations and data collectionBehavioral data were collected in situ by APM from May to November 2016 for G1 and from November 2017 to February 2018 for G2. All observational periods were also recorded using two waterproof cameras SJCAM SJ4000. Observations were conducted at the main pool between 8:00 a.m. and 11:00 a.m. Due to the schedules and dynamics of the zoo, we were unable to collect data outside this period. Dolphin social behavior was registered and videotaped for 30 min–2 h each day. Only data from sessions that lasted at least 30 min were included in the analysis. We did not collect any data during training or medical procedures and resumed the observational session a few minutes after the end of these events.We recorded all occurrences of affiliative and aggressive interactions, the identities of the involved individuals and the identity of the dolphin initiating the contact. Aggressive contacts were defined by the occurrence of chasing, biting, and hitting, as established in previous studies37,38,39,40,41. Affiliative contacts were defined as contact swimming, synchronous breathing and swimming (at least 30″ of continuous swimming) or flipper-rubbing, as established in previous studies37,39,40,41,43.To assess the strength of the affiliative bonds in both groups, we calculated the index of affiliative relationships (IA) between dolphins following the procedure described in Yamamoto et al. For calculating the IA we recorded the relative frequencies of synchronous swimming since it is a well-defined affiliative behavior in dolphins. Data of synchronous swimming were recorded using group 0–1 sampling44 at 3-min intervals. This method consists of the observation of individuals during short periods and the recording of the occurrence (assigning to that period a 1) or non-occurrence (assigning to that period a 0) of a well-defined behavior44. For calculating the IA for each couple, the number of sampling periods in which synchronous swimming between individuals A and B occurred (XAB) was divided by the number of sampling periods in which individuals A and B were observed (YAB): (IA=frac{{X}_{AB}}{{Y}_{AB}})39,45. Therefore, the IA reflects the level of affiliation for each dolphin dyad based on the pattern of synchronous swimming. This index served to construct the general affiliative social networks of both groups of dolphins.Temporal network constructionTemporal networks can provide insight into social events such as conflicts and post-conflict interactions in which the order of interactions and the timing is crucial. Furthermore, they allow us to calculate the probabilities of the different affiliative and aggressive interactions occurring in the group.We used behavioral observations to construct temporal networks for each group. Each dolphin was treated as a node (N) with their aggressive and affiliative interactions supplying the network links. We divided the daily observations into periods of 3 min. In each period, we assigned a positive (+ 1), negative (− 1) or neutral (0) interaction to each pair of dolphins. That is, if during the period a pair of dolphins displayed affiliative interactions, we assigned a + 1 to the link between that pair of nodes, if they were involved in a conflict, we assigned a − 1, and if the pair did not engage in any interaction, we assigned to that link a 0. If during the same period, the pair displayed both aggressive and affiliative interactions we considered the last observed interaction. Therefore, we obtained an adjacency matrix (an N × N matrix describing the links in the network) for each group of dolphins. Thus, for each day we had a series of different signed networks of the group, each network representing a 3-min period.Social network analysis: time-aggregated networks and network motifsWe collapsed the temporal networks of each day in time-aggregated networks. This procedure consists in aggregating the data collected over time within specific intervals to create weighted networks. The sign and the weight of the links characterize these networks, indicating the valence and duration of the interaction respectively. Thus, they are static representations of the social structure of the group of dolphins. To obtain these time-aggregated networks we proceeded as follows:First, for each day we aggregated the values of each interaction of the temporal networks until one link qualitatively changed. We considered a qualitative change if one interaction passed from being negative (− 1) to positive (+ 1) meaning that the pair of dolphins reconciled after the conflict or vice versa, or if a new affiliation (+ 1) or aggression (− 1) took place, that is the link changed from being neutral (0) to positive or negative. If a link changed from being negative or positive to being neutral, we did not consider that this interaction has changed qualitatively. For example, if dolphins interacted positively during two periods of time, then they ceased to interact (neutral) and finally they engaged in an aggressive interaction, the total weight of the interaction in the resulting time-aggregated network would be of + 2. Therefore, a conflict or an affiliation may extend over multiple periods containing several contacts, and is considered finished when the interaction changes its valence. In this way, we obtained a series of time-aggregated networks for each day, which retain the information on the duration, timing, and ordering of the affiliative and aggressive events in the group.We examined the local-scale structure of the affiliative-aggressive social networks using motif analysis. Thus, for each group, we analyzed the network motif representation of the temporal and time-aggregated networks, identifying and recording the number of occurrences of each motif.Model of affiliative and aggressive interactionsWe built two models (a simple and a complex one) that aim to simulate the dynamics of aggressive and affiliative interactions of a group of four dolphins. These models were created using the observed probabilities of each affiliative or aggressive interaction between individuals in group G1. We only used the data of G1 since we had more hours of video recordings and, thus, more statistics of the pattern of dolphins’ interactions. Both models return affiliative/aggressive temporal networks constituted by four nodes and different aggressive, affiliative, or neutral interactions between the six possible pairs of individuals in the network. We simulated data for 20 periods of 3 min per day for a total of 80 days to mimic the empirical data time structure. We obtained one temporal network for each period (1600 temporal networks in total) and ran 100 realizations of each model.Our models work as follows: At the beginning of the simulations, all the interactions between the four nodes are neutral (0). In each period, we select a pair of nodes randomly and assign to that link a positive (+ 1) or a negative (− 1) interaction with probability p (calculated previously for each type of interaction). These interactions correspond to spontaneous aggressions and affiliations. In the complex model, if in the previous period a conflict took place, before assessing spontaneous interactions we first evaluated the different possible post-conflict contacts that could occur (reconciliation, new aggressions, and affiliations). Therefore, for reconciliations, we change the valence of the interaction from negative to positive with a certain probability. Then, we also randomly choose a pair of nodes including one of the former opponents and assign to that link a positive or negative interaction with the observed probabilities to simulate the occurrence of new affiliations (third party-affiliation) or redirected aggressions arising from the previous conflict. We keep on doing this procedure period by period. Lastly, we obtained the time-aggregated networks for the two models.The simpler model only includes the probability of aggression and affiliation between group members, whereas the complex one also includes the patterns of conflict resolution previously observed. In this way, the complex model serves to assess the influence of post-conflict management mechanisms on the observed pattern of aggressive/affiliative networks. That is, the complex model also keeps track of past actions. Thus, depending on the interaction of the previous step, the probability of the following interaction changes based on the observed pattern of conflict resolution strategies.Calculation of the observed probabilities of affiliative and aggressive interactionsFor the simple model, we calculated the probability of general aggression and affiliation per day without distinguishing between types of positive and negative interactions. Thus, we obtained the number of periods in which an aggressive or affiliative contact took place per day and divided it by the total number of periods of that day (probability of general aggression or affiliation per 3-min period). With these probabilities, we calculated the mean probability of general aggression and affiliation per period.For the complex model, we calculated the probabilities of reconciliation, new affiliations/aggressions, and spontaneous affiliations/aggressions per day. That is, the probability that former opponents exchange affiliative contacts after an aggressive encounter (reconciliation), the probabilities that a conflict may promote new affiliations (third-party affiliation) or new conflicts (redirected aggression) between one of the opponents and a bystander in the same day, and the probability of affiliative or aggressive interactions not derived from a previous conflict (spontaneous interactions). To classify affiliations and aggressions in these categories we used the temporal networks, examining the interactions that took place after a conflict between opponents and between them and bystanders. If the opponents reconciled or affiliated with a bystander after a fight, we assumed that the following affiliative or aggressive interactions were spontaneous and were not a consequence of that conflict. Thus, to calculate the number of spontaneous affiliations, we subtracted the number of reconciliations and new affiliations from the total number of affiliations per day. For spontaneous aggressions, we subtracted the number of new aggressions to the total number of aggressions per day. Then, we obtained the probability of spontaneous affiliation and aggression per period.Using the previous probabilities, we obtained the rate (r) of reconciliation, new aggression and new affiliation per minute with the following formula:({p=1-e}^{-rDelta t}). Using the same formula, we finally calculated the probability of reconciliation, new aggression and affiliation per 3-min period used in the complex model (Supplementary Table 1 for details of probabilities calculation).Network-motif analysisWe also carried out a network-motif analysis. As we did not consider the identities or sex of the nodes in these models, we grouped the obtained motifs into equivalent categories considering the pattern of interactions between nodes. We also classified the motifs obtained from the real data of G1 into those equivalent categories. Finally, we compared the pattern of equivalent network motifs of the observed social network of dolphins and the ones of the two models. To do so we calculated the Spearman’s rank correlation coefficient (rs), defined as a nonparametric measure of the statistical dependence between the rankings of two variables: ({r}_{s}=frac{covleft({rg}_{X}{rg}_{Y}right)}{{sigma }_{{rg}_{X}}}{sigma }_{{rg}_{Y}}); rgX and rgY are the rank variables; cov (rgX rgY) is the covariance of the rank variables, and σrgX and σrgY are the standard deviations of the rank variables. Therefore, this coefficient allows us to assess the statistical dependence between the motif ranking of the real data and the one of each model.Computational implementationsAll the models, network construction, visualization and motif analysis were generated and implemented using MATLAB R2018b. More

  • in

    High-yield dairy cattle breeds improve farmer incomes, curtail greenhouse gas emissions and reduce dairy import dependency in Tanzania

    Meat, Milk and More: Policy Innovations to Shepherd Inclusive and Sustainable Livestock Systems in Africa (Malabo Montpellier Panel, 2020).Value of Agricultural Production (FAO, accessed August 25, 2022); https://www.fao.org/faostat/en/#data/QVJayne, T. & Sanchez, P. A. Agricultural productivity must improve in sub-Saharan Africa. Science 372, 1045–1047 (2021).Article 
    ADS 
    CAS 
    PubMed 

    Google Scholar 
    Dangal, S. R. S. et al. Methane emission from global livestock sector during 1890–2014: magnitude, trends and spatiotemporal patterns. Glob. Change Biol. 23, 4147–4161 (2017).Article 
    ADS 

    Google Scholar 
    Mottet, A. et al. Climate change mitigation and productivity gains in livestock supply chains: insights from regional case studies. Reg. Env. Change 17, 129–141 (2016).Article 

    Google Scholar 
    Valin, H. et al. Agricultural productivity and greenhouse gas emissions: trade-offs or synergies between mitigation and food security? Environ. Res. Lett. 8, 035019 (2013).Article 
    ADS 
    CAS 

    Google Scholar 
    González-Quintero, R. et al. Yield gap analysis to identify attainable milk and meat productivities and the potential for greenhouse gas emissions mitigation in cattle systems of Colombia. Agric. Syst. 195, 103303 (2022).Article 

    Google Scholar 
    Crops and Livestock Products (FAO, accessed August 17,2022); https://www.fao.org/faostat/en/#data/QCLLedo, J. et al. Persistent challenges in safety and hygiene control practices in emerging dairy chains: the case of Tanzania. Food Control 105, 164–173 (2019).Article 

    Google Scholar 
    Häsler, B. et al. Integrated food safety and nutrition assessments in the dairy cattle value chain in Tanzania. Glob. Food Sec. 18, 102–113 (2018).Article 

    Google Scholar 
    Supply Utilization Accounts (FAO, accessed August 26, 2022); https://www.fao.org/faostat/en/#data/SCLMichael, S. et al. Tanzania Livestock Master Plan (International Livestock Research Institute, 2018).Tanzania Livestock Sector Analysis (2016/2017–2030/2031) (United Republic of Tanzania Ministry of Livestock and Fisheries, 2017); https://www.mifugouvuvi.go.tz/uploads/projects/1553602287-LIVESTOCK%20SECTOR%20ANALYSIS.pdfNicholson, C. et al. Assessment of Investment Priorities for Tanzania’s Dairy Sector: Report on Activities and Accomplishments (International Livestock Research Institute, 2021).Chagunda, M. G. C., Romer, D. A. M. & Roberts, D. J. Effect of genotype and feeding regime on enteric methane, non-milk nitrogen and performance of dairy cows during the winter feeding period. Livest. Sci. 122, 323–332 (2009).Article 

    Google Scholar 
    Notenbaert, A. et al. Towards environmentally sound intensification pathways for dairy development in the Tanga region of Tanzania. Reg. Environ. Change 20, 138 (2020).Yesuf, G. A. et al. Embedding stakeholders’ priorities into the low-emission development of the East African dairy sector. Env. Res. Lett. 16, 064032 (2021).Article 
    CAS 

    Google Scholar 
    GLS (Greening Livestock Survey) (International Livestock Research Institute, 2019); https://data.ilri.org/portal/dataset/greeninglivestockIntended Nationally Determined Contributions (United Republic of Tanzania, 2021); https://unfccc.int/sites/default/files/NDC/2022-06/TANZANIA_NDC_SUBMISSION_30%20JULY%202021.pdfNdung’u, P. W. et al. Farm-level emission intensities of smallholder cattle (Bos indicus; B. indicus–B. taurus crosses) production systems in highlands and semi-arid regions. Animal 16, 100445 (2022).Article 
    PubMed 

    Google Scholar 
    Goopy, J. P. et al. Severe below-maintenance feed intake increases methane yield from enteric fermentation in cattle. Br. J. Nutr. 123, 1239–1246 (2020).Article 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Goopy, J. P. et al. A new approach for improving emission factors for enteric methane emissions of cattle in smallholder systems of East Africa—results for Nyando, Western Kenya. Agric. Syst. 161, 72–80 (2018).Article 

    Google Scholar 
    Supporting Low Emissions Development in the Tanzanian Dairy Cattle Sector—Reducing Enteric Methane for Food Security and Livelihoods (FAO, 2019).Gerssen-Gondelach, S. J. et al. Intensification pathways for beef and dairy cattle production systems: impacts on GHG emissions, land occupation and land use change. Agric. Ecosyst. Environ. 240, 135–147 (2017).Article 

    Google Scholar 
    Havlik, P. et al. Climate change mitigation through livestock system transitions. Proc. Natl Acad. Sci. USA 111, 3709–3714 (2014).Article 
    ADS 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Herrero, M. et al. Greenhouse gas mitigation potentials in the livestock sector. Nat. Clim. Change 6, 452–461 (2016).Article 
    ADS 

    Google Scholar 
    Dizyee, K., Baker, D. & Omore, A. Upgrading the smallholder dairy value chain: a system dynamics ex-ante impact assessment in Tanzania’s Kilosa district. J. Dairy Res. 86, 440–449 (2019).Article 
    CAS 
    PubMed 

    Google Scholar 
    Simões, A. R. P., Nicholson, C. F., Novakovicc, A. M. & Protil, R. M. Dynamic impacts of farm-level technology adoption on the Brazilian dairy supply chain. Int. Food Agribus. Manag. Rev. 23, 71–84 (2020).Article 

    Google Scholar 
    Rahimi, J. et al. Heat stress will detrimentally impact future livestock production in East Africa. Nat. Food. 2, 88–96 (2021).Article 

    Google Scholar 
    Mbululo, Y. & Nyihirani, F. Climate characteristics over southern highlands Tanzania. Atmos. Clim. Sci. 2, 454–463 (2012).
    Google Scholar 
    Kihoro, E. M., Schoneveld, G. C. & Crane, T. A. Pathways toward inclusive low-emission dairy development in Tanzania: producer heterogeneity and implications for intervention design. Agric. Syst. 190, 103073 (2021).Mruttu, H. et al. Animal Genetics Strategy and Vision for Tanzania (Tanzania Ministry of Agriculture, Livestock and Fisheries and ILRI, 2016).Agricultural Sample Survey 2018/19 Report on Livestock and Livestock Characteristics (Private Peasant Holdings) (Central Statistical Agency, 2019).2019/20 National Sample Census of Agriculture Main Report (Tanzania National Bureau of Statistics, 2022).Robinson, T. P. et al. Global Livestock Production Systems (FAO, 2011).Herrero, M. et al. Biomass use, production, feed efficiencies and greenhouse gas emissions from global livestock systems. Proc. Natl Acad. Sci. USA 110, 20888–20893 (2013).Article 
    ADS 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Baseline Study of the Tanzania Dairy Value Chain (United Republic of Tanzania Ministry of Agriculture, Livestock and Fisheries, 2016).Mbwambo, N., Nandonde, S., Ndomba, C. & Desta, S. Assessment of Animal Feed Resources in Tanzania (Tanzania Ministry of Agriculture, Livestock and Fisheries and ILRI, 2016).Hartung, C., Lerer, A., Anokwa, Y., Tseng, C., Brunette, W., & Borriello, G. Open data kit: tools to build information services for developing regions. Proc. 4th ACM/IEEE International Conference on Information and Communication Technologies and Development (Association for Computing Machinery, 2010).R Core Team. R: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, 2022).https://www.r-project.orgRufino, M. C. et al. Lifetime productivity of dairy cows in smallholder farming systems of the central highlands of Kenya. Animal 3, 1044–1056 (2009).Article 
    CAS 
    PubMed 

    Google Scholar 
    Hawkins, J. et al. Feeding efficiency gains can increase the greenhouse gas mitigation potential of the Tanzanian dairy sector. Sci. Rep. 11, 4190 (2021).Article 
    ADS 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Python Software Foundation (Python Software Foundation, 2019); https://www.python.org/psf/Kashoma, I. P. B. et al. Predicting body weight of Tanzania shorthorn zebu cattle using heart girth measurements. Livest. Res. Rural. Dev. 23, Table 1 (2011).Galukande, E. B., Mahadevan, P. & Black, J. G. Milk production in East African zebu cattle. Anim. Sci. 4, 329–336 (1962).Article 

    Google Scholar 
    Gillah, K. A., Kifaro, G. C. & Madsen, J. Effects of pre partum supplementation on milk yield, reproduction and milk quality of crossbred dairy cows raised in a peri urban farm of Morogoro town Tanzania. Livest. Res. Rural. Dev. 26 (2014).Njau, F. B. C., Lwelamira, J. & Hyandye, C. Ruminant livestock production and quality of pastures in the communal grazing land of semi-arid central Tanzania. Livest. Res. Rural. Dev. 8, Table 4 (2013).Mwambene, P. L. et al. Selecting indigenous cattle populations for improving dairy production in the Southern Highlands and Eastern Tanzania. Livest. Res. Rural. Dev. 26 (2014).Rege, J. E. O. et al. Cattle of Kenya: Uses, Performance, Farmer Preferences, Measures of Genetic Diversity and Options for Improved Use (International Livestock Research Institute, 2001).Beffa, L. M. Genotype × Environment Interaction in Afrikaner Cattle. PhD thesis, Univ. of the Free State (2005).Meaker, H. J., Coetsee, T. P. N. & Lishman, A. W. The effects of age at 1st calving on the productive and reproductive-performance of beef-cows. S. Afr. J. Anim. Sci. 10, 105–113 (1980).
    Google Scholar 
    Chenyambuga, S. W. & Mseleko, K. F. Reproductive and lactation performances of Ayrshire and Boran crossbred cattle kept in smallholder farms in Mufindi district, Tanzania. Livest. Res. Rural. Dev. 21, 100 (2009).
    Google Scholar 
    Ojango, J. M. K. et al. Dairy production systems and the adoption of genetic and breeding technologies in Tanzania, Kenya, India and Nicaragua. Anim. Genet. Resour. 59, 81–95 (2016).Article 

    Google Scholar 
    Feedipedia—Animal Feed Resources Information System (FAO, accessed 2021); https://www.feedipedia.org/Lukuyu, B. et al. (eds) Feeding Dairy Cattle in East Africa (East Africa Dairy Development Project, 2012).Rubanza, C. D. K. et al. Biomass production and nutritive potential of conserved forages in silvopastoral traditional fodder banks (Ngitiri) of Meatu District of Tanzania. Asian-Aust. J. Anim. Sci. 19, 978–983 (2006).Article 

    Google Scholar 
    Food Balances (2010-) (FAO, accessed September 29, 2021); http://www.fao.org/faostat/en/#data/FBSCrop Data for the United Republic of Tanzania (FAO, accessed September 22, 2021); http://www.fao.org/faost at/en/#data/QCGilbert, M. et al. Global distribution data for cattle, buffaloes, horses, sheep, goats, pigs, chickens and ducks in 2010. Sci. Data. 5, 180227 (2018).Article 
    PubMed 
    PubMed Central 

    Google Scholar 
    2014/15 Annual Agricultural Sample Survey Report (The United Republic of Tanzania, 2016).Basic Data for Livestock and Fisheries (The United Republic of Tanzania Ministry of Livestock and Fisheries, 2013).IPCC Guidelines for National Greenhouse Gas Inventories Vol. 4 Agriculture, Forestry and Other Land Use (IPCC, 2006).2019 Refinement to the IPCC Climate Change 2014: Synthesis Report (eds Core Writing Team, Pachauri, R. K. & Meyer L. A.) (IPCC, 2019).Fertilizers by Nutrient (FAO, accessed July 6, 2022); https://www.fao.org/faostat/en/#data/RFNHutton, M. O. et al. Toward a nitrogen footprint calculator for Tanzania. Env. Res. Lett. 12, 034016 (2017).Article 

    Google Scholar 
    Tanzania Fertilizer Assessment (International Fertilizer Development Center, 2012); http://tanzania.countrystat.org/fileadmin/user_upload/countrystat_fenix/congo/docs/Tanzania%20Fertilizer%20Assessment%202012.pdfA Common Carbon Footprint Approach for the Dairy Sector: The IDF Guide to Standard Life Cycle Methodology (International Dairy Federation, 2015); https://www.fil-idf.org/wp-content/uploads/2016/09/Bulletin479-2015_A-common-carbon-footprint-approach-for-the-dairy-sector.CAT.pdfBruzzone, L., Bovolo, F. & Arino, O. European Space Agency land cover climate change initiative. ESA LC CCI data: high resolution land cover data via Centre for Environmental Data Analysis; https://climate.esa.int/en/projects/high-resolution-land-cover/ (2021)Characteristics of Markets for Animal Feeds Raw Materials in the East African Community: Focus on Maize Bran and Sunflower Seed Cake (Kilimo Trust, 2017).Ngunga, D. & Mwendia, S. Forage Seed System in Tanzania: A Review Report (Alliance of Biodiversity and CIAT, 2020).Nkombe, B.M. Investigation of the Potential for Forage Species to Enhance the Sustainability of Degraded Rangeland and Cropland Soils. MSc thesis, Ohio State Univ. (2016).Producer Prices (FAO, accessed 2021); http://www.fao.org/faostat/en/#data/PP More