QTL for plant growth and morphology

Individuals and strains within most species exhibit heritable, and sometimes dramatic, differences in their growth and morphology. Advances in marker genotyping and statistical methods have increased the precision and sensitivity with which the quantitative trait loci (QTL) that are responsible for these differences can be mapped. This has resulted in both a more refined picture of the genetic architecture of many growth traits and the cloning of several of the genes that underlie plant QTL.     

Quantitative trait locus mapping reveals complex genetic architecture of quantitative virulence in the wheat pathogen Zymoseptoria tritici

We conducted a comprehensive analysis of virulence in the fungal wheat pathogen Zymoseptoria tritici using quantitative trait locus (QTL) mapping. High‐throughput phenotyping based on automated image analysis allowed the measurement of pathogen virulence on a scale and with a precision that was not previously possible. Across two mapping populations encompassing more than 520 progeny, 540 710 pycnidia were counted and their sizes and grey values were measured. A significant correlation was found between pycnidia size and both spore size and number. Precise measurements of percentage leaf area covered by lesions provided a quantitative measure of host damage. Combining these large and accurate phenotypic datasets with a dense panel of restriction site‐associated DNA sequencing (RADseq) genetic markers enabled us to genetically dissect pathogen virulence into components related to host damage and those related to pathogen reproduction. We showed that different components of virulence can be under separate genetic control. Large‐ and small‐effect QTLs were identified for all traits, with some QTLs specific to mapping populations, cultivars and traits and other QTLs shared among traits within the same mapping population. We associated the presence of four accessory chromosomes with small, but significant, increases in several virulence traits, providing the first evidence for a meaningful function associated with accessory chromosomes in this organism. A large‐effect QTL involved in host specialization was identified on chromosome 7, leading to the identification of candidate genes having a large effect on virulence.     

Tests of candidate genes in breed cross populations for QTL mapping in livestock

In recent years, several F₂ crosses between outbred lines of livestock have been developed to identify quantitative trait loci (QTL). These populations are valuable for further genetic analysis, including positional candidate gene loci (CGL). Analysis of CGL in F₂ populations is, however, hindered by extensive between-breed linkage disequilibrium (LD). The objectives here were to develop and evaluate three tests for CGL in simulated F₂ breed-cross populations. 1) A standard association test, based on the fixed effect of CGL genotype. This test was significant for CGL at considerable distances from the QTL. 2) A marker-assisted association test, based on a test at the CGL of the fixed effect of CGL genotype in a breed-cross QTL interval mapping model. This removed the impact of between-breed LD, but was not powerful in detecting CGL closely linked to the QTL, unless the CGL was the QTL. 3) An F-drop test, comparing F ratios for a QTL at the CGL with and without the CGL included as fixed effect. It had low power to distinguish close from distant CGL. Power to distinguish two CGL within 10 cM from the QTL was limited and little improved by including QTL effects associated with markers to remove between-breed LD, although the power was greater when one of the CGL was the causative mutation. Therefore, while we conclude that candidate gene tests in QTL mapping populations must be interpreted with caution, we now have a clearer picture of the value of candidate gene tests in these populations.     

QTL mapping and GWAS reveal candidate genes controlling capsaicinoid content in Capsicum

Capsaicinoids are unique compounds produced only in peppers (Capsicum spp.). Several studies using classical quantitative trait loci (QTLs) mapping and genomewide association studies (GWAS) have identified QTLs controlling capsaicinoid content in peppers; however, neither the QTLs common to each population nor the candidate genes underlying them have been identified due to the limitations of each approach used. Here, we performed QTL mapping and GWAS for capsaicinoid content in peppers using two recombinant inbred line (RIL) populations and one GWAS population. Whole‐genome resequencing and genotyping by sequencing (GBS) were used to construct high‐density single nucleotide polymorphism (SNP) maps. Five QTL regions on chromosomes 1, 2, 3, 4 and 10 were commonly identified in both RIL populations over multiple locations and years. Furthermore, a total of 109 610 SNPs derived from two GBS libraries were used to analyse the GWAS population consisting of 208 C. annuum‐clade accessions. A total of 69 QTL regions were identified from the GWAS, 10 of which were co‐located with the QTLs identified from the two biparental populations. Within these regions, we were able to identify five candidate genes known to be involved in capsaicinoid biosynthesis. Our results demonstrate that QTL mapping and GBS‐GWAS represent a powerful combined approach for the identification of loci controlling complex traits.     

