SNP discovery and genotyping using restriction‐site‐associated DNA sequencing in chickens

Single nucleotide polymorphisms (SNPs) are essential to the understanding of population genetic variation and diversity. Here, we performed restriction‐site‐associated DNA sequencing (RAD‐seq) on 72 individuals from 13 Chinese indigenous and three introduced chicken breeds. A total of 620 million reads were obtained using an Illumina Hiseq2000 sequencer. An average of 75 587 SNPs were identified from each individual. Further filtering strictly validated 28 895 SNPs candidates for all populations. When compared with the NCBI dbSNP (chicken_9031), 15 404 SNPs were new discoveries. In this study, RAD‐seq was performed for the first time on chickens, implicating the remarkable effectiveness and potential applications on genetic analysis and breeding technique for whole‐genome selection in chicken and other agricultural animals.     

Multiplex restriction amplicon sequencing: a novel next‐generation sequencing‐based marker platform for high‐throughput genotyping

To enable rapid selection of traits in marker‐assisted breeding, markers must be technically simple, low‐cost, high‐throughput and randomly distributed in a genome. We developed such a technology, designated as Multiplex Restriction Amplicon Sequencing (MRASeq), which reduces genome complexity by polymerase chain reaction (PCR) amplification of amplicons flanked by restriction sites. The first PCR primers contain restriction site sequences at 3’‐ends, preceded by 6‐10 bases of specific or degenerate nucleotide sequences and then by a unique M13‐tail sequence which serves as a binding site for a second PCR that adds sequencing primers and barcodes to allow sample multiplexing for sequencing. The sequences of restriction sites and adjacent nucleotides can be altered to suit different species. Physical mapping of MRASeq SNPs from a biparental population of allohexaploid wheat (Triticum aestivum L.) showed a random distribution of SNPs across the genome. MRASeq generated thousands of SNPs from a wheat biparental population and natural populations of wheat and barley (Hordeum vulgare L.). This novel, next‐generation sequencing‐based genotyping platform can be used for linkage mapping to screen quantitative trait loci (QTL), background selection in breeding and many other genetics and breeding applications of various species.     

Genetic dissection of a major Fusarium head blight QTL in tetraploid wheat

