SNP2CAPS: a SNP and INDEL analysis tool for CAPS marker development

Molecular selection of soybean towards adaptation to Central European agroclimatic conditions

Europe is highly dependent on soybean meal imports and anticipates an increase of domestic plant protein production. Ongoing climate change resulted in northward shift of plant hardiness zones, enabling spring-sowing of freezing-sensitive crops, including soybean. However, it requires efficient reselection of germplasm adapted to relatively short growing season and long-day photoperiod. In the present study, a PCR array has been implemented, targeting early maturity (E1-E4, E7, E9, and E10), pod shattering (qPHD1), and growth determination (Dt1) genes. This array was optimized for routine screening of soybean diversity panel (204 accessions), subjected to the 2018-2020 survey of phenology, morphology, and yield-related traits in a potential cultivation region in Poland. High broad-sense heritability (0.84-0.88) was observed for plant height, thousand grain weight, maturity date, and the first pod height. Significant positive correlations were identified between the number of seeds and pods per plant, between these two traits and seed yield per plant as well as between flowering, maturity, plant height, and first pod height. PCR array genotyping revealed high genetic diversity, yielding 98 allelic combinations. The most remarkable correlations were identified between flowering and E7 or E1, between maturity and E4 or E7 and between plant height and Dt1 or E4. The study demonstrated high applicability of this PCR array for molecular selection of soybean towards adaptation to Central Europe, designating recessive qPHD1 and dominant Dt1, E3, and E4 alleles as major targets to align soybean growth season requirements with the length of the frost-free period, improve plant performance, and increase yield.

Computational tools for plant genomics and breeding

Plant genomics and crop breeding are at the intersection of biotechnology and information technology. Driven by a combination of high-throughput sequencing, molecular biology and data science, great advances have been made in omics technologies at every step along the central dogma, especially in genome assembling, genome annotation, epigenomic profiling, and transcriptome profiling. These advances further revolutionized three directions of development. One is genetic dissection of complex traits in crops, along with genomic prediction and selection. The second is comparative genomics and evolution, which open up new opportunities to depict the evolutionary constraints of biological sequences for deleterious variant discovery. The third direction is the development of deep learning approaches for the rational design of biological sequences, especially proteins, for synthetic biology. All three directions of development serve as the foundation for a new era of crop breeding where agronomic traits are enhanced by genome design.

Identification of Chromosomal Regions and Candidate Genes for Round leaf Locus in Cucumis melo L

Leaf morphology plays a crucial role in plant classification and provides a significant model for studying plant diversity while directly impacting photosynthetic efficiency. In the case of melons, leaf shape not only influences production and classification but also represents a key genetic trait that requires further exploration. In this study, we utilized forward genetics to pinpoint a recessive locus, dubbed Cmrl (Round leaf), which is responsible for regulating melon leaf shape. Through bulked segregant analysis sequencing and extensive evaluation of a two-year F2 population, we successfully mapped the Cmrl locus to a 537.07 kb region on chromosome 8 of the melon genome. Subsequent genetic fine-mapping efforts, leveraging a larger F2 population encompassing 1322 plants and incorporating F2:3 phenotypic data, further refined the locus to an 80.27 kb interval housing five candidate genes. Promoter analysis and coding sequence cloning confirmed that one of these candidates, MELO3C019152.2 (Cmppr encoding a pentatricopeptide repeat-containing family protein, Cmppr), stands out as a strong candidate gene for the Cmrl locus. Notably, comparisons of Cmrl expressions across various stages of leaf development and different leaf regions suggest a pivotal role of Cmrl in the morphogenesis of melon leaves.

High-resolution mapping of Ryd4Hb, a major resistance gene to Barley yellow dwarf virus from Hordeum bulbosum

KEY MESSAGE

We mapped Ryd4Hb in a 66.5 kbp interval in barley and dissociated it from a sublethality factor. These results will enable a targeted selection of the resistance in barley breeding. Virus diseases are causing high yield losses in crops worldwide. The Barley yellow dwarf virus (BYDV) complex is responsible for one of the most widespread and economically important viral diseases of cereals. While no gene conferring complete resistance (immunity) has been uncovered in the primary gene pool of barley, sources of resistance were searched and identified in the wild relative Hordeum bulbosum, representing the secondary gene pool of barley. One such locus, Ryd4Hb, has been previously introgressed into barley, and was allocated to chromosome 3H, but is tightly linked to a sublethality factor that prevents the incorporation and utilization of Ryd4Hb in barley varieties. To solve this problem, we fine-mapped Ryd4Hb and separated it from this negative factor. We narrowed the Ryd4Hb locus to a corresponding 66.5 kbp physical interval in the barley 'Morex' reference genome. The region comprises a gene from the nucleotide-binding and leucine-rich repeat immune receptor family, typical of dominant virus resistance genes. The closest homolog to this Ryd4Hb candidate gene is the wheat Sr35 stem rust resistance gene. In addition to the fine mapping, we reduced the interval bearing the sublethality factor to 600 kbp in barley. Aphid feeding experiments demonstrated that Ryd4Hb provides a resistance to BYDV rather than to its vector. The presented results, including the high-throughput molecular markers, will permit a more targeted selection of the resistance in breeding, enabling the use of Ryd4Hb in barley varieties.

The Resistance of Soybean Variety Heinong 84 to Apple Latent Spherical Virus Is Controlled by Two Genetic Loci

Apple latent spherical virus (ALSV) is widely used as a virus-induced gene silencing (VIGS) vector for function genome study. However, the application of ALSV to soybeans is limited by the resistance of many varieties. In this study, the genetic locus linked to the resistance of a resistant soybean variety Heinong 84 was mapped by high-throughput sequencing-based bulk segregation analysis (HTS-BSA) using a hybrid population crossed from Heinong 84 and a susceptible variety, Zhonghuang 13. The results showed that the resistance of Heinong 84 to ALSV is controlled by two genetic loci located on chromosomes 2 and 11, respectively. Cleaved amplified polymorphic sequence (CAPS) markers were developed for identification and genotyping. Inheritance and biochemical analyses suggest that the resistance locus on chromosome 2 plays a dominant dose-dependent role, while the other locus contributes a secondary role in resisting ALSV. The resistance locus on chromosome 2 might encode a protein that can directly inhibit viral proliferation, while the secondary resistance locus on chromosome 11 may encode a host factor required for viral proliferation. Together, these data reveal novel insights on the resistance mechanism of Heinong 84 to ALSV, which will benefit the application of ALSV as a VIGS vector.

Molecular Mapping of Putative Genomic Regions Controlling Fruit and Seed Morphology of Watermelon

The genetic regulatory basis of qualitative and quantitative phenotypes of watermelon is being investigated in different types of molecular and genetic breeding studies around the world. In this study, biparental F2 mapping populations were developed over two experimental years, and the collected datasets of fruit and seed traits exhibited highly significant correlations. Whole-genome resequencing of comparative parental lines was performed and detected single nucleotide polymorphism (SNP) loci were converted into cleaved amplified polymorphic sequence (CAPS) markers. The screened polymorphic markers were genotyped in segregating populations and two genetic linkage maps were constructed, which covered a total of 2834.28 and 2721.45 centimorgan (cM) genetic lengths, respectively. A total of 22 quantitative trait loci (QTLs) for seven phenotypic traits were mapped; among them, five stable and major-effect QTLs (PC-8-1, SL-9-1, SWi-9-1, SSi-9-1, and SW-6-1) and four minor-effect QTLs (PC-2-1 and PC-2-2; PT-2-1 and PT-2-2; SL-6-1 and SSi-6-2; and SWi-6-1 and SWi-6-2) were observed with 3.77-38.98% PVE. The adjacent QTL markers showed a good fit marker-trait association, and a significant allele-specific contribution was also noticed for genetic inheritance of traits. Further, a total of four candidate genes (Cla97C09G179150, Cla97C09G179350, Cla97C09G180040, and Cla97C09G180100) were spotted in the stable colocalized QTLs of seed size linked traits (SL-9-1 and SWi-9-1) that showed non-synonymous type mutations. The gene expression trends indicated that the seed morphology had been formed in the early developmental stage and showed the genetic regulation of seed shape formation. Hence, we think that our identified QTLs and genes would provide powerful genetic insights for marker-assisted breeding aimed at improving the quality traits of watermelon.

Clorf Encodes Carotenoid Isomerase and Regulates Orange Flesh Color in Watermelon (Citrullus lanatus L.)

Flesh color is a significant characteristic of watermelon. Although various flesh-color genes have been identified, the inheritance and molecular basis of the orange flesh trait remain relatively unexplored. In the present study, the genetic analysis of six generations derived from W1-1 (red flesh) and W1-61 (orange flesh) revealed that the orange flesh color trait was regulated by a single recessive gene, Clorf (orange flesh). Bulk segregant analysis (BSA) locked the range to ∼4.66 Mb, and initial mapping situated the Clorf locus within a 688.35-kb region of watermelon chromosome 10. Another 1,026 F2 plants narrowed the Clorf locus to a 304.62-kb region containing 32 candidate genes. Subsequently, genome sequence variations in this 304.62-kb region were extracted for in silico BSA strategy among 11 resequenced lines (one orange flesh and ten nonorange flesh) and finally narrowed the Clorf locus into an 82.51-kb region containing nine candidate genes. Sequence variation analysis of coding regions and gene expression levels supports Cla97C10G200950 as the most possible candidate for Clorf, which encodes carotenoid isomerase (Crtiso). This study provides a genetic resource for investigating the orange flesh color of watermelon, with Clorf malfunction resulting in low lycopene accumulation and, thus, orange flesh.

