A local score approach improves GWAS resolution and detects minor QTL: application to Medicago truncatula quantitative disease resistance to multiple Aphanomyces euteiches isolates.

Quantitative trait loci (QTL) with small effects, which are pervasive in quantitative phenotypic variation, are difficult to detect in genome-wide association studies (GWAS). To improve their detection, we propose to use a local score approach that accounts for the surrounding signal due to linkage disequilibrium, by accumulating association signals from contiguous single markers. Simulations revealed that, in a GWAS context with high marker density, the local score approach outperforms single SNP p-value-based tests for detecting minor QTL (heritability of 5-10%) and is competitive with regard to alternative methods, which also aggregate p-values. Using more than five million SNPs, this approach was applied to identify loci involved in Quantitative Disease Resistance (QDR) to different isolates of the plant root rot pathogen Aphanomyces euteiches, from a GWAS performed on a collection of 174 accessions of the model legume Medicago truncatula. We refined the position of a previously reported major locus, underlying MYB/NB-ARC/tyrosine kinase candidate genes conferring resistance to two closely related A. euteiches isolates belonging to pea pathotype I. We also discovered a diversity of minor resistance QTL, not detected using p-value-based tests, some of which being putatively shared in response to pea (pathotype I and III) and/or alfalfa (race 1 and 2) isolates. Candidate genes underlying these QTL suggest pathogen effector recognition and plant proteasome as key functions associated with M. truncatula resistance to A. euteiches. GWAS on any organism can benefit from the local score approach to uncover many weak-effect QTL.

Accounting for linkage disequilibrium in genome scans for selection without individual genotypes: The local score approach.

Detecting genomic footprints of selection is an important step in the understanding of evolution. Accounting for linkage disequilibrium in genome scans increases detection power, but haplotype-based methods require individual genotypes and are not applicable on pool-sequenced samples. We propose to take advantage of the local score approach to account for linkage disequilibrium in genome scans for selection, cumulating (possibly small) signals from single markers over a genomic segment, to clearly pinpoint a selection signal. Using computer simulations, we demonstrate that this approach detects selection with higher power than several state-of-the-art single-marker, windowing or haplotype-based approaches. We illustrate this on two benchmark data sets including individual genotypes, for which we obtain similar results with the local score and one haplotype-based approach. Finally, we apply the local score approach to Pool-Seq data obtained from a divergent selection experiment on behaviour in quail and obtain precise and biologically coherent selection signals: while competing methods fail to highlight any clear selection signature, our method detects several regions involving genes known to act on social responsiveness or autistic traits. Although we focus here on the detection of positive selection from multiple population data, the local score approach is general and can be applied to other genome scans for selection or other genomewide analyses such as GWAS.

A medium density genetic map and QTL for behavioral and production traits in Japanese quail.

BACKGROUND: Behavioral traits such as sociability, emotional reactivity and aggressiveness are major factors in animal adaptation to breeding conditions. In order to investigate the genetic control of these traits as well as their relationships with production traits, a study was undertaken on a large second generation cross (F2) between two lines of Japanese Quail divergently selected on their social reinstatement behavior. All the birds were measured for several social behaviors (social reinstatement, response to social isolation, sexual motivation, aggression), behaviors measuring the emotional reactivity of the birds (reaction to an unknown object, tonic immobility reaction), and production traits (body weight and egg production). RESULTS: We report the results of the first genome-wide QTL detection based on a medium density SNP panel obtained from whole genome sequencing of a pool of individuals from each divergent line. A genetic map was constructed using 2145 markers among which 1479 could be positioned on 28 different linkage groups. The sex-averaged linkage map spanned a total of 3057 cM with an average marker spacing of 2.1 cM. With the exception of a few regions, the marker order was the same in Japanese Quail and the chicken, which confirmed a well conserved synteny between the two species. The linkage analyses performed using QTLMAP software revealed a total of 45 QTLs related either to behavioral (23) or production (22) traits. The most numerous QTLs (15) concerned social motivation traits. Interestingly, our results pinpointed putative pleiotropic regions which controlled emotional reactivity and body-weight of birds (on CJA5 and CJA8) or their social motivation and the onset of egg laying (on CJA19). CONCLUSION: This study identified several QTL regions for social and emotional behaviors in the Quail. Further research will be needed to refine the QTL and confirm or refute the role of candidate genes, which were suggested by bioinformatics analysis. It can be hoped that the identification of genes and polymorphisms related to behavioral traits in the quail will have further applications for other poultry species (especially the chicken) and will contribute to solving animal welfare issues in poultry production.

