Interleukin 10 polymorphisms in ankylosing spondylitis

Genetic polymorphisms of the IL10 promoter region have been implicated in many autoimmune diseases, including seronegative spondyloarthropathies. We studied three SNPs (IL10-1087,-824, and -597) and two microsafellites(IL10R and IL10G) lying within the promoter region of IL10 for association with susceptibility to and clinical manifestations of ankylosing spondylitis (AS), a common form of spondyloarthritis. Four hundred and sixty-eight individuals from 182 Finnish families affected with AS were studied. No association between individual IL10 promoter region polymorphisms or marker haplotype was observed with susceptibility to AS, but weak association was noted between the IL10-597 and -824 SNPs and age of disease onset (P= 0.01 for each SNP). The IL10.G4 allele was associated with BASFI (corrected for disease duration) (P= 0.03). We conclude that IL10 promoter polymorphisms have no significant effect on susceptibility to AS, but may play a minor role in determining age of disease onset and disease severity. © 2003 Nature Publishing Group All rights reserved.

Suggestive linkage of the parathyroid receptor type 1 to osteoporosis

We have investigated the role of 23 candidate genes in the control of bone mineral density (BMD) by linkage studies in families of probands with osteoporosis (lumbar spine [LS] or femoral neck [FN] BMD T score < -2.5) and low BMD relative to an age- and gender-matched cohort (Z score < -2.0). One hundred and fifteen probands (35 male, 80 female) and 499 of their first- or second-degree relatives (223 males and 276 females) were recruited for the study. BMD was measured at the LS and FN using dual-energy X-ray absorptiometry and expressed as age- and gender-matched Z scores corrected for body mass index. The candidate genes studied were the androgen receptor, type I collagen A1 (COLIA1), COLIA2, COLIIA1, vitamin D receptor (VDR), colony-stimulating factor 1, calcium-sensing receptor, epidermal growth factor (EGF), estrogen receptor 1 (ESR1), fibrillin type 1, insulin-like growth factor 1, interleukin-1 alpha (IL-1α), interleukin-4 (IL-4), interleukin-6 (IL-6), interleukin-11 (IL-11), osteopontin, parathyroid hormone (PTH), PTH-related peptide, PTH receptor type 1 (PTHR1), transforming growth factor-beta 1, and tumor necrosis factors alpha and beta. Sixty-four microsatellites lying close to or within these genes were investigated for linkage with BMD. Using the program MapMaker/Sibs there was suggestive evidence of linkage between BMD and PTHR1 (maximum LOD score obtained [MLS] 2.7-3.5). Moderate evidence of linkage was also observed with EGF (MLS 1.8), COLIA1 (MLS 1.7), COLIIA1/VDR (MLS 1.7), ESR1 (MLS 1.4), IL-1α (MLS 1.4), IL-4 (MLS 1.2), and IL-6 (MLS 1.2). Variance components analysis using the program ACT, correcting for proband-wise ascertainment, also showed evidence of linkage (p ≤0.05) at markers close to or within the candidate genes IL- 1α, PTHR1, IL-6, and COLIIA1/VDR. Further studies will be required to confirm these findings, to refine the location of gene responsible for the observed linkage, and to screen the candidate genes targeted at these loci for mutations.

A new data source for fisheries resource assessment: genetic estimates of the effective number of spawners. Final Report to the Fisheries Research and Development Corporation.

