Rapid assessment of genotype-by-environment interactions and heritability for growth rate in aquaculture species using in vitro fertilisation and DNA tagging

Commercial environments may receive only a fraction of expected genetic gains for growth rate as predicted from the selection environment This fraction is the result of undesirable genotype-by-environment interactions (G x E) and measured by the genetic correlation (r(g)) of growth between environments. Rapid estimates of genetic correlation achieved in one generation are notoriously difficult to estimate with precision. A new design is proposed where genetic correlations can be estimated by utilising artificial mating from cryopreserved semen and unfertilised eggs stripped from a single female. We compare a traditional phenotype analysis of growth to a threshold model where only the largest fish are genotyped for sire identification. The threshold model was robust to differences in family mortality differing up to 30%. The design is unique as it negates potential re-ranking of families caused by an interaction between common maternal environmental effects and growing environment. The design is suitable for rapid assessment of G x E over one generation with a true 0.70 genetic correlation yielding standard errors as low as 0.07. Different design scenarios were tested for bias and accuracy with a range of heritability values, number of half-sib families created, number of progeny within each full-sib family, number of fish genotyped, number of fish stocked, differing family survival rates and at various simulated genetic correlation levels

Generalised growth models for aquatic species with an application to blacklip abalone (Haliotis rubra)

This paper presents a maximum likelihood method for estimating growth parameters for an aquatic species that incorporates growth covariates, and takes into consideration multiple tag-recapture data. Individual variability in asymptotic length, age-at-tagging, and measurement error are also considered in the model structure. Using distribution theory, the log-likelihood function is derived under a generalised framework for the von Bertalanffy and Gompertz growth models. Due to the generality of the derivation, covariate effects can be included for both models with seasonality and tagging effects investigated. Method robustness is established via comparison with the Fabens, improved Fabens, James and a non-linear mixed-effects growth models, with the maximum likelihood method performing the best. The method is illustrated further with an application to blacklip abalone (Haliotis rubra) for which a strong growth-retarding tagging effect that persisted for several months was detected

Subsampling multi-species trawl catches from tropical northern Australia: Does it matter which part of the catch is sampled?

Subsampling is a common method for estimating the abundance of species in trawl catches. However, the accuracy of subsampling in representing the total catch has not been assessed. To estimate one possible source of bias due to subsampling, we tested whether the position on trawler sorting trays from which subsamples were taken affected their ability to represent species in catches. This was done by sorting catches into 10 kg subsamples and comparing subsamples taken from different positions on the sorting tray. Comparisons were made after species were grouped into three categories of abundance, either 'rare', 'common' or 'abundant'. A generalised linear model analysis showed that taking subsamples from different positions on the sorting tray had no major effect on estimating the total numbers or weights of fish or invertebrates, or the total number of fish or invertebrate taxa, recorded in each position. Some individual taxa showed differences between positions on the sorting tray (11.5% of taxa ina three-position design; 25% in a five-position design). But consistent and meaningful patterns in the position of these taxa on the sorting tray could only be seen for the pony fish Leiognathus moretoniensis and the saucer scallop Amusium pleuronectes. Because most bycatch laxa are well mixed throughout the catch, subsamples can be taken from any position on trawler sorting trays without introducing bias.

Spatial and temporal distribution of tropical phytoplankton species and biomass in the Gulf of Carpenteria, Australia

