Fine mapping of a QTL for ear size on porcine chromosome 5 and identification of high mobility group AT-hook 2 (HMGA2) as a positional candidate gene

Global dispersal and adaptive evolution of domestic cattle: a genomic perspective

Domestic cattle have spread across the globe and inhabit variable and unpredictable environments. They have been exposed to a plethora of selective pressures and have adapted to a variety of local ecological and management conditions, including UV exposure, diseases, and stall-feeding systems. These selective pressures have resulted in unique and important phenotypic and genetic differences among modern cattle breeds/populations. Ongoing efforts to sequence the genomes of local and commercial cattle breeds/populations, along with the growing availability of ancient bovid DNA data, have significantly advanced our understanding of the genomic architecture, recent evolution of complex traits, common diseases, and local adaptation in cattle. Here, we review the origin and spread of domestic cattle and illustrate the environmental adaptations of local cattle breeds/populations.

Whole genome variants across 57 pig breeds enable comprehensive identification of genetic signatures that underlie breed features

Background

A large number of pig breeds are distributed around the world, their features and characteristics vary among breeds, and they are valuable resources. Understanding the underlying genetic mechanisms that explain across-breed variation can help breeders develop improved pig breeds.

Results

In this study, we performed GWAS using a standard mixed linear model with three types of genome variants (SNP, InDel, and CNV) that were identified from public, whole-genome, sequencing data sets. We used 469 pigs of 57 breeds, and we identified and analyzed approximately 19 million SNPs, 1.8 million InDels, and 18,016 CNVs. We defined six biological phenotypes by the characteristics of breed features to identify the associated genome variants and candidate genes, which included coat color, ear shape, gradient zone, body weight, body length, and body height. A total of 37 candidate genes was identified, which included 27 that were reported previously (e.g., PLAG1 for body weight), but the other 10 were newly detected candidate genes (e.g., ADAMTS9 for coat color).

Conclusion

Our study indicated that using GWAS across a modest number of breeds with high density genome variants provided efficient mapping of complex traits.

Supplementary Information

Supplementary information accompanies this paper at 10.1186/s40104-020-00520-8.

Whole-genome sequencing of European autochthonous and commercial pig breeds allows the detection of signatures of selection for adaptation of genetic resources to different breeding and production systems

Background

Natural and artificial directional selection in cosmopolitan and autochthonous pig breeds and wild boars have shaped their genomes and resulted in a reservoir of animal genetic diversity. Signatures of selection are the result of these selection events that have contributed to the adaptation of breeds to different environments and production systems. In this study, we analysed the genome variability of 19 European autochthonous pig breeds (Alentejana, Bísara, Majorcan Black, Basque, Gascon, Apulo-Calabrese, Casertana, Cinta Senese, Mora Romagnola, Nero Siciliano, Sarda, Krškopolje pig, Black Slavonian, Turopolje, Moravka, Swallow-Bellied Mangalitsa, Schwäbisch-Hällisches Schwein, Lithuanian indigenous wattle and Lithuanian White old type) from nine countries, three European commercial breeds (Italian Large White, Italian Landrace and Italian Duroc), and European wild boars, by mining whole-genome sequencing data obtained by using a DNA-pool sequencing approach. Signatures of selection were identified by using a single-breed approach with two statistics [within-breed pooled heterozygosity (H_P) and fixation index (F_ST)] and group-based F_ST approaches, which compare groups of breeds defined according to external traits and use/specialization/type.

Results

We detected more than 22 million single nucleotide polymorphisms (SNPs) across the 23 compared populations and identified 359 chromosome regions showing signatures of selection. These regions harbour genes that are already known or new genes that are under selection and relevant for the domestication process in this species, and that affect several morphological and physiological traits (e.g. coat colours and patterns, body size, number of vertebrae and teats, ear size and conformation, reproductive traits, growth and fat deposition traits). Wild boar related signatures of selection were detected across all the genome of several autochthonous breeds, which suggests that crossbreeding (accidental or deliberate) occurred with wild boars.

