Novel Variants in the Putative Peroxisome Proliferator‐activated Receptor γ Promoter and Relationships with Obesity in Men

Yet unidentified variants within the peroxisome proliferator‐activated receptor γ (PPARγ) 2 promoter may explain the inconsistent reports on associations between variants in the coding region and obesity or diabetes. Thus, we examined the putative PPARγ2 promoter (?3371 to +43 bp) for variants in 83 subjects with obesity or type 2 diabetes. We identified eight variants, seven of which were novel, including ?792A>G, ?816C>T, ?882T>C, ?1505G>A, ?1881C>T, ?1884T>A, ?2604T>C, and ?2953A>G. The variants ?816C>T, ?1505G>A, ?1881C>T, and ?2604T>C were in total linkage disequilibrium, and there was a high degree of linkage disequilibrium between several of the novel variants and Pro12Ala. The novel variants were, together with Pro12Ala and 1431C>T, examined for relationships with obesity among 234 men with early‐onset obesity with a BMI at age ～20 years of 33.2 ± 2.5 kg/m² and 323 nonobese men with a BMI of 21.7 ± 2.5 kg/m², who were also reexamined after ～29 years. The prevalence of the identified variants was not significantly different between the two groups, and the variants did not affect changes in BMI over time. In conclusion, the identified novel variants in the PPARγ2 promoter region do not explain the reported discrepancies in the association of previously identified variants with obesity and type 2 diabetes.     

Association between single nucleotide polymorphisms of TPH1 and TPH2 genes,and depressive disorders

Tryptophan catabolites pathway disorders are observed in patients with depression. Moreover, single nucleotide polymorphisms of tryptophan hydroxylase genes may modulate the risk of depression occurrence. The objective of our study was to confirm the association between the presence of polymorphic variants of TPH1 and TPH2 genes, and the development of depressive disorders. Six polymorphisms were selected: c.804‐7C>A (rs10488682), c.‐1668T>A (rs623580), c.803+221C>A (rs1800532), c.‐173A>T (rs1799913)—TPH1, c.‐1449C>A (rs7963803), and c.‐844G>T (rs4570625)—TPH2. A total of 510 DNA samples (230 controls and 280 patients) were genotyped using TaqMan probes. Among the studied polymoorphisms, the G/G genotype and G allele of c.804‐7C>A—TPH1, the T/T homozygote of c.803+221C>A—TPH1, the A/A genotype and A allele of c.1668T>A—TPH1, the G/G homozygote and G allele of c.‐844G>T—TPH2, and the C/A heterozygote and A allele of c.‐1449C>A—TPH2 were associated with the occurrence of depression. However, the T/T homozygote of c.‐1668T>A—TPH1, the G/T heterozygote and T allele of c.‐844G>T—TPH2, and the C/C homozygote and C allele of c.‐1449C>A—TPH2 decreased the risk of development of depressive disorders . Each of the studied polymorphisms modulated the risk of depression for selected genotypes and alleles. These results support the hypothesis regarding the involvement of the pathway in the pathogenesis of depression.     

SNP genotypes of olfactory receptor genes associated with olfactory ability in German Shepherd dogs

To find out the relationship between SNP genotypes of canine olfactory receptor genes and olfactory ability, 28 males and 20 females from German Shepherd dogs in police service were scored by odor detection tests and analyzed using the Beckman GenomeLab SNPstream. The representative 22 SNP loci from the exonic regions of 12 olfactory receptor genes were investigated, and three kinds of odor (human, ice drug and trinitrotoluene) were detected. The results showed that the SNP genotypes at the OR10H1‐like:c.632C>T, OR10H1‐like:c.770A>T, OR2K2‐like:c.518G>A, OR4C11‐like:c.511T>G and OR4C11‐like:c.692G>A loci had a statistically significant effect on the scenting abilities (P < 0.001). The kind of odor influenced the performances of the dogs (P < 0.001). In addition, there were interactions between genotype and the kind of odor at the following loci: OR10H1‐like:c.632C>T, OR10H1‐like:c.770A>T, OR4C11‐like:c.511T>G and OR4C11‐like:c.692G>A (P < 0.001). The dogs with genotype CC at the OR10H1‐like:c.632C>T, genotype AA at the OR10H1‐like:c.770A>T, genotype TT at the OR4C11‐like:c.511T>G and genotype GG at the OR4C11‐like:c.692G>A loci did better at detecting the ice drug. We concluded that there was linkage between certain SNP genotypes and the olfactory ability of dogs and that SNP genotypes might be useful in determining dogs' scenting potential.     

