Novel Genes Affecting the Interaction between the Cabbage Whitefly and Arabidopsis Uncovered by Genome-Wide Association Mapping

Plants have evolved a variety of ways to defend themselves against biotic attackers. This has resulted in the presence of substantial variation in defense mechanisms among plants, even within a species. Genome-wide association (GWA) mapping is a useful tool to study the genetic architecture of traits, but has so far only had limited exploitation in studies of plant defense. Here, we study the genetic architecture of defense against the phloem-feeding insect cabbage whitefly (Aleyrodes proletella) in Arabidopsis thaliana. We determined whitefly performance, i.e. the survival and reproduction of whitefly females, on 360 worldwide selected natural accessions and subsequently performed GWA mapping using 214,051 SNPs. Substantial variation for whitefly adult survival and oviposition rate (number of eggs laid per female per day) was observed between the accessions. We identified 39 candidate SNPs for either whitefly adult survival or oviposition rate, all with relatively small effects, underpinning the complex architecture of defense traits. Among the corresponding candidate genes, i.e. genes in linkage disequilibrium (LD) with candidate SNPs, none have previously been identified as a gene playing a role in the interaction between plants and phloem-feeding insects. Whitefly performance on knock-out mutants of a number of candidate genes was significantly affected, validating the potential of GWA mapping for novel gene discovery in plant-insect interactions. Our results show that GWA analysis is a very useful tool to gain insight into the genetic architecture of plant defense against herbivorous insects, i.e. we identified and validated several genes affecting whitefly performance that have not previously been related to plant defense against herbivorous insects.     

Quantitative trait loci and candidate genes underlying genotype by environment interaction in the response of Arabidopsis thaliana to drought

Drought stress was imposed on two sets of Arabidopsis thaliana genotypes grown in sand under short‐day conditions and analysed for several shoot and root growth traits. The response to drought was assessed for quantitative trait locus (QTL) mapping in a genetically diverse set of Arabidopsis accessions using genome‐wide association (GWA) mapping, and conventional linkage analysis of a recombinant inbred line (RIL) population. Results showed significant genotype by environment interaction (G×E) for all traits in response to different watering regimes. For the RIL population, the observed G×E was reflected in 17 QTL by environment interactions (Q×E), while 17 additional QTLs were mapped not showing Q×E. GWA mapping identified 58 single nucleotide polymorphism (SNPs) associated with loci displaying Q×E and an additional 16 SNPs associated with loci not showing Q×E. Many candidate genes potentially underlying these loci were suggested. The genes for RPS3C and YLS7 were found to contain conserved amino acid differences when comparing Arabidopsis accessions with strongly contrasting drought response phenotypes, further supporting their candidacy. One of these candidate genes co‐located with a QTL mapped in the RIL population.     

Use of virus-induced gene silencing to characterize genes involved in modulating hypersensitive cell death in maize

Plant disease resistance proteins (R-proteins) detect specific pathogen-derived molecules, triggering a defence response often including a rapid localized cell death at the point of pathogen penetration called the hypersensitive response (HR). The maize Rp1-D21 gene encodes a protein that triggers a spontaneous HR causing spots on leaves in the absence of any pathogen. Previously, we used fine mapping and functional analysis in a Nicotiana benthamiana transient expression system to identify and characterize a number of genes associated with variation in Rp1-D21-induced HR. Here we describe a system for characterizing genes mediating HR, using virus-induced gene silencing (VIGS) in a maize line carrying Rp1-D21. We assess the roles of 12 candidate genes. Three of these genes, SGT1, RAR1, and HSP90, are required for HR induced by a number of R-proteins across several plant–pathogen systems. We confirmed that maize HSP90 was required for full Rp1-D21-induced HR. However, suppression of SGT1 expression unexpectedly increased the severity of Rp1-D21-induced HR while suppression of RAR1 expression had no measurable effect. We confirmed the effects on HR of two genes we had previously validated in the N. benthamiana system, hydroxycinnamoyltransferase and caffeoyl CoA O-methyltransferase. We further showed the suppression the expression of two previously uncharacterized, candidate genes, IQ calmodulin binding protein (IQM3) and vacuolar protein sorting protein 37, suppressed Rp1-D21-induced HR. This approach is an efficient way to characterize the roles of genes modulating the hypersensitive defence response and other dominant lesion phenotypes in maize.     

Confirmation and fine‐mapping of clinical mastitis and somatic cell score QTL in Nordic Holstein cattle

