GWAlpha: genome-wide estimation of additive effects (alpha) based on trait quantile distribution from pool-sequencing experiments

MOTIVATION: Sequencing pools of individuals (Pool-Seq) is a cost-effective way to gain insight into the genetics of complex traits, but as yet no parametric method has been developed to both test for genetic effects and estimate their magnitude. Here, we propose GWAlpha, a flexible method to obtain parametric estimates of genetic effects genome-wide from Pool-Seq experiments. RESULTS: We showed that GWAlpha powerfully replicates the results of Genome-Wide Association Studies (GWAS) from model organisms. We perform simulation studies that illustrate the effect on power of sample size and number of pools and test the method on different experimental data. AVAILABILITY AND IMPLEMENTATION: GWAlpha is implemented in python, designed to run on Linux operating system and tested on Mac OS. It is freely available at https://github.com/aflevel/GWAlpha CONTACT: [email protected] information: Supplementary data are available at Bioinformatics online

Construction of nested genetic core collections to optimize the exploitation of natural diversity in Vitis vinifera L. subsp. sativa

Background: The first high quality draft of the grape genome sequence has just been published. This is a critical step in accessing all the genes of this species and increases the chances of exploiting the natural genetic diversity through association genetics. However, our basic knowledge of the extent of allelic variation within the species is still not sufficient. Towards this goal, we constructed nested genetic core collections (G-cores) to capture the simple sequence repeat (SSR) diversity of the grape cultivated compartment (Vitis vinifera L. subsp. sativa) from the world's largest germplasm collection (Domaine de Vassal, INRA Hérault, France), containing 2262 unique genotypes. Results: Sub-samples of 12, 24, 48 and 92 varieties of V. vinifera L. were selected based on their genotypes for 20 SSR markers using the M-strategy. They represent respectively 58%, 73%, 83% and 100% of total SSR diversity. The capture of allelic diversity was analyzed by sequencing three genes scattered throughout the genome on 233 individuals: 41 single nucleotide polymorphisms (SNPs) were identified using the G-92 core (one SNP for every 49 nucleotides) while only 25 were observed using a larger sample of 141 individuals selected on the basis of 50 morphological traits, thus demonstrating the reliability of the approach. Conclusion: The G-12 and G-24 core-collections displayed respectively 78% and 88% of the SNPs respectively, and are therefore of great interest for SNP discovery studies. Furthermore, the nested genetic core collections satisfactorily reflected the geographic and the genetic diversity of grape, which are also of great interest for the study of gene evolution in this species

Structural Variants and Selective Sweep Foci Contribute to Insecticide Resistance in the Drosophila Genetic Reference Panel

Patterns of nucleotide polymorphism within populations of Drosophila melanogaster suggest that insecticides have been the selective agents driving the strongest recent bouts of positive selection. However, there is a need to explicitly link selective sweeps to the particular insecticide phenotypes that could plausibly account for the drastic selective responses that are observed in these non-target insects. Here, we screen the Drosophila Genetic Reference Panel with two common insecticides; malathion (an organophosphate) and permethrin (a pyrethroid). Genome-wide association studies map survival on malathion to the two of the largest sweeps in the D. melanogaster genome; Ace and Cyp6g1 Malathion survivorship also correlates with lines which have high levels of Cyp12d1, Jheh1 and Jheh2 transcript abundance. Permethrin phenotypes map to the largest cluster of P450 genes in the Drosophila genome, however in contrast to a selective sweep driven by insecticide use, the derived allele seems to be associated with susceptibility. These results underscore previous findings that highlight the importance of structural variation to insecticide phenotypes: Cyp6g1 exhibits copy number variation and transposable element insertions, Cyp12d1 is tandemly duplicated, the Jheh loci are associated with a Bari1 transposable element insertion, and a Cyp6a17 deletion is associated with susceptibility

GWAS on multiple traits identifies mitochondrial ACONITASE3 as important for acclimation to submergence stress

Mitochondrial ACONITASE3 is important for the acclimation to submergence stress by integrating carbon and nitrogen metabolism and impacting stress signaling pathways. Flooding causes severe crop losses in many parts of the world. Genetic variation in flooding tolerance exists in many species; however, there are few examples for the identification of tolerance genes and their underlying function. We conducted a genome-wide association study (GWAS) in 387 Arabidopsis (Arabidopsis thaliana) accessions. Plants were subjected to prolonged submergence followed by desubmergence, and seven traits (score, water content, Fv/Fm, and concentrations of nitrate, chlorophyll, protein, and starch) were quantified to characterize their acclimation responses. These traits showed substantial variation across the range of accessions. A total of 35 highly significant single-nucleotide polymorphisms (SNPs) were identified across the 20 GWA datasets, pointing to 22 candidate genes, with functions in TCA cycle, DNA modification, and cell division. Detailed functional characterization of one candidate gene, ACONITASE3 (ACO3), was performed. Chromatin immunoprecipitation followed by sequencing showed that a single nucleotide polymorphism in the ACO3 promoter co-located with the binding site of the master regulator of retrograde signaling ANAC017, while subcellular localization of an ACO3-YFP fusion protein confirmed a mitochondrial localization during submergence. Analysis of mutant and overexpression lines determined changes in trait parameters that correlated with altered submergence tolerance and were consistent with the GWAS results. Subsequent RNA-seq experiments suggested that impairing ACO3 function increases the sensitivity to submergence by altering ethylene signaling, whereas ACO3 overexpression leads to tolerance by metabolic priming. These results indicate that ACO3 impacts submergence tolerance through integration of carbon and nitrogen metabolism via the mitochondrial TCA cycle and impacts stress signaling during acclimation to stress.Peer reviewe

