Age,growth and reproductive life-history characteristics infer a high population productivity for the sustainably fished protogynous hermaphroditic yellowspotted rockcod (Epinephelus areolatus) in north-western Australia

The yellowspotted rockcod, Epinephelus areolatus, is a small-sized grouper that is widely distributed throughout the Indo-Pacific, where it forms a valuable component of the harvest derived from multispecies fisheries along continental and insular shelves. Samples of E. areolatus were collected from 2012 to 2018 from commercial catches and research surveys in the Kimberley, Pilbara and Gascoyne regions of north-western Australia to improve the understanding of the life history, inherent vulnerability and stock status of this species. Histological analysis of gonads (n = 1889) determined that E. areolatus was a monandric protogynous hermaphrodite. Non-functional spermatogenic crypts were dispersed within the ovaries of 23% of mature functioning females; nonetheless, these crypts were not observed during the immature female phase. The length and age at which 50% of females matured were 266 mm total length (L_T) and 2.7 years, respectively. The spawning period was protracted over 10–12 months of the year with biannual peaks at the start of spring and autumn (i.e., September and March) when the photoperiod was at its mid-range (i.e., 12.1 h). Estimates of the lengths and ages at which 50% of E. areolatus change sex from female to male were very similar (i.e., <5% difference) between the Kimberley and Pilbara regions, i.e., of 364 and 349 mm L_T and of 7.9 and 7.3 years, respectively. A maximum age of 19 years was observed in all three regions, but there was significant regional variation in growth. These variations in growth were not correlated with latitude; instead a parabolic relationship was evident, where the smallest mean length-at-age and fastest growth rates (k) occurred in the mid-latitudes of the Pilbara region. In the Kimberley and Pilbara regions, individuals were not fully selected by commercial fish traps until 5–6 years of age, hence, several years after reaching maturity. These life-history characteristics infer a high population productivity, which underpins the sustainable harvest of this species, despite comprising the largest catches of all epinephelids in the multispecies tropical fisheries across north-western Australia.     

Age and growth of the bull shark (Carcharhinus leucas) around Reunion Island,South West Indian Ocean

Sharks exhibit varied demographic strategies depending on both the species and the population location, which make them more or less vulnerable to fishing. Accurate evaluation of local age and growth parameters is therefore fundamental for the sustainable management of their stocks. Although demographic parameters have been assessed for bull shark (Carcharhinus leucas) populations in several locations of the world, this information is missing so far around Reunion Island, in the south-west Indian Ocean. To fill this gap of knowledge, age and growth data were gathered from the vertebrae of 140 individuals of C. leucas (77 females and 63 males, mostly adults) fished around the island between 2012 and 2019. After verification of the annual deposition of growth band pairs on these structures using relative marginal increment analysis on 40 individuals, band pairs were counted along the vertebral centrum for each individual. Thanks to this approach, growth was shown to significantly differ between male and female C. leucas around Reunion Island, with respective von Bertalanffy growth model equations of L_t =  and L_t = . Indeed, the females of the species fished in this area were significantly (P < 0.001) larger than local males, with an estimated difference in size of ~16.1 cm at 20 years old. They also apparently reach older ages, with an estimated maximum age of 33.50 years, against 29.75 years only for the males. The estimated size at birth around the island is larger than elsewhere in the world, varying from 92.30 to 100.00 cm depending on the method used. These results confirm that the population of C. leucas around Reunion Island exhibits a K-selected strategy, which makes it highly vulnerable to fishing pressure.     

Respirometer in a box: development and use of a portable field respirometer for estimating oxygen consumption of large-bodied fishes

This study developed a portable, low-cost field respirometer for measuring oxygen consumption rates of large-bodied fishes. The respirometer performed well in laboratory tests and was used to measure the oxygen consumption rates ( O₂) of bull sharks Carcharhinus leucas (mean: 249.21 ± 58.10 mg O₂ kg⁻¹ h⁻¹ at 27.05°C). Interspecific comparisons and assessments of oxygen degradation curves indicated that the respirometer provided reliable measurements of O₂. This system presents a field-based alternative to laboratory respirometers, opening opportunities for studies on species in remote localities, increasing the ability to validate physiological field studies.     