QTL and candidate genes for fecundity in sows

Fecundity in pigs is a trait of major economic interest but low heritability. For the improvement of fecundity, genetic markers for selection are desirable and therefore, several searches for genetic variation influencing fecundity have been performed. The aim of this review is to compare and to evaluate all published QTL analyses and candidate gene approaches concerning reproductive traits in sows. For this purpose, we present a comprehensive cytogenetic map comprising 54 QTL and 11 candidate genes with influence on reproductive traits. The evaluation and comparison of the results showed similarities, but also marked differences among studies. Reasons for different results are multicausal and are due to differences between resource populations, number of evaluated animals, mating systems, measured phenotypical traits and environmental influences. We could show that chromosome 8 and to a lower extend chromosome 7 are the most important chromosomes with regard to reproductive traits in pigs. For further research, fine mapping of the identified QTL regions is necessary in order to confirm and to narrow the most likely chromosomal intervals. Although difficult to perform, an advance would be a standardization of the experimental setup in particular, in respect to the collection of phenotypic data. Furthermore, we suggest to publish the information on further identified QTL and candidate genes as comprehensive and accurate as possible in order to allow a more transparent comparison and collation of the results.     

Small RNAs from the plant pathogenic fungus Sclerotinia sclerotiorum highlight host candidate genes associated with quantitative disease resistance

Fungal plant pathogens secrete effector proteins and metabolites to cause disease. Additionally, some species transfer small RNAs (sRNAs) into plant cells to silence host mRNAs through complementary base pairing and suppress plant immunity. The fungus Sclerotinia sclerotiorum infects over 600 plant species, but little is known about the molecular processes that govern interactions with its many hosts. In particular, evidence for the production of sRNAs by S. sclerotiorum during infection is lacking. We sequenced sRNAs produced by S. sclerotiorum in vitro and during infection of two host species, Arabidopsis thaliana and Phaseolus vulgaris. We found that S. sclerotiorum produces at least 374 distinct highly abundant sRNAs during infection, mostly originating from repeat-rich plastic genomic regions. We predicted the targets of these sRNAs in A. thaliana and found that these genes were significantly more down-regulated during infection than the rest of the genome. Predicted targets of S. sclerotiorum sRNAs in A. thaliana were enriched for functional domains associated with plant immunity and were more strongly associated with quantitative disease resistance in a genome-wide association study (GWAS) than the rest of the genome. Mutants in A. thaliana predicted sRNA target genes SERK2 and SNAK2 were more susceptible to S. sclerotiorum than wild-type, suggesting that S. sclerotiorum sRNAs may contribute to the silencing of immune components in plants. The prediction of fungal sRNA targets in plant genomes can be combined with other global approaches, such as GWAS, to assist in the identification of plant genes involved in quantitative disease resistance.     

Genetic mapping of sulfur assimilation genes reveals a QTL for onion bulb pungency

Onion exhibits wide genetic and environmental variation in bioactive organosulfur compounds that impart pungency and health benefits. A PCR-based molecular marker map that included candidate genes for sulfur assimilation was used to identify genomic regions affecting pungency in the cross 'W202A' × 'Texas Grano 438'. Linkage mapping revealed that genes encoding plastidic ferredoxin-sulfite reductase (SiR) and plastidic ATP sulfurylase (ATPS) are closely linked (1–2 cM) on chromosome 3. Inbred F₃ families derived from the F₂ population used to construct the genetic map were grown in replicated trials in two environments and bulb pungency was evaluated as pyruvic acid or lachrymatory factor. Broad-sense heritability of pungency was estimated to be 0.78–0.80. QTL analysis revealed significant associations of both pungency and bulb soluble solids content with marker intervals on chromosomes 3 and 5, which have previously been reported to condition pleiotropic effects on bulb carbohydrate composition. Highly significant associations (LOD 3.7–8.7) were observed between ATPS and SiR Loci and bulb pungency but not with bulb solids content. This association was confirmed in two larger, independently derived F₂ families from the same cross. Single-locus models suggested that the partially dominant locus associated with these candidate genes controls 30–50% of genetic variation in pungency in these pedigrees. These markers may provide a practical means to select for lower pungency without correlated selection for lowered solids.     

