Mapping, fine mapping, and molecular dissection of quantitative trait Loci in domestic animals

Engineering multi-degrading bacterial communities to bioremediate soils contaminated with pesticides residues

Parallel to the important use of pesticides in conventional agriculture there is a growing interest for green technologies to clear contaminated soil from pesticides and their degradation products. Bioaugmentation i. e. the inoculation of degrading micro-organisms in polluted soil, is a promising method still in needs of further developments. Specifically, improvements in the understanding of how degrading microorganisms must overcome abiotic filters and interact with the autochthonous microbial communities are needed in order to efficiently design bioremediation strategies. Here we designed a protocol aiming at studying the degradation of two herbicides, glyphosate (GLY) and isoproturon (IPU), via experimental modifications of two source bacterial communities. We used statistical methods stemming from genomic prediction to link community composition to herbicides degradation potentials. Our approach proved to be efficient with correlation estimates over 0.8 - between model predictions and measured pesticide degradation values. Multi-degrading bacterial communities were obtained by coalescing bacterial communities with high GLY or IPU degradation ability based on their community-level properties. Finally, we evaluated the efficiency of constructed multi-degrading communities to remove pesticide contamination in a different soil. While results are less clear in the case of GLY, we showed an efficient transfer of degrading capacities towards the receiving soil even at relatively low inoculation levels in the case of IPU. Altogether, we developed an innovative protocol for building multi-degrading simplified bacterial communities with the help of genomic prediction tools and coalescence, and proved their efficiency in a contaminated soil.

Selection on the promoter regions plays an important role in complex traits during duck domestication

BACKGROUND

Identifying the key factors that underlie complex traits during domestication is a great challenge for evolutionary and biological studies. In addition to the protein-coding region differences caused by variants, a large number of variants are located in the noncoding regions containing multiple types of regulatory elements. However, the roles of accumulated variants in gene regulatory elements during duck domestication and economic trait improvement are poorly understood.

RESULTS

We constructed a genomics, transcriptomics, and epigenomics map of the duck genome and assessed the evolutionary forces that have been in play across the whole genome during domestication. In total, 304 (42.94%) gene promoters have been specifically selected in Pekin duck among all selected genes. Joint multi-omics analysis reveals that 218 genes (72.01%) with selected promoters are located in open and active chromatin, and 267 genes (87.83%) with selected promoters were highly and differentially expressed in domestic trait-related tissues. One important candidate gene ELOVL3, with a strong signature of differentiation on the core promoter region, is known to regulate fatty acid elongation. Functional experiments showed that the nearly fixed variants in the top selected ELOVL3 promoter in Pekin duck decreased binding ability with HLF and increased gene expression, with the overexpression of ELOVL3 able to increase lipid deposition and unsaturated fatty acid enrichment.

CONCLUSIONS

This study presents genome resequencing, RNA-Seq, Hi-C, and ATAC-Seq data of mallard and Pekin duck, showing that selection of the gene promoter region plays an important role in gene expression and phenotypic changes during domestication and highlights that the variants of the ELOVL3 promoter may have multiple effects on fat and long-chain fatty acid content in ducks.

Genetic structure and characteristics of Tibetan chickens

Tibetan chicken is one of the most common and widely distributed highland breeds, and is often used as a model organism for understanding genetic adaptation to extreme environments in Tibet. Despite its apparent geographical diversity and large variations in plumage patterns, the genetic differences within breed were not accounted for in most studies and have not been systematically investigated. In order to reveal and genetically differentiate the current existing TBC sub-populations that might have major implications for genomic research in TBCs, we systematically evaluated the population structure and demography of current TBC populations. Based on 344 whole-genome sequenced birds including 115 Tibetan chickens that were mostly sampled from family-farms across Tibet, we revealed a clear separation of Tibetan chickens into 4 sub-populations that broadly aligns with their geographical distribution. Moreover, population structure, population size dynamics, and the extent of admixture jointly suggest complex demographic histories of these sub-populations, including possible multiple origins, inbreeding, and introgressions. While most of the candidate selected regions found between the TBC sub-populations and Red Jungle fowls were nonoverlapping, 2 genes RYR2 and CAMK2D were revealed as strong selection candidates in all 4 sub-populations. These 2 previously identified high altitude associated genes indicated that the sub-populations responded to similar selection pressures in an independent but functionally similar fashion. Our results demonstrate robust population structure in Tibetan chickens that will help inform future genetic analyses on chickens and other domestic animals alike in Tibet, recommending thoughtful experimental design.

Genome-Wide Association Study Identifies the Crucial Candidate Genes for Teat Number in Crossbred Commercial Pigs

The number of teats is a crucial reproductive trait with significant economic implications on maternal capacity and litter size. Consequently, improving this trait is essential to facilitate genetic selection for increased litter size. In this study, we performed a genome-wide association study (GWAS) of the number of teats in a three-way crossbred commercial Duroc × (Landrace × Yorkshire) (DLY) pig population comprising 1518 animals genotyped with the 50K BeadChip. Our analysis identified crucial quantitative trait loci (QTL) for the number of teats, containing the ABCD4 and VRTN genes on porcine chromosome 7. Our results establish SNP variants of ABCD4 and VRTN as new molecular markers for improving the number of teats in DLY pigs. Furthermore, the most significant noteworthy single nucleotide polymorphism (SNP) (7_97568284) was identified within the ABCD4 gene, exhibiting a significant association with the total teat number traits. This SNP accounted for a substantial proportion of the genetic variance, explaining 6.64% of the observed variation. These findings reveal a novel gene on SSC7 for the number of teats and provide a deeper understanding of the genetic mechanisms underlying reproductive traits.

Identification of candidate genes related to highland adaptation from multiple Chinese local chicken breeds by whole genome sequencing analysis

Understanding the genetic mechanism of highland adaptation is of great importance for breeding improvement of Tibetan chickens (TBC). Some studies of TBC have identified some candidate genes and pathways from multiple subgroups, but the related genetic mechanisms remain largely unknown. Different genetic backgrounds and the independent genetic basis of highland adaptation make it difficult to identity the selective region of highland adaptation with all TBC samples. In this study, we conducted pre-analysis in a large-scale population to select a TBC subgroup with the purest and highest level of highland-specific lineage for the further analysis. Finally, the 37 samples from a TBC subgroup and 19 Lahsa White chickens were used to represent the highland group for further analysis with 80 samples from five Chinese local lowland breeds as controls. Population structure analysis revealed that highland adaptation significantly affected population stratification in Chinese local chicken breeds. Genome-wide selection signal analysis identified 201 candidate genes associated with highland adaptation of TBC, and these genes were significantly enriched in calcium signaling, vascular smooth muscle contraction and the cellular response to oxidative stress pathways. Additionally, we identified a narrow 1.76 kb region containing an overlapping region between HBZ and an active enhancer, and our identified region showed a highly significant signal. The highland group selected the haplotype with high activity to improve the oxygen-carrying capacity, thus being adapted to a hypoxic environment. We also found that STX2 was significantly selected in the highland group, thus potentially reducing the oxidative stress caused by hypoxia, and that STX2 exhibited the opposite effects on highland adaptation and reproductive traits. Our findings advance our understanding of extreme environment adaptation of highland chickens, and provide some variants and genes beneficial to TBC genetic breeding improvement.

