Instituto de Biología Subtropical

The reconstruction of bacterial and archaeal genomes from shotgun metagenomes has enabled insights into the ecology and evolution of environmental and host-associated microbiomes. Here we applied this approach to >10,000 metagenomes collected from diverse habitats covering all of Earth's continents and oceans, including metagenomes from human and animal hosts, engineered environments, and natural and agricultural soils, to capture extant microbial, metabolic and functional potential. This comprehensive catalog includes 52,515 metagenome-assembled genomes representing 12,556 novel candidate species-level operational taxonomic units spanning 135 phyla. The catalog expands the known phylogenetic diversity of bacteria and archaea by 44% and is broadly available for streamlined comparative analyses, interactive exploration, metabolic modeling and bulk download. We demonstrate the utility of this collection for understanding secondary-metabolite biosynthetic potential and for resolving thousands of new host linkages to uncultivated viruses. This resource underscores the value of genome-centric approaches for revealing genomic properties of uncultivated microorganisms that affect ecosystem processes.

<p>Biomass is defined as organic matter from living organisms represented in all kingdoms. It is recognized to be an excellent source of proteins, polysaccharides and lipids and, as such, embodies a tailored feedstock for new products and processes to apply in green industries. The industrial processes focused on the valorization of terrestrial biomass are well established, but marine sources still represent an untapped resource. Oceans and seas occupy over 70% of the Earth’s surface and are used intensively in worldwide economies through the fishery industry, as logistical routes, for mining ores and exploitation of fossil fuels, among others. All these activities produce waste. The other source of unused biomass derives from the beach wrack or washed-ashore organic material, especially in highly eutrophicated marine ecosystems. The development of high-added-value products from these side streams has been given priority in recent years due to the detection of a broad range of biopolymers, multiple nutrients and functional compounds that could find applications for human consumption or use in livestock/pet food, pharmaceutical and other industries. This review comprises a broad thematic approach in marine waste valorization, addressing the main achievements in marine biotechnology for advancing the circular economy, ranging from bioremediation applications for pollution treatment to energy and valorization for biomedical applications. It also includes a broad overview of the valorization of side streams in three selected case study areas: Norway, Scotland, and the Baltic Sea.</p>

Engel, Michael S, Ceríaco, Luis M P, Daniel, Gimo M, Dellapé, Pablo M, Löbl, Ivan, Marinov, Milen, Reis, Roberto E, Young, Mark T, Dubois, Alain, Agarwal, Ishan, Lehmann A., Pablo, Alvarado, Mabel, Alvarez, Nadir, Andreone, Franco, Araujo-Vieira, Katyuscia, Ascher, John S, Baêta, Délio, Baldo, Diego, Bandeira, Suzana A, Barden, Phillip, Barrasso, Diego A, Bendifallah, Leila, Bockmann, Flávio A, Böhme, Wolfgang, Borkent, Art, Brandão, Carlos R F, Busack, Stephen D, Bybee, Seth M, Channing, Alan, Chatzimanolis, Stylianos, Christenhusz, Maarten J M, Crisci, Jorge V, D'elía, Guillermo, Da Costa, Luis M, Davis, Steven R, De Lucena, Carlos Alberto S, Deuve, Thierry, Fernandes Elizalde, Sara, Faivovich, Julián, Farooq, Harith, Ferguson, Adam W, Gippoliti, Spartaco, Gonçalves, Francisco M P, Gonzalez, Victor H, Greenbaum, Eli, Hinojosa-Díaz, Ismael A, Ineich, Ivan, Jiang, Jianping, Kahono, Sih, Kury, Adriano B, Lucinda, Paulo H F, Lynch, John D, Malécot, Valéry, Marques, Mariana P, Marris, John W M, Mckellar, Ryan C, Mendes, Luis F, Nihei, Silvio S, Nishikawa, Kanto, Ohler, Annemarie, Orrico, Victor G D, Ota, Hidetoshi, Paiva, Jorge, Parrinha, Diogo, Pauwels, Olivier S G, Pereyra, Martín O, Pestana, Lueji B, Pinheiro, Paulo D P, Prendini, Lorenzo, Prokop, Jakub, Rasmussen, Claus, Rödel, Mark-Oliver, Rodrigues, Miguel Trefaut, Rodríguez, Sara M, Salatnaya, Hearty, Sampaio, Íris, Sánchez-García, Alba, Shebl, Mohamed A, Santos, Bruna S, Solórzano-Kraemer, Mónica M, Sousa, Ana C A, Stoev, Pavel, Teta, Pablo, Trape, Jean-François, Dos Santos, Carmen Van-Dúnem, Vasudevan, Karthikeyan, Vink, Cor J, Vogel, Gernot, Wagner, Philipp, Wappler, Torsten, Ware, Jessica L, Wedmann, Sonja, Zacharie, Chifundera Kusamba (2021): EDITORIAL The taxonomic impediment: a shortage of taxonomists, not the lack of technical approaches. Zoological Journal of the Linnean Society 193 (2): 381-387, DOI: 10.1093/zoolinnean/zlab072, URL: https://academic.oup.com/zoolinnean/article/193/2/381/6374389