An image analysis method to quantify CFTR subcellular localization.

Aberrant protein subcellular localization caused by mutation is a prominent feature of many human diseases. In Cystic Fibrosis (CF), a recessive lethal disorder that results from dysfunction of the Cystic Fibrosis Transmembrane Conductance Regulator (CFTR), the most common mutation is a deletion of phenylalanine-508 (pF508del). Such mutation produces a misfolded protein that fails to reach the cell surface. To date, over 1900 mutations have been identified in CFTR gene, but only a minority has been analyzed at the protein level. To establish if a particular CFTR variant alters its subcellular distribution, it is necessary to quantitatively determine protein localization in the appropriate cellular context. To date, most quantitative studies on CFTR localization have been based on immunoprecipitation and western blot. In this work, we developed and validated a confocal microscopy-image analysis method to quantitatively examine CFTR at the apical membrane of epithelial cells. Polarized MDCK cells transiently transfected with EGFP-CFTR constructs and stained for an apical marker were used. EGFP-CFTR fluorescence intensity in a region defined by the apical marker was normalized to EGFP-CFTR whole cell fluorescence intensity, rendering "apical CFTR ratio". We obtained an apical CFTR ratio of 0.67 ± 0.05 for wtCFTR and 0.11 ± 0.02 for pF508del. In addition, this image analysis method was able to discriminate intermediate phenotypes: partial rescue of the pF508del by incubation at 27 °C rendered an apical CFTR ratio value of 0.23 ± 0.01. We concluded the method has a good sensitivity and accurately detects milder phenotypes. Improving axial resolution through deconvolution further increased the sensitivity of the system as rendered an apical CFTR ratio of 0.76 ± 0.03 for wild type and 0.05 ± 0.02 for pF508del. The presented procedure is faster and simpler when compared with other available methods and it is therefore suitable as a screening method to identify mutations that completely or mildly affect CFTR processing. Moreover, it could be extended to other studies on the biology underlying protein subcellular localization in health and disease.

Testing the existence of an unadmixed ancestor from a specific population t generations ago.

The ancestry of each locus of the genome can be estimated (local ancestry) based on sequencing or genotyping information together with reference panels of ancestral source populations. The length of those ancestry-specific genomic segments are commonly used to understand migration waves and admixture events. In short time scales, it is often of interest to determine the existence of the most recent unadmixed ancestor from a specific population t generations ago. We built a hypothesis test to determine if an individual has an ancestor belonging to a target ancestral population t generations ago based on these lengths of the ancestry-specific segments at an individual level. We applied this test on a data set that includes 20 Uruguayan admixed individuals to estimate for each one how many generations ago the most recent indigenous ancestor lived. As this method tests each individual separately, it is particularly suited to small sample sizes, such as our study or ancient genome samples.

Indigenous Ancestry and Admixture in the Uruguayan Population.

The Amerindian group known as the Charrúas inhabited Uruguay at the timing of European colonial contact. Even though they were extinguished as an ethnic group as a result of a genocide, Charrúan heritage is part of the Uruguayan identity both culturally and genetically. While mitochondrial DNA studies have shown evidence of Amerindian ancestry in living Uruguayans, here we undertake whole-genome sequencing of 10 Uruguayan individuals with self-declared Charruan heritage. We detect chromosomal segments of Amerindian ancestry supporting the presence of indigenous genetic ancestry in living descendants. Specific haplotypes were found to be enriched in "Charrúas" and rare in the rest of the Amerindian groups studied. Some of these we interpret as the result of positive selection, as we identified selection signatures and they were located mostly within genes related to the infectivity of specific viruses. Historical records describe contacts of the Charrúas with other Amerindians, such as Guaraní, and patterns of genomic similarity observed here concur with genomic similarity between these groups. Less expected, we found a high genomic similarity of the Charrúas to Diaguita from Argentinian and Chile, which could be explained by geographically proximity. Finally, by fitting admixture models of Amerindian and European ancestry for the Uruguayan population, we were able to estimate the timing of the first pulse of admixture between European and Uruguayan indigenous peoples in approximately 1658 and the second migration pulse in 1683. Both dates roughly concurring with the Franciscan missions in 1662 and the foundation of the city of Colonia in 1680 by the Spanish.