Association between METTL3 gene polymorphisms and neuroblastoma susceptibility: A nine‐centre case‐control study

Neuroblastoma ranks as the most commonly seen and deadly solid tumour in infancy. The aberrant activity of m⁶A‐RNA methyltransferase METTL3 is involved in human cancers. Therefore, functional genetic variants in the METTL3 gene may contribute to neuroblastoma risk. In the current nine‐centre case‐control study, we aimed to analyse the association between the METTL3 gene single nucleotide polymorphisms (SNPs) and neuroblastoma susceptibility. We genotyped four METTL3 gene SNPs (rs1061026 T>G, rs1061027 C>A, rs1139130 A>G, and rs1263801 G>C) in 968 neuroblastoma patients and 1814 controls in China. We found significant associations between these SNPs and neuroblastoma risk in neither single‐locus nor combined analyses. Interestingly, in the stratified analysis, we observed a significant risk association with rs1061027 AA in subgroups of children ≤ 18 months of age (adjusted OR = 1.87, 95% CI = 1.03‐3.41, P = .040) and females (adjusted OR = 1.86, 95% CI = 1.07‐3.24, P = .028). Overall, we identified a significant association between METTL3 gene rs1061027 C>A polymorphism and neuroblastoma risk in children ≤18 months of age and females. Our findings provide novel insights into the genetic determinants of neuroblastoma.     

Trans-activation-based risk assessment of <Emphasis Type="Italic">BRCA1</Emphasis> BRCT variants with unknown clinical significance

Background

Deleterious variants in the tumour suppressor BRCA1 are known to cause hereditary breast and ovarian cancer syndrome (HBOC). Missense variants in BRCA1 pose a challenge in clinical care, as their effect on protein functionality often remains unknown. Many of the pathogenic missense variants found in BRCA1 are located in the BRCA1 C-terminal (BRCT) domains, domains that are known to be vital for key functions such as homologous recombination repair, protein-protein interactions and trans-activation (TA). We investigated the TA activity of 12 BRCA1 variants of unknown clinical significance (VUSs) located in the BRCT domains to aid in the classification of these variants.

Results

Twelve BRCA1 VUSs were investigated using a modified version of the dual luciferase TA activity assay (TA assay) that yielded increased sensitivity and sample throughput. Variants were classified according to American College of Medical Genetics and Genomics (ACMG) criteria using TA assay results and available data. In combining our TA-assay results and available data, in accordance with the ACMG guidelines for variant classification, we proposed the following variant classifications: c.5100A>G, c.5326C>T, c.5348T>C and c.5477A>T as likely benign (class 2) variants. c.5075A>C, c.5116G>A and c.5513T>G were likely pathogenic (class 4), whereas c.5096G>A likely represents a likely pathogenic variant with moderate penetrance. Variants c.5123C>T, c.5125G>A, c.5131A>C and c.5504G>A remained classified as VUSs (class 3).

Conclusions

The modified TA assay provides efficient risk assessment of rare missense variants found in the BRCA1 BRCT-domains. We also report that increased post-transfection incubation time yielded a significant increase in TA assay sensitivity.

     

Haplotype analyses, mechanism and evolution of common double mutants in the human LDL receptor gene

Familial hypercholesterolemia (FH), an autosomal dominant inherited disorder resulting in increased levels of circulating plasma low-density lipoprotein (LDL), tendon xanthomas and premature coronary artery disease (CAD), is caused by defects in the LDL receptor gene (LDLR). Three widespread LDLR alterations not causing FH (c.1061-8T>C, c.2177C>T and c.829G>A) and one mutation (c.12G>A) with narrow geographical distribution and thought to cause disease were investigated. In an attempt to improve knowledge on their origin, spread and possible selective effects, estimations of the ages of these variants (t generations) and haplotype analysis were performed by genotyping 86 healthy individuals and 98 FH patients in Spain for five LDLR SNPs: c.81T>C, c.1413G>A, c.1725C>T, c.1959T>C, and c.2232G>A; most patients carried two of these LDLR variants simultaneously. It was found that both the c.1061-8T>C (t = 54) and c.2177C>T alterations (t = 62) arose at about the same time (54 and 62 generations ago, respectively) in the CGCTG haplotype, while the c.12G>A mutation (t = 70) appeared in a CGCCG haplotype carrying an earlier c.829G>A alteration (t = 83). The estimated ages of selectively neutral alterations could explain their distribution by migrations. The origin of the c.12G>A mutation could be in the Iberian Peninsula; despite its estimated age, a low selective pressure could explain its conservation in Spain from where it could have spread to China and Mexico, since the sixteenth century through the Spanish/Portuguese colonial expeditions.     

