Genotyping-by-sequencing Research Articles

Genome-wide association study revealed significant SNPs for anthracnose resistance, seed alkaloids and protein content in white lupin

White lupin (Lupinus albus L.) is a high-protein grain legume alternative to soybean in Central Europe, but its cultivation is risky due to the fungal disease anthracnose that can cause severe yield damage. In addition, management of seed alkaloids is critical for human nutrition and animal feed. We report on a white lupin collection of genebank accessions, advanced breeding lines and cultivars that was genotyped and phenotypically characterized for anthracnose resistance and seed alkaloids and protein levels. Using genotyping by sequencing (GBS), SeqSNP-targeted GBS, BiomarkX genotyping and Sanger sequencing, a genetic resource of genome-wide SNPs for white lupin was established. We determined anthracnose resistance in two years field trials at four locations with infection rows and measured seed alkaloids and protein levels by near-infrared spectroscopy (NIRS). Few white lupin breeding lines showed anthracnose resistance comparable or better than Celina and Frieda, currently the best commercial cultivars in Germany. NIRS estimates for seed alkaloids and protein levels revealed variation in the white lupin collection. Using genome-wide association studies (GWAS), we identified SNPs significantly associated with anthracnose resistance in the field representing known and new genomic regions. We confirmed the pauper locus and detected new SNP markers significantly associated with seed alkaloids. For the first time, we present loci associated with total grain protein content. Finally, we tested the potential of genomic prediction (GP) in predicting the phenotype of these three quantitative traits. Application of results and resources are discussed in the context of fostering breeding programs for white lupin.

Uncovering genes involved in pollinator-driven mating system shifts and selfing syndrome evolution in Brassica rapa.

Shifts in pollinator occurrence and their pollen transport effectiveness drive the evolution of mating systems in flowering plants. Understanding the genomic basis of these changes is essential for predicting the persistence of a species under environmental changes. We investigated the genomic changes in Brassica rapa over nine generations of pollination by hoverflies associated with rapid morphological evolution toward the selfing syndrome. We combined a genotyping-by-sequencing (GBS) approach with a genome-wide association study (GWAS) to identify candidate genes, and assessed their functional role in the observed morphological changes by studying mutations of orthologous genes in the model plant Arabidopsis thaliana. We found 31 candidate genes involved in a wide range of functions from DNA/RNA binding to transport. Our functional assessment of orthologous genes in A. thaliana revealed that two of the identified genes in B. rapa are involved in regulating the size of floral organs. We found a protein kinase superfamily protein involved in petal width, an important trait in plant attractiveness to pollinators. Moreover, we found a histone lysine methyltransferase (HKMT) associated with stamen length. Altogether, our study shows that hoverfly pollination leads to rapid evolution toward the selfing syndrome mediated by polygenic changes.

Fine Mapping of QTLs for Alkaline Tolerance in Crucian Carp (Carassius auratus) Using Genome-Wide SNP Markers.

Crucian carp (Carassius auratus) is widely distributed in the world and has become an economically freshwater fish. The population in Lake Dali Nur can tolerate the extreme alkaline environment with alkalinity over 50 mmol/L (pH 9.6), thus providing a special model for exploring alkali-tolerant molecular markers in an extremely alkaline environment. In this study, we constructed a high-density and high-resolution linkage map with 16,224 SNP markers based on genotyping-by-sequencing (GBS) consisting of 152 progenies and conducted QTL studies for alkali-tolerant traits. The total length of the linkage map was 3918.893 cM, with an average distance of 0.241 cM. Two QTLs for the ammonia-N-tolerant trait were detected on LG27 and LG45. A QTL for the urea-N-tolerant trait was detected on LG27. Interestingly, mapping the two QTLs on LG27 revealed that the mapped genes were both located in the intron of CDC42. GO functional annotation and KEGG enrichment analysis results indicated that the biological functions might be involved in the cell cycle, cellular senescence, MAPK, and Ras signaling pathways. These findings suggest that CDC42 may play an important role in the process of dealing with extremely alkaline environments.

Genotyping-by-Sequencing Analysis Reveals Associations between Agronomic and Oil Traits in Gamma Ray-Derived Mutant Rapeseed (Brassica napus L.).

