Projecting future expansion of invasive species: Comparing and improving methodologies for species distribution modeling

Modeling the distributions of species, especially of invasive species in non-native ranges, involves multiple challenges. Here, we developed some novel approaches to species distribution modeling aimed at reducing the influences of such challenges and improving the realism of projections. We estimated species-environment relationships with four modeling methods run with multiple scenarios of (1) sources of occurrences and geographically isolated background ranges for absences, (2) approaches to drawing background (absence) points, and (3) alternate sets of predictor variables. We further tested various quantitative metrics of model evaluation against biological insight. Model projections were very sensitive to the choice of training dataset. Model accuracy was much improved by using a global dataset for model training, rather than restricting data input to the species’ native range. AUC score was a poor metric for model evaluation and, if used alone, was not a useful criterion for assessing model performance. Projections away from the sampled space (i.e. into areas of potential future invasion) were very different depending on the modeling methods used, raising questions about the reliability of ensemble projections. Generalized linear models gave very unrealistic projections far away from the training region. Models that efficiently fit the dominant pattern, but exclude highly local patterns in the dataset and capture interactions as they appear in data (e.g. boosted regression trees), improved generalization of the models. Biological knowledge of the species and its distribution was important in refining choices about the best set of projections. A post-hoc test conducted on a new Partenium dataset from Nepal validated excellent predictive performance of our “best” model. We showed that vast stretches of currently uninvaded geographic areas on multiple continents harbor highly suitable habitats for Parthenium hysterophorus L. (Asteraceae; parthenium). However, discrepancies between model predictions and parthenium invasion in Australia indicate successful management for this globally significant weed. This article is protected by copyright. All rights reserved.

Projecting future expansion of invasive species: Comparing and improving methodologies for species distribution modeling

Modeling the distributions of species, especially of invasive species in non-native ranges, involves multiple challenges. Here, we developed some novel approaches to species distribution modeling aimed at reducing the influences of such challenges and improving the realism of projections. We estimated species-environment relationships with four modeling methods run with multiple scenarios of (1) sources of occurrences and geographically isolated background ranges for absences, (2) approaches to drawing background (absence) points, and (3) alternate sets of predictor variables. We further tested various quantitative metrics of model evaluation against biological insight. Model projections were very sensitive to the choice of training dataset. Model accuracy was much improved by using a global dataset for model training, rather than restricting data input to the species’ native range. AUC score was a poor metric for model evaluation and, if used alone, was not a useful criterion for assessing model performance. Projections away from the sampled space (i.e. into areas of potential future invasion) were very different depending on the modeling methods used, raising questions about the reliability of ensemble projections. Generalized linear models gave very unrealistic projections far away from the training region. Models that efficiently fit the dominant pattern, but exclude highly local patterns in the dataset and capture interactions as they appear in data (e.g. boosted regression trees), improved generalization of the models. Biological knowledge of the species and its distribution was important in refining choices about the best set of projections. A post-hoc test conducted on a new Partenium dataset from Nepal validated excellent predictive performance of our “best” model. We showed that vast stretches of currently uninvaded geographic areas on multiple continents harbor highly suitable habitats for Parthenium hysterophorus L. (Asteraceae; parthenium). However, discrepancies between model predictions and parthenium invasion in Australia indicate successful management for this globally significant weed. This article is protected by copyright. All rights reserved.

Influence of environmental and instrumental variables on the non-invasive prediction of Brix in pineapple using near infrared spectroscopy

