Genetic diversity of MHC class II <Emphasis Type="Italic">DRB</Emphasis> alleles in the continental and Japanese populations of the sable <Emphasis Type="Italic">Martes zibellina</Emphasis> (Mustelidae,Carnivora, Mammalia)

The sable (Martes zibellina) is a medium-sized mustelid inhabiting forest environments in Siberia, northern China, the Korean Peninsula, and Hokkaido Island, Japan. To further understand the molecular evolution of the major histocompatibility complex (MHC), we sequenced part of exon 2 in MHC class II DRB genes, including codons encoding the antigen binding site, from 33 individuals from continental Eurasia and Japan. We identified 16 MHC class II DRB alleles (Mazi-DRBs), some of which were geographically restricted and others broadly distributed, and eight putative pseudogenes. A single-breakpoint recombination analysis detected a recombination site in the middle of exon 2. A mixed effects model of evolution analysis identified five amino acid sites presumably under positive selection. These sites were all located in the region 3′ to the recombination site, suggesting that positive selection and recombination could be committed to the diversity of the M. zibellina DRB gene. In a Bayesian phylogenetic tree, all Mazi-DRBs and the presumed pseudogenes grouped within a Mustelidae clade. The Mazi-DRBs showed trans-species polymorphism, with some alleles most closely related to alleles from other mustelid species. This result suggests that the sable DRBs have evolved under long-lasting balancing selection.     

Population genetic diversity and geographical differentiation of MHC class II DAB genes in the vulnerable Chinese egret (<Emphasis Type="Italic">Egretta eulophotes</Emphasis>)

Major histocompatibility complex (MHC) genes are excellent markers for the study of adaptive genetic variation occurring over different geographical scales. The Chinese egret (Egretta eulophotes) is a vulnerable ardeid species with an estimated global population of 2600–3400 individuals. In this study, we sampled 172 individuals of this egret (approximately 6 % of the global population) from five natural populations that span the entire distribution range of this species in China. We examined their population genetic diversity and geographical differentiation at three MHC class II DAB genes by identifying eight exon 2 alleles at Egeu-DAB1, eight at Egeu-DAB2 and four at Egeu-DAB3. Allelic distributions at each of these three Egeu-DAB loci varied substantially within the five populations, while levels of genetic diversity varied slightly among the populations. Analysis of molecular variance showed low but significant genetic differentiation among five populations at all three Egeu-DAB loci (haplotype-based ?_ST: 0.029, 0.020 and 0.042; and distance-based ?_ST: 0.036, 0.027 and 0.043, respectively; all P < 0.01). The Mantel test suggested that this significant population genetic differentiation was likely due to an isolation-by-distance pattern of MHC evolution. However, the phylogenetic analyses and the Bayesian clustering analysis based on the three Egeu-DAB loci indicated that there was little geographical structuring of the genetic differentiation among five populations. These results provide fundamental population information for the conservation genetics of the vulnerable Chinese egret.     

Sequence diversity of MHC class-II <Emphasis Type="BoldItalic">DRB</Emphasis> gene in gazelles (<Emphasis Type="BoldItalic">Gazella subgutturosa</Emphasis>) raised in Sanliurfa of Turkey

In this study, we aimed to assess the sequence diversity of major histocompatibility complex (MHC) class-II DRB gene at exon 2 in gazelles raised in Sanliurfa Province of Turkey. Twenty DNA samples isolated from gazelles (Gazella subgutturosa) were used for sequencing exon 2 of MHC class-II DRB gene. Target region was amplified by polymerase chain reaction (PCR) and their products were directly sequenced. Nine of these 20 samples yielded unambiguously readable sequences. Three of the nine samples were homozygotes and each showed different sequences. A 262-bp sequence obtained from the three homozygote samples were submitted to GenBank (accession numbers: KC309405, KC309406 and KC309407). Using an allele specific PCR, we detected 10 additional haplotypes. Among 13 haplotypes, 45 nucleotide positions were polymorphic and most of the polymorphic nucleotide positions localized at peptide-binding region (PBR). Rates of nonsynonymous substitutions were significantly higher than synonymous substitutions at PBR. Phylogenetic analysis of the haplotypes showed that 10 haplotypes of the gazelles were clustered together while three were clustered with ovine and bovine haplotypes. The results indicated that at least 13 haplotypes at exon 2 of MHC class-II DRB gene were showing high degree of nucleotide and amino acid diversity, and certain haplotypes of G. subgutturosa were more similar to haplotypes from sheep or cattle than to each other. Rates of synonymous and nonsynonymous substitutions suggested that positive selection was a driving force for diversity at this locus in G. subgutturosa.     