Reproductive biology of female bigeye tuna Thunnus obesus in the western Pacific Ocean

The reproductive biology of female bigeye tuna Thunnus obesus was assessed by examining 888 fish (ranging from 84·9 to 174·4 cm fork length, L_F) caught by Taiwanese offshore longliners in the western Pacific Ocean from November 1997 to November 1998 and November to December 1999 and 258 gonad samples from these fish. The overall sex ratio of the catch during the sampling differed significantly from 0·5, but males were predominant in sizes >140 cm L_F. Reproductive activity (assessed by histology), a gonado‐somatic index, and the size‐frequency distributions of whole oocytes indicated that spawning occurred throughout the year and the major spawning season appeared to be from February to September. The estimated sizes at 50% maturity (L_F50) of females was 102·85 cm (95% c.i .: 90·79–110·21 cm) and the smallest mature female was 99·7 cm L_F. They are multiple spawners and oocytes develop asynchronously. The proportion of mature (0·63) and reproductively active (0·70) females with ovaries containing postovulatory follicles indicated that they spawn almost daily. Batch fecundity for 15 females with the most advanced oocytes (>730 µm) ranged from 0·84 to 8·56 million eggs (mean ± s.d . = 3·06 ± 2·09). The relationships between batch fecundity (F_B, in millions of eggs) and L_F (cm) and round mass (M_R, kg) were (r² = 0·84) and (r² = 0·80), respectively. The parameters estimated in this study are key information for stock assessments of T. obesus in the western Pacific Ocean and will contribute to the conservation and sustainable yield of this species.     

Developmental changes in oxygen consumption and hypoxia tolerance in the heat and hypoxia-adapted tabasco line of the Nile tilapia Oreochromis niloticus,with a survey of the metabolic literature for the genus Oreochromis

The genus Oreochromis is among the most popular of the tilapiine cichlid tribe for aquaculture. However, their temperature and hypoxia tolerance, if tested at all, is usually tested at temperatures of 20–25°C, rather than at the considerably higher temperatures of 30–35°C typical of tropical aquaculture. We hypothesized that both larvae and adults of the heat and hypoxia-adapted Tabasco-line of the Nile tilapia Oreochromis niloticus would be relatively hypoxia-tolerant. Oxygen consumption rate (), Q₁₀ and aquatic surface respiration (ASR) was measured using closed respirometry at 2 (c. 0.2 g), 30 (c. 2–5 g), 105 c. (10–15 g) and 240 (c. 250 g) days of development, at 25°C, 30°C and 35°C. at 30°C was inversely related to body mass: c. 90 μM O₂ g⁻¹/h in larvae down to c. 1 μM O₂ g⁻¹/h in young adults. Q₁₀ for was typical for fish over the range 25–35°C of 1.5–2.0. ASR was exhibited by 50% of the fish at pO₂ of 15–50 mmHg in a temperature-dependent fashion. However, the largest adults showed notable ASR only when pO₂ fell to below 10 mmHg. Remarkably, p_crit for was 12–17 mmHg at 25–30°C and still only 20–25 mmHg across development at 35°C. These values are among the lowest measured for teleost fishes. Noteworthy is that all fish maintain equilibrium, ventilated their gills and showed routine locomotor action for 10–20 min after ceased at near anoxia and when then returned to oxygenated waters, all fish survived, further indicating a remarkable hypoxic tolerance. Remarkably, data assembled for from >30 studies showed a > x2000 difference, which we attribute to calculation or conversion errors. Nonetheless, p_crit was very low for all Oreochromis sp. and lowest in the heat and hypoxia-adapted Tabasco line.     

Estimating the decision curve and its precision from three study designs

