Signatures of positive selection after the introduction of genomic selection in the Finnish Ayrshire population

The Finnish Ayrshire (FAY) belongs to the Nordic Red breeds and is characterized by high milk yield, high milk components, good fertility, and functional conformation. The FAY breeding program is based on genomic selection. Despite the benefits of selection on breeding values, autozygosity in the genome may increase due to selection, and increased autozygosity may cause inbreeding depression in selected traits. However, there is lack of studies concerning selection signatures in the FAY after genomic selection introduction. The aim of this study was to identify signatures of selection in FAY after the introduction of genomic selection. Genomic data included 45,834 SNPs. The genotyped animals were divided into 2 groups: animals born before genomic selection introduction (6,108 cows) and animals born after genomic selection introduction (47,361 cows). We identified the selection signatures using 3 complementary methods: 2 based on identification of selection signatures from runs of homozygosity (ROH) islands and one based on the decay of site-specific extended haplotype between populations at SNP sites (Rsb). In total, we identified 34 ROH islands on chromosomes 1, 3, 6, 8, 12-15, 17, 19, 22, and 26 in FAY animals born before genomic selection (between 1980 and 2011) and 30 ROH islands on chromosomes 1-3, 13-17, 22, and 25-26 in FAY animals born after genomic selection introduction (between 2015 and 2020). We additionally detected 22 ΔROH islands on chromosomes 2-3, 11, 13, 14, 16, 18, 20, and 25-26. Finally, a total of 31 Rsb regions on chromosomes 2, 3, 14, 18, 20, and 25 were identified. Based on the results, genomic selection has favored certain alleles and haplotypes on genomic regions related to traits relevant in the FAY breeding program: milk production, fertility, growth, beef production traits, and feed efficiency. Several genes related to these traits, e.g., PLA2G4A, MECR, CHUK, COX15, RICTOR, SHISA9, and SEMA4G overlapped or partially overlapped the observed selection signature regions. The association of genotypes within these regions and their effects on traits relevant in the FAY breeding program should be studied and genetic regions undergoing selection monitored in the FAY population.

Inheritance of feed intake-based resilience traits and their correlation with production traits in Finnish pig breeds

High resilience against diseases, changing environmental conditions, and other stress factors and the ability to efficiently recover to normal status, is becoming increasingly important in pig production. Finding new phenotypes that relate to resilience is a crucial step for improving the resilience of pigs through selection. The objective of this study was to extract resilience-related phenotypes based on fluctuations in daily feed intake (DFI, g) and time spent in feeding per day (TPD, min) and to estimate the heritability of these traits and genetic correlations with production traits (PT). Resilience-related traits with high enough heritability and with either favorable or neutral genetic correlation with PT could be used in the selection program to improve the productivity and welfare of pigs. In this study, we used data from 7,347 Finnish Yorkshire, Landrace, and crossbred pigs raised at the test station. Six pig-specific resilience-related phenotypes were extracted from the individual DFI and TPD: root mean square error (RMSE), quantile regression (QR), and coefficient of variation (CV). RMSE was calculated from the differences between the actual DFI (or TPD) and the pig-specific predicted values. QR was based on the number of days that a pig belonged to the group with the lowest 5% of pigs based on DFI (or TPD), and CV was calculated over the daily observations of DFI (or TPD). PT included average daily gain (ADG, g), backfat thickness (BF, mm), and feed conversion rate (FCR, g/g). The heritability estimates for resilience-related traits varied between 0.07 ± 0.02 (QRDFI) and 0.20 ± 0.03 (RMSETPD). The genetic correlations between resilience-related traits and PT were mostly neutral, but for example, RMSEDFI had a favorable genetic correlation with FCR and BF but an unfavorable correlation with ADG. Lastly, we observed that pigs belonging to the lowest 10% group based on their breeding value (BV) for QRTPD had a lower proportion (10% incidence) of sick days compared to the highest 10% BV group (30% incidence). Therefore, pigs exhibiting small TPD variation (related to high resilience) tend to be less susceptible to sickness than pigs with large TPD variation (related to low resilience). Given its moderate heritability, neutral genetic correlation with PT, and positive effect on health, QRTPD can be considered the most promising resilience-related trait in the Finnish production system.

Effect of carnosine synthesis precursors in the diet on jejunal metabolomic profiling and biochemical compounds in slow-growing Korat chicken

The slow-growing Korat chicken (KR) has been developed to provide an alternative breed for smallholder farmers in Thailand. Carnosine enrichment in the meat can distinguish KR from other chicken breeds. Therefore, our aim was to investigate the effect of enriched carnosine synthesis, obtained by the β-alanine and L-histidine precursor supplementation in the diet, on changes to metabolomic profiles and biochemical compounds in slow-growing KR jejunum tissue. Four hundred 21-day-old female KR chickens were divided into 4 experimental groups: a group with a basal diet, a group with a basal diet supplemented with 1.0% β-alanine, 0.5% L-histidine, and a mix of 1.0% β-alanine and 0.5% L-histidine. The feeding trial lasted 70 d. Ten randomly selected chickens from each group were slaughtered. Metabolic profiles were analyzed using proton nuclear magnetic resonance spectroscopy. In total, 28 metabolites were identified. Significant changes in the concentrations of these metabolites were detected between the groups. Partial least squares discriminant analysis was used to distinguish the metabolites between the experimental groups. Based on the discovered metabolites, 34 potential metabolic pathways showed differentiation between groups, and 8 pathways (with impact values higher than 0.05, P < 0.05, and FDR < 0.05) were affected by metabolite content. In addition, biochemical changes were monitored using synchrotron radiation-based Fourier transform infrared microspectroscopy. Supplementation of β-alanine alone in the diet increased the β-sheets and decreased the α-helix content in the amide I region, and supplementation of L-histidine alone in the diet also increased the β-sheets. Furthermore, the relationship between metabolite contents and biochemical compounds were confirmed using principal component analysis (PCA). Results from the PCA indicated that β-alanine and L-histidine precursor group was highly positively correlated with amide I, amide II, creatine, tyrosine, valine, isoleucine, and aspartate. These findings can help to understand the relationships and patterns between the spectral and metabolic processes related to carnosine synthesis.

Thigh muscle metabolic response is linked to feed efficiency and meat characteristics in slow-growing chicken

The Korat chicken (KR) is a slow-growing Thai chicken breed with relatively poor feed efficiency (FE) but very tasty meat with high protein and low fat contents, and a unique texture. To enhance the competitiveness of KR, its FE should be improved. However, selecting for FE has an unknown effect on meat characteristics. Thus, understanding the genetic basis underlying FE traits and meat characteristics is needed. In this study, 75 male KR birds were raised up to 10 wk of age. For each bird, the feed conversion ratio (FCR), residual feed intake (RFI), and physicochemical properties, flavor precursors, and biological compounds in the thigh meat were evaluated. At 10 wk of age, thigh muscle samples from 6 birds (3 with high FCR and 3 with low FCR values) were selected, and their proteomes were investigated using a label-free proteomic method. Weighted gene coexpression network analysis (WGCNA) was used to screen the key protein modules and pathways. The WGCNA results revealed that FE and meat characteristics significantly correlated with the same protein module. However, the correlation was unfavorable; improving FE may result in a decrease in meat quality through the alteration in biological processes including glycolysis/gluconeogenesis, metabolic pathway, carbon metabolism, biosynthesis of amino acids, pyruvate metabolism, and protein processing in the endoplasmic reticulum. The hub proteins of the significant module (TNNT1, TNNT3, TNNI2, TNNC2, MYLPF, MYH10, GADPH, PGK1, LDHA, and GPI) were also identified to be associated with energy metabolism, and muscle growth and development. Given that the same proteins and pathways are present in FE and meat characteristics but in opposite directions, selection practices for KR should simultaneously consider both trait groups to maintain the high meat quality of KR while improving FE.