Kinesin 1 regulates cilia length through an interaction with the Bardet-Biedl syndrome related protein CCDC28B.

Bardet-Biedl syndrome (BBS) is a ciliopathy characterized by retinal degeneration, obesity, polydactyly, renal disease and mental retardation. CCDC28B is a BBS-associated protein that we have previously shown plays a role in cilia length regulation whereby its depletion results in shortened cilia both in cells and Danio rerio (zebrafish). At least part of that role is achieved by its interaction with the mTORC2 component SIN1, but the mechanistic details of this interaction and/or additional functions that CCDC28B might play in the context of cilia remain poorly understood. Here we uncover a novel interaction between CCDC28B and the kinesin 1 molecular motor that is relevant to cilia. CCDC28B interacts with kinesin light chain 1 (KLC1) and the heavy chain KIF5B. Notably, depletion of these kinesin 1 components results in abnormally elongated cilia. Furthermore, through genetic interaction studies we demonstrate that kinesin 1 regulates ciliogenesis through CCDC28B. We show that kinesin 1 regulates the subcellular distribution of CCDC28B, unexpectedly, inhibiting its nuclear accumulation, and a ccdc28b mutant missing a nuclear localization motif fails to rescue the phenotype in zebrafish morphant embryos. Therefore, we uncover a previously unknown role of kinesin 1 in cilia length regulation that relies on the BBS related protein CCDC28B.

Selection signatures in worldwide sheep populations.

The diversity of populations in domestic species offers great opportunities to study genome response to selection. The recently published Sheep HapMap dataset is a great example of characterization of the world wide genetic diversity in sheep. In this study, we re-analyzed the Sheep HapMap dataset to identify selection signatures in worldwide sheep populations. Compared to previous analyses, we made use of statistical methods that (i) take account of the hierarchical structure of sheep populations, (ii) make use of linkage disequilibrium information and (iii) focus specifically on either recent or older selection signatures. We show that this allows pinpointing several new selection signatures in the sheep genome and distinguishing those related to modern breeding objectives and to earlier post-domestication constraints. The newly identified regions, together with the ones previously identified, reveal the extensive genome response to selection on morphology, color and adaptation to new environments.

Detecting signatures of selection through haplotype differentiation among hierarchically structured populations.

The detection of molecular signatures of selection is one of the major concerns of modern population genetics. A widely used strategy in this context is to compare samples from several populations and to look for genomic regions with outstanding genetic differentiation between these populations. Genetic differentiation is generally based on allele frequency differences between populations, which are measured by FST or related statistics. Here we introduce a new statistic, denoted hapFLK, which focuses instead on the differences of haplotype frequencies between populations. In contrast to most existing statistics, hapFLK accounts for the hierarchical structure of the sampled populations. Using computer simulations, we show that each of these two features-the use of haplotype information and of the hierarchical structure of populations-significantly improves the detection power of selected loci and that combining them in the hapFLK statistic provides even greater power. We also show that hapFLK is robust with respect to bottlenecks and migration and improves over existing approaches in many situations. Finally, we apply hapFLK to a set of six sheep breeds from Northern Europe and identify seven regions under selection, which include already reported regions but also several new ones. We propose a method to help identifying the population(s) under selection in a detected region, which reveals that in many of these regions selection most likely occurred in more than one population. Furthermore, several of the detected regions correspond to incomplete sweeps, where the favorable haplotype is only at intermediate frequency in the population(s) under selection.