Genotyping with RAPD and microsatellite markers resolves pathotype diversity in the ascochyta blight pathogen of chickpea

 The poor definition of variation in the ascochyta blight fungus (Ascochyta rabiei) has historically hindered breeding for resistance to the chickpea (Cicer arietinum L.) blight disease in West Asia and North Africa. We have employed 14 RAPD markers and an oligonucleotide probe complementary to the microsatellite sequence (GATA)₄ to construct a genotype-specific DNA fragment profile from periodically sampled Syrian field isolates of this fungus. By using conventional pathogenicity tests and genome analysis with RAPD and microsatellite markers, we demonstrated that the DNA markers distinguish variability within and among the major pathotypes of A. rabiei and resolved each pathotypes into several genotypes. The genetic diversity estimate based on DNA marker analysis within pathotypes was highest for the least-aggressive pathotype (pathotype I), followed by the aggressive (pathotype II) and the most-aggressive pathotype (pathotype III). The pair-wise genetic distance estimated for all the isolates varied from 0.00 to 0.39, indicating a range from a clonal to a diverse relationship. On the basis of genome analysis, and information on the spatial and temporal distribution of the pathogen, a general picture of A. rabiei evolution in Syria is proposed. Received: 10 January 1998 / Accepted: 23 January 1998     

Limber pine (Pinus flexilis James) genetic map constructed by exome‐seq provides insight into the evolution of disease resistance and a genomic resource for genomics‐based breeding

Limber pine ( Pinus flexilis ) is a keystone species of high‐elevation forest ecosystems of western North America, but some parts of the geographic range have high infection and mortality from the non‐native white pine blister rust caused by Cronartium ribicola . Genetic maps can provide essential knowledge for understanding genetic disease resistance as well as local adaptation to changing climates. Exome‐seq was performed to construct high‐density genetic maps in two seed families. Composite maps positioned 9612 unigenes across 12 linkage groups ( LG s). Syntenic analysis of genome structure revealed that the majority of orthologs were positional orthologous genes ( POG s) with localization on homologous LG s among conifer species. Gene ontology ( GO) enrichment analysis showed relatively fewer constraints for POG s with putative roles in adaptation to environments and relatively more conservation for POG s with roles in basic cell function and maintenance. The mapped genes included 639 nucleotide‐binding site leucine‐rich repeat genes ( NBS ‐ LRR s) , 290 receptor‐like protein kinase genes ( RLK s), and 1014 genes with potential roles in the defense response and induced systemic resistance to attack by pathogens. Orthologous loci for resistance to rust pathogens were identified and were co‐positioned with multiple members of the R gene family, revealing the evolutionary pressure acting upon them. This high‐density genetic map provides a genomic resource and practical tool for breeding and genetic conservation programs, with applications in genome‐wide association studies ( GWASs ), the characterization of functional genes underlying complex traits, and the sequencing and assembly of the full‐length genomes of limber pine and related Pinus species.     

QTL analysis for ascochyta blight resistance in an intraspecific population of chickpea (<Emphasis Type="Italic">Cicer arietinum</Emphasis> L.)

In both controlled environment and the field, six QTLs for ascochyta blight resistance were identified in three regions of the genome of an intraspecific population of chickpea using the IDS and AUDPC disease scoring systems. One QTL-region was detected from both environments, whereas the other two regions were detected from each environment. All the QTL-regions were significantly associated with ascochyta blight resistance using either of the disease scoring systems. The QTLs were verified by multiple interval mapping, and a two-QTL genetic model with considerable epistasis was established for both environments. The major QTLs generally showed additive gene action, as well as dominance inter-locus interaction in the multiple genetic model. All the QTLs were mapped near a RGA marker. The major QTLs were located on LG III, which was mapped with five different types of RGA markers. A CLRR-RGA marker and a STMS marker flanked QTL 6 for controlled environment resistance at 0.06 and 0.04 cM, respectively. Other STMS markers flanked QTL 1 for field resistance at a 5.6 cM interval. After validation, these flanking markers may be used in marker-assisted selection to breed for elite chickpea cultivars with durable resistance to ascochyta blight. The tight linkage of RGA markers to the major QTL on LG III will allow map-based cloning of the underlying resistance genes.Communicated by P. Langridge     

