Genetic Structure of Cochliobolus sativus Populations Sampled from Roots and Leaves of Barley and Wheat in North Dakota

Common root rot (CRR) and spot blotch, caused by Cochliobolus sativus (Ito and Kurib.) Drechsl. ex Dast., are important diseases of barley (Hordeum vulgare L.) and wheat (Triticum aestivum L.) worldwide. However, the population biology of C. sativus is still poorly understood. In this study, the genetic structure of three C. sativus populations, consisting of isolates sampled respectively from barley leaves (BL), barley roots (BR) and wheat roots (WR) in North Dakota, was analysed with amplified fragment length polymorphism (AFLP) markers. A total of 127 AFLP loci were generated among 208 C. sativus isolates analysed with three primer combinations. Gene diversity (H = 0.277-0.335) were high in all three populations. Genetic variation among C. sativus individuals within population accounted for 74%, whereas 26% of the genetic variation was explained among populations. Genetic differentiation was high (empty set PT = 0.261, corrected G ''(st)= 0.39), whereas gene flow (Nm) ranged from 1.27 to 1.56 among the three populations analysed. The multilocus linkage disequilibrium (LD) ((r) over bard = 0.0760.117) was moderate in C. sativus populations. Cluster analyses indicate that C. sativus populations differentiated according to the hosts (barley and wheat) and tissues (root and leaf) although generalists also exist in North Dakota. Crop breeding may benefit from combining genes for resistance against both specialists and generalists of C. sativus.

The Genetic Architecture of Climatic Adaptation of Tropical Cattle

Adaptation of global food systems to climate change is essential to feed the world. Tropical cattle production, a mainstay of profitability for farmers in the developing world, is dominated by heat, lack of water, poor quality feedstuffs, parasites, and tropical diseases. In these systems European cattle suffer significant stock loss, and the cross breeding of taurine x indicine cattle is unpredictable due to the dilution of adaptation to heat and tropical diseases. We explored the genetic architecture of ten traits of tropical cattle production using genome wide association studies of 4,662 animals varying from 0% to 100% indicine. We show that nine of the ten have genetic architectures that include genes of major effect, and in one case, a single location that accounted for more than 71% of the genetic variation. One genetic region in particular had effects on parasite resistance, yearling weight, body condition score, coat colour and penile sheath score. This region, extending 20 Mb on BTA5, appeared to be under genetic selection possibly through maintenance of haplotypes by breeders. We found that the amount of genetic variation and the genetic correlations between traits did not depend upon the degree of indicine content in the animals. Climate change is expected to expand some conditions of the tropics to more temperate environments, which may impact negatively on global livestock health and production. Our results point to several important genes that have large effects on adaptation that could be introduced into more temperate cattle without detrimental effects on productivity.

Causes and consequences of variation in nestling immune function

Defence against pathogens is a vital need of all living organisms that has led to the evolution of complex immune mechanisms. However, although immunocompetence the ability to resist pathogens and control infection has in recent decades become a focus for research in evolutionary ecology, the variation in immune function observed in natural populations is relatively little understood. This thesis examines sources of this variation (environmental, genetic and maternal effects) during the nestling stage and its fitness consequences in wild populations of passerines: the blue tit (Cyanistes caeruleus) and the collared flycatcher (Ficedula albicollis). A developing organism may face a dilemma as to whether to allocate limited resources to growth or to immune defences. The optimal level of investment in immunity is shaped inherently by specific requirements of the environment. If the probability of contracting infection is low, maintaining high growth rates even at the expense of immune function may be advantageous for nestlings, as body mass is usually a good predictor of post-fledging survival. In experiments with blue tits and haematophagous hen fleas (Ceratophyllus gallinae) using two methods, methionine supplementation (to manipulate nestlings resource allocation to cellular immune function) and food supplementation (to increase resource availability), I confirmed that there is a trade-off between growth and immunity and that the abundance of ectoparasites is an environmental factor affecting allocation of resources to immune function. A cross-fostering experiment also revealed that environmental heterogeneity in terms of abundance of ectoparasites may contribute to maintaining additive genetic variation in immunity and other traits. Animal model analysis of extensive data collected from the population of collared flycatchers on Gotland (Sweden) allowed examination of the narrow-sense heritability of PHA-response the most commonly used index of cellular immunocompetence in avian studies. PHA-response is not heritable in this population, but is subject to a non-heritable origin (presumably maternal) effect. However, experimental manipulation of yolk androgen levels indicates that the mechanism of the maternal effect in PHA-response is not in ovo deposition of androgens. The relationship between PHA-response and recruitment was studied for over 1300 collared flycatcher nestlings. Multivariate selection analysis shows that it is body mass, not PHA-response, that is under direct selection. PHA-response appears to be related to recruitment because of its positive relationship with body mass. These results imply that either PHA-response fails to capture the immune mechanisms that are relevant for defence against pathogens encountered by fledglings or that the selection pressure from parasites is not as strong as commonly assumed.