The decision curve plots the net benefit of a risk model for making decisions over a range of risk thresholds, corresponding to different ratios of misclassification costs. We discuss three methods to estimate the decision curve, together with corresponding methods of inference and methods to compare two risk models at a given risk threshold. One method uses risks (R) and a binary event indicator (Y) on the entire validation cohort. This method makes no assumptions on how well-calibrated the risk model is nor on the incidence of disease in the population and is comparatively robust to model miscalibration. If one assumes that the model is well-calibrated, one can compute a much more precise estimate of based on risks R alone. However, if the risk model is miscalibrated, serious bias can result. Case–control data can also be used to estimate if the incidence (or prevalence) of the event () is known. This strategy has comparable efficiency to using the full data, and its efficiency is only modestly less than that for the full data if the incidence is estimated from the mean of Y. We estimate variances using influence functions and propose a bootstrap procedure to obtain simultaneous confidence bands around the decision curve for a range of thresholds. The influence function approach to estimate variances can also be applied to cohorts derived from complex survey samples instead of simple random samples.     

Valid oxygen uptake measurements: using high r2 values with good intentions can bias upward the determination of standard metabolic rate

This analysis shows good intentions in the selection of valid and precise oxygen uptake (O₂) measurements by retaining only slopes of declining dissolved oxygen level in a respirometer that have very high values of the coefficient of determination, r², are not always successful at excluding nonlinear slopes. Much worse, by potentially removing linear slopes that have low r² only because of a low signal-to-noise ratio, this procedure can overestimate the calculation of standard metabolic rate (SMR) of the fish. To remedy this possibility, a few simple diagnostic tools are demonstrated to assess the appropriateness of a given minimum acceptable r², such as calculating the proportion of rejected O₂ determinations, producing a histogram of the r² values and a plot of r² as a function of O₂. The authors offer solutions for cases when many linear slopes have low r². The least satisfactory but easiest to implement is lowering the minimum acceptable r². More satisfactory solutions involve processing (smoothing) the raw signal of dissolved oxygen as a function of time to improve the signal-to-noise ratio and the r²s.     

An estimator of the Opportunity for Selection that is independent of mean fitness

Variation among individuals in number of offspring (fitness, k) sets an upper limit to the evolutionary response to selection. This constraint is quantified by Crow's Opportunity for Selection (I), which is the variance in relative fitness (I = σ²_k/(u_k)²). Crow's I has been widely used but remains controversial because it depends on mean offspring number in a sample (). Here, I used a generalized Wright-Fisher model that allows for unequal probabilities of producing offspring to evaluate behavior of Crow's I and related indices under a wide range of sampling scenarios. Analytical and numerical results are congruent and show that rescaling the sample variance (s²_k) to its expected value at a fixed removes dependence of I on mean offspring number, but the result still depends on choice of . A new index is introduced, Δ_I = Π– E(Î_drift) = Π– 1/, which makes Î independent of sample without the need for variance rescaling. Δ_I has a straightforward interpretation as the component of variance in relative fitness that exceeds that expected under a null model of random reproductive success. Δ_I can be used to directly compare estimates of the Opportunity for Selection for samples from different studies, different sexes, and different life stages.     

Excess post‐hypoxic oxygen consumption in Atlantic cod Gadus morhua

Atlantic cod Gadus morhua experienced oxygen deficit () when exposed to oxygen levels below their critical level (c. 73% of p_crit) and subsequent excess post‐hypoxic oxygen consumption (C_EPHO) upon return to normoxic conditions, indicative of an oxygen debt. The mean ± s.e . C_EPHO: was 6·9 ± 1·5, suggesting that resorting to anaerobic energy production in severe hypoxia is energetically expensive.     

Satiation meal and the effects of meal and body sizes on gastric evacuation rate in brook trout Salvelinus fontinalis fed commercial pellets

