Monitoring changes in the genetic structure of Brown Tsaiya duck selected for feeding efficiency by microsatellite markers

Animal care

The animals, management methods, and experimental protocols used in this study were approved by the Institutional Animal Care and Use Committee of the Ilan Branch of the Livestock Research Institute, Council of Agriculture, Executive Yuan, Taiwan (case numbers LRI IACUC 104-012, LRIIL IACUC107004, LRIIL IACUC108003, and LRIIL IACUC 109010).

Sample collection and DNA extraction

The experimental animals were from 10 populations of the second, fourth, sixth, seventh, and eighth generations of the selected and control lines. The two lines’ genotypes of the second and fourth generations were obtained from the previous study (i.e., the RFC selected line and control line of the 18th and 20th generations of Brown Tsaiya LRI 1) [13], and the other generations were genotyped in this study. To increase the representativeness of the samples, we selected two males and two females from each sire family of the two lines. Table 1 presents the number of samples in each population.

Table 1

The number of samples in each generation of RFC selected line and the control line used in this study

Microsatellite genotyping

For each duck in this study, 2 to 3 mL of blood were drawn from the superficial plantar metatarsal vein or wing vein. The genomic DNA of the ducks was extracted from the blood samples by removing plasma with a Gentra Puregene Blood Kit (Qiagen, Venlo, The Netherlands) under the instructions; we followed Bush et al [15] to adjust the amount of blood samples.

The microsatellite markers were screened from Brown Tsaiya ducks [16]. After testing several Tsaiya duck populations, we selected those with superior polymorphisms (Table 2), namely APT001, APT004, APT008, APT010, APT012, APT017, APT020, APT025, APT026, APT032, and APT033, for a total of 11 markers. For the microsatellite analysis, polymerase chain reaction (PCR) was performed using a Veriti 96 Well Fast Thermal Cycler (Applied Biosystems by Thermo Fisher Scientific, Waltham, MA, USA). A final volume of 10 μL solution containing 0.375 U of Taq DNA polymerase (TaKaRa, Kyoto, Japan), 1×PCR buffer (1.5 mM MgCl₂), 0.2 μM deoxyribose nucleotide triphosphate (dNTP), 0.2 μM forward and reverse primers with fluorescent-labeled in the forward primer, and 50 ng of genomic DNA. The PCR cycling program was as follows: 94°C for 10 min, 30 cycles at 94°C for 20 s, 60°C for 30 s, 72°C for 30 s, and final elongation at 72°C for 10 min. The PCR products were analyzed through gel electrophoresis in 1% agarose gels, visualized using ethidium bromide staining, and pooled based on fragment size and fluorescent color for subsequent capillary electrophoresis analysis. The amplified microsatellite PCR products were denatured using Hi-Di formamide and analyzed using a DNA analyzer (ABI PRISM 3730 DNA analyzer; Applied Biosystems by Thermo Fisher Scientific, USA). Fragment size was estimated using GeneScan 500 LIZ Size Standard (Applied Biosystems by Thermo Fisher Scientific, USA). The resulting fragment analysis data and size of the alleles were analyzed using Peak Scanner Software version 1.0 (Applied Biosystems by Thermo Fisher Scientific, USA).

Table 2

Primer sequences, annealing temperature, and orthologous microsatellites in the duck genome scaffold of the 11 microsatellite markers derived from Tsaiya ducks

Statistical analysis

The genetic variation data, which consist of the number of alleles (N_a), observed (H_O), and expected (H_E) heterozygosity [17], and the polymorphic information content (PIC) [18], were calculated using Cervus (version 3.0) [19]. POPGENE 1.32 [20] was used to calculate the number of effective alleles (N_e), and estimate Wright’s F-statistics for each locus, namely within-population inbreeding coefficient (F_IS), among-population genetic differentiation (F_ST), and population inbreeding coefficient (F_IT), and to perform the Hardy–Weinberg equilibrium test [21] and Ewens-Watterson neutrality test to evaluate the selective neutrality of the 11 microsatellite markers [22]. We also used FreeNA [23] to detect the null allele of the 11 microsatellite markers. Effective population size (N_e) was estimated using the Jordan–Ryman temporal method and plan I in NeEstimator (V2) [24] with consideration for the mating system and sampling. The temporal change in allele frequency was used to calculate the mean Ne at different sampling time points [25]. FSTAT (version 2.9.3) [26] was used to estimate the allele richness (AR) and pairwise F_ST between populations as described in Weir and Cockerham [27] to evaluate the degree of differentiation between populations. The analysis of molecular variance (AMOVA) and calculating the number of private alleles per locus (N_p) were performed using GenAlEx (v6.503) [28] to divide the total genetic variance into population- and individual-based components. These parameters were used to investigate the markers’ polymorphisms and the change of genetic variation in the selected and control lines through the generations.

