Deploying QTL-seq for rapid delineation of a potential candidate gene underlying major trait-associated QTL in chickpea

A rapid high-resolution genome-wide strategy for molecular mapping of major QTL(s)/gene(s) regulating important agronomic traits is vital for in-depth dissection of complex quantitative traits and genetic enhancement in chickpea. The present study for the first time employed a NGS-based whole-genome QTL-seq strategy to identify one major genomic region harbouring a robust 100-seed weight QTL using an intra-specific 221 chickpea mapping population (desi cv. ICC 7184 × desi cv. ICC 15061). The QTL-seq-derived major SW QTL (CaqSW1.1) was further validated by single-nucleotide polymorphism (SNP) and simple sequence repeat (SSR) marker-based traditional QTL mapping (47.6% R² at higher LOD >19). This reflects the reliability and efficacy of QTL-seq as a strategy for rapid genome-wide scanning and fine mapping of major trait regulatory QTLs in chickpea. The use of QTL-seq and classical QTL mapping in combination narrowed down the 1.37 Mb (comprising 177 genes) major SW QTL (CaqSW1.1) region into a 35 kb genomic interval on desi chickpea chromosome 1 containing six genes. One coding SNP (G/A)-carrying constitutive photomorphogenic9 (COP9) signalosome complex subunit 8 (CSN8) gene of these exhibited seed-specific expression, including pronounced differential up-/down-regulation in low and high seed weight mapping parents and homozygous individuals during seed development. The coding SNP mined in this potential seed weight-governing candidate CSN8 gene was found to be present exclusively in all cultivated species/genotypes, but not in any wild species/genotypes of primary, secondary and tertiary gene pools. This indicates the effect of strong artificial and/or natural selection pressure on target SW locus during chickpea domestication. The proposed QTL-seq-driven integrated genome-wide strategy has potential to delineate major candidate gene(s) harbouring a robust trait regulatory QTL rapidly with optimal use of resources. This will further assist us to extrapolate the molecular mechanism underlying complex quantitative traits at a genome-wide scale leading to fast-paced marker-assisted genetic improvement in diverse crop plants, including chickpea.     

Quantitative trait locus analysis of susceptibility to diet-induced atherosclerosis in recombinant inbred mice

Quantitative trait locus (QTL) analysis is a statistical method that can be applied to identify loci making a significant impact on a phenotype. For the phenotype of susceptibility to diet-induced atherosclerosis in the mouse, we have studied four quantitative traits: area of aortic fatty streaks and serum concentrations of high-density lipoprotein-bound cholesterol (HDL-cholesterol), apolipoprotein A-I, and apolipoprotein A-II (apo A-II). QTL analysis revealed a significant locus on chromosome 1 distal impacting serum apo A-II concentration on a high-fat diet and serum HDL-cholesterol concentration on a chow diet. This locus is presumablyApoa-2, the structural gene for apo A-II. QTL analysis of aortic fatty streaks failed to reveal a significant locus.     

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.     

Potential for detection of marker–trait associations in durum wheat using unbalanced, historical phenotypic datasets

Identification of marker?Ctrait associations in germplasm relevant to a breeding program can be an effective way to identify quantitative trait loci (QTL) useful for selection and is critical to the success of genome-wide selection strategies. This approach is most cost-effective if phenotypic data routinely collected by breeding programs is used, necessitating only addition of genotypic data. The objective of this work was to evaluate such an approach using unbalanced phenotypic data from durum wheat (Triticum turgidum L. var. durum) registration trials genotyped with diversity arrays technology (DArT) markers. Plant height, grain cadmium concentration and yellow colour loss during pasta manufacture were chosen as example traits because all are influenced by major genes associated with known QTL. A further evaluation was performed on semolina yellow pigment concentration, a more complexly-inherited trait, but with numerous QTL identified. In total, 870 informative DArT markers were used to detect marker?Ctrait associations. The genome coverage of markers was uneven, with low coverage of chromosomes 4B and 5A. The DArT coverage of chromosome 4B was too sparse to identify markers strongly associated with the semidwarf height locus Rht-B1 and the lipoxygenase locus Lpx-B1, both known to reside on 4B. The 20 DArT markers associated with pigment concentration localized to chromosomes 1B, 2A, 5B, 6A, 7A and 7B, linked to the trait in other studies. One DArT clone showed sequence identity to a single wheat expressed sequence tag that maps to the same deletion bin as Psy1-A1, a gene previously associated with yellow pigment concentration in durum wheat. Three markers were associated with grain cadmium and explained similar proportions of the phenotypic variance as the Xusw14 marker known to be physically linked to Cdu-B1, a major locus on 5B regulating cadmium accumulation. The sequences of these three DArT markers were 98?% identical, and were used to identify a single gene in rice that is physically linked to other rice genes that co-localize with Cdu-B1 in durum wheat. The results suggest that this historical phenotypic dataset is useful for QTL discovery and would potentially be a ??training population?? for genomic selection when a high-density, low-cost marker platform becomes available.     