The development of innovative methods of stock assessment is a priority for State and Commonwealth fisheries agencies. It is driven by the need to facilitate sustainable exploitation of naturally occurring fisheries resources for the current and future economic, social and environmental well being of Australia. This project was initiated in this context and took advantage of considerable recent achievements in genomics that are shaping our comprehension of the DNA of humans and animals. The basic idea behind this project was that genetic estimates of effective population size, which can be made from empirical measurements of genetic drift, were equivalent to estimates of the successful number of spawners that is an important parameter in process of fisheries stock assessment. The broad objectives of this study were to 1. Critically evaluate a variety of mathematical methods of calculating effective spawner numbers (Ne) by a. conducting comprehensive computer simulations, and by b. analysis of empirical data collected from the Moreton Bay population of tiger prawns (P. esculentus). 2. Lay the groundwork for the application of the technology in the northern prawn fishery (NPF). 3. Produce software for the calculation of Ne, and to make it widely available. The project pulled together a range of mathematical models for estimating current effective population size from diverse sources. Some of them had been recently implemented with the latest statistical methods (eg. Bayesian framework Berthier, Beaumont et al. 2002), while others had lower profiles (eg. Pudovkin, Zaykin et al. 1996; Rousset and Raymond 1995). Computer code and later software with a user-friendly interface (NeEstimator) was produced to implement the methods. This was used as a basis for simulation experiments to evaluate the performance of the methods with an individual-based model of a prawn population. Following the guidelines suggested by computer simulations, the tiger prawn population in Moreton Bay (south-east Queensland) was sampled for genetic analysis with eight microsatellite loci in three successive spring spawning seasons in 2001, 2002 and 2003. As predicted by the simulations, the estimates had non-infinite upper confidence limits, which is a major achievement for the application of the method to a naturally-occurring, short generation, highly fecund invertebrate species. The genetic estimate of the number of successful spawners was around 1000 individuals in two consecutive years. This contrasts with about 500,000 prawns participating in spawning. It is not possible to distinguish successful from non-successful spawners so we suggest a high level of protection for the entire spawning population. We interpret the difference in numbers between successful and non-successful spawners as a large variation in the number of offspring per family that survive – a large number of families have no surviving offspring, while a few have a large number. We explored various ways in which Ne can be useful in fisheries management. It can be a surrogate for spawning population size, assuming the ratio between Ne and spawning population size has been previously calculated for that species. Alternatively, it can be a surrogate for recruitment, again assuming that the ratio between Ne and recruitment has been previously determined. The number of species that can be analysed in this way, however, is likely to be small because of species-specific life history requirements that need to be satisfied for accuracy. The most universal approach would be to integrate Ne with spawning stock-recruitment models, so that these models are more accurate when applied to fisheries populations. A pathway to achieve this was established in this project, which we predict will significantly improve fisheries sustainability in the future. Regardless of the success of integrating Ne into spawning stock-recruitment models, Ne could be used as a fisheries monitoring tool. Declines in spawning stock size or increases in natural or harvest mortality would be reflected by a decline in Ne. This would be good for data-poor fisheries and provides fishery independent information, however, we suggest a species-by-species approach. Some species may be too numerous or experiencing too much migration for the method to work. During the project two important theoretical studies of the simultaneous estimation of effective population size and migration were published (Vitalis and Couvet 2001b; Wang and Whitlock 2003). These methods, combined with collection of preliminary genetic data from the tiger prawn population in southern Gulf of Carpentaria population and a computer simulation study that evaluated the effect of differing reproductive strategies on genetic estimates, suggest that this technology could make an important contribution to the stock assessment process in the northern prawn fishery (NPF). Advances in the genomics world are rapid and already a cheaper, more reliable substitute for microsatellite loci in this technology is available. Digital data from single nucleotide polymorphisms (SNPs) are likely to super cede ‘analogue’ microsatellite data, making it cheaper and easier to apply the method to species with large population sizes.

The genetic effective and adult census size of an Australian population of tiger prawns (Penaeus esculentus)

This study compares estimates of the census size of the spawning population with genetic estimates of effective current and long-term population size for an abundant and commercially important marine invertebrate, the brown tiger prawn (Penaeus esculentus). Our aim was to focus on the relationship between genetic effective and census size that may provide a source of information for viability analyses of naturally occurring populations. Samples were taken in 2001, 2002 and 2003 from a population on the east coast of Australia and temporal allelic variation was measured at eight polymorphic microsatellite loci. Moments-based and maximum-likelihood estimates of current genetic effective population size ranged from 797 to 1304. The mean long-term genetic effective population size was 9968. Although small for a large population, the effective population size estimates were above the threshold where genetic diversity is lost at neutral alleles through drift or inbreeding. Simulation studies correctly predicted that under these experimental conditions the genetic estimates would have non-infinite upper confidence limits and revealed they might be overestimates of the true size. We also show that estimates of mortality and variance in family size may be derived from data on average fecundity, current genetic effective and census spawning population size, assuming effective population size is equivalent to the number of breeders. This work confirms that it is feasible to obtain accurate estimates of current genetic effective population size for abundant Type III species using existing genetic marker technology.

Population genetic structure of the brown tiger prawn, Penaeus esculentus, in tropical northern Australia