The paradigm-changing opportunities of biologging sensors for ecological research, especially movement ecology, are vast, but the crucial questions of how best to match the most appropriate sensors and sensor combinations to specific biological questions and how to analyse complex biologging data, are mostly ignored. Here, we fill this gap by reviewing how to optimize the use of biologging techniques to answer questions in movement ecology and synthesize this into an Integrated Biologging Framework (IBF). We highlight that multisensor approaches are a new frontier in biologging, while identifying current limitations and avenues for future development in sensor technology. We focus on the importance of efficient data exploration, and more advanced multidimensional visualization methods, combined with appropriate archiving and sharing approaches, to tackle the big data issues presented by biologging. We also discuss the challenges and opportunities in matching the peculiarities of specific sensor data to the statistical models used, highlighting at the same time the large advances which will be required in the latter to properly analyse biologging data. Taking advantage of the biologging revolution will require a large improvement in the theoretical and mathematical foundations of movement ecology, to include the rich set of high-frequency multivariate data, which greatly expand the fundamentally limited and coarse data that could be collected using location-only technology such as GPS. Equally important will be the establishment of multidisciplinary collaborations to catalyse the opportunities offered by current and future biologging technology. If this is achieved, clear potential exists for developing a vastly improved mechanistic understanding of animal movements and their roles in ecological processes and for building realistic predictive models.

It is fundamentally important for many animal ecologists to quantify the costs of animal activities, although it is not straightforward to do so. The recording of triaxial acceleration by animal-attached devices has been proposed as a way forward for this, with the specific suggestion that dynamic body acceleration (DBA) be used as a proxy for movement-based power. Dynamic body acceleration has now been validated frequently, both in the laboratory and in the field, although the literature still shows that some aspects of DBA theory and practice are misunderstood. Here, we examine the theory behind DBA and employ modelling approaches to assess factors that affect the link between DBA and energy expenditure, from the deployment of the tag, through to the calibration of DBA with energy use in laboratory and field settings. Using data from a range of species and movement modes, we illustrate that vectorial and additive DBA metrics are proportional to each other. Either can be used as a proxy for energy and summed to estimate total energy expended over a given period, or divided by time to give a proxy for movement-related metabolic power. Nonetheless, we highlight how the ability of DBA to predict metabolic rate declines as the contribution of non-movement-related factors, such as heat production, increases. Overall, DBA seems to be a substantive proxy for movement-based power but consideration of other movement-related metrics, such as the static body acceleration and the rate of change of body pitch and roll, may enable researchers to refine movement-based metabolic costs, particularly in animals where movement is not characterized by marked changes in body acceleration.

