Selection Signature Analysis Research Articles

Assembly and annotation of the black spruce genome provide insights on spruce phylogeny and evolution of stress response.

Black spruce (Picea mariana [Mill.] B.S.P.) is a dominant conifer species in the North American boreal forest that plays important ecological and economic roles. Here, we present the first genome assembly of P. mariana with a reconstructed genome size of 18.3 Gbp and NG50 scaffold length of 36.0 kbp. A total of 66,332 protein-coding sequences were predicted in silico and annotated based on sequence homology. We analyzed the evolutionary relationships between P. mariana and 5 other spruces for which complete nuclear and organelle genome sequences were available. The phylogenetic tree estimated from mitochondrial genome sequences agrees with biogeography; specifically, P. mariana was strongly supported as a sister lineage to P. glauca and 3 other taxa found in western North America, followed by the European Picea abies. We obtained mixed topologies with weaker statistical support in phylogenetic trees estimated from nuclear and chloroplast genome sequences, indicative of ancient reticulate evolution affecting these 2 genomes. Clustering of protein-coding sequences from the 6 Picea taxa and 2 Pinus species resulted in 34,776 orthogroups, 560 of which appeared to be specific to P. mariana. Analysis of these specific orthogroups and dN/dS analysis of positive selection signatures for 497 single-copy orthogroups identified gene functions mostly related to plant development and stress response. The P. mariana genome assembly and annotation provides a valuable resource for forest genetics research and applications in this broadly distributed species, especially in relation to climate adaptation.

Comprehensive selection signature analyses in dairy cattle exploiting purebred and crossbred genomic data.

The main objective of the current research was to locate, annotate, and highlight specific areas of the bovine genome that are undergoing intense positive selection. Here, we are analyzing selection signatures in crossbred (Bos taurus X Bos indicus), taurine (Bos taurus), and indicine (Bos indicus) cattle breeds. Indicine cattle breeds found throughout India are known for their higher heat tolerance and disease resilience. More breeds and more methods can provide a better understanding of the selection signature. So, we have worked on nine distinct cattle breeds utilizing seven different summary statistics, which is a fairly extensive approach. In this study, we carried out a thorough genome-wide investigation of selection signatures using bovine50KSNP data. We have included the genotyped data of two taurine, two crossbreds, and five indicine cattle breeds, for a total of 320 animals. During the 1950s, these indicine (cebuine) cattle breeds were exported with the aim of enhancing the resilience of taurine breeds in Western countries. For this study, we employed seven summary statistics, including intra-population, i.e., Tajima's D, CLR, iHS, and ROH and inter-population statistics, i.e., FST, XP-EHH, and Rsb. The NCBI database, PANTHER 17.0, and CattleQTL database were used for annotation after finding the important areas under selection. Some genes, including EPHA6, CTNNA2, NPFFR2, HS6ST3, NPR3, KCNIP4, LIPK, SDCBP, CYP7A1, NSMAF, UBXN2B, UGDH, UBE2K, and DAB1, were shown to be shared by three or more different approaches. Therefore, it gives evidence of the most intense selection in these areas. These genes are mostly linked to milk production and adaptability traits. This study also reveals selection regions that contain genes which are crucial to numerous biological functions, including those associated with milkproduction, coat color, glucose metabolism, oxidative stress response, immunity and circadian rhythms.

Identification of SNPs in MITF associated with beak color of duck.

