Domestic-animal genomics: deciphering the genetics of complex traits

Beyond the hype: using AI, big data, wearable devices, and the internet of things for high-throughput livestock phenotyping

Phenotyping of animals is a routine task in agriculture which can provide large datasets for the functional annotation of genomes. Using the livestock farming sector to study complex traits enables genetics researchers to fully benefit from the digital transformation of society as economies of scale substantially reduces the cost of phenotyping animals on farms. In the agricultural sector genomics has transitioned towards a model of 'Genomics without the genes' as a large proportion of the genetic variation in animals can be modelled using the infinitesimal model for genomic breeding valuations. Combined with third generation sequencing creating pan-genomes for livestock the digital infrastructure for trait collection and precision farming provides a unique opportunity for high-throughput phenotyping and the study of complex traits in a controlled environment. The emphasis on cost efficient data collection mean that mobile phones and computers have become ubiquitous for cost-efficient large-scale data collection but that the majority of the recorded traits can still be recorded manually with limited training or tools. This is especially valuable in low- and middle income countries and in settings where indigenous breeds are kept at farms preserving more traditional farming methods. Digitalization is therefore an important enabler for high-throughput phenotyping for smaller livestock herds with limited technology investments as well as large-scale commercial operations. It is demanding and challenging for individual researchers to keep up with the opportunities created by the rapid advances in digitalization for livestock farming and how it can be used by researchers with or without a specialization in livestock. This review provides an overview of the current status of key enabling technologies for precision livestock farming applicable for the functional annotation of genomes.

Heritability of sudden cardiac death in Thoroughbred racehorses

BACKGROUND

Sudden cardiac death (SCD) is an important cause of exercise-associated fatalities in Thoroughbred racehorses. Equine deaths share similarities with fatalities in human athletes that result from inherited cardiac disease. Whilst genetic causes have been postulated in horses, these have not been confirmed and heritability of SCD has not previously been estimated in Thoroughbred racehorses.

OBJECTIVES

To determine the heritability of SCD in a sample population of Thoroughbred racehorses.

STUDY DESIGN

Retrospective case-control study.

METHODS

Steward and post-mortem reports of Thoroughbred racehorses in Australia between 2007 and 2020 were reviewed to identify horses with SCD. Control horses were randomly selected from races in which SCD occurred or from races on the date of the case fatality. A five-generation integrated pedigree chart was collected for each horse. Estimates of heritability were obtained using an animal model in the ASReml-R program with variance components estimated assuming SCD was normally distributed, and on the logit transformed scale. Inbreeding coefficients were calculated and the risk of producing SCD-affected progeny was calculated for stallions that sired ≥5 individuals in the case-control population.

RESULTS

Ninety-three horses with SCD and 465 control horses were identified. Heritability on the underlying scale was 0.15 ± 0.09 (logit animal) and 0.24 ± 0.12 (normal animal). Inbreeding coefficients were not significantly different between groups. Of the 16 first generation sires that appeared ≥5 times in the case-control data set, two sires more frequently produced affected progeny (OR 7.95-10.41).

MAIN LIMITATIONS

Challenges in definitively confirming SCD may lead to misclassification of some cases. Some control horses may have not been exposed to environmental influences of SCD. Case numbers are low and the studied population may not represent the entire Thoroughbred genetic pool.

CONCLUSION

The heritability of SCD in this population was relatively low. However, individual stallions appear more likely to produce affected progeny. Further studies are required to understand the genetic and environmental influences that contribute to disease expression.

Genetic characteristics and selection signatures between Southern Chinese local and commercial chickens

The introduction of exotic breeds and the cultivation of new lines by breeding companies have posed challenges to native chickens in South China, including loss of breed characteristics, decreased genetic diversity, and declining purity. Understanding the population genetic structure and genetic diversity of native chickens in South China is crucial for further advancements in breeding efforts. In this study, we analyzed the population genetic structure and genetic diversity of 321 individuals from 10 different breeds in South China. By comparing commercial chickens with native ones, we identified selection signatures occurring between local chickens and commercial breeds. The analysis of population genetic structure revealed that the native chicken populations in South China exhibited a considerable level of genetic diversity. Moreover, the commercial lines of Xiaobai chicken and Huangma chicken displayed even higher levels of genetic diversity, which distinguished them from other native varieties at the clustering level. However, certain individuals within these commercial varieties showed a discernible genetic relationship with the native populations. Notably, both commercial varieties also retained a significant degree of genetic similarity to their respective native counterparts. In order to investigate the genomic changes occurring during the commercialization of native chickens, we employed 4 methods (Fst, ROD, XPCLR, and XPEHH) to identify potential candidate regions displaying selective signatures in Southern Chinese native chicken population. A total of 168 (identified by Fst and ROD) and 86 (identified by XPCLR and XPEHH) overlapping genes were discovered. Functional annotation analysis revealed that these genes may be associated with reproduction and growth (SAMSN1, HYLS1, ROBO3, FGF14, PRSS23), musculoskeletal development (DNER, MYBPC1, DGKB, ORC1, KLF10), disease resistance and environmental adaptability (PUS3, CRB2, CALD1, USP15, SGCD, LTBP1), as well as egg production (ADGRB3, ACSF3). Overall, native chickens in South China harbor numerous selective sweep regions compared to commercial chickens, enriching valuable genomic resources for future genetic research and breeding conservation.

Construction of a high-density linkage map and QTL detection for growth traits in South African abalone (Haliotis midae)

Haliotis midae is one of the most important molluscs in South African commercial aquaculture. In this study, a high-resolution integrated linkage map was constructed, and QTL identified using 2b-RADseq for genotyping SNPs in three families. The final integrated linkage map was composed by merging the individual family maps, resulting in 3290 informative SNPs mapping to 18 linkage groups, conforming to the known haploid chromosome number for H. midae. The total map spanned 1798.25 cM with an average marker interval of 0.55 cM, representing a genome coverage of 98.76%. QTL analysis, across all three families, resulted in a total of five QTL identified for growth-related traits, shell width, shell length, and total body weight. For shell width and total body weight, one QTL was identified for each trait respectively, whilst three QTL were identified for shell length. The identified QTL respectively explained between 7.20% and 11.40% of the observed phenotypic variance. All three traits were significantly correlated (r = 0.862-0.970; p < 0.01) and shared overlapping QTL. The QTL for growth traits were mapped back to the H. midae draft genome and BLAST searches revealed the identity of candidate genes, such as egf-1, megf10, megf6, tnx, sevp1, kcp, notch1, and scube2 with possible functional roles in H. midae growth. The constructed high-density linkage map and mapped QTL have given valuable insights regarding the genetic architecture of growth-related traits and will be important genetic resources for marker-assisted selection. It remains, however, important to validate causal variants through linkage disequilibrium fine mapping in future.

Molecular genetic identification of Wuzhishan ant chicken, a newly discovered resource in China

Introduction

The Wuzhishan ant (MY) chicken exhibits significant differences from other chicken breeds. However, the molecular genetic relationship between the MY breed and other chicken breeds has not been assessed.

Methods

Whole-genome resequencing was used to compare genetic diversity, nucleotide diversity, the fixation index, the linkage disequilibrium coefficient, and phylogenetic tree relationships between the MY breed and the Wenchang (WC), Danzhou (DZ), Bawangling (BW), and Longsheng Feng (LF) breeds.

Results

A total of 21,586,378 singlenucleotide polymorphisms and 4,253,341 insertions/deletions were screened out among the five breeds. The MY breed had the second highest genomic genetic diversity and nucleotide diversity and the lowest LD coefficient among the five breeds. Moreover, the phylogenetic tree analysis showed that individual birds of each breed clustered together with those of their respective breeds.

Discussion

Our data indicated that the MY breed is distinct from the other breeds and can be considered a new genetic resource.

PC Gene Affects Milk Production Traits in Dairy Cattle

In previous work, we found that PC was differentially expressed in cows at different lactation stages. Thus, we deemed that PC may be a candidate gene affecting milk production traits in dairy cattle. In this study, we found the polymorphisms of PC by resequencing and verified their genetic associations with milk production traits by using an animal model in a cattle population. In total, we detected six single-nucleotide polymorphisms (SNPs) in PC. The single marker association analysis showed that all SNPs were significantly associated with the five milk production traits (p < 0.05). Additionally, we predicted that allele G of 29:g.44965658 in the 5' regulatory region created binding sites for TF GATA1 and verified that this allele inhibited the transcriptional activity of PC by the dual-luciferase reporter assay. In conclusion, we proved that PC had a prominent genetic effect on milk production traits, and six SNPs with prominent genetic effects could be used as markers for genomic selection (GS) in dairy cattle, which is beneficial for accelerating the improvement in milk yield and quality in Chinese Holstein cows.

Whole-genome resequencing reveals genomic variation and dynamics in Ethiopian indigenous goats