Abstract Home range estimation is routine practice in ecological research. While advances in animal tracking technology have increased our capacity to collect data to support home range analysis, these same advances have also resulted in increasingly autocorrelated data. Consequently, the question of which home range estimator to use on modern, highly autocorrelated tracking data remains open. This question is particularly relevant given that most estimators assume independently sampled data. Here, we provide a comprehensive evaluation of the effects of autocorrelation on home range estimation. We base our study on an extensive data set of GPS locations from 369 individuals representing 27 species distributed across five continents. We first assemble a broad array of home range estimators, including Kernel Density Estimation ( KDE ) with four bandwidth optimizers (Gaussian reference function, autocorrelated‐Gaussian reference function [ AKDE ], Silverman's rule of thumb, and least squares cross‐validation), Minimum Convex Polygon, and Local Convex Hull methods. Notably, all of these estimators except AKDE assume independent and identically distributed ( IID ) data. We then employ half‐sample cross‐validation to objectively quantify estimator performance, and the recently introduced effective sample size for home range area estimation () to quantify the information content of each data set. We found that AKDE 95% area estimates were larger than conventional IID ‐based estimates by a mean factor of 2. The median number of cross‐validated locations included in the hold‐out sets by AKDE 95% (or 50%) estimates was 95.3% (or 50.1%), confirming the larger AKDE ranges were appropriately selective at the specified quantile. Conversely, conventional estimates exhibited negative bias that increased with decreasing . To contextualize our empirical results, we performed a detailed simulation study to tease apart how sampling frequency, sampling duration, and the focal animal's movement conspire to affect range estimates. Paralleling our empirical results, the simulation study demonstrated that AKDE was generally more accurate than conventional methods, particularly for small . While 72% of the 369 empirical data sets had >1,000 total observations, only 4% had an >1,000, where 30% had an <30. In this frequently encountered scenario of small , AKDE was the only estimator capable of producing an accurate home range estimate on autocorrelated data.

The species-area relationship (SAR) gives a quantitative description of the increasing number of species in a community with increasing area of habitat. In conservation, SARs have been used to predict the number of extinctions when the area of habitat is reduced. Such predictions are most needed for landscapes rather than for individual habitat fragments, but SAR-based predictions of extinctions for landscapes with highly fragmented habitat are likely to be biased because SAR assumes contiguous habitat. In reality, habitat loss is typically accompanied by habitat fragmentation. To quantify the effect of fragmentation in addition to the effect of habitat loss on the number of species, we extend the power-law SAR to the species-fragmented area relationship. This model unites the single-species metapopulation theory with the multispecies SAR for communities. We demonstrate with a realistic simulation model and with empirical data for forest-inhabiting subtropical birds that the species-fragmented area relationship gives a far superior prediction than SAR of the number of species in fragmented landscapes. The results demonstrate that for communities of species that are not well adapted to live in fragmented landscapes, the conventional SAR underestimates the number of extinctions for landscapes in which little habitat remains and it is highly fragmented.

Aim: The introduction of aquatic non-indigenous species (ANS) has become a major driver for global changes in species biogeography. We examined spatial patterns and temporal trends of ANS detections since 1965 to inform conservation policy and management. Location: Global. Methods: We assembled an extensive dataset of first records of detection of ANS (1965-2015) across 49 aquatic ecosystems, including the (a) year of first collection, (b) population status and (c) potential pathway(s) of introduction. Data were analysed at global and regional levels to assess patterns of detection rate, richness and transport pathways. Results: ) primary detections of ANS occurred-one new detection every 8.4 days for 50 years. The global rate of detections was relatively stable during 1965-1995, but increased rapidly after this time, peaking at roughly 66 primary detections per year during 2005-2010 and then declining marginally. Detection rates were variable within and across regions through time. Arthropods, molluscs and fishes were the most frequently reported ANS. Most ANS were likely introduced as stowaways in ships' ballast water or biofouling, although direct evidence is typically absent. Main conclusions: This synthesis highlights the magnitude of recent ANS detections, yet almost certainly represents an underestimate as many ANS go unreported due to limited search effort and diminishing taxonomic expertise. Temporal rates of detection are also confounded by reporting lags, likely contributing to the lower detection rate observed in recent years. There is a critical need to implement standardized, repeated methods across regions and taxa to improve the quality of global-scale comparisons and sustain core measures over longer time-scales. It will be fundamental to fill in knowledge gaps given that invasion data representing broad regions of the world's oceans are not yet readily available and to maintain knowledge pipelines for adaptive management.

