Extremely Low Microsatellite Diversity but Distinct Population Structure in a Long-Lived Threatened Species, the Australian Lungfish Neoceratodus forsteri (Dipnoi)

The Australian lungfish is a unique living representative of an ancient dipnoan lineage, listed as ‘vulnerable’ to extinction under Australia’s Environment Protection and Biodiversity Conservation Act 1999. Historical accounts indicate this species occurred naturally in two adjacent river systems in Australia, the Burnett and Mary. Current day populations in other rivers are thought to have arisen by translocation from these source populations. Early genetic work detected very little variation and so had limited power to answer questions relevant for management including how genetic variation is partitioned within and among sub-populations. In this study, we use newly developed microsatellite markers to examine samples from the Burnett and Mary Rivers, as well as from two populations thought to be of translocated origin, Brisbane and North Pine. We test whether there is significant genetic structure among and within river drainages; assign putatively translocated populations to potential source populations; and estimate effective population sizes. Eleven polymorphic microsatellite loci genotyped in 218 individuals gave an average within-population heterozygosity of 0.39 which is low relative to other threatened taxa and for freshwater fishes in general. Based on F_ST values (average over loci = 0.11) and STRUCTURE analyses, we identify three distinct populations in the natural range, one in the Burnett and two distinct populations in the Mary. These analyses also support the hypothesis that the Mary River is the likely source of translocated populations in the Brisbane and North Pine rivers, which agrees with historical published records of a translocation event giving rise to these populations. We were unable to obtain bounded estimates of effective population size, as we have too few genotype combinations, although point estimates were low, ranging from 29 - 129. We recommend that, in order to preserve any local adaptation in the three distinct populations that they be managed separately.

Abundance indices for long-lived tropical snappers: estimating standardized catch rates from spatially and temporally coarse logbook data

The Northern Demersal Scalefish Fishery has historically comprised a small fleet (≤10 vessels year−1) operating over a relatively large area off the northwest coast of Australia. This multispecies fishery primarily harvests two species of snapper: goldband snapper, Pristipomoides multidens and red emperor, Lutjanus sebae. A key input to age-structured assessments of these stocks has been the annual time-series of the catch rate. We used an approach that combined Generalized Linear Models, spatio-temporal imputation, and computer-intensive methods to standardize the fishery catch rates and report uncertainty in the indices. These analyses, which represent one of the first attempts to standardize fish trap catch rates, were also augmented to gain additional insights into the effects of targeting, historical effort creep, and spatio-temporal resolution of catch and effort data on trap fishery dynamics. Results from monthly reported catches (i.e. 1993 on) were compared with those reported daily from more recently (i.e. 2008 on) enhanced catch and effort logbooks. Model effects of catches of one species on the catch rates of another became more conspicuous when the daily data were analysed and produced estimates with greater precision. The rate of putative effort creep estimated for standardized catch rates was much lower than estimated for nominal catch rates. These results therefore demonstrate how important additional insights into fishery and fish population dynamics can be elucidated from such “pre-assessment” analyses.

Pronounced genetic population structure in a potentially vagile fish species (Pristipomoides multidens, Teleostei; Perciformes; Lutjanidae) from the East Indies triangle

The East Indies triangle, bordered by the Phillipines, Malay Peninsula and New Guinea, has a high level of tropical marine species biodiversity. Pristipomoides multidens is a large, long-lived, fecund snapper species that is distributed throughout the East Indies and Indo-Pacific. Samples were analysed from central and eastern Indonesia and northern Australia to test for genetic discontinuities in population structure. Fish (n = 377) were collected from the Indonesian islands of Bali, Sumbawa, Flores, West Timor, Tanimbar and Tual along with 131 fish from two northern Australian locations (Arafura and Timor Seas) from a previous study. Genetic variation in the control region of the mitochondrial genome was assayed using restriction fragment length polymorphism and direct sequencing. Haplotype diversity was high (0.67-0.82), as was intraspecific sequence divergence (range 0-5.8%). FST between pairs of populations ranged from 0 to 0.2753. Genetic subdivision was apparent on a small spatial scale; FST was 0.16 over 191 km (Bali/Sumbawa) and 0.17 over 491 km (Bali/Flores). Constraints to dispersal that contribute to, and maintain, the observed degree of genetic subdivision are experienced presumably by all life history stages of this tropical marine finfish. The constraints may include (1) little or no movement of eggs or larvae, (2) little or no home range or migratory movement of adults and (3) loss of larval cohorts due to transport of larvae away from suitable habitat by prevailing currents