Introduction: Beak color-a pigment-related trait-is an important feature of duck breeds. Recently, little research has addressed genetic mechanism of the beak colors in poultry, whereas the process and the regulation factors of melanin deposition have been well described. Methods: To investigate the genetic mechanism of beak colors, we conducted an integrated analysis of genomic selection signatures to identify a candidate site associated with beak color. For this, we used black-billed (Yiyang I meat duck synthetic line H1, H2, H3&HF) and yellow-billed ducks (Cherry Valley ducks and white feather Putian black duck). Quantitative real-time PCR and genotyping approaches were used to verify the function of the candidate site. Results: We identified 3,895 windows containing 509 genes. After GO and KEGG enrichment analysis, nine genes were selected. Ultimately, MITF was selected by comparing the genomic differentiation (FST). After loci information selection, 41 extreme significantly different loci were selected, which are all located in intron regions of MITF and are in almost complete linkage disequilibrium. Subsequently, the site ASM874695v1:10:g.17814522T > A in MITF was selected as the marker site. Furthermore, we found that MITF expression is significantly higher in black-beaked ducks than in yellow-beaked ducks of the F2 generation (p < 0.01). After genotyping, most yellow-billed individuals are found with homozygous variant; at the same time, there are no birds with homozygous variant in black-billed populations, while the birds with homozygous and heterozygous variant share the same proportion. Conclusion: MITF plays a very critical role in the melanogenesis and melanin deposition of duck beaks, which can effectively affect the beak color. The MITF site, ASM874695v1:10:g.17814522T > A could be selected as a marker site for the duck beak color phenotype.

Landscape genomics analysis provides insights into future climate change-driven risk in rhesus macaque

Climate change significantly affects the suitability of wildlife habitats. Thus, understanding how animals adapt ecologically and genetically to climate change is important for targeted species protection. Rhesus macaques (Macaca mulatta) are widely distributed and multi-climatically adapted primates. This study explored how rhesus macaques adapt to climate change by integrating ecological and genetic methods and applying species distribution models (SDMs) and a gradient forest (GF) model. The findings suggested that temperature seasonality primarily affects habitat suitability and indicated that climate change will have a dramatic impact on macaque populations in the future. We also applied genotype-environment association (GEA) analyses and selection signature analyses to identify genes associated with climate change and provide possible explanations for the adaptation of rhesus macaques to climatic environments. The population genomics analyses suggested that the Taihang population has the highest genomic vulnerability with inbreeding and low heterozygosity. Combined with the higher ecological vulnerability, additional conservation strategies are required for this population under higher risk of climate change. Our work measured the impact of climate change and enabled the identification of populations that exhibit high vulnerability to severe climate change. Such information is useful for selecting populations of rhesus macaques as subject of long-term monitoring or evolutionary rescue under future climate change.

Exploring genome-wide differentiation and signatures of selection in Italian and North American Holstein populations

Among Italian dairy cattle, the Holstein is the most reared breed for the production of Parmigiano Reggiano protected designation of origin cheese, which represents one of the most renowned products in the entire Italian dairy industry. In this work, we used a medium-density genome-wide data set consisting of 79,464 imputed SNPs to study the genetic structure of Italian Holstein breed, including the population reared in the area of Parmigiano Reggiano cheese production, and assessing its distinctiveness from the North American population. Multidimensional scaling and ADMIXTURE approaches were used to explore the genetic structure among populations. We also investigated putative genomic regions under selection among these 3 populations by combining 4 different statistical methods based either on allele frequencies (single marker and window-based) or extended haplotype homozygosity (EHH; standardized log-ratio of integrated EHH and cross-population EHH). The genetic structure results allowed us to clearly distinguish the 3 Holstein populations; however, the most remarkable difference was observed between Italian and North American stock. Selection signature analyses identified several significant SNPs falling within or closer to genes with known roles in several traits such as milk quality, resistance to disease, and fertility. In particular, a total of 22 genes related to milk production have been identified using the 2 allele frequency approaches. Among these, a convergent signal has been found in the VPS8 gene which resulted to be involved in milk traits, whereas other genes (CYP7B1, KSR2, C4A, LIPE, DCDC1, GPR20, and ST3GAL1) resulted to be associated with quantitative trait loci related to milk yield and composition in terms of fat and protein percentage. In contrast, a total of 7 genomic regions were identified combining the results of standardized log-ratio of integrated EHH and cross-population EHH. In these regions candidate genes for milk traits were also identified. Moreover, this was also confirmed by the enrichment analyses in which we found that the majority of the significantly enriched quantitative trait loci were linked to milk traits, whereas the gene ontology and pathway enrichment analysis pointed to molecular functions and biological processes involved in AA transmembrane transport and methane metabolism pathway. This study provides information on the genetic structure of the examined populations, showing that they are distinguishable from each other. Furthermore, the selection signature analyses can be considered as a starting point for future studies in the identification of causal mutations and consequent implementation of more practical application.