Helping decision making for reliable and cost‐effective 2b‐RAD sequencing and genotyping analyses in non‐model species

High‐throughput sequencing has revolutionized population and conservation genetics. RAD sequencing methods, such as 2b‐RAD, can be used on species lacking a reference genome. However, transferring protocols across taxa can potentially lead to poor results. We tested two different IIB enzymes (AlfI and CspCI) on two species with different genome sizes (the loggerhead turtle Caretta caretta and the sharpsnout seabream Diplodus puntazzo) to build a set of guidelines to improve 2b‐RAD protocols on non‐model organisms while optimising costs. Good results were obtained even with degraded samples, showing the value of 2b‐RAD in studies with poor DNA quality. However, library quality was found to be a critical parameter on the number of reads and loci obtained for genotyping. Resampling analyses with different number of reads per individual showed a trade‐off between number of loci and number of reads per sample. The resulting accumulation curves can be used as a tool to calculate the number of sequences per individual needed to reach a mean depth ≥20 reads to acquire good genotyping results. Finally, we demonstrated that selective‐base ligation does not affect genomic differentiation between individuals, indicating that this technique can be used in species with large genome sizes to adjust the number of loci to the study scope, to reduce sequencing costs and to maintain suitable sequencing depth for a reliable genotyping without compromising the results. Here, we provide a set of guidelines to improve 2b‐RAD protocols on non‐model organisms with different genome sizes, helping decision‐making for a reliable and cost‐effective genotyping.     

Haplotype‐based genotyping‐by‐sequencing in oat genome research

In a de novo genotyping‐by‐sequencing (GBS) analysis of short, 64‐base tag‐level haplotypes in 4657 accessions of cultivated oat, we discovered 164741 tag‐level (TL) genetic variants containing 241224 SNPs. From this, the marker density of an oat consensus map was increased by the addition of more than 70000 loci. The mapped TL genotypes of a 635‐line diversity panel were used to infer chromosome‐level (CL) haplotype maps. These maps revealed differences in the number and size of haplotype blocks, as well as differences in haplotype diversity between chromosomes and subsets of the diversity panel. We then explored potential benefits of SNP vs. TL vs. CL GBS variants for mapping, high‐resolution genome analysis and genomic selection in oats. A combined genome‐wide association study (GWAS) of heading date from multiple locations using both TL haplotypes and individual SNP markers identified 184 significant associations. A comparative GWAS using TL haplotypes, CL haplotype blocks and their combinations demonstrated the superiority of using TL haplotype markers. Using a principal component‐based genome‐wide scan, genomic regions containing signatures of selection were identified. These regions may contain genes that are responsible for the local adaptation of oats to Northern American conditions. Genomic selection for heading date using TL haplotypes or SNP markers gave comparable and promising prediction accuracies of up to r = 0.74. Genomic selection carried out in an independent calibration and test population for heading date gave promising prediction accuracies that ranged between r = 0.42 and 0.67. In conclusion, TL haplotype GBS‐derived markers facilitate genome analysis and genomic selection in oat.     

SNP and INDEL detection in a QTL region on chicken chromosome 2 associated with muscle deposition