Genome-wide analysis reveals a complex pattern of genomic imprinting in mice

Parent-of-origin–dependent gene expression resulting from genomic imprinting plays an important role in modulating complex traits ranging from developmental processes to cognitive abilities and associated disorders. However, while gene-targeting techniques have allowed for the identification of imprinted loci, very little is known about the contribution of imprinting to quantitative variation in complex traits. Most studies, furthermore, assume a simple pattern of imprinting, resulting in either paternal or maternal gene expression; yet, more complex patterns of effects also exist. As a result, the distribution and number of different imprinting patterns across the genome remain largely unexplored. We address these unresolved issues using a genome-wide scan for imprinted quantitative trait loci (iQTL) affecting body weight and growth in mice using a novel three-generation design. We identified ten iQTL that display much more complex and diverse effect patterns than previously assumed, including four loci with effects similar to the callipyge mutation found in sheep. Three loci display a new phenotypic pattern that we refer to as bipolar dominance, where the two heterozygotes are different from each other while the two homozygotes are identical to each other. Our study furthermore detected a paternally expressed iQTL on Chromosome 7 in a region containing a known imprinting cluster with many paternally expressed genes. Surprisingly, the effects of the iQTL were mostly restricted to traits expressed after weaning. Our results imply that the quantitative effects of an imprinted allele at a locus depend both on its parent of origin and the allele it is paired with. Our findings also show that the imprinting pattern of a locus can be variable over ontogenetic time and, in contrast to current views, may often be stronger at later stages in life.     

Genome-wide Association Study and Meta-Analysis Identify ISL1 as Genome-wide Significant Susceptibility Gene for Bladder Exstrophy

The bladder exstrophy-epispadias complex (BEEC) represents the severe end of the uro-rectal malformation spectrum, and is thought to result from aberrant embryonic morphogenesis of the cloacal membrane and the urorectal septum. The most common form of BEEC is isolated classic bladder exstrophy (CBE). To identify susceptibility loci for CBE, we performed a genome-wide association study (GWAS) of 110 CBE patients and 1,177 controls of European origin. Here, an association was found with a region of approximately 220kb on chromosome 5q11.1. This region harbors the ISL1 (ISL LIM homeobox 1) gene. Multiple markers in this region showed evidence for association with CBE, including 84 markers with genome-wide significance. We then performed a meta-analysis using data from a previous GWAS by our group of 98 CBE patients and 526 controls of European origin. This meta-analysis also implicated the 5q11.1 locus in CBE risk. A total of 138 markers at this locus reached genome-wide significance in the meta-analysis, and the most significant marker (rs9291768) achieved a P value of 2.13 × 10⁻¹². No other locus in the meta-analysis achieved genome-wide significance. We then performed murine expression analyses to follow up this finding. Here, Isl1 expression was detected in the genital region within the critical time frame for human CBE development. Genital regions with Isl1 expression included the peri-cloacal mesenchyme and the urorectal septum. The present study identified the first genome-wide significant locus for CBE at chromosomal region 5q11.1, and provides strong evidence for the hypothesis that ISL1 is the responsible candidate gene in this region.     

New investigations around CYP11A1 and its possible involvement in an androstenone QTL characterised in Large White pigs

Background

Previously, in boars with extreme androstenone levels, differential expression of the CYP11A1 gene in the testes has been characterised. CYP11A1 is located in a region where a QTL influencing boar fat androstenone levels has been detected in a Large White pig population. Clarifying the role of CYP11A1 in boar taint is important because it catalyses the initial step of androstenone synthesis and also of steroid synthesis.

Results

