Comparison of selective genotyping strategies for prediction of breeding values in a population undergoing selection

Genomewide marker information can improve the reliability of breeding value predictions for young selection candidates in genomic selection. However, the cost of genotyping limits its use to elite animals, and how such selective genotyping affects predictive ability of genomic selection models is an open question. We performed a simulation study to evaluate the quality of breeding value predictions for selection candidates based on different selective genotyping strategies in a population undergoing selection. The genome consisted of 10 chromosomes of 100 cM each. After 5,000 generations of random mating with a population size of 100 (50 males and 50 females), generation G(0) (reference population) was produced via a full factorial mating between the 50 males and 50 females from generation 5,000. Different levels of selection intensities (animals with the largest yield deviation value) in G(0) or random sampling (no selection) were used to produce offspring of G(0) generation (G(1)). Five genotyping strategies were used to choose 500 animals in G(0) to be genotyped: 1) Random: randomly selected animals, 2) Top: animals with largest yield deviation values, 3) Bottom: animals with lowest yield deviations values, 4) Extreme: animals with the 250 largest and the 250 lowest yield deviations values, and 5) Less Related: less genetically related animals. The number of individuals in G(0) and G(1) was fixed at 2,500 each, and different levels of heritability were considered (0.10, 0.25, and 0.50). Additionally, all 5 selective genotyping strategies (Random, Top, Bottom, Extreme, and Less Related) were applied to an indicator trait in generation G(0), and the results were evaluated for the target trait in generation G(1), with the genetic correlation between the 2 traits set to 0.50. The 5 genotyping strategies applied to individuals in G(0) (reference population) were compared in terms of their ability to predict the genetic values of the animals in G(1) (selection candidates). Lower correlations between genomic-based estimates of breeding values (GEBV) and true breeding values (TBV) were obtained when using the Bottom strategy. For Random, Extreme, and Less Related strategies, the correlation between GEBV and TBV became slightly larger as selection intensity decreased and was largest when no selection occurred. These 3 strategies were better than the Top approach. In addition, the Extreme, Random, and Less Related strategies had smaller predictive mean squared errors (PMSE) followed by the Top and Bottom methods. Overall, the Extreme genotyping strategy led to the best predictive ability of breeding values, indicating that animals with extreme yield deviations values in a reference population are the most informative when training genomic selection models.

Mapping QTLs on chicken chromosome 1 for performance and carcass traits in a broiler x layer cross

With the objective of mapping quantitative trait loci (QTLs) for performance and carcass traits, an F-2 chicken population was developed by crossing broiler and layer lines. A total of 2063 F-2 chicks in 21 full-sib families were reared as broilers and slaughtered at 42 days of age. Seventeen performance and carcass traits were measured. Parental (F-0) and F-1 individuals were genotyped with 80 microsatellites from chicken chromosome 1 to select informative markers. Thirty-three informative markers were used for selective genotyping of F-2 individuals with extreme phenotypes for body weight at 42 days of age (BW42). Based on the regions identified by selective genotyping, seven full-sib families (649 F-2 chicks) were genotyped with 26 markers. Quantitative trait loci affecting body weight, feed intake, carcass weight, drums and thighs weight and abdominal fat weight were mapped to regions already identified in other populations. Quantitative trait loci for weights of gizzard, liver, lungs, heart and feet, as well as length of intestine, not previously described in the literature were mapped on chromosome 1. This F-2 population can be used to identify novel QTLs and constitutes a new resource for studies of genes related to growth and carcass traits in poultry.

Analysis of the pig genome for the identification of genomic regions affecting production traits

The aim of this work was to identify markers associated with production traits in the pig genome using different approaches. We focused the attention on Italian Large White pig breed using Genome Wide Association Studies (GWAS) and applying a selective genotyping approach to increase the power of the analyses. Furthermore, we searched the pig genome using Next Generation Sequencing (NSG) Ion Torrent Technology to combine selective genotyping approach and deep sequencing for SNP discovery. Other two studies were carried on with a different approach. Allele frequency changes for SNPs affecting candidate genes and at Genome Wide level were analysed to identify selection signatures driven by selection program during the last 20 years. This approach confirmed that a great number of markers may affect production traits and that they are captured by the classical selection programs. GWAS revealed 123 significant or suggestively significant SNP associated with Back Fat Thickenss and 229 associated with Average Daily Gain. 16 Copy Number Variant Regions resulted more frequent in lean or fat pigs and showed that different copies of those region could have a limited impact on fat. These often appear to be involved in food intake and behavior, beside affecting genes involved in metabolic pathways and their expression. By combining NGS sequencing with selective genotyping approach, new variants where discovered and at least 54 are worth to be analysed in association studies. The study of groups of pigs undergone to stringent selection showed that allele frequency of some loci can drastically change if they are close to traits that are interesting for selection schemes. These approaches could be, in future, integrated in genomic selection plans.

USE OF HIDDEN MARKOV MODELS FOR QTL MAPPING

