Impacts of traditional husbandry practices on exploitable levels of genetic diversity in cultured 'Tra' catfish (Pangasianodon hypophthalmus) in the Mekong Delta, Vietnam

Sutchi catfish (Pangasianodon hypophthalmus) – known more universally by the Vietnamese name ‘Tra’ is an economically important freshwater fish in the Mekong Delta in Vietnam that constitutes an important food resource. Artificial propagation technology for Tra catfish has only recently been developed along the main branches of the Mekong River where more than 60% of the local human population participate in fishing or aquaculture. Extensive support for catfish culture in general, and that of Tra (P. hypophthalmus) in particular, has been provided by the Vietnamese government to increase both the scale of production and to develop international export markets. In 2006, total Vietnamese catfish exports reached approximately 286,602 metric tons (MT) and were valued at 736.87 $M with a number of large new export destinations being developed. Total value of production from catfish culture has been predicted to increase to approximately USD 1 billion by 2020. While freshwater catfish culture in Vietnam has a promising future, concerns have been raised about long-term quality of fry and the effectiveness of current brood stock management practices, issues that have been largely neglected to date. In this study, four DNA markers (microsatellite loci: CB4, CB7, CB12 and CB13) that were developed specifically for Tra (P. hypophthalmus) in an earlier study were applied to examine the genetic quality of artificially propagated Tra fry in the Mekong Delta in Vietnam. The goals of the study were to assess: (i) how well available levels of genetic variation in Tra brood stock used for artificial propagation in the Mekong Delta of Vietnam (breeders from three private hatcheries and Research Institute of Aquaculture No2 (RIA2) founders) has been conserved; and (ii) whether or not genetic diversity had declined significantly over time in a stock improvement program for Tra catfish at RIA2. A secondary issue addressed was how genetic markers could best be used to assist industry development. DNA was extracted from fins of catfish collected from the two main branches of the Mekong River inf Vietnam, three private hatcheries and samples from the Tra improvement program at RIA2. Study outcomes: i) Genetic diversity estimates for Tra brood stock samples were similar to, and slightly higher than, wild reference samples. In addition, the relative contribution by breeders to fry in commercial private hatcheries strongly suggest that the true Ne is likely to be significantly less than the breeder numbers used; ii) in a stock improvement program for Tra catfish at RIA2, no significant differences were detected in gene frequencies among generations (FST=0.021, P=0.036>0.002 after Bonferroni correction); and only small differences were observed in alleles frequencies among sample populations. To date, genetic markers have not been applied in the Tra catfish industry, but in the current project they were used to evaluate the levels of genetic variation in the Tra catfish selective breeding program at RIA2 and to undertake genetic correlations between genetic marker and trait variation. While no associations were detected using only four loci, they analysis provided training in the practical applications of the use of molecular markers in aquaculture in general, and in Tra culture, in particular.

Intron-mediated improvement of a selectable marker gene for plant transformation using Agrobacterium tumefaciens

In binary vectors, the antibiotic resistance gene used for selection of transformed plant cells is also usually expressed in the transforming Agrobacterium cells. This expression gives the bacterium antibiotic resistance, an unnecessary advantage on selective medium containing the antibiotic. Insertion of a castor bean catalase-1 (CAT-1) gene intron or a Parasponia andersonii haemoglobin gene intron into the coding region of the selectable marker gene, hph, completely abolished the expression of the gene in Agrobacterium, rendering it susceptible to hygromycin B. Use of these modified binary vectors minimized the overgrowth of Agrobacterium during plant transformation. Both of the introns were correctly spliced in plant cells and significantly enhanced hph gene expression in transformed rice tissue. The presence of these introns in the hph coding sequence not only maintained the selection efficiency of the hph gene, but with the CAT-1 intron also substantially increased the frequency of rice transformation. Transgenic lines with an intron-hph gene generally contained fewer gene copies and produced substantially more mRNA of the predicted size. Our results also indicate that transgenic plants with many copies of the transgene were more likely to show gene silencing than plants with 1-3 copies.

