Polymorphic Single Nucleotide Polymorphism Research Articles

Genotyping of Soybean Cultivars With Medium-Density Array Reveals the Population Structure and QTNs Underlying Maturity and Seed Traits.

Soybean was domesticated about 5,000 to 6,000 years ago in China. Although genotyping technologies such as genotyping by sequencing (GBS) and high-density array are available, it is convenient and economical to genotype cultivars or populations using medium-density SNP array in genetic study as well as in molecular breeding. In this study, 235 cultivars, collected from China, Japan, USA, Canada and some other countries, were genotyped using SoySNP8k iSelect BeadChip with 7,189 single nucleotide polymorphisms (SNPs). In total, 4,471 polymorphic SNP markers were used to analyze population structure and perform genome-wide association study (GWAS). The most likely K value was 7, indicating this population can be divided into 7 subpopulations, which is well in accordance with the geographic origins of cultivars or accession studied. The LD decay rate was estimated at 184 kb, where r2 dropped to half of its maximum value (0.205). GWAS using FarmCPU detected a stable quantitative trait nucleotide (QTN) for hilum color and seed color, which is consistent with the known loci or genes. Although no universal QTNs for flowering time and maturity were identified across all environments, a total of 30 consistent QTNs were detected for flowering time (R1) or maturity (R7 and R8) on 16 chromosomes, most of them were corresponding to known E1 to E4 genes or QTL region reported in SoyBase (soybase.org). Of 16 consistent QTNs for protein and oil contents, 11 QTNs were detected having antagonistic effects on protein and oil content, while 4 QTNs soly for oil content, and one QTN soly for protein content. The information gained in this study demonstrated that the usefulness of the medium-density SNP array in genotyping for genetic study and molecular breeding.

Whole Genome Resequencing of Capsicum baccatum and Capsicum annuum to Discover Single Nucleotide Polymorphism Related to Powdery Mildew Resistance

The present study deals with genome wide identification of single-nucleotide polymorphism (SNP) markers related to powdery mildew (PM) resistance in two pepper varieties. Capsicum baccatum (PRH1- a PM resistant line) and Capsicum annuum (Saengryeg- a PM susceptible line), were resequenced to develop SNP markers. A total of 6,213,009 and 6,840,889 SNPs for PRH1 and Saengryeg respectively have been discovered. Among the SNPs, majority were classified as homozygous type SNPs, particularly in the resistant line. Moreover, the SNPs were differentially distributed among the chromosomes in both the resistant and susceptible lines. In total, 4,887,031 polymorphic SNP loci were identified between the two lines and 306,871 high-resolution melting (HRM) marker primer sets were designed. In order to understand the SNPs associated with the vital genes involved in diseases resistance and stress associated processes, chromosome-wise gene ontology analysis was performed. The results revealed the occurrence that SNPs related to diseases resistance genes were predominantly distributed in chromosome 4. In addition, 6281 SNPs associated with 46 resistance genes were identified. Among the lines, PRH1 consisted of maximum number of polymorphic SNPs related to NBS-LRR genes. The SNP markers were validated using HRM assay in 45 F4 populations and correlated with the phenotypic disease index.

High-density genetic mapping of a major QTL for resistance to multiple races of loose smut in a tetraploid wheat cross.

Loose smut, caused by Ustilago tritici (Pers.) Rostr., is a systemic disease of tetraploid durum wheat (Triticum turgidum L.). Loose smut can be economically controlled by growing resistant varieties, making it important to find and deploy new sources of resistance. Blackbird, a variety of T. turgidum L. subsp. carthlicum (Nevski) A. Love & D. Love, carries a high level of resistance to loose smut. Blackbird was crossed with the loose smut susceptible durum cultivar Strongfield to produce a doubled haploid (DH) mapping population. The parents and progenies were inoculated with U. tritici races T26, T32 and T33 individually and as a mixture at Swift Current, Canada in 2011 and 2012 and loose smut incidence (LSI) was assessed. Genotyping of the DH population and parents using an Infinium iSelect 90K single nucleotide polymorphism (SNP) array identified 12,952 polymorphic SNPs. The SNPs and 426 SSRs (previously genotyped in the same population) were mapped to 16 linkage groups spanning 3008.4 cM at an average inter-marker space of 0.2 cM in a high-density genetic map. Composite interval mapping analysis revealed three significant quantitative trait loci (QTL) for loose smut resistance on chromosomes 3A, 6B and 7A. The loose smut resistance QTL on 6B (QUt.spa-6B.2) and 7A (QUt.spa-7A.2) were derived from Blackbird. Strongfield contributed the minor QTL on 3A (QUt.spa-3A.2). The resistance on 6B was a stable major QTL effective against all individual races and the mixture of the three races; it explained up to 74% of the phenotypic variation. This study is the first attempt in durum wheat to identify and map loose smut resistance QTL using a high-density genetic map. The QTL QUt.spa-6B.2 would be an effective source for breeding resistance to multiple races of the loose smut pathogen because it provides near-complete broad resistance to the predominant virulence on the Canadian prairies.