Genome-wide association study reveals candidate genes relevant to body weight in female turkeys (Meleagris gallopavo)

The underlying genetic mechanisms affecting turkey growth traits have not been widely investigated. Genome-wide association studies (GWAS) is a powerful approach to identify candidate regions associated with complex phenotypes and diseases in livestock. In the present study, we performed GWAS to identify regions associated with 18-week body weight in a turkey population. The data included body weight observations for 24,989 female turkeys genotyped based on a 65K SNP panel. The analysis was carried out using a univariate mixed linear model with hatch-week-year and the 2 top principal components fitted as fixed effects and the accumulated polygenic effect of all markers captured by the genomic relationship matrix as random. Thirty-three significant markers were observed on 1, 2, 3, 4, 7 and 12 chromosomes, while 26 showed strong linkage disequilibrium extending up to 410 kb. These significant markers were mapped to 37 genes, of which 13 were novel. Interestingly, many of the investigated genes are known to be involved in growth and body weight. For instance, genes AKR1D1, PARP12, BOC, NCOA1, ADCY3 and CHCHD7 regulate growth, body weight, metabolism, digestion, bile acid biosynthetic and development of muscle cells. In summary, the results of our study revealed novel candidate genomic regions and candidate genes that could be managed within a turkey breeding program and adapted in fine mapping of quantitative trait loci to enhance genetic improvement in this species.

A 12 kb multi-allelic copy number variation encompassing a GC gene enhancer is associated with mastitis resistance in dairy cattle

Clinical mastitis (CM) is an inflammatory disease occurring in the mammary glands of lactating cows. CM is under genetic control, and a prominent CM resistance QTL located on chromosome 6 was reported in various dairy cattle breeds. Nevertheless, the biological mechanism underpinning this QTL has been lacking. Herein, we mapped, fine-mapped, and discovered the putative causal variant underlying this CM resistance QTL in the Dutch dairy cattle population. We identified a ~12 kb multi-allelic copy number variant (CNV), that is in perfect linkage disequilibrium with a lead SNP, as a promising candidate variant. By implementing a fine-mapping and through expression QTL mapping, we showed that the group-specific component gene (GC), a gene encoding a vitamin D binding protein, is an excellent candidate causal gene for the QTL. The multiplicated alleles are associated with increased GC expression and low CM resistance. Ample evidence from functional genomics data supports the presence of an enhancer within this CNV, which would exert cis-regulatory effect on GC. We observed that strong positive selection swept the region near the CNV, and haplotypes associated with the multiplicated allele were strongly selected for. Moreover, the multiplicated allele showed pleiotropic effects for increased milk yield and reduced fertility, hinting that a shared underlying biology for these effects may revolve around the vitamin D pathway. These findings together suggest a putative causal variant of a CM resistance QTL, where a cis-regulatory element located within a CNV can alter gene expression and affect multiple economically important traits.

Genome-wide association study of body size traits in Wenshang Barred chickens based on the specific-locus amplified fragment sequencing technology

Chicken body size (BS) is an economically important trait, which has been assessed in many studies for genetic selection. However, previous reports detected functional chromosome mutations or regions using gene chips. The present study used the specific-locus amplified fragment sequencing (SLAF-seq) technology to perform a genome-wide association study (GWAS) of purebred Wenshang Barred chickens. A total of 250 one-day-old male chickens were assessed in this study. Body size in individual birds was measured at 56 days. SLAF-seq was used to genotype and GWAS analysis was carried out using the general linear model (GLM) of the TASSEL program. A total of 1,286,715 single-nucleotide polymorphisms (SNPs) were detected, of which 175,211 were tested as candidate SNPs for genome-wide association analysis using the TASSEL general linear model. Three SNPs markers reached genome-wide significance. Of these, chrZ:81729634, chrZ:81841715, and chrZ:81954149 at 81,729,634, 81,841,715, and 81,954,149 bp of GGA Z were significantly associated with body diagonal length at 56 days (BDL56); and tibia length at 56 days (TL56). These SNPs were close to three genes, including ZCCHC7, PAX5, and MELK. These results open new horizons for studies on BS and should promote the use of Chinese chickens, especially Wenshang Barred chickens.

Analysis of four complete linkage sequence variants within a novel lncRNA located in a growth QTL on chromosome 1 related to growth traits in chickens

An increasing number of studies have shown that quantitative trait loci (QTLs) at the end of chromosome 1 identified in different chicken breeds and populations exert significant effects on growth traits in chickens. Nevertheless, the causal genes underlying the QTL effect remain poorly understood. Using an updated gene database, a novel lncRNA (named LncFAM) was found at the end of chromosome 1 and located in a growth and digestion QTL. This study showed that the expression level of LncFAM in pancreas tissues with a high weight was significantly higher than that in pancreas tissues with a low weight, which indicates that the expression level of LncFAM was positively correlated with various growth phenotype indexes, such as growth speed and body weight. A polymorphism screening identified four polymorphisms with strong linkage disequilibrium in LncFAM: a 5-bp indel in the second exon, an A/G base mutation, and 7-bp and 97-bp indels in the second intron. A study of a 97-bp insertion in the second intron using an F2 chicken resource population produced by Anka and Gushi chickens showed that the mutant individuals with genotype II had the highest values for body weight (BW) at 0 days and 2, 4, 6, 8, 10 and 12 weeks, shank girth (SG) at 4, 8 and 12 weeks, chest width (CW) at 4, 8 and 12 weeks, body slant length (BSL) at 8 and 12 weeks, and pelvic width (PW) at 4, 8 and 12 weeks, followed by ID and DD genotypes. The amplification and typing of 2,716 chickens from ten different breeds, namely, the F2 chicken resource population, dual-type chickens, including Xichuan black-bone chickens, Lushi green-shell layers, Dongxiang green-shell layers, Changshun green-shell layers, and Gushi chickens, and commercial broilers, including Ross 308, AA, Cobb and Hubbard broilers, revealed that II was the dominant genotype. Interestingly, only genotype II existed among the tested populations of commercial broilers. Moreover, the expression level in the pancreas tissue of Ross 308 chickens was significantly higher than that in the pancreas tissue of Gushi chickens (P < 0.001), which might be related to the conversion rates among different chickens. The prediction and verification of the target gene of LncFAM showed that LncFAM might regulate the expression of its target gene FAM48A through cis-expression. Our results provide useful information on the mutation of LncFAM, which can be used as a potential molecular breeding marker.

Population genomics of honey bees reveals a selection signature indispensable for royal jelly production