Markers of genetic susceptibility in human environmental hygiene and toxicology : the role of selected CYP, NAT and GST genes

Inherited genetic traits co-determine the susceptibility of an individual to a toxic chemical. Special emphasis has been put on individual responses to environmental and industrial carcinogens, but other chronic diseases are of increasing interest. Polymorphisms of relevant xenobiotic metabolising enzymes may be used as toxicological susceptibility markers. A growing number of genes encoding enzymes involved in biotransformation of toxicants and in cellular defence against toxicant-induced damage to the cells has been identified and cloned, leading to increased knowledge of allelic variants of genes and genetic defects that may result in a differential susceptibility toward environmental toxicants. "Low penetrating" polymorphisms in metabolism genes tend to be much more common in the population than allelic variants of "high penetrating" cancer genes, and are therefore of considerable importance from a public health point of view. Positive associations between cancer and CYP1A1 alleles, in particular the *2C I462V allele, were found for tissues following the aerodigestive tract. Again, in most cases, the effect of the variant CYP1A1 allele becomes apparent or clearer in connection with the GSTM1 null allele. The CYP1B1 codon 432 polymorphism (CYP1B1*3) has been identified as a susceptibility factor in smoking-related head-and-neck squameous cell cancer. The impact of this polymorphic variant of CYP1B1 on cancer risk was also reflected by an association with the frequency of somatic mutations of the p53 gene. Combined genotype analysis of CYP1B1 and the glutathione transferases GSTM1 or GSTT1 has also pointed to interactive effects. Of particular interest for the industrial and environmental field is the isozyme CYP2E1. Several genotypes of this isozyme have been characterised which seem to be associated with different levels of expression of enzyme activity. The acetylator status for NAT2 can be determined by genotyping or by phenotyping. In the pathogenesis of human bladder cancer due to occupational exposure to "classical" aromatic amines (benzidine, 4-aminodiphenyl, 1-naphthylamine) acetylation by NAT2 is regarded as a detoxication step. Interestingly, the underlying European findings of a higher susceptibility of slow acetylators towards aromatic amines are in contrast to findings in Chinese workers occupationally exposed to aromatic amines which points to different mechanisms of susceptibility between European and Chinese populations. Regarding human bladder cancer, the hypothesis has been put forward that genetic polymorphism of GSTM1 might be linked with the occurrence of this tumour type. This supports the hypothesis that exposure to PAH might causally be involved in urothelial cancers. The human polymorphic GST catalysing conjugation of halomethanes, dihalomethanes, ethylene oxide and a number of other industrial compounds could be characterised as a class theta enzyme (GSTT1) by means of molecular biology. "Conjugator" and "non-conjugator" phenotypes are coincident with the presence and absence of the GSTT1 gene. There are wide variations in the frequencies of GSTT1 deletion (GSTT1 *0/0) among different ethnicities. Human phenotyping is facilitated by the GST activity towards methyl bromide or ethylene oxide in erythrocytes which is representative of the metabolic GSTT1 competence of the entire organism. Inter-individual variations in xenobiotic metabolism capacities may be due to polymorphisms of the genes coding for the enzymes themselves or of the genes coding for the receptors or transcription factors which regulate the expression of the enzymes. Also, polymorphisms in several regions of genes may cause altered ligand affinity, transactivation activity or expression levels of the receptor subsequently influencing the expression of the downstream target genes. Studies of individual susceptibility to toxicants and gene-environment interaction are now emerging as an important component of molecular epidemiology.

The future for genetic studies in reproduction