The biomass and species composition of tropical phytoplankton in Albatross Bay, Gulf of Carpentaria, northern Australia, were examined monthly for 6 yr (1986 to 1992). Chlorophyll a (chl a) concentrations were highest (2 to 5.7 mu g l(-1)) in the wet season at inshore sites, usually coinciding with low salinities (30 to 33 ppt) and high temperatures (29 to 32 degrees C). At the offshore sites chi a concentrations were lower (0.2 to 2 mu g l(-1)) and did not vary seasonally. Nitrate and phosphate concentrations were generally low (0 to 3.68 mu M and 0.09 to 3 mu M for nitrate and phosphate respectively), whereas silicate was present in concentrations in the range 0.19 to 13 mu M. The phytoplankton community was dominated by diatoms, particularly at the inshore sites, as determined by a combination of microscopic and high-performance liquid chromatography (HPLC) pigment analyses. At the offshore sites the proportion of green flagellates increased. The cyanobacterium genus Trichodesmium and the diatom genera Chaetoceros, Rhizosolenia, Bacteriastrum and Thalassionema dominated the phytoplankton caught in 37 mu m mesh nets; however, in contrast to many other coastal areas studied worldwide there was no distinct species succession of the diatoms and only Trichodesmium showed seasonal changes in abundance. This reflects a stable phytoplankton community in waters without pulses of physical and chemical disturbances. These results are discussed in the context of the commercial prawn fishery in the Gulf of Carpentaria and the possible effect of phytoplankton on prawn larval growth and survival.

Inverse Correlation between Vitamin D and C-Reactive Protein in Newborns

Some studies suggested that adequate vitamin D might reduce inflammation in adults. However, little is known about this association in early life. We aimed to determine the relationship between cord blood 25-hydroxyvitamin D (25(OH)D) and C-reactive protein (CRP) in neonates. Cord blood levels of 25(OH)D and CRP were measured in 1491 neonates in Hefei, China. Potential confounders including maternal sociodemographic characteristics, perinatal health status, lifestyle, and birth outcomes were prospectively collected. The average values of cord blood 25(OH)D and CRP were 39.43 nmol/L (SD = 20.35) and 6.71 mg/L (SD = 3.07), respectively. Stratified by 25(OH)D levels, per 10 nmol/L increase in 25(OH)D, CRP decreased by 1.42 mg/L (95% CI: 0.90, 1.95) among neonates with 25(OH)D <25.0 nmol/L, and decreased by 0.49 mg/L (95% CI: 0.17, 0.80) among neonates with 25(OH)D between 25.0 nmol/L and 49.9 nmol/L, after adjusting for potential confounders. However, no significant association between 25(OH)D and CRP was observed among neonates with 25(OH)D ≥50 nmol/L. Cord blood 25(OH)D and CRP levels showed a significant seasonal trend with lower 25(OH)D and higher CRP during winter-spring than summer-autumn. Stratified by season, a significant linear association of 25(OH)D with CRP was observed in neonates born in winter-spring (adjusted β = −0.11, 95% CI: −0.13, −0.10), but not summer-autumn. Among neonates born in winter-spring, neonates with 25(OH)D <25 nmol/L had higher risk of CRP ≥10 mg/L (adjusted OR = 3.06, 95% CI: 2.00, 4.69), compared to neonates with 25(OH)D ≥25 nmol/L. Neonates with vitamin D deficiency had higher risk of exposure to elevated inflammation at birth.

Breastmilk-saliva interactions boost innate immunity by regulating the oral microbiome in early infancy

Introduction Xanthine oxidase (XO) is distributed in mammals largely in the liver and small intestine, but also is highly active in milk where it generates hydrogen peroxide (H2O2). Adult human saliva is low in hypoxanthine and xanthine, the substrates of XO, and high in the lactoperoxidase substrate thiocyanate, but saliva of neonates has not been examined. Results Median concentrations of hypoxanthine and xanthine in neonatal saliva (27 and 19 μM respectively) were ten-fold higher than in adult saliva (2.1 and 1.7 μM). Fresh breastmilk contained 27.3±12.2 μM H2O2 but mixing baby saliva with breastmilk additionally generated >40 μM H2O2, sufficient to inhibit growth of the opportunistic pathogens Staphylococcus aureus and Salmonella spp. Oral peroxidase activity in neonatal saliva was variable but low (median 7 U/L, range 2–449) compared to adults (620 U/L, 48–1348), while peroxidase substrate thiocyanate in neonatal saliva was surprisingly high. Baby but not adult saliva also contained nucleosides and nucleobases that encouraged growth of the commensal bacteria Lactobacillus, but inhibited opportunistic pathogens; these nucleosides/bases may also promote growth of immature gut cells. Transition from neonatal to adult saliva pattern occurred during the weaning period. A survey of saliva from domesticated mammals revealed wide variation in nucleoside/base patterns. Discussion and Conclusion During breast-feeding, baby saliva reacts with breastmilk to produce reactive oxygen species, while simultaneously providing growth-promoting nucleotide precursors. Milk thus plays more than a simply nutritional role in mammals, interacting with infant saliva to produce a potent combination of stimulatory and inhibitory metabolites that regulate early oral–and hence gut–microbiota. Consequently, milk-saliva mixing appears to represent unique biochemical synergism which boosts early innate immunity.