Applications of gene sequence polymorphisms in evolutionary genetic studies of Atlantic salmon (Salmo salar) and other teleost fish species

Evolutionary genetics incorporates traditional population genetics and studies of the origins of genetic variation by mutation and recombination, and the molecular evolution of genomes. Among the primary forces that have potential to affect the genetic variation within and among populations, including those that may lead to adaptation and speciation, are genetic drift, gene flow, mutations and natural selection. The main challenges in knowing the genetic basis of evolutionary changes is to distinguish the adaptive selection forces that cause existent DNA sequence variants and also to identify the nucleotide differences responsible for the observed phenotypic variation. To understand the effects of various forces, interpretation of gene sequence variation has been the principal basis of many evolutionary genetic studies. The main aim of this thesis was to assess different forms of teleost gene sequence polymorphisms in evolutionary genetic studies of Atlantic salmon (Salmo salar) and other species. Firstly, the level of Darwinian adaptive evolution affected coding regions of the growth hormone (GH) gene during the teleost evolution was investigated based on the sequence data existing in public databases. Secondly, a target gene approach was used to identify within population variation in the growth hormone 1 (GH1) gene in salmon. Then, a new strategy for single nucleotide polymorphisms (SNPs) discovery in salmonid fishes was introduced, and, finally, the usefulness of a limited number of SNP markers as molecular tools in several applications of population genetics in Atlantic salmon was assessed. This thesis showed that the gene sequences in databases can be utilized to perform comparative studies of molecular evolution, and some putative evidence of the existence of Darwinian selection during the teleost GH evolution was presented. In addition, existent sequence data was exploited to investigate GH1 gene variation within Atlantic salmon populations throughout its range. Purifying selection is suggested to be the predominant evolutionary force controlling the genetic variation of this gene in salmon, and some support for gene flow between continents was also observed. The novel approach to SNP discovery in species with duplicated genome fragments introduced here proved to be an effective method, and this may have several applications in evolutionary genetics with different species - e.g. when developing gene-targeted markers to investigate quantitative genetic variation. The thesis also demonstrated that only a few SNPs performed highly similar signals in some of the population genetic analyses when compared with the microsatellite markers. This may have useful applications when estimating genetic diversity in genes having a potential role in ecological and conservation issues, or when using hard biological samples in genetic studies as SNPs can be applied with relatively highly degraded DNA.

Associations of interleukin-1 (IL-1) gene variation and IL-1 receptor antagonist phenotypes with immunological responses, metabolic dysregulation and type 2 diabetes