Genetic factors contribute to risk of many common diseases affecting reproduction and fertility. In recent years, methods for genome-wide association studies(GWAS) have revolutionized gene discovery forcommontraits and diseases. Results of GWAS are documented in the Catalog of Published Genome-Wide Association Studies at the National Human Genome Research Institute and report over 70 publications for 32 traits and diseases associated with reproduction. These include endometriosis, uterine fibroids, age at menarche and age at menopause. Results that pass appropriate stringent levels of significance are generally well replicated in independent studies. Examples of genetic variation affecting twinning rate, infertility, endometriosis and age at menarche demonstrate that the spectrum of disease-related variants for reproductive traits is similar to most other common diseases.GWAS 'hits' provide novel insights into biological pathways and the translational value of these studies lies in discovery of novel gene targets for biomarkers, drug development and greater understanding of environmental factors contributing to disease risk. Results also show that genetic data can help define sub-types of disease and co-morbidity with other traits and diseases. To date, many studies on reproductive traits have used relatively small samples. Future genetic marker studies in large samples with detailed phenotypic and clinical information will yield new insights into disease risk, disease classification and co-morbidity for many diseases associated with reproduction and infertility.

Association of variants at 1q32 and STAT3 with ankylosing spondylitis suggests genetic overlap with Crohn's disease

Ankylosing spondylitis (AS) is a common inflammatory arthritic condition. Overt inflammatory bowel disease (IBD) occurs in about 10% of AS patients, and in addition 70% of AS cases may have subclinical terminal ileitis. Spondyloarthritis is also common in IBD patients. We therefore tested Crohn's disease susceptibility genes for association with AS, aiming to identify pleiotropic genetic associations with both diseases. Genotyping was carried out using Sequenom and Applied Biosystems TaqMan and OpenArray technologies on 53 markers selected from 30 Crohn's disease associated genomic regions. We tested genotypes in a population of unrelated individual cases (n = 2,773) and controls (n = 2,215) of white European ancestry for association with AS. Statistical analysis was carried out using a Cochran-Armitage test for trend in PLINK. Strong association was detected at chr1q32 near KIF21B (rs11584383, P = 1.66 x 10-10, odds ratio (OR) = 0.74, 95% CI:0.68-0.82). Association with disease was also detected for 2 variants within STAT3 (rs6503695, P = 4.6×10-4. OR = 0.86 (95% CI:0.79-0.93); rs744166, P = 2.6×10-5, OR = 0.84 (95% CI:0.77-0.91)). Association was confirmed for IL23R (rs11465804, P = 1.2×10-5, OR = 0.65 (95% CI:0.54-0.79)), and further associations were detected for IL12B (rs10045431, P = 5.261025, OR = 0.83 (95% CI:0.76-0.91)), CDKAL1 (rs6908425, P = 1.1×10-4, OR = 0.82 (95% CI:0.74-0.91)), LRRK2/MUC19 (rs11175593, P = 9.9×10-5, OR = 1.92 (95% CI: 1.38-2.67)), and chr13q14 (rs3764147, P = 5.9×10-4, OR = 1.19 (95% CI: 1.08-1.31)). Excluding cases with clinical IBD did not significantly affect these findings. This study identifies chr1q32 and STAT3 as ankylosing spondylitis susceptibility loci. It also further confirms association for IL23R and detects suggestive association with another 4 loci. STAT3 is a key signaling molecule within the Th17 lymphocyte differentiation pathway and further enhances the case for a major role of this T-lymphocyte subset in ankylosing spondylitis. Finally these findings suggest common aetiopathogenic pathways for AS and Crohn's disease and further highlight the involvement of common risk variants across multiple diseases.

Development of type I genetic markers from expressed sequence tags in highly polymorphic species

Expressed sequence tag (EST) databases provide a primary source of nuclear DNA sequences for genetic marker development in non-model organisms. To date, the process has been relatively inefficient for several reasons: - 1) priming site polymorphism in the template leads to inferior or erratic amplification; - 2) introns in the target amplicon are too large and/or numerous to allow effective amplification under standard screening conditions, and; - 3) at least occasionally, a PCR primer straddles an exon–intron junction and is unable to bind to genomic DNA template. The first is only a minor issue for species or strains with low heterozygosity but becomes a significant problem for species with high genomic variation, such as marine organisms with extremely large effective population sizes. Problems arising from unanticipated introns are unavoidable but are most pronounced in intron-rich species, such as vertebrates and lophotrochozoans. We present an approach to marker development in the Pacific oyster Crassostrea gigas, a highly polymorphic and intron-rich species, which minimizes these problems, and should be applicable to other non-model species for which EST databases are available. Placement of PCR primers in the 3′ end of coding sequence and 3′ UTR improved PCR success rate from 51% to 97%. Almost all (37 of 39) markers developed for the Pacific oyster were polymorphic in a small test panel of wild and domesticated oysters.