Extinction risk is inversely associated with maximum per capita population growth rate (r(max)). However, this parameter is not known for most threatened species, underscoring the value in identifying correlates of r(max) that, in the absence of demographic data, would indirectly allow one to identify species and populations at elevated risk of extinction and their associated recovery potential. We undertook a comparative life-history analysis of 199 species from three taxonomic classes: Chondrichthyes (e.g., sharks; n = 82), Actinopterygii (teleost or bony fishes; n = 47), and Mammalia (n = 70, including 16 marine species). Median r(max) was highest for (and similar between) terrestrial mammals (0.71) and teleosts (0.43), significantly lower among chondrichthyans (0.26), and lower still in marine mammals (0.07). Age at maturity was the primary (and negative) correlate of r(max). In contrast, although body size was negatively correlated with r(max) in chondrichthyans and mammals, evidence of an association in teleosts was equivocal, and fecundity was not related to r(max) in fishes, despite recurring assertions to the contrary. Our analyses suggest that age at maturity can serve as a universal predictor of extinction risk in fishes and mammals when r(max) itself is unknown. Moreover, in contrast to what is generally expected, the recovery potential of teleost fishes does not differ from that of terrestrial mammals. Our findings are supportive of the application of extinction-risk criteria that are based on generation time and that are independent of taxonomic affinity.

Estimates of abundance are extremely valuable for species conservation, yet determining abundance for elusive, wide-ranging, carnivores is difficult. We estimated density of pumas using remote cameras across study sites in Bolivia, Argentina, and Belize. We used obvious and subtle markings to identify individual pumas in photographs and conducted double-blind identifications to examine the degree of agreement among investigators. Average agreement on identification between pairs of investigators was nearly 80.0% and 3-way agreement was 72.9%. Identification of pumas as different individuals was uncommon (7.8% pairwise, 0.69% 3-way disagreement) with the remainder described as unidentifiable. Densities of pumas varied consistently from site to site regardless of investigator. Bolivian pumas moved the shortest distances between camera stations and Argentinean pumas the longest, but distances among cameras and area covered by surveys varied among sites. We applied a correction factor to the Bolivian data to account for the small area surveyed and found that, averaged across investigator, Bolivia had significantly more pumas per 100 km2 (mean ± SD; 6.80 ± 1.5) than Belize (3.42 ± 1.3) or Argentina (0.67 ± 0.2). Numbers of pumas in Argentina match those of low-density North American sites, and those for Belize are consistent with the Pantanal and high-density North American sites. Densities of pumas can be reliably estimated with remote cameras for these sites, and our work presents the 1st density estimates for Central America and for forested environments in South America.

Habitat fragmentation may disrupt original patterns of gene flow and lead to drift-induced differentiation among local population units. Top predators such as the jaguar may be particularly susceptible to this effect, given their low population densities, leading to small effective sizes in local fragments. On the other hand, the jaguar's high dispersal capabilities and relatively long generation time might counteract this process, slowing the effect of drift on local populations over the time frame of decades or centuries. In this study, we have addressed this issue by investigating the genetic structure of jaguars in a recently fragmented Atlantic Forest region, aiming to test whether loss of diversity and differentiation among local populations are detectable, and whether they can be attributed to the recent effect of drift. We used 13 microsatellite loci to characterize the genetic diversity present in four remnant populations, and observed marked differentiation among them, with evidence of recent allelic loss in local areas. Although some migrant and admixed individuals were identified, our results indicate that recent large-scale habitat removal and fragmentation among these areas has been sufficiently strong to promote differentiation induced by drift and loss of alleles at each site. Low estimated effective sizes supported the inference that genetic drift could have caused this effect within a short time frame. These results indicate that jaguars' ability to effectively disperse across the human-dominated landscapes that separate the fragments is currently very limited, and that each fragment contains a small, isolated population that is already suffering from the effects of genetic drift.