In the clustering analysis, the adegenet R package performed principal coordinate analysis (PCA) and plotting [29]. The clustering of individuals from multilocus genotypes was performed using the STRUCTURE (version 2.3) [30] admixture model. The number of clusters (K) was between 2 and 10. For each K, 100 independent runs were performed with 20,000 burn-in periods followed by 50,000 iterations. The optimal K was determined using STRUCTURE HARVESTER (v0.6.91) [31] following the Evanno method [32], and the results from STRUCTURE were integrated into CLUMPAK (version 1.1) [33] to visualize the long-term effects of selection on population differentiation.

Microsatellite marker polymorphism and F-statistics

Table 3 presents the polymorphisms and F-statistics of the 11 Tsaiya duck microsatellite markers in the selected and control lines. Fifty-four alleles were observed, and all 11 loci were polymorphic. The number of alleles (N_a) ranged from 2 to 8, with an average of 4.9 alleles per locus. The average number of effective alleles (N_e) per locus was 2.900 ranged from 1.911 (APT033) to 4.106 (APT004). The observed heterozygosity (H_O) among 11 microsatellite loci had a range from 0.227 (APT033) to 0.758 (APT026), with the average of 0.531. The expected heterozygosity (H_E) ranged from 0.477 (APT033) to 0.757 (APT004), with an average of 0.627. The PIC ranged from 0.363 (APT033) to 0.720 (APT026), with an average of 0.563. Among the 11 microsatellite markers, seven markers were highly polymorphic (PIC>0.5), and four markers were moderately polymorphic (0.5>PIC>0.25). Nevertheless, the Ewens-Watterson test indicates APT008 and APT026 are not selectively neutral (Supplementary Table S2). In addition, APT001, APT004, APT008, APT010, and APT017 deviated from the Hardy–Weinberg equilibrium (p<0.01). As for the F-statistics, the F_IS of the 11 microsatellite markers ranged from −0.097 (APT026) to 0.502 (APT033), with an average of 0.102. The F_IT ranged from −0.064 (APT026) to 0.538 (APT033), with an average of 0.154. The F_ST ranged from 0.021 (APT032) to 0.133 (APT010), with an average of 0.059. The results indicated that 5.9% of the overall genetic variation was due to differences between populations, whereas differences among individuals within each population caused 94.1% of total genetic variation.

Table 3

Characteristics of 11 Brown Tsaiya microsatellite markers used in five generations of RFC selected line and the control line

Inter- and intra-line comparison of genetic variation

Table 4 presents the average genetic variation of the 11 Tsaiya microsatellite markers in the various generations of the selected line (S) and the control line (C).

Table 4

Genetic variation analysis of the second, fourth, sixth, seventh, and eighth generations of RFC selected line and the control line with 11 Brown Tsaiya microsatellites markers

The inter-line comparison showed that N_a was higher in most of the generations of the selected line than in the control line; only S2 and S8 had slightly lower N_a values than did C2 and C8. N_e, H_O, H_E, and PIC in the control line were higher than those in the selected line, except for the eighth generation. The mean AR values estimated based on the minimum sample size 17 ranged from 3.189 (S8) to 3.727 in selected line, while 3.545 (C8) to 3.777 (C6) in control line. The selected line is slightly lower than the control line in each generation, especially in the eighth generation. The number of private alleles (N_p) ranged from 0.000 (S2) to 0.364 (S7) in selected line, and 0.000 (C4, C6, C7, and C8) to 0.091 (C2). Compared to the selected line, there was almost no private allele in the control line. However, except for Np no significant difference in genetic variation was observed between the lines (p>0.05).

Some markers deviated from the Hardy–Weinberg equilibrium in the sixth to eighth generations, namely APT001 and APT033 in the sixth generation of both lines, APT001 in S7, APT008 and APT033 in S8, and APT001 and APT008 in C8. F_IS was positive in all generations of both lines, but no difference between the lines was observed.