Genetic improvement is important for the poultry industry, contributing to increased efficiency of meat production and quality. Because breast muscle is the most valuable part of the chicken carcass, knowledge of polymorphisms influencing this trait can help breeding programs. Therefore, the complete genome of 18 chickens from two different experimental lines (broiler and layer) from EMBRAPA was sequenced, and SNPs and INDELs were detected in a QTL region for breast muscle deposition on chicken chromosome 2 between microsatellite markers MCW0185 and MCW0264 (105 849–112 649 kb). Initially, 94 674 unique SNPs and 10 448 unique INDELs were identified in the target region. After quality filtration, 77% of the SNPs (85 765) and 60% of the INDELs (7828) were retained. The studied region contains 66 genes, and functional annotation of the filtered variants identified 517 SNPs and three INDELs in exonic regions. Of these, 357 SNPs were classified as synonymous, 153 as non‐synonymous, three as stopgain, four INDELs as frameshift and three INDELs as non‐frameshift. These exonic mutations were identified in 37 of the 66 genes from the target region, three of which are related to muscle development (DTNA, RB1CC1 and MOS). Fifteen non‐tolerated SNPs were detected in several genes (MEP1B, PRKDC, NSMAF, TRAPPC8, SDR16C5, CHD7, ST18 and RB1CC1). These loss‐of‐function and exonic variants present in genes related to muscle development can be considered candidate variants for further studies in chickens. Further association studies should be performed with these candidate mutations as should validation in commercial populations to allow a better explanation of QTL effects.     

Functional screening of willow alleles in Arabidopsis combined with QTL mapping in willow (Salix) identifies SxMAX4 as a coppicing response gene

Willows (Salix spp.) are important biomass crops due to their ability to grow rapidly with low fertilizer inputs and ease of cultivation in short‐rotation coppice cycles. They are relatively undomesticated and highly diverse, but functional testing to identify useful allelic variation is time‐consuming in trees and transformation is not yet possible in willow. Arabidopsis is heralded as a model plant from which knowledge can be transferred to advance the improvement of less tractable species. Here, knowledge and methodologies from Arabidopsis were successfully used to identify a gene influencing stem number in coppiced willows, a complex trait of key biological and industrial relevance. The strigolactone‐related More AXillary growth (MAX) genes were considered candidates due to their role in shoot branching. We previously demonstrated that willow and Arabidopsis show similar response to strigolactone and that transformation rescue of Arabidopsis max mutants with willow genes could be used to detect allelic differences. Here, this approach was used to screen 45 SxMAX1, SxMAX2, SxMAX3 and SxMAX4 alleles cloned from 15 parents of 11 mapping populations varying in shoot‐branching traits. Single‐nucleotide polymorphism (SNP) frequencies were locus dependent, ranging from 29.2 to 74.3 polymorphic sites per kb. SxMAX alleles were 98%–99% conserved at the amino acid level, but different protein products varying in their ability to rescue Arabidopsis max mutants were identified. One poor rescuing allele, SxMAX4D, segregated in a willow mapping population where its presence was associated with increased shoot resprouting after coppicing and colocated with a QTL for this trait.     

Interaction of acid exudates in chickpea with biological activity of Bacillus thuringiensis towards Helicoverpa armigera

The gram pod borer, Helicoverpa armigera, is one of the most important constraints to chickpea production. High acidity of chickpea exudates is associated with resistance to pod borer, H. armigera; however, acidic exudates in chickpea might influence the biological activity of the bacterium, Bacillus thuringiensis (Bt), applied as a foliar spray or deployed in transgenic plants for controlling H. armigera. Therefore, studies were undertaken to evaluate the biological activity of Bt towards H. armigera on chickpea genotypes with different amounts of organic acids. Significantly lower leaf feeding, larval survival and larval weights were observed on ICC 506EB, followed by C 235, and ICCV 10 across Bt concentrations. Leaf feeding by the larvae and larval survival and weights decreased with an increase in Bt concentration. However, rate of decrease in leaf feeding and larval survival and weights with an increase in Bt concentration was greater on L 550 and ICCV 10 than on the resistant check, ICC 506EB, suggesting that factors in the resistant genotypes, particularly the acid exudates, resulted in lower levels of biological activity of Bt possibly because of antifeedant effects of the acid exudates. Antifeedant effects of acid exudates reduced food consumption and hence might reduce the efficacy of Bt sprays on insect‐resistant chickpea genotypes or Bt‐transgenic chickpeas, although the combined effect of plant resistance based on organic acids, and Bt had a greater effect on survival and development of H. armigera than Bt alone.     