Genetic Mapping and QTL Analysis of Fruit Traits in Melon (Cucumis melo L.)

Melon (Cucumis melo L.) is an important horticultural cash crop and its quality traits directly affect consumer choice and market price. These traits are controlled by genetic as well as environmental factors. In this study, a quantitative trait locus (QTL) mapping strategy was used to identify the potential genetic loci controlling quality traits of melons (i.e., exocarp and pericarp firmness and soluble solid content) based on newly derived whole-genome single nucleotide polymorphism-based cleaved amplified polymorphic sequence (SNP-CAPS) markers. Specifically, SNPs of two melon varieties, M4-5 and M1-15, as revealed by whole-genome sequencing, were converted to the CAPS markers, which were used to construct a genetic linkage map comprising 12 chromosomes with a total length of 1414.88 cM, in the F₂ population of M4-5 and M1-15. The six identified QTLs included: SSC6.1 and SSC11.1 related to soluble solid content; EF12.1 associated with exocarp firmness; and EPF3.1, EPF3.2 and EPF7.1 related to edible pericarp firmness. These genes were located on five chromosomes (3, 6, 7, 11, and 12) in the flanking regions of the CAPS markers. Moreover, the newly developed CAPS markers will be useful in guiding genetic engineering and molecular breeding in melon.

Integrated analysis of biparental and natural populations reveals CRIB domain-containing protein underlying seed coat crack trait in watermelon

The scc locus of the watermelon seed coat crack trait was fine mapped on chromosome 3. Cla97C03G056110 (annotated as CRIB domain-containing protein) was regarded as the most likely candidate gene Seed coat crack (scc) is a special characteristic of watermelon compared with other cucurbit crops. However, information regarding the genetic basis of this trait is limited. We conducted a genetic analysis of six generations derived from PI 192938 (scc) and Cream of Saskatchewan (COS) (non-scc) parental lines and found that the scc trait was regulated by a single recessive gene through two years. Bulk segregant analysis sequencing (BSA-seq) and initial mapping placed the scc locus into an 808.8 kb region on chromosome 3. Evaluation of another 1152 F₂ plants narrowed the scc locus to a 277.11 kb region containing 37 candidate genes. Due to the lack of molecular markers in the fine-mapping interval, we extracted the genome sequence variations in this 277.11 kb region with in silico BSA among seventeen re-sequenced lines (6 scc and 11 non-scc) and finally delimited the scc locus to an 8.34 kb region with only one candidate gene Cla97C03G056110 (CRIB domain-containing protein). Three single nucleotide polymorphism loci in the promoter region of Cla97C03G056110 altered cis-acting elements that were highly correlated with the nature watermelon panel. The expression of Cla97C03G056110 in seed coat tissue was higher in non-scc than in scc lines and was specifically expressed in seed coat compared with fruit flesh.

Identification of chromosome region and candidate genes for canary-yellow flesh (Cyf) locus in watermelon (Citrullus lanatus)

Genetic control of fruit flesh color in watermelon is complex, and significant knowledge gaps still exist. In the present study, we investigated the genetic basis of canary-yellow flesh color in watermelon inbred line PI 635597 using a segregating population derived from a cross between PI 635597 and another inbred line, Cream of Saskatchewan (pale yellow flesh color). We showed that a single dominant gene controls the canary-yellow flesh color for the Cyf (canary-yellow flesh) trait. Bulk segregant analysis (BSA) and fine genetic mapping narrowed down the Cyf locus to a 79.62-kb region on chromosome 6, which harbors 10 predicted genes. Sequence variation analysis in the promoter and coding regions and gene expression analysis in both parental lines and selected watermelon accessions with diverse fruit flesh colors support Cla97C06G122050 (unknown protein) and Cla97C06G122120 (pentatricopeptide repeat) as predicted candidate genes for the Cyf locus. Marker-assisted selection and sequence alignment showed that the Cyf locus could differentiate canary-yellow flesh and pale-yellow flesh. Our results indicate that the Cyf locus might be responsible for canary-yellow flesh color and carotenoid accumulation levels.

Primary mapping of quantitative trait loci regulating multivariate horticultural phenotypes of watermelon (Citrullus lanatus L.)

Watermelon fruits exhibit a remarkable diversity of important horticultural phenotypes. In this study, we initiated a primary quantitative trait loci (QTL) mapping to identify the candidate regions controlling the ovary, fruit, and seed phenotypes. Whole genome sequencing (WGS) was carried out for two differentiated watermelon lines, and 350 Mb (96%) and 354 Mb (97%) of re-sequenced reads covered the reference de novo genome assembly, individually. A total of 45.53% non-synonymous single nucleotide polymorphism (nsSNPs) and 54.47% synonymous SNPs (sSNPs) were spotted, which produced 210 sets of novel SNP-based cleaved amplified polymorphism sequence (CAPS) markers by depicting 46.25% co-dominant polymorphism among parent lines and offspring. A biparental F_2:3 mapping population comprised of 100 families was used for trait phenotyping and CAPS genotyping, respectively. The constructed genetic map spanned a total of 2,398.40 centimorgans (cM) in length and averaged 11.42 cM, with 95.99% genome collinearity. A total of 33 QTLs were identified at different genetic positions across the eight chromosomes of watermelon (Chr-01, Chr-02, Chr-04, Chr-05, Chr-06, Chr-07, Chr-10, and Chr-11); among them, eight QTLs of the ovary, sixteen QTLs of the fruit, and nine QTLs of the seed related phenotypes were classified with 5.32-25.99% phenotypic variance explained (PVE). However, twenty-four QTLs were identified as major-effect and nine QTLs were mapped as minor-effect QTLs across the flanking regions of CAPS markers. Some QTLs were exhibited as tightly localized across the nearby genetic regions and explained the pleiotropic effects of multigenic nature. The flanking QTL markers also depicted significant allele specific contributions and accountable genes were predicted for respective traits. Gene Ontology (GO) functional enrichment was categorized in molecular function (MF), cellular components (CC), and biological process (BP); however, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were classified into three main classes of metabolism, genetic information processing, and brite hierarchies. The principal component analysis (PCA) of multivariate phenotypes widely demonstrated the major variability, consistent with the identified QTL regions. In short, we assumed that our identified QTL regions provide valuable genetic insights regarding the watermelon phenotypes and fine genetic mapping could be used to confirm them.

Fine Mapping the Soybean Mosaic Virus Resistance Gene in Soybean Cultivar Heinong 84 and Development of CAPS Markers for Rapid Identification

Heinong 84 is one of the major soybean varieties growing in Northeast China, and is resistant to the infection of all strains of soybean mosaic virus (SMV) in the region including the most prevalent strain, N3. However, the resistance gene(s) in Heinong 84 and the resistant mechanism are still elusive. In this study, genetic and next-generation sequencing (NGS)-based bulk segregation analysis (BSA) were performed to map the resistance gene using a segregation population from the cross of Heinong 84 and a susceptible cultivar to strain N3, Zhonghuang 13. Results show that the resistance of Heinong 84 is controlled by a dominant gene on chromosome 13. Further analyses suggest that the resistance gene in Heinong 84 is probably an allele of Rsv1. Finally, two pairs of single-nucleotide-polymorphism (SNP)-based primers that are tightly cosegregated with the resistance gene were designed for rapidly identifying resistant progenies in breeding via the cleaved amplified polymorphic sequence (CAPS) assay.

Fine mapping of a stripe rust resistance gene YrZM175 in bread wheat

A stripe rust resistance gene YrZM175 in Chinese wheat cultivar Zhongmai 175 was mapped to a genomic interval of 636.4 kb on chromosome arm 2AL, and a candidate gene was predicted. Stripe rust, caused by Puccinia striiformis f. sp. tritici (PST), is a worldwide wheat disease that causes large losses in production. Fine mapping and cloning of resistance genes are important for accurate marker-assisted breeding. Here, we report the fine mapping and candidate gene analysis of stripe rust resistance gene YrZM175 in a Chinese wheat cultivar Zhongmai 175. Fifteen F₁, 7,325 F₂ plants and 117 F_2:3 lines derived from cross Avocet S/Zhongmai 175 were inoculated with PST race CYR32 at the seedling stage in a greenhouse, and F_2:3 lines were also evaluated for stripe rust reaction in the field using mixed PST races. Bulked segregant RNA-seq (BSR-seq) analyses revealed 13 SNPs in the region 762.50-768.52 Mb on chromosome arm 2AL. By genome mining, we identified SNPs and InDels between the parents and contrasting bulks and mapped YrZM175 to a 0.72-cM, 636.4-kb interval spanned by YrZM175-InD1 and YrZM175-InD2 (763,452,916-764,089,317 bp) including two putative disease resistance genes based on IWGSC RefSeq v1.0. Collinearity analysis indicated similar target genomic intervals in Chinese Spring, Aegilops tauschii (2D: 647.7-650.5 Mb), Triticum urartu (2A: 750.7-752.3 Mb), Triticum dicoccoides (2A: 771.0-774.5 Mb), Triticum turgidum (2B: 784.7-788.2 Mb), and Triticum aestivum cv. Aikang 58 (2A: 776.3-778.9 Mb) and Jagger (2A: 789.3-791.7 Mb). Through collinearity analysis, sequence alignments of resistant and susceptible parents and gene expression level analysis, we predicted TRITD2Bv1G264480 from Triticum turgidum to be a candidate gene for map-based cloning of YrZM175. A gene-specific marker for TRITD2Bv1G264480 co-segregated with the resistance gene. Molecular marker analysis and stripe rust response data revealed that YrZM175 was different from genes Yr1, Yr17, Yr32, and YrJ22 located on chromosome 2A. Fine mapping of YrZM175 lays a solid foundation for functional gene analysis and marker-assisted selection for improved stripe rust resistance in wheat.