Ethiopia has about 52 million indigenous goats with marked phenotypic variability, which is the outcome of natural and artificial selection. Here, we obtained whole-genome sequence data of three Ethiopian indigenous goat populations (Arab, Fellata, and Oromo) from northwestern Ethiopia and analyzed their genome-wide genetic diversity, population structure, and signatures of selection. We included genotype data from four other Ethiopian goat populations (Abergelle, Keffa, Gumuz, and Woyto-Guji) and goats from Asia; Europe; and eastern, southern, western, and northern Africa to investigate the genetic predisposition of the three Ethiopian populations and performed comparative genomic analysis. Genetic diversity analysis showed that Fellata goats exhibited the lowest heterozygosity values (Ho = 0.288 ± 0.005 and He = 0.334 ± 0.0001). The highest values were observed in Arab goats (Ho = 0.310 ± 0.010 and He = 0.347 ± 4.35e-05). A higher inbreeding coefficient (FROH = 0.137 ± 0.016) was recorded for Fellata goats than the 0.105 ± 0.030 recorded for Arab and the 0.112 ± 0.034 recorded for Oromo goats. This indicates that the Fellata goat population should be prioritized in future conservation activities. The three goat populations showed the majority (∼63%) of runs of homozygosity in the shorter (100-150 Kb) length category, illustrating ancient inbreeding and/or small founder effects. Population relationship and structure analysis separated the Ethiopian indigenous goats into two distinct genetic clusters lacking phylogeographic structure. Arab, Fellata, Oromo, Abergelle, and Keffa represented one genetic cluster. Gumuz and Woyto-Guji formed a separate cluster and shared a common genetic background with the Kenyan Boran goat. Genome-wide selection signature analysis identified nine strongest regions spanning 163 genes influencing adaptation to arid and semi-arid environments (HOXC12, HOXC13, HOXC4, HOXC6, and HOXC9, MAPK8IP2), immune response (IL18, TYK2, ICAM3, ADGRG1, and ADGRG3), and production and reproduction (RARG and DNMT1). Our results provide insights into a thorough understanding of genetic architecture underlying selection signatures in Ethiopian indigenous goats in a semi-arid tropical environment and deliver valuable information for goat genetic improvement, conservation strategy, genome-wide association study, and marker-assisted breeding.

A genome-wide association study on growth traits of Korean commercial pig breeds using Bayesian methods

OBJECTIVE

This study aims to identify the significant regions and candidate genes of growth-related traits (adjusted backfat thickness [ABF], average daily gain [ADG], and days to 90 kg [DAYS90]) in Korean commercial GGP pig (Duroc, Landrace, and Yorkshire) populations.

METHODS

A genome-wide association study (GWAS) was performed using single-nucleotide polymorphism (SNP) markers for imputation to Illumina PorcineSNP60. The BayesB method was applied to calculate thresholds for the significance of SNP markers. The identified windows were considered significant if they explained ≥1% genetic variance.

RESULTS

A total of 28 window regions were related to genetic growth effects. Bayesian GWAS revealed 28 significant genetic regions including 52 informative SNPs associated with growth traits (ABF, ADG, DAYS90) in Duroc, Landrace, and Yorkshire pigs, with genetic variance ranging from 1.00% to 5.46%. Additionally, 14 candidate genes with previous functional validation were identified for these traits.

CONCLUSION

The identified SNPs within these regions hold potential value for future markerassisted or genomic selection in pig breeding programs. Consequently, they contribute to an improved understanding of genetic architecture and our ability to genetically enhance pigs. SNPs within the identified regions could prove valuable for future marker-assisted or genomic selection in pig breeding programs.

Whole genome sequencing of mouse lines divergently selected for fatness (FLI) and leanness (FHI) revealed several genetic variants as candidates for novel obesity genes

BACKGROUND

Analysing genomes of animal model organisms is widely used for understanding the genetic basis of complex traits and diseases, such as obesity, for which only a few mouse models exist, however, without their lean counterparts.

OBJECTIVE

To analyse genetic differences in the unique mouse models of polygenic obesity (Fat line) and leanness (Lean line) originating from the same base population and established by divergent selection over more than 60 generations.

METHODS

Genetic variability was analysed using WGS. Variants were identified with GATK and annotated with Ensembl VEP. g.Profiler, WebGestalt, and KEGG were used for GO and pathway enrichment analysis. miRNA seed regions were obtained with miRPathDB 2.0, LncRRIsearch was used to predict targets of identified lncRNAs, and genes influencing adipose tissue amount were searched using the IMPC database.

RESULTS

WGS analysis revealed 6.3 million SNPs, 1.3 million were new. Thousands of potentially impactful SNPs were identified, including within 24 genes related to adipose tissue amount. SNP density was highest in pseudogenes and regulatory RNAs. The Lean line carries SNP rs248726381 in the seed region of mmu-miR-3086-3p, which may affect fatty acid metabolism. KEGG analysis showed deleterious missense variants in immune response and diabetes genes, with food perception pathways being most enriched. Gene prioritisation considering SNP GERP scores, variant consequences, and allele comparison with other mouse lines identified seven novel obesity candidate genes: 4930441H08Rik, Aff3, Fam237b, Gm36633, Pced1a, Tecrl, and Zfp536.

CONCLUSION

WGS revealed many genetic differences between the lines that accumulated over the selection period, including variants with potential negative impacts on gene function. Given the increasing availability of mouse strains and genetic polymorphism catalogues, the study is a valuable resource for researchers to study obesity.

Cis-eQTLs in seven duck tissues identify novel candidate genes for growth and carcass traits

BACKGROUND

Expression quantitative trait loci (eQTL) studies aim to understand the influence of genetic variants on gene expression. The colocalization of eQTL mapping and GWAS strategy could help identify essential candidate genes and causal DNA variants vital to complex traits in human and many farm animals. However, eQTL mapping has not been conducted in ducks. It is desirable to know whether eQTLs within GWAS signals contributed to duck economic traits.

RESULTS

In this study, we conducted an eQTL analysis using publicly available RNA sequencing data from 820 samples, focusing on liver, muscle, blood, adipose, ovary, spleen, and lung tissues. We identified 113,374 cis-eQTLs for 12,266 genes, a substantial fraction 39.1% of which were discovered in at least two tissues. The cis-eQTLs of blood were less conserved across tissues, while cis-eQTLs from any tissue exhibit a strong sharing pattern to liver tissue. Colocalization between cis-eQTLs and genome-wide association studies (GWAS) of 50 traits uncovered new associations between gene expression and potential loci influencing growth and carcass traits. SRSF4, GSS, and IGF2BP1 in liver, NDUFC2 in muscle, ELF3 in adipose, and RUNDC1 in blood could serve as the candidate genes for duck growth and carcass traits.

CONCLUSIONS

Our findings highlight substantial differences in genetic regulation of gene expression across duck primary tissues, shedding light on potential mechanisms through which candidate genes may impact growth and carcass traits. Furthermore, this availability of eQTL data offers a valuable resource for deciphering further genetic association signals that may arise from ongoing extensive endeavors aimed at enhancing duck production traits.

Convenient synthesis and delivery of a megabase-scale designer accessory chromosome empower biosynthetic capacity

Synthetic biology confers new functions to hosts by introducing exogenous genetic elements, yet rebuilding complex traits that are based on large-scale genetic information remains challenging. Here, we developed a CRISPR/Cas9-mediated haploidization method that bypasses the natural process of meiosis. Based on the programmed haploidization in yeast, we further developed an easy-to-use method designated HAnDy (Haploidization-based DNA Assembly and Delivery in yeast) that enables efficient assembly and delivery of large DNA, with no need for any fussy in vitro manipulations. Using HAnDy, a de novo designed 1.024 Mb synthetic accessory chromosome (synAC) encoding 542 exogenous genes was parallelly assembled and then directly transferred to six phylogenetically diverse yeasts. The synAC significantly promotes hosts' adaptations and increases the scope of the metabolic network, which allows the emergence of valuable compounds. Our approach should facilitate the assembly and delivery of large-scale DNA for expanding and deciphering complex biological functions.

Elucidation of the genetic determination of body weight and size in Chinese local chicken breeds by large-scale genomic analyses

BACKGROUND

Body weight and size are important economic traits in chickens. While many growth-related quantitative trait loci (QTLs) and candidate genes have been identified, further research is needed to confirm and characterize these findings. In this study, we investigate genetic and genomic markers associated with chicken body weight and size. This study provides new insights into potential markers for genomic selection and breeding strategies to improve meat production in chickens.

METHODS

We performed whole-genome resequencing of and Wenshang Barred (WB) chickens (n = 596) and three additional breeds with varying body sizes (Recessive White (RW), WB, and Luxi Mini (LM) chickens; (n = 50)). We then used selective sweeps of mutations coupled with genome-wide association study (GWAS) to identify genomic markers associated with body weight and size.

RESULTS

We identified over 9.4 million high-quality single nucleotide polymorphisms (SNPs) among three chicken breeds/lines. Among these breeds, 287 protein-coding genes exhibited positive selection in the RW and WB populations, while 241 protein-coding genes showed positive selection in the LM and WB populations. Genomic heritability estimates were calculated for 26 body weight and size traits, including body weight, chest breadth, chest depth, thoracic horn, body oblique length, keel length, pelvic width, shank length, and shank circumference in the WB breed. The estimates ranged from 0.04 to 0.67. Our analysis also identified a total of 2,522 genome-wide significant SNPs, with 2,474 SNPs clustered around two genomic regions. The first region, located on chromosome 4 (7.41-7.64 Mb), was linked to body weight after ten weeks and body size traits. LCORL, LDB2, and PPARGC1A were identified as candidate genes in this region. The other region, located on chromosome 1 (170.46-171.53 Mb), was associated with body weight from four to eighteen weeks and body size traits. This region contained CAB39L and WDFY2 as candidate genes. Notably, LCORL, LDB2, and PPARGC1A showed highly selective signatures among the three breeds of chicken with varying body sizes.