Rapeseed (Brassica napus L.) holds significant commercial value as one of the leading oil crops, with its agronomic features and oil quality being crucial determinants. In this investigation, 73,226 single nucleotide polymorphisms (SNPs) across 95 rapeseed mutant lines induced by gamma rays, alongside the original cultivar ('Tamra'), using genotyping-by-sequencing (GBS) analysis were examined. This study encompassed gene ontology (GO) analysis and a genomewide association study (GWAS), thereby concentrating on agronomic traits (e.g., plant height, ear length, thousand-seed weight, and seed yield) and oil traits (including fatty acid composition and crude fat content). The GO analysis unveiled a multitude of genes with SNP variations associated with cellular processes, intracellular anatomical structures, and organic cyclic compound binding. Through GWAS, we detected 320 significant SNPs linked to both agronomic (104 SNPs) and oil traits (216 SNPs). Notably, two novel candidate genes, Bna.A05p02350D (SFGH) and Bna.C02p22490D (MDN1), are implicated in thousand-seed weight regulation. Additionally, Bna.C03p14350D (EXO70) and Bna.A09p05630D (PI4Kα1) emerged as novel candidate genes associated with erucic acid and crude fat content, respectively. These findings carry implications for identifying superior genotypes for the development of new cultivars. Association studies offer a cost-effective means of screening mutants and selecting elite rapeseed breeding lines, thereby enhancing the commercial viability of this pivotal oil crop.

Characterization of Quantitative Trait Loci for Leaf Rust Resistance from CI 13227 in Three Winter Wheat Populations.

Leaf rust is a widespread foliar wheat disease causing substantial yield losses worldwide. Slow rusting is "adult plant" resistance that significantly slows epidemic development and thereby reduces yield loss. Wheat accession CI 13227 was previously characterized as having slow-rusting resistance. To validate the quantitative trait loci (QTLs) and develop diagnostic markers for slow rusting resistance in CI 13227, a new population of recombinant inbred lines of CI 13227 × Everest was evaluated for latent period, final severity, area under the disease progress curve, and infection type in greenhouses and genotyped using genotyping-by-sequencing. Four QTLs were identified on chromosome arms 2BL, 2DS, 3BS, and 7BL, explaining 6.82 to 28.45% of the phenotypic variance for these traits. Seven kompetitive allele-specific polymorphism markers previously reported to be linked to the QTLs in two other CI 13227 populations were validated. In addition, the previously reported QLr.hwwg-7AL was remapped to 2BL (renamed QLr.hwwg-2BL) after adding new markers in this study. Phenotypic data showed that the recombinant inbred lines harboring two or three of the QTLs had a significantly longer latent period. QLr.hwwg-2DS on 2DS showed a major effect on all rust resistance traits and was finely mapped to a 2.7-Mb interval by two newly developed flanking markers from exome capture. Three disease-resistance genes and two transporter genes were identified as the putative candidates for QLr.hwwg-2DS. The validated QTLs can be used as slow-rusting resistance resources, and the markers developed in this study will be useful for marker-assisted selection.

Relatedness of Luther Burbank’s Plum (Prunus sp.) Introductions Based on Genotyping by Sequencing

The renowned horticultural artist and plant breeder Luther Burbank worked with many species of plants. During his 50-year career, he introduced more than 800 cultivars, including more than 150 accessions of plums (Prunus spp.) in the late 1800s and early 1900s. Burbank preferred using wide, interspecific crosses to create a vast range of phenotypic variation and then artificially select from the extremes. Although a great artist, Burbank was a substandard scientist because he was derelict in pedigree note-taking. Although many of his introductions are extinct, hobbyists, enthusiasts, and international collections retain nearly a third of the economically viable cultivars he bred. For a century, many of his hybridizations remained inscrutable mysteries until modern genomic and computational tools developed their resolution and statistical power. Today, genotyping by sequencing (GBS) is a useful tool for pedigree reconstruction in the absence of reliable records. GBS can inform principal component analyses, identity by descent (IBD) kinship, and phylogenetic admixture, revealing complex relationships among taxa. In this study, whole genome sequencing was performed on 53 Prunus taxa used by Burbank in his breeding experiments in the most comprehensive genetic survey of his work to date. Exact parent–offspring relationships between this population may be impossible to discern due to years of back crossing, sibling mating, and open pollination. However, the proportion of genomic similarity among these taxa provides information on the relatedness of the genotypes in Burbank’s Prunus experiments, defining four primary lineages within his breeding population. These lineages comprised primarily P. salicina and P. simonii but also have influences from P. americana, P. cerasifera, P. domestica, and P. rivularis. The prevalence of P. simonii in Burbank’s Prunus introductions appears to have been vastly underreported, indicating that some of the seedstock founders of his breeding population could have been P. salicina × P. simonii hybrids at the inception of his career. This research has implications for pedigree reconstruction and prioritizing conservation in collections curation for future studies.