The devastating effect of Fusarium head blight (FHB) caused by Fusarium graminearum has led to significant financial losses across the Upper Midwest of the USA. These losses have spurred the need for research in biological, chemical, and genetic control methods for this disease. To date, most of the research on FHB resistance has concentrated on hexaploid wheat (Triticum aestivum L.) lines originating from China. Other sources of resistance to FHB would be desirable. One other source of resistance for both hexaploid wheat and tetraploid durum wheat (T. turgidum L. var. durum) is the wild tetraploid, T. turgidum L. var. dicoccoides (T. dicoccoides). Previous analysis of the `Langdon'-T. dicoccoides chromosome substitution lines, LDN(Dic), indicated that the chromosome 3A substitution line expresses moderate levels of resistance to FHB. LDN(Dic-3A) recombinant inbred chromosome lines (RICL) were used to generate a linkage map of chromosome 3A with 19 molecular markers spanning a distance of 155.2 cM. The individual RICL and controls were screened for their FHB phenotype in two greenhouse seasons. Analysis of 83 RICL identified a single major quantitative trait locus, Qfhs.ndsu-3AS, that explains 37% of the phenotypic or 55% of the genetic variation for FHB resistance. A microsatellite locus, Xgwm2, is tightly linked to the highest point of the QTL peak. A region of the LDN (Dic-3A) chromosome associated with the QTL for FHB resistance encompasses a 29.3 cM region from Xmwg14 to Xbcd828.     

High‐density SNP genotyping array for hexaploid wheat and its secondary and tertiary gene pool

In wheat, a lack of genetic diversity between breeding lines has been recognized as a significant block to future yield increases. Species belonging to bread wheat's secondary and tertiary gene pools harbour a much greater level of genetic variability, and are an important source of genes to broaden its genetic base. Introgression of novel genes from progenitors and related species has been widely employed to improve the agronomic characteristics of hexaploid wheat, but this approach has been hampered by a lack of markers that can be used to track introduced chromosome segments. Here, we describe the identification of a large number of single nucleotide polymorphisms that can be used to genotype hexaploid wheat and to identify and track introgressions from a variety of sources. We have validated these markers using an ultra‐high‐density Axiom^® genotyping array to characterize a range of diploid, tetraploid and hexaploid wheat accessions and wheat relatives. To facilitate the use of these, both the markers and the associated sequence and genotype information have been made available through an interactive web site.     

High‐density molecular characterization and association mapping in Ethiopian durum wheat landraces reveals high diversity and potential for wheat breeding

Durum wheat (Triticum turgidum subsp. durum) is a key crop worldwide, and yet, its improvement and adaptation to emerging environmental threats is made difficult by the limited amount of allelic variation included in its elite pool. New allelic diversity may provide novel loci to international crop breeding through quantitative trait loci (QTL) mapping in unexplored material. Here, we report the extensive molecular and phenotypic characterization of hundreds of Ethiopian durum wheat landraces and several Ethiopian improved lines. We test 81 587 markers scoring 30 155 single nucleotide polymorphisms and use them to survey the diversity, structure, and genome‐specific variation in the panel. We show the uniqueness of Ethiopian germplasm using a siding collection of Mediterranean durum wheat accessions. We phenotype the Ethiopian panel for ten agronomic traits in two highly diversified Ethiopian environments for two consecutive years and use this information to conduct a genome‐wide association study. We identify several loci underpinning agronomic traits of interest, both confirming loci already reported and describing new promising genomic regions. These loci may be efficiently targeted with molecular markers already available to conduct marker‐assisted selection in Ethiopian and international wheat. We show that Ethiopian durum wheat represents an important and mostly unexplored source of durum wheat diversity. The panel analysed in this study allows the accumulation of QTL mapping experiments, providing the initial step for a quantitative, methodical exploitation of untapped diversity in producing a better wheat.     

Development and validation of the Axiom®Apple480K SNP genotyping array

Cultivated apple (Malus × domestica Borkh.) is one of the most important fruit crops in temperate regions, and has great economic and cultural value. The apple genome is highly heterozygous and has undergone a recent duplication which, combined with a rapid linkage disequilibrium decay, makes it difficult to perform genome‐wide association (GWA) studies. Single nucleotide polymorphism arrays offer highly multiplexed assays at a relatively low cost per data point and can be a valid tool for the identification of the markers associated with traits of interest. Here, we describe the development and validation of a 487K SNP Affymetrix Axiom^® genotyping array for apple and discuss its potential applications. The array has been built from the high‐depth resequencing of 63 different cultivars covering most of the genetic diversity in cultivated apple. The SNPs were chosen by applying a focal points approach to enrich genic regions, but also to reach a uniform coverage of non‐genic regions. A total of 1324 apple accessions, including the 92 progenies of two mapping populations, have been genotyped with the Axiom^®Apple480K to assess the effectiveness of the array. A large majority of SNPs (359 994 or 74%) fell in the stringent class of poly high resolution polymorphisms. We also devised a filtering procedure to identify a subset of 275K very robust markers that can be safely used for germplasm surveys in apple. The Axiom^®Apple480K has now been commercially released both for public and proprietary use and will likely be a reference tool for GWA studies in apple.     

The xylanase inhibitor TAXI‐III counteracts the necrotic activity of a Fusarium graminearum xylanase in vitro and in durum wheat transgenic plants

The xylanase inhibitor TAXI‐III has been proven to delay Fusarium head blight (FHB) symptoms caused by Fusarium graminearum in transgenic durum wheat plants. To elucidate the molecular mechanism underlying the capacity of the TAXI‐III transgenic plants to limit FHB symptoms, we treated wheat tissues with the xylanase FGSG_03624, hitherto shown to induce cell death and hydrogen peroxide accumulation. Experiments performed on lemmas of flowering wheat spikes and wheat cell suspension cultures demonstrated that pre‐incubation of xylanase FGSG_03624 with TAXI‐III significantly decreased cell death. Most interestingly, a reduced cell death relative to control non‐transgenic plants was also obtained by treating, with the same xylanase, lemmas of TAXI‐III transgenic plants. Molecular modelling studies predicted an interaction between the TAXI‐III residue H395 and residues E122 and E214 belonging to the active site of xylanase FGSG_03624. These results provide, for the first time, clear indications in vitro and in planta that a xylanase inhibitor can prevent the necrotic activity of a xylanase, and suggest that the reduced FHB symptoms on transgenic TAXI‐III plants may be a result not only of the direct inhibition of xylanase activity secreted by the pathogen, but also of the capacity of TAXI‐III to avoid host cell death.     

SimRAD: an R package for simulation‐based prediction of the number of loci expected in RADseq and similar genotyping by sequencing approaches

Application of high‐throughput sequencing platforms in the field of ecology and evolutionary biology is developing quickly with the introduction of efficient methods to reduce genome complexity. Numerous approaches for genome complexity reduction have been developed using different combinations of restriction enzymes, library construction strategies and fragment size selection. As a result, the choice of which techniques to use may become cumbersome, because it is difficult to anticipate the number of loci resulting from each method. We developed SimRAD, an R package that performs in silico restriction enzyme digests and fragment size selection as implemented in most restriction site associated DNA polymorphism and genotyping by sequencing methods. In silico digestion is performed on a reference genome or on a randomly generated DNA sequence when no reference genome sequence is available. SimRAD accurately predicts the number of loci under alternative protocols when a reference genome sequence is available for the targeted species (or a close relative) but may be unreliable when no reference genome is available. SimRAD is also useful for fine‐tuning a given protocol to adjust the number of targeted loci. Here, we outline the functionality of SimRAD and provide an illustrative example of the use of the package (available on the CRAN at http://cran.r-project.org/web/packages/SimRAD ).     

Genotype‐specific SNP map based on whole chromosome 3B sequence information from wheat cultivars Arina and Forno

Agronomically important traits are frequently controlled by rare, genotype‐specific alleles. Such genes can only be mapped in a population derived from the donor genotype. This requires the development of a specific genetic map, which is difficult in wheat because of the low level of polymorphism among elite cultivars. The absence of sufficient polymorphism, the complexity of the hexaploid wheat genome as well as the lack of complete sequence information make the construction of genetic maps with a high density of reproducible and polymorphic markers challenging. We developed a genotype‐specific genetic map of chromosome 3B from winter wheat cultivars Arina and Forno. Chromosome 3B was isolated from the two cultivars and then sequenced to 10‐fold coverage. This resulted in a single‐nucleotide polymorphisms (SNP) database of the complete chromosome. Based on proposed synteny with the Brachypodium model genome and gene annotation, sequences close to coding regions were used for the development of 70 SNP‐based markers. They were mapped on a Arina × Forno Recombinant Inbred Lines population and found to be spread over the complete chromosome 3B. While overall synteny was well maintained, numerous exceptions and inversions of syntenic gene order were identified. Additionally, we found that the majority of recombination events occurred in distal parts of chromosome 3B, particularly in hot‐spot regions. Compared with the earlier map based on SSR and RFLP markers, the number of markers increased fourfold. The approach presented here allows fast development of genotype‐specific polymorphic markers that can be used for mapping and marker‐assisted selection.     

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.     

Next‐generation sequencing,FISH mapping and synteny‐based modeling reveal mechanisms of decreasing dysploidy in Cucumis

In the large Cucurbitaceae genus Cucumis, cucumber (C. sativus) is the only species with 2n = 2x = 14 chromosomes. The majority of the remaining species, including melon (C. melo) and the sister species of cucumber, C. hystrix, have 2n = 2x = 24 chromosomes, implying a reduction from n = 12 to n = 7. To understand the underlying mechanisms, we investigated chromosome synteny among cucumber, C. hystrix and melon using integrated and complementary approaches. We identified 14 inversions and a C. hystrix lineage‐specific reciprocal inversion between C. hystrix and melon. The results reveal the location and orientation of 53 C. hystrix syntenic blocks on the seven cucumber chromosomes, and allow us to infer at least 59 chromosome rearrangement events that led to the seven cucumber chromosomes, including five fusions, four translocations, and 50 inversions. The 12 inferred chromosomes (AK1–AK12) of an ancestor similar to melon and C. hystrix had strikingly different evolutionary fates, with cucumber chromosome C1 apparently resulting from insertion of chromosome AK12 into the centromeric region of translocated AK2/AK8, cucumber chromosome C3 originating from a Robertsonian‐like translocation between AK4 and AK6, and cucumber chromosome C5 originating from fusion of AK9 and AK10. Chromosomes C2, C4 and C6 were the result of complex reshuffling of syntenic blocks from three (AK3, AK5 and AK11), three (AK5, AK7 and AK8) and five (AK2, AK3, AK5, AK8 and AK11) ancestral chromosomes, respectively, through 33 fusion, translocation and inversion events. Previous results (Huang, S., Li, R., Zhang, Z. et al., 2009 , Nat. Genet. 41, 1275–1281; Li, D., Cuevas, H.E., Yang, L., Li, Y., Garcia‐Mas, J., Zalapa, J., Staub, J.E., Luan, F., Reddy, U., He, X., Gong, Z., Weng, Y. 2011a, BMC Genomics, 12, 396) showing that cucumber C7 stayed largely intact during the entire evolution of Cucumis are supported. Results from this study allow a fine‐scale understanding of the mechanisms of dysploid chromosome reduction that has not been achieved previously.     

High‐throughput SNP discovery in the rabbit (Oryctolagus cuniculus) genome by next‐generation semiconductor‐based sequencing

The European rabbit (Oryctolagus cuniculus) is a domesticated species with one of the broadest ranges of economic and scientific applications and fields of investigation. Rabbit genome information and assembly are available (oryCun2.0), but so far few studies have investigated its variability, and massive discovery of polymorphisms has not been published yet for this species. Here, we sequenced two reduced representation libraries (RRLs) to identify single nucleotide polymorphisms (SNPs) in the rabbit genome. Genomic DNA of 10 rabbits belonging to different breeds was pooled and digested with two restriction enzymes (HaeIII and RsaI) to create two RRLs which were sequenced using the Ion Torrent Personal Genome Machine. The two RRLs produced 2 917 879 and 4 046 871 reads, for a total of 280.51 Mb (248.49 Mb with quality >20) and 417.28 Mb (360.89 Mb with quality >20) respectively of sequenced DNA. About 90% and 91% respectively of the obtained reads were mapped on the rabbit genome, covering a total of 15.82% of the oryCun2.0 genome version. The mapping and ad hoc filtering procedures allowed to reliably call 62 491 SNPs. SNPs in a few genomic regions were validated by Sanger sequencing. The Variant Effect Predictor Web tool was used to map SNPs on the current version of the rabbit genome. The obtained results will be useful for many applied and basic research programs for this species and will contribute to the development of cost‐effective solutions for high‐throughput SNP genotyping in the rabbit.     

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.     

Towards a whole‐genome sequence for rye (Secale cereale L.)

We report on a whole‐genome draft sequence of rye (Secale cereale L.). Rye is a diploid Triticeae species closely related to wheat and barley, and an important crop for food and feed in Central and Eastern Europe. Through whole‐genome shotgun sequencing of the 7.9‐Gbp genome of the winter rye inbred line Lo7 we obtained a de novo assembly represented by 1.29 million scaffolds covering a total length of 2.8 Gbp. Our reference sequence represents nearly the entire low‐copy portion of the rye genome. This genome assembly was used to predict 27 784 rye gene models based on homology to sequenced grass genomes. Through resequencing of 10 rye inbred lines and one accession of the wild relative S. vavilovii, we discovered more than 90 million single nucleotide variants and short insertions/deletions in the rye genome. From these variants, we developed the high‐density Rye600k genotyping array with 600 843 markers, which enabled anchoring the sequence contigs along a high‐density genetic map and establishing a synteny‐based virtual gene order. Genotyping data were used to characterize the diversity of rye breeding pools and genetic resources, and to obtain a genome‐wide map of selection signals differentiating the divergent gene pools. This rye whole‐genome sequence closes a gap in Triticeae genome research, and will be highly valuable for comparative genomics, functional studies and genome‐based breeding in rye.     

A genome‐wide association study of 23 agronomic traits in Chinese wheat landraces

Uncovering the genetic basis of agronomic traits in wheat landraces is important for ensuring global food security via the development of improved varieties. Here, 723 wheat landraces from 10 Chinese agro‐ecological zones were evaluated for 23 agronomic traits in six environments. All accessions could be clustered into five subgroups based on phenotypic data via discriminant function analysis, which was highly consistent with genotypic classification. A genome‐wide association study was conducted for these traits using 52 303 DArT‐seq markers to identify marker‐trait associations and candidate genes. Using both the general linear model and the mixed linear model, 149 significant markers were identified for 21 agronomic traits based on best linear unbiased prediction values. Considering the linkage disequilibrium decay distance in this study, significant markers within 10 cM were combined as a quantitative trait locus (QTL), with a total of 29 QTL identified for 15 traits. Of these, five QTL for heading date, flag leaf width, peduncle length, and thousand kernel weight had been reported previously. Twenty‐five candidate genes associated with significant markers were identified. These included the known vernalization genes VRN‐B1 and vrn‐B3 and the photoperiod response genes Ppd and PRR. Overall, this study should be helpful in elucidating the underlying genetic mechanisms of complex agronomic traits and performing marker‐assisted selection in wheat.     

Identification of a large SNP dataset in Larimichthys crocea using specific‐locus amplified fragment sequencing

The large yellow croaker, Larimichthys crocea, is a commercially important drum fish (Family: Sciaenidae) native to the East and South China Sea. Habitat deterioration and overfishing have led to significant population decline and the collapse of its fishery over the past decades. Today, the market supply of L. crocea depends solely on stocks produced in hatcheries and farms. Common issues that occur in the culture of L. crocea include germplasm degradation, precocious puberty, elevated disease susceptibility and growth retardation. In this study, we employed SLAF‐seq (specific‐locus amplified fragment sequencing) technology to identify single nucleotide polymorphism (SNP) loci across the L. crocea genome. Sixty samples were selected for SLAF analysis out of 1000 progeny in the same cohort of a cultured stock. Our analysis obtained a total of 151 253 SLAFs, of which 65.88% (99 652) were identified to be polymorphic, scoring a total of 710 567 putative SNPs. Further filtration resulted in a final panel of 1782 SNP loci. The data derived from this work could be beneficial for understanding the genetics of complex phenotypic traits as well as for developing marker‐selection‐assisted breeding programs in L. crocea.     

Genome‐wide association study of six quality traits reveals the association of the TaRPP13L1 gene with flour colour in Chinese bread wheat

Flour colour, kernel hardness, grain protein content and wet gluten content are important quality properties that determine end use in bread wheat. Here, a wheat 90K genotyping assay was used for a genome‐wide association study (GWAS) of the six quality‐related traits in Chinese wheat cultivars in eight environments over four years. A total of 846 significant single nucleotide polymorphisms (SNPs) were identified, explaining approximately 30% of the phenotypic variation on average, and 103 multienvironment‐significant SNPs were detected in more than four environments. Quantitative trait loci (QTL) mapping in the biparent population confirmed some important SNP loci. Moreover, it was determined that some important genes were associated with the six quality traits, including some known functional genes and annotated unknown functional genes. Of the annotated unknown functional genes, it was verified that TaRPP13L1 was associated with flour colour. Wheat cultivars or lines with TaRPP13L1‐B1a showed extremely significantly higher flour redness and lower yellowness than those with TaRPP13L1‐B1b in the Chinese wheat natural population and the doubled haploid (DH) population. Two tetraploid wheat lines with premature stop codons of the TaRPP13L1 gene mutagenized by ethyl methanesulfonate (EMS) showed extremely significantly higher flour redness and lower yellowness than wild type. Our data suggest that the TaRPP13L1 gene plays an important role in modulating wheat flour colour. This study provides useful information for further dissection of the genetic basis of flour colour and also provides valuable genes or genetic loci for marker‐assisted selection to improve the process of breeding quality wheat in China.