Recovery and succession in a multi-species tropical seagrass meadow following experimental disturbance: the role of sexual and asexual reproduction

Recolonisation and succession in a multi-species tropical seagrass meadow was examined by creating gaps (50×50 cm) in the meadow and manipulating the supply of sexual and asexual propagules. Measurements of leaf shoot density and estimates of above-ground biomass were conducted monthly to measure recovery of gaps between September 1995 and November 1997. Measurements of the seeds stored in the sediment (seed bank) and horizontal rhizome growth of colonising species were also conducted to determine their role in the recovery process. Asexual colonisation through horizontal rhizome growth from the surrounding meadow was the main mechanism for colonisation of gaps created in the meadow. The seed bank played no role in recolonisation of cleared plots. Total shoot density and above-ground biomass (all species pooled) of cleared plots recovered asexually to the level of the undisturbed controls in 10 and 7 months, respectively. There was some sexual recruitment into cleared plots where asexual colonisation was prevented but seagrass abundance (shoot density and biomass) did not reach the level of unmanipulated controls. Seagrass species did not appear to form seed banks despite some species being capable of producing long-lived seeds. The species composition of cleared plots remained different to the undisturbed controls throughout the 26-month experiment. Syringodium isoetifolium was a rapid asexual coloniser of disturbed plots and remained at higher abundances than in the control treatments for the duration of the study. S. isoetifolium had the fastest horizontal rhizome growth of species asexually colonising cleared plots (6.9 mm day−1). Halophila ovalis was the most successful sexual coloniser but was displaced by asexually colonising species. H. ovalis was the only species observed to produce fruits during the study. Small disturbances in the meadow led to long-term (>2 years) changes in community composition. This study demonstrated that succession in tropical seagrass communities was not a deterministic process. Variations in recovery observed for different tropical seagrass communities highlighted the importance of understanding life history characteristics of species within individual communities to effectively predict their response to disturbance. A reproductive strategy involving clonal growth and production of long-lived, locally dispersed seeds is suggested which may provide an evolutionary advantage to plants growing in tropical environments subject to temporally unpredictable major disturbances such as cyclones

Effect of microwave radiation on seed mortality of rubber vine (Cryptostegia grandiflora R.Br.), parthenium (Parthenium hysterophorous L.) and bellyache bush (Jatropha gossypiifolia L.)

A trial was undertaken to evaluate the effect of microwaves on seed mortality of three weed species. Seeds of rubber vine (Cryptostegia grandiflora R.Br.), parthenium (Parthenium hysterophorous L.) and bellyache bush (Jatropha gossypiifolia L.) were buried at six depths (0, 2.5, 5, 10, 20 and 40 cm) in coarse sand maintained at one of two moisture levels, oven dry or wet (field capacity), and then subjected to one of five microwave radiation durations of (0, 2, 4, 8 and 16 min). Significant interactions between soil moisture level, microwave radiation duration, seed burial depth and species were detected for mortality of seeds of all three species. Maximum seed mortality of rubber vine (88%), parthenium (67%) and bellyache bush (94%) occurred in wet soil irradiated for 16 min. Maximum seed mortality of rubber vine and bellyache bush seeds occurred in seeds buried at 2.5 cm depth whereas that of parthenium occurred in seeds buried at 10 cm depth. Maximum soil temperatures of 114.1 and 87.5°C in dry and wet soil respectively occurred at 2.5 cm depth following 16 min irradiation. Irrespective of the greater soil temperatures recorded in dry soil, irradiating seeds in wet soil generally increased seed mortality 2.9-fold compared with dry soil. Moisture content of wet soil averaged 5.7% compared with 0.1% for dry soil. Results suggest that microwave radiation has the potential to kill seeds located in the soil seed bank. However, many factors, including weed species susceptibility, determine the effectiveness of microwave radiation on buried seeds. Microwave radiation may be an alternative to conventional methods at rapidly depleting soil seed banks in the field, particularly in relatively wet soils that contain long lived weed seeds.