Resolving cryptic species complexes of major tephritid pests

An FAO/IAEA Co-ordinated Research Project (CRP) on “Resolution of Cryptic Species Complexes of Tephritid Pests to Overcome Constraints to SIT Application and International Trade” was conducted from 2010 to 2015. As captured in the CRP title, the objective was to undertake targeted research into the systematics and diagnostics of taxonomically challenging fruit fly groups of economic importance. The scientific output was the accurate alignment of biological species with taxonomic names; which led to the applied outcome of assisting FAO and IAEA Member States in overcoming technical constraints to the application of the Sterile Insect Technique (SIT) against pest fruit flies and the facilitation of international agricultural trade. Close to 50 researchers from over 20 countries participated in the CRP, using coordinated, multidisciplinary research to address, within an integrative taxonomic framework, cryptic species complexes of major tephritid pests. The following progress was made for the four complexes selected and studied: Anastrepha fraterculus complex – Eight morphotypes and their geographic and ecological distributions in Latin America were defined. The morphotypes can be considered as distinct biological species on the basis of differences in karyotype, sexual incompatibility, post-mating isolation, cuticular hydrocarbon, pheromone, and molecular analyses. Discriminative taxonomic tools using linear and geometric morphometrics of both adult and larval morphology were developed for this complex. Bactrocera dorsalis complex – Based on genetic, cytogenetic, pheromonal, morphometric, and behavioural data, which showed no or only minor variation between the Asian/African pest fruit flies Bactrocera dorsalis, B. papayae, B. philippinensis and B. invadens, the latter three species were synonymized with B. dorsalis. Of the five target pest taxa studied, only B. dorsalis and B. carambolae remain as scientifically valid names. Molecular and pheromone markers are now available to distinguish B. dorsalis from B. carambolae. Ceratitis FAR Complex (C. fasciventris, C. anonae, C. rosa) – Morphology, morphometry, genetic, genomic, pheromone, cuticular hydrocarbon, ecology, behaviour, and developmental physiology data provide evidence for the existence of five different entities within this fruit fly complex from the African region. These are currently recognised as Ceratitis anonae, C. fasciventris (F1 and F2), C. rosa and a new species related to C. rosa (R2). The biological limits within C. fasciventris (i.e. F1 and F2) are not fully resolved. Microsatellites markers and morphological identification tools for the adult males of the five different FAR entities were developed based on male leg structures. Zeugodacus cucurbitae (formerly Bactrocera (Zeugodacus) cucurbitae) – Genetic variability was studied among melon fly populations throughout its geographic range in Africa and the Asia/Pacific region and found to be limited. Cross-mating studies indicated no incompatibility or sexual isolation. Host preference and genetic studies showed no evidence for the existence of host races. It was concluded that the melon fly does not represent a cryptic species complex, neither with regard to geographic distribution nor to host range. Nevertheless, the higher taxonomic classification under which this species had been placed, by the time the CRP was started, was found to be paraphyletic; as a result the subgenus Zeugodacus was elevated to genus level.

Resolution of Cryptic Species Complexes of Tephritid Pests to Enhance SIT Application and Facilitate International Trade