In order to interpret the molecular mechanisms that modulating the organism variations and selection signatures to drive adaptive evolutionary changes are indispensable goals in the new evolutionary ecological genetics. Here, we identified the gene locus associated to royal jelly production through whole-genome sequencing of the DNA from eight populations of honeybees. The analysis of the samples was composed of 120 individuals and each pointed extremely opposite trait values for a given phenotype. We identified functional single nucleotide polymorphisms (SNPs) candidate that might be essential in regulating the phenotypic traits of honeybee populations. Moreover, selection signatures were investigated using pooling sequencing of eight distinct honeybee populations, and the results provided the evidence of signatures of recent selection among populations under different selection objectives. Furthermore, gene ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses indicated that selected genes were potentially involved in several biological processes and molecular functioning, which could directly or indirectly influence the production of royal jelly. Our findings can be used to understand the genomic signatures, as well as implicate a profound glance on genomic regions that control the production trait of royal jelly in honey bees.

Unraveling genomic associations with feed efficiency and body weight traits in chickens through an integrative approach

BACKGROUND

Feed efficiency and growth rate have been targets for selection to improve chicken production. The incorporation of genomic tools may help to accelerate selection. We genotyped 529 individuals using a high-density SNP chip (600 K, Affymetrix®) to estimate genomic heritability of performance traits and to identify genomic regions and their positional candidate genes associated with performance traits in a Brazilian F2 Chicken Resource population. Regions exhibiting selection signatures and a SNP dataset from resequencing were integrated with the genomic regions identified using the chip to refine the list of positional candidate genes and identify potential causative mutations.

RESULTS

Feed intake (FI), feed conversion ratio (FC), feed efficiency (FE) and weight gain (WG) exhibited low genomic heritability values (i.e. from 0.0002 to 0.13), while body weight at hatch (BW1), 35 days-of-age (BW35), and 41 days-of-age (BW41) exhibited high genomic heritability values (i.e. from 0.60 to 0.73) in this F2 population. Twenty unique 1-Mb genomic windows were associated with BW1, BW35 or BW41, located on GGA1-4, 6-7, 10, 14, 24, 27 and 28. Thirty-eight positional candidate genes were identified within these windows, and three of them overlapped with selection signature regions. Thirteen predicted deleterious and three high impact sequence SNPs in these QTL regions were annotated in 11 positional candidate genes related to osteogenesis, skeletal muscle development, growth, energy metabolism and lipid metabolism, which may be associated with body weight in chickens.

CONCLUSIONS

The use of a high-density SNP array to identify QTL which were integrated with whole genome sequence signatures of selection allowed the identification of candidate genes and candidate causal variants. One novel QTL was detected providing additional information to understand the genetic architecture of body weight traits. We identified QTL for body weight traits, which were also associated with fatness in the same population. Our findings form a basis for further functional studies to elucidate the role of specific genes in regulating body weight and fat deposition in chickens, generating useful information for poultry breeding programs.

A functional regulatory variant of MYH3 influences muscle fiber-type composition and intramuscular fat content in pigs

Muscle development and lipid accumulation in muscle critically affect meat quality of livestock. However, the genetic factors underlying myofiber-type specification and intramuscular fat (IMF) accumulation remain to be elucidated. Using two independent intercrosses between Western commercial breeds and Korean native pigs (KNPs) and a joint linkage-linkage disequilibrium analysis, we identified a 488.1-kb region on porcine chromosome 12 that affects both reddish meat color (a*) and IMF. In this critical region, only the MYH3 gene, encoding myosin heavy chain 3, was found to be preferentially overexpressed in the skeletal muscle of KNPs. Subsequently, MYH3-transgenic mice demonstrated that this gene controls both myofiber-type specification and adipogenesis in skeletal muscle. We discovered a structural variant in the promotor/regulatory region of MYH3 for which Q allele carriers exhibited significantly higher values of a* and IMF than q allele carriers. Furthermore, chromatin immunoprecipitation and cotransfection assays showed that the structural variant in the 5'-flanking region of MYH3 abrogated the binding of the myogenic regulatory factors (MYF5, MYOD, MYOG, and MRF4). The allele distribution of MYH3 among pig populations worldwide indicated that the MYH3 Q allele is of Asian origin and likely predates domestication. In conclusion, we identified a functional regulatory sequence variant in porcine MYH3 that provides novel insights into the genetic basis of the regulation of myofiber type ratios and associated changes in IMF in pigs. The MYH3 variant can play an important role in improving pork quality in current breeding programs.

Band-collision gel electrophoresis

Electrophoretic mobility shift assays are widely used in gel electrophoresis to study binding interactions between different molecular species loaded into the same well. However, shift assays can access only a subset of reaction possibilities that could be otherwise seen if separate bands of reagent species might instead be collisionally reacted. Here, we adapt gel electrophoresis by fabricating two or more wells in the same lane, loading these wells with different reagent species, and applying an electric field, thereby producing collisional reactions between propagating pulse-like bands of these species, which we image optically. For certain pairs of anionic and cationic dyes, propagating bands pass through each other unperturbed; yet, for other pairs, we observe complexing and precipitation reactions, indicating strong attractive interactions. We generalize this band-collision gel electrophoresis (BCGE) approach to other reaction types, including acid-base, ligand exchange, and redox, as well as to colloidal species in passivated large-pore gels.

Electrophoretic mobility shift assays are widely used in gel electrophoresis to study binding interactions between different molecular species, but these assays access only a subset of reaction possibilities. Here, the authors develop a band-collision gel electrophoresis (BCGE) approach that demonstrates a much wider variety of reaction types.

Integration of genome wide association studies and whole genome sequencing provides novel insights into fat deposition in chicken

Excessive fat deposition is a negative factor for poultry production because it reduces feed efficiency, increases the cost of meat production and is a health concern for consumers. We genotyped 497 birds from a Brazilian F₂ Chicken Resource Population, using a high-density SNP array (600 K), to estimate the genomic heritability of fat deposition related traits and to identify genomic regions and positional candidate genes (PCGs) associated with these traits. Selection signature regions, haplotype blocks and SNP data from a previous whole genome sequencing study in the founders of this chicken F2 population were used to refine the list of PCGs and to identify potential causative SNPs. We obtained high genomic heritabilities (0.43-0.56) and identified 22 unique QTLs for abdominal fat and carcass fat content traits. These QTLs harbored 26 PCGs involved in biological processes such as fat cell differentiation, insulin and triglyceride levels, and lipid biosynthetic process. Three of these 26 PCGs were located within haplotype blocks there were associated with fat traits, five overlapped with selection signature regions, and 12 contained predicted deleterious variants. The identified QTLs, PCGs and potentially causative SNPs provide new insights into the genetic control of fat deposition and can lead to improved accuracy of selection to reduce excessive fat deposition in chickens.

Genome-wide association study of body weight in Wenshang Barred chicken based on the SLAF-seq technology