Fishery biology and management of Protonibea diacanthus (Sciaenidae) aggregations in far Northern Cape York Peninsula waters

The sciaenid Protonibea diacanthus is a large, long-lived predatory fish of inshore northern Australian waters, which forms annual aggregations that are fished extensively by traditional (subsistence) and recreational fishers. There are now widespread concerns that the resource is being overexploited. Indigenous fishers of the Cape York Northern Peninsula Area (NPA) relate that large adult fish (up to 1500 mm total length (TL)) made up the bulk of the catch from the sciaenid aggregations until about 1994. In contrast, sexually mature P. diacanthus comprised only a small component (12 fish out of 270=4.4%) examined in a 1999–2000 sampling programme that was biased towards the largest individuals available. At 790 mm TL, the minimum size at first maturity for female P. diacanthus in this study is much smaller than the 920 mm TL reported previously in Queensland waters. Developing ovaries were observed in specimens sampled from sciaenid aggregations which formed in NPA waters between May and September 2000. However, no fish with ripe or spent gonads were found in the study, so the current timing and location of the spawning season for P. diacanthus in the region remain unknown. Food items observed in the analysis of the diet of P. diacanthus from the NPA included a variety of teleosts and invertebrates. The range of animal taxa represented in the prey items support the description of an ‘opportunistic predator’ attributed to the species. In our sampling, the stomach contents of fish caught during the time of the aggregation events did not differ from those observed at other times of the year. A total of 114 P. diacanthus were tagged and released at aggregation sites during the study period, and 3 fish (2.6%) were subsequently recaptured. The low rate of tag returns from the wild stock tagging programme, both in this study (2.6%) and from recreational fisher tag/release programmes for the sciaenid elsewhere in Queensland (6.5%), were not explained by tag loss nor mortality, given the high retention rate of tags and the zero mortality seen in tank trials. In response to the biological findings from this study, indigenous community councils of the NPA imposed a 2-year fishing moratorium for P. diacanthus. Surveys at aggregation sites in 2002 and 2003 established that much larger fish (mean size 103.5 cm TL) were again present on the grounds, albeit in very low numbers. These recent preliminary results highlight the critical need for continued monitoring and management of the P. diacanthus fishery in the NPA, if prospects for resource recovery are to be realised. The NPA initiative has provided a rare opportunity to negotiate a co-management strategy, based on scientific data and traditional knowledge, for the recovery of a cultural and economically significant fished resource.

Stochastic demography and population dynamics in the red kangaroo Macropus rufus

1. Many organisms inhabit strongly fluctuating environments but their demography and population dynamics are often analysed using deterministic models and elasticity analysis, where elasticity is defined as the proportional change in population growth rate caused by a proportional change in a vital rate. Deterministic analyses may not necessarily be informative because large variation in a vital rate with a small deterministic elasticity may affect the population growth rate more than a small change in a less variable vital rate having high deterministic elasticity. 2. We analyse a stochastic environment model of the red kangaroo (Macropus rufus), a species inhabiting an environment characterized by unpredictable and highly variable rainfall, and calculate the elasticity of the stochastic growth rate with respect to the mean and variability in vital rates. 3. Juvenile survival is the most variable vital rate but a proportional change in the mean adult survival rate has a much stronger effect on the stochastic growth rate. 4. Even if changes in average rainfall have a larger impact on population growth rate, increased variability in rainfall may still be important also in long-lived species. The elasticity with respect to the standard deviation of rainfall is comparable to the mean elasticities of all vital rates but the survival in age class 3 because increased variation in rainfall affects both the mean and variability of vital rates. 5. Red kangaroos are harvested and, under the current rainfall pattern, an annual harvest fraction of c. 20% would yield a stochastic growth rate about unity. However, if average rainfall drops by more than c. 10%, any level of harvesting may be unsustainable, emphasizing the need for integrating climate change predictions in population management and increase our understanding of how environmental stochasticity translates into population growth rate.

