Exome Capture Data Research Articles

Analysis of a global wheat panel reveals a highly diverse introgression landscape and provides evidence for inter-homoeologue chromosomal recombination

Key messageThis study highlights the agronomic potential of rare introgressions, as demonstrated by a major QTL for powdery mildew resistance on chromosome 7D. It further shows evidence for inter-homoeologue recombination in wheat.Agriculturally important genes are often introgressed into crops from closely related donor species or landraces. The gene pool of hexaploid bread wheat (Triticum aestivum) is known to contain numerous such “alien” introgressions. Recently established high-quality reference genome sequences allow prediction of the size, frequency and identity of introgressed chromosome regions. Here, we characterise chromosomal introgressions in bread wheat using exome capture data from the WHEALBI collection. We identified 24,981 putative introgression segments of at least 2 Mb across 434 wheat accessions. Detailed study of the most frequent introgressions identified T. timopheevii or its close relatives as a frequent donor species. Importantly, 118 introgressions of at least 10 Mb were exclusive to single wheat accessions, revealing that large populations need to be studied to assess the total diversity of the wheat pangenome. In one case, a 14 Mb introgression in chromosome 7D, exclusive to cultivar Pamukale, was shown by QTL mapping to harbour a recessive powdery mildew resistance gene. We identified multiple events where distal chromosomal segments of one subgenome were duplicated in the genome and replaced the homoeologous segment in another subgenome. We propose that these examples are the results of inter-homoeologue recombination. Our study produced an extensive catalogue of the wheat introgression landscape, providing a resource for wheat breeding. Of note, the finding that the wheat gene pool contains numerous rare, but potentially important introgressions and chromosomal rearrangements has implications for future breeding.

Fine mapping of stem rust resistance derived from soft red winter wheat cultivar AGS2000 to an NLR gene cluster on chromosome 6D

The Puccinia graminis f. sp. tritici (Pgt) Ug99-emerging virulent races present a major challenge to global wheat production. To meet present and future needs, new sources of resistance must be found. Identification of markers that allow tracking of resistance genes is needed for deployment strategies to combat highly virulent pathogen races. Field evaluation of a DH population located a QTL for stem rust (Sr) resistance, QSr.nc-6D from the breeding line MD01W28-08-11 to the distal region of chromosome arm 6DS where Sr resistance genes Sr42, SrCad, and SrTmp have been identified. A locus for seedling resistance to Pgt race TTKSK was identified in a DH population and an RIL population derived from the cross AGS2000 × LA95135. The resistant cultivar AGS2000 is in the pedigree of MD01W28-08-11 and our results suggest that it is the source of Sr resistance in this breeding line. We exploited published markers and exome capture data to enrich marker density in a 10 Mb region flanking QSr.nc-6D. Our fine mapping in heterozygous inbred families identified three markers co-segregating with resistance and delimited QSr.nc-6D to a 1.3 Mb region. We further exploited information from other genome assemblies and identified collinear regions of 6DS harboring clusters of NLR genes. Evaluation of KASP assays corresponding to our co-segregating SNP suggests that they can be used to track this Sr resistance in breeding programs. However, our results also underscore the challenges posed in identifying genes underlying resistance in such complex regions in the absence of genome sequence from the resistant genotypes.

Between but Not Within-Species Variation in the Distribution of Fitness Effects.

New mutations provide the raw material for evolution and adaptation. The distribution of fitness effects (DFE) describes the spectrum of effects of new mutations that can occur along a genome, and is, therefore, of vital interest in evolutionary biology. Recent work has uncovered striking similarities in the DFE between closely related species, prompting us to ask whether there is variation in the DFE among populations of the same species, or among species with different degrees of divergence, that is whether there is variation in the DFE at different levels of evolution. Using exome capture data from six tree species sampled across Europe we characterized the DFE for multiple species, and for each species, multiple populations, and investigated the factors potentially influencing the DFE, such as demography, population divergence, and genetic background. We find statistical support for the presence of variation in the DFE at the species level, even among relatively closely related species. However, we find very little difference at the population level, suggesting that differences in the DFE are primarily driven by deep features of species biology, and those evolutionarily recent events, such as demographic changes and local adaptation, have little impact.

Clinal variation in PHY (PAS domain) and CRY (CCT domain)-Signs of local adaptation to light quality in Norway spruce.