The influence of habitat heterogeneity on patterns of connectivity among rabbit populations in southern Queensland

Patterns of connectivity among local populations influence the dynamics of regional systems, but most ecological models have concentrated on explaining the effect of connectivity on local population structure using dynamic processes covering short spatial and temporal scales. In this study, a model was developed in an extended spatial system to examine the hypothesis that long term connectivity levels among local populations are influenced by the spatial distribution of resources and other habitat factors. The habitat heterogeneity model was applied to local wild rabbit populations in the semi-arid Mitchell region of southern central Queensland (the Eastern system). Species' specific population parameters which were appropriate for the rabbit in this region were used. The model predicted a wide range of long term connectivity levels among sites, ranging from the extreme isolation of some sites to relatively high interaction probabilities for others. The validity of model assumptions was assessed by regressing model output against independent population genetic data, and explained over 80% of the variation in the highly structured genetic data set. Furthermore, the model was robust, explaining a significant proportion of the variation in the genetic data over a wide range of parameters. The performance of the habitat heterogeneity model was further assessed by simulating the widely reported recent range expansion of the wild rabbit into the Mitchell region from the adjacent, panmictic Western rabbit population system. The model explained well the independently determined genetic characteristics of the Eastern system at different hierarchic levels, from site specific differences (for example, fixation of a single allele in the population at one site), to differences between population systems (absence of an allele in the Eastern system which is present in all Western system sites). The model therefore explained the past and long term processes which have led to the formation and maintenance of the highly structured Eastern rabbit population system. Most animals exhibit sex biased dispersal which may influence long term connectivity levels among local populations, and thus the dynamics of regional systems. When appropriate sex specific dispersal characteristics were used, the habitat heterogeneity model predicted substantially different interaction patterns between female-only and combined male and female dispersal scenarios. In the latter case, model output was validated using data from a bi-parentally inherited genetic marker. Again, the model explained over 80% of the variation in the genetic data. The fact that such a large proportion of variability is explained in two genetic data sets provides very good evidence that habitat heterogeneity influences long term connectivity levels among local rabbit populations in the Mitchell region for both males and females. The habitat heterogeneity model thus provides a powerful approach for understanding the large scale processes that shape regional population systems in general. Therefore the model has the potential to be useful as a tool to aid in the management of those systems, whether it be for pest management or conservation purposes.

Inheritance of Resistance to Rice Tungro Spherical Virus in a Near-Isogenic Line Derived from Utri Merah and in Rice Cultivar TKM6

Resistance to rice virus diseases is an important requirement in many Southeast Asian rice breeding programs. Inheritance of resistance to rice tungro spherical virus (RTSV) in TW5, a near-isogenic line derived from Indonesian rice cultivar Utri Merah, was compared to that in TKM6, an Indian rice cultivar. Both TKM6 and Utri Merah are cultivars resistant to RTSV infections. Crosses were made between TKM6 and TN1, a susceptible cultivar, and between TW5 and TN1, and F3 lines were evaluated for their resistance to RTSV using two RTSV inoculum sources and a serological assay (ELISA). In TKM6, the resistance to the mixture of RTSV-V + RTBV inoculum source was controlled by a single recessive gene, whereas in TW5, the resistance was controlled by two recessive genes. A single recessive gene, however, controlled the resistance in TW5 when another RTSV variant, RTSV-VI, was used, suggesting that the resistance in TW5 depends on the nature of the RTSV inoculum used. RT-PCR, sequence, and phylogenetic analyses confirmed that RTSV-VI inoculum differs from RTSV-V inoculum and accurate phenotyping of the resistance to RTSV requires the use of a genetic marker.

Prediction of IBD based on population history for fine gene mapping