Detection of selection signatures of population‐specific genomic regions selected during domestication process in Jinhua pigs

Chinese pigs have been undergoing both natural and artificial selection for thousands of years. Jinhua pigs are of great importance, as they can be a valuable model for exploring the genetic mechanisms linked to meat quality and other traits such as disease resistance, reproduction and production. The purpose of this study was to identify distinctive footprints of selection between Jinhua pigs and other breeds utilizing genome‐wide SNP data. Genotyping by genome reducing and sequencing was implemented in order to perform cross‐population extended haplotype homozygosity to reveal strong signatures of selection for those economically important traits. This work was performed at a 2% genome level, which comprised 152 006 SNPs genotyped in a total of 517 individuals. Population‐specific footprints of selective sweeps were searched for in the genome of Jinhua pigs using six native breeds and three European breeds as reference groups. Several candidate genes associated with meat quality, health and reproduction, such as GH1, CRHR2, TRAF4 and CCK, were found to be overlapping with the significantly positive outliers. Additionally, the results revealed that some genomic regions associated with meat quality, immune response and reproduction in Jinhua pigs have evolved directionally under domestication and subsequent selections. The identified genes and biological pathways in Jinhua pigs showed different selection patterns in comparison with the Chinese and European breeds.     

Population and phylogenomic decomposition via genotyping‐by‐sequencing in Australian Pelargonium

Species delimitation has seen a paradigm shift as increasing accessibility of genomic‐scale data enables separation of lineages with convergent morphological traits and the merging of recently diverged ecotypes that have distinguishing characteristics. We inferred the process of lineage formation among Australian species in the widespread and highly variable genus Pelargonium by combining phylogenomic and population genomic analyses along with breeding system studies and character analysis. Phylogenomic analysis and population genetic clustering supported seven of the eight currently described species but provided little evidence for differences in genetic structure within the most widely distributed group that containing P. australe. In contrast, morphometric analysis detected three deep lineages within Australian Pelargonium; with P. australe consisting of five previously unrecognized entities occupying separate geographic ranges. The genomic approach enabled elucidation of parallel evolution in some traits formerly used to delineate species, as well as identification of ecotypic morphological differentiation within recognized species. Highly variable morphology and trait convergence each contribute to the discordance between phylogenomic relationships and morphological taxonomy. Data suggest that genetic divergence among species within the Australian Pelargonium may result from allopatric speciation while morphological differentiation within and among species may be more strongly driven by environmental differences.     

Taking advantage from phenotype variability in a local animal genetic resource: identification of genomic regions associated with the hairless phenotype in Casertana pigs

Casertana is an endangered autochthonous pig breed (raised in south‐central Italy) that is considered to be the descendant of the influential Neapolitan pig population that was used to improve British breeds in the 19th century. Casertana pigs are characterized by a typical, almost complete, hairless phenotype, even though a few Casertana pigs are normal haired. In this work, using Illumina PorcineSNP60 BeadChip data, we carried out a genome‐wide association study and an F_ST analysis with this breed by comparing animals showing the classical hairless phenotype (n = 81) versus pigs classified as haired (n = 15). Combining the results obtained with the two approaches, we identified two significant regions: one on porcine chromosome (SSC) 7 and one on SSC15. The SSC7 region contains the forkhead box N3 (FOXN3) gene, the most plausible candidate gene of this region, considering that mutations in another gene of the same family (forkhead box N1; Foxn1 or FOXN1) are responsible for the nude locus in rodents and alopecia in humans. Another potential candidate gene, rho guanine nucleotide exchange factor 10 (ARHGEF10), is located in the SSC15 region. FOXN3 and ARHGEF10 have been detected as differentially expressed in androgenetic and senescent alopecia respectively. This study on an autochthonous pig breed contributes to shed some light on novel genes potentially involved in hair development and growth and demonstrates that local animal breeds can be valuable genetic resources for disclosing genetic factors affecting unique traits, taking advantage of phenotype variability segregating in small populations.     