Kompetitive Allele Specific PCR (KASP™) genotyping of 48 polymorphisms at different caprine loci in French Alpine and Saanen goat breeds and their association with milk composition

Using a novel and fast genotyping method called Kompetitive Allele Specific PCR (KASP™), we carried out a pilot study on 48 single nucleotide polymorphisms (SNPs) belonging to 40 genes in French Alpine (n = 24) and Saanen (n = 25) goats reared in Romania. Furthermore, the associations of the 13 polymorphic genetic variants with milk production and composition were investigated. Thirty-five SNPs did not show polymorphism in the studied populations. Polymorphic SNPs were detected in the following genes: CAST, CLEC4E, DES, GHRHR, HSP90AA1, IL15RA, IL1RN, IL8, MITF, PPRC1, SOCS3, TNF and TNFSF13. The studied Alpine population was in Hardy-Weinberg disequilibrium at the g.62894878A>G locus (rs671391101) (P < 0.05). The results showed that four SNPs rs671391101 (GHRHR), rs640582069 (IL1RN) rs635583012 (SOCS3) and rs635969404 (IL15RA) out of the 13 polymorphic markers were significantly associated with milk production, protein, fat and lactose content in the Alpine breed. However, no significant effect was recorded in the Saanen population regarding milk yield or milk chemical composition. The current results provide new insights for the development of SNP marker-assisted selection technology in the goat industry and confirm the potential of using SNPs for the GHRHR, IL1RN, SOCS3, and IL15RA genes as candidate genes for selection, highlighting the direct implications of such genes for farm production outputs. The results from this study are relevant for future goat genomic studies and the inclusion of the associated traits into up-to-date selection schemes.

Development of Molecular Markers Linked to Powdery Mildew Resistance Gene Pm4b by Combining SNP Discovery from Transcriptome Sequencing Data with Bulked Segregant Analysis (BSR-Seq) in Wheat.

Powdery mildew resistance gene Pm4b, originating from Triticum persicum, is effective against the prevalent Blumeria graminis f. sp. tritici (Bgt) isolates from certain regions of wheat production in China. The lack of tightly linked molecular markers with the target gene prevents the precise identification of Pm4b during the application of molecular marker-assisted selection (MAS). The strategy that combines the RNA-Seq technique and the bulked segregant analysis (BSR-Seq) was applied in an F2:3 mapping population (237 families) derived from a pair of isogenic lines VPM1/7∗Bainong 3217 F4 (carrying Pm4b) and Bainong 3217 to develop more closely linked molecular markers. RNA-Seq analysis of the two phenotypically contrasting RNA bulks prepared from the representative F2:3 families generated 20,745,939 and 25,867,480 high-quality read pairs, and 82.8 and 80.2% of them were uniquely mapped to the wheat whole genome draft assembly for the resistant and susceptible RNA bulks, respectively. Variant calling identified 283,866 raw single nucleotide polymorphisms (SNPs) and InDels between the two bulks. The SNPs that were closely associated with the powdery mildew resistance were concentrated on chromosome 2AL. Among the 84 variants that were potentially associated with the disease resistance trait, 46 variants were enriched in an about 25 Mb region at the distal end of chromosome arm 2AL. Four Pm4b-linked SNP markers were developed from these variants. Based on the sequences of Chinese Spring where these polymorphic SNPs were located, 98 SSR primer pairs were designed to develop distal markers flanking the Pm4b gene. Three SSR markers, Xics13, Xics43, and Xics76, were incorporated in the new genetic linkage map, which located Pm4b in a 3.0 cM genetic interval spanning a 6.7 Mb physical genomic region. This region had a collinear relationship with Brachypodium distachyon chromosome 5, rice chromosome 4, and sorghum chromosome 6. Seven genes associated with disease resistance were predicted in this collinear genomic region, which included C2 domain protein, peroxidase activity protein, protein kinases of PKc_like super family, Mlo family protein, and catalytic domain of the serine/threonine kinases (STKc_IRAK like super family). The markers developed in the present study facilitate identification of Pm4b during its MAS practice.