Identifying Candidate Genes for Litter Size and Three Morphological Traits in Youzhou Dark Goats Based on Genome-Wide SNP Markers.

This study aimed to reveal the potential genetic basis for litter size, coat colour, black middorsal stripe and skin colour by combining genome-wide association analysis (GWAS) and selection signature analysis and ROH detection within the Youzhou dark (YZD) goat population (n = 206) using the Illumina GoatSNP54 BeadChip. In the GWAS, we identified one SNP (snp54094-scaffold824-899720) on chromosome 11 for litter size, two SNPs on chromosome 26 (snp11508-scaffold142-1990450, SORCS3) and chromosome 12 (snp55048-scaffold842-324525, LOC102187779) for coat colour and one SNP on chromosome 18 (snp56013-scaffold873-22716, TCF25) for the black middorsal stripe. In contrast, no SNPs were identified for skin colour. In selection signature analysis, 295 significant iHS genomic regions with a mean |iHS| score > 2.66, containing selection signatures encompassing 232 candidate genes were detected. In particular, 43 GO terms and one KEGG pathway were significantly enriched in the selected genes, which may contribute to the excellent environmental adaptability and characteristic trait formation during the domestication of YZD goats. In ROH detection, we identified 4446 ROH segments and 282 consensus ROH regions, among which nine common genes overlapped with those detected using the iHS method. Some known candidate genes for economic traits such as reproduction (TSHR, ANGPT4, CENPF, PIBF1, DACH1, DIS3, CHST1, COL4A1, PRKD1 and DNMT3B) and development and growth (TNPO2, IFT80, UCP2, UCP3, GHRHR, SIM1, CCM2L, CTNNA3 and CTNNA1) were revealed by iHS and ROH detection. Overall, this study is limited by the small population size, which affects the results of GWAS to a certain extent. Nevertheless, our findings could provide the first overview of the genetic mechanism underlying these important traits and provide novel insights into the future conservation and utilisation of Chinese goat germplasm resources.

Two mutations at KRT74 and EDAR synergistically drive the fine-wool production in Chinese sheep

IntroductionFine-wool sheep are the most common breed used by the wool industry worldwide. Fine-wool sheep have over a three-fold higher follicle density and a 50% smaller fiber diameter than coarse-wool sheep. ObjectivesThis study aims to clarify the underlying genetic basis for the denser and finer wool phenotype in fine-wool breeds. MethodWhole-genome sequences of 140 samples, Ovine HD630K SNP array data of 385 samples, including fine, semi-fine, and coarse wool sheep, as well as skin transcriptomes of nine samples were integrated for genomic selection signature analysis. ResultsTwo loci at keratin 74 (KRT74) and ectodysplasin receptor (EDAR) were revealed. Fine-scale analysis in 250 fine/semi-fine and 198 coarse wool sheep narrowed this association to one C/A missense variant of KRT74 (OAR3:133,486,008, P = 1.02E-67) and one T/C SNP in the regulatory region upstream of EDAR (OAR3:61,927,840, P = 2.50E-43). Cellular over-expression and ovine skin section staining assays confirmed that C-KRT74 activated the KRT74 protein and specifically enlarged cell size at the Huxley’s layer of the inner root sheath (P < 0.01). This structure enhancement shapes the growing hair shaft into the finer wool than the wild type. Luciferase assays validated that the C-to-T mutation upregulated EDAR mRNA expression via a newly created SOX2 binding site and potentially led to the formation of more hair placodes. ConclusionsTwo functional mutations driving finer and denser wool production were characterized and offered new targets for genetic breeding during wool sheep selection. This study not only provides a theoretical basis for future selection of fine wool sheep breeds but also contributes to improving the value of wool commodities.