A genome-wide association study located CYP11A1 at approximately 1300 kb upstream from SNP H3GA0021967, defining the centre of the region containing the QTL for androstenone variation. In this study, we partially sequenced the CYP11A1 gene and identified several new single nucleotide polymorphisms (SNP) within it. Characterisation of one animal, heterozygous for CYP11A1 testicular expression but homozygous for a haplotype of a large region containing CYP11A1, revealed that variation of CYP11A1 expression is probably regulated by a mutation located downstream from the SNP H3GA0021967. We analysed CYP11A1 expression in LW families according to haplotypes of the QTL region''s centre. Effects of haplotypes on CYP11A1 expression and on androstenone accumulation were not concordant.

Conclusion

This study shows that testicular expression of CYP11A1 is not solely responsible for the QTL influencing boar fat androstenone levels. As a conclusion, we propose to refute the hypothesis that a single mutation located near the centre of the QTL region could control androstenone accumulation in fat by regulating the CYP11A1 expression.     

The Calcitonin Receptor Gene Is a Candidate for Regulation of Susceptibility to Herpes simplex Type 1 Neuronal Infection Leading to Encephalitis in Rat

Herpes simplex encephalitis (HSE) is a fatal infection of the central nervous system (CNS) predominantly caused by Herpes simplex virus type 1. Factors regulating the susceptibility to HSE are still largely unknown. To identify host gene(s) regulating HSE susceptibility we performed a genome-wide linkage scan in an intercross between the susceptible DA and the resistant PVG rat. We found one major quantitative trait locus (QTL), Hse1, on rat chromosome 4 (confidence interval 24.3–31 Mb; LOD score 29.5) governing disease susceptibility. Fine mapping of Hse1 using recombinants, haplotype mapping and sequencing, as well as expression analysis of all genes in the interval identified the calcitonin receptor gene (Calcr) as the main candidate, which also is supported by functional studies. Thus, using unbiased genetic approach variability in Calcr was identified as potentially critical for infection and viral spread to the CNS and subsequent HSE development.     

A 12 kb multi-allelic copy number variation encompassing a GC gene enhancer is associated with mastitis resistance in dairy cattle

Clinical mastitis (CM) is an inflammatory disease occurring in the mammary glands of lactating cows. CM is under genetic control, and a prominent CM resistance QTL located on chromosome 6 was reported in various dairy cattle breeds. Nevertheless, the biological mechanism underpinning this QTL has been lacking. Herein, we mapped, fine-mapped, and discovered the putative causal variant underlying this CM resistance QTL in the Dutch dairy cattle population. We identified a ~12 kb multi-allelic copy number variant (CNV), that is in perfect linkage disequilibrium with a lead SNP, as a promising candidate variant. By implementing a fine-mapping and through expression QTL mapping, we showed that the group-specific component gene (GC), a gene encoding a vitamin D binding protein, is an excellent candidate causal gene for the QTL. The multiplicated alleles are associated with increased GC expression and low CM resistance. Ample evidence from functional genomics data supports the presence of an enhancer within this CNV, which would exert cis-regulatory effect on GC. We observed that strong positive selection swept the region near the CNV, and haplotypes associated with the multiplicated allele were strongly selected for. Moreover, the multiplicated allele showed pleiotropic effects for increased milk yield and reduced fertility, hinting that a shared underlying biology for these effects may revolve around the vitamin D pathway. These findings together suggest a putative causal variant of a CM resistance QTL, where a cis-regulatory element located within a CNV can alter gene expression and affect multiple economically important traits.     

Genome-wide association analyses reveal significant loci and strong candidate genes for growth and fatness traits in two pig populations

Background

Recently, genome-wide association studies (GWAS) have been reported on various pig traits. We performed a GWAS to analyze 22 traits related to growth and fatness on two pig populations: a White Duroc × Erhualian F₂ intercross population and a Chinese Sutai half-sib population.

Results

We identified 14 and 39 loci that displayed significant associations with growth and fatness traits at the genome-wide level and chromosome-wide level, respectively. The strongest association was between a 750 kb region on SSC7 (SSC for Sus scrofa) and backfat thickness at the first rib. This region had pleiotropic effects on both fatness and growth traits in F₂ animals and contained a promising candidate gene HMGA1 (high mobility group AT-hook 1). Unexpectedly, population genetic analysis revealed that the allele at this locus that reduces fatness and increases growth is derived from Chinese indigenous pigs and segregates in multiple Chinese breeds. The second strongest association was between the region around 82.85 Mb on SSC4 and average backfat thickness. PLAG1 (pleiomorphic adenoma gene 1), a gene under strong selection in European domestic pigs, is proximal to the top SNP and stands out as a strong candidate gene. On SSC2, a locus that significantly affects fatness traits mapped to the region around the IGF2 (insulin-like growth factor 2) gene but its non-imprinting inheritance excluded IGF2 as a candidate gene. A significant locus was also detected within a recombination cold spot that spans more than 30 Mb on SSCX, which hampered the identification of plausible candidate genes. Notably, no genome-wide significant locus was shared by the two experimental populations; different loci were observed that had both constant and time-specific effects on growth traits at different stages, which illustrates the complex genetic architecture of these traits.