Conclusions

Our findings provide a catalogue of genetic variants of many European pig populations and identify genome regions that can explain, at least in part, the phenotypic diversity of these genetic resources.

HMGA Genes and Proteins in Development and Evolution

HMGA (high mobility group A) (HMGA1 and HMGA2) are small non-histone proteins that can bind DNA and modify chromatin state, thus modulating the accessibility of regulatory factors to the DNA and contributing to the overall panorama of gene expression tuning. In general, they are abundantly expressed during embryogenesis, but are downregulated in the adult differentiated tissues. In the present review, we summarize some aspects of their role during development, also dealing with relevant studies that have shed light on their functioning in cell biology and with emerging possible involvement of HMGA1 and HMGA2 in evolutionary biology.

Genomic clues of the evolutionary history of Bos indicus cattle

Together with their sister subspecies Bos taurus, zebu cattle (Bos indicus) have contributed to important socioeconomic changes that have shaped modern civilizations. Zebu cattle were domesticated in the Indus Valley 8000 years before present (YBP). From the domestication site, they expanded to Africa, East Asia, southwestern Asia and Europe between 4000 and 1300 YBP, intercrossing with B. taurus to form clinal variations of zebu ancestry across the landmass of Afro-Eurasia. In the past 150 years, zebu cattle reached the Americas and Oceania, where they have contributed to the prosperity of emerging economies. The zebu genome is characterized by two mitochondrial haplogroups (I1 and I2), one Y chromosome haplogroup (Y3) and three major autosomal ancestral groups (Indian-Pakistani, African and Chinese). Phenotypically, zebu animals are recognized by their hump, large ears and excess skin. They are rustic, resilient to parasites and capable of bearing the hot and humid climates of the tropics. Many resources are available to study the zebu genome, including commercial arrays of SNP, reference assemblies and publicly available genotypes and whole-genome sequences. Nevertheless, many of these resources were initially developed to support research and subsidize industrial applications in B. taurus, and therefore they can produce bias in data analysis. The combination of genomics with precision agriculture holds great promise for the identification of genetic variants affecting economically important traits such as tick resistance and heat tolerance, which were naturally selected for millennia and played a major role in the evolution of B. indicus cattle.

Construction of a genetic linkage map in Pyropia yezoensis (Bangiales, Rhodophyta) and QTL analysis of several economic traits of blades

Pyropia yezoensis is an economically important seaweed but its molecular genetics is poorly understood. In the present study, we used a doubled haploid (DH) population that was established in our previous work to construct a genetic linkage map of P. yezoensis and analyze the quantitative trait loci (QTLs) of blades. The DH population was genotyped with fluorescent sequence-related amplified polymorphism (SRAP) markers. A chi-square test identified 301 loci with normal segregation (P ≥ 0.01) and 96 loci (24.18%) with low-level skewed segregation (0.001 ≤ P < 0.01). The genetic map was constructed after a total of 92 loci were assembled into three linkage groups (LGs). The map spanned 557.36 cM covering 93.71% of the estimated genome, with a mean interlocus space of 6.23 cM. Kolmogorov-Smirnov test (α = 5%) showed a uniform distribution of the markers along each LG. On the genetic map, 10 QTLs associated with five economic traits of blades were detected. One QTL was for length, one for width, two for fresh weight, two for specific growth rate of length and four for specific growth rate of fresh weight. These QTLs could explain 2.29–7.87% of the trait variations, indicating that their effects were all minor. The results may serve as a framework for future marker-assisted breeding in P. yezoensis.

Molecular cloning of WIF1 and HMGA2 reveals ear-preferential expression while uncovering a missense mutation associated with porcine ear size in WIF1

