Structural variations (SVs) are a major source of domestication and improvement traits. We present the first duck pan‐genome constructed using five genome assemblies capturing ∼40.98 Mb new sequences. This pan‐genome together with high‐depth sequencing data (∼46.5×) identified 101,041 SVs, of which substantial proportions were derived from transposable element (TE) activity. Many TE‐derived SVs anchoring in a gene body or regulatory region are linked to duck's domestication and improvement. By combining quantitative genetics with molecular experiments, we, for the first time, unraveled a 6945 bp Gypsy insertion as a functional mutation of the major gene IGF2BP1 associated with duck bodyweight. This Gypsy insertion, to our knowledge, explains the largest effect on bodyweight among avian species (27.61% of phenotypic variation). In addition, we also examined another 6634 bp Gypsy insertion in MITF intron, which triggers a novel transcript of MITF, thereby contributing to the development of white plumage. Our findings highlight the importance of using a pan‐genome as a reference in genomics studies and illuminate the impact of transposons in trait formation and livestock breeding.

Several geographically isolated populations of Tibetan pigs inhabit the high-altitude environment of the Tibetan Plateau. Their genetic relationships, contribution to the pool of genetic diversity, and their origin of domestication are unclear. In this study, whole-genome re-sequencing data from 10 geographically isolated Tibetan pig populations were collected and analyzed. Population genetic analyses revealed limited genetic differentiation among the Tibetan pig populations. Evidence from deleterious variant analysis indicated that population-specific deleterious variants were the major component of all mutational loci. Contribution to the meta-population was largest in the TT (Qinghai-Tibet Plateau) population, based on gene diversity or allelic diversity. Selective sweep analysis revealed numerous genes, including RXFP1, FZD1, OR1F1, TBX19, MSTN, ESR1, MC1R, HIF3A, and EGLN2 which are involved in lung development, hard palate development, coat color, hormone metabolism, facial appearance, and perception of smell. These findings increase our understanding of the origins and domestication of the Tibetan pig, and help optimize the strategy for their conservation.

Tibetan pigs from the Tibetan Plateau are characterized with a significant phenotypic difference relative to lowland pigs. In this study, a significant difference of the fatness and fatty acid composition traits was observed between the Tibetan and Yorkshire pigs. To uncover the involved mechanism, the expression profile of long noncoding RNAs (lncRNAs) and genes was compared between them. After serial filtered steps, 1,964 lncRNAs were obtained through our computational pipeline. In total, 63 and 715 lncRNAs and genes were identified to be differentially expressed. Evidence from cis- and trans-targeting analysis of lncRNAs demonstrated that some lncRNAs, such as MSTRG.14097 and MSTRG.8034, played important roles in the fatness and fatty acid composition traits. Bioinformatics analysis revealed that many candidate genes were responsible for the two traits. Of these, FASN, ACACA, SCD, ME3, PDHB, ACSS1, ACSS2, and ACLY were identified, which functioned in regulating the level of hexadecanoic acid, hexadecenoic acid, octadecenoic acid, and monounsaturated fatty acid. And LPGAT1, PDK4, ACAA1, and ADIPOQ were associated with the content of stearic acid, octadecadienoic acid, and polyunsaturated fatty acid. Candidate genes, which were responsible for fatness trait, consisted of FGF2, PLAG1, ADIPOQ, IRX3, MIF, IL-34, ADAM8, HMOX1, Vav1, and TLR8. In addition, association analysis also revealed that 34 and 57 genes significantly correlated to the fatness and fatty acid composition trait, respectively. Working out the mechanism caused by these lncRNAs and candidate genes is proven to be complicated but is invaluable to our understanding of fatness and fatty acid composition traits.

Litter size has increased and farrowing duration has also prolonged in recent years. The aim of this study was to analyze the effect of litter size and parity on farrowing duration (FAR) to estimate the possibility of selecting a short farrowing duration. We recorded 32,200 parturitions of 8420 Landrace × Yorkshire sows, determined farrowing duration, litter size, parity, gestation length. Results showed that total number of born (TNB) and parity obeyed a cubic (p = 0.0004, p = 0.004) relationship while number born alive (NBA) and number born dead (NBD) obeyed a linear (p = 0.0239, p = 0.0035) relationship with FAR. Gestation length obeyed a linear (p = 0.02) relationship with FAR. FAR of sows with stillbirth was longer than that of sows without stillbirth. Stillbirth rate increased rapidly from about 2% to 4%, especially when FAR was over 240 min. FAR gradually prolonged with the parities. FAR of 7th parity sows was longer than that of 1st~6th parity sows (p < 0.05), but different parity sows had little difference in the same FAR interval except for gilts. Results indicated it was possible and necessary to consider FAR into pig breeding without worrying about decreasing of live litter size or negative effect of parity if FAR was shorter than 300 min.

To achieve sustainable development of the poultry industry, the effective conservation of genetic resources has become increasingly important. In the present study, we systematically elucidated the population structure, conservation priority, and runs of homozygosity (ROH) patterns of Chinese native chicken breeds. We used a high-density genotyping dataset of 157 native chickens from eight breeds. The population structure showed different degrees of population stratification among the breeds. Chahua chicken was the most differentiated breed from the other breeds (Nei = 0.0813), and the Wannan three-yellow chicken (WanTy) showed the lowest degree of differentiation (Nei = 0.0438). On the basis of contribution priority, Xiaoshan chicken had the highest contribution to the total gene diversity (1.41%) and the maximum gene diversity of the synthetic population (31.1%). WanTy chicken showed the highest contribution to the total allelic diversity (1.31%) and the maximum allelic diversity of the syntenic population (17.0%). A total of 5242 ROH fragments and 5 ROH island regions were detected. The longest ROH fragment was 41.51 Mb. A comparison of the overlapping genomic regions between the ROH islands and QTLs in the quantitative trait loci (QTL) database showed that the annotated candidate genes were involved in crucial economic traits such as immunity, carcass weight, drumstick and leg muscle development, egg quality and egg production, abdominal fat precipitation, body weight, and feed intake. In conclusion, our findings revealed that Chahua, Xiaoshan, and WanTy should be the priority conservation breeds, which will help optimize the conservation and breeding programs for Chinese indigenous chicken breeds.