Chicken body weight (BW) is an economically important trait, and many studies have been conducted on genetic selection for BW. However, previous studies have detected functional chromosome mutations or regions using gene chips. The present study used the specific-locus amplified fragment sequencing (SLAF-seq) technology to perform a genome-wide association study (GWAS) on purebred Wengshang Barred chicken. A total of 1,286,715 single-nucleotide polymorphisms (SNPs) were detected, and 175,211 SNPs were selected as candidate SNPs for genome-wide association analysis using TASSEL general linear models. Six SNP markers reached genome-wide significance. Of these, rs732048524, rs735522839, rs738991545, and rs15837818 were significantly associated with body weight at 28 days (BW28), while rs314086457 and rs315694878 were significantly associated with BW120. These SNPs are close to seven genes (PRSS23, ME3, FAM181B, NABP1, SDPR, TSSK6L2, and RBBP8). Moreover, 24 BW-associated SNPs reached "suggestive" genome-wide significance. Of these, 6, 13, 1, and 4 SNPs were associated with BW28, BW56, BW80, and BW120, respectively. These results would enrich the studies on BW and promote the use of Chinese chicken, especially the Wenshang Barred chicken.

Invited review: A perspective on the future of genomic selection in dairy cattle

Genomic evaluation has been successfully implemented in the United States, Canada, Great Britain, Ireland, New Zealand, Australia, France, the Netherlands, Germany, and the Scandinavian countries. Adoption of this technology in the major dairy producing countries has led to significant changes in the worldwide dairy industry. Gradual elimination of the progeny test system has led to a reduction in the number of sires with daughter records and fewer genetic ties between years. As genotyping costs decrease, the number of cows genotyped will continue to increase, and these records will become the basic data used to compute genomic evaluations, most likely via application of "single-step" methodologies. Although genomic selection has been successful in increasing rates of genetic gain, we still know very little about the genetic architecture of quantitative variation. Apparently, a very large number of genes affect nearly all economic traits, in accordance with the infinitesimal model for quantitative traits. Less emphasis in selection goals will be placed on milk production traits, and more on health, reproduction, and efficiency traits and on environmentally friendly production with reduced waste and gas emission. Genetic variance for economic traits is maintained by the increase in frequency of rare alleles, new mutations, and changes in selection goals and management. Thus, it is unlikely that a selection plateau will be reached in the near future.

NGS-based reverse genetic screen for common embryonic lethal mutations compromising fertility in livestock

We herein report the result of a large-scale, next generation sequencing (NGS)-based screen for embryonic lethal (EL) mutations in Belgian beef and New Zealand dairy cattle. We estimated by simulation that cattle might carry, on average, ∼0.5 recessive EL mutations. We mined exome sequence data from >600 animals, and identified 1377 stop-gain, 3139 frame-shift, 1341 splice-site, 22,939 disruptive missense, 62,399 benign missense, and 92,163 synonymous variants. We show that cattle have a comparable load of loss-of-function (LoF) variants (defined as stop-gain, frame-shift, or splice-site variants) as humans despite having a more variable exome. We genotyped >40,000 animals for up to 296 LoF and 3483 disruptive missense, breed-specific variants. We identified candidate EL mutations based on the observation of a significant depletion in homozygotes. We estimated the proportion of EL mutations at 15% of tested LoF and 6% of tested disruptive missense variants. We confirmed the EL nature of nine candidate variants by genotyping 200 carrier × carrier trios, and demonstrating the absence of homozygous offspring. The nine identified EL mutations segregate at frequencies ranging from 1.2% to 6.6% in the studied populations and collectively account for the mortality of ∼0.6% of conceptuses. We show that EL mutations preferentially affect gene products fulfilling basic cellular functions. The resulting information will be useful to avoid at-risk matings, thereby improving fertility.

Fine mapping QTL for resistance to VNN disease using a high-density linkage map in Asian seabass

Asian seabass has suffered from viral nervous necrosis (VNN) disease. Our previous study has mapped quantitative trait loci (QTL) for resistance to VNN disease. To fine map these QTL and identify causative genes, we identified 6425 single nucleotide polymorphisms (SNPs) from 85 dead and 94 surviving individuals. Combined with 155 microsatellites, we constructed a genetic map consisting of 24 linkage groups (LGs) containing 3000 markers, with an average interval of 1.27 cM. We mapped one significant and three suggestive QTL with phenotypic variation explained (PVE) of 8.3 to 11.0%, two significant and two suggestive QTL with PVE of 7.8 to 10.9%, for resistance in three LGs and survival time in four LGs, respectively. Further analysis one QTL with the largest effect identified protocadherin alpha-C 2-like (Pcdhac2) as the possible candidate gene. Association study in 43 families with 1127 individuals revealed a 6 bp insertion-deletion was significantly associated with disease resistance. qRT-PCR showed the expression of Pcdhac2 was significantly induced in the brain, muscle and skin after nervous necrosis virus (NNV) infection. Our results could facilitate marker-assisted selection (MAS) for resistance to NNV in Asian seabass and set up the basis for functional analysis of the potential causative gene for resistance.

Variations at regulatory regions of the milk protein genes are associated with milk traits and coagulation properties in the Sarda sheep

Regulatory variation at the ovine casein genes could have important effects on the composition and coagulation properties of milk. Herewith, we have partially resequenced the promoters and the 3'-UTR of the four casein genes in 25 Sarda sheep. Alignment of these sequences allowed us to identify a total of 29 SNPs. This level of polymorphism (one SNP every 250 bp) is remarkably high if compared with SNP densities estimated in human genic regions (approximately one SNP per bp). The 29 SNPs identified in our resequencing experiment, plus three previously reported SNPs mapping to the lactalbumin, alpha (LALBA) and β-lactoglobulin (BLG, also known as progestagen-associated endometrial protein, PAEP) genes, were genotyped with a multiplex TaqMan Open Array Real-Time PCR assay in 760 Sarda sheep with records for milk composition and coagulation properties. Association analysis revealed the existence of significant associations of CSN1S2 and CSN3 genotypes with milk protein and casein contents. Moreover, genotypes at CSN1S1 were significantly associated with rennet coagulation time, curd firming time and curd firmness, whereas CSN2 was associated with curd firming time. These results suggest that SNPs mapping to the promoters and 3'-UTRs of ovine casein genes may exert regulatory effects on gene expression and that they could be used for improving sheep milk quality and technological traits at the population level through marker assisted selection.

Whole-genome resequencing of honeybee drones to detect genomic selection in a population managed for royal jelly