A genome‐wide association study of 2098 progeny‐tested Nordic Holstein bulls genotyped for 36 387 SNPs on 29 autosomes was conducted to confirm and fine‐map quantitative trait loci (QTL) for mastitis traits identified earlier using linkage analysis with sparse microsatellite markers in the same population. We used linear mixed model analysis where a polygenic genetic effect was fitted as a random effect and single SNPs were successively included as fixed effects in the model. We detected 143 SNP‐by‐trait significant associations (P < 0.0001) on 20 chromosomes affecting mastitis‐related traits. Among them, 21 SNP‐by‐trait combinations exceeded the genome‐wide significant threshold. For 12 chromosomes, both the present association study and the previous linkage study detected QTL, and of these, six were in the same chromosomal locations. Strong associations of SNPs with mastitis traits were observed on bovine autosomes 6, 13, 14 and 20. Possible candidate genes for these QTL were identified. Identification of SNPs in linkage disequilibrium with QTL will enable marker‐based selection for mastitis resistance. The candidate genes identified should be further studied to detect candidate polymorphisms underlying these QTL.     

Effector-Triggered Immune Response in Arabidopsis thaliana Is a Quantitative Trait

We identified loci responsible for natural variation in Arabidopsis thaliana (Arabidopsis) responses to a bacterial pathogen virulence factor, HopAM1. HopAM1 is a type III effector protein secreted by the virulent Pseudomonas syringae strain Pto DC3000. Delivery of HopAM1 from disarmed Pseudomonas strains leads to local cell death, meristem chlorosis, or both, with varying intensities in different Arabidopsis accessions. These phenotypes are not associated with differences in bacterial growth restriction. We treated the two phenotypes as quantitative traits to identify host loci controlling responses to HopAM1. Genome-wide association (GWA) of 64 Arabidopsis accessions identified independent variants highly correlated with response to each phenotype. Quantitative trait locus (QTL) mapping in a recombinant inbred population between Bur-0 and Col-0 accessions revealed genetic linkage to regions distinct from the top GWA hits. Two major QTL associated with HopAM1-induced cell death were also associated with HopAM1-induced chlorosis. HopAM1-induced changes in Arabidopsis gene expression showed that rapid HopAM1-dependent cell death in Bur-0 is correlated with effector-triggered immune responses. Studies of the effect of mutations in known plant immune system genes showed, surprisingly, that both cell death and chlorosis phenotypes are enhanced by loss of EDS1, a regulatory hub in the plant immune-signaling network. Our results reveal complex genetic architecture for response to this particular type III virulence effector, in contrast to the typical monogenic control of cell death and disease resistance triggered by most type III effectors.     

Genome-wide association study identifies Loci for body composition and structural soundness traits in pigs

Background

The recent completion of the swine genome sequencing project and development of a high density porcine SNP array has made genome-wide association (GWA) studies feasible in pigs.

Methodology/Principal Findings

Using Illumina''s PorcineSNP60 BeadChip, we performed a pilot GWA study in 820 commercial female pigs phenotyped for backfat, loin muscle area, body conformation in addition to feet and leg (FL) structural soundness traits. A total of 51,385 SNPs were jointly fitted using Bayesian techniques as random effects in a mixture model that assumed a known large proportion (99.5%) of SNPs had zero effect. SNP annotations were implemented through the Sus scrofa Build 9 available from pig Ensembl. We discovered a number of candidate chromosomal regions, and some of them corresponded to QTL regions previously reported. We not only have identified some well-known candidate genes for the traits of interest, such as MC4R (for backfat) and IGF2 (for loin muscle area), but also obtained novel promising genes, including CHCHD3 (for backfat), BMP2 (for loin muscle area, body size and several FL structure traits), and some HOXA family genes (for overall leg action). The candidate regions responsible for body conformation and FL structure soundness did not overlap greatly which implied that these traits were controlled by different genes. Functional clustering analyses classified the genes into categories related to bone and cartilage development, muscle growth and development or the insulin pathway suggesting the traits are regulated by common pathways or gene networks that exert roles at different spatial and temporal stages.

Conclusions/Significance

This study is one of the earliest GWA reports on important quantitative traits in pigs, and the findings will contribute to the further biological function analysis of the identified candidate genes and potential utilization of them in marker assisted selection.     

Association Mapping for Fruit,Plant and Leaf Morphology Traits in Eggplant