The differential expression of the blood group P1 -A4GALT and P2 -A4GALT alleles is stimulated by the transcription factor early growth response 1.

The P1 /P2 phenotypic polymorphism is one of the earliest blood groups discovered in humans. These blood groups have been connected to different levels of expression of the A4GALT gene in P1 and P2 red blood cells; however, the detailed molecular genetic mechanism that leads to these two phenotypes has not been established. After our previous identification of an association between the single-nucleotide polymorphisms (SNPs) rs2143918 and rs5751348 in A4GALT gene and the P1 /P2 phenotype, we conduct a survey of transcription factors that might connect these SNPs with the differential expression of the P1 -A4GALT and P2 -A4GALT alleles. An in silico analysis of potential transcription factor binding motifs within the polymorphic SNPs rs2143918 and rs5751348 genomic regions was performed, and this was followed by reporter assays examining the candidate transcription factors, gene expression profiling, electrophoretic mobility shift assays, and P1 -A4GALT and P2 -A4GALT allelic expression analysis. The results revealed that the differential binding of transcription factor early growth response 1 to the SNP rs5751348 genomic region with the different genotypes in the A4GALT gene leads to differential activation of P1 -A4GALT and P2 -A4GALT expression. The present investigation, together with our previous study (Lai et al., Transfusion 2014;54:3222-31), have elucidated the molecular genetic details associated with the P1 /P2 blood groups.

Genome-wide population structure and admixture analysis reveals weak differentiation among Ugandan goat breeds.

SummaryUganda has a large population of goats, predominantly from indigenous breeds reared in diverse production systems, whose existence is threatened by crossbreeding with exotic Boer goats. Knowledge about the genetic characteristics and relationships among these Ugandan goat breeds and the potential admixture with Boer goats is still limited. Using a medium‐density single nucleotide polymorphism (SNP) panel, we assessed the genetic diversity, population structure and admixture in six goat breeds in Uganda: Boer, Karamojong, Kigezi, Mubende, Small East African and Sebei. All the animals had genotypes for about 46 105 SNPs after quality control. We found high proportions of polymorphic SNPs ranging from 0.885 (Kigezi) to 0.928 (Sebei). The overall mean observed (HO) and expected (HE) heterozygosity across breeds was 0.355 ± 0.147 and 0.384 ± 0.143 respectively. Principal components, genetic distances and admixture analyses revealed weak population sub‐structuring among the breeds. Principal components separated Kigezi and weakly Small East African from other indigenous goats. Sebei and Karamojong were tightly entangled together, whereas Mubende occupied a more central position with high admixture from all other local breeds. The Boer breed showed a unique cluster from the Ugandan indigenous goat breeds. The results reflect common ancestry but also some level of geographical differentiation. admixture and f4 statistics revealed gene flow from Boer and varying levels of genetic admixture among the breeds. Generally, moderate to high levels of genetic variability were observed. Our findings provide useful insights into maintaining genetic diversity and designing appropriate breeding programs to exploit within‐breed diversity and heterozygote advantage in crossbreeding schemes.

Informativeness of Single Nucleotide Polymorphisms and Relationships among Onion Populations from Important World Production Regions

Single nucleotide polymorphisms (SNPs) were genotyped using a high-density array and DNAs from individual plants of important onion (Allium cepa L.) populations from major production regions and from the likely progenitor of onion, Allium vavilovii Popov et Vved. Genotypes at 1226 SNPs were used to estimate genetic relationships among these populations and revealed close associations among onions grown in Europe and those in North America, South America, and eastern Asia, supporting paths of introduction from Europe to the Americas and Asia. ‘Nasik Red’ is a population grown on the Indian subcontinent and was divergent from onions of European origin. Frequencies of SNPs among and within populations were used as a measure of informativeness, and 199 commonly polymorphic SNPs were identified distributed across the eight chromosomes of onion. These SNPs will be useful for estimations of relatedness among broader collections of onion populations, mapping of important phenotypes, fingerprinting of inbred lines and hybrids, and quality control of seed lots.