The prion‐related protein (testis‐specific) gene (PRNT) is highly polymorphic in Portuguese sheep

The objective of this study was to search for polymorphisms in the ovine prion‐related protein (testis‐specific) gene (PRNT). Sampling included 567 sheep from eight Portuguese breeds. The PRNT gene‐coding region was analyzed by single‐strand conformation polymorphism and sequencing, allowing the identification of the first ovine PRNT polymorphisms, in codons 6, 38, 43 and 48: c.17C>T (p.Ser6Phe, which disrupts a consensus arginine‐X‐X‐serine/threonine motif); c.112G>C (p.Gly38>Arg); c.129T>C and c.144A>G (synonymous) respectively. Polymorphisms in codons 6, 38 and 48 occur simultaneously in 50.6% of the animals, 38.8% presenting as heterozygous. To study the distribution of the polymorphism in codon 43, a restriction fragment length polymorphism analysis was performed. Polymorphic variant c.129C, identified in 89.8% of the animals with 32.8% presented as heterozygous, was considered the wild genotype in Portuguese sheep. Eight different haplotypes which have comparable distribution in all breeds were identified for the PRNT gene. In conclusion, the PRNT coding region is highly polymorphic in sheep, unlike the prion protein 2 dublet gene (PRND), in which we previously found only one synonymous substitution (c.78G>A), in codon 26. The absence or reduced number of PRND heterozygotes (c.78G>A) was significantly associated with three PRNT haplotypes (17C‐112G‐129T‐144A,17CT‐112GC‐129CT‐144AG and 17T‐112C‐129C‐144G), and the only three animals found homozygous at c.78A had the 17C‐112G‐129C‐144A PRNT haplotype. These results constitute evidence of an association between polymorphic variation in PRND and PRNT genes, as has already been observed for PRND and prion protein gene (PRNP).     

Novel Mutation in the <Emphasis Type="Italic">MECP2</Emphasis> Gene Identified in a Group of Rett Syndrome Patients from Ukraine

Mutations in the MECP2 gene are known to cause Rett syndrome (RTT)—a neurodevelopmental disorder, one of the most common causes of intellectual disability in females, with an incidence of 1 in 10000–15000. We have investigated exons 3 and 4 of the MECP2 gene, that coding MBD and TRD domains of the MeCP2 protein, in 21 RTT patients from Ukraine by PCR-DGGE analysis followed by Sanger sequencing of PCR fragments with abnormal migration profiles. In 13 of 21 (61.9%) patients 7 different mutations were identified one nonsense mutation—c. NC_000023.11:g.154031326G>A (MECP2:c.502C>T) and 4 missense mutation NC_000023.11:g.154031409G>T (MECP2:c.419C>T), NC_000023.11:g.154031355G>A (MECP2:c.473C>T), NC_000023.11:g.154031354A>C (MECP2:c.472A>C), NC_000023.11:g.154031431G>A (MECP2:c.397C>T) located in exon 4, a rare RTT-causing splice site mutation NC_000023.10:g.153296903T>G (MECP2:c.378-2A>C) in intron 3 and deletion NC_000023.10:g.1532 96079_153296122del44 in exon 4. The novel mutation MECP2:c.472A>C identified in our study in patients withclassic RTT phenotype leds to T158P substitution. It is one more confirmation of crucial role that 158 codon in MECP2 protein function.     

Identification of novel single nucleotide polymorphisms (SNPs) of the lipoprotein lipase (LPL) gene associated with fatty acid composition in Korean cattle

The lipoprotein lipase (LPL) gene can be considered a functional candidate gene that regulates fatty acid composition. In this study, genetic associations between fatty acid composition and exonic single nucleotide polymorphisms (SNPs) in the LPL gene were examined using 612 Korean cattle. We investigated the relationship between unsaturated fatty acids and five novel SNPs (c.322G>A, c.329A>T, c.527T>G, c.988C>T and c.1591G>A), and confirmed that three polymorphic SNPs (c.322G>A, c.329A>T and c.1591G>A) were associated with fatty acid composition. Korean cattle with an AA genotype of c.322G>A, c.329A>T, and GA genotype of c.1591G>A had higher levels of monounsaturated fatty acids and carcass traits (P < 0.05). Our findings confirmed that three novel SNPs we identified in the LPL gene can affect fatty acid composition and carcass traits. Therefore, selection for AA and GA genotypes should be recommended to genetically improve beef quality and flavor.     