Analysis of population structure and genetic diversity of Camellia tachangensis in Guizhou based on SNP markers

BackgroundCamellia tachangensis F. C. Zhang is a five-compartment species in the ovary of tea group plants, which represents the original germline of early differentiation of some tea group plants.Methods and resultsIn this study, we analyzed single-nucleotide polymorphisms (SNPs) at the genome level, constructed a phylogenetic tree, analyzed the genetic diversity, and further investigated the population structure of 100 C. tachangensis accessions using the genotyping-by-sequencing (GBS) method. A total of 91,959 high-quality SNPs were obtained. Population structure analysis showed that the 100 C. tachangensis accessions clustered into three groups: YQ-1 (Village Group), YQ-2 (Forest Group) and YQ-3 (Transition Group), which was further consistent with the results of phylogenetic analysis and principal component analyses (PCA). In addition, a comparative analysis of the genetic diversity among the three populations (Forest, Village, and Transition Groups) detected the highest genetic diversity in the Transition Group and the highest differentiation between Forest and Village Groups.ConclusionsC. tachangensis plants growing in the forest had different genetic backgrounds from those growing in villages. This study provides a basis for the effective protection and utilization of C. tachangensis populations and lays a foundation for future C. tachangensis breeding.

Genetic diversity and population structure of a rare flowering tree endemic to Appalachia, Stewartia ovata.

Stewartia ovata (cav.) Weatherby, commonly known as mountain stewartia, is an understory tree native to the southeastern United States (U.S.). This relatively rare species occurs in isolated populations in Virginia, Kentucky, Tennessee, North Carolina, South Carolina, Georgia, Alabama, and Mississippi. As a species, S. ovata has largely been overlooked, and limited information is available regarding its ecology, which presents obstacles to conservation efforts. Stewartia ovata has vibrant, large white flowers that bloom in summer with a variety of filament colors, suggesting potential horticultural traits prized by ornamental industry. However, S. ovata is relatively slow growing and, due to long seed dormancy, propagation is challenging with limited success rates. This has created a need to assess the present genetic diversity in S. ovata populations to inform potential conservation and restoration of the species. Here, we employ a genotyping-by-sequencing (GBS) approach to characterize the spatial distribution and genetic diversity of S. ovata in the southern Appalachia region of the eastern United States. A total of 4475 single nucleotide polymorphisms (SNPs) were identified across 147 individuals from 11 collection sites. Our results indicate low genetic diversity (He = 0.216), the presence of population structure (K = 2), limited differentiation (F ST = 0.039), and high gene flow (Nm = 6.16) between our subpopulations. Principal component analysis corroborated the findings of STRUCTURE, confirming the presence of two distinct S. ovata subpopulations. One subpopulation mainly contains genotypes from the Cumberland Plateau, Tennessee, while the other consists of genotypes present in the Great Smoky Mountain ranges in Tennessee, North Carolina, and portions of Nantahala, Chattahoochee-Oconee national forests in Georgia, highlighting that elevation likely plays a major role in its distribution. Our results further suggested low inbreeding coefficient (F IS = 0.070), which is expected with an outcrossing tree species. This research further provides necessary insight into extant subpopulations and has generated valuable resources needed for conservation efforts of S. ovata.