Estimating inbreeding rate and effective population size in the Finnish Ayrshire population in the era of genomic selection

Genomic selection has been applied in dairy cattle breeding over the last decade. Using genomic information may speed up genetic gain as breeding values can be predicted reasonably accurately directly after birth. However, genetic diversity may decrease if the inbreeding rate per generation increases and the effective population size decreases. Despite many positive qualities of the Finnish Ayrshire, for example, high average protein yield and fertility, over time the breed has lost its place as the most common dairy breed in Finland. Thus, maintaining the genetic variability of the breed is becoming more important. The aim of our research was to estimate the impact of genomic selection on inbreeding rate and effective population size using both pedigree and genomic data. The genomic data included 46,914 imputed single nucleotide polymorphism (SNP) variants from 75,038 individuals, and the pedigree data included 2,770,025 individuals. All animals in the data were born between 2000 and 2020. Genomic inbreeding coefficients were estimated as the proportion of SNPs in runs of homozygosity (ROH) out of the total number of SNPs. The inbreeding rate was estimated by regressing the mean genomic inbreeding coefficients on birth years. Effective population size was then estimated based on the inbreeding rate. Additionally, effective population size was estimated from the mean increase in individual inbreeding using pedigree data. Introduction of genomic selection was assumed to have taken place gradually; years 2012-2014 were treated as a transition period from the traditional phenotype-based breeding value estimation to genomic-based estimation. The median length of the identified homozygous segments was 5.5 Mbp, and a slight increase in the proportion of segments over 10 Mbp was observed after 2010. The inbreeding rate decreased from 2000 to 2011 and subsequently increased slightly. The pedigree- and genomic-based estimates of inbreeding rate were similar to each other. The estimates of effective population size based on the regression method were very sensitive to the number of years considered; thus, the estimates were not very reliable. The effective population size estimated from the mean increase in individual inbreeding reached its highest value of 160 in 2011 and decreased to 150 after that. In addition, the generation interval in the sire path has decreased from 5.5 years to 3.5 years after genomic selection was implemented. Based on our results, after the implementation of genomic selection, the proportion of long ROH stretches has increased, the generation interval in the sire path has decreased, the inbreeding rate has increased and the effective population size has decreased. However, the effective population size is still at a good level, allowing for an efficient selection scheme in the Finnish Ayrshire breed.

Comparative proteomics revealed duodenal metabolic function associated with feed efficiency in slow-growing chicken

The Korat chicken (KR), developed in Thailand, is a slow-growing breed developed as an alternative breed for Thai chicken producers. The growing interest in slow-growing chicken meat, due to its unique taste, distinct texture, health benefits, and higher broiler welfare have led to higher market demand for KR. However, its low feed efficiency (FE) has a significant negative impact on farm profitability. Understanding the molecular mechanism regulating FE allows for designing a suitable selection program and contributing to breeding more efficient chicken for poultry production. Thus, the objective of our study was to investigate the proteome differences and possible pathways associated with FE in male KR using a label-free quantitative proteomic approach. Seventy-five KR males were individually evaluated for FE, and duodenum samples from 6 animals (3 high-FE and 3 low-FE chickens) were collected at 10 wk of age for differential abundant proteins (DAPs), protein networks, functional enrichment, and pathway analyses. In this study, we found 40 DAPs significantly associated with FE pathways, including glycolysis/gluconeogenesis, peroxisome, oxidative phosphorylation, tight junction, and cysteine and methionine metabolism. Thus, variations in observed DAPs or genes related to DAPs could be interesting biomarker candidates for selection for higher feed utilization efficiency in chicken.

Single-step genomic evaluation of Russian dairy cattle using internal and external information

Genomic data are widely used in predicting the breeding values of dairy cattle. The accuracy of genomic prediction depends on the size of the reference population and how related the candidate animals are to it. For populations with limited numbers of progeny‐tested bulls, the reference populations must include cows and data from external populations. The aim of this study was to implement state‐of‐the‐art single‐step genomic evaluations for milk and fat yield in Holstein and Russian Black & White cattle in the Leningrad region (LR, Russia), using only a limited number of genotyped animals. We complemented internal information with external pseudo‐phenotypic and genotypic data of bulls from the neighbouring Danish, Finnish and Swedish Holstein (DFS) population. Three data scenarios were used to perform single‐step GBLUP predictions in the LR dairy cattle population. The first scenario was based on the original LR reference population, which constituted 1,080 genotyped cows and 427 genotyped bulls. In the second scenario, the genotypes of 414 bulls related to the LR from the DFS population were added to the reference population. In the third scenario, LR data were further augmented with pseudo‐phenotypic data from the DFS population. The inclusion of foreign information increased the validation reliability of the milk yield by up to 30%. Suboptimal data recording practices hindered the improvement of fat yield. We confirmed that the single‐step model is suitable for populations with a low number of genotyped animals, especially when external information is integrated into the evaluations. Genomic prediction in populations with a low number of progeny‐tested bulls can be based on data from genotyped cows and on the inclusion of genotypes and pseudo‐phenotypes from the external population. This approach increased the validation reliability of the implemented single‐step model in the milk yield, but shortcomings in the LR data recording scheme prevented improvements in fat yield.

Deletion of porcine BOLL is associated with defective acrosomes and subfertility in Yorkshire boars

A recessive sperm defect of Yorkshire boars was detected more than a decade ago. Affected boars produce ejaculates that contain spermatozoa with defective acrosomes, resulting in low fertility. The acrosome defect was mapped to porcine chromosome 15 but the causal mutation has not been identified. We re-analyzed microarray-derived genotypes of affected boars and confirmed that the acrosome defect maps to a 12.24 Mb segment of porcine chromosome 15. To detect the mutation causing defective acrosomes, we sequenced the genomes of two affected and three unaffected boars to an average coverage of 11-fold. Read depth analysis revealed a 55 kb deletion that is associated with the acrosome defect. The deletion encompasses the BOLL gene encoding the boule homolog, an RNA binding protein which is an evolutionarily conserved member of the DAZ (Deleted in AZoospermia) gene family. Lack of BOLL expression causes spermatogenic arrest and sperm maturation failure in many species. Boars that carry the deletion in the homozygous state produce sperm but their acrosomes are defective, suggesting that lack of porcine BOLL compromises acrosome formation. Our findings warrant further research to investigate the role of BOLL during spermatogenesis and sperm maturation in pigs.

Identification of runs of homozygosity affecting female fertility and milk production traits in Finnish Ayrshire cattle

Inbreeding gives rise to continuous lengths of homozygous genotypes called runs of homozygosity (ROH) that occur when identical haplotypes are inherited from both parents. ROHs are enriched for deleterious recessive alleles and can therefore be linked to inbreeding depression, defined as decreased phenotypic performance of the animals. However, not all ROHs within a region are expected to have harmful effects on the trait of interest. We aimed to identify ROHs that unfavourably affect female fertility and milk production traits in the Finnish Ayrshire population. The estimated effect of ROHs with the highest statistical significance varied between parities from 9 to 17 days longer intervals from calving to first insemination, from 13 to 38 days longer intervals from first to last insemination and from 0.3 to 1.0 more insemination per conception. Similarly, for milk production traits ROHs were associated with a reduction of 208 kg for milk yield, 7 kg for protein yield and 16 kg for fat yield. We also found regions where ROHs displayed unfavourable effects across multiple traits. Our findings can be exploited for more efficient control of inbreeding depression, for example by minimizing the occurrence of unfavourable haplotypes as homozygous state in breeding programmes.