A novel multiple regression method (RM) is developed to predict identity-by-descent probabilities at a locus L (IBDL), among individuals without pedigree, given information on surrounding markers and population history. These IBDL probabilities are a function of the increase in linkage disequilibrium (LD) generated by drift in a homogeneous population over generations. Three parameters are sufficient to describe population history: effective population size (Ne), number of generations since foundation (T), and marker allele frequencies among founders (p). IBD L are used in a simulation study to map a quantitative trait locus (QTL) via variance component estimation. RM is compared to a coalescent method (CM) in terms of power and robustness of QTL detection. Differences between RM and CM are small but significant. For example, RM is more powerful than CM in dioecious populations, but not in monoecious populations. Moreover, RM is more robust than CM when marker phases are unknown or when there is complete LD among founders or Ne is wrong, and less robust when p is wrong. CM utilises all marker haplotype information, whereas RM utilises information contained in each individual marker and all possible marker pairs but not in higher order interactions. RM consists of a family of models encompassing four different population structures, and two ways of using marker information, which contrasts with the single model that must cater for all possible evolutionary scenarios in CM.

Power of QTL detection using association tests with family controls

The power of testing for a population-wide association between a biallelic quantitative trait locus and a linked biallelic marker locus is predicted both empirically and deterministically for several tests. The tests were based on the analysis of variance (ANOVA) and on a number of transmission disequilibrium tests (TDT). Deterministic power predictions made use of family information, and were functions of population parameters including linkage disequilibrium, allele frequencies, and recombination rate. Deterministic power predictions were very close to the empirical power from simulations in all scenarios considered in this study. The different TDTs had very similar power, intermediate between one-way and nested ANOVAs. One-way ANOVA was the only test that was not robust against spurious disequilibrium. Our general framework for predicting power deterministically can be used to predict power in other association tests. Deterministic power calculations are a powerful tool for researchers to plan and evaluate experiments and obviate the need for elaborate simulation studies.

Candidate gene analysis for quantitative traits using the transmission disequilibrium test : the example of the melanocortin 4-receptor in pigs

Population-wide associations between loci due to linkage disequilibrium can be used to map quantitative trait loci (QTL) with high resolution. However, spurious associations between markers and QTL can also arise as a consequence of population stratification. Statistical methods that cannot differentiate between loci associations due to linkage disequilibria from those caused in other ways can render false-positive results. The transmission-disequilibrium test (TDT) is a robust test for detecting QTL. The TDT exploits within-family associations that are not affected by population stratification. However, some TDTs are formulated in a rigid-form, with reduced potential applications. In this study we generalize TDT using mixed linear models to allow greater statistical flexibility. Allelic effects are estimated with two independent parameters: one exploiting the robust within-family information and the other the potentially biased between-family information. A significant difference between these two parameters can be used as evidence for spurious association. This methodology was then used to test the effects of the fourth melanocortin receptor (MC4R) on production traits in the pig. The new analyses supported the previously reported results; i.e., the studied polymorphism is either causal of in very strong linkage disequilibrium with the causal mutation, and provided no evidence for spurious association.

Agrobacterium-mediated transformation of Australian rice cultivars Jarrah and Amaroo using modified promoters and selectable markers