The intra–line comparison showed that the selected line, N_a, N_e, H_O, H_E, and PIC slightly increased then decreased with each generation, while those continued decreasing slight in the control line. The change in AR showed a downward trend over time. F_IS ranged from 0.082 (S2) to 0.137 (S7) in the selected line and from 0.048 (C2) to 0.165 (C4) in the control line, and no consistent increase was observed in each line. In addition, no significant differences in genetic variation were observed between generations in each line (p>0.05).

The effective population size (Ne) in the second and eighth generations, estimated through the Jorde–Ryman temporal method, was 27.2 for the selected line (95% jackknife confidence interval: 17.3 to 63.4) and 23.9 for the control line (95% jackknife confidence interval: 16.6 to 42.3); Ne was lower in the control line than in the selected line. We also estimated Ne in two generations of each line; the closer the two generations were, the higher the Ne was. For example, the upper limits of the confidence intervals in the seventh and eighth generations of the two lines were infinity (Supplementary Table S3).

The result of FreeNA [23] (Supplementary Table S4) showed that the null allele frequencies of the APT001 and APT033 in 3 and 1 populations were greater than 0.2, respectively, and the null allele frequencies of the two markers in 3 and 6 populations were between 0.15 and 0.2. The null allele frequencies of APT008 in the fourth and eighth generations of both lines were higher than 0.2.

Population differentiation and genetic structure analysis

Table 5 presents a pairwise comparison of F_ST between generations of the selected line and its control line. The intra-line comparison revealed low F_ST values in the second and sixth generations: 0.0307 and 0.0217, respectively. Moderate F_ST values were observed in the fourth, seventh, and eighth generations: 0.0540, 0.1123, and 0.1271, respectively. F_ST gradually increased from the second to fourth generations, decreased in the sixth generation, and increased again to 0.1123 in the seventh generation and to 0.1271 in the eighth generation. Overall, differentiation tended to increase with each generation.

Table 5

F_ST values based on 11 Brown Tsaiya microsatellite markers for the second, fourth, sixth, seventh, and eighth generations of the RFC selected line and its control line

The intra-line comparison in the selected line showed that S8 was significantly differentiated from S2, S4, and S6, and the F_ST values range from 0.0074 to 0.0510. However, only moderate differentiation was observed between S2 and S8 (F_ST = 0.051). Significant differentiation was observed between C2 and C7 to C8 and between C8 and C4 to C6, with F_ST values ranging from 0.0240 to 0.0606. Only moderate differentiation was observed between C2 and C8 (F_ST = 0.0606). The results indicate that F_ST increased as the generation gap widened for both the selected and control lines.

The AMOVA results in Table 6 indicate that largest variation in the second and eighth generations were within individuals, accounting for 92% and 76%, respectively, followed by variation among individuals, accounting for 5% and 12%. The variation among populations accounted for 3% and 12%, respectively. The results also indicate that the proportion of genetic variation among individuals and among populations increased between the second and eighth generations and that the eighth generation had a higher proportion of genetic variation than the second.

Table 6

Summary of AMOVA results for second and eighth generations of the RFC selected line and its control line

Population PCA was performed using the adegenet R package (Supplementary Figure S1). The percentage of variation explained by the first, second, and third axes was 8.5%, 5.1%, and 4.4%, respectively. Because we focused on the cross-generation analysis of two lines, the proportion explained by PCA was low. Because of the low defined proportion, we plotted the first and second principal coordinates. Although the selected line and its control line were distributed into two groups in the PCA, numerous overlapping sections were observed. No significant differences in distribution range were observed among generations, indicating that the population structure was similar and that polymorphisms did not decrease.

Figure 1 presents the results of the STRUCTURE analysis visualized using CLUMPAK 1.1. The optimal K was 2, estimated through the Evanno method [32]. When K was 2, although the genetic structure of the second and fourth generations slightly differed, mixed clusters were observed. After the seventh generation, the two lines formed distinct clusters. When K was 3, although the populations formed mixed clusters, the genetic structure differed slightly between C2 and C7 and between S2 and S7. In the eighth generation, the difference from the second generation is more obvious than the seventh generation. The STRUCTURE results indicate that the differentiation between the selected and control lines increased with each generation, similar to the F_ST and AMOVA results.

Figure 1

STRUCTURE clustering of the residual feed consumption (RFC) selected and control lines. Numbers in parentheses indicate the identical solutions of the 100 runs at the 95% threshold, and colors correspond to genetic clusters. Each vertical bar represents an individual. S, RFC selected line; C, the control line.