Genome-wide SNP data unveils the globalization of domesticated pigs

BACKGROUND

Pigs were domesticated independently in Eastern and Western Eurasia early during the agricultural revolution, and have since been transported and traded across the globe. Here, we present a worldwide survey on 60K genome-wide single nucleotide polymorphism (SNP) data for 2093 pigs, including 1839 domestic pigs representing 122 local and commercial breeds, 215 wild boars, and 39 out-group suids, from Asia, Europe, America, Oceania and Africa. The aim of this study was to infer global patterns in pig domestication and diversity related to demography, migration, and selection.

RESULTS

A deep phylogeographic division reflects the dichotomy between early domestication centers. In the core Eastern and Western domestication regions, Chinese pigs show differentiation between breeds due to geographic isolation, whereas this is less pronounced in European pigs. The inferred European origin of pigs in the Americas, Africa, and Australia reflects European expansion during the sixteenth to nineteenth centuries. Human-mediated introgression, which is due, in particular, to importing Chinese pigs into the UK during the eighteenth and nineteenth centuries, played an important role in the formation of modern pig breeds. Inbreeding levels vary markedly between populations, from almost no runs of homozygosity (ROH) in a number of Asian wild boar populations, to up to 20% of the genome covered by ROH in a number of Southern European breeds. Commercial populations show moderate ROH statistics. For domesticated pigs and wild boars in Asia and Europe, we identified highly differentiated loci that include candidate genes related to muscle and body development, central nervous system, reproduction, and energy balance, which are putatively under artificial selection.

CONCLUSIONS

Key events related to domestication, dispersal, and mixing of pigs from different regions are reflected in the 60K SNP data, including the globalization that has recently become full circle since Chinese pig breeders in the past decades started selecting Western breeds to improve local Chinese pigs. Furthermore, signatures of ongoing and past selection, acting at different times and on different genetic backgrounds, enhance our insight in the mechanism of domestication and selection. The global diversity statistics presented here highlight concerns for maintaining agrodiversity, but also provide a necessary framework for directing genetic conservation.

Estimation of inbreeding depression on female fertility in the Finnish Ayrshire population

Single nucleotide polymorphism (SNP) data enable the estimation of inbreeding at the genome level. In this study, we estimated inbreeding levels for 19,075 Finnish Ayrshire cows genotyped with a low-density SNP panel (8K). The genotypes were imputed to 50K density, and after quality control, 39,144 SNPs remained for the analysis. Inbreeding coefficients were estimated for each animal based on the percentage of homozygous SNPs (F_PH ), runs of homozygosity (F_ROH ) and pedigree (F_PED ). Phenotypic records were available for 13,712 animals including non-return rate (NRR), number of inseminations (AIS) and interval from first to last insemination (IFL) for heifers and up to three parities for cows, as well as interval from calving to first insemination (ICF) for cows. Average F_PED was 0.02, F_ROH 0.06 and F_PH 0.63. A correlation of 0.71 was found between F_PED and F_ROH , 0.66 between F_PED and F_PH and 0.94 between F_ROH and F_PH . Pedigree-based inbreeding coefficients did not show inbreeding depression in any of the traits. However, when F_ROH or F_PH was used as a covariate, significant inbreeding depression was observed; a 10% increase in F_ROH was associated with 5 days longer IFL0 and IFL1, 2 weeks longer IFL3 and 3 days longer ICF2 compared to non-inbred cows.

Whole-genome association analysis of pork meat pH revealed three significant regions and several potential genes in Finnish Yorkshire pigs

Background

One of the most commonly used quality measurements of pork is pH measured 24 h after slaughter. The most probable mode of inheritance for this trait is oligogenic with several known major genes, such as PRKAG3. In this study, we used whole-genome SNP genotypes of over 700 AI boars; after a quality check, 42,385 SNPs remained for association analysis. All the boars were purebred Finnish Yorkshire. To account for relatedness of the animals, a pedigree-based relationship matrix was used in a mixed linear model to test the effect of SNPs on pH measured from loin. A bioinformatics analysis was performed to identify the most promising genes in the significant regions related to meat quality.

Results

Genome-wide association study (GWAS) revealed three significant chromosomal regions: one on chromosome 3 (39.9 Mb–40.1 Mb) and two on chromosome 15 (58.5 Mb–60.5 Mb and 132 Mb–135 Mb including PRKAG3). A conditional analysis with a significant SNP in the PRKAG3 region, MARC0083357, as a covariate in the model retained the significant SNPs on chromosome 3. Even though linkage disequilibrium was relatively high over a long distance between MARC0083357 and other significant SNPs on chromosome 15, some SNPs retained their significance in the conditional analysis, even in the vicinity of PRKAG3. The significant regions harbored several genes, including two genes involved in cyclic AMP (cAMP) signaling: ADCY9 and CREBBP. Based on functional and transcription factor-gene networks, the most promising candidate genes for meat pH are ADCY9, CREBBP, TRAP1, NRG1, PRKAG3, VIL1, TNS1, and IGFBP5, and the key transcription factors related to these genes are HNF4A, PPARG, and Nkx2-5.

Conclusions

Based on SNP association, pathway, and transcription factor analysis, we were able to identify several genes with potential to control muscle cell homeostasis and meat quality. The associated SNPs can be used in selection for better pork. We also showed that post-GWAS analysis reveals important information about the genes’ potential role on meat quality. The gained information can be used in later functional studies.

Heritability of hypothyroidism in the Finnish Hovawart population

BACKGROUND

The Hovawart is a working and companion dog breed of German origin. A few hundred Hovawart dogs are registered annually in Finland. The most common disease with a proposed genetic background in Hovawarts is hypothyroidism. The disease is usually caused by lymphocytic thyroiditis, an autoimmune disorder which destroys the thyroid gland. Hypothyroidism can be treated medically with hormone replacement. Its overall incidence could also be reduced through selection, provided that the trait shows an adequate genetic basis. The aim of this study was to estimate the heritability of hypothyroidism in the Finnish Hovawart population.

RESULTS

The pedigree data for the study were provided by the Finnish Kennel Club and the hypothyroidism data by the Finnish Hovawart Club. The data included 4953 dogs born between 1990 and 2010, of which 107 had hypothyroidism and 4846 were unaffected. Prior to the estimation of heritability, we studied the effects of gender, birth year, birth month, and inbreeding on susceptibility to hypothyroidism. Heritability was estimated with the probit model both via restricted maximum likelihood (REML) and Gibbs sampling, using litter and sire of the dog as random effects. None of the studied systematic effects or level of inbreeding had a significant effect on susceptibility to hypothyroidism. The estimated heritability of hypothyroidism varied from 0.47 (SE = 0.18) using REML to 0.62 (SD = 0.21) using Gibbs sampling.

CONCLUSIONS

Based on our analysis, the heritability of hypothyroidism is moderate to high, suggesting that its prevalence could be decreased through selection. Thus, breeders should notify the breed association of any affected dogs, and their use for breeding should be avoided.