Conclusions

We confirm several previously reported QTL and provide a list of novel loci for porcine growth and fatness traits in two experimental populations with Chinese Taihu and Western pigs as common founders. We showed that distinct loci exist for these traits in the two populations and identified HMGA1 and PLAG1 as strong candidate genes on SSC7 and SSC4, respectively.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-015-0089-5) contains supplementary material, which is available to authorized users.     

Genetic Regulation of Zfp30, CXCL1, and Neutrophilic Inflammation in Murine Lung

Allergic asthma is a complex disease characterized in part by granulocytic inflammation of the airways. In addition to eosinophils, neutrophils (PMN) are also present, particularly in cases of severe asthma. We sought to identify the genetic determinants of neutrophilic inflammation in a mouse model of house dust mite (HDM)-induced asthma. We applied an HDM model of allergic asthma to the eight founder strains of the Collaborative Cross (CC) and 151 incipient lines of the CC (preCC). Lung lavage fluid was analyzed for PMN count and the concentration of CXCL1, a hallmark PMN chemokine. PMN and CXCL1 were strongly correlated in preCC mice. We used quantitative trait locus (QTL) mapping to identify three variants affecting PMN, one of which colocalized with a QTL for CXCL1 on chromosome (Chr) 7. We used lung eQTL data to implicate a variant in the gene Zfp30 in the CXCL1/PMN response. This genetic variant regulates both CXCL1 and PMN by altering Zfp30 expression, and we model the relationships between the QTL and these three endophenotypes. We show that Zfp30 is expressed in airway epithelia in the normal mouse lung and that altering Zfp30 expression in vitro affects CXCL1 responses to an immune stimulus. Our results provide strong evidence that Zfp30 is a novel regulator of neutrophilic airway inflammation.     

Genetic mapping of a fusarium wilt resistance gene in Brassica oleracea

Fusarium wilt caused by Fusarium oxysporum f. sp. conglutinans is one of the most important diseases of Brassica crops, resulting in severe reductions in yield and quality. To characterize the inheritance pattern of fusarium resistance, a cross between a susceptible broccoli and a resistant cabbage was subjected to segregation analysis. Results indicated that resistance was controlled by a single dominant allele. This gene was named Foc-Bo1 and mapped to linkage group seven (O7) by both the segregation test and quantitative trait locus (QTL) analysis. The QTL on O7 was detected with a logarithm of odds score (LOD) of 19.5, which was above the threshold value with genome-wide 1% significance level (2.01). A minor QTL was also detected on O4 with a LOD score of 2.06. Inoculation tests indicated that stable expression of fusarium resistance at high temperatures required Foc-Bo1 homozygosity. The association between Foc-Bo1 and the closest simple sequence repeat marker (KBrS003O1N10) was analyzed in three F3 populations. Based on these studies, KBrS003O1N10 represents an effective marker-assisted selection (MAS) tool for breeding fusarium wilt resistance into Brassica oleracea crops. To our knowledge, this is the first paper to map the fusarium-resistance gene in Brassica species and to validate the effectiveness of MAS in improving fusarium resistance in these important plants.     

Powerful Identification of Cis-regulatory SNPs in Human Primary Monocytes Using Allele-Specific Gene Expression