Microsatellite marker polymorphism

So far, FAO has not offered recommended microsatellite markers for ducks, so this study used the 11 microsatellite markers from Brown Tsaiya ducks [16] with superior polymorphisms. The average N_e of Korean and Bangladeshi ducks is 2.8 to 3.7 [34], and that of local duck breeds in India is 2 to 2.9 [5]. In this study, the average N_e of each generation of the two lines was 2.4 to 3.0, which is not significantly different from those observed in other studies.

To confirm the polymorphism of the markers used in this study, we also tested the 11 microsatellite markers in Pekin duck, Kaiya, and several lines of Brown and White Tsaiya ducks to investigate polymorphisms in other breeds. The mean N_a per locus was 7.3, and the mean N_e per locus was 3.6. The mean H_O was 0.474, the mean H_E was 0.690, and the mean PIC was 0.646. All markers were highly polymorphic, except for APT001 and APT032, which were moderately polymorphic [35]. Compared to our study, the study with the newly screened markers from Brown Tsaiya duck, there was an average N_a was 11.3 per locus, the mean N_e was 5.4, the mean H_O was 0.591, the mean H_E was 0.747, and the mean PIC was 0.708 [14]; the genetic variation was slightly higher than in this study.

Through the Ewens-Watterson test, we found that APT008 and APT026 were non-neutral markers. However, Pampoulie et al [36] pointed out that both neutral and non-neutral genetic markers gave consistent results, the overall level of genetic differentiation was similar whether or not the non-neutral marker was removed, even combined with non-neutral markers can provide better understand the potential effects of selection [36]. In this study, the purpose is to investigate the genetic variation and genetic structure of RFC selected and control lines, so the non-neutral markers were not excluded.

All markers used in this study were four-nucleotide repeats to prevent genotyping errors than the dinucleotide and trinucleotide repeats. The duck genome has been gradually revealed since 2020. The microsatellite libraries of the 11 markers were queried in the RefSeq duck (A. platyrhynchos) genome database (GCF_015476345.1_ZJU1.0) through BLAST to find possible chromosomal locations. The 11 markers were distributed on chromosomes 1, 2, 3, 7, and Z (Table 2). However, all 11 markers were on different duck chromosome scaffolds when we queried them in the database’s previous version (GCF_000355885.1), suggesting a considerable distance between them. Therefore, we considered that the 11 Tsaiya microsatellite markers had a substantial-resolution in the genetic analysis.

However, in five microsatellite markers, the number of alleles identified is equal or less than four. The low N_a may affect the estimation of differentiation level, diversities and heterozygosity. For example, markers with high N_a are more likely to show differences in allele frequencies among different population. Taking biallelic single nucleotide polymorphisms (SNPs) as an example, it is necessary to increase the number of markers by 2 to 11 times to achieve the same resolution of microsatellite markers [37]. Therefore, we will add highly polymorphic markers from other studies and consider the repeat motifs and positions on chromosomes in the future.

Inter- and intra-line comparison of genetic variations

The results showed there were slightly decrease in the genetic variations, especially about the heterozygosity in the control line. Among them, the genetic variations of the second and fourth generations are slightly lower than those of the sixth generation. Hale et al [38] pointed out that when the sample size in a single population was less than 25, some alleles may not be detected, resulting in low N_a, and the allele frequencies and H_E estimated by samples may be quite different from the real population. However, the sample size of 25 to 30 is representative, and it is little benefit of further increasing the sample size per population. Since the sample size in the second and fourth generations of this study is less than 30, the genetic variation may be underestimated [38]. For subsequent genetic monitoring, at least 30 individuals per population should be sampled.

And there was also increase in F_IS in the both line, although no consistent increase in F_IS was observed in each line. It indicated there were heterozygosity deficits in both the selected and control lines. Studies have indicated that possible causes of heterozygosity deficits include selection, inbreeding, geographic isolation, sampling error, the Wahlund effect, and null alleles [39].