Abstract Aim The anthropogenic climate change and land use change are considered two of the main factors that are altering biodiversity at the global scale. An evaluation that combined both factors can be relevant to detect which species could be the most vulnerable and reveal the regions of highest stability or susceptibility to biodiversity. We aimed to: (a) assess the effect of climate change and land use on the distribution of Cerrado plant species for different countries where they occur, (b) evaluate the effectiveness of the current network of protected areas (PAs) to safeguards species under different greenhouse–gas (GHG) emissions and land use scenarios, and (c) estimate the vulnerability of species based on protection effectiveness and habitat loss. Location Bolivia, Brazil and Paraguay. Methods We modelled the distribution of 1,553 plant species of Cerrado and evaluated species range loss caused by present and future land use and two GHG for 2050 and 2080. We assessed species vulnerability combining the representativeness of species within conservation units with the loss of species’ ranges outside PAs. Results We found that climate change and land use will cause great damage to Cerrado flora by 2050 and 2080, even under optimistic conditions. The greatest impacts of land use will occur in the regions where the greatest richness will be harboured. The conservation of the species will be seriously affected since the PA network is not as effective in safeguarding them under current or future conditions. Main conclusions The low level of protection together with the losses caused by the advance of agricultural lands will lead most species being highly vulnerable. Due to the distinct impacts of climate and land use over the three countries, conservation strategies should be implemented at transboundary and national levels.

Although invasive alien species have long been recognized as a major threat to nature and people, until now there has been no comprehensive global review of the status, trends, drivers, impacts, management and governance challenges of biological invasions. The Intergovernmental Science-Policy Platform on Biodiversity and Ecosystem Services (IPBES) Thematic Assessment Report on Invasive Alien Species and Their Control (hereafter ‘IPBES invasive alien species assessment’) drew on more than 13,000 scientific publications and reports in 15 languages as well as Indigenous and local knowledge on all taxa, ecosystems and regions across the globe. Therefore, it provides unequivocal evidence of the major and growing threat of invasive alien species alongside ambitious but realistic approaches to manage biological invasions. The extent of the threat and impacts has been recognized by the 143 member states of IPBES who approved the summary for policymakers of this assessment. Here, the authors of the IPBES assessment outline the main findings of the IPBES invasive alien species assessment and highlight the urgency to act now. This Perspective highlights the global consensus on the urgency and growing threat of invasive alien species, and management needs, as found by the 2023 report on invasive alien species conducted by the Intergovernmental Science-Policy Platform on Biodiversity and Ecosystem Services (IPBES).

Sharks are one of the most threatened groups of marine animals, as high exploitation rates coupled with low resilience to fishing pressure have resulted in population declines worldwide. Designing conservation strategies for this group depends on basic knowledge of the geographic distribution and diversity of known species. So far, this information has been fragmented and incomplete. Here, we have synthesized the first global shark diversity pattern from a new database of published sources, including all 507 species described at present, and have identified hotspots of shark species richness, functional diversity and endemicity from these data. We have evaluated the congruence of these diversity measures and demonstrate their potential use in setting priority areas for shark conservation. Our results show that shark diversity across all species peaks on the continental shelves and at mid-latitudes (30-40 degrees N and S). Global hotspots of species richness, functional diversity and endemicity were found off Japan, Taiwan, the East and West coasts of Australia, Southeast Africa, Southeast Brazil and Southeast USA. Moreover, some areas with low to moderate species richness such as Southern Australia, Angola, North Chile and Western Continental Europe stood out as places of high functional diversity. Finally, species affected by shark finning showed different patterns of diversity, with peaks closer to the Equator and a more oceanic distribution overall. Our results show that the global pattern of shark diversity is uniquely different from land, and other well-studied marine taxa, and may provide guidance for spatial approaches to shark conservation. However, similar to terrestrial ecosystems, protected areas based on hotspots of diversity and endemism alone would provide insufficient means for safeguarding the diverse functional roles that sharks play in marine ecosystems.