CONCLUSION

Overall this study provides a comprehensive map of genomic variants associated with body weight and size in chickens. We propose two genomic regions, one on chromosome 1 and the other on chromosome 4, that could helpful for developing genome selection breeding strategies to enhance meat yield in chickens.

Genomic insights into Yak (Bos grunniens) adaptations for nutrient assimilation in high-altitudes

High-altitude environments present formidable challenges for survival and reproduction, with organisms facing limited oxygen availability and scarce nutrient resources. The yak (Bos grunniens), indigenous to the Tibetan Plateau, has notably adapted to these extreme conditions. This study delves into the genomic basis of the yak's adaptation, focusing on the positive selection acting on genes involved in nutrient assimilation pathways. Employing techniques in comparative genomics and molecular evolutionary analyses, we selected genes in the yak that show signs of positive selection associated with nutrient metabolism, absorption, and transport. Our findings reveal specific genetic adaptations related to nutrient metabolism in harsh climatic conditions. Notably, genes involved in energy metabolism, oxygen transport, and thermoregulation exhibited signs of positive selection, suggesting their crucial role in the yak's successful colonization of high-altitude regions. The study also sheds light on the yak's immune system adaptations, emphasizing genes involved in response to various stresses prevalent at elevated altitudes. Insights into the yak's genomic makeup provide valuable information for understanding the broader implications of high-altitude adaptations in mammalian evolution. They may contribute to efforts in enhancing livestock resilience to environmental challenges.

Managing genomic diversity in conservation programs of Chinese domestic chickens

BACKGROUND

Effective conservation and utilization of farm animals are fundamental for realizing sustainable increases in food production. In situ and ex situ conservation are the two main strategies that are currently used to protect the genetic integrity of Chinese domestic chicken breeds. However, genomic diversity and population structure have not been compared in these conserved populations.

RESULTS