The Genetic Diversity of Mink (Neovison vison) Populations in China.

The American mink (Neovison vison) is a semiaquatic species of Mustelid native to North America that is now widespread in China. However, the knowledge of genetic diversity of mink in China is still limited. In this study, we investigated the genetic diversity and identified significant single nucleotide polymorphisms (SNPs) in mink populations of five different color types in three different mink farms in China. Using double-digest restriction site-associated DNA sequencing, we identified a total of 1.3 million SNPs. After filtering the SNPs, phylogenetic tree, Fst, principal component, and population structure analyses were performed. The results demonstrated that red mink and black mink grouped, with separate clustering of all other color types. The population divergence index (Fst) study confirmed that different mink populations were distinct (K = 4). Two populations with different coat colors were subjected to the selection signature analysis, and 2300 genes were found to have a clear selection signature. The genes with a selection signature were subjected to Gene Ontology (GO) categorization and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis, the results revealed that the genes with a selection signature were enriched in the melanogenesis pathway. These study's findings have set the stage for improved breeding and conservation of genetic resources in real-world practical mink farming.

Genetic analysis of global faba bean diversity, agronomic traits and selection signatures

Key messageWe identified marker-trait associations for key faba bean agronomic traits and genomic signatures of selection within a global germplasm collection.Faba bean (Vicia faba L.) is a high-protein grain legume crop with great potential for sustainable protein production. However, little is known about the genetics underlying trait diversity. In this study, we used 21,345 high-quality SNP markers to genetically characterize 2678 faba bean genotypes. We performed genome-wide association studies of key agronomic traits using a seven-parent-MAGIC population and detected 238 significant marker-trait associations linked to 12 traits of agronomic importance. Sixty-five of these were stable across multiple environments. Using a non-redundant diversity panel of 685 accessions from 52 countries, we identified three subpopulations differentiated by geographical origin and 33 genomic regions subjected to strong diversifying selection between subpopulations. We found that SNP markers associated with the differentiation of northern and southern accessions explained a significant proportion of agronomic trait variance in the seven-parent-MAGIC population, suggesting that some of these traits were targets of selection during breeding. Our findings point to genomic regions associated with important agronomic traits and selection, facilitating faba bean genomics-based breeding.

Systematic Selection Signature Analysis of Chinese Gamecocks Based on Genomic and Transcriptomic Data.

Selection pressures driven by natural causes or human interference are key factors causing genome variants and signatures of selection in specific regions of the genome. Gamecocks were bred for cockfighting, presenting pea-combs, larger body sizes, stronger limbs, and higher levels of aggression than other chickens. In this study, we aimed to explore the genomic differences between Chinese gamecocks and commercial, indigenous, foreign, and cultivated breeds by detecting the regions or sites under natural or artificial selection using genome-wide association studies (GWAS), genome-wide selective sweeps based on the genetic differentiation index (FST), and transcriptome analyses. Ten genes were identified using GWAS and FST: gga-mir-6608-1, SOX5, DGKB, ISPD, IGF2BP1, AGMO, MEOX2, GIP, DLG5, and KCNMA1. The ten candidate genes were mainly associated with muscle and skeletal development, glucose metabolism, and the pea-comb phenotype. Enrichment analysis results showed that the differentially expressed genes between the Luxi (LX) gamecock and Rhode Island Red (RIR) chicken were mainly related to muscle development and neuroactive-related pathways. This study will help to understand the genetic basis and evolution of Chinese gamecocks and support the further use of gamecocks as an excellent breeding material from a genetic perspective.

Selection signature analysis reveals RDH5 performed key function in vision during sheep domestication process.