Considerable diversity exists in porcine ear size, which is an important morphological feature of pig breeds. Previously, we localized four crucial candidate genes-high mobility group AT-hook 2 (HMGA2), LEM domain-containing 3 (LEMD3), methionine sulfoxide reductase B3 (MSRB3) and Wnt inhibitory factor 1 (WIF1)-on Sus Scrofa chromosome 5 affecting porcine ear size, then cloned LEMD3 and MSBR3. In this study, we performed rapid amplification of cDNA ends to obtain full-length cDNA sequences of 2338-bp WIF1 and 2998-bp HMGA2. Using quantitative real-time PCR, we revealed that WIF1 expression was highest in ear cartilage of 60-day-old pigs and that this is therefore a better candidate gene for ear size than HMGA2. We further screened coding sequence variants in both genes and identified only one missense mutation (WIF1:c.1167C>G) in a conserved epidermal growth factor-like domain from the mammalian WIF1 protein. The protein-altering mutation was significantly associated with ear size across the Large White × Minzhu hybrid and Beijing Black pig populations. When WIF1:c.1167C>G was included as fixed effect in the model to re-run a genome-wide association study in the Large White × Minzhu intercross population the P-value of the peak SNP on SSC5 from re-running the genome-wide association study dropped from 2.45E-12 to 7.33E-05. Taken together, the WIF1:c.1167C>G could be an important mutation associated with ear size. Our findings provide helpful information for further studies of the molecular mechanisms controlling porcine ear size.

Copy number variation in the MSRB3 gene enlarges porcine ear size through a mechanism involving miR-584-5p

Background

The size and type of ears are important conformation characteristics that distinguish pig breeds. A significant quantitative trait locus (QTL) for ear size has been identified on SSC5 (SSC for Sus scrofa chromosome) but the underlying causative gene and mutation remain unknown. Thus, our aim was to identify the gene responsible for enlarged ears in pig.

Results

First, we narrowed down the QTL region on SSC5 to a 137.85-kb interval that harbors only the methionine sulfoxide reductase B3 (MSRB3) gene. Then, we identified a 38.7-kb copy number variation (CNV) that affects the last two exons of MSRB3 and could be the candidate causative mutation for this QTL. This CNV showed complete concordance with genotype at the QTL of the founder animals in a white Duroc × Erhualian F₂ intercross and was found only in pigs from six Chinese indigenous breeds with large ears and from the Landrace breed with half-floppy ears. Moreover, it accounted for the significant association with ear size on SSC5 across the five pig populations tested. eQTL mapping revealed that this CNV was significantly associated with the expression of the microRNA (miRNA) miR-584-5p, which interacts with MSRB3, one of its target genes. In vivo and in vitro experiments confirmed that miR-584-5p inhibits the translation of MSRB3 mRNA. Taken together, these results led us to conclude that presence of the 38.7-kb CNV in the genome of some pig breeds affects ear size by altering the expression of miR-584-5p, which consequently hinders the expression of one of its target genes (e.g. MSRB3).

Conclusions

Our findings shed insight into the underlying mechanism of development of external ears in mammals and contribute to a better understanding of how the presence of CNV can regulate gene expression.

Electronic supplementary material

The online version of this article (10.1186/s12711-018-0442-6) contains supplementary material, which is available to authorized users.

Association of Copy Number Variation at Intron 3 of HMGA2 With Navel Length in Bos indicus