Saturation mapping of QTL regions and identification of putative candidate genes for drought tolerance in rice

We have developed 85 new markers (50 RFLPs, 5 SSRs, 12 DD cDNAs, 9 ESTs, 8 HSP-encoding cDNAs and one BSA-derived AFLP marker) for saturation mapping of QTL regions for drought tolerance in rice, in our efforts to identify putative candidate genes. Thirteen of the markers were localized in the close vicinity of the targeted QTL regions. Fifteen of the additional markers mapped, respectively, inside one QTL region controlling osmotic adjustment on chromosome 3 ( oa3.1) and 14 regions that affect root traits on chromosomes 1, 2, 4, 5, 6, 7, 8, 9, 10 and 12. Differential display was used to identify more putative candidate genes and to saturate the QTL regions of the genetic map. Eleven of the isolated cDNA clones were found to be derived from drought-inducible genes. Two of them were unique and did not match any genes in the GenBank, while nine were highly similar to cDNAs encoding known proteins, including a DnaJ-related protein, a zinc-finger protein, a protease inhibitor, a glutathione-S-transferase, a DNA recombinase, and a protease. Twelve new cDNA fragments were mapped onto the genetic linkage map; seven of these mapped inside, or in close proximity to, the targeted QTL regions determining root thickness and osmotic adjustment capacity. The gene I12A1, which codes for a UDP-glucose 4-epimerase homolog, was identified as a putative target gene within the prt7.1/brt7.1 QTL region, as it is involved in the cell wall biogenesis pathway and hence may be implicated in modulating the ability of rice roots to penetrate further into the substratum when exposed to drought conditions. RNAs encoding elongation factor 1, a DnaJ-related protein, and a homolog of wheat zinc-finger protein were more prominently induced in the leaves of IR62266 (the lowland rice parent of the mapping materials used) than in those of CT9993 (the upland rice parent) under drought conditions. Homologs of 18S ribosomal RNA, and mRNAs for a multiple-stress induced zinc-finger protein, a protease inhibitor, and a glutathione-S-transferase were expressed at significantly higher levels in CT9993 than in IR62266. Thus several genes involved in the regulation of DNA structure and mRNA translation were found to be drought-regulated, and may be implicated in drought resistance.Communicated by R. Hagemann     

ATP prevents Woronin bodies from sealing septal pores in unwounded cells of the fungus Zymoseptoria tritici

Septa of filamentous ascomycetes are perforated by septal pores that allow communication between individual hyphal compartments. Upon injury, septal pores are plugged rapidly by Woronin bodies (WBs), thereby preventing extensive cytoplasmic bleeding. The mechanism by which WBs translocate into the pore is not known, but it has been suggested that wound‐induced cytoplasmic bleeding “flushes” WBs into the septal opening. Alternatively, contraction of septum‐associated tethering proteins may pull WBs into the septal pore. Here, we investigate WB dynamics in the wheat pathogen Zymoseptoria tritici. Ultrastructural studies showed that 3.4 ± 0.2 WBs reside on each side of a septum and that single WBs of 128.5 ± 3.6 nm in diameter seal the septal pore (41 ± 1.5 nm). Live cell imaging of green fluorescent ZtHex1, a major protein in WBs, and the integral plasma membrane protein ZtSso1 confirms WB translocation into the septal pore. This was associated with the occasional formation of a plasma membrane “balloon,” extruding into the dead cell, suggesting that the plasma membrane rapidly seals the wounded septal pore wound. Minor amounts of fluorescent ZtHex1‐enhanced green fluorescent protein (eGFP) appeared associated with the “ballooning” plasma membrane, indicating that cytoplasmic ZtHex1‐eGFP is recruited to the extending plasma membrane. Surprisingly, in ~15% of all cases, WBs moved from the ruptured cell into the septal pore. This translocation against the cytoplasmic flow suggests that an active mechanism drives WB plugging. Indeed, treatment of unwounded and intact cells with the respiration inhibitor carbonyl cyanide m‐chlorophenyl hydrazone induced WB translocation into the pores. Moreover, carbonyl cyanide m‐chlorophenyl hydrazone treatment recruited cytoplasmic ZtHex1‐eGFP to the lateral plasma membrane of the cells. Thus, keeping the WBs out of the septal pores, in Z. tritici, is an ATP‐dependent process.     