As one of the most successful domesticated animals in the Neolithic age, sheep gradually migrated all over the world with human activities. During the domestication process, remarkable changes have taken place in morphology, physiology, and behavior, resulting in different breeds with different characters via artificial and natural selection. However, the genetic background responsible for these phenotypic variations remains largely unclear. Here, we used whole genome resequencing technology to compare and analyze the genome differences between Asiatic mouflon wild sheep (Ovis orientalis) and Hu sheep (Ovis aries). A total of 755 genes were positively selected in the process of domestication and selection, and the genes related to sensory perception had directional evolution in the autosomal region, such as OPRL1, LEF1, TAS1R3, ATF6, VSX2, MYO1A, RDH5, and some novel genes. A missense mutation of c.T722C/p.M241T in exon 4 of RDH5 existing in sheep were found, and the T allele was completely fixed in Hu sheep. In addition, the mutation with the C allele reduced the retinol dehydrogenase activity encoding by RDH5, which can impair retinoic acid metabolism and further influenced the visual cycle. Overall, our results showed significant enrichment for positively selected genes involved in sensory perception development during sheep domestication; RDH5 and its variants may be related to the retinal degeneration in sheep. We infer that the wild sheep ancestors with weaker visual sensitivity were weeded out by humans, and the mutation was selective, swept by the dual pressures of natural and artificial selection.

Candidate Genes and Gene Networks Change with Age in Japanese Black Cattle by Blood Transcriptome Analysis.

Age is an important physiological factor that affects the metabolism and immune function of beef cattle. While there have been many studies using the blood transcriptome to study the effects of age on gene expression, few have been reported on beef cattle. To this end, we used the blood transcriptomes of Japanese black cattle at different ages as the study subjects and screened 1055, 345, and 1058 differential expressed genes (DEGs) in the calf vs. adult, adult vs. old, and calf vs. old comparison groups, respectively. The weighted co-expression network consisted of 1731 genes. Finally, blue, brown, and yellow age-specific modules were obtained, in which genes were enriched in signaling pathways related to growth and development and immune metabolic dysfunction, respectively. Protein-protein interaction (PPI) analysis showed gene interactions in each specific module, and 20 of the highest connectivity genes were chosen as potential hub genes. Finally, we identified 495, 244, and 1007 genes by exon-wide selection signature (EWSS) analysis of different comparison groups. Combining the results of hub genes, we found that VWF, PARVB, PRKCA, and TGFB1I1 could be used as candidate genes for growth and development stages of beef cattle. CORO2B and SDK1 could be used as candidate marker genes associated with aging. In conclusion, by comparing the blood transcriptome of calves, adult cattle, and old cattle, the candidate genes related to immunity and metabolism affected by age were identified, and the gene co-expression network of different age stages was constructed. It provides a data basis for exploring the growth, development, and aging of beef cattle.

Selection signature analysis and genome-wide divergence of South African Merino breeds from their founders.