An important aspect of the QTL mapping problem is the treatment of missing genotype data. If complete genotype data were available, QTL mapping would reduce to the problem of model selection in linear regression. However, in the consideration of loci in the intervals between the available genetic markers, genotype data is inherently missing. Even at the typed genetic markers, genotype data is seldom complete, as a result of failures in the genotyping assays or for the sake of economy (for example, in the case of selective genotyping, where only individuals with extreme phenotypes are genotyped). We discuss the use of algorithms developed for hidden Markov models (HMMs) to deal with the missing genotype data problem.

Sequential Quantitative Trait Locus Mapping in Experimental Crosses

The etiology of complex diseases is heterogeneous. The presence of risk alleles in one or more genetic loci affects the function of a variety of intermediate biological pathways, resulting in the overt expression of disease. Hence, there is an increasing focus on identifying the genetic basis of disease by sytematically studying phenotypic traits pertaining to the underlying biological functions. In this paper we focus on identifying genetic loci linked to quantitative phenotypic traits in experimental crosses. Such genetic mapping methods often use a one stage design by genotyping all the markers of interest on the available subjects. A genome scan based on single locus or multi-locus models is used to identify the putative loci. Since the number of quantitative trait loci (QTLs) is very likely to be small relative to the number of markers genotyped, a one-stage selective genotyping approach is commonly used to reduce the genotyping burden, whereby markers are genotyped solely on individuals with extreme trait values. This approach is powerful in the presence of a single quantitative trait locus (QTL) but may result in substantial loss of information in the presence of multiple QTLs. Here we investigate the efficiency of sequential two stage designs to identify QTLs in experimental populations. Our investigations for backcross and F2 crosses suggest that genotyping all the markers on 60% of the subjects in Stage 1 and genotyping the chromosomes significant at 20% level using additional subjects in Stage 2 and testing using all the subjects provides an efficient approach to identify the QTLs and utilizes only 70% of the genotyping burden relative to a one stage design, regardless of the heritability and genotyping density. Complex traits are a consequence of multiple QTLs conferring main effects as well as epistatic interactions. We propose a two-stage analytic approach where a single-locus genome scan is conducted in Stage 1 to identify promising chromosomes, and interactions are examined using the loci on these chromosomes in Stage 2. We examine settings under which the two-stage analytic approach provides sufficient power to detect the putative QTLs.

Genetic association studies under the population stratification, family pedigree and application to genome-wide association studies

This dissertation has three separate parts: the first part deals with the general pedigree association testing incorporating continuous covariates; the second part deals with the association tests under population stratification using the conditional likelihood tests; the third part deals with the genome-wide association studies based on the real rheumatoid arthritis (RA) disease data sets from Genetic Analysis Workshop 16 (GAW16) problem 1. Many statistical tests are developed to test the linkage and association using either case-control status or phenotype covariates for family data structure, separately. Those univariate analyses might not use all the information coming from the family members in practical studies. On the other hand, the human complex disease do not have a clear inheritance pattern, there might exist the gene interactions or act independently. In part I, the new proposed approach MPDT is focused on how to use both the case control information as well as the phenotype covariates. This approach can be applied to detect multiple marker effects. Based on the two existing popular statistics in family studies for case-control and quantitative traits respectively, the new approach could be used in the simple family structure data set as well as general pedigree structure. The combined statistics are calculated using the two statistics; A permutation procedure is applied for assessing the p-value with adjustment from the Bonferroni for the multiple markers. We use simulation studies to evaluate the type I error rates and the powers of the proposed approach. Our results show that the combined test using both case-control information and phenotype covariates not only has the correct type I error rates but also is more powerful than the other existing methods. For multiple marker interactions, our proposed method is also very powerful. Selective genotyping is an economical strategy in detecting and mapping quantitative trait loci in the genetic dissection of complex disease. When the samples arise from different ethnic groups or an admixture population, all the existing selective genotyping methods may result in spurious association due to different ancestry distributions. The problem can be more serious when the sample size is large, a general requirement to obtain sufficient power to detect modest genetic effects for most complex traits. In part II, I describe a useful strategy in selective genotyping while population stratification is present. Our procedure used a principal component based approach to eliminate any effect of population stratification. The paper evaluates the performance of our procedure using both simulated data from an early study data sets and also the HapMap data sets in a variety of population admixture models generated from empirical data. There are one binary trait and two continuous traits in the rheumatoid arthritis dataset of Problem 1 in the Genetic Analysis Workshop 16 (GAW16): RA status, AntiCCP and IgM. To allow multiple traits, we suggest a set of SNP-level F statistics by the concept of multiple-correlation to measure the genetic association between multiple trait values and SNP-specific genotypic scores and obtain their null distributions. Hereby, we perform 6 genome-wide association analyses using the novel one- and two-stage approaches which are based on single, double and triple traits. Incorporating all these 6 analyses, we successfully validate the SNPs which have been identified to be responsible for rheumatoid arthritis in the literature and detect more disease susceptibility SNPs for follow-up studies in the future. Except for chromosome 13 and 18, each of the others is found to harbour susceptible genetic regions for rheumatoid arthritis or related diseases, i.e., lupus erythematosus. This topic is discussed in part III.