The Brix content of pineapple fruit can be non-invasively predicted from the second derivative of near infrared reflectance spectra. Correlations obtained using a NIRSystems 6500 spectrophotometer through multiple linear regression and modified partial least squares analyses using a post-dispersive configuration were comparable with that from a pre-dispersive configuration in terms of accuracy (e.g. coefficient of determination, R2, 0.73; standard error of cross validation, SECV, 1.01°Brix). The effective depth of sample assessed was slightly greater using the post-dispersive technique (about 20 mm for pineapple fruit), as expected in relation to the higher incident light intensity, relative to the pre-dispersive configuration. The effect of such environmental variables as temperature, humidity and external light, and instrumental variables such as the number of scans averaged to form a spectrum, were considered with respect to the accuracy and precision of the measurement of absorbance at 876 nm, as a key term in the calibration for Brix, and predicted Brix. The application of post-dispersive near infrared technology to in-line assessment of intact fruit in a packing shed environment is discussed.

Non-invasive assessment of pineapple and mango fruit quality using near infra-red spectroscopy

The potential of near infra-red (NIR) spectroscopy for non-invasive measurement of fruit quality of pineapple (Ananas comosus var. Smooth Cayenne) and mango (Magnifera indica var. Kensington) fruit was assessed. A remote reflectance fibre optic probe, placed in contact with the fruit skin surface in a light-proof box, was used to deliver monochromatic light to the fruit, and to collect NIR reflectance spectra (760–2500 nm). The probe illuminated and collected reflected radiation from an area of about 16 cm2. The NIR spectral attributes were correlated with pineapple juice Brix and with mango flesh dry matter (DM) measured from fruit flesh directly underlying the scanned area. The highest correlations for both fruit were found using the second derivative of the spectra (d2 log 1/R) and an additive calibration equation. Multiple linear regression (MLR) on pineapple fruit spectra (n = 85) gave a calibration equation using d2 log 1/R at wavelengths of 866, 760, 1232 and 832 nm with a multiple coefficient of determination (R2) of 0.75, and a standard error of calibration (SEC) of 1.21 °Brix. Modified partial least squares (MPLS) regression analysis yielded a calibration equation with R2 = 0.91, SEC = 0.69, and a standard error of cross validation (SECV) of 1.09 oBrix. For mango, MLR gave a calibration equation using d2 log 1/R at 904, 872, 1660 and 1516 nm with R2 = 0.90, and SEC = 0.85% DM and a bias of 0.39. Using MPLS analysis, a calibration equation with R2 = 0.98, SEC = 0.54 and SECV = 1.19 was obtained. We conclude that NIR technology offers the potential to assess fruit sweetness in intact whole pineapple and DM in mango fruit, respectively, to within 1° Brix and 1% DM, and could be used for the grading of fruit in fruit packing sheds.

Phosphorus nutrition and tolerance of cotton to water stress: II. Water relations, free and bound water and leaf expansion rate

In previous experiments, increased leaf-Phosphorus (P) content with increasing P supply enhanced the individual leaf expansion and water content of fresh cotton leaves in a severely drying soil. In this paper, we report on the bulk water content of leaves and its components, free and bound water, along with other measures of plant water status, in expanding cotton leaves of various ages in a drying soil with different P concentrations. The bound water in living tissue is more likely to play a major role in tolerance to abiotic stresses by maintaining the structural integrity and/or cell wall extensibility of the leaves, whilst an increased amount of free water might be able to enhance solute accumulation, leading to better osmotic adjustment and tolerance to water stress, and maintenance of the volumes of sub-cellular compartments for expansive leaf growth. There were strong correlations between leaf-P%, leaf water (total, free and bound water) and leaf expansion rate (LER) under water stress conditions in a severely drying soil. Increased soil-P enhanced the uptake of P from a drying soil, leading to increased supply of osmotically active inorganic solutes to the cells in growing leaves. This appears to have led to the accumulation of free water and more bound water, ultimately leading to increased leaf expansion rates as compared to plants in low P soil under similar water stress conditions. The greater amount of bound and free water in the high-P plants was not necessarily associated with changes in cell turgor, and appears to have maintained the cell-wall properties and extensibility under water stressed conditions in soils that are nutritionally P-deficient.

Rapid Global Expansion of the Fungal Disease Chytridiomycosis into Declining and Healthy Amphibian Populations.