Tephritid fruit flies (Diptera: Tephritidae) are considered by far the most important group of horticultural pests worldwide. Female fruit flies lay eggs directly into ripening fruit, where the maggots feed causing fruit loss. Each and every continent is plagued by a number of fruit fly pests, both indigenous as well as invasive ones, causing tremendous economic losses. In addition to the direct losses through damage, they can negatively impact commodity trade through restrictions to market access. The quarantine and regulatory controls put in place to manage them are expensive, while the on-farm control costs and loss of crop affect the general well-being of growers. These constraints can have huge implications on loss in revenues and limitations to developing fruit and vegetable-based agroindustries in developing, emergent and developed nations. Because fruit flies are a global problem, the study of their biology and management requires significant international attention to overcome the hurdles they pose. The Joint Food and Agriculture Organisation / International Atomic Energy Agency (FAO/IAEA) Programme on Nuclear Techniques in Food and Agriculture has been on the foreground in assisting Member States in developing and validating environment-friendly fruit fly suppression systems to support viable fresh fruit and vegetable production and export industries. Such international attention has resulted in the successful development and validation of a Sterile Insect Technique (SIT) package for the Mediterranean fruit fly. Although demands for R&D support with respect to Mediterranean fruit fly are diminishing due to successful integration of this package into sustainable control programmes against this pest in many countries, there were increasing demands from Member States in Africa, Asia and Latin America, to address other major fruit fly pests and a related, but sometimes neglected issue of tephritid species complexes of economic importance. Any research, whether it is basic or applied, requires a taxonomic framework that provides reliable and universally recognized entities and names. Among the currently recognized major fruit fly pests, there are groups of species whose morphology is very similar or identical, but biologically they are distinct species. As such, some insect populations that are grouped taxonomically within the same pest species, display different biological and genetic traits and show reproductive isolation which suggest that they are different species. On the other hand, different species may have been taxonomically described, but there may be doubt as to whether they actually represent distinct biological species or merely geographical variants of the same species. This uncertain taxonomic status has practical implications on the effective development and use of the SIT against such complexes, particularly at the time of determining which species to mass-rear, and significantly affects international movement of fruit and vegetables through the establishment of trade barriers to important agricultural commodities which are hosts to these pest tephritid species...

Anticipatory and reactive response to falls: Muscle synergy activation of forearm muscles

We investigated the surface electromyogram response of six forearm muscles to falls onto the outstretched hand. The extensor carpi radialis longus, extensor carpi radialis brevis, extensor carpi ulnaris, abductor pollicis longus, flexor carpi radialis and flexor carpi ulnaris muscles were sampled from eight volunteers who underwent ten self-initiated falls. All muscles initiated prior to impact. Co-contraction is the most obvious surface electromyogram feature. The predominant response is in the radial deviators. The surface electromyogram timing we recorded would appear to be a complex anticipatory response to falling modified by the ef- fect on the forearm muscles following impact. The mitigation of the force of impact is probably more importantly through shoulder abduction and extension and elbow flexion rather than action of the forearm muscles.

Delineating environmental envelopes to improve mapping of species distributions, via a hurdle model with CART &/or MaxEnT