Adaptive divergence, historical population dynamics, and simulation of suitable distributions for Picea Meyeri and P. Mongolica at the whole-genome level

The taxonomic classification of Picea meyeri and P. mongolica has long been controversial. To investigate the genetic relatedness, evolutionary history, and population history dynamics of these species, genotyping-by-sequencing (GBS) technology was utilized to acquire whole-genome single nucleotide polymorphism (SNP) markers, which were subsequently used to assess population structure, population dynamics, and adaptive differentiation. Phylogenetic and population structural analyses at the genomic level indicated that although the ancestor of P. mongolica was a hybrid of P. meyeri and P. koraiensis, P. mongolica is an independent Picea species. Additionally, P. mongolica is more closely related to P. meyeri than to P. koraiensis, which is consistent with its geographic distribution. There were up to eight instances of interspecific and intraspecific gene flow between P. meyeri and P. mongolica. The P. meyeri and P. mongolica effective population sizes generally decreased, and Maxent modeling revealed that from the Last Glacial Maximum (LGM) to the present, their habitat areas decreased initially and then increased. However, under future climate scenarios, the habitat areas of both species were projected to decrease, especially under high-emission scenarios, which would place P. mongolica at risk of extinction and in urgent need of protection. Local adaptation has promoted differentiation between P. meyeri and P. mongolica. Genotype‒environment association analysis revealed 96,543 SNPs associated with environmental factors, mainly related to plant adaptations to moisture and temperature. Selective sweeps revealed that the selected genes among P. meyeri, P. mongolica and P. koraiensis are primarily associated in vascular plants with flowering, fruit development, and stress resistance. This research enhances our understanding of Picea species classification and provides a basis for future genetic improvement and species conservation efforts.

Characterization of Quantitative Trait Loci for Leaf Rust Resistance in the Uzbekistani Wheat Landrace Teremai Bugdai.

Leaf rust, caused by Puccinia triticina, is a major cause of wheat yield losses globally, and novel leaf rust resistance genes are needed to enhance wheat leaf rust resistance. Teremai Bugdai is a landrace from Uzebekistan that is highly resistant to many races of P. triticina in the United States. To unravel leaf rust resistance loci in Teremai Bugdai, a recombinant inbred line (RIL) population of Teremai Bugdai × TAM 110 was evaluated for response to P. triticina race Pt54-1 (TNBGJ) and genotyped using single nucleotide polymorphism (SNP) markers generated by genotyping-by-sequencing (GBS). Quantitative trait loci (QTL) analysis using 5,130 high-quality GBS-SNPs revealed three QTLs, QLr-Stars-2DS, QLr-Stars-6BL, and QLr.Stars-7BL, for leaf rust resistance in two experiments. QLr-Stars-2DS, which is either a new Lr2 allele or a new resistance locus, was delimited to an ∼19.47-Mb interval between 46.4 and 65.9 Mb on 2DS and explained 31.3 and 33.2% of the phenotypic variance in the two experiments. QLr-Stars-6BL was mapped in an ∼84.0-kb interval between 719.48 and 719.56 Mb on 6BL, accounting for 33 to 36.8% of the phenotypic variance in two experiments. QLr.Stars-7BL was placed in a 350-kb interval between 762.41 and 762.76 Mb on 7BL and explained 4.4 to 5.3% of the phenotypic variance. Nine GBS-SNPs flanking these QTLs were converted to kompetitive allele specific PCR (KASP) markers, and these markers can be used to facilitate their introgression into locally adapted wheat lines.

Biogeography of larches in eastern Siberia – using single nucleotide polymorphisms derived by genotyping by sequencing