Eight polymorphic microsatellite loci were analysed in six population samples from four locations of the Australian endemic brown tiger prawn, Penaeus esculentus. Tests of Hardy-Weinberg equilibrium were generally in accord with expectations, with only one locus, in two samples, showing significant deviations. Three samples were taken in different years from the Exmouth Gulf. These showed no significant heterogeneity, and it was concluded that they were from a single panmictic population. A sample from Shark Bay, also on the west coast of Australia, showed barely detectable differentiation from Exmouth Gulf (F (ST) = 0 to 0.0014). A northeast sample from the Gulf of Carpentaria showed low (F (ST) = 0.008) but significant differentiation from Moreton Bay, on the east coast. However, Exmouth Gulf/Shark Bay samples were well differentiated from the Gulf of Carpentaria/Moreton Bay (F (ST) = 0.047-0.063). The data do not fit a simple isolation by distance model. It is postulated that the east-west differentiation largely reflects the isolation of east and west coast populations that occurred at the last glacial maximum when there was a land bridge between north-eastern Australia and New Guinea.

Confirmation of QTL mapping and marker validation for partial seedling resistance to crown rot in wheat line '2-49'

We have tested the efficacy of putative microsatellite single sequence repeat (SSR) markers, previously identified in a 2-49 (Gluyas Early/Gala) × Janz doubled haploid wheat (Triticum aestivum) population, as being linked to partial seedling resistance to crown rot disease caused by Fusarium pseudograminearum. The quantitative trait loci (QTLs) delineated by these markers have been tested for linkage to resistance in an independent Gluyas Early × Janz doubled haploid population. The presence of a major QTL on chromosome 1DL (QCr.usq-1D1) and a minor QTL on chromosome 2BS (QCr.usq-2B1) was confirmed. However, a putative minor QTL on chromosome 2A was not confirmed. The QTL on 1D was inherited from Gluyas Early, a direct parent of 2-49, whereas the 2B QTL was inherited from Janz. Three other putative QTLs identified in 2-49 × Janz (on 1AL, 4BL, and 7BS) were inherited by 2-49 from Gala and were not able to be confirmed in this study. The screening of SSR markers on a small sample of elite wheat genotypes indicated that not all of the most tightly linked SSR markers flanking the major QTLs on 1D and 1A were polymorphic in all backgrounds, indicating the need for additional flanking markers when backcrossing into some elite pedigrees. Comparison of SSR haplotypes with those of other genotypes exhibiting partial crown rot resistance suggests that additional, novel sources of crown rot resistance are available.

Genetic variation within two sympatric spotted gum eucalypts exceeds between taxa variation

Population substructure and hybridization, among other factors, have the potential to cause erroneous associations in linkage disequilibrium (LD) mapping. Two closely related spotted gum eucalypts, Corymbia variegata and C. henryi (Myrtaceae) occur in sympatry in the east coast of Australia and potentially interbreed. They are morphologically similar but are distinguished as separate species based on capsule and foliage size. To determine whether they hybridize in nature and its implications for LD mapping, we investigated the level of molecular divergence between the two species at two sympatric locations separated by 300 kilometres. Very few individuals of intermediate morphology were identified, despite the two species occurring only metres apart. Analysis of genetic structure using 12 microsatellite loci showed that genetic differentiation between populations of the same species at different locations (FST = 0.07 for both species; p = 0.0001) was significantly higher than that observed between species at each location (mean FST = 0.02 and 0.04 for Cherry tree and Bunyaville respectively; p = 0.0001; all Mann-Whitney U-test p ≤ 0.01). No species-specific alleles or significant allele frequency differences were detected within a site, suggesting recurrent local gene flow between the two species. The lack of significant allele frequency differences implies no population stratification along taxonomic lines. This suggested that there is little concern for cryptic hybridization when sampling from sites of sympatry for LD mapping.

Heritability estimates for growth in the tropical abalone Haliotis asinina using microsatellites to assign parentage