Species distribution modelling (SDM) typically analyses species’ presence together with some form of absence information. Ideally absences comprise observations or are inferred from comprehensive sampling. When such information is not available, then pseudo-absences are often generated from the background locations within the study region of interest containing the presences, or else absence is implied through the comparison of presences to the whole study region, e.g. as is the case in Maximum Entropy (MaxEnt) or Poisson point process modelling. However, the choice of which absence information to include can be both challenging and highly influential on SDM predictions (e.g. Oksanen and Minchin, 2002). In practice, the use of pseudo- or implied absences often leads to an imbalance where absences far outnumber presences. This leaves analysis highly susceptible to ‘naughty-noughts’: absences that occur beyond the envelope of the species, which can exert strong influence on the model and its predictions (Austin and Meyers, 1996). Also known as ‘excess zeros’, naughty noughts can be estimated via an overall proportion in simple hurdle or mixture models (Martin et al., 2005). However, absences, especially those that occur beyond the species envelope, can often be more diverse than presences. Here we consider an extension to excess zero models. The two-staged approach first exploits the compartmentalisation provided by classification trees (CTs) (as in O’Leary, 2008) to identify multiple sources of naughty noughts and simultaneously delineate several species envelopes. Then SDMs can be fit separately within each envelope, and for this stage, we examine both CTs (as in Falk et al., 2014) and the popular MaxEnt (Elith et al., 2006). We introduce a wider range of model performance measures to improve treatment of naughty noughts in SDM. We retain an overall measure of model performance, the area under the curve (AUC) of the Receiver-Operating Curve (ROC), but focus on its constituent measures of false negative rate (FNR) and false positive rate (FPR), and how these relate to the threshold in the predicted probability of presence that delimits predicted presence from absence. We also propose error rates more relevant to users of predictions: false omission rate (FOR), the chance that a predicted absence corresponds to (and hence wastes) an observed presence, and the false discovery rate (FDR), reflecting those predicted (or potential) presences that correspond to absence. A high FDR may be desirable since it could help target future search efforts, whereas zero or low FOR is desirable since it indicates none of the (often valuable) presences have been ignored in the SDM. For illustration, we chose Bradypus variegatus, a species that has previously been published as an exemplar species for MaxEnt, proposed by Phillips et al. (2006). We used CTs to increasingly refine the species envelope, starting with the whole study region (E0), eliminating more and more potential naughty noughts (E1–E3). When combined with an SDM fit within the species envelope, the best CT SDM had similar AUC and FPR to the best MaxEnt SDM, but otherwise performed better. The FNR and FOR were greatly reduced, suggesting that CTs handle absences better. Interestingly, MaxEnt predictions showed low discriminatory performance, with the most common predicted probability of presence being in the same range (0.00-0.20) for both true absences and presences. In summary, this example shows that SDMs can be improved by introducing an initial hurdle to identify naughty noughts and partition the envelope before applying SDMs. This improvement was barely detectable via AUC and FPR yet visible in FOR, FNR, and the comparison of predicted probability of presence distribution for pres/absence.

Analytical model of reactive transport processes with spatially variable coefficients

Analytical solutions of partial differential equation (PDE) models describing reactive transport phenomena in saturated porous media are often used as screening tools to provide insight into contaminant fate and transport processes. While many practical modelling scenarios involve spatially variable coefficients, such as spatially variable flow velocity, v(x), or spatially variable decay rate, k(x), most analytical models deal with constant coefficients. Here we present a framework for constructing exact solutions of PDE models of reactive transport. Our approach is relevant for advection-dominant problems, and is based on a regular perturbation technique. We present a description of the solution technique for a range of one-dimensional scenarios involving constant and variable coefficients, and we show that the solutions compare well with numerical approximations. Our general approach applies to a range of initial conditions and various forms of v(x) and k(x). Instead of simply documenting specific solutions for particular cases, we present a symbolic worksheet, as supplementary material, which enables the solution to be evaluated for different choices of the initial condition, v(x) and k(x). We also discuss how the technique generalizes to apply to models of coupled multispecies reactive transport as well as higher dimensional problems.

Eliciting and encoding expert knowledge on variable selection into classical or Bayesian species distribution models

The quality of species distribution models (SDMs) relies to a large degree on the quality of the input data, from bioclimatic indices to environmental and habitat descriptors (Austin, 2002). Recent reviews of SDM techniques, have sought to optimize predictive performance e.g. Elith et al., 2006. In general SDMs employ one of three approaches to variable selection. The simplest approach relies on the expert to select the variables, as in environmental niche models Nix, 1986 or a generalized linear model without variable selection (Miller and Franklin, 2002). A second approach explicitly incorporates variable selection into model fitting, which allows examination of particular combinations of variables. Examples include generalized linear or additive models with variable selection (Hastie et al. 2002); or classification trees with complexity or model based pruning (Breiman et al., 1984, Zeileis, 2008). A third approach uses model averaging, to summarize the overall contribution of a variable, without considering particular combinations. Examples include neural networks, boosted or bagged regression trees and Maximum Entropy as compared in Elith et al. 2006. Typically, users of SDMs will either consider a small number of variable sets, via the first approach, or else supply all of the candidate variables (often numbering more than a hundred) to the second or third approaches. Bayesian SDMs exist, with several methods for eliciting and encoding priors on model parameters (see review in Low Choy et al. 2010). However few methods have been published for informative variable selection; one example is Bayesian trees (O’Leary 2008). Here we report an elicitation protocol that helps makes explicit a priori expert judgements on the quality of candidate variables. This protocol can be flexibly applied to any of the three approaches to variable selection, described above, Bayesian or otherwise. We demonstrate how this information can be obtained then used to guide variable selection in classical or machine learning SDMs, or to define priors within Bayesian SDMs.