Four main evolutionary lineages of A. mellifera have been described including eastern Europe (C) and western and northern Europe (M). Many apiculturists prefer bees from the C lineage due to their docility and high productivity. In France, the routine importation of bees from the C lineage has resulted in the widespread admixture of bees from the M lineage. The haplodiploid nature of the honeybee Apis mellifera, and its small genome size, permits affordable and extensive genomics studies. As a pilot study of a larger project to characterise French honeybee populations, we sequenced 60 drones sampled from two commercial populations managed for the production of honey and royal jelly. Results indicate a C lineage origin, whilst mitochondrial analysis suggests two drones originated from the O lineage. Analysis of heterozygous SNPs identified potential copy number variants near to genes encoding odorant binding proteins and several cytochrome P450 genes. Signatures of selection were detected using the hapFLK haplotype-based method, revealing several regions under putative selection for royal jelly production. The framework developed during this study will be applied to a broader sampling regime, allowing the genetic diversity of French honeybees to be characterised in detail.

Combatting African Animal Trypanosomiasis (AAT) in livestock: The potential role of trypanotolerance

African Animal Trypanosomiasis (AAT) is endemic in at least 37 of the 54 countries in Africa. It is estimated to cause direct and indirect losses to the livestock production industry in excess of US$ 4.5 billion per annum. A century of intervention has yielded limited success, owing largely to the extraordinary complexity of the host-parasite interaction. Trypanotolerance, which refers to the inherent ability of some African livestock breeds, notably Djallonke sheep, N'Dama cattle and West African Dwarf goats, to withstand a trypanosomiasis challenge and still remain productive without any form of therapy, is an economically sustainable option for combatting this disease. Yet trypanotolerance has not been adequately exploited in the fight against AAT. In this review, we describe new insights into the genetic basis of trypanotolerance and discuss the potential of exploring this phenomenon as an integral part of the solution for AAT, particularly, in the context of African animal production systems.

Mapping QTL for Resistance Against Viral Nervous Necrosis Disease in Asian Seabass

Viral nervous necrosis disease (VNN), caused by nervous necrosis virus (NNV), leads to mass mortality in mariculture. However, phenotypic selection for resistance against VNN is very difficult. To facilitate marker-assisted selection (MAS) for resistance against VNN and understanding of the genetic architecture underlying the resistance against this disease, we mapped quantitative trait loci (QTL) for resistance against VNN in Asian seabass. We challenged fingerlings at 37 days post-hatching (dph), from a single back-cross family, with NNV at a concentration of 9 × 10(6) TCID50/ml for 2 h. Daily mortalities were recorded and collected. A panel of 330 mortalities and 190 surviving fingerlings was genotyped using 149 microsatellites with 145 successfully mapped markers covering 24 linkage groups (LGs). Analysis of QTL for both resistance against VNN and survival time was conducted using interval mapping. Five significant QTL located in four LGs and eight suggestive QTL in seven LGs were identified for resistance. Another five significant QTL in three LGs and five suggestive QTL in three LGs were detected for survival time. One significant QTL, spanning 3 cM in LG20, was identified for both resistance and survival time. These QTL explained 2.2-4.1% of the phenotypic variance for resistance and 2.2-3.3% of the phenotypic variance for survival time, respectively. Our results suggest that VNN resistance in Asian seabass is controlled by many loci with small effects. Our data provide information for fine mapping of QTL and identification of candidate genes for a better understanding of the mechanism of disease resistance.

Targeted Exome Sequencing of Deafness Genes After Failure of Auditory Phenotype-Driven Candidate Gene Screening

OBJECTIVE

To demonstrate the efficacy and advantages of targeted exome sequencing (TES) of known deafness genes in cases with failed or misleading auditory phenotype-driven candidate gene screening.

STUDY DESIGN

Prospective cohort survey.

SETTING

Otolaryngology department of a tertiary referral hospital.

PATIENTS

Six hearing-impaired probands with seemingly non-syndromic features from six deaf families were enrolled in this study after failure of genetic diagnosis using auditory phenotype-driven candidate gene screening.

INTERVENTION

TES of known deafness genes was performed in the six probands, and a final causative variant was pursued using subsequent filtering steps.

MAIN OUTCOME MEASURE

Potential causative variants determined using TES were confirmed by previously introduced filtering steps.

RESULTS

We detected causative variants in three (50%) of six families, and these variants were in the COCH, PAX3, and GJB2 genes. Additionally, we also recapitulated the recent finding from other report arguing for the non-pathogenic potential of MYO1A variant.

CONCLUSIONS

TES of a deafness panel provides a comprehensive genetic screening tool that can be implemented without being misled by the audiogram configuration information and can complement incomplete clinical physical examinations. In addition, the secondary incidental finding obtained by TES contributes useful information regarding the deafness field.

Genome-association analysis of Korean Holstein milk traits using genomic estimated breeding value

Objective

Holsteins are known as the world’s highest-milk producing dairy cattle. The purpose of this study was to identify genetic regions strongly associated with milk traits (milk production, fat, and protein) using Korean Holstein data.

Methods

This study was performed using single nucleotide polymorphism (SNP) chip data (Illumina BovineSNP50 Beadchip) of 911 Korean Holstein individuals. We inferred each genomic estimated breeding values based on best linear unbiased prediction (BLUP) and ridge regression using BLUPF90 and R. We then performed a genome-wide association study and identified genetic regions related to milk traits.

Results

We identified 9, 6, and 17 significant genetic regions related to milk production, fat and protein, respectively. These genes are newly reported in the genetic association with milk traits of Holstein.

Conclusion

This study complements a recent Holstein genome-wide association studies that identified other SNPs and genes as the most significant variants. These results will help to expand the knowledge of the polygenic nature of milk production in Holsteins.

Association analysis for udder health based on SNP-panel and sequence data in Danish Holsteins

BACKGROUND

The sensitivity of genome-wide association studies for the detection of quantitative trait loci (QTL) depends on the density of markers examined and the statistical models used. This study compares the performance of three marker densities to refine six previously detected QTL regions for mastitis traits: 54 k markers of a medium-density SNP (single nucleotide polymorphism) chip (MD), imputed 777 k markers of a high-density SNP chip (HD), and imputed whole-genome sequencing data (SEQ). Each dataset contained data for 4496 Danish Holstein cattle. Comparisons were performed using a linear mixed model (LM) and a Bayesian variable selection model (BVS).

RESULTS

After quality control, 587, 7825, and 78 856 SNPs in the six targeted regions remained for MD, HD, and SEQ data, respectively. In general, the association patterns between SNPs and traits were similar for the three marker densities when tested using the same statistical model. With the LM model, 120 (MD), 967 (HD), and 7209 (SEQ) SNPs were significantly associated with mastitis, whereas with the BVS model, 43 (MD), 131 (HD), and 1052 (SEQ) significant SNPs (Bayes factor > 3.2) were observed. A total of 26 (MD), 75 (HD), and 465 (SEQ) significant SNPs were identified by both models. In addition, one, 16, and 33 QTL peaks for MD, HD, and SEQ data were detected according to the QTL intensity profile of SNP bins by post-analysis of the BVS model.

CONCLUSIONS