Gastric evacuation (GE) experiments were performed on brook trout Salvelinus fontinalis fed commercial food pellets. The experiments included small fish (36 g; 15 cm total length, L_T) fed meals of 0·2, 0·4 and 0·8 g and large fish (152 g; 23 cm) fed meals of 0·8, 2·0 and 4·0 g at temperatures ranging from 15·1 to 18·2° C. The stomach contents were thereafter sampled and weighed at 3 h intervals until the first empty stomach was observed. The course of GE was examined by use of a general power function of the data that revealed that the square‐root function described the GE rate (GER) by the current stomach content mass independently of original meal size. Using the square‐root function, the relationship between GER and fish size was described by a power function of fish length, whereas the effect of temperature was described by a simple exponential function. GER of the commercial pellets fed to S. fontinalis could thus be described by (g h^?1), where S_t is stomach mass (g) at time t (h), L is total fish length (cm) and T is temperature (° C). The result of this study should provide a useful tool for planning of feeding regimes in production of S. fontinalis by optimizing growth and minimizing food waste.     

Another chemolithotrophic metabolism missing in nature: sulfur comproportionation

Chemotrophic microorganisms gain energy for cellular functions by catalyzing oxidation–reduction (redox) reactions that are out of equilibrium. Calculations of the Gibbs energy ( ΔG _r ) can identify whether a reaction is thermodynamically favourable and quantify the accompanying energy yield at the temperature, pressure and chemical composition in the system of interest. Based on carefully calculated values of ΔG _r , we predict a novel microbial metabolism – sulfur comproportionation (3H₂S + + 2H⁺ ⇌ 4S⁰ + 4H₂O). We show that at elevated concentrations of sulfide and sulfate in acidic environments over a broad temperature range, this putative metabolism can be exergonic ( ΔG _r <0), yielding ~30–50 kJ mol⁻¹. We suggest that this may be sufficient energy to support a chemolithotrophic metabolism currently missing from the literature. Other versions of this metabolism, comproportionation to thiosulfate (H₂S + ⇌ + H₂O) and to sulfite (H₂S + 3 ⇌ 4 + 2H⁺), are only moderately exergonic or endergonic even at ideal geochemical conditions. Natural and impacted environments, including sulfidic karst systems, shallow-sea hydrothermal vents, sites of acid mine drainage, and acid–sulfate crater lakes, may be ideal hunting grounds for finding microbial sulfur comproportionators.     

NO and N2O transformations of diverse fungi in hypoxia: evidence for anaerobic respiration only in Fusarium strains

Fungal denitrification is claimed to produce non-negligible amounts of N₂O in soils, but few tested species have shown significant activity. We hypothesized that denitrifying fungi would be found among those with assimilatory nitrate reductase, and tested 20 such batch cultures for their respiratory metabolism, including two positive controls, Fusarium oxysporum and Fusarium lichenicola, throughout the transition from oxic to anoxic conditions in media supplemented with . Enzymatic reduction of (NIR) and NO (NOR) was assessed by correcting measured NO- and N₂O-kinetics for abiotic NO- and N₂O-production (sterile controls). Significant anaerobic respiration was only confirmed for the positive controls and for two of three Fusarium solani cultures. The NO kinetics in six cultures showed NIR but not NOR activity, observed through the accumulation of NO. Others had NOR but not NIR activity, thus reducing abiotically produced NO to N₂O. The presence of candidate genes (nirK and p450nor) was confirmed in the positive controls, but not in some of the NO or N₂O accumulating cultures. Based on our results, we conclude that only the Fusarium cultures were able to sustain anaerobic respiration and produced low amounts of N₂O as a response to an abiotic NO production from the medium.     

Optimal Strategies for Managing Wildlife Harvest Under System Change