The present distribution of Siberian boreal forests that are dominated by larches (Larix spp.) is influenced, to an unknown extent, by glacial history. Knowing the past treeline dynamics can improve our understanding of future treeline shifts under changing climate. Here, we study patterns in the genetic variability of Siberian Larix to help unravel biogeographic migration routes since the Last Glacial Maximum (LGM).We infer the spatial distribution and the postglacial demographic history of Larix using genome‐wide single nucleotide polymorphisms (SNPs) derived through genotyping by sequencing (GBS) from 130 individuals sampled across eastern Siberia.Our analysis gives statistical support for two or three clusters, spanning from western to eastern Siberia. These clusters reveal a genetic structure influenced by isolation resulting from geographical distance, barriers imposed by geographic features, and distinct glacial histories. Assuming three clusters, our demographic inference indicates that the common ancestor of the current Larix populations existed in northeast Siberia well before the LGM. This suggests that Larix persisted in the northern region throughout previous glacials.Our genetic studies suggest that Larix likely survived the cold LGM in northern refugia, enabling a fast colonization of Siberia. Instead of complete repopulation from southern areas postglacially, the northernmost Larix expansion during the Holocene seems to have benefitted from refugial populations ahead of the treeline. Present‐day migration is expected to be slow initially, due to the absence of current refugial populations in the far north, in contrast to the early‐Holocene situation.

Genotyping-by-sequencing and weighted gene co-expression network analysis of genes responsive against Potato virus Y in commercial potato cultivars.

Potato is considered a key component of the global food system and plays a vital role in strengthening world food security. A major constraint to potato production worldwide is the Potato Virus Y (PVY), belonging to the genus Potyvirus in the family of Potyviridae. Selective breeding of potato with resistance to PVY pathogens remains the best method to limit the impact of viral infections. Understanding the genetic diversity and population structure of potato germplasm is important for breeders to improve new cultivars for the sustainable use of genetic materials in potato breeding to PVY pathogens. While, genetic diversity improvement in modern potato breeding is facing increasingly narrow genetic basis and the decline of the genetic diversity. In this research, we performed genotyping-by-sequencing (GBS)-based diversity analysis on 10 commercial potato cultivars and weighted gene co-expression network analysis (WGCNA) to identify candidate genes related to PVY-resistance. WGCNA is a system biology technique that uses the WGCNA R software package to describe the correlation patterns between genes in multiple samples. In terms of consumption, these cultivars are a high rate among Iranian people. Using population structure analysis, the 10 cultivars were clustered into three groups based on the 118343 single nucleotide polymorphisms (SNPs) generated by GBS. Read depth ranged between 5 and 18. The average data size and Q30 of the reads were 145.98 Mb and 93.63%, respectively. Based on the WGCNA and gene expression analysis, the StDUF538, StGTF3C5, and StTMEM161A genes were associated with PVY resistance in the potato genome. Further, these three hub genes were significantly involved in defense mechanism where the StTMEM161A was involved in the regulation of alkalization apoplast, the StDUF538 was activated in the chloroplast degradation program, and the StGTF3C5 regulated the proteins increase related to defense in the PVY infected cells. In addition, in the genetic improvement programs, these hub genes can be used as genetic markers for screening commercial cultivars for PVY resistance. Our survey demonstrated that the combination of GBS-based genetic diversity germplasm analysis and WGCNA can assist breeders to select cultivars resistant to PVY as well as help design proper crossing schemes in potato breeding.

Multiple Pleistocene refugia for Arctic Bell‐Heather revealed with genomic analyses of modern and historic plants

AbstractAimArctic plants survived the Pleistocene glaciations in unglaciated refugia. The number, ages, and locations of these refugia are often unclear. We use high‐resolution genomic data from present‐day and Little‐Ice‐Age populations of Arctic Bell‐Heather to re‐evaluate the biogeography of this species and determine whether it had multiple independent refugia or a single refugium in Beringia.LocationCircumpolar Arctic and Coastal British Columbia (BC) alpine.Taxon Cassiope tetragona L., subspecies saximontana and tetragona, outgroup C. mertensiana (Ericaceae).MethodsWe built genotyping‐by‐sequencing (GBS) libraries using Cassiope tetragona tissue from 36 Arctic locations, including two ~250‐ to 500‐year‐old populations collected under glacial ice on Ellesmere Island, Canada. We assembled a de novo GBS reference to call variants. Population structure, genetic diversity and demography were inferred from PCA, ADMIXTURE, fastsimcoal2, SplitsTree, and several population genomics statistics.ResultsPopulation structure analyses identified 4–5 clusters that align with geographic locations. Nucleotide diversity was highest in Beringia and decreased eastwards across Canada. Demographic coalescent analyses dated the following splits with Alaska: BC subspecies saximontana (5 mya), Russia (~1.4 mya), Europe (>200–600 kya), and Greenland (~60 kya). Northern Canada populations appear to have formed during the current interglacial (7–9 kya). Admixture analyses show genetic variants from Alaska appear more frequently in present‐day than historic plants on Ellesmere Island.ConclusionsPopulation and demographic analyses support BC, Alaska, Russia, Europe and Greenland as all having had independent Pleistocene refugia. Northern Canadian populations appear to be founded during the current interglacial with genetic contributions from Alaska, Europe and Greenland. We found evidence, on Ellesmere Island, for continued recent gene flow in the last 250–500 years. These results suggest that a re‐analysis of other Arctic species with shallow population structure using higher resolution genomic markers and demographic analyses may help reveal deeper structure and other circumpolar glacial refugia.