Mapping‐by‐sequencing in complex polyploid genomes using genic sequence capture: a case study to map yellow rust resistance in hexaploid wheat

Previously we extended the utility of mapping‐by‐sequencing by combining it with sequence capture and mapping sequence data to pseudo‐chromosomes that were organized using wheat–Brachypodium synteny. This, with a bespoke haplotyping algorithm, enabled us to map the flowering time locus in the diploid wheat Triticum monococcum L. identifying a set of deleted genes (Gardiner et al., 2014). Here, we develop this combination of gene enrichment and sliding window mapping‐by‐synteny analysis to map the Yr6 locus for yellow stripe rust resistance in hexaploid wheat. A 110 MB NimbleGen capture probe set was used to enrich and sequence a doubled haploid mapping population of hexaploid wheat derived from an Avalon and Cadenza cross. The Yr6 locus was identified by mapping to the POPSEQ chromosomal pseudomolecules using a bespoke pipeline and algorithm (Chapman et al., 2015). Furthermore the same locus was identified using newly developed pseudo‐chromosome sequences as a mapping reference that are based on the genic sequence used for sequence enrichment. The pseudo‐chromosomes allow us to demonstrate the application of mapping‐by‐sequencing to even poorly defined polyploidy genomes where chromosomes are incomplete and sub‐genome assemblies are collapsed. This analysis uniquely enabled us to: compare wheat genome annotations; identify the Yr6 locus – defining a smaller genic region than was previously possible; associate the interval with one wheat sub‐genome and increase the density of SNP markers associated. Finally, we built the pipeline in iPlant, making it a user‐friendly community resource for phenotype mapping.     

Quantitative trait loci associated with resistance to Fusarium head blight and kernel discoloration in barley

Resistance to Fusarium head blight (FHB), deoxynivalenol (DON) accumulation, and kernel discoloration (KD) in barley are difficult traits to introgress into elite varieties because current screening methods are laborious and disease levels are strongly influenced by environment. To improve breeding strategies directed toward enhancing these traits, we identified genomic regions containing quantitative trait loci (QTLs) associated with resistance to FHB, DON accumulation, and KD in a breeding population of F_4:7 lines using restriction fragment length polymorphic (RFLP) markers. We evaluated 101 F_4:7 lines, derived from a cross between the cultivar Chevron and an elite breeding line, M69, for each of the traits in three or four environments. We used 94 previously mapped RFLP markers to create a linkage map. Using composite interval mapping, we identified 10, 11, and 4 QTLs associated with resistance to FHB, DON accumulation, and KD, respectively. Markers flanking these QTLs should be useful for introgressing resistance to FHB, DON accumulation, and KD into elite barley cultivars. Received: 8 November 1998 / Accepted: 8 January 1999     

Understanding sheath blight resistance in rice: the road behind and the road ahead

Rice sheath blight disease, caused by the basidiomycetous necrotroph Rhizoctonia solani, became one of the major threats to the rice cultivation worldwide, especially after the adoption of high‐yielding varieties. The pathogen is challenging to manage because of its extensively broad host range and high genetic variability and also due to the inability to find any satisfactory level of natural resistance from the available rice germplasm. It is high time to find remedies to combat the pathogen for reducing rice yield losses and subsequently to minimize the threat to global food security. The development of genetic resistance is one of the alternative means to avoid the use of hazardous chemical fungicides. This review mainly focuses on the effort of better understanding the host–pathogen relationship, finding the gene loci/markers imparting resistance response and modifying the host genome through transgenic development. The latest development and trend in the R. solani–rice pathosystem research with gap analysis are provided.