Wildlife populations are experiencing shifting dynamics due to climate and landscape change. Management policies that fail to account for non-stationary dynamics may fail to achieve management objectives. We establish a framework for understanding optimal strategies for managing a theoretical harvested population under non-stationarity. Building from harvest theory, we develop scenarios representing changes in population growth rate () or carrying capacity () and derive time-dependent optimal harvest policies using stochastic dynamic programming. We then evaluate the cost of falsely assuming stationarity by comparing the outcomes of forward projections in which either the optimal policy or a stationary policy is applied. When declines over time, the stationary policy leads to an underharvest of the population, resulting in less harvest over the short term but leaving the population in a higher-value state. When declines over time, the stationary policy leads to overharvest, resulting in greater harvest returns in the short term but leaving the population in a lower and potentially more vulnerable state. This work demonstrates the basic properties of time-dependent harvest management and provides a framework for evaluating the many outstanding questions about optimal management strategies under climate change. Published 2021. This article is a U.S. Government work and is in the public domain in the USA.     

The effects of hurricanes on the stochastic population growth of the endemic epiphytic orchid Broughtonia cubensis living in Cuba

We carried out a posthurricane evaluation of Broughtonia cubensis (Lindl.) Cogn., an endemic Cuban epiphytic orchid, after Hurricane Ivan (2004). We studied the transient responses in the stochastic dynamics of the species at three different sites over 13 successive years (2006–2019), monitored plot inventories (464 individuals in 10 transects) and built stochastic population models. The deterministic stochastic growth rate values () did not significantly differ (F = 2.76; p > 0.076) among the three sites over the 2006–2019 period. The long-term stochastic growth rate was 0.973 [0.932, 1.034]. The matrix elements that had the largest effect on were the transition to and stasis within the largest size class. Transient responses explained an average of 86% of the variation in the observed population growth rates , compared to 4% of the variation in the vital rates . Because transient dynamics are dependent on the population size composition, we ran extinction risk analyses under two scenarios: a population composed mainly of juveniles and another composed mainly of adults. There was little risk of falling below the quasi-extinction threshold before 25 year for both juveniles and adults. However, the risk of quasi-extinction was almost certain for both size classes by 80 year. We also simulated the effect of increasing the hurricane occurrence probability over 80 year on the population. There was little risk of extinction before 20 year in the baseline model, but there was a significant risk of extinction within 5 year when 90% of the individuals were affected by a new hurricane event.     

Age and growth of mangrove red snapper Lutjanus argentimaculatus at its cool‐water‐range limits

This study investigates the age and growth of Lutjanus argentimaculatus at its southern (cooler) range limits in eastern Australia. Specimens were collected from New South Wales and southern Queensland between November 2011 and December 2013. Fork lengths (L_F) ranged from 190 to 1019 mm, and ages ranged from 2+ to 57+ years. Growth was described by the von Bertalanffy growth function with coefficients L_∞ = 874·92 mm, K = 0·087 year^?1 and t₀ = ?2·76 years. Estimates of the instantaneous natural mortality rate (M) ranged from 0·072 to 0·25. The L_F (mm) and mass (W; g) relationship was represented by the equation: . The maximum age of 57+ years is the oldest reported for any lutjanid and comparisons with tropical studies suggest that the age‐based demography of L. argentimaculatus follows a latitudinal gradient. High maximum ages and low natural mortality rates indicate considerable vulnerability to overexploitation at the species' cool‐water‐range limits. These results demonstrate the need to identify underlying processes driving latitudinal gradients in fish demography.     

Genetic Characteristics of Restored Elk Populations in Kentucky

Translocations are a common management practice to restore or augment populations. Understanding the genetic consequences of translocation efforts is important for the long-term health of restored populations. The restoration of elk (Cervus canadensis) to Kentucky, USA, included source stocks from 6 western states, which were released at 8 sites in southeastern Kentucky during 1997–2002. We assessed genetic diversity in restored herds and compared genetic similarity to source stocks based on 15 microsatellite DNA loci. Genetic variation in the restored populations was comparable to source stocks ( allelic richness = 3.52 and 3.50; expected heterozygosity = 0.665 and 0.661 for restored and source, respectively). Genetic differentiation among all source and restored populations ranged from 0.000 to 0.065 for pairwise F_ST and 0.034 to 0.161 for pairwise Nei's D_A. Pairwise genetic differentiation and Bayesian clustering revealed that stocks from Utah and North Dakota, USA, contributed most to restored populations. Other western stocks appeared less successful and were not detected with our data, though our sampling was not exhaustive. We also inferred natural movements of elk among release sites by the presence of multiple genetic stocks. The success of the elk restoration effort in Kentucky may be due, in part, to the large number of elk (n = 1,548), repeated releases, and use of diverse source stocks. Future restoration efforts for elk in the eastern United States should consider the use of multiple stock sources and a large number of individuals. In addition, preservation of genetic samples of founder stock will enable detailed monitoring in the future. © 2020 The Authors. The Journal of Wildlife Management published by Wiley Periodicals, Inc. on behalf of The Wildlife Society.     