Navel injuries caused by friction against the pasture can promote infection, reproductive problems and costly treatments in beef cattle raised in extensive systems. A haplotype-based genome-wide association study (GWAS) was performed for visual scores of navel length at yearling in Nellore cattle (Bos indicus) using data from 2,016 animals and 503,088 single nucleotide polymorphism (SNP) markers. The strongest signal (p = 1.01 × 10^-9) was found on chromosome 5 spanning positions 47.9-48.2 Mbp. This region contains introns 3 and 4 and exons 4 and 5 of the high mobility group AT-hook 2 gene (HMGA2). Further inspection of the region with whole genome sequence data of 21 Nellore bulls revealed correlations between counts of the significant haplotype and copy number gains of a ∼6.2 kbp segment of intron 3 of HMGA2. Analysis of genome sequences from five African B. indicus and four European Bos taurus breeds revealed that the copy number variant (CNV) is indicine-specific. This intronic CNV was then validated through quantitative polymerase chain reaction (qPCR) using Angus animals as copy neutral controls. Importantly, the CNV was not detectable by means of conventional SNP-based GWAS or SNP probe intensity analyses. Given that HMGA2 affects the expression of the insulin-like growth factor 2 gene (IGF2) together with the pleomorphic adenoma gene 1 (PLAG1), and that the latter has been repeatedly shown to be associated with quantitative traits of economic importance in cattle, these findings highlight the emerging role of variants impacting the insulin-like growth factor pathway to cattle breeding.

Genome-wide association study reveals novel genes for the ear size in sheep (Ovis aries)

Variations in ear size can be observed in livestock such as sheep; however, the genetic basis of variable ear size in sheep is still poorly understood. To investigate causative genes associated with ear size in sheep, a genome-wide association study was performed in 115 adult Duolang sheep with different-sized floppy ears using the Ovine Infinium HD BeadChip. We found 38 significant SNPs at the genome-wide or chromosome-wise 5% significance level after Bonferroni correction. The most significant association (P = 1.61 × 10^-6 ) was found at SNP rs402740419, located in the DCC gene, which plays a critical role in ear development. Also, we observed two additional significant SNPs, rs407891215 in PTPRD and rs407769095 in SOX5, both of which are functionally associated with ear developmental processes. Our results are useful for future sheep breeding and provide insights into the genetic basis of ear size development in sheep and other livestock.

High mobility group A2 (HMGA2) deficiency in pigs leads to dwarfism, abnormal fetal resource allocation, and cryptorchidism

Expression of HMGA2 is strongly associated with body size and growth in mice and humans. In mice, inactivation of one or both alleles of Hmga2 results in body-size reductions of 20% and 60%, respectively. In humans, microdeletions involving the HMGA2 locus result in short stature, suggesting the function of the HMGA2 protein is conserved among mammals. To test this hypothesis, we generated HMGA2-deficient pigs via gene editing and somatic cell nuclear transfer (SCNT). Examination of growth parameters revealed that HMGA2^-/+ male and female pigs were on average 20% lighter and smaller than HMGA2^+/+ matched controls (P < 0.05). HMGA2^-/- boars showed significant size reduction ranging from 35 to 85% of controls depending on age (P < 0.05), and organ weights were also affected (P < 0.05). HMGA2^-/+ gilts and boars exhibited normal reproductive development and fertility, while HMGA2^-/- boars were sterile due to undescended testes (cryptorchidism). Crossbreeding HMGA2^-/+ boars and gilts produced litters lacking the HMGA2^-/- genotype. However, analysis of day (D) D40 and D78 pregnancies indicated that HMGA2^-/- fetuses were present at the expected Mendelian ratio, but placental abnormalities were seen in the D78 HMGA2^-/- concepti. Additionally, HMGA2^-/- embryos generated by gene editing and SCNT produced multiple pregnancies and viable offspring, indicating that lack of HMGA2 is not lethal per se. Overall, our results show that the effect of HMGA2 with respect to growth regulation is highly conserved among mammals and opens up the possibility of regulating body and organ size in a variety of mammalian species including food and companion animals.

Detecting selection signatures between Duroc and Duroc synthetic pig populations using high-density SNP chip