The fungal disease chytridiomycosis, caused by Batrachochytrium dendrobatidis, is enigmatic because it occurs globally in both declining and apparently healthy (non-declining) amphibian populations. This distribution has fueled debate concerning whether, in sites where it has recently been found, the pathogen was introduced or is endemic. In this study, we addressed the molecular population genetics of a global collection of fungal strains from both declining and healthy amphibian populations using DNA sequence variation from 17 nuclear loci and a large fragment from the mitochondrial genome. We found a low rate of DNA polymorphism, with only two sequence alleles detected at each locus, but a high diversity of diploid genotypes. Half of the loci displayed an excess of heterozygous genotypes, consistent with a primarily clonal mode of reproduction. Despite the absence of obvious sex, genotypic diversity was high (44 unique genotypes out of 59 strains). We provide evidence that the observed genotypic variation can be generated by loss of heterozygosity through mitotic recombination. One strain isolated from a bullfrog possessed as much allelic diversity as the entire global sample, suggesting the current epidemic can be traced back to the outbreak of a single clonal lineage. These data are consistent with the current chytridiomycosis epidemic resulting from a novel pathogen undergoing a rapid and recent range expansion. The widespread occurrence of the same lineage in both healthy and declining populations suggests that the outcome of the disease is contingent on environmental factors and host resistance.

Genetic effects and genetic relationships among shrunken (sh2) sweet corn lines and F1 hybrids

The objectives of this study were to quantify the components of genetic variance and the genetic effects, and to examine the genetic relationship of inbred lines extracted from various shrunken2(sh2) breeding populations. Ten diverse inbred lines developed from genetic background, were crossed in half diallel. Parents and their F1 hybrids were evaluated at three environments. The parents were genotyped using 20 polymorphic simple sequence repeats (SSR). Agronomic and quality traits were analysed by a mixed linear model according to additive-dominance genetic model. Genetic effects were estimated using an adjusted unbiased prediction method. Additive variance was more important than dominance variance in the expression of traits related to ear aspects (husk ratio and percentage of ear filled) and eating quality (flavour and total soluble solids). For agronomic traits, however, dominance variance was more important than additive variance. The additive genetic correlation between flavour and tenderness was strong (r = 0.84, P <0.01). Flavour, tenderness and kernel colour additive genetic effects were not correlated with yield related traits. Genetic distance (GD), estimated from SSR profiles on the basis of Jaccard's similarity coefficient varied from 0.10 to 0.77 with an average of 0.56. Cluster analysis classified parents according to their pedigree relationships. In most studied traits, F1 performance was not associated with GD.

Modelling temperature, photoperiod and vernalization responses of Brunonia australis (Goodeniaceae) and Calandrinia sp (Portulacaceae) to predict flowering time

Crop models for herbaceous ornamental species typically include functions for temperature and photoperiod responses, but very few incorporate vernalization, which is a requirement of many traditional crops. This study investigated the development of floriculture crop models, which describe temperature responses, plus photoperiod or vernalization requirements, using Australian native ephemerals Brunonia australis and Calandrinia sp. A novel approach involved the use of a field crop modelling tool, DEVEL2. This optimization program estimates the parameters of selected functions within the development rate models using an iterative process that minimizes sum of squares residual between estimated and observed days for the phenological event. Parameter profiling and jack-knifing are included in DEVEL2 to remove bias from parameter estimates and introduce rigour into the parameter selection process. Development rate of B. australis from planting to first visible floral bud (VFB) was predicted using a multiplicative approach with a curvilinear function to describe temperature responses and a broken linear function to explain photoperiod responses. A similar model was used to describe the development rate of Calandrinia sp., except the photoperiod function was replaced with an exponential vernalization function, which explained a facultative cold requirement and included a coefficient for determining the vernalization ceiling temperature. Temperature was the main environmental factor influencing development rate for VFB to anthesis of both species and was predicted using a linear model. The phenology models for B. australis and Calandrinia sp. described development rate from planting to VFB and from VFB to anthesis in response to temperature and photoperiod or vernalization and may assist modelling efforts of other herbaceous ornamental plants. In addition to crop management, the vernalization function could be used to identify plant communities most at risk from predicted increases in temperature due to global warming.