The F_IS of each population was similar to the one of the total population (Tables 3 and 4). The ratio of the markers deviated from the Hardy–Weinberg equilibrium after the sixth generation in the two lines were 9% to 18% (1/11 to 2/11), which were less than the one in Lai et al [14] (7/24, 29%). The study also indicated H_O is less than H_E in Taiwan’s duck breeds and lines, even in the germplasm-conserved Brown Tsaiya duck and White Tsaiya duck with rotational mating schemes and without selection. It’s just the degree of heterozygosity deficit in these two conserved populations was lower than that in other selected populations. Also one study genetically monitored five Italian populations with in situ conservation and a circular mating scheme for four years. They discovered that F_IS decreased in all but one slightly selected population, in which F_IS increased significantly [40]. In our study, the selected line has been selected for RFC for more than ten years. Although the control line has not been selected for any additional trait, its ancestor, Brown Tsaiya LRI 1, has been selected for egg-laying performance for a long time. So selection should be one of the reasons for heterozygosity deficit in this study.

From the zero to eighth generations, the RFC selected line in each generation comprised 65 to 86 drakes, with a selection rate of 14% to 28.6%, and 104 to 152 female ducks, with a selection rate of 20.8% to 82.6%. The control line comprised 64 to 84 drakes and 89 to 144 female duck, with a selection rate of 14.1% to 18.8% and 27.8% to 51.7%, respectively (Supplementary Table S1). Only the zero generation was not divided into two lines, and the third generation of the selected line may have been affected by the quality of the feed, resulting in a smaller population than other generations [13]. Because the original populations sizes were small and genetic drift, low-frequency alleles were unlikely to be observed (for example, the 201 bp allele of APT012 only appeared in C2) and this increased the differences in H_O, H_E, and F_IS and decreased heter-ozygosity.

Although there were increase in F_IS in the both line, however, the F_IS was not significantly different from 0 and did not increase sharply. This may have been due to attempting to prevent full- and half-sib mating in both lines. The FreeNA results (Supplementary Table S4) showed that the frequencies of null allele in APT001 and APT033 were higher or close to 0.2 in more than half of the populations. The results of BLAST to RefSeq indicated the two markers located on sex chromosomes, and may result in PCR failure in one of an allele. We think it may cause part of the null allele. After the markers with null allele frequencies higher than 0.2 were excluded, except for S6 (no markers have been removed), the average F_IS of each population of the selected line and its control line was less than 0.1.

We also pay attention to the changes of allele, the number of private alleles mainly appeared in the selected line, which corresponds to Petit et al [41] consider private allele to contribute to response to selection or have evolutionary significance. And it may also make us to see higher N_a on the selected line. In the other hand, the AR in our study showed a downward trend over time, especially in the S8. AR is considered sensitive to founder events, compared to heterozygosity, AR decrease while the loss of alleles during founder events [42]. Therefore, the reason for the decline in AR should be mainly related to RFC selection.

Finally, to understand the status of genetic diversity and conservation management, we estimated the effective population size of the selected and control line. Few temporal sampling studies have been conducted, and most studies on Ne have used the linkage disequilibrium (LD) method, which is not strongly affected by genetic drift, resulting in higher Ne values. The LD method is limited by random mating and lifetime monogamy [24], which are inapplicable to the population in this experiment. Therefore, we only used a temporal method for estimation. The estimated Ne values of both lines in this experiment were 27.2 in the selected line and 23.9 in the control line. The Ne in the control (23.9) line was lower than the selected line (27.2) may be result from lower actual population size since the sixth generation (Supplementary Table S1). However, the Ne in both line were higher than Hungarian chicken breed study using a temporal method (4.4 to 12.7) [11]. It is similar to Speckled Hungarian in the same study and a study on indigenous Chinese chickens estimated by genome-wide SNPs using the LD method [43].

In summary, the heterozygosity deficit and change in F_IS each generation may have been caused by selection for RFC, genetic drift, and null alleles. Also, RFC selection may result in the decline in AR. However, the consistent F_IS over time, N_e, and the performance monitoring results indicate no severe inbreeding depression. The follow-up genetic monitoring should be performed.

Population differentiation and genetic structure analysis

According to Wright [44], an F_ST lower than 0.05 indicates nearly no differentiation, 0.05 to 0.15 indicates moderate differentiation, 0.15 to 0.25 indicates high differentiation, and values higher than 0.25 indicate extremely highly differentiation. Except for the sixth generation, the second through eighth generations was increasingly differentiated between the two lines (Table 5). As for inter-line pairwise F_ST, the value between S6 and C6 indicated a decrease, followed by a gradual increase in the seventh and eighth generations. When we reproduced the fifth generation from the fourth generation, one of the drakes in the control line was mistakenly replaced with one from the selected line. The individuals with problematic pedigrees were gradually culled from the seventh and eighth generations. The temporary introduction of the selected line individual to the control line may have decreased the pairwise F_ST between S6 and C6 considerably, and have also caused the fluctuation in the F_IS. When we culled the problematic individuals from the seventh and eighth generations, F_ST increased as the gap between generations widened, which can be attributed to the RFC selection for the selected line; the control line was not selected for RFC.