A large number of genome-wide association studies have been performed during the past five years to identify associations between SNPs and human complex diseases and traits. The assignment of a functional role for the identified disease-associated SNP is not straight-forward. Genome-wide expression quantitative trait locus (eQTL) analysis is frequently used as the initial step to define a function while allele-specific gene expression (ASE) analysis has not yet gained a wide-spread use in disease mapping studies. We compared the power to identify cis-acting regulatory SNPs (cis-rSNPs) by genome-wide allele-specific gene expression (ASE) analysis with that of traditional expression quantitative trait locus (eQTL) mapping. Our study included 395 healthy blood donors for whom global gene expression profiles in circulating monocytes were determined by Illumina BeadArrays. ASE was assessed in a subset of these monocytes from 188 donors by quantitative genotyping of mRNA using a genome-wide panel of SNP markers. The performance of the two methods for detecting cis-rSNPs was evaluated by comparing associations between SNP genotypes and gene expression levels in sample sets of varying size. We found that up to 8-fold more samples are required for eQTL mapping to reach the same statistical power as that obtained by ASE analysis for the same rSNPs. The performance of ASE is insensitive to SNPs with low minor allele frequencies and detects a larger number of significantly associated rSNPs using the same sample size as eQTL mapping. An unequivocal conclusion from our comparison is that ASE analysis is more sensitive for detecting cis-rSNPs than standard eQTL mapping. Our study shows the potential of ASE mapping in tissue samples and primary cells which are difficult to obtain in large numbers.     

Genome-Wide Analyses Suggest Mechanisms Involving Early B-Cell Development in Canine IgA Deficiency

Immunoglobulin A deficiency (IgAD) is the most common primary immune deficiency disorder in both humans and dogs, characterized by recurrent mucosal tract infections and a predisposition for allergic and other immune mediated diseases. In several dog breeds, low IgA levels have been observed at a high frequency and with a clinical resemblance to human IgAD. In this study, we used genome-wide association studies (GWAS) to identify genomic regions associated with low IgA levels in dogs as a comparative model for human IgAD. We used a novel percentile groups-approach to establish breed-specific cut-offs and to perform analyses in a close to continuous manner. GWAS performed in four breeds prone to low IgA levels (German shepherd, Golden retriever, Labrador retriever and Shar-Pei) identified 35 genomic loci suggestively associated (p <0.0005) to IgA levels. In German shepherd, three genomic regions (candidate genes include KIRREL3 and SERPINA9) were genome-wide significantly associated (p <0.0002) with IgA levels. A ~20kb long haplotype on CFA28, significantly associated (p = 0.0005) to IgA levels in Shar-Pei, was positioned within the first intron of the gene SLIT1. Both KIRREL3 and SLIT1 are highly expressed in the central nervous system and in bone marrow and are potentially important during B-cell development. SERPINA9 expression is restricted to B-cells and peaks at the time-point when B-cells proliferate into antibody-producing plasma cells. The suggestively associated regions were enriched for genes in Gene Ontology gene sets involving inflammation and early immune cell development.     

A quantitative trait variant in Gabra2 underlies increased methamphetamine stimulant sensitivity

Psychostimulant (methamphetamine, cocaine) use disorders have a genetic component that remains mostly unknown. We conducted genome-wide quantitative trait locus (QTL) analysis of methamphetamine stimulant sensitivity. To facilitate gene identification, we employed a Reduced Complexity Cross between closely related C57BL/6 mouse substrains and examined maximum speed and distance traveled over 30 min following methamphetamine (2 mg/kg, i.p.). For maximum methamphetamine-induced speed following the second and third administration, we identified a single genome-wide significant QTL on chromosome 11 that peaked near the Cyfip2 locus (LOD = 3.5, 4.2; peak = 21 cM [36 Mb]). For methamphetamine-induced distance traveled following the first and second administration, we identified a genome-wide significant QTL on chromosome 5 that peaked near a functional intronic indel in Gabra2 coding for the alpha-2 subunit of the GABA-A receptor (LOD = 3.6–5.2; peak = 34–35 cM [66–67 Mb]). Striatal cis-expression QTL mapping corroborated Gabra2 as a functional candidate gene underlying methamphetamine-induced distance traveled. CRISPR/Cas9-mediated correction of the mutant intronic deletion on the C57BL/6J background to the wild-type C57BL/6NJ allele was sufficient to reduce methamphetamine-induced locomotor activity toward the wild-type C57BL/6NJ-like level, thus validating the quantitative trait variant (QTV). These studies show the power and efficiency of Reduced Complexity Crosses in identifying causal variants underlying complex traits. Functionally restoring Gabra2 expression decreased methamphetamine stimulant sensitivity and supports preclinical and human genetic studies implicating the GABA-A receptor in psychostimulant addiction-relevant traits. Importantly, our findings have major implications for studying psychostimulants in the C57BL/6J strain—the gold standard strain in biomedical research.