Association between XPG polymorphisms and stomach cancer susceptibility in a Chinese population

Xeroderma pigmentosum group G (XPG) protein plays an important role in the DNA repair process by cutting the damaged DNA at the 3′ terminus. Previous studies have indicated some polymorphisms in the XPG gene are associated with stomach cancer susceptibility. We performed this hospital‐based case–control study to evaluate the association of four potentially functional XPG polymorphisms (rs2094258 C>T, rs751402 C>T, rs2296147 T>C and rs873601G>A) with stomach cancer susceptibility. The four single nucleotide polymorphisms (SNPs) were genotyped in 692 stomach cancer cases and 771 healthy controls. Logistic regression analysis was conducted, and odds ratios (ORs) and 95% confidence intervals (CIs) were used to assess the association of interest. Of the studied SNPs, XPG rs873601G>A polymorphism was found to significantly associate with stomach cancer susceptibility (AA versus GG/AG: OR = 1.31, 95% CI = 1.03–1.66, P = 0.027). Combined analysis of all SNPs revealed that the individuals with two of risk genotypes had a significantly increased stomach cancer risk (OR = 1.52, 95% CI = 1.13–2.06). In the stratification analysis, the association between the rs873601AA genotype and stomach cancer risk was observed in older group (>59 year), as well as patients with non‐cardia stomach cancer. Further combined analysis indicated men, smokers, or non‐drinkers more than one risk genotypes had a significantly increased stomach cancer risk. Our results indicate that XPG rs873601G>A polymorphism may be associated with the risk of stomach cancer. Further prospective studies with different ethnicities and large sample sizes are needed to validate our findings.     

Mutational screening of <Emphasis Type="BoldItalic">PKD2</Emphasis> gene in the north Indian polycystic kidney disease patients revealed 28 genetic variations

Polycystic kidney disease (PKD) is a systemic disorder which adds majority of renal patients to end stage renal disease. Autosomal dominant polycystic kidney disease (ADPKD) is more prevalent and leading cause of dialysis and kidney transplant. Linkage analysis revealed some closely linked loci, two of which are identified as PKD1, PKD2 and an unidentified locus to ADPKD. This study was performed using PCR and automated DNA sequencing in 84 cases and 80 controls to test potential candidature of PKD2 as underlying cause of PKD by in silico and statistical analyses. Two associated symptoms, hypertension (19%) and liver cyst (31%) have major contribution to PKD. Gender-based analysis revealed that familial female patients (27%) and familial male patients (33%) are more hypertensive. Liver cyst, the second major contributing symptom presented by large percentage of sporadic males (46%). Genetic screening of all 15 exons of PKD2 revealed eight pathogenic (c.854_854delG, c.915C>A, c.973C>T, c.1050_1050delC, c.1604_1604delT, c.1790T>C, c.2182_2183delAG, c.2224C>T) and eight likely pathogenic (g.11732A>G, c.646T>C, c.1354A>G, g.39212G>C, c.1789C>A, c.1849C>A, c.2164G>T, c.2494A>G) DNA sequence variants. In our study, 27.38% (23/84) cases shown pathogenic / likely pathogenic variants in PKD2 gene. Some regions of PKD2 prone for genetic variation suggested to be linked with disease pathogenesis. This noticeable hot spot regions hold higher frequency (50%) of pathogenic / likely pathogenic genetic variants constituting single nucleotide variants than large deletion and insertion that actually represents only 41.08% of coding sequence of PKD2. Statistically significant association for IVS3-22AA genotype was observed with PKD, while association of IVS4+62C>T was found insignificant.     

Eight novel CARD15 variants detected by DNA sequence analysis of the CARD15 gene in 111 patients with inflammatory bowel disease