Estimating Hydrogen Cyanide in Forage Sorghum (Sorghum bicolor) by Near-Infrared Spectroscopy

Hydrogen cyanide (HCN) is a toxic chemical that can potentially cause mild to severe reactions in animals when grazing forage sorghum. Developing technologies to monitor the level of HCN in the growing crop would benefit graziers, so that they can move cattle into paddocks with acceptable levels of HCN. In this study, we developed near-infrared spectroscopy (MRS) calibrations to estimate HCN in forage sorghum and hay. The full spectral NIRS range (400-2498 nm) was used as well as specific spectral ranges within the full spectral range, i.e., visible (400-750 nm), shortwave (800-1100 nm) and near-infrared (NIR) (1100-2498 nm). Using the full spectrum approach and partial least-squares (PLS), the calibration produced a coefficient of determination (R-2) = 0.838 and standard error of cross-validation (SECV) = 0.040%, while the validation set had a R-2 = 0.824 with a low standard error of prediction (SEP = 0.047%). When using a multiple linear regression (MLR) approach, the best model (NIR spectra) produced a R-2 = 0.847 and standard error of calibration (SEC) = 0.050% and a R-2 = 0.829 and SEP = 0.057% for the validation set. The MLR models built from these spectral regions all used nine wavelengths. Two specific wavelengths 2034 and 2458 nm were of interest, with the former associated with C=O carbonyl stretch and the latter associated with C-N-C stretching. The most accurate PLS and MLR models produced a ratio of standard error of prediction to standard deviation of 3.4 and 3.0, respectively, suggesting that the calibrations could be used for screening breeding material. The results indicated that it should be feasible to develop calibrations using PLS or MLR models for a number of users, including breeding programs to screen for genotypes with low HCN, as well as graziers to monitor crop status to help with grazing efficiency.

Population Expansion and Genetic Structure in Carcharhinus brevipinna in the Southern Indo-Pacific

Background:Quantifying genetic diversity and metapopulation structure provides insights into the evolutionary history of a species and helps develop appropriate management strategies. We provide the first assessment of genetic structure in spinner sharks (Carcharhinus brevipinna), a large cosmopolitan carcharhinid, sampled from eastern and northern Australia and South Africa. Methods and Findings:Sequencing of the mitochondrial DNA NADH dehydrogenase subunit 4 gene for 430 individuals revealed 37 haplotypes and moderately high haplotype diversity (h = 0.6770 ±0.025). While two metrics of genetic divergence (ΦST and FST) revealed somewhat different results, subdivision was detected between South Africa and all Australian locations (pairwise ΦST, range 0.02717–0.03508, p values ≤ 0.0013; pairwise FST South Africa vs New South Wales = 0.04056, p = 0.0008). Evidence for fine-scale genetic structuring was also detected along Australia’s east coast (pairwise ΦST = 0.01328, p < 0.015), and between south-eastern and northern locations (pairwise ΦST = 0.00669, p < 0.04).Conclusions: The Indian Ocean represents a robust barrier to contemporary gene flow in C. brevipinna between Australia and South Africa. Gene flow also appears restricted along a continuous continental margin in this species, with data tentatively suggesting the delineation of two management units within Australian waters. Further sampling, however, is required for a more robust evaluation of the latter finding. Evidence indicates that all sampled populations were shaped by a substantial demographic expansion event, with the resultant high genetic diversity being cause for optimism when considering conservation of this commercially-targeted species in the southern Indo-Pacific.

The growth and productivity of ‘Festival’ strawberry plants growing in a subtropical environment