Detection of the genomic basis of local adaptation to environmental conditions is challenging in forest trees. Phytochromes (PHY) and cryptochromes (CRY) perceive the red (R)/far-red (FR)and blue lightrespectively, thus playing a fundamental role in regulating plant growth and development. PHYO and PHYP from conifers are the equivalents of PHYA/PHYCand PHYB in angiosperms, respectively. Norway spruce shows an adaptive latitudinal cline for shade (low R:FR or FR-enriched light) tolerance and requirement of FR light for its growth. We analyzed the exome capture data that included a uniquely large data set of 1654 Norway spruce trees sampled across many latitudes in Sweden to capture the natural clines for photoperiod and FR light exposure during the growth season. Statistically significant clinal variation was detected in allele and genotype frequencies of missense mutations in coding regions belonging to well-defined functional domains of PHYO (PAS-B),PHYP2 (PASfold-2), CRY1 (CCT1)and CRY2 (CCT2)that strongly correlates with the latitudinal gradient in response to variable light quality in Norway spruce. The missense SNP in PHYO resulting in Asn835Ser, displayed the steepest cline among all other polymorphisms. We propose that these variations in the photoreceptors represent signs of local adaptation to light quality.

QTL for induced resistance against leaf rust in barley.

Leaf rust caused by Puccinia hordei is one of the major diseases of barley (Hordeum vulgare L.) leading to yield losses up to 60%. Even though, resistance genes Rph1 to Rph28 are known, most of these are already overcome. In this context, priming may promote enhanced resistance to P. hordei. Several bacterial communities such as the soil bacterium Ensifer (syn. Sinorhizobium) meliloti are reported to induce resistance by priming. During quorum sensing in populations of gram negative bacteria, they produce N-acyl homoserine-lactones (AHL), which induce resistance in plants in a species- and genotype-specific manner. Therefore, the present study aims to detect genotypic differences in the response of barley to AHL, followed by the identification of genomic regions involved in priming efficiency of barley. A diverse set of 198 spring barley accessions was treated with a repaired E. meliloti natural mutant strain expR+ch producing a substantial amount of AHL and a transformed E. meliloti strain carrying the lactonase gene attM from Agrobacterium tumefaciens. For P. hordei resistance the diseased leaf area and the infection type were scored 12 dpi (days post-inoculation), and the corresponding relative infection and priming efficiency were calculated. Results revealed significant effects (p<0.001) of the bacterial treatment indicating a positive effect of priming on resistance to P. hordei. In a genome-wide association study (GWAS), based on the observed phenotypic differences and 493,846 filtered SNPs derived from the Illumina 9k iSelect chip, genotyping by sequencing (GBS), and exome capture data, 11 quantitative trait loci (QTL) were identified with a hot spot on the short arm of the barley chromosome 6H, associated to improved resistance to P. hordei after priming with E. meliloti expR+ch. Genes in these QTL regions represent promising candidates for future research on the mechanisms of plant-microbe interactions.

The evolutionary patterns of barley pericentromeric chromosome regions, as shaped by linkage disequilibrium and domestication.

SUMMARYThe distribution of recombination events along large cereal chromosomes is uneven and is generally restricted to gene‐rich telomeric ends. To understand how the lack of recombination affects diversity in the large pericentromeric regions, we analysed deep exome capture data from a final panel of 815 Hordeum vulgare (barley) cultivars, landraces and wild barleys, sampled from across their eco‐geographical ranges. We defined and compared variant data across the pericentromeric and non‐pericentromeric regions, observing a clear partitioning of diversity both within and between chromosomes and germplasm groups. Dramatically reduced diversity was found in the pericentromeres of both cultivars and landraces when compared with wild barley. We observed a mixture of completely and partially differentiated single‐nucleotide polymorphisms (SNPs) between domesticated and wild gene pools, suggesting that domesticated gene pools were derived from multiple wild ancestors. Patterns of genome‐wide linkage disequilibrium, haplotype block size and number, and variant frequency within blocks showed clear contrasts among individual chromosomes and between cultivars and wild barleys. Although most cultivar chromosomes shared a single major pericentromeric haplotype, chromosome 7H clearly differentiated the two‐row and six‐row types associated with different geographical origins. Within the pericentromeric regions we identified 22 387 non‐synonymous SNPs, 92 of which were fixed for alternative alleles in cultivar versus wild accessions. Surprisingly, only 29 SNPs found exclusively in the cultivars were predicted to be ‘highly deleterious’. Overall, our data reveal an unconventional pericentromeric genetic landscape among distinct barley gene pools, with different evolutionary processes driving domestication and diversification.