The tropical abalone Haliotis asinina is a wild-caught and cultured species throughout the Indo-Pacific as well as being an emerging model species for the study of haliotids. H. asinina has the fastest recorded natural growth rate of any abalone and reaches sexual maturity within one year. As such, it is a suitable abalone species for selective breeding for commercially important traits such as rapid growth. Estimating the amount of variation in size that is attributable to heritable genetic differences can assist the development of such a selective breeding program. Here we estimated heritability for growth-related traits at 12 months of age by creating a single cohort of 84 families in a full-factorial mating design consisting of 14 sires and 6 dams. Of 500 progeny sampled, 465 were successfully assigned to their parents based on shared alleles at 5 polymorphic microsatellite loci. Using an animal model, heritability estimates were 0.48 ± 0.15 for shell length, 0.38 ± 0.13 for shell width and 0.36 ± 0.13 for weight. Genetic correlations were > 0.98 between shell parameters and weight, indicating that breeding for weight gains could be successfully achieved by selecting for shell length.

Laboratory information management software for genotyping workflows: Applications in high throughput crop genotyping

Background: With the advances in DNA sequencer-based technologies, it has become possible to automate several steps of the genotyping process leading to increased throughput. To efficiently handle the large amounts of genotypic data generated and help with quality control, there is a strong need for a software system that can help with the tracking of samples and capture and management of data at different steps of the process. Such systems, while serving to manage the workflow precisely, also encourage good laboratory practice by standardizing protocols, recording and annotating data from every step of the workflow Results: A laboratory information management system (LIMS) has been designed and implemented at the International Crops Research Institute for the Semi-Arid Tropics (ICRISAT) that meets the requirements of a moderately high throughput molecular genotyping facility. The application is designed as modules and is simple to learn and use. The application leads the user through each step of the process from starting an experiment to the storing of output data from the genotype detection step with auto-binning of alleles; thus ensuring that every DNA sample is handled in an identical manner and all the necessary data are captured. The application keeps track of DNA samples and generated data. Data entry into the system is through the use of forms for file uploads. The LIMS provides functions to trace back to the electrophoresis gel files or sample source for any genotypic data and for repeating experiments. The LIMS is being presently used for the capture of high throughput SSR (simple-sequence repeat) genotyping data from the legume (chickpea, groundnut and pigeonpea) and cereal (sorghum and millets) crops of importance in the semi-arid tropics. Conclusions: A laboratory information management system is available that has been found useful in the management of microsatellite genotype data in a moderately high throughput genotyping laboratory. The application with source code is freely available for academic users and can be downloaded from http://www.icrisat.org/bt-software-d-lims.htm

ShadowBoxer and LocusEater: programs to optimise experimental design and multiplexing strategies for genetic mark-recapture

Genetic mark–recapture requires efficient methods of uniquely identifying individuals. 'Shadows' (individuals with the same genotype at the selected loci) become more likely with increasing sample size, and bias harvest rate estimates. Finding loci is costly, but better loci reduce analysis costs and improve power. Optimal microsatellite panels minimize shadows, but panel design is a complex optimization process. locuseater and shadowboxer permit power and cost analysis of this process and automate some aspects, by simulating the entire experiment from panel design to harvest rate estimation.

Identification of quantitative trait loci for resistance to two species of root-lesion nematode (Pratylenchus thornei and P. neglectus) in wheat

Pratylenchus thornei and P. neglectus are two species of root-lesion nematode that cause substantial yield losses in wheat. No commercially available wheat variety has resistance to both species. A doubled-haploid population developed from a cross between the synthetic hexaploid wheat line CPI133872 and the bread wheat Janz was used to locate and tag quantitative trait loci (QTLs) associated with resistance to both P. thornei and P. neglectus. Wheat plants were inoculated with both species of nematode in independent replicated glasshouse trials repeated over 2 years. Known locations of wheat microsatellite markers were used to construct a framework map. After an initial single-marker analysis to detect marker-trait linkages, chromosome regions associated with putative QTLs were targetted with microsatellite markers to increase map density in the chromosome regions of interest. In total, 148 wheat microsatellite markers and 21 amplified fragment length polymorphism markers were mapped. The codominant microsatellite marker Xbarc183 on the distal end of chromosome 6DS was allelic for resistance to both P. thornei and P. neglectus. The QTL were designated QRlnt.lrc-6D.1 and QRlnn.lrc-6D.1, for the 2 traits, respectively. The allele inherited from CPI133872 explained 22.0-24.2% of the phenotypic variation for P. thornei resistance, and the allele inherited from Janz accounted for 11.3-14.0% of the phenotypic variation for P. neglectus resistance. Composite interval mapping identified markers that flank a second major QTL on chromosome 6DL (QRlnt.lrc-6D.2) that explained 8.3-13.4% of the phenotypic variation for P. thornei resistance. An additional major QTL associated with P. neglectus resistance was detected on chromosome 4DS (QRlnn.lrc-4D.1) and explained a further 10.3-15.4% of the phenotypic variation. The identification and tagging of nematode resistance genes with molecular markers will allow appropriate allele combinations to be selected, which will aid the successful breeding of wheat with dual nematode resistance.