The intraspecies variability of <Emphasis Type="Italic">Spermophilus relictus</Emphasis> sensu lato Kaschkarov 1923 (Sciuridae,Rodentia)

The intraspecies variability of Spermophilus relictus sensu lato was studied based on 27 measurements of skulls from 67 specimens of relict ground squirrels (S. relictus) and 66 specimens of Tien Shan ground squirrels (S. ralli), as well as six specimens of the relict ground squirrel from the Gissar Ridge (Tien Shan). A colorimetric analysis of skins of relict ground squirrels (19 specimens, including three individuals from the Gissar Ridge) and Tien Shan ground squirrels (19 specimens) was made. Significant intraspecies variability was found in the relict and Issyk-Kul ground squirrels, whereas the interspecies differences were small, raising questions about the species independence of the Tien Shan ground squirrel (S. ralli).     

Effects of <Emphasis Type="Italic">HLA</Emphasis>-<Emphasis Type="Italic">DRB1</Emphasis>/<Emphasis Type="Italic">DQB1</Emphasis> Genetic Variants on Neuroimaging in Healthy,Mild Cognitive Impairment,and Alzheimer’s Disease Cohorts

Alzheimer’s disease (AD) is the most common form of dementia and exhibits a considerable level of heritability. Previous association studies gave evidence for the associations of HLA-DRB1/DQB1 alleles with AD. However, how and when the gene variants in HLA-DRB1/DQB1 function in AD pathogenesis has yet to be determined. Here, we firstly investigated the association of gene variants in HLA-DRB1/DQB1 alleles and AD related brain structure on magnetic resonance imaging (MRI) in a large sample from the Alzheimer’s Disease Neuroimaging Initiative (ADNI). We selected hippocampus, subregion, parahippocampus, posterior cingulate, precuneus, middle temporal, entorhinal cortex, and amygdala as regions of interest (ROIs). Twelve SNPs in HLA-DRB1/DQB1 were identified in the dataset following quality control measures. In the total group hybrid population analysis, our study (rs35445101, rs1130399, and rs28746809) were associated with the smaller baseline volume of the left posterior cingulate and rs2854275 was associated with the larger baseline volume of the left posterior cingulate. Furthermore, we detected the above four associations in mild cognitive impairment (MCI) sub-group analysis, and two risk loci (rs35445101 and rs1130399) were also the smaller baseline volume of the left posterior cingulate in (NC) sub-group analysis. Our study suggested that HLA-DRB1/DQB1 gene variants appeared to modulate the alteration of the left posterior cingulate volume, hence modulating the susceptibility of AD.     

Development and characterization of 24 chloroplast microsatellite markers for two species of <Emphasis Type="Italic">Eranthis</Emphasis> (Ranunculaceae)

Chloroplast microsatellites for two Korean endemic species, Eranthis byunsanensis and E. pungdoensis (Ranunculaceae), were isolated to address the questions of their distributional patterns and evolutionary relationships, using next-generation sequencing. Twenty-four polymorphic chloroplast microsatellite markers for these two species were developed, and then characterized in 65 individuals (55 individuals of E. byunsanensis and 10 individuals of E. pungdoensis). The number of alleles per locus ranged from 2 to 9; the average number of alleles across all the loci scored 4.792. The unbiased diversity per locus ranged from 0.089 to 0.880; the unbiased diversity averaged over all the loci was 0.646. The developed markers were successfully amplified for three congeneric species, E. stellata, E. pinnatifida, and E. longistipitata. The markers developed in this study can provide a valuable and important tool for understanding genetic variations, population structures, evolutionary histories and phylogeography of E. byunsanensis, E. pungdoensis, and related species.     

Isolation and Characterization of Microsatellite Markers for <Emphasis Type="Italic">Viola websteri</Emphasis> (Violaceae) and Cross-Species Amplification within the Genus <Emphasis Type="Italic">Viola</Emphasis>

We isolated and characterized microsatellite loci in Viola websteri (Violaceae), an endangered species from Korea and endemic to Northeast Asia. A total of 27 microsatellite loci were developed and tested in Korean and Chinese populations. The number of alleles per locus varied from two to eight. The observed and expected heterozygosities within two populations were 0.000 to 1.000 and 0.080 to 0.816, respectively. Korean and Chinese populations were clearly distinguished by the private alleles from 16 loci. A total of 21 loci out of the 27 developed loci were successfully cross-amplified in 39 other Viola species. We believe that these microsatellite loci will be useful for future studies on genetic diversity and population structure of V. websteri, as well as other Viola species.     