Fine mapping of a stripe rust resistance gene YrZM175 in bread wheat

A stripe rust resistance gene YrZM175 in Chinese wheat cultivar Zhongmai 175 was mapped to a genomic interval of 636.4 kb on chromosome arm 2AL, and a candidate gene was predicted. Stripe rust, caused by Puccinia striiformis f. sp. tritici (PST), is a worldwide wheat disease that causes large losses in production. Fine mapping and cloning of resistance genes are important for accurate marker-assisted breeding. Here, we report the fine mapping and candidate gene analysis of stripe rust resistance gene YrZM175 in a Chinese wheat cultivar Zhongmai 175. Fifteen F₁, 7,325 F₂ plants and 117 F_2:3 lines derived from cross Avocet S/Zhongmai 175 were inoculated with PST race CYR32 at the seedling stage in a greenhouse, and F_2:3 lines were also evaluated for stripe rust reaction in the field using mixed PST races. Bulked segregant RNA-seq (BSR-seq) analyses revealed 13 SNPs in the region 762.50-768.52 Mb on chromosome arm 2AL. By genome mining, we identified SNPs and InDels between the parents and contrasting bulks and mapped YrZM175 to a 0.72-cM, 636.4-kb interval spanned by YrZM175-InD1 and YrZM175-InD2 (763,452,916-764,089,317 bp) including two putative disease resistance genes based on IWGSC RefSeq v1.0. Collinearity analysis indicated similar target genomic intervals in Chinese Spring, Aegilops tauschii (2D: 647.7-650.5 Mb), Triticum urartu (2A: 750.7-752.3 Mb), Triticum dicoccoides (2A: 771.0-774.5 Mb), Triticum turgidum (2B: 784.7-788.2 Mb), and Triticum aestivum cv. Aikang 58 (2A: 776.3-778.9 Mb) and Jagger (2A: 789.3-791.7 Mb). Through collinearity analysis, sequence alignments of resistant and susceptible parents and gene expression level analysis, we predicted TRITD2Bv1G264480 from Triticum turgidum to be a candidate gene for map-based cloning of YrZM175. A gene-specific marker for TRITD2Bv1G264480 co-segregated with the resistance gene. Molecular marker analysis and stripe rust response data revealed that YrZM175 was different from genes Yr1, Yr17, Yr32, and YrJ22 located on chromosome 2A. Fine mapping of YrZM175 lays a solid foundation for functional gene analysis and marker-assisted selection for improved stripe rust resistance in wheat.

Biparental genetic mapping reveals that CmCLAVATA3 (CmCLV3) is responsible for the variation in carpel number in melon (Cucumis melo L.)

Genetic analysis revealed that CmCLV3 is a candidate gene for the variation in melon carpel number. Carpel number (CN) is an important trait in melon. Three-CN melon fruit is oval, while 5-CN melon fruit has a round or flat shape. Herein, a genetic analysis of a population in which the CN locus was segregated indicated that 3-CN is controlled by a major dominant effective gene. Bulked segregant analysis and initial linkage mapping placed the CN locus in a 6.67 Mb region on chromosome 12, and it was narrowed to 882.19 kb with molecular markers and recombinant plants. Fine mapping with a large F₂ population containing 1026 individuals further narrowed the locus to an 83.98 kb region harboring five annotated genes. Gene structure alignment between the parental lines revealed MELO3C035640.2 (annotated as CLAVATA3, CmCLV3) as the best candidate gene for the CN trait. CmCLV3 was more highly expressed in 3- than 5-CN lines and specifically expressed in terminal buds rather than in young leaves, hypocotyls, and roots. The CmCLV3 coding region was cloned from eight 3- or 5-CN melon accessions, and a nonsynonymous SNP site was highly correlated with CN variation. This SNP site was also related to CN variations among 40 melon lines according to their resequencing data, causing a helix alteration in the CmCLV3 protein. Promoter region sequence alignment and activity analysis showed that, unlike in cucumber and tomato, CmCLV3 promoter variation and activity were not the main reasons for CN alteration. Overall, this study provides a genetic resource for melon fruit development research and molecular breeding tools for melon CN improvement.

Development of Whole Genome SNP-CAPS Markers and Preliminary QTL Mapping of Fruit Pedicel Traits in Watermelon

Fruit pedicel (FP) is an important determinant of premium fruit quality that directly affects commercial market value. However, in-depth molecular and genetic basis of pedicel-related traits has not been identified in watermelon. Herein, a quantitative trait locus (QTL) mapping strategy was used to identify the potential genetic regions controlling FP traits based on newly derived whole-genome single nucleotide polymorphism based cleaved amplified polymorphism sequence (SNP-CAPS) markers. Next-generation sequencing based whole-genome re-sequencing of two watermelon parent lines revealed 98.30 and 98.40% of average coverage, 4,989,869 SNP variants, and 182,949 CAPS loci pairs across the reference genome, respectively. A total of 221 sets of codominant markers exhibited 46.42% polymorphism rate and were effectively genotyped within 100-F_2:3 derived mapping population. The developed linkage map covered a total of 2,630.49 cM genetic length with averaged 11.90 cM, and depicted a valid marker-trait association. In total, 6 QTLs (qFPL4.1, qFPW4.1, qFPD2.1, qFPD2.2, qFPD8.1, qFPD10.1) were mapped with five major effects and one minor effect between the whole genome adjacent markers positioned over distinct chromosomes (02, 04, 08, 10), based on the ICIM-ADD mapping approach. These significant QTLs were similarly mapped in delimited flanking regions of 675.10, 751.38, 859.24, 948.39, and 947.51 kb, which collectively explained 8.64-13.60% PVE, respectively. A highly significant and positive correlation was found among the observed variables. To our knowledge, we first time reported the mapped QTLs/genes affecting FP traits of watermelon, and our illustrated outcomes will deliver the potential insights for fine genetic mapping as well as functional gene analysis through MAS-based breeding approaches.

Development of Novel Markers for Yield in Hevea brasiliensis Muell. Arg. Based on Candidate Genes from Biosynthetic Pathways Associated with Latex Production

Scarcity of functional genetic markers associated with candidate genes (CGs) is a serious constraint for marker-assisted selection in the natural rubber producing tree, Hevea brasiliensis. In order to develop markers associated with rubber yield, five CGs involved in latex biosynthesis were characterized from 16 popular Hevea varieties. Novel SNPs and indels were identified and developed into markers using simple genotyping techniques like allele-specific PCR, CAPS, etc. A progeny population was genotyped using these markers to validate them, to understand their segregation pattern and to map them to a genetic linkage map. Parent-specific maps were constructed using pseudo-test cross strategy with the help of additional markers. The sequence structure information generated will be useful for future studies on gene mapping, functional relevance of coding SNPs and evolution of rubber biosynthesis genes in Hevea. Concurrently, the markers developed may serve as powerful tools for yield-based selection and for genetic diversity and pedigree studies in Hevea. Above all, the marker assays designed for genotyping could be economically carried out in any laboratory having basic molecular biology infrastructure and expertise.

Mapping of genetic loci controlling fruit linked morphological traits of melon using developed CAPS markers

BACKGROUND

Fruit morphology traits are important commercial traits that directly affect the market value. However, studying the genetic basis of these traits in un-explored botanical groups is a fundamental objective for crop genetic improvement through marker-assisted breeding.

METHODS AND RESULTS

In this study, a quantitative trait loci (QTLs) mapping strategy was used for dissecting the genomic regions of fruit linked morphological traits by single nucleotide polymorphism (SNP) based cleaved amplified polymorphism sequence (CAPS) molecular markers. Next-generation sequencing was done for the genomic sequencing of two contrasted melon lines (climacteric and non-climacteric), which revealed 97% and 96% of average coverage over the reference melon genome database, respectively. A total of 57.51% non-synonymous SNPs and 42.49% synonymous SNPs were found, which produced 149 sets of codominant markers with a 24% polymorphism rate. Total 138-F₂ derived plant populations were genotyped for linkage mapping and composite interval mapping based QTL mapping exposed 6 genetic loci, positioned over distinct chromosomes (02, 04, 08, 09, and 12) between the flanking intervals of CAPS markers, which explained an unlinked polygenic architecture in genome. Three minor QTLs of fruit weight (FWt2.1, FWt4.1, FWt9.1), one major QTL of fruit firmness (FrFir8.1), one major QTL of fruit length (FL12.1), and one major QTL of fruit shape (FS12.1) were determined and collectively explained the phenotypic variance from 5.64 to 15.64%. Fruit phenotypic correlation exhibited the significant relationship and principal component analysis also identified the potential variability. Multiple sequence alignments also indicated the significant base-mutations in the detected genetic loci, respectively.