A combination of two variants in PRKAG3 is needed for a positive effect on meat quality in pigs

BACKGROUND

Color and pH of meat measured 24 h post mortem are common selection objectives in pig breeding programs. Several amino acid substitutions in PRKAG3 have been associated with various meat quality traits. In our previous study ASGA0070625, a SNP next to PRKAG3, had the most significant association with meat quality traits in the Finnish Yorkshire. However, the known amino acid substitutions, including I199V, did not show any association. The aims of this study were to characterize further variation in PRKAG3 and its promoter region, and to test the association between these variants and the pH and color of pork meat.

RESULTS

The data comprised of 220 Finnish Landrace and 230 Finnish Yorkshire artificial insemination boars with progeny information. We sequenced the coding and promoter region of PRKAG3 in these and in three additional wild boars. Genotypes from our previous genome-wide scans were also included in the data. Association between SNPs or haplotypes and meat quality traits (deregressed estimates of breeding values from Finnish national breeding value estimation for pH, color lightness and redness measured from loin or ham) was tested using a linear regression model. Sequencing revealed several novel amino acid substitutions in PRKAG3, including K24E, I41V, K131R, and P134L. Linkage disequilibrium was strong among the novel variants, SNPs in the promoter region and ASGA0070625, especially for the Yorkshire. The strongest associations were observed between ASGA0070625 and the SNPs in the promoter region and pH measured from loin in the Yorkshire and between I199V and pH measured from ham in the Landrace. In contrast, ASGA0070625 was not significantly associated with meat quality traits in the Landrace and I199V not in the Yorkshire. Haplotype analysis showed a significant association between a haplotype consisting of 199I and 24E alleles (or g.-157C or g.-58A alleles in the promoter region) and pH measured from loin and ham in both breeds (P-values varied from 1.72 × 10⁻⁴ to 1.80 × 10⁻⁸).

CONCLUSIONS

We conclude that haplotype g.-157C - g.-58A - 24E - 199I in PRKAG3 has a positive effect on meat quality in pigs. Our results are readily applicable for marker-assisted selection in pigs.

Ectopic KIT copy number variation underlies impaired migration of primordial germ cells associated with gonadal hypoplasia in cattle (Bos taurus)

Impaired migration of primordial germ cells during embryonic development causes hereditary gonadal hypoplasia in both sexes of Northern Finncattle and Swedish Mountain cattle. The affected gonads exhibit a lack of or, in rare cases, a reduced number of germ cells. Most affected animals present left-sided gonadal hypoplasia. However, right-sided and bilateral cases are also found. This type of gonadal hypoplasia prevails in animals with white coat colour. Previous studies indicated that gonadal hypoplasia is inherited in an autosomal recessive fashion with incomplete penetrance. In order to identify genetic regions underlying gonadal hypoplasia, a genome-wide association study (GWAS) and a copy number variation (CNV) analysis were performed with 94 animals, including 21 affected animals, using bovine 777,962 SNP arrays. The GWAS and CNV results revealed two significantly associated regions on bovine chromosomes (BTA) 29 and 6, respectively (P=2.19 x 10^-13 and P=5.65 x 10^-6). Subsequent cytogenetic and PCR analyses demonstrated that homozygosity of a ~500 kb chromosomal segment translocated from BTA6 to BTA29 (Cs₂₉ allele) is the underlying genetic mechanism responsible for gonadal hypoplasia. The duplicated segment includes the KIT gene that is known to regulate the migration of germ cells and precursors of melanocytes. This duplication is also one of the two translocations associated with colour sidedness in various cattle breeds.

Evidence for three highly significant QTL for meat quality traits in the Finnish Yorkshire pig breed

Meat quality is important both to consumers and to the meat processing industry. Commonly used measures of porcine meat quality are the pH and color of the meat. The purpose of this study was to identify SNP associated with these meat quality traits in Finnish Yorkshire using the Illumina PorcineSNP60 BeadChip. The association of each SNP with the quality traits was tested with a weighted linear model. The relatedness of samples was accounted for by a random polygenic genetic effect with the accompanying full relationship matrix. The original EBV from single-trait evaluations were deregressed before analysis. The statistical significance of SNP was established using the Bonferroni correction to adjust for multiple testing. Three genomic regions were significant for the meat quality traits. The PRKAG3 region on chromosome 15 was significant for pH measured from loin and ham and for a* (redness) measured from loin. The smallest P-value in the region was obtained for pH measured from loin (ASGA0070634, P-value = 3.8 × 10(-13)). The allele substitution effect (-0.047) of the unfavorable allele A corresponds to 1 SD of the polygenic effect. The second significant region, on chromosome 2 at around 31 megabases (Mb), was associated with pH and L* (lightness) measured from loin. The most significant SNP (ASGA0009814, P-value = 3.89 × 10(-10)) had an allele substitution effect of 0.86, corresponding to 0.7 SD of the polygenic effect of L*. The third region, located on chromosome 6 at around 83 Mb, was significant for a* measured from ham. The P-value of the best SNP (ALGA0035896) was 8.71 × 10(-7) and the allele substitution effect -0.38, corresponding to 0.5 SD of the polygenic effect of a*. The significant association of PRKAG3 with pH was not due to the known AA substitutions. The candidate gene on chromosome 2 associated with color L* is RCN1, which has a high affinity Ca(2+)-binding motif, the EF hand. The significant region on chromosome 6 for color a* contains several genes, so more data are needed to identify the causative gene. Our results indicate that instead of the known AA substitutions of PRKAG3, some yet-unknown AA substitutions are causative for the pH variation in Finnish Yorkshire. Also, a new major QLT for L* was found on chromosome 2. The significant SNP identified in this study can be used in marker-assisted selection.

Genome-Wide Linkage Disequilibrium from 100,000 SNPs in the East Finland Founder Population

Information about linkage disequilibrium (LD) is important in understanding the genome structure and has its applications in association studies. Here we present the first genome-wide LD study based on a founder population (East Finland). The LD data consist of 118 unrelated individuals and around 480,000 SNP pairs genotyped with the Affymetrix 100K genotyping assay. Using the minor allele frequency (MAF) limit of .05, the squared correlation coefficient between two loci (r2) was .48, .37, .28, and .20 for distances of 5, 10, 20, and 40 kb respectively. MAF had a significant effect on the mean r2 so that the extent of useful LD (r2 > .3) varied from 17 kb to 80 kb depending on the limit set for the MAF. For D' the effect of MAF was smaller but reflected the possible age of the mutation: SNPs with high MAF had lower D' than those with low MAF. The X chromosome showed higher D' values than autosomes and the extent of useful LD (r2 > .3) was twice as long on the X chromosome than on the autosomes. Based on the results, LD varies across the genome and is correlated to local recombination rate between and within chromosomes. However, the recombination rate does not explain all the variation found in LD. We also report a number of long chromosomal regions where exceptionally high or low LD were detected.

An exonic insertion within Tex14 gene causes spermatogenic arrest in pigs

BACKGROUND

Male infertility is an increasing problem in all domestic species including man. Localization and identification of genes involved in defects causing male infertility provide valuable information of specific events in sperm development. Sperm development is a complex process, where diploid spermatogonia develop into haploid, highly specialized spermatozoa. Correct expression and function of various genes and their protein products are required for production of fertile sperm. We have identified an infertility defect in Finnish Yorkshire boars caused by spermatogenic arrest. The aim of this study was to locate the disease associated region using genome wide screen with the PorcineSNP60 Beadchip and identify the causal mutation by candidate gene approach.