Current state of the genetic polymorphism in spring barley (<Emphasis Type="Italic">Hordeum vulgare</Emphasis> L.) from Russia assessed by the alleles of hordein-coding loci

Starch gel electrophoresis was performed to study the polymorphism of hordeins encoded by the Hrd A, Hrd B, and Hrd F loci in 211 varieties of spring barley. For 41 of these varieties, the genetic formulas were established for the first time. In the two samples of varieties, the comparative analysis of allelic diversity and allele frequencies of hordein-coding loci was carried out. The first sample consisted of 101 spring barley varieties approved for the use on the territory of the Russian Federation in 1999, while the second sample included 160 spring barley varieties that were approved in 2014; 49 of these varieties were common for both samples. It is demonstrated that the current tendency to reduction of the proportion of heterogeneous spring barley varieties is mainly due to the introduction of foreign varieties homogeneous for the hordein-coding loci. At the same time, there is an increase in polymorphism of hordein-coding loci in modern spring barley varieties. The number of alleles for the Hrd A locus increased by five alleles, and for the Hrd B locus, by nine alleles. Along with the alleles recorded earlier in barley landrace populations and varieties bred in 20th century, three novel alleles of the Hrd A locus and four alleles of the Hrd B locus were identified. The number of alleles of the Hrd F locus remained unchanged (four), and the changes in their frequencies were small. At the same time, the changes in frequency observed for some alleles of the Hrd A and Hrd B loci were statistically significant. All newly identified alleles of hordein-coding loci were found with low frequencies (from 0.003 to 0.006), so despite the increased number of alleles, no statistically significant increase in genetic diversity in terms of μ and PIC indices was observed.     

Low molecular weight glutenin subunit gene composition at <Emphasis Type="Italic">Glu</Emphasis>-<Emphasis Type="Italic">D3</Emphasis> loci of <Emphasis Type="Italic">Aegilops tauschii</Emphasis> and common wheat and a further view of wheat evolution

Key message

A comprehensive comparison of LMW-GS genes between Ae. tauschii and its progeny common wheat.

Abstract

Low molecular weight glutenin subunits (LMW-GSs) are determinant of wheat flour processing quality. However, the LMW-GS gene composition in Aegilops tauschii, the wheat D genome progenitor, has not been comprehensively elucidated and the impact of allohexaploidization on the Glu-D3 locus remains elusive. In this work, using the LMW-GS gene molecular marker system and the full-length gene-cloning method, LMW-GS genes at the Glu-D3 loci of 218 Ae. tauschii and 173 common wheat (Triticum aestivum L.) were characterized. Each Ae. tauschii contained 11 LMW-GS genes, and the whole collection was divided into 25 haplotypes (AeH01–AeH25). The Glu-D3 locus in common wheat lacked the LMW-GS genes D3-417, D3-507 and D3-552, but shared eight genes of identical open reading frame (ORF) sequences when compared to that of Ae. tauschii. Therefore, the allohexaploidization induces deletions, but exerts no influence on LMW-GS gene coding sequences at the Glu-D3 locus. 92.17% Ae. tauschii had 7-9 LMW-GSs, more than the six subunits in common wheat. The haplotypes AeH16, AeH20 and AeH23 of Ae. tauschii ssp. strangulate distributed in southeastern Caspian Iran were the main putative D genome donor of common wheat. These results facilitate the utilization of the Ae. tauschii glutenin gene resources and the understanding of wheat evolution.

     

Acid phosphatase gene <Emphasis Type="Italic">GmHAD1</Emphasis> linked to low phosphorus tolerance in soybean,through fine mapping

Key message

Map-based cloning identified GmHAD1, a gene which encodes a HAD-like acid phosphatase, associated with soybean tolerance to low phosphorus stress.

Abstract