An eggplant (Solanum melongena) association panel of 191 accessions, comprising a mixture of breeding lines, old varieties and landrace selections was SNP genotyped and phenotyped for key breeding fruit and plant traits at two locations over two seasons. A genome-wide association (GWA) analysis was performed using the mixed linear model, which takes into account both a kinship matrix and the sub-population membership of the accessions. Overall, 194 phenotype/genotype associations were uncovered, relating to 30 of the 33 measured traits. These associations involved 79 SNP loci mapping to 39 distinct chromosomal regions distributed over all 12 eggplant chromosomes. A comparison of the map positions of these SNPs with those of loci derived from conventional linkage mapping showed that GWA analysis both validated many of the known controlling loci and detected a large number of new marker/trait associations. Exploiting established syntenic relationships between eggplant chromosomes and those of tomato and pepper recognized orthologous regions in ten eggplant chromosomes harbouring genes influencing breeders’ traits.     

Genome-wide association mapping of blood cell traits in mice

Genetic variations in blood cell parameters can impact clinical traits. We report here the mapping of blood cell traits in a panel of 100 inbred strains of mice of the Hybrid Mouse Diversity Panel (HMDP) using genome-wide association (GWA). We replicated a locus previously identified in using linkage analysis in several genetic crosses for mean corpuscular volume (MCV) and a number of other red blood cell traits on distal chromosome 7. Our peak for SNP association to MCV occurred in a linkage disequilibrium (LD) block spanning from 109.38 to 111.75 Mb that includes Hbb-b1, the likely causal gene. Altogether, we identified five loci controlling red blood cell traits (on chromosomes 1, 7, 11, 12, and 16), and four of these correspond to loci for red blood cell traits reported in a recent human GWA study. For white blood cells, including granulocytes, monocytes, and lymphocytes, a total of six significant loci were identified on chromosomes 1, 6, 8, 11, 12, and 15. An average of ten candidate genes were found at each locus and those were prioritized by examining functional variants in the HMDP such as missense and expression variants. These results provide intermediate phenotypes and candidate loci for genetic studies of atherosclerosis and cancer as well as inflammatory and immune disorders in mice.     

Use of Mutant-Assisted Gene Identification and Characterization (MAGIC) to identify novel genetic loci that modify the maize hypersensitive response

The partially dominant, autoactive maize disease resistance gene Rp1-D21 causes hypersensitive response (HR) lesions to form spontaneously on leaves and stems in the absence of pathogen recognition. The maize nested association mapping (NAM) population consists of 25 200-line subpopulations each derived from a cross between the maize line B73 and one of 25 diverse inbred lines. By crossing a line carrying the Rp1-D21 gene with lines from three of these subpopulations and assessing the F(1) progeny, we were able to map several novel loci that modify the maize HR, using both single-population quantitative trait locus (QTL) and joint analysis of all three populations. Joint analysis detected QTL in greater number and with greater confidence and precision than did single population analysis. In particular, QTL were detected in bins 1.02, 4.04, 9.03, and 10.03. We have previously termed this technique, in which a mutant phenotype is used as a "reporter" for a trait of interest, Mutant-Assisted Gene Identification and Characterization (MAGIC).     

Test for positional candidate genes for body composition on pig chromosome 6

One QTL affecting backfat thickness (BF), intramuscular fat content (IMF) and eye muscle area (MA) was previously localized on porcine chromosome 6 in an F₂ cross between Iberian and Landrace pigs. This work was done to study the effect of two positional candidate genes on these traits: H-FABP and LEPR genes. The QTL mapping analysis was repeated with a regression method using genotypes for seven microsatellites and two PCR-RFLPs in the H-FABP and LEPR genes. H-FABP and LEPR genes were located at 85.4 and 107 cM respectively, by linkage analysis. The effects of the candidate gene polymorphisms were analyzed in two ways. When an animal model was fitted, both genes showed significant effects on fatness traits, the H-FABP polymorphism showed significant effects on IMF and MA, and the LEPR polymorphism on BF and IMF. But when the candidate gene effect was included in a QTL regression analysis these associations were not observed, suggesting that they must not be the causal mutations responsible for the effects found. Differences in the results of both analyses showed the inadequacy of the animal model approach for the evaluation of positional candidate genes in populations with linkage disequilibrium, when the probabilities of the parental origin of the QTL alleles are not included in the model.     

Identification of Glu-B1-1 as a candidate gene for the quantity of high-molecular-weight glutenin in bread wheat (Triticum aestivum L.) by means of an association study