Identification of the resistance gene to powdery mildew in Chinese wheat landrace Baiyouyantiao

Powdery mildew, caused by Blumeria graminis f. sp. tritici (Bgt), is one of the most damaging diseases to wheat in the world. The cultivation of resistant varieties of wheat is essential for controlling the powdery mildew epidemic. Wheat landraces are important resources of resistance to many diseases. Mapping powdery mildew resistance genes from wheat landraces will promote the development of new varieties with disease resistance. The Chinese wheat landrace Baiyouyantiao possesses characteristic of disease resistance to powdery mildew. To identify the resistance gene in this landrace, Baiyouyantiao was crossed with the susceptible cultivar Jingshuang 16 and seedlings of parents and F1, BC1, F2, and F2:3 were tested with Bgt isolate E09. The genetic results showed that the resistance of Baiyouyantiao to E09 was controlled by a single recessive gene, tentatively designated PmBYYT. An Illumina wheat 90K single-nucleotide polymorphism (SNP) array was applied to screen polymorphisms between F2-resistant and F2-susceptible DNA bulks for identifying the chromosomal location of PmBYYT. A high percentage of polymorphic SNPs between the resistant and susceptible DNA bulks was found on chromosome 7B, indicating that PmBYYT may be located on this chromosome. A genetic linkage map of PmBYYT consisting of two simple sequence repeat markers and eight SNP markers was developed. The two flanking markers were SNP markers W7BL-8 and W7BL-15, with genetic distances of 3 and 2.9 cM, respectively. The results of this study demonstrated the rapid characterization of a wheat disease resistance gene and SNP marker development using the 90K SNP assay. The flanking markers of gene PmBYYT will benefit marker-assisted selection (MAS) and map-based cloning in breeding wheat cultivars with powdery mildew resistance.

QTL Mapping of Ineffective Nodulation and Nitrogen Utilization‐Related Traits in the IC‐1 Mutant of Cowpea

ABSTRACTBiological nitrogen fixation is a valuable component of sustainable agriculture. Improved understanding of nodulation pathways can facilitate enhancement of nitrogen utilization‐related traits. Previously, a cowpea [Vigna unguiculata (L.) Walp.] mutant, ‘IC‐1’, with ineffective nodules was identified. Heritability studies indicated that ineffective nodulation in IC‐1 was controlled by a single gene designated cpi. To map cpi, F2 mapping populations were developed from crosses between IC‐1 and two commercial cowpea cultivars. The F2 populations were evaluated for leaf and nodule color, nitrogen content, dry plant weight, nodule number, and nodule fresh weight. Nitrogen content, aboveground plant dry weight, and nodule fresh weight were each significantly lower in individuals with ineffective nodulation, whereas nodule number was significantly higher. Ninety individuals from one population were selected and genotyped with the Cowpea iSelect Consortium single‐nucleotide polymorphism (SNP) array. In total, 6083 polymorphic SNPs were used for linkage map construction and quantitative trait loci (QTL) analysis. A major QTL accounting for nearly 50% of the phenotypic variance was identified on linkage group (LG) 4 associated with all traits except nodule fresh weight, whereas a minor QTL for nodule fresh weight was mapped to LG 6. The major QTL was validated in a second F2 population with several polymerase chain reaction markers. Identification of these QTL provides a foundation for future identification of the cpi gene and potential loci for targeted improvement of nitrogen‐related traits in cowpea.

Combining Single Nucleotide Polymorphism Genotyping Array with Bulked Segregant Analysis to Map a Gene Controlling Adult Plant Resistance to Stripe Rust in Wheat Line 03031-1-5 H62.