Demography of three perennial grasses in a central Queensland eucalypt woodland.

The population dynamics of the palatable, perennial grasses Bothriochloa ewartiana (Domin) C.E.Hubb. (desert Mitchell grass), Chrysopogon fallax S.T.Blake (golden beard grass) and Heteropogon contortus (L.) P.Beauv. ex Roem. & Schult. (black speargrass), were studied in an extensive grazing study conducted in a eucalypt woodland within the Aristida-Bothriochloa pasture community in central Queensland between 1994 and 2000. Treatments were three grazing pressures based on light, medium and heavy utilisation of forage available at the end of summer and two timber treatments (trees intact and trees killed). Seasonal rainfall throughout this study was generally favourable for plant growth with no severe drought periods. Grazing pressure had a greater overall impact on plant dynamics than timber treatment, which had minimal impact. Grazing pressure had a large impact on H. contortus dynamics, an intermediate impact on B. ewartiana and no impact on C. fallax. Fluctuations in plant density of both B. ewartiana and C. fallax were small because both species were long lived with low levels of seedling recruitment and plant death, whereas fluctuations in H. contortus density were relatively high because of its relatively short life span and higher levels of both recruitment and death. Heavy grazing pressure increased the recruitment of B. ewartiana and H. contortus in some years but had no impact on that of C. fallax. Heavy grazing pressure reduced the survival of the original plants of both B. ewartiana and H.contortus but not of C. fallax. For H. contortus, the size of the original plants was larger where trees were killed than where trees were left intact and plants of the 1995 seedling cohort were larger in 1998 at heavy compared with those at light and medium grazing pressure. Grazing had a minor negative impact on the soil seed bank of H. contortus. Populations of all three species remained stable throughout this study, although the favourable seasonal rainfall experienced and the short duration of this study relative to the life span of these species may have masked longer term, deleterious impacts of heavy grazing pressure.

The genetic structure of Australasian green turtles (Chelonia mydas): Exploring the geographical scale of genetic exchange

Ecological and genetic studies of marine turtles generally support the hypothesis of natal homing, but leave open the question of the geographical scale of genetic exchange and the capacity of turtles to shift breeding sites. Here we combine analyses of mitochondrial DNA (mtDNA) variation and recapture data to assess the geographical scale of individual breeding populations and the distribution of such populations through Australasia. We conducted multiscale assessments of mtDNA variation among 714 samples from 27 green turtle rookeries and of adult female dispersal among nesting sites in eastern Australia. Many of these rookeries are on shelves that were flooded by rising sea levels less than 10 000 years (c. 450 generations) ago. Analyses of sequence variation among the mtDNA control region revealed 25 haplotypes, and their frequency distributions indicated 17 genetically distinct breeding stocks (Management Units) consisting either of individual rookeries or groups of rookeries in general that are separated by more than 500 km. The population structure inferred from mtDNA was consistent with the scale of movements observed in long-term mark-recapture studies of east Australian rookeries. Phylogenetic analysis of the haplotypes revealed five clades with significant partitioning of sequence diversity (Φ = 68.4) between Pacific Ocean and Southeast Asian/Indian Ocean rookeries. Isolation by distance was indicated for rookeries separated by up to 2000 km but explained only 12% of the genetic structure. The emerging general picture is one of dynamic population structure influenced by the capacity of females to relocate among proximal breeding sites, although this may be conditional on large population sizes as existed historically across this region.

Standard yet unusual mechanisms of long-distance dispersal: Seed dispersal of Corymbia torelliana by bees

