Mixture of time-dependent growth models with an application to blue swimmer crab length-frequency data

Understanding how aquatic species grow is fundamental in fisheries because stock assessment often relies on growth dependent statistical models. Length-frequency-based methods become important when more applicable data for growth model estimation are either not available or very expensive. In this article, we develop a new framework for growth estimation from length-frequency data using a generalized von Bertalanffy growth model (VBGM) framework that allows for time-dependent covariates to be incorporated. A finite mixture of normal distributions is used to model the length-frequency cohorts of each month with the means constrained to follow a VBGM. The variances of the finite mixture components are constrained to be a function of mean length, reducing the number of parameters and allowing for an estimate of the variance at any length. To optimize the likelihood, we use a minorization–maximization (MM) algorithm with a Nelder–Mead sub-step. This work was motivated by the decline in catches of the blue swimmer crab (BSC) (Portunus armatus) off the east coast of Queensland, Australia. We test the method with a simulation study and then apply it to the BSC fishery data.

Molecular tools for studying puberty in the grey mullet, Mugil cephalus: cloning of cDNAs that regulate reproductive function

The cDNAs coding for the brain GnRHs (AY373449-51), pituitary GH, SL and PRL, and liver IGFs (AY427954-5) were isolated. Partial cDNA sequences of the brain (Cyp19b) and gonadal (Cyp19a) aromatases have also been obtained. These tools would be utilized to study the endocrine regulation of puberty in the grey mullet.

Tree species diversity and ecosystem function: Can tropical multi-species plantations generate greater productivity?

Results from the humid tropics of Australia demonstrate that diverse plantations can achieve greater productivity than monocultures. We found that increases in both the observed species number and the effective species richness were significantly related to increased levels of productivity as measured by stand basal area or mean individual tree basal area. Four of five plantation species were more productive in mixtures with other species than in monocultures, offering on average, a 55% increase in mean tree basal area. A general linear model suggests that species richness had a significant effect on mean individual tree basal area when environmental variables were included in the model. As monoculture plantations are currently the preferred reforestation method throughout the tropics these results suggest that significant productivity and ecological gains could be made if multi-species plantations are more broadly pursued.

Modelling predicts positive and negative interactions between three Australian tropical tree species in monoculture and binary mixture