Using SPAdes, AUGUSTUS, and BLAST in an Automated Pipeline for Clustering Homologous Exome Sequences.

Cross-species exome sequencing approaches provide unprecedented avenues for obtaining genetic diversity, evolutionary relationships, and functional information from a variety of organisms including non-model species. These approaches offer cost-effective opportunities to study multiple individuals or species in parallel, but also create bioinformatics challenges in the application of multiple but powerful bioinformatics tools for the identification of homologous gene families across individual or species boundaries. Popular tools of this kind include SPAdes for sequence assembly, AUGUSTUS for ab initio gene prediction, and BLAST for building homologous sequence families. These tools can also be sophisticated in terms of installation and usage. Here, we present detailed steps on how to run these tools for the recovery and clustering of exon sequences from cross-species raw exome-capture data into homologous sequence families. We also present a utility pipeline, CODSEQCP, that automates these steps to cluster exon sequences, facilitating population genomics and evolutionary studies. © 2022 Wiley Periodicals LLC. Basic Protocol 1: Reads assembly using SPAdes Basic Protocol 2: Coding sequence extraction using AUGUSTUS Basic Protocol 3: Sequence clustering using BLAST Alternate Protocol: How to run CODSEQCP.

High-resolution mapping of Yr78, an adult plant resistance gene to wheat stripe rust.

Wheat stripe rust, caused by Puccinia striiformis f. sp. tritici (Pst), is responsible for significant yield losses worldwide, which can be minimized by the deployment of Pst resistance genes. Yr78 is an adult plant partial-resistance gene that has remained effective against the post-2000 virulent Pst races. In this study, we generated a high-resolution map of Yr78 based on 6,124 segregating chromosomes. We mapped Yr78 within a 0.05-cM interval on the short arm of chromosome 6B, which corresponds to an 11.16 Mb region between TraesCS6B02G116200 and TraesCS6B02G118000 in the 'Chinese Spring' Ref Seq. v1.1 genome. This interval is likely larger because it includes the unassembled NOR-B2 region, which may have contributed to the low recombination rate detected in this region. The Yr78 candidate region includes 15 genes that were prioritized for future functional studies based on their annotated function and polymorphisms between susceptible and resistant genotypes. Using exome capture data, we identified five major haplotypes in the candidate gene region, with the H1 haplotype associated with Yr78. The H1 haplotype was not detected in tetraploid wheat (Triticum turgidum L.) but was found in ∼30% of the common wheat cultivars (Triticum aestivum L.), suggesting that the associated resistance to stripe rust may have favored the selection of this haplotype. We developed two diagnostic molecular markers for the H1 haplotype that will facilitate the deployment of Yr78 in wheat breeding programs.

Application of the 3\u2032 mRNA-Seq using unique molecular identifiers in highly degraded RNA derived from formalin-fixed, paraffin-embedded tissue

BackgroundArchival formalin-fixed, paraffin-embedded (FFPE) tissue samples with clinical and histological data are a singularly valuable resource for developing new molecular biomarkers. However, transcriptome analysis remains challenging with standard mRNA-seq methods as FFPE derived-RNA samples are often highly modified and fragmented. The recently developed 3′ mRNA-seq method sequences the 3′ region of mRNA using unique molecular identifiers (UMI), thus generating gene expression data with minimal PCR bias. In this study, we evaluated the performance of 3′ mRNA-Seq using Lexogen QuantSeq 3′ mRNA-Seq Library Prep Kit FWD with UMI, comparing with TruSeq Stranded mRNA-Seq and RNA Exome Capture kit. The fresh-frozen (FF) and FFPE tissues yielded nucleotide sizes range from 13 to > 70% of DV200 values; input amounts ranged from 1 ng to 100 ng for validation.ResultsThe total mapped reads of QuantSeq 3′ mRNA-Seq to the reference genome ranged from 99 to 74% across all samples. After PCR bias correction, 3 to 56% of total sequenced reads were retained. QuantSeq 3′ mRNA-Seq data showed highly reproducible data across replicates in Universal Human Reference RNA (UHR, R > 0.94) at input amounts from 1 ng to 100 ng, and FF and FFPE paired samples (R = 0.92) at 10 ng. Severely degraded FFPE RNA with ≤30% of DV200 value showed good concordance (R > 0.87) with 100 ng input. A moderate correlation was observed when directly comparing QuantSeq 3′ mRNA-Seq data with TruSeq Stranded mRNA-Seq (R = 0.78) and RNA Exome Capture data (R > 0.67).ConclusionIn this study, QuantSeq 3′ mRNA-Seq with PCR bias correction using UMI is shown to be a suitable method for gene quantification in both FF and FFPE RNAs. 3′ mRNA-Seq with UMI may be applied to severely degraded RNA from FFPE tissues generating high-quality sequencing data.