Two sympatric spotted gum species are molecularly homogeneous.

Large fruited spotted gum eucalypt Corymbia henryi occurs sympatrically with small fruited spotted gum Corymbia citriodora subspecies variegata over a large portion of its range on the east coast of Australia. The two taxa are interfertile, have overlapping flowering times and share a common set of insect and vertebrate pollinators. Previous genetic analysis of both taxa from two geographically remote sites suggested that the two were morphotypes rather than genetically distinct species. In this study we further explore this hypothesis of genic species by expanding sampling broadly through their sympatric locations and examine local-scale spatial genetic structure in stands that differ in species and age composition. Delineation of populations at five microsatellite loci, using an individual-based approach and Bayesian modelling, as well as clustering of individuals based on allele frequencies showed the two species to be molecularly homogeneous. Genetic structure aligned largely with geographic areas of origin, and followed an isolation-by-distance model, where proximal populations were generally less differentiated than more distant ones. At the stand level, spotted gums also generally showed little structure consistent with the high levels of gene flow inferred across the species range. Disturbances in the uniformity of structuring were detected, however, and attributed to localised events giving rise to even aged stands, probably due to regeneration from a few individuals following fire.

Reproductive Biology of Corymbia citriodora subsp. variegata and Effective Pollination across its Native Range.

The spotted gum species complex represents a group of four eucalypt hardwood taxa that have a native range that spans the east coast of Australia, with a morphological cline from Victoria to northern Queensland. Of this group, Corymbia citriodora subsp. variegata (CCV) is widespread in south-eastern Queensland and northern New South Wales. It is currently the most commonly harvested native hardwood in Queensland. However, little basic knowledge of the reproductive biology of the species is available to inform genetic improvement and resource management programmes. Here we take an integrative approach, using both field and molecular data, to identify ecological factors important to mating patterns in native populations of CCV. Field observation of pollinator visitation and flowering phenology of 20 trees showed that foraging behaviour of pollinator guilds varies depending on flowering phenology and canopy structure. A positive effect of tree mean flowering effort was found on insect visitation, while bat visitation was predicted by tree height and by the number of trees simultaneously bearing flowers. Moreover, introduced honeybees were observed frequently, performing 73% of detected flower visits. Conversely, nectar-feeding birds and mammals were observed sporadically with lorikeets and honeyeaters each contributing to 11% of visits. Fruit bats, represented solely by the grey-headed flying fox, performed less than 2% of visits. Genotyping at six microsatellite markers in 301 seeds from 17 families sampled from four of Queensland's native forests showed that CCV displays a mixed-mating system that is mostly outcrossing (tm = 0.899 ± 0.021). Preferential effective pollination from near-neighbours was detected by means of maximum-likelihood paternity analysis with up to 16% of reproduction events resulting from selfing. Forty to 48% of fertilising pollen was also carried from longer distance (>60 m). Marked differences in foraging behaviour and visitation frequency between observed pollinator guilds suggests that the observed dichotomy of effective pollen movement in spotted gums may be due to frequent visit from introduced honeybees favouring geitonogamy and sporadic visits from honeyeaters and fruit bats resulting in potential long-distance pollinations.

Immunomodulatory effects of probiotic bacteria in healthy adults