Dissecting genetic architecture of grape proanthocyanidin composition through quantitative trait locus mapping

<p>Abstract</p> <p>Background</p> <p>Proanthocyanidins (PAs), or condensed tannins, are flavonoid polymers, widespread throughout the plant kingdom, which provide protection against herbivores while conferring organoleptic and nutritive values to plant-derived foods, such as wine. However, the genetic basis of qualitative and quantitative PA composition variation is still poorly understood. To elucidate the genetic architecture of the complex grape PA composition, we first carried out quantitative trait locus (QTL) analysis on a 191-individual pseudo-F1 progeny. Three categories of PA variables were assessed: total content, percentages of constitutive subunits and composite ratio variables. For nine functional candidate genes, among which eight co-located with QTLs, we performed association analyses using a diversity panel of 141 grapevine cultivars in order to identify causal SNPs.</p> <p>Results</p> <p>Multiple QTL analysis revealed a total of 103 and 43 QTLs, respectively for seed and skin PA variables. Loci were mainly of additive effect while some loci were primarily of dominant effect. Results also showed a large involvement of pairwise epistatic interactions in shaping PA composition. QTLs for PA variables in skin and seeds differed in number, position, involvement of epistatic interaction and allelic effect, thus revealing different genetic determinisms for grape PA composition in seeds and skin. Association results were consistent with QTL analyses in most cases: four out of nine tested candidate genes (<it>VvLAR1</it>, <it>VvMYBPA2</it>, <it>VvCHI1</it>, <it>VvMYBPA1</it>) showed at least one significant association with PA variables, especially <it>VvLAR1 </it>revealed as of great interest for further functional investigation. Some SNP-phenotype associations were observed only in the diversity panel.</p> <p>Conclusions</p> <p>This study presents the first QTL analysis on grape berry PA composition with a comparison between skin and seeds, together with an association study. Our results suggest a complex genetic control for PA traits and different genetic architectures for grape PA composition between berry skin and seeds. This work also uncovers novel genomic regions for further investigation in order to increase our knowledge of the genetic basis of PA composition.</p

Morphological features of distribution of branches of the ethmoid arteries on from the shape of the orbit

To study the topography of the lattice of the arteries of the labyrinth and their relationship to the form of the orbit was carried out anthropometric measurements by the method of V. N. Shevkunenko, A. M. Eselevich, T. V. Zolotareva, G. N. Toporov. To determine the shape of the orbit was measured by the following parameters: 1 - the height of the entrance to the eye socket - the distance between the upper and lower walls of the orbit; 2 - entrance width - the distance between the inner and outer walls of the orbit; 3 - the depth - the distance between the inner boundary of the entrance into the orbit to the optic canal

AtHMA4 drives natural variation in leaf Zn concentration of Arabidopsis thaliana

Zinc (Zn) is an essential element for plant growth and development, and Zn derived from crop plants in the diet is also important for human health. Here, we report that genetic variation in Heavy Metal-ATPase 4 (HMA4) controls natural variation in leaf Zn content. Investigation of the natural variation in leaf Zn content in a world-wide collection of 349 Arabidopsis thaliana wild collected accessions identified two accessions, Van-0 and Fab-2, which accumulate significantly lower Zn when compared with Col-0. Both quantitative trait loci (QTL) analysis and bulked segregant analysis (BSA) identified HMA4 as a strong candidate accounting for this variation in leaf Zn concentration. Genetic complementation experiments confirmed this hypothesis. Sequence analysis revealed that a 1-bp deletion in the third exon of HMA4 from Fab-2 is responsible for the lose of function of HMA4 driving the low Zn observed in Fab-2. Unlike in Fab-2 polymorphisms in the promoter region were found to be responsible for the weak function of HMA4 in Van-0. This is supported by both an expression analysis of HMA4 in Van-0 and through a series of T-DNA insertion mutants which generate truncated HMA4 promoters in the Col-0 background. In addition, we also observed that Fab-2, Van-0 and the hma4-2 null mutant in the Col-0 background show enhanced resistance to a combination of high Zn and high Cd in the growth medium, raising the possibility that variation at HMA4 may play a role in environmental adaptation