Phylogenomic Analysis Reveals Dispersal-Driven Speciation and Divergence with Gene Flow in Lesser Sunda Flying Lizards (Genus Draco).

The Lesser Sunda Archipelago offers exceptional potential as a model system for studying the dynamics of dispersal-driven diversification. The geographic proximity of the islands suggests the possibility for successful dispersal, but this is countered by the permanence of the marine barriers and extreme intervening currents that are expected to hinder gene flow. Phylogenetic and species delimitation analyses of flying lizards (genus Draco) using single mitochondrial genes, complete mitochondrial genomes, and exome-capture data sets identified 9-11 deeply divergent lineages including single-island endemics, lineages that span multiple islands, and parapatrically distributed nonsister lineages on the larger islands. Population clustering and PCA confirmed these genetic boundaries with isolation-by-distance playing a role in some islands or island sets. While gdi estimates place most candidate species comparisons in the ambiguous zone, migration estimates suggest 9 or 10 species exist with nuclear introgression detected across some intra-island contact zones. Initial entry of Draco into the archipelago occurred at 5.5-7.5 Ma, with most inter-island colonization events having occurred between 1-3 Ma. Biogeographical model testing favors scenarios integrating geographic distance and historical island connectivity, including an initial stepping-stone dispersal process from the Greater Sunda Shelf through the Sunda Arc as far eastward as Lembata Island. However, rather than reaching the adjacent island of Pantar by dispersing over the 15-km wide Alor Strait, Draco ultimately reached Pantar (and much of the rest of the archipelago) by way of a circuitous route involving at least five overwater dispersal events. These findings suggest that historical geological and oceanographic conditions heavily influenced dispersal pathways and gene flow, which in turn drove species formation and shaped species boundaries. [Biogeography; genomics, Indonesia; lizards; phylogeography; reptiles].

Genome-wide association study reveals the genetic complexity of fructan accumulation patterns in barley grain.

We profiled the grain oligosaccharide content of 154 two-row spring barley genotypes and quantified 27 compounds, mainly inulin- and neoseries-type fructans, showing differential abundance. Clustering revealed two profile groups where the 'high' set contained greater amounts of sugar monomers, sucrose, and overall fructans, but lower fructosylraffinose. A genome-wide association study (GWAS) identified a significant association for the variability of two fructan types: neoseries-DP7 and inulin-DP9, which showed increased strength when applying a novel compound ratio-GWAS approach. Gene models within this region included three known fructan biosynthesis genes (fructan:fructan 1-fructosyltransferase, sucrose:sucrose 1-fructosyltransferase, and sucrose:fructan 6-fructosyltransferase). Two other genes in this region, 6(G)-fructosyltransferase and vacuolar invertase1, have not previously been linked to fructan biosynthesis and showed expression patterns distinct from those of the other three genes, including exclusive expression of 6(G)-fructosyltransferase in outer grain tissues at the storage phase. From exome capture data, several single nucleotide polymorphisms related to inulin- and neoseries-type fructan variability were identified in fructan:fructan 1-fructosyltransferase and 6(G)-fructosyltransferase genes. Co-expression analyses uncovered potential regulators of fructan biosynthesis including transcription factors. Our results provide the first scientific evidence for the distinct biosynthesis of neoseries-type fructans during barley grain maturation and reveal novel gene candidates likely to be involved in the differential biosynthesis of various types of fructan in barley.

WGVD: an integrated web-database for wheat genome variation and selective signatures.