The upstream proinflammatory interleukin-1 (IL-1) cytokines, together with a naturally occurring IL-1 receptor antagonist (IL-1Ra), play a significant role in several diseases and physiologic conditions. The IL-1 proteins affect glucose homeostasis at multiple levels contributing to vascular injuries and metabolic dysregulations that precede diabetes. An association between IL-1 gene variations and IL-1Ra levels has been suggested, and genetic studies have reported associations with metabolic dysregulation and altered inflammatory responses. The principal aims of this study were to: 1) examine the associations of IL-1 gene variation and IL-1Ra expression in the development and persistence of thyroid antibodies in subacute thyroiditis; 2) investigate the associations of common variants in the IL-1 gene family with plasma glucose and insulin concentrations, glucose homeostasis measures and prevalent diabetes in a representative population sample; 3) investigate genetic and non-genetic determinants of IL-1Ra phenotypes in a cross-sectional setting in three independent study populations; 4) investigate in a prospective setting (a) whether variants of the IL-1 gene family are predictors for clinically incident diabetes in two population-based observational cohort studies; and (b) whether the IL-1Ra levels predict the progression of metabolic syndrome to overt diabetes during the median follow-up of 10.8 and 7.1 years. Results from on patients with subacte thyroiditis showed that the systemic IL-1Ra levels are elevated during a specific proinflammatory response and they correlated with C-reactive protein (CRP) levels. Genetic variation in the IL-1 family seemed to have an association with the appearance of thyroid peroxidase antibodies and persisting local autoimmune responses during the follow-up. Analysis of patients suffering from diabetes and metabolic traits suggested that genetic IL-1 variation and IL-1Ra play a role in glucose homeostasis and in the development of type 2 diabetes. The coding IL-1 beta SNP rs1143634 was associated with traits related to insulin resistance in cross-sectional analyses. Two haplotype variants of the IL-1 beta gene were associated with prevalent diabetes or incident diabetes in a prospective setting and both of these haplotypes were tagged by rs1143634. Three variants of the IL-1Ra gene and one of the IL-1 beta gene were consistently identified as significant, independent determinants of the IL-1Ra phenotype in two or three populations. The proportion of the phenotypic variation explained by the genetic factors was modest however, while obesity and other metabolic traits explained a larger part. Body mass index was the strongest predictor of systemic IL-1Ra concentration overall. Furthermore, the age-adjusted IL-1Ra concentrations were elevated in individuals with metabolic syndrome or diabetes when compared to those free of metabolic dysregulation. In prospective analyses the systemic IL-1Ra levels were found as independent predictors for the development of diabetes in people with metabolic syndrome even after adjustment for multiple other factors, including plasma glucose and CRP levels. The predictive power of IL-1Ra was better than that of CRP. The prospective results also provided some evidence for a role of common IL-1 alpha promoter SNP rs1800587 in the development of type 2 diabetes among men and suggested that the role may be gender specific. Likewise, common variations in the IL-1 beta coding region may have a gender specific association with diabetes development. Further research on the potential benefits of IL-1Ra measurements in identifying individuals at high risk for diabetes, who then could be targeted for specific treatment interventions, is warranted. It has been reported in the recent literature that IL-1Ra secreted from adipose tissue has beneficial effects on glucose homeostasis. Furthermore, treatment with recombinant human IL-1Ra has been shown to have a substantial therapeutic potential. The genetic results from the prospective analyses performed in this study remain inconclusive, but together with the cross-sectional analyses they suggest gender-specific effects of the IL-1 variants on the risk of diabetes. Larger studies with more extensive genotyping and resequencing may help to pinpoint the exact variants responsible and to further elucidate the biological mechanisms for the observed associations. This would improve our understanding of the pathways linking inflammation and obesity with glucose and insulin metabolism.

Association between molecular variation in the phosphoglucose isomerase (Pgi) locus and migration propensity in the Glanville fritillary butterfly

Habitat fragmentation produces patches of suitable habitat surrounded by unfavourable matrix habitat. A species may persist in such a fragmented landscape in an equilibrium between the extinctions and recolonizations of local populations, thus forming a metapopulation. Migration between local populations is necessary for the long-term persistence of a metapopulation. The Glanville fritillary butterfly (Melitaea cinxia) forms a metapopulation in the Åland islands in Finland. There is migration between the populations, the extent of which is affected by several environmental factors and variation in the phenotype of individual butterflies. Different allelic forms of the glycolytic enzyme phosphoglucose isomerase (Pgi) has been identified as a possible genetic factor influencing flight performance and migration rate in this species. The frequency of a certain Pgi allele, Pgi-f, follows the same pattern in relation to population age and connectivity as migration propensity. Furthermore, variation in flight metabolic performance, which is likely to affect migration propensity, has been linked to genetic variation in Pgi or a closely linked locus. The aim of this study was to investigate the association between Pgi genotype and the migration propensity in the Glanville fritillary both at the individual and population levels using a statistical modelling approach. A mark-release-recapture (MRR) study was conducted in a habitat patch network of M. cinxia in Åland to collect data on the movements of individual butterflies. Larval samples from the study area were also collected for population level examinations. Each butterfly and larva was genotyped at the Pgi locus. The MRR data was parameterised with two mathematical models of migration: the Virtual Migration Model (VM) and the spatially explicit diffusion model. VM model predicted and observed numbers of emigrants from populations with high and low frequencies of Pgi-f were compared. Posterior predictive data sets were simulated based on the parameters of the diffusion model. Lack-of-fit of observed values to the model predicted values of several descriptors of movements were detected, and the effect of Pgi genotype on the deviations was assessed by randomizations including the genotype information. This study revealed a possible difference in the effect of Pgi genotype on migration propensity between the two sexes in the Glanville fritillary. The females with and males without the Pgi-f allele moved more between habitat patches, which is probably related to differences in the function of flight in the two sexes. Females may use their high flight capacity to migrate between habitat patches to find suitable oviposition sites, whereas males may use it to acquire mates by keeping a territory and fighting off other intruding males, possibly causing them to emigrate. The results were consistent across different movement descriptors and at the individual and population levels. The effect of Pgi is likely to be dependent on the structure of the landscape and the prevailing environmental conditions.