Mellitochory, seed dispersal by bees, has been implicated in long-distance dispersal of the tropical rain forest tree, Corymbia torelliana (Myrtaceae). We examined natural and introduced populations of C. torelliana for 4 years to determine the species of bees that disperse seeds, and the extent and distance of seed dispersal. The mechanism of seed dispersal by bees was also investigated, including fruit traits that promote dispersal, foraging behaviour of bees at fruits, and the fate of seeds. The fruit structure of C. torelliana, with seed presented in a resin reward, is a unique trait that promotes seed dispersal by bees and often results in long-distance dispersal. We discovered that a guild of four species of stingless bees, Trigona carbonaria, T. clypearis, T. sapiens, and T. hockingsi, dispersed seeds of C. torelliana in its natural range. More than half of the nests found within 250 m of fruiting trees had evidence of seed transport. Seeds were transported minimum distances of 20-220 m by bees. Approximately 88% of seeds were dispersed by gravity but almost all fruits retained one or two seeds embedded in resin for bee dispersal. Bee foraging for resin peaked immediately after fruit opening and corresponded to a peak of seed dispersal at the hive. There were strong correlations between numbers of seeds brought in and taken out of each hive by bees (r = 0.753-0.992, P < 0.05), and germination rates were 95 ± 5%. These results showed that bee-transported seeds were effectively dispersed outside of the hive soon after release from fruits. Seed dispersal by bees is a non-standard dispersal mechanism for C. torelliana, as most seeds are dispersed by gravity before bees can enter fruits. However, many C. torelliana seeds are dispersed by bees, since seeds are retained in almost all fruits, and all of these are dispersed by bees.

The limit to wheat water-use efficiency in eastern Australia. II. Influence of rainfall patterns

We investigated the influence of rainfall patterns on the water-use efficiency of wheat in a transect between Horsham (36°S) and Emerald (23°S) in eastern Australia. Water-use efficiency was defined in terms of biomass and transpiration, WUEB/T, and grain yield and evapotranspiration, WUEY/ET. Our working hypothesis is that latitudinal trends in WUEY/ET of water-limited crops are the complex result of southward increasing WUEB/T and soil evaporation, and season-dependent trends in harvest index. Our approach included: (a) analysis of long-term records to establish latitudinal gradients of amount, seasonality, and size-structure of rainfall; and (b) modelling wheat development, growth, yield, water budget components, and derived variables including WUEB/T and WUEY/ET. Annual median rainfall declined from around 600 mm in northern locations to 380 mm in the south. Median seasonal rain (from sowing to harvest) doubled between Emerald and Horsham, whereas median off-season rainfall (harvest to sowing) ranged from 460 mm at Emerald to 156 mm at Horsham. The contribution of small events (≤ 5 mm) to seasonal rainfall was negligible at Emerald (median 15 mm) and substantial at Horsham (105 mm). Power law coefficients (τ), i.e. the slopes of the regression between size and number of events in a log-log scale, captured the latitudinal gradient characterised by an increasing dominance of small events from north to south during the growing season. Median modelled WUEB/T increased from 46 kg/ha.mm at Emerald to 73 kg/ha.mm at Horsham, in response to decreasing atmospheric demand. Median modelled soil evaporation during the growing season increased from 70 mm at Emerald to 172 mm at Horsham. This was explained by the size-structure of rainfall characterised with parameter τ, rather than by the total amount of rainfall. Median modelled harvest index ranged from 0.25 to 0.34 across locations, and had a season-dependent latitudinal pattern, i.e. it was greater in northern locations in dry seasons in association with wetter soil profiles at sowing. There was a season-dependent latitudinal pattern in modelled WUEY/ET. In drier seasons, high soil evaporation driven by a very strong dominance of small events, and lower harvest index override the putative advantage of low atmospheric demand and associated higher WUEB/T in southern locations, hence the significant southwards decrease in WUEY/ET. In wetter seasons, when large events contribute a significant proportion of seasonal rain, higher WUEB/T in southern locations may translate into high WUEY/ET. Linear boundary functions (French-Schultz type models) accounting for latitudinal gradients in its parameters, slope, and x-intercept, were fitted to scatter-plots of modelled yield v. evapotranspiration. The x-intercept of the model is re-interpreted in terms of rainfall size structure, and the slope or efficiency multiplier is described in terms of the radiation, temperature, and air humidity properties of the environment. Implications for crop management and breeding are discussed.

