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.

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.

Design, Synthesis And Biological Evaluation Of Hybrid Bioisoster Derivatives Of N-acylhydrazone And Furoxan Groups With Potential And Selective Anti-trypanosoma Cruzi Activity.

Hybrid bioisoster derivatives from N-acylhydrazones and furoxan groups were designed with the objective of obtaining at least a dual mechanism of action: cruzain inhibition and nitric oxide (NO) releasing activity. Fifteen designed compounds were synthesized varying the substitution in N-acylhydrazone and in furoxan group as well. They had its anti-Trypanosoma cruzi activity in amastigotes forms, NO releasing potential and inhibitory cruzain activity evaluated. The two most active compounds (6, 14) both in the parasite amastigotes and in the enzyme contain the nitro group in para position of the aromatic ring. The permeability screening in Caco-2 cell and cytotoxicity assay in human cells were performed for those most active compounds and both showed to be less cytotoxic than the reference drug, benznidazole. Compound 6 was the most promising, since besides activity it showed good permeability and selectivity index, higher than the reference drug. Thereby the compound 6 was considered as a possible candidate for additional studies.

Improved Selective Cholesterol Adsorption By Molecularly Imprinted Poly(methacrylic Acid)/silica (pmaa-sio₂) Hybrid Material Synthesized With Different Molar Ratios.

The present paper describes the synthesis of molecularly imprinted polymer - poly(methacrylic acid)/silica and reports its performance feasibility with desired adsorption capacity and selectivity for cholesterol extraction. Two imprinted hybrid materials were synthesized at different methacrylic acid (MAA)/tetraethoxysilane (TEOS) molar ratios (6:1 and 1:5) and characterized by FT-IR, TGA, SEM and textural data. Cholesterol adsorption on hybrid materials took place preferably in apolar solvent medium, especially in chloroform. From the kinetic data, the equilibrium time was reached quickly, being 12 and 20 min for the polymers synthesized at MAA/TEOS molar ratio of 6:1 and 1:5, respectively. The pseudo-second-order model provided the best fit for cholesterol adsorption on polymers, confirming the chemical nature of the adsorption process, while the dual-site Langmuir-Freundlich equation presented the best fit to the experimental data, suggesting the existence of two kinds of adsorption sites on both polymers. The maximum adsorption capacities obtained for the polymers synthesized at MAA/TEOS molar ratios of 6:1 and 1:5 were found to be 214.8 and 166.4 mg g(-1), respectively. The results from isotherm data also indicated higher adsorption capacity for both imprinted polymers regarding to corresponding non-imprinted polymers. Nevertheless, taking into account the retention parameters and selectivity of cholesterol in the presence of structurally analogue compounds (5-α-cholestane and 7-dehydrocholesterol), it was observed that the polymer synthesized at the MAA/TEOS molar ratio of 6:1 was much more selective for cholesterol than the one prepared at the ratio of 1:5, thus suggesting that selective binding sites ascribed to the carboxyl group from MAA play a central role in the imprinting effect created on MIP.

Genetic Diversity Of Giardia Duodenalis: Multilocus Genotyping Reveals Zoonotic Potential Between Clinical And Environmental Sources In A Metropolitan Region Of Brazil.

Giardia duodenalis is a flagellate protozoan that parasitizes humans and several other mammals. Protozoan contamination has been regularly documented at important environmental sites, although most of these studies were performed at the species level. There is a lack of studies that correlate environmental contamination and clinical infections in the same region. The aim of this study is to evaluate the genetic diversity of a set of clinical and environmental samples and to use the obtained data to characterize the genetic profile of the distribution of G. duodenalis and the potential for zoonotic transmission in a metropolitan region of Brazil. The genetic assemblages and subtypes of G. duodenalis isolates obtained from hospitals, a veterinary clinic, a day-care center and important environmental sites were determined via multilocus sequence-based genotyping using three unlinked gene loci. Cysts of Giardia were detected at all of the environmental sites. Mixed assemblages were detected in 25% of the total samples, and an elevated number of haplotypes was identified. The main haplotypes were shared among the groups, and new subtypes were identified at all loci. Ten multilocus genotypes were identified: 7 for assemblage A and 3 for assemblage B. There is persistent G. duodenalis contamination at important environmental sites in the city. The identified mixed assemblages likely represent mixed infections, suggesting high endemicity of Giardia in these hosts. Most Giardia isolates obtained in this study displayed zoonotic potential. The high degree of genetic diversity in the isolates obtained from both clinical and environmental samples suggests that multiple sources of infection are likely responsible for the detected contamination events. The finding that many multilocus genotypes (MLGs) and haplotypes are shared by different groups suggests that these sources of infection may be related and indicates that there is a notable risk of human infection caused by Giardia in this region.