Bread wheat is one of the most important crops worldwide. With the release of the complete wheat reference genome and the development of next-generation sequencing technology, a mass of genomic data from bread wheat and its progenitors has been yield and has provided genomic resources for wheat genetics research. To conveniently and effectively access and use these data, we established Wheat Genome Variation Database, an integrated web-database including genomic variations from whole-genome resequencing and exome-capture data for bread wheat and its progenitors, as well as selective signatures during the process of wheat domestication and improvement. In this version, WGVD contains 7 346 814 single nucleotide polymorphisms (SNPs) and 1 044 400 indels focusing on genic regions and upstream or downstream regions. We provide allele frequency distribution patterns of these variations for 5 ploidy wheat groups or 17 worldwide bread wheat groups, the annotation of the variant types and the genotypes of all individuals for 2 versions of bread wheat reference genome (IWGSC RefSeq v1.0 and IWGSC RefSeq v2.0). Selective footprints for Aegilops tauschii, wild emmer, domesticated emmer, bread wheat landrace and bread wheat variety are evaluated with two statistical tests (FST and Pi) based on SNPs from whole-genome resequencing data. In addition, we provide the Genome Browser to visualize the genomic variations, the selective footprints, the genotype patterns and the read coverage depth, and the alignment tool Blast to search the homologous regions between sequences. All of these features of WGVD will promote wheat functional studies and wheat breeding.Database URL http://animal.nwsuaf.edu.cn/code/index.php/Wheat

Gene editing of the wheat homologs of TONNEAU1-recruiting motif encoding gene affects grain shape and weight in wheat.

SummaryGrain size and weight are important components of a suite of yield‐related traits in crops. Here, we showed that the CRISPR‐Cas9 gene editing of TaGW7, a homolog of rice OsGW7 encoding a TONNEAU1‐recruiting motif (TRM) protein, affects grain shape and weight in allohexaploid wheat. By editing the TaGW7 homoeologs in the B and D genomes, we showed that mutations in either of the two or both genomes increased the grain width and weight but reduced the grain length. The effect sizes of mutations in the TaGW7 gene homoeologs coincided with the relative levels of their expression in the B and D genomes. The effects of gene editing on grain morphology and weight traits were dosage dependent with the double‐copy mutant showing larger effect than the respective single copy mutants. The TaGW7‐centered gene co‐expression network indicated that this gene is involved in the pathways regulating cell division and organ growth, also confirmed by the cellular co‐localization of TaGW7 with α‐ and β‐tubulin proteins, the building blocks of microtubule arrays. The analyses of exome capture data in tetraploid domesticated and wild emmer, and hexaploid wheat revealed the loss of diversity around TaGW7‐associated with domestication selection, suggesting that TaGW7 is likely to play an important role in the evolution of yield component traits in wheat. Our study showed how integrating CRISPR‐Cas9 system with cross‐species comparison can help to uncover the function of a gene fixed in wheat for allelic variants targeted by domestication selection and select targets for engineering new gene variants for crop improvement.

Identification of a Dominant Chlorosis Phenotype Through a Forward Screen of the Triticum turgidum cv. Kronos TILLING Population.

Durum wheat (Triticum turgidum) derives from a hybridization event approximately 400,000 years ago which led to the creation of an allotetraploid genome. The evolutionary recent origin of durum wheat means that its genome has not yet been fully diploidised. As a result, many of the genes present in the durum genome act in a redundant fashion, where loss-of-function mutations must be present in both gene copies to observe a phenotypic effect. Here, we use a novel set of induced variation within the cv. Kronos TILLING population to identify a locus controlling a dominant, environmentally dependent chlorosis phenotype. We carried out a forward screen of the sequenced cv. Kronos TILLING lines for senescence phenotypes and identified a line with a dominant early senescence and chlorosis phenotype. Mutant plants contained less chlorophyll throughout their development and displayed premature flag leaf senescence. A segregating population was classified into discrete phenotypic groups and subjected to bulked-segregant analysis using exome capture followed by next-generation sequencing. This allowed the identification of a single region on chromosome 3A, Yellow Early Senescence 1 (YES-1), which was associated with the mutant phenotype. While this phenotype was consistent across 4 years of field trials in the United Kingdom, the mutant phenotype was not observed when grown in Davis, CA (United States). To obtain further SNPs for fine-mapping, we isolated chromosome 3A using flow sorting and sequenced the entire chromosome. By mapping these reads against both the cv. Chinese Spring reference sequence and the cv. Kronos assembly, we could identify high-quality, novel EMS-induced SNPs in non-coding regions within YES-1 that were previously missed in the exome capture data. This allowed us to fine-map YES-1 to 4.3 Mb, containing 59 genes. Our study shows that populations containing induced variation can be sources of novel dominant variation in polyploid crop species, highlighting their importance in future genetic screens. We also demonstrate the value of using cultivar-specific genome assemblies alongside the gold-standard reference genomes particularly when working with non-coding regions of the genome. Further fine-mapping of the YES-1 locus will be pursued to identify the causal SNP underpinning this dominant, environmentally dependent phenotype.