CONCLUSION

In short, our illustrated genetic loci are expected to provide the reference insights for fine QTL mapping and candidate gene(s) mining through molecular genetic breeding approaches aimed at developing the new varieties.

Axenic interspecies and intraclonal hybrid formation in Leishmania: Successful crossings between visceral and cutaneous strains

Diseases caused by trypanosomatids are serious public health concerns in low-income endemic countries. Leishmaniasis is presented in two main clinical forms, visceral leishmaniasis-caused by L. infantum and L. donovani-and cutaneous leishmaniasis-caused by many species, including L. major, L. tropica and L. braziliensis. As for certain other trypanosomatids, sexual reproduction has been confirmed in these parasites, and formation of hybrids can contribute to virulence, drug resistance or adaptation to the host immune system. In the present work, the capability of intraclonal and interspecies genetic exchange has been investigated using three parental strains: L. donovani, L. tropica and L. major, which have been engineered to express different fluorescent proteins and antibiotic resistance markers in order to facilitate the phenotypic selection of hybrid parasites after mating events. Stationary and exponential-phase promastigotes of each species were used, in in vitro experiments, some of them containing LULO cells (an embryonic cell line derived from Lutzomyia longipalpis). Several intraclonal hybrids were obtained with L. tropica as crossing progenitor, but not with L. donovani or L. major. In interspecies crossings, three L. donovani x L. major hybrids and two L. donovani x L. tropica hybrids were isolated, thereby demonstrating the feasibility to obtain in vitro hybrids of parental lines causing different tropism of leishmaniasis. Ploidy analysis revealed an increase in DNA content in all hybrids compared to the parental strains, and nuclear analysis showed that interspecies hybrids are complete hybrids, i.e. each of them showing at least one chromosomal set from each parental. Regarding kDNA inheritance, discrepancies were observed between maxi and minicircle heritage. Finally, phenotypic studies showed either intermediate phenotypes in terms of growth profiles, or a decreased in vitro infection capacity compared to the parental cells. To the best of our knowledge, this is the first time that in vitro interspecies outcrossing has been demonstrated between Leishmania species with different tropism, thus contributing to shed light on the mechanisms underlying sexual reproduction in these parasites.

Nucleotide variation in the phytoene synthase (ClPsy1) gene contributes to golden flesh in watermelon (Citrullus lanatus L.)

A gene controlling golden flesh trait in watermelon was discovered and fine mapped to a 39.08 Kb region on chromosome 1 through a forward genetic strategy, and Cla97C01G008760 (annotated as phytoene synthase protein, ClPsy1 ) was recognized as the most likely candidate gene. Vitamin A deficiency is a worldwide public nutrition problem, and β-carotene is the precursor for vitamin A synthesis. Watermelon with golden flesh (gf, which occurs due to an accumulated abundance of β-carotene) is an important germplasm resource. In this study, a genetic analysis of segregated gf gene populations indicated that gf was controlled by a single recessive gene. BSA-seq (Bulked segregation analysis) and an initial linkage analysis placed the gf locus in a 290-Kb region on watermelon chromosome 1. Further fine mapping in a large population including over 1000 F₂ plants narrowed this region to 39.08 Kb harboring two genes, Cla97C01G008760 and Cla97C01G008770, which encode phytoene synthase (ClPsy1) and GATA zinc finger domain-containing protein, respectively. Gene sequence alignment and expression analysis between parental lines revealed Cla97C01G008760 as the best possible candidate gene for the gf trait. Nonsynonymous SNP mutations in the first exon of ClPsy1 between parental lines co-segregated with the gf trait only among individuals in the genetic population and were not related to flesh color in natural watermelon panels. Promoter sequence analysis of 26 watermelon accessions revealed two SNPs in the cis-acting element sequences corresponding to MYB and MYC2 transcription factors. RNA-seq data and qRT-PCR verification showed that two MYBs exhibited expression trends similar to that of ClPsy1 in the parental lines and may regulate the ClPsy1 expression. Further research findings indicate that the gf trait is determined not only by ClPsy1 but also by ClLCYB, ClCRTISO and ClNCED7, which play important roles in watermelon β-carotene accumulation.

VCF2CAPS-A high-throughput CAPS marker design from VCF files and its test-use on a genotyping-by-sequencing (GBS) dataset

Next-generation sequencing (NGS) is a powerful tool for massive detection of DNA sequence variants such as single nucleotide polymorphisms (SNPs), multi-nucleotide polymorphisms (MNPs) and insertions/deletions (indels). For routine screening of numerous samples, these variants are often converted into cleaved amplified polymorphic sequence (CAPS) markers which are based on the presence versus absence of restriction sites within PCR products. Current computational tools for SNP to CAPS conversion are limited and usually infeasible to use for large datasets as those generated with NGS. Moreover, there is no available tool for massive conversion of MNPs and indels into CAPS markers. Here, we present VCF2CAPS–a new software for identification of restriction endonucleases that recognize SNP/MNP/indel-containing sequences from NGS experiments. Additionally, the program contains filtration utilities not available in other SNP to CAPS converters–selection of markers with a single polymorphic cut site within a user-specified sequence length, and selection of markers that differentiate up to three user-defined groups of individuals from the analyzed population. Performance of VCF2CAPS was tested on a thoroughly analyzed dataset from a genotyping-by-sequencing (GBS) experiment. A selection of CAPS markers picked by the program was subjected to experimental verification. CAPS markers, also referred to as PCR-RFLPs, belong to basic tools exploited in plant, animal and human genetics. Our new software–VCF2CAPS–fills the gap in the current inventory of genetic software by high-throughput CAPS marker design from next-generation sequencing (NGS) data. The program should be of interest to geneticists involved in molecular diagnostics. In this paper we show a successful exemplary application of VCF2CAPS and we believe that its usefulness is guaranteed by the growing availability of NGS services.

Maize centromeric chromatin scales with changes in genome size

Centromeres are defined by the location of Centromeric Histone H3 (CENP-A/CENH3) which interacts with DNA to define the locations and sizes of functional centromeres. An analysis of 26 maize genomes including 110 fully assembled centromeric regions revealed positive relationships between centromere size and genome size. These effects are independent of variation in the amounts of the major centromeric satellite sequence CentC. We also backcrossed known centromeres into two different lines with larger genomes and observed consistent increases in functional centromere sizes for multiple centromeres. Although changes in centromere size involve changes in bound CENH3, we could not mimic the effect by overexpressing CENH3 by threefold. Literature from other fields demonstrate that changes in genome size affect protein levels, organelle size and cell size. Our data demonstrate that centromere size is among these scalable features, and that multiple limiting factors together contribute to a stable centromere size equilibrium.

Quantitative Control of Early Flowering in White Lupin (Lupinus albus L.)

White lupin (Lupinus albus L.) is a pulse annual plant cultivated from the tropics to temperate regions for its high-protein grain as well as a cover crop or green manure. Wild populations are typically late flowering and have high vernalization requirements. Nevertheless, some early flowering and thermoneutral accessions were found in the Mediterranean basin. Recently, quantitative trait loci (QTLs) explaining flowering time variance were identified in bi-parental population mapping, however, phenotypic and genotypic diversity in the world collection has not been addressed yet. In this study, a diverse set of white lupin accessions (n = 160) was phenotyped for time to flowering in a controlled environment and genotyped with PCR-based markers (n = 50) tagging major QTLs and selected homologs of photoperiod and vernalization pathway genes. This survey highlighted quantitative control of flowering time in white lupin, providing statistically significant associations for all major QTLs and numerous regulatory genes, including white lupin homologs of CONSTANS, FLOWERING LOCUS T, FY, MOTHER OF FT AND TFL1, PHYTOCHROME INTERACTING FACTOR 4, SKI-INTERACTING PROTEIN 1, and VERNALIZATION INDEPENDENCE 3. This revealed the complexity of flowering control in white lupin, dispersed among numerous loci localized on several chromosomes, provided economic justification for future genome-wide association studies or genomic selection rather than relying on simple marker-assisted selection.

Genetic variability and biodiversity of Ukrainian Gray cattle by the BoLA-DRB3 gene