Variation at Diabetes- and Obesity-Associated Loci May Mirror Neutral Patterns of Human Population Diversity and Diabetes Prevalence in India

South Asian populations harbor a high degree of genetic diversity, due in part to demographic history. Two studies on genome-wide variation in Indian populations have shown that most Indian populations show varying degrees of admixture between ancestral north Indian and ancestral south Indian components. As a result of this structure, genetic variation in India appears to follow a geographic cline. Similarly, Indian populations seem to show detectable differences in diabetes and obesity prevalence between different geographic regions of the country. We tested the hypothesis that genetic variation at diabetes-and obesity-associated loci may be potentially related to different genetic ancestries. We genotyped 2977 individuals from 61 populations across India for 18 SNPs in genes implicated in T2D and obesity. We examined patterns of variation in allele frequency across different geographical gradients and considered state of origin and language affiliation. Our results show that most of the 18 SNPs show no significant correlation with latitude, the geographic cline reported in previous studies, or by language family. Exceptions include KCNQ1 with latitude and THADA and JAK1 with language, which suggests that genetic variation at previously ascertained diabetes-associated loci may only partly mirror geographic patterns of genome-wide diversity in Indian populations.

Differences in host species relationships and biogeographic influences produce contrasting patterns of prevalence, community composition and genetic structure in two genera of avian malaria parasites in southern Melanesia

1. Host-parasite interactions have the potential to influence broadscale ecological and evolutionary processes, levels of endemism, divergence patterns and distributions in host populations. Understanding the mechanisms involved requires identification of the factors that shape parasite distribution and prevalence. 2. A lack of comparative information on community-level host-parasite associations limits our understanding of the role of parasites in host population divergence processes. Avian malaria (haemosporidian) parasites in bird communities offer a tractable model system to examine the potential for pathogens to influence evolutionary processes in natural host populations. 3. Using cytochrome b variation, we characterized phylogenetic diversity and prevalence of two genera of avian haemosporidian parasites, Plasmodium and Haemoproteus, and analysed biogeographic patterns of lineages across islands and avian hosts, in southern Melanesian bird communities to identify factors that explain patterns of infection. 4. Plasmodium spp. displayed isolation-by-distance effects, a significant amount of genetic variation distributed among islands but insignificant amounts among host species and families, and strong local island effects with respect to prevalence. Haemoproteus spp. did not display isolation-by-distance patterns, showed marked structuring of genetic variation among avian host species and families, and significant host species prevalence patterns. 5. These differences suggest that Plasmodium spp. infection patterns were shaped by geography and the abiotic environment, whereas Haemoproteus spp. infection patterns were shaped predominantly by host associations. Heterogeneity in the complement and prevalence of parasite lineages infecting local bird communities likely exposes host species to a mosaic of spatially divergent disease selection pressures across their naturally fragmented distributions in southern Melanesia. Host associations for Haemoproteus spp. indicate a capacity for the formation of locally co-adapted host-parasite relationships, a feature that may limit intraspecific gene flow or range expansions of closely related host species.

A Comparison of the SNP Variation in the Calcification PMCA Gene of Lottia gigantea between a Santa Barbara Population and a More Acidic Monterey Bay Population

As the atmospheric levels of CO2 rise from human activity, the carbonic acid levels of the ocean increase, causing ocean acidification. This increase in acidity breaks down the calcified bodies that many marine organisms depend upon. Upwelling regions such as Monterey Bay in California have pH levels that are not expected to reach the open ocean for a few decades. This study reviews one of the common intertidal animals of the California coast, the Owl Limpet Lottia gigantea, and its genetic variation of the plasma membrane Ca2+ ATPase (PMCA) in relation to the acidity of its environment. The PMCA protein functions in the calcification process of many organisms. Specifically in limpets, this gene functions to form its protective shell. Single-nucleotide polymorphisms (SNPs) were found among five sections of the gene to determine variation between the acidic environment population in Monterey, California and the non-acidic environment population in Santa Barbara, California. While some variation was determined, the Monterey Bay and Santa Barbara Lottia gigantea populations are not significantly distinct at the PMCA gene. Sections B, C, and D were found to be linked. Only one location in Section B was found to have an amino acid change within an exon. Section A has the strongest connection to the sampling location. Monterey individuals were seen to be more genetically recognizable, while Santa Barbara individuals showed slightly more variation. Understanding the trends of ocean acidification, upwelling region activities, and population genetics will assist in determining how the ocean environment will behave in the future.