Computer modelling promises to be an important tool for analysing and predicting interactions between trees within mixed species forest plantations. This study explored the use of an individual-based mechanistic model as a predictive tool for designing mixed species plantations of Australian tropical trees. The `spatially explicit individually based-forest simulator' (SeXI-FS) modelling system was used to describe the spatial interaction of individual tree crowns within a binary mixed-species experiment. The three-dimensional model was developed and verified with field data from three forest tree species grown in tropical Australia. The model predicted the interactions within monocultures and binary mixtures of Flindersia brayleyana, Eucalyptus pellita and Elaeocarpus grandis, accounting for an average of 42% of the growth variation exhibited by species in different treatments. The model requires only structural dimensions and shade tolerance as species parameters. By modelling interactions in existing tree mixtures, the model predicted both increases and reductions in the growth of mixtures (up to +/-50% of stem volume at 7 years) compared to monocultures. This modelling approach may be useful for designing mixed tree plantations.

Comparative productivity of irrigated short-term ryegrass (Lolium multiflorum) pasture receiving nitrogen, grown alone or in a mixture with white (Trifolium repens) and Persian (T. resupinatum) clovers

Dairy farms in subtropical Australia use irrigated, annually sown short-term ryegrass (Lolium multiflorum) or mixtures of short-term ryegrass and white (Trifolium repens) and Persian (shaftal) (T. resupinatum) clover during the winter-spring period in all-year-round milk production systems. A series of small plot cutting experiments was conducted in 3 dairying regions (tropical upland, north Queensland, and subtropical southeast Queensland and northern New South Wales) to determine the most effective rate and frequency of application of nitrogen (N) fertiliser. The experiments were not grazed, nor was harvested material returned to the plots, after sampling. Rates up to 100 kg N/ha.month (as urea or calcium ammonium nitrate) and up to 200 kg N/ha every 2 months (as urea) were applied to pure stands of ryegrass in 1991. In 1993 and 1994, urea, at rates up to 150 kg N/ha.month and to 200 kg N/ha every 2 months, was applied to pure stands of ryegrass; urea, at rates up to 50 kg N/ha.month, was also applied to ryegrass-clover mixtures. The results indicate that applications of 50-85 kg N/ha.month can be recommended for short-term ryegrass pastures throughout the subtropics and tropical uplands of eastern Australia, irrespective of soil type. At this rate, dry matter yields will reach about 90% of their potential, forage nitrogen concentration will be increased, there is minimal risk to stock from nitrate poisoning and there will be no substantial increase in soil N. The rate of N for ryegrass-clover pastures is slightly higher than for pure ryegrass but, at these rates, the clover component will be suppressed. However, increased ryegrass yields and higher forage nitrogen concentrations will compensate for the reduced clover component. At application rates up to 100 kg N/ha.month, build-up of NO3--N and NH4+-N in soil was generally restricted to the surface layers (0-20 cm) of the soil, but there was a substantial increase throughout the soil profile at 150 kg N/ha.month. The build-up of NO3--N and NH4+-N was greater and was found at lower rates on the lighter soil compared with heavy clays. Generally, most of the soil N was in the NO3--N form and most was in the top 20 cm.

The performance of irrigated mixtures of tall fescue, ryegrass and white clover in subtropical Australia. 1. The effects of sowing mixture combinations, nitrogen and oversowing on establishment, productivity, botanical composition and persistence.

In the subtropics of Australia, the ryegrass component of irrigated perennial ryegrass (Lolium perenne) - white clover (Trifolium repens) pastures declines by approximately 40% in the summer following establishment, being replaced by summer-active C4 grasses. Tall fescue (Festuca arundinacea) is more persistent than perennial ryegrass and might resist this invasion, although tall fescue does not compete vigorously as a seedling. This series of experiments investigated the influence of ryegrass and tall fescue genotype, sowing time and sowing mixture as a means of improving tall fescue establishment and the productivity and persistence of tall fescue, ryegrass and white clover-based mixtures in a subtropical environment. Tall fescue frequency at the end of the establishment year decreased as the number of companion species sown in the mixture increased. Neither sowing mixture combinations nor sowing rates influenced overall pasture yield (of around 14 t/ha) in the establishment year but had a significant effect on botanical composition and component yields. Perennial ryegrass was less competitive than short-rotation ryegrass, increasing first-year yields of tall fescue by 40% in one experiment and by 10% in another but total yield was unaffected. The higher establishment-year yield (3.5 t/ha) allowed Dovey tall fescue to compete more successfully with the remaining pasture components than Vulcan (1.4 t/ha). Sowing 2 ryegrass cultivars in the mixture reduced tall fescue yields by 30% compared with a single ryegrass (1.6 t/ha), although tall fescue alone achieved higher yields (7.1 t/ha). Component sowing rate had little influence on composition or yield. Oversowing the ryegrass component into a 6-week-old sward of tall fescue and white clover improved tall fescue, white clover and overall yields in the establishment year by 83, 17 and 11%, respectively, but reduced ryegrass yields by 40%. The inclusion of red (T. pratense) and Persian (T. resupinatum) clovers and chicory (Cichorium intybus) increased first-year yields by 25% but suppressed perennial grass and clover components. Yields were generally maintained at around 12 t/ha/yr in the second and third years, with tall fescue becoming dominant in all 3 experiments. The lower tall fescue seeding rate used in the first experiment resulted in tall fescue dominance in the second year following establishment, whereas in Experiments 2 and 3 dominance occurred by the end of the first year. Invasion by the C4 grasses was relatively minor (<10%) even in the third year. As ryegrass plants died, tall fescue and, to a lesser extent, white clover increased as a proportion of the total sward. Treatment effects continued into the second, but rarely the third, year and mostly affected the yield of one of the components rather than total cumulative yield. Once tall fescue became dominant, it was difficult to re-introduce other pasture components, even following removal of foliage and moderate renovation. Severe renovation (reducing the tall fescue population by at least 30%) seems a possible option for redressing this situation.

Study of the structure and function of the sheath of tropically adapted bulls and risk factors for preputial eversion and preputial prolapse

Preputial prolapse is an obvious condition affecting bulls from many breeds. Unfortunately, the losses in production and welfare concerns associated with preputial prolapse can remain undetected for long periods of time in the extensive beef areas of northern Australia where the bulls are not inspected regularly. Thus, there is a critical need to identify the structural factors predisposing to preputial prolapse in young bulls so that they can be culled early. Despite there being no firm scientific evidence of an association between preputial eversion and preputial prolapse, it seems logical that the increased exposure of the sensitive prepuce as a consequence of preputial eversion may increase the risk of bulls developing preputial pathology, in particular preputial prolapse. This may be particularly relevant in Bos indicus bulls as they have a more pendulous sheath and thus eversion of the prepuce may be associated with a greater risk of injury to the prepuce compared to that in Bos taurus bulls. Further, studies of preputial eversion in Bos taurus bulls have concluded that there is an association between polledness and increased prevalence and severity (length of everted prepuce and duration of eversion) of preputial eversion due primarily to the absence or poor development of the caudal preputial muscles. No similar definitive work in Bos indicus bulls has been conducted and thus anatomical studies reported in this thesis were conducted to determine if a similar association occurred in Bos indicus bulls. A survey of a sample of large beef breeding herds in northern Australia found that preputial prolapse is a significant problem in Bos indicus and Bos indicus derived bulls and affected both young and older bulls. The importance of preputial prolapse confirmed the value of further research into the causes of this problem. A series of anatomical studies confirmed that preputial eversion in Bos indicus derived bulls was not more prevalent in polled bulls than horned bulls and was not associated with deficiency of the caudal preputial muscles as was established in Bos taurus bulls. An anatomical study of Bos indicus derived bulls with preputial prolapse found that preputial prolapse occurred in horned bulls of varying ages and these bulls did not have any evidence of deficiency in the caudal preputial muscles. However, preputial prolapse was observed in young polled bulls that had poorly developed or absent caudal preputial muscles. It was concluded that deficiency of the caudal preputial muscles in polled Bos indicus derived bulls may predispose to preputial prolapse at an early age, but no predisposing anatomical factors were found for horned Bos indicus derived bulls. In these studies, preputial eversion and preputial prolapse were found in horned Bos indicus derived bulls that did not have any preputial muscle deficiency and it was noted that preputial eversion was not related to the length of the prepuce. Further studies confirmed that preputial eversion was linearly and consistently associated with position of the glans penis within the sheath in Bos indicus derived bulls, and movement of the glans penis towards the preputial orifice consistently resulted in preputial eversion in these bulls. A method to objectively measure the relationship between movement of the glans penis within the sheath and preputial eversion was developed. Studies in humans have linked function of some abdominal muscles to function of the pelvic organs. This relationship was investigated in Bos indicus derived bulls to determine whether the function of specific abdominal muscles affected position of the penis in the sheath. Using the method developed to objectively measure the relationship between penis movement and preputial eversion, the abdominal muscles that potentially were associated with movement of the glans penis or preputial eversion were examined but no significant relationships were observed. In the anatomical study of Bos indicus derived bulls not affected with preputial prolapse a more pendulous sheath was associated with increased prevalence of preputial eversion. This relationship was confirmed for horned and polled bulls in the penis movement studies. Bos indicus derived bulls with more pendulous sheaths evert their prepuces more than bulls with less pendulous sheaths thus increasing the risk of damage to the prepuce either from the environment, other bulls, or from them inadvertently stepping on the everted prepuce when they get to their feet. Culling Bos indicus derived bulls with more pendulous sheaths should reduce the incidence of preputial eversion and possibly preputial prolapse. The anatomical study of Bos indicus derived bulls that did not have preputial prolapse demonstrates that there are herds of bulls where the polled bulls do not have any evidence of deficiency of the caudal preputial iv muscles. There is a need to develop a practical and cost effective test to identify polled Bos indicus bulls that have a deficiency in their caudal preputial muscles.

Animal-Level Factors Affecting Ovarian Function in Bos indicus Heifers Treated to Synchronize Ovulation with Intravaginal Progesterone-Releasing Devices and Oestradiol Benzoate

The primary objective of this study was to investigate the impact of animal-level factors including energy balance and environmental/management stress, on the ovarian function of Bos indicus heifers treated to synchronize ovulation. Two-year-old Brahman (BN) (n = 30) and BN-cross (n = 34) heifers were randomly allocated to three intravaginal progesterone-releasing device (IPRD) treatment groups: (i) standard-dose IPRD [Cue-Mate (R) (CM) 1.56 g; n = 17]; (ii) half-dose IPRD [0.78 g progesterone (P4); CM 0.78 g; n = 15]; (iii) half-dose IPRD + 300 IU equine chorionic gonadotrophin at IPRD removal (CM 0.78 g + G; n = 14); (iv) and a control group, 2x PGF2a [500 mu g prostaglandin F2a (PGF2a)] on Day -16 and -2 (n = 18). Intravaginal progesterone-releasing device-treated heifers received 250 mu g PGF2a at IPRD insertion (Day -10) and IPRD removal (Day -2) and 1 mg oestradiol benzoate on Day -10 and -1. Heifers were managed in a small feedlot and fed a defined ration. Ovarian function was evaluated by ultrasonography and plasma P4 throughout the synchronized and return cycles. Energy balance was evaluated using plasma insulin-like growth factor 1 (IGF-I) and glucose concentrations. The impact of environmental stressors was evaluated using plasma cortisol concentration. Heifers that had normal ovarian function had significantly higher IGF-I concentrations at commencement of the experiment (p = 0.008) and significantly higher plasma glucose concentrations at Day -2 (p = 0.040) and Day 4 (p = 0.043), than heifers with abnormal ovarian function. There was no difference between the mean pre-ovulatory cortisol concentrations of heifers that ovulated or did not ovulate. However, heifers that ovulated had higher cortisol concentrations at Day 4 (p = 0.056) and 6 (p = 0.026) after ovulation than heifers that did not ovulate.

Mixture Effects of Benzene, Toluene, Ethylbenzene and Xylenes (BTEX) to Lung Carcinoma Cells via a Hanging Drop Air Exposure System

A recently developed hanging drop air exposure system for toxicity studies of volatile chemicals was applied to evaluate the cell viability of lung carcinoma A549 cells after 1 h and 24 h of exposure to benzene, toluene, ethylbenzene and xylenes (BTEX) as individual compounds and mixtures of 4 or 6 components. The cellular chemical concentrations causing 50% reduction of cell viability (EC50) were calculated use a mass balance model and came to 17, 12, 11, 9, 4 and 4 mmol/kg cell dry weight for benzene, toluene, ethylbenzene, m-xylene, o-xylene and p-xylene respectively after 1 h of exposure. The EC50 decreased by a factor of four after 24 h of exposure. All mixture effects were best described by the mixture toxicity model of concentration addition, which is valid for chemicals with the same mode of action. Good agreement with the model predictions were found for benzene, toluene, ethylbenzene and m-xylene at four different representative fixed concentration ratios after 1 h of exposure but lower agreement to mixture prediction was obtained after 24 h of exposure. A recreated car exhaust mixture, which involved the contribution of the more toxic p-xylene and o-xylene, yielded an acceptable but lower quality prediction as well.

Identifying the Function of Sorghum's Drought Tolerance Stay-Green QTL

The goal of this research is to understand the function of allelic variation of genes underpinning the stay-green drought adaptation trait in sorghum in order to enhance yield in water-limited environments. Stay-green, a delayed leaf senescence phenotype in sorghum, is primarily an emergent consequence of the improved balance between the supply and demand of water. Positional and functional fine-mapping of candidate genes associated with stay-green in sorghum is the focus of an international research partnership between Australian (UQ/DAFFQ) and US (Texas A&M University) scientists. Stay-green was initially mapped to four chromosomal regions (Stg1, Stg2, Stg3, and Stg4) by a number of research groups in the US and Australia. Physiological dissection of near-isolines containing single introgressions of Stg QTL (Stg1-4) indicate that these QTL reduce water demand before flowering by constricting the size of the canopy, thereby increasing water availability during grain filling and, ultimately, grain yield. Stg and root angle QTL are also co-located and, together with crop water use data, suggest the role of roots in the stay-green phenomenon. Candidate genes have been identified in Stg1-4, including genes from the PIN family of auxin efflux carriers in Stg1 and Stg2, with 10 of 11 PIN genes in sorghum co-locating with Stg QTL. Modified gene expression in some of these PIN candidates in the stay-green compared with the senescent types has been found in preliminary RNA expression profiling studies. Further proof-of-function studies are underway, including comparative genomics, SNP analysis to assess diversity at candidate genes, reverse genetics and transformation.

Hazard Function Models to Estimate Mortality Rates Affecting Fish Populations with Application to the Sea Mullet (Mugil cephalus) Fishery on the Queensland Coast (Australia)

Fisheries management agencies around the world collect age data for the purpose of assessing the status of natural resources in their jurisdiction. Estimates of mortality rates represent a key information to assess the sustainability of fish stocks exploitation. Contrary to medical research or manufacturing where survival analysis is routinely applied to estimate failure rates, survival analysis has seldom been applied in fisheries stock assessment despite similar purposes between these fields of applied statistics. In this paper, we developed hazard functions to model the dynamic of an exploited fish population. These functions were used to estimate all parameters necessary for stock assessment (including natural and fishing mortality rates as well as gear selectivity) by maximum likelihood using age data from a sample of catch. This novel application of survival analysis to fisheries stock assessment was tested by Monte Carlo simulations to assert that it provided unbiased estimations of relevant quantities. The method was applied to the data from the Queensland (Australia) sea mullet (Mugil cephalus) commercial fishery collected between 2007 and 2014. It provided, for the first time, an estimate of natural mortality affecting this stock: 0.22±0.08 year −1 .