At the current stage of genetic studies of cattle, more and more attention is being drawn to autochthonous breeds. Native cattle have a number of prominent phenotypic traits and have preserved unique genes and their combinations lost by modern commercial breeds, which would be valuable to use in selective programs. We surveyed polymorphism of the Ukrainian autochthonous Gray breed according to alleles of exon 2 of the BoLA-DRB3 gene. The uniqueness of the gene lies in the broad variability of its allele variants. Significant informativeness at DNA level is quite important for genetic studies. We surveyed allele polymorphism using the PCR-RLFP method on DNA isolated from 88 samples of blood of cows and 5 samples of sperm. We identified 28 alleles, of which 23 variants were nomenclature ones and 5 (jba, *jab, *jbb, *nad and *nda) were “without established nomenclature”, their share accounting for 8.9%. Four alleles *06, *12, *16 and *jba had a frequency above 5% and occupied 69.9% of the breed’s allele fund overall. The commonest allele was BoLA-DRB3.2*16 (44.1%). In total, we found 40 genotypes. Considering the significant dominance of variant *16, as expected, 5 genotypes with its inclusion occurred: *16/*16, *12/*16, *06/*16, *16/*24 and *jba/*16. It was present in the genotype of two out three studied animals. Parameters of heterozygosity, effective number of alleles, Shannon and Pielou indices indicate that Ukrainian Gray cattle are characterized by lowest level of genetic variability and biodiversity according to the BoLA-DRB3 gene compared with other breeds. Due to significant dominance of allele *16, the breed has no inbred motifs. We noted deviation toward increase in homozygosity without deviations from the norm of the distribution according to Hardy-Weinberg equilibrium. The obtained results will be used for genetic-populational programs with the purpose of improving the genetic potential of cattle breeds in terms of economically beneficial traits and diseases of cattle.

New data for systematics of the genus Zeravschania (Apiaceae) and related taxa based on the nrDNA (ITS) sequences

All 13 species of the genus Zeravschania were studied in molecular and morphological analyzes. The molecular study showed that Zeravschania is a clearly paraphyletic taxon including monotypic genera Demavendia and Haussknechtia. Further phylogenetic, anatomical and morphological studies are needed to confirm taxonomic position of these genera. Zeravschania afghanica confirms its placement in the genus Zeravschania, not in Cephalopodum or Selinum. Dichoropetalum knappii confirms its placement in the genus Dichoropetalum, not in Zeravschania.

Mutation of the ALBOSTRIANS Ohnologous Gene HvCMF3 Impairs Chloroplast Development and Thylakoid Architecture in Barley

Gene pairs resulting from whole genome duplication (WGD), so-called ohnologous genes, are retained if at least one member of the pair undergoes neo- or sub-functionalization. Phylogenetic analyses of the ohnologous genes ALBOSTRIANS (HvAST/HvCMF7) and ALBOSTRIANS-LIKE (HvASL/HvCMF3) of barley (Hordeum vulgare) revealed them as members of a subfamily of genes coding for CCT motif (CONSTANS, CONSTANS-LIKE and TIMING OF CAB1) proteins characterized by a single CCT domain and a putative N-terminal chloroplast transit peptide. Recently, we showed that HvCMF7 is needed for chloroplast ribosome biogenesis. Here we demonstrate that mutations in HvCMF3 lead to seedlings delayed in development. They exhibit a yellowish/light green - xantha - phenotype and successively develop pale green leaves. Compared to wild type, plastids of mutant seedlings show a decreased PSII efficiency, impaired processing and reduced amounts of ribosomal RNAs; they contain less thylakoids and grana with a higher number of more loosely stacked thylakoid membranes. Site-directed mutagenesis of HvCMF3 identified a previously unknown functional domain, which is highly conserved within this subfamily of CCT domain containing proteins. HvCMF3:GFP fusion constructs were localized to plastids and nucleus. Hvcmf3Hvcmf7 double mutants exhibited a xantha-albino or albino phenotype depending on the strength of molecular lesion of the HvCMF7 allele. The chloroplast ribosome deficiency is discussed as the primary observed defect of the Hvcmf3 mutants. Based on our observations, the genes HvCMF3 and HvCMF7 have similar but not identical functions in chloroplast development of barley supporting our hypothesis of neo-/sub-functionalization between both ohnologous genes.

Identification of the 'Haryejosaeng' mandarin cultivar by multiplex PCR-based SNP genotyping

Most satsuma mandarin (Citrus unshiu Marc.) cultivars are difficult to identify in the seedling stage based only on morphological traits. Therefore, simple polymerase chain reaction (PCR)-based single-nucleotide polymorphism (SNP) markers were developed to specifically and rapidly distinguish the 'Haryejosaeng' cultivar, which is generally supplied to breeders of other satsuma mandarin cultivars. SNP markers were verified using high-resolution melt (HRM)-specific primers. PCR was performed to distinguish 'Haryejosaeng' from eight other satsuma mandarin cultivars using six SNP markers (P1-P6) specific for 'Haryejosaeng', with one negative control SNP primer pair. The best results were obtained using three SNP markers (P1, P2, and P5). In the multiplex PCR, markers P1, P2, and P5 yielded 165-, 150-, and 526-base pair amplicons, respectively, in 'Haryejosaeng', distinguishing it from other satsuma mandarin cultivars. The selected SNP markers were validated by HRM with HRM-specific primers. The multiplex PCR with P1/P5 and P2/P5 also identified 'Haryejosaeng' obtained from a farm growing 17 different cultivars of satsuma mandarin. Specific SNP molecular markers were determined for accurately identifying the 'Haryejosaeng' cultivar by multiplex PCR to save the time and costs associated with its supply to breeders of satsuma mandarin.

Innovative transcriptome-based genotyping highlights environmentally responsive genes for phenology, growth and yield in a non-model grain legume

The narrow-leafed lupin, Lupinus angustifolius L., is a grain legume crop, cultivated both as a green manure and as a source of protein for animal feed and human food production. During its domestication process, numerous agronomic traits were improved, however, only two trait-related genes were identified hitherto, both by linkage mapping. Genome-wide association studies (GWAS), exploiting genomic sequencing, did not select any novel candidate gene. In the present study, an innovative method of 3'-end reduced representation transcriptomic profiling, a massive analysis of cDNA ends, has been used for genotyping of 126 L. angustifolius lines surveyed by field phenotyping. Significant genotype × environment interactions were identified for all phenology and yield traits analysed. Principal component analysis of population structure evidenced European domestication bottlenecks, visualized by clustering of breeding materials and cultivars. GWAS provided contribution towards deciphering vernalization pathway in legumes, and, apart from highlighting known domestication loci (Ku/Julius and mol), designated novel candidate genes for L. angustifolius traits. Early phenology was associated with genes from vernalization, cold-responsiveness and phosphatidylinositol signalling pathways whereas high yield with genes controlling photosynthesis performance and abiotic stress (drought or heat) tolerance. PCR-based toolbox was developed and validated to enable tracking desired alleles in marker-assisted selection. Narrow-leafed lupin was genotyped with an innovative method of transcriptome profiling and phenotyped for phenology, growth and yield traits in field. Early phenology was found associated with genes from cold-response, vernalization and phosphatidylinositol signalling pathways, whereas high yield with genes running photosystem II and drought or heat stress response. Key loci were supplied with PCR-based toolbox for marker-assisted selection.

A Genome-Wide Analysis of the Pentatricopeptide Repeat (PPR) Gene Family and PPR-Derived Markers for Flesh Color in Watermelon (Citrullus lanatus)

Watermelon (Citrullus lanatus) is an economically important fruit crop grown for consumption of its large edible fruit flesh. Pentatricopeptide-repeat (PPR) encoding genes, one of the large gene families in plants, are important RNA-binding proteins involved in the regulation of plant growth and development by influencing the expression of organellar mRNA transcripts. However, systematic information regarding the PPR gene family in watermelon remains largely unknown. In this comprehensive study, we identified and characterized a total of 422 C. lanatus PPR (ClaPPR) genes in the watermelon genome. Most ClaPPRs were intronless and were mapped across 12 chromosomes. Phylogenetic analysis showed that ClaPPR proteins could be divided into P and PLS subfamilies. Gene duplication analysis suggested that 11 pairs of segmentally duplicated genes existed. In-silico expression pattern analysis demonstrated that ClaPPRs may participate in the regulation of fruit development and ripening processes. Genotyping of 70 lines using 4 single nucleotide polymorphisms (SNPs) from 4 ClaPPRs resulted in match rates of over 0.87 for each validated SNPs in correlation with the unique phenotypes of flesh color, and could be used in differentiating red, yellow, or orange watermelons in breeding programs. Our results provide significant insights for a comprehensive understanding of PPR genes and recommend further studies on their roles in watermelon fruit growth and ripening, which could be utilized for cultivar development of watermelon.

Development of PCR-based markers and whole-genome selection model for anthracnose resistance in white lupin (Lupinus albus L.)

White lupin (Lupinus albus L.) is a high-protein grain legume crop, grown since ancient Greece and Rome. Despite long domestication history, its cultivation remains limited, partly because of susceptibility to anthracnose. Only some late-flowering, bitter, low-yielding landraces from Ethiopian mountains displayed resistance to this devastating disease. The resistance is controlled by various genes, thereby complicating the breeding efforts. The objective of this study was developing tools for molecular tracking of Ethiopian resistance genes based on genotyping-by-sequencing (GBS) data, envisaging linkage mapping and genomic selection approaches. Twenty GBS markers from two major quantitative trait loci (QTLs), antr04_1/antr05_1 and antr04_2/antr05_2, were converted to PCR-based markers using assigned transcriptome sequences. Newly developed markers improved mapping resolution around both anthracnose resistance loci, providing more precise QTL estimation. PCR-based screening of diversified domesticated and primitive germplasm revealed the high specificity of two markers for the antr04_1/antr05_1 locus (TP222136 and TP47110) and one for the antr04_2/antr05_2 locus (TP338761), highlighted by simple matching coefficients of 0.96 and 0.89, respectively. Moreover, a genomic selection approach based on GBS data of a recombinant inbred line mapping population was assessed, providing an average predictive ability of 0.56. These tools can be used for preselection of candidate white lupin germplasm for anthracnose resistance assays.

Role of NGS and SNP genotyping methods in sugarcane improvement programs

Sugarcane (Saccharum spp.) is one of the most economically significant crops because of its high sucrose content and it is a promising biomass feedstock for biofuel production. Sugarcane genome sequencing and analysis is a difficult task due to its heterozygosity and polyploidy. Long sequence read technologies, PacBio Single-Molecule Real-Time (SMRT) sequencing, the Illumina TruSeq, and the Oxford Nanopore sequencing could solve the problem of genome assembly. On the applications side, next generation sequencing (NGS) technologies played a major role in the discovery of single nucleotide polymorphism (SNP) and the development of low to high throughput genotyping platforms. The two mainstream high throughput genotyping platforms are the SNP microarray and genotyping by sequencing (GBS). This paper reviews the NGS in sugarcane genomics, genotyping methodologies, and the choice of these methods. Array-based SNP genotyping is robust, provides consistent SNPs, and relatively easier downstream data analysis. The GBS method identifies large scale SNPs across the germplasm. A combination of targeted GBS and array-based genotyping methods should be used to increase the accuracy of genomic selection and marker-assisted breeding.

Genome-based breeding approaches in major vegetable crops

Vegetable crops are major nutrient sources for humanity and have been well-cultivated since thousands of years of domestication. With the rapid development of next-generation sequencing and high-throughput genotyping technologies, the reference genome of more than 20 vegetables have been well-assembled and published. Resequencing approaches on large-scale germplasm resources have clarified the domestication and improvement of vegetable crops by human selection; its application on genetic mapping and quantitative trait locus analysis has led to the discovery of key genes and molecular markers linked to important traits in vegetables. Moreover, genome-based breeding has been utilized in many vegetable crops, including Solanaceae, Cucurbitaceae, Cruciferae, and other families, thereby promoting molecular breeding at a single-nucleotide level. Thus, genome-wide SNP markers have been widely used, and high-throughput genotyping techniques have become one of the most essential methods in vegetable breeding. With the popularization of gene editing technology research on vegetable crops, breeding efficiency can be rapidly increased, especially by combining the genomic and variomic information of vegetable crops. This review outlines the present genome-based breeding approaches used for major vegetable crops to provide insights into next-generation molecular breeding for the increasing global population.

High Resolution Mapping of a Hordeum bulbosum-Derived Powdery Mildew Resistance Locus in Barley Using Distinct Homologous Introgression Lines

Powdery mildew caused by Blumeria graminis f. sp. hordei (Bgh) is one of the main foliar diseases in barley (Hordeum vulgare L.; Hv). Naturally occurring resistance genes used in barley breeding are a cost effective and environmentally sustainable strategy to minimize the impact of pathogens, however, the primary gene pool of H. vulgare contains limited diversity owing to recent domestication bottlenecks. To ensure durable resistance against this pathogen, more genes are required that could be unraveled by investigation of secondary barley gene-pool. A large set of Hordeum bulbosum (Hb) introgression lines (ILs) harboring a diverse set of desirable resistance traits have been developed and are being routinely used as source of novel diversity in gene mapping studies. Nevertheless, this strategy is often compromised by a lack of recombination between the introgressed fragment and the orthologous chromosome of the barley genome. In this study, we fine-mapped a Hb gene conferring resistance to barley powdery mildew. The initial genotyping of two Hb ILs mapping populations with differently sized 2HS introgressions revealed severely reduced interspecific recombination in the region of the introgressed segment, preventing precise localization of the gene. To overcome this difficulty, we developed an alternative strategy, exploiting intraspecific recombination by crossing two Hv/Hb ILs with collinear Hb introgressions, one of which carries a powdery mildew resistance gene, while the other doesn't. The intraspecific recombination rate in the Hb-introgressed fragment of 2HS was approximately 20 times higher than it was in the initial simple ILs mapping populations. Using high-throughput genotyping-by-sequencing (GBS), we allocated the resistance gene to a 1.4 Mb interval, based on an estimate using the Hv genome as reference, in populations of only 103 and 146 individuals, respectively, similar to what is expected at this locus in barley. The most likely candidate resistance gene within this interval is part of the coiled-coil nucleotide-binding-site leucine-rich-repeat (CC-NBS-LLR) gene family, which is over-represented among genes conferring strong dominant resistance to pathogens. The reported strategy can be applied as a general strategic approach for identifying genes underlying traits of interest in crop wild relatives.

Genetic mapping of green curd gene Gr in cauliflower

Gr5.1 is the major locus for cauliflower green curd color and mapped to an interval of 236 Kbp with four most likely candidate genes. Cauliflower with colored curd enhances not only the visual appeal but also the nutritional value of the crop. Green cauliflower results from ectopic development of chloroplasts in the normal white curd. However, the underlying genetic basis is unknown. In this study, we employed QTL-seq analysis to identify the loci that were associated with green curd phenotype in cauliflower. A F2 population was generated following a cross between a white curd (Stovepipe) and a green curd (ACX800) cauliflower plants. By whole-genome resequencing and SNP analysis of green and white F2 bulks, two QTLs were detected on chromosomes 5 (Gr5.1) and 7 (Gr7.1). Validation by traditional genetic mapping with CAPS markers suggested that Gr5.1 represented a major QTL, whereas Gr7.1 had a minor effect. Subsequent high-resolution mapping of Gr5.1 in the second large F2 population with additional CAPS markers narrowed down the target region to a genetic and physical distance of 0.3 cM and 236 Kbp, respectively. This region contained 35 genes with four of them representing the best candidates for the green curd phenotype in cauliflower. They are LOC106295953, LOC106343833, LOC106345143, and LOC106295954, which encode UMP kinase, DEAD-box RNA helicase 51-like, glutathione S-transferase T3-like, and protein MKS1, respectively. These findings lay a solid foundation for the isolation of the Gr gene and provide a potential for marker-assisted selection of the green curd trait in cauliflower breeding. The eventual isolation of Gr will also facilitate better understanding of chloroplast biogenesis and development in plants.

Linkage Mapping and Comparative Transcriptome Analysis of Firmness in Watermelon (Citrullus lanatus)

Watermelon fruit texture and quality are determined by flesh firmness. As a quality trait, flesh firmness is controlled by multigenes. Defining the key regulatory factors of watermelon flesh firmness is of great significance for watermelon genetic breeding. In this study, the hard-flesh egusi seed watermelon PI186490 was used as the male parent, the soft-flesh cultivated watermelon W1-1 was used as the female parent, and 175 F₂ generations were obtained from selfing F₁. Primary mapping of the major genes controlling center flesh firmness was achieved by bulked-segregant analysis (BSA)-Seq analysis and molecular marker technology. Finally, major genes were delimited in the physical interval between 6,210,787 and 7,742,559 bp on chromosome 2 and between 207,553 and 403,137 bp on chromosome 8. The content of each cell wall component and hormone was measured, and comparative transcriptome analysis was performed during fruit development in watermelon. The protopectin, cellulose, hemicellulose, indole-3-acetic acid (IAA) and abscisic acid (ABA) contents were measured, and paraffin sections were made during the three fruit developmental stages. The results revealed that protopectin, celluloses, and hemicelluloses exhibited similar trends for flesh firmness, while the IAA and ABA concentrations continued to decrease with fruit ripening. Paraffin sections showed that PI186490 cells were more numerous, were more tightly packed, had clearer cell wall edges and had thicker cell walls than W1-1 cells at every developmental stage. Comparative transcriptome analysis was conducted on RNA samples of flesh during fruit development and ripening in W1-1 and PI186490. The results from the localization interval transcriptome analysis showed that Cla016033 (DUF579 family member), which may influence the cell wall component contents to adjust the flesh firmness in watermelon fruit, was different in W1-1 and PI186490 and that Cla012507 (MADS-box transcription factor) may be involved in the regulation of fruit ripening and affect the hardness of watermelon fruit.

Candidate Domestication-Related Genes Revealed by Expression Quantitative Trait Loci Mapping of Narrow-Leafed Lupin (Lupinus angustifolius L.)

The last century has witnessed rapid domestication of the narrow-leafed lupin (Lupinus angustifolius L.) as a grain legume crop, exploiting discovered alleles conferring low-alkaloid content (iucundus), vernalization independence (Ku and Julius), and reduced pod shattering (lentus and tardus). In this study, a L. angustifolius mapping population was subjected to massive analysis of cDNA ends (MACE). The MACE yielded 4185 single nucleotide polymorphism (SNP) markers for linkage map improvement and 30,595 transcriptomic profiles for expression quantitative trait loci (eQTL) mapping. The eQTL highlighted a high number of cis- and trans-regulated alkaloid biosynthesis genes with gene expression orchestrated by a regulatory agent localized at iucundus locus, supporting the concept that ETHYLENE RESPONSIVE TRANSCRIPTION FACTOR RAP2-7 may control low-alkaloid phenotype. The analysis of Ku shed light on the vernalization response via FLOWERING LOCUS T and FD regulon in L. angustifolius, providing transcriptomic evidence for the contribution of several genes acting in C-repeat binding factor (CBF) cold responsiveness and in UDP-glycosyltransferases pathways. Research on lentus selected a DUF1218 domain protein as a candidate gene controlling the orientation of the sclerified endocarp and a homolog of DETOXIFICATION14 for purplish hue of young pods. An ABCG transporter was identified as a hypothetical contributor to sclerenchyma fortification underlying tardus phenotype.

Rapid identification and characterization of genetic loci for defective kernel in bread wheat

BACKGROUND

Wheat is a momentous crop and feeds billions of people in the world. The improvement of wheat yield is very important to ensure world food security. Normal development of grain is the essential guarantee for wheat yield formation. The genetic study of grain phenotype and identification of key genes for grain filling are of great significance upon dissecting the molecular mechanism of wheat grain morphogenesis and yield potential.

RESULTS

Here we identified a pair of defective kernel (Dek) isogenic lines, BL31 and BL33, with plump and shrunken mature grains, respectively, and constructed a genetic population from the BL31/BL33 cross. Ten chromosomes had higher frequency of polymorphic single nucleotide polymorphism (SNP) markers between BL31 and BL33 using Wheat660K chip. Totally 783 simple sequence repeat (SSR) markers were chosen from the above chromosomes and 15 of these were integrated into two linkage groups using the genetic population. Genetic mapping identified three QTL, QDek.caas-3BS.1, QDek.caas-3BS.2 and QDek.caas-4AL, explaining 14.78-18.17%, 16.61-21.83% and 19.08-28.19% of phenotypic variances, respectively. Additionally, five polymorphic SNPs from Wheat660K were successfully converted into cleaved amplified polymorphic sequence (CAPS) markers and enriched the target regions of the above QTL. Biochemical analyses revealed that BL33 has significantly higher grain sucrose contents at filling stages and lower mature grain starch contents than BL31, indicating that the Dek QTL may be involved in carbohydrate metabolism. As such, the candidate genes for each QTL were predicated according to International Wheat Genome Sequence Consortium (IWGSC) RefSeq v1.0.

CONCLUSIONS

Three major QTL for Dek were identified and their causal genes were predicted, laying a foundation to conduct fine mapping and dissect the regulatory mechanism underlying Dek trait in wheat.

Whole genome resequencing of watermelons to identify single nucleotide polymorphisms related to flesh color and lycopene content

Cultivated watermelon (Citrullus lanatus) is one of the most important food crops in the Cucurbitaceae family. Diversification after domestication has led cultivated watermelons to exhibit diverse fruit flesh colors, including red, yellow, and orange. Recently, there has been increased interest in red-fleshed watermelons because they contain the antioxidant cis-isomeric lycopene. We performed whole genome resequencing (WGRS) of 24 watermelons with different flesh colors to identify single-nucleotide polymorphisms (SNPs) related to high lycopene content. The resequencing data revealed 203,894-279,412 SNPs from read mapping between inbred lines and the 97103 reference genome. In total, 295,065 filtered SNPs were identified, which had an average polymorphism information content of 0.297. Most of these SNPs were intergenic (90.1%) and possessed a transversion (Tv) rate of 31.64%. Overall, 2,369 SNPs were chosen at 0.5 Mb physical intervals to analyze genetic diversity across the 24 inbred lines. A neighbor-joining dendrogram and principal coordinate analysis (PCA) based on the 2,369 SNPs revealed that the 24 inbred lines could be grouped into high and low lycopene-type watermelons. In addition, we analyzed SNPs that could discriminate high lycopene content, red-fleshed watermelon from low lycopene, yellow or orange watermelon inbred lines. For validation, 19 SNPs (designated as WMHL1-19) were chosen randomly, and cleavage amplified polymorphic sequence (CAPS) markers were designed. Genotyping of the above 24 lines and 12 additional commercial cultivars using WMHL1-19 CAPS markers resulted in match rates of over 0.92 for most validated markers in correlation with the flesh color phenotypes. Our results provide valuable genomic information regarding the high lycopene content phenotype of red-fleshed cultivated watermelons, and the identified SNPs will be useful for the development of molecular markers in the marker-assisted breeding of watermelons with high lycopene content.

Development of gene-based molecular markers tagging low alkaloid pauper locus in white lupin (Lupinus albus L.)

White lupin (Lupinus albus L.) is a legume grain crop cultivated since ancient Greece and Egypt. Modern white lupin cultivars are appreciated as a source of protein with positive nutraceutical impact. However, white lupins produce anti-nutritional compounds, quinolizidine alkaloids, which provide bitter taste and have a negative influence on human health. During domestication of this species, several recessive alleles at unlinked loci controlling low alkaloid content were selected. One of these loci, pauper, was exploited worldwide providing numerous low-alkaloid cultivars. However, molecular tracking of pauper has been hampered due to the lack of diagnostic markers. In the present study, the synteny-based approach was harnessed to target pauper locus. Single-nucleotide polymorphisms flanking pauper locus on white lupin linkage map as well as candidate gene sequences elucidated from the narrow-leafed lupin (L. angustifolius L.) chromosome segment syntenic to the pauper linkage group region were transformed to PCR-based molecular markers. These markers were analyzed both in the mapping population and world germplasm collection. From fourteen markers screened, eleven were localized at a distance below 1.5 cM from this locus, including five co-segregating with pauper. The linkage of these markers was confirmed by high LOD values (up to 58.4). Validation performed in the set of 127 bitter and 23 sweet accessions evidenced high applicability of one marker, LAGI01_35805_F1_R1, for pauper locus selection, highlighted by the low ratio of false-positive scores (2.5%). LAGI01_35805 represents a homolog of L. angustifolius acyltransferase-like (LaAT) gene which might hypothetically participate in the alkaloid biosynthesis process in lupins.

Leaf Variegation and Impaired Chloroplast Development Caused by a Truncated CCT Domain Gene in albostrians Barley

Chloroplasts fuel plant development and growth by converting solar energy into chemical energy. They mature from proplastids through the concerted action of genes in both the organellar and the nuclear genome. Defects in such genes impair chloroplast development and may lead to pigment-deficient seedlings or seedlings with variegated leaves. Such mutants are instrumental as tools for dissecting genetic factors underlying the mechanisms involved in chloroplast biogenesis. Characterization of the green-white variegated albostrians mutant of barley (Hordeum vulgare) has greatly broadened the field of chloroplast biology, including the discovery of retrograde signaling. Here, we report identification of the ALBOSTRIANS gene HvAST (also known as Hordeum vulgare CCT Motif Family gene 7, HvCMF7) by positional cloning as well as its functional validation based on independently induced mutants by Targeting Induced Local Lesions in Genomes (TILLING) and RNA-guided clustered regularly interspaced short palindromic repeats-associated protein 9 endonuclease-mediated gene editing. The phenotypes of the independent HvAST mutants imply residual activity of HvCMF7 in the original albostrians allele conferring an imperfect penetrance of the variegated phenotype even at homozygous state of the mutation. HvCMF7 is a homolog of the Arabidopsis (Arabidopsis thaliana) CONSTANS, CO-like, and TOC1 (CCT) Motif transcription factor gene CHLOROPLAST IMPORT APPARATUS2, which was reported to be involved in the expression of nuclear genes essential for chloroplast biogenesis. Notably, in barley we localized HvCMF7 to the chloroplast, without any clear evidence for nuclear localization.

Identification and characterization of a new dwarf locus DS-4 encoding an Aux/IAA7 protein in Brassica napus

A dominant dwarfing gene, ds - 4 , encodes an Aux/IAA protein that negatively regulates plant stature through an auxin signaling pathway. Dwarfism is an important agronomic trait affecting yield in many crop species. The molecular mechanisms underlying dwarfism in oilseed rape (Brassica napus) are poorly understood, restricting the progress of breeding dwarf varieties in this species. Here, we identified and characterized a new dwarf locus, DS-4, in B. napus. Next-generation sequencing-assisted genetic mapping and candidate gene analysis found that DS-4 encodes a nucleus-targeted auxin/indole-3-acetic acid (Aux/IAA) protein. A substitution (P87L) was found in the highly conserved degron motif of the Aux/IAA7 protein in the ds-4 mutant. This mutation co-segregated with the phenotype of individuals in the BC₁F₂ population. The P87L substitution was confirmed as the cause of the extreme dwarf phenotype by ectopic expression of the mutant allele BnaC05.iaa7 (equivalent to ds-4) in Arabidopsis. The P87L substitution blocked the interaction of BnaC05.iaa7 with TRANSPORT INHIBITOR RESPONSE 1 in the presence of auxin. The BnaC05.IAA7 gene is highly expressed in the cotyledons, hypocotyls, stems and leaves, but weakly in the roots and seeds of B. napus. Our findings provide new insights into the molecular mechanisms underlying dominant (gain-of-function) dwarfism in B. napus. Our identification of a distinct gene locus controlling plant height may help to improve lodging resistance in oilseed rape.

Characterisation of barley resistance to rhynchosporium on chromosome 6HS

Major resistance gene to rhynchosporium, Rrs18, maps close to the telomere on the short arm of chromosome 6H in barley. Rhynchosporium or barley scald caused by a fungal pathogen Rhynchosporium commune is one of the most destructive and economically important diseases of barley in the world. Testing of Steptoe × Morex and CIho 3515 × Alexis doubled haploid populations has revealed a large effect QTL for resistance to R. commune close to the telomere on the short arm of chromosome 6H, present in both populations. Mapping markers flanking the QTL from both populations onto the 2017 Morex genome assembly revealed a rhynchosporium resistance locus independent of Rrs13 that we named Rrs18. The causal gene was fine mapped to an interval of 660 Kb using Steptoe × Morex backcross 1 S₂ and S₃ lines with molecular markers developed from Steptoe exome capture variant calling. Sequencing RNA from CIho 3515 and Alexis revealed that only 4 genes within the Rrs18 interval were transcribed in leaf tissue with a serine/threonine protein kinase being the most likely candidate for Rrs18.

High Resolution Genetic and Physical Mapping of a Major Powdery Mildew Resistance Locus in Barley

Powdery mildew caused by Blumeria graminis f. sp. hordei is a foliar disease with highly negative impact on yield and grain quality in barley. Thus, breeding for powdery mildew resistance is an important goal and requires constantly the discovery of new sources of natural resistance. Here, we report the high resolution genetic and physical mapping of a dominant race-specific powdery mildew resistance locus, originating from an Ethiopian spring barley accession 'HOR2573,' conferring resistance to several modern mildew isolates. High-resolution genetic mapping narrowed down the interval containing the resistance locus to a physical span of 850 kb. Four candidate genes with homology to known disease resistance gene families were identified. The mapped resistance locus coincides with a previously reported resistance locus from Hordeum laevigatum, suggesting allelism at the same locus in two different barley lines. Therefore, we named the newly mapped resistance locus from HOR2573 as MlLa-H. The reported co-segregating and flanking markers may provide new tools for marker-assisted selection of this resistance locus in barley breeding.

Characterization of viral RNA splicing using whole-transcriptome datasets from host species

RNA alternative splicing (AS) is an important post-transcriptional mechanism enabling single genes to produce multiple proteins. It has been well demonstrated that viruses deploy host AS machinery for viral protein productions. However, knowledge on viral AS is limited to a few disease-causing viruses in model species. Here we report a novel approach to characterizing viral AS using whole transcriptome dataset from host species. Two insect transcriptomes (Acheta domesticus and Planococcus citri) generated in the 1,000 Insect Transcriptome Evolution (1KITE) project were used as a proof of concept using the new pipeline. Two closely related densoviruses (Acheta domesticus densovirus, AdDNV, and Planococcus citri densovirus, PcDNV, Ambidensovirus, Densovirinae, Parvoviridae) were detected and analyzed for AS patterns. The results suggested that although the two viruses shared major AS features, dramatic AS divergences were observed. Detailed analysis of the splicing junctions showed clusters of AS events occurred in two regions of the virus genome, demonstrating that transcriptome analysis could gain valuable insights into viral splicing. When applied to large-scale transcriptomics projects with diverse taxonomic sampling, our new method is expected to rapidly expand our knowledge on RNA splicing mechanisms for a wide range of viruses.

Unravelling the Carbohydrate-Binding Preferences of the Carbohydrate-Binding Modules of AMP-Activated Protein Kinase

The β subunit of adenosine monophosphate (AMP)-activated protein kinase (AMPK), which exists as two isoforms (β1 and β2) in humans, has a carbohydrate-binding module (CBM) that interacts with glycogen. Although the β1- and β2-CBMs are structurally similar, with strictly conserved ligand-contact residues, they show different carbohydrate affinities. β2-CBM shows the strongest affinity for both branched and unbranched oligosaccharides and it has recently been shown that a Thr insertion into β2-CBM (Thr101) forms a pocket to accommodate branches. This insertion does not explain why β2-CBM binds all carbohydrates with stronger affinity. Herein, it is shown that residue 134 (Val for β2 and Thr for β1), which does not come into contact with a carbohydrate, appears to account for the affinity difference. Characterisation by NMR spectroscopy, however, suggests that mutant β2-Thr101Δ/Val134Thr differs from that of β1-CBM, and mutant β1-Thr101ins/Thr134Val differs from that of β2-CBM. Furthermore, these mutants are less stable to chemical denaturation, relative to that of wild-type β-CBMs, which confounds the affinity analyses. To support the importance of Thr101 and Val134, the ancestral CBM has been constructed. This CBM retains Thr101 and Val134, which suggests that the extant β1-CBM has a modest loss of function in carbohydrate binding. Because the ancestor bound carbohydrate with equal affinity to that of β2-CBM, it is concluded that residue 134 plays an indirect role in carbohydrate binding.

Development and Applications of a High Throughput Genotyping Tool for Polyploid Crops: Single Nucleotide Polymorphism (SNP) Array

Polypoid species play significant roles in agriculture and food production. Many crop species are polyploid, such as potato, wheat, strawberry, and sugarcane. Genotyping has been a daunting task for genetic studies of polyploid crops, which lags far behind the diploid crop species. Single nucleotide polymorphism (SNP) array is considered to be one of, high-throughput, relatively cost-efficient and automated genotyping approaches. However, there are significant challenges for SNP identification in complex, polyploid genomes, which has seriously slowed SNP discovery and array development in polyploid species. Ploidy is a significant factor impacting SNP qualities and validation rates of SNP markers in SNP arrays, which has been proven to be a very important tool for genetic studies and molecular breeding. In this review, we (1) discussed the pros and cons of SNP array in general for high throughput genotyping, (2) presented the challenges of and solutions to SNP calling in polyploid species, (3) summarized the SNP selection criteria and considerations of SNP array design for polyploid species, (4) illustrated SNP array applications in several different polyploid crop species, then (5) discussed challenges, available software, and their accuracy comparisons for genotype calling based on SNP array data in polyploids, and finally (6) provided a series of SNP array design and genotype calling recommendations. This review presents a complete overview of SNP array development and applications in polypoid crops, which will benefit the research in molecular breeding and genetics of crops with complex genomes.

The genome sequence of sweet cherry (Prunus avium) for use in genomics-assisted breeding

We determined the genome sequence of sweet cherry (Prunus avium) using next-generation sequencing technology. The total length of the assembled sequences was 272.4 Mb, consisting of 10,148 scaffold sequences with an N50 length of 219.6 kb. The sequences covered 77.8% of the 352.9 Mb sweet cherry genome, as estimated by k-mer analysis, and included >96.0% of the core eukaryotic genes. We predicted 43,349 complete and partial protein-encoding genes. A high-density consensus map with 2,382 loci was constructed using double-digest restriction site–associated DNA sequencing. Comparing the genetic maps of sweet cherry and peach revealed high synteny between the two genomes; thus the scaffolds were integrated into pseudomolecules using map- and synteny-based strategies. Whole-genome resequencing of six modern cultivars found 1,016,866 SNPs and 162,402 insertions/deletions, out of which 0.7% were deleterious. The sequence variants, as well as simple sequence repeats, can be used as DNA markers. The genomic information helps us to identify agronomically important genes and will accelerate genetic studies and breeding programs for sweet cherries. Further information on the genomic sequences and DNA markers is available in DBcherry (http://cherry.kazusa.or.jp (8 May 2017, date last accessed)).

Candidate genes and molecular markers associated with heat tolerance in colonial Bentgrass

Elevated temperature is a major abiotic stress limiting the growth of cool-season grasses during the summer months. The objectives of this study were to determine the genetic variation in the expression patterns of selected genes involved in several major metabolic pathways regulating heat tolerance for two genotypes contrasting in heat tolerance to confirm their status as potential candidate genes, and to identify PCR-based markers associated with candidate genes related to heat tolerance in a colonial (Agrostis capillaris L.) x creeping bentgrass (Agrostis stolonifera L.) hybrid backcross population. Plants were subjected to heat stress in controlled-environmental growth chambers for phenotypic evaluation and determination of genetic variation in candidate gene expression. Molecular markers were developed for genes involved in protein degradation (cysteine protease), antioxidant defense (catalase and glutathione-S-transferase), energy metabolism (glyceraldehyde-3-phosphate dehydrogenase), cell expansion (expansin), and stress protection (heat shock proteins HSP26, HSP70, and HSP101). Kruskal-Wallis analysis, a commonly used non-parametric test used to compare population individuals with or without the gene marker, found the physiological traits of chlorophyll content, electrolyte leakage, normalized difference vegetative index, and turf quality were associated with all candidate gene markers with the exception of HSP101. Differential gene expression was frequently found for the tested candidate genes. The development of candidate gene markers for important heat tolerance genes may allow for the development of new cultivars with increased abiotic stress tolerance using marker-assisted selection.

PCR-based assays for validation of single nucleotide polymorphism markers in rice and mungbean

BACKGROUND

Single nucleotide polymorphism (SNP) markers are the method of choice for genetic analyses including diversity and quantitative trait loci (QTL) studies. Marker validation is essential for QTL studies, but the cost and workload are considerable when large numbers of markers need to be verified. Marker systems with low development costs would be most suitable for this task.

RESULTS

We have tested allele specific polymerase chain reaction (PCR), tetra markers and a genotyping tool based on the single strand specific nuclease CEL-I to verify randomly selected SNP markers identified previously either with a SNP array or by genotyping by sequencing in rice and mungbean, respectively. The genotyping capacity of allele-specific PCR and tetra markers was affected by the sequence context surrounding the SNP; SNPs located in repeated sequences and in GC-rich stretches could not be correctly identified. In contrast, CEL-I digestion of mixed fragments produced from test and reference DNA reliably pinpointed the correct genotypes, yet scoring of the genotypes became complicated when multiple SNPs were present in the PCR fragments. A cost analysis showed that as long the sample number remains small, CEL-I genotyping is more cost-effective than tetra markers.

CONCLUSIONS

CEL-I genotyping performed better in terms of genotyping accuracy and costs than tetra markers. The method is highly useful for validating SNPs in small to medium size germplasm panels.