Temporal genetic structure patterns in tropical maize populations under reciprocal recurrent selection.

In order to investigate the effect of long term recurrent selection on the pattern of gene diversity, thirty randomly-selected individuals from the progenitors (p) and four selection cycles (C0, C3, C6 and C11) were sampled for DNA analysis from the tropical maize (Zea mays L.) breeding populations, Atherton 1 (AT1) and Atherton 2 (AT2). Fifteen polymorphic Simple Sequence Repeat markers amplified a total of 284 and 257 alleles in AT1 and AT2 populations, respectively. Reductions in the number of alleles were observed at advanced selection cycles. About 11 and 12% of the alleles in AT1 and AT2 populations respectively, were near to fixation. However, a higher number of alleles (37% in AT1 and 33% in AT2) were close to extinction. Fisher's exact test and analysis of molecular variance (AMOVA) showed significant population differentiations. Gene diversity estimates and AMOVA revealed increased genetic differentiations at the expense of loss of heterozygosity. Population differentiations were mainly due to fixation of complementary alleles at a locus in the two breeding populations. The estimates of effective population at an advanced selection cycle were close to the population size predicted by the breeding method.

Stock structure of exploited shark species in north-eastern Australia

The project has provided management and other stakeholders with information necessary to make informed decisions about the management of four of the key exploited shark species caught in the Queensland inshore net fishery and northern New South Wales line fishery. The project has determined that spatial management of milk sharks within Queensland, and scalloped hammerhead, common black tip and Australian black tip sharks within Queensland and New South Wales is appropriate. The project has determined that both black tip shark species are likely to require co-operative management arrangements between Queensland and New South Wales. For scalloped hammerheads separate stocks between the two jurisdictions were identified from the fisheriesdependent samples, however genetic exchange across borders is likely to be facilitated by movement of adult females and perhaps larger males to a lesser extent. This information will greatly assist compliance with the Commonwealth Environment Protection and Biodiversity Conservation Act (1999) for shark fisheries in north-eastern Australia by providing the necessary basis for robust assessment of the status of stocks of the study species, thereby helping to deliver their sustainable harvest. It also helps to achieve objectives of the Australian National Shark Plan. The project provides the appropriate spatial framework for future monitoring and assessment of the study species. This is at a time when shark fisheries are receiving close attention from all sectors and when monitoring programs are being implemented, aimed at better assessment of stock status. This project has provided the crucial information for developing an appropriate monitoring design as well as the necessary basis for making statements about stock status. The project has addressed research priorities identified by the Queensland Fisheries Research Advisory Board, Great Barrier Reef Marine Park Authority and Queensland Fisheries. Previously management has assumed a single stock for each species on the east coast of Queensland, and management of shark fisheries in New South Wales (NSW) and Queensland has been independent of one another. The project has been able to enhance and develop links between research, management and industry. Strong positive relationships with commercial fishers were crucial in the collection of samples throughout the study area and fisheries managers were part of the project team throughout the study period. During the project the study area was extended to include both Queensland and NSW waters, creating mutualistic and positive links between the States’ research and management agencies. Extension of project results included management representatives from NSW and Queensland, as well as the Northern Territory where similar shark fisheries operate and similar species are targeted. The project was able to provide significant human capital development opportunities providing considerable value to the project outcomes. Use of vertebral microchemistry and life history characteristics as stock determination methods provided material for two PhD students based at James Cook University: Ron Schroeder, vertebral chemistry; and Alastair Harry, life history characteristic. The project has developed novel research methods that have great capacity for future application, including: • Development of a simple and rapid genetic diagnostic tool (RT-HRM-PCR assay) for differentiating among the black tip shark species, for which no simple morphological identifier exists; and • Development of laser ablation inductively coupled plasma mass spectrometry (LA-ICPMS) methods for analysing and interpreting microchemical composition of shark vertebrae. The study has provided further confirmation of the effectiveness of using a holistic approach in stock structure studies and justifies investment into such studies.

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.