The power to detect significant associations increased with increasing marker density. The BVS model resulted in clearer boundaries between linked QTL than the LM model. Using SEQ data, the six targeted regions were refined to 33 candidate QTL regions for udder health. The comparison between these candidate QTL regions and known genes suggested that NPFFR2, SLC4A4, DCK, LIFR, and EDN3 may be considered as candidate genes for mastitis susceptibility.

Genome-wide copy number profiling using high-density SNP array in chickens

Phenotypic diversity is a direct consequence resulting mainly from the impact of underlying genetic variation, and recent studies have shown that copy number variation (CNV) is emerging as an important contributor to both phenotypic variability and disease susceptibility. Herein, we performed a genome-wide CNV scan in 96 chickens from 12 diversified breeds, benefiting from the high-density Affymetrix 600 K SNP arrays. We identified a total of 231 autosomal CNV regions (CNVRs) encompassing 5.41 Mb of the chicken genome and corresponding to 0.59% of the autosomal sequence. The length of these CNVRs ranged from 2.6 to 586.2 kb with an average of 23.4 kb, including 130 gain, 93 loss and eight both gain and loss events. These CNVRs, especially deletions, had lower GC content and were located particularly in gene deserts. In particular, 102 CNVRs harbored 128 chicken genes, most of which were enriched in immune responses. We obtained 221 autosomal CNVRs after converting probe coordinates to Galgal3, and comparative analysis with previous studies illustrated that 153 of these CNVRs were regarded as novel events. Furthermore, qPCR assays were designed for 11 novel CNVRs, and eight (72.73%) were validated successfully. In this study, we demonstrated that the high-density 600 K SNP array can capture CNVs with higher efficiency and accuracy and highlighted the necessity of integrating multiple technologies and algorithms. Our findings provide a pioneering exploration of chicken CNVs based on a high-density SNP array, which contributes to a more comprehensive understanding of genetic variation in the chicken genome and is beneficial to unearthing potential CNVs underlying important traits of chickens.

Combined QTL and selective sweep mappings with coding SNP annotation and cis-eQTL analysis revealed PARK2 and JAG2 as new candidate genes for adiposity regulation

Very few causal genes have been identified by quantitative trait loci (QTL) mapping because of the large size of QTL, and most of them were identified thanks to functional links already known with the targeted phenotype. Here, we propose to combine selection signature detection, coding SNP annotation, and cis-expression QTL analyses to identify potential causal genes underlying QTL identified in divergent line designs. As a model, we chose experimental chicken lines divergently selected for only one trait, the abdominal fat weight, in which several QTL were previously mapped. Using new haplotype-based statistics exploiting the very high SNP density generated through whole-genome resequencing, we found 129 significant selective sweeps. Most of the QTL colocalized with at least one sweep, which markedly narrowed candidate region size. Some of those sweeps contained only one gene, therefore making them strong positional causal candidates with no presupposed function. We then focused on two of these QTL/sweeps. The absence of nonsynonymous SNPs in their coding regions strongly suggests the existence of causal mutations acting in cis on their expression, confirmed by cis-eQTL identification using either allele-specific expression or genetic mapping analyses. Additional expression analyses of those two genes in the chicken and mice contrasted for adiposity reinforces their link with this phenotype. This study shows for the first time the interest of combining selective sweeps mapping, coding SNP annotation and cis-eQTL analyses for identifying causative genes for a complex trait, in the context of divergent lines selected for this specific trait. Moreover, it highlights two genes, JAG2 and PARK2, as new potential negative and positive key regulators of adiposity in chicken and mice.

A splice mutation in the PHKG1 gene causes high glycogen content and low meat quality in pig skeletal muscle

Glycolytic potential (GP) in skeletal muscle is economically important in the pig industry because of its effect on pork processing yield. We have previously mapped a major quantitative trait loci (QTL) for GP on chromosome 3 in a White Duroc × Erhualian F₂ intercross. We herein performed a systems genetic analysis to identify the causal variant underlying the phenotype QTL (pQTL). We first conducted genome-wide association analyses in the F₂ intercross and an F₁₉ Sutai pig population. The QTL was then refined to an 180-kb interval based on the 2-LOD drop method. We then performed expression QTL (eQTL) mapping using muscle transcriptome data from 497 F₂ animals. Within the QTL interval, only one gene (PHKG1) has a cis-eQTL that was colocolizated with pQTL peaked at the same SNP. The PHKG1 gene encodes a catalytic subunit of the phosphorylase kinase (PhK), which functions in the cascade activation of glycogen breakdown. Deep sequencing of PHKG1 revealed a point mutation (C>A) in a splice acceptor site of intron 9, resulting in a 32-bp deletion in the open reading frame and generating a premature stop codon. The aberrant transcript induces nonsense-mediated decay, leading to lower protein level and weaker enzymatic activity in affected animals. The mutation causes an increase of 43% in GP and a decrease of>20% in water-holding capacity of pork. These effects were consistent across the F₂ and Sutai populations, as well as Duroc × (Landrace × Yorkshire) hybrid pigs. The unfavorable allele exists predominantly in Duroc-derived pigs. The findings provide new insights into understanding risk factors affecting glucose metabolism, and would greatly contribute to the genetic improvement of meat quality in Duroc related pigs.

Glycogen storage diseases (GSD) are a group of inherited disorders characterized by storage of excess glycogen, which are mainly caused by the abnormality of a particular enzyme essential for releasing glucose from glycogen. GSD-like conditions have been described in a wide variety of species. Pigs are a valuable model for the study of human GSD. Moreover, pigs affected by GSD usually produce inferior pork with a lower ultimate pH (so-called “acid meat”) and less processing yield due to post-mortem degradation of the excess glycogen. So far, only one causal variant, PRKAG3 R225Q, has been identified for GSD in pigs. Here we reported a loss-of-function mutation in the PHKG1 gene that causes the deficiency of the glycogen breakdown, consequently leading to GSD and acid meat in Duroc-sired pigs. Eliminating the undesirable mutation from the breeding stock by a diagnostic DNA test will greatly reduce the incidence of GSD and significantly improve pork quality and productivity in the pig.

Re-sequencing data for refining candidate genes and polymorphisms in QTL regions affecting adiposity in chicken

In this study, we propose an approach aiming at fine-mapping adiposity QTL in chicken, integrating whole genome re-sequencing data. First, two QTL regions for adiposity were identified by performing a classical linkage analysis on 1362 offspring in 11 sire families obtained by crossing two meat-type chicken lines divergently selected for abdominal fat weight. Those regions, located on chromosome 7 and 19, contained a total of 77 and 84 genes, respectively. Then, SNPs and indels in these regions were identified by re-sequencing sires. Considering issues related to polymorphism annotations for regulatory regions, we focused on the 120 and 104 polymorphisms having an impact on protein sequence, and located in coding regions of 35 and 42 genes situated in the two QTL regions. Subsequently, a filter was applied on SNPs considering their potential impact on the protein function based on conservation criteria. For the two regions, we identified 42 and 34 functional polymorphisms carried by 18 and 24 genes, and likely to deeply impact protein, including 3 coding indels and 4 nonsense SNPs. Finally, using gene functional annotation, a short list of 17 and 4 polymorphisms in 6 and 4 functional genes has been defined. Even if we cannot exclude that the causal polymorphisms may be located in regulatory regions, this strategy gives a complete overview of the candidate polymorphisms in coding regions and prioritize them on conservation- and functional-based arguments.