Three hundred and sixty-one individuals from three Chinese domestic chicken breeds were collected from populations conserved in situ and ex situ and genotyped using genotyping-by-sequencing (GBS). First, we used different parameters based on heterozygosity, genomic inbreeding, and linkage disequilibrium to estimate the genomic diversity of these populations, and applied principal component analysis (PCA), neighbor-joining tree, and ADMIXTURE to analyze population structure. We found that the small ex situ conserved populations, which have been maintained in controlled environments, retained less genetic diversity than the in situ conserved populations. In addition, genetic differentiation was detected between the in situ and ex situ conserved populations of the same breed. Next, we analyzed signatures of selection using three statistical methods (fixation index (FST), nucleotide diversity (Pi), and cross-population extended haplotype homozygosity (XP-EHH) to study the genetic footprints that underlie the differentiation between in situ and ex situ conserved populations. We concluded that, in these small populations, differentiation might be caused by genetic drift or by mutations from the original populations. The differentiation observed in the population of Beijing You chicken probably reflects adaptation to environmental changes in temperature and humidity that the animals faced when they were moved from their place of origin to the new site for ex situ conservation.

CONCLUSIONS

Conservation programs of three Chinese domestic chicken breeds have maintained their genomic diversity to a sustainable degree. The small ex situ conserved populations, which are maintained in controlled environments, retain less genetic diversity than populations conserved in situ. In addition, the transfer of populations from their place of origin to another site for conservation purposes results in genetic differentiation, which may be caused by genetic drift or adaptation. This study provides a basis for further optimization of in situ and ex situ conservation programs for domestic chicken breeds in China.

Exploring the mechanism of artificial selection signature in Chinese indigenous pigs by leveraging multiple bioinformatics database tools

BACKGROUND

Chinese indigenous pigs in Yunnan exhibit considerable phenotypic diversity, but their population structure and the biological interpretation of signatures of artificial selection require further investigation. To uncover population genetic diversity, migration events, and artificial selection signatures in Chinese domestic pigs, we sampled 111 Yunnan pigs from four breeds in Yunnan which is considered to be one of the centres of livestock domestication in China, and genotyped them using Illumina Porcine SNP60K BeadChip. We then leveraged multiple bioinformatics database tools to further investigate the signatures and associated complex traits.

RESULTS

Population structure and migration analyses showed that Diannanxiaoer pigs had different genetic backgrounds from other Yunnan pigs, and Gaoligongshan may undergone the migration events from Baoshan and Saba pigs. Intriguingly, we identified a possible common target of sharing artificial selection on a 265.09 kb region on chromosome 5 in Yunnan indigenous pigs, and the genes on this region were associated with cardiovascular and immune systems. We also detected several candidate genes correlated with dietary adaptation, body size (e.g., PASCIN1, GRM4, ITPR2), and reproductive performance. In addition, the breed-sharing gene MMP16 was identified to be a human-mediated gene. Multiple lines of evidence at the mammalian genome, transcriptome, and phenome levels further supported the evidence for the causality between MMP16 variants and the metabolic diseases, brain development, and cartilage tissues in Chinese pigs. Our results suggested that the suppression of MMP16 would directly lead to inactivity and insensitivity of neuronal activity and skeletal development in Chinese indigenous pigs.

CONCLUSION

In this study, the population genetic analyses and identification of artificial selection signatures of Yunnan indigenous pigs help to build an understanding of the effect of human-mediated selection mechanisms on phenotypic traits in Chinese indigenous pigs. Further studies are needed to fully characterize the process of human-mediated genes and biological mechanisms.

Genetic variation in the OPN gene affects milk composition in Chinese Holstein cows

The aim of this study was to investigate the association between genotypes and haplotypes of OPN, and milk composition in dairy cows. A total of 317 Chinese Holstein cows were genotyped via DNA sequencing in this study. Three single nucleotide polymorphisms (SNPs), g.2916G > A, g.58675C > T and g.58899C > A, and eight haplotypes were identified. Of the eight possible haplotypes, four haplotypes i.e., Hap2 (ACC; 55.30%), Hap6 (GCC, 15.6%), Hap1 (ACA, 13.6%) and Hap4 (ATC, 5.70%), were considered to be major with a cumulative estimated frequency of >90%. Single markers (g.2916G > A and g.58899C > A) and Haplotype Hap6/4 were found to be associated with an increase in butter-fat percentage (p < 0.05). Taken together, our results provided evidence that polymorphisms in OPN are associated with milk composition, and could potentially be used for marker-assisted selection in Chinese Holstein cows.

Anthropogenic hybridization and its influence on the adaptive potential of the Sardinian wild boar (Sus scrofa meridionalis)

The wild boar (Sus scrofa meridionalis) arrived in Sardinia with the first human settlers in the early Neolithic with the potential to hybridize with the domestic pig (S. s. domesticus) throughout its evolution on the island. In this paper, we investigated the possible microevolutionary effects of such introgressive hybridization on the present wild boar population, comparing Sardinian wild specimens with several commercial pig breeds and Sardinian local pigs, along with a putatively unadmixed wild boar population from Central Italy, all genotyped with a medium density SNP chip. We first aimed at identifying hybrids in the population using different approaches, then examined genomic regions enriched for domestic alleles in the hybrid group, and finally we applied two methods to find regions under positive selection to possibly highlight instances of domestic adaptive introgression into a wild population. We found three hybrids within the Sardinian sample (3.1% out of the whole dataset). We reported 11 significant windows under positive selection with a method that looks for overly differentiated loci in the target population, compared with other two populations. We also identified 82 genomic regions with signs of selection in the domestic pig but not in the wild boar, two of which overlapped with genomic regions enriched for domestic alleles in the hybrid pool. Genes in these regions can be linked with reproductive success. Given our results, domestic introgression does not seem to be pervasive in the Sardinian wild boar. Nevertheless, we suggest monitoring the possible spread of advantageous domestic alleles in the coming years.

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.

Comparative neurogenetics of dog behavior complements efforts towards human neuropsychiatric genetics

Domestic dogs display a wide array of heritable behaviors that have intermediate genetic complexity thanks to a long history of human-influenced selection. Comparative genetics in dogs could address the scarcity of non-human neurogenetic systems relevant to human neuropsychiatric disorders, which are characterized by mental, emotional, and behavioral symptoms and involve vastly complex genetic and non-genetic risk factors. Our review describes the diverse behavioral "phenome" of domestic dogs, past and ongoing sources of behavioral selection, and the state of canine behavioral genetics. We highlight two naturally disordered behavioral domains that illustrate how dogs may prove useful as a comparative, forward neurogenetic system: canine age-related cognitive dysfunction, which can be examined more rapidly given the attenuated lifespan of dogs, and compulsive disorders, which may have genetic roots in purpose-bred behaviors. Growing community science initiatives aimed at the companion dog population will be well suited to investigating such complex behavioral phenotypes and offer a comparative resource that parallels human genomic initiatives in scale and dimensionality.

Unveiling the Mysteries of Non-Mendelian Heredity in Plant Breeding

Mendelian heredity is the cornerstone of plant breeding and has been used to develop new varieties of plants since the 19th century. However, there are several breeding cases, such as cytoplasmic inheritance, methylation, epigenetics, hybrid vigor, and loss of heterozygosity (LOH), where Mendelian heredity is not applicable, known as non-Mendelian heredity. This type of inheritance can be influenced by several factors besides the genetic architecture of the plant and its breeding potential. Therefore, exploring various non-Mendelian heredity mechanisms, their prevalence in plants, and the implications for plant breeding is of paramount importance to accelerate the pace of crop improvement. In this review, we examine the current understanding of non-Mendelian heredity in plants, including the mechanisms, inheritance patterns, and applications in plant breeding, provide an overview of the various forms of non-Mendelian inheritance (including epigenetic inheritance, cytoplasmic inheritance, hybrid vigor, and LOH), explore insight into the implications of non-Mendelian heredity in plant breeding, and the potential it holds for future research.

Formation, Application, and Significance of Chicken Primordial Germ Cells: A Review

Chicken is one of the most widely consumed sources of protein globally. Primordial germ cells (PGCs) are the precursors for ova and sperm. One of the early embryogenesis events in most animals is the segregation of the somatic and germ lineages. PGC cultures occur in the germline, and PGCs are less studied in many species. It is relatively challenging to separate, cultivate, and genetically alter chicken without mutating the basic germline. The present study aims to gather previous research about chicken PGCs and provide a customized review of studies and developments in the field of PGCs, especially for avian species. Furthermore, we show that the propagation of chicken PGCs into embryonic germ cells that contribute to somatic tissues may be produced in vitro. Primordial germ cells offer an ideal system in developmental biology, as these cells play a vital role in the genetic modification and treatment of infertility. Cryopreservation helps to maintain genetic resources and sustainable production in the poultry industry. Keeping in mind the significance of cryopreservation for storage and gametogenesis, we discuss its role in the preservation of primordial germ cells. Transgenesis and genetic modifications in chicken lead to the development of various medicinal chicken varieties and aid in improving their production and quality for consumption purposes. Additionally, these characteristics open up new possibilities for modifying the chicken genome for agricultural and medical purposes.

Genomics of Dwarfism in Italian Local Chicken Breeds

The identification of the dwarf phenotype in chicken is based on body weight, height, and shank length, leaving the differentiation between dwarf and small breeds ambiguous. The aims of the present study were to characterize the sequence variations associated with the dwarf phenotype in three Italian chicken breeds and to investigate the genes associated with their phenotype. Five hundred and forty-one chickens from 23 local breeds (from 20 to 24 animals per breed) were sampled. All animals were genotyped with the 600 K chicken SNP array. Three breeds were described as “dwarf”, namely, Mericanel della Brianza (MERI), Mugellese (MUG), and Pepoi (PPP). We compared MERI, MUG, and PPP with the four heaviest breeds in the dataset by performing genome-wide association studies. Results showed significant SNPs associated with dwarfism in the MERI and MUG breeds, which shared a candidate genomic region on chromosome 1. Due to this similarity, MERI and MUG were analyzed together as a meta-population, observing significant SNPs in the LEMD3 and HMGA2 genes, which were previously reported as being responsible for dwarfism in different species. In conclusion, MERI and MUG breeds seem to share a genetic basis of dwarfism, which differentiates them from the small PPP breed.

Selection signatures for fiber production in commercial species: A review

Natural fibers derived from diverse animal species have gained increased attention in recent years due to their favorable environmental effects, long-term sustainability benefits, and remarkable physical and mechanical properties that make them valuable raw materials used for textile and non-textile production. Domestication and selective breeding for the economically significant fiber traits play an imperative role in shaping the genomes and, thus, positively impact the overall productivity of the various fiber-producing species. These selection pressures leave unique footprints on the genome due to alteration in the allelic frequencies at specific loci, characterizing selective sweeps. Recent advances in genomics have enabled the discovery of selection signatures across the genome using a variety of methods. The increased demand for 'green products' manufactured from natural fibers necessitates a detailed investigation of the genomes of the various fiber-producing plant and animal species to identify the candidate genes associated with important fiber attributes such as fiber diameter/fineness, color, length, and strength, among others. The objective of this review is to present a comprehensive overview of the concept of selection signature and selective sweeps, discuss the main methods used for its detection, and address the selection signature studies conducted so far in the diverse fiber-producing animal species.

Morphostructural Differences between the Historical Genetic Lines of the Spanish Merino Sheep

The Merino breed, which originates from Spain, is the most emblematic livestock breed in the world, since it is the first with a worldwide extension and has had an important impact on the genetic origin of several of the main current sheep populations. For this reason, it is of vital importance to typify the historical genetic lines of the original Spanish Merino breed and thereby ensure the conservation of its variability. In the present study, we used 337 purebred animals (males and females) registered in the Genealogical Book of the Native Merino Breed. All the animals were descendants of herds from six ancestral genetic lines (Maesso, Egea, Granda, López-Montenegro, Hidalgo, and Donoso). Significant differences were found in all the morphometric traits and indexes between the different genetic lines. Using discriminant analysis, 84% of the animals were classified correctly into their historical genetic lines. Furthermore, the distances between the lines, calculated by a cluster test, showed that Hidalgo, Maesso, and Donoso had the most clearly defined lines, while the Granda, López-Montenegro, and Egea lines were more similar to each other. All this demonstrates the rich genetic variability existing in the genuine gene pool of the Merino sheep breed.

Common protein-coding variants influence the racing phenotype in galloping racehorse breeds

Selection for system-wide morphological, physiological, and metabolic adaptations has led to extreme athletic phenotypes among geographically diverse horse breeds. Here, we identify genes contributing to exercise adaptation in racehorses by applying genomics approaches for racing performance, an end-point athletic phenotype. Using an integrative genomics strategy to first combine population genomics results with skeletal muscle exercise and training transcriptomic data, followed by whole-genome resequencing of Asian horses, we identify protein-coding variants in genes of interest in galloping racehorse breeds (Arabian, Mongolian and Thoroughbred). A core set of genes, G6PC2, HDAC9, KTN1, MYLK2, NTM, SLC16A1 and SYNDIG1, with central roles in muscle, metabolism, and neurobiology, are key drivers of the racing phenotype. Although racing potential is a multifactorial trait, the genomic architecture shaping the common athletic phenotype in horse populations bred for racing provides evidence for the influence of protein-coding variants in fundamental exercise-relevant genes. Variation in these genes may therefore be exploited for genetic improvement of horse populations towards specific types of racing.

Three-Dimensional Organization of Chicken Genome Provides Insights into Genetic Adaptation to Extreme Environments

The high-throughput chromosome conformation capture (Hi-C) technique is widely used to study the functional roles of the three-dimensional (3D) architecture of genomes. However, the knowledge of the 3D genome structure and its dynamics during extreme environmental adaptations remains poor. Here, we characterized 3D genome architectures using the Hi-C technique for chicken liver cells. Upon comparing Lindian chicken (LDC) liver cells with Wenchang chicken (WCC) liver cells, we discovered that environmental adaptation contributed to the switching of A/B compartments, the reorganization of topologically associated domains (TADs), and TAD boundaries in both liver cells. In addition, the analysis of the switching of A/B compartments revealed that the switched compartmental genes (SCGs) were strongly associated with extreme environment adaption-related pathways, including tight junction, notch signaling pathway, vascular smooth muscle contraction, and the RIG-I-like receptor signaling pathway. The findings of this study advanced our understanding of the evolutionary role of chicken 3D genome architecture and its significance in genome activity and transcriptional regulation.

Low-coverage sequencing in a deep intercross of the Virginia body weight lines provides insight to the polygenic genetic architecture of growth: novel loci revealed by increased power and improved genome-coverage

Genetic dissection of highly polygenic traits is a challenge, in part due to the power necessary to confidently identify loci with minor effects. Experimental crosses are valuable resources for mapping such traits. Traditionally, genome-wide analyses of experimental crosses have targeted major loci using data from a single generation (often the F₂) with individuals from later generations being generated for replication and fine-mapping. Here, we aim to confidently identify minor-effect loci contributing to the highly polygenic basis of the long-term, bi-directional selection responses for 56-d body weight in the Virginia body weight chicken lines. To achieve this, a strategy was developed to make use of data from all generations (F₂-F₁₈) of the advanced intercross line, developed by crossing the low and high selected lines after 40 generations of selection. A cost-efficient low-coverage sequencing based approach was used to obtain high-confidence genotypes in 1Mb bins across 99.3% of the chicken genome for >3,300 intercross individuals. In total, 12 genome-wide significant, and 30 additional suggestive QTL reaching a 10% FDR threshold, were mapped for 56-d body weight. Only 2 of these QTL reached genome-wide significance in earlier analyses of the F₂ generation. The minor-effect QTL mapped here were generally due to an overall increase in power by integrating data across generations, with contributions from increased genome-coverage and improved marker information content. The 12 significant QTL explain >37% of the difference between the parental lines, three times more than 2 previously reported significant QTL. The 42 significant and suggestive QTL together explain >80%. Making integrated use of all available samples from multiple generations in experimental crosses are economically feasible using the low-cost, sequencing-based genotyping strategies outlined here. Our empirical results illustrate the value of this strategy for mapping novel minor-effect loci contributing to complex traits to provide a more confident, comprehensive view of the individual loci that form the genetic basis of the highly polygenic, long-term selection responses for 56-d body weight in the Virginia body weight chicken lines.

A composite strategy of genome-wide association study and copy number variation analysis for carcass traits in a Duroc pig population

BACKGROUND

Carcass traits are important in pig breeding programs for improving pork production. Understanding the genetic variants underlies complex phenotypes can help explain trait variation in pigs. In this study, we integrated a weighted single-step genome-wide association study (wssGWAS) and copy number variation (CNV) analyses to map genetic variations and genes associated with loin muscle area (LMA), loin muscle depth (LMD) and lean meat percentage (LMP) in Duroc pigs.

RESULTS

Firstly, we performed a genome-wide analysis for CNV detection using GeneSeek Porcine SNP50 Bead chip data of 3770 pigs. A total of 11,100 CNVs were detected, which were aggregated by overlapping 695 CNV regions (CNVRs). Next, we investigated CNVs of pigs from the same population by whole-genome resequencing. A genome-wide analysis of 21 pigs revealed 23,856 CNVRs that were further divided into three categories (851 gain, 22,279 loss, and 726 mixed), which covered 190.8 Mb (~ 8.42%) of the pig autosomal genome. Further, the identified CNVRs were used to determine an overall validation rate of 68.5% for the CNV detection accuracy of chip data. CNVR association analyses identified one CNVR associated with LMA, one with LMD and eight with LMP after applying stringent Bonferroni correction. The wssGWAS identified eight, six and five regions explaining more than 1% of the additive genetic variance for LMA, LMD and LMP, respectively. The CNVR analyses and wssGWAS identified five common regions, of which three regions were associated with LMA and two with LMP. Four genes (DOK7, ARAP1, ELMO2 and SLC13A3) were highlighted as promising candidates according to their function.

CONCLUSIONS

We determined an overall validation rate for the CNV detection accuracy of low-density chip data and constructed a genomic CNV map for Duroc pigs using resequencing, thereby proving a value genetic variation resource for pig genome research. Furthermore, our study utilized a composite genetic strategy for complex traits in pigs, which will contribute to the study for elucidating the genetic architecture that may be influenced and regulated by multiple forms of variations.

Trajectory of livestock genomics in South Asia: A comprehensive review

Livestock plays a central role in sustaining human livelihood in South Asia. There are numerous and distinct livestock species in South Asian countries. Several of them have experienced genetic development in recent years due to the application of genomic technologies and effective breeding programs. This review discusses genomic studies on cattle, buffalo, sheep, goat, pig, horse, camel, yak, mithun, and poultry. The frontiers covered in this review are genetic diversity, admixture studies, selection signature research, QTL discovery, genome-wide association studies (GWAS), and genomic selection. The review concludes with recommendations for South Asian livestock systems to increasingly leverage genomic technologies, based on the lessons learned from the numerous case studies. This paper aims to present a comprehensive analysis of the dichotomy in the South Asian livestock sector and argues that a realistic approach to genomics in livestock can ensure long-term genetic advancements.

Single Nucleotide Polymorphisms of ALDH18A1 and MAT2A Genes and Their Genetic Associations with Milk Production Traits of Chinese Holstein Cows

Our preliminary work had suggested two genes, aldehyde dehydrogenase 18 family member A1 (ALDH18A1) and methionine adenosyltransferase 2A (MAT2A), related to amino acid synthesis and metabolism as candidates affecting milk traits by analyzing the liver transcriptome and proteome of dairy cows at different lactation stages. In this study, the single nucleotide polymorphisms (SNPs) of ALDH18A1 and MAT2A genes were identified and their genetic effects and underlying causative mechanisms on milk production traits in dairy cattle were analyzed, with the aim of providing effective genetic information for the molecular breeding of dairy cows. By resequencing the entire coding and partial flanking regions of ALDH18A1 and MAT2A, we found eight SNPs located in ALDH18A1 and two in MAT2A. Single-SNP association analysis showed that most of the 10 SNPs of these two genes were significantly associated with the milk yield traits, 305-day milk yield, fat yield, and protein yield in the first and second lactations (corrected p ≤ 0.0488). Using Haploview 4.2, we found that the seven SNPs of ALDH18A1 formed two haplotype blocks; subsequently, the haplotype-based association analysis showed that both haplotypes were significantly associated with 305-day milk yield, fat yield, and protein yield (corrected p ≤ 0.014). Furthermore, by Jaspar and Genomatix software, we found that 26:g.17130318 C>A and 11:g.49472723G>C, respectively, in the 5′ flanking region of ALDH18A1 and MAT2A genes changed the transcription factor binding sites (TFBSs), which might regulate the expression of corresponding genes to affect the phenotypes of milk production traits. Therefore, these two SNPs were considered as potential functional mutations, but they also require further verification. In summary, ALDH18A1 and MAT2A were proved to probably have genetic effects on milk production traits, and their valuable SNPs might be used as candidate genetic markers for dairy cattle’s genomic selection (GS).

Whole Genome Analyses Reveal Novel Genes Associated with Chicken Adaptation to Tropical and Frigid Environments

INTRODUCTION

Investigating the genetic footprints of historical temperature selection can get insights to the local adaptation and feasible influences of climate change on long-term population dynamics.

OBJECT

Chicken is a significative species to study genetic adaptation on account of its similar domestication track related to human activity with the most diversified varieties. Yet, few studies have demonstrated the genetic signatures of its adaptation to naturally tropical and frigid environments.

METHOD

Here, we generated whole genome resequencing of 119 domesticated chickens in China including the following breeds which are in order of breeding environmental temperature from more tropical to more frigid: Wenchang chicken (WCC), green-shell chicken (GSC), Tibetan chicken (TBC), and Lindian chicken (LDC).

RESULTS

Our results showed WCC branched off earlier than LDC with an evident genetic admixture between WCC and LDC, suggesting their closer genetic relationship. Further comparative genomic analyses solute carrier family 33 member 1 (SLC33A1) and thyroid stimulating hormone receptor (TSHR) genes exhibited stronger signatures for positive selection in the genome of the more tropical WCC. Furthermore, genotype data from about 3,000 African local ecotypes confirmed that allele frequencies of single nucleotide polymorphisms (SNPs) in these 2 genes appeared strongly associated with tropical environment adaptation. In addition, the NADH:ubiquinone oxidoreductase subunit S4 (NDUFS4) gene exhibited a strong signature for positive selection in the LDC genome, and SNPs with marked allele frequency differences indicated a significant relationship with frigid environment adaptation.

CONCLUSION

Our findings partially clarify how selection footprints from environmental temperature stress can lead to advantageous genomic adaptions to tropical and frigid environments in poultry and provide a valuable resource for selective breeding of chickens.

Camel Proteins and Enzymes: A Growing Resource for Functional Evolution and Environmental Adaptation

In less agroecological parts of the Asian, Arabian, and African deserts, Camelus dromedarius play an important role in human survival. For many years, camels have been employed as a source of food, a tool of transportation, and a means of defense. They are becoming increasingly important as viable livestock animals in many desert climates. With the help of camel genetics, genomics and proteomics known so far, this review article will summarize camel enzymes and proteins, which allow them to thrive under varied harsh environmental situations. An in-depth study of the dromedary genome revealed the existence of protein-coding and fast-developing genes that govern a variety of metabolic responses including lipid and protein metabolism, glucoamylase, flavin-containing monooxygenase and guanidinoacetate methyltransferase are other metabolic enzymes found in the small intestine, liver, pancreas, and spleen. In addition, we will discuss the handling of common medications by camel liver cytochrome p 450, which are different from human enzymes. Moreover, camels developed several paths to get optimum levels of trace elements like copper, zinc, selenium, etc., which have key importance in their body for normal regulation of metabolic events. Insulin tolerance, carbohydrate and energy metabolism, xenobiotics metabolizing enzymes, vimentin functions, behavior during the rutting season, resistance to starvation and changes in blood composition and resistance to water loss were among the attractive aspects of camel enzymes and proteins peculiarities in the camels. Resolving the enigma of the method of adaptation and the molecular processes linked with camel life is still a developing repository full of mysteries that need additional exploration.

Effect of genetic variants in the SMAD1 and SMAD5 genes promoter on growth and beef quality traits in cattle

The SMAD1 and SMAD5 genes belong to mothers against decapentaplegic proteins family, which participate in the BMP pathway to control skeletal myogenesis and growth. In the present study, we analyzed the associations between polymorphisms of SMAD1 and SMAD5 genes promoter and important economical traits in Qinchuan cattle. Four SNPs in the SMAD1 gene promoter and three SNPs in the SMAD5 promoter were identified by sequencing of 448 Qinchuan cattles. Allelic and frequency analyses of these SNPs resulted in eight haplotypes both in the promoters of the two genes promoter and identified potential cis-regulatory transcription factor (TF) components. In addition, correlation analysis showed that cattle SMAD1 promoter activity of individuals with Hap4 (P < 0.01) was stronger than that of individuals with Hap2. while the transcriptional activity of individuals with Hap3 within SMAD5 gene promoter was significantly (P < 0.01) higher followed by H2. Uniformly, diplotypes H4-H6 of SMAD1 gene and H1-H3 of SMAD5 gene performed significant (P < 0.01) associations with body measurement and improved carcass quality traits. All these results have indicated that polymorphisms in SMAD1 and SMAD5 genes promoter could impact the transcriptional regulation and then affect muscle content in beef cattle. Moreover, both the SMAD1 and SMAD5 genes were expressed ubiquitously in 10 tissues and had higher expression in the longissimus thoracis tissue from 6-month-old and 12-month-old cattle than in cattle of other ages. We can conclude that SMAD1 and SMAD5 genes may play an important role in muscle growth and development, and the variants mapped within SMAD1 and SMAD5 genes can be utilized in molecular marker-assisted selection for cattle carcass quality and body measurement traits in breed improvement programs of Qinchuan cattle.

Whole-Genome Analyses Reveal Genomic Characteristics and Selection Signatures of Lincang Humped Cattle at the China-Myanmar Border

The location on the Yunnan border with Myanmar and its unique cultural landscape has shaped Lincang humped cattle over time. In the current study, we investigated the genetic characteristics of 22 Lincang humped cattle using whole-genome resequencing data. We found that Lincang humped cattle derived from both Indian indicine and Chinese indicine cattle depicted higher levels of genomic diversity. Based on genome-wide scans, candidate genomic regions were identified that were potentially involved in local thermal and humid environmental adaptions, including genes associated with the body size (TCF12, SENP2, KIF1C, and PFN1), immunity (LIPH, IRAK3, GZMM, and ELANE), and heat tolerance (MED16, DNAJC8, HSPA4, FILIP1L, HELB, BCL2L1, and TPX2). Missense mutations were detected in candidate genes IRAK3, HSPA4, and HELB. Interestingly, eight missense mutations observed in the HELB gene were specific to the indicine cattle pedigree. These mutations may reveal differences between indicine and taurine cattle adapted to variable climatic conditions. Our research provides new insights into the genetic characteristics of Lincang humped cattle representing Lincang and Pu'er areas as an important channel for the migration of Indian indicine from domestication centers toward southwestern China.

A large deletion at the cortex locus eliminates butterfly wing patterning

As the genetic basis of natural and domesticated variation has been described in recent years, a number of hotspot genes have been repeatedly identified as the targets of selection, Heliconius butterflies display a spectacular diversity of pattern variants in the wild and the genetic basis of these patterns has been well-described. Here, we sought to identify the mechanism behind an unusual pattern variant that is instead found in captivity, the ivory mutant, in which all scales on both the wings and body become white or yellow. Using a combination of autozygosity mapping and coverage analysis from 37 captive individuals, we identify a 78-kb deletion at the cortex wing patterning locus, a gene which has been associated with wing pattern evolution in H. melpomene and 10 divergent lepidopteran species. This deletion is undetected among 458 wild Heliconius genomes samples, and its dosage explains both homozygous and heterozygous ivory phenotypes found in captivity. The deletion spans a large 5′ region of the cortex gene that includes a facultative 5′UTR exon detected in larval wing disk transcriptomes. CRISPR mutagenesis of this exon replicates the wing phenotypes from coding knock-outs of cortex, consistent with a functional role of ivory-deleted elements in establishing scale color fate. Population demographics reveal that the stock giving rise to the ivory mutant has a mixed origin from across the wild range of H. melpomene, and supports a scenario where the ivory mutation occurred after the introduction of cortex haplotypes from Ecuador. Homozygotes for the ivory deletion are inviable while heterozygotes are the targets of artificial selection, joining 40 other examples of allelic variants that provide heterozygous advantage in animal populations under artificial selection by fanciers and breeders. Finally, our results highlight the promise of autozygosity and association mapping for identifying the genetic basis of aberrant mutations in captive insect populations.

Identification of important genomic footprints using eight different selection signature statistics in domestic cattle breeds

In the present study, the population genomic data of different cattle breeds were explored to decipher the genomic regions affected due to selective events and reflected in the productive, reproductive, thermo-tolerance, and health-related traits. To find out these genomic deviations due to selective sweeps, we used eight different statistical tools (Tajima's D, Fu & Li's D*, CLR, ROH, iHS, F_ST, FLK, and hapFLK) on seven indigenous and five exotic cattle breeds. We further performed composite analysis by comparing their covariance matrix. Several candidate genes were found to be related to milk production (ADARB, WDR70, and CA8), reproductive (PARN, FAM134B2, and ZBTB20), and health-related traits (SP110, CXCL2, CLXCL3, CXCL5, IRF8, and MYOM1). The outcome of this investigation provides a basis for detecting selective sweeps that explain the genetic variation of traits. They may possess functional importance for multiple cattle breeds in different subcontinents. However, further studies are required to improve the findings using high-density arrays or whole-genome sequencing with higher resolution and greater sample sizes.

[Opportunities and risks of the use of genetic resistances to infectious diseases in pigs - an overview]

Demands for health, performance and welfare in pigs, as well as the desire for consumer protection and reduced antibiotic use, require optimal measures in advance of disease development. This includes, in principle, the use of genetically more resistant lines and breeding animals, whose existence has been proven for a wide range of pathogen-host interactions. In addition, attempts are being made to identify the gene variants responsible for disease resistance in order to force the selection of suitable populations, also using modern biotechnical technics. The present work is intended to provide an overview of the research status achieved in this context and to highlight opportunities and risks for the future.The evaluation of the international literature shows that genetic disease resistance exist in many areas of swine diseases. However, polygenic inheritance, lack of animal models and the influence of environmental factors during evaluation render their implementation in practical breeding programs demanding. This is where modern molecular genetic methods, such as Gene Editing, come into play. Both approaches possess their pros and cons, which are discussed in this paper. The most important infectious diseases in pigs, including general diseases and epizootics, diseases of the respiratory and digestive tract and diseases of the immune system are taken into account.

Assessing Bos taurus introgression in the UOA Bos indicus assembly

Background

Reference genomes are essential in the analysis of genomic data. As the cost of sequencing decreases, multiple reference genomes are being produced within species to alleviate problems such as low mapping accuracy and reference allele bias in variant calling that can be associated with the alignment of divergent samples to a single reference individual. The latest reference sequence adopted by the scientific community for the analysis of cattle data is ARS_UCD1.2, built from the DNA of a Hereford cow (Bos taurus taurus—B. taurus). A complementary genome assembly, UOA_Brahman_1, was recently built to represent the other cattle subspecies (Bos taurus indicus—B. indicus) from a Brahman cow haplotype to further support analysis of B. indicus data. In this study, we aligned the sequence data of 15 B. taurus and B. indicus breeds to each of these references.

Results

The alignment of B. taurus individuals against UOA_Brahman_1 detected up to five million more single-nucleotide variants (SNVs) compared to that against ARS_UCD1.2. Similarly, the alignment of B. indicus individuals against ARS_UCD1.2 resulted in one and a half million more SNVs than that against UOA_Brahman_1. The number of SNVs with nearly fixed alternative alleles also increased in the alignments with cross-subspecies. Interestingly, the alignment of B. taurus cattle against UOA_Brahman_1 revealed regions with a smaller than expected number of counts of SNVs with nearly fixed alternative alleles. Since B. taurus introgression represents on average 10% of the genome of Brahman cattle, we suggest that these regions comprise taurine DNA as opposed to indicine DNA in the UOA_Brahman_1 reference genome. Principal component and admixture analyses using genotypes inferred from this region support these taurine-introgressed loci. Overall, the flagged taurine segments represent 13.7% of the UOA_Brahman_1 assembly. The genes located within these segments were previously reported to be under positive selection in Brahman cattle, and include functional candidate genes implicated in feed efficiency, development and immunity.

Conclusions

We report a list of taurine segments that are in the UOA_Brahman_1 assembly, which will be useful for the interpretation of interesting genomic features (e.g., signatures of selection, runs of homozygosity, increased mutation rate, etc.) that could appear in future re-sequencing analysis of indicine cattle.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00688-1.

Detecting the differential genomic variants using cross-population phenotype-associated variant (XP-PAV) of the Landrace and Yorkshire pigs in Korea

Although there have been many genome-wide association studies (GWAS) and selective sweep analyses to understand pig genomic regions related to growth performance, these methods considered only the gene effect and selection signal, respectively. In this study, we suggest the cross-population phenotype associated variant (XP-PAV) analysis as a novel method to determine the genomic variants with different effects between the two populations. XP-PAV analysis could reveal the differential genetic variants between the two populations by considering the gene effect and selection signal simultaneously. In this study, we used daily weight gain (DWG) and back fat thickness (BF) as phenotypes and the Landrace and Yorkshire populations were used for XP-PAV analysis. The main aim was to reveal the differential selection by considering the gene effect between Landrace and Yorkshire pigs. In the gene ontology analysis of XP-PAV results, differential selective genes in DWG analysis were involved in the regulation of interleukin-2 production and cell cycle G2/M transition. The protein modification and glycerophospholipid biosynthetic processes were the most enriched terms in the BF analysis. Therefore, we could identify genetic differences for immune and several metabolic pathways between Landrace and Yorkshire breeds using the XP-PAV analysis. In this study, we expect that XP-PAV analysis will play a role in determining useful selective variants with gene effects and provide a new interpretation of the genetic differences between the two populations.

An Integrative Analysis of Transcriptome and GWAS Data to Identify Potential Candidate Genes Influencing Meat Quality Traits in Pigs

Understanding the genetic factors behind meat quality traits is of great significance to animal breeding and production. We previously conducted a genome-wide association study (GWAS) for meat quality traits in a White Duroc × Erhualian F2 pig population using Illumina porcine 60K SNP data. Here, we further investigate the functional candidate genes and their network modules associated with meat quality traits by integrating transcriptomics and GWAS information. Quantitative trait transcript (QTT) analysis, gene expression QTL (eQTL) mapping, and weighted gene co-expression network analysis (WGCNA) were performed using the digital gene expression (DGE) data from 493 F2 pig's muscle and liver samples. Among the quantified 20,108 liver and 23,728 muscle transcripts, 535 liver and 1,014 muscle QTTs corresponding to 416 and 721 genes, respectively, were found to be significantly (p < 5 × 10^-4) correlated with 22 meat quality traits measured on longissiums dorsi muscle (LM) or semimembranosus muscle (SM). Transcripts associated with muscle glycolytic potential (GP) and pH values were enriched for genes involved in metabolic process. There were 42 QTTs (for 32 genes) shared by liver and muscle tissues, of which 10 QTTs represent GP- and/or pH-related genes, such as JUNB, ATF3, and PPP1R3B. Furthermore, a genome-wide eQTL mapping revealed a total of 3,054 eQTLs for all annotated transcripts in muscle (p < 2.08 × 10^-5), including 1,283 cis-eQTLs and 1771 trans-eQTLs. In addition, WGCNA identified five modules relevant to glycogen metabolism pathway and highlighted the connections between variations in meat quality traits and genes involved in energy process. Integrative analysis of GWAS loci, eQTL, and QTT demonstrated GALNT15/GALNTL2 and HTATIP2 as strong candidate genes for drip loss and pH drop from postmortem 45 min to 24 h, respectively. Our findings provide valuable insights into the genetic basis of meat quality traits and greatly expand the number of candidate genes that may be valuable for future functional analysis and genetic improvement of meat quality.

Genetic prediction of complex traits with polygenic scores: a statistical review

Accurate genetic prediction of complex traits can facilitate disease screening, improve early intervention, and aid in the development of personalized medicine. Genetic prediction of complex traits requires the development of statistical methods that can properly model polygenic architecture and construct a polygenic score (PGS). We present a comprehensive review of 46 methods for PGS construction. We connect the majority of these methods through a multiple linear regression framework which can be instrumental for understanding their prediction performance for traits with distinct genetic architectures. We discuss the practical considerations of PGS analysis as well as challenges and future directions of PGS method development. We hope our review serves as a useful reference both for statistical geneticists who develop PGS methods and for data analysts who perform PGS analysis.

Revealing the coexistence of differentiation and communication in an endemic hare, Lepus yarkandensis (Mammalia, Leporidae) using specific-length amplified fragment sequencing

BACKGROUND

The Yarkand hare (Lepus yarkandensis Günther, 1875) is endemic to oasis and desert areas around the Tarim Basin in the Xinjiang Uyghur Autonomous Region of northwest China; however, genome-wide information for this species remains limited. Moreover, the genetic variation, genetic structure, and phylogenetic relationships of Yarkand hare from the plateau mountain regions have not been reported. Thus, we used specific-length amplified fragment sequencing (SLAF-seq) technology to evaluate the genetic diversity of 76 Yarkand hares from seven geographic populations in the northern and southwestern parts of the Tarim Basin to investigate single-nucleotide polymorphism (SNP) marker-based population differentiation and evolutionary processes. Selective sweep analysis was conducted to identify genetic differences between populations.

RESULTS

Using SLAF-seq, a total of 1,835,504 SNPs were initially obtained, of which 308,942 high-confidence SNPs were selected for further analysis. Yarkand hares exhibited a relatively high degree of genetic diversity at the SNP level. Based on pairwise FST estimates, the north and southwest groups showed a moderate level of genetic differentiation. Phylogenetic tree and population structure analyses demonstrated evident systematic phylogeographical structure patterns consistent with the geographical distribution of the hares. Hierarchical analysis of molecular variation further indicated that genetic variation was mainly observed within populations. Low to moderate genetic differentiation also occurred among populations despite a common genomic background, likely due to geographical barriers, genetic drift, and differential selection pressure of distinct environments. Nevertheless, the observed lineage-mixing pattern, as indicated by the evolutionary tree, principal component analysis, population structure, and TreeMix analyses, suggests a certain degree of gene flow between the north and southwest groups. This may be related to the migration of hares to high-altitude water sources southwest of the basin during glacial climatic oscillations, as well as river re-diffusion and oasis restoration in the basin following the glacial period. We also identified candidate genes, and their associated gene ontology terms and pathways, related to the adaptation of Yarkand hares to different environmental habitats.

CONCLUSIONS

The identified genome-wide SNPs, genetic diversity, and population structure of Yarkand hares expand our understanding of the genetic background of this endemic species and provide valuable insights into its environmental adaptation, allowing for further exploration of the underlying mechanisms.

Average semivariance yields accurate estimates of the fraction of marker-associated genetic variance and heritability in complex trait analyses

The development of genome-informed methods for identifying quantitative trait loci (QTL) and studying the genetic basis of quantitative variation in natural and experimental populations has been driven by advances in high-throughput genotyping. For many complex traits, the underlying genetic variation is caused by the segregation of one or more ‘large-effect’ loci, in addition to an unknown number of loci with effects below the threshold of statistical detection. The large-effect loci segregating in populations are often necessary but not sufficient for predicting quantitative phenotypes. They are, nevertheless, important enough to warrant deeper study and direct modelling in genomic prediction problems. We explored the accuracy of statistical methods for estimating the fraction of marker-associated genetic variance (p) and heritability (HM2) for large-effect loci underlying complex phenotypes. We found that commonly used statistical methods overestimate p and HM2. The source of the upward bias was traced to inequalities between the expected values of variance components in the numerators and denominators of these parameters. Algebraic solutions for bias-correcting estimates of p and HM2 were found that only depend on the degrees of freedom and are constant for a given study design. We discovered that average semivariance methods, which have heretofore not been used in complex trait analyses, yielded unbiased estimates of p and HM2, in addition to best linear unbiased predictors of the additive and dominance effects of the underlying loci. The cryptic bias problem described here is unrelated to selection bias, although both cause the overestimation of p and HM2. The solutions we described are predicted to more accurately describe the contributions of large-effect loci to the genetic variation underlying complex traits of medical, biological, and agricultural importance.

The contributions of individual genes to the phenotypic variation observed for genetically complex traits has been an ongoing and important challenge in biology, medicine, and agriculture. While many genes have statistically undetectable effects, those with large effects often warrant in-depth study and can be important predictors of complex phenotypes such as disease risk in humans or disease resistance in domesticated plants and animals. The genes identified through associations with genetic markers in complex trait analyses typically account for a fraction of the heritable variation, a genetic parameter we called ‘marker heritability’. We discovered that textbook statistical methods systematically overestimate marker heritability and thus overestimate the contributions of specific genes to the phenotypic variation observed for complex traits in natural and experimental populations. We describe the source of the upward bias, validate our findings through computer simulation, describe methods for bias-correcting estimates of marker heritability, and illustrate their application through empirical examples. The statistical methods we describe supply investigators with more accurate estimates of the contributions of specific genes or networks of interacting genes to the heritable variation observed in complex trait studies.

Integration of selection signatures and multi-trait GWAS reveals polygenic genetic architecture of carcass traits in beef cattle

Carcass merits are widely considered as economically important traits affecting beef production in the beef cattle industry. However, the genetic basis of carcass traits remains to be well understood. Here, we applied multiple methods, including the Composite of Likelihood Ratio (CLR) and Genome-wide Association Study (GWAS), to explore the selection signatures and candidate variants affecting carcass traits. We identified 11,600 selected regions overlapping with 2214 candidate genes, and most of those were enriched in binding and gene regulation. Notably, we identified 66 and 110 potential variants significantly associated with carcass traits using single-trait and multi-traits analyses, respectively. By integrating selection signatures with single and multi-traits associations, we identified 12 and 27 putative genes, respectively. Several highly conserved missense variants were identified in OR5M13D, NCAPG, and TEX2. Our study supported polygenic genetic architecture of carcass traits and provided novel insights into the genetic basis of complex traits in beef cattle.

Genome-Wide Scanning for Signatures of Selection Revealed the Putative Genomic Regions and Candidate Genes Controlling Milk Composition and Coat Color Traits in Sahiwal Cattle

Background

In the evolutionary time scale, selection shapes the genetic variation and alters the architecture of genome in the organisms. Selection leaves detectable signatures at the genomic coordinates that provide clues about the protein-coding regions. Sahiwal is a valuable indicine cattle adapted to tropical environments with desirable milk attributes. Insights into the genomic regions under putative selection may reveal the molecular mechanisms affecting the quantitative and other important traits. To understand this, the present investigation was undertaken to explore signatures of selection in the genome of Sahiwal cattle using a medium-density genotyping INDUS chip.

Result

De-correlated composite of multiple selection signals (DCMS), which combines five different univariate statistics, was computed in the dataset to detect the signatures of selection in the Sahiwal genome. Gene annotations, Quantitative Trait Loci (QTL) enrichment, and functional analyses were carried out for the identification of significant genomic regions. A total of 117 genes were identified, which affect a number of important economic traits. The QTL enrichment analysis highlighted 14 significant [False Discovery Rate (FDR)-corrected p-value ≤ 0.05] regions on chromosomes BTA 1, 3, 6, 11, 20, and 21. The top three enriched QTLs were found on BTA 6, 20, and 23, which are associated with exterior, health, milk production, and reproduction traits. The present study on selection signatures revealed some key genes related with coat color (PDGFRA, KIT, and KDR), facial pigmentation (LEF), milk fat percent (MAP3K1, HADH, CYP2U1, and SGMS2), sperm membrane integrity (OSTC), lactation persistency (MRPS30, NNT, CCL28, HMGCS1, NIM1K, ZNF131, and CCDC152), milk yield (GHR and ZNF469), reproduction (NKX2-1 and DENND1A), and bovine tuberculosis susceptibility (RNF144B and PAPSS1). Further analysis of candidate gene prioritization identified four hub genes, viz., KIT, KDR, MAP3K1, and LEF, which play a role in coat color, facial pigmentation, and milk fat percentage in cattle. Gene enrichment analysis revealed significant Gene ontology (GO) terms related to breed-specific coat color and milk fat percent.

Conclusion

The key candidate genes and putative genomic regions associated with economic traits were identified in Sahiwal using single nucleotide polymorphism data and the DCMS method. It revealed selection for milk production, coat color, and adaptability to tropical climate. The knowledge about signatures of selection and candidate genes affecting phenotypes have provided a background information that can be further utilized to understand the underlying mechanism involved in these traits in Sahiwal cattle.

Ecology and population structure of some indigenous geese breeds and the impact of four GH and Pit-1 SNPs on their body weights

This study aims to determine the genetic correlation using nine microsatellite markers to reconstruct the history of some indigenous geese populations, along with the use of four single nucleotides polymorphisms (SNPs) to investigate their correlation with the geese body weight. Microsatellite markers are mainly used to provide updated information on changes in the population structure of geese breeds. The eight goose populations reported 24% private alleles specific for each population. Expected heterozygosity (He) ranged from 0.46 to 0.70. Three breeds were reported highly polymorphic. Inbreeding coefficient (Fis) revealed that three breeds were in a minimum level of extinction danger, while one breed was in a potential endangered situation. Phylogenetic tree, principal component analysis (PCA), and self-organizing map (SOM) were constructed using MATLAB to study the population distribution and relationship among these breeds. Four SNPs were detected, two SNPs at GH gene exon (C123T and C158T), and two SNPs at Pit-1 gene exons (G161A and T282G). Four SNP loci were reported to have a significant effect on geese body weight. They were CT genotype for C123T locus, TT genotype for C158T locus, GG genotype for G161A locus, and GG genotype for T282G locus.

Genome-wide scan reveals genetic divergence in Italian Holstein cows bred within PDO cheese production chains

Dairy cattle breeds have been exposed to intense artificial selection for milk production traits over the last fifty years. In Italy, where over 80% of milk is processed into cheese, selection has also focused on cheese-making traits. Due to a deep-rooted tradition in cheese-making, currently fifty Italian cheeses are marked with the Protected Designation of Origin (PDO) label as they proved traditional land of origin and procedures for milk transformation. This study aimed to explore from a genetic point of view if the presence of such diverse productive contexts in Italy have shaped in a different manner the genome of animals originally belonging to a same breed. We analyzed high density genotype data from 1000 Italian Holstein cows born between 2014 and 2018. Those animals were either farmed in one of four Italian PDO consortia or used for drinkable milk production only. Runs of Homozygosity, Bayesian Information Criterion and Discriminant Analysis of Principal Components were used to evaluate potential signs of genetic divergence within the breed. We showed that the analyzed Italian Holstein cows have genomic inbreeding level above 5% in all subgroups, reflecting the presence of ongoing artificial selection in the breed. Our study provided a comprehensive representation of the genetic structure of the Italian Holstein breed, highlighting the presence of potential genetic subgroups due to divergent dairy farming systems. This study can be used to further investigate genetic variants underlying adaptation traits in these subgroups, which in turn might be used to design more specialized breeding programs.

Association of SNPs in AKIRIN2, TTN, EDG1 and MYBPC1 genes with growth and carcass traits in Qinchuan cattle

Growth and carcass traits are the main breeding objectives in beef cattle. The aim of this study was to confirm genetic effects of the c.*188G>A SNP of AKIRIN2, the g.231054C>T SNP of TTN, the g.1471620G>T SNP of EDG1 and the g.70014208A>G SNP of MYBPC1 on growth and carcass traits in Chinese Qinchuan (QC) cattle, as well as to compare the frequencies of the well-characterized alleles of these SNPs among six Chinese cattle populations, three Japanese cattle populations, two European cattle populations and one Korean cattle population. In this study, a total of 665 cattle samples were genotyped using MassARRAY and PCR-RFLP. Association analysis explored effects of four SNPs on growth and carcass traits including body length, wither height, hip height, hip width, rump length, chest depth, chest circumference, back fat thickness, ultrasound longissimus muscle area and ultrasound longissimus muscle depth in QC (P<0.05 to P<0.001). The well-characterized A (c.*188G>A), T (g.231054C>T) and T (g.1471620G>T) alleles in Japanese Black cattle were significantly higher than Chinese cattle breeds, on the contrary, the G allele (g.70014208A>G) was markedly higher in Chinese cattle breeds than other cattle breeds. These results suggest that the four SNPs might be useful as a molecular marker for growth-related traits in Chinese QC cattle.

Detection of selection signatures in farmed coho salmon (Oncorhynchus kisutch) using dense genome-wide information

Animal domestication and artificial selection give rise to gradual changes at the genomic level in populations. Subsequent footprints of selection, known as selection signatures or selective sweeps, have been traced in the genomes of many animal livestock species by exploiting variation in linkage disequilibrium patterns and/or reduction of genetic diversity. Domestication of most aquatic species is recent in comparison with land animals, and salmonids are one of the most important fish species in aquaculture. Coho salmon (Oncorhynchus kisutch), cultivated primarily in Chile, has been subjected to breeding programs to improve growth, disease resistance traits, and flesh color. This study aimed to identify selection signatures that may be involved in adaptation to culture conditions and traits of productive interest. To do so, individuals of two domestic populations cultured in Chile were genotyped with 200 thousand SNPs, and analyses were conducted using iHS, XP-EHH and CLR. Several signatures of selection on different chromosomal regions were detected across both populations. Some of the identified regions under selection contained genes such anapc2, alad, chp2 and myn, which have been previously associated with body weight in Atlantic salmon, or sec24d and robo1, which have been associated with resistance to Piscirickettsia salmonis in coho salmon. Findings in our study can contribute to an integrated genome-wide map of selection signatures, to help identify the genetic mechanisms of phenotypic diversity in coho salmon.

Genome-wide insights on gastrointestinal nematode resistance in autochthonous Tunisian sheep

Gastrointestinal nematode (GIN) infections have negative impacts on animal health, welfare and production. Information from molecular studies can highlight the underlying genetic mechanisms that enhance host resistance to GIN. However, such information often lacks for traditionally managed indigenous livestock. Here, we analysed 600 K single nucleotide polymorphism genotypes of GIN infected and non-infected traditionally managed autochthonous Tunisian sheep grazing communal natural pastures. Population structure analysis did not find genetic differentiation that is consistent with infection status. However, by contrasting the infected versus non-infected cohorts using ROH, LR-GWAS, F_ST and XP-EHH, we identified 35 candidate regions that overlapped between at least two methods. Nineteen regions harboured QTLs for parasite resistance, immune capacity and disease susceptibility and, ten regions harboured QTLs for production (growth) and meat and carcass (fatness and anatomy) traits. The analysis also revealed candidate regions spanning genes enhancing innate immune defence (SLC22A4, SLC22A5, IL-4, IL-13), intestinal wound healing/repair (IL-4, VIL1, CXCR1, CXCR2) and GIN expulsion (IL-4, IL-13). Our results suggest that traditionally managed indigenous sheep have evolved multiple strategies that evoke and enhance GIN resistance and developmental stability. They confirm the importance of obtaining information from indigenous sheep to investigate genomic regions of functional significance in understanding the architecture of GIN resistance.