A previous study in wheat (Triticum aestivum L.) identified two candidate genes controlling a quantitative trait locus (QTL) for high-molecular-weight glutenin subunit (HMW-GS) GluBx. These candidates were Glu-B1-1, the structural gene coding for Glu1Bx, and the B homoeologous gene coding for SPA (spa-B), a seed storage protein activator. The goal of this study was to identify the best candidate gene for this QTL. Single nucleotide polymorphisms (SNPs) are an abundant source of DNA polymorphisms that have been successfully used to identify loci associated with particular phenotypes. As no linkage disequilibrium was detected between Glu-B1-1 and spa-B, we performed an association study to identify the individual gene responsible for the QTL. Six SNPs, three located in Glu-B1-1 and three in spa-B, were genotyped by mass spectrometry in a collection of 113 bread wheat lines. These lines were also evaluated for protein content as well as the total quantity of HMW-GSs and of each HMW-GS in seed samples from two harvest years. Significant associations were detected only between Glu-B1-1 polymorphism and most of the traits evaluated. Spa-B was unambiguously discarded as a candidate. To our knowledge, this is the first report on an association study that was successfully used to discriminate between two candidate genes.     

Molecular and Functional Analyses of a Maize Autoactive NB-LRR Protein Identify Precise Structural Requirements for Activity

Plant disease resistance is often mediated by nucleotide binding-leucine rich repeat (NLR) proteins which remain auto-inhibited until recognition of specific pathogen-derived molecules causes their activation, triggering a rapid, localized cell death called a hypersensitive response (HR). Three domains are recognized in one of the major classes of NLR proteins: a coiled-coil (CC), a nucleotide binding (NB-ARC) and a leucine rich repeat (LRR) domains. The maize NLR gene Rp1-D21 derives from an intergenic recombination event between two NLR genes, Rp1-D and Rp1-dp2 and confers an autoactive HR. We report systematic structural and functional analyses of Rp1 proteins in maize and N. benthamiana to characterize the molecular mechanism of NLR activation/auto-inhibition. We derive a model comprising the following three main features: Rp1 proteins appear to self-associate to become competent for activity. The CC domain is signaling-competent and is sufficient to induce HR. This can be suppressed by the NB-ARC domain through direct interaction. In autoactive proteins, the interaction of the LRR domain with the NB-ARC domain causes de-repression and thus disrupts the inhibition of HR. Further, we identify specific amino acids and combinations thereof that are important for the auto-inhibition/activity of Rp1 proteins. We also provide evidence for the function of MHD2, a previously uncharacterized, though widely conserved NLR motif. This work reports several novel insights into the precise structural requirement for NLR function and informs efforts towards utilizing these proteins for engineering disease resistance.     

GWA Analysis for Milk Production Traits in Dairy Sheep and Genetic Support for a QTN Influencing Milk Protein Percentage in the LALBA Gene

In this study, we used the Illumina OvineSNP50 BeadChip to conduct a genome-wide association (GWA) analysis for milk production traits in dairy sheep by analyzing a commercial population of Spanish Churra sheep. The studied population consisted of a total of 1,681 Churra ewes belonging to 16 half-sib families with available records for milk yield (MY), milk protein and fat yields (PY and FY) and milk protein and fat contents (PP and FP). The most significant association identified reached experiment-wise significance for PP and FP and was located on chromosome 3 (OAR3). These results confirm the population-level segregation of a previously reported QTL affecting PP and suggest that this QTL has a significant pleiotropic effect on FP. Further associations were detected at the chromosome-wise significance level on 14 other chromosomal regions. The marker on OAR3 showing the highest significant association was located at the third intron of the alpha-lactalbumin (LALBA) gene, which is a functional and positional candidate underlying this association. Sequencing this gene in the 16 Churra rams of the studied resource population identified additional polymorphisms. One out of the 31 polymorphisms identified was located within the coding gene sequence (LALBA_g.242T>C) and was predicted to cause an amino acid change in the protein (Val27Ala). Different approaches, including GWA analysis, a combined linkage and linkage disequilibrium study and a concordance test with the QTL segregating status of the sires, were utilized to assess the role of this mutation as a putative QTN for the genetic effects detected on OAR3. Our results strongly support the polymorphism LALBA_g.242T>C as the most likely causal mutation of the studied OAR3 QTL affecting PP and FP, although we cannot rule out the possibility that this SNP is in perfect linkage disequilibrium with the true causal polymorphism.