Insights into the genetic architecture of Phytophthora capsici root rot resistance in chile pepper (Capsicum spp.) from multi-locus genome-wide association study

BackgroundPhytophthora root rot, a major constraint in chile pepper production worldwide, is caused by the soil-borne oomycete, Phytophthora capsici. This study aimed to detect significant regions in the Capsicum genome linked to Phytophthora root rot resistance using a panel consisting of 157 Capsicum spp. genotypes. Multi-locus genome wide association study (GWAS) was conducted using single nucleotide polymorphism (SNP) markers derived from genotyping-by-sequencing (GBS). Individual plants were separately inoculated with P. capsici isolates, ‘PWB-185’, ‘PWB-186’, and ‘6347’, at the 4–8 leaf stage and were scored for disease symptoms up to 14-days post-inoculation. Disease scores were used to calculate disease parameters including disease severity index percentage, percent of resistant plants, area under disease progress curve, and estimated marginal means for each genotype.ResultsMost of the genotypes displayed root rot symptoms, whereas five accessions were completely resistant to all the isolates and displayed no symptoms of infection. A total of 55,117 SNP markers derived from GBS were used to perform multi-locus GWAS which identified 330 significant SNP markers associated with disease resistance. Of these, 56 SNP markers distributed across all the 12 chromosomes were common across the isolates, indicating association with more durable resistance. Candidate genes including nucleotide-binding site leucine-rich repeat (NBS-LRR), systemic acquired resistance (SAR8.2), and receptor-like kinase (RLKs), were identified within 0.5 Mb of the associated markers.ConclusionsResults will be used to improve resistance to Phytophthora root rot in chile pepper by the development of Kompetitive allele-specific markers (KASP®) for marker validation, genomewide selection, and marker-assisted breeding.

Genetic tapestry of Capsicum fruit colors: a comparative analysis of four cultivated species.

Genome-wide association study of color spaces across the four cultivated Capsicum spp. revealed a shared set of genes influencing fruit color, suggesting mechanisms and pathways across Capsicum species are conserved during the speciation. Notably, Cytochrome P450 of the carotenoid pathway, MYB transcription factor, and pentatricopeptide repeat-containing protein are the major genes responsible for fruit color variation across the Capsicum species. Peppers (Capsicum spp.) rank among the most widely consumed spices globally. Fruit color, serving as a determinant for use in food colorants and cosmeceuticals and an indicator of nutritional contents, significantly influences market quality and price. Cultivated Capsicum species display extensive phenotypic diversity, especially in fruit coloration. Our study leveraged the genetic variance within four Capsicum species (Capsicum baccatum, Capsicum chinense, Capsicum frutescens, and Capsicum annuum) to elucidate the genetic mechanisms driving color variation in peppers and related Solanaceae species. We analyzed color metrics and chromatic attributes (Red, Green, Blue, L*, a*, b*, Luminosity, Hue, and Chroma) on samples cultivated over six years (2015-2021). We resolved genomic regions associated with fruit color diversity through the sets of SNPs obtained from Genotyping by Sequencing (GBS) and genome-wide association study (GWAS) with a Multi-Locus Mixed Linear Model (MLMM). Significant SNPs with FDR correction were identified, within the Cytochrome P450, MYB-related genes, Pentatricopeptide repeat proteins, and ABC transporter family were the most common among the four species, indicating comparative evolution of fruit colors. We further validated the role of a pentatricopeptide repeat-containing protein (Chr01:31,205,460) and a cytochrome P450 enzyme (Chr08:45,351,919) via competitive allele-specific PCR (KASP) genotyping. Our findings advance the understanding of the genetic underpinnings of Capsicum fruit coloration, with developed KASP assays holding potential for applications in crop breeding and aligning with consumer preferences. This study provides a cornerstone for future research into exploiting Capsicum's diverse fruit color variation.