Merino sheep are a breed of choice across the world, popularly kept for their wool and mutton value. They are often reared as a pure breed or used in crossbreeding and are a common component in synthetic breed development. This study evaluated genetic diversity, population structure, and breed divergence in 279 animals of Merino and Merino-based sheep breeds in South Africa using the Illumina Ovine SNP 50K BeadChip. The sheep breeds analysed included the three Merino-derived breeds of Dohne Merino (n = 50); Meatmaster (n = 47); and Afrino (n = 52) and five presumed ancestral populations of Merinos (Merino (n = 46); South African Merino (n = 10); and South African Mutton Merino (n = 8)); and the non-Merino founding breeds of Damara (n = 20); Ronderib Afrikaner (n = 17); and Nguni (n = 29). Highest genetic diversity values were observed in the Dohne Merino (DM), with H o = 0.39 ± 0.01, followed by the Meatmaster and South African Merino (SAM), with H o = 0.37 ± 0.03. The level of inbreeding ranged from 0.0 ± 0.02 (DM) to 0.27 ± 0.05 (Nguni). Analysis of molecular variance (AMOVA) showed high within-population variance (>80%) across all population categories. The first principal component (PC1) separated the Merino, South African Mutton Merino (SAMM), DM, and Afrino (AFR) from the Meatmaster, Damara, Nguni, and Ronderib Afrikaner (RDA). PC2 aligned each Merino-derived breed with its presumed ancestors and separated the SAMM from the Merino and SAM. The iHS analysis yielded selection sweeps across the AFR (12 sweeps), Meatmaster (four sweeps), and DM (29 sweeps). Hair/wool trait genes such as FGF12; metabolic genes of ICA1, NXPH1, and GPR171; and immune response genes of IL22, IL26, IFNAR1, and IL10RB were reported. Other genes include HMGA, which was observed as selection signatures in other populations; WNT5A, important in the development of the skeleton and mammary glands; ANTXR2, associated with adaptation to variation in climatic conditions; and BMP2, which has been reported as strongly selected in both fat-tailed and thin-tailed sheep. The DM vs. SAMM shared all six sweep regions on chromosomes 1, 10, and 11 with AFR vs. SAMM. Genes such as FGF12 on OAR 1:191.3-194.7Mb and MAP2K4 on OAR 11:28.6-31.3Mb were observed. The selection sweep on chromosome 10 region 28.6-30.3Mb harbouring the RXFP2 for polledness was shared between the DM vs. Merino, the Meatmaster vs. Merino, and the Meatmaster vs. Nguni. The DM vs. Merino and the Meatmaster vs. Merino also shared an Rsb-based selection sweep on chromosome 1 region 268.5-269.9Mb associated with the Calpain gene, CAPN7. The study demonstrated some genetic similarities between the Merino and Merino-derived breeds emanating from common founding populations and some divergence driven by breed-specific selection goals. Overall, information regarding the evolution of these composite breeds from their founding population will guide future breed improvement programs and management and conservation efforts.

Selection Signature Analysis for Important Economic Traits in Sheep (Ovis aries)

Background: Many genomic imprinting left during long-term selective evolution may drive the formation of important economical traits of introduced and indigenous sheep breeds. Methods: In this study, OvineSNP50 BeadChip data of the introduced breeds GMM and DOR and the indigenous breed Sunit were downloaded from public databases. Selective sweep base on Pi and Fst strategy initially screened for a number of candidate genes that could explain the genetic mechanism elucidated the genetic-molecular mechanisms underlying the differentiation of growth development and reproductive traits between introduced and indigenous breeds. Result: In the DOR_VS_GMM group, 32 genes were identified. In the DOR_VS_SUN group, 26 genes were identified. In the GMM_VS_SUN group, 25 genes were identified. GO, KEGG enrichment analysis showed that these genes are mainly involved in the pathways of muscle development and reproductive regulation. We speculate that EHBP1, CSRP1, TNNI1, MBNL1, AADACL3, RDH13, TNNI3, BMP2, Wfdc8, PPP1CC2 and MOV10L1 may have left deep marks on the genome during long-term domestication and evolution and may be responsible for the productive and reproductive performance in the Merino and Sunit sheep. These genes have played an important role in the evolution of selection in the three breeds. It is expected that the selected candidate genes will provide some insights for future sheep germplasm mining and foreign introductions.

Whole-Genome Sequence Analysis Reveals the Origin of the Chakouyi Horse.