Intracellular effects of atmospheric-pressure plasmas on melanoma cancer cells

Gas discharge plasmas formed at atmospheric pressure and near room temperature have recently been shown as a promising tool for cancer treatment. The mechanism of the plasma action is attributed to generation of reactive oxygen and nitrogen species, electric fields, charges, and photons. The relative importance of different modes of action of atmospheric-pressure plasmas depends on the process parameters and specific treatment objects. Hence, an in-depth understanding of biological mechanisms that underpin plasma-induced death in cancer cells is required to optimise plasma processing conditions. Here, the intracellular factors involved in the observed anti-cancer activity in melanoma Mel007 cells are studied, focusing on the effect of the plasma treatment dose on the expression of tumour suppressor protein TP73. Over-expression of TP73 causes cell growth arrest and/or apoptosis, and hence can potentially be targeted to enhance killing efficacy and selectivity of the plasma treatment. It is shown that the plasma treatment induces dose-dependent up-regulation of TP73 gene expression, resulting in significantly elevated levels of TP73 RNA and protein in plasma-treated melanoma cells. Silencing of TP73 expression by means of RNA interference inhibited the anticancer effects of the plasma, similar to the effect of caspase inhibitor z-VAD or ROS scavenger N-acetyl cysteine. These results confirm the role of TP73 protein in dose-dependent regulation of anticancer activity of atmospheric-pressure plasmas.

Pro-apoptotic NOXA is implicated in atmospheric-pressure plasma-induced melanoma cell death

Atmospheric-pressure plasma (APP) has been successfully used to treat several types of cancers in vivo and in vitro, with the effect being primarily attributed to the generation of reactive oxygen species (ROS). However, the mechanisms by which APP induces apoptosis in cancer cells require further elucidation. In this study, the effects of APP on the expression of 500 genes in melanoma Mel007 cancer cells were examined. Pro-apoptotic phorbol-12-myristate-13-acetate-induced protein (PMAIP1), also known as NOXA, was highly expressed as a result of APP treatment in a dose-dependent manner. Blocking of ROS using scavenger NAC or silencing of NOXA gene by RNA interference inhibited the APP-induced NOXA genes upregulation and impaired caspases 3/7 mediated apoptosis, confirming the important role plasma-generated ROS species and pro-apoptotic NOXA play in APP-induced cancer cell death.

Multi-species sequence comparison reveals conservation of ghrelin gene-derived splice variants encoding a truncated ghrelin peptide

The peptide hormone ghrelin is a potent orexigen produced predominantly in the stomach. It has a number of other biological actions, including roles in appetite stimulation, energy balance, the stimulation of growth hormone release and the regulation of cell proliferation. Recently, several ghrelin gene splice variants have been described. Here, we attempted to identify conserved alternative splicing of the ghrelin gene by cross-species sequence comparisons. We identified a novel human exon 2-deleted variant and provide preliminary evidence that this splice variant and in1-ghrelin encode a C-terminally truncated form of the ghrelin peptide, termed minighrelin. These variants are expressed in humans and mice, demonstrating conservation of alternative splicing spanning 90 million years. Minighrelin appears to have similar actions to full-length ghrelin, as treatment with exogenous minighrelin peptide stimulates appetite and feeding in mice. Forced expression of the exon 2-deleted preproghrelin variant mirrors the effect of the canonical preproghrelin, stimulating cell proliferation and migration in the PC3 prostate cancer cell line. This is the first study to characterise an exon 2-deleted preproghrelin variant and to demonstrate sequence conservation of ghrelin gene-derived splice variants that encode a truncated ghrelin peptide. This adds further impetus for studies into the alternative splicing of the ghrelin gene and the function of novel ghrelin peptides in vertebrates.