Stripe rust, caused by Puccinia striiformis f. sp. tritici, is one of the most devastating diseases of wheat worldwide. Growing resistant cultivars is considered the best approach to manage this disease. In order to identify the resistance gene(s) in wheat line 03031-1-5 H62, which displayed high resistance to stripe rust at adult plant stage, a cross was made between 03031-1-5 H62 and susceptible cultivar Avocet S. The mapping population was tested with Chinese P. striiformis f. sp. tritici race CYR32 through artificial inoculation in a field in Yangling, Shaanxi Province and under natural infection in Tianshui, Gansu Province. The segregation ratios indicated that the resistance was conferred by a single dominant gene, temporarily designated as YrH62. A combination of bulked segregant analysis (BSA) with wheat 90K single nucleotide polymorphism (SNP) array was used to identify molecular markers linked to YrH62. A total of 376 polymorphic SNP loci identified from the BSA analysis were located on chromosome 1B, from which 35 kompetitive allele-specific PCR (KASP) markers selected together with 84 simple sequence repeat (SSR) markers on 1B were used to screen polymorphism and a chromosome region associated with rust resistance was identified. To saturate the chromosomal region covering the YrH62 locus, a 660K SNP array was used to identify more SNP markers. To develop tightly linked markers for marker-assisted selection of YrH62 in wheat breeding, 18 SNPs were converted into KASP markers. A final linkage map consisting of 15 KASP and 3 SSR markers was constructed with KASP markers AX-109352427 and AX-109862469 flanking the YrH62 locus in a 1.0 cM interval. YrH62 explained 63.8 and 69.3% of the phenotypic variation for disease severity and infection type, respectively. YrH62 was located near the centromeric region of chromosome 1BS based on the positions of the SSR markers in 1B deletion bins. Based on the origin, responses to P. striiformis f. sp. tritici races, and marker distances, YrH62 is likely different from the other reported stripe rust resistance genes/quantitative trait loci on 1B. The gene and tightly linked KASP markers will be useful for breeding wheat cultivars with resistance to stripe rust.

Mapping QTLs using a novel source of salinity tolerance from Hasawi and their interaction with environments in rice

BackgroundSalinity is one of the most severe and widespread abiotic stresses that affect rice production. The identification of major-effect quantitative trait loci (QTLs) for traits related to salinity tolerance and understanding of QTL × environment interactions (QEIs) can help in more precise and faster development of salinity-tolerant rice varieties through marker-assisted breeding. Recombinant inbred lines (RILs) derived from IR29/Hasawi (a novel source of salinity) were screened for salinity tolerance in the IRRI phytotron in the Philippines (E1) and in two other diverse environments in Senegal (E2) and Tanzania (E3). QTLs were mapped for traits related to salinity tolerance at the seedling stage.ResultsThe RILs were genotyped using 194 polymorphic SNPs (single nucleotide polymorphisms). After removing segregation distortion markers (SDM), a total of 145 and 135 SNPs were used to construct a genetic linkage map with a length of 1655 and 1662 cM, with an average marker density of 11.4 cM in E1 and 12.3 cM in E2 and E3, respectively. A total of 34 QTLs were identified on 10 chromosomes for five traits using ICIM-ADD and segregation distortion locus (SDL) mapping (IM-ADD) under salinity stress across environments. Eight major genomic regions on chromosome 1 between 170 and 175 cM (qSES1.3, qSES1.4, qSL1.2, qSL1.3, qRL1.1, qRL1.2, qFWsht1.2, qDWsht1.2), chromosome 4 at 32 cM (qSES4.1, qFWsht4.2, qDWsht4.2), chromosome 6 at 115 cM (qFWsht6.1, qDWsht6.1), chromosome 8 at 105 cM (qFWsht8.1, qDWsht8.1), and chromosome 12 at 78 cM (qFWsht12.1, qDWsht12.1) have co-localized QTLs for the multiple traits that might be governing seedling stage salinity tolerance through multiple traits in different phenotyping environments, thus suggesting these as hot spots for tolerance of salinity. Forty-nine and 30 significant pair-wise epistatic interactions were detected between QTL-linked and QTL-unlinked regions using single-environment and multi-environment analyses.ConclusionsThe identification of genomic regions for salinity tolerance in the RILs showed that Hasawi possesses alleles that are novel for salinity tolerance. The common regions for the multiple QTLs across environments as co-localized regions on chromosomes 1, 4, 6, 8, and 12 could be due to linkage or pleiotropic effect, which might be helpful for multiple QTL introgression for marker-assisted breeding programs to improve the salinity tolerance of adaptive and popular but otherwise salinity-sensitive rice varieties.

GWAS for plant growth stages and yield components in spring wheat (Triticum aestivum L.) harvested in three regions of Kazakhstan