Short-Term effect of COX-2 selective inhibitor as an adjunct for the treatment of periodontal disease: a clinical double-blind study in humans

Adjunctive therapeutic strategies that modulate the inflammatory mediators can play a significant role in periodontal therapy. In this double-blind, placebo-controlled study, 60 subjects diagnosed as periodontitis patients were evaluated for 28 days after periodontal treatment combined with selective cyclooxygenase-2 (COX-2) inhibitor. The experimental group received scaling and root planning (SRP) combined with the Loxoprofen antiinflammatory drug (SRP+Loxoprofen). The control group received SRP combined with placebo (SRP+placebo). Plaque index (PI), probing pocket depth (PD) and bleeding on probing (BOP) were monitored with an electronic probe at baseline and after 14 and 28 days. Both groups displayed clinical improvement in PD, PI and BOP. They also showed statistically similar values (p>0.05) of PD reduction on day 14 (0.4 mm) and on day 28 (0.6 mm). At the baseline, few deeper sites (>7 mm) from SRP+Loxoprofen group were responsible and most PD reduction was observed after 14 days (p<0.05). The percentage of remaining deep pockets (>7 mm) after 14 days in the SRP+Loxoprofen group was significantly lower (p<0.05) than in the SRP+placebo group. Loxoprofen presents potential effect as an adjunct of periodontal disease treatment, but long-term clinical trials are necessary to confirm its efficacy.

Molecular markers and genetic diversity of Plasmodium vivax

Enhanced understanding of the transmission dynamics and population genetics for Plasmodium vivax is crucial in predicting the emergence and spread of novel parasite phenotypes with major public health implications, such as new relapsing patterns, drug resistance and increased virulence. Suitable molecular markers are required for these population genetic studies. Here, we focus on two groups of molecular markers that are commonly used to analyse natural populations of P. vivax. We use markers under selective pressure, for instance, antigen-coding polymorphic genes, and markers that are not under strong natural selection, such as most minisatellite and microsatellite loci. First, we review data obtained using genes encoding for P. vivax antigens: circumsporozoite protein, merozoite surface proteins 1 and 3α, apical membrane antigen 1 and Duffy binding antigen. We next address neutral or nearly neutral molecular markers, especially microsatellite loci, providing a complete list of markers that have already been used in P. vivax populations studies. We also analyse the microsatellite loci identified in the P. vivax genome project. Finally, we discuss some practical uses for P. vivax genotyping, for example, detecting multiple-clone infections and tracking the geographic origin of isolates.

Performance of transport and selective media for swine Bordetella bronchiseptica recovery and it comparison to polymerase chain reaction detection

Three comparative assays were performed seeking to improve the sensitivity of the diagnosis of Bordetella bronchiseptica infection analyzing swine nasal swabs. An initial assay compared the recovery of B. bronchiseptica from swabs simultaneously inoculated with B. bronchiseptica and some interfering bacteria, immersed into three transport formulations (Amies with charcoal, trypticase soy broth and phosphate buffer according to Soerensen supplemented with 5% of bovine fetal serum) and submitted to different temperatures (10ºC and 27ºC) and periods of incubation (24, 72 and 120 hours). A subsequent assay compared three selective media (MacConkey agar, modified selective medium G20G and a ceftiofur medium) for their recovery capabilities from clinical specimens. One last assay compared the polymerase chain reaction to the three selective media. In the first assay, the recovery of B. bronchiseptica from transport systems was better at 27ºC and the three formulations had good performances at this temperature, but the collection of qualitative and quantitative analysis indicated the advantage of Amies medium for nasal swabs transportation. The second assay indicated that MacConkey agar and modified G20G had similar results and were superior to the ceftiofur medium. In the final assay, polymerase chain reaction presented superior capability of B. bronchiseptica detection to culture procedures.