Although C6 was affected by the misuse of the selected line drakes, the increase in pairwise F_ST between generations in the control line may have been caused by the smaller population size and genetic drift. Similar to the genetic structure study of local chicken breeds in Hungary in 2002 and 2017 [11], the pairwise F_ST of the same breed at different times ranged from 0.03 to 0.08; this indicates slight differentiation. Because the populations in a Hungarian study were small, approximately 200 individuals, the differentiation may have also been caused by genetic drift.

The pairwise F_ST values of the seventh and eighth generations between the selected and control lines (i.e., Brown Tsaiya LRI 1) in this study were 0.1123 and 0.1271, respectively. Similar to the genetic analysis of Taiwan’s duck breeds and lines [14] indicated that the range of F_ST between Brown Tsaiya duck lines (preserved Brown Tsaiya duck, Brown Tsaiya LRI 1, LRI 2, and LRI 3) separated by more than ten generations was 0.092 to 0.167. The result indicated that the differentiation between the selected and control lines has reached the level of different lines. And the pairwise F_ST between lines from 0.0307 (C2 and S2) increased to 0.1271 (C8 and S8), the result is similar to Hungarian study. The pairwise F_ST between breeds ranged from 0.15 to 0.28 in 2002 and slightly increased to 0.21 to 0.30 after 15 years, the differentiation increased, and genetic similarity between breeds decreased over time [11].

F_IT (0.154) was mainly affected by F_IS (0.102) rather than by F_ST (Table 3). The AMOVA results for the second and eighth generations (Table 6) indicate that the primary genetic variation was between individuals. The variance ratio in the eighth generation is similar to that of other study on Brown Tsaiya ducks [14]. However, because the study included several lines of Brown Tsaiya ducks, the variance among the populations was slightly higher than that in this study. In this study, the variance among populations increased from 3% to 12% between the second and eighth generations, and the variance among individuals increased from 5% to 12%. The results are similar to those of the Hungarian study; the variance among populations increased from 23% to 25%, and the variance among individuals increased from 1% to 6% between 2002 and 2017 [11]. Therefore, the change in AMOVA results may also have been caused by genetic drift or selection.

As for PCA (Supplementary Figure S1), the first and second axes do not explain a considerable portion of the variation, and figure 1 indicates that each population did not form independent clusters. Our results are similar to those of Palinkas-Bodzsar et al [11], the distributions of the same breed sampled at an interval of 15 years were still very close, even though the first and second coordinates explained 9.13% and 17.47% of variations, higher than our study.

The STRUCTURE results are consistent with the pairwise F_ST results. When K was 2, the two lines formed independent clusters by the seventh generation. However, the proportions of inferred clusters on the left one-third of the samples in C6 were highly similar to the S6, and those in the other two-thirds of the samples were similar to C4. This may have been caused by the temporary introduced the selected line to the control line. Overall, the differentiation between the selected and control lines increased with each generation.

With the strengthening effects of climate change on the livestock industry, promoting “Better Feed Efficiency Brown Tsaiya” can help farmers reduce feed costs by 10% and mitigate the negative effects of extreme weather. However, maintaining large breeding populations is difficult with a limited workforce and management costs. Smaller poultry populations also often result in changes in allele frequency. Therefore, changes in heterozygosity over time must be monitored. Few studies have genetically monitored chickens over time, and no such study has been conducted on ducks. To ensure the sustainable management of key duck breeds, we used microsatellite markers to monitor ducks over time genetically. The results indicate that the markers reflect changes in heterozygosity and differentiation between the selected line and its control line. Although selection, genetic drift, and null alleles may have slightly decreased the heterozygosity and may result in the loss of allele in the two lines, heterozygosity remained, and F_IS did not increase considerably, indicating that the breeding and population management strategies were effective. To preserve genetic heterozygosity, rotational mating schemes, or increasing the number of drake families under the same population can be used if necessary. The results of this study indicate that genetic variation and structure did not change considerably after six generations with the current breeding and mating strategy and can be used as a reference for other poultry conservation research. Genetic monitoring should be performed regularly to conserve the germplasm, but the intervals between analyses can be longer to reduce costs.