Frequent intra- and inter-species introgression shapes the landscape of genetic variation in bread wheat

BackgroundBread wheat is one of the most important and broadly studied crops. However, due to the complexity of its genome and incomplete genome collection of wild populations, the bread wheat genome landscape and domestication history remain elusive.ResultsBy investigating the whole-genome resequencing data of 93 accessions from worldwide populations of bread wheat and its diploid and tetraploid progenitors, together with 90 published exome-capture data, we find that the B subgenome has more variations than A and D subgenomes, including SNPs and deletions. Population genetics analyses support a monophyletic origin of domesticated wheat from wild emmer in northern Levant, with substantial introgressed genomic fragments from southern Levant. Southern Levant contributes more than 676 Mb in AB subgenomes and enriched in the pericentromeric regions. The AB subgenome introgression happens at the early stage of wheat speciation and partially contributes to their greater genetic diversity. Furthermore, we detect massive alien introgressions that originated from distant species through natural and artificial hybridizations, resulting in the reintroduction of ~ 709 Mb and ~ 1577 Mb sequences into bread wheat landraces and varieties, respectively. A large fraction of these intra- and inter-introgression fragments are associated with quantitative trait loci of important traits, and selection events are also identified.ConclusionWe reveal the significance of multiple introgressions from distant wild populations and alien species in shaping the genetic components of bread wheat, and provide important resources and new perspectives for future wheat breeding.

A High-Resolution Map of Wheat QYr.ucw-1BL, an Adult Plant Stripe Rust Resistance Locus in the Same Chromosomal Region as Yr29.

The appearance of highly virulent and more aggressive races of f. sp. () during the last two decades has led to stripe rust epidemics worldwide and to the rapid erosion of effective resistance genes. In this study, we mapped an adult-plant resistance locus from the Argentinean wheat ( L.) cultivar Klein Chajá, which is effective against these new races. By using wheat exome capture data and a large population of 2480 segregating plants (4960 gametes), we mapped within a 0.24-cM region [332 kb in International Wheat Genome Sequencing Consortium (IWGSC) RefSeq version 1.0] on chromosome arm 1BL. This region overlaps with current maps of the adult-plant resistance gene , which has remained effective for more than 60 yr. An allelism test failed to find recombination between and and yielded similar resistance phenotypes for the two loci. These results, together with similar haplotypes in the candidate region, suggested that and might represent the same gene. However, we cannot rule out the possibility of tightly linked but different genes because most of the 13 genes in the candidate region are annotated with functions associated with disease resistance. To evaluate their potential as candidate genes, we characterized their polymorphisms between resistant and susceptible haplotypes. Finally, we used these polymorphisms to develop high-throughput markers to accelerate the deployment of these resistance loci in wheat breeding programs.

Development and Evaluation of a Barley 50k iSelect SNP Array.

High-throughput genotyping arrays continue to be an attractive, cost-effective alternative to sequencing based approaches. We have developed a new 50k Illumina Infinium iSelect genotyping array for barley, a cereal crop species of major international importance. The majority of SNPs on the array have been extracted from variants called in exome capture data of a wide range of European barley germplasm. We used the recently published barley pseudomolecule assembly to map the exome capture data, which allowed us to generate markers with accurate physical positions and detailed gene annotation. Markers from an existing and widely used barley 9k Infinium iSelect array were carried over onto the 50k chip for backward compatibility. The array design featured 49,267 SNP markers that converted into 44,040 working assays, of which 43,461 were scorable in GenomeStudio. Of the working assays, 6,251 are from the 9k iSelect platform. We validated the SNPs by comparing the genotype calls from the new array to legacy datasets. Rates of agreement averaged 98.1 and 93.9% respectively for the legacy 9k iSelect SNP set (Comadran et al., 2012) and the exome capture SNPs. To test the utility of the 50k chip for genetic mapping, we genotyped a segregating population derived from a Golden Promise × Morex cross (Liu et al., 2014) and mapped over 14,000 SNPs to genetic positions which showed a near exact correspondence to their known physical positions. Manual adjustment of the cluster files used by the interpreting software for genotype scoring improved results substantially, but migration of cluster files between sites led to a deterioration of results, suggesting that local adjustment of cluster files is required on a site-per-site basis. Information relating to the markers on the chip is available online at https://ics.hutton.ac.uk/50k.