A molecular genetic investigation of the population structure of Atlantic menhaden (Brevoortia tyrannus)

Atlantic menhaden (Brevoortia tyrannus), through landings, support one of the largest commercial fisheries in the United States. Recent consolidation of the once coast-wide reduction fishery to waters within and around Chesapeake Bay has raised concerns over the possibility of the loss of unique genetic variation resulting from concentrated fishing pressure. To address this question, we surveyed variation at the mitochondrial cytochrome c oxidase subunit I (COI) gene region and seven nuclear microsatellite loci to evaluate stock structure of Atlantic menhaden. Samples were collected from up to three cohorts of Atlantic menhaden at four geographic locations along the U.S. Atlantic coast in 2006 and 2007, and from the closely related Gulf menhaden (B. patronus) in the Gulf of Mexico. Genetic divergence between Atlantic menhaden and Gulf menhaden, based on the COI gene region sequences and microsatellite loci, was more characteristic of conspecific populations than separate species. Hierarchical analyses of molecular variance indicated a homogeneous distribution of genetic variation within Atlantic menhaden. No significant variation was found between young-of-the-year menhaden (YOY) collected early and late in the season within Chesapeake Bay, between young-of-the-year and yearling menhaden collected in the Chesapeake Bay during the same year, between YOY and yearling menhaden taken in Chesapeake Bay in successive years, or among combined YOY and yearling Atlantic menhaden collected in both years from the four geographic locations. The genetic connectivity between the regional collections indicates that the concentration of fishing pressure in and around Chesapeake Bay will not result in a significant loss of unique genetic variation.

Exploring Genetic Factors Involved in Huntington Disease Age of Onset: E2F2 as a New Potential Modifier Gene

Age of onset (AO) of Huntington disease (HD) is mainly determined by the length of the CAG repeat expansion (CAGexp) in exon 1 of the HTT gene. Additional genetic variation has been suggested to contribute to AO, although the mechanism by which it could affect AO is presently unknown. The aim of this study is to explore the contribution of candidate genetic factors to HD AO in order to gain insight into the pathogenic mechanisms underlying this disorder. For that purpose, two AO definitions were used: the earliest age with unequivocal signs of HD (earliest AO or eAO), and the first motor symptoms age (motor AO or mAO). Multiple linear regression analyses were performed between genetic variation within 20 candidate genes and eAO or mAO, using DNA and clinical information of 253 HD patients from REGISTRY project. Gene expression analyses were carried out by RT-qPCR with an independent sample of 35 HD patients from Basque Country Hospitals. We found suggestive association signals between HD eAO and/or mAO and genetic variation within the E2F2, ATF7IP, GRIN2A, GRIN2B, LINC01559, HIP1 and GRIK2 genes. Among them, the most significant was the association between eAO and rs2742976, mapping to the promoter region of E2F2 transcription factor. Furthermore, rs2742976 T allele patient carriers exhibited significantly lower lymphocyte E2F2 gene expression, suggesting a possible implication of E2F2-dependent transcriptional activity in HD pathogenesis. Thus, E2F2 emerges as a new potential HD AO modifier factor.

Genetic evidence of two sibling species within the Contracoecum ogmorhini Johnson & Mawson 1941 complex (Nematoda; Anisakidae) from otariid seals in boreal and austral regions

Genetic variation of Contracaecum ogmorhini (sensu lato) populations from different otariid seals of the northern and southern hemisphere was studied on the basis of 18 enzyme loci as well as preliminary sequence analysis of the mitochondrial cyt b gene (260 bp). Samples were collected from Zalophus californianus in the boreal region and from Arctocephalus pusillus pusillus, A. pusillus doriferus and A. australis from the austral region. Marked genetic heterogeneity was found between C. ogmorhini (sensu lato) samples from the boreal and austral region, respectively. Two loci (Mdh-2 and NADHdh) showed fixed differences and a further three loci (Iddh, Mdh-1 and 6Pgdh) were highly differentiated between boreal and austral samples. Their average genetic distance was DNei = 0.36 at isozyme level. At mitochondrial DNA level, an average proportion of nucleotide substitution of 3.7% was observed. These findings support the existence of two distinct sibling species, for which the names C. ogmorhini (sensu stricto) and C. margolisi n. sp., respectively, for the austral and boreal taxon, are proposed. A description for C. margolisi n. sp. is provided. No diagnostic morphological characters have so far been detected; on the other hand, two enzyme loci, Mdh-2 and NADHdh, fully diagnostic between the two species, can be used for the routine identification of males, females and larval stages. Mirounga leonina was found to host C. ogmorhini (s.s.) inmixed infections with C. osculatum (s.l.) (of which C. ogmorhini (s.l.) was in the past considered to be a synonym) and C. miroungae; no hybrid genotypes were found,confirming the reproductive isolation of these three anisakid species. The hosts and geographical range so far recorded for C. margolisi n. sp. and C. ogmorhini (s.s.) are given.