Probiooteilla kantakohtaisia vaikutuksia ihmisen immuunijärjestelmään terveillä aikuisilla Probiooteilla on kantakohtaisia tulehduksen välittäjäaineita vähentäviä vaikutuksia ja probioottien yhdistelmien vaikutukset eroavat yksittäisten kantojen vaikutuksista selviää TtM Riina Kekkosen tuoreesta väitöstutkimuksesta. TtM Riina Kekkonen on selvittänyt väitöskirjassaan eri probioottikantojen vaikutuksia immuunivasteeseen valkosolumallissa sekä terveillä aikuisilla lumekontrolloiduissa kliinisissä tutkimuksissa. Aikaisemmin probioottien vaikutuksia on tutkittu lähinnä allergian ja erilaisten vatsavaivojen ehkäisyssä ja hoidossa. Probiootteja sisältäviä tuotteita käyttävät kuluttajat ovat kuitenkin useimmiten terveitä aikuisia, ja probioottien vaikutus terveiden aikuisten immuunijärjestelmään on ollut puutteellisesti selvitettyä. Valkosolumallissa probioottikantojen havaittiin poikkeavan toisistaan niiden kyvyssä aktivoida immuunivasteen välittäjäaineiden, sytokiinien, tuotantoa. Anti-inflammatorisia, eli tulehdusta lievittäviä vaikutuksia nähtiin lähinnä Bifidobacterium ja Propionibacterium sukuihin kuuluvilla kannoilla. Streptococcus ja Leuconostoc sukuihin kuuluvat kannat puolestaan aktivoivat Th1 tyyppistä, soluvälitteistä immuunivastetta. Eri probioottien kombinaatiot eivät saaneet aikaan voimakkaampaa aktivaatiota yksittäisiin kantoihin verrattuna, joka viittaa probioottien keskinäiseen kilpailuun niiden ollessa kontaktissa ihmisen solujen kanssa. Probioottikantojen valinta kliinisiin tutkimuksiin tehtiin niiden anti-inflammatoristen ominaisuuksien perusteella. Parhaita anti-inflammatorisia kantoja olivat B. lactis ssp. animalis Bb12 ja P. freudenreichii ssp. shermanii JS, joiden lisäksi tutkimuksiin valittiin myös L. rhamnosus GG (LGG) hyvin tutkittuna referenssikantana. Solutöiden tulokset eivät olleet täysin verrannollisia kliinisen työn tuloksiin, koska LGG näytti omaavan parhaat anti-inflammatoriset ominaisuudet kliinisissä tutkimuksissa vaikka solutyössä sen aikaansaamat vasteet olivat melko vaimeita. Kolmen viikon kliinisessä tutkimuksessa terveillä aikuisilla LGG alensi mm. tulehdusta kuvaavan C-reaktiivisen proteiinin ja inflammatoristen sytokiinien määrää. Pidemmässä kolmen kuukauden pituisessa kliinisessä tutkimuksessa LGG:llä ei ollut vaikutusta terveiden aikuisten infektiosairastavuuteen, mutta LGG lyhensi vatsavaivojen kestoa. Probioottien vaikutukset immuunijärjestelmään näyttävät olevan kantakohtaisia ja erityisesti Lactobacillus rhamnosus GG:llä havaittiin anti-inflammatorisia vaikutuksia. Valkosolumallia ei tulisi käyttää ainoana probioottikantojen skriinausmenetelmänä niiden immunologisia vaikutuksia selvitettäessä, koska solutöiden tulokset eivät olleet täysin verrannollisia kliinisten tutkimusten tuloksiin. Sen sijaan veren perifeeristen lymfosyyttien eristäminen ja niiden aktivoitumisen selvittäminen lyhytaikaisessa kliinisessä tutkimuksessa voisi toimia suhteellisen helppona skiinausmenetelmänä.

Population structure, effective population size and adverse effects of stocking in the endangered Australian eastern freshwater cod Maccullochella ikei

Microsatellite markers were used to examine spatio-temporal genetic variation in the endangered eastern freshwater cod Maccullochella ikei in the Clarence River system, eastern Australia. High levels of population structure were detected. A model-based clustering analysis of multilocus genotypes identified four populations that were highly differentiated by F-statistics (FST = 0· 09 − 0· 49; P < 0· 05), suggesting fragmentation and restricted dispersal particularly among upstream sites. Hatchery breeding programmes were used to re-establish locally extirpated populations and to supplement remnant populations. Bayesian and frequency-based analyses of hatchery fingerling samples provided evidence for population admixture in the hatchery, with the majority of parental stock sourced from distinct upstream sites. Comparison between historical and contemporary wild-caught samples showed a significant loss of heterozygosity (21%) and allelic richness (24%) in the Mann and Nymboida Rivers since the commencement of stocking. Fragmentation may have been a causative factor; however, temporal shifts in allele frequencies suggest swamping with hatchery-produced M. ikei has contributed to the genetic decline in the largest wild population. This study demonstrates the importance of using information on genetic variation and population structure in the management of breeding and stocking programmes, particularly for threatened species.