The jaguar is the top predator of the Atlantic Forest (AF), which is a highly threatened biodiversity hotspot that occurs in Brazil, Paraguay and Argentina. By combining data sets from 14 research groups across the region, we determine the population status of the jaguar and propose a spatial prioritization for conservation actions. About 85% of the jaguar's habitat in the AF has been lost and only 7% remains in good condition. Jaguars persist in around 2.8% of the region, and live in very low densities in most of the areas. The population of jaguars in the AF is probably lower than 300 individuals scattered in small sub-populations. We identified seven Jaguar Conservation Units (JCUs) and seven potential JCUs, and only three of these areas may have ≥50 individuals. A connectivity analysis shows that most of the JCUs are isolated. Habitat loss and fragmentation were the major causes for jaguar decline, but human induced mortality is the main threat for the remaining population. We classified areas according to their contribution to jaguar conservation and we recommend management actions for each of them. The methodology in this study could be used for conservation planning of other carnivore species.

Abstract Aim Trait data are widely used in ecological and evolutionary phylogenetic comparative studies, but often values are not available for all species of interest. Traditionally, researchers have excluded species without data from analyses, but estimation of missing values using imputation has been proposed as a better approach. However, imputation methods have largely been designed for randomly missing data, whereas trait data are often not missing at random (e.g., more data for bigger species). Here, we evaluate the performance of approaches for handling missing values when considering biased datasets. Location Any. Time period Any. Major taxa studied Any. Methods We simulated continuous traits and separate response variables to test the performance of nine imputation methods and complete‐case analysis (excluding missing values from the dataset) under biased missing data scenarios. We characterized performance by estimating the error in imputed trait values (deviation from the true value) and inferred trait–response relationships (deviation from the true relationship between a trait and response). Results Generally, Rphylopars imputation produced the most accurate estimate of missing values and best preserved the response–trait slope. However, estimates of missing data were still inaccurate, even with only 5% of values missing. Under severe biases, errors were high with every approach. Imputation was not always the best option, with complete‐case analysis frequently outperforming Mice imputation and, to a lesser degree, BHPMF imputation. Mice , a popular approach, performed poorly when the response variable was excluded from the imputation model. Main conclusions Imputation can handle missing data effectively in some conditions but is not always the best solution. None of the methods we tested could deal effectively with severe biases, which can be common in trait datasets. We recommend rigorous data checking for biases before and after imputation and propose variables that can assist researchers working with incomplete datasets to detect data biases and minimize errors.

Accurately estimating home range and understanding movement behavior can provide important information on ecological processes. Advances in data collection and analysis have improved our ability to estimate home range and movement parameters, both of which have the potential to impact species conservation. Fitting continuous-time movement model to data and incorporating the autocorrelated kernel density estimator (AKDE), we investigated range residency of forty-four jaguars fit with GPS collars across five biomes in Brazil and Argentina. We assessed home range and movement parameters of range resident animals and compared AKDE estimates with kernel density estimates (KDE). We accounted for differential space use and movement among individuals, sex, region, and habitat quality. Thirty-three (80%) of collared jaguars were range resident. Home range estimates using AKDE were 1.02 to 4.80 times larger than KDE estimates that did not consider autocorrelation. Males exhibited larger home ranges, more directional movement paths, and a trend towards larger distances traveled per day. Jaguars with the largest home ranges occupied the Atlantic Forest, a biome with high levels of deforestation and high human population density. Our results fill a gap in the knowledge of the species' ecology with an aim towards better conservation of this endangered/critically endangered carnivore-the top predator in the Neotropics.