Genetic differentiation among Atlantic cod (Gadus morhua) from Browns Bank, Georges Bank, and Nantucket Shoals

This study examines genetic variation at five microsatellite loci and at the vesicle membrane protein locus, pantophysin, of Atlantic cod (Gadus morhua) from Browns Bank, Georges Bank, and Nantucket Shoals. The Nantucket Shoals sample represents the first time cod south of Georges Bank have been genetically evaluated. Heterogeneity of allelic distribution was not observed (P>0.05) between two temporally separated Georges Bank samples indicating potential genetic stability of Georges Bank cod. When Bonferroni corrections (α=0.05, P<0.017) were applied to pairwise measures of population differentiation and estimates of FST, significance was observed between Nantucket Shoals and Georges Bank cod and also between Nantucket Shoals and Browns Bank cod. However, neither significant differentiation nor significant estimates of FST were observed between Georges Bank and the Browns Bank cod. Our research suggests that the cod spawning on Nantucket Shoals are genetically differentiated from cod spawning on Browns Bank and Georges Bank. Managers may wish to consider Nantucket Shoals cod a separate stock for assessment and management purposes in the future.

Protein polymorphism and genetic divergence in slow loris (Genus nycticebus)

In this study, protein electrophoresis was assayed to detect genetic variation in Genus Nycticebus. A total of 29 samples (2 N. coucang and 27 N. pygmaeus) were analyzed for 42 genetic loci. In the 27 samples of N. pygmaeus, 4 loci were observed to be polymorphic. Therefore, the estimated P value (proportion of polymorphic loci) is 0.095, the A value (average number of alleles each locus) is 1.045, and the H value (mean individual heterozygosity) is 0.040. After comparing the H of N. pygmaeus with those of other primates reported, we found that the protein variation in N. pygmaeus is slightly lower than the average level. Additionally, we also observed obvious allele difference between N. pygmaeus and N. coucang. There are no shared alleles between these two species in eight loci. The NEI's genetic distance between them was calculated as 0.2541, which falls in the spectrum of genetic difference between species in primates.

Allozyme variation in four hatchery populations of Thai pangas, Pangasius hypophthalmus in Bogra, Bangladesh

Genetic variation of four hatchery stocks of Thai pangas, Pangasius hypophthalmu [sic] of Bogra region, Bangladesh was studied from 1 January 2002 to 31 December 2003. Muscle samples were collected for allozyme analysis from four (Bhai-Bhai, Jahangir, Belal and Bhai-Bon) different hatchery populations. For allozyme electrophoresis, eight enzymes were used and 11 loci viz. Adh-1*, Est-1*, GJpdh-1*, Gpi-1*, Gpi-2*, Jdhp-1*, Ldh-1*, Ldh-2*, Mdh-1*, Mdh-2*and Pgm* were identified, of which three loci (Est-1*, Gpi-2*, G3pdh-1 and Pgm*) were polymorphic in all the four populations. The mean proportion of polymorphic loci per population and the mean proportion of heterozygous loci per individual was 36.36% and 13.33, respectively for all the population studied. The highest variability measured by the mean number of alleles per locus was 1.545 in Bhai-Bon hatchery population. Based on Nei's (1972) genetic distance, the dendrogram (UPGMA) shows that four populations have made two clusters by D-value (D=0.043). Bhai-Bhai and Jahangir hatchery populations have made cluster-I, and Belal Uddin and Bhai-Bon hatchery populations formed cluster-II. Among the four populations, BhaiBhai and Jahangir hatchery populations were differentiated from each other by the D-value of 0.013, and Belal Uddin and Bhai-Bon populations were differentiated from each other by the D-value of 0.002, which suggests that the four populations may be fallen into the local population or race.