Identification of loci that cause phenotypic variation in diverse species with the reciprocal hemizygosity test

The reciprocal hemizygosity test is a straightforward genetic test that can positively identify genes that have evolved to contribute to a phenotypic difference between strains or between species. The test involves a comparison between hybrids that are genetically identical throughout the genome except at the test locus, which is rendered hemizygous for alternative alleles from the two parental strains. If the two reciprocal hemizygotes display different phenotypes, then the two parental alleles must have evolved. New methods for targeted mutagenesis will allow application of the reciprocal hemizygosity test in many organisms. This review discusses the principles, advantages, and limitations of the test.

Genome-wide association study for wool production traits in a Chinese Merino sheep population

Genome-wide association studies (GWAS) provide a powerful approach for identifying quantitative trait loci without prior knowledge of location or function. To identify loci associated with wool production traits, we performed a genome-wide association study on a total of 765 Chinese Merino sheep (JunKen type) genotyped with 50 K single nucleotide polymorphisms (SNPs). In the present study, five wool production traits were examined: fiber diameter, fiber diameter coefficient of variation, fineness dispersion, staple length and crimp. We detected 28 genome-wide significant SNPs for fiber diameter, fiber diameter coefficient of variation, fineness dispersion, and crimp trait in the Chinese Merino sheep. About 43% of the significant SNP markers were located within known or predicted genes, including YWHAZ, KRTCAP3, TSPEAR, PIK3R4, KIF16B, PTPN3, GPRC5A, DDX47, TCF9, TPTE2, EPHA5 and NBEA genes. Our results not only confirm the results of previous reports, but also provide a suite of novel SNP markers and candidate genes associated with wool traits. Our findings will be useful for exploring the genetic control of wool traits in sheep.

The Application of Genomic Technologies to Investigate the Inheritance of Economically Important Traits in Goats

Goat genomics has evolved at a low pace because of a lack of molecular tools and sufficient investment. Whilst thousands and hundreds of quantitative trait loci (QTL) have been identified in cattle and sheep, respectively, about nine genome scans have been performed in goats dealing with traits as conformation, growth, fiber quality, resistance to nematodes, and milk yield and composition. In contrast, a great effort has been devoted to the characterization of candidate genes and their association with milk, meat, and reproduction phenotypes. In this regard, causal mutations have been identified in the αS1-casein gene that has a strong effect on milk composition and the PIS locus that is linked to intersexuality and polledness. In recent times, the development of massive parallel sequencing technologies has allowed to build a reference genome for goats as well as to monitor the expression of mRNAs and microRNAs in a broad array of tissues and experimental conditions. Besides, the recent design of a 52K SNP chip is expected to have a broad impact in the analysis of the genetic architecture of traits of economic interest as well as in the study of the population structure of goats at a worldwide scale.

The X-Linked Hypothesis of Brain Disorders

In this article, we propose an X-linked hypothesis of brain disorders that postulates a neuronal origin of those neurodegenerative and psychiatric disorders with a greater male prevalence. The hypothesis is based on the accumulated genetics and genomic evidence linking X chromosome genes and transcripts to neuronal cells. The behavioral genetics literature has long pointed to the link between postsynaptic protein complexes coded on chromosome X and mental retardation. More recently, novel genomic evidence has emerged of X-linked mRNA overexpression of neuronal source in the human brain. We review the evidence for this hypothesis and its consistency with the distribution across genders of brain disorders of known aetiology. We then provide examples of the utilization of this hypothesis in the investigation of the pathophysiology of complex brain disorders in both the stratification of disease cohorts and the development of realistic preclinical models. We conclude by providing a general framework for testing its validity, which will be exploited in future studies, and provide future directions for research.

Porcine ubiquitin-like 5 (UBL5) gene: genomic organization, polymorphisms, mRNA cloning, splicing variants and association study

Ubiquitin-like 5 (UBL5), which is supposed to be involved in regulation of feed intake, energy metabolism, obesity and type 2 diabetes, is located at position 62.1 cM on the pig chromosome 2 region harbouring quantitative trait loci for carcass and meat quality. The 4,354 bp genomic sequence (FR798948) of the porcine gene encompassing the promoter and entire gene was cloned by polymerase chain reaction. Comparative sequencing revealed 13 polymorphisms in noncoding regions. Synthesis of full-length cDNA sequences using rapid amplification of 5' and 3' ends showed three splice variants. Variants 1 and 2 differ in transcription length for the untranslated part of exon 1 with deduced protein of 73 amino acid (aa) residues and 100 % identities between human, mouse and other species. Variant 3, with 4 bp deletion at the 3' end of exon 2, encodes a truncated protein with 28 aa residues. In a Wild boar×Meishan F2 population (n = 334) with 47 recorded traits, loci FR798948:g.2788G>A and FR798948:g.2141T>C were associated at nominal P < 0.05 with fat deposition, growth and fattening and muscling but after adjustment for multiple testing (Benjamini and Hochberg, J R Stat Soc B 57:289-300, 1995) only eight fat deposition traits showed suggestive association with FR798948:g.2788G>A at adjusted P < 0.10. In a Meishan×Large White (MLW) cross (n = 562) with six trait records available, FR798948:g.2141T>C showed suggestive association with growth (adjusted P = 0.0690). As association mapping conducted in the outbred MLW population is more precise than in the three generation F2 population the UBL5 gene tends to be associated with growth rather than with fat accretion.

Genome-wide patterns of genetic variation in two domestic chickens

Domestic chickens are excellent models for investigating the genetic basis of phenotypic diversity, as numerous phenotypic changes in physiology, morphology, and behavior in chickens have been artificially selected. Genomic study is required to study genome-wide patterns of DNA variation for dissecting the genetic basis of phenotypic traits. We sequenced the genomes of the Silkie and the Taiwanese native chicken L2 at ∼23- and 25-fold average coverage depth, respectively, using Illumina sequencing. The reads were mapped onto the chicken reference genome (including 5.1% Ns) to 92.32% genome coverage for the two breeds. Using a stringent filter, we identified ∼7.6 million single-nucleotide polymorphisms (SNPs) and 8,839 copy number variations (CNVs) in the mapped regions; 42% of the SNPs have not found in other chickens before. Among the 68,906 SNPs annotated in the chicken sequence assembly, 27,852 were nonsynonymous SNPs located in 13,537 genes. We also identified hundreds of shared and divergent structural and copy number variants in intronic and intergenic regions and in coding regions in the two breeds. Functional enrichments of identified genetic variants were discussed. Radical nsSNP-containing immunity genes were enriched in the QTL regions associated with some economic traits for both breeds. Moreover, genetic changes involved in selective sweeps were detected. From the selective sweeps identified in our two breeds, several genes associated with growth, appetite, and metabolic regulation were identified. Our study provides a framework for genetic and genomic research of domestic chickens and facilitates the domestic chicken as an avian model for genomic, biomedical, and evolutionary studies.