We performed a limited DNA sequence analysis of the CARD15 gene in 89 patients with Crohn’s disease (CD), 19 patients with ulcerative colitis (UC), and three patients with indeterminate colitis (IC), who were heterozygous carriers of one of the common CARD15 mutations [c.2104C>T (p.R702W), c.2722G>C (p.G908R), or c.3019_3020insC (p.Leu1007fsX1008)], the c.2462+10A>C variant, or of a new amino acid substitution in the 3′-end of exon 4. CARD15 exons 4, 5, 6, 8, and 11 were amplified by PCR and completely sequenced, thereby theoretically covering 73.9% of the described CARD15 variants and 96.6% of the mutated alleles. Using this approach, eight novel amino acid substitutions [c.1171C>T (p.R391C), c.1387C>G (p.P463A), c.2138G>A (p.R713H), c.2278C>T (p.R760C), c.2368C>T (p.R790W), c.2371C>T (p.R791W), c.2475C>G (p.N825K), and c.2546C>T (p.A849V)] were detected in six CD and two IC patients, and one UC patient. A severe disease phenotype was observed especially in patients who are compound-heterozygous for a common and a novel CARD15 mutation.Schnitzler and Brand contributed equally     

Effect of FASN gene on milk yield and milk composition in the Chinese Holstein dairy population

Fatty acid synthase (FASN) is a multifunctional protein that catalyzes de novo synthesis of fatty acids in cells. It plays a key role in the lipid biosynthesis as well as in the general metabolism of all living animals. We herein investigated polymorphisms of FASN. As a result, six single nucleotide polymorphisms (SNPs) were found and then genotyped in 752 Chinese Holstein cows. It was found that g.17924A>G was non‐synonymous, g.13965 C>T, g.16907 T>C and g.18663T>C were synonymous mutations and two other two SNPs, g.8948 C>T (ss491228481) and g.14439T>C (rs133498277), were in intronic sequences of the gene. All such identified SNPs were found to be associated with milk yield and composition traits (P = 0.0441 to <0.0001). Significant additive and allele substitution effects were observed for three yield traits at all six loci as well (P < 0.05 to <0.01). Complete linkage disequilibrium among the five SNPs, with the exception of g.8948 C>T, was observed.     

Novel polymorphisms in PDLIM3 and PDLIM5 gene encoding Z‐line proteins increase risk of idiopathic dilated cardiomyopathy

Idiopathic dilated cardiomyopathy (IDCM), characterized by ventricular dilation and impaired systolic function, is a primary cardiomyopathy resulting in heart failure. During heart contraction, the Z‐line is responsible for transmitting force between sarcomeres and is also a hot spot for muscle cell signalling. Mutations in Z‐line proteins have been linked to cardiomyopathies in both humans and mice. Actinin‐associated LIM protein (ALP) and enigma homolog protein (ENH), encoded by PDLIM3 and PDLIM5, are components of the muscle cytoskeleton and localize to the Z‐line. A PDLIM3 or PDLIM5 deficiency in mice leads to dilated cardiomyopathy. Since PDLIM3 and PDLIM5 are candidate IDCM susceptibility genes, the current study aims to investigate whether polymorphisms within PDLIM3 and PDLIM5 could be correlated with IDCM. We designed a case‐control study, and exons of the PDLIM3 and PDLIM5 were amplified by polymerase chain reactions in 111 IDCM patients and 137 healthy controls. We found that five synonymous polymorphisms had statistical distribution differences between IDCM patients and controls, including rs4861669, rs4862543, c.731 + 131 T > G, c.1789‐3 C > T and rs7690296, according to genotype and allele distribution. Haplotype G‐C‐C‐C and A‐T‐C‐T (rs2306705, rs10866276, rs12644280 and rs4635850 synthesized) were regarded as risk factors for IDCM patients when compared with carriers of other haplotypes (all P < .05). Furthermore, IDCM patients with two novel polymorphisms (c.731 + 131 T > G and c.1789‐3 C > T) had lower systolic blood pressure. In conclusion, these five synonymous polymorphisms might constitute a genetic background that increases the risk of the development of IDCM in the Chinese Han population.     

CMAH genotyping survey for blood types A,B and C (AB) in purpose‐bred cats