The Chakouyi horse is an ancient Chinese indigenous horse breed distributed in Gansu Province in northwestern China, and is also one of the key breeds protected by the government. However, the origin of the Chakouyi horse remains unclear. As it is distributed in a key region of the Silk Road, it was speculated that the origin of the Chakouyi horse might involve the foreign horse breeds found along this ancient commercial artery. In this study, whole-genome resequencing data of 12 horse breeds, including both indigenous and foreign horses, were applied to reveal the genetic relationships between the Chakouyi horse and other breeds, as well as the ancestry of this ancient breed. An analysis of the population structure and admixture showed that there is no close genetic affinity between the Chakouyi horse and the foreign horses while Chinese indigenous horse populations were grouped together in accordance with their geographic locations, and the Chakouyi horse showed a closer relationship with Kazak horses, Mongolian horses, and Tibetan horses. The results from the ancestral composition prediction indicated that the Kazak horse and the Mongolian horse might be two ancestors of the Chakouyi horse. Furthermore, the genome-wide selection signature analysis revealed that the DMRT3 gene was positively selected in the Chakouyi horse and related to the gait trait of the breed. Our results provide insights into the native origin of the Chakouyi horse and indicate that Kazak and Mongolian horses played important roles in the formation of the Chakouyi horse. Genetic communication between the Chakouyi horse and other horse populations could be attributed, at least partially, to population migrations and trade activities along the ancient commercial routes.

Bayesian model and selection signature analyses reveal risk factors for canine atopic dermatitis.

Canine atopic dermatitis is an inflammatory skin disease with clinical similarities to human atopic dermatitis. Several dog breeds are at increased risk for developing this disease but previous genetic associations are poorly defined. To identify additional genetic risk factors for canine atopic dermatitis, we here apply a Bayesian mixture model adapted for mapping complex traits and a cross-population extended haplotype test to search for disease-associated loci and selective sweeps in four dog breeds at risk for atopic dermatitis. We define 15 associated loci and eight candidate regions under selection by comparing cases with controls. One associated locus is syntenic to the major genetic risk locus (Filaggrin locus) in human atopic dermatitis. One selection signal in common type Labrador retriever cases positions across the TBC1D1 gene (body weight) and one signal of selection in working type German shepherd controls overlaps the LRP1B gene (brain), near the KYNU gene (psoriasis). In conclusion, we identify candidate genes, including genes belonging to the same biological pathways across multiple loci, with potential relevance to the pathogenesis of canine atopic dermatitis. The results show genetic similarities between dog and human atopic dermatitis, and future across-species genetic comparisons are hereby further motivated.

Population genetic structure analysis and identification of backfat thickness loci of Chinese synthetic Yunan pigs.

Yunan is a crossed lean meat pig breed in China. Backfat thickness is the gold standard for carcass quality grading. However, over 14years after breed registration, the backfat of Yunan thickened and the consistency of backfat thickness decreased. Meanwhile, no genetic study has been ever performed on Yunan population. So, in this study we collected all the 120 nucleus individuals of Yunan and recorded six backfat traits of them, carried out population genetic structure analysis, selection signals analysis and genome-wide association study of Yunan pigs with the help of their founder population Duroc and Chinese native Huainan pigs, to determine the genomic loci on backfat of Yunan. Genetic diversity indexes suggested Yunan pigs had no inbreeding risk while population genetic structure showed they had few molecular pedigrees and were stratified. A total of 71 common selection signals affecting growth and fat deposition were detected by F ST and XP-CLR methods. 34 significant loci associated with six backfat traits were detected, among which a 1.40Mb region on SSC4 (20.03-21.43Mb) were outstanding as the strong region underlying backfat. This region was common with the results of selection signature analysis, former reported QTLs for backfat and was common for different kinds of backfat traits at different development stage. ENPP2, EXT1 and SLC30A8 genes around were fat deposition related genes and were of Huainan pig's origin, among which Type 2 diabetes related gene SLC30A8 was the most reasonable for being in a 193.21Kb haplotype block of the 1.40Mb region. Our results had application value for conservation, mating and breeding improvement of backfat thickness of Yunan pigs and provided evidence for a human function gene might be reproduced in pigs.

A ∼ 4.1 kb deletion in IRX1 gene upstream is completely associated with rumplessness in Piao chicken