BackgroundSpring wheat is the largest agricultural crop grown in Kazakhstan with an annual sowing area of 12 million hectares in 2016. Annually, the country harvests around 15 million tons of high quality grain. Despite environmental stress factors it is predicted that the use of new technologies may lead to increases in productivity from current levels of 1.5 to up to 3 tons per hectare. One way of improving wheat productivity is by the application of new genomic oriented approaches in plant breeding projects. Genome wide association studies (GWAS) are emerging as powerful tools for the understanding of the inheritance of complex traits via utilization of high throughput genotyping technologies and phenotypic assessments of plant collections. In this study, phenotyping and genotyping data on 194 spring wheat accessions from Kazakhstan, Russia, Europe, and CIMMYT were assessed for the identification of marker-trait associations (MTA) of agronomic traits by using GWAS.ResultsField trials in Northern, Central and Southern regions of Kazakhstan using 194 spring wheat accessions revealed strong correlations of yield with booting date, plant height, biomass, number of spikes per plant, and number of kernels per spike. The accessions from Europe and CIMMYT showed high breeding potential for Southern and Central regions of the country in comparison with the performance of the local varieties. The GGE biplot method, using average yield per plant, suggested a clear separation of accessions into their three breeding origins in relationship to the three environments in which they were evaluated. The genetic variation in the three groups of accessions was further studied using 3245 polymorphic SNP (single nucleotide polymorphism) markers. The application of Principal Coordinate analysis clearly grouped the 194 accessions into three clades according to their breeding origins. GWAS on data from nine field trials allowed the identification of 114 MTAs for 12 different agronomic traits.ConclusionsField evaluation of foreign germplasm revealed its poor yield performance in Northern Kazakhstan, which is the main wheat growing region in the country. However, it was found that EU and CIMMYT germplasm has high breeding potential to improve yield performance in Central and Southern regions. The use of Principal Coordinate analysis clearly separated the panel into three distinct groups according to their breeding origin. GWAS based on use of the TASSEL 5.0 package allowed the identification of 114 MTAs for twelve agronomic traits. The study identifies a network of key genes for improvement of yield productivity in wheat growing regions of Kazakhstan.

The MGMT promoter single-nucleotide polymorphism rs1625649 had prognostic impact on patients with MGMT methylated glioblastoma

Promoter methylation is the most significant mechanism to regulate O6-methylguanine-DNA-methyltransferase (MGMT) expression. Single-nucleotide polymorphisms (SNPs) in the MGMT promoter region may also play a role. The aim of this study was to evaluate the clinical significance of SNPs in the MGMT promoter region of glioblastoma. Genomic DNAs from 118 glioblastomas were collected for polymerase chain reaction (PCR) amplification. Sanger sequencing was used to sequence the MGMT promoter region to detect SNPs. The results were correlated with MGMT status and patient survival. Rs1625649 was the only polymorphic SNP located at the MGMT promoter region in 37.5% of glioblastomas. Homozygous rs1625649 (AA genotype) was correlated with a higher MGMT methylation level and a lower protein expression, but the result was not statistically significant. In patients with MGMT methylated glioblastoma, cases with homozygous rs1625649 (AA genotype) were significantly associated with a lack of MGMT protein expression and a better progression-free survival (PFS) than the cases with wild type rs1625649 (CC genotype) or heterozygous rs1625649 (CA genotype). The survival impact was significant in multivariate analyses. In conclusion, the MGMT promoter homozygous rs1625649 (AA genotype) was found to correlate with a better PFS in patients with MGMT methylated glioblastoma.

High Density Linkage Map Construction and QTL Detection for Three Silique-Related Traits in Orychophragmus violaceus Derived Brassica napus Population.