Genomic analysis of differentiation between soil types reveals candidate genes for local adaptation in Arabidopsis lyrata

Serpentine soil, which is naturally high in heavy metal content and has low calcium to magnesium ratios, comprises a difficult environment for most plants. An impressive number of species are endemic to serpentine, and a wide range of non-endemic plant taxa have been shown to be locally adapted to these soils. Locating genomic polymorphisms which are differentiated between serpentine and non-serpentine populations would provide candidate loci for serpentine adaptation. We have used the Arabidopsis thaliana tiling array, which has 2.85 million probes throughout the genome, to measure genetic differentiation between populations of Arabidopsis lyrata growing on granitic soils and those growing on serpentinic soils. The significant overrepresentation of genes involved in ion transport and other functions provides a starting point for investigating the molecular basis of adaptation to soil ion content, water retention, and other ecologically and economically important variables. One gene in particular, calcium-exchanger 7, appears to be an excellent candidate gene for adaptation to low CaratioMg ratio in A. lyrata.     

Quantitative trait locus (QTL) mapping reveals a role for unstudied genes in Aspergillus virulence

Infections caused by the fungus Aspergillus are a major cause of morbidity and mortality in immunocompromised populations. To identify genes required for virulence that could be used as targets for novel treatments, we mapped quantitative trait loci (QTL) affecting virulence in the progeny of a cross between two strains of A. nidulans (FGSC strains A4 and A91). We genotyped 61 progeny at 739 single nucleotide polymorphisms (SNP) spread throughout the genome, and constructed a linkage map that was largely consistent with the genomic sequence, with the exception of one potential inversion of ～527 kb on Chromosome V. The estimated genome size was 3705 cM and the average intermarker spacing was 5.0 cM. The average ratio of physical distance to genetic distance was 8.1 kb/cM, which is similar to previous estimates, and variation in recombination rate was significantly positively correlated with GC content, a pattern seen in other taxa. To map QTL affecting virulence, we measured the ability of each progeny strain to kill model hosts, larvae of the wax moth Galleria mellonella. We detected three QTL affecting in vivo virulence that were distinct from QTL affecting in vitro growth, and mapped the virulence QTL to regions containing 7-24 genes, excluding genes with no sequence variation between the parental strains and genes with only synonymous SNPs. None of the genes in our QTL target regions have been previously associated with virulence in Aspergillus, and almost half of these genes are currently annotated as "hypothetical". This study is the first to map QTL affecting the virulence of a fungal pathogen in an animal host, and our results illustrate the power of this approach to identify a short list of unknown genes for further investigation.     

Genome wide QTL mapping to identify candidate genes for carcass traits in Hanwoo (Korean Cattle)

Meat quality traits are the most economically important traits affecting the beef industry in Korea. We performed a whole genome quantitative trait locus (QTL) mapping study of carcass data in Hanwoo Korean cattle. Two hundred sixty-six Hanwoo steers from 65 sires were genotyped using a 10K Affymetrix SNP chip. The average SNP interval across the bovine genome was 1.5Mb. Associations between each individual SNP and four carcass traits [carcass weight (CWT), eye muscle area (EMA), back fat thickness (BFT), and marbling (MAR)] were assessed using a linear mixed model of each trait. Combined linkage and linkage disequilibrium analysis (LDLA) detected six potential QTL on BTA04, 06, 13, 16, 17, and 23 at the chromosome-wise level (P<0.05). Two MAR QTL were detected at 52.2 cM of BTA06 and 46.04 cM of BTA17. We identified three genes (ARAP2, LOC539460, and LOC511424) in the QTL region of BTA06 and seven genes (RPS14, SCARB1, LOC782103, BRI3BP, AACS, DHX37, and UBC) in the QTL region of BTA17. One significant QTL for CWT was detected at 100 cM on BTA04 and the corresponding QTL region spanned 1.7 cM from 99.7 to 101.4 cM. For EMA QTL, one significant QTL was detected at 3.9 cM of BTA23 and the most likely QTL interval was 1.4 cM, placing 15 candidate genes in the marker bracket. Finally, two QTL for BFT were identified at 68 cM on BTA13 and 24 cM on BTA16. The LPIN3 gene, which is functionally associated with lipodystrophy in humans, is located in the BFT QTL on BTA13. Thus, two potential candidate genes, acetoacetyl-CoA synthetase (AACS) and lipin (LPIN), were detected in QTL regions on BTA17 for MAR and BTA13 for BFT, respectively. In conclusion, LDLA analysis can be used to detect chromosome regions harboring QTL and candidate genes with a low density SNP panel, yielding relatively narrow confidence intervals regarding location.     