The development of high throughput genotyping techniques has facilitated the identification of selection signatures of pigs. The detection of genomic selection signals in a population subjected to differential selection pressures may provide insights into the genes associated with economically and biologically important traits. To identify genomic regions under selection, we genotyped 488 Duroc (D) pigs and 155 D × Korean native pigs (DKNPs) using the Porcine SNP70K BeadChip. By applying the F_ST and extended haplotype homozygosity (EHH-Rsb) methods, we detected genes under directional selection associated with growth/stature (DOCK7, PLCB4, HS2ST1, FBP2 and TG), carcass and meat quality (TG, COL14A1, FBXO5, NR3C1, SNX7, ARHGAP26 and DPYD), number of teats (LOC100153159 and LRRC1), pigmentation (MME) and ear morphology (SOX5), which are all mostly near or at fixation. These results could be a basis for investigating the underlying mutations associated with observed phenotypic variation. Validation using genome-wide association analysis would also facilitate the inclusion of some of these markers in genetic evaluation programs.

Porcine methionine sulfoxide reductase B3: molecular cloning, tissue-specific expression profiles, and polymorphisms associated with ear size in Sus scrofa

Background

In Sus scrofa, methionine sulfoxide reductase B3 (MSRB3) is a crucial candidate gene for ear size, and an important conformational trait of pig breeds. However, challenges in MSRB3 cDNA amplification have prevented further identification of MSRB3 allelic variants influencing pig ear size.

Results

We cloned a full-length cDNA sequence of porcine MSRB3 by rapid-amplification of cDNA ends. The 3,765-bp gene contained a 5’-untranslated region (UTR) (190 bp), a coding region (552 bp), and a 3’-UTR (3,016 bp) and shared 84 %, 84 %, 87 %, 86 %, and 70 % sequence identities with human, orangutan, mouse, chicken, and zebrafish, respectively. The gene encoded a 183-amino acid protein, which shared 88 %, 91 %, 89 %, 86 %, and 67 % identities with human, orangutan, mouse, chicken, and zebrafish, respectively. Tissue expression analysis using qRT-PCR revealed that MSRB3 was expressed in the heart, liver, lung, kidney, spleen, ear, muscle, fat, lymph, skeletal, and hypothalamic tissues. Three single nucleotide polymorphisms (SNPs) were identified in MSRB3: c.-735C > T in the 5’ flanking region, c.2571 T > C in the 3’-UTR, and a synonymous mutation of c.484 T > C in the coding region. The SNPs c.-735C > T and c.2571 T > C were significantly associated with ear size in a Large White × Minzhu F2 population other than in Beijing Black pigs. Subsequently, at SNP c.-735C > T, the mRNA of MSRB3 was significantly higher expressed in ears of individuals with the TT genotype (Minzhu) than those with CC (Large White).

Conclusions

The porcine MSRB3 owned a 3,765-bp full-length cDNA sequence and was detected to express in ear tissue. Two SNPs of this gene were shown to be significantly associated with ear size in a Large White × Minzhu intercross population instead of Beijing Black pig population. What’s more, the individuals with higher mRNA expression of MSRB3 have larger ear sizes. These results provide useful information for further functional analyses of MSRB3 influencing ear size in pigs.

Electronic supplementary material

The online version of this article (doi:10.1186/s40104-015-0060-x) contains supplementary material, which is available to authorized users.

Genome-wide scan reveals LEMD3 and WIF1 on SSC5 as the candidates for porcine ear size

The quantitative trait loci (QTL) for porcine ear size was previously reported to mainly focus on SSC5 and SSC7. Recently, a missense mutation, G32E, in PPARD in the QTL interval on SSC7 was identified as the causative mutation for ear size. However, on account of the large interval of QTL, the responsible gene on SSC5 has not been identified. In this study, an intercross population was constructed from the large-eared Minzhu, an indigenous Chinese pig breed, and the Western commercial Large White pig to examine the genetic basis of ear size diversity. A GWAS was performed to detect SNPs significantly associated with ear size. Thirty-five significant SNPs defined a 10.78-Mb (30.14–40.92 Mb) region on SSC5. Further, combining linkage disequilibrium and haplotype sharing analysis, a reduced region of 3.07-Mb was obtained. Finally, by using a selective sweep analysis, a critical region of about 450-kb interval containing two annotated genes LEMD3 and WIF1 was refined in this work. Functional analysis indicated that both represent biological candidates for porcine ear size, with potential application in breeding programs. The two genes could also be used as novel references for further study of the mechanism underlying human microtia.