RESULTS

In the Finnish Yorkshire pig population the spermatogenic arrest (SA) defect appears to be of genetic origin and causes severe degeneration of germ cells and total absence of spermatozoa. Genome wide scan with the PorcineSNP60 Beadchip localized the SA defect to porcine chromosome 12 in a 2 Mbp region. Sequencing of a candidate gene Tex14 revealed a 51 bp insertion within exon 27, which caused differential splicing of the exon and created a premature translation stop codon. The expression of Tex14 was markedly down regulated in the testis of a SA affected boar compared to control boars and no protein product was identified by Western blotting. The SA insertion sequence was also found within intron 27 in all analyzed animals, thus the insertion appears to be a possible duplication event.

CONCLUSION

In this study we report the identification of a causal mutation for infertility caused by spermatogenic arrest at an early meiotic phase. Our results highlight the role of TEX14 specifically in spermatogenesis and the importance of specific genomic remodeling events as causes for inherited defects.

Whole-genome SNP association analysis of reproduction traits in the Finnish Landrace pig breed

Background

Good genetic progress for pig reproduction traits has been achieved using a quantitative genetics-based multi-trait BLUP evaluation system. At present, whole-genome single nucleotide polymorphisms (SNP) panels provide a new tool for pig selection. The purpose of this study was to identify SNP associated with reproduction traits in the Finnish Landrace pig breed using the Illumina PorcineSNP60 BeadChip.

Methods

Association of each SNP with different traits was tested with a weighted linear model, using SNP genotype as a covariate and animal as a random variable. Deregressed estimated breeding values of the progeny tested boars were used as the dependent variable and weights were based on their reliabilities. Statistical significance of the associations was based on Bonferroni-corrected P-values.

Results

Deregressed estimated breeding values were available for 328 genotyped boars. Of the 62 163 SNP in the chip, 57 868 SNP had a call rate > 0.9 and 7 632 SNP were monomorphic. Statistically significant results (P-value < 2.0E-06) were obtained for total number of piglets born in first and later parities and piglet mortality between birth and weaning in later parity, and suggestive associations (P-value < 4.0E-06) for piglet mortality between birth and weaning in first parity, number of stillborn piglets in later parity, first farrowing interval and second farrowing interval. Two of the statistically significant regions for total number of piglets born in first and later parities are located on chromosome 9 around 95 and 79 Mb. The estimated SNP effect in these regions was approximately one piglet between the two homozygote classes. By combining the two most significant SNP in these regions, favourable double homozygote animals are expected to have 1.3 piglets (P-value = 1.69E-08) more than unfavourable double homozygote animals. A region on chromosome 9 (66 Mb) was statistically significant for piglet mortality between birth and weaning in later parity (0.44 piglets between homozygotes, P-value = 6.94E-08).

Conclusions

Three separate regions on chromosome 9 gave significant results for litter size and pig mortality. The frequencies of favourable alleles of the significant SNP are moderate in the Finnish Landrace population and these SNP are thus valuable candidates for possible marker-assisted selection.

L1 insertion within SPEF2 gene is associated with increased litter size in the Finnish Yorkshire population

Immotile, short-tail sperm defect (ISTS) expanded in the Finnish Yorkshire population in the end of 1990s. The causal mutation for this defect is a recent L1 insertion within the SPEF2 gene in chromosome 16. Even though all homozygous boars are eliminated from the population because of infertility, the amount of affected boars increased rapidly until marker-assisted selection against the defect was established. To elucidate the associated effects of the ISTS defect on production traits, we have investigated the association of the L1 insertion and PRLR haplotype with reproduction traits in the Finnish Yorkshire population. Two data sets including 357 sows and 491 AI-boars were genotyped for the presence of the L1 insertion and analysed for association with reproduction traits. A Proc Mixed procedure (SAS Inc) and a software package for analysing multivariate mixed models (DMU) were used to study the effect of polymorphisms on reproduction traits. The L1-insertion within SPEF2 gene was associated with litter size in the first parity. The SPEF2 gene is located adjacent to a candidate gene for litter size in the pig, PRLR. Haplotypes within PRLR exon 10 were analysed in data set of 93 AI-boars for the association with reproduction traits. However, no associations were detected within the analysed data set indicating that PRLR sequence variants are not the causal cause for the identified effect on litter size.

Mapping of QTL affecting incidence of blood and meat inclusions in egg layers

BACKGROUND

Occurrence of blood and meat inclusions is an internal egg quality defect. Mass candling reveals most of the spots, but because brown eggshell hampers selection in brown chicken lines it has not been possible to eliminate the defect by selection. Estimated frequency of blood and meat inclusions in brown layers is about 18% whereas it is 0.5% in white egg layers. Several factors are known to increase the incidence of this fault: genetic background, low level of vitamin A and/or D, stress or infections, for instance. To study the genetic background of the defect, a mapping population of 1599 F2 hens from a cross of White Rock and Rhode Island Red lines was set up.

RESULTS

Our histopathological analyses show that blood spots consist of mainly erythrocytes and that meat spots are accumulations of necrotic material. Linkage analysis of 27 chromosomes with 162 microsatellite markers revealed one significant quantitative trait locus (QTL) affecting blood spot and meat spot frequency. We sequenced a fragment of a candidate gene within the region, ZO-2, coding for a tight junction protein. Nine polymorphisms were detected and two of them were included in fine-mapping and association analysis. Fine-mapping defined the QTL result. To further verify the QTL, association analyses were carried out in two independent commercial breeding lines with the marker MCW241 and surrounding SNPs. Association was found mainly in a 0.8 Mb-wide chromosomal area on GGAZ.

CONCLUSIONS

There was good agreement between the location of the QTL region on chromosome Z and the association results in the commercial breeds analyzed. Variations found in tight junction protein ZO-2 and microRNA gga-mir-1556 may predispose egg layers to blood and meat spot defects. This paper describes the first results of detailed QTL analyses of the blood and meat spots trait(s) in chickens.

Knobbed acrosome defect is associated with a region containing the genes STK17b and HECW2 on porcine chromosome 15

BACKGROUND

Male infertility is an increasing problem in all domestic species including man. Localization and identification of genes involved in defects causing male infertility provide valuable information of specific events in sperm development. Correct condensation of the sperm head and development of the acrosome are required for fertile sperm. In the Finnish Yorkshire pig population a knobbed acrosome defect (KAD) has been reported which appears to be of genetic origin. In previous studies we have shown that a large number of affected spermatozoa have a cystic swelling anterior to the apical part of the acrosome.

RESULTS

Characterization of the knobbed acrosome affected sperm revealed that both the acrosomal granules and chromatin are affected. This type of KAD appears to be a previously unknown and serious form of the defect. A genome wide scan with PorcineSNP60 Genotyping BeadChip defined the KAD associated region within 0.7 Mbp on porcine chromosome 15. Two genes, STK17b and HECW2, located within this region were sequenced. The expression of these genes appeared comparable in KA-affected and control boars. The known function of HECW2 in acrosome development highlighted this gene as a good candidate responsible for the KAD. One nonsynonymous SNP was identified within the HECW2 gene. However, as this mutation was found in homozygous state in individuals with normal sperm, this is not likely to be the causal mutation.

CONCLUSIONS

In this study we identified two candidate genes for a severe defect affecting both the sperm acrosome and chromatin that causes infertility. One of these genes, HECW2, plays an important role in ubiquitination, a prerequisite for chromatin remodelling and acrosome formation, highlighting the involvement of this gene in the knobbed acrosome defect and male infertility.