Piao chicken, a Chinese indigenous rumpless chicken breed, lacks pygostyle, caudal vertebra, uropygial gland and tail feathers. The rumplessness in Piao chicken presents an autosomal dominant inheritance pattern. However, the molecular genetic mechanisms underlying the rumplessness in Piao chicken remains unclear. In this study, whole-genome resequencing was performed for 146 individuals from 10 chicken breeds, including 9 tailed chicken breeds and Piao rumpless breed. Tailbone CT scan for Piao chickens and WL chickens, revealed that some Piao chicken tails were normal in number, and for a few Piao chickens tail length and tail bone numbers were between the rumpless and the normal tailed chickens. The results showed that the rumpless phenotype has not been completely fixed in Piao chicken breed. Using selection signature analysis and structural variation detection, we found a 4174 bp deletion located in the upstream region of IRX1 gene on chromosome 2 related to rumpless phenotype. Structural variation genotyping showed that the deletion was present in all 32 rumpless Piao chickens (del/del, wild/del) and absent from all 112 tailed chickens included in the dataset for the other 9 breeds and 2 tailed Piao chickens (wild/wild). In summary, all rumpless Piao chickens tested here carry this deletion mutation, to show a complete linkage association with rumplessness trait. We suggested that the 4174 bp deletion could be causative for rumpless phenotype in Piao chicken since this is the only mutation to show the complete linkage disequilibrium with rumplessness on whole genome level across all of 146 chickens from the 10 breeds. This study could facilitate a better understanding of the genetic characteristics of Piao chicken.

Candidate genes for litter size in Xinjiang sheep identified by Specific Locus Amplified Fragment (SLAF) sequencing

The aim of this study was to investigate the selection signatures at a genome-wide level in ‘Pishan’ sheep using Specific Locus Amplified Fragment (SLAF)-seq. Blood samples from 126 ewes were sequenced using SLAF tags, and the ovarian tissues from 8 ewes (Bashbay sheep, a single litter size group (SG group); ‘Pishan’ sheep, double litter size group (DG group)) were collected to detect expression levels by quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR). Selection signature analysis was performed using global fixation index (Fst) and nucleotide diversity (π) ratio. A total of 1,192,168 high-quality SLAFs were identified. Notably, 2380 candidate regions under selection using two approaches were identified. A total of 2069 genes were identified, which were involved in dopaminergic synapses, thyroid hormone synthesis, ovarian steroidogenesis and thyroid hormone signalling pathways. Furthermore, Growth Differentiation Factor 9 (GDF9), Period Circadian Regulator 2 (PER2), Thyroid Stimulating Hormone Receptor (TSHR), and Nuclear Receptor Coactivator 1 (NCOA1) reside within these regions and pathways. The expression levels of GDF9 and PER2 genes in sheep tissue of the DG group were significantly higher than those in the SG group. These genes are interesting candidates for litter size and provide a starting point for further identification of conservation strategies for ‘Pishan’ sheep.

Blood Transcriptome Analysis of Beef Cow with Different Parity Revealed Candidate Genes and Gene Networks Regulating the Postpartum Diseases.

Maternal parity is an important physiological factor influencing beef cow reproductive performance. However, there are few studies on the influence of different calving periods on early growth and postpartum diseases. Here, we conducted blood transcriptomic analysis on cows of different parities for gene discovery. We used Short Time Series Expression Miner (STEM) analysis to determine gene expression levels in cows of various parities and divided multiple parities into three main periods (nulliparous, primiparous, and multiparous) for subsequent analysis. Furthermore, the top 15,000 genes with the lowest median absolute deviation (MAD) were used to build a co-expression network using weighted correlation network analysis (WGCNA), and six independent modules were identified. Combing with Exon Wide Selection Signature (EWSS) and protein-protein interaction (PPI) analysis revealed that TPCN2, KIF22, MICAL3, RUNX2, PDE4A, TESK2, GPM6A, POLR1A, and KLHL6 involved in early growth and postpartum diseases. The GO and KEGG enrichment showed that the Parathyroid hormone synthesis, secretion, and action pathway and stem cell differentiation function-related pathways were enriched. Collectively, our study revealed candidate genes and gene networks regulating the early growth and postpartum diseases and provided new insights into the potential mechanism of reproduction advantages of different parity selection.