Major chromosome rearrangements in intergeneric wheat × rye hybrids in compatible and incompatible crosses detected by GBS read coverage analysis

The presence of incompatibility alleles in primary amphidiploids constitutes a reproductive barrier in newly synthesized wheat-rye hybrids. To overcome this barrier, the genome stabilization process includes large-scale chromosome rearrangements. In incompatible crosses resulting in fertile amphidiploids, the elimination of one of the incompatible alleles Eml-A1 or Eml-R1b can occur already in the somatic tissue of the wheat × rye hybrid embryo. We observed that the interaction of incompatible loci Eml-A1 of wheat and Eml-R1b of rye after overcoming embryo lethality leads to hybrid sterility in primary triticale. During subsequent seed reproductions (R1, R2 or R3) most of the chromosomes of A, B, D and R subgenomes undergo rearrangement or eliminations to increase the fertility of the amphidiploid by natural selection. Genotyping-by-sequencing (GBS) coverage analysis showed that improved fertility is associated with the elimination of entire and partial chromosomes carrying factors that either cause the disruption of plant development in hybrid plants or lead to the restoration of the euploid number of chromosomes (2n = 56) in the absence of one of the incompatible alleles. Highly fertile offspring obtained in compatible and incompatible crosses can be successfully adapted for the production of triticale pre-breeding stocks.

Genetic Diversity and Genome-Wide Association Analysis of the Hulled/Naked Trait in a Barley Collection from Shanghai Agricultural Gene Bank.

Barley is one of the most important cereal crops in the world, and its value as a food is constantly being revealed, so the research into and the use of barley germplasm are very important for global food security. Although a large number of barley germplasm samples have been collected globally, their specific genetic compositions are not well understood, and in many cases their origins are even disputed. In this study, 183 barley germplasm samples from the Shanghai Agricultural Gene Bank were genotyped using genotyping-by-sequencing (GBS) technology, SNPs were identified and their genetic parameters were estimated, principal component analysis (PCA) was preformed, and the phylogenetic tree and population structure of the samples were also analyzed. In addition, a genome-wide association study (GWAS) was carried out for the hulled/naked grain trait, and a KASP marker was developed using an associated SNP. The results showed that a total of 181,906 SNPs were identified, and these barley germplasm samples could be roughly divided into three categories according to the phylogenetic analysis, which was generally consistent with the classification of the traits of row type and hulled/naked grain. Population structure analysis showed that the whole barley population could be divided into four sub-populations (SPs), the main difference from previous classifications being that the two-rowed and the hulled genotypes were sub-divided into two SPs. The GWAS analysis of the hulled/naked trait showed that many associated loci were unrelated to the Nud/nud locus, indicating that there might be new loci controlling the trait. A KASP marker was developed for one exon-type SNP on chromosome 7. Genotyping based on the KASP assay was consistent with that based on SNPs, indicating that the gene of this locus might be associated with the hulled/naked trait. The above work not only lays a good foundation for the future utilization of this barley germplasm population but it provides new loci and candidate genes for the hulled/naked trait.

Cutting-Edge of Genotyping by Sequencing (GBS) for Improving Drought and Salinity Stress Tolerance in Watermelon (Citrullus lanatus L.): A Review