Phosphorus (P) deficiency in soils is a major limiting factor for crop growth worldwide. Plants may adapt to low phosphorus (LP) conditions via changes to root morphology, including the number, length, orientation, and branching of the principal root classes. To elucidate the genetic mechanisms for LP tolerance in soybean, quantitative trait loci (QTL) related to root morphology responses to LP were identified via hydroponic experiments. In total, we identified 14 major loci associated with these traits in a RIL population. The log-likelihood scores ranged from 2.81 to 7.43, explaining 4.23–13.98% of phenotypic variance. A major locus on chromosome 08, named qP8-2, was co-localized with an important P efficiency QTL (qPE8), containing phosphatase genes GmACP1 and GmACP2. Another major locus on chromosome 10 named qP10-2 explained 4.80–13.98% of the total phenotypic variance in root morphology. The qP10-2 contains GmHAD1, a gene which encodes an acid phosphatase. In the transgenic soybean hairy roots, GmHAD1 overexpression increased P efficiency by 8.4–16.5% relative to the control. Transgenic Arabidopsis plants had higher biomass than wild-type plants across both short- and long-term P reduction. These results suggest that GmHAD1, an acid phosphatase gene, improved the utilization of organic phosphate by soybean and Arabidopsis plants.

     

Low MHC variation in the endangered Galápagos penguin (<Emphasis Type="Italic">Spheniscus mendiculus</Emphasis>)

The major histocompatibility complex (MHC) is one of the most polymorphic regions of the genome, likely due to balancing selection acting to maintain alleles over time. Lack of MHC variability has been attributed to factors such as genetic drift in small populations and relaxed selection pressure. The Galápagos penguin (Spheniscus mendiculus), endemic to the Galápagos Islands, is the only penguin that occurs on the equator. It relies upon cold, nutrient-rich upwellings and experiences severe population declines when ocean temperatures rise during El Niño events. These bottlenecks, occurring in an already small population, have likely resulted in reduced genetic diversity in this species. In this study, we used MHC class II exon 2 sequence data from a DRB1-like gene to characterize the amount of genetic variation at the MHC in 30 Galápagos penguins, as well as one Magellanic penguin (S. magellanicus) and two king penguins (Aptenodytes patagonicus), and compared it to that in five other penguin species for which published data exist. We found that the Galápagos penguin had the lowest MHC diversity (as measured by number of polymorphic sites and average divergence among alleles) of the eight penguin species studied. A phylogenetic analysis showed that Galápagos penguin MHC sequences are most closely related to Humboldt penguin (Spheniscus humboldti) sequences, its putative sister species based on other loci. An excess of non-synonymous mutations and a pattern of trans-specific evolution in the neighbor-joining tree suggest that selection is acting on the penguin MHC.     

Integration of sudden death syndrome resistance loci in the soybean genome

Key message

Complexity and inconsistencies in resistance mapping publications of soybean sudden death syndrome (SDS) result in interpretation difficulty. This review integrates SDS mapping literature and proposes a new nomenclature system for reproducible SDS resistance loci.

Abstract

Soybean resistance to sudden death syndrome (SDS) is composed of foliar resistance to phytotoxins and root resistance to pathogen invasion. There are more than 80 quantitative trait loci (QTL) and dozens of single nucleotide polymorphisms (SNPs) associated with soybean resistance to SDS. The validity of these QTL and SNPs is questionable because of the complexity in phenotyping methodologies, the disease synergism between SDS and soybean cyst nematode (SCN), the variability from the interactions between soybean genotypes and environments, and the inconsistencies in the QTL nomenclature. This review organizes SDS mapping results and proposes the Rfv (resistance to Fusarium virguliforme) nomenclature based on supporting criteria described in the text. Among ten reproducible loci receiving our Rfv nomenclature, Rfv18-01 is mostly supported by field studies and it co-localizes to the SCN resistance locus rhg1. The possibility that Rfv18-01 is a pleiotropic resistance locus and the concern about Rfv18-01 being confounded with Rhg1 is discussed. On the other hand, Rfv06-01, Rfv06-02, Rfv09-01, Rfv13-01, and Rfv16-01 were identified both by screening soybean leaves against phytotoxic culture filtrates and by evaluating SDS severity in fields. Future phenotyping using leaf- and root-specific resistance screening methodologies may improve the precision of SDS resistance, and advanced genetic studies may further clarify the interactions among soybean genotypes, F. virguliforme, SCN, and environments. The review provides a summary of the SDS resistance literature and proposes a framework for communicating SDS resistance loci for future research considering molecular interactions and genetic breeding for soybean SDS resistance.

     

Characterization of an IAA-glucose hydrolase gene <Emphasis Type="Italic">TaTGW6</Emphasis> associated with grain weight in common wheat (<Emphasis Type="Italic">Triticum aestivum</Emphasis> L.)