The mammary gland in domestic ruminants: a systems biology perspective

Milk and dairy products are central elements in the human diet. It is estimated that 108kg of milk per year are consumed per person worldwide. Therefore, dairy production represents a relevant fraction of the economies of many countries, being cattle, sheep, goat, water buffalo, and other ruminants the main species used worldwide. An adequate management of dairy farming cannot be achieved without the knowledge on the biological mechanisms behind lactation in ruminants. Thus, understanding the morphology, development and regulation of the mammary gland in health, disease and production is crucial. Presently, innovative and high-throughput technologies such as genomics, transcriptomics, proteomics and metabolomics allow a much broader and detailed knowledge on such issues. Additionally, the application of a systems biology approach to animal science is vastly growing, as new advances in one field of specialization or animal species lead to new lines of research in other areas or/and are expanded to other species. This article addresses how modern research approaches may help us understand long-known issues in mammary development, lactation biology and dairy production.

BIOLOGICAL SIGNIFICANCE

Dairy production depends upon the knowledge of the morphology and regulation of the mammary gland and lactation. High-throughput technologies allow a much broader and detailed knowledge on the biology of the mammary gland. This paper reviews the major contributions that genomics, transcriptomics, metabolomics and proteomics approaches have provided to understand the regulation of the mammary gland in health, disease and production. In the context of mammary gland "omics"-based research, the integration of results using a Systems Biology Approach is of key importance.

The G32E functional variant reduces activity of PPARD by nuclear export and post-translational modification in pigs

Peroxisome proliferator-activated receptor beta/delta (PPARD) is a crucial and multifaceted determinant of diverse biological functions including lipid metabolism, embryonic development, inflammatory response, wound healing and cancer. Recently, we proposed a novel function of porcine PPARD (sPPARD) in external ear development. A missense mutation (G32E) in an evolutionary conservative domain of sPPARD remarkably increases external ear size in pigs. Here, we investigated the underlying molecular mechanism of the causal mutation at the cellular level. Using a luciferase reporter system, we showed that the G32E substitution reduced transcription activity of sPPARD in a ligand-dependent manner. By comparison of the subcellular localization of wild-type and mutated sPPARD in both PK-15 cells and pinna cartilage-derived primary chondrocytes, we found that the G32E substitution promoted CRM-1 mediated nuclear exportation of sPPARD. With the surface plasmon resonance technology, we further revealed that the G32E substitution had negligible effect on its ligand binding affinity. Finally, we used co-immunoprecipitation and luciferase reporter assays to show that the G32E substitution greatly reduced ubiquitination level by blocking ubiquitination of the crucial A/B domain and consequently decreased transcription activity of sPPARD. Taken together, our findings strongly support that G32E is a functional variant that plays a key role in biological activity of sPPARD, which advances our understanding of the underlying mechanism of sPPARD G32E for ear size in pigs.

Identification of quantitative trait loci underlying milk traits in Spanish dairy sheep using linkage plus combined linkage disequilibrium and linkage analysis approaches

In this study, 2 procedures were used to analyze a data set from a whole-genome scan, one based on linkage analysis information and the other combing linkage disequilibrium and linkage analysis (LDLA), to determine the quantitative trait loci (QTL) influencing milk production traits in sheep. A total of 1,696 animals from 16 half-sib families were genotyped using the OvineSNP50 BeadChip (Illumina Inc., San Diego, CA) and analysis was performed using a daughter design. Moreover, the same data set has been previously investigated through a genome-wide association (GWA) analysis and a comparison of results from the 3 methods has been possible. The linkage analysis and LDLA methodologies yielded different results, although some significantly associated regions were common to both procedures. The linkage analysis detected 3 overlapping genome-wise significant QTL on sheep chromosome (OAR) 2 influencing milk yield, protein yield, and fat yield, whereas 34 genome-wise significant QTL regions were detected using the LDLA approach. The most significant QTL for protein and fat percentages was detected on OAR3, which was reported in a previous GWA analysis. Both the linkage analysis and LDLA identified many other chromosome-wise significant associations across different sheep autosomes. Additional analyses were performed on OAR2 and OAR3 to determine the possible causality of the most significant polymorphisms identified for these genetic effects by the previously reported GWA analysis. For OAR3, the analyses demonstrated additional genetic proof of the causality previously suggested by our group for a single nucleotide polymorphism located in the α-lactalbumin gene (LALBA). In summary, although the results shown here suggest that in commercial dairy populations, the LDLA method exhibits a higher efficiency to map QTL than the simple linkage analysis or linkage disequilibrium methods, we believe that comparing the 3 analysis methods is the best approach to obtain a global picture of all identifiable QTL segregating in the population at both family-based and population-based levels.

Association mapping of quantitative trait loci for carcass and meat quality traits at the central part of chromosome 2 in Italian Large White pigs

Association mapping of the central part of porcine chromosome 2 harboring QTLs for carcass and meat quality traits was performed with 17 gene-tagged SNPs located between 44.0 and 77.5 Mb on a physical map (Sscrofa10.2) in Italian Large White pigs. For the analyzed animals records of estimated breeding values for average daily gain, back fat thickness, lean cuts, ham weight, feed conversion ratio, pH1, pHu, CIE L*, CIE a*, CIE b* and drip loss were available. A significant QTL for fat deposition (adjusted P=0.0081) and pH1 (adjusted P=0.0972) to MYOD1 at position 44.4 Mb and a QTL for growth and meatiness (adjusted P=0.0238-0.0601) to UBL5 at position 68.9 Mb were mapped. These results from association mapping are much more accurate than those from linkage mapping and facilitate further search for position candidate genes and causative mutations needed for application of markers through marker assisted selection.

Genomics of complex traits

The analysis of complex genetic traits, including mapping and identification of causative genes, has long been an enigma of genetic biology, whether in the animal sciences or in medical sciences. Traits of agricultural interest and traits of medical interest are often under the influence of both environmental factors and multiple genes, each with modest contributions to the total variance in the trait. Although the number of known mutations underlying complex traits is still relatively small, advances in genomics have greatly enhanced traditional pathways to their analysis and gene mining. The candidate gene approach, linkage analysis, and association studies are all significantly more powerful with recent advances in genome mapping, sequencing, and analysis of individual variation. Avenues to gene discovery are discussed with emphasis on genome wide association studies (GWAS) and the use of single nucleotide polymorphisms (SNPs) as revealed by increasingly powerful commercially available microarrays.