Complex patterns of local adaptation in teosinte

Populations of widely distributed species often encounter and adapt to specific environmental conditions. However, comprehensive characterization of the genetic basis of adaptation is demanding, requiring genome-wide genotype data, multiple sampled populations, and a good understanding of population structure. We have used environmental and high-density genotype data to describe the genetic basis of local adaptation in 21 populations of teosinte, the wild ancestor of maize. We found that altitude, dispersal events and admixture among subspecies formed a complex hierarchical genetic structure within teosinte. Patterns of linkage disequilibrium revealed four mega-base scale inversions that segregated among populations and had altitudinal clines. Based on patterns of differentiation and correlation with environmental variation, inversions and nongenic regions play an important role in local adaptation of teosinte. Further, we note that strongly differentiated individual populations can bias the identification of adaptive loci. The role of inversions in local adaptation has been predicted by theory and requires attention as genome-wide data become available for additional plant species. These results also suggest a potentially important role for noncoding variation, especially in large plant genomes in which the gene space represents a fraction of the entire genome

Repeated local emergence of carbapenem-resistant Acinetobacter baumannii in a single hospital ward.

We recently reported a dramatic increase in the prevalence of carbapenem-resistant Acinetobacter baumannii infections in the intensive care unit (ICU) of a Vietnamese hospital. This upsurge was associated with a specific oxa23-positive clone that was identified by multilocus VNTR analysis. Here, we used whole-genome sequence analysis to dissect the emergence of carbapenem-resistant A. baumannii causing ventilator-associated pneumonia (VAP) in the ICU during 2009-2012. To provide historical context and distinguish microevolution from strain introduction, we compared these genomes with those of A. baumannii asymptomatic carriage and VAP isolates from this same ICU collected during 2003-2007. We identified diverse lineages co-circulating over many years. Carbapenem resistance was associated with the presence of oxa23, oxa40, oxa58 and ndm1 genes in multiple lineages. The majority of resistant isolates were oxa23-positive global clone GC2; fine-scale phylogenomic analysis revealed five distinct GC2 sublineages within the ICU that had evolved locally via independent chromosomal insertions of oxa23 transposons. The increase in infections caused by carbapenem-resistant A. baumannii was associated with transposon-mediated transmission of a carbapenemase gene, rather than clonal expansion or spread of a carbapenemase-harbouring plasmid. Additionally, we found evidence of homologous recombination creating diversity within the local GC2 population, including several events resulting in replacement of the capsule locus. We identified likely donors of the imported capsule locus sequences amongst the A. baumannii isolated on the same ward, suggesting that diversification was largely facilitated via reassortment and sharing of genetic material within the localized A. baumannii population

Genetic diversity, linkage disequilibrium and power of a large grapevine (Vitis vinifera L) diversity panel newly designed for association studies

UMR-AGAP Equipe DAVV (Diversité, adaptation et amélioration de la vigne) ; équipe ID (Intégration de Données)International audienceAbstractBackgroundAs for many crops, new high-quality grapevine varieties requiring less pesticide and adapted to climate change are needed. In perennial species, breeding is a long process which can be speeded up by gaining knowledge about quantitative trait loci linked to agronomic traits variation. However, due to the long juvenile period of these species, establishing numerous highly recombinant populations for high resolution mapping is both costly and time-consuming. Genome wide association studies in germplasm panels is an alternative method of choice, since it allows identifying the main quantitative trait loci with high resolution by exploiting past recombination events between cultivars. Such studies require adequate panel design to represent most of the available genetic and phenotypic diversity. Assessing linkage disequilibrium extent and panel power is also needed to determine the marker density required for association studies.ResultsStarting from the largest grapevine collection worldwide maintained in Vassal (France), we designed a diversity panel of 279 cultivars with limited relatedness, reflecting the low structuration in three genetic pools resulting from different uses (table vs wine) and geographical origin (East vs West), and including the major founders of modern cultivars. With 20 simple sequence repeat markers and five quantitative traits, we showed that our panel adequately captured most of the genetic and phenotypic diversity existing within the entire Vassal collection. To assess linkage disequilibrium extent and panel power, we genotyped single nucleotide polymorphisms: 372 over four genomic regions and 129 distributed over the whole genome. Linkage disequilibrium, measured by correlation corrected for kinship, reached 0.2 for a physical distance between 9 and 458 Kb depending on genetic pool and genomic region, with varying size of linkage disequilibrium blocks. This panel achieved reasonable power to detect associations between traits with high broad-sense heritability (> 0.7) and causal loci with intermediate allelic frequency and strong effect (explaining > 10 % of total variance).ConclusionsOur association panel constitutes a new, highly valuable resource for genetic association studies in grapevine, and deserves dissemination to diverse field and greenhouse trials to gain more insight into the genetic control of many agronomic traits and their interaction with the environment