Precision of the age–length increments of three cyprinids: effects of fish number and sub‐sampling strategy

The effects of number of fish that are aged and scale sub‐sampling strategies on the precision of estimates of mean age–length increments from populations of Rutilus rutilus, Leuciscus leuciscus and Leuciscus cephalus were tested. Analyses used data derived from river fish communities in eastern England, U.K.. Regarding the number of fishes analysed in each age group, for each species and mean fork‐length increment at age, significant relationships were detected between sample size (n) and the coefficient of variation of the mean (Z) and mean length increment and measured variance (s²). This enabled calculation of the number of scales for producing a mean length increment at age according to . Outputs indicated that the number of scales requiring ageing increased substantially as precision increased, but with little variation between species per age category. Ageing between seven and 12 scales per age group would thus provide estimates at 10% precision. As the ages of fishes are not known in advance of scale ageing, the effect of scale sub‐sampling regime on precision was also tested using randomized strategies of 10 fish per 5 mm, five per 5 mm, three per 5 mm, 10 per 10 mm, five per 10 mm and three per 10 mm. These were applied to the datasets and the consequences of their reduction in the number of scales for precision were determined using . When compared to no sub‐sampling, three per 10 mm always significantly reduced data precision, whereas 10 per 5 mm never significantly reduced precision. These outputs can thus be applied to the design of fish sampling protocols where age and growth estimates are required, with the randomized sub‐sampling likely to be the most useful strategy.     

The design heatmap: A simple visualization of -optimality design problems

Optimal experimental designs are often formal and specific, and not intuitively plausible to practical experimenters. However, even in theory, there often are many different possible design points providing identical or nearly identical information compared to the design points of a strictly optimal design. In practical applications, this can be used to find designs that are a compromise between mathematical optimality and practical requirements, including preferences of experimenters. For this purpose, we propose a derivative-based two-dimensional graphical representation of the design space that, given any optimal design is already known, will show which areas of the design space are relevant for good designs and how these areas relate to each other. While existing equivalence theorems already allow such an illustration in regard to the relevance of design points only, our approach also shows whether different design points contribute the same kind of information, and thus allows tweaking of designs for practical applications, especially in regard to the splitting and combining of design points. We demonstrate the approach on a toxicological trial where a -optimal design for a dose–response experiment modeled by a four-parameter log-logistic function was requested. As these designs require a prior estimate of the relevant parameters, which is difficult to obtain in a practical situation, we also discuss an adaption of our representations to the criterion of Bayesian -optimality. While we focus on -optimality, the approach is in principle applicable to different optimality criteria as well. However, much of the computational and graphical simplicity will be lost.     

Body shape and robustness response to water flow during development of brown trout Salmo trutta parr

Domesticated brown trout Salmo trutta parr were subjected to increased, variable flow under controlled experimental conditions. Using geometric morphometric analyses, (a mass–length index) and caudal fin area–body length ratio, this study assessed morphological responses in lateral body depth, growth and robustness and propulsive potential, respectively, of parr over the course of 32 weeks. Geometric morphometric analyses did not reveal an effect of exercise on either lateral body depth or caudal fin area. However, improved overall robustness and growth trajectories in exercised parr showed a positive adaptive response to the enriched habitat. Exercise and habitat heterogeneity thus have the potential to improve survivability of domesticated salmonids in the wild.