Extent of linkage disequilibrium and effective population size in Finnish Landrace and Finnish Yorkshire pig breeds

The extent of linkage disequilibrium (LD) and effective population size in Finnish Landrace and Finnish Yorkshire pig populations were studied using a whole genome SNP panel (Illumina PorcineSNP60 BeadChip) and pedigree data. Genotypic data included 86 Finnish Landrace and 32 Finnish Yorkshire boars. Pedigree data included 608,138 Finnish Landrace 554,237 and Finnish Yorkshire pigs, and on average 15 ancestral generations were known for the reference animals, born in 2005 to 2009. The breeding animals of the 2 populations have been kept separate in the breeding programs. Based on the pedigree data, the current effective population size for Finnish Landrace is 91 and for Finnish Yorkshire 61. Linkage disequilibrium measures (D' and r(2)) were estimated for over 1.5 million pairs of SNP. Average r(2) for SNP 30 kb apart was 0.47 and 0.49 and for SNP 5 Mb apart 0.09 and 0.12 for Finnish Landrace and Finnish Yorkshire, respectively. Average LD (r(2)) between adjacent SNP in the Illumina PorcineSNP60 BeadChip was 0.43 (57% of the adjacent SNP pairs had r(2) > 0.2) for Finnish Landrace and 0.46 (60% of the adjacent SNP pairs had r(2) > 0.2) for Finnish Yorkshire, and average r(2) > 0.2 extended to 1.0 and 1.5 Mb for Finnish Landrace and Finnish Yorkshire, respectively. Effective population size estimates based on the decay of r(2) with distance were similar to those based on the pedigree data: 80 and 55 for Finnish Landrace and Finnish Yorkshire, respectively. Thus, the results indicate that the effective population size of Finnish Yorkshire is smaller than of Finnish Landrace and has a clear effect on the extent of LD. The current effective population size of both breeds is above the recommended minimum of 50 but may get smaller than that in the near future, if no action is taken to balance the inbreeding rate and selection response. Because a moderate level of LD extends over a long distance, selection based on whole genome SNP markers (genomic selection) is expected to be efficient for both breeds.

Type 2 diabetes whole-genome association study in four populations: the DiaGen consortium

Type 2 diabetes (T2D) is a common, polygenic chronic disease with high heritability. The purpose of this whole-genome association study was to discover novel T2D-associated genes. We genotyped 500 familial cases and 497 controls with >300,000 HapMap-derived tagging single-nucleotide-polymorphism (SNP) markers. When a stringent statistical correction for multiple testing was used, the only significant SNP was at TCF7L2, which has already been discovered and confirmed as a T2D-susceptibility gene. For a replication study, we selected 10 SNPs in six chromosomal regions with the strongest association (singly or as part of a haplotype) for retesting in an independent case-control set including 2,573 T2D cases and 2,776 controls. The most significant replicated result was found at the AHI1-LOC441171 gene region.

Genome-wide linkage disequilibrium from 100,000 SNPs in the East Finland founder population

Information about linkage disequilibrium (LD) is important in understanding the genome structure and has its applications in association studies. Here we present the first genome-wide LD study based on a founder population (East Finland). The LD data consist of 118 unrelated individuals and around 480,000 SNP pairs genotyped with the Affymetrix 100K genotyping assay. Using the minor allele frequency (MAF) limit of .05, the squared correlation coefficient between two loci (r(2)) was .48, .37, .28, and .20 for distances of 5, 10, 20, and 40 kb respectively. MAF had a significant effect on the mean r(2) so that the extent of useful LD (r(2) > .3) varied from 17 kb to 80 kb depending on the limit set for the MAF. For D' the effect of MAF was smaller but reflected the possible age of the mutation: SNPs with high MAF had lower D' than those with low MAF. The X chromosome showed higher D' values than autosomes and the extent of useful LD (r(2) > .3) was twice as long on the X chromosome than on the autosomes. Based on the results, LD varies across the genome and is correlated to local recombination rate between and within chromosomes. However, the recombination rate does not explain all the variation found in LD. We also report a number of long chromosomal regions where exceptionally high or low LD were detected.

Haplotypic association of DDAH1 with susceptibility to pre-eclampsia

Association between pre-eclampsia (PEE1) and the dimethylarginine dimethylaminohydrolase (DDAH) 1 and 2 genes, which play a role in the regulation of nitric oxide synthesis and release, was studied. In a case-control study design single nucleotide polymorphisms (SNPs) were determined at eight sites in the DDAH1 gene and at one site (Pro231Pro) in the DDAH2 gene from 132 women with pre-eclampsia and 112 healthy controls. Three SNPs in the DDAH1 gene were associated with pre-eclampsia, showing complete linkage disequilibrium with each other, but none of the associations in the allele or genotype data reached statistical significance in either of the genes after the correction for multiple testing. Haplotype frequencies were estimated using a population based on a maximum likelihood method (EM algorithm). Four common DDAH1 haplotypes were present and a significant association of haplotypes H2 and H3 with pre-eclampsia (P=0.03) was found. The risk of pre-eclampsia was greatest in individuals (odds ratio: 3.93; 95% confidence interval: 1.54-9.99) who had two copies of the high-risk haplotypes (H2 or H3). The observed haplotypic association provides the first evidence of the importance of DDAH1 polymorphisms in pre-eclampsia susceptibility.

New paraoxonase 1 polymorphism I102V and the risk of prostate cancer in Finnish men

BACKGROUND

Human serum paraoxonase eliminates carcinogenic lipid-soluble radicals. Because expression of the main human paraoxonase gene PON1 varies widely in humans, certain PON1 polymorphisms might be associated with increased risks of cancer. We sought new functional mutations in PON1 and determined whether known or new PON1 mutations were associated with the risk for prostate cancer in a prospective, random, population-based sample of Finnish men and in a case-control study.

METHODS

Serum paraoxonase activity was measured in 835 healthy men in the Kuopio Ischaemic Heart Disease Risk Factor Study. PON1 mutations were identified by hierarchical phenotype-targeted sequencing in DNAs from the 100 men with the lowest paraoxonase activity in this cohort, and 1595 men in the cohort were genotyped for PON1 mutations by restriction fragment length polymorphism. Multivariable analysis was used to investigate the association of known and new PON1 mutations with incident prostate cancer in 1569 cancer-free men in the cohort followed for 9-14 years. In a case-control study of Finnish men, the association of prostate cancer with the PON1 mutation identified in the cohort study was investigated in 69 case patients with familial prostate cancer and 69 unmatched healthy control subjects.

RESULTS

We identified a new single-nucleotide PON1 polymorphism associated with decreased serum paraoxonase activity that caused an isoleucine-->valine change at codon 102 in exon 4 (I102V). Of the 1569 men cancer-free at baseline, 56 (3.6%) were carriers of the I102V mutation. After adjusting for age and cholesterol-lowering medications, the relative risk for developing prostate cancer during follow-up was 6.3 (95% confidence interval [CI] = 2.1 to 19.2) among 102V allele carriers compared with noncarriers. Other PON1 alleles were not statistically significantly associated with prostate cancer. In the case-control study, patients with familial prostate cancer were more likely to be carriers of the PON1 I102V mutation than control subjects (odds ratio = 4.3, 95% CI = 0.9 to 21.5).

CONCLUSION

The PON1 102V allele appears to be associated with an increased risk for prostate cancer.

Chromosome 19q13 and multiple sclerosis susceptibility in Finland: a linkage and two-stage association study

Several studies have previously provided some albeit weak evidence for linkage or association between chromosome 19q13 and multiple sclerosis (MS) susceptibility. We performed a two-stage association analysis with 19 markers spanning 7 Mb/5.5 cM of 19q13. In stage 1 analysis (135 MS families) allelic and haplotypic associations were found with markers within or close to the ApoE-ApoC subregion. These observations were taken as a hypothesis, which was tested in stage 2 in 125 families. However, none of the initial associations were replicated suggesting that they were most likely due to chance. Linkage analysis was performed in 27 Finnish multiplex families using 10 microsatellites spanning 23 Mb/24 cM of 19q13. DNA was available from 72 MS patients and 150 unaffected relatives. Parametric and non-parametric linkage analyses did not provide evidence for linkage when all families were tested. After stratifying the families according to HLA-DR15 there was weak evidence for linkage to the 19q13.1 subregion in DR15 negative families (LOD(max)=1.8). Taken together these results do not support a major role of chromosome 19q13.2-q13.3 in MS susceptibility among Finnish MS patients, whereas conclusions on the 19q13.1 subregion are less clear and this region requires further study.

Bayesian association mapping for quantitative traits in a mixture of two populations

We introduce a novel Bayesian approach to estimate and account for population structure simultaneously with association mapping of multiple quantitative trait loci. The method is designed for an analysis of unrelated individuals from a mixture of two populations (no admixture), where the individual population memberships are unknown. In our approach, the population structure is estimated and accounted for by using data on additional "grouping" markers which are assumed to be in Hardy-Weinberg equilibrium within the populations but have different allele frequencies between the populations. We use Bayesian hierarchical modeling and Markov chain Monte Carlo estimation, where we allow both population stratification and genetic heterogeneity. In our model the number of quantitative trait loci and their positions are treated as random variables, and we obtain their posterior distributions. Here we select the candidate and the grouping markers based on results from a preliminary SOLAR analysis.

Multiple QTL mapping in related plant populations via a pedigree-analysis approach

QTL mapping experiments in plant breeding may involve multiple populations or pedigrees that are related through their ancestors. These known relationships have often been ignored for the sake of statistical analysis, despite their potential increase in power of mapping. We describe here a Bayesian method for QTL mapping in complex plant populations and reported the results from its application to a (previously analysed) potato data set. This Bayesian method was originally developed for human genetics data, and we have proved that it is useful for complex plant populations as well, based on a sensitivity analysis that was performed here. The method accommodates robustness to complex structures in pedigree data, full flexibility in the estimation of the number of QTL across multiple chromosomes, thereby accounting for uncertainties in the transmission of QTL and marker alleles due to incomplete marker information, and the simultaneous inclusion of non-genetic factors affecting the quantitative trait.

Mapping of an immotile short tail sperm defect in the Finnish Yorkshire on porcine Chromosome 16

An immotile short tail sperm defect has recently been identified as a hereditary disorder present within the Finnish Yorkshire pig population. The syndrome is inherited as an autosomal recessive disease exclusively expressed in male individuals as shorter sperm tail length and immotile spermatozoa. Based on the assumption of a recent common origin of the disease-causing mutation, a genome-wide search was performed with 228 evenly spaced microsatellites by homozygosity mapping of affected and unaffected DNA pools. One locus, SW2411 on Chr 16, demonstrated a significantly skewed allele distribution between the two pools. Linkage analysis of five markers in this region mapped the disease-causing gene within a 6-cM confidence interval region with a highest LOD score of 7.7 at marker SW419. It appears that three-marker haplotypes can be used for marker-assisted selection within analyzed pedigrees. Furthermore, future fine mapping may reveal a more precise population-wide associated haplotype and facilitate identification of a new gene affecting sperm tail development.

Bayesian oligogenic analysis of quantitative and qualitative traits in general pedigrees

A Bayesian method for multipoint oligogenic analysis of quantitative and qualitative traits is presented. This method can be applied to general pedigrees, which do not necessarily have to be "peelable" and can have large numbers of markers. The number of quantitative/qualitative trait loci (QTL), their map positions in the genome, and phenotypic effects (mode of inheritances) are all estimated simultaneously within the same framework. The summaries of the estimated parameters are based on the marginal posterior distributions that are obtained through Markov chain Monte Carlo (MCMC) methods. The method uses founder alleles together with segregation indicators in order to determine the genotypes of the trait loci of all individuals in the pedigree. To improve mixing properties of the sampler, we propose (1) joint sampling of map position and segregation indicators, (2) omitting data augmentation for untyped or uninformative markers (homozygous parent), and (3) updating several markers jointly within a single block. The performance of the method was tested with two replicate GAW10 data sets (considering two levels of available marker information). The results were concordant and similar to those presented earlier with other methods. These analyses clearly illustrate the utility and wide applicability of the method.

Respiratory distress syndrome: evaluation of genetic susceptibility and protection by transmission disequilibrium test

Based on epidemiological data and genetic association studies, neonatal respiratory distress syndrome (RDS) is a complex disease with a multigenic background. The genes coding for surfactant proteins (SP) A and B have been assigned as the most likely genes in the etiology of RDS. The major factor predisposing to RDS is prematurity, and thus the phenotype of a very premature newborn infant that does not develop the disease can be regarded as hypernormal. Altogether 107 father-mother-offspring trios were divided into two sets according to the proband's phenotype, to evaluate familial segregation of candidate gene polymorphisms by the transmission disequilibrium test. A set of 76 trios were analyzed for transmission disequilibrium from parents to affected offspring. Another set of 31 trios were studied for allele transmission from parents to hypernormal offspring born very prematurely before the gestational age of 32 weeks. SP-A1-A2 haplotype 6A(2)-1A(0) showed significant excess transmission to affected infants and SP-A1 allele 6A(2) decreased transmission to the hypernormals. The present family study provides strong support for a direct or indirect role of the SP-A alleles as genetic predisposers to RDS in premature infants. The inclusion of parent-hypernormal offspring trios in transmission disequilibrium test is a useful approach to test for genetic protection against a disease.

Candidate gene regions and genetic heterogeneity in gluten sensitivity

BACKGROUND

Gluten sensitivity is a common multifactorial disorder, manifested in the small intestine or on the skin as typical coeliac disease or dermatitis herpetiformis, respectively. The only established genetic risk factor is HLA DQ2.

AIMS

We tested genetic linkage of previously reported chromosomal loci 5q and 11q in Finnish families with gluten sensitivity. We also tested if genetic linkage to candidate loci on 5q, 11q, 2q33, and HLA DQ differed with respect to clinical manifestations or sex.

SUBJECTS

We studied 102 Finnish families with affected sibpairs. For heterogeneity analysis, families were divided into subgroups according to sex and the presence of dermatitis herpetiformis, the skin manifestation of gluten sensitivity.

METHODS

Non-parametric linkage between microsatellite markers and disease was tested. Linkage heterogeneity between subgroups was tested using the M test. The transmission/disequilibrium test and association analysis were performed.

RESULTS

Evidence of linkage to 11q (MLS 1.37), but not to 5q, was found in the entire dataset of 102 families. Heterogeneity between subgroups was suggested: families with only the intestinal disease showed linkage mainly to 2q33 whereas families with dermatitis herpetiformis showed linkage to 11q and 5q, but not to 2q33. Linkage in all three non-HLA loci was strongest in families with predominantly male patients. HLA DQ2 conferred much stronger susceptibility to females than males.

CONCLUSIONS

Independent evidence for the suggested genetic linkage between 11q and gluten sensitivity was obtained. The possible linkage heterogeneity suggests genetic differences between intestinal and skin manifestations, and the gender dependent effect of HLA DQ2.

A Bayesian Markov chain Monte Carlo approach to map disease genes in simulated GAW11 data

A Bayesian method for multipoint mapping of disease genes based on Markov chain Monte Carlo algorithms was applied to the simulated GAW11 data (Study 2). The method is based on repeated Gibbs and more general Metropolis-Hastings steps. For simplicity we assumed a single disease locus model with two alleles. A normal distribution for the underlying latent variable of the qualitative phenotype was assumed. Based on a single replicate of the data no clear evidence of any of the genes underlying the simulated disease was found. However, when three replicates were combined the method was able to locate the locus C correctly on chromosome 3.

The effect of using different culling regimens on genetic response with two-trait, two-stage selection in a nucleus broiler stock

Stochastic simulation was used to study the effect on genetic response and inbreeding of various two-stage two-trait culling strategies. Four different parameter sets were considered for the two traits, BW and egg number. Selection of replacement animals was based on animal model best linear unbiased prediction (BLUP) to obtain estimated breeding values (EBV) at the second stage. Culling at Stage 1 was based on either animal model BLUP or phenotypes, and information from culled animals was either available or not available for calculation of second stage EBV. Besides founder individuals, six discrete generations were considered. Culling based on BLUP of two traits at Stage 1 produced higher response than culling on phenotypic evaluations. It was found that culling based on phenotypic evaluation and not carrying information to the second stage reduce rates of response by 9 to 17% and produced inbreeding higher than or close to that of BLUP selection. This study clearly shows that a double penalty of less response and higher inbreeding is generally paid for not using all information. Optimum selection schemes will depend on relative costs and benefits of collecting and processing the extra information required for full BLUP selection schemes.

Mapping quantitative trait loci for milk production and health of dairy cattle in a large outbred pedigree

Quantitative trait loci (QTL) affecting milk production and health of dairy cattle were mapped in a very large Holstein granddaughter design. The analysis included 1794 sons of 14 sires and 206 genetic markers distributed across all 29 autosomes and flanking an estimated 2497 autosomal cM using Kosambi's mapping function. All families were analyzed jointly with least-squares (LS) and variance components (VC) methods. A total of 6 QTL exceeding approximate experiment-wise significance thresholds, 24 QTL exceeding suggestive thresholds, and 34 QTL exceeding chromosome-wise thresholds were identified. Significance thresholds were determined via data permutation (for LS analysis) and chi-square distribution (for VC analysis). The average bootstrap confidence interval for the experiment-wise significant QTL was 48 cM. Some chromosomes harbored QTL affecting several traits, and these were always in coupling phase, defined by consistency with genetic correlations among traits. Chromosome 17 likely harbors 2 QTL affecting milk yield, and some other chromosomes showed some evidence for 2 linked QTL affecting the same trait. In each of these cases, the 2 QTL were in repulsion phase in those families appearing to be heterozygous for both QTL, a finding which supports the build-up of linkage disequilibrium due to selection.

Advances in statistical methods to map quantitative trait loci in outbred populations

Statistical methods to map quantitative trait loci (QTL) in outbred populations are reviewed, extensions and applications to human and plant genetic data are indicated, and areas for further research are identified. Simple and computationally inexpensive methods include (multiple) linear regression of phenotype on marker genotypes and regression of squared phenotypic differences among relative pairs on estimated proportions of identity-by-descent at a locus. These methods are less suited for genetic parameter estimation in outbred populations but allow the determination of test statistic distributions via simulation or data permutation; however, further inferences including confidence intervals of QTL location require the use of Monte Carlo or bootstrap sampling techniques. A method which is intermediate in computational requirements is residual maximum likelihood (REML) with a covariance matrix of random QTL effects conditional on information from multiple linked markers. Testing for the number of QTLs on a chromosome is difficult in a classical framework. The computationally most demanding methods are maximum likelihood and Bayesian analysis, which take account of the distribution of multilocus marker-QTL genotypes on a pedigree and permit investigators to fit different models of variation at the QTL. The Bayesian analysis includes the number of QTLs on a chromosome as an unknown.

Mapping-linked quantitative trait loci using Bayesian analysis and Markov chain Monte Carlo algorithms

A Bayesian method for mapping linked quantitative trait loci (QTL) using multiple linked genetic markers is presented. Parameter estimation and hypothesis testing was implemented via Markov chain Monte Carlo (MCMC) algorithms. Parameters included were allele frequencies and substitution effects for two biallelic QTL, map positions of the QTL, and markers, allele frequencies of the markers, and polygenic and residual variances. Missing data were polygenic effects and multi-locus marker-QTL genotypes. Three different MCMC schemes for testing the presence of a single or two linked QTL on the chromosome were compared. The first approach includes a model indicator variable representing two unlinked QTL, affecting the trait, one linked and one unlinked QTL, or both QTL linked with the markers. The second approach incorporates an indicator variable for each QTL into the model for phenotype, allowing or not allowing for a substitution effect of a QTL, on phenotype, and the third approach is based on model determination by reversible jump MCMC. Methods were evaluated empirically by analyzing simulated granddaughter designs. All methods identified correctly a second, linked QTL and did not reject the one-QTL model when there was only a single QTL, and no additional or an unlinked QTL.

Sensitivity of segregation analysis to data structure and transformation: a case study of trypanotolerance in mice

Sensitivity of segregation analysis for data structure and data transformation was studied using data from two trials in which mice were challenged at three months of age with a cloned isolate of Trypanosoma congolense and survival time was recorded. Data included records from three inbred strains (C57BL/6 (tolerant), A/J, and BALB/c (both susceptible)) and their crosses. Data were standardized and normalized using a modified power transformation. Segregation analysis was applied to both untransformed and transformed data to determine the genetic inheritance of trypanotolerance in these mice. Data from the two trials were analysed separately and combined. Four genetic models were compared; a one locus model, a polygenic model, a mixed model with common variance, and a mixed model with different variances for each major genotype. Even though the separate data sets and the combined data set all supported the hypothesis of a major gene (or a tightly linked cluster of genes) with different variances within each genotype, parameter estimates were highly sensitive to data transformation and several sets of parameter estimates gave similar likelihood values because of high dependency between parameters. Based on the results segregation analysis can be very sensitive to data structure in a crossbreeding design and to data transformation. Interpretation of the results can be misleading if the entire parameter space is not studied carefully.

The use of multiple markers in a Bayesian method for mapping quantitative trait loci

Information on multiple linked genetic markers was used in a Bayesian method for the statistical mapping of quantitative trait loci (QTL). Bayesian parameter estimation and hypothesis testing were implemented via Markov chain Monte Carlo algorithms. Variables sampled were the augmented data (marker-QTL genotypes, polygenic effects), an indicator variable for linkage or nonlinkage, and the parameters. The parameter vector included allele frequencies at the markers and the QTL, map distances of the markers and the QTL, QTL substitution effect, and polygenic and residual variances. The criterion for QTL detection was the marginal posterior probability of a QTL being located on the chromosome carrying the markers. The method was evaluated empirically by analyzing simulated granddaughter designs consisting of 2000 sons, 20 related sires, and their ancestors.