Simultaneous mapping of epistatic QTL in chickens reveals clusters of QTL pairs with similar genetic effects on growth

We used simultaneous mapping of interacting quantitative trait locus (QTL) pairs to study various growth traits in a chicken F2 intercross. The method was shown to increase the number of detected QTLs by 30 % compared with a traditional method detecting QTLs by their marginal genetic effects. Epistasis was shown to be an important contributor to the genetic variance of growth, with the largest impact on early growth (before 6 weeks of age). There is also evidence for a discrete set of interacting loci involved in early growth, supporting the previous findings of different genetic regulation of early and late growth in chicken. The genotype-phenotype relationship was evaluated for all interacting QTL pairs and 17 of the 21 evaluated QTL pairs could be assigned to one of four clusters in which the pairs in a cluster have very similar genetic effects on growth. The genetic effects of the pairs indicate commonly occurring dominance-by-dominance, heterosis and multiplicative interactions. The results from this study clearly illustrate the increase in power obtained by using this novel method for simultaneous detection of epistatic QTL, and also how visualization of genotype-phenotype relationships for epistatic QTL pairs provides new insights to biological mechanisms underlying complex traits.     

QTL mapping and the genetic basis of adaptation: recent developments

Quantitative trait loci (QTL) mapping has been used in a number of evolutionary studies to study the genetic basis of adaptation by mapping individual QTL that explain the differences between differentiated populations and also estimating their effects and interaction in the mapping population. This analysis can provide clues about the evolutionary history of populations and causes of the population differentiation. QTL mapping analysis methods and associated computer programs provide us tools for such an inference on the genetic basis and architecture of quantitative trait variation in a mapping population. Current methods have the capability to separate and localize multiple QTL and estimate their effects and interaction on a quantitative trait. More recent methods have been targeted to provide a comprehensive inference on the overall genetic architecture of multiple traits in a number of environments. This development is important for evolutionary studies on the genetic basis of multiple trait variation, genotype by environment interaction, host–parasite interaction, and also microarray gene expression QTL analysis.     

Evolutionary analyses of the avirulence effector AvrStb6 in global populations of Zymoseptoria tritici identify candidate amino acids involved in recognition

We analysed the population genetic diversity of AvrStb6, the first avirulence gene cloned from the wheat pathogen Zymoseptoria tritici, using 142 Z. tritici strains sampled from four wheat fields growing on three continents. Although AvrStb6 was located in a recombination hotspot, it was found in every strain, with 71 polymorphic sites that produced 41 distinct DNA haplotypes encoding 30 AvrStb6 protein isoforms. An AvrStb6 homologue was found in the closest known relative, Z. pseudotritici, but not in three other closely related Zymoseptoria species, indicating that this gene has emerged in Zymoseptoria quite recently. Two AvrStb6 homologues with nucleotide similarities greater than 70% were identified on chromosome 10 in all Z. tritici isolates, suggesting that AvrStb6 belongs to a multigene family of candidate effectors that has expanded recently through gene duplication. The AvrStb6 sequences exhibited strong evidence for non‐neutral evolution, including a large number of non‐synonymous mutations, with significant positive diversifying selection operating on nine of the 82 codons. It appears that balancing selection is operating across the entire gene in natural field populations. There was also evidence for co‐evolving codons within the gene that may reflect compensatory mutations associated with the evasion of recognition by Stb6. Intragenic recombination also appears to have affected the diversity of AvrStb6.