We report the first successful Agrobacterium-mediated transformation of Australian elite rice cultivars, Jarrah and Amaroo, using binary vectors with our improved promoters and selectable markers. Calli derived from mature embryos were used as target tissues. The binary vectors contained hph (encoding hygromycin resistance) or bar (encoding herbicide resistance) as the selectable marker gene and uidA (gus) or sgfpS65T as the reporter gene driven by different promoters. Use of Agrobacterium strain AGL1 carrying derivatives of an improved binary vector pWBVec8, wherein the CaMV35S driven hph gene is interrupted by the castor bean catalase 1 intron, produced a 4-fold higher number of independent transgenic lines compared to that produced with the use of strain EHA101 carrying the binary vector pIG121-Hm wherein the CaMV35S driven hph is intronless. The Ubiquitin promoter produced 30-fold higher β-glucuronidase (GUS) activity (derivatives of binary vector pWBVec8) in transgenic plants than the CaMV35S promoter (pIG121-Hm). The two modified SCSV promoters produced GUS activity comparable to that produced by the Ubiquitin promoter. Progeny analysis (R1) for hygromycin resistance and GUS activity with selected lines showed both Mendelian and non-Mendelian segregation. Lines showing very high levels of GUS activity in T0 showed a reduced level of GUS activity in their T1 progeny, while lines with moderate levels of GUS activity showed increased levels in T1 progeny. Stable heritable green fluorescent protein (GFP) expression was also observed in few transgenic plants produced with the binary vector pTO134 which had the CaMV35S promoter-driven selectable marker gene bar and a modified CaMV35S promoter-driven reporter gene sgfpS65T.

A multi-genome analysis approach enables tracking of the invasion of a single Russian wheat aphid (Diuraphis noxia) clone throughout the New World

This study investigated the population genetics, demographic history and pathway of invasion of the Russian wheat aphid (RWA) from its native range in Central Asia, the Middle East and Europe to South Africa and the Americas. We screened microsatellite markers, mitochondrial DNA and endosymbiont genes in 504 RWA clones from nineteen populations worldwide. Following pathway analyses of microsatellite and endosymbiont data, we postulate that Turkey and Syria were the most likely sources of invasion to Kenya and South Africa, respectively. Furthermore, we found that one clone transferred between South Africa and the Americas was most likely responsible for the New World invasion. Finally, endosymbiont DNA was found to be a high resolution population genetic marker, extremely useful for studies of invasion over a relatively short evolutionary history time frame. This study has provided valuable insights into the factors that may have facilitated the recent global invasion by this damaging pest.

Genome-wide SNP detection, validation, and development of an 8K SNP array for apple

As high-throughput genetic marker screening systems are essential for a range of genetics studies and plant breeding applications, the International RosBREED SNP Consortium (IRSC) has utilized the Illumina Infinium® II system to develop a medium- to high-throughput SNP screening tool for genome-wide evaluation of allelic variation in apple (Malus×domestica) breeding germplasm. For genome-wide SNP discovery, 27 apple cultivars were chosen to represent worldwide breeding germplasm and re-sequenced at low coverage with the Illumina Genome Analyzer II. Following alignment of these sequences to the whole genome sequence of 'Golden Delicious', SNPs were identified using SoapSNP. A total of 2,113,120 SNPs were detected, corresponding to one SNP to every 288 bp of the genome. The Illumina GoldenGate® assay was then used to validate a subset of 144 SNPs with a range of characteristics, using a set of 160 apple accessions. This validation assay enabled fine-tuning of the final subset of SNPs for the Illumina Infinium® II system. The set of stringent filtering criteria developed allowed choice of a set of SNPs that not only exhibited an even distribution across the apple genome and a range of minor allele frequencies to ensure utility across germplasm, but also were located in putative exonic regions to maximize genotyping success rate. A total of 7867 apple SNPs was established for the IRSC apple 8K SNP array v1, of which 5554 were polymorphic after evaluation in segregating families and a germplasm collection. This publicly available genomics resource will provide an unprecedented resolution of SNP haplotypes, which will enable marker-locus-trait association discovery, description of the genetic architecture of quantitative traits, investigation of genetic variation (neutral and functional), and genomic selection in apple.

The interleukin 1 gene cluster contains a major susceptibility locus for ankylosing spondylitis

Ankylosing spondylitis (AS) is a common and highly heritable inflammatory arthropathy. Although the gene HLA-B27 is almost essential for the inheritance of the condition, it alone is not sufficient to explain the pattern of familial recurrence of the disease. We have previously demonstrated suggestive linkage of AS to chromosome 2q13, a region containing the interleukin 1 (IL-1) family gene cluster, which includes several strong candidates for involvement in the disease. In the current study, we describe strong association and transmission of IL-1 family gene cluster single-nucleotide polymorphisms and haplotypes with AS.