In domestic cats, the AB blood group system consists of the three types A, B and C (also called AB). Mismatches can cause acute hemolytic transfusion reactions and hemolysis of the newborn (neonatal isoerythrolysis, NI). As blood types B and C are inherited recessively to A, breeders need to know the genotype to predict blood types in offspring and avoid NI. Several CMAH variants have been described as being associated with the b and a^c alleles, and different genotyping schemes exist. Here, we genotyped 2145 cats with the original SNV panel, including SNVs c.142G>A and ?‐53, and our new scheme, with SNVs c.179G>T, c.268T>A and c.1322delT, to differentiate types A and B and added the SNV for the common a^c (c.364C>T). Based upon the new scheme, all samples were assigned the correct genotype. No discordances appeared for the A allele, and new breed‐specific SNVs (c.179G>T, c.1322delT) for the b allele were discovered. Furthermore, the genotypes A/a^c (type A), a^c/a^c (C) and a^c/b (C) could be detected. We found the variant c.179G>T in additional breeds: Ragdoll, Siberian, Scottish Fold, Chartreux, Neva Masquerade, British Shorthair and Highlander. Also, the variant c.364C>T was detected in additional breeds: Bengal, British Shorthair, Maine Coon, and Scottish Fold. We conclude that our new SNV panel is superior in genotyping cats than the original SNV panel and assures correct assignments of types A, B and C to assist veterinary clinicians and breeders to recognize, confirm and avoid blood incompatibilities such as acute hemolytic transfusion reactions and NI.     

Identification of novel SNPs in the Sarda breed goats POU1F1 gene and their association with milk productive performance

The aim of the study was to detect polymorphism in the POU1F1 gene in Sarda breed goat, as well as to establish if SNPs could be associated with milk productive traits. The research was conducted on 129 Sarda breed goats from 4 to 5 years old, multiparous, lactating and in their third to fifth lactation. We report nine exonic and seven non-coding regions SNPs within the Sarda goat POU1F1 gene, namely, Ex 1 61 G>C; Ex 1 108 G>A; Ex 3 C>T; Ex 3 92 C>T; Ex 4 110 A>G; Ex 5 34 G>A resulting in Arg213Lys change; IVS4 641 G>A, IVS4 643 A>C, IVS4 659 G>A, IVS4 677 A>C, IVS4 G699Del, IVS4 709 C>G, Ex 6 17 G>A resulting in Arg228Ser change, Ex 6 58 G>T, Ex 6 172 T>C, 3′UTR 110 T>C. A statistically significant association was found between genotype TT, in position 17 of the exon 6 (3.1 % of frequency), and increased milk yield (P < 0.01) while genotype GT (25.6 % of frequency) was associated with a higher fat content. Genotype TT in position 58 of the exon 6 (3.9 % of frequency) was found to be associated with a higher fat (P < 0.01) and protein content (P < 0.05). Twenty-eight haplotypes were detected, but no significant association between the haplotypes and the milk production traits have been found. Our data, as well as providing new SNPs extending the POU1F1 gene characterization, evidence a relationship between polymorphism and milk production traits in Sarda goat breed.     

Significant association between SNPs in the superoxide dismutase 3, extracellular (SOD3) gene and resistance to Aeromonas hydrophila in the freshwater mussel Hyriopsis cumingii

Extracellular superoxide dismutase (SOD3) is a major antioxidant enzyme that protects organs from damage by reactive oxygen species (ROS). In this study, the SOD3 gene was identified and characterized from the freshwater mussel Hyriopsis cumingii (Hc‐SOD3). The cDNA sequence consists of 763 bp, encoding a protein of 208 amino acids. The amino acid sequence possesses two CuZnSOD signature sequences, and amino acids required for binding of Cu (His‐93, ‐95, ‐110 and ‐169) and Zn (His‐110, ‐118, ‐129 and Asp‐132) were conserved in Hc‐SOD3. The Hc‐SOD3 genomic sequence was 9165 bp in length, containing four exons and three introns. Eighteen single nucleotide polymorphisms were detected in the Hc‐SOD3 gene from resistant stock (RS) and susceptible stock (SS) of H. cumingii to Aeromonas hydrophila. The genotype and allele distribution were examined in resistant and susceptible stocks. Among them, a C/G substitution at the g.7994C>G locus and G/C substitution at the g.8087G>C locus were significantly associated with resistance/susceptibility of H. cumingii to A. hydrophila, both in genotype (P = 0.017, P = 0.004 respectively) and allele frequency (P = 0.021, P = 0.006 respectively). Linkage disequilibrium analysis revealed that g.7994C>G, g.8001A>G, g.8035G>A, g.8087G>C and g.8191T>A were in linkage disequilibrium. The results suggest that the two polymorphic loci, g.7994C>G and g.8087G>C, could be potential genetic markers for future molecular selection of strains that are resistant to diseases.