Seeds per silique (SS), seed weight (SW), and silique length (SL) are important determinant traits of seed yield potential in rapeseed (Brassica napus L.), and are controlled by naturally occurring quantitative trait loci (QTLs). Mapping QTLs to narrow chromosomal regions provides an effective means of characterizing the genetic basis of these complex traits. Orychophragmus violaceus is a crucifer with long siliques, many SS, and heavy seeds. A novel B. napus introgression line with many SS was previously selected from multiple crosses (B. rapa ssp. chinesis × O. violaceus) × B. napus. In present study, a doubled haploid (DH) population with 167 lines was established from a cross between the introgression line and a line with far fewer SS, in order to detect QTLs for silique-related traits. By screening with a Brassica 60K single nucleotide polymorphism (SNP) array, a high-density linkage map consisting of 1,153 bins and spanning a cumulative length of 2,209.1 cM was constructed, using 12,602 high-quality polymorphic SNPs in the DH population. The average recombination bin densities of the A and C subgenomes were 1.7 and 2.4 cM, respectively. 45 QTLs were identified for the three traits in all, which explained 4.0–34.4% of the total phenotypic variation; 20 of them were integrated into three unique QTLs by meta-analysis. These unique QTLs revealed a significant positive correlation between SS and SL and a significant negative correlation between SW and SS, and were mapped onto the linkage groups A05, C08, and C09. A trait-by-trait meta-analysis revealed eight, four, and seven consensus QTLs for SS, SW, and SL, respectively, and five major QTLs (cqSS.A09b, cqSS.C09, cqSW.A05, cqSW.C09, and cqSL.C09) were identified. Five, three, and four QTLs for SS, SW, and SL, respectively, might be novel QTLs because of the existence of alien genetic loci for these traits in the alien introgression. Thirty-eight candidate genes underlying nine QTLs for silique-related traits were identified.

Genotyping-by-sequencing enables linkage mapping in three octoploid cultivated strawberry families.

Genotyping-by-sequencing (GBS) was used to survey genome-wide single-nucleotide polymorphisms (SNPs) in three biparental strawberry (Fragaria × ananassa) populations with the goal of evaluating this technique in a species with a complex octoploid genome. GBS sequence data were aligned to the F. vesca ‘Fvb’ reference genome in order to call SNPs. Numbers of polymorphic SNPs per population ranged from 1,163 to 3,190. Linkage maps consisting of 30–65 linkage groups were produced from the SNP sets derived from each parent. The linkage groups covered 99% of the Fvb reference genome, with three to seven linkage groups from a given parent aligned to any particular chromosome. A phylogenetic analysis performed using the POLiMAPS pipeline revealed linkage groups that were most similar to ancestral species F. vesca for each chromosome. Linkage groups that were most similar to a second ancestral species, F. iinumae, were only resolved for Fvb 4. The quantity of missing data and heterogeneity in genome coverage inherent in GBS complicated the analysis, but POLiMAPS resolved F. × ananassa chromosomal regions derived from diploid ancestor F. vesca.

Fine Mapping of a Clubroot Resistance Gene in Chinese Cabbage Using SNP Markers Identified from Bulked Segregant RNA Sequencing.

Clubroot, caused by Plasmodiophora brassicae, is an important disease of canola (Brassica napus) in western Canada and worldwide. In this study, a clubroot resistance gene (Rcr2) was identified and fine mapped in Chinese cabbage cv. “Jazz” using single-nucleotide polymorphisms (SNP) markers identified from bulked segregant RNA sequencing (BSR-Seq) and molecular markers were developed for use in marker assisted selection. In total, 203.9 million raw reads were generated from one pooled resistant (R) and one pooled susceptible (S) sample, and >173,000 polymorphic SNP sites were identified between the R and S samples. One significant peak was observed between 22 and 26 Mb of chromosome A03, which had been predicted by BSR-Seq to contain the causal gene Rcr2. There were 490 polymorphic SNP sites identified in the region. A segregating population consisting of 675 plants was analyzed with 15 SNP sites in the region using the Kompetitive Allele Specific PCR method, and Rcr2 was fine mapped between two SNP markers, SNP_A03_32 and SNP_A03_67 with 0.1 and 0.3 cM from Rcr2, respectively. Five SNP markers co-segregated with Rcr2 in this region. Variants were identified in 14 of 36 genes annotated in the Rcr2 target region. The numbers of poly variants differed among the genes. Four genes encode TIR-NBS-LRR proteins and two of them Bra019410 and Bra019413, had high numbers of polymorphic variants and so are the most likely candidates of Rcr2.

Toward Identification of Black Lemma and Pericarp Gene Blp1 in Barley Combining Bulked Segregant Analysis and Specific-Locus Amplified Fragment Sequencing.