AbstractGenotyping by sequencing (GBS) is a very helpful approach and one of the most useful techniques for examining and analysing the genetic variation of various lines and varieties. GBS technology was used to identify single nucleotide polymorphisms (SNPs) and assess genetic variation in several watermelon accessions. GBS application in watermelon breeding programs has recently become a popular technique among many breeders. Watermelon (Citrullus lanatusL.) is a warm-season crop that is widely cultivated for its delicious fruit. And it is one of the most economically significant crops in the world. However, watermelon cultivation is frequently hampered by abiotic stressors such as drought and salinity. Recently, there has been a growing body of research on the mechanisms that allow watermelon to tolerate these stresses and improve crop yield. Generally, cucurbits are beneficial to human health, they provide necessary minerals, fibre, and nutrient components. Therefore, this review demonstrates the cutting edge of using GBS technology to identify genetic design of several features in watermelon to improve abiotic stresses (drought and saline). The application of the GBS technique has provided a distinct advantage in watermelon breeding studies. Based on GBS approach, many new candidate genes in watermelon lines control a variety of traits including saline and drought tolerance, fruit rind color, disease tolerance, nutrient components, size, and fruit shape were discovered. Modern breeding techniques are being used to develop economically viable vegetable crops that will meet customer preferences and needs. Further research is needed to enhance watermelon production.

Genetic Diversity and Genome-Wide Association in Cowpeas (Vigna unguiculata L. Walp)

Cowpea is one of the most popular dry-land legumes cultivated for food and forage in arid and semi-arid areas. Genetic diversity for global germplasm can be organized into core collections providing optimum resources to serve breeding requirements. Here, we present diversity analysis and genome-wide association study (GWAS) results for part of the cowpea core collection of the United States Department of Agriculture (USDA) along with breeding line controls. Included in the analysis were a total of 373 accessions analyzed with 6880 Single Nucleotide Polymorphism (SNP) markers from Genotyping by Sequencing (GBS). Population structure differentiated accessions into two groups irrespective of geographical origin and formed three clusters based on taxa upon phylogenetic analysis. A total of 56 SNPs were significantly associated to nine traits including pod length (25 Quantitative Trait Nucleotides, QTNs), seed anti-oxidant content (7 QTNs), dry pod color (7 QTNs), plant maturity (5 QTNs), flower color (5 QTNs), seed weight (4 QTNs), tolerance to low phosphate (1 QTN), growth habit (1 QTN), and response to rock phosphate (1 QTN) using Bayesian-information, Linkage-disequilibrium Iteratively Nested Keyway (BLINK), and Fixed and random model Circulating Probability Unification (FarmCPU) association models. Key genes related to all significant SNPs were identified based on annotations of the cowpea reference genome, including a flavonoid gene controlling flower color (Vigun08g040200.1), a root nodulation regulator for tolerance to low phosphate (Vigun11g168000.1), and numerous genes involved in signaling, biosynthesis, metabolite transport, and abiotic stress. Our results highlight the importance of maintaining public phenotyping databases at USDA and strengthening collaborations for data collection in cowpea to maximize research impacts.

A genomic investigation on the origins of the Korean brown planthopper, Nilaparvata lugens (Stål) (Hemiptera: Delphasidae)

AbstractNilaparvata lugens, the brown planthopper (BPH), is a serious pest species. BPHs cause significant damage to rice plants in Korea as well as other countries in East and Southeast Asia. As BPHs cannot survive winter in Korea, they annually migrate into Korea from foreign countries. The BPHs found in Korea are believed to originate from China, but most BPHs in China are also known to originate from Southeast Asia. To understand the origin of Korean BPHs, we conducted a genotyping‐by‐sequencing (GBS) study. We sampled BPHs from five locations in Korea and five countries in Southeast Asia and analyzed the GBS and sequencing results using various methods based on the f statistics and admixture graph analyses. We confirmed that the domestic BPHs shared a greater genetic drift than the BPHs from Southeast Asian populations over several years, implying that a continuous genetic substratum of Korea exists. This genetic substratum is genetically closer to BPHs from the southern part (Thailand and Cambodia) of Southeast Asia than to BPHs from the northern part (Bhutan, Myanmar and Laos) of Southeast Asia. In addition, direct gene flows from Southeast Asia seem possible, so Korea is considered a hotspot where diverse BPH populations mix. Therefore, the origin of Korean BPHs extends beyond China and as far as southern Southeast Asia. This result will help to understand and control the population dynamics of the Korean BPH population.