In rice, the TGW6 gene determines grain weight and encodes a protein with indole-3-acetic acid (IAA)-glucose hydrolase activity. Its homolog in wheat, TaTGW6, is considered as a candidate gene related to grain development. To amplify this gene, we designed primers based on a homologous conserved domain of the rice TGW6 gene. Sequence analysis indicated that TaTGW6 comprises only one exon, with 1656 bp in total and an open reading frame of 1035 bp. Three alleles at TaTGW6 locus detected by the primer pair TG23 were designated as TaTGW6-a, TaTGW6-b and TaTGW6-c, respectively. Compared with TaTGW6-a, TaTGW6-b had a 6-bp InDel at the position 170 downstream of initiation codon, and TaTGW6-c was a null mutant. Both TaTGW6-b and TaTGW6-c could significantly increase grain size and weight other than TaTGW6-a; however, the former two alleles showed a low frequency distribution in modern varieties. TaTGW6 was located on chromosome 4AL using a recombinant inbred line population and a set of Chinese Spring nullisomic-tetrasomic lines. It was linked to the SSR locus Xbarc1047 with a genetic distance of 6.62 cM and explained 15.8–21.0 % of phenotypic variation of grain weight in four environments. Association analysis using a natural population and Chinese wheat mini-core collections additionally validated the relationship of TaTGW6 with grain weight; the gene could explain 7.7–12.4 % of phenotypic variation in three environments. Quantitative real-time PCR revealed that TaTGW6-b showed relatively lower expression than TaTGW6-a in immature grain at 20 and 30 days post-anthesis and in mature grain. The low expression of TaTGW6 generally associated with low IAA content, but with high grain weight. The novel functional marker, designated as TG23, can be used for marker-assisted selection to improve grain weight in wheat and also provides insights into the regulatory mechanism underlying grain weight.     

Nucleotide polymorphisms of candidate genes of adaptive significance in the ural populations of <Emphasis Type="Italic">Larix sibirica</Emphasis> Ledeb.

The objectives of conservation and sustainable forest management require in depth study of genomes of woody plants and definition of their intraspecific genetic diversity. In recent years, an approach was developed based on the study of “candidate genes” that can potentially be involved in the formation of adaptive traits. In this study, we investigated nucleotide polymorphism of several adaptive candidate genes in the populations of Siberian larch (Larix sibirica Ledeb.) in the Urals. Representatives of this genus are among the most valuable and widely distributed forest tree species in Russia. From ten selected gene loci in the genome of L. sibirica, we isolated and investigated three loci, one of which (ABA-WDS) was sequenced in L. sibirica for the first time. The total length of the analyzed sequence in each individual amounted to 2865 bp. The length of locus alignment was from 360 bp to 1395 bp. In total, we identified 200 polymorphic positions. The most conservative is locus 4CL1-363, and the most polymorphic is locus sSPcDFD040B03103-274. The studied populations of L. sibirica are characterized by a high level of nucleotide polymorphism in comparison with other species and genuses (Picea, Pinus, Pseudotsuga, Abies) conifers plants (Hd = 0.896; π = 0.007; θ_W = 0.015). The most selectively neutral polymorphism (D _T =–0.997) was attributed to locus 4CL1-363, and polymorphism with high probability of adaptability (D _T =–1.807) was determined for the ABA-WDS locus. We identified 54 SNP markers, only five of which were nonsynonymous (9.26%) replacements. The average frequency of SNPs in the three studied loci of L. sibirica was one SNP in 53 bp. We detected unique SNP markers for eight populations, which could potentially be used to identify populations. Populations that are characterized by the highest number of unique SNP markers can be recommended for selection in order to preserve the gene pool of the species.     

Genetic structure of the populations of <Emphasis Type="Italic">Dactylorhiza ochroleuca</Emphasis> and <Emphasis Type="Italic">D. incarnata</Emphasis> (Orchidaceae) in the area of their joint growth in Russia and Belarus