Standardised terminology in science is important for clarity of interpretation and communication. In invasion science - a dynamic and rapidly evolving discipline - the proliferation of technical terminology has lacked a standardised framework for its development. The result is a convoluted and inconsistent usage of terminology, with various discrepancies in descriptions of damage and interventions. A standardised framework is therefore needed for a clear, universally applicable, and consistent terminology to promote more effective communication across researchers, stakeholders, and policymakers. Inconsistencies in terminology stem from the exponential increase in scientific publications on the patterns and processes of biological invasions authored by experts from various disciplines and countries since the 1990s, as well as publications by legislators and policymakers focusing on practical applications, regulations, and management of resources. Aligning and standardising terminology across stakeholders remains a challenge in invasion science. Here, we review and evaluate the multiple terms used in invasion science (e.g. 'non-native', 'alien', 'invasive' or 'invader', 'exotic', 'non-indigenous', 'naturalised', 'pest') to propose a more simplified and standardised terminology. The streamlined framework we propose and translate into 28 other languages is based on the terms (i) 'non-native', denoting species transported beyond their natural biogeographic range, (ii) 'established non-native', i.e. those non-native species that have established self-sustaining populations in their new location(s) in the wild, and (iii) 'invasive non-native' - populations of established non-native species that have recently spread or are spreading rapidly in their invaded range actively or passively with or without human mediation. We also highlight the importance of conceptualising 'spread' for classifying invasiveness and 'impact' for management. Finally, we propose a protocol for classifying populations based on (i) dispersal mechanism, (ii) species origin, (iii) population status, and (iv) impact. Collectively and without introducing new terminology, the framework that we present aims to facilitate effective communication and collaboration in invasion science and management of non-native species.

Abstract Species distribution models (SDM) are widely used in diverse research areas because of their simple data requirements and application versatility. However, SDM outcomes are sensitive to data input and methodological choices. Such sensitivity and diverse applications mean that flexibility is necessary to create SDMs with tailored protocols for a given set of data and model use. We introduce the r package flexsdm for supporting flexible species distribution modelling workflows. flexsdm functions and their arguments serve as building blocks to construct a specific modelling protocol for user's needs. The main flexsdm features are modelling flexibility, integration with other modelling tools, simplicity of the objects returned and function speed. As an illustration, we used flexsdm to define a complete workflow for California red fir Abies magnifica . This package provides modelling flexibility by incorporating comprehensive tools structured in three steps: (a) The Pre‐modelling functions that prepare input, for example, sampling bias correction, sampling pseudo‐absences and background points, data partitioning, and reducing collinearity in predictors. (b) The Modelling functions allow fitting and evaluating different modelling approaches, including individual algorithms, tuned models, ensembles of small models and ensemble models. (c) The Post‐modelling functions include tools related to models' predictions, interpolation and overprediction correction. Because flexsdm comprises a large part of the SDM process, from outlier detection to overprediction correction, flexsdm users can delineate partial or complete workflows based on the combination functions to meet specific modelling needs.

Wood biophysical properties and the dynamics of water storage discharge and refilling were studied in the trunk of canopy tree species with diverse life history and functional traits in subtropical forests of northeast Argentina. Multiple techniques assessing capacitance and storage capacity were used simultaneously to improve our understanding of the functional significance of internal water sources in trunks of large trees. Sapwood capacitances of 10 tree species were characterized using pressure-volume relationships of sapwood samples obtained from the trunk. Frequency domain reflectometry was used to continuously monitor the volumetric water content in the main stems. Simultaneous sap flow measurements on branches and at the base of the tree trunk, as well as diurnal variations in trunk contraction and expansion, were used as additional measures of stem water storage use and refilling dynamics. All evidence indicates that tree trunk internal water storage contributes from 6 to 28% of the daily water budget of large trees depending on the species. The contribution of stored water in stems of trees to total daily transpiration was greater for deciduous species, which exhibited higher capacitance and lower sapwood density. A linear relationship across species was observed between wood density and growth rates with the higher wood density species (mostly evergreen) associated with lower growth rates and the lower wood density species (mostly deciduous) associated with higher growth rates. The large sapwood capacitance in deciduous species may help to avoid catastrophic embolism in xylem conduits. This may be a low-cost adaptation to avoid water deficits during peak water use at midday and under temporary drought periods and will contribute to higher growth rates in deciduous tree species compared with evergreen ones. Large capacitance appears to have a central role in the rapid growth patterns of deciduous species facilitating rapid canopy access as these species are less shade tolerant than evergreen species.