Black barley is caused by phytomelanin synthesized in lemma and/or pericarp and the trait is controlled by one dominant gene Blp1. The gene is mapped on chromosome 1H by molecular markers, but it is yet to be isolated. Specific-locus amplified fragment sequencing (SLAF-seq) is an effective method for large-scale de novo single nucleotide polymorphism (SNP) discovery and genotyping. In the present study, SLAF-seq with bulked segregant analysis (BSA) was employed to obtain sufficient markers to fine mapping Blp1 gene in an F2 population derived from Hatiexi No.1 × Zhe5819. Based on SNP screening criteria, a total of 77,542 polymorphic SNPs met the requirements for association analysis. Combining two association analysis methods, the overlapped region with a size of 32.41 Mb on chromosome 1H was obtained as the candidate region of Blp1 gene. According to SLAF-seq data, markers were developed in the target region and were used for mapping the Blp1 gene. Linkage analysis showed that Blp1 co-segregated with HZSNP34 and HZSNP36, and was delimited by two markers (HZSNP35 and HZSNP39) spanning 8.1 cM in 172 homozygous yellow grain F2 plants of Hatiexi No.1 × Zhe5819. More polymorphic markers were screened in the reduced target region and were used to genotype the population. As a result, Blp1 was delimited within a 1.66 Mb on chromosome 1H by the upstream marker HZSNP63 and the downstream marker HZSNP59. Our results demonstrated the utility of SLAF-seq-BSA approach to identify the candidate region and discover polymorphic markers at the specific targeted genomic region.

Genome-wide association mapping reveals a rich genetic architecture of stripe rust resistance loci in emmer wheat (Triticum turgidum ssp. dicoccum)

Key messageSNP-based genome scanning in worldwide domesticated emmer germplasm showed high genetic diversity, rapid linkage disequilibrium decay and 51 loci for stripe rust resistance, a large proportion of which were novel.Cultivated emmer wheat (Triticum turgidum ssp. dicoccum), one of the oldest domesticated crops in the world, is a potentially rich reservoir of variation for improvement of resistance/tolerance to biotic and abiotic stresses in wheat. Resistance to stripe rust (Puccinia striiformis f. sp. tritici) in emmer wheat has been under-investigated. Here, we employed genome-wide association (GWAS) mapping with a mixed linear model to dissect effective stripe rust resistance loci in a worldwide collection of 176 cultivated emmer wheat accessions. Adult plants were tested in six environments and seedlings were evaluated with five races from the United States and one from Italy under greenhouse conditions. Five accessions were resistant across all experiments. The panel was genotyped with the wheat 90,000 Illumina iSelect single nucleotide polymorphism (SNP) array and 5106 polymorphic SNP markers with mapped positions were obtained. A high level of genetic diversity and fast linkage disequilibrium decay were observed. In total, we identified 14 loci associated with field resistance in multiple environments. Thirty-seven loci were significantly associated with all-stage (seedling) resistance and six of them were effective against multiple races. Of the 51 total loci, 29 were mapped distantly from previously reported stripe rust resistance genes or quantitative trait loci and represent newly discovered resistance loci. Our results suggest that GWAS is an effective method for characterizing genes in cultivated emmer wheat and confirm that emmer wheat is a rich source of stripe rust resistance loci that can be used for wheat improvement.

SNP discovery from next-generation transcriptome sequencing data and their validation using KASP assay in wheat (Triticum aestivum L.)

Single nucleotide polymorphisms (SNPs) are becoming the most amenable form of DNA-based molecular markers for genetic analysis. In hexaploid bread wheat (Triticum aestivum L.), it is difficult to discern true polymorphic SNPs due to homoeologous and paralogous genes. Two serial analysis of gene expression (SAGE) libraries were developed utilizing leaves from resistant plants carrying leaf rust resistance gene Lr28; one library was derived from leaves that were mock inoculated and the other was derived from leaves inoculated with the urediniospores of the leaf rust pathogen Puccinia triticina. Next-generation sequencing reads, after quality trimming and removal of fungal sequences, were mapped to wheat reference sequences at Ensembl Plants. CLC Genomics Workbench and Freebayes softwares were employed for SNP calling. A total of 611 SNPs were predicted to be common by both softwares, of which 207 varietal SNPs were identified by ConservedPrimer software. A subset of 100 SNPs was used for validation across 47 wheat genotypes using Kompetitive Allele Specific PCR (KASP) assay; 83 SNPs could be successfully validated. These SNPs were positioned on wheat subgenomes and chromosome arms. When functionally annotated, many sequences harboring SNPs showed homology to resistance and resistance-like genes listed in Plant Resistance Gene database (PRGdb) as well as pathogenesis-related (PR) and stress-responsive genes. The results of the present study involving discovery of SNPs associated with resistance to leaf rust, a major threat to wheat production worldwide, will be valuable for molecular breeding for rust resistance.