We carried out an allozyme analysis to investigate polymorphism and genetic structure of the populations of D. incarnata and D. ochroleuca in regions of their joint growth in Russia and Belarus. We found that D. ochroleuca individuals in the populations of the Urals and Siberia, which are distant fragments from the main range of the species, do not differ significantly from individuals within the main part of the area (Belarus) on the basis of the allelic composition of eight gene loci. We revealed that D. ochroleuca and D. incarnata are differentiated by different alleles of the GDH locus. Thus, we established a genetic marker suitable to distinguish these closely related taxa. In addition to the GDH locus, D. ochroleuca and D. incarnata in the places of their joint growth, differ in the allelic structure of the PGI and NADHD loci. D. incarnata from the Urals and Siberia were polymorphic for both loci, and individuals from Belarus were polymorphic for one locus (PGI). In contrast, all D. ochroleuca individuals growing in sympatric populations with polymorphic D. incarnata were homozygous for the same alleles. Thus, comparison of the genetic structure of D. ochroleuca and D. incarnata points to the existence of a genetic isolation and a functioning isolation mechanism even under conditions of their joint growth. We found that the GDH locus in D. incarnata is polymorphic only in populations which grow together with D. ochroleuca, with exception a few examples. Thus, we conclude that variability of the GDH locus in D. incarnata is associated with hybridization with D. ochroleuca.     

Genome-wide association study of heading and flowering dates and construction of its prediction equation in Chinese common wheat

Heading date is one of the most important traits in wheat breeding as it affects adaptation and yield potential. A genome-wide association study (GWAS) using the 90 K iSelect SNP genotyping assay indicated that a total of 306 loci were significantly associated with heading and flowering dates in 13 environments in Chinese common wheat from the Yellow and Huai wheat region. Of these, 105 loci were significantly correlated with both heading and flowering dates and were found in clusters on chromosomes 2, 5, 6, and 7. Based on differences in distribution of the vernalization and photoperiod genes among chromosomes, arms, or block regions, 13 novel, environmentally stable genetic loci were associated with heading and flowering dates, including RAC875_c41145_189 on 1DS, RAC875_c50422_299 on 2BL, and RAC875_c48703_148 on 2DS, that accounted for more than 20% phenotypic variance explained (PVE) of the heading/flowering date in at least four environments. GWAS and t test of a combination of SNPs and vernalization and photoperiod alleles indicated that the Vrn-B1, Vrn-D1, and Ppd-D1 genes significantly affect heading and flowering dates in Chinese common wheat. Based on the association of heading and flowering dates with the vernalization and photoperiod alleles at seven loci and three significant SNPs, optimal linear regression equations were established, which show that of the seven loci, the Ppd-D1 gene plays the most important role in modulating heading and flowering dates in Chinese wheat, followed by Vrn-B1 and Vrn-D1. Additionally, three novel genetic loci (RAC875_c41145_189, Excalibur_c60164_137, and RAC875_c50422_299) also show important effect on heading and flowering dates. Therefore, Ppd-D1, Vrn-B1, Vrn-D1, and the novel genetic loci should be further investigated in terms of improving heading and flowering dates in Chinese wheat. Further quantitative analysis of an F₁₀ recombinant inbred lines population identified a major QTL that controls heading and flowering dates within the Ppd-D1 locus with PVEs of 28.4% and 34.0%, respectively; this QTL was also significantly associated with spike length, peduncle length, fertile spikelets number, cold resistance, and tiller number.     

Genetic markers of adaptive processes in the Far Eastern pink salmon <Emphasis Type="Italic">Oncorhynchus gorbuscha</Emphasis>: Allelic diversity at the locus of major histocompatibility complex MHC I-<Emphasis Type="Italic">A1</Emphasis>

To clarify allelic diversity at the locus of major histocompatibility complex MHC class I-A1 in the Far Eastern pink salmon Oncorhynchus gorbuscha, sequencing of the electrophoretic alleles isolated from the gel (DGGE alleles) was performed. In 47 individuals, the genotypes of which consisted of ten DGGE alleles, 18 MHC I-A1 nucleotide sequences were revealed, and thus, eight cryptic alleles not detected by electrophoresis were identified. Eleven of these alleles were identified earlier in pink salmon from Hokkaido, Alaska, and British Columbia, and seven, possibly, were unique to the populations from some Far Eastern regions. Six of the previously determined DGGE alleles corresponded to more than one nucleotide sequence. However, the sequences attributed to the same DGGE allele differed on average by less than 1 nucleotide. These findings point to sufficient sensitivity of the DGGE method, although the genetic diversity and differentiation estimates obtained with it will obviously be somewhat underestimated. Considerable predominance of nonsynonymous substitutions over the synonymous ones in the codons of the MHC I-A1 antigen-binding site confirms the presence of positive selection aimed at providing the population resistance to local spectrum of pathogens. Refinement of the allelic composition of the adaptively important MHC genetic marker will contribute to more complete understanding of the adaptive genetic structure of pink salmon as an important element of the overall population structure of the species.