Application of selection mapping to identify genomic regions associated with dairy production in sheep

In Europe, especially in Mediterranean areas, the sheep has been traditionally exploited as a dual purpose species, with income from both meat and milk. Modernization of husbandry methods and the establishment of breeding schemes focused on milk production have led to the development of "dairy breeds." This study investigated selective sweeps specifically related to dairy production in sheep by searching for regions commonly identified in different European dairy breeds. With this aim, genotypes from 44,545 SNP markers covering the sheep autosomes were analysed in both European dairy and non-dairy sheep breeds using two approaches: (i) identification of genomic regions showing extreme genetic differentiation between each dairy breed and a closely related non-dairy breed, and (ii) identification of regions with reduced variation (heterozygosity) in the dairy breeds using two methods. Regions detected in at least two breeds (breed pairs) by the two approaches (genetic differentiation and at least one of the heterozygosity-based analyses) were labeled as core candidate convergence regions and further investigated for candidate genes. Following this approach six regions were detected. For some of them, strong candidate genes have been proposed (e.g. ABCG2, SPP1), whereas some other genes designated as candidates based on their association with sheep and cattle dairy traits (e.g. LALBA, DGAT1A) were not associated with a detectable sweep signal. Few of the identified regions were coincident with QTL previously reported in sheep, although many of them corresponded to orthologous regions in cattle where QTL for dairy traits have been identified. Due to the limited number of QTL studies reported in sheep compared with cattle, the results illustrate the potential value of selection mapping to identify genomic regions associated with dairy traits in sheep.

Signatures of diversifying selection in European pig breeds

Following domestication, livestock breeds have experienced intense selection pressures for the development of desirable traits. This has resulted in a large diversity of breeds that display variation in many phenotypic traits, such as coat colour, muscle composition, early maturity, growth rate, body size, reproduction, and behaviour. To better understand the relationship between genomic composition and phenotypic diversity arising from breed development, the genomes of 13 traditional and commercial European pig breeds were scanned for signatures of diversifying selection using the Porcine60K SNP chip, applying a between-population (differentiation) approach. Signatures of diversifying selection between breeds were found in genomic regions associated with traits related to breed standard criteria, such as coat colour and ear morphology. Amino acid differences in the EDNRB gene appear to be associated with one of these signatures, and variation in the KITLG gene may be associated with another. Other selection signals were found in genomic regions including QTLs and genes associated with production traits such as reproduction, growth, and fat deposition. Some selection signatures were associated with regions showing evidence of introgression from Asian breeds. When the European breeds were compared with wild boar, genomic regions with high levels of differentiation harboured genes related to bone formation, growth, and fat deposition.

The domestic pig, an important source of protein worldwide, was domesticated from the ancestral wild boar in multiple locations throughout the world. In Europe, local types were developed following domestication, but phenotypically distinct breeds only arose in the eighteenth century with the advent of systematic breeding. Recently developed molecular tools for pigs (as well as other livestock species) now allow a genetic characterisation of breed histories, including identification of regions of the genome that have been under selection in the establishment of breeds. We have applied these tools to identify genomic regions associated with breed development in a set of commercial and traditional pig breeds. We found strong evidence of genetic differentiation between breeds near genes associated with traits that are used to define breed standards, such as ear morphology and coat colour, as well as in regions of the genome that are associated with pork production traits. It is well documented that crosses with Asian pigs have been used to modify European breeds. We have found evidence of genetic influence from Asian pigs in European breeds, again in regions of the genome associated with breed standard characteristics, including ear shape and coat colour, as well as production traits.