Potential applications of reproductive and molecular genetic technologies in the selective breeding of aquaculture species

The use of reproductive and genetic technologies can increase the efficiency of selective breeding programs for aquaculture species. Four technologies are considered, namely: marker-assisted selection, DNA fingerprinting, in-vitro fertilization, and cryopreservation. Marker-assisted selection can result in greater genetic gain, particularly for traits difficult or expensive to measure, than conventional selection methods, but its application is currently limited by lack of high density linkage maps and by the high cost of genotyping. DNA fingerprinting is most useful for genetic tagging and parentage verification. Both in-vitro fertilization and cryopreservation techniques can increase the accuracy of selection while controlling accumulation of inbreeding in long-term selection programs. Currently, the cost associated with the utilization of reproductive and genetic techniques is possibly the most important factor limiting their use in genetic improvement programs for aquatic species.

Fluorescent amplified fragment length polymorphism genotyping of human and animal Staphylococcus aureus isolates from dairy farms with manual milking

Fluorescence amplified fragment length polymorphism (fAFLP) was used to assess the genetic relatedness of 40 Staphylococcus aureus strains isolated from human and animal skin samples in seven dairy farms with manual milking. S. aureus was isolated from 11 out of 30 (36%) human skin samples and from 29 out of 100 (29%) teat skin samples from apparently healthy cows. Genomic DNA from each isolate was double-digested with EcoRI and MseI and complementary oligonucleotide adaptors were ligated to the restriction fragments. Pre-selective and selective, amplification reactions were performed, the amplified fragments were separated by electrophoresis in an ABI377 sequencer and analysed using GeneScan 3.1 and Genotyper 2.5. Three single isolates (a-c), a predominant cluster with 35 isolates (d) and another cluster with two isolates (e) were identified. Both clusters d and e included human and animal isolates genetically related, because the profiles had 90-100% homology. Since no cluster was comprised uniquely of human or animal isolates and given the close genetic relatedness among human and animal samples in the farms, the present findings support the. hypothesis that dairy workers can spread S. aureus through manual milking. (C) 2005 Elsevier B.V. All rights reserved.

Genomic divergence of zebu and taurine cattle identified through high-density SNP genotyping

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Linkage & Selective Mapping and QTL Analysis in a Harosoy x Clark population of 300 RILs

Four of the 12 major Glycine max ancestors of all modern elite U.S.A. soybean cultivars were the grandparents of Harosoy and Clark, so a Harosoy x Clark population would include some of that genetic diversity. A mating of eight Harosoy and eight Clark plants generated eight F1 plants. The eight F1:2 families were advanced via a plant-to-row selfing method to produce 300 F6-derived RILs that were genotyped with 266 SSR, 481 SNP, and 4 classical markers. SNPs were genotyped with the Illumina 1536-SNP assay. Three linkage maps, SSR, SNP, and SSR-SNP, were constructed with a genotyping error of < 1 %. Each map was compared with the published soybean consensus map. The best subset of 94 RILs for a high-resolution framework (joint) map was selected based on the expected bin length statistic computed with MapPop. The QTLs of seven traits measured in a 2-year replicated performance trial of the 300 RILs were identified using composite interval mapping (CIM) and multiple-interval mapping (MIM). QTL x Year effects in multiple trait analysis were compared with results of multiple-interval mapping. QTL x QTL effects were identified in MIM.

Campylobacter jejuni and campylobacter coli genotyping by high-resolution melting analysis of a flaA fragment

The highly variable flagellin-encoding flaA gene has long been used for genotyping Campylobacter jejuni and Campylobacter coli. High-resolution melting (HRM) analysis is emerging as an efficient and robust method for discriminating DNA sequence variants. The objective of this study was to apply HRM analysis to flaA-based genotyping. The initial aim was to identify a suitable flaA fragment. It was found that the PCR primers commonly used to amplify the flaA short variable repeat (SVR) yielded a mixed PCR product unsuitable for HRM analysis. However, a PCR primer set composed of the upstream primer used to amplify the fragment used for flaA restriction fragment length polymorphism (RFLP) analysis and the downstream primer used for flaA SVR amplification generated a very pure PCR product, and this primer set was used for the remainder of the study. Eighty-seven C. jejuni and 15 C. coli isolates were analyzed by flaA HRM and also partial flaA sequencing. There were 47 flaA sequence variants, and all were resolved by HRM analysis. The isolates used had previously also been genotyped using single-nucleotide polymorphisms (SNPs), binary markers, CRISPR HRM, and flaA RFLP. flaAHRManalysis provided resolving power multiplicative to the SNPs, binary markers, and CRISPR HRM and largely concordant with the flaA RFLP. It was concluded that HRM analysis is a promising approach to genotyping based on highly variable genes.