Experiments were conducted over 5 years to understand the seasonal phenology of bare-rooted ?Festival? strawberry plants (Fragaria ?ananassa) growing at Nambour in southeastern Queensland, Australia. Yields ranged from 661 to 966 g/plant, and average seasonal fruit fresh weight ranged from 15 to 18 g. The growth of the leaves, crowns, roots, flowers and fruit over time followed a linear or sigmoid pattern. Maximum values of leaf, crown and root dry weight towards the end of the growing season about 190 days after planting were 30, 15 and 7 g/plant, respectively. The rates of leaf and crown growth were lower than those achieved in California under a Mediterranean climate. There were strong relationships between the allocation of dry matter to the leaves, crowns and roots and plant dry weight. Allocation to the leaves, and especially to the crowns and roots, declined as the plants grew. The number of fruit/plant increased initially over time with a decline later in the season. Average fruit fresh weight was generally higher early in the season and then declined as fruit production increased. There were strong relationships between the growth of the whole plant and the growth of the flowers and immature fruit, and leaf expansion, across the growing season and across the 5 different years. These results indicate that seasonal growth and potential productivity were strongly linked to the expansion of the leaves in this environment.

Large-scale expansion of no-take closures within the Great Barrier Reef has not enhanced fishery production

A rare opportunity to test hypotheses about potential fishery benefits of large-scale closures was initiated in July 2004 when an additional 28.4% of the 348 000 km2 Great Barrier Reef (GBR) region of Queensland, Australia was closed to all fishing. Advice to the Australian and Queensland governments that supported this initiative predicted these additional closures would generate minimal (10%) initial reductions in both catch and landed value within the GBR area, with recovery of catches becoming apparent after three years. To test these predictions, commercial fisheries data from the GBR area and from the two adjacent (non-GBR) areas of Queensland were compared for the periods immediately before and after the closures were implemented. The observed means for total annual catch and value within the GBR declined from pre-closure (2000–2003) levels of 12 780 Mg and Australian $160 million, to initial post-closure (2005–2008) levels of 8143 Mg and $102 million; decreases of 35% and 36% respectively. Because the reference areas in the non-GBR had minimal changes in catch and value, the beyond-BACI (before, after, control, impact) analyses estimated initial net reductions within the GBR of 35% for both total catch and value. There was no evidence of recovery in total catch levels or any comparative improvement in catch rates within the GBR nine years after implementation. These results are not consistent with the advice to governments that the closures would have minimal initial impacts and rapidly generate benefits to fisheries in the GBR through increased juvenile recruitment and adult spillovers. Instead, the absence of evidence of recovery in catches to date currently supports an alternative hypothesis that where there is already effective fisheries management, the closing of areas to all fishing will generate reductions in overall catches similar to the percentage of the fished area that is closed.

Sirex woodwasp range expansion in Australia: performance and parasitism on two commercial pine species

Sirex woodwasp was detected in Queensland in 2009 and rapidly established in softwood plantations (Pinus radiata and P. taeda) in southern border regions. Biocontrol inoculations of Deladenus siricidicola began soon after, and adults were monitored to assess the success of the programme. Wasp size, sex ratios, emergence phenology and nematode parasitism rates were recorded, along with the assessment of wild-caught females. Patterns varied within and among seasons, but overall, P. taeda appeared to be a less suitable host than P. radiata, producing smaller adults, lower fat body content and fewer females. Sirex emerging from P. taeda also showed lower levels of nematode parasitism, possibly due to interactions with the more abundant blue-stain fungus in this host. Sirex adults generally emerged between November and March, with distinct peaks in January and March, separated by a marked drop in emergence in early February. Temperature provided the best correlate of seasonal emergence, with fortnights with higher mean minimum temperatures having higher numbers of Sirex emerging. This has implications for the anticipated northward spread of Sirex into sub-tropical coastal plantation regions. Following four seasons of inundative release of nematodes in Queensland, parasitism rates remain low and have resulted in only partial sterilization of infected females.