Haplotype Analysis of the Pre-harvest Sprouting Resistance Locus Phs-A1 Reveals a Causal Role of TaMKK3-A in Global Germplasm.

Pre-harvest sprouting (PHS) is an important cause of quality loss in many cereal crops and is particularly prevalent and damaging in wheat. Resistance to PHS is therefore a valuable target trait in many breeding programs. The Phs-A1 locus on wheat chromosome arm 4AL has been consistently shown to account for a significant proportion of natural variation to PHS in diverse mapping populations. However, the deployment of sprouting resistance is confounded by the fact that different candidate genes, including the tandem duplicated Plasma Membrane 19 (PM19) genes and the mitogen-activated protein kinase kinase 3 (TaMKK3-A) gene, have been proposed to underlie Phs-A1. To further define the Phs-A1 locus, we constructed a physical map across this interval in hexaploid and tetraploid wheat. We established close proximity of the proposed candidate genes which are located within a 1.2 Mb interval. Genetic characterization of diverse germplasm used in previous genetic mapping studies suggests that TaMKK3-A, and not PM19, is the major gene underlying the Phs-A1 effect in European, North American, Australian and Asian germplasm. We identified the non-dormant TaMKK3-A allele at low frequencies within the A-genome diploid progenitor Triticum urartu genepool, and show an increase in the allele frequency in modern varieties. In United Kingdom varieties, the frequency of the dormant TaMKK3-A allele was significantly higher in bread-making quality varieties compared to feed and biscuit-making cultivars. Analysis of exome capture data from 58 diverse hexaploid wheat accessions identified fourteen haplotypes across the extended Phs-A1 locus and four haplotypes for TaMKK3-A. Analysis of these haplotypes in a collection of United Kingdom and Australian cultivars revealed distinct major dormant and non-dormant Phs-A1 haplotypes in each country, which were either rare or absent in the opposing germplasm set. The diagnostic markers and haplotype information reported in the study will help inform the choice of germplasm and breeding strategies for the deployment of Phs-A1 resistance into breeding germplasm.

Recombination Rate Variation, Hitchhiking, and Demographic History Shape Deleterious Load in Poplar.

Deleterious alleles are expected to be purged by purifying selection or maintained at low frequency. However, many additional evolutionary forces may shape the pattern of deleterious mutations across the genome and among populations, including selection, hitchhiking, recombination, and demographic history. We used exome capture data to estimate the genome-wide distribution of deleterious alleles across natural populations of the model tree black cottonwood (Populus trichocarpa). Although deleterious alleles were on average present at low frequency suggesting purifying selection, they were preferentially enriched both within genomic regions of low-recombination and in regions showing evidence of positive selection. The demographic history of this species also appeared to play a role in the distribution of deleterious alleles among populations, with peripheral populations having higher rates of deleterious homozygosity. This be due to less efficient selection arising from smaller effective population sizes at the edges of the range, and possibly also due to recent bottlenecks associated with postglacial recolonization. Finally, correlations between deleterious homozygosity and plant growth suggest a significant effect of deleterious load on fitness. Our results show that both genomic context and historical demography play a role in shaping the distribution of deleterious alleles across the genome and range of P. trichocarpa.

Somatic mutation detection using ensemble of flexible neural tree model

The advances on next-generation sequencing technology (NGS) have enabled researchers to detect somatic mutations. Much effect has been devoted to improve accuracy of discovering somatic mutations from tumour/normal NGS data. In this study, flexible neural tree model (FNT) is proposed to detect somatic mutations in tumour-normal paired sequencing data. To improve the classification accuracy further, a new classification ensemble approach based on Radial Basis Function (RBF) neural networks as nonlinear combination function is proposed. The proposed method is applied to real biological dataset from exome capture data and the whole genome shotgun data. Results show that the obtained FNT model has a fewer number of variables with reduced number of input features and with significant improvement in the detection accuracy using the proposed ensemble learning method. Our method also selects 10 import features for somatic mutation detection, which could be used to analyze NGS mutations further.